English

Recent Results on Domain Decomposition Preconditioning for the High-frequency Helmholtz Equation using Absorption

Numerical Analysis 2016-06-24 v1

Abstract

In this paper we present an overview of recent progress on the development and analysis of domain decomposition preconditioners for discretised Helmholtz problems, where the preconditioner is constructed from the corresponding problem with added absorption. Our preconditioners incorporate local subproblems that can have various boundary conditions, and include the possibility of a global coarse mesh. While the rigorous analysis describes preconditioners for the Helmholtz problem with added absorption, this theory also informs the development of efficient multilevel solvers for the "pure" Helmholtz problem without absorption. For this case, 2D experiments for problems containing up to about 5050 wavelengths are presented. The experiments show iteration counts of order about O(n0.2)\mathcal{O}(n^{0.2}) and times (on a serial machine) of order about O(nα)\mathcal{O}(n^{\alpha}), { with α[1.3,1.4]\alpha \in [1.3,1.4]} for solving systems of dimension nn. This holds both in the pollution-free case corresponding to meshes with grid size O(k3/2)\mathcal{O}(k^{-3/2}) (as the wavenumber kk increases), and also for discretisations with a fixed number of grid points per wavelength, commonly used in applications. Parallelisation of the algorithms is also briefly discussed.

Keywords

Cite

@article{arxiv.1606.07172,
  title  = {Recent Results on Domain Decomposition Preconditioning for the High-frequency Helmholtz Equation using Absorption},
  author = {I. G. Graham and E. A. Spence and E. Vainikko},
  journal= {arXiv preprint arXiv:1606.07172},
  year   = {2016}
}

Comments

To appear in: Modern Solvers for Helmholtz problems edited by D. Lahaye, J. Tang and C. Vuik, Springer Geosystems Mathematics series, 2016

R2 v1 2026-06-22T14:32:17.089Z