Abstract
An iterative solution method, in the form of a preconditioner for a Krylov subspace method, is presented for the Helmholtz equation. The preconditioner is based on a Helmholtz-type differential operator with a complex term. A multigrid iteration is used for approximately inverting the preconditioner. The choice of multigrid components for the corresponding preconditioning matrix with a complex diagonal is validated with Fourier analysis. Multigrid analysis results are verified by numerical experiments. High wavenumber Helmholtz problems in heterogeneous media are solved indicating the performance of the preconditioner.
Original language | English |
---|---|
Pages (from-to) | 1471-1492 |
Number of pages | 22 |
Journal | SIAM Journal on Scientific Computing |
Volume | 27 |
Issue number | 4 |
DOIs | |
Publication status | Published - 2006 |
Externally published | Yes |
Keywords
- Complex multigrid preconditioner
- Fourier analysis
- Helmholtz equation
- Nonconstant high wavenumber