Computing the Stable Manifold of a Saddle Slow Manifold (original) (raw)
Related papers
Spike Adding and the Role of Saddle Slow Manifolds
2017
In this thesis, we study spike adding in slow-fast systems, i.e., systems in which some variables evolve significantly faster than others. In slow-fast systems, there exist so-called slow manifolds, along which the flow is much slower than the dynamics far from them. Slow manifolds, together with the invariant manifolds of saddle equilibria and periodic orbits organise the dynamics of slow-fast systems. The slow manifolds can be attracting, repelling or of saddle type, depending on the stability of their associated set of equilibria in the so-called fast subsystem when the slow variables are treated as parameters. In this thesis, we focus on slow manifolds of saddle type and their associated (un)stable manifolds in three-dimensional systems with one slow and two fast variables. In this case, the saddle slow manifold (SSM) exists as a (nonunique) one-dimensional curve and the associated stable and unstable manifolds are (nonunique) two-dimensional surfaces. Numerical methods for comp...
The Dynamics of Slow Manifolds
Journal of the Indonesian Mathematical Society, 2007
After reviewing a number of results from geometric singular perturbation theory, we discuss several approaches to obtain periodic solutions in a slow manifold. Regarding nonhyperbolic transitions we consider relaxation oscillations and canard-like solutions. The results are illustrated by prey-predator systems.
Discrete and Continuous Dynamical Systems, 2016
The present work deals with numerical methods for computing slow stable invariant manifolds as well as their invariant stable and unstable normal bundles. The slow manifolds studied here are sub-manifolds of the stable manifold of a hyperbolic equilibrium point. Our approach is based on studying certain partial differential equations equations whose solutions parameterize the invariant manifolds/bundles. Formal solutions of the partial differential equations are obtained via power series arguments, and truncating the formal series provides an explicit polynomial representation for the desired chart maps. The coefficients of the formal series are given by recursion relations which are amenable to computer calculations. The parameterizations conjugate the dynamics on the invariant manifolds and bundles to a prescribed linear dynamical systems. Hence in addition to providing accurate representation of the invariant manifolds and bundles our methods describe the dynamics on these objects as well. Example computations are given for vector fields which arise as Galerkin projections of a partial differential equation. As an application we illustrate the use of the parameterized slow manifolds and their linear bundles in the computation of heteroclinic orbits.
An Approach for approximating slow and fast dynamics of multi-scale dynamical systems
Journal of Computational Physics, 2011
In this paper, we introduce a fictitious dynamics for describing the only fast relaxation of a stiff ordinary differential equation (ODE) system towards a stable low-dimensional invariant manifold in the phase-space (slow invariant manifold – SIM). As a result, the demanding problem of constructing SIM of any dimensions is recast into the remarkably simpler task of solving a properly devised ODE system by stiff numerical schemes available in the literature. In the same spirit, a set of equations is elaborated for local construction of the fast subspace, and possible initialization procedures for the above equations are discussed. The implementation to a detailed mechanism for combustion of hydrogen and air has been carried out, while a model with the exact Chapman–Enskog solution of the invariance equation is utilized as a benchmark.
Rigorous Enclosures of a Slow Manifold
SIAM Journal on Applied Dynamical Systems, 2012
Slow-fast dynamical systems have two time scales and an explicit parameter representing the ratio of these time scales. Locally invariant slow manifolds along which motion occurs on the slow time scale are a prominent feature of slow-fast systems. This paper introduces a rigorous numerical method to compute enclosures of the slow manifold of a slow-fast system with one fast and two slow variables. A triangulated first order approximation to the two dimensional invariant manifold is computed "algebraically". Two translations of the computed manifold in the fast direction that are transverse to the vector field are computed as the boundaries of an initial enclosure. The enclosures are refined to bring them closer to each other by moving vertices of the enclosure boundaries one at a time. As an application we use it to prove the existence of tangencies of invariant manifolds in the problem of singular Hopf bifurcation and to give bounds on the location of one such tangency.
2 Rigorous Enclosures of a Slow Manifold
2016
Slow-fast dynamical systems have two time scales and an explicit parameter representing the ratio of these time scales. Locally invariant slow manifolds along which motion occurs on the slow time scale are a prominent feature of slow-fast systems. This paper introduces a rigorous numerical method to compute enclosures of the slow manifold of a slow-fast system with one fast and two slow variables. A triangulated first order approximation to the two dimensional invariant manifold is computed "algebraically". Two translations of the computed manifold in the fast direction that are transverse to the vector field are computed as the boundaries of an initial enclosure. The enclosures are refined to bring them closer to each other by moving vertices of the enclosure boundaries one at a time. As an application we use it to prove the existence of tangencies of invariant manifolds in the problem of singular Hopf bifurcation and to give bounds on the location of one such tangency.
Approximating stable and unstable manifolds in experiments
Physical Review E, 2003
We introduce a procedure to reveal invariant stable and unstable manifolds, given only experimental data. We assume a model is not available and show how coordinate delay embedding coupled with invariant phase space regions can be used to construct stable and unstable manifolds of an embedded saddle. We show that the method is able to capture the fine structure of the manifold, is independent of dimension, and is efficient relative to previous techniques.
Calculation of Slow Invariant Manifolds for Reactive Systems
47th AIAA Aerospace Sciences Meeting including The New Horizons Forum and Aerospace Exposition, 2009
One-dimensional slow invariant manifolds for dynamical systems arising from modeling unsteady, isothermal, isochoric, spatially homogeneous, closed reactive systems are calculated. The technique is based on global analysis of the composition space of the reactive system. The identification of all the system's finite and infinite critical points plays a major role in calculating the system's slow invariant manifold. The slow invariant manifolds are constructed by calculating heteroclinic orbits which connect appropriate critical points to the critical point which corresponds to the unique stable physical critical point of chemical equilibrium. The technique is applied to small and large detailed kinetics mechanisms for hydrogen combustion.
Analysis of the accuracy and convergence of equation-free projection to a slow manifold
ESAIM: Mathematical Modelling and Numerical Analysis, 2009
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