A numerical approach for a general class of the spatial segregation of reactiondiffusion systems arising in population dynamics
Abstract.
In the current work we consider the numerical solutions of equations of stationary states for a general class of the spatial segregation of reactiondiffusion systems with population densities. We introduce a discrete multiphase minimization problem related to the segregation problem, which allows to prove the existence and uniqueness of the corresponding finite difference scheme. Based on that scheme, we suggest an iterative algorithm and show its consistency and stability. For the special case we show that the problem gives rise to the generalized version of the socalled twophase obstacle problem. In this particular case we introduce the notion of viscosity solutions and prove convergence of the difference scheme to the unique viscosity solution. At the end of the paper we present computational tests, for different internal dynamics, and discuss numerical results.
Key words and phrases:
Free boundary, Twophase obstacle problem, Reactiondiffusion systems, Finite difference, Viscosity solution2000 Mathematics Subject Classification:
35R35, 65N06, 65N22, 92D251. Introduction and known results
1.1. The setting of the problem
In recent years there have been intense studies of spatial segregation for reactiondiffusion systems. The existence of spatially inhomogeneous solutions for competition models of LotkaVolterra type in the case of two and more competing densities have been considered in [11, 12, 13, 15, 21, 20, 25]. The aim of this paper is to study the numerical solutions for a certain class of the spatial segregation of reactiondiffusion system with population densities.
Let be a connected and bounded domain with smooth boundary and be a fixed integer. We consider the steadystates of competing species coexisting in the same area . Let denotes the population density of the component with the internal dynamic prescribed by .
We call the tuple a segregated state if
The problem amounts to
(1) 
over the set
where for and on the boundary .
We assume that
where is Lipschitz continuous in uniformly continuous in and .
Remark 1.
Functions ’s are defined only for non negative values of s (recall that our densities ’s are assumed non negative); thus we can arbitrarily define such functions on the negative semiaxis. For the sake of convenience, when , we will let . This extension preserves the continuity due to the conditions on defined above. In the same way, each is extended as an even function.
Remark 2.
We emphasize that for the case the assumption is that for all the functions are nonnegative and uniformly continuous in . Also for simplicity, throughout the paper we shall call both and as internal dynamics.
Remark 3.
We would like to point out that the only difference between our minimization problem (1) and the problem discussed by Conti,Terrachini and Verzini [12], is the sign in front of the internal dynamics . In our case, the plus sign of allows to get rid of some additional conditions, which are imposed in [12, Section ]. Those conditions are important to provide coercivity of a minimizing functional in [12]. But in our case the above given conditions together with convexity assumption on with respect to the variable are enough to conclude which in turn implies coercivity of a functional (1).
In order to speak on the local properties of the population densities, let us introduce the notion of multiplicity of a point in .
Definition 1.
The multiplicity of the point is defined by:
where .
For the local properties of the same results as in [12] with the opposite sign in front of the internal dynamics hold. Below, for the sake of clarity, we write down those results from [12] with appropriate changes.
Lemma 1.
(Proposition in [12])
Assume that then the following holds:

If then there exists such that for every ;

If then there are and such that in

