Next Article in Journal
New Applications of Sălăgean and Ruscheweyh Operators for Obtaining Fuzzy Differential Subordinations
Next Article in Special Issue
Existence and Uniqueness of Nontrivial Periodic Solutions to a Discrete Switching Model
Previous Article in Journal
Two Gaussian Bridge Processes for Mapping Continuous Trait Evolution along Phylogenetic Trees
Previous Article in Special Issue
Existence Results for p1(x,·) and p2(x,·) Fractional Choquard–Kirchhoff Type Equations with Variable s(x,·)-Order
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Traveling Waves Solutions for Delayed Temporally Discrete Non-Local Reaction-Diffusion Equation

School of Mathematics and Information Science, Guangzhou University, Guangzhou 510006, China
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(16), 1999; https://doi.org/10.3390/math9161999
Submission received: 9 July 2021 / Revised: 16 August 2021 / Accepted: 18 August 2021 / Published: 20 August 2021
(This article belongs to the Special Issue Functional Differential Equations and Epidemiological Modelling)

Abstract

:
This paper deals with the existence of traveling wave solutions to a delayed temporally discrete non-local reaction diffusion equation model, which has been derived recently for a single species with age structure. When the birth function satisfies monotonic condition, we obtained the traveling wavefront by using upper and lower solution methods together with monotonic iteration techniques. Otherwise, without the monotonicity assumption for birth function, we constructed two auxiliary equations. By means of the traveling wavefronts of the auxiliary equations, using the Schauder’ fixed point theorem, we proved the existence of a traveling wave solution to the equation under consideration with speed c > c * , where c * > 0 is some constant. We found that the delayed temporally discrete non-local reaction diffusion equation possesses the dynamical consistency with its time continuous counterpart at least in the sense of the existence of traveling wave solutions.

1. Introduction

Very recently, we derived in [1] the following delayed temporally discrete non-local reaction-diffusion equation
w ( n + 1 , x ) = D Δ x w ( n + 1 , x ) + ( 1 d ) w ( n , x ) + η + K D i ( r 1 , x y ) b ( w ( n + 1 r , y ) ) d y ,
where w ( n , x ) denotes the density of the total matured population of a single species at time n N , location x R , N represents the set of non-negative integers. D > 0 is the diffusion rate of matured individuals. Δ x is the Laplacian operator with respect to x R , r 2 is a positive integer, which denotes the maturation time for the species. 0 < d < 1 represents the death rate of matured individuals. 0 < η < 1 is the survival rate of an individual from birth to maturation. b ( · ) is the birth function. The integral kernel function K D i reads as
K D i ( r 1 , x ) = 1 D i e | x | D i m = 0 r 2 ( r 1 ) · r ( r + m 2 ) m ! ( r 2 m ) ! 2 r + m 1 | x | D i r 2 m ,
where D i is the diffusion rate of immatured individuals. In particular, we use the convention that if m = 0 then m ! = 1 and ( r 1 ) r ( r + m 2 ) = 1 .
For the case r = 1 , K D i ( r 1 , x ) = K D i ( 0 , x ) = δ ( x ) , i.e., the Dulac δ function. Therefore, (1) is reduced to
w ( n + 1 , x ) = D Δ x w ( n + 1 , x ) + ( 1 d ) w ( n , x ) + η b ( w ( n , x ) ) ,
Equation (1) reflects the changes of matured population of a single species in an unbounded habitat with dimension one. It can also be viewed as a non-standard discretization of the following well-known non-local reaction-diffusion equation model with delay,
w t = D 2 w x 2 d w + ε b ( w ( t r , y ) ) f α ( x y ) d y ,
which was derived in 2001 by So, Wu and Zou [2]. The model (4) describes the adult population’s evolution of a single species that have two age classes in an unbounded spatial domain R 1 . In above Equation (4), D > 0 is the diffusion rate of the adult population and d > 0 is the death rate of the adult population. ε > 0 reflects the impact of the death rate of the immature on the matured population. b ( · ) is the birth function. r > 0 is the maturation time of the species. α 0 also reflect the impact of the dispersal rate of the immature on the matured population. The integral kernel function is f α ( x ) = 1 4 π α e x 2 4 α . When α approaches to 0, that is, the immature do not disperse, (4) is reduced to
w t = D 2 w x 2 d w + ε b ( w ( t r , x ) ) ,
and the non-local effect disappears.
In our previous paper [1], we have already given a detailed derivation of Equation (1). Especially, when D i = 0 (the immature individuals do not disperse), (1) becomes
w ( n + 1 , x ) = D Δ x w ( n + 1 , x ) + ( 1 d ) w ( n , x ) + η b ( w ( n + 1 r , x ) ) .
Above equation is a non-standard temporally discrete version of (5).
It is well known that the continuous-time reaction-diffusion equation has been widely used to describe diffusive phenomena in physics, engineering, chemistry, biology, and so on (see, for example, [3,4,5,6,7,8,9,10]). In general, the dynamical behaviors of solutions to a non-linear reaction-diffusion equation are very complicated, and it is often very difficult to find exact solutions. For the sake of understanding the properties of the solutions numerically, we need to study its discrete analogue. However, the basic principle of constructing appropriate discretization of differential equations is to preserve the properties of the corresponding differential equations. Since many classical (standard) discretizations cannot achieve dynamic consistency, non-standard discretizations are usually used to ensure it. Thus, Mickens first introduced the concept of dynamical consistency in [11,12] for ordinary differential equations and since then, some dynamical consistent discrete schemes have been constructed. See, for example, Refs. [13,14,15,16,17], and references therein.
In [1], we established the existence of traveling wavefronts of Equation (6) by using upper and lower solution methods and iterate techniques. We found that (6) possesses the dynamical consistency with its time continuous counterpart (5) at least in the sense of the existence of traveling wave solutions. In the sense of propagation, Equation (6) is also a good approximation of corresponding continuous time model (5).
From the perspective of mathematical biological modeling, almost all of the data collected are discrete in time because observations are always discontinuous. For example, satellite photographs used for scientific research are usually taken periodically, but the spatial distribution can be seen as continuous. Temporally discrete and spatially continuous diffusion model will be more suitable than its corresponding time continuous diffusion model to study the dynamic behavior of a single species that living in a spatially continuous habitat in population ecology. In 2002, Weinberg, Lewis, and Li in [18] gave some reasons for studying discrete-time models rather than just reaction diffusion models. They also pointed out the advantage of a discrete-time model over a reaction-diffusion model. We note that in [18], the authors studied the discrete-time recursion system
u n + 1 = Q [ u n ] ,
where u n ( x ) denotes the population distributions of species and Q is an operator that models the growth, migration, and interaction of the species.
Although Equations (1) and (6) are different from (7) formally, we have proved in [1] that the following general non-linear equation
w ( n + 1 , x ) = D Δ x w ( n + 1 , x ) + ( 1 d ) w ( n , x ) + f ( n , w ( n , x ) ) , x R .
is equivalent to an integral-difference equation as below,
w ( n + 1 , x ) = + k ( x y ) [ ( 1 d ) w ( n , y ) + f ( n , w ( n , y ) ) ] d y ,
where k ( x ) = 1 2 D e | x | D . Clearly, (9) is a special case of (7) with
Q ( w ( n , x ) ) = + k ( x y ) [ ( 1 d ) w ( n , y ) + f ( n , w ( n , y ) ) ] d y .
As for the general integro-difference equations,
u n + 1 ( x ) = Ω k ( n , x , y ) f ( n , u n ( y ) ) d y ,
Kot and Schaffer [19] are the first to apply it modeling temporally discrete and spatially continuous dispersal phenomena, and studying the dispersal of a single species with non-overlapping generations. They showed that, the above model will has more complex dynamic behavior than its corresponding time-continuous one. Moreover, even chaos could occur.
Since the selection of kernel function k ( n , x , y ) plays a key role in the dynamical behavior of (10), using such a model to describe some biological phenomena will have some uncertainty. Especially, when we discuss dynamical behaviors of populations of some species living in a bounded domain, the choice of suitable integral kernels is very difficult because we may cope with various different boundary value problems.
In contrast with integral difference equations, we found that there is no such problems for temporally discrete reaction diffusion equations like (8) (or (1)). Furthermore, from the point of view of the mathematical modeling, (8) (or (1)) has the same biological explanations as those for integral difference equations. In fact, in our previous paper [1], the life cycle of individuals of the population is divided into relatively sedentary and dispersal stages. This coincides with the explanations by Kot and Schaffer [19] in establishing Equation (7). To distinguish the difference of dynamical behaviors which occurs at different stages, we assume that the evolution (the relatively sedentary stage) occurs at time n and dispersal occurs at time n + 1 .
For temporally discrete reaction diffusion models, there are only a few results in the literature. In 2006, Lin and Li [20] studied following equation with delay:
u n ( x ) u n 1 ( x ) = d Δ u n ( x ) + f ( u n ( x ) , u n τ ( x ) ) , n N , x R .
They established the existence of traveling wavefronts and showed that (11) is a good approximation of its continuous time model in the sense of propagation. For more researches on this topic see [4,21,22]. However, we note that in the existing research literature, researchers simply assumed that the non-standard discretizations preserve the dynamical consistency of the continuous-time reaction-diffusion equations, but they do not provide reasonable biological explanations for the modeling process.
Although (1) has been derived in [1], the existence of traveling wave solutions is proved without non-local effect, and monotonic condition for birth function is assumed. In order to better understand the dynamical behaviors of (1) with non-local diffusion caused by immature individuals dispersion, we will study traveling wave solutions whenever the birth function is monotonic or non-monotonic, respectively.
The rest of this paper is organized as follow. In Section 2, by using the upper-lower solutions and monotone iteration technique, we studied the existence of traveling wavefronts of (1) for the cases that the birth function b ( w ) is increasing in [ 0 , w * ] , where w * > 0 is the unique solution to the equation η b ( w ) = d w . As for the case that birth function b ( w ) is non-monotone in [ 0 , w * ] , the theory of monotone dynamical system cannot be directly used. By using a similar idea as Ma in [23,24], we establish the existence of traveling waves of (1) in Section 3. Finally, we give a short discussion in Section 4.

