Efficient and stable perfectly matched layer for CEM

Kenneth Duru, Gunilla Kreiss

Research output: Contribution to journalArticlepeer-review

9 Scopus citations

Abstract

An 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.
Original languageEnglish (US)
Pages (from-to)34-47
Number of pages14
JournalApplied Numerical Mathematics
Volume76
DOIs
StatePublished - Feb 2014
Externally publishedYes

Bibliographical note

KAUST Repository Item: Exported on 2020-10-01
Acknowledgements: This 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.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.

Fingerprint

Dive into the research topics of 'Efficient and stable perfectly matched layer for CEM'. Together they form a unique fingerprint.

Cite this