Abstract

The aim of this paper is to study the dynamics of a reaction-diffusion SIR epidemic model with specific nonlinear incidence rate. The global existence, positivity, and boundedness of solutions for a reaction-diffusion system with homogeneous Neumann boundary conditions are proved. The local stability of the disease-free equilibrium and endemic equilibrium is obtained via characteristic equations. By means of Lyapunov functional, the global stability of both equilibria is investigated. More precisely, our results show that the disease-free equilibrium is globally asymptotically stable if the basic reproduction number is less than or equal to unity, which leads to the eradication of disease from population. When the basic reproduction number is greater than unity, then disease-free equilibrium becomes unstable and the endemic equilibrium is globally asymptotically stable; in this case the disease persists in the population. Numerical simulations are presented to illustrate our theoretical results.

1. Introduction

In this paper, we consider the following SIR epidemic model with a specific nonlinear incidence rate described by where , , and are susceptible, infectious, and recovered classes, respectively. is the recruitment rate of the population, is the natural death rate of the population, is the death rate due to disease, is the recovery rate of the infective individuals, is the infection coefficient, and is the incidence rate, where are constants. It is very important to note that this incidence rate becomes the bilinear incidence rate if , the saturated incidence rate if or , the Beddington-DeAngelis functional response introduced in [1, 2] and used in [3] when , and Crowley-Martin functional response presented in [46] if . Moreover, the function satisfies the hypotheses , , and of general incidence rate presented by Hattaf et al. in [7]. From the biological point of view, the transmission rate of infectious diseases remains unknown in detail and may be different from one disease to another. In the classical epidemic models, this rate was assumed to be linear with respect to the numbers of susceptible and infected individuals. This assumption is based on the law of mass action which is more appropriate for communicable diseases such as influenza but not for sexually transmitted diseases such as HIV/AIDS. For one reason, the transmission rate in system (1) is assumed to be nonlinear and has the form that measures the saturation effect which represents that the number of individual contacts reaches a certain maximum value due to social or spatial distribution of the population. For more details on the choice of the nonlinearity of the incidence rate, we refer the reader to the book of Capasso [8].

On the other hand, the spatial content of the environment has been ignored in the model (1). However, due to the large mobility of people within a country or even worldwide, spatially uniform models are not sufficient to give a realistic picture of disease diffusion. For this reason, the spatial effects cannot be neglected in studying the spread of epidemics.

Therefore, we consider the following SIR epidemic model with specific nonlinear incidence rate and spatial diffusion: where , , and represent the numbers of susceptible, infected, and removed individuals at location and time , respectively. The positive constants , , and denote the corresponding diffusion rates for these three classes of individuals.

The aim of this work is to investigate the global dynamics of the reaction-diffusion system (2). Note that does not appear in the first two equations; this allows us to study the system with homogeneous Neumann boundary conditions and initial conditions Here, is a bounded domain in with smooth boundary . and are, respectively, the normal derivatives of and on .

The rest of paper is organized as follows. The next section deals with the global existence, positivity, and boundedness of solutions of problem (3)–(5). In Section 3, we discuss the stability analysis of equilibria. In Section 4, we present the numerical simulation to illustrate our result. Finally, the conclusion of our paper is in Section 5.

2. Global Existence, Positivity, and Boundedness of Solutions

In this section, we establish the global existence, positivity, and boundedness of solutions of problem (3)–(5) because this model describes the population. Hence, the population should remain nonnegative and bounded.

Proposition 1. For any given initial data satisfying the condition (5), there exists a unique solution of problem (3)–(5) defined on and this solution remains nonnegative and bounded for all .

Proof. System (3)–(5) can be written abstractly in the Banach space of the form where , , , and It is clear that is locally Lipschitz in . From [9], we deduce that system (6) admits a unique local solution on , where is the maximal existence time for solution of system (6).
In addition, system (3) can be written in the form It is easy to see that the functions and are continuously differentiable satisfying and for all . Since initial data of system (3) are nonnegative, we deduce the positivity of the local solution (see the book of Smoller [10]).
Now, we show the boundedness of solution. From (3)–(5) we have By the comparison principle [11], we have , where is the solution of the problem Since for , we have that
From Theorem 2 given by Alikakos in [12], to establish the uniform boundedness of , it is sufficient to show the uniform boundedness of .
Since and , we get Hence, which implies that . Using [12, Theorem 3.1], we deduce that there exists a positive constant that depends on and on such that
From the above, we have proved that and are bounded on . Therefore, it follows from the standard theory for semilinear parabolic systems (see [13]) that . This completes the proof of the proposition.

3. Qualitative Analysis of the Spatial Model

