An upwinding boundary condition capturing method for Maxwell's equations in media with material interfaces

By using ghost points on either side of the interfaces, a global second-order accurate upwinding boundary condition capturing method for time-domain Maxwell's equations in media with material interfaces is proposed. The equations are discretized on a uniform Cartesian grid and the interfaces are allowed to intersect the grid in an arbitrary fashion. The method is then obtained by combining central finite difference schemes with applicable nodes being replaced by the ghost points and upwinding technique with jump conditions across the interfaces being captured in a manner that the upwind property is always satisfied. The resulting discretization has the desirable property that the allowed time step size is independent of the locations and the shapes of the interfaces. Numerical examples are then given to demonstrate the second-order accuracy as well as the stability of the method, where it is used to study wave equations with various types of material interfaces, including electromagnetic scattering of a plane incident wave by a dielectric circular cylinder.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

[15]  T. Aslam.  A partial differential equation approach to multidimensional extrapolation , 2004 .

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

[17]  David Gottlieb,et al.  On the construction and analysis of absorbing layers in CEM , 1998 .

[18]  J. Strikwerda.  Finite Difference Schemes and Partial Differential Equations , 1989 .

[19]  H. Kreiss.  Stability theory for difference approximations of mixed initial boundary value problems. I , 1968 .

[20]  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 .

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

[22]  J. Hesthaven.  High-Order Accurate Methods in Time-Domain Computational Electromagnetics , 2002 .

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

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

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

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

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

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

[29]  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 .

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

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

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

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

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

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

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

[37]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

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

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

[40]  J. Hesthaven,et al.  Modeling Dielectric Interfaces in the FDTD-Method: A Comparative Study , 2000 .

[41]  H. Kreiss,et al.  Stability Theory of Difference Approximations for Mixed Initial Boundary Value Problems. II , 1972 .

[42]  Richard W. Ziolkowski,et al.  Time-derivative Lorentz material model-based absorbing boundary condition , 1997 .