A Fourth-Order Upwinding Embedded Boundary Method (UEBM) for Maxwell's Equations in Media with Material Interfaces: Part I

In this paper, we present a new fourth-order upwinding embedded boundary method (UEBM) over Cartesian grids, originally proposed in the Journal of Computational Physics (190, pp. 159, 2003) as a second-order method for treating material interfaces for Maxwell's equations. In addition to the idea of the UEBM to evolve solutions at interfaces, we utilize the ghost fluid method to construct finite dierence approximation of spatial deriva- tives at Cartesian grid points near the material interfaces. As a result, Runge-Kutta type time discretization can be used for the semi-discretized system to yield an overall fourth-order method, in contrast to the original second-order UEBM based on a Lax-Wendro type dier- ence. The final scheme allows time step sizes independent of the interface locations. Numerical examples are given to demonstrate the fourth-order accuracy as well as the stability of the method. We tested the scheme for several wave problems with various material interface lo- cations, including electromagnetic scattering of a plane wave incident on a planar boundary and a two-dimensional electromagnetic application with an interface parallel to the y-axis.

[1]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

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

[3]  S.,et al.  Numerical Solution of Initial Boundary Value Problems Involving Maxwell’s Equations in Isotropic Media , 1966 .

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

[5]  R. Fedkiw,et al.  A Boundary Condition Capturing Method for Poisson's Equation on Irregular Domains , 2000 .

[6]  Ronald Fedkiw,et al.  A Boundary Condition Capturing Method for Multiphase Incompressible Flow , 2000, J. Sci. Comput..

[7]  B. Gustafsson.  The convergence rate for difference approximations to mixed initial boundary value problems , 1975 .

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

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

[10]  Heinz-Otto Kreiss,et al.  Difference Approximations for the Second Order Wave Equation , 2002, SIAM J. Numer. Anal..

[11]  Joël Piraux,et al.  Numerical treatment of two-dimensional interfaces for acoustic and elastic waves , 2004 .

[12]  M.Y. Hussaini,et al.  Low-Dissipation and Low-Dispersion Runge-Kutta Schemes for Computational Acoustics , 1994 .

[13]  Endre Süli,et al.  Error estimates for Yee's method on non-uniform grids , 1994 .

[14]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

[15]  Xiangyu Hu,et al.  An interface interaction method for compressible multifluids , 2004 .

[16]  Bertil Gustafsson,et al.  The convergence rate for difference approximations to general mixed initial boundary value problems , 1981 .

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

[18]  Heinz-Otto Kreiss,et al.  Difference Approximations of the Neumann Problem for the Second Order Wave Equation , 2004, SIAM J. Numer. Anal..

[19]  Randall J. LeVeque,et al.  H-Box Methods for the Approximation of Hyperbolic Conservation Laws on Irregular Grids , 2003, SIAM J. Numer. Anal..

[20]  Heinz-Otto Kreiss,et al.  An Embedded Boundary Method for the Wave Equation with Discontinuous Coefficients , 2005, SIAM J. Sci. Comput..

[21]  R. Fedkiw,et al.  A fourth order accurate discretization for the Laplace and heat equations on arbitrary domains, with applications to the Stefan problem , 2005 .

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

[23]  F. Xiao,et al.  Three-dimensional numerical simulation of flows with complex geometries in a regular Cartesian grid and its application to blood flow in cerebral artery with multiple aneurysms , 2005 .

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

[25]  HEINZ-OTTO KREISS,et al.  A Second Order Accurate Embedded Boundary Method for the Wave Equation with Dirichlet Data , 2005, SIAM J. Sci. Comput..

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

[27]  Li-Tien Cheng,et al.  A second-order-accurate symmetric discretization of the Poisson equation on irregular domains , 2002 .

[28]  R. Fedkiw,et al.  A boundary condition capturing method for incompressible flame discontinuities , 2001 .

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

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

[31]  Ronald Fedkiw,et al.  Regular Article: The Ghost Fluid Method for Deflagration and Detonation Discontinuities , 1999 .

[32]  John B. Bell,et al.  Cartesian grid method for unsteady compressible flow in irregular regions , 1995 .

[33]  Paul Lorrain,et al.  Electromagnetic fields and waves : including electric circuits , 1988 .