The role of slow system dynamics in predicting the degeneracy of slow invariant manifolds: The case of vdP relaxation-oscillations

Abstract The development of dissipative time scales, which are much faster than the rest, force the solution of a system of 1st order ODEs to land and then evolve on a slow invariant manifold. The slow evolution on this manifold can be approximated by a reduced system, which is constructed with asymptotic analysis and is free of the fast scales. In this work, the dynamics of the reduced system are compared to the fast and slow dynamics of the original system. The comparison is based on the van der Pol system, which is typical of those exhibiting relaxation–oscillation behavior and is characterized by successive dissipative and explosive phases. Special attention is given to the dissipative phase, during which the solution evolves on a slow invariant manifold and the fast/slow time scale gap tends to contract. It is shown that one cannot predict from the dynamics of the van der Pol system the approaching explosive phase and the associated degeneracy of the manifold. Such a prediction is possible only through the dynamics of the reduced system, which involves an explosive component that is absent in the original system. This explosive component accelerates as higher order effects are incorporated in the reduced system, resulting in a smaller gap between the fast time scale of the original system and the fast time scale of the reduced one. This feature sets a limit on the accuracy singular perturbation techniques can deliver, even in the presence of a reasonably large fast/slow time scale gap in the dynamics of the original system. Given a fast/slow system, predicting forthcoming changes in the dimensions of the manifold is very important when it is desired to identify the underlying major physical mechanisms or to control the process.

[1]  Christopher Jones,et al.  Geometric singular perturbation theory , 1995 .

[2]  S. Lam,et al.  The CSP method for simplifying kinetics , 1994 .

[3]  Vladimir Gol'dshtein,et al.  Comparative analysis of two asymptotic approaches based on integral manifolds , 2004 .

[4]  Habib N. Najm,et al.  A CSP and tabulation-based adaptive chemistry model , 2007 .

[5]  Ioannis G. Kevrekidis,et al.  Constraint-Defined Manifolds: a Legacy Code Approach to Low-Dimensional Computation , 2005, J. Sci. Comput..

[6]  Ulrich Maas,et al.  Simplifying chemical kinetics: Intrinsic low-dimensional manifolds in composition space , 1992 .

[7]  Dimitris A. Goussis,et al.  Quasi steady state and partial equilibrium approximations: their relation and their validity , 2012 .

[8]  Dimitris A. Goussis,et al.  Physical understanding of complex multiscale biochemical models via algorithmic simplification: Glycolysis in Saccharomyces cerevisiae , 2010 .

[9]  Mauro Valorani,et al.  An efficient iterative algorithm for the approximation of the fast and slow dynamics of stiff systems , 2006, J. Comput. Phys..

[10]  Hans G. Kaper,et al.  Two perspectives on reduction of ordinary differential equations , 2005 .

[11]  S. H. Lam,et al.  Understanding complex chemical kinetics with computational singular perturbation , 1989 .

[12]  C. Kuehn.  A mathematical framework for critical transitions: Bifurcations, fast–slow systems and stochastic dynamics , 2011, 1101.2899.

[13]  N. Levinson,et al.  Small Periodic Pertubations of an Autonomous System with a Stable Orbit , 1950 .

[14]  Peter Szmolyan,et al.  Extending Geometric Singular Perturbation Theory to Nonhyperbolic Points - Fold and Canard Points in Two Dimensions , 2001, SIAM J. Math. Anal..

[15]  S. H. Lam,et al.  Conventional asymptotics and computational singular perturbation for simplified kinetics modelling , 1991 .

[16]  O. Rössler.  An equation for continuous chaos , 1976 .

[17]  H. Weitzner,et al.  Perturbation Methods in Applied Mathematics , 1969 .

[18]  Habib N. Najm,et al.  Higher order corrections in the approximation of low-dimensional manifolds and the construction of simplified problems with the CSP method , 2005 .

[19]  M. Krupa,et al.  Relaxation Oscillation and Canard Explosion , 2001 .

[20]  S. H. Lam,et al.  Using CSP to Understand Complex Chemical Kinetics , 1993 .

[21]  N. Semenoff,et al.  Zur Theorie des Verbrennungsprozesses , 1928 .

[22]  Ferdinand Verhulst,et al.  Periodic Solutions and Slow Manifolds , 2007, Int. J. Bifurc. Chaos.

[23]  Habib N. Najm,et al.  CSP analysis of a transient flame-vortex interaction: time scales and manifolds , 2003 .

[24]  H. Najm,et al.  Reactive and reactive-diffusive time scales in stiff reaction-diffusion systems , 2005 .

[25]  Frank W. J. Olver,et al.  Introduction to Asymptotics and Special Functions , 1974 .

[26]  P. A. Lagerstrom,et al.  Matched Asymptotic Expansions , 1988 .

[27]  L. A. Segel,et al.  The Quasi-Steady-State Assumption: A Case Study in Perturbation , 1989, SIAM Rev..

[28]  Carroll,et al.  Synchronization in chaotic systems. , 1990, Physical review letters.

[29]  S. Orszag,et al.  Advanced Mathematical Methods For Scientists And Engineers , 1979 .

[30]  Prodromos Daoutidis,et al.  Model reduction and control of multi-scale reaction-convection processes , 2008 .

[31]  Habib N. Najm,et al.  Analysis of methane–air edge flame structure , 2010 .

[32]  A. B. Vasil’eva.  ASYMPTOTIC BEHAVIOUR OF SOLUTIONS TO CERTAIN PROBLEMS INVOLVING NON-LINEAR DIFFERENTIAL EQUATIONS CONTAINING A?SMALL PARAMETER MULTIPLYING THE HIGHEST DERIVATIVES , 1963 .

[33]  W. D. Evans,et al.  PARTIAL DIFFERENTIAL EQUATIONS , 1941 .

[34]  Iliya V. Karlin,et al.  Method of invariant manifold for chemical kinetics , 2003 .

[35]  S. H. Lam,et al.  A study of homogeneous methanol oxidation kinetics using CSP , 1992 .

[36]  Habib N. Najm,et al.  Model Reduction and Physical Understanding of Slowly Oscillating Processes: The Circadian Cycle , 2006, Multiscale Model. Simul..

[37]  Neil Fenichel.  Geometric singular perturbation theory for ordinary differential equations , 1979 .

[38]  Ulrich Maas,et al.  Model Reduction for Combustion Chemistry , 2011 .

[39]  N. K. Rozov,et al.  Differential Equations with Small Parameters and Relaxation Oscillations , 1980 .

[40]  P. Holmes,et al.  Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields , 1983, Applied Mathematical Sciences.

[41]  Hans G. Kaper,et al.  Fast and Slow Dynamics for the Computational Singular Perturbation Method , 2004, Multiscale Model. Simul..

[42]  Terese Løvås,et al.  Model Reduction Techniques for Chemical Mechanisms , 2012 .

[43]  F. Verhulst.  Methods and applications of singular perturbations , 2005 .

[44]  V. N. Bogaevski,et al.  Algebraic methods in nonlinear perturbation theory , 1991 .

[45]  J. Grasman.  Asymptotic Methods for Relaxation Oscillations and Applications , 1987 .

[46]  Marc R. Roussel,et al.  Geometry of the steady-state approximation: Perturbation and accelerated convergence methods , 1990 .

[47]  Kaspar Nipp,et al.  Invariant manifolds of singularly perturbed ordinary differential equations , 1985 .

[48]  Dimitris A. Goussis,et al.  Asymptotic Solution of Stiff PDEs with the CSP Method: The Reaction Diffusion Equation , 1998, SIAM J. Sci. Comput..