## Advanced Nonlinear Dynamics of Population Biology and Epidemiology

View this Special IssueResearch Article | Open Access

# Existence of Traveling Wave Solutions for Cholera Model

**Academic Editor:**Kaifa Wang

#### Abstract

To investigate the spreading speed of cholera, CodeÃ§oâ€™s cholera model (2001) is developed by a reaction-diffusion model that incorporates both indirect environment-to-human and direct human-to-human transmissions and the pathogen diffusion. The two transmission incidences are supposed to be saturated with infective density and pathogen density. The basic reproduction number is defined and the formula for minimal wave speed is given. It is proved by shooting method that there exists a traveling wave solution with speed for cholera model if and only if .

#### 1. Introduction

Cholera has been a serious threat to human health in the past and at present, which is an acute, diarrheal illness caused by infection of the intestine with the bacterium* Vibrio cholera*. An estimated 3â€“5 million cases and over 100,000 deaths occur each year around the world [1]. The cholera bacterium is usually found in water or food sources that have been contaminated by feces from a person infected with cholera. Cholera is most likely to be found and to spread in places with inadequate water treatment, poor sanitation, and inadequate hygiene. Therefore, cholera outbreaks have occurred in developing countries, for example, Iraq (2007-2008), Guinea Bissau (2008), Zimbabwe (2008-2009), Haiti (2010), Democratic Republic of Congo (2011-2012), and Sierra Leone (2012) [2].

To understand the propagation mechanism of cholera, many mathematical models were proposed, whose earlier one was established by Capasso and Paveri-Fontana [3] to study the 1973 cholera epidemic in the Mediterranean region as follows:
where and denote the concentrations of the pathogen and the infective populations, respectively. In addition, CodeÃ§o [4] investigated the role of the aquatic pathogen in dynamics of cholera through the following susceptible-infective-pathogen model:
where is the susceptible individuals. In this model, human is divided into two groups: the susceptible group and the infective group. As pointed out in [4â€“8], bacterium* Vibrio cholera* can spread by direct human-to-human and indirect environment-to-human modes. To understand the complex dynamics of cholera, model (2) is extended by [8â€“15] and so forth.

In all previous models the influences of space distribution of human on the transmission of cholera are omitted. Cholera usually spreads in spatial wave [16]. Cholera bacteria live in rivers and interact with the plankton on the surface of the water [17]. When individuals drink contaminated water and are infected, they will release cholera bacteria through excretion [18]. Capasso et al. [19â€“23] developed model (1) by incorporating the bacterium diffusion in a bounded area and studied the existence and stability of solutions. To deeply investigate the interaction of transmission modes and bacterium diffusion, Bertuzzo et al. [24, 25] incorporated patchy structure into model (2) and supposed that pathogen in water could diffuse among these patches. Furthermore, Mari et al. [26] studied the influence of diffusion of both human and pathogen on cholera dynamics through a patchy model.

Infectious case is usually found firstly at some location and then spreads to other areas. Consequently, the most important question for cholera is what the spreading speed of cholera is. However, the above spatial models mainly focus on the stability of solutions not the spreading speed. Traveling wave solution is an important tool used to study the spreading speed of infectious diseases [27â€“29]. Based on Capassoâ€™s model (1), Zhao and Wang [30], Xu and Zhao [31], Jin and Zhao [32], and Hsu and Yang [33] studied the influences of pathogen diffusion on the spread speed of cholera.

The studies of traveling wave solutions of Capassoâ€™s model (1) incorporating pathogen diffusion provide insight into the spreading speed of cholera. However, some pieces of information are omitted, such as the interaction of direct human-to-human and indirect environment-to-human transmissions. In this paper, a reaction-diffusion model with pathogen diffusion and both transmission paths is proposed by developing CodeÃ§oâ€™s model (2). Based on model (2) and ignoring the disease-related death, a general diffusive cholera model can be formulated as the following reaction-diffusion system:
where and denote the concentrations of susceptible and infected individuals, respectively, and is the concentration of the infectious agents. is the total human population, stands for the natural birth and death rate, denotes the contribution of each infected person to the concentration of cholera, and is the net death rate of* vibrio cholera*. and are the human-to-human and environment-to-human transmission incidences, respectively. Similar to [10], we assume that and satisfy (A1), , ;(A2), , , , and is strictly monotonously increasing in .It is easy to conclude that , , and and are nonincreasing. Obviously, hypotheses (A1) and (A2) imply that the two transmission paths are saturated. In Tian and Wang [10], and have the following expressions:
Obviously, as a special case, such selections satisfy (A1) and (A2).

