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

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.

References in corpus (1)

• Asymptotic analysis of two reduction methods for systems of chemical reactions

Cited by in corpus (25)

• Implicit Methods for Equation-Free Analysis: Convergence Results and Analysis of Emergent Waves in Microscopic Traffic Models
• Minimal Curvature Trajectories: Riemannian Geometry Concepts for Model Reduction in Chemical Kinetics
• Data-driven Control of Agent-based Models: an Equation/Variable-free Machine Learning Approach
• Slow manifold reduction for plasma science
• A variational principle for computing slow invariant manifolds in dissipative dynamical systems
• Slow Invariant Manifolds of Slow-Fast Dynamical Systems
• Guiding center dynamics as motion on a formal slow manifold in loop space
• Slow Invariant Manifolds of Singularly Perturbed Systems via Physics-Informed Machine Learning
• A physics-informed neural network method for the approximation of slow invariant manifolds for the general class of stiff systems of ODEs
• Two novel families of multiscale staggered patch schemes efficiently simulate large-scale, weakly damped, linear waves
• Exponential estimates of symplectic slow manifolds
• A Numerical Slow Manifold Approach to Model Reduction for Optimal Control of Multiple Time Scale ODE
• On fundamental unifying concepts for trajectory-based slow invariant attracting manifold computation in multiscale models of chemical kinetics
• Towards differential geometric characterization of slow invariant manifolds in extended state space: Sectional Curvature and Flow Invariance
• Simulating Stochastic Inertial Manifolds by a Backward-Forward Approach
• Convergence of equation-free methods in the case of finite time scale separation with application to deterministic and stochastic systems
• On convergence of higher order schemes for the projective integration method for stiff ordinary differential equations
• Characterization of Separatrices in Holomorphic Dynamical Systems
• A Numerical Method for the Parametrization of Stable and Unstable Manifolds of Microscopic Simulators
• An iterative method for the approximation of fibers in slow-fast systems
• Towards an efficient multiscale modeling of low-dimensional reactive systems: study of numerical closure procedures
• A micro/macro algorithm to accelerate Monte Carlo simulation of stochastic differential equations
• Computation of saddle type slow manifolds using iterative methods
• An online slow manifold approach for efficient optimal control of multiple time-scale kinetics
• Considering Slow Manifold Based Model Reduction for Multiscale Chemical Optimal Control Problems