# Permanence of Diffusive Models for Three Competing Species in Heterogeneous Environments

Authors: Benlong Xu; Zhenzhang Ni
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

We address the question of the long-term coexistence of three competing species whose dynamics are governed by the partial differential equations. We obtain criteria for permanent coexistence in a Lotka-Volterra system modeling the interaction of three competing species in a bounded habitat whose exterior is lethal to each species. It is also proved that if the intercompeting strength is very weak, the system is always permanent, provided that each single one of the three species can survive in the absence of the two  other species.


---

## Body

## 1. Introduction

Reaction-diffusion systems are some of the most widely used models for population dynamics in situations where spatial dispersal plays a significant role. Generally speaking, the dynamical behaviors of the systems with two competing species are relatively simple and have been widely investigated in the past years. However, the investigation and knowledge about the systems with three competing species are very limited; see [1–5] for some related results. It is well known, in general, that the dynamical behavior of systems with three competing species may be extremely complex, even in nondiffusive case, namely, ODE case; see, for example, [6–9] and the references therein.The object of this paper is to study the problem of coexistence for three competing species dispersing through a spatially heterogeneous region. We model the population dynamics of the species with a system of three diffusive Lotka-Volterra equations.The model we consider has the general form(1)
∂
u
1
∂
t
=
μ
1
Δ
u
1
+
u
1
(
a
1
(
x
)
-
b
11
(
x
)
u
1
-
b
12
(
x
)
u
2
-
b
13
(
x
)
u
3
)
in
Ω
,
∂
u
2
∂
t
=
μ
2
Δ
u
2
+
u
2
(
a
2
(
x
)
-
b
21
(
x
)
u
1
-
b
22
(
x
)
u
2
-
b
23
(
x
)
u
3
)
in
Ω
,
∂
u
3
∂
t
=
μ
3
Δ
u
3
+
u
3
(
a
3
(
x
)
-
b
31
(
x
)
u
1
-
b
32
(
x
)
u
2
-
b
33
(
x
)
u
3
)
in
Ω
with u
i
=
0 on ∂
u
i
×
R
+, where Ω
⊆
R
N is a bounded domain and a
i
(
x
)
>
0, b
i
i
(
x
)
>
0, b
i
j
(
x
)
≥
0
(
i
,
j
=
1,2
,
3
,
i
≠
j
) are all continuous in Ω
¯. The variables u
i represent population densities of the competing species. The boundary condition u
i
=
0 describes a situation that the boundary of Ω is lethal to the species.We will use the criterion of permanence to characterize coexistence. A system is said to be permanent if any solutions with all components positive initially must ultimately enter and remain within a fixed set of positive states that are strictly bounded away from zero in each component. For some investigations of coexistence characterized by permanence, mostly in the case for two interacting species, see [10–12], and also some parts of the books [13, 14].This paper is organized as follows. In Section2, we introduce some necessary material about semiflow and state V. Hutson’s average Lyapunov theorem which will be the basic tool in the proof of our main theorem. We also discuss the dissipative property of our systems, which is necessary for using V. Hutson’s average Lyapunov function theorem. In Section 3, we establish our main permanent result, Theorem 11. We also prove that if the intercompeting strength is “very weak,” the system (1) is always permanent, provided that either single one of the three species can survive with the absence of the two other species. In Section 4, we give some summery and discussion of the results of this paper.


## 2. Semiflows and Dissipativity

In this section, we will make some preparations for establishing our main results in the next section. Firstly, we will introduce some terminologies and results about semiflows (semigroups) for readers’ convenience. The materials can be found in [15] or [13, 14].Let(
Y
,
d
) be a metric space, with points in Y being denoted by u
,
v
,
… and subsets Y by U
,
V
,
…. The following two unsymmetric distances of sets will be used:
(2)
d
¯
(
U
,
V
)
=
sup
⁡
u
∈
U
d
(
u
,
V
)
,
d
_
(
U
,
V
)
=
inf
⁡
u
∈
U
d
(
u
,
V
)
.
The triple (
Y
,
π
,
R
+
) is said to be a semiflow (or semigroup), if π
:
Y
×
R
+
→
Y is continuous and satisfies:(i)
π
(
u
,
0
)
=
u,
(ii)
π
(
π
(
u
,
t
)
,
s
)
=
π
(
u
,
t
+
s
)
(
s
,
t
∈
R
+
),for allu
∈
Y. For convenience, we often write π
(
u
,
t
)
=
u
·
t. The symbols γ
+
(
u
) and ω
(
u
) denote the semiorbit through u and the omega limit set of u, respectively, and the equivalent expressions for sets are defined by taking unions.A solutionϕ through u is a continuous map ϕ
:
R
→
Y such that ϕ
(
0
)
=
u and π
(
ϕ
(
τ
)
,
t
)
=
ϕ
(
t
+
τ
) for t
∈
R
+
,
τ
∈
R. The range of ϕ is denoted by γ
(
u
) and is called an orbit through u.A setU is said to be forward invariant if γ
+
(
u
)
⊂
U and invariant if γ
(
u
)
⊂
U. The semiflows is said to be dissipative if there is a bounded set U such that lim
⁡
t
→
∞
⁡
d
(
u
·
t
,
U
)
=
0 for all u
∈
Y. U is said to be a global attractor of the semiflow if it is compact invariant and for all bounded V, lim
⁡
t
→
∞
d
¯
(
V
·
t
,
U
)
=
0, where V
·
t
=
{
v
·
t
,
v
∈
V
}.Theorem 1.
LetY be complete and suppose that the semiflow isdissipative. Assume that there is a t
0
≥
0 such that π
(
·
,
t
) is compact for t
>
t
0; then there is a nonempty global attractor.Consider next the concept of permanence in the abstract semiflow context. We suppose thatY
=
Y
0
∪
∂
Y
0, where Y
0 is open, and assume that Y
0, ∂
Y
0 are forward invariant. In relation to the remarks in the introduction, ∂
Y
0 will consist of the states with at least one species identically zero.Definition 2.
The semiflow is said to be permanent if there exists a bounded setU with d
_
(
U
,
∂
Y
0
)
>
0 such that lim
⁡
t
→
∞
d
(
v
·
t
,
U
)
=
0 for all v
∈
Y
0.A setU
⊂
Y
0 is said to be strongly bounded if it is bounded and d
_
(
U
,
∂
Y
0
)
>
0. A
0 is said to be a global attractor relative to strongly bounded sets if it is a compact invariant subset of Y
0 and lim
⁡
t
→
∞
d
¯
(
U
·
t
,
A
0
)
=
0 for all strongly bounded U.Permanence is obviously an asymptotic property. It can thus be studied by examining the semiflow restricted to a forward invariant set derived from anɛ-neighbourhood B
(
A
,
ɛ
) of the global attractor A of Theorem 1. Set then X
=
c
l
π
(
B
(
A
,
ɛ
)
,
[
1
,
∞
)
), the closure of π
(
B
(
A
,
ɛ
)
,
[
1
,
∞
)
), and take S
=
X
∩
∂
Y
0. The following Hutson’s theorem on average Lyapunov functions is the basic tool for establishing our main theorem in Section 3.Theorem 3 (see [15]).
Assume that the conditions of Theorem1 hold, and let X, S be as defined above. Suppose that P
:
X
∖
S
→
R
+ is continuous, strictly positive, and bounded, and for u
∈
S define
(3)
α
(
t
,
u
)
=
liminf
⁡
υ
→
u
,
υ
∈
X
∖
S
(
P
(
υ
·
t
)
P
(
υ
)
)
.
Then the semiflow is permanent if
(4)
sup
⁡
t
>
0
α
(
t
,
u
)
>
1
u
∈
ω
(
S
)
,
sup
⁡
t
>
0
α
(
t
,
u
)
>
0
u
∈
S
.For our application of the last theorem, we will cast the system (1) to the abstract frame of the semiflow. It is well known that if (1) is simply viewed as a parabolic system, it actually generated a semiflow in C
1
(
Ω
¯
), and solutions which belong to C
1
(
Ω
¯
) must have Hölder-continuous second derivatives on Ω (see [13]). Thus, we may use maximum principles to obtain a priori bounds, even though we will ultimately want to view our semiflow as acting on C
1
(
Ω
¯
). Then by maximum principle (or comparison theorem), any solution (
u
1
,
u
2
,
u
3
) of (1) must satisfy a uniform bound of the form 0
≤
u
i
≤
C
i for some constants C
i, i
=
1,2
,
3, after finite time. So the semiflow generated by (1) is dissipative. We state it as follows.Theorem 4.
The semiflow generated by the system (1) is dissipative.


