The Scientific World Journal

Volume 2014 (2014), Article ID 721403, 10 pages

http://dx.doi.org/10.1155/2014/721403

Research Article

Dynamics of a Diffusive Predator-Prey Model with General Nonlinear Functional Response

Wensheng Yang

School of Mathematics and Computer Science, Fujian Normal University, Fuzhou, Fujian 350007, China

Received 7 August 2013; Accepted 13 November 2013; Published 3 February 2014

Academic Editors: M. Bellassoued, J. Fernández, and J.-L. Liu

Copyright © 2014 Wensheng Yang. 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 diffusive predator-prey model with nonconstant death rate and general nonlinear functional response. Firstly, stability analysis of the equilibrium for reduced ODE system is discussed. Secondly, sufficient and necessary conditions which guarantee the predator and the prey species to be permanent are obtained. Furthermore, sufficient conditions for the global asymptotical stability of the unique positive equilibrium of the system are derived by using the method of Lyapunov function. Finally, we show that there are no nontrivial steady state solutions for certain parameter configuration.

1. Introduction

One of the dominant themes in both ecology and mathematical ecology is the dynamic relationship between predators and their prey due to its universal existence and importance in population dynamics. The investigations on predator-prey models are developed during these thirty years, and more realistic models are derived in view of laboratory experiments and observations (see [1–6] and the references therein).

In [7], Duque and Lizana considered the following reaction-diffusion predator-prey model with Holling type II functional response:

where and represent the population density of prey and predator at and at time , respectively. The parameter is the specific growth rate of the prey in the absence of predation and without environment limitation; in the absence of predators, the prey population grows logistically to a carrying capacity ; the functional response of the predator is a Holling type II function ; , , and are satiation coefficients or conversion rates. The specific mortality of predators in the absence of prey

depends on the quantity of predator, is the mortality at low density, and is the maximal mortality with the natural assumption . are constants, , while denotes the Laplace operator in with bounded and connected. The advantage of the present model over the more often used models is that here the predator mortality is neither a constant nor an unbounded function, and still it is increasing with quantity.

In population dynamics, a functional response of the predator to the prey density refers to the change in the density of prey attached per unit time per predator as the prey density changes. The simplest functional response is a Lotka-Volterra function which is described as

which is also called a Holling type I function. Michaelis and Menten proposed the response function

in studying enzymatic reactions. It is now referred to as a Michaelis-Menten function or a Holling type II function. However, in Holling type I function and Holling type II function is monotonic in the first quadrant. It implies that, as the prey population increases, the consumption rate of prey per predator increases. But some experiments and observations indicate that a nonmonotonic response occurs at this level: when the nutrient concentration reaches a high level, an inhibitory effect on the specific growth rate may occur. To model such an inhibitory effect, Andrews [8] suggested a function

called the Monod-Haldane function and also called a Holling type IV function. Sokol and Howell [9] proposed a simplified Holling type IV function of the form

Motivated by the above question a functional response may be monotonic or nonmonotonic; we consider the following diffusive predator-prey model with general nonlinear functional response:

where is a general functional response, and the function is assumed to satisfy the following assumptions:

(i) is of class , ;(ii) is decreasing on and there exists a such that ;(iii)there exists a such that for ;(iv)there exists a such that for ;(v)there exists a such that is increasing for any , where .

Note that hypotheses (i)–(v) are satisfied if function represents Holling type II functional response or Holling type IV functional response; that is, or .

If , system (7) reduces to the system (1). Duque and Lizana [7] obtained necessary and sufficient conditions under which the system (1) is permanent. To the best of the author knowledge, for the predator-prey system (7) with general functional response, whether one could obtain the sufficient and necessary conditions which insure the permanence of the system (7) or not is still an open problem.

The aim of this paper is, by further developing the analysis technique of Duque and Lizana [7], to obtain sufficient and necessary conditions which ensure the permanence of the system (7) and describe the global dynamic of the system (7). More concretely, we firstly discuss stability analysis of the equilibrium for reduced ODE system. Secondly, we obtain necessary and sufficient conditions under which the system is dissipative and permanent. Furthermore, we study the global stability of the nontrivial equilibrium when it is unique. Finally, we show that system (7) has no nontrivial steady state solutions for certain parameter configuration.

