SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research Article

Global Behavior for a Diffusive Predator-Prey Model with Stage Structure and Nonlinear Density Restriction-I: The Case in

Rui Zhang1, Ling Guo2 and Shengmao Fu2*

Author Affiliations

1 Department of Mathematics, Lanzhou Jiaotong University, Lanzhou 730070, China

2 Department of Mathematics, Northwest Normal University, Lanzhou 730070, China

For all author emails, please log on.

Boundary Value Problems 2009, 2009:378763  doi:10.1155/2009/378763


The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2009/1/378763


Received:2 April 2009
Accepted:31 August 2009
Published:27 September 2009

© 2009 The Author(s)

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.

This paper deals with a Holling type III diffusive predator-prey model with stage structure and nonlinear density restriction in the space . We first consider the asymptotical stability of equilibrium points for the model of ODE type. Then, the existence and uniform boundedness of global solutions and stability of the equilibrium points for the model of weakly coupled reaction-diffusion type are discussed. Finally, the global existence and the convergence of solutions for the model of cross-diffusion type are investigated when the space dimension is less than 6.

1. Introduction

Population models with stage structure have been investigated by many researchers, and various methods and techniques have been used to study the existence and qualitative properties of solutions [19]. However, most of the discussions in these works are devoted to either systems of ODE or weakly coupled systems of reaction-diffusion equations. In this paper we investigate the global existence and convergence of solutions for a strongly coupled cross-diffusion predator-prey model with stage structure and nonlinear density restriction. Nonlinear problems of this kind are quite difficult to deal with since the usual idea to apply maximum principle arguments to get priori estimates cannot be used here [10].

Consider the following predator-prey model with stage-structure:

(11)

where , denote the density of the immature and mature population of the prey, respectively, is the density of the predator. For the prey, the immature population is nonlinear density restriction. is assumed to consume with Holling type III functional response and contributes to its growth with rate . For more details on the backgrounds of this model see references [11, 12].

Using the scaling and redenoting by , we can reduce the system (1.1) to

(12)

where

To take into account the natural tendency of each species to diffuse, we are led to the following PDE system of reaction-diffusion type:

(13)

where is a bounded domain in with smooth boundary , is the outward unit normal vector on , and . are nonnegative smooth functions on . The diffusion coefficients are positive constants. The homogeneous Neumann boundary condition indicates that system (1.3) is self-contained with zero population flux across the boundary. The knowledge for system (1.3) is limited (see [1317]).

In the recent years there has been considerable interest to investigate the global behavior for models of interacting populations with linear density restriction by taking into account the effect of self-as well as cross-diffusion [1826]. In this paper we are led to the following cross-diffusion system:

(14)

where are the diffusion rates of the three species, respectively. are referred as self-diffusion pressures, and is cross-diffusion pressure. The term self-diffusion implies the movement of individuals from a higher to a lower concentration region. Cross-diffusion expresses the population fluxes of one species due to the presence of the other species. The value of the cross-diffusion coefficient may be positive, negative, or zero. The term positive cross-diffusion coefficient denotes the movement of the species in the direction of lower concentration of another species and negative cross-diffusion coefficient denotes that one species tends to diffuse in the direction of higher concentration of another species [27]. For , problem (1.4) becomes strongly coupled with a full diffusion matrix. As far as the authors are aware, very few results are known for cross-diffusion systems with stage-structure.

The main purpose of this paper is to study the asymptotic behavior of the solutions for the reaction-diffusion system (1.3), the global existence, and the convergence of solutions for the cross-diffusion system (1.4). The paper will be organized as follows. In Section 2 a linear stability analysis of equilibrium points for the ODE system (1.2) is given. In Section 3 the uniform bound of the solution and stability of the equilibrium points to the weakly coupled system (1.3) are proved. Section 4 deals with the existence and the convergence of global solutions for the strongly coupled system (1.4).

2. Global Stability for System (1.2)

Let . If , then (1.2) has semitrivial equilibria , where . To discuss the existence of the positive equilibrium point of (1.2), we give the following assumptions:

(21)

where . Let one curve : , and the other curve : . Obviously, passes the point . Noting that attains its maximum at , thus when , . has the asymptote and passes the point . In this case, and have unique intersection , as shown in Figure 1. is the unique positive equilibrium point of (1.2), where , , . In addition, the restriction of the existence of the positive equilibrium can be removed, if .

