A two-level shifted Laplace preconditioner for Helmholtz problems: Field-of-values analysis and wavenumber-independent convergence

arXiv2006.08750MaRDI QIDQ6342995FDOQ6342995

Luis García Ramos, Reinhard Nabben

Publication date: 15 June 2020

Abstract: One of the main tools for solving linear systems arising from the discretization of the Helmholtz equation is the shifted Laplace preconditioner, which results from the discretization of a perturbed Helmholtz problem $-Deltau-(k^2+ivarepsilon)u=f$ where $0eqvarepsiloninmathbbR$ is an absorption parameter. In this work we revisit the idea of combining the shifted Laplace preconditioner with two-level deflation and apply it to Helmholtz problems discretized with linear finite elements. We use the convergence theory of GMRES based on the field of values to prove that GMRES applied to the two-level preconditioned system with a shift parameter $varepsilonsim{k}^2$ converges in a number of iterations independent of the wavenumber $k$,provided that the coarse mesh size $H$ satisfies a condition of the form $H{k}^{2}leqC$ for some constant $C$ depending on the domain but independent of the wavenumber $k$. This behaviour is sharply different to the standalone shifted Laplacian, for which wavenumber-independent GMRES convergence has been established only under the condition that $varepsilonsimk$ by [M.J. Gander, I.G. Graham and E.A. Spence, Numer. Math., 131 (2015), 567-614]. Finally, we present numerical evidence that wavenumber-independent convergence of GMRES also holds for pollution-free meshes, where the coarse mesh size satisfies $H{k}^{3/2}leqC$, and inexact coarse grid solves.

Mathematics Subject Classification ID

Preconditioners for iterative methods (65F08) Iterative numerical methods for linear systems (65F10) Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30)

This page was built for publication: A two-level shifted Laplace preconditioner for Helmholtz problems: Field-of-values analysis and wavenumber-independent convergence

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6342995)