Development and analysis of both finite element and fourth-order in space finite difference methods for an equivalent Berenger's PML model

Abstract

This paper deals with an equivalent Berenger's Perfectly Matched Layer (PML) model. We first develop a finite element scheme using edge elements to solve this model. We prove a discrete stability of this method, which inherits the stability obtained in the continuous case. Then we propose a fourth-order in space finite difference scheme for solving this PML model. Numerical stability similar to the continuous stability and the optimal error estimate are established for the difference scheme. Here only second order time discretizations are considered for both schemes. Finally, numerical results are presented to justify our analysis and demonstrate the effectiveness of this PML model for absorbing impinging waves.