A high-order accurate scheme for Maxwell's equations with a generalized dispersive material model

Abstract A high-order accurate scheme for solving the time-domain Maxwell's equations with a generalized dispersive material model is described. The equations for the electric field are solved in second-order form, and a general dispersion model is treated with the addition of one or more polarization vectors which obey a set of auxiliary differential equations (ADE). Numerical methods are developed for both second-order and fourth-order accuracy in space and time. The equations are discretized using finite-differences, and advanced in time with a single-stage, three-level, space–time scheme which remains stable up to the usual explicit CFL restriction, as proven using mode analysis. Because the equations are treated in their second-order form, there is no need for grid staggering, and instead a collocated grid is used. Composite overlapping grids are used to treat complex geometries with boundary-conforming grids, and a high-order upwind dissipation is added to ensure robust and stable approximations on overlapping grids. Numerical results in two and three space dimensions confirm the accuracy and stability of the new schemes.

[1]  T. Sakurai,et al.  A projection method for generalized eigenvalue problems using numerical integration , 2003 .

[2]  K. Kunz,et al.  A frequency-dependent finite-difference time-domain formulation for transient propagation in plasma , 1991 .

[3]  Jeffrey W. Banks,et al.  Upwind schemes for the wave equation in second-order form , 2012, J. Comput. Phys..

[4]  Raymond J. Luebbers,et al.  FDTD for Nth-order dispersive media , 1992 .

[5]  A. Kildishev,et al.  Fast Eigensolver for Plasmonic Metasurfaces , 2014 .

[6]  Bernard Jecko,et al.  Application of fractional derivatives to the FDTD modeling of pulse propagation in a Cole–Cole dispersive medium , 1996 .

[7]  R. W. Lau,et al.  The dielectric properties of biological tissues: III. Parametric models for the dielectric spectrum of tissues. , 1996, Physics in medicine and biology.

[8]  Ichiro Fukai,et al.  A finite-difference time-domain formulation for transient propagation in dispersive media associated with Cole-Cole's circular arc law , 1990 .

[9]  J. Nédélec.  Mixed finite elements in ℝ3 , 1980 .

[10]  A. Robert.  Dielectric permittivity of concrete between 50 Mhz and 1 Ghz and GPR measurements for building materials evaluation , 1998 .

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

[12]  A. Kildishev,et al.  Optical Dispersion Models for Time-Domain Modeling of Metal-Dielectric Nanostructures , 2011, IEEE Transactions on Magnetics.

[13]  T. Tsiboukis,et al.  HIGHER-ORDER SPATIAL FDTD SCHEMES FOR EM PROPAGATION IN DISPERSIVE MEDIA ∗ , 2005 .

[14]  Jeffrey L. Young,et al.  Toward the construction of a fourth-order difference scheme for transient EM wave simulation: staggered grid approach , 1997 .

[15]  F. Gantmacher,et al.  Applications of the theory of matrices , 1960 .

[16]  K. Cole,et al.  Dispersion and Absorption in Dielectrics I. Alternating Current Characteristics , 1941 .

[17]  Jeffrey W. Banks,et al.  High-order accurate FDTD schemes for dispersive Maxwell's equations in second-order form using recursive convolutions , 2017, J. Comput. Appl. Math..

[18]  S. Havriliak,et al.  A complex plane representation of dielectric and mechanical relaxation processes in some polymers , 1967 .

[19]  M. Centini,et al.  Second- and third-harmonic generation in metal-based structures , 2010 .

[20]  S.-H. Oh,et al.  A hybridizable discontinuous Galerkin method for computing nonlocal electromagnetic effects in three-dimensional metallic nanostructures , 2017, J. Comput. Phys..

[21]  D. Sullivan,et al.  On the dispersion errors related to (FD)/sup 2/TD type schemes , 1995 .

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

[23]  Jeffrey W. Banks,et al.  High-order upwind schemes for the wave equation on overlapping grids: Maxwell's equations in second-order form , 2018, J. Comput. Phys..

