A multigrid-based preconditioned Krylov subspace method for the Helmholtz equation with PML (Q2275513)

The authors consider the Helmholtz equation where the absorbing boundary conditions are often used to truncate the infinite domain into a finite domain; here they use the perfectly matched layer (PML) to truncate the domain. PML is an absorbing boundary condition with the property of generating almost no reflection at the PML interface between the interior medium and the artificial absorbing medium. To discretize Helmholtz-PML equation, the authors use a rotated 9-point difference scheme which uses fewer grids per weblength than the conventional 5-point difference scheme while maintaining comparable accuracy. Spectral analysis of the linear system is carried out and a new matrix-based interpolation operator is proposed for the multigrid method, which is used to approximately invert the preconditioner. Numerical experiments are presented to illustrate the efficiency of the preconditioned Bi-CGSTAB method with the multigrid based on the new interpolation operator.