Using the results presented by Hattaf et al. in [7], it is easy to get that the basic reproduction number of disease in the absence of spatial dependence is given by which describes the average number of secondary infections produced by a single infectious individual during the entire infectious period.

It is not hard to show that the system (3) is always a disease-free equilibrium of the form . Further, if , the system (3) has an endemic stationary state where with and .

The objective of this section is to discuss the local and global stability of the equilibria.

3.1. Local Stability of the Equilibria

First, we linearize the dynamical system (3) around arbitrary spatially homogeneous fixed point for small space- and time-dependent fluctuations and expand them in Fourier space. For this, let where and ; and are the wavenumber vector and frequency, respectively. Then we can obtain the corresponding characteristic equation as follows: where is the identity matrix, is the diffusion matrix, and is the Jacobian matrix of (3) without diffusion atwhich is given by The characterization of the local stability of disease-free equilibrium is given by the following result.

Theorem 2. The disease-free equilibrium is locally asymptotically stable if and it is unstable if .

Proof. Evaluating (18) at , we have Clearly, the roots of (20) are and . Note that is negative if for all . Hence is locally asymptotically stable if . If , is positive. So is unstable.

Next, we focus on the local stability of the endemic equilibrium .

Theorem 3. The endemic equilibrium is locally asymptotically stable if .

Proof. Evaluating (18) at , we have where We have and ; then is locally asymptotically stable.

3.2. Global Stability of the Equilibria

The purpose of this subsection is to determine the global stability for reaction-diffusion equations (3)–(5) by constructing Lyapunov functionals. These Lyapunov functionals are obtained from those for ordinary differential equations (1) by applying the method of Hattaf and Yousfi presented in [14].

The system (1) is particular case of the model proposed by Hattaf et al. [7] with . To study the global stability of for (1), the authors Hattaf et al. [7] proposed the following Lyapunov functional: where .

From [14], we construct the Lyapunov functional for system (3)–(5) at as follows: Calculating the time derivative of along the solution of system (3)–(5), we have Since , we have . Thus, the disease-free equilibrium is stable, and if and only if and . We discuss two cases as follows.(i)If , then .(ii)If , from and the first equation of (3), we have Then, . Since and , then . By the above discussion, we deduce that the largest compact invariant set in is just the singleton . From LaSalle invariance principle [15], we conclude that is globally asymptotically stable.

Using same technique, we construct a Lyapunov functional for system (3)–(5) at from the Lyapunov functional defined by Hattaf et al. in [7]. It is easy to show that verifies the condition (15) given in [14]. Hence, it follows from [14, Proposition 2.1] that is a Lyapunov functional for the reaction-diffusion system (3)–(5) at when . We summarize the above in the following result.

Theorem 4. (i) If , the disease-free equilibrium of (3)–(5) is globally asymptotically stable for all diffusion coefficients.
(ii) If , the endemic equilibrium of (3)–(5) is globally asymptotically stable for all diffusion coefficients.

4. Numerical Simulations

In this section, we present the numerical simulations to illustrate our theoretical results. To simplify, we consider system (3) under Neumann boundary conditions and initial conditions

In Figure 1, we show that the maximum value of the numbers of susceptibles and infectious individuals is concentrated at the middle of the interval and these numbers decrease linearly to zero at the boundaries and .

Now, we choose the following data set of system (3): , , , , , , , , and . By calculation, we have . In this case, system (3) has a disease-free equilibrium . Hence, by Theorem 4(i), is globally asymptotically stable. Numerical simulation illustrates our result (see Figure 2).

In Figure 3, we choose and do not change the other parameter values. By calculation, we have which satisfy Theorem 4(ii); then the disease-free equilibrium is still present and the system (2) has a unique endemic equilibrium . Therefore, by Theorem 2 and Theorem 4(ii), is unstable, while is globally asymptotically stable. Numerical simulation illustrates our result (see Figure 3).

5. Conclusion

In this paper, we investigated the dynamics of a reaction-diffusion epidemic model with specific nonlinear incidence rate. This specific nonlinear incidence rate includes the traditional bilinear incidence rate, the saturated incidence rate, the Beddington-DeAngelis functional response, and Crowley-Martin functional response. The global dynamics of the model are completely determined by the basic reproduction number . We proved that the disease-free equilibrium is globally asymptotically stable if ,, which leads to the eradication of disease from population. When then disease-free equilibrium becomes unstable and a unique endemic equilibrium exists and is globally asymptotically stable, which means that the disease persists in the population.

From our theoretical and numerical results, we conclude that the spatial diffusion has no effect on the stability behavior of equilibria in the case of Neumann conditions and spatially constant coefficients.

Conflict of Interests

The authors declare that they have no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank the editor and the anonymous referees for their valuable remarks and comments which have led to the improvement of the quality of their study.