A Source Transfer Domain Decomposition Method for Helmholtz Equations in Unbounded Domain

Abstract

We propose and study a domain decomposition method for solving the truncated perfectly matched layer (PML) approximation in bounded domain of Helmholtz scattering problems. The method is based on the decomposition of the domain into nonoverlapping layers and the idea of source transfer which transfers the sources equivalently layer by layer so that the solution in the final layer can be solved using a PML method defined locally outside the last two layers. The convergence of the method is proved for the case of constant wave number based on the analysis of the fundamental solution of the PML equation. The method can be used as an efficient preconditioner in the preconditioned GMRES method for solving discrete Helmholtz equations with constant and heterogeneous wave numbers. Numerical examples are included.