Time-domain matched interface and boundary (MIB) modeling of Debye dispersive media with curved interfaces

A new finite-difference time-domain (FDTD) method is introduced for solving transverse magnetic Maxwell's equations in Debye dispersive media with complex interfaces and discontinuous wave solutions. Based on the auxiliary differential equation approach, a hybrid Maxwell-Debye system is constructed, which couples the wave equation for the electric component with Maxwell's equations for the magnetic components. This hybrid formulation enables the calculation of the time dependent parts of the interface jump conditions, so that one can track the transient changes in the regularities of the electromagnetic fields across a dispersive interface. Effective matched interface and boundary (MIB) treatments are proposed to rigorously impose the physical jump conditions which are not only time dependent, but also couple both Cartesian directions and both magnetic field components. Based on a staggered Yee lattice, the proposed MIB scheme can deal with arbitrarily curved interfaces and nonsmooth interfaces with sharped edges. Second order convergences are numerically achieved in solving dispersive interface problems with constant curvatures, general curvatures, and nonsmooth corners.

[1]  S. S. Stuchly,et al.  Propagation of transients in dispersive dielectric media , 1991 .

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

[3]  Shan Zhao,et al.  High order FDTD methods for transverse magnetic modes with dispersive interfaces , 2014, Appl. Math. Comput..

[4]  Doina Cioranescu,et al.  Homogenization of Periodically Varying Coefficients in Electromagnetic Materials , 2006, J. Sci. Comput..

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

[6]  Pingwen Zhang,et al.  Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions , 2004 .

[7]  D. Jiao,et al.  Time-domain finite-element modeling of dispersive media , 2001, IEEE Microwave and Wireless Components Letters.

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

[9]  Dongwei Shi,et al.  Superconvergence analysis of nonconforming mixed finite element methods for time-dependent Maxwell's equations in isotropic cold plasma media , 2013, Appl. Math. Comput..

[10]  Yu-Qiang Zhang,et al.  A Unified FDTD Approach for Electromagnetic Analysis of Dispersive Objects , 2009 .

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

[12]  Jichun Li.  Error analysis of finite element methods for 3-D Maxwell's equations in dispersive media , 2006 .

[13]  Shan Zhao,et al.  Author ' s personal copy High order matched interface and boundary methods for the Helmholtz equation in media with arbitrarily curved interfaces , 2010 .

[14]  Harvey Thomas Banks,et al.  Analysis of stability and dispersion in a finite element method for Debye and Lorentz dispersive media , 2009 .

[15]  Zhimin Zhang,et al.  Error analysis of a discontinuous Galerkin method for Maxwell equations in dispersive media , 2010, J. Comput. Phys..

[16]  Yang Hao,et al.  Finite-Difference Time-Domain Study of Guided Modes in Nano-Plasmonic Waveguides , 2007, IEEE Transactions on Antennas and Propagation.

[17]  Shan Zhao,et al.  A unified discontinuous Galerkin framework for time integration , 2014, Mathematical methods in the applied sciences.

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

[19]  Glenn S. Smith,et al.  A fully three-dimensional simulation of a ground-penetrating radar: FDTD theory compared with experiment , 1996, IEEE Trans. Geosci. Remote. Sens..

[20]  B.D. Van Veen,et al.  An overview of ultra-wideband microwave imaging via space-time beamforming for early-stage breast-cancer detection , 2005, IEEE Antennas and Propagation Magazine.

[21]  R. Luebbers,et al.  Debye Function Expansions of Complex Permittivity Using a Hybrid Particle Swarm-Least Squares Optimization Approach , 2007, IEEE Transactions on Antennas and Propagation.

[22]  Shan Zhao,et al.  High order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources , 2006, J. Comput. Phys..

[23]  Shan Zhao,et al.  Matched interface and boundary (MIB) for the implementation of boundary conditions in high‐order central finite differences , 2009, International journal for numerical methods in engineering.

[24]  C. Balanis.  Advanced Engineering Electromagnetics , 1989 .

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

[26]  Shan Zhao,et al.  High-order FDTD methods for transverse electromagnetic systems in dispersive inhomogeneous media. , 2011, Optics letters.

[27]  Om P. Gandhi,et al.  A frequency-dependent finite-difference time-domain formulation for general dispersive media , 1993 .

[28]  Ahmad Mohammadi,et al.  Dispersive contour-path algorithm for the two-dimensional finite-difference time-domain method. , 2008, Optics express.

[29]  Yunqing Huang,et al.  Interior Penalty Discontinuous Galerkin Method for Maxwell’s Equations in Cold Plasma , 2009, J. Sci. Comput..

[30]  A Taflove,et al.  Direct time integration of Maxwell's equations in linear dispersive media with absorption for scattering and propagation of femtosecond electromagnetic pulses. , 1991, Optics letters.

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

[32]  A. Deinega,et al.  Subpixel smoothing for conductive and dispersive media in the finite-difference time-domain method. , 2007, Optics Letters.