right-parenthesis w Superscript n Baseline greater-than-or-equal-to 0 comma n equals 0 comma ellipsis comma upper N minus 1 comma w Superscript 0 Baseline greater-than-or-equal-to 0"/>
implies that
Based on this definition, we see that the implicit Euler scheme is always positive while the explicit Euler scheme is positive if the term:
is positive. Thus, if the function
Definition 2.2 A difference scheme is stable if its solution is based in much the same way as the solution of the continuous problem (2.1) (see Theorem 2.1), that is:
(2.19)
where:
and:
Based on the fact that a scheme is stable and consistent (see Dahlquist and Björck (1974)), we can state in general that the error between the exact and discrete solutions is bounded by some polynomial power of the step-size
(2.20)
where
(2.21)
Thus, we see that the Box method is second-order accurate and is better than the implicit Euler scheme, which is only first-order accurate.
2.4 SPECIAL SCHEMES
We introduce exponentially fitted schemes that are used for boundary layer problems (for example, convection-dominated PDEs) and the extrapolation method to increase the accuracy of finite difference schemes. We shall see how to apply these techniques to more complex problems in later chapters. We also discuss predictor-corrector methods.
2.4.1 Exponential Fitting
We now introduce a special class of schemes with desirable properties. These are schemes that are suitable for problems with rapidly increasing or decreasing solutions. In the literature these are called stiff or singular perturbation problems (see Duffy (1980)). We can motivate these schemes in the present context. Let us take the problem (2.1) when
If a is large then the derivatives of
(2.23)
