A singular perturbation analysis for the Brusselator

Maximilian Engel,Guillermo Olicón-Méndez
2023-12-18
Abstract:In this work we study the Brusselator - a prototypical model for chemical oscillations - under the assumption that the bifurcation parameter is of order $O(1/\epsilon)$ for positive $\epsilon\ll 1$. The dynamics of this mathematical model exhibits a time scale separation visible via fast and slow regimes along its unique attracting limit cycle. Noticeably this limit cycle accumulates at infinity as $\epsilon\rightarrow 0$, so that in polar coordinates $(\theta,r)$, and by doing a further change of variable $r\mapsto r^{-1}$, we analyse the dynamics near the line at infinity, corresponding to the set $\{r=0\}$. This object becomes a nonhyperbolic invariant manifold for which we use a desingularising rescaling, in order to study the closeby dynamics. Further use of geometric singular perturbation techniques allows us to give a decomposition of the Brusselator limit cycle in terms of four different fully quantified time scales.
Dynamical Systems
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the time - scale separation phenomenon in the system dynamics when the value of the bifurcation parameter \(b\) is much larger than \(a\), that is, \(b\gg a\) in the Brusselator model. Specifically, the paper focuses on the dynamic behavior of the Brusselator model under the conditions of \(b = \frac{a}{\epsilon}\) and \(\epsilon\ll1\), especially the properties of its attractor limit cycle. ### Main research questions: 1. **Time - scale separation**: The paper explores how the dynamics of the Brusselator model shows the separation of fast and slow dynamics under the above - mentioned parameter settings. This separation leads to relaxation oscillations of the system, that is, the system alternately performs fast and slow motions on different time scales. 2. **Properties of the limit cycle**: The paper analyzes in detail the limit cycle in the Brusselator model, especially its time scales in different regions. Through the Geometric Singular Perturbation Theory (GSPT) and the rescaling method, the author decomposes the four different time scales of the limit cycle and gives the exact expressions of the time scales of each branch. 3. **Non - hyperbolic invariant manifolds**: The paper also studies the non - hyperbolic invariant manifolds at infinity, especially the dynamic behavior at \(r = 0\). Through appropriate coordinate transformations and time rescaling, the author reveals the dynamic structures near these non - hyperbolic points. ### Specific questions: - **Approximate curves of the limit cycle**: The paper provides approximate curves of the limit cycle in different regions and analyzes the properties of these curves. - **Quantification of time scales**: The paper gives the quantified expressions of the time scales of the system in different regions, such as \(O(\epsilon^3)\), \(O(\epsilon^{- 1})\), \(O(1)\) and \(O(\epsilon^{-3/2})\). - **Treatment of non - hyperbolic points**: The paper uses the rescaling and "blow - up" methods to deal with non - hyperbolic points, thereby revealing the dynamic behavior near these points. ### Conclusion: Through detailed mathematical analysis and geometric methods, the paper successfully describes the dynamic behavior of the Brusselator model under large bifurcation parameters, especially the properties of its limit cycle and the time - scale separation phenomenon. These results not only deepen the understanding of the Brusselator model, but also provide an important reference for studying other chemical reaction systems with multiple time scales.