Table of Contents Author Guidelines Submit a Manuscript
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

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

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.

Abstract

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 [13]. 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 [68] 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 [911] and references wherein. In this paper, we are interested in the following model having Holling-type II functional response with disease in the prey:

There are two species in the model: prey and predator . In (1a), is susceptible prey, and 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 , and is called the transmission coefficient. The disease is not genetically inherited. The infected populations do not recover or become immune. Infected 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 Michaelis-Menten functional response function; that is, Infected prey has a death rate and is a death rate of predators. Prey and predators are inhabited in a bounded region 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 time-dependent 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:which imply that predator and prey die out on the boundary and initial valuesBesides, we first attempt to investigate existence and nonexistence of inhomogeneous positive solutions of (1a), (1b), and (1c) in steady state as follows:

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.

2. Preliminary

In this section, some notations and basic well-known results are introduced. For any , the linear eigenvalue problem [13, 14],has an infinite sequence of eigenvalues. Let be the principal eigenvalue of (9) which is simple and its eigenfunction does not change sign in . We denote by for simplicity and the eigenfunction corresponding to is denoted by . with being not identically equal to implies .

Now, consider the nonlinear boundary value problemIf 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 For , we define We always assume that . Let be a compact linear operator on . We say that has property on if there exists and such that . Let be a compact operator with a fixed point and is Fréchet differentiable at Let be the Fréchet derivative of at Then maps into itself. We denote by the fixed point index of at relative to

Proposition 1 (see [1618]). Assume that is invertible on . Then (i)if has property on , then ;(ii)if does not have property on , then , where is the sum of multiplicities of all eigenvalues of which is greater than

For a linear operator , we denote by the spectral radius of .

Proposition 2 (see [16, 18]). Let and be a positive constant such that on . Let be the first eigenvalue of the problemWe have the following conclusions:(a),(b),(c)

Now, we introduce the degree calculations, which was introduced by Dancer and Du in [19]. Assume that and are ordered Banach spaces with positive cones and , respectively. Let and Then is an ordered Banach space with positive cone Let be an open set in containing and be completely continuous operators, Denote by a general element in with and Let be defined by Define

Theorem 3 (see [15]). Suppose that is relatively open and bounded and Suppose that extends to a continuously differentiable mapping of a neighborhood of into , is dense in , and Then the following are true: (i) for small, if, for any , the spectral radius and is not an eigenvalue of corresponding to a positive eigenvector;(ii) for small, if, for any , the spectral radius

3. Existence of a Global Positive Solution of (1a), (1b), and (1c)

For convenience, we introduce some notations. For , we denote , , , and and extend these notations to real-valued functions. If is partially ordered vector space, we denote its positive cone by . For , , let . Assume and its norm is defined as .

Hence, system (1a), (1b), and (1c) can be rewritten as an abstract differential equation:where and ,with

Next, the existence of the local solution for (9) will be proved. First, we give out the following lemma.

Lemma 4. For every , Cauchy problem (9) has a unique maximal local solutionwhich satisfies the following Duhamel formula for :where is a operator and the closure of in . Moreover, if then .

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 in [21], we complete this proof.

Lemma 5. For initial-boundary problem (1a), (1b), and (1c), the components of its solution , , are nonnegative.

Proof. To prove for and , we consider the following system:with boundary and initial conditions After multiplying (14a), (14b), and (14c) by , , and , respectively, and integrating by parts on the domain , we obtain Hence, for , Consequently, Now, we can know that is a solution of (1a), (1b), and (1c), and according to Lemma 4, we finally get that

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.

Theorem 6. Let be the solution of (1a), (1b), and (1c) with initial value (1c); then for , where , are defined by

Proof. In view of Lemma 5, we only prove the upper bound of (1a), (1b), and (1c) for . From (1a), we haveWe get , by maximum principle, where
Let . From (1b) we have Hence, by maximum principle we have that By using the same method, we consider ; then from (1c) we have that Thus, the maximum principle shows that

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 , which implies that the solution of (1a), (1b), and (1c) is a global solution.

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

Theorem 8. System (1a), (1b), and (1c) has a unique, nonnegative, and bounded solution , such that

4. Extinction of (1a), (1b), and (1c)

Theorem 9. Let be a positive solution of (1a), (1b), and (1c), then we have the following: (i)if and then as ;(ii)if and then as

Proof. (i) First, we observe that solution of (1a), (1b), and (1c) satisfiesLet be a sufficiently small positive constant such that Since , we see that uniformly as by using comparison principle for elliptic problems. Then there exists such that for all Therefore, we havewhich concludes that uniformly as . Similarly, there exists such that for all . Further, since is monotone increasing with respect to the variable , we havewhich concludes that as
(ii) Let be a sufficiently small positive constant such that Since , we get from (28) and comparison theorem thatThen, there exists such that for all Thus, we haveSince satisfies (a) and , we conclude that uniformly as Similarly, we can see that uniformly as
Besides, there exists such that for all Hence, we haveSince satisfies (b), we get thatThus, 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 of (2) has a priori bounds:where

Proof. From the first equation of (2), we haveand we get by maximum principle. Let with . From (1b), we have Thus, by maximum principle, we get We set with , and then we have By using maximum principle again, we obtain Therefore, we complete this proof.

Now, we introduce some notations:, where ,, where ,, according to Theorem 10.Taking we can define a positive and compact operator byIt 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 of (1a), (1b), and (1c) relative to

Lemma 11. Assume that , then one has(i),(ii)

