FEM and CIP-FEM for Helmholtz Equation with High Wave Number and PML truncation

Abstract: The Helmholtz scattering problem with high wave number is truncated by the perfectly matched layer (PML) technique and then discretized by the linear continuous interior penalty finite element method (CIP-FEM). It is proved that the truncated PML problem satisfies the inf--sup condition with inf--sup constant of order $O(k^{-1})$. Stability and convergence of the truncated PML problem are discussed. In particular, the convergence rate is twice of the previous result. The preasymptotic error estimates in the energy norm of the linear CIP-FEM as well as FEM are proved to be $C_1kh+C_2k^3h^2$ under the mesh condition that $k^3h^2$ is sufficiently small. Numerical tests are provided to illustrate the preasymptotic error estimates and show that the penalty parameter in the CIP-FEM may be tuned to reduce greatly the pollution error.