Application of the Complex Envelope Vectorization to a Boundary Element Formulation

Author(s):  
Oliviero Giannini ◽  
Aldo Sestieri

The complex envelope vectorization (CEV) is a recent method that has been successfully applied to structural and internal acoustic problems. Unlike other methods proposed in the last two decades to solve high frequency problems, CEV is not an energy method, although it shares with all the other techniques a variable transformation of the field variable. By such transformation involving a Hilbert transform, CEV allows the representation of a fast oscillating signal through a set of low oscillating signals. Thanks to such transformation it is possible to solve a high frequency dynamic problem at a computational cost that is lower than that required by finite elements. In fact, by using finite elements, a high frequency problem usually implies large matrices. On the contrary the CEV formulation is obtained by solving a set of linear problems of highly reduced dimensions. Although it was proved that CEV is in general a successful procedure, it was shown that it is particularly appropriate when the modes of the system have a negligible role on the solution. Moreover, the numerical advantage of the CEV formulation is much more pronounced when full matrices are used. Thus, for the first time it is applied to a boundary element formulation (BEM). Both external and internal acoustic fields of increasing complexity are considered: the internal and external field generated by a pulsating sphere; the external field of a forced box, where the velocity field is determined by finite elements; a set of 4 plates that form an open cavity. The results are compared with those obtained by a BEM procedure (SYSNOISE), highlighting the good quality of the proposed approach. An estimate of the computational advantage is also provided. Finally it is worthwhile to point out that the reduction of the BE matrices allows for an in-core solution even for large problems.

1999 ◽  
Vol 223 (2) ◽  
pp. 273-296 ◽  
Author(s):  
M. Viktorovitch ◽  
F. Thouverez ◽  
L. Jezequel

Author(s):  
A. Gibbs ◽  
D. P. Hewett ◽  
D. Huybrechs ◽  
E. Parolin

Abstract We present a hybrid numerical-asymptotic (HNA) boundary element method (BEM) for high frequency scattering by two-dimensional screens and apertures, whose computational cost to achieve any prescribed accuracy remains bounded with increasing frequency. Our method is a collocation implementation of the high order hp HNA approximation space of Hewett et al. (IMA J Numer Anal 35:1698–1728, 2015), where a Galerkin implementation was studied. An advantage of the current collocation scheme is that the one-dimensional highly oscillatory singular integrals appearing in the BEM matrix entries are significantly easier to evaluate than the two-dimensional integrals appearing in the Galerkin case, which leads to much faster computation times. Here we compute the required integrals at frequency-independent cost using the numerical method of steepest descent, which involves complex contour deformation. The change from Galerkin to collocation is nontrivial because naive collocation implementations based on square linear systems suffer from severe numerical instabilities associated with the numerical redundancy of the HNA basis, which produces highly ill-conditioned BEM matrices. In this paper we show how these instabilities can be removed by oversampling, and solving the resulting overdetermined collocation system in a weighted least-squares sense using a truncated singular value decomposition. On the basis of our numerical experiments, the amount of oversampling required to stabilise the method is modest (around 25% typically suffices), and independent of frequency. As an application of our method we present numerical results for high frequency scattering by prefractal approximations to the middle-third Cantor set.


Sign in / Sign up

Export Citation Format

Share Document