- 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 216913, 13 pages
Traveling Wave Solutions in a Reaction-Diffusion Epidemic Model
1College of Mathematics and Information Science, Wenzhou University, Wenzhou 325035, China
2College of Physics and Electronic Information Engineering, Wenzhou University, Wenzhou 325035, China
Received 4 February 2013; Revised 17 March 2013; Accepted 17 March 2013
Academic Editor: Anke Meyer-Baese
Copyright © 2013 Sheng Wang 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.
We investigate the traveling wave solutions in a reaction-diffusion epidemic model. The existence of the wave solutions is derived through monotone iteration of a pair of classical upper and lower solutions. The traveling wave solutions are shown to be unique and strictly monotonic. Furthermore, we determine the critical minimal wave speed.
Recently, great attention has been paid to the study of the traveling wave solutions in reaction-diffusion models [1–17]. In the sense of epidemiology, the traveling wave solutions describe the transition from a disease-free equilibrium to an endemic equilibrium; the existence and nonexistence of nontrivial traveling wave solutions indicate whether or not the disease can spread . The results contribute to predicting the developing tendency of infectious diseases, to determining the key factors of the spread of infectious disease, and to seeking the optimum strategies of preventing and controlling the spread of the infectious diseases [18–21].
Some methods have been used to derive the existence of traveling wave solutions in reaction-diffusion models, and the monotone iteration method has been proved to be an effective one. Such a method reduces the existence of traveling wave solutions to that of an ordered pair of upper-lower solutions [6, 7, 9, 10, 14, 15].
In , Berezovsky and coworkers introduced a simple epidemic model through the incorporation of variable population, disease-induced mortality, and emigration into the classical model of Kermack and McKendrick . The total population is divided into two groups of susceptible and infectious ; that is to say, . The model describing the relations between the state variables is where the reproduction of susceptible follows a logistic equation with the intrinsic growth rate and the carrying capacity , denotes the contact transmission rate (the infection rate constant), is the natural mortality; denotes the disease-induced mortality, and is the per-capita emigration rate of uninfected.
For model (1), the epidemic threshold, the so-called basic reproduction number , is then computed as . The disease will successfully invade when but will die out if . is usually a threshold whether the disease goes to extinction or goes to an endemic. Large values of may indicate the possibility of a major epidemic . In addition, the basic demographic reproductive number is given by . It can be shown that if the population grows, while implies that the population does not survive .
For simplicity, rescaling model (1) by letting , , and leads to the following model: where is defined by the ratio of the average life span of susceptibles to that of infectious.
In this paper, we are interested in the existence of traveling wave solutions in the following reaction-diffusion epidemic model : where are all positive constants, is the diffusion coefficient, and .
This paper is arranged as follows. In Section 2, we construct a pair of ordered upper-lower solutions of model (3) and establish the uniqueness and strict monotonicity of the traveling wave solutions.
2. Existence and Asymptotic Decay Rates
In this section, we will establish the existence of traveling wave solutions of model (3) by constructing a pair of ordered upper-lower solutions. The definition of the upper solution and the lower solution is standard. We assume that the inequality between two vectors throughout this paper is componentwise.
Setting then model (3) can be written as
For model (3), the equilibria are and , where and for model (7), the equilibria are and , where Obviously, For simplicity, we define the following functions and constants: And we will always assume the following hypotheses throughout the rest of this paper: [H1] [H2]
Then we can obtain the following.
Lemma 2. For model (7), if [H1] holds, then is unstable, and is stable.
For the sake of convenience, let . For simplicity, we still use the variables , , and instead of , , and , respectively, then model (7) could be rewritten as
Following the definition of quasi-monotonicity , we can obtain the following results.
Lemma 3. Model (14) is a quasi-monotone decreasing system in .
From , we can know that the functions and are said to possess a quasi-monotone nonincreasing system, if the sign of and are both nonpositive.
Since Then, Let then obviously, is the unique real root of .
Since , consider , then we can get And hence, has two positive roots.
Since , thus .
According to conditions and , we can get Then, . Hence, .
That is to say, model (14) is a quasi-monotone system in .
Obviously, we can know the following.
Remark 4. Model (24) is also a quasi-monotone system in .
Now we establish the existence of traveling wave solutions of model (24) through monotone iteration of a pair of smooth upper and lower solutions. Following , we give the definitions of the upper and lower solutions of model (24) as follows, respectively.
Definition 5. A smooth function () is an upper solution of model (24) if its derivatives and are continuous on , and satisfies with the following boundary value conditions
Definition 6. A smooth function () is a lower solution of model (24) if its derivatives and are continuous on , and satisfies with the following boundary value conditions
The construction of the smooth upper-lower solution pair is based on the solution of the following KPP equation: where and in the open interval with , , and . First, let us recall the following result.
Lemma 7 (see [1, 15]). Corresponding to every , model (29) has a unique (up to a translation of the origin) monotonically increasing traveling wave solution for . The traveling wave solution has the following asymptotic behaviors.(i) For the wave solution with noncritical speed , one has where and are positive constants.(ii) For the wave with critical speed , one has where the constant is negative, is positive, and .
For constructing the upper solution of the model (24), we start with the following model:
Define , , one can verify that all of the following conditions are satisfied:(i);(ii)for all and ;(iii), .
Define then we can get the following result.
Proof. On the boundary,
As for the component, we have
As for the component, since , then . And(i)if , then ;(ii)if , then .
Thus we can get: Hence, forms a smooth upper solution for model (24).
For constructing the lower solution of the model (24), we start with the following model:
Define , . One can easily verify that all of the following conditions hold:(i);(ii), for all and ;(iii), .
Define then we have the following result:
Proof. On the boundary, As for the component, we have As for the component, we have Thus forms a smooth lower solution for model (24).
Proof. Our proof is only for , and the proof for the case of is similar to it.
First, we derive the asymptotic behaviors of the upper solution and the lower solution at infinities.
According to Lemma 7, when , we can obtain:
And let , , when , we can get where, , , , are all positive constants.
Since for any , is also a solution of model (32). Thus, is an upper solution of model (24). So, according to Lemma 7, when , we can get:
Since , we can choose a large enough number , such that hence, there exists a large number , such that
By using a similar argument as above, there exists a large enough number , such that
Second, we show that
We deal with such two possible cases:
Case 1. If then, the proof is completed.
Case 2. If there exists a point , such that satisfying or .
In this case, we use the Sliding Domain method .
Step 1. we shift to the left by increasing the number until finding a new number such that on the smaller interval .
Step 2. we shift back to the right by decreasing to a smaller number such that one of the branches of the upper solution touches its counterpart of the lower solution at some point in the interval . On the endpoints of the interval , we still have .
Let and , where
For , we get that where , . Since the above model is monotone and the cube is convex, thus we can deduce by Maximum Principle that for . So does not exist and we can decrease further to . It is calculated that the point does not exist either. The proof of this lemma is completed.
To ease the burden of notations, we still use to denote the shifted upper solution as given in Lemma 8. Let
With such constructed ordered upper-lower solution pair, we can get the following.
Theorem 11. For , model (24) has a unique (up to a translation of the origin) traveling wave solution. The traveling wave solution is strictly increasing and has the following asymptotic properties:(i) if , when , when , and if , then while : where, , , , , , , and are all positive constants.(ii) if , when , when , and if , then while , where , , , and , , , are all positive constants.
Proof. From Lemma 3 and Remark 4, we know that model (24) is a quasi-monotone nonincreasing system in , and by using the monotone iteration scheme given in [3, 13], we can obtain the existence of the solution to the first two equations in model (24) for every , which satisfies
According to the above inequality, we can get that, on the boundary, the solution tends to as and as .
To derive the asymptotic decay rate of the traveling wave solutions as , we just let and be the traveling wave solution of model (24) generated form the monotone iteration, since the case of (ii) is similar to it.
We differentiate model (24) with respect to , and note that satisfies where
Now, we study the exponential decay rate of the traveling wave solution as . The asymptotic model of model (62) as is where
The second equation of model (64) has two independent solutions with the following form:
Relating the second equation of model (62) with the second equation of model (64), we can deduce that has the following property as : for some constants and . Thus, we can obtain that where So we obtain that Thus, .
Now, we consider the first equation of model (64). We rewrite it as
One can verify that is not a characteristic of
The above equation has two independent solutions of the following form: Thus, when , has the following property: for some constants ; . Since , thus . So, when , we have the following formula:
Then, we study the exponential decay rate of the traveling wave solution as . The asymptotic model of model (62) as is
By setting , , we rewrite model (76) as a first order model of ordinary differential equation in the four components :
In the case of (i) , we can obtain that the solution of model (77) has the following form: where and are the eigenvectors of the constant matrix with as the corresponding eigenvalues, are arbitrary constants. Since thus