Next Article in Journal
Factoring Continuous Characters Defined on Subgroups of Products of Topological Groups
Next Article in Special Issue
Bounded Perturbation Resilience of Two Modified Relaxed CQ Algorithms for the Multiple-Sets Split Feasibility Problem
Previous Article in Journal
A Novel Numerical Method for Solving Fractional Diffusion-Wave and Nonlinear Fredholm and Volterra Integral Equations with Zero Absolute Error
Previous Article in Special Issue
Delay-Dependent Stability, Integrability and Boundedeness Criteria for Delay Differential Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Stability of a Lotka-Volterra Competition-Diffusion-Advection System with Different Positive Diffusion Distributions

College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, China
*
Author to whom correspondence should be addressed.
Submission received: 22 June 2021 / Revised: 25 July 2021 / Accepted: 26 July 2021 / Published: 28 July 2021
(This article belongs to the Special Issue Special Issue in Honor of the 60th Birthday of Professor Hong-Kun Xu)

Abstract

:
In this paper, the problem of a Lotka–Volterra competition–diffusion–advection system between two competing biological organisms in a spatially heterogeneous environments is investigated. When two biological organisms are competing for different fundamental resources, and their advection and diffusion strategies follow different positive diffusion distributions, the functions of specific competition ability are variable. By virtue of the Lyapunov functional method, we discuss the global stability of a non-homogeneous steady-state. Furthermore, the global stability result is also obtained when one of the two organisms has no diffusion ability and is not affected by advection.

1. Introduction

For researchers from the fields of biology and mathematics, advancing the exploration of dynamic systems is a long-term challenge (see [1,2,3]). The competitive system of two diffusive organisms is often used to simulate population dynamics in biomathematics; for an example, see [1,2,4]. The key to spatial heterogeneity has been discussed in a lot of work, such as [2,5] and its references. In 2020, by proposing a new Lyapunov functional, Ni et al. [6] first studied and proved the global stability of a diffusive, competitive two-organism system, and then extended it to multiple organisms.
Since various methods in the reaction–diffusion–convection system cannot continue to work well, the global dynamics is far from being fully understood. In competitive diffusion advection systems, some progress has been made in [7,8,9,10,11]. Li et al. introduced the weighted Lyapunov functional related to the advection term to study global stability results in 2020 (see [12]), and studied the stability and bifurcation analysis of the model with the time delay term in 2021 (see [11]). Similarly, in 2021, Ma et al. described the overlapping characteristics of bifurcation solutions and studied the influence of advection on the stability of bifurcation solutions. Their results showed that the advection term may change its stability (see [13]). In 2021, Zhou et al. studied the global dynamics of a parabolic system using the competition coefficient (see [14]).
Motivated by the efforts of the aforementioned papers, we will investigate the global stability of a non-homogeneous steady-state solution of a Lotka–Volterra model between two organisms in heterogeneous environments, where two competing organisms have different intrinsic growth rates, advection and diffusion strategies, and follow different positive diffusion distributions.
Hence, we discuss the following advection system:
U t = · [ μ 1 ( x ) ( U ρ 1 ( x ) ) R 1 ( x ) U ρ 1 ( x ) B 1 ( x ) ] + U [ λ 1 ( x ) ϖ 11 ( x ) U ϖ 12 ( x ) V ] , i n Ω × R + , V t = · [ μ 2 ( x ) ( V ρ 2 ( x ) ) R 2 ( x ) V ρ 2 ( x ) B 2 ( x ) ] + V [ λ 2 ( x ) ϖ 21 ( x ) U ϖ 22 ( x ) V ] , i n Ω × R + , μ 1 ( x ) n ( U ρ 1 ) R 1 ( x ) U ρ 1 B 1 ( x ) n = 0 , o n Ω × R + , μ 2 ( x ) n ( V ρ 2 ) R 2 ( x ) V ρ 2 B 2 ( x ) n = 0 , o n Ω × R + , U ( x , 0 ) = U 0 ( x ) , 0 , V ( x , 0 ) = V 0 ( x ) , 0 , i n Ω ,
Here, U ( x , t ) and V ( x , t ) are the population densities of biological organisms, location x Ω , time t > 0 , which are supposed to be nonnegative. μ 1 ( x ) , μ 2 ( x ) > 0 correspond to the dispersal rates of two competing biological organisms, respectively. R 1 ( x ) , R 2 ( x ) > 0 correspond to the advection rates of two competing biological organisms, and B 1 ( x ) , B 2 ( x ) C 2 ( Ω ¯ ) are the nonconstant functions and represent the advective directions. Two bounded functions λ 1 ( x ) and λ 2 ( x ) are the intrinsic growth rates of competing organisms, ρ 1 ( x ) , ρ 2 ( x ) C 2 ( Ω ¯ ) are two positive diffusion distributions, respectively. ϖ i j ( x ) > 0 , i = 1 , 2 , j = 1 , 2 show the strength of competition ability. The spatial habitat Ω R N is a bounded smooth domain, 1 N Z ; n denotes the outward unit normal vector on the boundary Ω . No one can enter or leave the habitat boundary.
The following are our basic assumptions:
Hypothesis 1.
0 < μ i ( x ) , R i ( x ) C 1 + ϱ ( Ω ¯ ) , 0 < λ i ( x ) , ϖ i j ( x ) C ϱ ( Ω ¯ ) , ϱ ( 0 , 1 ) .
Hypothesis 2.
μ 1 ( x ) R 1 ( x ) = : c 1 > 0 , μ 2 ( x ) R 2 ( x ) = : c 2 > 0 , x Ω ¯ , w h e r e c 1 a n d c 2 a r e c o n s t a n t s .
To simplify the calculation, by letting u = e c 1 B 1 ( x ) U ρ 1 ( x ) , v = e c 2 B 2 ( x ) V ρ 2 ( x ) , the system (1) converts into the following coupled system
u t = e c 1 B 1 ( x ) ρ 1 ( x ) [ μ 1 ( x ) e c 1 B 1 ( x ) u ] + u [ λ 1 ( x ) ϖ 11 ( x ) u e c 1 B 1 ( x ) ρ 1 ( x ) ϖ 12 ( x ) v e c 2 B 2 ( x ) ρ 2 ( x ) ] , i n Ω × R + , v t = e c 2 B 2 ( x ) ρ 2 ( x ) [ μ 2 ( x ) e c 2 B 2 ( x ) v ] + v [ λ 2 ( x ) ϖ 21 ( x ) u e c 1 B 1 ( x ) ρ 1 ( x ) ϖ 22 ( x ) v e c 2 B 2 ( x ) ρ 2 ( x ) ] , i n Ω × R + , u n = v n = 0 , o n Ω × R + , u ( x , 0 ) = e c 1 B 1 ( x ) U 0 ( x ) ρ 1 ( x ) , 0 , v ( x , 0 ) = e c 2 B 2 ( x ) V 0 ( x ) ρ 2 ( x ) , 0 , i n Ω ,
when c 1 = c 2 = 0 , ρ 1 ( x ) = ρ 2 ( x ) = 1 , the model (2) has been studied in Ni et al. [6]. c 1 = c 2 , B 1 ( x ) = B 2 ( x ) , ρ 1 ( x ) = ρ 2 ( x ) = 1 , the model (2) has been studied in Li et al. [12].
The rest of this article is arranged as follows. In Section 2, we carry out some preparatory work and give four lemmas, where some related properties of the system (1) are deduced from the properties of a single organism model (4). Using the Lyapunov functional method, we will provide and prove our main results in Section 3. In Section 4, one example is given to explain our conclusions.