The Jacobian matrix of the equilibrium is

(22)

The characteristic equation of () is . is a saddle for . In addition, the dimensions of the local unstable and stable manifold of are 1 and 2, respectively. is locally asymptotically stable for .

The Jacobian matrix of the equilibrium is

(23)

where . The characteristic equation of () is , where

(24)

According to Routh-Hurwitz criterion, is locally asymptotically stable for and , that is, and .

The Jacobian matrix of the equilibrium is

(25)

where

(26)

The characteristic equation of () is , where

(27)

According to Routh-Hurwitz criterion, is locally asymptotically stable for . Obviously, can be checked by (2.1).

Now we discuss the global stability of equilibrium points for (1.2).

Theorem 2.1.

(i) Assume that (2.1),

(28)

hold, then the equilibrium point of (1.2) is globally asymptotically stable.

(ii) Assume that , and hold, then the equilibrium point of (1.2) is globally asymptotically stable.

(iii) Assume that holds, then the equilibrium point of (1.2) is globally asymptotically stable.

Proof.

(i) Define the Lyapunov function

(29)

Calculating the derivative of along the positive solution of (1.2), we have

(210)

When , the minimum of and is and 0, respectively; the maximum of is are and , respectively. Thus, when (2.8) hold, According to the Lyapunov-LaSalle invariance principle [28], is globally asymptotically stable if (2.1)–(2.3) hold.

(ii) Let

(211)

Then

(212)

Noting that the maximum of is , and , we find . Therefore,

(iii) Let

(213)

then

(214)

Thus, for . This completes the proof of Theorem 2.1.

3. Global Behavior of System (1.3)

In this section we discuss the existence, uniform boundedness of global solutions, and the stability of constant equilibrium solutions for the weakly coupled reaction-diffusion system (1.3). In particular, the unstability results in Section 2 also hold for system (1.3) because solutions of (1.2) are also solutions of (1.3).

Theorem 3.1.

Let be nonnegative smooth functions on . Then system (1.3) has a unique nonnegative solution , and

(31)

on . In particular, if , then for all .

Proof.

It is easily seen that is sufficiently smooth in and possesses a mixed quasimonotone property in . In addition, and are a pair of lower-upper solutions of problem (1.3) (cf. in (3.1)). From [29, Theorem  5.3.4], we conclude that (1.3) exists a unique classical solution satisfying (3.1). According to strong maximum principle, it follows that . So the proof of the Theorem is completed.

Remark 3.2.

When (namely ), system (1.3) reduces to a system in which the immature population of the prey is linear density restriction. Similar to the proof of Theorem 3.1, we have

(32)

Now we show the local and global stability of constant equilibrium solutions for (1.3), respectively.

Theorem 3.3.

(i) Assume that (2.1) holds, then the equilibrium point of (1.3) is locally asymptotically stable.

(ii) Assume that , , and hold, then the equilibrium point of (1.3) is locally asymptotically stable.

(iii) Assume that holds, then the equilibrium point of (1.3) is locally asymptotically stable.

Proof.

Let be the eigenvalues of the operator on with Neumann boundary condition, and let be the eigenspace corresponding to in . Let

(33)

where is an orthonormal basis of , then

(34)

Let , , where

(35)

The linearization of (1.3) is at . For each , is invariant under the operator L, and is an eigenvalue of L on , if and only if is an eigenvalue of the matrix . The characteristic equation is , where

(36)

From Routh-Hurwitz criterion, we can see that three eigenvalues (denoted by , , ) all have negative real parts if and only if . Noting that , we must have . It is easy to check that if (see Section 2).

We can conclude that there exists a positive constant , such that

(37)

In fact, let , then

(38)

Since as , it follows that

(39)

Clearly, has the three roots . Let . By continuity, there exists such that the three roots of satisfy

(310)

Let , then . Let , then (3.7) holds. According to [30, Theorem  5.1.1], we have the locally asymptotically stability of .

(ii) The linearization of (1.4) is at , where , and

(311)

The characteristic equation of is , where

(312)

The three roots of all have negative real parts for and . Namely, is the locally asymptotically stable, if and .

