Applying GMRES to the Helmholtz equation with shifted Laplacian preconditioning: what is the largest shift for which wavenumber-independent convergence is guaranteed?

• Published in: Numerische Mathematik Vol. 131; no. 3; pp. 567 - 614
• Main Authors: Gander, M. J., Graham, I. G., Spence, E. A.
• Format: Journal Article
• Language: English
• Published: Berlin/Heidelberg Springer Berlin Heidelberg 01.11.2015
• Subjects: Mathematical and Computational Engineering; Mathematical and Computational Physics; Mathematical Methods in Physics; Mathematics; Mathematics and Statistics; Numerical Analysis; Numerical and Computational Physics; Simulation; Theoretical; 78A45; 35J05; 65F10; 65F08; 65N30
• ISSN: 0029-599X; 0945-3245
• DOI: 10.1007/s00211-015-0700-2

There has been much recent research on preconditioning discretisations of the Helmholtz operator Δ + k 2 (subject to suitable boundary conditions) using a discrete version of the so-called “shifted Laplacian” Δ + ( k 2 + i ε ) for some ε > 0 . This is motivated by the fact that, as ε increases, the shifted problem becomes easier to solve iteratively. Despite many numerical investigations, there has been no rigorous analysis of how to chose the shift. In this paper, we focus on the question of how large ε can be so that the shifted problem provides a preconditioner that leads to k -independent convergence of GMRES, and our main result is a sufficient condition on ε for this property to hold. This result holds for finite element discretisations of both the interior impedance problem and the sound-soft scattering problem (with the radiation condition in the latter problem imposed as a far-field impedance boundary condition). Note that we do not address the important question of how large ε should be so that the preconditioner can easily be inverted by standard iterative methods.