Issue |
ESAIM: M2AN
Volume 43, Number 4, July-August 2009
Special issue on Numerical ODEs today
|
|
---|---|---|
Page(s) | 757 - 784 | |
DOI | https://doi.org/10.1051/m2an/2009026 | |
Published online | 08 July 2009 |
Analysis of the accuracy and convergence of equation-free projection to a slow manifold
1
Department of Mathematics,
University of Amsterdam,
Amsterdam, The Netherlands.
2
Modeling, Analysis and Simulation,
Centrum Wiskunde & Informatica,
Amsterdam, The Netherlands.
3
Department of Chemical Engineering,
Princeton University,
Princeton, NJ 08544, USA.
4
NEC Laboratories USA,
retired.
5
Department of Mathematics and Center for BioDynamics,
Boston University,
Boston, MA 02215, USA. tasso@math.bu.edu
6
Program in Applied and Computational Mathematics,
Princeton University,
Princeton, NJ 08544, USA.
Received:
10
July
2008
In [C.W. Gear, T.J. Kaper, I.G. Kevrekidis and A. Zagaris,
SIAM J. Appl. Dyn. Syst. 4 (2005) 711–732],
we developed
a class of iterative algorithms
within the context
of equation-free methods
to approximate
low-dimensional,
attracting,
slow manifolds
in systems
of differential equations
with multiple time scales.
For user-specified values
of a finite number
of the observables,
the mth member
of the class
of algorithms
()
finds iteratively
an approximation
of the appropriate zero
of the (m+1)st time derivative
of the remaining variables
and
uses this root
to approximate the location
of the point
on the slow manifold
corresponding to these values
of the observables.
This article
is the first
of two articles
in which
the accuracy and convergence
of the iterative algorithms
are analyzed.
Here,
we work directly
with fast-slow systems,
in which there is
an explicit small parameter, ε,
measuring the separation
of time scales.
We show that,
for each
,
the fixed point
of the iterative algorithm
approximates the slow manifold
up to and including
terms of
.
Moreover,
for each m,
we identify explicitly
the conditions
under which
the mth iterative algorithm
converges to this fixed point.
Finally,
we show that
when
the iteration
is unstable
(or
converges slowly)
it may be stabilized
(or
its convergence
may be accelerated)
by application
of the Recursive Projection Method.
Alternatively,
the Newton-Krylov
Generalized Minimal Residual Method
may be used.
In the subsequent article,
we will consider
the accuracy and convergence
of the iterative algorithms
for a broader class
of systems – in which
there need not be
an explicit
small parameter – to which
the algorithms also apply.
Mathematics Subject Classification: 35B25 / 35B42 / 37M99 / 65L20 / 65P99
Key words: Iterative initialization / DAEs / singular perturbations / legacy codes / inertial manifolds.
© EDP Sciences, SMAI, 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.