Pseudospectrum Analysis on the Computation of Resonances in Open Cavities with Perfectly Matched Layers

Abstract eng:
During the last decade, several authors have addressed that the PML technique [1] can be used not only for the computation of the near-field in time-dependent and time-harmonic scattering problems, but also to compute numerically the resonances in open cavities (see for instance [2]). Despite such complex resonances are not natural eigen-frequencies of the physical system, the numerical determination of this kind of eigenvalues provides information about the model, what can be used in further applications. For instance, the computation of leaky modes in an open waveguide is used as the spectral basis for a modal decomposition technique [3]. In other contexts, such computations allow us to obtain the localization of maximums and minimums for the frequency response function of a mechanical model [2] stated in an unbounded domain. In the framework of the PML technique, to ensure that the complex eigenvalues are independent of the particular choices of the so-called PML absorption coefficient, the theoretical background is related with the complex scaling technique developed in the seventies by Aguilar and Combes. However, the theoretical results do not guarantee any convergence result on the PML error, which comes from the spurious reflections originated by the truncation of the PML layer to a finite thickness. Moreover, the numerical computation of these complex resonances have been revealed highly unstable (see, for instance [2]) for standard finite element methods (FEM). The numerical results depend strongly not only on the thickness of the PML layer but also on the position of its inner PML boundary with respect to physical domain of the problem. All these numerical difficulties can be interpreted by using the concept of pseudospetrum [4] associated to the PML problem and its nonnormal behavior. In addition, a comparison between a standard FEM and a Chebyshev collocation discretization is done to illustrate the computational cost required to achieve accurate results. References [1] J. P. Bérenger, A perfectly matched layer for the absorption of electromagnetic waves, Journal of Computational Physics, 114, 185-200, 1994. [2] S. Hein, T. Hohage, W. Koch, J. Schöberl, Acoustic resonances in a high-lift configuration, Journal of Fluid Mechanics, 582, 179-202, 2007. [3] A.-S. Bonnet-Ben Dhia, B. Goursaud, C. Hazard, and A. Prieto. Finite element computation of leaky modes in stratified waveguides. In Ultrasonic Wave Propagation in Non Homogeneous Media, vol. 128, 73–86, Springer, 2009. [4] L. N. Trefethen, Spectra and pseudospectra, Princeton University Press, 2005.

• Export as: BibTeX | MARC | MARCXML | DC | EndNote | NLM | RefWorks
• Add to your basket