2. Preliminaries

In order to describe our main results, we present the following uniform estimates for the parabolic equation:
w t = ϖ i j ( x ) D i j w + β j ( x ) D j w + λ ( x ) w + H ( x , t , w ) , i n Ω × R + , w n = 0 , o n Ω × R + , w ( x , 0 ) = w 0 ( x ) , 0 , i n Ω ,
where Ω R N is bounded and Ω C 2 + ϱ ( ϱ ( 0 , 1 ) ) is a smooth boundary. The initial condition w 0 ( x ) W 2 , p ( Ω ) , p > 1 + N 2 .
Setting the following assumptions:
( A 1 ) Let ϖ i j , β j , λ C ( Ω ¯ ) , χ 1 , χ 2 > 0 , such that
χ 1 | y | 2 1 i , j N ϖ i j ( x ) y i y j χ 2 | y | 2 , | β j ( x ) | , | λ ( x ) | χ 2 , f o r a l l x Ω , y R N .
( A 2 ) Let Λ > 0 be a constant, such that
ϖ i j C ϱ ( Ω ¯ ) , β j C ϱ ( Ω ¯ ) , λ C ϱ ( Ω ¯ ) Λ .
( A 3 ) H L ( Ω × [ 0 , ) × [ τ 1 , τ 2 ] ) for some τ 1 < τ 2 and there is Λ ( τ 1 , τ 2 ) > 0 such that
| H ( x , t , w 1 ) H ( x , t , w 2 ) | Λ ( τ 1 , τ 2 ) | w 1 w 2 | , f o r a l l ( x , t ) Ω × [ 0 , ) , w 1 , w 2 [ τ 1 , τ 2 ] ,
and there exists Λ > 0 , satisfying
| H ( x 1 , t 1 , w ) H ( x 2 , t 2 , w ) | Λ ( | x 1 x 2 | ϱ + | t 1 t 2 | ϱ 2 ) f o r a l l ( x 1 , t 1 ) , ( x 2 , t 2 ) Ω ¯ × [ d , d + 3 ] , u [ τ 1 , τ 2 ] , d 0 .
The following lemma (see [15,16]) is the boundedness result of the solution w ( x , t ) in (3).
Lemma 1.
Let w ( x , t ) be a solution of (3) with τ 1 < w < τ 2 , τ 1 , τ 2 R . Suppose that f, ϖ i j , β j , λ satisfy the assumptions ( A 1 ) ( A 3 ) , then for any κ 1 , there is a constant Λ ( κ ) > 0 such that
max x Ω ¯ w t ( x , · ) C ϱ 2 ( [ κ , + ) ) + max t κ w t ( · , t ) C ( Ω ¯ ) + max t κ w ( · , t ) C 2 + ϱ ( Ω ¯ ) Λ ( κ ) .
In the proof of global stability, the following calculus theory and integral inequality are very important. For details, see [6,17].
Lemma 2
([17]). Let β , λ > 0 be constants, φ ( t ) 0 in [ β , ) . Assume that ϕ C 1 ( [ β , ) ) has lower bound, ϕ ( t ) λ φ ( t ) in [ β , ) . If one of the following alternatives holds:
  • φ C 1 ( [ β , ) ) and φ ( t ) P in [ β , ) for P > 0 ,
  • φ C ϱ ( [ β , ) ) and φ C ϱ ( [ β , ) ) P for 0 < m < 1 and P > 0 ,
where P and m are constants, then lim t φ ( t ) = 0 .
Lemma 3
([6]). Let α , α C 2 ( Ω ¯ ) with α , α > 0 and m C 1 ( Ω ¯ ) , b C 2 ( Ω ¯ ) with m , b 0 , α , α , m , b are functions. If the following conditions holds:
  • q 1 is a constant, the function h C 0 , 1 ( Ω × [ 0 , ) ) , x Ω , h ( x , K ) K is a non-increasing function for K [ 0 , ) ,
  • ( b ( x ) α ) ν = h ( x , α ) , ( b ( x ) α ) ν = h ( x , α ) on Ω ,
then
Ω b ( x ) α [ α q α q ] α q { m ( x ) [ b ( x ) α ] } α α { m ( x ) [ b ( x ) α ] } d x Ω q m b 2 α 2 ( α α ) q 1 | α α | 2 d x 0 .
Next, we consider the following scalar evolution eqution
u t = e c B ( x ) ρ ( x ) [ μ ( x ) e c B ( x ) u ] + u [ λ ( x ) ϖ ( x ) u e c B ( x ) ρ ( x ) ] , i n Ω × R + , u n = 0 , o n Ω × R + , u ( x , 0 ) = e c B ( x ) U 0 ( x ) ρ ( x ) , 0 , i n Ω ,
where μ ( x ) , c , ϖ ( x ) , λ ( x ) satisty
0 < μ ( x ) , R ( x ) C 1 + ϱ ( Ω ¯ ) , 0 < λ ( x ) , ϖ ( x ) C ϱ ( Ω ¯ ) , ϱ ( 0 , 1 ) , μ ( x ) R ( x ) = c , w h e r e c i s a c o n s t a n t .
Now we see the following useful lemma.
Lemma 4
([1]). Assume that 0 < μ ( x ) , λ ( x ) , ρ ( x ) , ϖ ( x ) on Ω ¯ , then the elliptic problem:
e c B ( x ) ρ ( x ) [ μ ( x ) e c B ( x ) u ] + u [ λ ( x ) ϖ ( x ) u e c B ( x ) ρ ( x ) ] = 0 , i n Ω , u n = 0 , o n Ω ,
has a unique positive solution, denoted by u θ .