Shooting method is very important in proving the existence of traveling wave solutions, which was proposed by Dunbar [34, 35] and was applied to many models (e.g., [36â€“40]). In this paper, the existence of traveling wave solutions of system (3) will be proved by shooting method and the formula for minimal wave speed will be given.

This paper is organized as follows. In next section, the main theorem and the formula for minimal wave speed will be given. In Section 3, the nonexistence of the traveling wave solutions for is proved by geometric method. Section 4 is devoted to shooting arguments and the construction of Wazewski set. In Section 5, we prove the existence of traveling wave solutions for and then give the existence of traveling wave solution for by limit arguments. The final section is devoted to the simulations.

#### 2. Main Results

For convenience, we introduce dimensionless variables and parameters. By setting model (3) has the form where and .

Denote , which is the basic reproduction number of (6). Then hypotheses (A1) and (A2) imply that system (6) has two nonnegative constant solutions and if and only if , where is the only one positive root of equation and . Biologically, corresponds to disease-free equilibrium and corresponds to endemic equilibrium. To study the spreading wave of cholera, it is assumed that holds in this paper; that is

A traveling wave solution of system (6) with speed is a nonnegative solution of the form

Substituting traveling profile into system (6) yields the following equations: where denotes . To investigate invasion question by cholera, we will study the positive solutions of (10) such that

Before giving the main theorem, we introduce the equation for minimal wave speed where

Theorem 1. *There exists a constant which is the greatest positive root of (12). When , system (6) has a traveling wave solution satisfying boundary condition (11). When , system (6) has no traveling wave solutions satisfying boundary condition (11).*

#### 3. Nonexistence of Traveling Wave Solutions for

From (10), we have

Consequently, if , then

Hence, the traveling profile with boundary condition (11) must satisfy

Therefore, to study traveling wave solutions we assume (16) satisfies. Setting in system (10) and noticing (16), it follows

If , then by system (10). Therefore, we suppose for any ; that is, .

Obviously, system (17) has two equilibria and . A profile solution of (10) which satisfies boundary condition (11) corresponds to the positive solution of system (17) which satisfies where . Therefore, to study the solutions of (10), it is sufficient to study those of system (17) satisfying boundary condition (18).

Firstly, we investigate the dynamics near . Simple calculations show that the characteristic equation of the linearization of system (17) at is where

Because (19) has a negative real root, which is denoted by . Let and be the other two eigenvalues of (19) and suppose that . To investigate the distribution of roots of (19), denote and introduce the following lemma [41].

Lemma 2. * â€‰ â€‰If , (19) has one real root and two nonreal complex conjugate roots.**â€‰ â€‰If , (19) has a multiple root and all its roots are real.**â€‰ â€‰If , (19) has three distinct real roots.**Direct calculations show that , where is defined by (12).*

Lemma 3. *â€‰â€‰The real parts of and are positive.**â€‰â€‰Assume . Then, there exists which is the only positive root of . When , , and are real. When , , and are complex and nonreal.**â€‰â€‰Assume that . Then, there exist two positive constants which are all positive roots of . and are complex and nonreal if and only if . If , then ; if , then , where .**â€‰â€‰ if and only if or .*

*Proof. *Suppose is the root of (19). Substituting into (19) and comparing real and imaginary parts show that and . Since , then . However, it is impossible that and by the expressions of and . Therefore, the real parts of and are not zero. Furthermore, since it is impossible that and , Routh-Hurwitz theorem implies that it is impossible that the real parts of both and are negative. Consequently, there are two cases: (i) and are complex conjugate roots with positive real parts; (ii) and are real and at least one is positive. However, Descartesâ€™ rule of signs shows that the number of positive roots of (19) is zero or two. Thus, if case (ii) is true, both of and are real and positive. Therefore, (a) is proved.

