Abstract:Reaction–diffusion equations have been widely used to describe biological pattern formation. Nonuniform steady states of reaction–diffusion models correspond to stationary spatial patterns supported by these models. Frequently these steady states are not unique and correspond to various spatial patterns observed in biology. Traditionally, time-marching methods or steady state solvers based on Newton's method were used to compute such solutions. However, the solutions that these methods converge to highly depend on the initial conditions or guesses. In this paper, we present a systematic method to compute multiple nonuniform steady states for reaction–diffusion models and determine their dependence on model parameters. The method is based on homotopy continuation techniques and involves mesh refinement, which significantly reduces computational cost. The method generates one-parameter steady state bifurcation diagrams that may contain multiple unconnected components, as well as two-parameter solution maps that divide the parameter space into different regions according to the number of steady states. We applied the method to two classic reaction–diffusion models and compared our results with available theoretical analysis in the literature. The first is the Schnakenberg model which has been used to describe biological pattern formation due to diffusion-driven instability. The second is the Gray–Scott model which was proposed in the 1980s to describe autocatalytic glycolysis reactions. In each case, the method uncovers many, if not all, nonuniform steady states and their stabilities.

Analysis on Bifurcation and Stability of a Generalized Gray–Scott Chemical Reaction Model

Application Of Grey Theory To Fault Prediction Of Chemical Equipment System

Existence of Traveling Waves of General Gray-Scott Models

Bogdanov–Takens bifurcation of an enzyme-catalyzed reaction model

Steady-State Bifurcation and Spatial Patterns of a Chemical Reaction System

On a reversible Gray-Scott type system from energetic variational approach and its irreversible limit

Analysis and Simulation of a Nonlocal Gray–Scott Model

Composite Molecules and Decoupling in Reaction Diffusion Models

Dynamic analysis and pattern formation in a Generalized Klausmeier-Gray-Scott Model

Structural Bifurcation Analysis in Chemical Reaction Networks

Stochastic analysis of chemical reactions in multi-component interacting systems at criticality

Unified framework for localized patterns in reaction–diffusion systems; the Gray–Scott and Gierer–Meinhardt cases

Turing Instability and Dynamic Bifurcation for the One‐Dimensional Gray–Scott Model

Mathematical Analysis of Chemical Reaction Systems

Fractional Gray–Scott model: Well-posedness, discretization, and simulations

Stability analysis of the stochastic Grey‐Scott model using spectral method

A high-precision numerical approach to solving space fractional Gray-Scott model

Modeling of chemical reaction systems with detailed balance using gradient structures

Spatial pattern formation in reaction–diffusion models: a computational approach

Numerical Solutions of Certain New Models of the Time-Fractional Gray-Scott

Numerical simulation of the time fractional Gray-Scott model on 2D space domains using radial basis functions