3. Main Results

In this section, firstly, by utilizing the Lyapunov function method, the global stability of the model (5) is obtained, and we can see that the non-constant steady-state for (5) is equivalent to the solution u θ of (7).
Theorem 1.
Assume that u 0 ( x ) 0 . If μ , ρ , c , λ , ϖ satisfy (6), then Equation (5) has a unique solution u ( x , t ) > 0 with lim t u ( x , t ) = u θ in C 2 ( Ω ) .
Proof. 
According to the upper–lower solutions method [1,18], we obtain (5) with a unique solution u ( x , t ) > 0 . Let M be a upper solution of (5), we have 0 < u ( x , t ) < M , ( x , t ) Ω ¯ × ( 0 , ) .
By applying Lemma 1, we can obtain that there exists a constant Λ > 0 such that
max t 1 u t ( · , t ) C ( Ω ¯ ) + max t 1 u ( · , t ) C 2 + ϱ ( Ω ¯ ) Λ .
Then, define a function Φ : [ 0 , ) R by
Φ ( t ) = Ω ρ u θ e c B ( u u θ u θ ln u u θ ) d x .
Then, Φ ( t ) 0 , t 0 . By (2) and (4), we have
Φ ( t ) = Ω ρ u θ e c B ( 1 u θ u ) u t d x = Ω ρ u θ e c B ( 1 u θ u ) [ e c B ρ ( μ e c B u ) + u ( λ ϖ u e c B ρ ) ] d x = Ω ρ u θ e c B ( 1 u θ u ) [ e c B ρ ( μ e c B u ) u e c B u θ ρ ( μ e c B u θ ) ] d x + Ω ρ u θ e c B ( 1 u θ u ) [ u ( λ ϖ u e c B ρ ) u u θ u θ ( λ ϖ u θ e c B ρ ) ] d x Ω μ e c B u 2 | u θ u | 2 d x Ω ρ 2 u θ e 2 c B ϖ ( u u θ ) 2 d x .
We get
Φ ( t ) Ω ρ 2 u θ e 2 c B ϖ ( u u θ ) 2 d x = : φ ( t ) 0 .
By virtue of (8), we get | φ ( t ) | Λ in [ 1 , ) for some Λ > 0 . From Lemma 2, it follows that
lim t φ ( t ) = lim t Ω ρ 2 u θ e 2 c B ϖ ( u u θ ) 2 d x = 0 .
Applying (8) again, { u ( · , t ) : t 1 } is relatively compact in C 2 ( Ω ¯ ) . It can be found that there exists some function u ( x ) C 2 ( Ω ¯ ) such that
u ( · , t s ) u C 2 ( Ω ¯ ) 0 a s t s .
Combining with (12), we get u ( x ) = u θ ( x ) where x Ω . Hence, we deduce
lim t u ( x , t ) = u θ ( x ) i n C 2 ( Ω ¯ ) .
In addition, taking advantage of Lyapunov function method, the global stability results of (2) are obtained.
Theorem 2.
Suppose that u 0 ( x ) , v 0 ( x ) , 0 , ( H 1 ) and ( H 2 ) hold, the system (2) admits a non-homogeneous steady-state ( u θ ˜ ( x ) , v θ ˜ ( x ) ) > 0 and there exists
η 1 > 0 , η 2 > 0 s u c h t h a t η 1 u θ ˜ ( x ) v θ ˜ ( x ) η 2 , x Ω ¯ .
Suppose that
η 2 η 1 < min Ω ¯ ϖ 11 ϖ 22 ϖ 12 ϖ 21 .
Then, the system (2) admits a solution ( u ( x , t ) , v ( x , t ) ) that satisfies
lim t u ( x , t ) = u θ ˜ ( x ) , lim t v ( x , t ) = v θ ˜ ( x ) i n C 2 ( Ω ¯ ) .
Proof. 
Assume that the inequality (15) holds, let Φ : [ 0 , + ) R defined by
Φ ( t ) = Ω ρ 1 u θ ˜ e c 1 B 1 ( u u θ ˜ u θ ˜ ln u u θ ˜ ) d x + Ω ξ ρ 2 v θ ˜ e c 2 B 2 ( v v θ ˜ v θ ˜ ln v v θ ˜ ) d x ,
where 0 < ξ ( x ) : = ϖ 12 η 1 η 2 ϖ 21 . Clearly, Φ ( t ) 0 . By (2) and (4), we have
Φ ( t ) = Ω [ ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) u t + ξ ρ 2 v θ ˜ e c 2 B 2 ( 1 v θ ˜ v ) v t ] d x = Ω ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) [ e c 1 B 1 ρ 1 ( μ 1 e c 1 B 1 u ) + u ( λ 1 ϖ 11 u e c 1 B 1 ρ 1 ϖ 12 v e c 2 B 2 ρ 2 ) ] d x + Ω ξ ρ 2 v θ ˜ e c 2 B 2 ( 1 v θ ˜ v ) [ e c 2 B 2 ρ 2 ( μ 2 e c 2 B 2 v ) + v ( λ 2 ϖ 21 u e c 1 B 1 ρ 1 ϖ 22 v e c 2 B 2 ρ 2 ) ] d x = Ω ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) [ e c 1 B 1 ρ 1 ( μ 1 e c 1 B 1 u ) u e c 1 B 1 u θ ˜ ρ 1 ( μ 1 e c 1 B 1 u θ ˜ ) ] d x + Ω ξ ρ 2 v θ ˜ e c 2 B 2 ( 1 v θ ˜ v ) [ e c 2 B 2 ρ 2 ( μ 2 e c 2 B 2 v ) v e c 2 B 2 v θ ˜ ρ 2 ( μ 2 e c 2 B 2 v θ ˜ ) ] d x + Ω ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) u ( λ 1 ϖ 11 u e c 1 B 1 ρ 1 ϖ 12 v e c 2 B 2 ρ 2 ) d x Ω ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) u u θ ˜ u θ ˜ ( λ 1 ϖ 11 u θ ˜ e c 1 B 1 ρ 1 ϖ 12 v θ ˜ e c 2 B 2 ρ 2 ) ] d x + Ω ξ ρ 2 v θ ˜ e c 2 B 2 ( 1 v θ ˜ v ) v ( λ 2 ϖ 21 u e c 1 B 1 ρ 1 ϖ 22 v e c 2 B 2 ρ 2 ) d x Ω ξ ρ 2 v θ ˜ e c 2 B 2 ( 1 v θ ˜ v ) v v θ ˜ v θ ˜ ( λ 2 ϖ 21 u θ ˜ e c 1 B 1 ρ 1 ϖ 22 v θ ˜ e c 2 B 2 ρ 2 ) d x Ω μ 1 e c 1 B 1 u 2 | u θ ˜ u | 2 d x Ω μ 2 e c 2 B 2 v 2 | v θ ˜ v | 2 d x Ω ρ 1 2 u θ ˜ e 2 c 1 B 1 ϖ 11 ( u u θ ˜ ) 2 d x Ω ρ 1 ρ 2 e c 1 B 1 + c 2 B 2 ( ϖ 12 u θ ˜ + ξ ϖ 21 v θ ˜ ) ( u u θ ˜ ) ( v v θ ˜ ) d x Ω ξ ρ 2 2 v θ ˜ e 2 c 2 B 2 ϖ 22 ( v v θ ˜ ) 2 d x .
Note that (14) and (15) give rise to
2 ρ 1 2 u θ ˜ e 2 c 1 B 1 ϖ 11 ξ ρ 2 2 v θ ˜ e 2 c 2 B 2 ϖ 22 ρ 1 ρ 2 e c 1 B 1 + c 2 B 2 ( ϖ 12 u θ ˜ + ξ ϖ 21 v θ ˜ ) = 2 ρ 1 ρ 2 e c 1 B 1 + c 2 B 2 ξ u θ ˜ v θ ˜ ϖ 11 ϖ 22 ρ 1 ρ 2 e c 1 B 1 + c 2 B 2 ( ϖ 12 u θ ˜ + ξ ϖ 21 v θ ˜ ) = ρ 1 ρ 2 e c 1 B 1 + c 2 B 2 ( 2 ξ u θ ˜ v θ ˜ ϖ 11 ϖ 22 ξ u θ ˜ v θ ˜ ( ϖ 12 u θ ˜ ξ v θ ˜ + ϖ 21 ξ v θ ˜ u θ ˜ ) ) = ρ 1 ρ 2 e c 1 B 1 + c 2 B 2 ξ u θ ˜ v θ ˜ ( 2 ϖ 11 ϖ 22 ( ϖ 12 η 2 ξ + ϖ 21 ξ η 1 ) ) ρ 1 ρ 2 e c 1 B 1 + c 2 B 2 ξ u θ ˜ v θ ˜ ( 2 ϖ 11 ϖ 22 2 ϖ 12 ϖ 21 η 2 η 1 ) > 0 .
Choosing 0 < ε 1 , we have
2 ρ 1 2 u θ ˜ e 2 c 1 B 1 ( ϖ 11 ε ) ξ ρ 2 2 v θ ˜ e 2 c 2 B 2 ( ϖ 22 ε ) ρ 1 ρ 2 e c 1 B 1 + c 2 B 2 ( ϖ 12 u θ ˜ + ξ ϖ 21 v θ ˜ ) > 0 .
Combining with (17), we can deduce
Φ ( t ) Ω [ ρ 1 2 u θ ˜ e 2 c 1 B 1 ε ( u u θ ˜ ) 2 + ξ ρ 2 2 v θ ˜ e 2 c 2 B 2 ε ( v v θ ˜ ) 2 ] d x = : φ ( t ) 0 .
From (13), it follows that
lim t u ( x , t ) = u θ ˜ ( x ) , lim t v ( x , t ) = v θ ˜ ( x ) i n C 2 ( Ω ¯ ) .
Finally, we consider that if one of the two organisms has no diffusion ability and is not affected by advection, the Lyapunov function method can also deduce the following global stability results in (2).
Theorem 3.
If u 0 , v 0 C ( Ω ¯ ) satisfy u 0 ( x ) , 0 and v 0 ( x ) > 0 on Ω ¯ . Let μ 1 ( x ) R 1 ( x ) = : c 1 , μ 2 ( x ) = R 2 ( x ) = 0 for x Ω ¯ , and
ϖ 12 ( x ) ϖ 21 ( x ) < ϖ 11 ( x ) ϖ 22 ( x ) , x Ω ¯ .
( i ) If
ϖ 22 ( x ) λ 1 ( x ) ϖ 12 ( x ) λ 2 ( x ) > 0 , x Ω ¯ ,
and
min Ω ¯ λ 2 ( x ) ρ 1 ( x ) ϖ 21 ( x ) e c 1 B 1 ( x ) > max Ω ¯ ϖ 22 ( x ) λ 1 ( x ) ϖ 12 ( x ) λ 2 ( x ) ρ 1 ( x ) e c 1 B 1 ( x ) ( ϖ 11 ( x ) ϖ 22 ( x ) ϖ 12 ( x ) ϖ 21 ( x ) ) ,
then there is a unique non-homogeneous steady-state ( u θ ˜ ( x ) , v θ ˜ ( x ) ) > 0 for the model (2) such that
lim t ( u ( x , t ) , v ( x , t ) ) = ( u θ ˜ ( x ) , v θ ˜ ( x ) ) i n C 1 ( Ω ¯ ) × L 2 ( Ω ) .
( i i ) If
λ 2 ( x ) ρ 1 ( x ) ϖ 21 ( x ) e c 1 B 1 ( x ) u θ ˜ ( x ) , x Ω ¯ ,
then there exists a semi-trivial steady-state ( u θ ˜ ( x ) , 0 ) for the model (2) such that
lim t ( u ( x , t ) , v ( x , t ) ) = ( u θ ˜ ( x ) , 0 ) i n C 1 ( Ω ¯ ) × L 2 ( Ω ¯ ) .
( i i i ) Let
ϖ 22 ( x ) ϖ 12 ( x ) λ 2 ( x ) λ 1 ( x ) , x Ω ¯ ,
then the model (2) has a semi-trivial steady-state ( 0 , v θ ˜ ( x ) ) ,
lim t ( u ( x , t ) , v ( x , t ) ) = ( 0 , v θ ˜ ( x ) ) i n C 1 ( Ω ¯ ) × L 2 ( Ω ¯ ) ,
where v θ ˜ ( x ) = λ 2 ( x ) ρ 2 ( x ) ϖ 22 ( x ) e c 2 B 2 ( x ) .
Proof. 
( i ) When μ 2 ( x ) = R 2 ( x ) = 0 , x Ω ¯ , ( u θ ˜ ( x ) , v θ ˜ ( x ) ) of the model (2) satisfies
e c 1 B 1 ( x ) ρ 1 ( x ) [ μ 1 ( x ) e c 1 B 1 ( x ) u ] + u [ λ 1 ( x ) ϖ 12 ( x ) ϖ 22 ( x ) λ 2 ( x ) ρ 1 ( x ) u e c 1 B 1 ( x ) ( ϖ 11 ( x ) ϖ 12 ( x ) ϖ 21 ( x ) ϖ 22 ( x ) ) ] = 0 , x Ω , u n = 0 , x Ω ,
and v θ ˜ = λ 2 ϖ 21 ρ 1 u θ ˜ e c 1 B 1 ϖ 22 ρ 2 e c 2 B 2 .
If (18) and (19) hold, we see μ 1 , λ 1 ϖ 12 ϖ 22 λ 2 , ρ 1 e c 1 B 1 ( ϖ 11 ϖ 12 ϖ 21 ϖ 22 ) > 0 , then by Lemma 4, the problem (23) has a unique solution u θ ˜ ( x ) > 0 . By using the maximum principle in elliptic equation, we infer
u θ ˜ < max Ω ¯ ϖ 22 λ 1 ϖ 12 λ 2 ρ 1 e c 1 B 1 ( ϖ 11 ϖ 22 ϖ 12 ϖ 21 ) .
According to (20), we can get v θ ˜ = λ 2 ϖ 21 ρ 1 u θ ˜ e c 1 B 1 ϖ 22 ρ 2 e c 2 B 2 > 0 , hence there exists a unique steady-state for (2), ( u θ ˜ ( x ) , v θ ˜ ( x ) ) > 0 .
Let us define a function Φ : [ 0 , ) R ,
Φ ( t ) = Ω ρ 1 u θ ˜ e c 1 B 1 ( u u θ ˜ u θ ˜ ln u u θ ˜ ) d x + Ω ξ ρ 2 e c 2 B 2 ( v v θ ˜ v θ ˜ ln v v θ ˜ ) d x ,
where ξ ( x ) = ϖ 12 ( x ) u θ ˜ ( x ) ϖ 21 ( x ) > 0 . Clearly, Φ ( t ) 0 . From (2) and (4), we get
Φ ( t ) = Ω [ ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) u t + ξ ρ 2 e c 2 B 2 ( 1 v θ ˜ v ) v t ] d x = Ω ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) [ e c 1 B 1 ρ 1 ( μ 1 e c 1 B 1 u ) + u ( λ 1 ϖ 11 u e c 1 B 1 ρ 1 ϖ 12 v e c 2 B 2 ρ 2 ) ] d x + Ω ξ ρ 2 v θ ˜ e c 2 B 2 ( 1 v θ ˜ v ) [ v ( λ 2 ϖ 21 u e c 1 B 1 ρ 1 ϖ 22 v e c 2 B 2 ρ 2 ) ] d x = Ω ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) [ e c 1 B 1 ρ 1 ( μ 1 e c 1 B 1 u ) u e c 1 B 1 u θ ˜ ρ 1 ( μ 1 e c 1 B 1 u θ ˜ ) ] d x + Ω ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) u ( λ 1 ϖ 11 u e c 1 B 1 ρ 1 ϖ 12 v e c 2 B 2 ρ 2 ) d x Ω ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) u u θ ˜ u θ ˜ ( λ 1 ϖ 11 u θ ˜ e c 1 B 1 ρ 1 ϖ 12 v θ ˜ e c 2 B 2 ρ 2 ) ] d x + Ω ξ ρ 2 e c 2 B 2 ( 1 v θ ˜ v ) v ( λ 2 ϖ 21 u e c 1 B 1 ρ 1 ϖ 22 v e c 2 B 2 ρ 2 ) d x Ω ξ ρ 2 e c 2 B 2 ( 1 v θ ˜ v ) v v θ ˜ v θ ˜ ( λ 2 ϖ 21 u θ ˜ e c 1 B 1 ρ 1 ϖ 22 v θ ˜ e c 2 B 2 ρ 2 ) d x Ω μ 1 e c 1 B 1 u 2 | u θ ˜ u | 2 d x Ω ρ 1 2 u θ ˜ e 2 c 1 B 1 ϖ 11 ( u u θ ˜ ) 2 d x Ω ρ 1 ρ 2 e c 1 B 1 + c 2 B 2 ( ϖ 12 u θ ˜ + ξ ϖ 21 ) ( u u θ ˜ ) ( v v θ ˜ ) d x Ω ξ ρ 2 2 e 2 c 2 B 2 ϖ 22 ( v v θ ˜ ) 2 d x .
We can choose 0 < ε 1 and use (18), such that
2 ρ 1 2 u θ ˜ e 2 c 1 B 1 ( ϖ 11 ε ) ξ ρ 2 2 e 2 c 2 B 2 ( ϖ 22 ε ) ρ 1 ρ 2 e c 1 B 1 + c 2 B 2 ( ϖ 12 u θ ˜ + ξ ϖ 21 ) > 0 .
Combining this with (24), we can deduce
Φ ( t ) Ω [ ρ 1 2 u θ ˜ e 2 c 1 B 1 ε ( u u θ ˜ ) 2 + ξ ρ 2 2 e 2 c 2 B 2 ε ( v v θ ˜ ) 2 ] d x = : φ ( t ) 0 .
Applying the Lemma 1 and Sobolev embedding theorem, we deduce that u and v are bounded in Ω × [ 0 , ) and there is a constant Λ > 0 such that
max t 1 u ( · , t ) C 1 + ϱ ( Ω ¯ ) Λ f o r s o m e 0 < ϱ < 1 .
Combining with (2) and | φ ( t ) | < Λ 1 in [ 1 , ) for some Λ 1 > 0 , and making use of Lemma 2, we get lim t φ ( t ) = 0 and we deduce that
lim t u ( x , t ) = u θ ˜ ( x ) , lim t v ( x , t ) = v θ ˜ ( x ) i n L 2 ( Ω ) .
Applying Theorem 2, we get lim t u ( x , t ) = u θ ˜ ( x ) in C 1 ( Ω ¯ ) .
( i i ) Let’s define a function Φ : [ 0 , ) R ,
Φ ( t ) = Ω ρ 1 u θ ˜ e c 1 B 1 ( u u θ ˜ u θ ˜ ln u u θ ˜ ) d x + Ω ξ ρ 2 e c 2 B 2 v d x ,
where ξ ( x ) = ϖ 12 ( x ) u θ ˜ ( x ) ϖ 21 ( x ) > 0 . From (4) and (21), we have
Φ ( t ) = Ω ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) [ e c 1 B 1 ρ 1 ( μ 1 e c 1 B 1 u ) u e c 1 B 1 u θ ˜ ρ 1 ( μ 1 e c 1 B 1 u θ ˜ ) ] d x + Ω ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) u ( λ 1 ϖ 11 u e c 1 B 1 ρ 1 ϖ 12 v e c 2 B 2 ρ 2 ) d x Ω ρ 1 u θ ˜ e c 1 B 1 ( 1 u θ ˜ u ) u u θ ˜ u θ ˜ ( λ 1 ϖ 11 u θ ˜ e c 1 B 1 ρ 1 ) ] d x + Ω ξ ρ 2 e c 2 B 2 v ( λ 2 ϖ 21 u e c 1 B 1 ρ 1 ϖ 22 v e c 2 B 2 ρ 2 ) d x Ω μ 1 e c 1 B 1 u 2 | u θ ˜ u | 2 d x Ω ρ 1 u θ ˜ e c 1 B 1 ( u u θ ˜ ) ( ϖ 11 e c 1 B 1 ρ 1 ( u u θ ˜ ) ϖ 12 v e c 2 B 2 ρ 2 ) d x + Ω ξ ρ 2 e c 2 B 2 v [ ( λ 2 ϖ 21 u θ ˜ e c 1 B 1 ρ 1 ) ϖ 21 e c 1 B 1 ρ 1 ( u u θ ˜ ) ϖ 22 v e c 2 B 2 ρ 2 ] d x Ω ρ 1 2 u θ ˜ e 2 c 2 B 2 ϖ 11 ( u u θ ˜ ) 2 d x Ω ξ ρ 2 2 e 2 c 2 B 2 ϖ 22 v 2 d x Ω ρ 1 ρ 2 e c 1 B 1 + c 2 B 2 ( ϖ 12 u θ ˜ + ξ ϖ 21 ) ( u u θ ˜ ) v d x .
The following discussion will refer to the part ( i ) , then we will not repeat it.
( i i i ) Clearly, (2) has a semi-trivial steady-state ( 0 , λ 2 ( x ) ρ 2 ( x ) ϖ 22 ( x ) e c 2 B 2 ( x ) ) . Let us define a function Φ : [ 0 , ) R ,
Φ ( t ) = Ω ρ 1 e c 1 B 1 u d x + Ω ξ ρ 2 e c 2 B 2 ( v v θ ˜ v θ ˜ ln v v θ ˜ ) d x ,
where ξ ( x ) = ϖ 12 ( x ) ϖ 21 ( x ) > 0 and v θ ˜ ( x ) = λ 2 ( x ) ρ 2 ( x ) ϖ 22 ( x ) e c 2 B 2 ( x ) . From (22), we have
Φ ( t ) = Ω ρ 1 e c 1 B 1 u ( λ 1 ϖ 11 u e c 1 B 1 ρ 1 ϖ 12 v e c 2 B 2 ρ 2 ) d x + Ω ξ ρ 2 e c 2 B 2 v ( λ 2 ϖ 21 u e c 1 B 1 ρ 1 ϖ 22 v e c 2 B 2 ρ 2 ) d x Ω ξ ρ 2 e c 2 B 2 v v θ ˜ v θ ˜ ( λ 2 ϖ 22 v θ ˜ e c 2 B 2 ρ 2 ) ] d x = Ω ρ 1 e c 1 B 1 u [ ( λ 1 ϖ 12 v θ ˜ e c 2 B 2 ρ 2 ) ϖ 11 u e c 1 B 1 ρ 1 ϖ 12 e c 2 B 2 ρ 2 ( v v θ ˜ ) ] d x Ω ξ ρ 2 e 2 c 2 B 2 ( v v θ ˜ ) [ ϖ 21 u e c 1 B 1 ρ 1 ϖ 22 e c 2 B 2 ρ 2 ( v v θ ˜ ) ] d x Ω ρ 1 2 e 2 c 2 B 2 ϖ 11 u 2 d x Ω ξ ρ 2 2 e 2 c 2 B 2 ϖ 22 ( v v θ ˜ ) 2 d x Ω ρ 1 ρ 2 e c 1 B 1 + c 2 B 2 ( ϖ 12 + ξ ϖ 21 ) ( v v θ ˜ ) u d x .
The following discussion is similar to the part ( i ) , so we omit it. □

