Perfectly matched layer for the time domain finite element method

A new perfectly matched layer (PML) formulation for the time domain finite element method is described and tested for Maxwell's equations. In particular, we focus on the time integration scheme which is based on Galerkin's method with a temporally piecewise linear expansion of the electric field. The time stepping scheme 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 without PML. Extensive numerical tests on scattering from infinitely long metal cylinders in two dimensions show good accuracy and no signs of instabilities. For a circular cylinder, the proposed scheme indicates the expected second order convergence toward the analytic solution and gives less than 2% root-mean-square error in the bistatic radar cross section (RCS) for resolutions with more than 10 points per wavelength. An ogival cylinder, which has sharp corners supporting field singularities, shows similar accuracy in the monostatic RCS.

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

[2]  Analysis of exponential time-differencing for FDTD in lossy dielectrics , 1997 .

[3]  Fernando L. Teixeira,et al.  On causality and dynamic stability of perfectly matched layers for FDTD simulations , 1999 .

[4]  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).

[5]  W. Chew,et al.  Systematic derivation of anisotropic PML absorbing media in cylindrical and spherical coordinates , 1997 .

[6]  Jian-Ming Jin,et al.  Simple moment method program for computing scattering from complex cylindrical obstacles , 1989 .

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

[8]  A. Bondeson,et al.  Tunable Integration Scheme for the Finite-Element Method , 1991 .

[9]  Jian-Ming Jin,et al.  Complex coordinate stretching as a generalized absorbing boundary condition , 1997 .

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

[11]  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.

[12]  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, IMS 1996.

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

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

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

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

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

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

[19]  Jin-Fa Lee,et al.  Time-domain finite-element methods , 1997 .

[20]  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.

[21]  Accurate extrapolation to zero cell size by Padé approximation , 2003 .

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

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

[24]  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 .

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

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

[27]  C. Balanis.  Advanced Engineering Electromagnetics , 1989 .

[28]  Thomas Rylander,et al.  Application of stable FEM-FDTD hybrid to scattering problems , 2002 .

[29]  Patrick Joly,et al.  Stability of perfectly matched layers, group velocities and anisotropic waves , 2003 .

[30]  E. Turkel,et al.  Absorbing PML boundary layers for wave-like equations , 1998 .

[31]  Boualem Boashash,et al.  Estimating and interpreting the instantaneous frequency of a signal. I. Fundamentals , 1992, Proc. IEEE.

[32]  J. Nédélec.  Mixed finite elements in ℝ3 , 1980 .

[33]  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 .

[34]  Weng Cho Chew,et al.  PML-FDTD in cylindrical and spherical grids , 1997 .