- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Abstract and Applied Analysis

Volume 2013 (2013), Article ID 921401, 7 pages

http://dx.doi.org/10.1155/2013/921401

## The Existence of Positive Nonconstant Steady States in a Reaction: Diffusion Epidemic Model

^{1}College of Mathematics and Information Science, Wenzhou University, Wenzhou 325035, China^{2}Department of Mathematics, Hubei Minzu University, Enshi, Hubei 445000, China

Received 19 November 2013; Accepted 1 December 2013

Academic Editor: Kaifa Wang

Copyright © 2013 Yuan Yuan 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 investigate the disease’s dynamics of a reaction-diffusion epidemic model. We first give a priori estimates of upper and lower bounds for positive solutions to model and then give the conditions of the existence and nonexistence of the positive nonconstant steady states, which guarantees the existence of the stationary patterns.

#### 1. Introduction

Infectious diseases are the second leading cause of death worldwide, after heart disease, and are responsible for more deaths annually than cancer [1]. Since the pioneer work of Kermark and McKendrick [2], mathematical models have been contributing to improve our understanding of infectious disease dynamics and help us develop preventive measures to control infection spread qualitatively and quantitatively.

Many studies indicate that spatial epidemiology with self-diffusion has become a principal scientific discipline aiming at understanding the causes and consequences of spatial heterogeneity in disease transmission [3]. In these studies, reaction-diffusion equations have been intensively used to describe spatiotemporal dynamics. In particular, the spatial spread of infections has been studied by analyzing traveling wave solutions and calculating spread rates [4–10].

Besides, there has been some research on pattern formation in the spatial epidemic model, starting with Turing’s seminal paper [11]. Turing’s revolutionary idea was that the passive diffusion could interact with chemical reaction in such a way that even if the reaction by itself has no symmetry-breaking capabilities, diffusion can destabilize the symmetric solutions with the result that the system with diffusion has them [12]. In these studies [3, 13–20], via standard linear analysis, the authors obtained the conditions of Turing instability, and, via numerical simulation, they showed the pattern formation induced by self-diffusion or cross-diffusion and found that model dynamics exhibits a diffusion controlled formation growth to stripes, spots, and coexistence or chaos pattern replication.

Recently, the researchers are interested in research on the stationary patterns due to the existence and nonexistence nonconstant solutions of the reaction-diffusion model [21–29]. But the research on the existence and nonexistence nonconstant solutions of reaction-diffusion epidemic model, seems rare [3].

In this paper, we will focus on the disease’s dynamics through studying the existence of the constant and nonconstant steady states of a simple reaction-diffusion epidemic model.

The rest of this paper is organized as follows. In Section 2, we derive a reaction-diffusion epidemic model. In Section 3, we give a priori estimates of upper and lower bounds for positive solutions to model. In Section 4, we give the main results on the existence and nonexistence of positive nonconstant steady states of the model. The paper ends with a brief discussion in Section 5.

#### 2. Basic Model

In [30], Berezovsky and coworkers introduced a simple epidemic model through the incorporation of variable population, disease induced mortality, and emigration into the classic model of Kermark and McKendrick [2]. The total population () is divided into two groups susceptible () and infectious (); that is, . The model describing the relations between the state variables is where the birth process incorporates density dependent effects via a logistic equation with the intrinsic growth rate and the carrying capacity ; , represent population densities of susceptible and infected population, respectively; denotes the transmission rate (the infection rate constant); is the natural mortality; denotes the disease-induced mortality; is the per-capita emigration rate of noninfective.

For model (1), the epidemic threshold of basic reproduction number is then computed as

The basic demographic reproductive number is given by

For simplicity, rescalling the model (1) by letting , , and leads to the following model: where defined by the ratio of the average life-span of susceptibles to that of infections and .

See [30] for more details.

Assume that the habitat () is a bounded domain with smooth boundary (when ), and is the outward unit normal vector on . We consider the following reaction-diffusion epidemic model: where the nonnegative constants and are the diffusion coefficients of and , respectively. The symbol is the Laplacian operator. The homogeneous Neumann boundary condition implies that the above model is self-contained and there is no infection across the boundary.

The corresponding kinetic model (5) with has been investigated by Wang et al. [20].

In this paper, we concentrated on the steady states of model (5) which satisfy

Throughout this paper, the positive solution satisfying model (6) refers to a classical one with , on . Clearly, model (6) has a unique positive constant solution (endemic equilibrium) if and , where

#### 3. A Priori Estimates for Positive Solutions to Model (6)

