Page 340 -
P. 340

7.5. The Maximum Principle for the One-Step-Theta Method  323


        According to the above considerations, it is nonexpansive only if λτ ≤ 1
        holds. For large numbers λ, this is a very restrictive step size condition; see
        also the discussion of (7.85) to (7.87).
          Due to their better stability properties, implicit methods such as the
        Crank–Nicolson and the implicit Euler methods do not have such step
        size restrictions. Nevertheless, the application of implicit methods is not
        free from surprises. For example, in the case of large numbers λ, an order
        reduction can occur.



        Exercises

         7.10 Determine the corresponding domain of stability S R of the one-step-
                                                              1
        theta method for the following values of the parameter Θ : 0, , 1.
                                                              2
         7.11 Show the L-stability of the implicit Euler method.

         7.12  (a) Show that the discretization
                         n
                        ξ = ξ n−2  +2τf(t n−1 ,ξ n−1 ) ,  n =2,...N
             (midpoint rule), applied to the model equation ξ = f(t, ξ)with

             f(t, ξ)= −λξ and λ> 0 leads, for a sufficiently small step size
             τ> 0, to a general solution that can be additively decomposed into a
             decaying and an increasing (by absolute value) oscillating component.
         (b) Show that the oscillating component can be damped if additionally
             the quantity ξ N  is computed (modified midpoint rule):
                         ∗
                                                        '
                                                      N
                           ξ  N  =  1 & ξ N  + ξ  N−1  + τf(t N ,ξ ) .
                                 2
                            ∗
         7.13 Let m ∈ N be given. Find a polynomial R m (z)=1 + z +

          m     j
          j=2  γ j z (γ j ∈ R) such that the corresponding domain of absolute sta-
        bility for R(z):= R m (z) contains an interval of the negative real axis that
        is as large as possible.

        7.5 The Maximum Principle for the
               One-Step-Theta Method

        In Section 1.4 we have seen that for a discrete problem of the form (1.31)
        there is a hierarchy of properties ranging from a comparison principle to
        a strong maximum principle, which is in turn applied by a hierarchy of
        conditions, partly summarized as (1.32) or (1.32) . To remind the reader,
                                                   ∗
        we regroup these conditions accordingly:
   335   336   337   338   339   340   341   342   343   344   345