Toward the construction of a fourth-order difference scheme for transient EM wave simulation: staggered grid approach

A compact central-difference approximation in conjunction with the Yee (1966) grid is used to compute the spatial derivatives in Maxwell's equations. To advance the semi-discrete equations, the four-stage Runge-Kutta (RK) integrator is invoked. This combination of spatial and temporal differencing leads to a scheme that is fourth-order accurate, conditionally stable, and highly efficient. Moreover, the use of compact differencing allows one to apply the compact operator in the vicinity of a perfect conductor-an attribute not found in other higher order methods. Results are provided that quantify the spectral properties of the method. Simulations are conducted on problem spaces that span one and three dimensions and whose domains are of the open and closed type. Results from these simulations are compared with exact closed-form solutions; the agreement between these results is consistent with numerical analysis.

[1]  D. Fyfe,et al.  Economical evaluation of Runge-Kutta formulae , 1966 .

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

[3]  J. Williamson Low-storage Runge-Kutta schemes , 1980 .

[4]  D. Gaitonde,et al.  Practical aspects of high-order accurate finite-volume schemes for electromagnetics , 1997 .

[5]  D. Gaitonde,et al.  High-order finite-volume schemes in wave propagation phenomena , 1996 .

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

[7]  A. Cangellaris,et al.  Numerical stability and numerical dispersion of a compact 2-D/FDTD method used for the dispersion analysis of waveguides , 1993, IEEE Microwave and Guided Wave Letters.

[8]  D. Gottlieb,et al.  The Stability of Numerical Boundary Treatments for Compact High-Order Finite-Difference Schemes , 1993 .

[9]  B. Gustafsson The convergence rate for difference approximations to mixed initial boundary value problems , 1975 .

[10]  J. S. Shang,et al.  Characteristic-based algorithms for solving the Maxwell equations in the time domain , 1995 .

[11]  James R. Wait Introduction to antennas & propagation , 1986 .

[12]  G. Mur Absorbing Boundary Conditions for the Finite-Difference Approximation of the Time-Domain Electromagnetic-Field Equations , 1981, IEEE Transactions on Electromagnetic Compatibility.

[13]  Om P. Gandhi,et al.  Use of the Finite-Difference Time-Domain Method in Calculating EM Absorption in Human Tissues , 1987, IEEE Transactions on Biomedical Engineering.

[14]  A. Taflove,et al.  Radar Cross Section of General Three-Dimensional Scatterers , 1983, IEEE Transactions on Electromagnetic Compatibility.

[15]  William H. Press,et al.  Numerical recipes , 1990 .

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

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

[18]  Andrew F. Peterson,et al.  Relative accuracy of several finite-difference time-domain methods in two and three dimensions , 1993 .