Pattern formations in two-dimensional Gray–Scott model: existence of single-spot solutions and their stability
Introduction
Recently, there is a great interest in the study of self-replicating patterns observed in the irreversible Gray–Scott model [14], [15], [16] which governs the chemical reactions and in a gel reactor, where and are two chemical species, catalyzes its own reaction with and an inert product. Letting U=U(x,t) and V=V(x,t) denote the concentrations of and , respectively, the pair of coupled reaction–diffusion equations governing these reactions iswhere F denotes the rate at which is fed from the reservoir into the reactor, the concentration of in the reservoir is assumed to be zero, and k the rate at which is converted to an inert product. Here the diffusivities, DU and DV, of and , respectively, can be any chemically relevant positive numbers.
In 1993, Pearson [37] did some numerical simulations on the Gray–Scott model in a square of size 2.5 in R2 with periodic boundary conditions. By choosing DU=2×10−5, DV=10−5 and varying the parameters F and k, several interesting patterns were discovered. One of them is that the spot may self-replicate in a self-sustaining fashion and develop into a variety of time-dependent and time-independent asymptotic states. Lin et al. [28] reported their experiments in a ferro-cyanide–iodate–sulfite reaction which showed strong qualitative agreement with the self-replication regimes in simulations of [37]. Moreover, those same experiments led to the discovery of other new patterns, such as annular patterns emerging from circular spots. See [27] for more details on the set-up. (However, no stable single-spot was observed.) In 1D, numerical simulations were done by Reynolds et al. [38], [39], independently by Petrov et al. [36], and again self-replication phenomena were observed. However, in 1D, self-replication patterns were observed when DU=1, DV=δ2=0.01. Some formal asymptotics and dynamics in 1D were contained in [36], [38]. Recent numerical simulations of [9] in 1D and [29], [33] in 2D show that the single spot may be stable in some very narrow parameter regimes.
The first rigorous result in constructing a single-spot solution was due to Doelman et al. [9] in 1997. In [9], by using Mel’nikov method, Doelman et al. rigorously constructed single- and multiple-pulse solutions for (1.1) in the case N=1, DU=1, DV=δ2≪1. In their paper [9], it is assumed that F∼δ2, F+k∼δ2α/3, where . In this case, they showed that U=O(δα), V=O(δ−α/3). Later, the stability of single- and multi-pulse solutions in 1D are obtained in [6], [7], [8]. Hale et al. [21], [22] studied the case DU=DV in 1D and the existence of single- and multiple-pulse solutions are established. Nishiura and Ueyama [33] proposed a skeleton structure of self-replicating dynamics. Some related results on the existence and stability of solutions to the Gray–Scott model in 1D can be found in [11], [12], [38].
In higher dimension, as far as the author knows, there are no rigorous result on the existence or stability of spotty solutions for (1.1). Since self-replicating patterns were first observed in R2 [37], it is important and interesting to have a rigorous study on spotty solutions in R2. We remark that in R2 and R3, Muratov and Osipov [29] have given some formal asymptotic analysis on the construction and stability of spot solution. In [47], the author studied (1.1) in a bounded domain for the shadow system case, namely, DU≫1, DV≪1 and F=O(1), F+k=O(1). The shadow system can be reduced to a single equation. For spot solutions for single equations, see [1], [2], [4], [13], [17], [19], [20], [24], [25], [30], [31], [32], [34], [35], [40], [41], [42], [43], [44], [45], [46], [48], [49], and the references therein.
In this paper, we rigorously establish the existence and stability of single-spot solution in 2D. Let us first transform the system (1.1). We follow the notations in [29]. SetLet us drop the bar from now on. It is easy to see that (1.1) is equivalent to the following:Note that there are three parameters (α,ϵ,A) in Eq. (1.2). Throughout this paper, we always assume thatBy definition, we must haveTo study (1.2), we first consider the stationary equation of (1.2):From (1.1), it is also required that u and v are positive.
Before we state our results, let us first gain some insights into Eq. (1.5) and therefore introduce an important quantity. We are interested in radially symmetric solutions to (1.5). By (1.3), u is a slow variable and v a fast variable. Hence near 0, formally one can assume that u(x)∼u(0). Substituting that into equation for v, we have that v satisfiesBy the following scaling:we see that v̂ satisfiesProblem (1.7) has been studied in [18], [23]. It was proved that there exists a unique radially symmetric solution to the following problem: Remark 1.1 Unlike 1D case, the solution to (1.8) cannot be written explicitly. Let us denote the unique solution of (1.8) by w(y)=w(∣y∣). Going back to the equation for v, we see thatFrom the equation for u, we have thatwhere K(∣z∣) is the fundamental solution of −Δ+1 in R2. It is well known thatBy , , we obtainwhich suggests that the following quantity play a very important role:Thus formally, we have by the definition of L in (1.13)Eq. (1.14) admits a positive solution for u(0) if and only if . It is easy to see that for each , (1.14) admits two solutionsThus formally, we should have at least two solutions — one is small and another is large. This is exactly what we will prove in this paper.
Section snippets
Main results: existence and stability
We now state our main results in this paper. The existence result will be stated in terms of the quantity L introduced in (1.13) and the stability result will involve both L and α. Note that both L and α may depend on ϵ.
LetOur first result is on the existence of single-spot solution to (1.5). Theorem 2.1 Assume thatThen problem (1.5) admits two solutions (vϵi,uϵi), i=1,2 with the following properties: vϵi,uϵi are radially symmetric functions. vϵi=(1/Auϵi(0))(1+o(1))w(∣x∣/ϵ),i=1,2,
Outline of proofs and two technical lemmas
We prove Theorem 2.1 by the following steps. Step 3.1 We rewrite the system (1.5) as a single equation with a nonlocal term (rescaled):where T[v] solves the following equation:where ξϵ=ξϵ1 or ξϵ2. We will show that T is a well-defined operator. Step 3.2 We construct solutions of (3.1) of the following form:where φϵ is suitably small in some Sobolev spaces. To this end, it is important to study the linearized operatorAs ϵ
Proof of Theorem 2.1
In this section, we prove Theorem 2.1. We construct two single-spot solutions in R2 for the following Gray–Scott system:First we rescale Eq. (4.1). Recall that there are two solutions ξϵ1<ξϵ2 to Eq. (1.14). Moreover, it is important to note that by (2.2),From now on, we take ξϵ≔ξϵ1. The other case can be treated similarly. We rescaleThen will satisfy
Spectrum analysis I: reduction to NLEP and large eigenvalue estimates
Let be the two solutions constructed in Section 2. Without loss of generality, we use (vϵ,uϵ) to denote any one of the solutions. We now study the eigenvalue problem associated with (vϵ,uϵ). We assume that α=α0ϵγ for some 0≤γ<2 and positive number α0.
We need to analyze the following eigenvalue problem (letting x=ϵy):Let , . Problem (5.1) becomes
Spectrum analysis II: study of small eigenvalues
In this section, we discuss Case II: and finish Proof of Theorem 2.2(3). In this case, we need to show that and the multiplicity is exactly 2. We only need to consider (vϵ,uϵ)=(vϵ1,uϵ1).
It is easy to see that are solutions of (5.1) with λϵ=0. (Here we use x=(x1,x2)∈R2, y=x/ϵ=(y1,y2)∈R2.) Furthermore, since (vϵ,uϵ) are radially symmetric functions, we have that Φ1⊥Φ2 in L2(R2)⊕L2(R2). Here, we equip L2(R2)⊕L2(R2) with the
Some discussions
There are several interesting problems left.
- 1.
When ϵ2≪α≪O(1) andfrom Proof of Theorem 2.2, we can actually prove that all unstable eigenvalues for (vϵ1,uϵ1) must satisfy that (see (6.24)). That is, there are only small eigenvalues. (For large eigenvalues, by , , . So they are all stable by Theorem 5.1(1).) We can also decompose (Φϵ,Ψϵ) as in , . However, we cannot show that (Φϵ⊥,Ψϵ⊥)=(0,0) (one can see this from (6.24)). This means that the unstable modes, if
Acknowledgements
This research is supported by an earmarked research grant from RGC of Hong Kong. The author would like to thank Prof. S.-I. Ei, Prof. Y. Nishiura, and Prof. M.J. Ward for useful discussions. I thank the referees for their valuable remarks and suggestions which greatly improved the presentation of the paper.
References (50)
- et al.
Stability of singular patterns in the 1D Gray–Scott model: a matched asymptotic approach
Physica D
(1998) - et al.
Autocatalytic reactions in the isothermal, continuous stirred tank reactor: isolas and other forms of multistability
Chem. Eng. Sci.
(1983) - et al.
Autocatalytic reactions in the isothermal, continuous stirred tank reactor: oscillations and instabilities in the system A+2B→3B,B→C
Chem. Eng. Sci.
(1984) - et al.
Multiple interior peak solutions for some singularly perturbed Neumann problems
J. Diff. Eqs.
(1999) - et al.
Multiple boundary peak solutions for some singularly perturbed Neumann problems
Ann. Inst. H. Poincaré Anal. Nonlinéaire
(2000) - et al.
Stability and instability of Gray–Scott model: the case of equal diffusivities
Appl. Math. Lett.
(1999) - et al.
A skeleton structure of self-replicating dynamics
Physica D
(1999) On the construction of single-peaked solutions to a singularly perturbed semilinear Dirichlet problem
J. Diff. Eqs.
(1996)On the boundary spike layer solutions of singularly perturbed semilinear Neumann problem
J. Diff. Eqs.
(1997)- et al.
Stationary solutions for the Cahn–Hilliard equation
Ann. Inst. H. Poincaré Anal. Nonlinéaire
(1998)