B-methods for the numerical solution of evolution problems with blow-up solutions Part I: Variation of the constant

Mélanie Beck, Martin J. Gander, Wing Hong Felix KWOK

Research output: Contribution to journalArticlepeer-review

2 Citations (Scopus)

Abstract

In the last two decades, the field of geometric numerical integration and structurepreserving algorithms has focused on the design of numerical methods that preserve properties of Hamiltonian systems, evolution problems on manifolds, and problems with highly oscillatory solutions. In this paper, we show that a different geometric property, namely, the blow up of solutions in finite time, can also be taken into account in the numerical integrator, giving rise to geometric methods we call B-methods. We give a first systematic approach for deriving such methods for scalar and systems of semi- and quasi-linear parabolic and hyperbolic partial differential equations. We show both analytically and numerically that B-methods have substantially better approximation properties than standard numerical integrators as the solution approaches the blow-up time.

Original languageEnglish
Pages (from-to)A2998-A3029
JournalSIAM Journal of Scientific Computing
Volume37
Issue number6
DOIs
Publication statusPublished - 2015

Scopus Subject Areas

  • Computational Mathematics
  • Applied Mathematics

User-Defined Keywords

  • Blow-up solutions
  • Geometric integration
  • Nonlinear partial differential equations
  • Nonlinear systems of equations

Fingerprint

Dive into the research topics of 'B-methods for the numerical solution of evolution problems with blow-up solutions Part I: Variation of the constant'. Together they form a unique fingerprint.

Cite this