WOLFRAM|DEMONSTRATIONS PROJECT

Evolution in a Cellular-Automaton Model of Gray-Scott Reaction-Diffusion System

​
field size
20
30
40
pattern
spot I
spot II
reset
show field
U
V
run/stop simulation
This Demonstration shows the spatiotemporal evolution of a common spot pattern emerging from the numerical simulation of a simple reaction-diffusion system described by the Gray–Scott model.
The Gray–Scott model is defined by the following two irreversible reactions between two substances that give an inert final product. Let
U
and
V
be the concentration fields of the two substances and
P
that of the final product.

U+2V→3V
V→P
The system of coupled reaction-diffusion equations is given by

∂U
∂t
=
D
U
2
∇
U-U
2
V
+F(1-U)
∂V
∂t
=
D
V
2
∇
V+U
2
V
-(F+k)V
.
Here
D
U
and
D
V
are the respective diffusion coefficients in dimensionless units,
F
is the dimensionless feed rate, and
k
is the dimensionless rate constant entering the second equation (see [1]).
This Demonstration is an implementation of a cellular automaton simulation of the forward Euler integration of the coupled partial differential equations after stencil-type discretization of the Laplacian appearing in the diffusion operators. Von Neumann boundary conditions are implemented according to [2]. At the beginning of the simulation, the entire system is placed in the trivial state with
(U=1,V=0)
.
A square mesh grid is used with the center of the grid perturbed to
(U=0.5,V=0.25)
. Hence, the lattice for the
V
field appears as a central red square on a blue background, according to the chosen color legend. The grid values are perturbed with random noise in order to break the square symmetry of the lattice. The system is then integrated for several time steps until a spatial pattern emerges for the
V
field.
With the parameter set used in the initialization code of this simulation
(
D
U
=0.18,
D
V
=0.08,F=0.035,k=0.065)
,
you can observe the time evolution of the central spot, which initially increases in size, growing radially outward, and then divides into two spots. Then these two spots get further apart, elongate and further split in a way visually similar to cell division. The process goes on until the field is filled with spots.
This spot multiplication occurs in response to the finite-amplitude perturbations applied to a nonlinear system. It is an example of Turing instability patterns, which often occur in natural systems such as animal coats (e.g. zebra stripes). Some of these can be found in [1, Figure 2], which may require hundreds of thousands of time integration steps to be observed.
In order to see the spot patterns in this Demonstration, first choose the field size of the cellular automaton and a pattern (spot type I or II) by selecting the appropriate control at the top of the Manipulate panel. Then click the Play button to start running the simulation. You can visualize either the
U
or
V
field during the animation, but changes to the field size or pattern type are disabled.
To restart the simulation or change pattern type or field size, first stop the running simulation with the Pause control and check the box to reset the system to the initial state configuration.