(iii) The linearization of (1.3) is at , where , and

(313)

Similar to (i), is locally asymptotically stable, when .

Remark 3.4.

When denote . If , then (1.3) has the semitrivial equilibrium point , where . If , then (1.3) has a unique positive equilibrium point . Similar as Theorem 3.3, we have the following.

(i)If , , and (namely, , , ), then is locally asymptotically stable.

(ii)If and , then is locally asymptotically stable.

(iii)If , then is locally asymptotically stable.

Before discussing the global stability, we give an important lemma which has been proved in [31, Lemma ] or in [32, Lemma ].

Lemma 3.5.

Let be positive constants. Assume that , , and is bounded from below. If and for some positive constant , then

Theorem 3.6.

(i) Assume that (2.1),

(314)

hold, then the equilibrium point of system (1.3) is globally asymptotically stable.

(ii) Assume that , and hold, then the equilibrium point of system (1.3) is globally asymptotically stable.

(iii)Assume that and hold, then the equilibrium point of system (1.3) is globally asymptotically stable.

Proof.

Let be the unique positive solution of (1.3). By Theorem 3.1, there exists a positive constant C which is independent of and such that , for . By [33, Theorem ],

(315)

(i) Define the Lyapunov function

(316)

By Theorem 3.1, is defined well for all solutions of (1.3) with the initial functions . It is easily see that and if and only if .

Calculating the derivative of along positive solution of (1.3) by integration by parts and the Cauchy inequality, we have

(317)

It is not hard to verify that

(318)

if (3.14) hold. Applying Lemma 3.5, we can obtain

(319)

Recomputing , we find

(320)

From (3.15), we can see that is bounded in , . It follows from Lemma 3.5 and (3.15) that as . Namely,

(321)

Using the Pioncaré inequality, we have

(322)

where Noting that

(323)

according to (3.19) and (3.22), we can see

(324)

Thus, there exists as . Applying the boundness of , there exists a subsequence of , denoted still by , such that On the one hand

(325)

On the other hand

(326)

According to (3.19) to compute the limit of the previous equation and using the uniqueness of the limit, we have , and

(327)

It follows from (3.15) that there exists a subsequence of , denoted still by , and nonnegative functions , such that

(328)

Applying (3.19)–(3.27), we obtain that , and

(329)

In view of Theorem 3.3, we can conclude that is globally asymptotically stable.

(ii) Let

(330)

Then

(331)

Therefore, It follows that the equilibrium point of (1.3) is globally asymptotically stable.

(iii) Define

(332)

Then

(333)

When ,

(334)

The following proof is similar to (i).

Remark 3.7.

When , Theorem 3.6 shows the following.

(i) Assume that ,

(335)

hold, then the equilibrium point of (1.3) is globally asymptotically stable.

(ii) Assume that hold, then the equilibrium point of (1.3) is globally asymptotically stable.

(iii) Assume that and hold, then the equilibrium point of (1.3) is globally asymptotically stable.

Example 3.8.

Consider the following system:

(336)

Using the software Matlab, one can obtain , . It is easy to see that the previous system satisfies the all conditions of Theorem 3.6(i). So the positive equilibrium point (0.5637,0.5637,0.1199) of the previous system is globally asymptotically stable.

4. Global Existence and Stability of Solutions for the System (1.4)

By [3436], we have the following result.

Theorem 4.1.

If , then (1.4) has a unique nonnegative solution , where is the maximal existence time of the solution. If the solution satisfies the estimate

(41)

then . If, in addition, , then

In this section, we consider the existence and the convergence of global solutions to the system (1.4).

Theorem 4.2.

Let and the space dimension . Suppose that are nonnegative functions and satisfy zero Neumann boundary conditions. Then (1.4) has a unique nonnegative solution

In order to prove Theorem 4.2, some preparations are collected firstly.

Lemma 4.3.

Let be a solution of (1.4). Then

(42)

where .

Proof.

From the maximum principle for parabolic equations, it is not hard to verify that and is bounded.

Multiplying the second equation of (1.4) by , adding up the first equation of (1.4), and integrating the result over , we obtain

(43)

Using Young inequality and Hlder inequality, we have

(44)

where It follows from (4.3) and (4.4) that

(45)

