Back to Search Start Over

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

Authors :
Gander, M.
Graham, I.
Spence, E.
Source :
Numerische Mathematik; Nov2015, Vol. 131 Issue 3, p567-614, 48p
Publication Year :
2015

Abstract

There has been much recent research on preconditioning discretisations of the Helmholtz operator $$\Delta + k^2 $$ (subject to suitable boundary conditions) using a discrete version of the so-called 'shifted Laplacian' $$\Delta + (k^2+ \mathrm{i}\varepsilon )$$ for some $$\varepsilon >0$$ . This is motivated by the fact that, as $$\varepsilon $$ 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 $$\varepsilon $$ 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 $$\varepsilon $$ 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 $$\varepsilon $$ should be so that the preconditioner can easily be inverted by standard iterative methods. [ABSTRACT FROM AUTHOR]

Details

Language :
English
ISSN :
0029599X
Volume :
131
Issue :
3
Database :
Complementary Index
Journal :
Numerische Mathematik
Publication Type :
Academic Journal
Accession number :
109345420
Full Text :
https://doi.org/10.1007/s00211-015-0700-2