The main purpose of this section is to give a priori upper and lower positive bounds for positive solution of model (6). To this aim, we first cite two known results. The first is due to Lin et al. [31] and the second to Lou and Ni [32]. In the following, let us denote the constants , , and collectively by . The positive constants , , , , and so forth will depend only on the domains and .

Lemma 1 (Harnack inequality [31]). *Let be a positive solution to , where , satisfying the homogeneous Neumann boundary conditions. Then there exists a positive constant , such that
*

Lemma 2 (maximum principle [32]). *Let be a bounded Lipschitz domain in and .*(a)*Assume that and satisfies
* *If , then .*(b)*Assume that and satisfies
* *If , then .*

Theorem 3. *If , then the positive solution of model (6) satisfies
*

*Proof. *Assume that is a positive solution of model (6). We set
By applying Lemma 2, we have
and . This clearly gives

Theorem 4. *Assume that and . Let and be fixed positive constants. Then there exists a positive constant such that, if , every positive solution of model (6) satisfies
*

*Proof. *Let
In view of Theorem 3, there exists a positive constant such that , provided that . As and satisfy
It follows from Lemma 1 that there exists a positive constant such that
for .

Now, on the contrary, suppose that (15) is not true, then there exist sequences , with and the positive solution of model (6) corresponding to , such that
It follows from Lemma 1 that
satisfies
Integrating by parts, we obtain that, for ,
By the regularity theory for elliptic equations [33], we see that there exist a subsequence of , which we will still denote by , and two nonnegative functions , such that in as . By (20), we have that or .

Letting in (22) we obtain that
*Case 1 ( **, ** or **, **).* Since satisfies the second inequality of (18), on . Therefore, on as . Hence, for sufficiently large which contradicts the second integral identity of (22).*Case 2 ( **, **).* As above, on . It follows from the first integral identity of (23) that
This fact combines with yielding to , which implies that uniformly on as , since uniformly on . As , this contradicts the second integral identity of (23) and the fact that . This completes the proof.

#### 4. Existence and Nonexistence of Positive Nonconstant Steady States

In this section, we provide some sufficient conditions for the existence and nonexistence of nonconstant positive solution of model (6) by using the Leray-Schauder degree theory [34]. From now on, we denote by the eigenvalues of the operator on with the zero-flux boundary conditions.

##### 4.1. Nonexistence for Positive Nonconstant Steady States to Model (6)

This section is devoted to the consideration of the nonexistence for the nonconstant positive solutions of model (6), and, in the following results, the diffusion coefficients do play a significant role.

Theorem 5. *Assume that . Let be a fixed positive constant with . Then there exists a positive constant such that model (6) has no positive nonconstant solution provided that and .*

*Proof. *Let be any positive solution of model (6) and denote . Then, multiplying the first equation of model (6) by , integrating over , by virtue of Theorem 3, we have that
where depend only on . In a similar manner, we multiply the second equation in model (6) by to have
It follows from (26), (27) and the -Young inequality that
where . It follows from the well-known Poincaré inequality that
Since from the assumption, we can find a sufficiently small such that . Finally, by taking one can conclude that and , which asserts our results.

##### 4.2. Existence for Positive Nonconstant Steady States to Model (6)

In this section, we discuss the global existence of nonconstant positive classical solutions to model (6), which guarantees the existence of the *stationary patterns* [21, 24, 26, 27].

Unless otherwise specified, in this section, we always require that and , which guarantees that model (6) has one positive constant steady state . From now on, let us denote

Let and . Then we write model (6) in the form where

Define a compact operator by where is the inverse operator of subject to the zero-flux boundary condition. Then is a positive solution of model (31) if and only if satisfies

To apply the index theory, we investigate the eigenvalue of the problem where and with is an eigenvalue of (35) if and only if is an eigenvalue of the matrix for any . Therefore, is invertible if and only if, for any , the matrix is invertible. A straightforward computation yields where . For the sake of convenience, we denote Then .

If , then has two real roots given by

Set , , and the multiplicity of .

To compute , we can assert the following conclusion by Pang and Wang [22].

Lemma 6 (see [22]). *Suppose for all . Then
**
where
**
In particular, if for all , then .*

From Lemma 6, we see that to calculate the index of , the key step is to determine the range of for which .

Theorem 7. *Assume that . If , for some , and is odd, then there exists a positive constant such that model (6) has at least one nonconstant solution if .*

*Proof. *Since , equivalently, , it follows that if is large enough, then and . Furthermore,
Since for some , there exists such that

By Theorem 5, we know that there exists such that model (6) with diffusion coefficients and has no nonconstant solutions. Moreover, we can choose so large that . It follows that there exists such that
We shall prove that, for any , model (6) has at least one nonconstant positive solution. On the contrary, suppose that this assertion is not true for some . In the following, we will derive a contradiction by using a homotopy argument.

