Unformatted text preview: implied that stability restrictions
were necessary for explicit nitedi erence approximations of the heat equation. The
answer to the dilemma in this case is linked to a loss of consistency unless t ! 0 at
a faster rate than x (again, cf. Problem 1 at the end of this section). Should the
computational mesh be re ned with t= x = (a constant), then the du FortFrankel
scheme (4.3.8) is consistent with the hyperbolic equation ut ; uxx + 2u tt = 0: Thus, the time step of the du FortFrankel scheme must still be restricted to be O( x2)
however, now for reasons of consistency rather than for stability.
Multistep schemes can be started with a onestep method. We must be careful to
preserve accuracy of the multistep scheme when the onestep method has a lower order.
Suppose, for example, that the time step of a secondorder multistep method is tM ,
then the time step of the rstorder starting method should be O( t2 ) in order to keep
M
the local errors in balance. This, however, is not the case when the forward timecentered
space scheme (4.1.2) is used as a starting method for the du FortFrankel (4.3.8) scheme.
Although the order of (4.1.2) is O( t) + O( x2 ) and that of (4.3.8) is O( t2) + O( x2 ),
both schemes restrict t to O( x2 ). Thus, both schemes have O( x2) accuracy. 4.3.1 Matrix Stability Analysis
When using matrix methods to analyze the stability of a multilevel scheme it is convenient
to write the scheme as an equivalent onelevel scheme. This is easily done and we'll
illustrate it for the du FortFrankel scheme (4.3.8). For simplicity, consider an initialboundary value problem for (4.3.8) with homogeneous Dirichlet boundary conditions,
then the vector form of (4.3.8) is Un+1 = AUn;1 + BUn
where
2r
A = 1 ; 2r I
1+ 2
01
61 0 1
2
B = 1 +r2r 6
6
...
4
10 3
7
7
7
5 (4.3.9a)
2
6
Un = 6
6
4 U1n
U2n
...
n
UJ ;1 3
7
7:
7
5 (4.3.9b) 4.3. Multilevel Schemes
Rewrite (4.3.9a) as 21 Un+1 = B A
Un
I0 Un :
Un;1 In this form, the multilevel du FortFrankel scheme looks like the onelevel scheme Wn+1 = L Wn (4.3.9c) with Un
Wn = Un;1 B
L= AA
0 (4.3.9d) Thus, matrix stability analysis and techniques described in Section 3.3 are directly applicable to multilevel di erence methods using the matrix L of (4.3.9d). 4.3.2 The Method of Lines
Excellent software is available for solving ordinary di erential equations. Modern multistep methods use codes containing either Adams or backward di erence formulas 5,
12, 13] to integrate the system in time. Adams methods are preferred for nonsti problems while backward di erence methods are useful for sti problems. In order to use the
ordinary di erential equations software to solve partial di erential equations of the form ut = Lu (4.3.10a) we introduce a spatial grid (as shown in Figure 4.3.2) and replace all of the spatial
derivatives appearing in the operator L by nite di erence approximations to obtain
_
V = L xV (4.3.10b) where L x is the discrete approximation of L and (_) denotes time di erentiation. The
elements Vj (t) of V(t) are approximations of u(xj t). Let us illustrate the idea for a heat
conduction problem.
Example 4.3.1. The spatial operator for the heat conduction equation (4.1.1a) is
Lu = uxx. Discretizing the second spatial derivative using centered di erences gives
(t
_
Vj (t) = Vj;1(t) ; 2Vjx2 ) + Vj+1(t) : 22 Parabolic PDEs
t x, j
0 1 2 j J1 J Figure 4.3.2: Spatial discretization to be used with the method of lines.
Consider a problem with homogeneous Dirichlet boundary conditions and let V(t) =
V1(t) V2 (t) : : : Vj;1(t)]T , then the matrix form of this problem is
2
3
;2 1
6 1 ;2 1
7
7:
L x = x2 6
6
7
...
4
5
1 ;2
The above arguments and example emphasize that the partial di erential equation has
been \reduced," by spatial discretization, to a system of ordinary di erential equations.
The semidiscrete system (4.3.10b) can be integrated in time by most ordinary di erential
equations software. This software automatically selects time steps to satisfy a prescribed
(local) temporal error tolerance and to maintain stability. Most multistep codes also
adjust the order of accuracy of the formulas in order to improve performance. A user of
the software would only have to provide a temporal error tolerance, initial values V(0),
and a procedure for de ning the spatiallydiscrete operator L xV.
For heat conduction problems, the small divisor ( x2 ) present in L x suggests that
the ordinary di erential equations (4.3.10b) will be sti . These notes are not the appropriate place to initiate an extended discussion of \sti ness." Let us simply state that sti
systems have more stringent stability restrictions than nonsti systems. The accepted 4.3. Multilevel Schemes 23 S1 2
3
0 1 2/3 6/11
0 1 4/3 18/11
1/3 9/11
1
2/11
2
3
4
5 4
12/25
48/25
36/25
16/25
3/25 5
60/137
300/137
300/137
200/137
75/137
12/137 6
60/147
360/147
450/147
400/147
225/147
72/147
10/147 Table 4.3.1: Coe cients of backward di erence formuas (4.3.11) 11].
remedy is to use implicit ordinary di erential equations software and the backward difference codes seem to be adequate for most problems 5, 11, 13].
The simplest backward di erence method is the implicit or backward Euler method
_
Vn+1 = V...
View
Full Document
 Spring '14
 JosephE.Flaherty
 The Land, Boundary value problem, Partial differential equation, Dirichlet boundary condition, Neumann boundary condition, di erences

Click to edit the document details