Proof. (i) It is easy to see that has no fixed point on , so is well defined. For , we define a positive and compact operator by and then For each , a fixed point is a solution of the following problem:Form Theorem 10 has no fixed point and is well defined. Hence, Note that (46) has only the trivial solution when . Set Let for some Obviously,
Hence, has no nontrivial kernel in . Since , it is easy to see that 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 . Let ; then Assume that for some , and thenSince and , from (51) and (52), we know that If with , then from (50), which contradicts with the assumption. Thus So has nontrivial kernel in
For , by Proposition 2, we see that , and is the principle eigenvalue of with a corresponding eigenfunction . Since , we see that Set , and then Thus, has the property By Proposition 1 we have

Now, we give out the index of the semitrivial solution of (2).

Lemma 12. Assume that , and then one has(i)if , then ;(ii)if , then

Proof. Let ; then Let for some ThenSince , and . So we get that Hence, from the first equation of (54) we haveIf , then by Proposition 2. On the other hand, , and we get a contradiction. Hence, Then has no nontrivial kernel in .
(i) Since , by Proposition 1, we have which is an eigenvalue of with a corresponding eigenfunction Since , we have Set , and thenHence, has the property on resulting in according to Proposition 1.
(ii) First, we prove that does not have the property on
Since , we have On the contrary, we suppose that has the property on Then there exist and such that SoThen is an eigenvalue of the operator , which contradicts So does not have the property on By Proposition 2, we have , where is the sum of the multiplicities of all eigenvalues of . Next, we will prove that
Assume that is an eigenvalue of with corresponding eigenfunction , and then ; that is,Since and , it follows from the third equation (58) that
We suppose that , and it follows from the second equation of (58) and Proposition 2 that which contradicts , and so
If , it follows from the first equation of (58) that Thus, we get the contradiction and this implies that Hence, , which implies that has no eigenvalue being greater than ; that is, , and We complete the proof.

Next, we consider the other semitrivial solutions of (2). Let us consider the following three subsystems:It is obvious that (63) has trivial solution and that (62) only has a trivial solution and a semitrivial solution if and only if 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 with if and only if and Moreover, if and , then the positive solution satisfies and , where is the unique positive solution of the following equations:

Let be the set of all positive solutions of (61), denoted by , and It is easy to see that (61) has semitrivial solution from Theorem 13. Next, we give the degree of in relative to

Lemma 14. Assume , and then one has the following:(i)if , for any , then ;(ii)if , for any , then .

Proof. Let Define Then and we are in the framework to use Theorem 3. We choose a neighborhood of such that
Now, with if and only if By the results in [22], we have For any , a direct calculation shows that Notice that, by our choice of , so it follows from the maximum principle that is -positive in the sense of [23] with Hence, and is the only eigenvalue corresponding to a positive eigenvector. By a simple calculation, we can easily show that if and only if and if and only if . Therefore, by Theorem 3 and the discussion above, we have Again since the degree does not depend on the particular choice of and , means and . Thus, we complete this proof.

From the discussion above, we can obtain the existence of positive solutions of (2).

Theorem 15. Assume Equation (2) has at least one positive solution if and then (2) has at least one positive solution.

Proof. Since , we obtain and from Lemma 11. Hence, we haveSince , we have from Lemma 12. Moreover, Theorem 13 implies Let Since and by Theorem 13, we have and thus from Lemma 14. Therefore, 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);(ii) and

Proof. Assume that (2) has positive solution , and then satisfy the following three equations, respectively:(i) Since in , we get from (72) that By Proposition 2 and comparison property of principle eigenvalue, this contradicts Then (73) and (74) have and
(ii) Since in , we get from (73) that by Proposition 2 and comparison property of principle eigenvalue, we get a contradiction

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. View at MathSciNet
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  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. View at Publisher · View at Google Scholar · View at MathSciNet
  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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  9. H. Yin, X. Xiao, and X. Wen, “Turing patterns in a predator-prey system with self-diffusion,” Abstract and Applied Analysis, vol. 2013, Article ID 891738, 10 pages, 2013. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  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. View at Publisher · View at Google Scholar · View at Scopus
  13. J. Zhou and J. Shi, “The existence, bifurcation and stability of positive stationary solutions of a diffusive Leslie-Gower predator-prey model with Holling-type II functional responses,” Journal of Mathematical Analysis and Applications, vol. 405, no. 2, pp. 618–630, 2013. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  14. J. Zhou, “Positive solutions of a diffusive predator-prey model with modified Leslie-Gower and Holling-type II schemes,” Journal of Mathematical Analysis and Applications, vol. 389, no. 2, pp. 1380–1393, 2012. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  15. 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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  16. 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. View at Publisher · View at Google Scholar · View at MathSciNet
  17. 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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  18. 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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  19. E. N. Dancer and Y. Du, “Positive solutions for a three-species competition system with diffusion—I. General existence results,” Nonlinear Analysis: Theory, Methods & Applications, vol. 24, no. 3, pp. 337–357, 1995. View at Publisher · View at Google Scholar · View at Scopus
  20. A. Pazy, Semigroups of Linear Operators and Applications to Partial Differential Equations, vol. 44 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1983. View at Publisher · View at Google Scholar · View at MathSciNet
  21. J. W. Cholewa and T. Dlotko, Global Attractors in Abstract Parabolic Problems, vol. 278, Cambridge University Press, 2000. View at Publisher · View at Google Scholar · View at MathSciNet
  22. 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. View at Publisher · View at Google Scholar · View at MathSciNet
  23. M. A. Krasnoselskij, Positive Solutions of Operator Equations, Noordhoff, Groningen, The Netherlands, 1964.