FDTD Schemes for Maxwell’s Equations with Embedded Perfect Electric Conductors Based on the Correction Function Method

In this work, we propose staggered FDTD schemes based on the correction function method (CFM) to discretize Maxwell's equations with embedded perfect electric conductor (PEC) boundary conditions. The CFM uses a minimization procedure to compute a correction to a given FD scheme in the vicinity of the embedded boundary to retain its order. The minimization problem associated with CFM approaches is analyzed in the context of Maxwell's equations with embedded boundaries. In order to obtain a well-posed problem, we propose fictitious interface conditions to fulfill the lack of information, namely the surface current and charge density, on the embedded boundary. Fictitious interfaces can induce some issues for long time simulations and therefore the penalization coefficient associated with fictitious interface conditions must be chosen small enough. We introduce CFM-FDTD schemes based on the well-known Yee scheme and a fourth-order staggered FDTD scheme. Long time simulations and convergence studies are performed in 2-D for various geometries of the embedded boundary. CFM-FDTD schemes have shown high-order convergence.

[1]  Zhiming Chen,et al.  The adaptive immersed interface finite element method for elliptic and Maxwell interface problems , 2009, J. Comput. Phys..

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

[3]  Guo-Wei Wei,et al.  Matched interface and boundary (MIB) method for elliptic problems with sharp-edged interfaces , 2007, J. Comput. Phys..

[4]  B. Griffith,et al.  An immersed boundary method for rigid bodies , 2014, 1505.07865.

[5]  Gary Cohen,et al.  A spatial high-order hexahedral discontinuous Galerkin method to solve Maxwell's equations in time domain , 2006, J. Comput. Phys..

[6]  Jean-Christophe Nave,et al.  Imposing jump conditions on nonconforming interfaces for the Correction Function Method: A least squares approach , 2019, J. Comput. Phys..

[7]  Chi-Wang Shu,et al.  Locally divergence-free discontinuous Galerkin methods for the Maxwell equations , 2004, Journal of Computational Physics.

[8]  Jean-Christophe Nave,et al.  Treatment of Complex Interfaces for Maxwell’s Equations with Continuous Coefficients Using the Correction Function Method , 2018, J. Sci. Comput..

[9]  Jan S. Hesthaven,et al.  Spectral Simulations of Electromagnetic Wave Scattering , 1997 .

[10]  Jean-Christophe Nave,et al.  A correction function method for the wave equation with interface jump conditions , 2016, J. Comput. Phys..

[11]  Jean-Christophe Nave,et al.  A Correction Function Method for Poisson problems with interface jump conditions , 2010, J. Comput. Phys..

[12]  Qing Huo Liu,et al.  Multidomain pseudospectral time-domain simulations of scattering by objects buried in lossy media , 2002, IEEE Trans. Geosci. Remote. Sens..

[13]  Jan S. Hesthaven,et al.  High-order accurate methods in time-domain computational electromagnetics: A review , 2003 .

[14]  H. Kreiss,et al.  Time-Dependent Problems and Difference Methods , 1996 .

[15]  Jean-Christophe Nave,et al.  A Fourier penalty method for solving the time-dependent Maxwell's equations in domains with curved boundaries , 2016, J. Comput. Phys..

[16]  Jean-Christophe Nave,et al.  High order solution of Poisson problems with piecewise constant coefficients and interface jumps , 2014, J. Comput. Phys..

[17]  Shan Zhao,et al.  A fourth order finite difference method for waveguides with curved perfectly conducting boundaries , 2010 .

[18]  Tobin A. Driscoll,et al.  Staggered Time Integrators for Wave Equations , 2000, SIAM J. Numer. Anal..

[19]  Shan Zhao,et al.  High-order FDTD methods via derivative matching for Maxwell's equations with material interfaces , 2004 .

[20]  J. Hesthaven,et al.  Convergent Cartesian Grid Methods for Maxwell's Equations in Complex Geometries , 2001 .

[21]  J. Hesthaven,et al.  Nodal high-order methods on unstructured grids , 2002 .

[22]  Shan Zhao,et al.  Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves , 2016 .

[23]  Robert D. Guy,et al.  Immersed boundary smooth extension: A high-order method for solving PDE on arbitrary smooth domains using Fourier spectral methods , 2015, J. Comput. Phys..