Stiffness and Multiste Methods Chapter 26刚度和多步方法章26.pptVIP

by Lale Yurttas, Texas AM University Chapter 26 Stiffness and Multistep MethodsChapter 26 Two areas are covered: Stiff ODEs will be described - ODEs that have both fast and slow components to their solution. Implicit solution technique and multistep methods will be described. Stiffness A stiff system is the one involving rapidly changing components together with slowly changing ones. Both individual and systems of ODEs can be stiff: If y(0)=0, the analytical solution is developed as: Figure 26.1 Insight into the step size required for stability of such a solution can be gained by examining the homogeneous part of the ODE: The solution starts at y(0)=y0 and asymptotically approaches zero. If Euler’s method is used to solve the problem numerically: The stability of this formula depends on the step size h: Thus, for transient part of the equation, the step size must be 2/1000=0.002 to maintain stability. While this criterion maintains stability, an even smaller step size would be required to obtain an accurate solution. Rather than using explicit approaches, implicit methods offer an alternative remedy. An implicit form of Euler’s method can be developed by evaluating the derivative at a future time. Figure 26.2 Figure 26.3 Multistep Methods The Non-Self-Starting Heun Method/ Huen method uses Euler’s method as a predictor and trapezoidal rule as a corrector. Predictor is the weak link in the method because it has the greatest error, O(h2). One way to improve Heun’s method is to develop a predictor that has a local error of O(h3). Figure 26.4 Step-Size Control/ Constant Step Size. A value for h must be chosen prior to computation. It must be small enough to yield a sufficiently small truncation error. It should also be as large as possible to minimize run time cost and round-off error. Variable Step Size. If the corrector error is greater than some specified error, the step size is decreased. A step size is chosen so that the convergence cr