## 3. Permanent Coexistence Results

In this section, we establish criteria of permanent coexistence of the system (1). The main tool is Hutson’s average Lyapunov function theorem stated in Section 2, and the key step is to set up a suitable average Lyapunov function. To construct the average Lyapunov functions, we must have a detailed knowledge of the ω-limit set of the generated semiflow in the boundary of the positive cone.The maximum principle implies that solutions of (1) with nonnegative nonzero initial data for a given component will have that component strictly positive in Ω for all t
>
0. In the case of Dirichlet conditions, such solutions will have normal derivatives on ∂
Ω which are bounded above by a negative constant. Hence, the only trajectories which remain in the boundary of the positive cone have one or both components identically zero. Thus, to determine the ω-limit set of the semiflow generated by (1) on the boundary, we need only consider the steady state solutions (equilibrium points) of subsystems of (1).Letλ
1
>
0 denote the principal eigenvalue for problem
(5)
Δ
ϕ
+
λ
ϕ
=
0
in
Ω
,
ϕ
=
0
on
⁡
∂
Ω
.In the rest of this paper, we denote(6)
M
i
=
max
⁡
x
∈
Ω
¯
a
i
(
x
)
b
i
i
(
x
)
.The point(
0,0
,
0
) is always a steady state, which means that there are no species in the domain. With the absence of two species and only the other one species left, system (1) becomes a scalar equation, and we have the following well-known result (see [14]).Lemma 5.
Suppose thata
(
x
)
>
μ
i
λ
1 for x
∈
Ω
¯. Then the following problem
(7)
u
i
t
=
μ
i
Δ
u
i
+
u
i
(
a
i
(
x
)
-
b
i
i
(
x
)
u
i
)
i
n
Ω
×
(
0
,
∞
)
,
u
i
=
0
o
n
∂
Ω
×
(
0
,
∞
)
,
for =
1,2
,
3, has a unique positive steady state solution u
¯
i
(
x
) which is globally approximately stable.Remark 6.
Lemma5 tells us that if a
(
x
)
>
μ
i
λ
1 for x
∈
Ω
¯, i
=
1,2
,
3, then either single one of the three species can survive in the absence of the other two species.Now we consider the subsystem of (1) with only one species to be absent. There are three cases.Case 1.
One hasu
3
≡
0,
(8)
u
1
t
=
μ
1
Δ
u
1
+
u
1
(
a
1
(
x
)
-
b
11
(
x
)
u
1
-
b
12
(
x
)
u
2
)
hhhhhhhhhhhhhhhhhhhhhhhhhhh
in
Ω
×
(
0
,
∞
)
,
u
2
t
=
μ
2
Δ
u
2
+
u
2
(
a
2
(
x
)
-
b
21
(
x
)
u
1
-
b
22
(
x
)
u
2
)
hhhhhhhhhhhhhhhhhhhhhhhhhhh
in
Ω
×
(
0
,
∞
)
,
u
1
=
u
2
=
0
on
⁡
∂
Ω
×
(
0
,
∞
)
.Case 2.
One hasu
2
≡
0,
(9)
u
1
t
=
μ
1
Δ
u
1
+
u
1
(
a
1
(
x
)
-
b
11
(
x
)
u
1
-
b
13
(
x
)
u
3
)
hhhhhhhhhhhhhhhhhhhhhhhhhhh
in
Ω
×
(
0
,
∞
)
,
u
3
t
=
μ
3
Δ
u
3
+
u
3
(
a
3
(
x
)
-
b
31
(
x
)
u
1
-
b
33
(
x
)
u
3
)
hhhhhhhhhhhhhhhhhhhhhhhhhhh
in
Ω
×
(
0
,
∞
)
,
u
1
=
u
3
=
0
on
⁡
∂
Ω
×
(
0
,
∞
)
.Case 3.
One hasu
1
≡
0,
(10)
u
2
t
=
μ
2
Δ
u
2
+
u
2
(
a
2
(
x
)
-
b
22
(
x
)
u
2
-
b
23
(
x
)
u
3
)
hhhhhhhhhhhhhhhhhhhhhhhhhhh
in
Ω
×
(
0
,
∞
)
,
u
3
t
=
μ
3
Δ
u
3
+
u
3
(
a
3
(
x
)
-
b
32
(
x
)
u
2
-
b
33
(
x
)
u
3
)
hhhhhhhhhhhhhhhhhhhhhhhhhhh
in
Ω
×
(
0
,
∞
)
,
u
1
=
u
3
=
0
on
⁡
∂
Ω
×
(
0
,
∞
)
.
With some modification of lower-supper solution methods of parabolic systems (see [14]), it is easy to prove the following results.Theorem 7.
WithM
i
,
i
=
1,2
,
3, as defined in (6),(1)
if(11)
a
1
(
x
)
>
μ
1
λ
1
+
M
2
b
12
(
x
)
,
a
2
(
x
)
>
μ
2
λ
1
+
M
1
b
21
(
x
)
,
∀
x
∈
Ω
¯
,
the system (8) has a steady state solution (
u
~
(
12
)
(
x
)
,
u
~
(
21
)
(
x
)
) with u
~
(
12
)
(
x
)
>
0
,
u
~
(
21
)
(
x
)
>
0 for x
∈
Ω; similarly,
(2)
if(12)
a
1
(
x
)
>
μ
1
λ
1
+
M
3
b
13
(
x
)
,
a
3
(
x
)
>
μ
3
λ
1
+
M
1
b
31
(
x
)
,
∀
x
∈
Ω
¯
,
the system (9) has a steady state solution(
u
~
(
13
)
(
x
)
,
u
~
(
31
)
(
x
)
) with u
~
(
13
)
(
x
)
>
0, u
~
(
31
)
(
x
)
>
0 for x
∈
Ω;
(3)
if(13)
a
2
(
x
)
>
μ
2
λ
1
+
M
3
b
23
(
x
)
,
a
3
(
x
)
>
μ
3
λ
1
+
M
2
b
32
(
x
)
,
∀
x
∈
Ω
¯
,
the system (10) has a steady state solution (
u
~
(
23
)
(
x
)
,
u
~
(
32
)
(
x
)
) with u
~
(
23
)
(
x
)
>
0, u
~
(
32
)
(
x
)
>
0 for x
∈
Ω.Proof.
We only prove the part (1); the proof of part (2) and part (3) are similar.
Letw
(
x
) be a positive eigenfunction of the principal eigenvalue λ
1 for the eigenvalue problem (5). Choose r
1
>
0 sufficiently small,
(14)
μ
1
Δ
(
r
1
w
(
x
)
)
+
r
1
w
(
x
)
×
[
a
1
(
x
)
-
b
11
(
x
)
r
1
w
(
x
)
-
b
12
(
x
)
u
2
]
≥
r
1
w
(
x
)
[
a
1
(
x
)
-
μ
1
λ
1
hhhhhhhhhhh
-
b
11
(
x
)
r
1
w
(
x
)
-
b
12
(
x
)
M
2
]
≥
0
,
μ
1
Δ
(
M
1
)
+
M
1
[
(
a
1
(
x
)
-
b
11
(
x
)
M
1
-
b
12
(
x
)
u
2
]
=
M
1
[
(
a
1
(
x
)
-
b
11
(
x
)
M
1
-
b
12
(
x
)
u
2
]
≤
M
1
[
(
a
1
(
x
)
-
b
11
(
x
)
M
1
]
≤
0
for all 0
≤
u
2
≤
M
2. So, (
M
1
,
r
1
w
(
x
)
) is a set of upper and lower solutions for u
1 in (8).
Similarly, chooser
2
>
0 sufficiently small; (
M
2
,
r
2
w
(
x
)
) is a set of upper and lower solutions for u
2 in (8).
By coupled upper and lower theorem (see [16, Theorem 1.4-2]), the system (8) has a steady state solution (
u
~
(
12
)
(
x
)
,
u
~
(
21
)
(
x
)
) with u
~
(
12
)
(
x
)
>
0, u
~
(
21
)
(
x
)
>
0 for x
∈
Ω.Remark 8.
By the comparison principle, it is easy to see that, inΩ, we have u
~
(
12
), u
~
(
13
)
<
u
¯
1
<
M
1, u
~
(
21
), u
~
(
23
)
<
u
¯
2
<
M
2, and u
~
(
31
), u
~
(
32
)
<
u
¯
3
<
M
3.Now, we consider the unique problem of the steady state solutions above.Letθ
μ
,
A
(
x
) denote the unique positive solution of
(15)
μ
i
Δ
θ
+
θ
(
A
(
x
)
-
θ
)
=
0
in
Ω
,
θ
=
0
on
⁡
∂
Ω
,
for any A
(
x
)
>
μ
i
λ
1. For brevity, we denote
(16)
u
¯
(
1
)
=
θ
μ
1
,
a
1
(
x
)
,
u
¯
(
2
)
=
θ
μ
2
,
a
2
(
x
)
,
u
¯
(
3
)
=
θ
μ
3
,
a
3
(
x
)
,
u
_
(
12
)
=
θ
μ
1
,
(
a
1
(
x
)
-
M
2
b
12
(
x
)
)
,
u
_
(
21
)
=
θ
μ
2
,
(
a
2
(
x
)
-
M
1
b
21
(
x
)
)
,
u
_
(
13
)
=
θ
μ
1
,
(
a
1
(
x
)
-
M
3
b
13
(
x
)
)
,
u
_
(
31
)
=
θ
μ
3
,
(
a
3
(
x
)
-
M
1
b
31
(
x
)
)
,
u
_
(
23
)
=
θ
μ
2
,
(
a
2
(
x
)
-
M
3
b
23
(
x
)
)
,
u
_
(
32
)
=
θ
μ
3
,
(
a
3
(
x
)
-
M
2
b
32
(
x
)
)
.
They are all positive functions in Ω.Theorem 9.
Assume that all the hypotheses of Theorem7 are satisfied.(1)
If(17)
b
12
2
(
x
)
u
¯
(
1
)
b
11
2
(
x
)
u
_
(
21
)
+
2
b
12
(
x
)
b
21
(
x
)
b
11
(
x
)
b
22
(
x
)
+
b
21
2
(
x
)
u
¯
(
2
)
b
22
2
(
x
)
u
_
(
12
)
<
4
,
then the strictly positive steady state solution(
u
~
(
12
)
(
x
)
,
u
~
(
21
)
(
x
)
) of system (8) is unique; similarly
(2)
if(18)
b
13
2
(
x
)
u
¯
(
1
)
b
11
2
(
x
)
u
_
(
31
)
+
2
b
13
(
x
)
b
31
(
x
)
b
11
(
x
)
b
33
(
x
)
+
b
31
2
(
x
)
u
¯
(
3
)
b
33
2
(
x
)
u
_
(
13
)
<
4
,
then the steady state solution(
u
~
(
13
)
(
x
)
,
u
~
(
31
)
(
x
)
) of system (9) is unique;
(3)
if(19)
b
23
2
(
x
)
u
¯
(
2
)
b
22
2
(
x
)
u
_
(
32
)
+
2
b
23
(
x
)
b
32
(
x
)
b
22
(
x
)
b
33
(
x
)
+
b
32
2
(
x
)
u
¯
(
3
)
b
33
2
(
x
)
u
_
(
23
)
<
4
,
then the steady state solution(
u
~
(
23
)
(
x
)
,
u
~
(
32
)
(
x
)
) of system (10) is unique.Remark 10.
For fixed functionsa
i
(
x
)
>
0, b
i
i
(
x
)
>
0, (
i
=
1,2
,
3
), hypothesis (17) will be satisfied for b
12
(
x
), b
21
(
x
)
≥
0 sufficiently small. This is true because of Hopf’s strong maximal value theorem, and also because u
_
(
12
)(resp., u
_
(
21
)) increases as b
12
(
x
) (resp., b
21
(
x
)) decreases for x
∈
Ω. Thus (
b
21
2
(
x
)
u
¯
(
2
)
/
b
22
2
(
x
)
u
_
(
12
)
) (resp., b
12
2
(
x
)
u
¯
(
1
)
/
b
11
2
(
x
)
u
_
(
21
)) decreases as b
12
(
x
) (resp., b
21
(
x
)) decreases. Similarly, hypotheses (18) and (19) will be satisfied if b
13
(
x
), b
31
(
x
), b
23
(
x
), and b
32
(
x
) are sufficiently small.Proof.
We only give the proof of part (1) in Theorem9, since the arguments of part (2) and part (3) are similar.
Assume that(
u
~
(
12
)
(
x
)
,
u
~
(
21
)
(
x
)
), (
v
~
(
12
)
(
x
)
,
v
~
(
21
)
(
x
)
) are two strictly positive steady state solutions of system (8) in Ω.
Let(20)
p
(
x
)
=
u
~
(
12
)
(
x
)
-
v
~
(
12
)
(
x
)
,
q
(
x
)
=
u
~
(
21
)
(
x
)
-
v
~
(
21
)
(
x
)
I
12
=
b
11
(
x
)
v
~
(
12
)
(
x
)
p
(
x
)
+
b
12
(
x
)
v
~
(
12
)
(
x
)
q
(
x
)
I
21
=
b
21
(
x
)
u
~
(
21
)
(
x
)
p
(
x
)
+
b
22
(
x
)
u
~
(
21
)
(
x
)
q
(
x
)
;
then
(21)
μ
1
Δ
p
(
x
)
+
[
a
1
(
x
)
-
b
11
(
x
)
u
~
(
12
)
(
x
)
h
h
h
h
h
h
h
h
h
h
-
b
12
(
x
)
u
~
(
21
)
(
x
)
]
p
(
x
)
-
I
12
=
0
hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh
in
Ω
,
μ
2
Δ
q
(
x
)
+
[
a
2
(
x
)
-
b
21
(
x
)
v
~
(
12
)
(
x
)
h
h
h
h
h
h
h
h
h
h
-
b
22
(
x
)
v
~
(
21
)
(
x
)
]
q
(
x
)
-
I
21
=
0
hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh
in
Ω
,
p
(
x
)
=
q
(
x
)
=
0
on
⁡
∂
Ω
.
Sinceu
~
(
12
)
(
x
) is a strictly positive solution of(22)
μ
1
Δ
ψ
+
[
a
1
(
x
)
-
b
11
(
x
)
u
~
(
12
)
(
x
)
hhhhhhhh
-
b
12
(
x
)
u
~
(
21
)
(
x
)
]
ψ
+
α
ψ
=
0
hhhhhhhhhhhhhhhhhhhhhhhhhhh
in
Ω
,
ψ
=
0
on
⁡
∂
Ω
,
with α
=
0, the number α
=
0 must be the smallest eigenvalue of the above problem. Moreover, by variational properties, we have
(23)
∫
Ω
z
(
-
μ
1
Δ
z
-
[
a
1
(
x
)
-
b
11
(
x
)
u
~
(
12
)
(
x
)
hhhhhhhhhhhh
-
b
12
(
x
)
u
~
(
21
)
(
x
)
]
z
)
d
x
≥
0
,
for any z
∈
C
2
(
Ω
¯
) which vanishes on ∂
Ω. Similarly, since v
~
(
21
)
(
x
) is strictly positive solution of
(24)
μ
2
Δ
ψ
+
[
a
2
(
x
)
-
b
21
(
x
)
v
~
(
12
)
(
x
)
hhhhhhh
-
b
22
(
x
)
v
~
(
21
)
(
x
)
]
ψ
+
α
ψ
=
0
hhhhhhhhhhhhhhhhhhhhhhhhhh
in
Ω
,
ψ
=
0
on
⁡
∂
Ω
,
with α
=
0, the number α
=
0 must be the smallest eigenvalue of the above problem. Moreover,
(25)
∫
Ω
z
(
-
μ
2
Δ
z
-
[
a
2
(
x
)
-
b
21
(
x
)
v
~
(
12
)
(
x
)
-
b
22
(
x
)
v
~
(
21
)
(
x
)
]
z
)
d
x
≥
0
,
for any z
∈
C
2
(
Ω
¯
) which vanishes on ∂
Ω. Multiplying the first equation of (21) by -
p
(
x
), the second by -
q
(
x
), integrating over Ω, and adding, we deduce from (23) and (25) that
(26)
∫
Ω
[
b
11
(
x
)
v
~
(
12
)
(
x
)
p
2
(
x
)
+
(
b
12
(
x
)
v
~
(
12
)
(
x
)
+
b
21
(
x
)
u
~
(
21
)
(
x
)
)
p
(
x
)
q
(
x
)
+
b
22
(
x
)
u
~
(
21
)
(
x
)
q
2
(
x
)
]
d
x
≤
0
.
By comparison of scalar equations using upper and lower solutions we can readily obtain for i
=
1,2
,
x
∈
Ω,
(27)
(
1
b
11
(
x
)
)
u
_
(
12
)
≤
u
~
(
12
)
,
v
~
(
12
)
≤
(
1
b
11
(
x
)
)
u
¯
(
1
)
,
(
1
b
22
(
x
)
)
u
_
(
21
)
≤
u
~
(
21
)
,
v
~
(
21
)
≤
(
1
b
22
(
x
)
)
u
¯
(
2
)
.
From (27), we have
(28)
b
22
(
x
)
b
12
2
(
x
)
u
¯
(
1
)
b
11
(
x
)
u
_
(
21
)
+
2
b
12
(
x
)
b
21
(
x
)
+
b
11
(
x
)
b
21
2
(
x
)
u
¯
(
2
)
b
22
(
x
)
u
_
(
12
)
>
b
12
2
(
x
)
v
~
(
12
)
(
x
)
u
~
(
21
)
(
x
)
+
2
b
12
(
x
)
b
21
(
x
)
+
b
21
2
(
x
)
u
~
(
21
)
(
x
)
v
~
(
12
)
(
x
)
in Ω. It follows from (17) that
(29)
b
12
2
(
x
)
v
~
(
12
)
(
x
)
u
~
(
21
)
(
x
)
+
2
b
12
(
x
)
b
21
(
x
)
+
b
21
2
(
x
)
u
~
(
21
)
(
x
)
v
~
(
12
)
(
x
)
<
4
b
11
(
x
)
b
22
(
x
)
,
in
Ω
.
Then it is easy to see that the quadratic expression in the integrand of (26) is positive definite for each x
∈
Ω. Consequently, we must have p
(
x
) and q
(
x
) identically equal to zero in Ω. That is (
u
~
(
12
)
(
x
)
,
u
~
(
21
)
(
x
)
)
≡
(
v
~
(
12
)
(
x
)
,
v
~
(
21
)
(
x
)
) in Ω.Now we are ready to develop the permanent coexistence results of the system (1). To construct the average Lyapunov functions, we will use the positive eigenfunctions ψ
1, ψ
2, and ψ
3 corresponding to the principal eigenvalues σ
1, σ
2, and σ
3 of
(30)
μ
1
Δ
ψ
1
+
ψ
1
(
a
1
(
x
)
-
b
12
(
x
)
u
~
(
23
)
-
b
13
(
x
)
u
~
(
32
)
)
=
σ
ψ
1
in
Ω
,
ψ
1
=
0
on
⁡
∂
Ω
,
(31)
μ
2
Δ
ψ
2
+
ψ
2
(
a
2
(
x
)
-
b
21
(
x
)
u
~
(
13
)
-
b
23
(
x
)
u
~
(
31
)
)
=
σ
ψ
2
in
Ω
,
ψ
2
=
0
on
⁡
∂
Ω
,
(32)
μ
3
Δ
ψ
3
+
ψ
3
(
a
3
(
x
)
-
b
31
(
x
)
u
~
(
12
)
-
b
32
(
x
)
u
~
(
21
)
)
=
σ
ψ
3
in
Ω
,
ψ
3
=
0
on
⁡
∂
Ω
,
respectively.Theorem 11.
Suppose that (11)–(13), (17)–(19) hold, and the principal eigenvalues σ
1, σ
2, and σ
3 of (30), (31), and (32), respectively, are all positive. Assume also that
(33)
σ
1
+
b
12
(
x
)
u
~
(
23
)
+
b
13
(
x
)
u
~
(
32
)
>
max
⁡
{
b
12
(
x
)
u
¯
2
,
b
13
(
x
)
u
¯
3
}
,
σ
2
+
b
21
(
x
)
u
~
(
13
)
+
b
23
(
x
)
u
~
(
31
)
>
max
⁡
{
b
21
(
x
)
u
¯
1
,
b
23
(
x
)
u
¯
3
}
,
σ
3
+
b
31
(
x
)
u
~
(
12
)
+
b
32
(
x
)
u
~
(
21
)
>
max
⁡
{
b
31
(
x
)
u
¯
1
,
b
32
(
x
)
u
¯
2
}
,
for all x
∈
Ω
¯; then the semiflow on [
C
0
1
(
Ω
¯
)
]
3 generated by (1) under homogeneous Dirichlet boundary conditions is permanent.Proof.
The hypotheses imply that (1) is dissipative and that the ω-limit set in the boundary of the positive cone consists of the following seven cases: (
0,0
,
0
), (
u
¯
1
,
0,0
), (
0
,
u
¯
2
,
0
), (
0,0
,
u
¯
3
), (
u
~
(
12
)
,
u
~
(
21
)
,
0
), (
u
~
(
13
)
,
0
,
u
~
(
31
)
), and (
0
,
u
~
(
23
)
,
u
~
(
32
)
). Thus we need only find a suitable average Lyapunov function to conclude permanence. Since the semiflow is dissipative, we may restrict our attention to a bounded absorbing subset X of the positive cone in [
C
0
1
(
Ω
¯
)
]
3 that contains the global attractor whose existence is asserted in Theorem 1. Let S denote the intersection of X with the boundary of the positive cone. Choose eigenfunctions ψ
1, ψ
2, ψ
3
>
0 for (30), (31), and (32), respectively, and define
(34)
P
(
(
v
1
,
v
2
,
v
3
)
)
=
(
∫
Ω
‍
ψ
1
v
1
d
x
)
(
∫
Ω
‍
ψ
2
v
2
d
x
)
(
∫
Ω
‍
ψ
3
v
3
d
x
)
=
exp
⁡
[
log
⁡
∫
Ω
‍
ψ
1
v
1
d
x
+
log
⁡
∫
Ω
‍
ψ
2
v
2
d
x
h
h
h
h
h
h
+
log
⁡
∫
Ω
‍
ψ
3
v
3
d
x
]
.
We have for (
u
1
,
u
2
,
u
3
)
∈
S that
(35)
α
(
t
,
(
u
1
,
u
2
,
u
3
)
)
=
liminf
⁡
(
v
1
,
v
2
,
v
3
)
∈
X
∖
S
(
v
1
,
v
2
,
v
3
)
→
(
u
1
,
u
2
,
u
3
)
P
(
(
v
1
,
v
2
,
v
3
)
·
t
)
P
(
(
v
1
,
v
2
,
v
3
)
)
we need sup
⁡
t
>
0
⁡
α
(
t
,
(
u
1
,
u
2
,
u
3
)
)
>
0 for (
u
1
,
u
2
,
u
3
)
∈
S and sup
⁡
t
>
0
⁡
α
(
t
,
(
u
1
,
u
2
,
u
3
)
)
>
1 for (
u
1
,
u
2
,
u
3
)
∈
ω
(
S
). Let (
w
1
(
t
)
,
w
2
(
t
)
,
w
3
(
t
)
)
=
(
(
v
1
,
v
2
,
v
3
)
·
t
). Computation using (34) yields
(36)
P
(
(
v
1
,
v
2
,
v
3
)
·
t
)
P
(
(
v
1
,
v
2
,
v
3
)
)
=
exp
⁡
[
(
log
⁡
∫
Ω
‍
ψ
1
w
1
d
x
|
t
-
log
⁡
∫
Ω
‍
ψ
1
w
1
d
x
|
0
)
+
(
log
⁡
∫
Ω
‍
ψ
2
w
2
d
x
|
t
-
log
⁡
∫
Ω
‍
ψ
2
w
2
d
x
|
0
)
+
(
log
⁡
∫
Ω
‍
ψ
3
w
3
d
x
|
t
-
log
⁡
∫
Ω
‍
ψ
3
w
3
d
x
|
0
)
]
=
exp
⁡
[
∫
0
t
‍
(
∫
Ω
‍
ψ
1
w
1
t
d
x
∫
Ω
‍
ψ
1
w
1
d
x
)
+
∫
0
t
‍
(
∫
Ω
‍
ψ
2
w
2
t
d
x
∫
Ω
‍
ψ
2
w
2
d
x
)
+
∫
0
t
‍
(
∫
Ω
‍
ψ
3
w
3
t
d
x
∫
Ω
‍
ψ
3
w
3
d
x
)
]
.
We have fort
>
0
(37)
∫
Ω
‍
ψ
1
w
1
t
d
x
=
∫
Ω
‍
ψ
1
[
μ
1
Δ
w
1
+
w
1
(
a
1
(
x
)
-
b
11
(
x
)
w
1
hhhhhhhhhhhhhhhhhhh
-
b
12
(
x
)
w
2
-
b
13
(
x
)
w
3
)
]
d
x
=
∫
Ω
‍
[
(
μ
1
Δ
ψ
1
)
w
1
+
(
a
1
(
x
)
-
b
11
(
x
)
w
1
hhhhhhhhhhhhhhhhhh
-
b
12
(
x
)
w
2
-
b
13
(
x
)
w
3
)
ψ
1
w
1
]
d
x
=
∫
Ω
‍
[
σ
1
-
(
a
1
(
x
)
-
b
12
(
x
)
u
~
(
23
)
-
b
13
(
x
)
u
~
(
32
)
)
hhhhhh
+
(
a
1
(
x
)
-
b
11
(
x
)
w
1
-
b
12
(
x
)
w
2
hhhhhhhhh
-
b
13
(
x
)
w
3
)
]
ψ
1
w
1
d
x
=
∫
Ω
‍
[
σ
1
-
b
11
(
x
)
w
1
-
b
12
(
x
)
(
w
2
-
u
~
(
23
)
)
hhhhhhh
-
b
13
(
x
)
(
w
3
-
u
~
(
32
)
)
]
ψ
1
w
1
d
x
.
Similarly,(38)
∫
Ω
‍
ψ
2
w
2
t
d
x
=
∫
Ω
‍
[
σ
2
-
(
a
2
(
x
)
-
b
21
(
x
)
u
~
(
13
)
-
b
23
(
x
)
u
~
(
31
)
)
hhhhhh
+
(
a
2
(
x
)
-
b
22
(
x
)
w
2
-
b
21
(
x
)
w
1
hhhhhhhhh
-
b
23
(
x
)
w
3
)
(
a
2
(
x
)
-
b
21
(
x
)
u
~
(
13
)
-
b
23
(
x
)
u
~
(
31
)
)
]
ψ
2
w
2
d
x
=
∫
Ω
‍
[
σ
2
-
b
22
(
x
)
w
2
-
b
21
(
x
)
(
w
1
-
u
~
(
13
)
)
hhhhhhh
-
b
23
(
x
)
(
w
3
-
u
~
(
31
)
)
]
ψ
2
w
2
d
x
,
(39)
∫
Ω
‍
ψ
3
w
3
t
d
x
=
∫
Ω
‍
[
σ
3
-
(
a
3
(
x
)
-
b
31
(
x
)
u
~
(
12
)
-
b
32
(
x
)
u
~
(
21
)
)
hhhhhhhh
+
(
a
3
(
x
)
-
b
33
(
x
)
w
3
-
b
31
(
x
)
w
1
hhhhhhhh
-
b
32
(
x
)
w
2
)
b
31
(
x
)
u
~
(
12
)
]
ψ
2
w
2
d
x
=
∫
Ω
‍
[
σ
3
-
b
33
(
x
)
w
3
-
b
31
(
x
)
(
w
1
-
u
~
(
12
)
)
hhhhhhh
-
b
32
(
x
)
(
w
2
-
u
~
(
21
)
)
]
ψ
3
w
3
d
x
.
SinceX is bounded,
(40)
σ
1
-
b
11
(
x
)
w
1
-
b
12
(
x
)
(
w
2
-
u
~
(
23
)
)
-
b
13
(
x
)
(
w
3
-
u
~
(
32
)
)
,
σ
2
-
b
22
(
x
)
w
2
-
b
21
(
x
)
(
w
1
-
u
~
(
13
)
)
-
b
23
(
x
)
(
w
3
-
u
~
(
31
)
)
,
σ
3
-
b
33
(
x
)
w
3
-
b
31
(
x
)
(
w
1
-
u
~
(
12
)
)
-
b
32
(
x
)
(
w
2
-
u
~
(
21
)
)
are bounded below on X
∖
S, so that the ratios (
∫
Ω
‍
ψ
i
w
i
t
d
x
/
∫
Ω
‍
ψ
i
w
i
d
x
) are bounded below. It follows that P
(
(
v
1
,
v
2
,
v
3
)
·
t
)
/
P
(
(
v
1
,
v
2
,
v
3
)
) has a strictly positive lower bounded since the expression inside the exponential in (36) is bounded away from -
∞. To see what happens as (
v
1
,
v
2
,
v
3
)
→
ω
(
S
), we must examine how
(41)
∫
Ω
‍
[
σ
1
-
b
11
(
x
)
w
1
-
b
12
(
x
)
(
w
2
-
u
~
(
23
)
)
-
b
13
(
x
)
(
w
3
-
u
~
(
32
)
)
]
ψ
1
w
1
d
x
×
(
∫
Ω
‍
ψ
1
w
1
d
x
)
-
1
+
∫
Ω
‍
[
σ
2
-
b
22
(
x
)
w
2
-
b
21
(
x
)
(
w
1
-
u
~
(
13
)
)
-
b
23
(
x
)
(
w
3
-
u
~
(
31
)
)
]
ψ
2
w
2
d
x
×
(
∫
Ω
‍
ψ
2
w
2
d
x
)
-
1
+
∫
Ω
‍
[
σ
3
-
b
33
(
x
)
w
3
-
b
31
(
x
)
(
w
1
-
u
~
(
12
)
)
-
b
32
(
x
)
(
w
2
-
u
~
(
21
)
)
]
ψ
3
w
3
d
x
×
(
∫
Ω
‍
ψ
3
w
3
d
x
)
-
1
behaves as (
v
1
,
v
2
,
v
3
)
→
(
u
1
,
u
2
,
u
3
)
∈
ω
(
S
), (
v
1
,
v
2
,
v
3
)
∈
X
∖
S. By the continuity of the semiflow π, w
1, w
2, and w
3 will be uniformly close to v
1, v
2, and v
3, respectively, for t
>
0 sufficiently small. Hence, if we can show that the expression
(42)
σ
*
(
v
1
,
v
2
,
v
3
)
=
∫
Ω
‍
[
σ
1
-
b
11
(
x
)
v
1
-
b
12
(
x
)
(
v
2
-
u
~
(
23
)
)
-
b
13
(
x
)
(
v
3
-
u
~
(
32
)
)
]
ψ
1
v
1
d
x
×
(
∫
Ω
‍
ψ
1
v
1
d
x
)
-
1
+
∫
Ω
‍
[
σ
2
-
b
22
(
x
)
v
2
-
b
21
(
x
)
(
v
1
-
u
~
(
13
)
)
-
b
23
(
x
)
(
v
3
-
u
~
(
31
)
)
]
ψ
2
v
2
d
x
×
(
∫
Ω
‍
ψ
2
v
2
d
x
)
-
1
+
∫
Ω
‍
[
σ
3
-
b
33
(
x
)
v
3
-
b
31
(
x
)
(
v
1
-
u
~
(
12
)
)
-
b
32
(
x
)
(
v
2
-
u
~
(
21
)
)
]
ψ
3
v
3
d
x
×
(
∫
Ω
‍
ψ
3
v
3
d
x
)
-
1
always has a positive liminf
⁡ as (
v
1
,
v
2
,
v
3
)
→
(
u
1
,
u
2
,
u
3
)
∈
ω
(
S
), then by (35)–(39), we have α
(
t
,
(
u
1
,
u
2
,
u
3
)
)
>
1 for (
u
1
,
u
2
,
u
3
)
∈
ω
(
S
) and t sufficiently small.
If we let(
v
1
,
v
2
,
v
3
)
→
(
0,0
,
0
), then
(43)
σ
1
-
b
11
(
x
)
v
1
-
b
12
(
x
)
(
v
2
-
u
~
(
23
)
)
-
b
13
(
x
)
(
v
3
-
u
~
(
32
)
)
⟶
σ
1
+
b
12
(
x
)
u
~
(
23
)
+
b
13
(
x
)
u
~
(
32
)
≥
σ
1
>
0
,
σ
2
-
b
22
(
x
)
v
2
-
b
21
(
x
)
(
v
1
-
u
~
(
13
)
)
-
b
23
(
x
)
(
v
3
-
u
~
(
31
)
)
⟶
σ
2
+
b
21
(
x
)
u
~
(
13
)
+
b
23
(
x
)
u
~
(
31
)
≥
σ
2
>
0
,
σ
3
-
b
33
(
x
)
v
3
-
b
31
(
x
)
(
v
1
-
u
~
(
12
)
)
-
b
32
(
x
)
(
v
2
-
u
~
(
21
)
)
⟶
σ
3
+
b
31
(
x
)
u
~
(
12
)
+
b
32
(
x
)
u
~
(
21
)
≥
σ
3
>
0
.
Then(44)
liminf
⁡
(
v
1
,
v
2
,
v
3
)
∈
X
∖
S
(
v
1
,
v
2
,
v
3
)
→
(
0,0
,
0
)
σ
*
(
v
1
,
v
2
,
v
3
)
≥
σ
1
+
σ
2
+
σ
3
>
0
.
As (
v
1
,
v
2
,
v
3
)
→
(
u
¯
1
,
0,0
),
(45)
∫
Ω
‍
[
σ
1
-
(
a
1
(
x
)
-
b
12
(
x
)
u
~
(
23
)
-
b
13
(
x
)
u
~
(
32
)
)
+
(
a
1
(
x
)
-
b
11
(
x
)
v
1
-
b
12
(
x
)
v
2
-
b
13
(
x
)
v
3
)
b
12
(
x
)
u
~
(
23
)
]
ψ
1
v
1
d
x
⟶
∫
Ω
‍
[
σ
1
-
(
a
1
(
x
)
-
b
12
(
x
)
u
~
(
23
)
-
b
13
(
x
)
u
~
(
32
)
)
hhhhhhh
+
(
a
1
(
x
)
-
b
11
(
x
)
u
¯
1
)
b
13
(
x
)
u
~
(
32
)
]
ψ
1
u
¯
1
d
x
.
Also
(46)
∫
Ω
‍
(
a
1
(
x
)
-
b
11
(
x
)
u
¯
1
)
ψ
1
u
¯
1
d
x
=
-
∫
Ω
‍
ψ
1
μ
1
Δ
u
¯
1
d
x
=
-
∫
Ω
‍
u
¯
1
μ
1
Δ
ψ
1
d
x
=
-
∫
Ω
‍
u
¯
1
ψ
1
[
σ
1
-
(
a
1
(
x
)
-
b
12
(
x
)
u
~
(
23
)
hhhhhhhhhhhhhhhhh
-
b
13
(
x
)
u
~
(
32
)
)
]
d
x
so that as (
v
1
,
v
2
,
v
3
)
→
(
u
¯
1
,
0,0
),
(47)
∫
Ω
[
σ
1
-
b
11
(
x
)
v
1
-
b
12
(
x
)
(
v
2
-
u
~
(
23
)
)
‍
hhhh
-
b
13
(
x
)
(
v
3
-
u
~
(
32
)
)
]
ψ
1
v
1
d
x
⟶
0
.
Since ∫
Ω
‍
ψ
1
v
1
d
x
→
∫
Ω
‍
ψ
1
u
¯
1
d
x
>
0, the first ratio of integrals in (42) has limit 0 as (
v
1
,
v
2
,
v
3
)
→
(
u
¯
1
,
0,0
). For the latter two ratios of integrals, note that as (
v
1
,
v
2
,
v
3
)
→
(
u
¯
1
,
0,0
).
By the assumption (33), we have
(48)
σ
2
-
b
21
(
x
)
(
u
¯
1
-
u
~
(
13
)
)
+
b
23
(
x
)
u
~
(
31
)
>
0
,
σ
3
-
b
31
(
x
)
(
u
¯
1
-
u
~
(
12
)
)
+
b
32
(
x
)
u
~
(
21
)
>
0
.
Hence we have
(49)
liminf
⁡
(
v
1
,
v
2
,
v
3
)
∈
X
∖
S
(
v
1
,
v
2
,
v
3
)
→
(
u
¯
1
,
0,0
)
σ
*
(
v
1
,
v
2
,
v
3
)
>
0
.
As (
v
1
,
v
2
,
v
3
)
→
(
0
,
u
¯
2
,
0
), computations similar to (45) and (46) show that
(50)
∫
Ω
[
σ
2
-
b
22
(
x
)
v
2
-
b
21
(
x
)
(
v
1
-
u
~
(
13
)
)
‍
-
b
23
(
x
)
(
v
3
-
u
~
(
31
)
)
]
ψ
2
v
2
d
x
×
(
∫
Ω
‍
ψ
3
v
3
d
x
)
-
1
⟶
0
.
By assumption (33), we have
(51)
σ
1
-
b
12
(
x
)
(
u
¯
2
-
u
~
(
23
)
)
+
b
13
(
x
)
u
~
(
32
)
>
0
,
σ
3
+
b
31
(
x
)
u
~
(
12
)
-
b
32
(
x
)
(
u
¯
2
-
u
~
(
21
)
)
>
0
.
Hence we have(52)
liminf
⁡
(
v
1
,
v
2
,
v
3
)
∈
X
∖
S
(
v
1
,
v
2
,
v
3
)
→
(
0
,
u
¯
2
,
0
)
σ
*
(
v
1
,
v
2
,
v
3
)
>
0
.
As (
v
1
,
v
2
,
v
3
)
→
(
0,0
,
u
¯
3
), computations similar to (45) and (46) show that
(53)
∫
Ω
[
σ
3
-
b
33
(
x
)
v
3
-
b
31
(
x
)
(
v
1
-
u
~
(
12
)
)
‍
-
b
32
(
x
)
(
v
2
-
u
~
(
21
)
)
]
ψ
3
v
3
d
x
×
(
∫
Ω
‍
ψ
3
v
3
d
x
)
-
1
⟶
0
.
By assumption (33), we get
(54)
σ
1
+
b
12
(
x
)
u
~
(
23
)
-
b
13
(
x
)
(
u
¯
3
-
u
~
(
32
)
)
>
0
,
σ
2
+
b
21
(
x
)
u
~
(
13
)
-
b
23
(
x
)
(
u
¯
3
-
u
~
(
31
)
)
>
0
.
Hence we have
(55)
liminf
⁡
(
v
1
,
v
2
,
v
3
)
∈
X
∖
S
(
v
1
,
v
2
,
v
3
)
→
(
0,0
,
u
¯
3
)
σ
*
(
v
1
,
v
2
,
v
3
)
>
0
.
As (
v
1
,
v
2
,
v
3
)
→
(
u
~
(
12
)
,
u
~
(
21
)
,
0
),
(56)
∫
Ω
‍
[
σ
1
-
(
a
1
(
x
)
-
b
12
(
x
)
u
~
(
23
)
-
b
13
(
x
)
u
~
(
32
)
)
+
(
a
1
(
x
)
-
b
11
(
x
)
v
1
-
b
12
(
x
)
v
2
-
b
13
(
x
)
v
3
)
(
a
1
(
x
)
-
b
12
(
x
)
u
~
(
23
)
-
b
13
(
x
)
u
~
(
32
)
)
]
ψ
1
v
1
d
x
⟶
∫
Ω
‍
[
σ
1
-
(
a
1
(
x
)
-
b
12
(
x
)
u
~
(
23
)
-
b
13
(
x
)
u
~
(
32
)
)
+
(
a
1
(
x
)
-
b
11
(
x
)
u
~
(
12
)
-
b
12
(
x
)
u
~
(
21
)
)
]
ψ
1
u
~
(
12
)
d
x
.
Also
(57)
∫
Ω
‍
[
a
1
(
x
)
-
b
11
(
x
)
u
~
(
12
)
-
b
12
(
x
)
u
~
(
21
)
]
ψ
1
u
~
(
12
)
d
x
=
-
∫
Ω
‍
ψ
1
μ
1
Δ
u
~
(
12
)
d
x
=
-
∫
Ω
‍
u
~
(
12
)
μ
1
Δ
ψ
1
d
x
=
-
∫
Ω
‍
u
~
(
12
)
ψ
1
[
σ
1
-
(
a
1
(
x
)
-
b
12
(
x
)
u
~
(
23
)
-
b
13
(
x
)
u
~
(
32
)
)
]
d
x
so that as (
v
1
,
v
2
,
v
3
)
→
(
u
~
(
12
)
,
u
~
(
21
)
,
0
),
(58)
∫
Ω
[
σ
1
-
b
11
(
x
)
v
1
-
b
12
(
x
)
(
v
2
-
u
~
(
23
)
)
‍
-
b
13
(
x
)
(
v
3
-
u
~
(
32
)
)
]
ψ
1
v
1
d
x
×
(
∫
Ω
‍
ψ
1
v
1
d
x
)
-
1
⟶
0
.
For the second ratio of integrals,
(59)
∫
Ω
‍
[
σ
2
-
(
a
2
(
x
)
-
b
21
(
x
)
u
~
(
13
)
-
b
23
(
x
)
u
~
(
31
)
)
+
(
a
2
(
x
)
-
b
21
(
x
)
v
1
-
b
22
(
x
)
v
2
-
b
23
(
x
)
v
3
)
]
ψ
2
v
2
d
x
⟶
∫
Ω
‍
[
σ
2
-
(
a
2
(
x
)
-
b
21
(
x
)
u
~
(
13
)
-
b
23
(
x
)
u
~
(
31
)
)
hhhhhhh
+
(
a
2
(
x
)
-
b
21
(
x
)
u
~
(
12
)
-
b
22
(
x
)
u
~
(
21
)
)
]
ψ
2
u
~
(
21
)
d
x
.
Also
(60)
∫
Ω
[
a
2
(
x
)
-
b
21
(
x
)
u
~
(
12
)
-
b
22
(
x
)
u
~
(
21
)
]
ψ
2
u
~
(
12
)
d
x
=
-
∫
Ω
‍
ψ
2
μ
2
Δ
u
~
(
21
)
d
x
=
-
∫
Ω
‍
u
~
(
21
)
μ
2
Δ
ψ
2
d
x
=
-
∫
Ω
‍
u
~
(
21
)
ψ
2
[
σ
2
-
(
a
2
(
x
)
-
b
21
(
x
)
u
~
(
13
)
-
b
23
(
x
)
u
~
(
31
)
)
]
d
x
so that as (
v
1
,
v
2
,
v
3
)
→
(
u
~
(
12
)
,
u
~
(
21
)
,
0
),
(61)
∫
Ω
‍
[
σ
2
-
b
22
(
x
)
v
2
-
b
21
(
x
)
(
v
1
-
u
~
(
13
)
)
-
b
23
(
x
)
(
v
3
-
u
~
(
31
)
)
]
ψ
2
v
2
d
x
×
(
∫
Ω
‍
ψ
2
v
2
d
x
)
-
1
⟶
0
.
For the third ratio of integrals, as (
v
1
,
v
2
,
v
3
)
→
(
u
~
(
12
)
,
u
~
(
21
)
,
0
),
(62)
σ
3
-
b
33
(
x
)
v
3
-
b
31
(
x
)
(
v
1
-
u
~
(
12
)
)
-
b
32
(
x
)
(
v
2
-
u
~
(
21
)
)
⟶
σ
3
>
0
.
Hence we have
(63)
liminf
⁡
(
v
1
,
v
2
,
v
3
)
∈
X
∖
S
(
v
1
,
v
2
,
v
3
)
→
(
u
~
(
12
)
,
u
~
(
12
)
,
0
)
σ
*
(
v
1
,
v
2
,
v
3
)
≥
σ
3
>
0
.
Using the same method, we can get the conclusion about the other two cases.
As(
v
1
,
v
2
,
v
3
)
→
(
u
~
(
13
)
,
0
,
u
~
(
31
)
),
(64)
liminf
⁡
(
v
1
,
v
2
,
v
3
)
∈
X
∖
S
(
v
1
,
v
2
,
v
3
)
→
(
u
~
(
13
)
,
0
,
u
~
(
31
)
)
σ
*
(
v
1
,
v
2
,
v
3
)
≥
σ
2
>
0
.
and as (
v
1
,
v
2
,
v
3
)
→
(
0
,
u
~
(
23
)
,
u
~
(
32
)
),
(65)
liminf
⁡
(
v
1
,
v
2
,
v
3
)
∈
X
∖
S
(
v
1
,
v
2
,
v
3
)
→
(
0
,
u
~
(
23
)
,
u
~
(
32
)
)
σ
*
(
v
1
,
v
2
,
v
3
)
≥
σ
1
>
0
.
In conclusion
(66)
liminf
⁡
(
v
1
,
v
2
,
v
3
)
∈
X
∖
S
(
v
1
,
v
2
,
v
3
)
→
ω
(
S
)
σ
*
(
v
1
,
v
2
,
v
3
)
>
0
so that α
(
t
,
(
u
1
,
u
2
,
u
3
)
)
>
1 for (
u
1
,
u
2
,
u
3
)
∈
ω
(
S
) and t sufficiently small, and permanence follows from the abstract results of the previous section.Now we are ready to prove an interesting result, namely, permanent existence always holds among three “very weak” competing species. The result may be known, but to our knowledge, it does not appear in the published literature.Theorem 12.
Suppose thata
i
(
x
)
>
μ
i
λ
1, b
i
i
(
x
)
>
0 for x
∈
Ω
¯, i
=
1,2
,
3. Then there is a constant δ
>
0, such that if b
i
j
<
δ for i, j
=
1,2
,
3 with i
≠
j, the semiflow on [
C
0
1
(
Ω
¯
)
]
3 generated by (1) under homogeneous Dirichlet boundary conditions is permanent.Proof.
For fixeda
i
>
μ
i
λ
1 and b
i
i
(
x
)
>
0 in Ω
¯, i
=
1,2
,
3, when b
i
j
≥
0 (i, j
=
1,2
,
3, i
≠
j) are sufficiently small, conditions (11)–(13) and (17)–(19) are satisfied, and by the variational expression of principal eigenvalues (see [13]) and the boundedness of u
~
(
i
j
)
(
i
≠
j
) (see Remark 8), the principal eigenvalues σ
1, σ
2, and σ
3 of (30), (31), and (32), respectively, are all positive. And then it is easy to see that condition (33) is satisfied when b
i
j (i
≠
j) are even smaller, if necessary. The proof is completed by Theorem 11.


## 4. Summary and Discussion

The important aspect of the approach of our main result, Theorem11, is that it produces conditions that in terms of coefficients describing the strength of interactions. The criterions depend only on the behavior near the species equilibria of the two-dimensional subsystems, which is relatively easy to examine experimentally. The eigenvalues depending on the coefficients can be computed explicitly in simple cases, and there is a huge literature available on their estimation, numerical approximation, qualitative properties, and so on. Eigenvalue-based criteria for coexistence can then be used to reach biological conclusions that have been difficult to deduce through other methods.Theorem12 tell us that if the intercompeting coefficients are sufficiently small, the system (1) is always permanent, provided that a
i
>
μ
i
λ
1 in Ω
¯, which is, by Lemma 5, every single species can survive when the other two species become extinct. This sounds reasonable, and we guess that this is true for systems of any finite multiple competing species. A related topic on the systems of finite multiple competing species with constant coefficients and Neumann boundary condition can be found in [17], Theorem 9.The significant feature of our conditions for permanence is that they do not require any special assumptions the uniqueness of coexistence states, the existence of a global Lyapunov function, or other such properties, and so the methods used here can be used to other more general situations.



---

*Source: 376919-2014-08-12.xml*