W. Auzinger, R. Frank:
"Asymptotic expansions of the global discretization error for stiff problems";
SIAM Journal on Scientific Computing,
We discuss the existence of asymptotic expansions of the global discretization error for a general class of nonlinear stiff differential equations y'(t) = A(t)y(t) + \phi(t,y(t)) where A(t) has a `stiff spectrum' characterized by a small parameter \epsilon and where \phi(t,y) is smooth. The following methods are considered: Implicit Euler, implicit midpoint and trapezoidal rules.
In strongly stiff situations (\epsilon significantly smaller than the stepsize h) the implicit Euler scheme admits a full asymptotic expansion; the same is true for the midpoint rule and for the trapezoidal rule under certain coupling conditions. In those strongly stiff cases where a full expansion does not exist for the midpoint or trapezoidal rule, the remainder term is of a reduced order but shows a regular, oscillating behavior which we describe in detail. In mildly stiff situations, order reductions of the remainder term inevitably occur in any case after the start or after the change of stepsize but - as can be shown by discrete singular perturbation techniques - these order reductions are rapidly damped out as the integration proceeds. Our results are illustrated by various numerical examples; in particular, numerical experience with extrapolation and Defect Correction is reported.
"Offizielle" elektronische Version der Publikation (entsprechend ihrem Digital Object Identifier - DOI)
Erstellt aus der Publikationsdatenbank der Technischen Universitšt Wien.