In this paragraph, we consider the case . Firstly, suppose that . Obviously, and . By the expression of , we have
since . Now, assume ; that is, . Then. , , , and . Then, if , Descartes' rule of signs shows that there exists which is the only positive root of , where for and for . Using Lemma 2 completes the proof of (b).

Suppose that in this paragraph and, thus, . Calculations show that
and that has two roots and , where
and . By letting and using trivial calculations, we get (see Figure 1)

Therefore, if , then . Since is the only minimum-value point of , and then for any and both of and are not real. Lemma 2 shows that . Thus, since and , there exist two positive roots for equation such that . Then, using (25) and Lemma 2 completes the proof of (c) and (d).

**(a)**

**(b)**

Direct calculations show that corresponding eigenvectors of eigenvalue are where . Since and thus

Then, we have the following lemma.

Lemma 4. *If , there exist no traveling wave solutions which satisfy boundary condition (11).*

*Proof. *Assume that and . Then, (b) of Lemma 3 implies that and are complex conjugate eigenvalues and there exits locally unstable manifold and locally stable manifold . If a solution of (17) tends to when , then it will be spiral on . By the structures of and , at some time , which shows that there exist no traveling wave solutions departing from .

Suppose that . If , (c) of Lemma 3 shows that and are complex conjugate eigenvalues and similar arguments to that of previous paragraph finish the proof. If , (c) of Lemma 3 shows that and are real; however, . If a solution of (17) tends to when , structures of and indicate that there is an such that . The proof is completed.

From Section 4 to Section 5.2, we suppose that , which implies .

#### 4. Shooting Method and Wazewski Set

To prove the existence of traveling wave, shooting method developed by Dunbar [34] is used. Firstly, we give the shooting arguments.

Consider the differential equation where from to satisfies Lipschitz condition about . Let denote the unique solution of (29) with initial value . It is convenient to give the notations and . To describe the shooting method (or Wazewski theorem), some definitions are necessary.

*Definition 5. *(a) For , define immediate exit set of as

(b) For , let .

(c) Given , define exit time of by

Then, Wazewski theorem is formulated as follows.

Lemma 6 (see [34]). *Suppose that *(1)*if and , then .*(2)*If , and , then there exists an open set about disjoint from .*(3)*If , is compact and intersects a trajectory of (29) only once.**Then, the mapping is a homeomorphism from to its image on .*

A set satisfying conditions (1) and (2) of Lemma 6 is called a Wazewski set. In the following, we first construct the Wazewski set . Fundamental idea to construct a Wazewski set is that the characteristic vectors corresponding eigenvalues with positive real parts should be removed from and that those characteristic vectors corresponding eigenvalues with negative real parts should be included. Therefore, we set where

It is obvious that . Firstly, we give the construction of , which is described in Figure 2.

Lemma 7. *The construction of is as follows:
**
where .*

*Proof. *It is enough to analyze the behavior of solution on . We only study and omit the proof of since the analysis of is similar to that of and is simpler. In the process of this proof, we use some notations to simplify the proof. Set

From hypotheses (A1) and (A2), we find that and are monotonously decreasing, is strictly monotonously decreasing for , and is the only positive root of . The set is classified into two cases according to variable . â€‰(a) Case . This case is classified as follows.â€‰(1) Case . Then and the solution of (17) will enter .â€‰(2) Case . Then
â€‰The solution of (17) will enter .â€‰(3) Case . Then
â€‰The solution of (17) will enter .â€‰(4) Case . Then and the solution of (17) will enter .â€‰(5) Case and . The solution of (17) will enter .â€‰(6) Case . Then ,
â€‰and . Therefore, the solution of (17) will enter .â€‰(b)â€‰â€‰Case . This case is classified as follows.â€‰(1) Case . Then and the solution of (17) will enter .â€‰(2) Case . Then , . The solution of (17) will enter .â€‰(3) Case . Then , , and . The solution of (17) will enter .â€‰(4) Case . Then and . The solution of (17) will enter .â€‰(5) Case . In this case, is equilibrium and is constant.â€‰(6) Case and . Then and . The solution of (17) will enter .â€‰(7) Case . Then is equilibrium and is constant.

