Traveling Wavefronts of a Diffusive Hematopoiesis Model with Time Delay ()
1. Introduction
It is well know that the traveling wave theory was initiated in 1937 by Kolmogorov, Petrovskii, Piskunov [1] and Fisher [2] . Now, the theory of traveling wave solutions to reaction-diffusion equations is one of the fast developing areas of modern mathematics and has attracted much attention due to its significance in biology, chemistry, epidemiology and physics, see [3] [4] and the reference cited therein. In recent years, the traveling wave problem for reaction-diffusion systems with delay has been widely studied. For example, Gomez and Trofimchuk [5] considered the Fisher-KPP equation and their results showed that each monotone traveling wave could be found via an iteration procedure by using the special montone integral operators. Schaaf [6] systematically studied two scalar reaction-diffusion equations with a single discrete delay by using the phase plane technique, the maximum principle for parabolic functional differential equations and the general theory of ordinary differential equations. The degree theory has been adopted in [7] [8] .
In this paper we consider the following reaction-diffusion equation with a discrete time delay:
(1.1)
When is independent on the spatial variable, the above equation reduces to the following ordinary differential equation
(1.2)
which was first proposed by Mackey and Glass [9] to describe the dynamics of blood cell production. Here, denotes the density of mature stem cells in blood circulation and is the time delay between the production of immature stem cells in bone marrow and their maturation for release in the circulating blood stream; and are positive constants that represent some specific meanings in blood circula-
tion. For instance, is the lost rate of the cells from the circulation. The term shows that the
flux of the cells into the circulation from the stem cell compartment depends on the number of cells at time. For more details about Hematopoiesis model, we refer the readers to the articles of Mackey [9] -[11] and the references given in them.
Equation (1.2) has been studied by many authors such as in [12] -[14] . Weng and Dai [13] proved that the positive equilibrium to Equation (1.2) could be a global attractor under some conditions. Wu, Li and Zhou [14] derived a sufficient and necessary condition that guarantees the existence of positive periodic solutions of Equation (1.2) with periodic coefficients.
Equation (1.2) can be generalized as the following functional differential equation
(1.3)
Wang [15] investigated the generalized equation with Neumann boundary condition and obtained the oscillatory behavior of solutions about the positive equilibrium of (1.3). Further, they derived the sufficient and necessary conditions for global attractivity of the zero solution. In addition, global attractivity of the positive equilibrium of (1.3) was investigated by Gopalsamy and Kulenvic [16] . Cheng and Zhang [17] and Jiang et al. [18] (-dimentional case) instead investigated the existence of positive periodic solutions of Equation (1.3) by using the Krasnosel skii fixed point theorem.
The aim of this paper is to consider the existence of traveling wave front solutions for (1.1) in the case of one dimensional space.
This paper is outlined as follows. The next section, we will introduce the technique of upper and lower solutions developed by Wu and Zou [19] . The conditions for establishing the positive equilibria and obtaining the existence of traveling waves are derived in Section 2.
To investigate the existence of traveling wave fronts of (1.1), we describe briefly the technique of upper and lower solutions developed by Wu and Zhou [19] .
Consider a scalar reaction-diffusion equation with time delay:
(1.4)
where, and is the diffusion coefficient. The function is continuous and is an element in parameterized by and given by
Looking for traveling wave solutions of the form leads to a second-order functional differential Equation
(1.5)
where is define by
Now we assume that
A1. There exists such that and for, where denotes the constant function taking the value on.
A2. There exists such that
for with.
If for some, (1.5) has a monotone solution defined on such that
(1.6)
then is called a traveling wave front of (1.4) with a wave speed.
Define a profile set for traveling wave fronts of (1.1) by
The upper and lower solution for (2.1) are defined as follows:
Definition 1 The piecewise smooth functions and in are called upper and lower solution of (1.4) if and if
and satisfies the above differential inequalities in reversed order.
Now we are in the position to state a scalar version of [19] (Theorem 3.6).
Theorem 1 If the conditions (A1) and (A2) hold, suppose that (1.5) has an upper solution in and a lower solution (which is not necessarily in) with and in, then the problem (1.4) admits a traveling wave front.
2. Existence of Traveling Wave Fronts
Assume that, and we can get two equilibria of (1.1)
We will tackle the existence of solutions of (3.1) with the asymptotic boundary condition
which corresponds to the traveling wave fronts of (1.2) connecting and.
Substituting into (1.2), and denoting the movingvariable still by, the resulting wave equation becomes
(1.7)
Define the function
Lemma 1 If, then, and, for any, where denotes the constant function taking the value on.
Next we show that satisfies quasi-monotonicity condition with some assumptions.
Lemma 2 If, , then satisfies the following quasi-monotonicity condition:
Take, and we have
for all with,.
Proof 1 Consider the function, and it is obvious that
If, then
since.
It demonstrates that the function is increasing on. A direct computation shows that
and then
Therefore, if choosing, we have
.
This completes the proof.
Remark 1 If, the function is increasing for all. When the condition (A2)
is hold.
Define the profile set
Next we will discuss the existence problem by using the method of upper and lower solutions that are defined as follows:
Definition 2 The piecewise smooth functions and in are called upper and lower solution of (3.1) if and if
and satisfies the above differential inequalities in reversed order.
Define
,
then we have the following lemmas.
Lemma 3 There exists such that for, has two positive real roots, and
Since the proof of this lemma is similar to that of Claim 2.3 of [19] , we omit it. Next we first construct the upper solution of (3.1).
Lemma 4 Assume, then is an upper solution of (3.1) and.
Proof 2 It is easy to verify that. We show that is an upper solution of (3.1).
Let be such that.
i) For, , , ,. thus
ii) For, and, thus
According to the discussion above, we know that is an upper solution of (1.5). This completes the proof.
We now give the lower solution to (1.5). Let and be the same as those given in Lemma 3.2. Take such that. Define, where the constant is to be determined.
Lemma 5 For, is a lower solution for Equation (3.1).
Proof 3 Let, then for and
i) For, , and. Hence
ii) For, we have
It is easy to check that. It follows that. Therefore
Note that
Since, we have. Thus
If we choose
then. Thus is a lower solution of (3.1).
It is clear that. Summarizing the above conclusions, we give our main result of this paper below.
Theorem 2 If, , or, then for every, the problem (1.2) has a traveling wave front which connects the equilibria and.
NOTES
*Corresponding author.