4. Example

See the following parabolic problem:
u t = e c B ( x ) ρ ( x ) [ μ 1 ( x ) e c B ( x ) u ] + u [ λ ¯ 1 φ ( x ) e c B ( x ) ρ ( x ) + ε 1 g 1 ( x ) ϖ ¯ 11 φ ( x ) u e c B ( x ) ρ ( x ) ϖ ¯ 12 φ ( x ) v e c B ( x ) ρ ( x ) ] , i n Ω × R + , v t = e c B ( x ) ρ ( x ) [ μ 2 ( x ) e c B ( x ) v ] + v [ λ ¯ 2 φ ( x ) e c B ( x ) ρ ( x ) + ε 2 g 2 ( x ) ϖ ¯ 21 φ ( x ) u e c B ( x ) ρ ( x ) ϖ ¯ 22 φ ( x ) v e c B ( x ) ρ ( x ) ] , i n Ω × R + , u n = v n = 0 , o n Ω × R + , u ( x , 0 ) = e c B ( x ) U 0 ( x ) ρ ( x ) , 0 , v ( x , 0 ) = e c B ( x ) V 0 ( x ) ρ ( x ) , 0 , i n Ω ,
where λ ¯ i , ϖ ¯ i j , ε i are all positive constants, B , ρ C 2 ( Ω ¯ ) , μ i C 1 + ϱ ( Ω ¯ ) , φ , g i C ϱ ( Ω ¯ ) and φ ( x ) , μ i ( x ) > 0 on Ω ¯ .
Proposition 1.
If 0 ε i 1 and ϖ ¯ 21 ϖ ¯ 11 < λ ¯ 2 λ ¯ 1 < ϖ ¯ 22 ϖ ¯ 12 , ϖ ¯ 11 ϖ ¯ 22 ϖ ¯ 12 ϖ ¯ 21 > 1 , then there exists η 1 > 0 , η 2 > 0 such that
ϖ ¯ 11 ϖ ¯ 22 ϖ ¯ 12 ϖ ¯ 21 > η 2 η 1
and the system (25) admits a positive non-homogeneous steady-state ( u θ ˜ ( x ) , v θ ˜ ( x ) ) , which satisfies η 1 u θ ˜ ( x ) v θ ˜ ( x ) η 2 .
Proof. 
The steady-state of (25) satisfies the following elliptic problem
e c B ( x ) ρ ( x ) [ μ 1 ( x ) e c B ( x ) u ] + u [ λ ¯ 1 φ ( x ) e c B ( x ) ρ ( x ) + ε 1 g 1 ( x ) ϖ ¯ 11 φ ( x ) u e c B ( x ) ρ ( x ) ϖ ¯ 12 φ ( x ) v e c B ( x ) ρ ( x ) ] = 0 , i n Ω , e c B ( x ) ρ ( x ) [ μ 2 ( x ) e c B ( x ) v ] + v [ λ ¯ 2 φ ( x ) e c B ( x ) ρ ( x ) + ε 2 g 2 ( x ) ϖ ¯ 21 φ ( x ) u e c B ( x ) ρ ( x ) ϖ ¯ 22 φ ( x ) v e c B ( x ) ρ ( x ) ] = 0 , i n Ω , u n = v n = 0 , o n Ω .
Set k ¯ i = max Ω ¯ g i ( x ) φ ( x ) e c B ( x ) ρ ( x ) , k ̲ i = min Ω ¯ g i ( x ) φ ( x ) e c B ( x ) ρ ( x ) for i = 1 , 2 . Applying 0 < ε i 1 and ϖ ¯ 21 ϖ ¯ 11 < λ ¯ 2 λ ¯ 1 < ϖ ¯ 22 ϖ ¯ 12 , we have the linear system
λ ¯ 1 + ε 1 k ̲ 1 ϖ ¯ 11 u ̲ ϖ ¯ 12 v ¯ = 0 , λ ¯ 2 + ε 2 k ¯ 2 ϖ ¯ 21 u ̲ ϖ ¯ 22 v ¯ = 0 , λ ¯ 1 + ε 1 k ¯ 1 ϖ ¯ 11 u ¯ ϖ ¯ 12 v ̲ = 0 , λ ¯ 2 + ε 2 k ̲ 2 ϖ ¯ 21 u ¯ ϖ ¯ 22 v ̲ = 0 .
Then
u ¯ = ϖ ¯ 22 ( λ ¯ 1 + ε 1 k ¯ 1 ) ϖ ¯ 12 ( λ ¯ 2 + ε 2 k ̲ 2 ) ϖ ¯ 11 ϖ ¯ 22 ϖ ¯ 12 ϖ ¯ 21 , u ̲ = ϖ ¯ 22 ( λ ¯ 1 + ε 1 k ̲ 1 ) ϖ ¯ 12 ( λ ¯ 2 + ε 2 k ¯ 2 ) ϖ ¯ 11 ϖ ¯ 22 ϖ ¯ 12 ϖ ¯ 21 ,
v ¯ = ϖ ¯ 11 ( λ ¯ 2 + ε 2 k ¯ 2 ) ϖ ¯ 21 ( λ ¯ 1 + ε 1 k ̲ 1 ) ϖ ¯ 11 ϖ ¯ 22 ϖ ¯ 12 ϖ ¯ 21 , v ̲ = ϖ ¯ 11 ( λ ¯ 2 + ε 2 k ̲ 2 ) ϖ ¯ 21 ( λ ¯ 1 + ε 1 k ¯ 1 ) ϖ ¯ 11 ϖ ¯ 22 ϖ ¯ 12 ϖ ¯ 21 .
Hence, the system (25) has a positive non-homogeneous steady-state ( u θ ˜ ( x ) , v θ ˜ ( x ) ) and 0 < u ̲ < u θ ˜ ( x ) < u ¯ and 0 < v ̲ < v θ ˜ ( x ) < v ¯ . Let
η 1 = u ̲ v ¯ , η 2 = u ¯ v ̲ .
we have η 1 u θ ˜ ( x ) v θ ˜ ( x ) η 2 . Applying (28), we get lim ε 1 , ε 2 0 η 2 η 1 = 1 . Hence, for 0 < ε i 1 ,
min Ω ¯ ϖ ¯ 11 φ ( x ) ϖ ¯ 22 φ ( x ) ϖ ¯ 12 φ ( x ) ϖ ¯ 21 φ ( x ) = ϖ ¯ 11 ϖ ¯ 22 ϖ ¯ 12 ϖ ¯ 21 > η 2 η 1 .
The proof is completed. □
Example 1.
In the above (25), let c = 2 , B ( x ) = x , ρ ( x ) = e x , μ 1 ( x ) = μ 2 ( x ) = e x , R 1 ( x ) = R 2 ( x ) = 1 2 e x , φ ( x ) = e x , g 1 ( x ) = g 2 ( x ) = 1 + cos ( π 2 x ) , λ ¯ 1 = 1 , λ ¯ 2 = 2 , ϖ ¯ 11 = ϖ ¯ 12 = ϖ ¯ 21 = 1 , ϖ ¯ 22 = 3 , and ε 1 = ε 2 = 1 3 , x Ω = [ 0 , 10 ] . Then the problem (25) becomes the following model
u t = e x [ e x u ] + u [ 1 + 1 3 ( 1 + cos ( π 2 x ) ) u v ] , i n Ω × R + , v t = e x [ e x v ] + v [ 2 + 1 3 ( 1 + cos ( π 2 x ) ) u 3 v ] , i n Ω × R + , u n = v n = 0 , o n Ω × R + , u ( x , 0 ) = e x ( 2 + cos ( π x ) ) , 0 , v ( x , 0 ) = e x ( 2 + cos ( π x ) ) , 0 , i n Ω ,
where u 0 ( x ) , v 0 ( x ) , 0 . It is not difficult to verify that ( H 1 ) and ( H 2 ) hold. We can find η 1 = 1 > 0 , η 2 = 7 3 > 0 , s u c h t h a t η 1 u θ ˜ ( x ) v θ ˜ ( x ) η 2 a n d η 2 η 1 < min Ω ¯ ϖ 11 ϖ 22 ϖ 12 ϖ 21 . According to Theorem 2, the model (29) admits a solution ( u ( x , t ) , v ( x , t ) ) that satisfies
lim t u ( x , t ) = u θ ˜ ( x ) , lim t v ( x , t ) = v θ ˜ ( x ) i n C 2 ( Ω ¯ ) .
Indeed, the steady-state of (29) satisfies the following elliptic problem
e x [ e x u ] + u [ 1 + 1 3 ( 1 + cos ( π 2 x ) ) u v ] = 0 , i n Ω , e x [ e x v ] + v [ 2 + 1 3 ( 1 + cos ( π 2 x ) ) u 3 v ] = 0 , i n Ω , u n = v n = 0 , o n Ω .
It is not difficult to see that k ¯ 1 = k ¯ 2 = 2 , k ̲ 1 = k ̲ 2 = 0 . By calculation, we can obtain
1 u ̲ v ¯ = 0 , 2 + 2 ε 2 u ̲ 3 v ¯ = 0 , 1 + 2 ε 1 u ¯ v ̲ = 0 , 2 u ¯ 3 v ̲ = 0 .
Then
u ¯ = 1 + 6 ε 1 2 = 3 2 > 0 , u ̲ = 1 2 ε 2 2 = 1 6 > 0 ,
v ¯ = 1 + 2 ε 2 2 = 5 6 > 0 , v ̲ = 1 2 ε 1 2 = 1 6 > 0 .
Hence, 0 < u ̲ < u θ ˜ ( x ) < u ¯ and 0 < v ̲ < v θ ˜ ( x ) < v ¯ , which yield that there exists a positive non-homogeneous steady-state ( u θ ˜ ( x ) , v θ ˜ ( x ) ) of (29).

