1
EE256 Numerical Electromagnetics
H. O. #6
Marshall
24 June 2008
Summer 2008
LECTURE 3
3.1
FINITE DIFFERENCE SOLUTIONS OF THE CONVECTION EQUATION
As the simplest example of a PDE, we consider ﬁnite difference solutions of the convection equation [2.20].
As mentioned before, this equation constitutes the simplest prototype for hyperbolic equations, and is also
highly relevant for our FDTD analyses of Maxwell’s equations, since it has the same form as the component
curl equations that we shall solve in FDTD solutions of electromagnetic problems. Consider the convection
equation written in terms of voltage on a transmission line:
∂
V
∂t
+
v
p
∂
V
∂z
= 0
[3.1]
where
v
p
’
1
ft(ns)

1
, i.e., the speed of light in free space. As mentioned in Lecture #2, the exact
solution of this convection equation is:
V
(
z,t
) =
V
+
(
z

v
p
t
)
[3.2]
where
V
+
(
·
)
is an arbitrary function. A particular solution is determined by means of the initial condition
at time
t
= 0
, i.e.,
V
+
(
z,
0) = Φ(
z
)
, where
Φ(
z
)
is the initial voltage distribution on the line. The solution
at any later time is then given by:
V
(
z,t
) = Φ(
z

v
p
t
)
[3.3]
In other words, the initial voltage distribution simply is convected (or propagated) to the right at the
velocity
v
p
, preserving its magnitude and shape.
We shall now consider FDTD solutions of [3.1] for a particular initial voltage distribution:
V
(
z,
0) = Φ(
z
) =
12
z
V
0
≤
z
≤
0
.
5ft
,

12(1

z
)
V
0
.
5ft
≤
z
≤
1
.
0ft
,
0
V
Elsewhere
[3.4]
The exact solution is simply the triangular pulse moving to the right at the speed of
v
p
’
1
ft(ns)

1
, as
shown in Figure 3.1.
We now proceed to numerically solve for
V
(
z,t
)
, using two different methods, namely (i) the forward
time centeredspace method, and (ii) the leapfrog method.