Spectral Deferred Correction Methods
Anita Layton and Michael Minion
|
Many physical systems with dynamics that involve
two or more processes with widely-differing characteristic
time scales are of interest to researchers
in the physical and biological sciences.
Well-known examples include combustion,
the transport of air pollutants, and the
movement of contaminants or microorganisms
in ground water systems. In each of these
examples, the mathematical models used to describe
the dynamics consist of systems of partial differential
equations (PDEs), which specify the advection,
diffusion, and reaction of chemical species
within a moving medium.
Hence, the construction of accurate, stable, and
efficient numerical methods for the solution of
advection-diffusion-reaction (A-D-R) equations
is a topic of significant current interest.
In this study, high-order multi-implicit
spectral deferred correction (MISDC) methods are presented for solving
A-D-R equations. MISDC methods are similar to operating-splitting methods
both types of methods give rise
to multiple implicit equations, for which different time steps may be used.
In theory, however, the temporal order of accuracy of MISDC methods can be
arbitrarily high because both the integration and splitting errors
are eliminated during the deferred correction process.
MISDC methods, which are a generalization of semi-implicit spectral deferred
correction
(SISDC) methods introduced by Minion (which are in turn modifications
of the explicit and implicit spectral deferred correction methods previously
developed by Dutt et al.),
use a low-order numerical method to compute a high-order approximation.
This is achieved by using the low-order numerical method to solve a series of
correction equations, each of which increases the order of accuracy
of the approximation by one.
The accuracy and stability of MISDC methods for A-D-R equations
and its efficiency relative to SISDC methods are investigated in this study.
The above figure shows results of an efficiency comparison between MISDC
methods and semi-implicit additive Runge-Kutta (ARK) methods
developed by Kennedy et al.,
using a simple, one-dimensional model of flamelets, with a large reaction
coefficient, as a test problem.
These results indicate that for problems with stiff reactions,
MISDC methods can be constructed to be competitive with ARK methods.
|
Related Publications
- Anne Bourlioux, Anita T. Layton, and Michael L. Minion,
``
High-Order Multi-Implicit Spectral Deferred Correction Methods
for Problems of Reactive Flows,''
Journal of Computational Physics, Vol. 189, No. 2, pp. 651-675, 2003.
- Anita T. Layton and Michael L. Minion,
``
Conservative multi-implicit spectral deferred correction methods
for reacting gas dynamics,''
Journal of Computational Physics, in press.
- Michael L. Minion,
``
Spectral deferred correction and projection method,''
to appear in proceedings for Workshop for Innovative Time Integrators for PDEs.
- Michael L. Minion,
``
Semi-implicit spectral deferred correction methods for ordinary
differential equations,''
Communications in Mathematical Sciences, in press, 2003.
- Michael L. Minion,
``
Higher-order semi-implicit projection methods,''
Lawrence Livermore National Laboratory,
Technical Report, UCRL-JC-145295, 2001.
|
Last updated: September 29, 2003.
Mail comments to Anita Layton.
|
|