View source for Finite-difference time-domain method

{{Short description|Numerical analysis technique}}
[[File:Yee cell.png|thumb|250px|In finite-difference time-domain method, "Yee lattice" is used to discretize [[Maxwell's equations]] in space. This scheme involves the placement of [[Electric field|electric]] and [[magnetic fields]] on a staggered grid.]]
'''Finite-difference time-domain''' ('''FDTD''') or '''Yee's method''' (named after the Chinese American applied mathematician [[Kane S. Yee]], born 1934) is a [[numerical analysis]] technique used for modeling [[computational electrodynamics]] (finding approximate solutions to the associated system of [[differential equation]]s).  Since it is a [[time domain|time-domain]] method, FDTD solutions can cover a wide [[frequency]] range with a single [[computer simulation|simulation]] run, and treat nonlinear material properties in a natural way.

The FDTD method belongs in the general class of [[Discretization|grid]]-based differential numerical modeling methods ([[finite difference methods]]). The time-dependent [[Maxwell's equations]] (in [[Partial differential equation|partial differential]] form) are discretized using [[central difference|central-difference]] approximations to the space and time [[partial derivative]]s. The resulting [[finite difference method|finite-difference]] equations are solved in either software or hardware in a [[leapfrog integration|leapfrog]] manner:  the [[electric field]] [[vector component]]s in a volume of space are solved at a given instant in time; then the [[magnetic field]] vector components in the same spatial volume are solved at the next instant in time; and the process is repeated over and over again until the desired transient or steady-state electromagnetic field behavior is fully evolved.

== History ==
Finite difference schemes for time-dependent [[partial differential equation]]s (PDEs) have been employed for many years in [[computational fluid dynamics]] problems,<ref name="vonneumann49" /> including the idea of using centered finite difference operators on staggered grids in space and time to achieve second-order accuracy.<ref name="vonneumann49" />
The novelty of Kane Yee's FDTD scheme, presented in his seminal 1966 paper,<ref name="yee66" /> was to apply centered finite difference operators on staggered grids in space and time for each electric and magnetic vector field component in Maxwell's curl equations.
The descriptor "Finite-difference time-domain" and its corresponding "FDTD" acronym were originated by [[Allen Taflove]] in 1980.<ref name="taflove80" />
Since about 1990, FDTD techniques have emerged as primary means to computationally model many scientific and engineering problems dealing with [[electromagnetic wave]] interactions with material structures. Current FDTD modeling applications range from near-[[Direct current|DC]] (ultralow-frequency [[geophysics]] involving the entire Earth-[[ionosphere]] waveguide) through [[microwaves]] (radar signature technology, [[Antenna (radio)|antennas]], wireless communications devices, digital interconnects, biomedical imaging/treatment) to [[visible light]] ([[photonic crystal]]s, nano[[plasmon]]ics, [[soliton]]s, and [[biophotonics]]).<ref name="taflove05" />  In 2006, an estimated 2,000 FDTD-related publications appeared in the science and engineering literature (see [[#Popularity|Popularity]]). As of 2013, there are at least 25 commercial/proprietary FDTD software vendors; 13 free-software/[[Open source|open-source]]-software FDTD projects; and 2 freeware/closed-source FDTD projects, some not for commercial use (see [[#External links|External links]]).

=== Development of FDTD and Maxwell's equations===
{{Copypaste|section|date=January 2021}}<!-- Contents of the chronology, despite being referenced with the original articles, appears to be largely taken in verbatim from Taflove and Hagness's book. (Chapter 1) -->
An appreciation of the basis, technical development, and possible future of FDTD numerical techniques for Maxwell's equations can be developed by first considering their history.  The following lists some of the key publications in this area.

{| class="wikitable" width="90%" style="text-align:left"
! colspan="2" | Partial chronology of FDTD techniques and applications for Maxwell's equations.<ref>Adapted with permission from Taflove and Hagness (2005).</ref>
|-
! width="10%"  | year
! width="80%"  | event
|-
| 1928 || Courant, Friedrichs, and Lewy (CFL) publish seminal paper with the discovery of conditional stability of explicit time-dependent finite difference schemes, as well as the classic FD scheme for solving second-order wave equation in 1-D and 2-D.<ref name="courant1928" />
|-
| 1950 || First appearance of von Neumann's method of stability analysis for implicit/explicit time-dependent finite difference methods.<ref name="obrien1950" />
|-
| 1966 || Yee described the FDTD numerical technique for solving Maxwell's curl equations on grids staggered in space and time.<ref name="yee66" />
|-
| 1969 || Lam reported the correct numerical CFL stability condition for Yee's algorithm by employing von Neumann stability analysis.<ref name="lam69" />
|-
| 1975 || Taflove and Brodwin reported the first sinusoidal steady-state FDTD solutions of two- and three-dimensional electromagnetic wave interactions with material structures;<ref name="taflove75a" /> and the first bioelectromagnetics models.<ref name="taflove75b" />
|-
| 1977 || Holland and Kunz & Lee applied Yee's algorithm to EMP problems.<ref name="holland77" /><ref name="kunz77" />
|-
| 1980 || Taflove coined the FDTD acronym and published the first validated FDTD models of sinusoidal steady-state electromagnetic wave penetration into a three-dimensional metal cavity.<ref name="taflove80" />
|- 
| 1981 || Mur published the first numerically stable, second-order accurate, absorbing boundary condition (ABC) for Yee's grid.<ref name="mur81" />
|-
| 1982–83 || Taflove and Umashankar developed the first FDTD electromagnetic wave scattering models computing sinusoidal steady-state near-fields, far-fields, and radar cross-section for two- and three-dimensional structures.<ref name="umashankar82" /><ref name="taflove83" />
|-
| 1984 || Liao ''et al'' reported an improved ABC based upon space-time extrapolation of the field adjacent to the outer grid boundary.<ref name="liao84" />
|-
| 1985 || Gwarek introduced the lumped equivalent circuit formulation of FDTD.<ref name="gwarek85" />
|-
| 1986 || Choi and Hoefer published the first FDTD simulation of waveguide structures.<ref name="choi86" />
|-
| 1987–88 || Kriegsmann ''et al'' and Moore ''et al'' published the first articles on ABC theory in ''IEEE  Transactions on Antennas and Propagation''.<ref name="kriegsmann87" /><ref name="moore88" />
|-
| 1987–88, 1992  || Contour-path subcell techniques were introduced by Umashankar ''et al'' to permit FDTD modeling of thin wires and wire bundles,<ref name="umashankar87" /> by Taflove ''et al'' to model penetration through cracks in conducting screens,<ref name="taflove88" /> and by Jurgens ''et al'' to conformally model the surface of a smoothly curved scatterer.<ref name="jurgens92" />
|-
| 1988 || Sullivan ''et al'' published the first 3-D FDTD model of sinusoidal steady-state electromagnetic wave absorption by a complete human body.<ref name="sullivan88" />
|-
| 1988 || FDTD modeling of microstrips was introduced by Zhang ''et al''.<ref name="zhang88" />
|-
| 1990–91 || FDTD modeling of frequency-dependent dielectric permittivity was introduced by Kashiwa and Fukai,<ref name="kashiwa90" /> Luebbers ''et al'',<ref name="luebbers90" /> and Joseph ''et al''.<ref name="joseph91" />
|-
| 1990–91 || FDTD modeling of antennas was introduced by Maloney ''et al'',<ref name="maloney90" /> Katz ''et al'',<ref name="katz91" /> and Tirkas and Balanis.<ref name="tirkas91" />
|-
| 1990 || FDTD modeling of picosecond optoelectronic switches was introduced by Sano and Shibata,<ref name="sano90" /> and El-Ghazaly ''et al''.<ref name="el-ghazaly90" />
|-
| 1992–94 || FDTD modeling of the propagation of optical pulses in nonlinear dispersive media was introduced, including the first temporal solitons in one dimension by Goorjian and Taflove;<ref name="goorjian92" />  beam self-focusing by Ziolkowski and Judkins;<ref name="ziolkowski93" /> the first temporal solitons in two dimensions by Joseph ''et al'';<ref name="joseph93" />  and the first spatial solitons in two dimensions by Joseph and Taflove.<ref name="joseph94" />
|-
| 1992 || FDTD modeling of lumped electronic circuit elements was introduced by Sui ''et al''.<ref name="sui92" />
|-
| 1993 || Toland ''et al'' published the first FDTD models of gain devices (tunnel diodes and Gunn diodes) exciting cavities and antennas.<ref name="toland93" />
|-
| 1993 || Aoyagi ''et al'' present a hybrid Yee algorithm/scalar-wave equation and demonstrate equivalence of Yee scheme to finite difference scheme for [[electromagnetic wave equation]].<ref name="aoyagi93" />
|-
| 1994 || Thomas ''et al'' introduced a Norton's equivalent circuit for the FDTD space lattice, which permits the SPICE circuit analysis tool to implement accurate subgrid models of nonlinear electronic components or complete circuits embedded within the lattice.<ref name="thomas94" />
|-
| 1994 || Berenger introduced the highly effective, perfectly matched layer (PML) ABC for two-dimensional FDTD grids,<ref name="berenger94" /> which was extended to non-orthogonal meshes by Navarro ''et al'',<ref name="navarro94" /> and three dimensions by Katz ''et al'',<ref name="katz94" />  and to dispersive waveguide terminations by Reuter ''et al''.<ref name="reuter94" />
|-
| 1994 || Chew and Weedon introduced the coordinate stretching PML that is easily extended to three dimensions, other coordinate systems and other physical equations.<ref name="chewweedon94" />
|-
| 1995–96 || Sacks ''et al'' and Gedney introduced a physically realizable, uniaxial perfectly matched layer (UPML) ABC.<ref name="gedney96" /><ref name="sacks95" />
|-
| 1997 || Liu introduced the pseudospectral time-domain (PSTD) method, which permits extremely coarse spatial sampling of the electromagnetic field at the Nyquist limit.<ref name="liu97" />
|-
| 1997 || Ramahi introduced the complementary operators method (COM) to implement highly effective analytical ABCs.<ref name="ramahi97" />
|-
| 1998 || Maloney and Kesler introduced several novel means to analyze periodic structures in the FDTD space lattice.<ref name="maloney98" />
|-
| 1998 || Nagra and York introduced a hybrid FDTD-quantum mechanics model of electromagnetic wave interactions with materials having electrons transitioning between multiple energy levels.<ref name="nagra98" />
|-
| 1998 || Hagness ''et al'' introduced FDTD modeling of the detection of breast cancer using ultrawideband radar techniques.<ref name="hagness98" />
|-
| 1999 || Schneider and Wagner introduced a comprehensive analysis of FDTD grid dispersion based upon complex wavenumbers.<ref name="schneider99" />
|-
| 2000–01 || Zheng, Chen, and Zhang introduced the first three-dimensional alternating-direction implicit (ADI) FDTD algorithm with provable unconditional numerical stability.<ref name="zhen00" /><ref name="zheng01" />
|-
| 2000 || Roden and Gedney introduced the advanced convolutional PML (CPML) ABC.<ref name="roden00" />
|-
| 2000 || Rylander and Bondeson introduced a provably stable FDTD - finite-element time-domain hybrid technique.<ref name="rylander00" />
|-
| 2002 || Hayakawa ''et al'' and Simpson and Taflove independently introduced FDTD modeling of the global Earth-ionosphere waveguide for extremely low-frequency geophysical phenomena.<ref name="hayakawa02" /><ref name="simpson02" />
|-
| 2003 || DeRaedt introduced the unconditionally stable, “one-step” FDTD technique.<ref name="de_raedt03" />
|-
| 2004 || Soriano and Navarro derived the stability condition for Quantum FDTD technique.<ref name="SorianoNavarro2004" />
|-
| 2008 || Ahmed, Chua, Li and Chen introduced the three-dimensional locally one-dimensional (LOD)FDTD method and proved unconditional numerical stability.<ref name="Ahmed2008" />
|-
| 2008 || Taniguchi, Baba, Nagaoka and Ametani introduced a Thin Wire Representation for FDTD Computations for conductive media<ref name="baba08" />
|-
| 2009 || Oliveira and Sobrinho applied the FDTD method for simulating lightning strokes in a power substation<ref name="oliveira09" />
|-
| 2012 || Moxley ''et al'' developed a generalized finite-difference time-domain quantum method for the N-body interacting Hamiltonian.<ref name="Moxley2012" />
|-
| 2013 || Moxley ''et al'' developed a generalized finite-difference time-domain scheme for solving nonlinear Schrödinger equations.<ref name="Moxley2013" />
|-
| 2014 || Moxley ''et al'' developed an implicit generalized finite-difference time-domain scheme for solving nonlinear Schrödinger equations.<ref name="Moxley2014" />
|-
| 2021 || Oliveira and Paiva developed the Least Squares Finite-Difference Time-Domain method (LS-FDTD) for using time steps beyond FDTD CFL limit.<ref name="oliveira2021" />
|-
|}

== FDTD models and methods ==
When [[Maxwell's differential equations]] are examined, it can be seen that the change in the E-field in time (the time derivative) is dependent on the change in the H-field across space (the [[curl (mathematics)|curl]]).  This results in the basic FDTD time-stepping relation that, at any point in space, the updated value of the E-field in time is dependent on the stored value of the E-field and the numerical curl of the local distribution of the H-field in space.<ref name="yee66" />

The H-field is time-stepped in a similar manner. At any point in space, the updated value of the H-field in time is dependent on the stored value of the H-field and the numerical curl of the local distribution of the E-field in space. Iterating the E-field and H-field updates results in a marching-in-time process wherein sampled-data analogs of the continuous electromagnetic waves under consideration propagate in a numerical grid stored in the computer memory.

[[Image:FDTD Yee grid 2d-3d.svg|thumb|right|450px|Illustration of a standard Cartesian Yee cell used for FDTD, about which electric and magnetic field vector components are distributed.<ref name="yee66" /> Visualized as a cubic [[voxel]], the electric field components form the edges of the cube, and the magnetic field components form the normals to the faces of the cube. A three-dimensional space lattice consists of a multiplicity of such Yee cells. An electromagnetic wave interaction structure is mapped into the space lattice by assigning appropriate values of permittivity to each electric field component, and permeability to each magnetic field component.]]
This description holds true for 1-D, 2-D, and 3-D FDTD techniques. When multiple dimensions are considered, calculating the numerical curl can become complicated. Kane Yee's seminal 1966 paper proposed spatially staggering the vector components of the E-field and H-field about rectangular unit cells of a Cartesian computational grid so that each E-field vector component is located midway between a pair of H-field vector components, and conversely.<ref name="yee66" />  This scheme, now known as a '''Yee lattice''', has proven to be very robust, and remains at the core of many current FDTD software constructs.

Furthermore, Yee proposed a leapfrog scheme for marching in time wherein the E-field and H-field updates are staggered so that E-field updates are conducted midway during each time-step between successive H-field updates, and conversely.<ref name="yee66" /> On the plus side, this explicit time-stepping scheme avoids the need to solve simultaneous equations, and furthermore yields dissipation-free numerical wave propagation. On the minus side, this scheme mandates an upper bound on the time-step to ensure numerical stability.<ref name="taflove75a" /> As a result, certain classes of simulations can require many thousands of time-steps for completion.

=== Using the FDTD method ===
To implement an FDTD solution of Maxwell's equations, a computational domain must first be established. The computational domain is simply the physical region over which the simulation will be performed. The E and H fields are determined at every point in space within that computational domain. The material of each cell within the computational domain must be specified. Typically, the material is either free-space (air), [[metal]], or [[dielectric]].  Any material can be used as long as the [[Permeability (electromagnetism)|permeability]], [[permittivity]], and [[electrical conductivity|conductivity]] are specified.

The permittivity of dispersive materials in tabular form cannot be directly substituted into the FDTD scheme.
Instead, it can be approximated using multiple Debye, Drude, Lorentz or critical point terms.
This approximation can be obtained using open fitting programs<ref name="fitting" /> and does not necessarily have physical meaning.

Once the computational domain and the grid materials are established, a source is specified. The source can be current on a wire,  applied electric field or impinging plane wave.
In the last case FDTD can be used to simulate light scattering from arbitrary shaped objects, planar periodic structures at various incident angles,<ref name="obl_it" /><ref name="obl_sfdtd" /> and photonic band structure of infinite periodic structures.<ref name="TMatrix" /><ref name="Hao" />

Since the E and H fields are determined directly, the output of the simulation is usually the E or H field at a point or a series of points within the computational domain. The simulation evolves the E and H fields forward in time.

Processing may be done on the E and H fields returned by the simulation. Data processing may also occur while the simulation is ongoing.

While the FDTD technique computes electromagnetic fields within a compact spatial region, scattered and/or radiated far fields can be obtained via near-to-far-field transformations.<ref name="umashankar82" />

=== Strengths of FDTD modeling ===
Every modeling technique has strengths and weaknesses, and the FDTD method is no different.

* FDTD is a versatile modeling technique used to solve Maxwell's equations. It is intuitive, so users can easily understand how to use it and know what to expect from a given model.
* FDTD is a time-domain technique, and when a broadband pulse (such as a Gaussian pulse) is used as the source, then the response of the system over a wide range of frequencies can be obtained with a single simulation. This is useful in applications where resonant frequencies are not exactly known, or anytime that a broadband result is desired.
* Since FDTD calculates the E and H fields everywhere in the computational domain as they evolve in time, it lends itself to providing animated displays of the electromagnetic field movement through the model. This type of display is useful in understanding what is going on in the model, and to help ensure that the model is working correctly.
* The FDTD technique allows the user to specify the material at all points within the computational domain. A wide variety of linear and nonlinear dielectric and magnetic materials can be naturally and easily modeled.
* FDTD allows the effects of apertures to be determined directly. Shielding effects can be found, and the fields both inside and outside a structure can be found directly or indirectly.
* FDTD uses the E and H fields directly. Since most EMI/EMC modeling applications are interested in the E and H fields, it is convenient that no conversions must be made after the simulation has run to get these values.

=== Weaknesses of FDTD modeling===
[[File:Numerical dispersion of a pulse signal in 1D FDTD.ogg|thumb|right|290px|[[Numerical dispersion]] of a square pulse signal in a simple one-dimensional FDTD scheme. [[Ringing artifacts]] around the edges of the pulse are heavily accentuated ([[Gibbs phenomenon]]) and the signal distorts as it propagates, even in the absence of a [[Dispersion (optics)|dispersive medium]]. This artifact is a direct result of the discretization scheme.<ref name="taflove05"/>]]
* Since FDTD requires that the entire computational domain be gridded, and the grid spatial discretization must be sufficiently fine to resolve both the smallest electromagnetic wavelength and the smallest geometrical feature in the model, very large computational domains can be developed, which results in very long solution times. Models with long, thin features, (like wires) are difficult to model in FDTD because of the excessively large computational domain required. Methods such as [[eigenmode expansion]] can offer a more efficient alternative as they do not require a fine grid along the z-direction.<ref name="phot_cad" />
* There is no way to determine unique values for permittivity and permeability at a material interface.
* Space and time steps must satisfy the [[Courant–Friedrichs–Lewy condition|CFL condition]], or the [[leapfrog integration]] used to solve the partial differential equation is likely to become unstable.
* FDTD finds the E/H fields directly everywhere in the computational domain. If the field values at some distance are desired, it is likely that this distance will force the computational domain to be excessively large. Far-field extensions are available for FDTD, but require some amount of postprocessing.<ref name="taflove05" />
* Since FDTD simulations calculate the E and H fields at all points within the computational domain, the computational domain must be finite to permit its residence in the computer memory. In many cases this is achieved by inserting artificial boundaries into the simulation space. Care must be taken to minimize errors introduced by such boundaries. There are a number of available highly effective absorbing boundary conditions (ABCs) to simulate an infinite unbounded computational domain.<ref name="taflove05" />  Most modern FDTD implementations instead use a special absorbing "material", called a [[perfectly matched layer]] (PML) to implement absorbing boundaries.<ref name="berenger94" /><ref name="gedney96" />
* Because FDTD is solved by propagating the fields forward in the time domain, the electromagnetic time response of the medium must be modeled explicitly.  For an arbitrary response, this involves a computationally expensive time convolution, although in most cases the time response of the medium (or [[Dispersion (optics)]]) can be adequately and simply modeled using either the recursive convolution (RC) technique, the auxiliary differential equation (ADE) technique, or the Z-transform technique. An alternative way of solving [[Maxwell's equations]] that can treat arbitrary dispersion easily is the [[Computational electrodynamics#Pseudo-spectral spatial domain .28PSSD.29|pseudo-spectral spatial domain (PSSD)]], which instead propagates the fields forward in space.

===Grid truncation techniques===
The most commonly used grid truncation techniques for open-region FDTD modeling problems are the Mur absorbing boundary condition (ABC),<ref name="mur81" /> the Liao ABC,<ref name="liao84" /> and various [[perfectly matched layer]] (PML) formulations.<ref name="taflove05" /><ref name="navarro94" /><ref name="berenger94" /><ref name="gedney96" /> The Mur and Liao techniques are simpler than PML. However, PML (which is technically an absorbing region rather than a boundary condition ''per se'') can provide orders-of-magnitude lower reflections. The PML concept was introduced by J.-P. Berenger in a seminal 1994 paper in the Journal of Computational Physics.<ref name="berenger94" />  Since 1994, Berenger's original split-field implementation has been modified and extended to the uniaxial PML (UPML), the convolutional PML (CPML), and the higher-order PML. The latter two PML formulations have increased ability to absorb evanescent waves, and therefore can in principle be placed closer to a simulated scattering or radiating structure than Berenger's original formulation.

To reduce undesired numerical reflection from the PML additional back absorbing layers technique can be used.<ref name="back_pml" />

== Popularity ==
{{Original research|section|date=August 2013}}
{{Copypaste|section|date=January 2021}}<!-- Contents of the section, despite being referenced with the original source, the content of the section appears to be largely taken in verbatim from Taflove and Hagness's book. (Chapter 1) -->

<!-- The following text is from Computational Electrodynamics:  The
Finite-Difference Time-Domain Method, by Taflove.  Taflove and Artech House granted permission to use this text under the CC-A-SA-3.0 license; the permission was sent by the publisher to permissions-en@wikimedia.org on 21 February 2012. -->

Notwithstanding both the general increase in academic publication
throughput during the same period and the overall expansion of interest
in all Computational electromagnetics (CEM)  techniques, there are
seven primary reasons for the tremendous expansion of interest in FDTD
computational solution approaches for Maxwell's equations:

# FDTD does not require a matrix inversion.  Being a fully explicit computation, FDTD avoids the difficulties with matrix inversions that limit the size of frequency-domain integral-equation and finite-element electromagnetics models to generally fewer than 10<sup>9</sup> electromagnetic field unknowns.<ref name="taflove05" /> FDTD models with as many as 10<sup>9</sup> field unknowns have been run; there is no intrinsic upper bound to this number.<ref name="taflove05" />
# FDTD is accurate and robust.  The sources of error in FDTD calculations are well understood, and can be bounded to permit accurate models for a very large variety of electromagnetic wave interaction problems.<ref name="taflove05" />
# FDTD treats impulsive behavior naturally.  Being a time-domain technique, FDTD directly calculates the impulse response of an electromagnetic system.  Therefore, a single FDTD simulation can provide either ultrawideband temporal waveforms or the sinusoidal steady-state response at any frequency within the excitation spectrum.<ref name="taflove05" />
# FDTD treats nonlinear behavior naturally.  Being a time-domain technique, FDTD directly calculates the nonlinear response of an electromagnetic system.  This allows natural hybriding of FDTD with sets of auxiliary differential equations that describe nonlinearities from either the classical or semi-classical standpoint.<ref name="taflove05" />  One research frontier is the development of hybrid algorithms which join FDTD classical electrodynamics models with phenomena arising from quantum electrodynamics, especially vacuum fluctuations, such as the [[Casimir effect]].<ref name="taflove05" /><ref>S. G. Johnson, "[https://arxiv.org/abs/arXiv:1007.0966 Numerical methods for computing Casimir interactions]," in Casimir Physics (D. Dalvit, [[P. Milonni]], D. Roberts, and F. da Rosa, eds.), vol. 834 of ''Lecture Notes in Physics'', ch. 6, pp. 175–218, Berlin: Springer, June 2011.</ref>
# FDTD is a systematic approach. With FDTD, specifying a new structure to be modeled is reduced to a problem of mesh generation rather than the potentially complex reformulation of an integral equation.  For example, FDTD requires no calculation of structure-dependent Green functions.<ref name="taflove05" />
# Parallel-processing computer architectures have come to dominate supercomputing.  FDTD scales with high efficiency on parallel-processing CPU-based computers, and extremely well on recently developed GPU-based accelerator technology.<ref name="taflove05" />
# Computer visualization capabilities are increasing rapidly. While this trend positively influences all numerical techniques, it is of particular advantage to FDTD methods, which generate time-marched arrays of field quantities suitable for use in color videos to illustrate the field dynamics.<ref name="taflove05" />
 
Taflove has argued that these factors combine to suggest that FDTD will remain one of
the dominant computational electrodynamics techniques (as well as potentially other [[multi-physics|multiphysics]] problems).<ref name="taflove05" />

==See also==
* [[Computational electromagnetics]]
* [[Eigenmode expansion]]
* [[Beam propagation method]]
* [[Finite-difference frequency-domain]]
* [[Finite element method]]
* [[Scattering-matrix method]]
* [[Discrete dipole approximation]]

==References==
{{Reflist|3|refs=
<ref name="courant1928">
{{cite journal
  |author1=Richard Courant |author2=Kurt Otto Friedrichs |author3=Hans Lewy | title = Über die partiellen Differenzengleichungen der mathematischen Physik
  | journal = [[Mathematische Annalen]]
  | volume = 100
  | issue = 1
  | pages = 32–74
  | year = 1928
  | language = de
  | url = http://resolver.sub.uni-goettingen.de/purl?GDZPPN002272636
  | doi = 10.1007/BF01448839
  | jfm = 54.0486.01
  | mr = 1512478
 |bibcode = 1928MatAn.100...32C |s2cid=120760331 }}</ref>

<ref name="obrien1950">
{{cite journal
  | author = G. G. O’Brien, M. A Hyman, and S. Kaplan
  | title = A study of the numerical solution of partial differential equations
  | journal = Journal of Mathematical Physics
  | volume = 29
  | issue = 1
  | pages = 223–251
  | year = 1950
  | mr = 0040805
 | doi = 10.1002/sapm1950291223
  }}</ref>

<ref name="Moxley2012">
{{cite journal
  |author1=F. I. Moxley III |author2=T. Byrnes |author3=F. Fujiwara |author4=W. Dai | title = A generalized finite-difference time-domain quantum method for the N-body interacting Hamiltonian
  | journal = Computer Physics Communications
  | volume = 183
  | issue = 11
  | pages = 2434–2440
  | year = 2012
  | doi=10.1016/j.cpc.2012.06.012
 |bibcode = 2012CoPhC.183.2434M }}</ref>

<ref name="Moxley2014">{{cite book 
  | author=Frederick Moxley 
  | display-authors=etal 
  | title=Contemporary Mathematics: Mathematics of Continuous and Discrete Dynamical Systems 
  | publisher=American Mathematical Society 
  | year=2014 
  | isbn=978-0-8218-9862-8 
  | url=https://www.ams.org/bookstore-getitem?item=CONM-618| author-link=Frederick Moxley 
 }}
</ref>

<ref name="Moxley2013">
{{cite journal
  |author1=F. I. Moxley III |author2=D. T. Chuss |author3=W. Dai | title = A generalized finite-difference time-domain scheme for solving nonlinear Schrödinger equations
  | journal = Computer Physics Communications
  | volume = 184
  | issue = 8
  | pages = 1834–1841
  | year = 2013
  | doi=10.1016/j.cpc.2013.03.006
 |bibcode = 2013CoPhC.184.1834M }}</ref>

<ref name="vonneumann49">
{{cite journal
  |author1=J. von Neumann |author2=RD Richtmyer | title = A method for the numerical calculation of hydrodynamic shocks
  | journal = Journal of Applied Physics
  | volume = 21
  |issue=3 | pages = 232–237
  |date=March 1950
  | doi=10.1063/1.1699639
 |bibcode = 1950JAP....21..232V }}</ref>

<ref name="phot_cad">
{{cite journal
| author= D. Gallagher
| url=http://www.photond.com/files/docs/leos_newsletter_feb08_article.pdf
| title=Photonics CAD Matures
| journal=LEOS Newsletter
| year=2008
}}</ref>

<ref name="back_pml">
{{cite journal
 |author1=A. Deinega |author2=I. Valuev |title=Long-time behavior of PML absorbing boundaries for layered periodic structures
 |journal=Comput. Phys. Commun.
 |volume= 182
 |issue=1 |pages= 149–151
 |year=2011
 |doi=10.1016/j.cpc.2010.06.006
|bibcode = 2011CoPhC.182..149D }}</ref>

<ref name="navarro94">
{{cite journal
 |author1=E.A. Navarro |author2=C. Wu |author3=P.Y. Chung |author4=J. Litva
 |title=Application of PML superabsorbing boundary condition to non-orthogonal FDTD method
 |journal=Electronics Letters
 |volume= 30
 |issue=20 |pages= 1654–1656
 |year=1994
 |doi=10.1049/el:19941139
 |bibcode=1994ElL....30.1654N }}</ref>

<ref name="SorianoNavarro2004">
{{cite journal
 |author1=A. Soriano |author2=E.A. Navarro |author3=J. Portí |author4=V. Such
 |title=Analysis of the finite difference time domain technique to solve the Schrödinger equation for quantum devices
 |journal=Journal of Applied Physics
 |volume= 95
 |issue=12 |pages= 8011–8018
 |year=2004
 |doi= 10.1063/1.1753661
 |bibcode=2004JAP....95.8011S |hdl=10550/12837
 |hdl-access=free
 }}</ref>

<ref name="berenger94">{{cite journal 
  | author= J. Berenger 
  | title= A perfectly matched layer for the absorption of electromagnetic waves 
  | journal= Journal of Computational Physics 
  | year= 1994 
  | volume= 114 
  | pages= 185–200 
  | url=http://www.stanford.edu/class/ee256/Berenger1994.pdf 
  | doi= 10.1006/jcph.1994.1159|bibcode = 1994JCoPh.114..185B 
  | issue= 2 }}
</ref>

<ref name="choi86">{{cite journal 
  |author1=D. H. Choi |author2=W. J. Hoefer | title= The finite-difference time-domain method and its application to eigenvalue problems 
  | journal= IEEE Transactions on Microwave Theory and Techniques
  | year= 1986 
  | volume= 34 
  | pages= 1464–1470  
  | doi= 10.1109/TMTT.1986.1133564|bibcode = 1986ITMTT..34.1464C 
  | issue= 12 }}  
</ref>

<ref name="de_raedt03">{{cite journal 
  |author1=H. De Raedt |author2=K. Michielsen |author3=J. S. Kole |author4=M. T. Figge | title= Solving the Maxwell equations by the Chebyshev method: A one-step finite difference time-domain algorithm 
  | journal= IEEE Transactions on Antennas and Propagation 
  | year= 2003 
  | volume= 51 
  | pages= 3155–3160  
  | doi=10.1109/TAP.2003.818809|arxiv = physics/0208060 |bibcode = 2003ITAP...51.3155D 
  | issue= 11 |s2cid=119095479 }}
</ref>

<ref name="el-ghazaly90">{{cite journal 
  |author1=S. M. El-Ghazaly |author2=R. P. Joshi |author3=R. O. Grondin | title= Electromagnetic and transport considerations in subpicosecond photoconductive switch modeling 
  | journal= IEEE Transactions on Microwave Theory and Techniques 
  | year= 1990 
  | volume= 38 
  | pages= 629–637  
  | doi=10.1109/22.54932|bibcode = 1990ITMTT..38..629E 
  | issue= 5 }} 
</ref>

.<ref name="fitting">
{{cite web
 |url=http://fdtd.kintechlab.com/en/fitting
 |title=Fitting of dielectric function
}}</ref>

<ref name="gedney96">{{cite journal 
  | author= S. D. Gedney 
  | title= An anisotropic perfectly matched layer absorbing media for the truncation of FDTD lattices 
  | journal= IEEE Transactions on Antennas and Propagation
  | year= 1996 
  | volume= 44 
  | pages= 1630–1639  
  | doi=10.1109/8.546249|bibcode = 1996ITAP...44.1630G 
  | issue= 12 }} 
</ref>

<ref name="goorjian92">{{cite journal 
  |author1=P. M. Goorjian |author2=A. Taflove | title= Direct time integration of Maxwell's equations in nonlinear dispersive media for propagation and scattering of femtosecond electromagnetic solitons 
  | journal= Optics Letters 
  | year= 1992 
  | volume= 17 
  | pages= 180–182 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper37.pdf 
  | doi= 10.1364/OL.17.000180|pmid=19784268 |bibcode = 1992OptL...17..180G 
  | issue= 3 }} 
</ref>

<ref name="gwarek85">{{cite journal 
  | author= W. Gwarek 
  | title= Analysis of an arbitrarily shaped planar circuit — A time-domain approach 
  | journal= IEEE Transactions on Microwave Theory and Techniques
  | year= 1985 
  | volume= 33 
  | pages= 1067–1072 
  | doi= 10.1109/TMTT.1985.1133170|bibcode = 1985ITMTT..33.1067G 
  | issue= 10 }}  
</ref>

<ref name="hagness98">{{cite journal 
  |author1=S. C. Hagness |author2=A. Taflove |author3=J. E. Bridges | title= Two-dimensional FDTD analysis of a pulsed microwave confocal system for breast cancer detection: Fixed-focus and antenna-array sensors 
  | journal= IEEE Transactions on Biomedical Engineering
  | year= 1998 
  | volume= 45 
  | pages= 1470–1479 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper66.pdf 
  | doi= 10.1109/10.730440 
  | pmid= 9835195 
  | issue= 12|s2cid=6169784 }} 
</ref>

<ref name="Hao">{{cite book 
  |author1=Y. Hao |author2=R. Mittra | title=FDTD Modeling of Metamaterials: Theory and Applications
  | publisher=Artech House Publishers 
  | year=2009
  | url=http://www.artechhouse.com/Detail.aspx?strIsbn=978-1-59693-160-2}}
</ref>

<ref name="hayakawa02">{{cite journal 
  |author1=M. Hayakawa |author2=T. Otsuyama | title= FDTD analysis of ELF wave propagation in inhomogeneous subionospheric waveguide models 
  | journal= ACES Journal  
  | year= 2002 
  | volume= 17 
  | pages= 239–244 
  | url=http://handle.dtic.mil/100.2/ADP013476| archive-url=https://web.archive.org/web/20120527162538/http://handle.dtic.mil/100.2/ADP013476| url-status=dead| archive-date=May 27, 2012}}
</ref>

<ref name="simpson02">{{cite journal
 |doi         = 10.1109/LAWP.2002.805123
 |author1     = J. J. Simpson
 |author2     = A. Taflove
 |title       = Two-dimensional FDTD model of antipodal ELF propagation and Schumann resonance of the Earth
 |journal     = IEEE Antennas and Wireless Propagation Letters
 |year        = 2002
 |volume      = 1
 |pages       = 53–56
 |issue       = 2
 |url         = http://www.ece.unm.edu/~simpson/Paper1.pdf
 |bibcode     = 2002IAWPL...1...53S
 |url-status     = dead
 |archive-url  = https://web.archive.org/web/20100617131759/http://www.ece.unm.edu/~simpson/Paper1.pdf
 |archive-date = 2010-06-17
|citeseerx     = 10.1.1.694.4837
 |s2cid = 368964
 }}
</ref>

<ref name="holland77">{{cite journal 
  | author= R. Holland 
  | title= Threde: A free-field EMP coupling and scattering code 
  | journal= IEEE Transactions on Nuclear Science
  | year= 1977 
  | volume= 24 
  | pages= 2416–2421  
  | doi= 10.1109/TNS.1977.4329229|bibcode = 1977ITNS...24.2416H 
  | issue= 6 | s2cid= 35395821 
 }} 
</ref>

<ref name="joseph91">{{cite journal 
  |author1=R. M. Joseph |author2=S. C. Hagness |author3=A. Taflove | title= Direct time integration of Maxwell's equations in linear dispersive media with absorption for scattering and propagation of femtosecond electromagnetic pulses 
  | journal= Optics Letters 
  | year= 1991 
  | volume= 16 
  | pages= 1412–4 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper35.pdf 
  | doi= 10.1364/OL.16.001412 
  | pmid= 19776986 
  | issue= 18|bibcode = 1991OptL...16.1412J }} 
</ref>

<ref name="joseph93">{{cite journal 
  |author1=R. M. Joseph |author2=P. M. Goorjian |author3=A. Taflove | title= Direct time integration of Maxwell's equations in 2-D dielectric waveguides for propagation and scattering of femtosecond electromagnetic solitons 
  | journal= Optics Letters 
  | year= 1993 
  | volume= 18 
  | pages= 491–3 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper44.pdf 
  | doi= 10.1364/OL.18.000491 
  | pmid= 19802177 
  | issue= 7|bibcode = 1993OptL...18..491J }} 
</ref>

<ref name="joseph94">{{cite journal 
  |author1=R. M. Joseph |author2=A. Taflove | title= Spatial soliton deflection mechanism indicated by FDTD Maxwell's equations modeling 
  | journal= IEEE Photonics Technology Letters
  | year= 1994 
  | volume= 2 
  |issue=10 | pages= 1251–1254 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper54.pdf
  | doi=10.1109/68.329654|bibcode = 1994IPTL....6.1251J |s2cid=46710331 }} 
</ref>

<ref name="jurgens92">{{cite journal 
  |author1=T. G. Jurgens |author2=A. Taflove |author3=K. R. Umashankar |author4=T. G. Moore | title= Finite-difference time-domain modeling of curved surfaces 
  | journal= IEEE Transactions on Antennas and Propagation 
  | year= 1992 
  | volume= 40 
  | pages= 357–366 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper39.pdf 
  | doi= 10.1109/8.138836|bibcode = 1992ITAP...40..357J 
  | issue= 4 }} 
</ref>

<ref name="kashiwa90">{{cite journal 
  |author1=T. Kashiwa |author2=I. Fukai | title= A treatment by FDTD method of dispersive characteristics associated with electronic polarization 
  | journal= [[Microwave and Optical Technology Letters]] 
  | year= 1990 
  | volume= 3 
  | pages= 203–205 
  | doi= 10.1002/mop.4650030606 
  | issue= 6}} 
</ref>

<ref name="katz91">{{cite journal 
  |author1=D. S. Katz |author2=A. Taflove|author2-link=Allen Taflove |author3=M. J. Piket-May|authorlink3=Melinda Piket-May |author4=K. R. Umashankar | title= FDTD analysis of electromagnetic wave radiation from systems containing horn antennas 
  | journal= IEEE Transactions on Antennas and Propagation 
  | year= 1991 
  | volume= 39 
  | pages= 1203–1212 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper34.pdf 
  | doi= 10.1109/8.97356|bibcode = 1991ITAP...39.1203K 
  | issue= 8 }} 
</ref>

<ref name="katz94">{{cite journal 
  |author1=D. S. Katz |author2=E. T. Thiele |author3=A. Taflove | title= Validation and extension to three dimensions of the Berenger PML absorbing boundary condition for FDTD meshes 
  | journal= IEEE Microwave and Guided Wave Letters
  | year= 1994 
  | volume= 4 
  | pages= 268–270 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper51.pdf 
  | doi= 10.1109/75.311494 
  | issue= 8|s2cid=10156811 }} 
</ref>

<ref name="kriegsmann87">{{cite journal 
  |author1=G. A. Kriegsmann |author2=A. Taflove |author3=K. R. Umashankar | title= A new formulation of electromagnetic wave scattering using an on-surface radiation boundary condition approach 
  | journal= IEEE Transactions on Antennas and Propagation
  | year= 1987 
  | volume= 35 
  | pages= 153–161 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper14.pdf 
  | doi= 10.1109/TAP.1987.1144062|bibcode = 1987ITAP...35..153K 
  | issue= 2 }} 
</ref>

<ref name="kunz77">{{cite journal 
  |author1=K. S. Kunz |author2=K. M. Lee | title= A three-dimensional finite-difference solution of the external response of an aircraft to a complex transient EM environment 
  | journal= [[IEEE Transactions on Electromagnetic Compatibility|IEEE Trans. Electromagn. Compat.]]
  | year= 1978 
  | volume= 20 
  | pages= 333–341 
  | doi=10.1109/TEMC.1978.303727 
  | issue= 2|s2cid=31666283 }}  
</ref>

<ref name="liao84">{{cite journal 
  |author1=Z. P. Liao |author2=H. L. Wong |author3=B. P. Yang |author4=Y. F. Yuan | title= A transmitting boundary for transient wave analysis
  | journal= Scientia Sinica, Series A 
  | year= 1984 
  | volume= 27 
  | pages= 1063–1076}}
</ref>

<ref name="liu97">{{cite book 
  | author= Q. H. Liu 
  | title= IEEE Antennas and Propagation Society International Symposium 1997. Digest 
 | chapter= The pseudospectral time-domain (PSTD) method: A new algorithm for solutions of Maxwell's equations 
  | year= 1997 
  | volume= 1 
  | pages= 122–125  
  | doi=10.1109/APS.1997.630102 
  | isbn= 978-0-7803-4178-4| s2cid= 21345353 
 }} 
</ref>

<ref name="luebbers90">{{cite journal 
  |author1=R. Luebbers |author2=F. Hunsberger |author3=K. Kunz |author4=R. Standler |author5=M. Schneider | title= A frequency-dependent finite-difference time-domain formulation for dispersive materials 
  | journal= [[IEEE Transactions on Electromagnetic Compatibility|IEEE Trans. Electromagn. Compat.]]
  | year= 1990 
  | volume= 32 
  | pages= 222–227  
  | doi=10.1109/15.57116 
  | issue= 3}} 
</ref>

<ref name="maloney90">{{cite journal 
  |author1=J. G. Maloney |author2=G. S. Smith |author3=W. R. Scott Jr. | title= Accurate computation of the radiation from simple antennas using the finite-difference time-domain method 
  | journal= IEEE Transactions on Antennas and Propagation 
  | year= 1990 
  | volume= 38 
  | pages= 1059–1068  
  | doi=10.1109/8.55618|bibcode = 1990ITAP...38.1059M 
  | issue= 7 |s2cid=31583883 }} 
</ref>

<ref name="maloney98">{{cite journal 
  |author1=J. G. Maloney |author2=M. P. Kesler | title= Analysis of Periodic Structures 
  | journal= Chap. 6 in Advances in Computational Electrodynamics: The Finite-Difference Time-Domain Method, A. Taflove, Ed., Artech House, Publishers 
  | year= 1998}} 
</ref>

<ref name="moore88">{{cite journal 
  |author1=T. G. Moore |author2=J. G. Blaschak |author3=A. Taflove |author4=G. A. Kriegsmann | title= Theory and application of radiation boundary operators 
  | journal= IEEE Transactions on Antennas and Propagation
  | year= 1988 
  | volume= 36 
  | pages= 1797–1812 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper20.pdf 
  | doi= 10.1109/8.14402|bibcode = 1988ITAP...36.1797M 
  | issue= 12 }} 
</ref>

<ref name="mur81">{{cite journal 
  | author= G. Mur 
  | title= Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic field equations 
  | journal= [[IEEE Transactions on Electromagnetic Compatibility|IEEE Trans. Electromagn. Compat.]]
  | year= 1981 
  | volume= 23 
  | pages= 377–382 
  | doi= 10.1109/TEMC.1981.303970 
  | issue= 4| s2cid= 25768246 
 }}
</ref>

<ref name="nagra98">{{cite journal 
  |author1=A. S. Nagra |author2=R. A. York | title= FDTD analysis of wave propagation in nonlinear absorbing and gain media 
  | journal= IEEE Transactions on Antennas and Propagation
  | year= 1998 
  | volume= 46 
  | pages= 334–340  
  | doi=10.1109/8.662652|bibcode = 1998ITAP...46..334N 
  | issue= 3 }} 
</ref>

<ref name="obl_it">
{{cite journal
 |author1=I. Valuev |author2=A. Deinega |author3=S. Belousov |title=Iterative technique for analysis of periodic structures at oblique incidence in the finite-difference time-domain method
 |journal=Opt. Lett.
 |volume= 33
 |issue=13 |pages= 1491–3
 |year=2008
 |doi=10.1364/ol.33.001491
|pmid=18594675 |bibcode = 2008OptL...33.1491V }}</ref>

<ref name="obl_sfdtd">
{{cite journal
 |author1=A. Aminian |author2=Y. Rahmat-Samii |title=Spectral FDTD: a novel technique for the analysis of oblique incident plane wave on periodic structures
 |journal=IEEE Transactions on Antennas and Propagation
 |volume= 54
 |issue=6 |pages= 1818–1825
 |year=2006
 |doi=10.1109/tap.2006.875484
|bibcode = 2006ITAP...54.1818A |s2cid=25120679 }}</ref>

<ref name="ramahi97">{{cite journal 
  | author= O. M. Ramahi 
  | title= The complementary operators method in FDTD simulations 
  | journal= IEEE Antennas and Propagation Magazine
  | year= 1997 
  | volume= 39 
  | pages= 33–45  
  | doi=10.1109/74.646801|bibcode = 1997IAPM...39...33R 
  | issue= 6 }} 
</ref>

<ref name="reuter94">{{cite journal 
  |author1=C. E. Reuter |author2=R. M. Joseph |author3=E. T. Thiele |author4=D. S. Katz |author5=A. Taflove | title= Ultrawideband absorbing boundary condition for termination of waveguiding structures in FDTD simulations 
  | journal= IEEE Microwave and Guided Wave Letters
  | year= 1994 
  | volume= 4 
  | pages= 344–346 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper53.pdf 
  | doi= 10.1109/75.324711 
  | issue= 10|s2cid=24572883 }} 
</ref>

<ref name="chewweedon94">{{cite journal 
  |author1=W.C. Chew |author2=W.H. Weedon | title= A 3D perfectly matched medium from modified Maxwell's equations with stretched coordinates
  | journal= Microwave and Optical Technology Letters 
  | year= 1994 
  | volume= 7 
  | pages= 599–604
  | doi= 10.1002/mop.4650071304
  | issue= 13|bibcode=1994MiOTL...7..599C }} 
</ref>

<ref name="roden00">{{cite journal 
  |author1=J. A. Roden |author2=S. D. Gedney | title= Convolution PML (CPML): An efficient FDTD implementation of the CFS-PML for arbitrary media 
  | journal= [[Microwave and Optical Technology Letters]]  
  | year= 2000 
  | volume= 27 
  | pages= 334–339 
  | url= http://www3.interscience.wiley.com/journal/73504513/abstract
  | archive-url= https://archive.today/20130105125334/http://www3.interscience.wiley.com/journal/73504513/abstract
  | url-status= dead
  | archive-date= 2013-01-05
  | doi=10.1002/1098-2760(20001205)27:5<334::AID-MOP14>3.0.CO;2-A 
  | issue= 5}}
</ref>

<ref name="rylander00">{{cite journal 
  |author1=T. Rylander |author2=A. Bondeson | title= Stable FDTD-FEM hybrid method for Maxwell's equations 
  | journal= Computer Physics Communications  
  | year= 2000 
  | volume= 125 
  |issue=1–3 | pages= 75–82 
  | doi=10.1016/S0010-4655(99)00463-4}}
</ref>

<ref name="sacks95">{{cite journal 
  |author1=Z. S. Sacks |author2=D. M. Kingsland |author3=R. Lee |author4=J. F. Lee | title= A perfectly matched anisotropic absorber for use as an absorbing boundary condition 
  | journal= IEEE Transactions on Antennas and Propagation
  | year= 1995 
  | volume= 43 
  | pages= 1460–1463 
  | doi=10.1109/8.477075|bibcode = 1995ITAP...43.1460S 
  | issue= 12 }} 
</ref>

<ref name="sano90">{{cite journal 
  |author1=E. Sano |author2=T. Shibata | title= Fullwave analysis of picosecond photoconductive switches 
  | journal= IEEE Journal of Quantum Electronics 
  | year= 1990 
  | volume= 26 
  | pages= 372–377  
  | doi=10.1109/3.44970|bibcode = 1990IJQE...26..372S 
  | issue= 2 }} 
</ref>

<ref name="schneider99">{{cite journal 
  |author1=J. B. Schneider |author2=C. L. Wagner | title= FDTD dispersion revisited: Faster-than-light propagation 
  | journal= IEEE Microwave and Guided Wave Letters
  | year= 1999 
  | volume= 9 
  | pages= 54–56  
  | doi=10.1109/75.755044 
  | issue= 2|citeseerx=10.1.1.77.9132 }} 
</ref>

<ref name="sui92">{{cite journal 
  |author1=W. Sui |author2=D. A. Christensen |author3=C. H. Durney | title= Extending the two-dimensional FDTD method to hybrid electromagnetic systems with active and passive lumped elements 
  | journal= IEEE Transactions on Microwave Theory and Techniques 
  | year= 1992 
  | volume= 40 
  | pages= 724–730  
  | doi=10.1109/22.127522|bibcode = 1992ITMTT..40..724S 
  | issue= 4 }} 
</ref>

<ref name="sullivan88">{{cite journal 
  |author1=D. M. Sullivan |author2=O. P. Gandhi |author3=A. Taflove | title= Use of the finite-difference time-domain method in calculating EM absorption in man models 
  | journal= IEEE Transactions on Biomedical Engineering
  | year= 1988 
  | volume= 35 
  | pages= 179–186 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper18.pdf 
  | doi= 10.1109/10.1360 
  | pmid= 3350546 
  | issue= 3|s2cid=20350396 }} 
</ref>

<ref name="taflove75a">{{cite journal 
  |author1=A. Taflove |author2=M. E. Brodwin | title= Numerical solution of steady-state electromagnetic scattering problems using the time-dependent Maxwell's equations 
  | journal= IEEE Transactions on Microwave Theory and Techniques
  | year= 1975 
  | volume= 23 
  | pages= 623–630
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper2.pdf 
  | doi= 10.1109/TMTT.1975.1128640 |bibcode = 1975ITMTT..23..623T 
  | issue= 8 }}
</ref>

<ref name="taflove75b">{{cite journal 
  |author1=A. Taflove |author2=M. E. Brodwin | title= Computation of the electromagnetic fields and induced temperatures within a model of the microwave-irradiated human eye 
  | journal= IEEE Transactions on Microwave Theory and Techniques
  | year= 1975 
  | volume= 23 
  | pages= 888–896
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper3.pdf 
  | doi= 10.1109/TMTT.1975.1128708 |bibcode = 1975ITMTT..23..888T 
  | issue= 11 }} 
</ref>

<ref name="taflove80">{{cite journal 
  | author= A. Taflove 
  | title= Application of the finite-difference time-domain method to sinusoidal steady state electromagnetic penetration problems 
  | journal= [[IEEE Transactions on Electromagnetic Compatibility|IEEE Trans. Electromagn. Compat.]]
  | year= 1980 
  | volume= 22 
  | pages= 191–202
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper7.pdf 
  | doi= 10.1109/TEMC.1980.303879 
  | issue= 3 | bibcode= 1980ITElC..22..191T 
 | s2cid= 39236486 
 }}
</ref>

<ref name="taflove83">{{cite journal 
  |author1=A. Taflove |author2=K. R. Umashankar | title= Radar cross section of general three-dimensional scatterers 
  | journal= [[IEEE Transactions on Electromagnetic Compatibility|IEEE Trans. Electromagn. Compat.]]
  | year= 1983 
  | volume= 25 
  | pages= 433–440 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper10.pdf 
  | doi=10.1109/TEMC.1983.304133 
  | issue= 4 |s2cid=40419955 }} 
</ref>

<ref name="taflove88">{{cite journal 
  |author1=A. Taflove |author2=K. R. Umashankar |author3=B. Beker |author4=F. A. Harfoush |author5=K. S. Yee | title= Detailed FDTD analysis of electromagnetic fields penetrating narrow slots and lapped joints in thick conducting screens 
  | journal= IEEE Transactions on Antennas and Propagation 
  | year= 1988 
  | volume= 36 
  | pages= 247–257 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper17.pdf 
  | doi= 10.1109/8.1102|bibcode = 1988ITAP...36..247T 
  | issue= 2 }} 
</ref>

<ref name="taflove05">{{cite book 
  | author=[[Allen Taflove]] and [[Susan Hagness|Susan C. Hagness]] 
  | title=Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed. 
  | publisher=Artech House Publishers 
  | year=2005 
  | isbn=978-1-58053-832-9 
  | url=http://www.artechhouse.com/Detail.aspx?strBookId=1123}}
</ref>

<ref name="thomas94">{{cite journal 
  |author1=V. A. Thomas |author2=M. E. Jones |author3=M. J. Piket-May |authorlink3=Melinda Piket-May|author4=A. Taflove |author5=E. Harrigan | title= The use of SPICE lumped circuits as sub-grid models for FDTD high-speed electronic circuit design 
  | journal= IEEE Microwave and Guided Wave Letters
  | year= 1994 
  | volume= 4 
  |issue=5 | pages= 141–143 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper49.pdf
  | doi=10.1109/75.289516|s2cid=32905331 }} 
</ref>

<ref name="tirkas91">{{cite book 
  |author1=P. A. Tirkas |author2=C. A. Balanis |title=Antennas and Propagation Society Symposium 1991 Digest |chapter=Finite-difference time-domain technique for radiation by horn antennas | year= 1991 
  | volume= 3 
  | pages= 1750–1753  
  | doi=10.1109/APS.1991.175196 
  | isbn= 978-0-7803-0144-3|s2cid=122038624 }} 
</ref>

<ref name="TMatrix">
{{cite journal
 |author1=A. Deinega |author2=S. Belousov |author3=I. Valuev |s2cid=27742034 |title=Hybrid transfer-matrix FDTD method for layered periodic structures
 |journal=Opt. Lett. 
 |volume= 34
 |issue=6 |pages= 860–2
 |year=2009
 |doi=10.1364/ol.34.000860
|pmid=19282957 |bibcode = 2009OptL...34..860D }}</ref>

<ref name="toland93">{{cite journal 
  |author1=B. Toland |author2=B. Houshmand |author3=T. Itoh | title= Modeling of nonlinear active regions with the FDTD method 
  | journal= IEEE Microwave and Guided Wave Letters 
  | year= 1993 
  | volume= 3 
  | pages= 333–335  
  | doi=10.1109/75.244870 
  | issue= 9|s2cid=27549555 }} 
</ref>

<ref name="umashankar82">{{cite journal 
  |author1=K. R. Umashankar |author2=A. Taflove | title= A novel method to analyze electromagnetic scattering of complex objects 
  | journal= [[IEEE Transactions on Electromagnetic Compatibility|IEEE Trans. Electromagn. Compat.]]
  | year= 1982 
  | volume= 24 
  | pages= 397–405
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper9.pdf 
  | doi= 10.1109/TEMC.1982.304054 
  | issue= 4 |bibcode=1982ITElC..24..397U |s2cid=37962500 }}
</ref>

<ref name="umashankar87">{{cite journal 
  |author1=K. R. Umashankar |author2=A. Taflove |author3=B. Beker | title= Calculation and experimental validation of induced currents on coupled wires in an arbitrary shaped cavity 
  | journal= IEEE Transactions on Antennas and Propagation 
  | year= 1987 
  | volume= 35 
  | pages= 1248–1257 
  | url=http://www.ece.northwestern.edu/ecefaculty/taflove/Paper16.pdf 
  | doi= 10.1109/TAP.1987.1144000|bibcode = 1987ITAP...35.1248U 
  | issue= 11 }} 
</ref>

<ref name="yee66">{{cite journal 
  | author= Kane Yee 
  | title= Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media 
  | journal= IEEE Transactions on Antennas and Propagation
  | year= 1966 
  | volume= 14 
  | pages= 302–307 
  | doi= 10.1109/TAP.1966.1138693 
  |bibcode = 1966ITAP...14..302Y 
  | issue= 3 | s2cid= 122712881 
 }}
</ref>

<ref name="lam69">{{cite journal 
  | author= Dong-Hoa Lam
  | title= Finite Difference Methods for Electromagnetic Scattering Problems
  | journal= Mississippi State University, Interaction Notes
  | volume= 44 
  | year= 1969 
  | url= http://ece-research.unm.edu/summa/notes/In/0044.pdf }}
</ref>

<ref name="zhang88">{{cite journal 
  |author1=X. Zhang |author2=J. Fang |author3=K. K. Mei |author4=Y. Liu | title= Calculation of the dispersive characteristics of microstrips by the time-domain finite-difference method 
  | journal= IEEE Transactions on Microwave Theory and Techniques
  | year= 1988 
  | volume= 36 
  | pages= 263–267  
  | doi=10.1109/22.3514|bibcode = 1988ITMTT..36..263Z 
  | issue= 2 }}  
</ref>

<ref name="zhen00">{{cite journal 
  |author1=F. Zhen |author2=Z. Chen |author3=J. Zhang | title= Toward the development of a three-dimensional unconditionally stable finite-difference time-domain method 
  | journal= IEEE Transactions on Microwave Theory and Techniques  
  | year= 2000 
  | volume= 48 
  | pages= 1550–1558  
  | doi=10.1109/22.869007|bibcode = 2000ITMTT..48.1550Z 
  | issue= 9 }}
</ref>

<ref name="zheng01">{{cite journal 
  |author1=F. Zheng |author2=Z. Chen | title= Numerical dispersion analysis of the unconditionally stable 3-D ADI-FDTD method 
  | journal= IEEE Transactions on Microwave Theory and Techniques  
  | year= 2001 
  | volume= 49 
  | pages= 1006–1009  
  | doi=10.1109/22.920165|bibcode = 2001ITMTT..49.1006Z 
  | issue= 5 }}
</ref>

<ref name="ziolkowski93">{{cite journal 
  |author1=R. W. Ziolkowski |author2=J. B. Judkins | title= Full-wave vector Maxwell's equations modeling of self-focusing of ultra-short optical pulses in a nonlinear Kerr medium exhibiting a finite response time 
  | journal= Journal of the Optical Society of America B
  | year= 1993 
  | volume= 10 
  | pages= 186–198 
  | doi= 10.1364/JOSAB.10.000186|bibcode = 1993JOSAB..10..186Z 
  | issue= 2 }} 
</ref>

<ref name="Ahmed2008">{{cite journal 
  |author1=I. Ahmed |author2=E. K. Chua |author3=E. P. Li |author4=Z. Chen | title= Development of the three-dimensional unconditionally stable LOD-FDTD method
  | journal= IEEE Transactions on Antennas and Propagation
  | year= 2008 
  | volume= 56 
  | pages= 3596–3600 
  | doi= 10.1109/TAP.2008.2005544|bibcode = 2008ITAP...56.3596A 
  | issue= 11 |s2cid=31351974 }} 
</ref>

<ref name="aoyagi93">{{cite journal 
  | title=A hybrid Yee algorithm/scalar-wave equation approach
  | journal=IEEE Transactions on Microwave Theory and Techniques
  |author=Aoyagi, P.H. |author2=Lee, J.F. |author3=Mittra, R.
  | volume=41
  | number=9
  | pages=1593–1600
  | year=1993
  |bibcode = 1993ITMTT..41.1593A |doi = 10.1109/22.245683 }} 
</ref>

<ref name="oliveira09">{{cite journal 
  | title=Computational Environment for Simulating Lightning Strokes in a Power Substation by Finite-Difference Time-Domain Method
  |author=R. M. S. de Oliveira |author2=C. L. S. S. Sobrinho
  | journal=IEEE Transactions on Electromagnetic Compatibility
  | volume=51
  | issue=4
  | pages=995–1000
  | year=2009
  | doi = 10.1109/TEMC.2009.2028879 }} 
</ref>

<ref name="baba08">{{cite journal 
  | title=An Improved Thin Wire Representation for FDTD Computations
  |author= Taniguchi, Y. |author2=Baba, Y.  |author3=N. Nagaoka |author4=A. Ametani
  | journal=IEEE Transactions on Antennas and Propagation
  | volume=56
  | number=10
  | pages=3248–3252
  | year=2008
  | doi = 10.1109/TAP.2008.929447 |bibcode = 2008ITAP...56.3248T |s2cid= 29617214 }}
</ref>

<ref name="oliveira2021">{{cite journal 
  | title=Least Squares Finite-Difference Time-Domain
  | author=R. M. S. de Oliveira | author2=R. R. Paiva
  | journal=IEEE Transactions on Antennas and Propagation
  | year=2021
  | volume=69 | issue=9 | pages=6111–6115 | doi = 10.1109/TAP.2021.3069576 | bibcode=2021ITAP...69.6111D | s2cid=234307029 }} 
</ref>

}}

== Further reading ==
{{Refbegin}}
The following article in ''Nature Milestones: Photons'' illustrates the historical significance of the FDTD method as related to Maxwell's equations:
*{{cite journal
  | url=http://www.nature.com/milestones/milephotons/full/milephotons02.html
  | journal=Nature Milestones: Photons 
  | title=Milestone 2 (1861) Maxwell's equations
  | author=David Pile
  |date=May 2010
  | doi=10.1038/nmat2639
  | access-date=17 June 2010
| doi-access=free
  }}

Allen Taflove's interview, "Numerical Solution," in the January 2015 focus issue of ''Nature Photonics'' honoring the 150th anniversary of the publication of Maxwell's equations.  This interview touches on how the development of FDTD ties into the century and one-half history of Maxwell's theory of electrodynamics:
*[http://www.nature.com/nphoton/focus/maxwell-anniversary/index.html ''Nature Photonics interview'']

The following university-level textbooks provide a good general introduction to the FDTD method:

*{{cite book 
  |author=Karl S. Kunz |author2=Raymond J. Luebbers 
  | title=The Finite Difference Time Domain Method for Electromagnetics 
  | publisher=CRC Press 
  | year=1993 
  | isbn=978-0-8493-8657-2 
  | url=http://www.crcpress.com/shopping_cart/products/product_detail.asp?sku=8657&af=W1129 
  | access-date=2006-08-05 
  | archive-url=https://web.archive.org/web/20071210045441/http://www.crcpress.com/shopping_cart/products/product_detail.asp?sku=8657&af=W1129 
  | archive-date=2007-12-10 
  | url-status=dead 
  |author2-link=Raymond Luebbers 
 }}
 
*{{cite book
  |author=Allen Taflove |author2=Susan C. Hagness 
  | title=Computational Electrodynamics: The Finite-Difference Time-Domain Method, 3rd ed. 
  | publisher=Artech House Publishers 
  | year=2005 
  | isbn=978-1-58053-832-9 
  | url=http://www.artechhouse.com/Detail.aspx?strBookId=1123|author-link=Allen Taflove 
 }}

*{{cite book
  |author1=Wenhua Yu |author2=Raj Mittra |author3=Tao Su |author4=Yongjun Liu |author5=Xiaoling Yang | title=Parallel Finite-Difference Time-Domain Method 
  | publisher=Artech House Publishers 
  | year=2006 
  | isbn=978-1-59693-085-8
  | url=http://www.artechhouse.com/default.asp?frame=book.asp&book=1-59693-085-3&Country=US&Continent=NO&State=}}

*{{cite book
  | author=John B. Schneider 
  | title=Understanding the FDTD Method
  | publisher=available online
  | year=2010
  | url=http://eecs.wsu.edu/~schneidj/ufdtd/index.php}}
{{Refend}}
*

== External links ==
{{commons category}}
[[Free software]]/[[Open-source software]] FDTD projects:
* [http://www.fdtdxx.com FDTD++]: advanced, fully featured FDTD software, along with sophisticated material models and predefined fits as well as discussion/support forums and email support
* [http://openEMS.de openEMS] (Fully 3D Cartesian & Cylindrical graded mesh EC-FDTD Solver, written in C++, using a [[Matlab]]/[[GNU Octave|Octave]]-Interface)
* [https://web.archive.org/web/20110517102321/http://www.its.caltech.edu/~seheon/FDTD.html pFDTD] (3D C++ FDTD codes developed by Se-Heon Kim)
* [https://web.archive.org/web/20090626051810/http://www.thecomputationalphysicist.com/ JFDTD] (2D/3D C++ FDTD codes developed for nanophotonics by Jeffrey M. McMahon)
* [http://www.ece.ncsu.edu/oleg/wiki/WOLFSIM WOLFSIM] {{Webarchive|url=https://web.archive.org/web/20080702190617/http://www.ece.ncsu.edu/oleg/wiki/WOLFSIM |date=2008-07-02 }} (NCSU) (2-D)
* [http://ab-initio.mit.edu/meep/ Meep] ([[Massachusetts Institute of Technology|MIT]], 2D/3D/cylindrical parallel FDTD)
* [http://freshmeat.net/projects/radarfdtd/ (Geo-) Radar FDTD]
* [http://sourceforge.net/projects/bigboy bigboy] (unmaintained, no release files. must get source from cvs)
* [http://sourceforge.net/projects/pfdtd/files/ Parallel (MPI&OpenMP) FDTD codes in C++] (developed by Zs. Szabó)
* [https://archive.today/20121217222254/http://cs.tu-berlin.de/~peutetre/sfdtd/ FDTD code in Fortran 90]
* [http://code.google.com/p/emwave2d/ FDTD code in C for 2D EM Wave simulation]
* [http://angorafdtd.org Angora] (3D parallel FDTD software package, maintained by Ilker R. Capoglu)
* [http://gsvit.net/ GSvit] (3D FDTD solver with graphics card computing support, written in C, graphical user interface XSvit available)
*[http://www.gprmax.com gprMax] (Open Source (GPLv3), 3D/2D FDTD modelling code in Python/Cython developed for GPR but can be used for general EM modelling.)
[[Freeware]]/[[Closed source]] FDTD projects (some not for commercial use):
* [http://fdtd.kintechlab.com/en/start EMTL (Electromagnetic Template Library)] (Free С++ library for electromagnetic simulations. The current version implements mainly the FDTD).

{{Numerical PDE}}

{{DEFAULTSORT:Finite-Difference Time-Domain Method}}
[[Category:Numerical software]]
[[Category:Simulation software]]
[[Category:Electromagnetic radiation]]
[[Category:Numerical differential equations]]
[[Category:Computational science]]
[[Category:Computational electromagnetics]]
[[Category:Electromagnetism]]
[[Category:Electrodynamics]]
[[Category:Scattering, absorption and radiative transfer (optics)]]