Non-Stationary Electromagnetics, 2nd Edition

The calculation algorithm is similar to the integration algorithm for the quadratic mesh described in the previous section. Because the odd time step points are spatially shifted onto the half-cell from the points on the even time steps, the calculation routines for even and odd time iterations differ slightly.

(vii) Validation of the algorithm

Let us consider the case when the layer permittivity changes abruptly at τ = 0 from ε to ε₁ and the initial field is a plane wave E₀ = cos(2πη (τ − ξ)). The normalized frequency η = ωα/(3πv) = α/λ shows the number of full wavelengths λ per layer width a where λ is the wavelength in the background medium. The exact solution for this case is constructed in the previous Section 3.4 and for the central spatial ...