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:
Here, and are the population densities of biological organisms, location time , which are supposed to be nonnegative. correspond to the dispersal rates of two competing biological organisms, respectively. correspond to the advection rates of two competing biological organisms, and are the nonconstant functions and represent the advective directions. Two bounded functions and are the intrinsic growth rates of competing organisms, , are two positive diffusion distributions, respectively. show the strength of competition ability. The spatial habitat is a bounded smooth domain, ; 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. , , .
Hypothesis 2. To simplify the calculation, by letting
,
, the system (
1) converts into the following coupled system
when
, the model (
2) has been studied in Ni et al. [
6].
, 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:
where
is bounded and
is a smooth boundary. The initial condition
.
Setting the following assumptions:
Let
,
, such that
Let
be a constant, such that
for some
and there is
such that
and there exists
, satisfying
The following lemma (see [
15,
16]) is the boundedness result of the solution
in (
3).
Lemma 1. Let be a solution of (3) with , Suppose that f, satisfy the assumptions , then for any , there is a constant such that 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 be constants, in . Assume that has lower bound, in . If one of the following alternatives holds: and in for ,
and for and ,
where P and m are constants, then .
Lemma 3 ([
6]).
Let with and , with , are functions. If the following conditions holds: is a constant, the function , , is a non-increasing function for ,
, on ,
Next, we consider the following scalar evolution eqution
where
satisty
Now we see the following useful lemma.
Lemma 4 ([
1]).
Assume that on , then the elliptic problem:has a unique positive solution, denoted by . 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
of (
7).
Theorem 1. Assume that . If satisfy (6), then Equation (5) has a unique solution with in . Proof. According to the upper–lower solutions method [
1,
18], we obtain (
5) with a unique solution
. Let
M be a upper solution of (
5), we have
,
.
By applying Lemma 1, we can obtain that there exists a constant
such that
Then, define a function
by
Then,
,
. By (
2) and (
4), we have
We get
By virtue of (
8), we get
in
for some
. From Lemma 2, it follows that
Applying (
8) again,
is relatively compact in
. It can be found that there exists some function
such that
Combining with (
12), we get
where
. Hence, we deduce
□
In addition, taking advantage of Lyapunov function method, the global stability results of (
2) are obtained.
Theorem 2. Suppose that , and hold, the system (2) admits a non-homogeneous steady-state and there exists Then, the system (2) admits a solution that satisfies Proof. Assume that the inequality (
15) holds, let
defined by
where
. Clearly,
. By (
2) and (
4), we have
Note that (
14) and (
15) give rise to
Choosing
, we have
Combining with (
17), we can deduce
From (
13), it follows that
□
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 satisfy and on . Let for , and Ifandthen there is a unique non-homogeneous steady-state for the model (2) such that Ifthen there exists a semi-trivial steady-state for the model (2) such that Letthen the model (2) has a semi-trivial steady-state ,where . Proof. When
,
,
of the model (
2) satisfies
and
.
If (
18) and (
19) hold, we see
, then by Lemma 4, the problem (
23) has a unique solution
. By using the maximum principle in elliptic equation, we infer
According to (
20), we can get
, hence there exists a unique steady-state for (
2),
.
Let us define a function
,
where
. Clearly,
. From (
2) and (
4), we get
We can choose
and use (18), such that
Combining this with (24), we can deduce
Applying the Lemma 1 and Sobolev embedding theorem, we deduce that
u and
v are bounded in
and there is a constant
such that
Combining with (
2) and
in
for some
, and making use of Lemma 2, we get
and we deduce that
Applying Theorem 2, we get in .
Let’s define a function
,
where
. From (
4) and (
21), we have
The following discussion will refer to the part , then we will not repeat it.
Clearly, (
2) has a semi-trivial steady-state
. Let us define a function
,
where
and
. From (
22), we have
The following discussion is similar to the part , so we omit it. □
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.