31 right a von neumann stability analysis reveals

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: implied that stability restrictions were necessary for explicit nite-di 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 Fort-Frankel scheme (4.3.8) is consistent with the hyperbolic equation ut ; uxx + 2u tt = 0: Thus, the time step of the du Fort-Frankel 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 one-step method. We must be careful to preserve accuracy of the multistep scheme when the one-step method has a lower order. Suppose, for example, that the time step of a second-order multistep method is tM , then the time step of the rst-order 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 time-centered space scheme (4.1.2) is used as a starting method for the du Fort-Frankel (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 one-level scheme. This is easily done and we'll illustrate it for the du Fort-Frankel 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 Fort-Frankel scheme looks like the one-level 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 J-1 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 semi-discrete 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 spatially-discrete 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

Ask a homework question - tutors are online