By virtue of Theorems 3 and 4, there exists a positive constant such that the positive solution of model (6) satisfies .

Set
and define
by
where

It is clear that finding the positive solution of model (31) is equivalent to finding the fixed point of in . Further, by virtue of the definition of , we have that has no fixed point in for all .

Since is compact, the Leray-Schauder topological degree is well defined. From the invariance of Leray-Schauder degree at the homotopy, we deduce

In view of and , we have . Clearly, . Thus, if model (6) has no other solutions except the constant one , then Lemma 6 shows that

On the contrary, by the choice of and , we have that and is the only fixed point of . It therefore follows from Lemma 6 that
From (50)–(52), we get a contradiction. Therefore, there exists a nonconstant solution of model (6). The proof is completed.

#### 5. Discussion

In this paper, we investigate the disease’s dynamics through studying the existence and nonexistence positive constant steady states of a reaction-diffusion epidemic model. We give a priori estimates for positive solutions to model and show that the nonconstant positive steady states exist due to the emergence of diffusion, which demonstrates that *stationary patterns* can be found as a result of diffusion. The numerical results about the *stationary patterns* for model (5) can be found in [20].

On the other hand, there are plenty of papers which focus on the pattern formation of reaction-diffusion population models via standard linear analysis method and numerical simulations. But there is little literature analytically concerning the existence of a stationary patterns via theory and methods of partial differential equations infrequently. The methods and results in the present paper may enrich the research of pattern formation in the spatial epidemic model.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors thank the anonymous referee for very helpful suggestions and comments which led to improvement of our original paper. This research was supported by the National Science Foundation of China (61373005) and Zhejiang Provincial Natural Science Foundation (LY12A01014).

#### References

