×

Convergence analysis of Yee-FDTD schemes for 3D Maxwell’s equations in linear dispersive media. (English) Zbl 1499.65428

Summary: In this paper, we develop and analyze finite difference methods for the 3D Maxwell’s equations in the time domain in three different types of linear dispersive media described as Debye, Lorentz and cold plasma. These methods are constructed by extending the Yee-Finite Difference Time Domain (FDTD) method to linear dispersive materials. We analyze the stability criterion for the FDTD schemes by using the energy method. Based on energy identities for the continuous models, we derive discrete energy estimates for the FDTD schemes for the three dispersive models. We also prove the convergence of the FDTD schemes with perfect electric conducting boundary conditions, which describes the second order accuracy of the methods in both time and space. The discrete divergence-free conditions of the FDTD schemes are studied. Lastly, numerical examples are given to demonstrate and confirm our results.

MSC:

65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
65M12 Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs
78M20 Finite difference methods applied to problems in optics and electromagnetic theory
PDFBibTeX XMLCite
Full Text: Link

References:

[1] Aziz, M. M. Sub-nanosecond electromagnetic-micromagnetic dynamic simulations using the finite-difference time-domain method. Progress In Electromagnetics Research, 15 (2009), 1-29.
[2] Banks, H. T., Buksas, M. W., and Lin, T. Electromagnetic Material Interrogation Using Conductive Interfaces and Acoustic Wavefronts, vol. 21. Society for Industrial Mathematics (SIAM), 2000. · Zbl 1008.78500
[3] Bécache, E., and Joly, P. On the analysis of Bérenger’s perfectly matched layers for Maxwell’s equations. Mathematical Modelling and Numerical Analysis, 36, 1 (2002), 87-119. · Zbl 0992.78032
[4] Bidégaray-Fesquet, B. Stability of FD-TD schemes for Maxwell-Debye and Maxwell-Lorentz equations. SIAM Journal on Numerical Analysis, 46, 5 (2008), 2551-2566. · Zbl 1177.78056
[5] Bokil, V. A., Cheng, Y., Jiang, Y., Li, F., and Sakkaplangkul, P. High spatial order energy stable FDTD methods for Maxwell’s equations in nonlinear optical media in one dimension. Journal of Scientific Computing, 77 (2018), 330-371. · Zbl 1407.65094
[6] Bokil, V. A., and Gibson, N. L. Analysis of Spatial High-Order Finite Difference Methods for Maxwell’s Equations in Dispersive Media. IMA Journal of Numerical Analysis 32, 3 (2012), 926-956. · Zbl 1253.78048
[7] Bokil, V. A., and Gibson, N. L. Convergence analysis of Yee schemes for Maxwell’s equations in Debye and Lorentz dispersive media. International Journal of Numerical Analysis and Modeling, 11, 4 (2014), 657-687. · Zbl 1482.65137
[8] Bokil, V. A., Keefer, O., and Leung, A.-Y. Operator splitting methods for Maxwell’s equa-tions in dispersive media with orientational polarization. Journal of Computational and Applied Mathematics, 263 (2014), 160-188. · Zbl 1301.65091
[9] Bokil, V. A., and Sakkaplangkul, P. Construction and analysis of weighted sequential split-ting FDTD methods for the 3D Maxwell’s equations. International Journal of Numerical Analysis and Modeling, 15, 6 (2018), 747-784. · Zbl 1412.65067
[10] Causley, M. F., Petropoulos, P. G., and Jiang, S. Incorporating the Havriliak-Negami di-electric model in the FD-TD method. Journal of Computational Physics, 230, 10 (2011), 3884-3899. · Zbl 1219.78136
[11] Chen, W., Li, X., and Liang, D. Energy-conserved splitting FDTD methods for Maxwell’s equations. Numerische Mathematik, 108, 3 (2008), 445-485. · Zbl 1185.78020
[12] Chen, W., Li, X., and Liang, D. Energy-conserved splitting finite-difference time-domain methods for Maxwell’s equations in three dimensions. SIAM Journal on Numerical Analysis, 48, 4 (2010), 1530-1554. · Zbl 1220.78116
[13] Cole, K. S., and Cole, R. H. Dispersion and absorption in dielectrics i. alternating current characteristics. Journal of Chemical Physics, 9 (1941), 341-351.
[14] Cummer, S. A. An analysis of new and existing FDTD methods for isotropic cold plasma and a method for improving their accuracy. IEEE Transactions on Antennas and Propagation, 45, 3 (1997), 392-400.
[15] Da Silva, F., Pinto, M. C., Després, B., and Heuraux, S. Stable explicit coupling of the yee scheme with a linear current model in fluctuating magnetized plasmas. Journal of Compu-tational Physics, 295 (2015), 24-45. · Zbl 1349.76563
[16] Debye, P. J. W. Polar Molecules. Dover New York, 1929. · JFM 55.1179.04
[17] Gansen, A., El Hachemi, M., Belouettar, S., Hassan, O., and Morgan, K. A 3D Unstructured Mesh FDTD Scheme for EM Modelling. Archives of Computational Methods in Engineering, (2020), 1-33.
[18] Gao, L., and Zhang, B. Stability and superconvergence analysis of the FDTD scheme for the 2D Maxwell equations in a lossy medium. Science China Mathematics, 54, 12 (2011), 2693-2712. · Zbl 1268.78023
[19] Gedney, S. D. Introduction to the finite-difference time-domain (fdtd) method for electro-magnetics. Synthesis Lectures on Computational Electromagnetics, 6, 1 (2011), 1-250. · Zbl 1271.78001
[20] Gilles, L., Hagness, S., and Vázquez, L. Comparison between staggered and unstaggered finite-difference time-domain grids for few-cycle temporal optical soliton propagation. Jour-nal of Computational Physics, 161, 2 (2000), 379-400. · Zbl 0960.78015
[21] Hao, Y., and Mittra, R. FDTD modeling of metamaterials: Theory and applications. Artech house, 2008.
[22] Huang, Y., Li, J., and Yang, W. Interior penalty DG methods for Maxwell’s equations in dispersive media. Journal of Computational Physics, 230, 12 (2011), 4559-4570. · Zbl 1220.78015
[23] Jackson, J. D. Classical electrodynamics. John Wiley & Sons, 2007. · Zbl 0114.42903
[24] James, R. Wait, electromagnetic waves in stratified media. Oxford: Pergmon Press, 2 (1962), 17. · Zbl 0122.45005
[25] Jiang, Y., Sakkaplangkul, P., Bokil, V. A., Cheng, Y., and Li, F. Dispersion analysis of finite difference and discontinuous Galerkin schemes for Maxwell’s equations in linear Lorentz media. Journal of Computational Physics, 394 (2019), 100-135. · Zbl 1452.78032
[26] Joseph, R. M., Hagness, S. C., and Taflove, A. Direct time integration of Maxwell’s equations in linear dispersive media with absorption for scattering and propagation of femtosecond electromagnetic pulses. Optics Letters, 16, 18 (1991), 1412-1414.
[27] Kashiwa, T., and Fukai, I. A treatment by the FD-TD method of the dispersive character-istics associated with electronic polarization. Microwave and Optical Technology Letters, 3, 6 (1990), 203-205.
[28] Kashiwa, T., Yoshida, N., and Fukai, I. A treatment by the finite-difference time domain method of the dispersive characteristics associated with orientational polarization. IEEE Trans. IEICE, 73, 8 (1990), 1326-1328.
[29] Lanteri, S., and Scheid, C. Convergence of a discontinuous Galerkin scheme for the mixed time-domain Maxwell’s equations in dispersive media. IMA Journal of Numerical Analysis, 33, 2 (2013), 432-459. · Zbl 1277.78035
[30] Li, J. Numerical convergence and physical fidelity analysis for Maxwell’s equations in meta-materials. Computer Methods in Applied Mechanics and Engineering, 198, 37 (2009), 3161-3172. · Zbl 1229.78024
[31] Li, J. Unified analysis of leap-frog methods for solving time-domain Maxwell’s equations in dispersive media. Journal of Scientific Computing, 47, 1 (2011), 1-26. · Zbl 1230.78030
[32] Li, J., and Huang, Y. Time-domain finite element methods for Maxwell’s equations in meta-materials, vol. 43. Springer, 2013. · Zbl 1304.78002
[33] Li, J., Huang, Y., and Lin, Y. Developing finite element methods for Maxwell’s equations in a Cole-Cole dispersive medium. SIAM Journal on Scientific Computing, 33, 6 (2011), 3153-3174. · Zbl 1238.65097
[34] Li, J., and Shields, S. Superconvergence analysis of yee scheme for metamaterial maxwells equations on non-uniform rectangular meshes. Numerische Mathematik, 134, 4 (2016), 741-781. · Zbl 1359.78025
[35] Li, W., and Liang, D. The energy conservative splitting FDTD scheme and its energy iden-tities for metamaterial electromagnetic lorentz system. Computer Physics Communications, 239 (2019), 94-111. · Zbl 07684938
[36] Liang, D., and Yuan, Q. The spatial fourth-order energy-conserved S-FDTD scheme for Maxwells equations. Journal of Computational Physics, 243 (2013), 344-364. · Zbl 1349.78094
[37] Lu, T., Zhang, P., and Cai, W. Discontinuous Galerkin methods for dispersive and lossy Maxwell’s equations and PML boundary conditions. Journal of Computational Physics, 200, 2 (2004), 549-580. · Zbl 1115.78330
[38] Luebbers, R., Hunsberger, F. P., Kunz, K. S., Standler, R. B., and Schneider, M. A frequency-dependent finite-difference time-domain formulation for dispersive materials. IEEE Transactions on Electromagnetic Compatibility, 32, 3 (1990), 222-227.
[39] Luebbers, R. J., Hunsberger, F., and Kunz, K. S. A frequency-dependent finite-difference time-domain formulation for transient propagation in plasma. IEEE Transactions on Anten-nas and Propagation, 39, 1 (1991), 29-34.
[40] Monk, P., and Süli, E. A convergence analysis of Yee’s scheme on nonuniform grids. SIAM Journal on Numerical Analysis, 31, 2 (1994), 393-412. · Zbl 0805.65121
[41] Monk, P., and Suli, E. Error estimates for yee’s method on non-uniform grids. IEEE Trans-actions on Magnetics, 30, 5 (1994), 3200-3203.
[42] Oughstun, K. E., and Sherman, G. C. Uniform asymptotic description of electromagnetic pulse propagation in a linear dispersive medium with absorption (the Lorentz medium). Journal of the Optical Society of America A, 6, 9 (1989), 1394-1420.
[43] Petropoulos, P. G. Stability and phase error analysis of FDTD in dispersive dielectrics. Antennas and Propagation, IEEE Transactions on, 42, 1 (1994), 62-69.
[44] Stix, T. Waves in Plasmas. American Institute of Physics, 1992.
[45] Sun, J., Collino, F., Monk, P., and Wang, L. An eddy-current and micromagnetism model with applications to disk write heads. International Journal for Numerical Methods in En-gineering, 60, 10 (2004), 1673-1698. · Zbl 1109.78325
[46] Swanson, D. G. Plasma waves. Elsevier, 2012.
[47] Taflove, A., and Hagness, S. C. Computational electrodynamics: The Finite-Difference Time-Domain method, 3rd ed. Artech House, Norwood, MA, 2005.
[48] Wang, B., Xie, Z., and Zhang, Z. Error analysis of a discontinuous Galerkin method for Maxwell equations in dispersive media. Journal of Computational Physics, 229, 22 (2010), 8552-8563. · Zbl 1203.78052
[49] Xie, J., Liang, D., and Zhang, Z. Energy-preserving local mesh-refined splitting FDTD schemes for two dimensional Maxwell’s equations. Journal of Computational Physics, 425 (2021), 109896. · Zbl 07508494
[50] Yee, K. Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media. Antennas and Propagation, IEEE Transactions on, 14, 3 (1966), 302-307. · Zbl 1155.78304
[51] Young, J. A full finite difference time domain implementation for radio wave propagation in a plasma. Radio Science, 29, 6 (1994), 1513-1522.
[52] Young, J. A higher order FDTD method for EM propagation in a collisionless cold plasma. IEEE Antennas and Propagation Magazine, 44, 9 (1996), 1283-1289.
[53] Young, J. L., and Nelson, R. O. A Summary and Systematic Analysis of FDTD Algorithms for Linearly Dispersive Media. IEEE Antennas and Propagation Magazine, 43 (2001), 61-77.
[54] Department of Mathematics, Faculty of Science, King Mongkut’s Institute of Technology Lad-krabang, Ladkrabang, Bangkok 10520, Thailand. E-mail: puttha.sa@kmitl.ac.th
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.