An Auxiliary Differential Equation Formulation for the Complex-Frequency Shifted PML

An efficient auxiliary-differential equation (ADE) form of the complex frequency shifted perfectly matched layer (CPML) absorbing media derived from a stretched coordinate PML formulation is presented. It is shown that a unit step response of the ADE-CPML equations leads to a discrete form that is identical to Roden's convolutional PML method for FDTD implementations. The derivation of discrete difference operators for the ADE-CPML equations for FDTD is also presented. The ADE-CPML method is also extended in a compact form to a multiple-pole PML formulation. The advantage of the ADE-CPML method is that it provides a more flexible representation that can be extended to higher-order methods. In this paper, it is applied to the discontinuous Galerkin finite element time-domain (DGFETD) method. It is demonstrated that the ADE-CPML maintains the exponential convergence of the DGFETD method.

[1]  Jian-Ming Jin,et al.  On the development of a higher-order PML , 2005, SBMO/IEEE MTT-S International Conference on Microwave and Optoelectronics, 2005..

[2]  S. Gedney,et al.  Scaled CFS-PML: it is more robust, more accurate, more efficient, and simple to implement. Why aren't you using it? , 2005, 2005 IEEE Antennas and Propagation Society International Symposium.

[3]  Stephen D. Gedney,et al.  Convolution PML (CPML): An efficient FDTD implementation of the CFS–PML for arbitrary media , 2000 .

[4]  S. Gedney,et al.  The Discontinuous Galerkin Finite Element Time Domain Method (DGFETD): A High Order, Globally-Explicit Method for Parallel Computation , 2007, 2007 IEEE International Symposium on Electromagnetic Compatibility.

[5]  Qing Huo Liu,et al.  Three‐dimensional unstructured‐grid discontinuous Galerkin method for Maxwell's equations with well‐posed perfectly matched layer , 2005 .

[6]  J. P. Webb.  Hierarchal vector basis functions of arbitrary order for triangular and tetrahedral finite elements , 1999 .

[7]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .

[8]  S. Gedney.  An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices , 1996 .

[9]  Qing Huo Liu,et al.  The PSTD algorithm: A time-domain method requiring only two cells per wavelength , 1997 .

[10]  An Ping Zhao,et al.  Application of the material-independent PML absorbers to the FDTD analysis of electromagnetic waves in nonlinear media , 1998 .

[11]  J. Hesthaven,et al.  High-Order Accurate Methods for Time-domain Electromagnetics , 2004 .

[12]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[13]  Raj Mittra,et al.  Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers , 1996 .

[14]  Jin-Fa Lee,et al.  A perfectly matched anisotropic absorber for use as an absorbing boundary condition , 1995 .

[15]  M.W. Chevalier,et al.  A PML using a convolutional curl operator and a numerical reflection coefficient for general linear media , 2004, IEEE Transactions on Antennas and Propagation.

[16]  Stephen D. Gedney,et al.  An Anisotropic PML Absorbing Media for the FDTD Simulation of Fields in Lossy and Dispersive Media , 1996 .

[17]  Stephen D. Gedney,et al.  The Discontinuous Galerkin Finite-Element Time-Domain Method Solution of Maxwell ’ s Equations , 2009 .

[18]  Chang-Hong Liang,et al.  A new implementation of CFS‐PML for ADI‐FDTD method , 2006 .

[19]  J.-P. Wrenger,et al.  Numerical reflection from FDTD-PMLs: a comparison of the split PML with the unsplit and CFS PMLs , 2002 .

[20]  Davi Correia,et al.  Performance of regular PML, CFS‐PML, and second‐order PML for waveguide problems , 2006 .

[21]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .