×

Investigation of numerical time-integrations of Maxwell’s equations using the staggered grid spatial discretization. (English) Zbl 1108.78016

The authors have considered 4 numerical methods for solving the Maxwell equations: Yee-method, Namiki-Zheng-Chen-Zhang alternating direction implicit method (NZCZ), the Kole-Figge-de Raedt method (KFR) and a Krylov-space method. They have showed how these methods are connected with the approximation of matrix exponential. The authors proved the unconditional stability of the NZCZ method based on this approach. The methods are compared in many numerical experiments.

MSC:

78M20 Finite difference methods applied to problems in optics and electromagnetic theory
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs

Software:

Expokit
PDFBibTeX XMLCite
Full Text: DOI Link

References:

[1] Yee, Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media, IEEE Transactions on Antennas and Propagation 14 (3) pp 302– (1966) · Zbl 1155.78304
[2] Namiki, A new FDTD algorithm based on alternating-direction implicit method, IEEE Transactions on Microwave Theory and Techniques 47 (10) pp 2003– (1999)
[3] Zheng, Toward the development of a three-dimensional unconditionally stable finite-difference time-domain method, IEEE Transactions on Microwave Theory and Techniques 48 (9) pp 1550– (2000)
[4] Kole, Unconditionally stable algorithms to solve the time-dependent Maxwell equations, Physical Review E 64 pp 066705– (2001)
[5] Taflove, Computational Electrodynamics: The Finite-Difference Time-Domain Method (2000) · Zbl 0963.78001
[6] Holland, Implicit three-dimensional finite differencing of Maxwell’s equations, IEEE Transactions on Nuclear Science NS-31 pp 1322– (1984)
[7] Holland, Technical Report to Harry Diamond Labs (1986)
[8] Darms, Asymmetry effects in the ADI-FDTD algorithm, IEEE Microwave Guided Wave Letters 12 pp 491– (2002)
[9] Fornberg, Technical Report #472 (2001)
[10] Fornberg, Topics in Computational Wave Propagation pp 265– (2003)
[11] Hochbruck, On Krylov-Subspace approximations to the matrix exponential operator, SIAM Journal on Numerical Analysis 34 (5) pp 1911– (1997) · Zbl 0888.65032
[12] Remis, A modified lanczos algorithm for the computation of transient electromagnetic wavefields, IEEE Transactions on Microwave Theory and Techniques 45 (12) pp 2139– (1997)
[13] Sidje RB EXPOKIT: software package for computing matrix exponentials 1998
[14] Sidje, A numerical study of large sparse matrix exponentials arising in Markov chains, Computational Statistics and Data Analysis 29 pp 345– (1999) · Zbl 1042.65508
[15] van der Vorst, An iterative solution method for solving f(A)x=b, using Krylov subspace information obtained for the symmetric positive definite matrix A, Journal of Computational and Applied Mathematics 18 pp 249– (1987) · Zbl 0621.65022
[16] Golub, Matrix Computations (1996)
[17] Richmyer, Difference Methods for Initial Value Problems (1967)
[18] Yoshida, Construction of higher order symplectic integrators, Physics Letters A 150 pp 262– (1990)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.