Positive Solutions of a Diffusive Predator-Prey System including

0 downloads 0 Views 2MB Size Report
Sep 14, 2016 - Proposition 1 (see [16–18]). Assume that ... Theorem 3 (see [15]). Suppose ...... [18] M. Wang and Q. Wu, “Positive solutions of a prey-predator.
Missing:
Hindawi Publishing Corporation Discrete Dynamics in Nature and Society Volume 2016, Article ID 2323752, 10 pages http://dx.doi.org/10.1155/2016/2323752

Research Article Positive Solutions of a Diffusive Predator-Prey System including Disease for Prey and Equipped with Dirichlet Boundary Condition Xiaoqing Wen,1 Yue Chen,2 and Hongwei Yin1,3 1

School of Sciences, Nanchang University, Nanchang 330031, China School of Civil Engineering and Architecture, Nanchang University, Nanchang 330031, China 3 Numerical Simulation and High-Performance Computing Laboratory, School of Sciences, Nanchang University, Nanchang 330031, China 2

Correspondence should be addressed to Hongwei Yin; [email protected] Received 30 June 2016; Accepted 14 September 2016 Academic Editor: Sanja Konjik Copyright © 2016 Xiaoqing Wen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. We study a three-dimensional system of a diffusive predator-prey model including disease spread for prey and with Dirichlet boundary condition and Michaelis-Menten functional response. By semigroup method, we are able to achieve existence of a global solution of this system. Extinction of this system is established by spectral method. By using bifurcation theory and fixed point index theory, we obtain existence and nonexistence of inhomogeneous positive solutions of this system in steady state.

1. Introduction In recent decade, predator-prey systems are significant in mathematical biology and applications [1–3]. These systems are usually described by ordinary differential equation (ODEs). Besides, epidemic models have been deeply studied by differential equations in [4, 5], based on the pioneer work of the classical SIR model by Kermack and Mckendrick. However, epidemic disease can spread within a species except for interactions among species. For this case, some ecoepidemiological models were introduced in [6–8] to study how disease affects the dynamic of predator-prey systems. Since we live in the real world, spatial factor plays an important role in populations dynamic. Moreover, space could induce some significant phenomena different from the corresponding ODE systems. At present, a lot of literature has showed that predator-prey systems with spatial diffusion can produce a variety of pattern formation by Turing bifurcation [9–11] and references wherein. In this paper, we are interested

in the following model having Holling-type II functional response with disease in the prey: 𝑢𝑡 − Δ𝑢 = 𝑢 (𝑎 − 𝑢 − V) − 𝛽𝑢V,

(𝑥, 𝑡) ∈ Ω × R+ ,

𝑐𝑤V , (𝑥, 𝑡) ∈ Ω × R+ , 𝑚𝑤 + V 𝑒𝑤V 𝑤𝑡 − Δ𝑤 = −𝑑𝑤 + , (𝑥, 𝑡) ∈ Ω × R+ . 𝑚𝑤 + V V𝑡 − ΔV = 𝛽𝑢V − 𝑏V −

(1a)

There are two species in the model: prey (𝑢+V) and predator (𝑤). In (1a), 𝑢 is susceptible prey, and V denotes infected prey. Only susceptible prey 𝑢 is capable of reproducing and its reproducing rate is 𝑎. However, the infected prey does possess the capability of reproducing and still contributes to population growth of susceptible prey 𝑢 according to the logistic growth law. The disease is transmitted only within prey in the form of 𝛽𝑢V, and 𝛽 ∈ (0, 1) is called the transmission coefficient. The disease is not genetically inherited. The infected populations do not recover or become immune. Infected

2

Discrete Dynamics in Nature and Society

individuals are less active and can be caught more easily by predator resulting in increasing survival of susceptible prey. Catching rate of predators is a ratio-dependent MichaelisMenten functional response function; that is, 𝑐𝑤V/(𝑚𝑤 + V). Infected prey has a death rate 𝑏 > 0 and 𝑑 > 0 is a death rate of predators. Prey and predators are inhabited in a bounded region Ω ⊆ R𝑁 with the smooth boundary 𝜕Ω, where 𝑁 is a spatial dimension number. Δ is the Laplace operator and represents diffusion of individuals from a high density region to low density one. Zhang et al. [12] have studied (1a) under the homogeneous Neumann boundary condition and obtain permanence and stability. However, few author, as far as we know, investigated mathematical property of (1a) equipped with Dirichlet boundary condition, such as existence of global timedependent positive solution and inhomogeneous positive solutions in steady state. To this end, in this paper we make the first attempt to fill this gap and study existence of a global positive solution (1a) with homogeneous Dirichlet boundary conditions: 𝑢 (𝑥, 𝑡) = V (𝑥, 𝑡) = 𝑤 (𝑥, 𝑡) = 0,

(𝑥, 𝑡) ∈ 𝜕Ω × R+ ,

(1b)

which imply that predator and prey die out on the boundary 𝜕Ω and initial values 𝑢 (𝑥, 0) = 𝑢0 (𝑥) ≥ 0, V (𝑥, 0) = V0 (𝑥) ≥ 0, 𝑤 (𝑥, 0) = 𝑤0 (𝑥) ≥ 0,

(1c)

𝑥 ∈ Ω.

has an infinite sequence of eigenvalues. Let 𝜆 1 (𝑞) be the principal eigenvalue of (9) which is simple and its eigenfunction does not change sign in Ω. We denote 𝜆 1 (0) by 𝜆 1 for simplicity and the eigenfunction corresponding to 𝜆 1 is denoted by 𝜑1 (𝑥). 𝑞1 (𝑥) ≤ 𝑞2 (𝑥) with 𝑞1 being not identically equal to 𝑞2 implies 𝜆 1 (𝑞1 ) ≤ 𝜆 1 (𝑞2 ). Now, consider the nonlinear boundary value problem −Δ𝑢 = 𝑢 (𝑎 − 𝑢) , 𝑥 ∈ Ω, 𝑢 = 0,

𝑥 ∈ 𝜕Ω.

(4)

If 𝑎 > 𝜆 1 then (4) has a unique positive solution (cf. [15, 16]). We denote this unique positive solution by 𝑢𝑎 and 𝑢𝑎 < 𝑎. Next, the fixed point index theory, which is used later, is introduced. Let 𝐸 be a real Banach space. 𝑊 is called a wedge if 𝑊 is a closed convex set and 𝛽𝑊 ⊂ 𝑊 for all 𝛽 ≥ 0. For 𝑦 ∈ 𝑊, we define 𝑊𝑦 = {𝑥 ∈ 𝐸 : ∃𝑟 = 𝑟 (𝑥) > 0, s.t. 𝑦 + 𝑟𝑥 ∈ 𝑊} , 𝑆𝑦 = {𝑥 ∈ 𝑊𝑦 : −𝑥 ∈ 𝑊𝑦 } .

(5)

We always assume that 𝐸 = 𝑊 − 𝑊. Let 𝑇 : 𝑊𝑦 → 𝑊𝑦 be a compact linear operator on 𝐸. We say that 𝑇 has property 𝛼 on 𝑊𝑦 if there exists 𝑡 ∈ (0, 1) and 𝑤 ∈ 𝑊𝑦 \ 𝑆𝑦 such that (1 − 𝑡𝑇)𝑤 ∈ 𝑆𝑦 . Let 𝐴 : 𝑊 → 𝑊 be a compact operator with a fixed point 𝑦 ∈ 𝑊 and 𝐴 is Fr´echet differentiable at 𝑦. Let 𝐿 = 𝐴󸀠 (𝑦) be the Fr´echet derivative of 𝐴 at 𝑦. Then 𝐿 maps 𝑊𝑦 into itself. We denote by index𝑊(𝐴, 𝑦) the fixed point index of 𝐴 at 𝑦 relative to 𝑊.

Besides, we first attempt to investigate existence and nonexistence of inhomogeneous positive solutions of (1a), (1b), and (1c) in steady state as follows:

Proposition 1 (see [16–18]). Assume that 𝐼 − 𝐿 is invertible on 𝑊𝑦 . Then

−Δ𝑢 = 𝑢 (𝑎 − 𝑢 − V) − 𝛽𝑢V, 𝑥 ∈ Ω,

(i) if 𝐿 has property 𝛼 on 𝑊𝑦 , then index𝑊(𝐴, 𝑦) = 0;

𝑐𝑤V , 𝑥 ∈ Ω, 𝑚𝑤 + V 𝑒𝑤V −Δ𝑤 = −𝑑𝑤 + , 𝑥 ∈ Ω, 𝑚𝑤 + V