If then there are and such that for every and , we have and in
Next, we state the following uniqueness Theorem due to Conti, Terrachini and Verzini.
Theorem 1 (Theorem in [12]).
This theorem will play a crucial role in studying the difference scheme, especially for the case where we will reformulate it as a generalized twophase obstacle problem. Note that in this case, the problem will be reduced to:
(2) 
over the set
Here with property on the boundary .
1.2. Known results
In last years there has been much interest given to study the numerical approximation of reactiondiffusion type equations. For instance, the equations arising in the study of population ecology when high competitive interactions between different species occurs.
We refer the reader to [12, 16, 17, 18, 19, 20, 21] for models involving Dirichlet boundary data. A complete analysis of the stationary case has been studied in [12]. Also numerical simulation for the spatial segregation limit of two diffusive LotkaVolterra models in presence of strong competition and inhomogeneous Dirichlet boundary conditions is provided in [26].
In the work [6] Bozorgnia proposed two numerical algorithms for the problem (1) with the internal dynamics . The finite element approximation is based on the local properties of the solution. In this case the author was able to provide the convergence of the method. The second approach is a finite difference method, but lack of its analysis in [6]. This finite difference method has been generalized in [9] for the case of nonnegative . In [9] the authors present a numerical consistent variational system with strong interaction, and provide disjointness condition of populations during the iteration of the scheme. In this case the proposed algorithm is lack of deep analysis, especially for the case of three and more competing populations.
The present work concerns to close that gap and provides theoretical results for finite difference scheme, with competing populations and general internal dynamics satisfying certain suitable conditions. We introduce the discrete analogue of minimization problem and prove the existence and uniqueness of the difference scheme. Moreover for the special case we introduce viscosity solution and prove the convergence of corresponding difference scheme.
1.3. Notations
We will make the notations for the onedimensional and twodimensional cases parallely. For the sake of simplicity, we will assume that in onedimensional case and in twodimensional case in the rest of the paper, keeping in mind that the method works also for more complicated domains.
Let be a positive integer, and
We use the notation and (or simply , where is one or twodimensional index) for the finite difference scheme approximation to and , respectively. Concerning the boundary functions we assume they are extended to be zero everywhere outside the boundary for all . The discrete approximation for these functions will be and respectively (or simply , where is one or twodimensional index). Note that for the case with two population densities, we additionally will use the following notations:
and
in one and twodimensional cases, respectively.
In this paper we will use also notations , (not to be confused with functions ).
Denote
in one and two dimensional cases, respectively, and
In onedimensional case, we consider the following approximation for Laplace operator: for any ,
and for twodimensional case we introduce the following 5point stencil approximation for Laplacian:
for any .
2. Segregation problem with population densities
2.1. The minimizing functional and existence of difference scheme
In this section, we introduce the discrete counterpart of the spatial segregation problem for the general case. As in the previous section we assume that are convex in the variable for all and satisfy the properties stated in introduction.
In the rest of the paper the following notation
for elements , will play a crucial role. We focus on the following functional
(3) 
defined over the set
(4) 
where
Here for and , , the inner product is defined by
Theorem 2.
Proof.
We are going to prove that functional (3) is coercive in the set . Due to the standard arguments of calculus of variations coercivity and lower semicontinuity will imply the existence of minimizers over the closed set . To this end, we observe that
For every fixed and such that we have . Thus,
due to the simple fact that implies , which in turn yields . Therefore
for some constant . Recalling that we finally obtain that the functional (3) is coercive. ∎
Proposition 1.
If an element solves the following minimization problem:
then for every and
(6) 
Proof.
First of all, it is easy to verify that, if two vectors and belong to the set then for arbitrary we have
and
Let a vector be a minimizer to our discrete minimization problem (5). If we obtain:
for every . Thus,
We choose the vector such that for every and the following condition holds
This implies the following identity:
Hence,
(7) 
(8) 
Next, dividing both sides in (8) by and letting we arrive at:
(9) 
Let . We set as follows:
(10) 
It is easy to see, that the chosen vector satisfies . Therefore, the inequality (9) holds for this vector, which means it can be substituted into (9). We clearly obtain
and
Hence,
(11) 
In the same way, for every and we apparently have:
In this case we choose the vector such that for every and the following condition holds
This implies the following identity:
After proceeding the same steps as above, we obtain
(12) 
Lemma 2.
Let an element solves the minimization problem (5), and moreover assume that for some and . Then
Proof.
Assume that for some fixed and . There are two possibilities: either or . Let then we take a vector as follows:
Next theorem shows the existence of the corresponding finite difference scheme for two or more population densities. Note that for the particular case we obtain the difference scheme for the socalled multiphase obstacle problem.
Theorem 3.
If an element is a minimizer to (5), then the vector where solves the discrete system:
(16) 
for every and . Here for a given uniform mesh on we define
as the average of for all neighbour points of .
Proof.
Remark 4.
We would like to emphasize that the theorem statement holds also for more general uniformly discretized schemes without restriction on the number of stencil (neighbour) points corresponding to a discrete Laplacian.
2.2. Uniqueness of difference scheme
In this section our goal is to show the uniqueness of the difference scheme, which solves the discrete system (16). To this aim, we need two auxiliary lemmas.
For the sake of convenience we denote by the set of all closest neighbour points corresponding to a mesh point .
Lemma 3.
Let the functions be nondecreasing with respect to the variable . We take any two elements and in . If and are satisfying the discrete system (16), then the following equation holds:
for all .
Proof.
We argue by contradiction. Suppose for some we have
(17) 
Then taking into account the following simple chain of inclusions
(18) 
we obviously see that implies . On the other hand, the discrete system (16) gives us
Therefore
Thus,
which implies that