Numerical continuation of boundaries in parameter space between stable and unstable periodic travelling wave (wavetrain) solutions of partial differential equations

Research output: Contribution to journalArticle

16 Citations (Scopus)

Abstract

A variety of numerical methods are available for determining the stability of a given solution of a partial differential equation. However for a family of solutions, calculation of boundaries in parameter space between stable and unstable solutions remains a major challenge. This paper describes an algorithm for the calculation of such stability boundaries, for the case of periodic travelling wave solutions of spatially extended local dynamical systems. The algorithm is based on numerical continuation of the spectrum. It is implemented in a fully automated way by the software package wavetrain, and two examples of its use are presented. One example is the Klausmeier model for banded vegetation in semi-arid environments, for which the change in stability is of Eckhaus (sideband) type; the other is the two-component Oregonator model for the photosensitive Belousov-Zhabotinskii reaction, for which the change in stability is of Hopf type.

Original languageEnglish
Pages (from-to)175-192
Number of pages18
JournalAdvances in Computational Mathematics
Volume39
Issue number1
DOIs
Publication statusPublished - 1 Jul 2013

Keywords

  • Numerical continuation
  • Periodic traveling wave
  • Wavetrain
  • Auto
  • Eckhaus
  • Hopf
  • Spectral stability
  • GINZBURG-LANDAU EQUATION
  • SEMIARID ENVIRONMENTS
  • BANDED VEGETATION
  • LINEAR-OPERATORS
  • PATTERN-FORMATION
  • KLAUSMEIER MODEL
  • DYNAMICS
  • SPECTRA
  • STABILITY
  • SYSTEMS

Cite this