A staggered fourth-order accurate explicit finite difference scheme for the time-domain Maxwell's equations

Abstract We consider a model explicit fourth-order staggered finite-difference method for the hyperbolic Maxwell's equations. Appropriate fourth-order accurate extrapolation and one-sided difference operators are derived in order to complete the scheme near metal boundaries and dielectric interfaces. An eigenvalue analysis of the overall scheme provides a necessary, but not sufficient, stability condition and indicates long-time stability. Numerical results verify both the stability analysis, and the scheme's fourth-order convergence rate over complex domains that include dielectric interfaces and perfectly conducting surfaces. For a fixed error level, we find the fourth-order scheme is computationally cheaper in comparison to the Yee scheme by more than an order of magnitude. Some open problems encountered in the application of such high-order schemes are also discussed.

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

[2]  David W. Zingg,et al.  High-Accuracy Finite-Difference Schemes for Linear Wave Propagation , 1996, SIAM J. Sci. Comput..

[3]  J. S. Shang,et al.  High-Order Compact-Difference Schemes for Time-Dependent Maxwell Equations , 1999 .

[4]  Arje Nachman.  A Brief Perspective on Computational Electromagnetics , 1996 .

[5]  Peter G. Petropoulos,et al.  Phase error control for FD-TD methods of second and fourth order accuracy , 1994 .

[6]  Kenneth M. Mitzner,et al.  How Scattering Increases as an Edge Is Blunted: The Case of an Electric Field Parallel to the Edge , 1990 .

[7]  T. Senior,et al.  Electromagnetic and Acoustic Scattering by Simple Shapes , 1969 .

[8]  Eli Turkel,et al.  On the construction of a high order difference scheme for complex domains in a Cartesian grid , 2000 .

[9]  Yen Liu,et al.  Fourier Analysis of Numerical Algorithms for the Maxwell Equations , 1993 .

[10]  R. LeVeque,et al.  The immersed interface method for acoustic wave equations with discontinuous coefficients , 1997 .

[11]  Eli Turkel,et al.  Fourth order compact implicit method for the Maxwell equations with discontinuous coefficients , 2000 .

[12]  Tobin A. Driscoll,et al.  Staggered Time Integrators for Wave Equations , 2000, SIAM J. Numer. Anal..

[13]  J. Hesthaven,et al.  Staircase-free finite-difference time-domain formulation for general materials in complex geometries , 2001 .

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

[15]  J. Hesthaven,et al.  Convergent Cartesian Grid Methods for Maxwell's Equations in Complex Geometries , 2001 .

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

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