Higher order corrections in the approximation of low-dimensional manifolds and the construction of simplified problems with the CSP method

Authors:
Mauro Valorani;Dimitris A. Goussis;Francesco Creta;Habib N. Najm
Affiliations:
Dipartimento di Meccanica e Aeronautica, University of Rome La Sapienza, Via Eudossiana 18, 00184 Rome, Italy;Agiou Georgiou 49, 26500 Patra, Greece;Dipartimento di Meccanica e Aeronautica, University of Rome La Sapienza, Via Eudossiana 18, 00184 Rome, Italy;Sandia National Laboratories, Livermore, CA 94551, United States
Venue:
Journal of Computational Physics
Year:
2005

Citing 2
Cited 5

Asymptotic Solution of Stiff PDEs with the CSP Method: The Reaction Diffusion Equation

SIAM Journal on Scientific Computing
Explicit time-scale splitting algorithm for stiff problems: auto-ignition of gaseous mixtures behind a steady shock

Journal of Computational Physics

An efficient iterative algorithm for the approximation of the fast and slow dynamics of stiff systems

Journal of Computational Physics
The G-Scheme: A framework for multi-scale adaptive model reduction

Journal of Computational Physics
Eigenvalues of the Jacobian of a Galerkin-Projected Uncertain ODE System

SIAM Journal on Scientific Computing
Approximation of slow and fast dynamics in multiscale dynamical systems by the linearized Relaxation Redistribution Method

Journal of Computational Physics
Multiscale Stochastic Preconditioners in Non-intrusive Spectral Projection

Journal of Scientific Computing

Quantified Score

Hi-index	31.46

Visualization

Abstract

In systems of stiff Ordinary Differential Equations (ODEs) both fast and slow time scales are encountered. The fast time scales are responsible for the development of low-dimensional manifolds on which the solution moves according to the slow time scales. In this paper, methodologies for constructing highly accurate (i) expressions describing the manifold, and (ii) simplified non-stiff equations governing the slow evolution of the solution on the manifold are developed, according to an iterative procedure proposed in the Computational Singular Perturbation (CSP) method. It is shown that the increasing accuracy achieved with each iteration is directly related to the time rates of change of the CSP vectors spanning the manifold along the solution trajectory. Here, an algorithm is presented which implements these calculations and is validated on the basis of two simple examples.