(ii) if 𝐿 does not have property 𝛼 on 𝑊𝑦 , then index𝑊(𝐴, 𝑦) = (−1)𝜎 , where 𝜎 is the sum of multiplicities of all eigenvalues of 𝐿 which is greater than 1.

−ΔV = 𝛽𝑢V − 𝑏V −

𝑢 = V = 𝑤 = 0,

(2)

For a linear operator 𝐿, we denote by 𝑟(𝐿) the spectral radius of 𝐿.

𝑥 ∈ 𝜕Ω.

The remainder of the paper is organized as follows: in Section 2, some necessary notations and theoretical results are introduced; in Section 3, we examine existence of a global positive solution to (1a), (1b), and (1c); in Section 4 we discuss extinction of system (1a), (1b), and (1c); existence and nonexistence of inhomogeneous positive solutions are discussed by bifurcation theory and degree theory in Section 5.

Proposition 2 (see [16, 18]). Let 𝑞(𝑥) ∈ 𝐶(Ω) and 𝑀 be a positive constant such that 𝑀 − 𝑞(𝑥) > 0 on Ω. Let 𝜆 1 (𝑞) be the first eigenvalue of the problem −Δ𝑢 + 𝑞 (𝑥) 𝑢 = 𝜆𝑢, 𝑥 ∈ Ω, 𝑢 = 0,

𝑥 ∈ 𝜕Ω.

(6)

We have the following conclusions:

2. Preliminary In this section, some notations and basic well-known results are introduced. For any 𝑞 ∈ 𝐶(Ω), the linear eigenvalue problem [13, 14], −Δ𝑢 + 𝑞 (𝑥) 𝑢 = 𝜆𝑢, 𝑥 ∈ Ω, 𝑢 = 0, 𝑥 ∈ 𝜕Ω,

(3)

(a) 𝜆 1 (𝑞) < 0 ⇒ 𝑟[(𝑀 − Δ)−1 (𝑀 − 𝑞(𝑥))] > 1, (b) 𝜆 1 (𝑞) > 0 ⇒ 𝑟[(𝑀 − Δ)−1 (𝑀 − 𝑞(𝑥))] < 1, (c) 𝜆 1 (𝑞) = 0 ⇒ 𝑟[(𝑀 − Δ)−1 (𝑀 − 𝑞(𝑥))] = 1. Now, we introduce the degree calculations, which was introduced by Dancer and Du in [19]. Assume that 𝐸1 and 𝐸2

Discrete Dynamics in Nature and Society

3

are ordered Banach spaces with positive cones 𝑊1 and 𝑊2 , respectively. Let 𝐸 = 𝐸1 ⊕ 𝐸2 and 𝑊 = 𝑊1 ⊕ 𝑊2 . Then 𝐸 is an ordered Banach space with positive cone 𝑊. Let 𝐷 be an open set in 𝑊 containing 0 and 𝐹𝑖 : 𝐷 → 𝑊𝑖 be completely continuous operators, 𝑖 = 1, 2. Denote by (𝑢, V) a general element in 𝑊 with 𝑢 ∈ 𝑊1 and V ∈ 𝑊2 . Let 𝐹 : 𝐷 → 𝑊 be defined by 𝐹 (𝑢, V) = (𝐹1 (𝑢, V) , 𝐹2 (𝑢, V)) .

(7)

Define 𝑊2 (𝜖) = {V ∈ 𝑊2 : ‖V‖𝐸2 < 𝜖}. Theorem 3 (see [15]). Suppose that 𝑈 ⊂ 𝑊1 ∩ 𝐷 is relatively open and bounded and 𝐹1 (𝑢, 0) ≠ 𝑢

𝑓𝑜𝑟 𝑢 ∈ 𝜕𝑈;

𝐹2 (𝑢, 0) = 0

𝑓𝑜𝑟 𝑢 ∈ 𝑈.

(8)

(i) deg𝑊(𝐼 − 𝐹, 𝑈 × 𝑊2 (𝜀), 0) = 0 for 𝜀 > 0 small, if, for any 𝑢 ∈ Φ, the spectral radius 𝑟(𝐹2󸀠 (𝑢, 0)|𝑊2 ) > 1 and 1 is not an eigenvalue of 𝐹2󸀠 (𝑢, 0)|𝑊2 corresponding to a positive eigenvector; (ii) deg𝑊(𝐼 − 𝐹, 𝑈 × 𝑊2 (𝜀), 0) = deg𝑊(𝐼 − 𝐹1 |𝑊1 , 𝑈, 0) for 𝜀 > 0 small, if, for any 𝑢 ∈ Φ, the spectral radius 𝑟(𝐹2󸀠 (𝑢, 0)|𝑊2 ) < 1.

3. Existence of a Global Positive Solution of (1a), (1b), and (1c) For convenience, we introduce some notations. For 𝑥, 𝑦 ∈ R, we denote 𝑥 ∧ 𝑦 fl min{𝑥, 𝑦}, 𝑥 ∨ 𝑦 fl max{𝑥, 𝑦}, 𝑥+ fl 𝑥 ∨ 0, and 𝑥− fl (−𝑥) ∨ 0 and extend these notations to real-valued functions. If (𝑉, ≥) is partially ordered vector space, we denote its positive cone by 𝑉+ fl {V ∈ 𝑉 : V ≥ 0}. For 𝑝 ∈ [1, +∞), 𝑝 > 𝑁/2, let 𝑋 = {𝐿𝑝 (Ω)}3 . Assume 𝑧 = (𝑢, V, 𝑤) ∈ 𝑋 and its norm is defined as ‖𝑧‖ = ‖𝑢‖𝑝 + ‖V‖𝑝 + ‖𝑤‖𝑝 . Hence, system (1a), (1b), and (1c) can be rewritten as an abstract differential equation: 𝑡 ∈ R+ ,

(9)

𝑧 (0) = 𝑧0 , 𝑡 = 0, where 𝑧 = (𝑢(𝑥, 𝑡), V(𝑥, 𝑡), 𝑤(𝑥, 𝑡))𝑇 and 𝑧0 = (𝑢0 , V0 , 𝑤0 )𝑇 ,

0 0 Δ

= 0} , 𝑢 (𝑎 − 𝑢 − V) − 𝛽𝑢V 𝑐𝑤V 𝐻 (𝑧) = (𝛽𝑢V − 𝑏V − 𝑚𝑤 + V ) . 𝑒𝑤V −𝑑𝑤 + 𝑚𝑤 + V

(10)

(11)

Next, the existence of the local solution for (9) will be proved. First, we give out the following lemma. Lemma 4. For every 𝑧0 ∈ 𝑋+ , Cauchy problem (9) has a unique maximal local solution

∩ 𝐶 ((0, 𝑇max ) ; 𝐷 (𝐴 𝑝 )) ,

(12)

which satisfies the following Duhamel formula for 𝑡 ∈ [0, 𝑇max ): 𝑡

𝑧 (𝑡) = 𝑒𝑡𝐴 𝑝 𝑧0 + ∫ 𝑒(𝑡−𝑠)𝐴 𝑝 𝐻 [𝑧 (𝑠)] 𝑑𝑠, 0

(13)

where 𝐴 𝑝 is a operator and the closure of 𝐴 in 𝑋. Moreover, if 𝑇max < ∞ then lim sup𝑇max − ‖𝑧(𝑡)‖ = ∞. Proof. Since the operator 𝐴 𝑝 is the closure of 𝐴 in 𝑋, it generates an analytic, condensed, strong continuous operator semigroup 𝑒𝑡𝐴 𝑝 . Moreover, we observe that 𝐻 : 𝐷(𝐴 𝑝 ) → 𝑋 is locally Lipschitz on a bounded set [20]. Via Theorem 7.2.1 in [21], we complete this proof. Lemma 5. For initial-boundary problem (1a), (1b), and (1c), the components of its solution 𝑢(𝑥, 𝑡), V(𝑥, 𝑡), 𝑤(𝑥, 𝑡) are nonnegative. Proof. To prove 𝑢(𝑥, 𝑡), V(𝑥, 𝑡), 𝑤(𝑥, 𝑡) ≥ 0 for 𝑡 ∈ [0, 𝑇max ) and 𝑥 ∈ Ω, we consider the following system: 𝜕𝑢󸀠 = Δ𝑢󸀠 + 𝑢+󸀠 (𝑎 − 𝑢+󸀠 − V+󸀠 ) − 𝛽𝑢+󸀠 V+󸀠 , 𝜕𝑡

(14a)

𝑐𝑤+󸀠 V+󸀠 𝜕V󸀠 , = ΔV󸀠 + 𝛽𝑢+󸀠 V+󸀠 − 𝑏V+󸀠 − 𝜕𝑡 𝑚𝑤+󸀠 + V+󸀠