Thus,

(46)

where depends on and coefficients of (1.4). In addition, there exists a positive constant , such that

(47)

Integrating the first equation of (1.4) over , we have

(48)

Integrating (4.8) from to , we have

(49)

According to (4.7), there exists a positive constant , such that

(410)

Multiplying the second equation of (1.4) by and integrating it over , we obtain

(411)

Integrating the previous inequation from to , we have

(412)

Lemma 4.4.

Let be a solution of (1.4), , and . Then there exists a positive constant depending on and , such that

(413)

Furthermore and

Proof.

satisfies the equation

(414)

where are functions of and so are bounded because of Lemma 4.3.

Multiply the second equation of (1.4) by and integrate it over to obtain

(415)

Then

(416)

and . From a disposal similar to the proof of Lemma  2.2 in [23], we have . Using a standard embedding result, we obtain

Lemma 4.5 (see [23, Lemmas  2.3 and 2.4]).

Let , , and let be any number which may depend on . Then there is a constant depending on , and such that

(417)

for any with for all .

To obtain -estimates of , we establish -estimates of in the following lemma.

Lemma 4.6.

Let , , then there exist positive constants and , such that

(418)

Proof.

Multiply the first equation of (1.4) by for and integrate by parts over to obtain

(419)

Integrating (4.19) from 0 to , we have

(420)

Then substitution of , into (4.20) leads to

(421)

It follows from Hlder inequality and Lemma 4.3 that

(422)

Note that , and for . From Hlder inequality, Young inequality, and Lemma 4.4, we have

(423)

Substitution of (4.22) and (4.23) into (4.21) leads to

(424)

where is arbitrary and .

Choose such that

(425)

then it follows from (4.24) that

(426)

Let

(427)

Then for

(428)

According to Lemma 4.5 and the definition of , we can see

(429)

Combining (4.26) and (4.29), we have

(430)

where . Therefore is bounded from (4.30).

From (4.29), we have . Namely, , . Combining (4.28), we have , where .

Setting in (4.20) (it is easily checked that , i.e., ), we have .

Multiplying the second equation of (1.4) by and integrating it over , we have

(431)

The result of can be obtained from an analogue of the previous proof of 's.

Lemma 4.7.

Let , then there exists a positive constant such that

(432)

Proof.

We will prove this lemma by [37, Theorem  7.1, page 181]. At first, we rewrite the first two equations of (1.4) as

(433)

where , , , is symbol. It follows from Lemma 4.6 that , .

By the third equation of (1.4), we have

(434)

It follows from Lemma 4.3 that is bounded in . Applying Theorem [37, Page 204] to (4.34), we have

(435)

Recall that satisfy (4.14) in Lemma 4.4, that is,

(436)

where is bounded. Since by (4.35), applying Theorem [37, page 341-342] to (4.36), we have

(437)

It follows from [37, Lemma , page 80] that and so . Recall from Lemma 4.6 that , so that by applying Theorem [37, Page 181] to (4.33).

Proof of Theorem 4.2.

Firstly, Theorem 4.2 can be proved in a similar way as Theorem in [21, 25] when the space dimension .

Secondly, for , applying Lemma [37, Page 80] to (4.36), we have

(438)

Since , we obtain

(439)

The first two equations can be written in the divergence form as

(440)

where . It follows from Lemmas 4.1, 4.5, and (4.39) that are bounded. Thus applying Theorem [37, Page 204] to (4.40) leads to

(441)

We rewrite the third equation of (1.4) as

(442)

where . Applying Schauder estimate [29, Theorem , page 114] to (4.42) gives

(443)

Let

(444)

then

(445)

where , . From (4.41), we have . It follows from (4.41) and (4.43) that . Applying Schauder estimate to (4.45) gives

(446)

Solving equations (4.44) for , respectively, we have

(447)

In particular, to conclude , we need to repeat the above bootstrap technique. Since is arbitrary, so the classical solution of (1.4) exists globally in time.

Now we discuss the global stability of the positive equilibrium (see Section 2) for (1.4).

Theorem 4.8.

Assume that the all conditions in Theorem 4.2, (2.1), and

(448)

hold. Let be the unique positive equilibrium point of (1.4), and let be a positive solution for (1.4). Then

(449)

