Developing and analyzing an explicit unconditionally stable finite element scheme for an equivalent B´erenger’s PML model *

Abstract

The original B´erenger’s perfectly matched layer (PML) was quite effective in simulating wave propagation problem in unbounded domains. But its stability is very challenging to prove. Later, some equivalent PML models were developed by B´ecache and Joly [4] and their stabilities were established. Hence studying and developing efficicent numerical methods for solving those equivalent PML models are needed and interesting. Here we propose a novel explicit unconditionally stable finite element scheme to solve an equivalent B´erenger’s PML model. Both the stability and convergence analysis are proved for the proposed scheme. Numerical results justifying the theoretical analysis are presented. We also demonstrate the effectiveness of this PML in simulating wave propagation in the free space. To our best knowledge, this is the first explicit unconditionally stable finite element scheme developed for this PML model.