(14b)

𝑒𝑤+󸀠 V+󸀠 𝜕𝑤󸀠 = Δ𝑤󸀠 − 𝑑𝑤+󸀠 + 𝜕𝑡 𝑚𝑤+󸀠 + V+󸀠

(14c)

with boundary and initial conditions 󵄨 󵄨 󵄨 𝑢󸀠 󵄨󵄨󵄨󵄨𝜕Ω = V󸀠 󵄨󵄨󵄨󵄨𝜕Ω = 𝑤󸀠 󵄨󵄨󵄨󵄨𝜕Ω = 0, 𝑢󸀠 (𝑥, 0) = 𝑢0 ,

Δ 0 0 𝐴 = (0 Δ 0)

3

𝐷 (𝐴) = {(𝑢, V, 𝑤) ∈ {𝐶2 (Ω)} , 𝑢|𝜕Ω = V|𝜕Ω = 𝑤|𝜕Ω

𝑧 ∈ 𝐶 ([0, 𝑇max ) ; 𝑋) ∩ 𝐶1 ((0, 𝑇max ) ; 𝑋)

Suppose that 𝐹2 : 𝐷 → 𝑊2 extends to a continuously differentiable mapping of a neighborhood of 𝐷 into 𝐸2 , 𝑊2 − 𝑊2 is dense in 𝐸2 , and Φ = {𝑢 ∈ 𝑈 : 𝑢 = 𝐹1 (𝑢, 0)}. Then the following are true:

𝑑𝑧 = 𝐴𝑧 + 𝐻 (𝑧) , 𝑑𝑡

with

V󸀠 (𝑥, 0) = V0 , 𝑤󸀠 (𝑥, 0) = 𝑤0 .

(15)

4

Discrete Dynamics in Nature and Society

After multiplying (14a), (14b), and (14c) by 𝑢−󸀠 , V−󸀠 , and 𝑤−󸀠 , respectively, and integrating by parts on the domain Ω, we obtain −

1 𝑑 󵄨 󵄨2 󵄨 󵄨2 ∫ 󵄨󵄨󵄨𝑢󸀠 󵄨󵄨󵄨 𝑑𝑥 − ∫ 󵄨󵄨󵄨󵄨∇𝑢−󸀠 󵄨󵄨󵄨󵄨 𝑑𝑥 = 0, 2 𝑑𝑡 Ω 󵄨 − 󵄨 Ω

1 𝑑 󵄨 󵄨2 󵄨 󵄨2 − ∫ 󵄨󵄨󵄨󵄨V−󸀠 󵄨󵄨󵄨󵄨 𝑑𝑥 − ∫ 󵄨󵄨󵄨󵄨∇V−󸀠 󵄨󵄨󵄨󵄨 𝑑𝑥 = 0, 2 𝑑𝑡 Ω Ω −

(𝑎 + 𝑑) 𝑒𝛽𝑐1 𝑐𝑑 (1 + 𝛽) 󵄩 󵄩 𝑒𝑐 𝑒 𝛽 󵄩󵄩𝑢0 󵄩󵄩 󵄩 󵄩 󵄩 󵄩 + 2 , ( 󵄩 󵄩∞ + 󵄩󵄩󵄩V0 󵄩󵄩󵄩∞ ) + 󵄩󵄩󵄩𝑤0 󵄩󵄩󵄩∞ } . 𝑐 𝑐 1+𝛽

𝑐3 fl max {

(21) (16)

1 𝑑 󵄨 󵄨2 󵄨 󵄨2 ∫ 󵄨󵄨󵄨𝑤󸀠 󵄨󵄨󵄨 𝑑𝑥 − ∫ 󵄨󵄨󵄨󵄨∇𝑤−󸀠 󵄨󵄨󵄨󵄨 𝑑𝑥 = 0. 2 𝑑𝑡 Ω 󵄨 − 󵄨 Ω

Proof. In view of Lemma 5, we only prove the upper bound of (1a), (1b), and (1c) for (𝑥, 𝑡) ∈ Ω × [0, 𝑇max ). From (1a), we have 𝑢𝑡 ≤ Δ𝑢 + 𝑢 (𝑎 − 𝑢) , in Ω × R+ , 𝑢 (𝑥, 𝑡) = 0,

Hence, for 𝑡 ∈ (0, 𝑇max ), 󵄩2 󵄩 󵄩2 󵄩 󵄩2 󵄩 󵄩2 󵄩󵄩 󸀠 󵄩󵄩𝑢− (⋅, 𝑡)󵄩󵄩󵄩 + 󵄩󵄩󵄩V−󸀠 (⋅, 𝑡)󵄩󵄩󵄩 ≤ 󵄩󵄩󵄩𝑢−󸀠 (⋅, 0)󵄩󵄩󵄩 + 󵄩󵄩󵄩V−󸀠 (⋅, 0)󵄩󵄩󵄩 󵄩2 󵄩 󵄩2 󵄩 󵄩2 󵄩 󵄩2 󵄩

(17)

Consequently,

We get 𝑢(𝑥, 𝑡) ≤ 𝑐1 , (𝑥, 𝑡) ∈ Ω × (0, 𝑇max ) by maximum principle, where 𝑐1 = max{𝑎, ‖𝑢0 ‖∞ }. Let V1 = 𝛽𝑢/(1 + 𝛽) + V. From (1b) we have (V1 )𝑡 − ΔV1 =

󸀠

𝑢 (𝑥, 𝑡) ≥ 0, (18)

𝑤󸀠 (𝑥, 𝑡) ≥ 0. Now, we can know that (𝑢󸀠 (𝑥, 𝑡), V󸀠 (𝑥, 𝑡), 𝑤󸀠 (𝑥, 𝑡)) is a solution of (1a), (1b), and (1c), and according to Lemma 4, we finally get that

(19)

𝑤 (𝑥, 𝑡) = 𝑤󸀠 (𝑥, 𝑡) ,

Next, we will prove the global existence of the positive solution for (1a), (1b), and (1c). Via Lemmas 4 and 5, we only need to show that the solution of (1a), (1b), and (1c) is uniformly bounded, that is, dissipation.

0 ≤ 𝑢 (𝑥, 𝑡) ≤ 𝑐1 ,

𝛽 (𝑎 + 𝑏) 𝑐1 − 𝑏V1 . 1+𝛽

(24)

fl 𝑐2 , 𝑡 ∈ [0, 𝑇) .



(20)

for (𝑥, 𝑡) ∈ Ω × [0, 𝑇max ), where 𝑐𝑖 , 𝑖 = 1, 2, 3 are defined by

𝑒𝛽 𝑏𝑒 𝑢 (𝑎 − 𝑢) − V − 𝑑𝑤 𝑐 𝑐 (1 + 𝛽) 𝑎𝑒𝛽 𝑢 𝑐 (1 + 𝛽) 𝑒𝛽𝑢 𝑒V − ] − 𝑑 [𝑤1 − 𝑐 𝑐 (1 + 𝛽)



0 ≤ 𝑤 (𝑥, 𝑡) ≤ 𝑐3

󵄩 󵄩 𝛽 (𝑎 + 𝑏) 𝑐1 𝛽 󵄩󵄩󵄩𝑢0 󵄩󵄩󵄩∞ 󵄩󵄩 󵄩󵄩 , + 󵄩󵄩V0 󵄩󵄩∞ } , 1+𝛽 𝑏 (1 + 𝛽)



(23)

V (𝑥, 𝑡) ≤ V1 (𝑥, 𝑡)

Theorem 6. Let (𝑢(𝑥, 𝑡), V(𝑥, 𝑡), 𝑤(𝑥, 𝑡)) be the solution of (1a), (1b), and (1c) with initial value (1c); then