2. The Model without Diffusion

For the reaction-diffusion predator-prey system (7), the reduced system is an ordinary differential equation of the form

For simplicity, let us rewrite (8) as

where

The equilibria of (9) consist of two trivial critical points and on the boundary of and a set of nontrivial critical points obtained as the intersection of the curves

It is well known that persistence implies the existence of a positive equilibrium; see, for instance, [10]. From Theorem 8 in Section 3, one can easily know that system (9) has a positive equilibrium point , if , where

The Jacobian matrix of system (9) at is

and then it is clear that the Jacobian matrix has two eigenvalues , . Hence, the equilibrium is a saddle point and the stable and unstable manifolds lie on the -axis and -axis, respectively.

The Jacobian matrix of system (9) at is

and then it is clear that the Jacobian matrix has two eigenvalues , . Hence the equilibrium is a saddle point, if , and it is locally asymptotically stable when .

The Jacobian matrix of system (9) at is

and the eigenvalues of the matrix are the roots of the polynomial

where

The following result holds.

Proposition 1. If , then the equilibrium is locally asymptotically stable with respect to system (9).

Proof. Taking into account that , the assumptions (i)–(v), and (2), it follows that

which in turn implies that . Moreover, since and , it follows that , which completes the proof.

3. Extinction of the Predator and Permanence

In this section, we will show that any nonnegative solution of (7) lies in a certain bounded region as for all . We first show a well-known conclusion on the Logistic equation.

Lemma 2 (see [6]). Assume that is defined by

then .

Theorem 3. All the solutions of (7) are nonnegative and defined for all . Furthermore, the nonnegative solution of (7) satisfies

Proof. The nonnegativity of the solutions of (7) is clear since the initial value is nonnegative. We only consider the latter of the theorem.

From the first equation of system (7), we have

then from comparison principle of parabolic equations and Lemma 2, for an arbitrary , there exists such that, for any ,

This implies

Now, let us prove the boundedness of . Taking into account the second equation of system (7) and the fact that , we get

Let us show that there cannot exist a such that

where

Indeed, assume that (25) is true. Choosing and taking into account that

we get that

where . Then, from (28) and the first equation of system (7), we obtain

This immediately implies that

where is the solution of the equation

Since

Having in mind that (30) and (33) hold, we get that

From (34), we can easily obtain that there exists such that

Now, taking into account the second equation of system (7) and (35), we get

By using the comparison principle, we obtain that , which is a contradiction.

Let us define the function , and . We have to consider two possibilities. If there is some such that for any , we say that the zeroes of are bounded. If this is not true we say that the zeroes are unbounded; in this case, for a sequence of tending to as .

If the zeroes of are bounded, then for any .

If the zeroes are unbounded, the -axis is divided by the zeroes into a sequence of intervals . In each, such interval is of constant sign. Let us assume that for . Since on the intervals , arguing as in the proof that (25) is not possible, we obtain that

for and which is large enough, and . So, by the comparison principle, for , where is the solution of the equation , with a suitable initial condition. Thus,

and which is large enough. This contradiction completes the proof of our claim.

Remark 4. An immediate consequence of the proof of the former result is that for a given ,

is an absorbing set for system (7). Hence, system (7) is dissipative.

Now, we will summarize some facts contained in Hale and Waltman studies, see [11], about the permanence for abstract dynamical systems.