provided that is large enough.

Proof.

Define the Lyapunov function

(450)

Let be a positive solution of (1.4), Then

(451)

The first integrand in the right hand of the previous inequality is positive definite if

(452)

Therefore, when the all conditions in Theorem 4.8 hold, there exists a positive constant such that

(453)

This implies that . So the proof of Theorem 4.8 is completed.

Acknowledgments

This work has been partially supported by the China National Natural Science Foundation (no. 10871160), the NSF of Gansu Province (no. 096RJZA118), the Scientific Research Fund of Gansu Provincial Education Department, and NWNU-KJCXGC-03-47 Foundation.

References

  1. Aiello, WG, Freedman, HI: A time-delay model of single-species growth with stage structure. Mathematical Biosciences. 101(2), 139–153 (1990). PubMed Abstract | Publisher Full Text OpenURL

  2. Zhang, X, Chen, L, Neumann, AU: The stage-structured predator-prey model and optimal harvesting policy. Mathematical Biosciences. 168(2), 201–210 (2000). PubMed Abstract | Publisher Full Text OpenURL

  3. Liu, S, Chen, L, Liu, Z: Extinction and permanence in nonautonomous competitive system with stage structure. Journal of Mathematical Analysis and Applications. 274(2), 667–684 (2002). Publisher Full Text OpenURL

  4. Lin, Z: Time delayed parabolic system in a two-species competitive model with stage structure. Journal of Mathematical Analysis and Applications. 315(1), 202–215 (2006). Publisher Full Text OpenURL

  5. Xu, R: A reaction-diffusion predator-prey model with stage structure and nonlocal delay. Applied Mathematics and Computation. 175(2), 984–1006 (2006). Publisher Full Text OpenURL

  6. Xu, R, Chaplain, MAJ, Davidson, FA: Global convergence of a reaction-diffusion predator-prey model with stage structure for the predator. Applied Mathematics and Computation. 176(1), 388–401 (2006). Publisher Full Text OpenURL

  7. Xu, R, Chaplain, MAJ, Davidson, FA: Global convergence of a reaction-diffusion predator-prey model with stage structure and nonlocal delays. Computers & Mathematics with Applications. 53(5), 770–788 (2007). PubMed Abstract | Publisher Full Text OpenURL

  8. Wang, M: Stability and Hopf bifurcation for a prey-predator model with prey-stage structure and diffusion. Mathematical Biosciences. 212(2), 149–160 (2008). PubMed Abstract | Publisher Full Text OpenURL

  9. Wang, Z, Wu, J: Qualitative analysis for a ratio-dependent predator-prey model with stage structure and diffusion. Nonlinear Analysis: Real World Applications. 9(5), 2270–2287 (2008). Publisher Full Text OpenURL

  10. Galiano, G, Garzón, ML, Jüngel, A: Semi-discretization in time and numerical convergence of solutions of a nonlinear cross-diffusion population model. Numerische Mathematik. 93(4), 655–673 (2003). Publisher Full Text OpenURL

  11. Chen, L: Mathematical Models and Methods in Ecology, Science Press, Beijing, China (1988)

  12. Chen, LJ, Sun, JH: The uniqueness of a limit cycle for a class of Holling models with functional responses. Acta Mathematica Sinica. 45(2), 383–388 (2002)

  13. Li, W-T, Wu, S-L: Traveling waves in a diffusive predator-prey model with Holling type-III functional response. Chaos, Solitons & Fractals. 37(2), 476–486 (2008). PubMed Abstract | Publisher Full Text OpenURL

  14. Ko, W, Ryu, K: Qualitative analysis of a predator-prey model with Holling type II functional response incorporating a prey refuge. Journal of Differential Equations. 231(2), 534–550 (2006). Publisher Full Text OpenURL

  15. Zhang, H, Georgescu, P, Chen, L: An impulsive predator-prey system with Beddington-DeAngelis functional response and time delay. International Journal of Biomathematics. 1(1), 1–17 (2008). Publisher Full Text OpenURL

  16. Fan, Y, Wang, L, Wang, M: Notes on multiple bifurcations in a delayed predator-prey model with nonmonotonic functional response. International Journal of Biomathematics. 2(2), 129–138 (2009). Publisher Full Text OpenURL

  17. Wang, F, An, Y: Existence of nontrivial solution for a nonlocal elliptic equation with nonlinear boundary condition. Boundary Value Problems. 2009, (2009)

  18. Lou, Y, Ni, W-M: Diffusion, self-diffusion and cross-diffusion. Journal of Differential Equations. 131(1), 79–131 (1996). Publisher Full Text OpenURL

  19. Lou, Y, Ni, W-M, Wu, Y: On the global existence of a cross-diffusion system. Discrete and Continuous Dynamical Systems. 4(2), 193–203 (1998)

  20. Shim, S-A: Uniform boundedness and convergence of solutions to cross-diffusion systems. Journal of Differential Equations. 185(1), 281–305 (2002). Publisher Full Text OpenURL

  21. Shim, S-A: Uniform boundedness and convergence of solutions to the systems with cross-diffusions dominated by self-diffusions. Nonlinear Analysis: Real World Applications. 4(1), 65–86 (2003). Publisher Full Text OpenURL

  22. Choi, YS, Lui, R, Yamada, Y: Existence of global solutions for the Shigesada-Kawasaki-Teramoto model with weak cross-diffusion. Discrete and Continuous Dynamical Systems. 9(5), 1193–1200 (2003)

  23. Choi, YS, Lui, R, Yamada, Y: Existence of global solutions for the Shigesada-Kawasaki-Teramoto model with strongly coupled cross-diffusion. Discrete and Continuous Dynamical Systems. 10(3), 719–730 (2004)

  24. Pang, PYH, Wang, MX: Existence of global solutions for a three-species predator-prey model with cross-diffusion. Mathematische Nachrichten. 281(4), 555–560 (2008). Publisher Full Text OpenURL

  25. Fu, S, Wen, Z, Cui, S: Uniform boundedness and stability of global solutions in a strongly coupled three-species cooperating model. Nonlinear Analysis: Real World Applications. 9(2), 272–289 (2008). Publisher Full Text OpenURL

  26. Yang, F, Fu, S: Global solutions for a tritrophic food chain model with diffusion. The Rocky Mountain Journal of Mathematics. 38(5), 1785–1812 (2008). Publisher Full Text OpenURL

  27. Dubey, B, Das, B, Hussain, J: A predator-prey interaction model with self and cross-diffusion. Ecological Modelling. 141(1–3), 67–76 (2001)

  28. Hale, JK: Ordinary Differential Equations, Krieger, Malabar, Fla, USA (1980)

  29. Ye, Q, Li, Z: Introduction to Reaction-Diffusion Equations, Science Press, Beijing, China (1999)

  30. Henry, D: Geometric Theory of Semilinear Parabolic Equations, Lecture Notes in Mathematics, Springer, Berlin, Germany (1993)

  31. Lin, Z, Pedersen, M: Stability in a diffusive food-chain model with Michaelis-Menten functional response. Nonlinear Analysis: Theory, Methods & Applications. 57(3), 421–433 (2004). PubMed Abstract | Publisher Full Text OpenURL

  32. Wang, M: Nonliear Parabolic Equation of Parabolic Type, Science Press, Beijing, China (1993)

  33. Brown, KJ, Dunne, PC, Gardner, RA: A semilinear parabolic system arising in the theory of superconductivity. Journal of Differential Equations. 40(2), 232–252 (1981). Publisher Full Text OpenURL

  34. Amann, H: Dynamic theory of quasilinear parabolic equations. I. Abstract evolution equations. Nonlinear Analysis: Theory, Methods & Applications. 12(9), 895–919 (1988). PubMed Abstract | Publisher Full Text OpenURL

  35. Amann, H: Dynamic theory of quasilinear parabolic equations. II. Reaction-diffusion systems. Differential and Integral Equations. 3(1), 13–75 (1990)

  36. Amann, H: Dynamic theory of quasilinear parabolic systems. III. Global existence. Mathematische Zeitschrift. 202(2), 219–250 (1989). Publisher Full Text OpenURL

  37. Ladyženskaja, OA, Solonnikov, VA, Ural'ceva, NN: Linear and Quasilinear Equations of Parabolic Type, Translations of Mathematical Monographs,p. xi+648. American Mathematical Society, Providence, RI, USA (1967)