Numerical Challenges for Resolving Spike Dynamics for Two One‐Dimensional Reaction‐Diffusion Systems

WT Sun,T Tang,MJ Ward,JC Wei
DOI: https://doi.org/10.1111/1467-9590.t01-1-00227
2003-01-01
Studies in Applied Mathematics
Abstract:Asymptotic and numerical methods are used to highlight different types of dynamical behaviors that occur for the motion of a localized spike‐type solution to the singularly perturbed Gierer–Meinhardt and Schnakenberg reaction‐diffusion models in a one‐dimensional spatial domain. Depending on the parameter range in these models, there can either be a slow evolution of a spike toward the midpoint of the domain, a sudden oscillatory instability triggered by a Hopf bifurcation leading to an intricate temporal oscillation in the height of the spike, or a pulse‐splitting instability leading to the creation of new spikes in the domain. Criteria for the onset of these oscillatory and pulse‐splitting instabilities are obtained through asymptotic and numerical techniques. A moving‐mesh numerical method is introduced to compute these different behaviors numerically, and results are compared with corresponding results computed using a method of lines based software package.
What problem does this paper attempt to address?