- A. Lahrouz and A. Settati, “Asymptotic properties of switching diffusion epidemic model with varying population size,”
*Applied Mathematics and Computation*, vol. 219, no. 24, pp. 11134–11148, 2013. View at Publisher · View at Google Scholar - M. D. Kermark and A. G. Mckendrick, “Contributions to the mathematical theory of epidemics,”
*Proceedings of the Royal Society of London A*, vol. 115, pp. 700–721, 1927. View at Google Scholar - Y. Cai, D. Chi, W. Liu, and W. Wang, “Stationary patterns of a cross-diffusion epidemic model,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 852698, 10 pages, 2013. View at Publisher · View at Google Scholar - E. E. Holmes, M. A. Lewis, J. E. Banks, and R. R. Veit, “Partial differential equations in ecology: spatial interactions and population dynamics,”
*Ecology*, vol. 75, no. 1, pp. 17–29, 1994. View at Publisher · View at Google Scholar · View at Scopus - Y. Hosono and B. Ilyas, “Traveling waves for a simple diffusive epidemic model,”
*Mathematical Models and Methods in Applied Sciences*, vol. 5, no. 7, pp. 935–966, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - B. T. Grenfell, O. N. Bjørnstad, and J. Kappey, “Travelling waves and spatial hierarchies in measles epidemics,”
*Nature*, vol. 414, no. 6865, pp. 716–723, 2001. View at Publisher · View at Google Scholar · View at Scopus - A. L. Lloyd and V. A. A. Jansen, “Spatiotemporal dynamics of epidemics: synchrony in metapopulation models,”
*Mathematical Biosciences*, vol. 188, no. 1-2, pp. 1–16, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - K. Wang, W. Wang, and S. Song, “Dynamics of an HBV model with diffusion and delay,”
*Journal of Theoretical Biology*, vol. 253, no. 1, pp. 36–44, 2008. View at Publisher · View at Google Scholar · View at Scopus - R. Xu and Z. Ma, “An HBV model with diffusion and time delay,”
*Journal of Theoretical Biology*, vol. 257, no. 3, pp. 499–509, 2009. View at Publisher · View at Google Scholar · View at Scopus - S. Wang, W. Liu, Z. Guo, and W. Wang, “Traveling wave solutions in a reaction-diffusion epidemic model,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 216913, 13 pages, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. M. Turing, “The chemical basis of morphogenesis,”
*Philosophical Transactions of the Royal Society of London B*, vol. 237, no. 641, pp. 37–72, 1952. View at Publisher · View at Google Scholar - N. F. Britton,
*Essential Mathematical Biology*, Springer, 2003. - G. Sun, Z. Jin, Q.-X. Liu, and L. Li, “Pattern formation in a spatial S-I model with non-linear incidence rates,”
*Journal of Statistical Mechanics*, vol. 2007, no. 11, Article ID P11011, 2007. View at Publisher · View at Google Scholar · View at Scopus - M. Bendahmane and M. Langlais, “A reaction-diffusion system with cross-diffusion modeling the spread of an epidemic disease,”
*Journal of Evolution Equations*, vol. 10, no. 4, pp. 883–904, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - Y. Cai and W. Wang, “Spatiotemporal dynamics of a reaction-diffusion epidemic model with nonlinear incidence rate,”
*Journal of Statistical Mechanics*, vol. 2011, no. 2, Article ID P02025, 2011. View at Publisher · View at Google Scholar · View at Scopus - W. Wang, Y. Lin, H. Wang, H. Liu, and Y. Tan, “Pattern selection in an epidemic model with self and cross diffusion,”
*Journal of Biological Systems*, vol. 19, no. 1, pp. 19–31, 2011. View at Publisher · View at Google Scholar · View at Scopus - W.-M. Wang, H.-Y. Liu, Y.-L. Cai, and Z.-Q. Li, “Turing pattern selection in a reaction-diffusion epidemic model,”
*Chinese Physics B*, vol. 20, no. 7, Article ID 074702, 2011. View at Publisher · View at Google Scholar · View at Scopus - S. Berres and R. Ruiz-Baier, “A fully adaptive numerical approximation for a two-dimensional epidemic model with nonlinear cross-diffusion,”
*Nonlinear Analysis*, vol. 12, no. 5, pp. 2888–2903, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - F. Rao, W. Wang, and Z. Li, “Stability analysis of an epidemic model with diffusion and stochastic perturbation,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 17, no. 6, pp. 2551–2563, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - W. Wang, Y. Cai, M. Wu, K. Wang, and Z. Li, “Complex dynamics of a reaction-diffusion epidemic model,”
*Nonlinear Analysis*, vol. 13, no. 5, pp. 2240–2258, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - M. Wang, “Non-constant positive steady states of the Sel'kov model,”
*Journal of Differential Equations*, vol. 190, no. 2, pp. 600–620, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - P. Y. H. Pang and M. Wang, “Qualitative analysis of a ratio-dependent predator-prey system with diffusion,”
*Royal Society of Edinburgh A*, vol. 133, no. 4, pp. 919–942, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - P. Y. H. Pang and M. Wang, “Strategy and stationary pattern in a three-species predator-prey model,”
*Journal of Differential Equations*, vol. 200, no. 2, pp. 245–273, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - M. Wang, “Stationary patterns for a prey-predator model with prey-dependent and ratio-dependent functional responses and diffusion,”
*Physica D*, vol. 196, no. 1-2, pp. 172–192, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - M. Wang, “Stationary patterns caused by cross-diffusion for a three-species prey-predator model,”
*Computers and Mathematics with Applications*, vol. 52, no. 5, pp. 707–720, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - R. Peng, J. Shi, and M. Wang, “Stationary pattern of a ratio-dependent food chain model with diffusion,”
*SIAM Journal on Applied Mathematics*, vol. 67, no. 5, pp. 1479–1503, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - R. Peng, M. Wang, and G. Yang, “Stationary patterns of the holling-tanner prey-predator model with diffusion and cross-diffusion,”
*Applied Mathematics and Computation*, vol. 196, no. 2, pp. 570–577, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - R. Peng and J. Shi, “Non-existence of non-constant positive steady states of two Holling type-II predator-prey systems: strong interaction case,”
*Journal of Differential Equations*, vol. 247, no. 3, pp. 866–886, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - C. Bianca, “Existence of stationary solutions in kinetic models with gaussian thermostats,”
*Mathematical Methods in the Applied Sciences*, vol. 36, no. 13, pp. 1768–1775, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - F. Berezovsky, G. Karev, B. Song, and C. Castillo-Chavez, “A simple epidemic model with surprising dynamics,”
*Mathematical Biosciences and Engineering*, vol. 2, no. 1, pp. 133–152, 2005. View at Google Scholar · View at Zentralblatt MATH - C.-S. Lin, W.-M. Ni, and I. Takagi, “Large amplitude stationary solutions to a chemotaxis system,”
*Journal of Differential Equations*, vol. 72, no. 1, pp. 1–27, 1988. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - Y. Lou and W.-M. Ni, “Diffusion, self-diffusion and cross-diffusion,”
*Journal of Differential Equations*, vol. 131, no. 1, pp. 79–131, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - D. Gilbarg and N. S. Trudinger,
*Elliptic Partial Differential Equations of Second Order*, Springer, 1983. - L. Nirenberg,
*Topics in Nonlinear Functional Analysis*, AMS Bookstore, 2001.