Suppose that is a complete metric space with for an open set , where is the boundary of the set . We will typically choose to be the positive cone in an ordered Banach space. A flow or semiflow on under which and are forward invariant is said to be permanent if it is dissipative and if there is a number such that any trajectory starting in will be at least a distance from for all sufficiently large . To state a theorem implying permanence, we need a few definitions. An invariant set for the flow or semiflow is said to be isolated if it has a neighborhood such that is the maximal invariant subset of . Let denote the union of the sets over . This differs from the standard definition of the -limit set of a set but it is more convenient for our purposes; see [10] for a discussion. The set is said to be isolated if it has a covering of pairwise disjoint, both sets which are isolated and invariant with respect to the flow or semiflow both on and on . The covering is then called an isolated covering. Suppose and are isolated invariant sets (not necessarily distinct). The set is said to be chained to (denoted by ) if there exists with . As usual, and denote the unstable and stable manifolds, respectively. A finite sequence of isolated invariant sets is a chain if . This is possible for if . The chain is called a cycle if . The set is said to be acyclic if there exists an isolated covering such that no subset of is a cycle. We now state a lemma that can be used to establish permanence.

Lemma 5 (see [11]). Suppose that is a complete metric space with , where is open. Suppose that a semiflow on leaves both and forward invariant, maps bounded sets in to precompact sets for , and is dissipative. If in addition

(1) is isolated and acyclic,(2) for all , where is the isolated covering used in the definition of acyclicity of ,

then the semiflow is permanent; that is, there exists such that any trajectory with initial data in will be bounded away from by a distance greater than for which is sufficiently large.

Now, we show that system (7) is permanent. From the viewpoint of biology, this implies that the two species of prey and predator will always coexist at any time and any location of the inhabit domain, no matter what their diffusion coefficients are. First, we identify a less interesting situation where the predator cannot survive.

Proposition 6. If , then is globally asymptotically stable with respect to nonnegative initial functions.

Proof. Let us assume that , we can easily obtain that, for small enough ,

From Theorem 3, it follows that, for the above , there exists a , such that , for all .

Therefore, from the second equation of system (7), we obtain that

for all . Invoking again the comparison principle, we get that for any , where is the solution of the equation , with , where

This implies that as . Henceforth, we may assume that, under the hypothesis of our claim, for any , , for all . Now, from the first equation of system (7), we obtain

By the comparison principle, we get that for any , where is the solution of the equation . Finally, our assertion follows from the arbitrariness of and the fact that the solutions with positive initial data of the equation tend exponentially to .

From the above discussion and Theorem 3, we obtain that is globally asymptotically stable.

Remark 7. If , system (7) reduces to system (1). So, it is shown that our result supplements and covers the Proposition 4.1 of Duque and Lizana's paper [7].

Theorem 8. System (7) is permanent if and only if

Proof. Let us set , , and , where

Henceforth, considering also an initial condition, system (7) can be rewritten as

Denote , where . Assume that operator is compact for , which is defined by , where is the classical solution of system (7). Moreover, from Theorem 3, it follows that the semiflow is pointwise dissipative.

A direct application of the maximum principle shows that is positively invariant on , which in turn implies that is positively invariant on .

Now, let us show that . First, we consider

There is no loss of generality assuming that for all . Indeed, if for all with and applying the strong maximum principle, we obtain that for all and . One could change the origin of time to a positive time and choose .

Let and be the solutions of the equation such that and , respectively. From the comparison principle, we obtain

Taking into account that attracts any positive solution, we conclude that as .

Let us now consider

Since , we get that

by the comparison principle, we obtain that as . Henceforth, we conclude that and it is isolated and acyclic. By choosing and , then is the covering required by Lemma 2.

Taking into account that is positively invariant and the fact that the stable and unstable manifolds of the equilibrium lie on , we obtain that . In order to show that , let us linearize (7) about , obtaining the following system:

From (44) and by using comparison techniques, we obtain the unstable manifold to point into the region . Moreover, from the above performed computations we know that the stable manifold corresponding to the equilibrium lies on . This implies that .

Since all hypotheses of Lemma 5 are fulfilled, we may conclude that system (7) is permanent.

Permanence implies that condition (44) is an immediate consequence of Theorem 3. This concludes the proof of our claim.

Remark 9. From Theorem 8, we obtain that system (7) with nonmonotonic functional response is also permanent if and only if (44) holds. That is, there exist positive constants , , , and such that , as time converges to infinity. If , system (7) reduces to system (1). So, it is shown that our result supplements and covers Theorem 4.1 of Duque and Lizana's paper [7].

