2. 1D Wave Equation
The one-dimensional (1D) wave equation governs basic wave propagation.
The wave equation is easily discretized by using the central finite difference model.
We can numerically solve the discrete wave equation and understand wave phenomena.
Standard FDTD Algorithm
The one-dimensional wave equation is defined by

. . . (1)
where ∂t=∂/∂t, (∂x is analogously defined to ∂t), and v is the wave speed.
Discretizing x = iΔx, t = nΔt (i, n = integer), we simply write

. . . (2)
Using the central finite difference model, we find

. . . (3)
where

. . . (4)
Expanding the temporal difference, dt2ψin, we obtain a
"1D standard finite difference time domain (FDTD) algorithm,"

. . . (5)
where

. . . (6)
Nonstandard FDTD Algorithm
The central finite difference model is given by

. . . (7)
whereas we introduce an exact expression,

. . . (8)
where s(Δx) is a correction function. Here we suppose a monochromatic wave,

. . . (9)
where k = wavenumber and ω = angular frequency.
Note that any electromagnetic wave is a superposition of monochromatic waves, we find

. . . (10)
Putting s(Δx) into s(k, Δx), we obtain

. . . (11)
Although this approach is not always valid because the following definition must be satisfied,

. . . (12)
(12) is fortunately satisfied. Thus, (8) is an exact central finite difference model what called
nonstandard finite difference model" [1, 2].
Similarly for the temporal derivative, we find

. . . (13)
Expanding the temporal difference, dt2ψin, we obtain a
"1D nonstandard FDTD algorithm,"

. . . (14)
where

. . . (15)
where we used v = ω / k.
Bibliography:
[1] R. E. Mickens, "Nonstandard Finite Difference Models of Differential Equation," World Scientific (1994).
[2] J. B. Cole, "A high accuracy fdtd algorithm to solve microwave propagation and scattering problems on a coarse grid," IEEE Trans. Microwave Theory and Tech, 43, 9, pp. 2053-2058 (1995).
Copyright (C) 2011 Naoki Okada, All Rights Reserved.