AbstractAn efficient unsplit perfectly matched layer for numerical simulation of electromagnetic waves in unbounded domains is derived via a complex change of variables. In order to surround a Cartesian grid with the PML, the time-dependent PML requires only one (scalar) auxiliary variable in two space dimensions and six (scalar) auxiliary variables in three space dimensions. It is therefore cheap and straightforward to implement. We use Fourier and energy methods to prove the stability of the PML. We extend the stability result to a semi-discrete PML approximated by central finite differences of arbitrary order of accuracy and to a fully discrete problem for the 'Leap-Frog' schemes. This makes precise the usefulness of the derived PML model for longtime simulations. Numerical experiments are presented, illustrating the accuracy and stability of the PML. © 2013 IMACS.
CitationDuru K, Kreiss G (2014) Efficient and stable perfectly matched layer for CEM. Applied Numerical Mathematics 76: 34–47. Available: http://dx.doi.org/10.1016/j.apnum.2013.09.005.
SponsorsThis project was completed during the author's postdoctoral program at the Geophysics Department, Stanford University, California. This work was supported by King Abdullah University of Science and Technology (KAUST) through a joint KAUST Academic Excellence Alliance (AEA) grant with Stanford.
JournalApplied Numerical Mathematics