Remark 10. We can easily obtain that in case the condition (44) is not satisfied, system (7) has no nontrivial equilibria. This fact is not surprising, it is well-known that persistence implies the existence of a positive equilibrium; see, for instance [10].

Remark 11. From Theorem 8, it is shown that diffusion has no influence on the permanence of system (7).

4. Global Stability of the Nontrivial Equilibrium

Hereafter in this paper, we restrict our attention to the case , in this case system (9) has a nontrivial equilibrium point. It is obvious that the equilibria of system (9) are solutions of (7). We will focus our attention on the nontrivial equilibrium of the system (9). Concretely, in this section, we will analyze the local and the global stability of nontrivial equilibrium of (7).

Let us first briefly study the local stability of the homogeneous stationary solution of (7) via the linearized stability analysis (see [12, pp. 68–70]). Setting and recalling that

the linearized system of the reaction-diffusion equation (7) about is given by

The trivial solution is asymptotically stable if and only if every solution of (53) decays to zero as . Following [12, pp. 68–70], we obtain that the trivial solution of (53) is asymptotically stable if and only if each solution of the system

tends to zero as for each . denotes the th eigenvalue corresponding to the th eigenfunction of the operator on with no-flux boundary conditions. Let us recall that the eigenvalues satisfy the following inequality:

where

See, for instance, [13, pp. 205–208].

Henceforth, the asymptotic stability of is equivalent to that each matrix has two eigenvalues with negative real parts for all . And the eigenvalues of the matrix are the roots of the polynomial

where

The following result holds.

Proposition 12. If , then the equilibrium is locally asymptotically stable with respect to system (7).

Proof. Taking into account that , the assumptions (i)–(v), and (2), it follows that

which in turn implies that . Moreover, since and , it follows that , which completes the proof.

Remark 13. From Propositions 1 and 12, we can easily obtain that diffusion has no influence on the stability of the equilibrium .

Now, we will focus on the global stability of . Biologically, our statement of the global stability of means that no matter the two species diffuse, they will be spatially homogeneously distributed as time converges to infinity.

Theorem 14. Assume that

hold; then is globally asymptotically stable.

Proof. Let us define the following function:

It is easy to show that for all . Differentiating along the solutions of the system (7), we obtain

From (60), we can easily obtain that

Thus, satisfies Lyapunov's asymptotic stability theorem, and the positive equilibrium of system (7) is globally stable. This completes the proof.

5. Nonexistence of Nonconstant Positive Steady States

The main purpose of this section is to prove that, under certain restriction on the coefficients, (7) has no nontrivial positive steady state solutions. To accomplish this goal, first, we will obtain a priori bounds for the positive steady states of system (7). The corresponding steady state problem of (7) is the elliptic system

Now let us state the following result that we will need in what follows.

Lemma 15 (see [14]). Suppose that ; then the following results hold.

(1)If satisfies for , for , and , then .(2)If satisfies for , for , and , then .

The following result gives an upper estimate for classical solutions of (64), where classical solution means solution in . Let us point out that the same upper bounds that we obtained for the dissipativity are valid for the elliptic problem. But, by using Lemma 15, we will obtain the sharpest bounds.

Proposition 16. Let such that ; then for any positive classical solution of system (64), the following estimates hold:

Proof. Let us set . From the first equation of system (64), we obtain that

Let be such that . Then by Lemma 15 and positiveness of , we get

Hence, . This in turn implies that .

Analogously, by setting , we get from the second equation of system (64) that

Let be such that . Then, by Lemma 15 and positiveness of , we get

That is,

Hence,

So, . This in turn implies that

Theorem 17. Let be the smallest positive eigenvalue of the operator on with the homogeneous Neumann boundary condition. If

where and , then (64) has no nonconstant positive classical solution.

Proof. Assume that is a nonconstant positive classical solution of (64). By (65), there exists a positive constant such that . Let and . Multiplying the first equation of (64) by and integrating over by parts, we have

