On the numerical continuation of isolas of equilibria

D. Avitabile*, M. Desroches, S. Rodrigues

*Corresponding author for this work

Research output: Contribution to JournalArticleAcademicpeer-review

Abstract

We present a numerical strategy to compute one-parameter families of isolas of equilibrium solutions in ODEs. Isolas are solution branches closed in parameter space. Numerical continuation is required to compute one single isola since it contains at least one unstable segment. We show how to use pseudo-arclength predictor-corrector schemes in order to follow an entire isola in parameter space, as an individual object, by posing a suitable algebraic problem. We continue isolas of equilibria in a two-dimensional dynamical system, the so-called continuous stirred tank reactor model, and also in a three-dimensional model related to plasma physics. We then construct a toy model and follow a family of isolas past a fold and illustrate how to initiate the computation close to a formation center, using approximate ellipses in a model inspired by the van der Pol equation. We also show how to introduce node adaptivity in the discretization of the isola, so as to concentrate on nodes in region with higher curvature. We conclude by commenting on the extension of the proposed numerical strategy to the case of isolas of periodic orbits.

Original languageEnglish
Article number1250277
JournalInternational Journal of Bifurcation and Chaos
Volume22
Issue number11
DOIs
Publication statusPublished - 1 Jan 2012
Externally publishedYes

Keywords

  • bifurcation theory
  • isolas
  • Numerical continuation

Fingerprint Dive into the research topics of 'On the numerical continuation of isolas of equilibria'. Together they form a unique fingerprint.

Cite this