2. Traveling Wavefronts for the Monotone Case

In this section, we will consider the existence of traveling wavefronts to equation (1) for monotone case. We are interested in finding traveling waves w ( n , x ) = φ ( x + c n ) of following equation
w ( n + 1 , x ) = D Δ w ( n + 1 , x ) + ( 1 d ) w ( n , x ) + η + K D i ( r 1 , x y ) b ( w ( n + 1 r , y ) ) d y .
For this purpose, we will find a solution φ ( ξ ) to (12) with ξ = x + c n . Clearly, φ ( ξ ) satisfies the following wave profile equation
φ ( ξ + c ) = D φ ( ξ + c ) + ( 1 d ) φ ( ξ ) + η + b ( φ ( ξ + c y c r ) ) K D i ( r 1 , y ) d y .
Let ξ = ξ + c and still denote it by ξ . Then, the above equation becomes
D φ ( ξ ) + φ ( ξ ) = ( 1 d ) φ ( ξ c ) + η + b ( φ ( ξ y c r ) ) K D i ( r 1 , y ) d y
Throughout this section, we always assume that
Hypothesis 1 (H1).
b ( w ) = w g ( w ) , where g ( w ) is a continuously differentiable function and satisfying g ( w ) > 0 , g ( w ) < 0 for w 0 and lim w g ( w ) 0 ;
Hypothesis 2 (H2).
b ( w ) and b ( w ) are bounded;
Hypothesis 3 (H3).
d < b ( 0 ) η = g ( 0 ) η .
From (H1) and (H3), Equation (12) has only two constant equilibria w = 0 and w = w * , where w * is the unique solution of the equation η g ( w ) = d .
We will study the existence of non-decreasing solutions to Equation (14) subject to the boundary value conditions
lim ξ φ ( ξ ) = 0 , lim ξ + φ ( ξ ) = w * .
Our approach is similar to that in [1], which is based on the monotonic iteration techniques combined with upper and lower solution methods that was developed in [2]. To this end, we further assume that
Hypothesis 4 (H4).
b ( w ) is increasing in [ 0 , w * ] .
To proceed further, for readers’ convenience, we introduce some results on the following temporally discrete reaction-diffusion equation which will be used later,
u ( n + 1 , x ) u ( n , x ) = D Δ x u ( n + 1 , x ) + f ( u ( n , x ) ) ,
where f ( · ) is assumed to be a continuous function. For detailed proofs on these results, readers can refer to [1].
Lemma 1
([1]). The initial value problem
u ( n + 1 , x ) u ( n , x ) = D Δ x u ( n + 1 , x ) , x R , n 0 . u ( 0 , x ) = φ ( x )
has a unique solution u ( n , x ) and
u ( n , x ) = K ( n , x y ) φ ( y ) d y , n 1
where φ is a continuous and absolutely integrable function on R 1 satisfying lim | x | φ ( x ) = 0 .
K ( n , x ) = e | x | D 1 D m = 0 n 1 n ( n + 1 ) ( n + m 1 ) m ! ( n 1 m ) ! 2 n + m | x | D n 1 m .
In particular, we use the convention that if m = 0 , then m ! = 1 and n ( n + 1 ) ( n + m 1 ) = 1 . The integral kernel function K ( n , x ) satisfies the following properties.
Proposition 1
([1]). 
1.
+ K ( n , x ) d x = 1 , for n > 0 ;
2.
K ( 0 , x ) = δ ( x ) , for n = 0 ;
3.
For x 0 , n > 0 , K ( n , x ) is the solution of equation
u ( n + 1 , x ) u ( n , x ) = D Δ x u ( n + 1 , x ) .
Lemma 2
([1]). The function u ( n , x ) is a solution of (15) if, and only if, it is a solution of following non-linear integro-difference equation
u ( n + 1 , x ) = + 1 2 D e | x y | D [ u ( n , y ) + f ( u ( n , y ) ) ] d y , n 0 .
For φ C ( R , R ) , define
H ( φ ) ( ξ ) = ( 1 d ) φ ( ξ c ) + η + b ( φ ( ξ y c r ) ) K D i ( r 1 , y ) d y .
By (H2) and Proposition 1, H ( φ ) ( ξ ) is well defined. Then, Equation (14) becomes
D φ ( ξ ) + φ ( ξ ) = H ( φ ) ( ξ ) .
Let Γ = { φ C ( R , R ) | φ is non decerasing in ξ R , lim ξ φ ( ξ ) = 0 , lim ξ + φ ( ξ ) = w * } .
For φ , ψ Γ , for ξ R , we say ψ φ , if ψ ( ξ ) φ ( ξ ) . Then from the definition of H and Proposition 1, we can easily obtained the following result.
Lemma 3.
Assume that (H4) holds. Then, for every φ Γ , H ( φ ) Γ , and H ( φ ) is increasing on Γ, that is, H ( ψ ) H ( φ ) for φ , ψ Γ with ψ φ .
By direct computations, we have
Lemma 4.
For h Γ , the equation
D φ ( ξ ) + φ ( ξ ) = h ( ξ ) .
has a unique solution φ Γ satisfying
φ ( ξ ) = 1 2 D + e | ξ ζ | D h ( ζ ) d ζ = + K ( 1 , ξ ζ ) h ( ζ ) d ζ ,
where K ( 1 , x ) is defined in (18).
Further define the mapping F : Γ Γ ,
F ( ψ ) ( ξ ) = + K ( 1 , ξ ζ ) H ( ψ ) ( ζ ) d ζ , ψ Γ
The following theorem is a direct consequence of Lemmas 3 and 4.
Theorem 1.
φ Γ is a solution of Equation (14) if, and only if, it satisfies
φ ( ξ ) = + K ( 1 , ξ ζ ) ( 1 d ) φ ( ζ c ) + η + b ( φ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ .
In other words, φ is a fixed point of mapping F in Γ.
Above theorem shows that the existence of traveling wavefronts of Equation (12) is equivalent to the existence of fixed points to mapping F . Therefore, in what follows, we will construct upper and lower solutions of (22) to prove that the mapping F has a unique fixed point in Γ , as long as c is greater than a certain constant c * .
The linearized equation of (14) at φ = 0 is as below.
D φ ( ξ ) + φ ( ξ ) = ( 1 d ) φ ( ξ c ) + η + g ( 0 ) φ ( ξ y c r ) ) K D i ( r 1 , y ) d y .
Then, its characteristic equation is given by
D λ 2 + 1 = ( 1 d ) e λ c + g ( 0 ) η + K D i ( r 1 , y ) e λ ( y + c r ) d y .
Define
Δ ( c , λ ) = ( 1 d ) e λ c + D λ 2 1 + g ( 0 ) η + K D i ( r 1 , y ) e λ ( y + c r ) d y .
Note that for | λ | < 1 D i ,
+ K D i ( r 1 , y ) e λ y d y = + 1 D i e | y | D i m = 0 r 2 ( r 1 ) · r ( r + m 2 ) m ! ( r 2 m ) ! 2 r + m 1 | y | D i r 2 m e λ y d y = m = 0 r 2 ( r 1 ) · r ( r + m 2 ) m ! ( r 2 m ) ! 2 r + m 1 0 + 1 D i e ( λ + 1 D i ) y y D i r 2 m d y + m = 0 r 2 ( r 1 ) · r ( r + m 2 ) m ! ( r 2 m ) ! 2 r + m 1 0 1 D i e ( λ 1 D i ) y y D i r 2 m d y = m = 0 r 2 ( r 1 ) · r ( r + m 2 ) m ! ( r 2 m ) ! 2 r + m 1 1 ( 1 + λ D i ) r 1 m 0 + e z z r 2 m d z + m = 0 r 2 ( r 1 ) · r ( r + m 2 ) m ! ( r 2 m ) ! 2 r + m 1 1 ( 1 λ D i ) r 1 m 0 + e z z r 2 m d z .
Since 0 + e z z r 2 m d z = ( r 2 m ) ! , we have
+ K D i ( r 1 , y ) e λ y d y = m = 0 r 2 ( r 1 ) · r ( r + m 2 ) m ! 2 r + m 1 1 ( 1 + λ D i ) r 1 m + 1 ( 1 λ D i ) r 1 m = m = 0 r 2 ( r 1 ) · r ( r + m 2 ) m ! 2 r + m 1 · ( 1 + λ D i ) r 1 m + ( 1 λ D i ) r 1 m ( 1 λ 2 D i ) r 1 m .
By direct computations, we have
m = 0 r 2 ( r 1 ) · r ( r + m 2 ) m ! 2 r + m 1 · ( 1 + λ D i ) r 1 m + ( 1 λ D i ) r 1 m ( 1 λ 2 D i ) r 1 m = 1 ( 1 λ 2 D i ) r 1 .
That is,
+ K D i ( r 1 , y ) e λ y d y = 1 ( 1 λ 2 D i ) r 1 .
In fact, one can give another proof of the equality as follows. By Lemma 1, for | λ | < 1 D i , the initial value problem
u ( n + 1 , x ) u ( n , x ) = D i Δ x u ( n + 1 , x ) , x R , n 0 . u ( 0 , x ) = e λ x
has a unique solution u ( n , x ) and
u ( n , x ) = K D i ( n , x y ) e λ y d y , n 1
On the other hand, one can easily verified that
u ( n , x ) = e λ x ( 1 λ 2 D i ) n
is a solution to initial value problem (26). Therefore,
K D i ( n , x y ) e λ y d y = e λ x ( 1 λ 2 D i ) n ,
which leads to the equality (25). Then
Δ ( c , λ ) = ( 1 d ) e λ c + D λ 2 1 + g ( 0 ) η e λ c r ( 1 λ 2 D i ) r 1 ,
and
Δ λ ( c , λ ) = c ( 1 d ) e λ c + 2 D λ + g ( 0 ) η e λ c r [ 2 λ ( r 1 ) D i c r ( 1 λ 2 D i ) ] ( 1 λ 2 D i ) r ,
Δ λ ( c , λ ) = c 2 ( 1 d ) e λ c + 2 D + g ( 0 ) η e λ c r ( 1 λ 2 D i ) r + 1 2 λ ( r 1 ) D i c r ( 1 λ 2 D i ) 2 + 2 ( r 1 ) D i ( 1 λ 2 D i ) + 4 λ 2 D i 2 ( r 1 ) .
It is easy to see that for any given c 0 and | λ | < 1 D i , Δ λ ( c , λ ) > 0 . Then, for any given c 0 , Δ λ ( c , λ ) strictly increasing in λ ( 0 , 1 D i ) . Additionally, from
Δ λ ( c , 0 ) = c ( 1 d ) g ( 0 ) η c r < 0 , c > 0
and lim λ 1 D i 0 Δ λ ( c , λ ) = + , for any c > 0 , there exists unique λ = λ ( c ) , 0 < λ ( c ) < 1 D i , such that Δ λ ( c , λ ( c ) ) = 0 and
Δ m i n = Δ ( c , λ ( c ) ) = ( 1 d ) e λ ( c ) c + D λ 2 ( c ) 1 + g ( 0 ) η e λ ( c ) c r ( 1 λ 2 ( c ) D i ) r 1 .
If Δ ( c , λ ( c ) ) > 0 , the characteristic Equation (24) has no positive real roots; If Δ ( c , λ ( c ) ) = 0 , the characteristic Equation (24) has only one positive real root; If Δ ( c , λ ( c ) ) < 0 , the characteristic Equation (24) has exactly two positive real roots.
Lemma 5.
There exist a positive constant c * > 0 , such that
(i) 
If c > c * , then the characteristic Equation (24) has two positive real roots, 0 < λ 1 ( c ) < λ 2 ( c ) ;
(ii) 
If c < c * , then the characteristic Equation (24) has no positive real roots;
(iii) 
If c = c * , the characteristic Equation (24) has only one positive real roots.
Proof. 
Consider the function Δ ( c , λ ) = ( 1 d ) e λ c + D λ 2 1 + g ( 0 ) η e λ c r ( 1 λ 2 D i ) r 1 . When c = 0 and 0 λ < 1 D i , due to g ( 0 ) η > d , we have
Δ ( 0 , λ ) = ( 1 d ) + D λ 2 1 + g ( 0 ) η ( 1 λ 2 D i ) r 1 > D λ 2 d + g ( 0 ) η > g ( 0 ) η d > 0 ,
which implies that (24) will has no positive real roots.
By the continuity of Δ ( c , λ ) in c, for sufficiently small c > 0 , (24) will also have no positive real roots. Otherwise, there will be sequences { c n } and λ n with c n > 0 , lim n c n = 0 and 0 < λ n < 1 D i , such that Δ ( c n , λ n ) = 0 , i.e.,
( 1 d ) e λ n c n + D λ n 2 1 + g ( 0 ) η e λ n c n r ( 1 λ n 2 D i ) r 1 = 0 .
Clearly, there exists a convergent subsequence { λ n k } , and let lim k λ n k = λ * . If λ * [ 0 , 1 D i ) , then,
0 = ( 1 d ) + D λ * 2 1 + g ( 0 ) η ( 1 λ * 2 D i ) r 1 > D λ * 2 + g ( 0 ) η d > 0 .
Contradiction. If λ * = 1 D i , then
0 = ( 1 d ) + D D i 1 + lim k g ( 0 ) η ( 1 λ n k 2 D i ) r 1 = .
It is also a contradiction.
Next, we consider Δ c ( c , λ ) for any given λ ( 0 , min { 1 D , 1 D i } ) ,
Δ c ( c , λ ) = ( 1 d ) λ e λ c λ r g ( 0 ) η e λ c r ( 1 λ 2 D i ) r 1 < 0 .
Since Δ ( 0 , λ ) > 0 , lim c + Δ ( c , λ ) = D λ 2 1 < 0 , then there exists unique c = c ( λ ) , such that Δ ( c ( λ ) , λ ) = 0 . Denote c * = min { c ( λ ) | λ ( 0 , min { 1 D , 1 D i } ) } . Obviously c * > 0 . Moreover, Δ ( c * , λ ) = 0 has exactly one positive real root and Δ ( c * , λ ) attains its minimum at the point λ = λ ( c * ) . This completes the proof. □
In what follows, we will consider the case where c > c * . In this case, the characteristic Equation (24) has two positive real roots λ 1 ( c ) and λ 2 ( c ) satisfying 0 < λ 1 ( c ) < λ 2 ( c ) < min { 1 D , 1 D i } and Δ ( c , λ ) < 0 for λ ( λ 1 ( c ) , λ 2 ( c ) ) .
Definition 1.
A continuous bounded function ϕ : R [ 0 , w * ] is called an upper solution of (22) if it satisfies
ϕ ( ξ ) + K ( 1 , ξ ζ ) ( 1 d ) ϕ ( ζ c ) + η + b ( ϕ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ ,
A lower solution of (22) is defined in a similar way by reversing the inequality in (29).
Now fixed c > c * and let 0 < λ 1 ( c ) < λ 2 ( c ) < min { 1 D , 1 D i } . Choose sufficiently large constant q > 1 to be determined later, define the functions ϕ ¯ and ϕ ̲ by
ϕ ¯ ( ξ ) = min { w * , w * e λ 1 ( c ) ξ } , ϕ ̲ ( ξ ) = max { 0 , w * [ e λ 1 ( c ) ξ q e p λ 1 ( c ) ξ ] } ,
where p ( 1 , min { 2 , λ 2 ( c ) λ 1 ( c ) } ) , w * satisfies η g ( w * ) = d .
Proposition 2.
ϕ ¯ ( ξ ) is an upper solution of (22).
Proof. 
By the definition of ϕ ¯ ( ξ ) , for any ξ R , ϕ ¯ ( ξ ) w * . Thus, by (H4) and Proposition 1, we achieve
+ K ( 1 , ξ ζ ) ( 1 d ) ϕ ¯ ( ζ c ) + η + b ( ϕ ¯ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ + K ( 1 , ξ ζ ) [ ( 1 d ) w * + η b ( w * ) ] d ζ w * + K ( 1 , ξ ζ ) d ζ = w * .
For ξ R , if ϕ ¯ ( ξ ) = w * , clearly, we have,
+ K ( 1 , ξ ζ ) ( 1 d ) ϕ ¯ ( ζ c ) + η + b ( ϕ ¯ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ ϕ ¯ ( ξ ) .
If ϕ ¯ ( ξ ) = w * e λ 1 ( c ) ξ , that is, ξ 0 , by (H1) and the definition of ϕ ¯ ( ξ ) ,
+ K ( 1 , ξ ζ ) ( 1 d ) ϕ ¯ ( ζ c ) + η + b ( ϕ ¯ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ + K ( 1 , ξ ζ ) ( 1 d ) ϕ ¯ ( ζ c ) + η g ( 0 ) + ϕ ¯ ( ζ y c r ) K D i ( r 1 , y ) d y d ζ + K ( 1 , ξ ζ ) ( 1 d ) w * e λ 1 ( c ) ( ζ c ) + η g ( 0 ) + w * e λ 1 ( c ) ( ζ y c r ) K D i ( r 1 , y ) d y d ζ = w * + K ( 1 , ξ ζ ) e λ 1 ( c ) ζ ( 1 d ) e λ 1 ( c ) c + η g ( 0 ) + e λ 1 ( c ) ( y + c r ) K D i ( r 1 , y ) d y d ζ = w * [ 1 D λ 1 2 ( c ) ] + K ( 1 , ξ ζ ) e λ 1 ( c ) ζ d ζ .
The last equality follows from the fact that λ 1 ( c ) is a root of Δ ( c , λ ) . Calculate above integral directly,
+ K ( 1 , ξ ζ ) e λ 1 ( c ) ζ d ζ = 1 2 D + e | ξ ζ | D e λ 1 ( c ) ζ d ζ = 1 2 D ξ + e ζ ξ D e λ 1 ( c ) ζ d ζ + 1 2 D ξ e ξ ζ D e λ 1 ( c ) ζ d ζ = e ξ D 2 D [ λ 1 ( c ) 1 D ] e [ λ 1 ( c ) 1 D ] ξ + e ξ D 2 D [ λ 1 ( c ) + 1 D ] e [ λ 1 ( c ) + 1 D ] ξ = 1 1 D λ 1 2 ( c ) e λ 1 ( c ) ξ .
It follows that,
+ K ( 1 , ξ ζ ) ( 1 d ) ϕ ¯ ( ζ c ) + η + b ( ϕ ¯ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ ϕ ¯ ( ξ ) .
This proves that ϕ ¯ ( ξ ) is an upper solution of (22). □
Proposition 3.
For any given p ( 1 , min { 2 , λ 2 ( c ) λ 1 ( c ) } ) , there exist a sufficiently large positive constant q > 1 , such that ϕ ̲ ( ξ ) is a lower solution of (22).
Proof. 
Let
κ ( ξ ) = w * [ e λ 1 ( c ) ξ q e p λ 1 ( c ) ξ ] .
Then
κ ( ξ ) = w * λ 1 ( c ) [ e λ 1 ( c ) ξ q p e p λ 1 ( c ) ξ ] .
There must exist ξ 0 = ln q λ 1 ( c ) ( 1 p ) < 0 , such that κ ( ξ 0 ) = 0 and ξ 1 = ln p q λ 1 ( c ) ( 1 p ) < ξ 0 , such that κ ( ξ 1 ) = 0 . We know that κ is increasing on ( , ξ 1 ) and decreasing on ( ξ 1 , + ) , and
κ ( ξ 1 ) = w * p 1 p 1 ( p q ) 1 / ( p 1 ) ( 0 , w * ) .
Consequently, 0 ϕ ̲ ( ξ ) < w * , for ξ R .
If ϕ ̲ ( ξ ) = 0 ,
+ K ( 1 , ξ ζ ) ( 1 d ) ϕ ̲ ( ζ c ) + η + b ( ϕ ̲ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ 0 = ϕ ̲ ( ξ ) .
If ϕ ̲ ( ξ ) = w * [ e λ 1 ( c ) ξ q e p λ 1 ( c ) ξ ] ,
+ K ( 1 , ξ ζ ) ( 1 d ) ϕ ̲ ( ζ c ) + η + b ( ϕ ̲ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ = + K ( 1 , ξ ζ ) ( 1 d ) ϕ ̲ ( ζ c ) + η + ϕ ̲ ( ζ y c r ) g ( ϕ ̲ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ = + K ( 1 , ξ ζ ) ( 1 d ) ϕ ̲ ( ζ c ) + η + ϕ ̲ ( ζ y c r ) g ( 0 ) K D i ( r 1 , y ) d y d ζ + η + K ( 1 , ξ ζ ) + ϕ ̲ ( ζ y c r ) [ g ( ϕ ̲ ( ζ y c r ) ) g ( 0 ) ] K D i ( r 1 , y ) d y d ζ : = I 1 + I 2 ,
where
I 1 = + K ( 1 , ξ ζ ) ( 1 d ) ϕ ̲ ( ζ c ) + η + ϕ ̲ ( ζ y c r ) g ( 0 ) K D i ( r 1 , y ) d y d ζ ,
I 2 = η + K ( 1 , ξ ζ ) + ϕ ̲ ( ζ y c r ) [ g ( ϕ ̲ ( ζ y c r ) ) g ( 0 ) ] K D i ( r 1 , y ) d y d ζ .
For I 1 , we have
I 1 + K ( 1 , ξ ζ ) ( 1 d ) w * [ e λ 1 ( c ) ( ζ c ) q e p λ 1 ( c ) ( ζ c ) ] + η g ( 0 ) + w * [ e λ 1 ( c ) ( ζ y c r ) q e p λ 1 ( c ) ( ζ y c r ) ] K D i ( r 1 , y ) d y d ζ = w * + K ( 1 , ξ ζ ) ( 1 d ) e λ 1 ( c ) ( ζ c ) + g ( 0 ) η + e λ 1 ( c ) ( ζ y c r ) K D i ( r 1 , y ) d y d ζ w * q + K ( 1 , ξ ζ ) ( 1 d ) e p λ 1 ( c ) ( ζ c ) + g ( 0 ) η + e p λ 1 ( c ) ( ζ y c r ) K D i ( r 1 , y ) d y d ζ = w * ( 1 d ) e λ 1 ( c ) c + g ( 0 ) η + e λ 1 ( c ) ( y + c r ) K D i ( r 1 , y ) d y + K ( 1 , ξ ζ ) e λ 1 ( c ) ζ d ζ w * q ( 1 d ) e p λ 1 ( c ) c + g ( 0 ) η + e p λ 1 ( c ) ( y + c r ) K D i ( r 1 , y ) d y + K ( 1 , ξ ζ ) e p λ 1 ( c ) ζ d ζ .
Note that,
Δ ( c , λ 1 ( c ) ) = ( 1 d ) e λ 1 ( c ) c + g ( 0 ) η + e λ 1 ( c ) ( y + c r ) K D i ( r 1 , y ) d y + D λ 1 2 ( c ) 1 = 0 ,
Δ ( c , p λ 1 ( c ) ) = ( 1 d ) e p λ 1 ( c ) c + g ( 0 ) η + e p λ 1 ( c ) ( y + c r ) K D i ( r 1 , y ) d y + D p 2 λ 1 2 ( c ) 1 .
Therefore,
I 1 w * [ 1 D λ 1 2 ( c ) ] + K ( 1 , ξ ζ ) e λ 1 ( c ) ζ d ζ w * q [ Δ ( c , p λ 1 ( c ) ) + 1 D p 2 λ 1 2 ( c ) ] + K ( 1 , ξ ζ ) e p λ 1 ( c ) ζ d ζ .
From the proof of Proposition 2, we know
+ K ( 1 , ξ ζ ) e λ 1 ( c ) ζ d ζ = 1 1 D λ 1 2 ( c ) e λ 1 ( c ) ξ ,
and similar arguments have
+ K ( 1 , ξ ζ ) e p λ 1 ( c ) ζ d ζ = 1 1 D p 2 λ 1 2 ( c ) e p λ 1 ( c ) ξ .
Therefore,
I 1 w * [ e λ 1 ( c ) ξ q e p λ 1 ( c ) ξ ] w * q Δ ( c , p λ 1 ( c ) ) 1 1 D p 2 λ 1 2 ( c ) e p λ 1 ( c ) ξ ϕ ̲ ( ξ ) w * q Δ ( c , p λ 1 ( c ) ) 1 D p 2 λ 1 2 ( c ) e p λ 1 ( c ) ξ .
Since λ 1 ( c ) < p λ 1 ( c ) < λ 2 ( c ) < 1 D , Δ ( c , p λ 1 ( c ) ) < 0 , we have
w * q Δ ( c , p λ 1 ( c ) ) 1 D p 2 λ 1 2 ( c ) > 0 .
For I 2 , by (H1), we know g ( w ) 0 . Let L g ( w ) 0 , w [ 0 , w * ] . Then
I 2 = η + K ( 1 , ξ ζ ) + ϕ ̲ ( ζ y c r ) [ g ( ϕ ̲ ( ζ y c r ) ) g ( 0 ) ] K D i ( r 1 , y ) d y d ζ = η + K ( 1 , ξ ζ ) + ϕ ̲ 2 ( ζ y c r ) g ( θ ϕ ̲ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ ( θ ( 0 , 1 ) ) L η + K ( 1 , ξ ζ ) + ϕ ̲ 2 ( ζ y c r ) K D i ( r 1 , y ) d y d ζ .
According to the definition of ξ 0 ,
+ K ( 1 , ξ ζ ) + ϕ ̲ 2 ( ζ y c r ) K D i ( r 1 , y ) d y d ζ = ( w * ) 2 + K ( 1 , ξ ζ ) ζ c r ξ 0 + [ e λ 1 ( c ) ( ζ y c r ) q e p λ 1 ( c ) ( ζ y c r ) ] 2 K D i ( r 1 , y ) d y d ζ ( w * ) 2 + K ( 1 , ξ ζ ) ζ c r ξ 0 + e 2 λ 1 ( c ) ( ζ y c r ) K D i ( r 1 , y ) d y d ζ ( w * ) 2 + K ( 1 , ξ ζ ) ζ c r ξ 0 + e p λ 1 ( c ) ( ζ y c r ) K D i ( r 1 , y ) d y d ζ = ( w * ) 2 + K ( 1 , ξ ζ ) e p λ 1 ( c ) ( ζ c r ) ζ c r ξ 0 + e p λ 1 ( c ) y K D i ( r 1 , y ) d y d ζ ( w * ) 2 + K ( 1 , ξ ζ ) e p λ 1 ( c ) ( ζ c r ) + e p λ 1 ( c ) y K D i ( r 1 , y ) d y d ζ .
Due to (25),
+ e p λ 1 ( c ) y K D i ( r 1 , y ) d y = 1 ( 1 p 2 λ 1 2 ( c ) D i ) r 1 .
Then
+ K ( 1 , ξ ζ ) + ϕ ̲ 2 ( ζ y c r ) K D i ( r 1 , y ) d y d ζ ( w * ) 2 ( 1 p 2 λ 1 2 ( c ) D i ) r 1 + K ( 1 , ξ ζ ) e p λ 1 ( c ) ( ζ c r ) d ζ ( w * ) 2 ( 1 p 2 λ 1 2 ( c ) D i ) r 1 e p λ 1 ( c ) c r + K ( 1 , ξ ζ ) e p λ 1 ( c ) ζ d ζ = ( w * ) 2 ( 1 p 2 λ 1 2 ( c ) D i ) r 1 e p λ 1 ( c ) c r 1 1 D p 2 λ 1 2 ( c ) e p λ 1 ( c ) ξ .
Thus,
I 2 L η ( w * ) 2 e p λ 1 ( c ) ( ξ c r ) 1 p 2 λ 1 2 ( c ) D i r 1 1 D p 2 λ 1 2 ( c ) .
Based on the above results,
I 1 + I 2 ϕ ̲ ( ξ ) w * q Δ ( c , p λ 1 ( c ) ) 1 D p 2 λ 1 2 ( c ) e p λ 1 ( c ) ξ L η ( w * ) 2 e p λ 1 ( c ) ( ξ c r ) 1 p 2 λ 1 2 ( c ) D i r 1 1 D p 2 λ 1 2 ( c ) = ϕ ̲ ( ξ ) + w * e p λ 1 ( c ) ξ 1 D p 2 λ 1 2 ( c ) q Δ ( c , p λ 1 ( c ) ) L η w * e p λ 1 ( c ) c r 1 p 2 λ 1 2 ( c ) D i r 1 .
Since Δ ( c , p λ 1 ( c ) ) < 0 , when q > 1 sufficiently large,
q Δ ( c , p λ 1 ( c ) ) L η w * e p λ 1 ( c ) c r 1 p 2 λ 1 2 ( c ) D i r 1 > 0 .
Hence, for such q,
I 1 + I 2 ϕ ̲ ( ξ ) .
That is, ϕ ̲ ( ξ ) is a lower solution of (22). This completes the proof. □
We have already obtained an upper and a lower solution of (22). Using the classical upper and lower solution method together with iteration techniques, we find the following existence result.
Theorem 2.
Assume that (H1)–(H4) hold. Then for any c > c * , (12) admits a traveling wavefront solution w ( n , x ) = φ ( x + c n ) connecting the equilibrium 0 and w * .

3. Traveling Wave Solution for the Non-Monotone Case

This section is devoted to the existence of traveling wave solutions to (1) when the birth function b ( w ) is not increasing in [ 0 , w * ] . Our approach is to construct two auxiliary temporally discrete diffusion equations with birth functions satisfying monotonic conditions. Then, by using a similar method developed in [23] and applying the results obtained in the Section 1 to these auxiliary equations, we can prove that (1) possesses a traveling wave solution connecting equilibrium 0 and w * .
Throughout of this section, we suppose that b ( w ) satisfies assumptions (H1)–(H3).
By (H3) and the continuity of g ( w ) at zero, there exists a small constant δ > 0 , such that for any w ( 0 , δ ) , g ( w ) > d η . Then, for any w ( 0 , δ ) , b ( w ) d η w = w ( g ( w ) d η ) > 0 . For such w, max { b ( z ) | 0 z w } > d η w .
On the other hand, when w > w * , g ( w ) < d η . This means that for w > w * , b ( w ) d η w < 0 . Then
max { b ( z ) | w * < z w } < d η w .
Let W * = η d max { b ( w ) | 0 w w * } . Then W * w * . For any w ( w * , W * ] , b ( w ) d η w d η W * . Consequently,
d W * η max { b ( w ) | 0 w W * }
Let W * = η d min { b ( w ) | w * w W * } . Then W * > 0 . Clearly, 0 < W * w * and b ( w ) > d w η for all w ( 0 , W * ) .
We note that, if b ( w * ) = max { b ( w ) | 0 w w * } , then W * = W * = w * . For sufficiently small ε < W * , we define two auxiliary functions b * ( w ) and b ϵ ( w ) as follows.
b * ( w ) = b ( 0 ) w , for w [ 0 , d b ( 0 ) η W * ] , d W * η , for w ( d b ( 0 ) η W * , W * ] , max { d W * η , b ( w ) } , for w > W * ,
and
b ε ( w ) = inf s [ w , W * ] { b ( s ) , d ( W * ε ) η } , for w [ 0 , W * ] , min { b ( w ) , d ( W * ε ) η } , for w > W * .
Lemma 6.
The following statements hold true:
(i) 
b * and b ε are continuous on [ 0 , + ) and non-decreasing on [ 0 , W * ] ;
(ii) 
b ( 0 ) w b * ( w ) b ( w ) b ε ( w ) > 0 for all w > 0 ;
(iii) 
b * ( 0 ) = 0 , b * ( W * ) = d η W * and b * ( w ) > d η w for all w ( 0 , W * ) ;
(iv) 
b ε ( 0 ) = 0 , b ε ( W * ε ) = d η ( W * ε ) and b ε ( w ) > d η w for all w ( 0 , W * ε ) ;
(v) 
( b * ) ( 0 ) = b ε ( 0 ) = b ( 0 ) .
The proof of Lemma is a direct verification, we omit it.
Now we consider the following two auxiliary temporally discrete diffusion equations
w ( n + 1 , x ) = D Δ w ( n + 1 , x ) + ( 1 d ) w ( n , x ) + η + K D i ( r 1 , x y ) b * ( w ( n + 1 r , y ) ) d y
and
w ( n + 1 , x ) = D Δ w ( n + 1 , x ) + ( 1 d ) w ( n , x ) + η + K D i ( r 1 , x y ) b ε ( w ( n + 1 r , y ) ) d y .
Then, the corresponding wave equations of (30) and (31) read as
D φ ( ξ ) + φ ( ξ ) = ( 1 d ) φ ( ξ c ) + η + K D i ( r 1 , y ) b * ( φ ( ξ y c r ) ) d y
and
D φ ( ξ ) + φ ( ξ ) = ( 1 d ) φ ( ξ c ) + η + K D i ( r 1 , y ) b ε ( φ ( ξ y c r ) ) d y ,
respectively. They are equivalent to the following two integral equations:
φ ( ξ ) = + K ( 1 , ξ ζ ) ( 1 d ) φ ( ζ c ) + η + b * ( φ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ ,
φ ( ξ ) = + K ( 1 , ξ ζ ) ( 1 d ) φ ( ζ c ) + η + b ε ( φ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ .
Moreover, the traveling wavefronts of (30) and (31) are fixed points of operators F * and F ε , where
F * ( φ ) ( ξ ) = + K ( 1 , ξ ζ ) ( 1 d ) φ ( ζ c ) + η + b * ( φ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ ,
F ε ( φ ) ( ξ ) = + K ( 1 , ξ ζ ) ( 1 d ) φ ( ζ c ) + η + b ε ( φ ( ζ y c r ) ) K D i ( r 1 , y ) d y d ζ .
By Lemma 6, it is easily to verify that b * and b ε satisfy all the assumptions (H1)–(H4). Therefore, the results below follow from Theorem 2.
Lemma 7.
Assume that (H1)–(H3) hold. Then for each c > c * , there exist traveling wavefronts φ * ( x + c n ) and φ ε ( x + c n ) of (30) and (31), respectively, which satisfy φ * ( + ) = W * , φ * ( ) = 0 and φ ε ( + ) = W * ε , φ ε ( ) = 0 .
In the sequel, we will always assume c > c * . Therefore, there exist two positive roots λ 1 ( c ) and λ 2 ( c ) with 0 < λ 1 ( c ) < λ 2 ( c ) < min { 1 D i , 1 D } . In order to proceed further, we also need the following assumptions.
Hypothesis 5 (H5).
b ( w ) η < d ( 2 w * w ) , w [ W * , w * ) , b ( w ) η > d ( 2 w * w ) , w ( w * , W * ] .
Let us first define
C [ 0 , W * ] = { φ C ( R , R ) : 0 φ ( ξ ) W * , ξ R } ,
and for given λ ( 0 , λ 1 ( c ) ) , denote
X λ = { φ C ( R , R ) : sup ξ R | φ ( ξ ) | e λ ξ < } ,
φ λ = sup ξ R | φ ( ξ ) | e λ ξ .
Clearly, ( X λ , · λ ) is a Banach space. Then, define
H * ( φ ) ( ξ ) = ( 1 d ) φ ( ξ c ) + η + K D i ( r 1 , y ) b * ( φ ( ξ y c r ) ) d y , ξ R ,
and
H ε ( φ ) ( ξ ) = ( 1 d ) φ ( ξ c ) + η + K D i ( r 1 , y ) b ε ( φ ( ξ y c r ) ) d y , ξ R .
For any φ , ψ ( R , [ 0 , W * ] ) with φ ( ξ ) ψ ( ξ ) , ξ R , we have
H * ( φ ) ( ξ ) H * ( ψ ) ( ξ ) and H ε ( φ ) ( ξ ) H ε ( ψ ) ( ξ ) for all ξ R .
Set
Γ * = φ C ( R , [ 0 , W * ] ) ( i ) φ ϵ ( ξ ) φ ( ξ ) φ * ( ξ ) for all ξ R ( ii ) | φ ( ξ 1 ) φ ( ξ 2 ) | 2 W * D | ξ 1 ξ 2 | for all ξ 1 , ξ 2 R .
Lemma 8.
Γ * X λ is nonempty, convex and compact in X λ .
Proof. 
Firstly, we note that φ * ( ξ ) and φ ε ( ξ ) are traveling wavefronts of (32) and (33), respectively, and they are obtained by iteration procedures. Set
ϕ * ¯ ( ξ ) = min { W * , W * e λ 1 ( c ) ξ } , ϕ ε ¯ ( ξ ) = min { W * ε , ( W * ε ) e λ 1 ( c ) ξ } .
From the proof of Theorem 2, ϕ * ¯ ( ξ ) and ϕ ε ¯ ( ξ ) are upper solutions to (34) and (35), respectively. Let
ϕ 0 ( ξ ) = ϕ * ¯ ( ξ ) , ϕ n + 1 ( ξ ) = F * ( ϕ n ) ( ξ ) , n = 0 , 1 , ,
ψ 0 ( ξ ) = ϕ ε ¯ ( ξ ) , ψ n + 1 ( ξ ) = F ε ( ψ n ) ( ξ ) , n = 0 , 1 , .
Since ϕ 0 ( ξ ) ψ 0 ( ξ ) , b * ( ξ ) b ε ( ξ ) , by the monotonicity of b * and b ε on [ 0 , W * ] , we have
ϕ 1 ( ξ ) = F * ( ϕ 0 ) ( ξ ) F * ( ψ 0 ) ( ξ ) F ε ( ψ 0 ) ( ξ ) = ψ 1 ( ξ ) ,
and inductively,
ϕ n ( ξ ) ψ n ( ξ ) , n = 1 , 2 , .
Moreover,
ϕ 0 ( ξ ) ϕ 1 ( ξ ) ϕ n ( ξ ) , ψ 0 ( ξ ) ψ 1 ( ξ ) ψ n ( ξ ) ,
ϕ 0 ( ξ ) φ * ( ξ ) = lim n ϕ n ( ξ ) lim n ψ n ( ξ ) = φ ε ( ξ ) 0 .
Since ϕ 0 ( ξ ) X λ for λ ( 0 , λ 1 ( c ) ) , we obtain φ ε X λ and φ * X λ . Hence, Γ * X λ .
To prove Γ * is non-empty, it suffices to show that φ * Γ * .
Note that
φ * ( ξ ) = F * ( φ * ) ( ξ ) = + K ( 1 , ξ ζ ) H * ( φ * ( ζ ) ) d ζ , ξ R ,
| H * ( φ ( ξ ) ) | ( 1 d ) | φ * ( ξ c ) | + η + K D i ( r 1 , y ) | b * ( φ * ( ξ y c r ) ) | d y ( 1 d ) W * + η b ( W * ) + K D i ( r 1 , y ) d y = W * .
For any ξ 1 , ξ 2 R , ξ 1 < ξ 2 ,
| φ * ( ξ 1 ) φ * ( ξ 2 ) | = | F * ( φ * ) ( ξ 1 ) F * ( φ * ) ( ξ 2 ) | + | K ( 1 , ξ 1 ζ ) K ( 1 , ξ 2 ζ ) | | H * ( φ * ( ζ ) ) | d ζ , W * + | K ( 1 , ξ 1 ζ ) K ( 1 , ξ 2 ζ ) | d ζ , W * 2 D ξ 1 | e 1 D ( ξ 1 ζ ) e 1 D ( ξ 2 ζ ) | d ζ + W * 2 D ξ 2 + | e 1 D ( ξ 1 ζ ) e 1 D ( ξ 2 ζ ) | d ζ + W * 2 D ξ 1 ξ 2 e 1 D ( ξ 2 ζ ) d ζ + W * 2 D ξ 1 ξ 2 e 1 D ( ξ 1 ζ ) d ζ 2 W * | 1 e 1 D ( ξ 2 ξ 1 ) | 2 W * D | ξ 2 ξ 1 | .
Therefore, φ * Γ * . Similarly, we can prove φ ε Γ * . Next, we proof that Γ * is convex in X λ . For any φ , ψ Γ , 0 λ 1 ,
λ φ ε ( ξ ) λ φ ( ξ ) λ φ * ( ξ ) , ( 1 λ ) ψ ε ( ξ ) ( 1 λ ) ψ ( ξ ) ( 1 λ ) ψ * ( ξ ) .
Thus, we have
λ φ ε ( ξ ) + ( 1 λ ) ψ ε ( ξ ) λ φ ( ξ ) + ( 1 λ ) ψ ( ξ ) λ φ * ( ξ ) + ( 1 λ ) ψ * ( ξ ) ,
| λ φ ( ξ 1 ) + ( 1 λ ) ψ ( ξ 1 ) [ λ φ ( ξ 2 ) + ( 1 λ ) ψ ( ξ 2 ) ] | 2 W * D | ξ 2 ξ 1 | .
Therefore, λ φ + ( 1 λ ) ψ Γ * . Proofs of the other conclusions of Lemma 8 are easy, we omit them. This complete the proof. □
Now, we are in position to state our main result.
Theorem 3.
Assume that (H1)–(H3) hold true. Then, there exists c * > 0 , such that for every c > c * , (14) admits a traveling wave solution φ ( x + c n ) satisfying
lim ξ φ ( ξ ) = 0 , W * lim inf ξ + φ ( ξ ) lim sup ξ + φ ( ξ ) W * .
Furthermore, if (H5) hold, then φ connects the equilibrium 0 and w * , i.e., lim ξ φ ( ξ ) = 0 , lim ξ + φ ( ξ ) = w * .
Proof. 
We will use Schauder fixed point theorem to prove that (14) has a traveling wave solution.
Recall that
F ( φ ) ( ξ ) = + K ( 1 , ξ ζ ) H ( φ ) ( ζ ) d ζ
where
H ( φ ) ( ξ ) = ( 1 d ) φ ( ξ c ) + η + K D i ( r 1 , y ) b ( φ ( ξ y c r ) ) d y , ξ R .
Obviously, for any φ Γ * ,
0 H ε ( φ ) ( ξ ) H ( φ ) ( ξ ) H * ( φ ) ( ξ ) W * , ξ R .
It follows that
0 F ( φ ) ( ξ ) W * ,
and, hence, F : Γ * C ( R , [ 0 , W * ] ) is well defined.
For any φ , ψ X λ , by the continuity of b ( w ) on [ 0 , W * ] , we have
| H ( φ ) ( ξ ) H ( ψ ) ( ξ ) | e λ ξ ( 1 d ) | φ ( ξ c ) ψ ( ξ c ) | e λ ξ + η + K D i ( r 1 , y ) | b ( φ ( ξ y c r ) ) b ( ψ ( ξ y c r ) ) | e λ ξ d y ( 1 d ) φ ψ λ + g ( 0 ) η e c r λ + K D i ( r 1 , y ) e λ y φ ψ λ d y ( 1 d ) + g ( 0 ) η e c r λ ( 1 D i λ 2 ) r 1 φ ψ λ = M φ ψ λ ,
where M = ( 1 d ) + g ( 0 ) η e c r λ ( 1 D i λ 2 ) r 1 . Therefore, we have
| F ( φ ) ( ξ ) F ( ψ ) ( ξ ) | e λ ξ + K ( 1 , ξ ζ ) | H ( φ ) ( ζ ) H ( ψ ) ( ζ ) | e λ ξ d ζ = + K ( 1 , ξ ζ ) e λ ( ξ ζ ) | H ( φ ) ( ζ ) H ( ψ ) ( ζ ) | e λ ζ d ζ M 1 λ 2 D φ ψ λ ,
which yields
F ( φ ) ( ξ ) F ( ψ ) ( ξ ) λ M 0 φ ψ λ .
where M 0 = M 1 λ 2 D . Therefore, F ( Γ * ) C ( R , [ 0 , W * ] ) is continuous.
In what follows, we verify that F ( Γ * ) Γ * . Note that φ ε ( ξ ) is a solution of (35), i.e.,
φ ϵ ( ξ ) = + K ( 1 , ξ ζ ) H ϵ ( φ ε ) ( ζ ) d ζ .
Then, for any φ Γ * ,
F ( φ ) ( ξ ) = + K ( 1 , ξ ζ ) H ( φ ) ( ζ ) d ζ + K ( 1 , ξ ζ ) H ε ( φ ) ( ζ ) d ζ + K ( 1 , ξ ζ ) H ε ( φ ε ) ( ζ ) d ζ = φ ε ( ξ )
Similarly, F ( φ ) ( ξ ) φ * ( ξ ) for all ξ R .
For any φ Γ * and ξ 1 , ξ 2 R with ξ 1 < ξ 2 , by the same argument as in the proof Lemma 8, we have
| F ( φ ) ( ξ 1 ) F ( φ ) ( ξ 2 ) | 2 W * D | ξ 1 ξ 2 | .
Therefore, we conclude that F ( φ ) Γ * for all φ Γ * .
Using Schauder’s fixed point theorem, we obtain that F has a fixed point φ in Γ * X λ , which satisfies
φ ( ξ ) = + K ( 1 , ξ ζ ) ( 1 d ) φ ( ζ c ) + η + K D i ( r 1 , y ) b ( φ ( ζ y c r ) ) d y d ζ .
and
φ ϵ ( ξ ) φ ( ξ ) φ * ( ξ ) for all ξ R .
Taking limits ξ and ξ + in (43), respectively, we have φ ( ) = 0 and
W * ε lim inf ξ + φ ( ξ ) lim sup ξ + φ ( ξ ) W * .
Then, letting ε 0 + in the above inequality, we obtain
W * lim inf ξ + φ ( ξ ) lim sup ξ + φ ( ξ ) W * .
This finishes the first part of the proof of Theorem 3.
Next, we denote
α = lim inf ξ + φ ( ξ ) , β = lim sup ξ + φ ( ξ ) .
Then W * α β W * . We will verify that α = β .
If this is invalid, then α < β . It is easy to see that if there exists a large number M > 0 , such that φ 0 or φ 0 on [ M , + ) , then lim ξ + φ ( ξ ) exists and α = β , which leads to a contradiction. So there must be a sequence { ξ j } j N with ξ j + as j + , such that φ ( ξ j ) = 0 , φ ( ξ j ) 0 and φ ( ξ j ) β as j + . It follows from (14) that
0 = D φ ( ξ j ) φ ( ξ j ) + ( 1 d ) φ ( ξ j c ) + η + K D i ( r 1 , y ) b ( φ ( ξ j y c r ) ) d y φ ( ξ j ) + ( 1 d ) φ ( ξ j c ) + η + K D i ( r 1 , y ) b ( φ ( ξ j y c r ) ) d y ,
thus
φ ( ξ j ) ( 1 d ) φ ( ξ j c ) η + K D i ( r 1 , y ) b ( φ ( ξ j y c r ) ) d y .
For any ε > 0 , there exists a sufficiently large constant N > 0 , such that
d W * / η | y | > N K D i ( r 1 , y ) d y < ε 2 .
Due to the continuity of b, we can choose δ > 0 with δ < ε satisfying
max { b ( w ) | w [ α δ , β + δ ] } < max { b ( w ) | α w β } + ε 2 η .
For such δ > 0 , take N 1 > 0 , such that
φ ( ξ ) [ α δ , β + δ ] , for all ξ N 1 .
Then choose J 0 > 0 , such that
ξ j N 1 + N + c r , for all j J 0 .
Therefore, for j J 0 , we have
φ ( ξ j ) ( 1 d ) φ ( ξ j c ) η | y | N K D i ( r 1 , y ) b ( φ ( ξ j y c r ) ) d y + η | y | > N K D i ( r 1 , y ) b ( φ ( ξ j y c r ) ) d y η max { b ( w ) | w [ α δ , β + δ ] } + d W * / η | y | > N K D i ( r 1 , y ) d y η max { b ( w ) | α w β } + ε .
Let j + , we find
d β η max { b ( w ) | α w β } + ε ,
It follows from ε 0 + that
d β η max { b ( w ) | α w β } .
Similarly, we have
d α η min { b ( w ) | α w β } .
If α < β w * , then by (47), we have d α η min { b ( w ) | α w β } > d α , a contradiction. If w * α < β , then (46) implies that d β η max { b ( w ) | α w β } < d β , also a contradiction. Therefore, we must have α < w * < β .
Let w 1 , w 2 [ α , β ] , such that b ( w 1 ) = max { b ( w ) | α w β } and b ( w 2 ) = min { b ( w ) | α w β } . We have the following three cases:
Case (i). w * w 1 β . If w 1 = β , then d β η b ( β ) follows from (46), which is impossible since β > w * . Therefore, we have w 1 < β and, hence,
d β η max { b ( w ) | α w β } = b ( w 1 ) η d w 1 < d β ,
which is a contradiction;
Case (ii). α w 2 w * . Using the similar argument that used in (i), we find α < w 2 and
d α η min { b ( w ) | α w β } = b ( w 2 ) η d w 2 > d α ,
which is also a contradiction;
Case (iii). w 1 < w * < w 2 . In this case, we have w 1 = α and w 2 = β . Otherwise, we have
d β d α η b ( w 1 ) η b ( w 2 ) < d ( 2 w * w 1 ) d ( 2 w * w 2 ) = d ( w 2 w 1 ) < d β d α
which is impossible.
Thus, α = β , and hence lim ξ + φ ( ξ ) = α [ W * , W * ] exists. Using the Lebesgue’s dominated convergence theorem and taking the limit as ξ + in (42), we have
α = lim ξ + + K ( 1 , ξ ζ ) ( 1 d ) φ ( ζ c ) + η + K D i ( r 1 , y ) b ( φ ( ζ y c r ) ) d y d ζ = ( 1 d ) α + η b ( α ) ,
which yields d α = η b ( α ) , hence lim ξ + φ ( ξ ) = α = w * . □

4. Discussion

In this paper, we studied the traveling wave solutions of the temporally discrete reaction-diffusion equation with monotone and non-monotone birth functions, respectively. Now we compare our results with the counterparts for continuous model. In [2], the authors studied the existence of traveling wavefronts of following equation,
w t = D 2 w x 2 d w + ε p + b ( w ( t r , y ) ) f α ( x y ) d y ,
where they took a particular birth function b ( w ) = p w e a w , p > 0 and a > 0 are parameters. They showed that if 1 < ε p d e , then there exists c ¯ * > 0 such that for every c > c ¯ * , (48) admits a traveling wavefront solution connecting the trivial equilibrium w 1 = 0 and positive equilibrium w 2 = 1 a l n ε p d . In [23], the author established a general existence result of traveling wave solutions for non-local reaction diffusion equation. As a special case, he proved that if e < ε p d e 2 , (48) still possess a traveling wave solution connecting w 1 and w 2 .
Using the non-standard discretization as the form of (1), the corresponding temporally discrete reaction diffusion equation reads as
w ( n + 1 , x ) = D Δ w ( n + 1 , x ) + ( 1 d ) w ( n , x ) + ε p + K D i ( r 1 , x y ) w ( n + 1 r , y ) e a w ( n + 1 r , y ) d y .
If 1 < ε p d e , then the assumptions (H1)–(H4) can be easily verified. By Theorem 2, there exists c * > 0 , such that, for every c > c * , (49) admits a traveling wavefront with speed c. When e < ε p d e 2 , by Theorem 3, we can prove that (49) still possess a traveling wave solution connecting w 1 and w 2 . This implies that in the sense of the existence of traveling wavefronts, (49) is a dynamically persistent discretization.
Now that both (48) and (49) have critical propagation speeds c ¯ * and c * , respectively. What is the relationship between the two propagation speeds c ¯ * and c * ? To answer this question, we can consider following two functions,
Δ ¯ ( c , λ ) = D λ 2 d c λ + ε p e α λ 2 λ c r ,
Δ ( c , λ ) = D λ 2 1 + ( 1 d ) e λ c + ε p e λ c r ( 1 D i λ 2 ) r 1 .
where α = ( r 1 ) D i . In fact, Δ ¯ ( c , λ ) = 0 and Δ ( c , λ ) = 0 are characteristic equations of linearizations at zero of wave profile equations corresponding to (48) and (49).
By direct computations, we find
d c λ < 1 + ( 1 d ) e λ c , e α λ 2 < 1 ( 1 D i λ 2 ) r 1 , λ ( 0 , 1 D i ) ,
which implies that for any given c > 0 , Δ ( c , λ ) > Δ ¯ ( c , λ ) for λ ( 0 , 1 D i ) .
By the definition of c ¯ * , there exists a unique λ = λ ( c ¯ * ) > 0 , such that Δ ¯ ( c ¯ * , λ ( c ¯ * ) ) = 0 . Moreover, Δ ¯ ( c ¯ * , λ ) > 0 for 0 < λ λ ( c ¯ * ) . Consequently, Δ ( c ¯ * , λ ) > Δ ¯ ( c ¯ * , λ ) 0 for any 0 λ 1 D i . It follows that c * > c ¯ * .

Author Contributions

Writing—original draft preparation, H.G.; writing—review and editing, Z.G. and H.G. Both authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by the National Natural Science Foundation of China (No. 11771104), Science and Technology Planning Project of Guangdong Province of China (2020A1414010106), The Natural Science Foundation of Guangdong Province of China (2021A1515010310) and Graduate Basic Innovation Project of Guangzhou University (No. 2019-GDJC-D03).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guo, Z.; Guo, H.; Chen, Y. Traveling wavefronts of a delayed temporally discrete reaction-diffusion equation. J. Math. Anal. Appl. 2021, 496, 124787. [Google Scholar] [CrossRef]
  2. So, J.; Wu, J.; Zou, X. A reaction-diffusion model for a single species with age structure. I Travelling wavefronts on unbounded domains. Proc. R. Soc. Lond. A. 2001, 457, 1841–1853. [Google Scholar] [CrossRef]
  3. Britton, N. Reaction-Diffusion Equations and Their Applications to Biology; Academic Press: New York, NY, USA, 1986. [Google Scholar]
  4. Bian, Q.; Zhang, W.; Yu, Z. Temporally discrete three-species Lotka-Volterra competitive systems with time delays. Taiwan J. Math. 2016, 20, 49–75. [Google Scholar] [CrossRef]
  5. Fife, P. Mathematical Aspects of Reaction and Diffusing Systems; Lecture Notes in Biomathematics; Springer: Berlin, Germany, 1979; Volume 28. [Google Scholar]
  6. Freedman, H. Deterministic Mathematical Models in Population Ecology; Marcel Dekker: New York, NY, USA, 1980. [Google Scholar]
  7. Goh, B. Management and Analysis of Biological Populations; Elsevier: Amsterdam, The Netherlands, 1980. [Google Scholar]
  8. Martin, R.H.; Smith, H.L. Abstract functional differential equations and reaction-diffusion systems. Trans. Am. Math. Soc. 1990, 321, 1–44. [Google Scholar]
  9. Murray, J.D. Mathematial Biology; Biomathematics; Springer: Berlin, Germany, 1993; Volume 19. [Google Scholar]
  10. Weinberger, H.F. Long-time behavior of a class of biological models. SIAM J. Math. Anal. 1982, 13, 353–396. [Google Scholar] [CrossRef]
  11. Mickens, R.E. Nonstandard Finite Difference Models of Differential Equations; World Scientific: Singapore, 1994. [Google Scholar]
  12. Mickens, R.E. Dynamic consistency: A fundamental principle for constructing nonstandard finite difference dchemes for differential equations. J. Differ. Equ. Appl. 2005, 11, 645–653. [Google Scholar] [CrossRef]
  13. Anguelov, R.; Dumont, Y.; Lubuma, J.M.-S.; Shillor, M. Dynamically consistent nonstandard finite difference schemes for epidemiological models. J. Comput. Appl. Math. 2014, 255, 161–182. [Google Scholar] [CrossRef]
  14. Dang, Q.A.; Hoang, M.T. Dynamically consistent discrete metapopulation model. J. Differ. Equ. Appl. 2016, 22, 1325–1349. [Google Scholar] [CrossRef]
  15. Dimitrov, D.T.; Kojouharov, H.V. Dynamically consistent numerical methods for general productive-destructive systems. J. Differ. Equ. Appl. 2011, 17, 1721–1736. [Google Scholar] [CrossRef]
  16. Roeger, L.-I.W. Dynamically consistent discrete-time Lotka-Volterra competition models. Discrete Contin. Dyn. Syst. 2009, 2009, 650–658. [Google Scholar]
  17. Suryanto, A.; Kusumawinahyu, W.M.; Darti, I.; Yanti, I. Dynamically consistent discrete epidemic model with modified saturated incidence rate. Comput. Appl. Math. 2013, 32, 373–383. [Google Scholar] [CrossRef]
  18. Weinberger, H.F.; Lewis, M.A.; Li, B. Analysis of linear determinacy for spread in cooperative models. J. Math. Biol. 2002, 45, 183–218. [Google Scholar] [CrossRef] [PubMed]
  19. Kot, M.; Schaffer, W.M. Discrete-time growth-dispersal models. J. Math. Biol. 1986, 80, 109–136. [Google Scholar] [CrossRef]
  20. Lin, G.; Li, W. Traveling wavefronts in temporally discrete reaction-diffusion equations with delay. Nonlinear Anal. Real World Appl. 2008, 9, 197–205. [Google Scholar] [CrossRef]
  21. Xue, H.; Huang, J.; Yu, Z. Invasion traveling wave solutions in temporally discrete random-diffusion systems with delays. Discret. Contin. Dyn. Syst. Ser. S 2017, 10, 1107–1131. [Google Scholar] [CrossRef] [Green Version]
  22. Xue, H.; Huang, J.; Yu, Z. Existence and asymptotic behavior of invasion wave solutions in temporally discrete diffusion systems with delays. Int. J. Biomath. 2018, 11, 1850016. [Google Scholar] [CrossRef]
  23. Ma, S. Traveling waves for non-local delayed diffusion equations via auxiliary equations. J. Differ. Equ. 2007, 237, 259–277. [Google Scholar] [CrossRef] [Green Version]
  24. Ma, S.; Weng, P.; Zou, X. Asymptotic speed of propagation and traveling wavefronts in a non-local delayed lattice differential equation. Nonlinear Anal. 2006, 65, 1858–1890. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, H.; Guo, Z. Traveling Waves Solutions for Delayed Temporally Discrete Non-Local Reaction-Diffusion Equation. Mathematics 2021, 9, 1999. https://doi.org/10.3390/math9161999

AMA Style

Guo H, Guo Z. Traveling Waves Solutions for Delayed Temporally Discrete Non-Local Reaction-Diffusion Equation. Mathematics. 2021; 9(16):1999. https://doi.org/10.3390/math9161999

Chicago/Turabian Style

Guo, Hongpeng, and Zhiming Guo. 2021. "Traveling Waves Solutions for Delayed Temporally Discrete Non-Local Reaction-Diffusion Equation" Mathematics 9, no. 16: 1999. https://doi.org/10.3390/math9161999

APA Style

Guo, H., & Guo, Z. (2021). Traveling Waves Solutions for Delayed Temporally Discrete Non-Local Reaction-Diffusion Equation. Mathematics, 9(16), 1999. https://doi.org/10.3390/math9161999

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop