Mathematical Problems in Engineering

Volume 2015 (2015), Article ID 235935, 10 pages

http://dx.doi.org/10.1155/2015/235935

## The Impact of Awareness Programs with Recruitment and Delay on the Spread of an Epidemic

Department of Mathematics, North University of China, Taiyuan, Shanxi 030051, China

Received 2 April 2014; Revised 18 July 2014; Accepted 1 August 2014

Academic Editor: Yuxin Zhao

Copyright © 2015 Lixia Zuo 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

A compartment epidemic model with delay is given to discuss the impact of awareness programs on the spread and control of infectious diseases in a given region. It is assumed that there is a constant recruitment rate in the cumulative density of awareness programs, and further it is assumed that awareness programs can influence the susceptible to a limited extent. The system exhibits two equilibria: the disease-free equilibrium is stable if the basic reproduction number is less than unity for any delay and the unique endemic equilibrium exhibits Hopf-bifurcation under certain conditions. Numerical simulations prove the results of analysis and the significance of awareness programs in preventing and controlling the diseases, by investigating the relationship between the proportion of the infective and the dissemination rate and the implementation rate, respectively.

#### 1. Introduction

Plenty of evidence shows that awareness programs, which can influence the susceptible to a limited extent due to some objective factors, play an important role in the spread and control of infectious diseases. For example, during the outbreak of SARS, H1N1 influenza pandemic, and HIV epidemic, public media had massive reports on the number of the infections and deaths per day, which had a great impact on the diseases control [1–4]. That is due to the fact that the spread of diseases is often accompanied by a rise in awareness of those in the social vicinity of infected individuals and a subsequent change in behavior, such as keeping social distancing, wearing protective masks, and vaccination [5, 6]. Such reactions can manifest themselves in lower susceptibility as people try to prevent themselves from catching the disease, but also in lower infectivity because of self-imposed quarantine or better hygiene, shorter durations of infectiousness, or longer immunity [7]. And once the infective are cured, they will be aware of the disease [8].

Recently some scholars used mathematical models to discuss the impact of awareness programs on the diseases spreading and controlling in a given region [6, 8–11]. Joshi et al. [9] formulated a model to investigate the effect on the HIV epidemic in Uganda and compared their model with three types of the susceptible to a standard SIR model and then pointed out that the awareness programs in Uganda are successful in combating diseases. Li and Cui [1] analyzed a SIS epidemic model incorporating media coverage under constant and pulse vaccination; then they obtained the exact periodic infection-free solution which is globally asymptotically stable under some conditions. In order to better describe populations mixed condition, some authors studied infectious diseases models on different networks [12–14]. In [12], the authors introduced a constant to represent the density level of media coverage from other regions. Funk et al. [14] formulated and analyzed a mathematical model in a host population; then they put forward that awareness programs can result in a lower size of the outbreak but do not affect the epidemic threshold, and if the behavioral response is treated as a local effect arising in the proximity of an outbreak, it can completely stop a disease from spreading, although only if the infection rate is below a threshold.

Some scholars consider some other factors in awareness program models. Liu and Wang et al. took into account the random perturbation [15, 16]. In [16], the authors extended a deterministic SIRS epidemic model to a stochastic differential equation, and then they discussed the exponential* p*-stability and global stability of unique positive solution. Some scholars focused their attention at the contact rate, and most of them assumed that awareness programs will aid in modifying the contact rate between the susceptible and the infective [17–29]. Liu and Cui [17] used the contact rate and found that both of the two equilibria are asymptotically stable. Mirsa et al. [20] used the constant with the contact rate to limit the effect of awareness programs on the susceptible and sought out the conditions that Hopf-bifurcation occurs. Tchuenche and Bauch [23] used an exponentially decreasing function to affect the force of infection; then numerical results showed the potential short-term beneficial effect of awareness programs. Pang and Cui [26] used to show the contact rate after awareness programs alerts and found that though it is not a determined fact to eradicate the infection of the diseases, the effective awareness programs can postpone the arrival of the infection peak. Elenbaas et al. [28] introduced a segmented function to describe the media impact when they formulate an epidemic model. A Filippov epidemic model was proposed to describe the real characteristics of media impact on the spread of infectious diseases by incorporating a piecewise continuous transmission rate in [29]. Mathematical and bifurcation analyses with regard to the local and global stability of equilibria and local sliding bifurcations are performed.

But what we regard as unreasonable is that most of the articles assume that the cured infective become unaware of the disease. In fact people have a certain consciousness about the disease once they get sick. Therefore we propose a delayed mathematical model for predicting the future course of any epidemic by considering some of the infect join the aware susceptible after recovery. We hold the opinion that awareness programs can influence the susceptible to a limited extent for some objective factors and then consider the interaction between the susceptible and awareness programs as Holling type-II functional response. In addition we make a constant to represent the density level of media coverage from other regions with the disease because other regions can also effect the region that we consider. In fact the results about global stability of other delayed systems could be further utilized for other related problems [30–32].

The rest of this paper is organized as follows. In the next section, a mathematical model with delay has been proposed to capture the dynamics of the effect of awareness programs. It is assumed that diseases spread due to the contact between the susceptible and the infective only. Then we analyze the conditions of the stability of equilibria and the existence of Hopf-bifurcation in Section 3. Furthermore, in Section 4 we perform some numerical examples to validate the analysis in Section 3 and then introduce the importance of the dissemination rate, implementation rate, and the delay in disease control. In Section 5 we discuss the above contents.

#### 2. Mathematical Model and Equilibrium Analysis

In the region under consideration the rate of immigration of the susceptible is . It is assumed that the disease spreads due to the direct contact between the susceptible and the infective only and due to awareness programs the susceptible avoid being in contact with the infective and form a different class with a proportion , named the aware susceptible. So the total population is divided into three classes: the susceptible, the aware susceptible, and the infective, the proportions of which at time in the total population are , , and . It is assumed that the aware susceptible may lose awareness with passage of time and become susceptible with the proportion of again. The infective can be cured with a proportion and a fraction of recovered people will become aware and join the aware susceptible class whereas the remaining fraction will join the susceptible.

Consider that the cumulative density of awareness programs driven by media in that region at time is , which is related to the infective. We make the constant represent the executed rate of awareness programs. As the time passes, some campaigns lose their impact on people and lead to the diminution of the awareness programs, so we introduce to denote the rate of their depletion. Moreover represents the density level of media coverage on the disease from other regions. Using the fact that there is a limited extent of the awareness programs influence on the susceptible due to some objective factors, we introduce to limit the effect and consider the interaction between them as Holling type-II functional response. It is plausible that the policy makers need some time to gather the cases of the infective generally; then we introduce and consider that the cumulative density of awareness programs at time being executed will be in accordance with the infected cases reported at time . Keeping the above facts in mind, the dynamics of model is governed by the following system of nonlinear delay differential equations: Here , , and for and .

In the above model, the constants and , respectively, represent the contact rate of the unaware susceptible with the infective and the natural death rate. All the constants in the system are assumed to be positive.

Using the fact that , system (1) is reduced to the following system: Now it is sufficient to study system (2) in detail rather than system (1).

For the analysis of system (2), we need the region of attraction which is given by the set , and it attracts all solutions initiating in the interior of the positive orthant.

The system (2) has two equilibria.(i)Disease-free equilibrium .(ii)Endemic equilibrium .

Define the basic reproduction number . The existence of equilibrium is trivial; then we prove the existence of in detail. When , in the equilibrium the values of , are with satisfying the equation where

Solving (4) we get We obtain , and when , then we get for .

*Remark 1. *From the expression of , it is easy to note that and , which shows that the equilibrium proportion of the infective decreases as the rate of dissemination of awareness among the susceptible and the implementation rate of awareness programs increases.

#### 3. Stability Analysis

In this section we present the local stability of , and explore the conditions of Hopf-bifurcation by taking delay as a bifurcation parameter.

##### 3.1. Stability of Equilibria without Delay ()

Theorem 2. *The disease-free equilibrium is locally asymptotically stable if .*

* Proof. *The Jacobian matrix corresponding to the system (2) when is given below:The characteristic equation at is of the form where is the eigenvalue. We get so when . Thus is locally asymptotically stable if .

The form of characteristic equation at is where

Theorem 3. *When the endemic equilibrium exists, it is locally asymptotically stable provided .*

*Proof. *For the characteristic equation it is easy to show And where so . According to Hurwitz criterion, we can claim that all the eigenvalues will be either negative or having negative real part. Thus, the equilibrium is locally asymptotically stable.

##### 3.2. Stability of Equilibria with Delay ()

Theorem 4. *The disease-free equilibrium is locally asymptotically stable if .*

Similar to the proof of Theorem 2, this proof is omitted.

We linearize system (2) about to study the stability of the endemic equilibrium and get whereIn the above, , , and are small perturbations around .

Then the form of the characteristic equation of the system is where is the eigenvalue.

To show the Hopf-bifurcation, we need to show that (19) has a pair of purely imaginary roots. For this, substituting into (19) and separating real and imaginary parts, we get the following transcendental equations: Squaring and adding the above equations and substituting , we get where , , and .

If the coefficients in satisfy the conditions of Routh-Hurwitz criterion, then is locally asymptotically stable for all delays , provided it is stable in absence of delay. In the following we consider that the values of do not satisfy the Routh-Hurwitz criterion.

Lemma 5. *If (21) will have at least one positive root .*

* Proof. *For , it is obvious that Then we have and ; thus (21) has at least one positive root .

Denote ; then (21) has a pair of purely imaginary roots .

Now we turn to the bifurcation analysis. We use the delay as the bifurcation parameter. We view the solutions of (21) as a function of the bifurcation parameter , and let be the eigenvalue of (21) such that, for the initial value of the bifurcation parameter , we have and (). To establish the Hopf bifurcation at , we need to show that .

Lemma 6. *One has the following transversality condition: *

* Proof. *Differentiating with respect to from (21), we get So This proves the lemma. Now we have the following theorem.

*Theorem 7. The endemic equilibrium of the system is locally asymptotically stable when and becomes unstable for provided When , a Hopf-bifurcation occurs, leading a family of periodic solutions bifurcating from as passes through the critical value , where *

*4. Numerical Simulations and Results*

*4. Numerical Simulations and Results*

*To check the feasibility of our analysis in Section 3 about , we present some numerical computations in this section using MatLab by choosing the following set of parameter values: , , , , , , , , , and (see Figure 1), and then we get Figure 2 with different initial values, from which we can know that is stable and all the trajectories approach it. The numerical simulations support the analysis given.*