On the immersed interface method for solving time-domain Maxwell's equations in materials with curved dielectric interfaces

This paper deals with accurate numerical simulation of two-dimensional time-domain Maxwell's equations in materials with curved dielectric interfaces. The proposed fully second-order scheme is a hybridization between the immersed interface method (IIM), introduced to take into account curved geometries in structured schemes, and the Lax-Wendroff scheme, usually used to improve order of approximations in time for partial differential equations. In particular, the IIM proposed for two-dimensional acoustic wave equations with piecewise constant coefficients [C. Zhang, R.J. LeVeque, The immersed interface method for acoustic wave equations with discontinuous coefficients, Wave Motion 25 (1997) 237-263] is extended through a simple least squares procedure to such Maxwell's equations. Numerical results from the simulation of electromagnetic scattering of a plane incident wave by a dielectric circular cylinder appear to indicate that, compared to the original IIM for the acoustic wave equations, the augmented IIM with the proposed least squares fitting greatly improves the long-time stability of the time-domain solution. Semi-discrete finite difference schemes using the IIM for spatial discretization are also discussed and numerically tested in the paper.

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

[2]  Kendall E. Atkinson.  An introduction to numerical analysis , 1978 .

[3]  Zhilin Li.  AN OVERVIEW OF THE IMMERSED INTERFACE METHOD AND ITS APPLICATIONS , 2003 .

[4]  Qing Huo Liu,et al.  A Staggered Upwind Embedded Boundary (SUEB) method to eliminate the FDTD staircasing error , 2004 .

[5]  Bo Zhang,et al.  An explicit fourth-order staggered finite-difference time-domain method for Maxwell's equations , 2002 .

[6]  Chi Hou Chan,et al.  An explicit fourth-order orthogonal curvilinear staggered-grid FDTD method for Maxwell's equations , 2002 .

[7]  Zhilin Li,et al.  The Immersed Interface Method: Numerical Solutions of PDEs Involving Interfaces and Irregular Domains (Frontiers in Applied Mathematics) , 2006 .

[8]  Shaozhong Deng,et al.  An upwinding boundary condition capturing method for Maxwell's equations in media with material interfaces , 2007, J. Comput. Phys..

[9]  Kazufumi Ito,et al.  Maximum Principle Preserving Schemes for Interface Problems with Discontinuous Coefficients , 2001, SIAM J. Sci. Comput..

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

[11]  R. LeVeque,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[12]  Eli Turkel,et al.  Fourth order compact implicit method for the Maxwell equations with discontinuous coefficients , 2000 .

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

[14]  Amir Yefet,et al.  A staggered fourth-order accurate explicit finite difference scheme for the time-domain Maxwell's equations , 2001 .

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

[16]  R. LeVeque,et al.  The immersed interface method for acoustic wave equations with discontinuous coefficients , 1997 .

[17]  Wei Cai,et al.  An upwinding embedded boundary method for Maxwell's equations in media with material interfaces: 2D case , 2003 .