ACCURACY OF THREE UNCONDITIONALLY-STABLE FDTD SCHEMES

作者

  • Guilin Sun Electromagnetic Compatibility Laboratory, Concordia University, 7141 Sherbrooke Street West, Montreal, Quebec, Canada H4B 1R6
  • Christopher W. Trueman Electromagnetic Compatibility Laboratory, Concordia University, 7141 Sherbrooke Street West, Montreal, Quebec, Canada H4B 1R6

关键词:

ACCURACY OF THREE UNCONDITIONALLY-STABLE FDTD SCHEMES

摘要

This paper discusses accuracy limitations due to numerical dispersion and time step size for three implicit unconditionally-stable FDTD methods: Alternate-Direction-Implicit (ADI), Crank-Nicolson (CN) and Crank-Nicolson-Douglas-Gunn (CNDG). It is shown that for a uniform mesh, the three methods have the same numerical phase velocity along the axes, but have large differences along the diagonals. The ADI method has two orders-of-magnitude larger anisotropy than that of CN and CNDG. CNDG has no anisotropy at certain Courant numbers and mesh densities. At the limit of zero spatial mesh size, the three methods have different “intrinsic temporal dispersion” for a given time step size: CN has no anisotropy; ADI has positive anisotropy and CNDG has negative anisotropy, which is much smaller than ADI. The Nyquist sampling theorem provides a fundamental upper bound on the time step size for all three methods. It is shown that for ADI and CN the practical upper bound is close to the Nyquist limit, but for CNDG is half the Nyquist limit.

##plugins.generic.usageStats.downloads##

##plugins.generic.usageStats.noStats##

参考

K.S., Yee, “Numerical solution of initial boundary

value problems involving Maxwell’s equations in

isotropic media,” IEEE Trans. Antennas and

Propagation, Vol.14, May, 1966, pp. 302-307.

ACES JOURNAL, VOL. 18, NO. 4, NOVEMBER 2003

A.,Taflove, and S.C., Hagness, Computational

electrodynamics—the finite-difference time-domain

method, 2 nd ed., Artech House, Boston, 2000.

J.Stoer, R.Bulirsch, “Introduction to numerical

analysis,” 3 rd ed, Translated by R.Bartels, W. Gautschi

and C.Witzgall, Springer-Verlag, New York, 2002.

R.Holland, “Implicit three-dimensional finite

differencing of Maxwell’s equations for computational

electromagnetism,” IEEE Trans. Nucl. Sci., Vol.NS-

, 1984, pp. 1322-1326.

J. S. Shang, “A fractional-step method for solving 3D,

time-domain Maxwell Equations,” J. of. Comput.

Phys., Vol. 118, 1195, pp. 109-119.

A.Fijany, M. A. Jensen, Y. Rahmat-Samii and J.

Barhen, “A massively parallel computation strategy for

FDTD: time and space parallelism applied to

electromagnetics problems,” IEEE Trans. Antenna and

Propagation, Vol.43, Dec., 1995, pp. 1441-1449.

T. Namiki, “A new FDTD algorithm based on

alternating-direction implicit method,” IEEE Trans.

Microwave Theory Tech., Vol. 47, Oct., 1999, pp.

-2007.

F. Zhen, Z. Chen and J. Zhang, “Toward the

development of a three-dimensional unconditionally

stable finite-difference time-domain method,” IEEE

Trans. Microwave Theory Tech., Vol.48, Sept., 2000,

pp. 1550-1558.

J.H. Beggs and W. R. Briley, “An implicit characteristic

based method for electromagnetics,” NASA/TM-2001-

, May, 2001

G. Sun and C.W. Trueman, “Unconditionally stable

Crank-Nicolson scheme for solving the two-

dimensional Maxwell’s Equations,” IEE Electronics

Lett., Vol. 39, April, 2003, pp. 595-597.

G. Sun, C. W. Trueman, “Analysis and Numerical

Experiments on the Numerical Dispersion of Two-

Dimensional ADI-FDTD”, IEEE Antennas and

Wireless Propaga. Letters, to be published.

J. C. Tannehill, D. A. Anderson, R.H. Pletcher,

“Computational fluid mechanics and heat transfer,”

nd Ed., Taylor&Francis, Washington, DC, 1997.

Y. Liu, “Fourier analysis of numerical algorithms for

the Maxwell’s Equations,” J. Comput. Phys., Vol. 124,

March, 1996, pp. 396-416

##submission.downloads##

已出版

2022-06-18

栏目

General Submission