Error analysis of fully discrete mixed finite element schemes for 3-D Maxwell’s equations in dispersive media

Abstract We consider the time dependent Maxwell’s equations in dispersive media in a bounded three-dimensional (3-D) domain. Fully discrete mixed finite element methods are developed for three most popular dispersive media models: i.e., the isotropic cold plasma, the one-pole Debye medium and the two-pole Lorentz medium. Optimal error estimates are proved for those three dispersive media under the assumption that the solutions are smooth enough. Numerical results justifying our analysis are provided. We believe that this is the first error analysis carried out for fully discrete mixed finite element methods for solving 3-D Maxwell’s equations in dispersive media.

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

[2]  Jun Zou,et al.  Fully discrete finite element approaches for time-dependent Maxwell's equations , 1999, Numerische Mathematik.

[3]  Bruce Archambeault,et al.  The Finite-Difference Time-Domain Method , 1998 .

[4]  Niel K. Madsen,et al.  A mixed finite element formulation for Maxwell's equations in the time domain , 1990 .

[5]  P. Neittaanmäki,et al.  Mathematical and Numerical Aspects of Wave Propagation WAVES 2003 , 2003 .

[6]  Jichun Li,et al.  Analysis of a time-domain finite element method for 3-D Maxwell’s equations in dispersive media , 2006 .

[7]  Jichun Li.  Error analysis of finite element methods for 3-D Maxwell's equations in dispersive media , 2006 .

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

[9]  J. Hesthaven,et al.  Nodal high-order methods on unstructured grids , 2002 .

[10]  R. Luebbers,et al.  The Finite Difference Time Domain Method for Electromagnetics , 1993 .

[11]  Cynthia Furse,et al.  Currents induced in the human body for exposure to ultrawideband electromagnetic pulses , 1997 .

[12]  L. Demkowicz,et al.  hp-adaptive finite elements in electromagnetics , 1999 .

[13]  Chi-Wang Shu,et al.  Locally divergence-free discontinuous Galerkin methods for the Maxwell equations , 2004, Journal of Computational Physics.

[14]  Steven A. Cummer,et al.  An analysis of new and existing FDTD methods for isotropic cold plasma and a method for improving their accuracy , 1997 .

[15]  Peter Monk.  A Comparison of Three Mixed Methods for the Time-Dependent Maxwell's Equations , 1992, SIAM J. Sci. Comput..

[16]  Pingwen Zhang,et al.  Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions , 2004 .

[17]  Qiang Chen,et al.  An FDTD formulation for dispersive media using a current density , 1998 .

[18]  Jian-Ming Jin,et al.  Time-domain finite element modeling of dispersive media , 2001, IEEE Antennas and Propagation Society International Symposium. 2001 Digest. Held in conjunction with: USNC/URSI National Radio Science Meeting (Cat. No.01CH37229).

[19]  Daniel A. White,et al.  A high order mixed vector finite element method for solving the time dependent Maxwell equations on unstructured grids , 2005 .

[20]  D. Boffi,et al.  Computational Models of Electromagnetic Resonators: Analysis of Edge Element Approximation , 1999 .

[21]  V. Girault,et al.  Vector potentials in three-dimensional non-smooth domains , 1998 .

[22]  S. S. Stuchly,et al.  Propagation of transients in dispersive dielectric media , 1991 .

[23]  L. Demkowicz,et al.  Maxwell eigenvalues and discrete compactness in two dimensions , 2000 .

[24]  L. Demkowicz,et al.  De Rham diagram for hp finite element spaces , 2000 .

[25]  Propagation of ultra-wide-band electromagnetic pulses through dispersive media , 1995 .

[26]  Ilaria Perugia,et al.  Energy norm a posteriori error estimation for mixed discontinuous Galerkin approximations of the Maxwell operator , 2005 .

[27]  J. Young.  A higher order FDTD method for EM propagation in a collisionless cold plasma , 1996 .

[28]  Dennis M. Sullivan,et al.  A frequency-dependent FDTD method for biological applications , 1992 .

[29]  D. Sullivan,et al.  On the dispersion errors related to (FD)/sup 2/TD type schemes , 1995 .

[30]  Joseph E. Pasciak,et al.  Analysis of a Multigrid Algorithm for Time Harmonic Maxwell Equations , 2004, SIAM J. Numer. Anal..

[31]  A. Quarteroni,et al.  Numerical Approximation of Partial Differential Equations , 2008 .

[32]  R. Hiptmair.  Finite elements in computational electromagnetism , 2002, Acta Numerica.

[33]  Ralf Hiptmair,et al.  Mixed finite element approximation of eddy current problems , 2004 .

[34]  Dongwoo Sheen,et al.  A NONCONFORMING MIXED FINITE ELEMENT METHOD FOR MAXWELL'S EQUATIONS , 2000 .

[35]  P. Raviart,et al.  Study of an implicit scheme for integrating Maxwell's equations , 1980 .

[36]  Peter Monk,et al.  A mixed method for approximating Maxwell's equations , 1991 .

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

[38]  Ralf Hiptmair,et al.  Boundary Element Methods for Maxwell Transmission Problems in Lipschitz Domains , 2003, Numerische Mathematik.

[39]  Qing Huo Liu,et al.  Simulations of GPR in dispersive media using a frequency-dependent PSTD algorithm , 1999, IEEE Trans. Geosci. Remote. Sens..

[40]  Jian-Ming Jin,et al.  Frequency-domain and time-domain finite-element solution of Maxwell's equations using spectral Lanczos decomposition method , 1999 .

[41]  J. Tinsley Oden,et al.  Finite Elements, Mathematical Aspects. , 1986 .

[42]  Jianming Jin,et al.  Efficient time-domain and frequency-domain finite-element solution of Maxwell's equations using spectral Lanczos decomposition method , 1998 .