An efficient Newton-Krylov implementation of the constrained runs scheme for initializing on a slow manifold

Christophe Vandekerckhove, Ioannis Kevrekidis, Dirk Roose

Research output: Contribution to journalArticlepeer-review

14 Scopus citations


The long-term dynamic behavior of many dynamical systems evolves on a low-dimensional, attracting, invariant slow manifold, which can be parameterized by only a few variables ("observables"). The explicit derivation of such a slow manifold (and thus, the reduction of the long-term system dynamics) is often extremely difficult or practically impossible. For this class of problems, the equation-free framework has been developed to enable performing coarse-grained computations, based on short full model simulations. Each full model simulation should be initialized so that the full model state is consistent with the values of the observables and close to the slow manifold. To compute such an initial full model state, a class of constrained runs functional iterations was proposed (Gear and Kevrekidis, J. Sci. Comput. 25(1), 17-28, 2005; Gear et al., SIAM J. Appl. Dyn. Syst. 4(3), 711-732, 2005). The schemes in this class only use the full model simulator and converge, under certain conditions, to an approximation of the desired state on the slow manifold. In this article, we develop an implementation of the constrained runs scheme that is based on a (preconditioned) Newton-Krylov method rather than on a simple functional iteration. The functional iteration and the Newton-Krylov method are compared in detail using a lattice Boltzmann model for one-dimensional reaction-diffusion as the full model simulator. Depending on the parameters of the lattice Boltzmann model, the functional iteration may converge slowly or even diverge. We show that both issues are largely resolved by using the Newton-Krylov method, especially when a coarse grid correction preconditioner is incorporated.

Original languageEnglish (US)
Pages (from-to)167-188
Number of pages22
JournalJournal of Scientific Computing
Issue number2
StatePublished - May 2009

All Science Journal Classification (ASJC) codes

  • Software
  • General Engineering
  • Computational Mathematics
  • Theoretical Computer Science
  • Applied Mathematics
  • Numerical Analysis
  • Computational Theory and Mathematics


  • Constrained runs
  • Equation-free computing
  • Initialization
  • Lattice Boltzmann model
  • Newton-Krylov method
  • Slow manifold


Dive into the research topics of 'An efficient Newton-Krylov implementation of the constrained runs scheme for initializing on a slow manifold'. Together they form a unique fingerprint.

Cite this