A multiparameter singular perturbation analysis of the Robertson model

Abstract

The Robertson model describing a chemical reaction involving three reactants is one of the classical examples of stiffness in ODEs. The stiffness is caused by the occurrence of three reaction rates 𝑘₁, 𝑘₂, and 𝑘₃, with largely differing orders of magnitude, acting as parameters. The model has been widely used as a numerical test problem. Surprisingly, no asymptotic analysis of this multiscale problem seems to exist. In this paper, we provide a full asymptotic analysis of the Robertson model under the assumption 𝑘₁, 𝑘₃ ≪ 𝑘₂. We rewrite the equations as a two-parameter singular perturbation problem in the rescaled small parameters (𝜀₁, 𝜀₂) ∶= (𝑘₁∕𝑘₂, 𝑘₃∕𝑘₂), which we then analyze using geometric singular perturbation theory (GSPT). To deal with the multiparameter singular structure, we perform blowups in parameter- and variable space. We identify four distinct regimes in a neighborhood of the singular limit (𝜀₁, 𝜀₂) = (0, 0). Within these four regimes, we use GSPT and additional blowups to analyze the dynamics and the structure of solutions. Our asymptotic results are in excellent qualitative and quantitative agreement with the numerics.