Stability and Super Convergence Analysis of ADI-FDTD for the 2D Maxwell Equations in a Lossy Medium

Abstract Several new energy identities of the two dimensional(2D) Maxwell equations in a lossy medium in the case of the perfectly electric conducting boundary conditions are proposed and proved. These identities show a new kind of energy conservation in the Maxwell system and provide a new energy method to analyze the alternating direction implicit finite difference time domain method for the 2D Maxwell equations (2D-ADI-FDTD). It is proved that 2D-ADI-FDTD is approximately energy conserved, unconditionally stable and second order convergent in the discrete L2 and H1 norms, which implies that 2D-ADI-FDTD is super convergent. By this super convergence, it is simply proved that the error of the divergence of the solution of 2D-ADI-FDTD is second order accurate. It is also proved that the difference scheme of 2D-ADI-FDTD with respect to time t is second order convergent in the discrete H1 norm. Experimental results to confirm the theoretical analysis on stability, convergence and energy conservation are presented.

[1]  E. Süli,et al.  A convergence analysis of Yee's scheme on nonuniform grids , 1994 .

[2]  Jingjing Zhang,et al.  Splitting multisymplectic integrators for Maxwell's equations , 2010, J. Comput. Phys..

[3]  T. Namiki,et al.  A new FDTD algorithm based on alternating-direction implicit method , 1999 .

[4]  Tae-Woo Lee,et al.  On the accuracy of the ADI-FDTD method , 2002, IEEE Antennas and Wireless Propagation Letters.

[5]  Bo Zhang,et al.  Optimal error estimates and modified energy conservation identities of the ADI-FDTD scheme on staggered grids for 3D Maxwell’s equations , 2013 .

[6]  John R. Cary,et al.  Divergence preservation in the ADI algorithms for electromagnetics , 2009, J. Comput. Phys..

[7]  H. H. Rachford,et al.  The Numerical Solution of Parabolic and Elliptic Differential Equations , 1955 .

[8]  Qun Lin,et al.  Superconvergence analysis for Maxwell's equations in dispersive media , 2007, Math. Comput..

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

[10]  Roy A. Nicolaides,et al.  Convergence analysis of a covolume scheme for Maxwell's equations in three dimensions , 1998, Math. Comput..

[11]  B. D. Welfert Analysis of iterated ADI-FDTD schemes for Maxwell curl equations , 2007, J. Comput. Phys..

[12]  Dong Liang,et al.  Energy-Conserved Splitting Finite-Difference Time-Domain Methods for Maxwell's Equations in Three Dimensions , 2010, SIAM J. Numer. Anal..

[13]  H. H. Rachford,et al.  On the numerical solution of heat conduction problems in two and three space variables , 1956 .

[14]  An Ping Zhao,et al.  Analysis of the numerical dispersion of the 2D alternating-direction implicit FDTD method , 2002 .

[15]  Qun Lin,et al.  Global superconvergence for Maxwell's equations , 2000, Math. Comput..

[16]  Fenghua Zhen,et al.  Toward the development of a three-dimensional unconditionally stable finite-difference time-domain method , 2000 .

[17]  Gang Liu,et al.  Perfectly matched layer media with CFS for an unconditionally stable ADI-FDTD method , 2001 .

[18]  Fang Nengsheng,et al.  Analysis of FDTD to UPML for Maxwell equations in polar coordinates , 2011 .

[19]  Tie Zhang,et al.  Superconvergence of discontinuous Galerkin methods for hyperbolic systems , 2009 .

[20]  Y. Hao,et al.  ADI-FDTD algorithm in curvilinear co-ordinates , 2005 .

[21]  Nengsheng Fang,et al.  Stability analysis of FDTD to UPML for time dependent Maxwell equations , 2009 .