The proof is completed.

#### 5. Existence of Traveling Wave Solution for

In this section, we prove the existence of traveling wave solution for . Firstly, we study the behaviors of solutions near .

##### 5.1. Behaviors of Solutions Near

Lemma 8. *Suppose is a solution of (17) satisfying initial conditions
**
where . Then, for every , we have
*

*Proof. *From Lemma 3, we have . To finish the proof, it is sufficient to prove that the set
is positively invariant. It is obvious that
where

Suppose that . Then, and

The last inequality is given since . Suppose that . If , then

If , we have
Consequently, the solution of system (17) departing from cannot intersect . If , then . Since is equilibrium, in summary, is positive invariant.

Since , stable manifold theorem implies that there exists a one-dimensional strong unstable manifold tangent to at such that the point on near can be expressed by Furthermore, there is a two-dimensional unstable manifold tangent to span at such that near can be expressed by

Lemma 9. *Suppose that is a solution of (17) such that for small . Then, will leave and enter .*

*Proof. *Obviously, satisfies initial condition (39) by the structure of , and Lemma 8 implies ( means that and , ) for every .

Furthermore, Lemma 8 shows that , implying . Since , it follows . Suppose that for every . Then
for large since and are strictly monotonous increasing with respect to and , respectively. Thus, we have that , contradicting for any . Therefore, there exists such that . Without losing generality, let . Obviously, we have . If , thenwhich is a contradiction. Therefore, and . Then, the construction of shows that will leave and enter .

Let be a small circle on centered at . Then, points on can be expressed in terms of local coordinate by where , , and is chosen such that lies on with . Then, stable manifold theorem shows that when . Denote .

Lemma 10. *There exists a such that
**
and that
**
for .*

*Proof. *From (51), we have
where , , and since . Therefore, and imply that . Obviously, for any . However,

Then, equality , together with the last of (55), reveals ; that is, . For , the first and second equalities of (55) imply that where since and .

Let

By Lemma 10, is an arc of circle, , and the solution of (17) with initial value being the endpoint will enter since . From Lemma 9, the solution of (17) with initial value being the endpoint will enter .

##### 5.2. Traveling Wave Solution for

Lemma 11. *Let be a solution of (17) such that . If for any , then for any , where
**
and .*

*Proof. *Set . Suppose the conclusion is false; that is, . Obviously, and where

In Figure 3, we find , , , , , , and .

Since , thus . If , we have because for and . However, which is a contradiction. Therefore, . If , thencontradicting . If , then
since , contradicting . If , then which is a contradiction. In conclusion, . If , then and for any . Hence, for any , which implies that . From this contradiction we find . Because is a constant solution, we get . In summary, and . The proof is completed.

Lemma 12. *There exists a point such that the solution of (17) with initial value being will stay in for any .*

*Proof. *It is sufficient to prove . Suppose that . Firstly, we verify Conditions and of Lemma 6. Condition of Lemma 6 is valid since is closed.

Suppose , and . Then, and since . The structure of implies that , , and . By the proof of Lemma 11, we have that for . Therefore, and . Condition of Lemma 6 holds.

Lemma 6 shows that is homeomorphic to . Since
and is disconnected, we have that is disconnected, contradicting the connection of . Thus, and the proof is completed.

Lemma 13. *Let . Then, there exists a positive solution of (17) such that
*

*Proof. *By Lemma 12 there exists a point such that the solution of (17) with initial value being will stay in for any . Furthermore, Lemma 11 shows for any . Stable manifold theorem implies that for any and . Therefore, is a positive solution.

To complete the proof, it is sufficient to show that . By Lemma 11, we know that for any since remains in for all . Because , then the limit of exists; that is, and . Suppose that . The first equation of (17) shows that
for large , which implies that there is an