Finite-Difference Time-Domain Solution of Maxwell's Equations

For over 100 years after the publication of Maxwell's equations in 1865, essentially all solution techniques for electromagnetic fields and waves were based on Fourier-domain concepts, assuming a priori a time-harmonic (sinusoidal steady-state) field variation and possibly the existence of a particular Green's function or a set of spatial modes. In 1966, Kane Yee's seminal paper introduced a complete paradigm shift in how to solve Maxwell's equations, reporting a field evolution-in-time technique that subsequently evolved into the finite-difference time-domain (FDTD) method. In the decades since the publication of Yee's paper, there has been an explosion of interest in FDTD and related grid-based time-marched solutions of Maxwell's equations among scientists and engineers. During this period, FDTD modeling has evolved to an advanced stage enabling large-scale simulations of full-wave time-domain electromagnetic wave interactions with volumetrically complex structures over large frequency ranges, spatial scales, and timescales. Currently, FDTD modeling spans the electromagnetic spectrum from ultralow frequencies to visible light. FDTD modeling is routinely conducted as an invaluable virtual laboratory bench in scientific inquiry and exploration in electrodynamics; as an integral part of the electromagnetic engineering design and optimization process; and as a powerful forward solver in imaging and sensing inverse problems. This article reviews the technical basis of the key features of FDTD solution techniques for Maxwell's equations and provides 18 modeling examples spanning the electromagnetic spectrum to illustrate the power, flexibility, and robust nature of FDTD computational electrodynamics simulations. Keywords: computational electrodynamics; finite-difference time-domain; FDTD; Maxwell's equations

[1]  S. Gedney.  An anisotropic perfectly matched layer-absorbing medium for the truncation of FDTD lattices , 1996 .

[2]  T. Krauss,et al.  Flatband slow light in photonic crystals featuring spatial pulse compression and terahertz bandwidth. , 2007, Optics express.

[3]  J. Simpson,et al.  An Efficient 3-D FDTD Model of Electromagnetic Wave Propagation in Magnetized Plasma , 2015, IEEE Transactions on Antennas and Propagation.

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

[5]  Shuiping Luo,et al.  Generalization of the Finite-Difference-Based Time-Domain Methods Using the Method of Moments , 2006, IEEE Transactions on Antennas and Propagation.

[6]  A. Taflove,et al.  Three-dimensional FDTD modeling of impulsive ELF propagation about the earth-sphere , 2004, IEEE Transactions on Antennas and Propagation.

[7]  Stephen D. Gedney,et al.  Convolution PML (CPML): An efficient FDTD implementation of the CFS–PML for arbitrary media , 2000 .

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

[9]  R. Grondin,et al.  Electromagnetic and transport considerations in subpicosecond photoconductive switch modeling , 1990 .

[10]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .

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

[12]  S. Hagness,et al.  Terahertz conductivity of doped silicon calculated using the ensemble Monte Carlo/finite-difference time-domain simulation technique , 2010 .

[13]  Carey M. Rappaport,et al.  Perfectly matched absorbing boundary conditions based on anisotropic lossy mapping of space , 1995 .

[14]  Richard Holland,et al.  Total-Field versus Scattered-Field Finite-Difference Codes: A Comparative Assessment , 1983, IEEE Transactions on Nuclear Science.

[15]  Niels Kuster,et al.  Suitability of FDTD-based TCAD tools RF design of mobile phones , 2003 .

[16]  Seng-Tiong Ho,et al.  FDTD microcavity simulations: design and experimental realization of waveguide-coupled single-mode ring and whispering-gallery-mode disk resonators , 1997 .

[17]  Allen Taflove,et al.  Application of the Finite-Difference Time-Domain Method to Sinusoidal Steady-State Electromagnetic-Penetration Problems , 1980, IEEE Transactions on Electromagnetic Compatibility.

[18]  Qing Huo Liu,et al.  Direct $Z$-Transform Implementation of the CFS-PML Based on Memory-Minimized Method , 2015, IEEE Transactions on Microwave Theory and Techniques.

[19]  J. Yamauchi,et al.  Efficient implicit FDTD algorithm based on locally one-dimensional scheme , 2005 .

[20]  Zhizhang Chen,et al.  Numerical dispersion analysis of the unconditionally stable 3-D ADI-FDTD method , 2001 .

[21]  Jin-Fa Lee,et al.  A perfectly matched anisotropic absorber for use as an absorbing boundary condition , 1995 .

[22]  K. Kunz,et al.  Finite-Difference Analysis of EMP Coupling to Lossy Dielectric Structures , 1980, IEEE Transactions on Electromagnetic Compatibility.

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

[24]  Tae-Woo Lee,et al.  On the accuracy of the ADI-FDTD method , 2002, IEEE Antennas and Wireless Propagation Letters.

