Computational singular perturbation with non-parametric tabulation of slow manifolds for time integration of stiff chemical kinetics

Bert J. Debusschere, Youssef M. Marzouk, Habib N. Najm, Blane Rhoads, Dimitris A. Goussis, Mauro Valorani

Research output: Contribution to journalArticlepeer-review

17 Scopus citations


This paper presents a novel tabulation strategy for the adaptive numerical integration of chemical kinetics using the computational singular perturbation (CSP) method. The strategy stores and reuses CSP quantities required to filter out fast dissipative processes, resulting in a non-stiff chemical source term. In particular, non-parametric regression on low-dimensional slow invariant manifolds (SIMs) in the chemical state space is used to approximate the CSP vectors spanning the fast chemical subspace and the associated fast chemical time-scales. The relevant manifold and its dimension varies depending on the local number of exhausted modes at every location in the chemical state space. Multiple manifolds are therefore tabulated, corresponding to different numbers of exhausted modes (dimensions) and associated radical species. Non-parametric representations are inherently adaptive, and rely on efficient approximate-nearest-neighbor queries. As the CSP information is only a function of the non-radical species in the system and has relatively small gradients in the chemical state space, tabulation occurs in a lower-dimensional state space and at a relatively coarse level, thereby improving scalability to larger chemical mechanisms. The approach is demonstrated on the simulation of homogeneous constant pressure H 2-air and CH 4-air ignition, over a range of initial conditions. For CH 4-air, results are shown that outperform direct implicit integration of the stiff chemical kinetics while maintaining good accuracy.

Original languageBritish English
Pages (from-to)173-198
Number of pages26
JournalCombustion Theory and Modelling
Issue number1
StatePublished - Feb 2012


  • chemical kinetics
  • computational singular perturbation
  • kd-trees
  • nearest neighbors
  • non-parametric regression
  • slow manifold


Dive into the research topics of 'Computational singular perturbation with non-parametric tabulation of slow manifolds for time integration of stiff chemical kinetics'. Together they form a unique fingerprint.

Cite this