Perfectly matched layer in three dimensions for the time-domain finite element method applied to radiation problems

A new perfectly matched layer (PML) formulation for the time-domain finite element method is described and tested for problems in three-dimensional space. An integral part of this PML algorithm is a novel time integration scheme based on Galerkin's method with a piecewise linear temporal expansion of the electric field. Our time integration procedure is constructed by forming a linear combination of exact and trapezoidal integration applied to the temporal weak form, which reduces to the well-known Newmark scheme in the case where the PML is not present. This technique allows for a consistent time discretization of the electric field vector wave equation augmented with the PML, which includes terms that involve the electric field itself, time derivatives thereof and exponential functions temporally convolved with the electric field. The auxiliary variables that are associated with such convolutions are efficiently updated in an explicit manner by matrix-vector products. The new PML formulation is successfully tested on canonical problems. The PML can very efficiently absorb propagating electromagnetic waves, which is demonstrated by reflections coefficients as low as -100 dB. The PML formulation is capable of giving highly accurate results for very broad frequency bands which is confirmed by tests that show an error smaller than 2% over more than six octaves in frequency. Finally, we apply the PML formulation to a bow-tie antenna and comparisons with measurements show good agreement.

[1]  Jian-Ming Jin,et al.  Time-domain finite-element simulation of three-dimensional scattering and radiation problems using perfectly matched layers , 2003 .

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

[3]  P. Monk,et al.  Gauss Point Mass Lumping Schemes for Maxwell's Equations , 1998 .

[4]  Nathan M. Newmark,et al.  A Method of Computation for Structural Dynamics , 1959 .

[5]  Jian-Ming Jin,et al.  Perfectly matched layer for the time domain finite element method , 2004 .

[6]  V. Mathis.  An anisotropic perfectly matched layer-absorbing medium in finite element time domain method for Maxwell's equations , 1997, IEEE Antennas and Propagation Society International Symposium 1997. Digest.

[7]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

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

[9]  Jianming Jin,et al.  A novel efficient algorithm for scattering from a complex BOR using mixed finite elements and cylindrical PML , 1999 .

[10]  Jin-Fa Lee,et al.  A comparison of anisotropic PML to Berenger's PML and its application to the finite-element method for EM scattering , 1997 .

[11]  R. Cools,et al.  Monomial cubature rules since “Stroud”: a compilation , 1993 .

[12]  O. Picon,et al.  A finite element method based on Whitney forms to solve Maxwell equations in the time domain , 1995 .

[13]  S. Eisenstat.  Efficient Implementation of a Class of Preconditioned Conjugate Gradient Methods , 1981 .

[14]  Boualem Boashash,et al.  Estimating and interpreting the instantaneous frequency of a signal. II. A/lgorithms and applications , 1992, Proc. IEEE.

[15]  Raj Mittra,et al.  Investigation of nonplanar perfectly matched absorbers for finite-element mesh truncation , 1997 .

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

[17]  Gary Cohen.  Higher-Order Numerical Methods for Transient Wave Equations , 2001 .

[18]  Jian-Ming Jin,et al.  The Finite Element Method in Electromagnetics , 1993 .

[19]  R. Cools.  Monomial cubature rules since “Stroud”: a compilation—part 2 , 1999 .

[20]  Tatsuo Itoh,et al.  An unconditionally stable extended (USE) finite-element time-domain solution of active nonlinear microwave circuits using perfectly matched layers , 2001, IMS 2001.

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

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

[23]  Andreas C. Cangellaris,et al.  GT-PML: generalized theory of perfectly matched layers and its application to the reflectionless truncation of finite-difference time-domain grids , 1996, 1996 IEEE MTT-S International Microwave Symposium Digest.

[24]  Jian-Ming Jin,et al.  An effective algorithm for implementing perfectly matched layers in time-domain finite-element simulation of open-region EM problems , 2002 .

[25]  Carey M. Rappaport.  Interpreting and improving the PML absorbing boundary condition using anisotropic lossy mapping of space , 1996 .

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

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

[28]  Jian-Ming Jin,et al.  Finite element analysis of complex axisymmetric radiating structures , 1999, IEEE Antennas and Propagation Society International Symposium. 1999 Digest. Held in conjunction with: USNC/URSI National Radio Science Meeting (Cat. No.99CH37010).