[25]  Chen Bin,et al.  An Improved Paralleling-in-Order Solving Scheme for AH-FDTD Method Using Eigenvalue Transformation , 2015, IEEE Transactions on Antennas and Propagation.

[26]  T. Weiland.  Time Domain Electromagnetic Field Computation with Finite Difference Methods , 1996 .

[27]  Er-Ping Li,et al.  Development of the Three-Dimensional Unconditionally Stable LOD-FDTD Method , 2008, IEEE Transactions on Antennas and Propagation.

[28]  A Taflove,et al.  Interferometric spectroscopy of scattered light can quantify the statistics of subdiffractional refractive-index fluctuations. , 2013, Physical review letters.

[29]  Raj Mittra,et al.  Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers , 1996 .

[30]  Dennis M. Sullivan,et al.  Electromagnetic Simulation Using the FDTD Method: Sullivan/Electromagnetic Simulation Using the FDTD Method , 2013 .

[31]  R. Holland.  THREDE: A Free-Field EMP Coupling and Scattering Code , 1977, IEEE Transactions on Nuclear Science.

[32]  S. Hagness,et al.  Global modeling of carrier-field dynamics in semiconductors using EMC–FDTD , 2009 .

[33]  T. Namiki,et al.  A new FDTD algorithm based on alternating-direction implicit method , 1999 .

[34]  S. Hagness,et al.  Clustered impurities and carrier transport in supported graphene , 2014, 1405.4938.

[35]  S. Hagness,et al.  Multiphysics simulation of high-frequency carrier dynamics in conductive materials , 2011 .

[36]  John B. Schneider,et al.  Divergent fields, charge, and capacitance in FDTD simulations , 1998 .

[37]  S. Gedney,et al.  On the long-time behavior of unsplit perfectly matched layers , 2004, IEEE Transactions on Antennas and Propagation.

[38]  Soon-Hong Kwon,et al.  Electrically Driven Single-Cell Photonic Crystal Laser , 2004, Science.

[39]  A. Weisshaar,et al.  CAD-oriented equivalent-circuit modeling of on-chip interconnects on lossy silicon substrate , 2000 .

[40]  Yen Liu,et al.  Fourier Analysis of Numerical Algorithms for the Maxwell Equations , 1993 .

[41]  Vadim Backman,et al.  High-speed spectral nanocytology for early cancer screening , 2013, Journal of biomedical optics.

[42]  Ali Akbarzadeh-Sharbaf,et al.  On the Development of Nonoverlapping and Stable Hybrid FETD-FDTD Formulations , 2014, IEEE Transactions on Antennas and Propagation.

[43]  S. Hagness,et al.  Terahertz-frequency electronic transport in graphene , 2014 .

[44]  A. Taflove,et al.  Reinventing electromagnetics: emerging applications for FD-TD computation , 1995 .

[45]  Allen Taflove,et al.  A Novel Method to Analyze Electromagnetic Scattering of Complex Objects , 1982, IEEE Transactions on Electromagnetic Compatibility.

[46]  Peter B. Johns,et al.  Numerical solution of 2-dimensional scattering problems using a transmission-line matrix , 1971 .

[47]  Steven G. Johnson,et al.  Casimir forces in the time domain: Applications , 2009, 0906.5170.

[48]  Allen Taflove,et al.  The determination of the effective radius of a filamentary source in the FDTD mesh , 2000 .

[49]  Mark A. Ratner,et al.  Classical Electrodynamics Coupled to Quantum Mechanics for Calculation of Molecular Optical Properties: a RT-TDDFT/FDTD Approach , 2010 .

[50]  Tapan K. Sarkar,et al.  An unconditionally stable scheme for the finite-difference time-domain method , 2003 .

[51]  Dennis M. Sullivan,et al.  Z-transform theory and the FDTD method , 1996 .

[52]  A. Taflove,et al.  Computational and experimental study of a microwave electromagnetic bandgap structure with waveguiding defect for potential use as a bandpass wireless interconnect , 2004, IEEE Microwave and Wireless Components Letters.

[53]  R. Luebbers,et al.  A finite-difference time-domain near zone to far zone transformation (electromagnetic scattering) , 1991 .

[54]  M. Potter,et al.  FDTD Discrete Planewave (FDTD-DPW) Formulation for a Perfectly Matched Source in TFSF Simulations , 2010, IEEE Transactions on Antennas and Propagation.

[55]  Bin Chen,et al.  A New Unconditionally Stable Scheme for FDTD Method Using Associated Hermite Orthogonal Functions , 2014, IEEE Transactions on Antennas and Propagation.

[56]  S. Hagness,et al.  EMC/FDTD/MD simulation of carrier transport and electrodynamics in two-dimensional electron systems , 2013, 1405.4968.

[57]  D. Katz,et al.  Validation and extension to three dimensions of the Berenger PML absorbing boundary condition for FD-TD meshes , 1994, IEEE Microwave and Guided Wave Letters.