[24]  David R. Smith,et al.  Hydrodynamic model for plasmonics: a macroscopic approach to a microscopic problem. , 2013, Chemphyschem : a European journal of chemical physics and physical chemistry.

[25]  Qing Huo Liu,et al.  Higher-Order Numerical Methods for Transient Wave Equations , 2003 .

[26]  Thomas A. Klar,et al.  Numerical modeling of active plasmonic metamaterials , 2011, Optical Systems Design.

[27]  Stephan W. Koch,et al.  Classical theory for second-harmonic generation from metallic nanoparticles. Phys Rev B 79:235109 , 2008, 0807.3575.

[28]  Jian-Ming Jin,et al.  Time-domain finite element modeling of dispersive media , 2001, IEEE Antennas and Propagation Society International Symposium. 2001 Digest. Held in conjunction with: USNC/URSI National Radio Science Meeting (Cat. No.01CH37229).

[29]  Theodoros D. Tsiboukis,et al.  An FDTD Algorithm for Wave Propagation in Dispersive Media using Higher-Order Schemes , 2004 .

[30]  Jeffrey W. Banks,et al.  Numerical methods for solid mechanics on overlapping grids: Linear elasticity , 2010, J. Comput. Phys..

[31]  J. Young.  A higher order FDTD method for EM propagation in a collisionless cold plasma , 1996 .

[32]  S. Havriliak,et al.  A complex plane analysis of α‐dispersions in some polymer systems , 2007 .

[33]  Traianos V. Yioultsis,et al.  Approximation of Grünwald–Letnikov Fractional Derivative for FDTD Modeling of Cole–Cole Media , 2014, IEEE Transactions on Magnetics.

[34]  R. W. Christy,et al.  Optical Constants of the Noble Metals , 1972 .

[35]  Youngjoo Chung,et al.  PLRC AND ADE IMPLEMENTATIONS OF DRUDE- CRITICAL POINT DISPERSIVE MODEL FOR THE FDTD METHOD , 2013 .

[36]  R. B. Standler,et al.  A frequency-dependent finite-difference time-domain formulation for dispersive materials , 1990 .

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

[38]  Ichiro Fukai,et al.  A treatment by the FD‐TD method of the dispersive characteristics associated with electronic polarization , 1990 .

[39]  Alexander V. Kildishev,et al.  Efficient time-domain model of the graphene dielectric function , 2013, Optics & Photonics - NanoScience + Engineering.

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

[41]  Diego Caratelli,et al.  Fractional Derivative Based FDTD Modeling of Transient Wave Propagation in Havriliak–Negami Media , 2014, IEEE Transactions on Microwave Theory and Techniques.

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

[43]  R. Cole,et al.  Dielectric Relaxation in Glycerol, Propylene Glycol, and n‐Propanol , 1951 .

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

[45]  S. D. Gedney,et al.  A Discontinuous Galerkin Finite Element Time-Domain Method Modeling of Dispersive Media , 2012, IEEE Transactions on Antennas and Propagation.

[46]  Norinobu Yoshida,et al.  A Treatment by the Finite-Difference Time-Domain Method of the Dispersive Characteristics Associated with Orientation Polarization , 1990 .

[47]  William D. Henshaw,et al.  A High-Order Accurate Parallel Solver for Maxwell's Equations on Overlapping Grids , 2005, SIAM J. Sci. Comput..

[48]  John E. Sipe,et al.  Analysis of second-harmonic generation at metal surfaces , 1980 .

[49]  Yong Zeng,et al.  Generalization of the FDTD algorithm for simulations of hydrodynamic nonlinear Drude model , 2010, J. Comput. Phys..

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

[51]  Alexey V. Krasavin,et al.  Nonperturbative hydrodynamic model for multiple harmonics generation in metallic nanostructures , 2014, 1405.4903.

[52]  R. J. Luebbers,et al.  Piecewise linear recursive convolution for dispersive media using FDTD , 1996 .

[53]  P. Etchegoin,et al.  An analytic model for the optical properties of gold. , 2006, The Journal of chemical physics.

[54]  W. Henshaw,et al.  Composite overlapping meshes for the solution of partial differential equations , 1990 .

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