Spatial Segregation Limit of a Quasilinear Competition-Diffusion System ()
1. Introduction
In this paper, we study the spatial and temporal behavior of interacting biological species. Assuming the reaction rates of competition follow the Holling-Tanner interaction mechanism, the quasilinear reaction-diffusion model under consideration can be given by
(1)
here, , , where is a bounded domain in. are all positive constants. and stand for their population densities of the competing species at the time t and at the habitat. is the respective intrinsic growth rates, and represent the intra-specific competition rates, whereas and represent the inter-spe- cific competition rates. The boundary condition models the fact that species have no-flux near the boundary, where is the outward normal unit vector to. may not be equal to from an ecological point of view, but for the convenience of presentation, we may assume here.
Quasilinear parabolic equations have received a great attention in recent years. We can refer to [1] -[6] and the references therein for more details. However, the main concerns in above works are for the existence of a global solution, a weak solution, periodic solutions, the existence-uniqueness of positive solutions, blow-up property of the solution, and the qualitative property of the solution including finite time extinction and large time behavior of the solution.
Our main interest is different from those of the above works, we mainly consider the spatial segregation limit of (1) when only the interspecific competition rates and are very large. To study this case, it is convenient to rewrite (1) as the following equivalent form:
(2)
where and k are positive constants derived from, and k is the only parameter which is large. For similar studies, here we refer [7] -[15] to the interested readers for more information. A striking difference between (2) and above relevant works is that the diffusion term in (2) is quasilinear. When and, the system (2) is reduced to the classical Volterra-Lotka competition model, which has been studied in [9] , where Dancer et al. showed that the two competition species spatially segregate as k tends to infinity. Moreover, they proved that, for any, there exist subsequences and of the k-dependent non- negative solutions converging weakly in to the positive and negative parts respectively of a limit function w satisfying a scalar equation of the form
(P)
where, , and they also showed that the limit problem (P) turns out to be an explicit Stefan-like type free boundary problem.
Motivated by [9] , our main purpose of this paper is to extend most of results of [9] to systems (2) with quasilinear diffusion terms. In addition, we will get the convergence results for the further improvement. Specifically, we have strong convergence in.
Note that the study of strong-competition limits in corresponding elliptic of parabolic systems is of interest not only for questions of spatial segregation and coexistence in population dynamics, as here and in [7] [9] [13] [16] -[19] but also is key to the understanding of phase separation in Hartree-Fock type approximations of systems of modelling Bose-Einstein condensates, see [10] [20] [21] [23] , and reference therein.
To conclude, we observe that a couple of problems addressed and solved for family of solutions to (2) remains for further study in our general context: firstly, to develop a common regularity theory for the solutions of the system, which is independent of the competition rate, as in [16] -[18] [21] [22] ; secondly, to study the regularity of the class of limiting profiles, both in terms of the densities and in terms of the emerging free boundary problem, as in [10] [16] [23] [24] ; thirdly, the precise description of the singular set in the emerging free boundary problem, as in [25] [26] . These will be object of future investigation.
The outline of this paper is arranged as follows. In Section 2, we give some a prior estimates and some convergence results for solutions of problem (2). Section 3 is focused on the limit problem as. In Section 4, we get the further convergence results in the special case of. Concluding remarks are given in the last section.
2. Preliminaries
In order to study the limit case as, we rewrite problem (2) as
(3)
Throughout this paper, we let and suppose the initial functions and satisfy
(4)
We say a pair is a solution of (3) in the sense that and satisfy (3). We now prove some basic facts of solutions for problem (3), which will be used later.
Lemma 1. The solution of problem (3) exists and is unique. Moreover, there exist constants and such that
Proof. The existence and uniqueness of solutions of (3) are followed from the standard parabolic equations theory [4] .
By using the maximum principle, the solution is positive for and. For the upper bound, it follows from the comparison principle that for, where
which is the solution of the problem
(5)
Thus we have
Similarly, there exists a constant such that
Lemma 2. Let be the solution of problem (3), then
(6)
where is a constant which is independent of k.
Proof. Integrating the equation for in (3) over and using Green’s formula yield
(7)
By Lemma 1 and noting that the right side of (7) is independent of k, we get (6).
Lemma 3. Let be the solution of problem (3), then
(8)
where is a constant which is independent of k.
Proof. Multiplying the equation for in (3) by, integrating over and applying Green’s formula, we yield
which leads to
where we have used Lemma 1. To get the first estimate of (8), we simply integrate the above inequality from 0 to T. The second inequality of (8) can be derived similarly.
In order to derive a free boundary problem, we also need to introduce a new function
which is related with. Then satisfies the scalar problem
(9)
(10)
(11)
The following result yields uniform boundedness of.
Lemma 4. The sequence is bounded in uniformly with respect to k.
Proof. Multiplying the Equation (9) with, and integrating over using integration by parts, we get
where is the duality product between the space and. By Lemmas 1 and 3, we then have
where M is a positive constant which is independent of k or. This implies
With the above discussion, below we study some convergence properties. It follows from Lemmas 1 and 3 that and are uniformly bounded in. Hence, there exist subsequences of and (still denoted by and), and two functions such that
(12)
and
(13)
as. Furthermore, by Lemma 2, we have
Lemma 5. a.e. in.
Below we manage to build the relations between u, v and w.
Lemma 6. The subsequences and are such that
(14)
as, where and. Moreover,
(15)
Proof. Let be such that
In order to prove the theorem, we need to divide our proof into three cases:
In case, according to the definition of limit, there exists a positive constant such that
then we have
Due to Lemma 2, above inequality implies that
Next we consider case. We proceed as in the proof of case (a), then there exists a positive constant such that
Recalling Lemma 2, we claim that
For the last case. We claim that
Otherwise, if there is a subsequence of, which we still denote by, such that, it follows that, consequently, which contradicts the fact. Similarly, it is impossible to have that.
From the boundedness of and, it is easy to achieve convergence in. To the end, we get (15) from (14).
3. The Limit Problem as
Lemma 6 illustrates that and weakly in as. We set
(16)
and
(17)
In this section, we mainly consider the scalar equation
(18)
First, we show that problem (18) has a weak solution, which are defined as follows:
Definition 3.1 We say that a function is a weak solution of (16) if it satisfies
(19)
for all and any test function with.
Theorem 1. The function defined by (15) is the unique weak solution of problem (18). Moreover, for some.
Proof. From Lemmas 1 and 3, we easily have, and Lemma 4 yields. is derived from by a standard regularity result (see for example [27] , Theorem 3, p.287).
Multiplying (9) by a test function with, and using integration by parts, we deduce
Let along the sequence for which (12) holds. By the dominated convergence theorem and Lemma 1, we have
(20)
Note that (16) and (17) yield
and
With (20), we then have that z satisfies
for all and any test function with. Namely, z satisfies the differential equation in (18) as well as the homogeneous Neumann boundary condition in the sense of distributions, and the initial condition
This follows easily that z is the weak solution of problem (18).
It is clear from [2] that the weak solution of problem (18) is unique. Last, for the regularity of z, we refer to Theorems 1.1 and 1.3 in [28] .
According to the above discussion, there exists a family of closed hypersurfaces , which separates the two strongly competing species. That is
and
We denote
Finally, as in [9] , we rewrite a strong form of the limit problem (18), where the equations can be described a classical two-phase Stefan-like free boundary problem.
Theorem 2. Let z be a weak solution of limit problem (18), if is smooth enough, and if the functions
are smooth up to, then u and v satisfy
(21)
where we suppose that.
4. Further Convergence Results
In this section, we prove that the subsequences and of k-dependent non-negative solutions to (3) converge strongly in. For the convenience of presentation, we consider the special case of.
Theorem 3. If and a.e. in, then up to a subsequence,
(22)
and hence in, in, in.
Proof. By arguments as in the proof of Theorem 1, we first obtain
Thus
This implies
by Lemma 7.6 in [29] . Hence Also by Lemma 7.7 in [29] and Lemma 5, we get
(23)
Now, multiplying the second equation in (3) by the limit u and integrating it over, , we
have
Integrating by parts gives
(24)
Integrating (24) with respect to over gives
With (4), (12) and Lemma 5, as, we obtain
and
Since is bounded in and a.e. in, we may apply Fubini theorem to obtain
(25)
as. Similarly, by (12), (23) and Lemma 5, we have
Therefore, (24) yields
This implies that
(26)
Next, multiplying the first equation in (3) by the limit u and integration it over, we have
(27)
Integrating above equation in and passing to the limit as yield
(28)
by using (4), (12) and (26).
Finally, multiplying the first equation in (1) again by and integrating it over, we deduce
This concludes that
(29)
by (28). It follows from (12) and weak lower semi-continuity that
By Fatou’s lemma, we have
which together with (29) implies that there exists a subsequence, which we denote again by such that for a.e.
In other words,
Hence in. Similarly, we claim that in. The rest of the conclusions in this theorem follow consequently.
5. Concluding Remarks
The study of spatial behavior of the interacting species has been attracting much attention in population ecology, in particular, in the case when the interactions are large and of competitive type. Many different models based on partial differential equations can be successfully employed to investigate the phenomenon of coexistence and exclusions of competing species. In this paper, we have attempted to study a class of quasilinear parabolic system (3) describing a Holling-Tanner’s competitive interaction of two species. We prove that if inter-specific competition rates tend to infinity, then spatial segregation of the densities and a scalar limit problem (21) are given. In particular, we have obtained the strong convergence results in in the special case of. Ecologically, our results show that competition leads to segregation.
Finally, we want to mention that there are still many interesting questions to do for this kind of problem. First of all, noting that the diffusion term of the first equation in (2) can be written as, the term describes the “self-diffusion”. Naturally to ask whether our results can be extended to parabolic systems with “cross-diffusion”? Moreover, as mentioned in the introduction, we have seen that limit profiles of solutions to (2) are segregated configurations, it is then natural to define the free boundary as the nodal set. The regularity of the nodal set remains a challenge, and it will be the object of a forthcoming paper.
Acknowledgements
We thank the Editor and the referee for their comments. This work is partially supported by PRC grant NSFC 11501494 and NSF of the Higher Education Institutions of Jiangsu Province (12KJD110008). This support is greatly appreciated.
NOTES
*Corresponding author.