Analysis of the accuracy and convergence of equation-free projection to a slow manifold

TL;DR

This paper analyzes the accuracy and convergence of iterative algorithms designed to approximate slow manifolds in explicit fast-slow systems, demonstrating their effectiveness and stability under certain conditions.

Contribution

It provides explicit error bounds and convergence conditions for a class of iterative algorithms in explicit fast-slow systems, enhancing understanding of their accuracy and stability.

Findings

Fixed points approximate the slow manifold up to O(ε^m)

Conditions for convergence of the iterative algorithms are explicitly identified

Stability and convergence can be improved using the Recursive Projection Method or Newton-Krylov methods

Abstract

In [C.W. Gear, T.J. Kaper, I.G. Kevrekidis, and A. Zagaris, Projecting to a Slow Manifold: Singularly Perturbed Systems and Legacy Codes, 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 m-th member of the class of algorithms (m = 0, 1, ...) 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 explicit fast--slow systems, in which there is an explicit small parameter, epsilon, measuring the separation of time scales. We show that, for each m = 0, 1, ..., the fixed point of the iterative algorithm approximates the slow manifold up to and including terms of O(epsilon^m). Moreover, for each m, we identify explicitly the conditions under which the m-th 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.