5. Discussion

In this paper, by using the Lyapunov functional method, we mainly analyzed the global stability of non-homogeneous steady-state for the Lotka–Volterra competition–diffusion–advection system between two competing biological organisms in heterogeneous environments, where two biological organisms are competing for different fundamental resources, their advection and diffusion strategies follow different positive diffusion distributions, and the functions of specific competition ability are variable. Moreover, we also obtained the global stability result when one of the two organisms has no diffusion ability and is not affected by advection.
At the end of this section, we propose an interesting research problem. To the best of our knowledge, for the Lotka–Volterra competition–diffusion–advection system between two competing biological organisms in heterogeneous environments, we did not obtain any results under the condition of cross-diffusion, such as the existence and stability of nontrivial positive steady state. We leave this challenge to future investigations.

Author Contributions

All authors contributed equally and significantly in writing this article. Conceptualization, L.C.; Formal analysis, Y.Z.; Funding acquisition, L.C.; Writing—original draft, S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Shandong Provincial Natural Science Foundation under grant ZR2020MA006 and the Introduction and Cultivation Project of Young and Innovative Talents in Universities of Shandong Province.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to express our thanks to the anonymous referees and the editor for their constructive comments and suggestions, which greatly improved this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cantrell, R.; Cosner, C. Spatial Ecology via Reaction-Diffusion Equations; Wiley: New York, NY, USA, 2003. [Google Scholar]
  2. He, X.; Ni, W. Global dynamics of the Lotka-Volterra competition-diffusion system: Diffusion and spatial heterogeneity I. Commun. Pure Appl. Math. 2016, 69, 981–1014. [Google Scholar] [CrossRef]
  3. Korobenko, L.; Braverman, E. On evolutionary stability of carrying capacity driven dispersal in competition with regularly diffusing populations. J. Math. Biol. 2014, 69, 1181–1206. [Google Scholar] [CrossRef] [PubMed]
  4. Hsu, S.; Waltman, P. On a system of reaction-diffusion equations arising from competition in an unstirred chemostat. SIAM J. Appl. Math. 1993, 53, 1026–1044. [Google Scholar] [CrossRef] [Green Version]
  5. Lam, K.; Ni, W. Uniqueness and complete dynamics in heterogeneous competition–diffusion systems. SIAM J. Appl. Math. 2012, 72, 1695–1712. [Google Scholar] [CrossRef] [Green Version]
  6. Ni, W.; Shi, J.; Wang, M. Global stability of nonhomogeneous equilibrium solution for the diffusive Lotka-Volterra competition model. Calc. Var. Partial Differ. Equ. 2020, 59, 267–270. [Google Scholar] [CrossRef]
  7. Cosner, C. Reaction-diffusion-advection models for the effects and evolution of dispersal. Discret. Contin. Dyn. Syst. 2014, 34, 1701–1745. [Google Scholar] [CrossRef]
  8. Lou, Y.; Zhao, X.; Zhou, P. Global dynamics of a Lotka-Volterra competition-diffusion-advection system in heterogeneous environments. J. Math. Pures Appl. 2019, 121, 47–82. [Google Scholar] [CrossRef]
  9. Zhou, P.; Xiao, D. Global dynamics of a classical Lotka-Volterra competition-diffusion-advection system. J. Funct. Anal. 2018, 275, 356–380. [Google Scholar] [CrossRef]
  10. Wang, Q. On a Lotka-Volterra competition-diffusion-advection model in general heterogeneous environments. J. Math. Anal. Appl. 2020, 489, 124127. [Google Scholar] [CrossRef]
  11. Li, Z.; Dai, B. Stability and Hopf bifurcation analysis in a Lotka–Volterra competition–diffusion–advection model with time delay effect. Nonlinearity 2021, 34, 3271–3313. [Google Scholar] [CrossRef]
  12. Li, Z.; Dai, B.; Dong, X. Global stability of nonhomogeneous steady-state solution in a Lotka-Volterra competition-diffusion-advection model. Appl. Math. Lett. 2020, 107, 106480. [Google Scholar] [CrossRef]
  13. Ma, L.; Guo, S. Bifurcation and stability of a two-species reaction-diffusion-advection competition model. Nonlinear Anal-Real 2021, 59, 103241. [Google Scholar] [CrossRef]
  14. Zhou, P.; Tang, D.; Xiao, D. On Lotka-Volterra competitive parabolic systems: Exclusion, coexistence and bistability. J. Differ. Equ. 2021, 282, 596–625. [Google Scholar] [CrossRef]
  15. Wang, M. A diffusive logistic equation with a free boundary and sign-changing coefficient in time-periodic environment. J. Funct. Anal. 2016, 270, 483–508. [Google Scholar] [CrossRef]
  16. Wang, M.; Zhang, Y. Dynamics for a diffusive prey-predator model with different free boundaries. J. Differ. Equ. 2018, 264, 3527–3558. [Google Scholar] [CrossRef] [Green Version]
  17. Wang, M. Note on the Lyapunov functional method. Appl. Math. Lett. 2018, 75, 102–107. [Google Scholar] [CrossRef]
  18. Pao, C. Nonlinear Parabolic and Elliptic Equations; Plenum Press: New York, NY, USA, 1992. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, L.; Lin, S.; Zhao, Y. Global Stability of a Lotka-Volterra Competition-Diffusion-Advection System with Different Positive Diffusion Distributions. Axioms 2021, 10, 166. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms10030166

AMA Style

Chen L, Lin S, Zhao Y. Global Stability of a Lotka-Volterra Competition-Diffusion-Advection System with Different Positive Diffusion Distributions. Axioms. 2021; 10(3):166. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms10030166

Chicago/Turabian Style

Chen, Lili, Shilei Lin, and Yanfeng Zhao. 2021. "Global Stability of a Lotka-Volterra Competition-Diffusion-Advection System with Different Positive Diffusion Distributions" Axioms 10, no. 3: 166. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms10030166

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