ACFD Lecture 3
ACFD Lecture 3
Control-volume approach:
Eq. (3.96) is somewhat different from Eq. (3.95),
which followed from the most obvious application
of the Taylor-series method to approximate the
mathematical statement of the BC.
Methodology
Taylor-series method: provided difference
approximations to derivatives and the
representation for the complete PDE was mad up
from the addition of serval such representations.
Finite-volume method: employs the conservation
statement corresponding to the entire PDE. The
distinctive characteristic of the FV approach is
that a “ balance” of some physical quantity is
made on the CV in the neighborhood of a grid
point.
It would appear that the discretization developed
by the FV approach would almost certainly have
the conservative property.
Notes
It is difficult to appreciate the subtle differences
that may occur in the difference representations
obtained by the same PDE by using the four
different methods discussed without working a
large number of examples.
In many cases, and especially for simple, linear
equations, the resulting difference equations can
be identical.
There is no guarantee that difference equations
developed by any of the methods will be
numerically stable, so that the same difference
scheme developed by all four methods could turn
out to be worthless.
Stability Considerations
A finite-difference approximation to a PDE may
be consistent, but the solution will not
necessarily converge to the solution of the PDE.
The Lax equivalence theorem states that a
stable numerical method must also be used.
Consider a marching problem in which initial
values at time level n are known and values of
the unknown at time level n+1 are required.
1-D Heat Equation: simple explicit scheme
errors.
The exact solution D and the error H must both satisfy the
same difference equation.
This means that the numerical error and the exact
numerical solution both possess the same growth property
in time and either could be used to examine stability.
We assume the error H(x,t) can be written as a
¦
series of the form: H (x, t) b m (t)e
m ik x
(3.104)
m
where the period of the fundamental frequency
(m=1) is assumed to be 2L.
The wave number km=2Sm/(2L)=S/'x
m=0,1,2,…M(=L/'x)
If M=2 (interval of length 2L is subdivided using
five points) the corresponding frequencies are
fm=km/(2S)=m/(2L) m=0,1,2
f0=0 m=0, f1=1/2L m=1, f2=1/L m=2
Since the difference equation is linear, superposition
equal to 1, a general component of the error will
not grow form one time step to the next.
_c't/'x_d1.
stability requirement is the usual CFL condition:
Notes
It should be noted that the stability analysis
presented above does not include the effect of BCs
even though a matrix notation for the system is used.
The influence of BCs is easily included for systems of
difference equations.
Equation (3.116) shows that the stability of a finite-
difference operator is related to the amplification
matrix. We may also write Eq. (3.116) as
requires that for some positive W, the matrices
The stability condition (Richtmyer & Morton, 1967)
BACK