Figure 1
Spatiotemporal pattern of reactant concentration
for
,
(left), and
(right). A concentration of
is represented in white (black), and time increases from top to bottom in each panel. Reactant
is initially distributed homogeneously,
, over the entire domain of length
. Species
is randomly seeded, indicated by the white rectangles in the top left and right panels, with a seed length of
, an initial concentration of
, and the constraint that at least two seeds are separated by a minimum length of
. Both realizations of initial conditions have similar extinction times,
(left) and
(right). The upper panels are calculated for a total simulation time of
, where the vertical bars represent
time units and the horizontal bars
space units. The extinction process during the last
of the simulation is shown in the corresponding lower panels, with a fivefold increase in temporal resolution, with the vertical bars representing
time units. An explicit Euler method (
), with a three-point approximation of the Laplacian operator (
), and periodic boundary conditions were used for the numerical integration of the one-dimensional array.
Reuse & Permissions