On the Accuracy and Stability of Several Widely Used FDTD Approaches for Modeling Lorentz Dielectrics

A rigorous and comparative study on the approximation accuracy and stability limits of several widely used finite-difference time-domain (FDTD) approaches, namely the auxiliary differential equation (ADE) approach, the bilinear transform (BT) approach, the Z-transform approach (ZT) and the piecewise linear recursive convolution (PLRC) approach, for modeling dispersive Lorentz dielectrics is presented following the given updating equations between the electric flux density and the electric field intensity. We find the ZT approach with modified material parameters is much more accurate than the original ZT approach and the other three approaches for modeling Lorentz dielectrics. The stability limits of the ADE, ZT and PLRC approaches in simulating Lorentz dielectrics are also shown to be a bit more stringent than that of BT approach which preserves the Courant stability limit as previously reported.

[1]  Angel Vegas,et al.  Analyzing the stability of the FDTD technique by combining the von Neumann method with the Routh-Hurwitz criterion , 2001 .

[2]  Andre Knoesen,et al.  DISPERSIVE MODELS FOR THE FINITE-DIFFERENCE TIME-DOMAIN METHOD : DESIGN, ANALYSIS, AND IMPLEMENTATION , 1994 .

[3]  Ichiro Fukai,et al.  A treatment by the FD‐TD method of the dispersive characteristics associated with electronic polarization , 1990 .

[4]  Dennis M. Sullivan,et al.  Frequency-dependent FDTD methods using Z transforms , 1992 .

[5]  M. A. Stuchly,et al.  Simple treatment of multi-term dispersion in FDTD , 1997 .

[6]  Juan Chen,et al.  Error between unconditionally stable FDTD methods and conventional FDTD method , 2006 .

[7]  Jeffrey L. Young,et al.  A summary and systematic analysis of FDTD algorithms for linearly dispersive media , 2001 .

[8]  R. J. Luebbers,et al.  Piecewise linear recursive convolution for dispersive media using FDTD , 1996 .

[9]  Joe LoVetri,et al.  A comparison of numerical techniques for modeling electromagnetic dispersive media , 1995 .

[10]  A. Grande,et al.  On the Equivalence of Several FDTD Formulations for Modeling Electromagnetic Wave Propagation in Double-Negative Metamaterials , 2007, IEEE Antennas and Wireless Propagation Letters.

[11]  Jianguo Wang,et al.  PEC condition implementation for the ADI‐FDTD method , 2007 .

[12]  A Taflove,et al.  Direct time integration of Maxwell's equations in linear dispersive media with absorption for scattering and propagation of femtosecond electromagnetic pulses. , 1991, Optics letters.

[13]  P. Bevelacqua,et al.  Finite-difference and pseudospectral time-domain methods applied to backward-wave metamaterials , 2004, IEEE Transactions on Antennas and Propagation.

[14]  Nicholas Chako,et al.  Wave propagation and group velocity , 1960 .

[15]  Y. Takayama,et al.  Reinterpretation of the auxiliary differential equation method for FDTD , 2002, IEEE Microwave and Wireless Components Letters.

[16]  Dennis M. Sullivan.  Digital filtering techniques for use with the FDTD method , 1999 .

[17]  An Ping Zhao.  A novel implementation for two‐dimensional unconditionally stable FDTD method , 2003 .