𝑐2 fl max {

𝑎𝛽𝑢 𝛽𝑢 − 𝑏 (V1 − ) 1+𝛽 1+𝛽

Hence, by maximum principle we have that

(𝑤1 )𝑡 − Δ𝑤1 =

󵄩 󵄩 𝑐1 fl max {𝑎, 󵄩󵄩󵄩𝑢0 󵄩󵄩󵄩∞ } ,



By using the same method, we consider 𝑤1 = (𝑒/𝑐)(𝛽𝑢/(1 + 𝛽) + V) + 𝑤; then from (1c) we have that

𝑡 ∈ (0, 𝑇max ) .

0 ≤ V (𝑥, 𝑡) ≤ 𝑐2 ,

𝛽 𝑐𝑤V 𝑢 (𝑎 − 𝑢) − 𝑏V − 1+𝛽 𝑚𝑤 + V

󵄩 󵄩 𝛽 (𝑎 + 𝑏) 𝑐1 𝛽 󵄩󵄩󵄩𝑢0 󵄩󵄩󵄩∞ 󵄩󵄩 󵄩󵄩 ≤ max { , + 󵄩󵄩V0 󵄩󵄩∞ } 1+𝛽 𝑏 (1 + 𝛽)

𝑢 (𝑥, 𝑡) = 𝑢󸀠 (𝑥, 𝑡) ≥ 0, V (𝑥, 𝑡) = V󸀠 (𝑥, 𝑡) ≥ 0,

(22)

𝑢 (𝑥, 0) = 𝑢0 , in Ω.

= 0.

V󸀠 (𝑥, 𝑡) ≥ 0,

on 𝜕Ω × R+ ,

(25)

(𝑎 + 𝑑) 𝑒𝛽𝑐1 𝑑𝑒𝑐2 + − 𝑑𝑤1 . 𝑐 𝑐 (1 + 𝛽)

Thus, the maximum principle shows that 𝑤 (𝑥, 𝑡) ≤ 𝑤1 (𝑥, 𝑡) ≤ max {

(𝑎 + 𝑑) 𝑒𝛽𝑐1 𝑐𝑑 (1 + 𝛽)

󵄩 󵄩 𝑒𝑐2 𝑒 𝛽 󵄩󵄩󵄩𝑢0 󵄩󵄩󵄩∞ 󵄩󵄩 󵄩󵄩 󵄩 󵄩 + , ( + 󵄩󵄩V0 󵄩󵄩∞ ) + 󵄩󵄩󵄩𝑤0 󵄩󵄩󵄩∞ } fl 𝑐3 . 𝑐 𝑐 1+𝛽

(26)

Discrete Dynamics in Nature and Society

5

Remark 7. From Theorem 6, we know that the upper bound of the solution of (1a), (1b), and (1c) is independent of the maximal time 𝑇max , which implies that the solution of (1a), (1b), and (1c) is a global solution.

Since 𝑎 > 𝜆 1 , we get from (28) and comparison theorem that

By Lemmas 4 and 5 and Theorem 6, we can get the following theorem.

Then, there exists 𝑇𝜀 > 0 such that 𝑢(𝑥, 𝑡) ≤ 𝑢𝑎 + 𝜀 for all 𝑡 > 𝑇𝜀 . Thus, we have

Theorem 8. System (1a), (1b), and (1c) has a unique, nonnegative, and bounded solution 𝑧(𝑥, 𝑡) = (𝑢(𝑥, 𝑡), V(𝑥, 𝑡), 𝑤(𝑥, 𝑡)), such that 𝑧 (𝑥, 𝑡) ∈ 𝐶 ([0, ∞) , 𝑋) ∩ 𝐶1 ((0, ∞) , 𝑋) ∩ 𝐶 ((𝑥, ∞) , 𝐷 (𝐴)) .

(27)

4. Extinction of (1a), (1b), and (1c) Theorem 9. Let (𝑢(𝑥, 𝑡), V(𝑥, 𝑡), 𝑤(𝑥, 𝑡)) be a positive solution of (1a), (1b), and (1c), then we have the following: (i) if 𝑎 ≤ 𝜆 1 and 𝑒 < 𝑑 then (𝑢, V, 𝑤) → (0, 0, 0) as 𝑡 → ∞; (ii) if 𝑎 < 𝜆 1 and −𝑏 < 𝜆 1 (−𝛽𝑢𝑎 ) then (𝑢, V, 𝑤) → (𝑢𝑎 , 0, 0) as 𝑡 → ∞. Proof. (i) First, we observe that solution (𝑢, V, 𝑤) of (1a), (1b), and (1c) satisfies 𝑢𝑡 − Δ𝑢 ≤ 𝑢 (𝑎 − 𝑢) , 𝑢 = 0,

(𝑥, 𝑡) ∈ Ω × R+ ,

(𝑥, 𝑡) ∈ 𝜕Ω × R+ .

(28)

Let 𝜀 be a sufficiently small positive constant such that 𝜀 ≤ 𝑏/𝛽. Since 𝑎 ≤ 𝜆 1 , we see that 𝑢(𝑥, 𝑡) → 0 uniformly as 𝑡 → ∞ by using comparison principle for elliptic problems. Then there exists 𝑇𝜀 ≥ 0 such that 𝑢(𝑥, 𝑡) ≤ 𝜀 for all 𝑡 ≥ 𝑇𝜀 . Therefore, we have V𝑡 − ΔV ≤ V (𝛽𝜀 − 𝑏) , V = 0,

(𝑥, 𝑡) ∈ Ω × (𝑇𝜀 , ∞) ,

(𝑥, 𝑡) ∈ 𝜕Ω × (𝑇𝜀 , ∞) ,

(29)

which concludes that V(𝑥, 𝑡) → 0 uniformly as 𝑡 → ∞. Similarly, there exists 𝑇𝜀󸀠 > 0 such that V(𝑥, 𝑡) ≤ 𝜀 for all 𝑡 > 𝑇𝜀󸀠 . Further, since 𝑒𝑤V/(𝑚𝑤 + V) is monotone increasing with respect to the variable V, we have

lim sup 𝑢 (𝑥, 𝑡) ≤ 𝑢𝑎 . 𝑡→∞

V𝑡 − ΔV ≤ V [𝛽 (𝑢𝑎 + 𝜀) − 𝑏] , V = 0,

(𝑥, 𝑡) ∈ 𝜕Ω × (𝑇𝜀 , ∞) .

𝑢𝑡 − Δ𝑢 ≥ 𝑢 (𝑎 − 𝑢 − (𝛽 + 1) 𝜀) , (𝑥, 𝑡) ∈ Ω × (𝑇𝜀󸀠 , ∞) , (34) 𝑢 = 0,

(𝑥, 𝑡) ∈ 𝜕Ω × (𝑇𝜀󸀠 , ∞) .

Since 𝜀 satisfies (b), we get that lim inf 𝑢 (𝑥, 𝑡) ≥ 𝑢(𝑎−(𝛽+1)𝜀) . 𝑡→∞

Thus, from (32), (35), and the arbitrariness of 𝜀, we can conclude that 𝑢(𝑥, 𝑡) → 𝑢𝑎 as 𝑡 → ∞. We finish this proof.

5. Existence and Nonexistence of Positive Solutions of (2) To prove the existence of positive solution of (2) by using fixed point index theory, we need a priori estimate for nonnegative solutions of (1a), (1b), and (1c). So we first give the following theorem. Theorem 10. Any nonnegative solution (𝑢(𝑥), V(𝑥), 𝑤(𝑥)) of (2) has a priori bounds: 𝑢 (𝑥) ≤ 𝑎, V (𝑥) ≤ 𝑐4 ,

which concludes that 𝑤(𝑥, 𝑡) → 0 as 𝑡 → ∞. (ii) Let 𝜀 be a sufficiently small positive constant such that

𝑎 − 𝜆1 . (b) 𝜀 < 𝛽+1

(36)

where 𝑐4 fl

(𝑥, 𝑡) ∈ 𝜕Ω × (𝑇𝜀󸀠 , ∞) ,

𝜆 1 (𝑏 − 𝛽𝑢𝑎 ) , 𝛽

(35)

𝑤 (𝑥) ≤ 𝑐5 ,

(𝑥, 𝑡) ∈ Ω × (𝑇𝜀󸀠 , ∞) , (30)

(a) 𝜀
𝑇𝜀󸀠 . Hence, we have

𝑒𝜀 𝑤𝑡 − Δ𝑤 ≤ 𝑤 ( − 𝑑) ≤ 𝑤 (𝑒 − 𝑑) , 𝑚𝑤 + 𝜀

𝑤 = 0,

(𝑥, 𝑡) ∈ Ω × (𝑇𝜀 , ∞) ,

(32)

𝛽 (𝑎 + 𝑏) 𝑎 , 𝑏 (1 + 𝛽)

𝑎 (𝑎 + 𝑑) 𝑒𝛽 𝑒𝑐4 𝑐5 fl + . 𝑐 𝑐𝑑 (1 + 𝛽)

(37)

Proof. From the first equation of (2), we have (31)

−Δ𝑢 ≤ 𝑢 (𝑎 − 𝑢) , 𝑥 ∈ Ω, 𝑢 = 0,

𝑥 ∈ 𝜕Ω,

(38)

6

Discrete Dynamics in Nature and Society

and we get 𝑢(𝑥) ≤ 𝑎 by maximum principle. Let V1 = 𝛽𝑢/(1 + 𝛽) + V with V1 |𝜕Ω = 0. From (1b), we have −ΔV1 = ≤

𝛽 𝑐𝑤V 𝑢 (𝑎 − 𝑢) − 𝑏V − 1+𝛽 𝑚𝑤 + V 𝑎𝛽𝑢 𝛽𝑢 𝛽 (𝑎 + 𝑏) 𝑎 − 𝑏 (V1 − )≤ − 𝑏V1 . 1+𝛽 1+𝛽 1+𝛽

(39)

Thus, by maximum principle, we get V (𝑥) ≤ V1 (𝑥) ≤

𝛽 (𝑎 + 𝑏) 𝑎 = 𝑐4 . 𝑏 (1 + 𝛽)

𝑒𝛽 𝑏𝑒 𝑢 (𝑎 − 𝑢) − V − 𝑑𝑤 𝑐 𝑐 (1 + 𝛽)

𝑒𝛽𝑢 𝑎𝑒𝛽 𝑒V ≤ 𝑢 − 𝑑 [𝑤1 − − ] 𝑐 𝑐 (1 + 𝛽) 𝑐 (1 + 𝛽) ≤

(i) deg𝑊(1 − 𝐹, 𝐷) = 1, (ii) index𝑊(𝐹, (0, 0, 0)) = 0. Proof. (i) It is easy to see that 𝐹 has no fixed point on 𝜕Ω, so deg𝑊(𝐼 − 𝐹, 𝐷) is well defined. For 𝑠 ∈ [0, 1], we define a positive and compact operator 𝐹𝑠 : 𝐸 → 𝐸 by 𝐹𝑠 (𝑢, V)

(40)

We set 𝑤1 = (𝑒/𝑐)(𝛽𝑢/(1 + 𝛽) + V) + 𝑤 with 𝑤1 |𝜕Ω = 0, and then we have −Δ𝑤1 =

Lemma 11. Assume that 𝑎 > 𝜆 1 , then one has

(41)

𝑠𝑢 (𝑎 − 𝑢 − V) − 𝑠𝛽𝑢V + 𝑀𝑢 𝑠𝑐𝑤V (45) = (−Δ + 𝑀)−1 (𝑠𝛽𝑢V − 𝑠𝑏V − 𝑚𝑤 + V + 𝑀V) , 𝑠𝑒𝑤V −𝑑𝑠𝑤 + + 𝑀𝑤 𝑚𝑤 + V and then 𝐹1 = 𝐹. For each 𝑠, a fixed point 𝐹𝑠 is a solution of the following problem: −Δ𝑢 = 𝑠𝑢 (𝑎 − 𝑢 − V) − 𝑠𝛽𝑢V,

𝑎 (𝑎 + 𝑑) 𝑒𝛽 𝑑𝑒𝑐4 + − 𝑑𝑤1 . 𝑐 𝑐 (1 + 𝛽)

𝑠𝑐𝑤V , 𝑥 ∈ Ω, 𝑚𝑤 + V 𝑠𝑒𝑤V −Δ𝑤 = −𝑑𝑠𝑤 + , 𝑥 ∈ Ω, 𝑚𝑤 + V −ΔV = 𝑠𝛽𝑢V − 𝑠𝑏V −

By using maximum principle again, we obtain 𝑤 ≤ 𝑤1 ≤

𝑎 (𝑎 + 𝑑) 𝑒𝛽 𝑒𝑐4 + = 𝑐5 . 𝑐 𝑐𝑑 (1 + 𝛽)

(42)

Therefore, we complete this proof.

𝑢 = V = 0,

(46)

𝑥 ∈ 𝜕Ω.

Form Theorem 10 𝐹𝑠 has no fixed point 𝜕Ω and deg𝑊(1 − 𝐹𝑠 , 𝐷) is well defined. Hence,

Now, we introduce some notations: 𝐸 = {𝐶01 (Ω)}3 , where 𝐶01 (Ω) = {𝜙 ∈ 𝐶1 (Ω) : 𝜙 = 0, on 𝜕Ω}, 3

𝑥 ∈ Ω,

deg (𝐼 − 𝐹, 𝐷) = deg𝑊 (𝐼 − 𝐹1 , 𝐷)

1

𝑊 = {𝐾} , where 𝐾 = {𝜙 ∈ 𝐶 (Ω) : 𝜙 ≥ 0 in Ω and 𝜙 = 0 on 𝜕Ω}, 𝐷 = {(𝑢, V, 𝑤) ∈ 𝑊 : 𝑢 < 𝑎 + 1, V < 𝑐4 + 1, 𝑤 < 𝑐5 + 1}, according to Theorem 10.

Note that (46) has only the trivial solution (0, 0, 0) when 𝑠 = 0. Set

Taking

𝑀 0 𝑐 𝑀 = max {𝑎2 + 𝑎𝑐4 + 𝛽𝑎𝑐4 , 𝑏𝑐4 + , 𝑑𝑐5 } , 𝑚

(43)

𝐿=

𝐹0󸀠

𝑢 (𝑎 − 𝑢 − V) − 𝛽𝑢V + 𝑀𝑢 𝑐𝑤V = (−Δ + 𝑀)−1 (𝛽𝑢V − 𝑏V − 𝑚𝑤 + V + 𝑀V) . 𝑒𝑤V −𝑑𝑤 + + 𝑀𝑤 𝑚𝑤 + V

(44)

It is obvious that 𝐹 having a positive fixed point in interior of 𝐷 implies that (2) has a positive solution. Next, we give the degree of 𝐼 − 𝐹 in 𝐷 relative to 𝐾 and the fixed point index of 𝐹 at the trivial solution (0, 0, 0) of (1a), (1b), and (1c) relative to 𝐾.

0

−1

(0, 0, 0) = (−Δ + 𝑀) ( 0 𝑀 0 ) . 0

we can define a positive and compact operator 𝐹 : 𝐸 → 𝐸 by 𝐹 (𝑢, V, 𝑤)

(47)

= deg𝑊 (𝐼 − 𝐹0 , 𝐷) .

(48)

0 𝑀

Let 𝐿(𝜉1 , 𝜉2 , 𝜉3 ) = (𝜉1 , 𝜉2 , 𝜉3 ) for some (𝜉1 , 𝜉2 , 𝜉3 ) ∈ 𝑊(0,0,0) = 𝐾 × 𝐾 × 𝐾. Obviously, (𝜉1 , 𝜉2 , 𝜉3 ) = (0, 0, 0). Hence, 𝐼 − 𝐿 has no nontrivial kernel in 𝑊(0,0,0) . Since 𝜆 1 > 0, it is easy to see that 𝑟((−Δ + 𝑀)−1 (𝑀)) < 1 by Proposition 2. It implies that 𝐿 has no property 𝛼. Therefore, by Proposition 1, we get the first result of Lemma 11. (ii) It is obvious that 𝐹(0, 0, 0) = (0, 0, 0). Let 𝐿 = 𝐹󸀠 (0, 0, 0); then 𝑎+M −1

𝐿 = (−Δ + 𝑀) ( 0 0

0

0

−𝑏 + 𝑀

0

0

−𝑑 + 𝑀

).

(49)

Discrete Dynamics in Nature and Society

7

Assume that 𝐿(𝜉1 , 𝜉2 , 𝜉3 ) = (𝜉1 , 𝜉2 , 𝜉3 ) for some (𝜉1 , 𝜉2 , 𝜉3 ) ∈ 𝑊(0,0,0) , and then −Δ𝜉1 = 𝑎𝜉1 , 𝑥 ∈ Ω, 𝜉1 = 0,

(50)

𝑥 ∈ 𝜕Ω,

−Δ𝜉2 = −𝑏𝜉2 , 𝑥 ∈ Ω, 𝜉2 = 0,

(51)

𝑥 ∈ 𝜕Ω,

−Δ𝜉3 = −𝑑𝜉3 , 𝑥 ∈ Ω, 𝜉3 = 0,

0

(52)

𝑥 ∈ 𝜕Ω.

(𝐼 − 𝑠𝐿) ( 𝜙0 )

Since 𝑏 > 0 and 𝑑 > 0, from (51) and (52), we know that 𝜉2 = 𝜉3 = 0. If 𝜉1 ∈ 𝐾 with 𝜉1 ≢ 0, then 𝑎 = 𝜆 1 from (50), which contradicts with the assumption. Thus (𝜉1 , 𝜉2 , 𝜉3 ) = (0, 0, 0). So 𝐼 − 𝐿 has nontrivial kernel in 𝑊(0,0,0) . For 𝑎 > 𝜆 1 , by Proposition 2, we see that 𝑟0 = 𝑟((−Δ + 𝑀)−1 (𝑎 + 𝑀)) > 1, and 𝑟0 is the principle eigenvalue of (−Δ + 𝑀)−1 (𝑎 + 𝑀) with a corresponding eigenfunction 𝜙0 > 0. Since 𝑆(0,0,0) = {(0, 0, 0)}, we see that (𝜙0 , (0, 0, 0)) ∈ 𝑊(0,0,0) \ 𝑆(0,0,0) . Set 𝑠 ∈ 1/𝑟0 ∈ (0, 1), and then (𝐼 − 𝑠𝐿)(𝜙0 , (0, 0, 0)) = (0, 0, 0) ∈ 𝑆(0,0,0) . Thus, 𝐿 has the property 𝛼. By Proposition 1 we have index𝑊(𝐹, (0, 0, 0)) = 0. Now, we give out the index of the semitrivial solution (𝑢𝑎 , 0, 0) of (2). Lemma 12. Assume that 𝑎 > 𝜆 1 , and then one has (i) if −𝑏 > 𝜆 1 (−𝛽𝑢𝑎 ), then index𝑊(𝐹, (𝑢𝑎 , 0, 0)) = 0; (ii) if −𝑏 < 𝜆 1 (−𝛽𝑢𝑎 ), then index𝑊(𝐹, (𝑢𝑎 , 0, 0)) = 1. Proof. Let 𝐿 = 𝐹󸀠 (𝑢𝑎 , 0, 0); then

0 0

− (1 + 𝛽) 𝑢𝑎 𝜙0 −1

= ( 𝜙0 ) − 𝑠 (−Δ + 𝑀) ((𝛽𝑢𝑎 − 𝑏 + 𝑀) 𝜙0 ) (56) 0

0

𝑠 (−Δ + 𝑀)−1 (1 + 𝛽) 𝑢𝑎 𝜙0 =(

) ∈ 𝑆(𝑢𝑎 ,0,0) .

0 0

Hence, 𝐿 has the property 𝛼 on 𝑊(𝑢𝑎 ,0,0) resulting in index𝑊(𝐹, (𝑢𝑎 , 0, 0)) = 0 according to Proposition 1. (ii) First, we prove that 𝐿 does not have the property 𝛼 on 𝑊(𝑢𝑎 ,0,0) . Since −𝑏 < 𝜆 1 (−𝛽𝑢𝑎 ), we have 𝑟((−Δ + 𝑀)−1 (−𝛽𝑢𝑎 )) < 1. On the contrary, we suppose that 𝐿 has the property 𝛼 on 𝑊(𝑢𝑎 ,0,0) . Then there exist 𝑠 ∈ (0, 1) and (𝜉1 , 𝜉2 , 𝜉3 ) ∈ 𝑊(𝑢𝑎 ,0,0) such that (𝐼 − 𝑠𝐿)(𝜉1 , 𝜉2 , 𝜉3 ) ∈ 𝑆(𝑢𝑎 ,0,0) . So (−Δ + 𝑀)−1 (𝛽𝑢𝑎 − 𝑏 + 𝑀) 𝜉2 =

−1

𝐿 = (−Δ + 𝑀)

𝑎 − 2𝑢𝑎 + 𝑀 − (1 + 𝛽) 𝑢𝑎 ⋅(

If 𝜉1 ≠ 0, then 𝜆 1 (2𝑢𝑎 − 𝑎) = 0 by Proposition 2. On the other hand, 𝜆 1 (2𝑢𝑎 −𝑎) > 𝜆 1 (𝑢𝑎 −𝑎) = 0, and we get a contradiction. Hence, (𝜉1 , 𝜉2 , 𝜉3 ) = (0, 0, 0). Then 𝐼 − 𝐿 has no nontrivial kernel in 𝑊(𝑢𝑎 ,0,0) . (i) Since −𝑏 > 𝜆 1 (−𝛽𝑢𝑎 ), by Proposition 1, we have 𝑟0 = 𝑟((−Δ + 𝑀)−1 (𝛽𝑢𝑎 − 𝑏 + 𝑀)) > 1 which is an eigenvalue of (−Δ + 𝑀)−1 (𝛽𝑢𝑎 − 𝑏 + 𝑀) with a corresponding eigenfunction 𝜙0 > 0. Since 𝑆(𝑢𝑎 ,0,0) = 𝐶01 (Ω) × {0} × {0}, we have (0, 𝜙0 , 0) ∈ 𝑊(𝑢𝑎 ,0,0) \ 𝑆(𝑢𝑎 ,0,0) . Set 𝑠 ∈ 1/𝑟0 ∈ (0, 1), and then

0

0

𝛽𝑢𝑎 − 𝑏 + 𝑀

0

0

0

−𝑑 + 𝑀

).

(53)

Let 𝐿(𝜉1 , 𝜉2 , 𝜉3 ) = (𝜉1 , 𝜉2 , 𝜉3 ) for some (𝜉1 , 𝜉2 , 𝜉3 ) ∈ 𝑊(𝑢𝑎 ,0,0) = 𝐶01 (Ω) × 𝐾 × 𝐾. Then −Δ𝜉1 + (2𝑢𝑎 − 𝑎) 𝜉1 = − (1 + 𝛽) 𝑢𝑎 𝜉2 , 𝑥 ∈ Ω, −Δ𝜉2 + (𝑏 − 𝛽𝑢𝑎 ) 𝜉2 = 0,

𝑥 ∈ Ω,

−Δ𝜉3 + 𝑑𝜉3 = 0,

𝑥 ∈ Ω,

𝜉1 = 𝜉2 = 𝜉3 = 0,

(54)

(57)

Then 1/𝑠 is an eigenvalue of the operator (−Δ + 𝑀)−1 (𝛽𝑢𝑎 − 𝑏 + 𝑀), which contradicts 𝑟((−Δ + 𝑀)−1 (𝛽𝑢𝑎 − 𝑏 + 𝑀)) < 1. So 𝐿 does not have the property 𝛼 on 𝑊(𝑢𝑎 ,0,0) . By Proposition 2, we have index𝑊(𝐹, (𝑢𝑎 , 0, 0)) = (−1)𝜎 , where 𝜎 is the sum of the multiplicities of all eigenvalues of 𝐿. Next, we will prove that 𝜎 = 0. Assume that 1/𝜌 > 1 is an eigenvalue of 𝐿 with corresponding eigenfunction (𝜉1 , 𝜉2 , 𝜉3 ), and then 𝐿(𝜉1 , 𝜉2 , 𝜉3 ) = 1/𝜌(𝜉1 , 𝜉2 , 𝜉3 ); that is, −Δ𝜉1 + 𝑀𝜉1 = 𝜌 ((𝑎 − 2𝑢𝑎 + 𝑀) 𝜉1 − (1 + 𝛽) 𝑢𝑎 𝜉2 ) , 𝑥 ∈ Ω,

𝑥 ∈ 𝜕Ω.

−Δ𝜉2 + 𝑀𝜉2 = 𝜌 (𝛽𝑢𝑎 − 𝑏 + 𝑀) 𝜉2 ,

Since 𝜉2 , 𝜉3 ∈ 𝐾, −𝑏 ≠ 𝜆 1 (−𝛽𝑢𝑎 ) and 𝑑 > 0. So we get that 𝜉2 = 𝜉3 = 0. Hence, from the first equation of (54) we have −Δ𝜉1 + (2𝑢𝑎 − 𝑎) 𝜉1 = − (1 + 𝛽) 𝑢𝑎 𝜉2 , 𝜉1 = 0, 𝑥 ∈ 𝜕Ω.

𝜉2 . 𝑠

𝑥 ∈ Ω,

(55)

−Δ𝜉3 + 𝑀𝜉3 = 𝜌 (−𝑑 + 𝑀) 𝜉3 ,

𝑥 ∈ Ω,

(58)

𝑥 ∈ Ω,

𝜉1 = 𝜉2 = 𝜉3 = 0, 𝑥 ∈ 𝜕Ω. Since 0 < 𝜌 < 1 and 𝑑 > 0, it follows from the third equation (58) that 𝜉3 = 0.

8

Discrete Dynamics in Nature and Society

We suppose that 𝜉2 ≠ 0, and it follows from the second equation of (58) and Proposition 2 that 0 = 𝜆 1 (𝑀 (1 − 𝜌) − 𝜌 (𝛽𝑢𝑎 − 𝑏)) > 𝜆 1 (−𝛽𝑢𝑎 ) + 𝑏,

(59)

which contradicts −𝑏 < 𝜆 1 (−𝛽𝑢𝑎 ), and so 𝜉2 = 0. If 𝜉1 ≠ 0, it follows from the first equation of (58) that 0 = 𝜆 1 (𝑀 (1 − 𝜌) − 𝜌 (𝑎 − 2𝑢𝑎 )) > 𝜆 1 (2𝑢𝑎 − 𝑎)

Lemma 14. Assume Ψ ≠ 0, and then one has the following: (i) if −𝑑 > 𝜆 1 − 𝑒, for any (𝑢∗ , V∗ , 0) ∈ Ψ, then deg𝑊(𝐼 − 𝐹, Ψ) = 0; (ii) if −𝑑 < 𝜆 1 − 𝑒, for any (𝑢∗ , V∗ , 0) ∈ Ψ, then deg𝑊(𝐼 − 𝐹, Ψ) = 1. Proof. Let 2

𝐸1 = {𝐶01 (Ω)} ,

(60)

> 𝜆 1 (𝑢𝑎 − 𝑎) = 0.

𝐸2 = 𝐶01 (Ω) ,

Thus, we get the contradiction and this implies that 𝜉1 = 0. Hence, (𝜉1 , 𝜉2 , 𝜉3 ) = (0, 0, 0), which implies that 𝐿 has no eigenvalue being greater than 1; that is, 𝜎 = 0, and index𝑊(𝐹, (𝑢𝑎 , 0, 0)) = 1. We complete the proof. Next, we consider the other semitrivial solutions of (2). Let us consider the following three subsystems: −Δ𝑢 = 𝑢 (𝑎 − 𝑢 − V) − 𝛽𝑢V, 𝑥 ∈ Ω, −ΔV = 𝛽𝑢V − 𝑏V,

𝑥 ∈ Ω,

(61)

𝑢 = V = 0, 𝑥 ∈ 𝜕Ω.

𝑢 = 𝑤 = 0,

𝑥 ∈ Ω,

(62)

𝑥 ∈ 𝜕Ω.

𝑐𝑤V , 𝑥 ∈ Ω, 𝑚𝑤 + V 𝑒𝑤V −Δ𝑤 = −𝑑𝑤 + , 𝑥 ∈ Ω, 𝑚𝑤 + V −ΔV = −𝑏V −

V = 𝑤 = 0,

(63)

𝑥 ∈ 𝜕Ω.

It is obvious that (63) has trivial solution (V, 𝑤) = (0, 0) and that (62) only has a trivial solution (𝑢, 𝑤) = (0, 0) and a semitrivial solution (𝑢, 𝑤) = (𝑢𝑎 , 0) if and only if 𝑎 > 𝜆 1 . For subsystem (61), a lot of papers [15, 22] and references wherein have discussed systems in this case, so we directly give out some results about (61). ∗



Theorem 13. Equation (61) has a positive solution (𝑢 , V ) with 𝑢∗ < 𝑢𝑎 if and only if 𝑎 > 𝜆 1 and −𝑏 > 𝜆 1 (−𝛽𝑢𝑎 ). Moreover, if 𝑎 − 𝛽 − 1 > 𝜆 1 and −𝑏 > 𝜆 1 (−𝛽𝑢(𝑎−𝛽−1) ), then the positive solution (𝑢∗ , V∗ ) satisfies 𝑢(𝑎−𝛽−1) ≤ 𝑢∗ and V ≤ V∗ , where V is the unique positive solution of the following equations: −ΔV = V (𝛽𝑢(𝑎−𝛽−1) − 𝑏) , 𝑥 ∈ Ω, V = 0,

𝑥 ∈ 𝜕Ω.

𝑊2 = 𝐾. Define 𝐹1 (𝑢, V, 𝑤) = (−Δ + 𝑀)−1 (

𝑢 (𝑎 − 𝑢) − 𝛽𝑢V + 𝑀𝑢 ), 𝑐V𝑤 𝛽𝑢V − 𝑏V − + 𝑀V 𝑚𝑤 + V

(66)

𝑒V𝑤 − 𝑑𝑤) . 𝑚𝑤 + V Then 𝐹 = (𝐹1 , 𝐹2 ) and we are in the framework to use Theorem 3. We choose a neighborhood 𝑈 ⊂ 𝑊1 ∩𝐷 of Ψ∩𝑊1 such that (𝑢𝑎 , 0) ∉ 𝑈. Now, 𝐹1 (𝑢, V, 0) = (𝑢, V) with (𝑢, V) ∈ 𝑈 if and only if (𝑢, V, 0) ∈ Ψ. By the results in [22], we have 󵄨 deg𝑊1 (𝐼 − 𝐹1 󵄨󵄨󵄨𝑊1 , 𝑈, 0) 𝐹2 (𝑢, V, 𝑤) = (−Δ + 𝑀)−1 (

−Δ𝑢 = 𝑢 (𝑎 − 𝑢) , 𝑥 ∈ Ω, −Δ𝑤 = −𝑑𝑤,

(65)

𝑊1 = {𝐾}2 ,

(64)

Let {(𝑢, V)} be the set of all positive solutions of (61), denoted by Φ, and Ψ = {(𝑢, V, 0), (𝑢, V) ∈ Φ}. It is easy to see that (61) has semitrivial solution (𝑢∗ , V∗ , 0) ∈ Ψ from Theorem 13. Next, we give the degree of (𝐼 − 𝐹) in Ψ relative to 𝑊.

1, if 𝑎 > 𝜆 1 , −𝑏 > 𝜆 1 (−𝛽𝑢𝑎 ) , { { { { = {−1, if 𝑎 < 𝜆 1 , −𝑏 < 𝜆 1 (−𝛽𝑢𝑎 ) , { { { if (𝑎 − 𝜆 1 ) ⋅ [−𝑏 − 𝜆 1 (−𝛽𝑢𝑎 )] < 0. {0,

(67)

For any (𝑢∗ , V∗ , 0) ∈ Ψ, a direct calculation shows that 󵄨 𝐹2󸀠 (𝑢∗ , V∗ , 0)󵄨󵄨󵄨󵄨𝑊 𝑤 = (−Δ + 𝑀)−1 (𝑀 − 𝑏 + 𝑒) 𝑤. (68) 2 Notice that, by our choice of 𝑀, 𝑀 − 𝑏 + 𝑒 > 0, so it follows from the maximum principle that 𝐹2󸀠 (𝑢∗ , V∗ , 0)|𝑊2 is 𝑢0 -positive in the sense of [23] with 𝑢0 = (−Δ)−1 . Hence, 𝑟(𝐹2󸀠 (𝑢∗ , V∗ , 0)|𝑊2 ) > 0 and is the only eigenvalue corresponding to a positive eigenvector. By a simple calculation, we can easily show that 𝑟(𝐹2󸀠 (𝑢∗ , V∗ , 0)|𝑊2 ) > 1 if and only if −𝑑 > 𝜆 1 − 𝑒 and 𝑟(𝐹2󸀠 (𝑢∗ , V∗ , 0)|𝑊2 ) < 1 if and only if −𝑑 < 𝜆 1 − 𝑒. Therefore, by Theorem 3 and the discussion above, we have deg𝑊 (𝐼 − 𝐹, 𝑈 × 𝑊2 (𝜀) , 0) if − 𝑑 > 𝜆 1 − 𝑒, (69) {0, ={ 󵄨 deg (𝐼 − 𝐹1 󵄨󵄨󵄨𝑊1 , 𝑈, 0) , if − 𝑑 < 𝜆 1 − 𝑒. { 𝑊 Again since the degree does not depend on the particular choice of 𝑈 and 𝜀, Ψ ≠ 0 means 𝑎 > 𝜆 1 and −𝑏 > 𝜆 1 (−𝛽𝑢𝑎 ). Thus, we complete this proof.

Discrete Dynamics in Nature and Society

9

From the discussion above, we can obtain the existence of positive solutions of (2). Theorem 15. Assume 𝑎 > 𝜆 1 . Equation (2) has at least one positive solution if 𝑎 > 𝜆 1 + 𝛽 + 1, −𝑏 > 𝜆 1 (−𝛽𝑢(𝑎−𝛽−1) ) ,

(70)

−𝑑 > 𝜆 1 − 𝑒, and then (2) has at least one positive solution. Proof. Since 𝑎 > 𝜆 1 , we obtain deg𝑊(𝐼 − 𝐹, 𝐷) = 1 and index𝑊(𝐹, (0, 0, 0)) = 0 from Lemma 11. Hence, we have index𝑊 (𝐹, (𝑢𝑎 , 0, 0)) + deg𝑊 (𝐼 − 𝐹, 𝐷) = 1.

(71)

Since −𝑏 > 𝜆 1 (−𝛽𝑢(𝑎−𝛽−1) ) > 𝜆 1 (−𝛽𝑢𝑎 ), we have index𝑊(𝐹, (𝑢𝑎 , 0, 0)) = 0 from Lemma 12. Moreover, Theorem 13 implies Ψ ≠ 0. Let (𝑢∗ , V∗ , 0) ∈ Ψ. Since 𝑢∗ ≥ 𝑢(𝑎−𝛽) and V∗ ≥ V by Theorem 13, we have −𝑑 > 𝜆 1 − 𝑒 and thus deg𝑊(𝐼 − 𝐹, Ψ) = 0 from Lemma 14. Therefore, index𝑊(𝐹, (𝑢𝑎 , 0, 0)) + deg𝑊(𝐼 − 𝐹, Ψ) = 0. So (71) holds. We complete this proof. Next, we give out the result of nonexistence of positive solutions. Theorem 16. If any one of the following conditions holds, then (2) has no positive solution: (i) 𝑎 < 𝜆 1 ; (ii) 𝑎 > 𝜆 1 and −𝑏 ≤ 𝜆 1 (−𝛽𝑢𝑎 ). Proof. Assume that (2) has positive solution (𝑢, V, 𝑤), and then 𝑢, V, 𝑤 satisfy the following three equations, respectively: −Δ𝑢 = 𝑢 (𝑎 − 𝑢 − V) − 𝛽𝑢V, 𝑥 ∈ Ω, 𝑢 = 0, 𝑥 ∈ 𝜕Ω, −ΔV = 𝛽𝑢V − 𝑏V −

𝑐𝑤V , 𝑥 ∈ Ω, 𝑚𝑤 + V

(72)

(73)

V = 0, 𝑥 ∈ 𝜕Ω, −Δ𝑤 =

𝑒𝑤V − 𝑑𝑤, 𝑥 ∈ Ω, 𝑚𝑤 + V

(74)

𝑤 = 0, 𝑥 ∈ 𝜕Ω. (i) Since 𝑢 > 0 in Ω, we get from (72) that 𝑎 = 𝜆 1 (𝑢 + (𝛽 + 1) V) > 𝜆 1 .

(75)

By Proposition 2 and comparison property of principle eigenvalue, this contradicts 𝑎 < 𝜆 1 . Then (73) and (74) have V = 0 and 𝑤 = 0. (ii) Since V > 0 in Ω, we get from (73) that −𝑏 = 𝜆 1 (

𝑐𝑤 − 𝛽𝑢) > 𝜆 1 (−𝛽𝑢𝑎 ) , 𝑚𝑤 + V

(76)

by Proposition 2 and comparison property of principle eigenvalue, we get a contradiction −𝑏 ≤ 𝜆 1 (−𝛽𝑢𝑎 ).

Competing Interests The authors declare that they have no competing interests.

Acknowledgments Thanks are due to reviewers’ comments and suggestions for improving this paper. This work is supported by Grants 61563033, 11563005, and 71363043 from the National Natural Science Foundation of China and 20161BAB201010, 20151BAB212011, and 20151BAB201021 from the Natural Science Foundation of Jiangxi Province.

References [1] J. D. Murray, Mathematical Biology. I. An Introduction, vol. 17 of Interdisciplinary Applied Mathematics, Springer, New York, NY, USA, 2002. [2] Y. Li and D. Xiao, “Bifurcations of a predator-prey system of Holling and Leslie types,” Chaos, Solitons & Fractals, vol. 34, no. 2, pp. 606–620, 2007. [3] M. Liao, X. Tang, and C. Xu, “Stability and instability analysis for a ratio-dependent predator-prey system with diffusion effect,” Nonlinear Analysis: Real World Applications, vol. 12, no. 3, pp. 1616–1626, 2011. [4] R. M. Anderson and R. M. May, Infectious Diseases of Humans, Volume 1, Oxford University Press, Oxford, UK, 1991. [5] L. Rass and J. Radcliffe, Spatial Deterministic Epidemics, vol. 102, American Mathematical Society, 2003. [6] K. P. Hadeler and H. I. Freedman, “Predator-prey populations with parasitic infection,” Journal of Mathematical Biology, vol. 27, no. 6, pp. 609–631, 1989. [7] Y. Xiao and L. Chen, “Modeling and analysis of a predator-prey model with disease in the prey,” Mathematical Biosciences, vol. 171, no. 1, pp. 59–82, 2001. [8] Y. Xiao and L. Chen, “A ratio-dependent predator-prey model with disease in the prey,” Applied Mathematics and Computation, vol. 131, no. 2-3, pp. 397–414, 2002. [9] H. Yin, X. Xiao, and X. Wen, “Turing patterns in a predatorprey system with self-diffusion,” Abstract and Applied Analysis, vol. 2013, Article ID 891738, 10 pages, 2013. [10] H. Yin, J. Zhou, X. Xiao, and X. Wen, “Analysis of a diffusive Leslie-Gower predator-prey model with nonmonotonic functional response,” Chaos, Solitons & Fractals, vol. 65, pp. 51–61, 2014. [11] H. Yin, X. Xiao, X. Wen, and K. Liu, “Pattern analysis of a modified Leslie-Gower predator-prey model with Crowley-Martin functional response and diffusion,” Computers & Mathematics with Applications, vol. 67, no. 8, pp. 1607–1621, 2014. [12] X. Zhang, Y. Huang, and P. Weng, “Permanence and stability of a diffusive predator-prey model with disease in the prey,” Computers and Mathematics with Applications, vol. 68, no. 10, pp. 1431–1445, 2014. [13] J. Zhou and J. Shi, “The existence, bifurcation and stability of positive stationary solutions of a diffusive LeslieGower predator-prey model with Holling-type II functional responses,” Journal of Mathematical Analysis and Applications, vol. 405, no. 2, pp. 618–630, 2013. [14] J. Zhou, “Positive solutions of a diffusive predator-prey model with modified Leslie-Gower and Holling-type II schemes,”

10

[15]

[16]

[17]

[18]

[19]

[20]

[21] [22]

[23]

Discrete Dynamics in Nature and Society Journal of Mathematical Analysis and Applications, vol. 389, no. 2, pp. 1380–1393, 2012. J. Zhou and C. Mu, “Positive solutions for a three-trophic food chain model with diffusion and Beddington-Deangelis functional response,” Nonlinear Analysis: Real World Applications, vol. 12, no. 2, pp. 902–917, 2011. L. Li, “Coexistence theorems of steady states for predator-prey interacting systems,” Transactions of the American Mathematical Society, vol. 305, no. 1, pp. 143–166, 1988. E. N. Dancer, “On the indices of fixed points of mappings in cones and applications,” Journal of Mathematical Analysis and Applications, vol. 91, no. 1, pp. 131–151, 1983. M. Wang and Q. Wu, “Positive solutions of a prey-predator model with predator saturation and competition,” Journal of Mathematical Analysis and Applications, vol. 345, no. 2, pp. 708– 718, 2008. E. N. Dancer and Y. Du, “Positive solutions for a threespecies competition system with diffusion—I. General existence results,” Nonlinear Analysis: Theory, Methods & Applications, vol. 24, no. 3, pp. 337–357, 1995. A. Pazy, Semigroups of Linear Operators and Applications to Partial Differential Equations, vol. 44 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1983. J. W. Cholewa and T. Dlotko, Global Attractors in Abstract Parabolic Problems, vol. 278, Cambridge University Press, 2000. E. N. Dancer, “On positive solutions of some pairs of differential equations,” Transactions of the American Mathematical Society, vol. 284, no. 2, pp. 729–743, 1984. M. A. Krasnoselskij, Positive Solutions of Operator Equations, Noordhoff, Groningen, The Netherlands, 1964.

Advances in

Operations Research Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Advances in

Decision Sciences Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Applied Mathematics

Algebra

Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Probability and Statistics Volume 2014

The Scientific World Journal Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of

Differential Equations Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

Submit your manuscripts at http://www.hindawi.com International Journal of

Advances in

Combinatorics Hindawi Publishing Corporation http://www.hindawi.com

Mathematical Physics Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Journal of

Complex Analysis Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of Mathematics and Mathematical Sciences

Mathematical Problems in Engineering

Journal of

Mathematics Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Discrete Mathematics

Journal of

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Discrete Dynamics in Nature and Society

Journal of

Function Spaces Hindawi Publishing Corporation http://www.hindawi.com

Abstract and Applied Analysis

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

International Journal of

Journal of

Stochastic Analysis

Optimization

Hindawi Publishing Corporation http://www.hindawi.com

Hindawi Publishing Corporation http://www.hindawi.com

Volume 2014

Volume 2014