Low-dispersion algorithms based on the higher order (2,4) FDTD method

This paper discusses the enhancement of numerical dispersion characteristics in the context of the finite-difference time-domain method based on a (2,4) computational stencil. Rather than implementing the conventional approach-based on Taylor analysis-for the determination of the finite-difference operators, two alternative procedures that result in numerical schemes with diverse wide-band behavior are proposed. First, an algorithm that performs better than the standard counterpart over all frequencies is constructed by requiring the mutual cancellation of terms with equal order in the corresponding dispersion relation. In addition, a second method is derived, which is founded on the separate optimization of the spatial and temporal derivatives. In this case, analysis proves that significant error compensation is accomplished around a specific design frequency, while reduced errors are obtained for higher frequencies, thus enabling the reliable execution of wide-band simulations as well. The quality and efficiency of the proposed techniques, which exhibit the same computational requirements as the standard (2,4) approach, are investigated theoretically, and subsequently, validated by means of numerical experimentation.

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

[2]  Zhongxiang Shen,et al.  A generalized higher order finite-difference time-domain method and its application in guided-wave problems , 2003 .

[3]  Shlomo Ta'asan,et al.  Finite difference schemes for long-time integration , 1994 .

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

[5]  W. Chew,et al.  A time-domain method with isotropic dispersion and increased stability on an overlapped lattice , 2002 .

[6]  Melinda Piket-May,et al.  A modified FDTD (2, 4) scheme for modeling electrically large structures with high-phase accuracy , 1997 .

[7]  John B. Schneider,et al.  Comparison of the dispersion properties of several low-dispersion finite-difference time-domain algorithms , 2003 .

[8]  C. Balanis,et al.  Higher-order finite-difference schemes for electromagnetic radiation, scattering, and penetration .2. Applications , 2002 .

[9]  Fernando L. Teixeira,et al.  A three-dimensional angle-optimized finite-difference time-domain algorithm , 2003 .

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

[11]  Jeffrey L. Young,et al.  Toward the construction of a fourth-order difference scheme for transient EM wave simulation: staggered grid approach , 1997 .

[12]  C. Tam,et al.  Dispersion-relation-preserving finite difference schemes for computational acoustics , 1993 .

[13]  Stavros V. Georgakopoulos,et al.  Higher-Order Finite-Difference Schemes for Electromagnetic Radiation , Scattering , and Penetration , Part I : Theory , 2004 .

[14]  Nikolaos V. Kantartzis,et al.  A higher-order FDTD technique for the implementation of enhanced dispersionless perfectly matched layers combined with efficient absorbing boundary conditions , 1998 .

[15]  Fernando L. Teixeira,et al.  Dispersion-relation-preserving FDTD algorithms for large-scale three-dimensional problems , 2003 .

[16]  J. M. Sanz-Serna,et al.  High-order methods , 1994 .

[17]  J. B. Cole.  A high-accuracy realization of the Yee algorithm using non-standard finite differences , 1997 .

[18]  Amir Yefet,et al.  A staggered fourth-order accurate explicit finite difference scheme for the time-domain Maxwell's equations , 2001 .

[19]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[20]  Robert J. Lee,et al.  Reducing the phase error for finite-difference methods without increasing the order , 1998 .