Similarly, we get

Adding (74) and (75), we obtain

By using the Poincare inequality (Theorem 11.11, page 112, in [15]), it follows that

Taking into account (73), we obtain

From (77) and (78), we get that and , which is a contradiction. So, system (64) has no nonconstant positive classical solution. This completes the proof.

6. Conclusion and Discussion

In this paper, we study a diffusive predator-prey model with nonconstant death rate and general nonlinear functional response. Firstly, stability analysis of the equilibrium for reduced ODE system is discussed. Secondly, sufficient and necessary conditions which guarantee the predator and the prey species to be permanent are obtained. Furthermore, sufficient conditions for the global asymptotical stability of the unique positive equilibrium of the system are derived by using the method of Lyapunov function. Finally, we show that there are no nontrivial steady state solutions for certain parameter configuration.

Note that hypotheses (i)–(v) in Section 1 are satisfied if function represents Holling type II functional response or Holling type IV functional response; that is, or . Hence, we can obtain the sufficient and necessary conditions which guaranting the predator and the prey species of the system with monotonic or nonmonotonic functional response to be permanent. If , system (7) reduces to system (1) that was studied by Duque and Lizana in [7]. It is shown that our result supplements and covers the main results of Duque and Lizana's paper [7].

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the Natural Science Foundation of China (11102041, 11201072), the Natural Science Foundation of Fujian Province (2011J01002, 2012J01002), and the Foundation of Fujian Education Bureau (JB12030).

References

S. Chen and J. Shi, “Global stability in a diffusive Holling-Tanner predatorprey model,” Applied Mathematics Letters, vol. 25, no. 3, pp. 614–618, 2012.

W. Ko and K. Ryu, “Qualitative analysis of a predator-prey model with Holling type II functional response incorporating a prey refuge,” Journal of Differential Equations, vol. 231, no. 2, pp. 534–550, 2006.

C. V. Pao, “On nonlinear reaction-diffusion systems,” Journal of Mathematical Analysis and Applications, vol. 87, no. 1, pp. 165–198, 1982.

E. Sáez and E. González-Olivares, “Dynamics of a predator-prey model,” SIAM Journal on Applied Mathematics, vol. 59, no. 5, pp. 1867–1878, 1999.

J. T. Tanner, “The stability and the intrinsic growth rates of prey and predator populations,” Ecology, vol. 56, pp. 855–867, 1975.

Q. X. Ye and Z. Y. Li, Introduction to Reaction-Diffusion Equations, Science Press, Beijing, China, 1990.

C. Duque and M. Lizana, “On the dynamics of a predatorprey model with nonconstant death rate and diffusion,” Nonlinear Analysis: Real World Applications, vol. 12, no. 4, pp. 2198–2210, 2011.

J. F. Andrews, “A mathematical model for the continuous culture of microorganisms utilizing inhibitory substrates,” Biotechnology and Bioengineering, vol. 10, no. 6, pp. 707–723, 1968.

W. Sokol and J. A. Howell, “Kinetics of phenol oxidation by washed cells,” Biotechnology and Bioengineering, vol. 23, pp. 2039–2049, 1980.

V. Hutson and K. Schmitt, “Permanence and the dynamics of biological systems,” Mathematical Biosciences, vol. 111, no. 1, pp. 1–17, 1992.

J. K. Hale and P. Waltman, “Persistence in infinite-dimensional systems,” SIAM Journal on Mathematical Analysis, vol. 20, no. 2, pp. 388–395, 1989.

P. Grindrod, Patterns and Waves. The Theory and Applications of Reaction-Difusion Equations, Clarendon Press, Oxford, UK, 1991.

G. B. Folland, Introduction to Partial Differential Equations, Princeton University Press, 1976.

Y. Lou and W.-M. Ni, “Diffusion, self-diffusion and cross-diffusion,” Journal of Differential Equations, vol. 131, no. 1, pp. 79–131, 1996.

J. Smoller, Shock Waves and Reaction-Diffusion Equations, Springer, New York, NY, USA, 2nd edition, 1994.