Energy-preserving local mesh-refined splitting FDTD schemes for two dimensional Maxwell's equations

Abstract In this paper, we develop and analyze two types of new energy-preserving local mesh-refined splitting finite difference time-domain (EP-LMR-S-FDTD) schemes for two-dimensional Maxwell's equations. For the local mesh refinements, it is challenging to define the suitable local interface schemes which can preserve energy and guarantee high accuracy. The important feature of the work is that based on energy analysis, we propose the efficient local interface schemes on the interfaces of coarse and fine grids that ensure the energy conservation property, keep spatial high accuracy and avoid oscillations and meanwhile, we propose a fast implementation of the EP-LMR-S-FDTD schemes, which overcomes the difficulty in solving unknowns on the “trifuecate structure” of refinement by first solving the values of course mesh unknowns and the average values of fine mesh unknowns on a line-structure and then solving the values of fine mesh unknowns and the course mesh unknown on an inverted “U-form” structure for each loop. The EP-LMR-S-FDTD schemes can be solved in a series of tridiagonal linear systems of unknowns which can be efficiently implemented at each time step. We prove the EP-LMR-S-FDTD schemes to be energy preserving and unconditionally stable. We further prove the convergence of the schemes and obtain the error estimates. Numerical experiments are given to show the performance of the EP-LMR-S-FDTD schemes which confirm theoretical results.

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

[2]  M. Brio,et al.  FDTD based second-order accurate local mesh refinement method for Maxwell's equations in two space dimensions , 2004 .

[3]  Francesca Rapetti,et al.  MORTARING THE TWO-DIMENSIONAL EDGE FINITE ELEMENTS FOR THE DISCRETIZATION OF SOME ELECTROMAGNETIC MODELS , 2004 .

[4]  P. Moin,et al.  Fully Conservative Higher Order Finite Difference Schemes for Incompressible Flow , 1998 .

[5]  Tomi,et al.  THE USE OF PLANE WAVES TO APPROXIMATE WAVE PROPAGATION IN ANISOTROPIC MEDIA , 2007 .

[6]  I. S. Kim,et al.  A local mesh refinement algorithm for the time domain-finite difference method using Maxwell's curl equations , 1990 .

[7]  Dong Liang,et al.  The spatial fourth-order energy-conserved S-FDTD scheme for Maxwell's equations , 2013, J. Comput. Phys..

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

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

[10]  A. B. Strong,et al.  A Fully Conservative Second-Order Finite Difference Scheme for Incompressible Flow on Nonuniform Grids , 2002 .

[11]  R. LeVeque,et al.  Adaptive Mesh Refinement Using Wave-Propagation Algorithms for Hyperbolic Systems , 1998 .

[13]  Patrick Joly,et al.  A Conservative Space-time Mesh Refinement Method for the 1-D Wave Equation. Part II: Analysis , 2003, Numerische Mathematik.

[14]  Nathan Louis Gibson,et al.  Analysis of spatial high-order finite difference methods for Maxwell's equations in dispersive media , 2012 .

[15]  Faker Ben Belgacem,et al.  The Mortar finite element method with Lagrange multipliers , 1999, Numerische Mathematik.

[16]  Zhizhang Chen,et al.  A hybrid ADI-FDTD scheme and its comparisons with the FDTD method , 2003, IEEE Antennas and Propagation Society International Symposium. Digest. Held in conjunction with: USNC/CNC/URSI North American Radio Sci. Meeting (Cat. No.03CH37450).

[17]  Ping Lin,et al.  Mass conservative and energy stable finite difference methods for the quasi-incompressible Navier–Stokes–Cahn–Hilliard system: Primitive variable and projection-type schemes , 2017, 1703.06606.

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

[19]  Wenbin Chen,et al.  A Multigrid Method for the Mortar-Type Morley Element Approximation of a Plate Bending Problem , 2001, SIAM J. Numer. Anal..

[20]  P. Joly,et al.  Space-time mesh refinement for elastodynamics. Numerical results , 2005 .

[21]  F. BEN BELGACEM,et al.  The Mortar Finite Element Method for 3D Maxwell Equations: First Results , 2001, SIAM J. Numer. Anal..

[22]  Oleg V. Vasilyev High Order Finite Difference Schemes on Non-uniform Meshes with Good Conservation Properties , 2000 .

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

[24]  Dong Liang,et al.  Energy-conserved splitting FDTD methods for Maxwell’s equations , 2007, Numerische Mathematik.

[25]  Patrick Joly,et al.  An Error Analysis of Conservative Space-Time Mesh Refinement Methods for the One-Dimensional Wave Equation , 2005, SIAM J. Numer. Anal..

[26]  Wenbin Chen,et al.  Symmetric Energy-Conserved Splitting FDTD Scheme for the Maxwell's Equations , 2009 .