In this paper, a nonlinear diffusion SI epidemic model with a general incidence rate in heterogeneous environment is studied. Global behavior of classical solutions under certain restrictions on the coefficients is considered. We first establish the global existence of classical solutions of the system under heterogeneous environment by energy estimate and maximum principles. Based on such estimates, we then study the large-time behavior of the solution of system under homogeneous environment. The model and mathematical results in [M. Kirane, S. Kouachi, Global solutions to a system of strongly coupled reaction-diffusion equations, Nonlinear Anal., 26 (1996), 1387-1396.] are generalized.
Citation: Shenghu Xu, Xiaojuan Li. Global behavior of solutions to an SI epidemic model with nonlinear diffusion in heterogeneous environment[J]. AIMS Mathematics, 2022, 7(4): 6779-6791. doi: 10.3934/math.2022377
[1] | Dazhao Chen . Endpoint estimates for multilinear fractional singular integral operators on Herz and Herz type Hardy spaces. AIMS Mathematics, 2021, 6(5): 4989-4999. doi: 10.3934/math.2021293 |
[2] | Anas A. Hijab, Laith K. Shaakir, Sarah Aljohani, Nabil Mlaiki . Double composed metric-like spaces via some fixed point theorems. AIMS Mathematics, 2024, 9(10): 27205-27219. doi: 10.3934/math.20241322 |
[3] | Xianyong Huang, Shanhe Wu, Bicheng Yang . A Hardy-Hilbert-type inequality involving modified weight coefficients and partial sums. AIMS Mathematics, 2022, 7(4): 6294-6310. doi: 10.3934/math.2022350 |
[4] | M. Zakarya, Ghada AlNemer, A. I. Saied, H. M. Rezk . Novel generalized inequalities involving a general Hardy operator with multiple variables and general kernels on time scales. AIMS Mathematics, 2024, 9(8): 21414-21432. doi: 10.3934/math.20241040 |
[5] | Limin Yang, Ruiyun Yang . Some new Hardy-Hilbert-type inequalities with multiparameters. AIMS Mathematics, 2022, 7(1): 840-854. doi: 10.3934/math.2022050 |
[6] | Suriyakamol Thongjob, Kamsing Nonlaopon, Sortiris K. Ntouyas . Some (p, q)-Hardy type inequalities for (p, q)-integrable functions. AIMS Mathematics, 2021, 6(1): 77-89. doi: 10.3934/math.2021006 |
[7] | Chunyan Luo, Yuping Yu, Tingsong Du . Estimates of bounds on the weighted Simpson type inequality and their applications. AIMS Mathematics, 2020, 5(5): 4644-4661. doi: 10.3934/math.2020298 |
[8] | Sabir Hussain, Javairiya Khalid, Yu Ming Chu . Some generalized fractional integral Simpson’s type inequalities with applications. AIMS Mathematics, 2020, 5(6): 5859-5883. doi: 10.3934/math.2020375 |
[9] | Irshad Ayoob, Ng Zhen Chuan, Nabil Mlaiki . Hardy-Rogers type contraction in double controlled metric-like spaces. AIMS Mathematics, 2023, 8(6): 13623-13636. doi: 10.3934/math.2023691 |
[10] | Qing Yang, Chuanzhi Bai . Fixed point theorem for orthogonal contraction of Hardy-Rogers-type mapping on O-complete metric spaces. AIMS Mathematics, 2020, 5(6): 5734-5742. doi: 10.3934/math.2020368 |
In this paper, a nonlinear diffusion SI epidemic model with a general incidence rate in heterogeneous environment is studied. Global behavior of classical solutions under certain restrictions on the coefficients is considered. We first establish the global existence of classical solutions of the system under heterogeneous environment by energy estimate and maximum principles. Based on such estimates, we then study the large-time behavior of the solution of system under homogeneous environment. The model and mathematical results in [M. Kirane, S. Kouachi, Global solutions to a system of strongly coupled reaction-diffusion equations, Nonlinear Anal., 26 (1996), 1387-1396.] are generalized.
Hardy type inequalities play crucial roles in analysis, probability, and partial differential equations. We first recall the classical L2- Hardy inequality:
∫Rn|∇ϕ|2dx≥(n−22)2∫Rn|ϕ|2|x|2dx, |
for n≥3 and ϕ∈C∞0(Rn∖{0}). The constant (n−22)2 is sharp and is never attained by nontrivial functions. The Rellich inequality is a natural generalization of the above Hardy inequality. In Rn, it reads as follows:
∫Rn|Δϕ|2dx≥(n(n−4)4)2∫Rn|ϕ|2|x|4dx, |
where n≥5 and ϕ∈C∞0(Rn∖{0}).
In recent years, there have been many results in the literature on the Hardy and Rellich type inequalities in the context of a complete Riemannian manifold. In particular, the following Hardy inequality has been first established on Riemannian manifold (M,g) by Carron [1]:
∫Mρα|∇ϕ|2dvg≥(C+α−1)24∫Mρα−2ϕ2dvg, |
where α∈R, C+α−1>0, ϕ∈C∞0(M∖ρ−1{0}). The weight function ρ is nonnegative and it satisfies |∇ρ|=1 and Δρ≥Cρ in the sense of distribution. Here, dvg, ∇, and Δ denote the volume element, the gradient, and the Laplace operator on M. Under the same geometric assumptions on the weight function ρ, Kombe and Özaydin [2] established the weighted Lp-Hardy inequality:
∫Mρα|∇ϕ|pdvg≥(C+1+α−pp)p∫Mρα−pϕpdvg, |
where 1≤p<∞, C+1+α−p>0, and ϕ∈C∞0(M∖ρ−1{0}). Kombe and Özaydin [3] also proved a new weighted Hardy-Poincaré inequality. They showed that if M is a complete non-compact Riemannian manifold of dimension n>1 and ρ is a nonnegative function on M such that |∇ρ|=1 and Δρ≥Cρ in the sense of distribution, where C>0, the following inequality holds:
∫Mρα+p|∇ρ⋅∇ϕ|pdvg≥(C+1+αp)p∫Mρα|ϕ|pdvg. |
Xia [4] proved the following Hardy type inequality on a complete non-compact Riemannian manifold. Let M be an n-dimensional complete non-compact Riemannian manifold, where n≥2, and let ρ be a nonnegative function on M such that |∇ρ|=1 and Δρ≥Cρ+H in the sense of distribution, where H is a continuous function on M and C is a constant. The result is the following: For any p,q∈R with q>1+C and any compactly supported smooth function ϕ∈C∞0(M∖ρ−1{0}), the following inequality holds:
(q−C−1)∫M|ϕ|pρqdvg≤|p|∫M|ϕ|p−1ρq−1|∇ϕ|dvg+∫MH|ϕ|pρq−1dvg. |
Huang and Ye [5] considered the first order Hardy inequalities using simple identities. This basic setting not only permits to derive quickly many well-known Hardy inequalities with optimal constants, but also supplies improved or new estimates in miscellaneous situations. We also refer the interested reader to [6,7,8,9,10,11,12,13,14,15,16], which are excellent monographs on the topic.
In this paper, we are interested in proving some Hardy type identities and inequalities on the smooth metric measure spaces related to the divergence type operator Lf,V. Before that, we would like to briefly introduce the smooth metric measure spaces and the divergence type operator Lf,V.
A smooth metric measure space (Mn,g,dμ) is a Riemannian manifold (Mn,g) equipped with a conformal Riemannian volume dμ=1Ve−fdvg, where dvg denotes the Riemannian volume measure on M, V is a positive twice differentiable function on M, and f is a real-valued smooth function on M.
On a smooth metric measure space (Mn,g,dμ), we can define the weighted Ricci curvature ^RicVf,m [17] given by
^RicVf,m=ΔfVVg−1V∇2V+Ricmf, |
where
Δf⋅:=efdiv(e−f∇⋅)=Δ⋅−∇f∇, | (1.1) |
is the f-Laplacian (see [18,19,20,21]),
Ricmf=Ric+∇2f−1m−n∇f⊗∇f, |
is the m-Bakry-Émery Ricci curvature, and m is a real constant. Here, ∇, Δ, ∇2, and div are the gradient, Laplacian, Hessian, and divergence operator on M with respect to g, respectively.
We introduce the divergence type operator Lf,V on M as follows (see [17,22] for detail):
Lf,Vφ=efdiv[e−fV2∇(φV)], | (1.2) |
where φ is a smooth function and V is a positive twice differentiable function defined on M. In general, the divergence type operator Lf,V is not self-adjoint with respect to the standard L2-inner product because of the first order term, but it is self-adjoint with respect to the weighted measure dμ=1Ve−fdvg. That is, for any smooth functions u,v∈C∞0(M), we have
∫MuLf,Vvdμ=−∫MV3∇uV∇vVdμ=∫MvLf,Vudμ. | (1.3) |
We observe that the divergence type operator Lf,V encompasses, as very special cases, many differential linear operators. In the case that V=1, i.e., Lf,1φ=Δfφ. In the case that f is constant, i.e., Lf,Vφ=LVφ=VΔφ−φΔV, we call LV by the generalized Schrödinger operator. From the viewpoint of geometry, the generalized Schrödinger operator LV plays an important role in the geometric understanding of the sub-static manifolds; see [23].
Du and Mao [24] proved that some Hardy and Rellich type inequalities on the smooth metric measure spaces related to the f-Laplace. Also, Li, Abolarinwa, Alkhaldi, and Ali [25] generalized some integral inequalities of Hardy type to the setting of the smooth metric measure spaces. These studies are without the curvature conditions. On the other hand, the case of additional curvature conditions has been studied. Kolesnikov and Milman [26] proved the Hardy-Poincaré inequality under the curvature condition Ricmf>0 on M:
mm−1∫Mφ2e−fdvg≤∫Me−f⟨(Ricmf)−1∇φ,∇φ⟩dvg, |
where 1m∈(−∞,1n] and φ∈C1(M). Huang and Zhu [22] gave the following Hardy-Poincaré inequality: Let M be a compact Riemannian manifold and ^RicVf,m>0, then
mm−1∫MVφ2e−fdvg≤∫Me−fV⟨(^RicVf,m)−1∇φ,∇φ⟩dvg, |
where φ∈C∞(M). Huang and Zhu [22] studied weighted L2-Hardy-Poincaré inequalities on a smooth metric measure space related to the divergence type operator Lf,V under the curvature condition.
The primary objectives of the present paper are twofold: First, we want to establish some L2 Hardy type identities and Lp Hardy type inequalities related to the divergence type operator Lf,V on the smooth metric measure spaces without curvature conditions. Second, as application we would like to show some Rellich type inequalities for the divergence type operator Lf,V.
The remainder of the paper is as follows: In Section 2, we will prove Hardy type identity related to the divergence type operator Lf,V and several related corollaries. Then, we will prove some Lp Hardy type inequalities in Section 3. In the last section, we will prove a Rellich type inequality related to the divergence type operator Lf,V.
In this section, we will prove some Hardy type identities related to the divergence type operator Lf,V and several related corollaries.
Theorem 2.1. Let (Mn,g,dμ) be an n-dimensional (n≥2) complete non-compact smooth metric measure space. Suppose that V is a positive twice differentiable function and f is a real-valued smooth function on M and W∈C1(M). Then, for any compactly supported smooth functions uV∈C10(M) and φV∈C20(M), the following identity holds:
∫MWV3|∇(uV)|2dμ=∫MWV3|∇(uV)−(uV)(φV)−1∇(φV)|2dμ−∫MefV(uV)2(φV)−1div(e−fV2W∇(φV))dμ, | (2.1) |
where dμ=1Ve−fdvg and dvg is the Riemannian volume element related to g.
Proof. As a consequence of the integration by parts (1.3), we have
∫MWV3|∇(uV)−(uV)(φV)−1∇(φV)|2dμ=∫Me−fV2W|∇(uV)|2dvg−2∫Me−fV2W⟨∇(uV),(uV)(φV)−1∇(φV)⟩dvg+∫Me−fV2W(uV)2(φV)−2|∇(φV)|2dvg=∫Me−fV2W|∇(uV)|2dvg−∫Me−fV2W⟨(φV)−1∇(uV)2,∇(φV)⟩dvg+∫Me−fV2W(uV)2(φV)−2|∇(φV)|2dvg=∫Me−fV2W|∇(uV)|2dvg+∫M(uV)2div(e−fV2W(φV)−1∇(φV))dvg+∫Me−fV2W(uV)2(φV)−2|∇(φV)|2dvg=∫Me−fV2W|∇(uV)|2dvg+∫M(φV)−1(uV)2div(e−fV2W∇(φV))dvg−∫M(φV)−2(uV)2⟨∇(φV),e−fV2W∇(φV)⟩dvg+∫Me−fV2W(uV)2(φV)−2|∇(φV)|2dvg=∫MWV3|∇(uV)|2dμ+∫MefV(uV)2(φV)−1div(e−fV2W∇(φV))dμ. |
Then, we can get
∫MWV3|∇(uV)|2dμ=∫MWV3|∇(uV)−(uV)(φV)−1∇(φV)|2dμ−∫MefV(uV)2(φV)−1div(e−fV2W∇(φV))dμ. |
This completes the proof of Theorem 2.1.
As a special case for W=1 in Theorem 2.1, we have the following.
Corollary 2.2. Let (Mn,g,dμ) be an n-dimensional (n≥2) complete non-compact smooth metric measure space. Suppose that V is a positive twice differentiable function and f is a real-valued smooth function on M. Then, for any uV∈C10(M) and φV∈C20(M), the following identity holds:
∫MV3|∇(uV)|2dμ=∫MV3|∇(uV)−(uV)(φV)−1∇(φV)|2dμ−∫MV(uV)2(φV)−1(Lf,Vφ)dμ, | (2.2) |
where dμ=1Ve−fdvg and dvg is the Riemannian volume element related to g.
By using Corollary 2.2, we have the following.
Corollary 2.3. Let (Mn,g,dμ) be an n-dimensional (n≥2) complete non-compact smooth metric measure space. Suppose that V is a positive twice differentiable function and f is a real-valued smooth function on M. Then, for any uV∈C10(M) and φV∈C20(M), the following inequality holds:
∫MV3|∇(uV)|2dμ≥−∫MV(uV)2(φV)−1(Lf,Vφ)dμ, | (2.3) |
where dμ=1Ve−fdvg and dvg is the Riemannian volume element related to g.
As the special case that φV=(ρV)1−C2 in the Corollary 2.3, we have the following.
Corollary 2.4. Let (Mn,g,dμ) be an n-dimensional (n≥2) complete non-compact smooth metric measure space. Suppose that V is a positive twice differentiable function and f is a real-valued smooth function on M. Let ρ be a nonnegative function on M such that |∇(ρV)|=1 and Lf,Vρ≥C(ρV)−1V2 in the sense of distribution, where C is a constant and C>1. Then, for any uV∈C10(M∖(ρV)−1{0}), the following inequality holds:
∫MV3|∇(uV)|2dμ≥(C−1)24∫MV3(uV)2(ρV)−2dμ, |
where dμ=1Ve−fdvg and dvg is the Riemannian volume element related to g.
Proof. Taking φV=(ρV)1−C2, we have
−(φV)−1(Lf,Vφ)=−ef(ρV)C−12div(e−fV2∇(ρV)1−C2)=−(ρV)C−12efdiv((1−C2)e−fV2(ρV)−C−12∇(ρV))=−(1−C2)(ρV)C−12ef⟨∇(ρV)−C−12,e−fV2∇(ρV)⟩−(1−C2)(ρV)C−12(ρV)−C−12efdiv(e−fV2∇(ρV))=1−C24V2(ρV)−2|∇(ρV)|2+C−12(ρV)−1Lf,Vρ. |
Then, we have
−(φV)−1(Lf,Vφ)=1−C24V2(ρV)−2|∇(ρV)|2+C−12(ρV)−1Lf,Vρ≥1−C24V2(ρV)−2+C2−C2V2(ρV)−2=C2−2C+14V2(ρV)−2=(C−1)24V2(ρV)−2, |
where we use the assumption |∇(ρV)|=1 and Lf,Vρ≥C(ρV)−1V2. As a result, we can get
−(φV)−1(Lf,Vφ)≥(C−1)24V2(ρV)−2. | (2.4) |
Substituting (2.4) into (2.3), we can get
∫MV3|∇(uV)|2dμ≥(C−1)24∫MV3(uV)2(ρV)−2dμ. |
This completes the proof of Corollary 2.4.
Remark 2.5. In the special case that f is a constant and V=1, Corollary 2.4 reduces to Carron's result in [1, Proposition 2.1].
In this section, we will prove some Lp-Hardy type inequalities on smooth metric measure spaces. Our first result is the following.
Theorem 3.1. Let (Mn,g,dμ) be an n-dimensional (n≥2) complete non-compact smooth metric measure space. Suppose that V is a positive twice differentiable function and f is a real-valued smooth function on M. Let ρ be a nonnegative function on M such that |∇(ρV)|=1 and Lf,Vρ≤(H+C(ρV)−1)V2 in the sense of distribution, where H is a continuous function on M and C is a constant. Then, for any p,q∈R,q>1+C, and any compactly supported smooth functions ϕV∈C∞0(M∖(ρV)−1{0}), the following inequality holds:
(q−C−1)∫MV3|ϕV|p(ρV)−qdμ≤|p|∫MV3|ϕV|p−1(ρV)1−q|∇(ϕV)|dμ+∫MHV3|ϕV|p(ρV)1−qdμ. | (3.1) |
Proof. For a vector field X on M, we denote by div X the divergence of X. Note that
div(e−fV2|ϕV|p(ρV)1−q∇(ρV))=pe−fV2|ϕV|p−1⟨∇|ϕV|,(ρV)1−q∇(ρV)⟩+|ϕV|pdiv(e−fV2(ρV)1−q∇(ρV))=pe−fV2|ϕV|p−1⟨∇|ϕV|,(ρV)1−q∇(ρV)⟩+|ϕV|p(ρV)1−qdiv(e−fV2∇(ρV))+(1−q)e−fV2|ϕV|p⟨(ρV)−q∇(ρV),∇(ρV)⟩=pe−fV2|ϕV|p−1⟨∇|ϕV|,(ρV)1−q∇(ρV)⟩+e−f|ϕV|p(ρV)1−q(Lf,Vρ)+(1−q)e−fV2|ϕV|p(ρV)−q|∇(ρV)|2. |
By direct computation, we have
div(e−fV2|ϕV|p(ρV)1−q∇(ρV))=pe−fV2|ϕV|p−1⟨∇|ϕV|,(ρV)1−q∇(ρV)⟩+e−f|ϕV|p(ρV)1−q(Lf,Vρ)+(1−q)e−fV2|ϕV|p(ρV)−q|∇(ρV)|2≤pe−fV2|ϕV|p−1⟨∇|ϕV|,(ρV)1−q∇(ρV)⟩+e−fV2|ϕV|p(ρV)1−q(H+C(ρV)−1)+(1−q)e−fV2|ϕV|p(ρV)−q=pe−fV2|ϕV|p−1⟨∇|ϕV|,(ρV)1−q∇(ρV)⟩+(1+C−q)e−fV2|ϕV|p(ρV)−q+e−fV2H|ϕV|p(ρV)1−q, | (3.2) |
where we use the assumption |∇(ρV)|=1 and Lf,Vρ≤(H+C(ρV)−1)V2.
Since q>1+C and |∇(ρV)|=1, we have
(q−C−1)∫Me−fV2|ϕV|p(ρV)−qdvg≤p∫Me−fV2|ϕV|p−1⟨∇|ϕV|,(ρV)1−q∇(ρV)⟩dvg+∫Me−fV2H|ϕV|p(ρV)1−qdvg≤|p∫Me−fV2|ϕV|p−1⟨∇|ϕV|,(ρV)1−q∇(ρV)⟩dvg|+∫Me−fV2H|ϕV|p(ρV)1−qdvg≤|p|∫Me−fV2|ϕV|p−1(ρV)1−q|∇(ϕV)||∇(ρV)|dvg+∫Me−fV2H|ϕV|p(ρV)1−qdvg=|p|∫Me−fV2|ϕV|p−1(ρV)1−q|∇(ϕV)|dvg+∫Me−fV2H|ϕV|p(ρV)1−qdvg. |
Then, we have
(q−C−1)∫MV3|ϕV|p(ρV)−qdμ≤|p|∫MV3|ϕV|p−1(ρV)1−q|∇(ϕV)|dμ+∫MHV3|ϕV|p(ρV)1−qdμ. | (3.3) |
This completes the proof of Theorem 3.1.
Remark 3.2. In the special case that V=1 and f is a constant, the inequality (3.1) reduces to Xia's result in [4, Theorem 2.1].
Theorem 3.3. Let (Mn,g,dμ) be an n-dimensional (n≥2) complete non-compact smooth metric measure space. Suppose that V is a positive twice differentiable function and f is a real-valued smooth function on M. Let ρ be a nonnegative function on M such that |∇(ρV)|=1 in the sense of distributions. Then, for any p,q∈R,1<p<+∞,0≤q≤p, and any compactly supported smooth function ϕV∈C∞0(M∖(ρV)−1{0}), we have
(i) When Lf,Vρ≤C(ρV)−1V2 in the sense of distributions, where C<q−1 is a constant, the following inequality holds:
∫MV3|ϕV|p(ρV)−qdμ≤(pq−C−1)q(∫MV3|ϕV|pdμ)p−qp(∫MV3|∇(ϕV)|pdμ)qp. | (3.4) |
(ii) When Lf,Vρ≥C(ρV)−1V2 in the sense of distributions, where C>q−1 is a constant, the following inequality holds:
∫MV3|ϕV|p(ρV)−qdμ≥(pC+1−q)q(∫MV3|ϕV|pdμ)p−qp(∫MV3|∇(ϕV)|pdμ)qp. | (3.5) |
Proof. (ⅰ) In the special case that H=0 in (3.1), we can get
(q−C−1)∫MV3|ϕV|p(ρV)−qdμ≤|p|∫MV3|ϕV|p−1(ρV)1−q|∇(ϕV)|dμ. | (3.6) |
It follows from the Hölder inequality that
∫Me−fV2|ϕV|p−1(ρV)1−q|∇(ϕV)|dvg=∫M(e−fV2)1p(e−fV2)p−1p|ϕV|p−1(ρV)1−q|∇(ϕV)|dvg≤(∫Me−fV2|∇(ϕV)|pdvg)1p(∫Me−fV2|ϕV|p(ρV)p(1−q)p−1dvg)p−1p≤(∫Me−fV2|∇(ϕV)|pdvg)1p(∫Me−fV2|ϕV|p(ρV)−qdvg)q−1q(∫Me−fV2|ϕV|pdvg)(1−p(q−1)q(p−1))p−1p. | (3.7) |
Substituting (3.7) into (3.6), we can get
(q−C−1)∫Me−fV2|ϕV|p(ρV)−qdvg≤|p|(∫Me−fV2|∇(ϕV)|pdvg)1p(∫Me−fV2|ϕV|p(ρV)−qdvg)q−1q(∫Me−fV2|ϕV|pdvg)(1−p(q−1)q(p−1))p−1p≤|p|(∫Me−fV2|∇(ϕV)|pdvg)1p(∫Me−fV2|ϕV|p(ρV)−qdvg)q−1q(∫Me−fV2|ϕV|pdvg)p−qpq. |
Then,
∫MV3|ϕV|p(ρV)−qdμ≤(pq−C−1)q(∫MV3|ϕV|pdμ)p−qp(∫MV3|∇(ϕV)|pdμ)qp. |
This completes the proof of (ⅰ).
(ⅱ) Similar to the derivation of (3.4) above, the inequality (3.5) can be obtained without any difficulty. This completes the proof of Theorem 3.3.
Remark 3.4. In the special case that V=1, Theorem 3.3 reduces to the result of Du and Mao in [24, Theorem 2.1].
Then, we prove the following Hardy-Poincaré type inequalities related to the divergence type operator Lf,V.
Theorem 3.5. Let (Mn,g,dμ) be an n-dimensional (n≥2) complete non-compact smooth metric measure space. Suppose that V is a positive twice differentiable function and f is a real-valued smooth function on M. Let ρ be a nonnegative function on M such that |∇(ρV)|=1 and Lf,Vρ≥(C(ρV)−1+G)V2 in the sense of distribution, where C>0 is a constant and G is a continuous function. Then, for any p,q,α∈R, φ∈C∞0(M∖(ρV)−1{0}), p∈(1,∞), Aα=(C+α+1)p with C+α+1>0, we have the following inequality:
∫MV3(ρV)α+p|⟨∇(ρV),∇(φV)⟩|pdμ≥Apα∫MV3(ρV)α|φV|pdμ+Ap−1α∫MV3G(ρV)α+1|φV|pdμ. | (3.8) |
Proof. It follows from |∇(ρV)|=1 and Lf,Vρ≥(C(ρV)−1+G)V2 that
efdiv(e−fV2(ρV)∇(ρV))=⟨∇(ρV),V2∇(ρV)⟩+ef(ρV)div(e−fV2∇(ρV))=V2|∇(ρV)|2+(ρV)Lf,Vρ≥V2(1+C+G(ρV)). | (3.9) |
Multiplying (3.9) by (ρV)α|φV|p and integrating both sides over M gives
(1+C)∫Me−fV2(ρV)α|φV|pdvg+∫Me−fGV2(ρV)α+1|φV|pdvg≤∫M(ρV)α|φV|pdiv(e−fV2(ρV)∇(ρV))dvg=−∫M⟨(ρV)∇(ρV),e−fV2∇((ρV)α|φV|p)⟩dvg=−α∫Me−fV2(ρV)α−1|φV|p(ρV)⟨∇(ρV),∇(ρV)⟩dvg−p∫Me−fV2|φV|p−2(φV)(ρV)α(ρV)⟨∇(φV),∇(ρV)⟩dvg=−α∫Me−fV2(ρV)α|φV|pdvg−p∫Me−fV2(ρV)α+1|φV|p−2(φV)⟨∇(ρV),∇(φV)⟩dvg, | (3.10) |
which implies
(1+C+α)∫Me−fV2(ρV)α|φV|pdvg+∫Me−fV2G(ρV)α+1|φV|pdvg≤−p∫Me−fV2(ρV)α+1|φV|p−2(φV)⟨∇(ρV),∇(φV)⟩dvg. | (3.11) |
It follows from the Hölder inequality that
−p∫Me−fV2(ρV)α+1|φV|p−2(φV)⟨∇(ρV),∇(φV)⟩dvg≤|p∫Me−fV2(ρV)α+1|φV|p−2(φV)⟨∇(ρV),∇(φV)⟩dvg|≤p∫Me−fV2(ρV)α+1|φV|p−1|⟨∇(ρV),∇(φV)⟩|dvg=p∫M(e−fV2(ρV)α|φV|p)p−1p(e−fV2(ρV)α+p|⟨∇(ρV),∇(φV)⟩|p)1pdvg≤p(∫Me−fV2(ρV)α|φV|pdvg)p−1p(∫Me−fV2(ρV)α+p|⟨∇(ρV),∇(φV)⟩|pdvg)1p. |
Thus,
−p∫Me−fV2(ρV)α+1|φV|p−2(φV)⟨∇(ρV),∇(φV)⟩dvg≤p(∫Me−fV2(ρV)α|φV|pdvg)p−1p(∫Me−fV2(ρV)α+p|⟨∇(ρV),∇(φV)⟩|pdvg)1p. | (3.12) |
Then, we use the Young inequality in (3.12), which is described as follows: Denoting
Φ=:(∫Me−fV2(ρV)α+p|⟨∇(ρV),∇(φV)⟩|pdvg)1pandΨ=:(∫Me−fV2(ρV)α|φV|pdvg)p−1p, |
then for any ϵ>0,
ΦΨ=ϵΦΨϵ≤1p(ϵΦ)p+1q(Ψϵ)qwithq=pp−1relating top, |
we can get
(∫Me−fV2(ρV)α|φV|pdvg)p−1p(∫Me−fV2(ρV)α+p|⟨∇(ρV),∇(φV)⟩|pdvg)1p≤(p−1)pϵpp−1∫Me−fV2(ρV)α|φV|pdvg+1pϵp∫Me−fV2(ρV)α+p|⟨∇(ρV),∇(φV)⟩|pdvg, |
namely,
−p∫Me−fV2(ρV)α+1|φV|p−2(φV)⟨∇(ρV),∇(φV)⟩dvg≤(p−1)ϵpp−1∫Me−fV2(ρV)α|φV|pdvg+ϵp∫Me−fV2(ρV)α+p|⟨∇(ρV),∇(φV)⟩|pdvg. | (3.13) |
Hence, putting (3.11) and (3.13) together, we can get
∫Me−fV2(ρV)α+p|⟨∇(ρV),∇(φV)⟩|pdvg≥ϵ−p(1+C+α−p−1ϵpp−1)∫Me−fV2(ρV)α|φV|pdvg+ϵ−p∫Me−fV2G(ρV)α+1|φV|pdvg. | (3.14) |
We take Y(ϵ) to be the function Y(ϵ)=ϵ−p(1+C+α−p−1ϵpp−1). By direct computation, we can conclude that Y(ϵ) reaches its maximum value when ϵ=(p1+C+α)p−1p(for details, see [25, Theorem 3]).
Finally, the required inequality can be determined by substituting ϵ=(p1+C+α)p−1p into (3.14) as follows:
∫MV3(ρV)α+p|⟨∇(ρV),∇(φV)⟩|pdμ≥Apα∫MV3(ρV)α|φV|pdμ+Ap−1α∫MV3G(ρV)α+1|φV|pdμ. |
This completes the proof of Theorem 3.5.
Using Theorem 3.5, we can obtain the following applications:
Theorem 3.6. Let (Mn,g,dμ) be an n-dimensional complete noncompact smooth metric measure space. Suppose that V is a positive twice differentiable function and f is a real-valued smooth function on M. Let ρ be a nonnegative function on M such that |∇(ρV)|=1 in the sense of distributions. Then for any p,q∈R,1<p<+∞,0≤q≤p, and any compactly supported smooth function ϕV∈C∞0(M∖(ρV)−1{0}), we have
(i) When Lf,Vρ≤C(ρV)−1V2 in the sense of distributions, where C>0 is a constant and C+α<−1, the following inequality holds:
∫MV3(ρV)α+p|⟨∇(ρV),∇(ϕV)⟩|pdμ≥(|C+α+1|p)p∫MV3(ρV)α|ϕV|pdμ. | (3.15) |
(ii) When Lf,Vρ≥C(ρV)−1V2 in the sense of distributions, where C>0 is a constant and C+α>−1, the following inequality holds:
∫MV3(ρV)α+p|⟨∇(ρV),∇(ϕV)⟩|pdμ≥(C+α+1p)p∫MV3(ρV)α|ϕV|pdμ. | (3.16) |
Proof. (ⅰ) It follows from |∇(ρV)|=1 and Lf,Vρ≥C(ρV)−1V2 that
div(e−fV2(ρV)∇(ρV))=e−fV2|∇(ρV)|2+(ρV)div(e−fV2∇(ρV))≥(1+C)e−fV2. |
Thus,
div(e−fV2(ρV)∇(ρV))≥(1+C)e−fV2. | (3.17) |
Multiplying (3.17) by (ρV)α|ϕV|pand integrating over M yields
(1+C)∫Me−fV2(ρV)α|ϕV|pdvg≤∫M(ρV)α|ϕV|pdiv(e−fV2(ρV)∇(ρV))dvg=−∫M⟨(ρV)∇(ρV),e−fV2∇((ρV)α|ϕV|p)⟩dvg=−α∫Me−fV2(ρV)α|ϕV|pdvg−p∫Me−fV2|ϕV|p−2(ϕV)(ρV)α+1⟨∇(ρV),∇(ϕV)⟩dvg. | (3.18) |
Since C+α+1<0, by using the Hölder inequality, we can infer from (3.18) that
|C+α+1|∫Me−fV2(ρV)α|ϕV|pdvg≤−p∫Me−fV2|ϕV|p−2(ϕV)(ρV)α+1⟨∇(ρV),∇(ϕV)⟩dvg≤|p∫Me−fV2|ϕV|p−2(ϕV)(ρV)α+1⟨∇(ρV),∇(ϕV)⟩dvg|≤p∫Me−fV2|ϕV|p−1(ρV)α+1|⟨∇(ρV),∇(ϕV)⟩|dvg≤p∫M(e−fV2(ρV)α|ϕV|p)p−1p(e−fV2(ρV)α+p|⟨∇(ρV),∇(ϕV)⟩|p)1pdvg≤p(∫Me−fV2(ρV)α|ϕV|pdvg)p−1p(∫Me−fV2(ρV)α+p|⟨∇(ρV),∇(ϕV)⟩|pdvg)1p. | (3.19) |
It follows from the Young inequality that
|C+α+1|∫Me−fV2(ρV)α|ϕV|pdvg≤p(∫Me−fV2(ρV)α|ϕV|pdvg)p−1p(∫Me−fV2(ρV)α+p|⟨∇(ρV),∇(ϕV)⟩|pdvg)1p≤(p−1)[ϵ−1(∫Me−fV2(ρV)α|ϕV|pdvg)p−1p]pp−1+[ϵ(∫Me−fV2(ρV)α+p|⟨∇(ρV),∇(ϕV)⟩|pdvg)1p]p=(p−1)ϵ−pp−1∫Me−fV2(ρV)α|ϕV|pdvg+ϵp∫Me−fV2(ρV)α+p|⟨∇(ρV),∇(ϕV)⟩|pdvg. | (3.20) |
Thus, for any ϵ>0, from (3.20), we have
∫Me−fV2(ρV)α+p|⟨∇(ρV),∇(ϕV)⟩|pdvg≥ϵ−p(|C+α+1|−(p−1)ϵ−pp−1)∫Me−fV2(ρV)α|ϕV|pdvg. |
Taking
ϵ=(p|C+α+1|)p−1p, |
in the above inequality, we can get
∫MV3(ρV)α+p|⟨∇(ρV),∇(ϕV)⟩|pdμ≥(|C+α+1|p)p∫MV3(ρV)α|ϕV|pdμ. |
(ⅱ) Similar to the proof of (3.15) above, the inequality (3.16) can be obtained without any difficulty. This completes the proof of Theorem 3.6.
Remark 3.7. In the special case that V=1, Theorem 3.6 reduces to the result of Du and Mao in [24, Theorem 4.1].
Then, we will prove the weighted Lp-Hardy type inequality.
Theorem 3.8. Let (Mn,g,dμ) be an n-dimensional (n≥2) complete non-compact smooth metric measure space. Suppose that V is a positive twice differentiable function and f is a real-valued smooth function on M. Let ρ be a nonnegative function on M such that (ρV)Lf,Vρ≥CV2|∇(ρV)|2 in the sense of distribution, where C>0 is a constant. Then, the following inequality holds for any ϕV∈C∞0(M∖ρ−1{0}):
∫MV3(ρV)α(sinhβ(ρV))|∇(ρV)|2−p|∇(ϕV)|pdμ≥(C+1+α+β−pp)p∫MV3(ρV)α−p(sinhβ(ρV))|ϕV|p|∇(ρV)|2dμ, | (3.21) |
where p,α,β∈R,1≤p<∞, and C+1+α+β−p>0.
Proof. By direct computation, we have
∫Me−fV2(ρV)α−psinhβ(ρV)|ϕV|p|∇(ρV)|2dvg≤1C∫Me−f(ρV)α−p+1sinhβ(ρV)|ϕV|p(Lf,Vρ)dvg=1C∫M(ρV)α−p+1|ϕV|psinhβ(ρV)div(e−fV2∇(ρV))dvg=−1C∫M⟨∇(sinhβ(ρV)|ϕV|p(ρV)α−p+1),e−fV2∇(ρV)⟩dvg=−∫M(α−p+1+β(ρV)coth(ρV))Ce−fV2(ρV)α−psinhβ(ρV)|ϕV|p|∇(ρV)|2dvg−pC∫Me−fV2(ρV)α−p+1sinhβ(ρV)|ϕV|p−2(ϕV)⟨∇(ρV),∇(ϕV)⟩dvg, | (3.22) |
where we have used the assumption (ρV)Lf,Vρ≥CV2|∇(ρV)|2 in the first line of the inequality.
Observing that (ρV)coth(ρV)≥1, we can rewrite the inequality (3.22) as
∫Me−fV2(ρV)α−psinhβ(ρV)|ϕV|p|∇(ρV)|2dvg≤−(α−p+1+β)C∫Me−fV2(ρV)α−psinhβ(ρV)|ϕV|p|∇(ρV)|2dvg−pC∫Me−fV2(ρV)α−p+1sinhβ(ρV)|ϕV|p−2(ϕV)⟨∇(ρV),∇(ϕV)⟩dvg. |
Then,
(C+α−p+1+β)∫Me−fV2(ρV)α−psinhβ(ρV)|ϕV|p|∇(ρV)|2dvg≤p∫Me−fV2(ρV)α−p+1sinhβ(ρV)|ϕV|p−1⟨∇(ρV),∇(ϕV)⟩dvg. | (3.23) |
It follows from the Hölder inequality that
(C+α−p+1+β)∫Me−fV2(ρV)α−psinhβ(ρV)|ϕV|p|∇(ρV)|2dvg≤p∫Me−fV2(ρV)α−p+1sinhβ(ρV)|ϕV|p−1⟨∇(ρV),∇(ϕV)⟩dvg≤|p∫Me−fV2(ρV)α−p+1sinhβ(ρV)|ϕV|p−1⟨∇(ρV),∇(ϕV)⟩dvg|≤p∫M(e−fV2(ρV)α−psinhβ(ρV)|ϕV|p|∇(ρV)|2)p−1p(e−fV2(ρV)αsinhβ(ρV)|∇(ρV)|2−p|∇(ϕV)|p)1pdvg≤p(∫Me−fV2(ρV)α−psinhβ(ρV)|ϕV|p|∇(ρV)|2dvg)p−1p(∫Me−fV2(ρV)αsinhβ(ρV)|∇(ρV)|2−p|∇(ϕV)|pdvg)1p. | (3.24) |
namely,
∫MV3(ρV)αsinhβ(ρV)|∇(ρV)|2−p|∇(ϕV)|pdμ≥(C+1+α+β−pp)p∫MV3(ρV)α−psinhβ(ρV)|ϕV|p|∇(ρV)|2dμ. |
This completes the proof of Theorem 3.8.
In this section, by applying Theorem 3.8 of Section 3, we can give the following Rellich type inequality.
Theorem 4.1. Let (Mn,g,dμ) be an n-dimensional (n≥2) complete non-compact smooth metric measure space. Suppose that V is a positive twice differentiable function and f is a real-valued smooth function on M. Let ρ be a nonnegative function on M such that |∇(ρV)|=1 and (ρV)Lf,Vρ≥CV2 in the sense of distributions, where C>0 is a constant and max{3−C,(7−C)/3}≤α≤2. Then, the following inequality holds for ϕV∈C∞0(M∖(ρV)−1{0}):
∫MV3(ρV)α|Lf,Vϕ|2dμ≥(C+1−α2)2∫MV3(ρV)α−2|∇(ϕV)|2dμ, | (4.1) |
where dμ=1Ve−fdvg and dvg is the Riemannian volume element related to g.
Proof. As a consequence of integration by parts (1.3), we compute and estimate the righthand side:
∫Me−fV2(ρV)α−2|∇(ϕV)|2dvg=∫Me−fV2(ρV)α−2⟨∇(ϕV),∇(ϕV)⟩dvg=∫M⟨∇(ϕV),e−fV2(ρV)α−2∇(ϕV)⟩dvg=−∫M(ϕV)div(e−fV2(ρV)α−2∇(ϕV))dvg=−∫M(ϕV)(e−fV2⟨∇(ρV)α−2,∇(ϕV)⟩+(ρV)α−2div(e−fV2∇(ϕV)))dvg=−∫Me−fV2(ϕV)⟨∇(ρV)α−2,∇(ϕV)⟩dvg−∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg=−(α−2)∫Me−fV2(ϕV)(ρV)α−3⟨∇(ρV),∇(ϕV)⟩dvg−∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg=−(α−2)2∫Me−fV2(ρV)α−3⟨∇(ρV),∇(ϕV)2⟩dvg−∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg=−(α−2)2∫M(ρV)α−3⟨e−fV2∇(ρV),∇(ϕV)2⟩dvg−∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg=(α−2)2∫M(ϕV)2div(e−fV2(ρV)α−3∇(ρV))dvg−∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg. | (4.2) |
It follows from |∇(ρV)|=1 and (ρV)Lf,Vρ≥CV2 that
∫Me−fV2(ρV)α−2|∇(ϕV)|2dvg=(α−2)2∫M(ϕV)2div(e−fV2(ρV)α−3∇(ρV))dvg−∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg=−(2−α)2∫M((α−3)e−fV2(ρV)α−4(ϕV)2|∇(ρV)|2+(ρV)α−3(ϕV)2div(e−fV2∇(ρV)))dvg−∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg=−(2−α)2∫M((α−3)e−fV2(ρV)α−4(ϕV)2+e−f(ρV)α−3(ϕV)2(Lf,Vρ))dvg−∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg≤−(2−α)(C+α−3)2∫Me−fV2(ρV)α−4(ϕV)2dvg−∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg. |
Thus,
∫Me−fV2(ρV)α−2|∇(ϕV)|2dvg≤−(2−α)(C+α−3)2∫Me−fV2(ρV)α−4(ϕV)2dvg−∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg. | (4.3) |
Applying the Hölder inequality to the above, we have
−∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg≤|∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg|≤(∫Me−fV2(ρV)α|Lf,Vϕ|2dvg)12(∫Me−fV−2(ϕV)2(ρV)α−4dvg)12. |
Thus,
−∫Me−f(ρV)α−2(ϕV)(Lf,Vϕ)dvg≤(∫Me−fV2(ρV)α|Lf,Vϕ|2dvg)12(∫Me−fV−2(ϕV)2(ρV)α−4dvg)12. | (4.4) |
Then, substituting (4.4) into (4.3),
(∫Me−fV2(ρV)α|Lf,Vϕ|2dvg)12(∫Me−fV−2(ϕV)2(ρV)α−4dvg)12≥∫Me−fV2(ρV)α−2|∇(ϕV)|2dvg+(2−α)(C+α−3)2∫Me−fV2(ρV)α−4(ϕV)2dvg. |
We denote
A=∫Me−fV2(ρV)α−2|∇(ϕV)|2dvg,B=∫Me−fV2(ρV)α−4(ϕV)2dvg,D=∫Me−fV2(ρV)α|Lf,Vϕ|2dvg, |
and then can get
√D≥√A(√AB+(2−α)(C+α−3)2√BA). |
Denoting further S=AB,
√D≥√A(√S+(2−α)(C+α−3)2√1S):=√Af(S). |
Now, we compute a lower bound for the function f(S)=√S+((2−α)(C+α−3)/2)√1/S. To begin, we use Theorem 3.8 (in the special case thatβ=0,p=2,|∇(ρV)|=1) and then can get
A=∫Me−fV2(ρV)α−2|∇(ϕV)|2dvg≥(C+α−32)2∫Me−fV2(ρV)α−4(ϕV)2dvg=(C+α−32)2B. |
We have the following condition on S:
S=AB≥(C+α−32)2. |
Next, the function f(S) is increasing for S∈(0,+∞) and it attains its minimum at S=(2−α)(C+α−3)/2. However, (2−α)(C+α−3)/2≤((C+α−3)/2)2 when α≥(7−C)/3, so f(S) attains its minimum at S=((C+α−3)/2)2, and this minimum is equal to ((C+1−α)/2)2. Finally, we obtain the following inequality:
D=∫MV3(ρV)α|Lf,Vϕ|2dμ≥(C+1−α2)2A=(C+1−α2)2∫MV3(ρV)α−2|∇(ϕV)|2dμ. |
This completes the proof of Theorem 4.1.
In this paper, we have established some Hardy type identities and inequalities for the divergence type operator Lf,V on smooth metric measure spaces. First, we have established some L2 Hardy type identities. As their corollary, we have obtained a L2 Hardy type inequality. Second, we have established some Lp Hardy type inequalities. As their corollary, we have obtained a Lp Rellich type inequality. From the proof of the above results, we see that our method does not work for us to obtain the sharp constants. Hence, we shall further pursue sharp Hardy type inequalities for the divergence type operator Lf,V on smooth metric measure spaces in the subsequent papers.
Pengyan Wang: Methodology, writing-original draft; Jiahao Wang: validation, writing-review and editing. All authors have read and agreed to the published version of the manuscript.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The research of authors is supported by NSFC (No.12101530), the Science and Technology Project of Henan Province (No.232102310321), Nanhu Scholars Program for Young Scholars of XYNU (No.2023), and Xinyang Normal University Graduate Student Innovation Fund Project (2024KYJJ060).
The authors declare no conflict of interest.
[1] |
N. T. J. Bailey, The mathematical theory of infectious diseases and its applications, J. R. Stat. Soc. C-Appl., 26 (1977), 85–87. https://doi.org/10.2307/2346882 doi: 10.2307/2346882
![]() |
[2] |
C. Holing, The functional response of predators to prey density and its role in mimicry and population regulation, Mem. Entomol. Soc. Can., 45 (1965), 1–65. https://doi.org/10.4039/entm9745fv doi: 10.4039/entm9745fv
![]() |
[3] |
R. Arditi, L. R. Ginzburg, Coupling in predator-prey dynamics: Ratio dependence, J. Theor. Biol., 139 (1989), 311–326. https://doi.org/10.1016/S0022-5193(89)80211-5 doi: 10.1016/S0022-5193(89)80211-5
![]() |
[4] | J. R. Beddington, Mutual interference between parasites or predators and its effect on searhing efficiency, J. Anim. Ecodogy, 1975,331–340. |
[5] |
D. L. DeAngelis, R. A. Goldstein, R. V. A. O'Neill, A model for tropic interaction, Ecology, 56 (1975), 881–892. https://doi.org/10.2307/1936298 doi: 10.2307/1936298
![]() |
[6] |
P. Crowley, E. Martin, Functional responses and interference within and between year classes of a dragonfly population, J. North. Am. Benth. Soc., 8 (1989), 211–221. https://doi.org/10.2307/1467324 doi: 10.2307/1467324
![]() |
[7] |
M. Kirane, S. Kouachi, Global solutions to a system of strongly coupled reaction-diffusion equations, Nonlinear Anal., 26 (1996), 1387–1396. https://doi.org/10.1016/0362-546x(94)00337-h doi: 10.1016/0362-546x(94)00337-h
![]() |
[8] |
Y. Choi, R. Lui, Y. Yamada, Existence of global solutions for the Shigesada-Kawasaki-Teramoto model with strongly coupled cross-diffusion, Discrete Con. Dyn., 10 (2004), 719–730. https://doi.org/10.3934/dcds.2004.10.719 doi: 10.3934/dcds.2004.10.719
![]() |
[9] |
P. Tuoc, On global existence of solutions to a cross-diffusion system, J. Math. Anal. Appl., 343 (2008), 826–834. https://doi.org/10.1016/j.jmaa.2008.01.089 doi: 10.1016/j.jmaa.2008.01.089
![]() |
[10] |
S. Fu, S. Cui, Global existence and stability of solution of a reaction-diffusion model for cancer invasion, Nonlinear Anal. RWA, 10 (2009), 1362–1369. https://doi.org/10.1016/j.nonrwa.2008.01.011 doi: 10.1016/j.nonrwa.2008.01.011
![]() |
[11] | H. Xu, S. Xu, Existence and stability of global solutions for a cross-diffusion predator-prey model with sex-structure, Nonlinear Anal. RWA, 13 (2012) 999–1009. https://doi.org/10.1016/j.nonrwa.2010.08.029 |
[12] | S. Fu, L. Zhang, P. Hu, Global behavior of solutions in a Lotka-Volterra predator-prey model with prey-stage structure, Nonlinear Anal. RWA, 14 (2013) 2027–2045. https://doi.org/10.1016/j.nonrwa.2013.02.007 |
[13] |
X. Li, S. Fu, Global stability of the virus dynamics model with intracellular delay and Crowley-Martin functional response, Math. Meth. Appl. Sci., 37 (2014), 1405–1411. https://doi.org/10.1002/mma.2895 doi: 10.1002/mma.2895
![]() |
[14] |
X. Li, S. Fu, Global stability of a virus dynamics model with intracellular delay and CTL immune response, Math. Meth. Appl. Sci., 38 (2015), 420–430. https://doi.org/10.1002/mma.3078 doi: 10.1002/mma.3078
![]() |
[15] | H. Amann, Dynamic theory of quasilinear parabolic equations: Reaction-diffusion, Differ. Integral Equ., 3 (1990), 13–75. https://doi.org/10.1007/BF02571246 |
[16] | H. Amann, Dynamic theory of quasilinear parabolic equations: Global existence, Math. Z., 202 (1989), 219–250. https://doi.org/10.1007/BF01215256 |
[17] | O. A. Ladyzenskaja, V. A. Solonnikov, N. N. Ural'ceva, Linear and quasilinear equations of parabolic type, Translations of Mathematical Monographs, AMS, 1968. |
[18] | A. Haraux, M. Kirane, Estimations C1 pour des problemes paraboliques semi-linéaires, Annales de la Faculte des Sciences de Toulouse, University Paul Sabatier, 1983. https://doi.org/10.5802/afst.598 |
[19] | M. Protter, H. Weinerger, Maximum principles in differential equations, (Section Editions), Springer-Verlag, New York, 1984. https://doi.org/10.1007/978-1-4612-5282-5 |
[20] |
P. Pang, M. Wang, Strategy and stationary pattern in a three-species predator-prey model, J. Differ. Eq., 200 (2004), 245–273. https://doi.org/10.1016/j.jde.2004.01.004 doi: 10.1016/j.jde.2004.01.004
![]() |
[21] |
D. Jiang, Z. C. Wang, L. Zhang, A reaction-diffusion-advection SIS epidemic model in a spatially-temporally heterogeneous environment, Discrete Con. Dyn. Sys.-B, 23 (2018), 4557–4578. https://doi.org/10.3934/dcdsb.2018176 doi: 10.3934/dcdsb.2018176
![]() |
[22] | H. Li, R. Peng, T. Xiang, Dynamics and asymptotic profiles of endemic equilibrium for two frequency-dependent SIS epidemic models with cross-diffusion, Eur. J. Appl. Math., 2017, 1–31. https://doi.org/10.1017/S0956792518000463 |
[23] |
F. Y. Yang, W. T. Li, Dynamics of a nonlocal dispersal SIS epidemic model, Commun. Pure Appl. Anal., 16 (2017), 781–798. https://doi.org/10.3934/cpaa.2017037 doi: 10.3934/cpaa.2017037
![]() |
[24] |
D. Hnaien, K. Ferdaous, L. Rafika, Asymptotic behavior of global solutions of an anomalous diffusion system, J. Math. Anal. Appl., 421 (2015), 1519–1530. https://doi.org/10.1016/j.jmaa.2014.07.083 doi: 10.1016/j.jmaa.2014.07.083
![]() |
[25] |
Y. Wang, J. Wang, L. Zhang, Cross diffusion-induced pattern in an SI model, Appl. Math. Comput., 217 (2010), 1965–1970. https://doi.org/10.1016/j.amc.2010.06.052 doi: 10.1016/j.amc.2010.06.052
![]() |
[26] |
Y. Fan, Pattern formation of an epidemic model with cross diffusion, Appl. Math. Comput., 228 (2014), 311–319. https://doi.org/10.1016/j.amc.2013.11.090 doi: 10.1016/j.amc.2013.11.090
![]() |
[27] | C. S. Ruoja, C. Surulescu, A. Zhigun, On a model for epidemic spread with interpopulation contact and repellent taxis, arXiv: 1902.02171, 2019. |
[28] | D. Le, L. V. Nguyen, T. T. Nguyen, Shigesada-Kawasaki-Teramoto model on higher dimen-sional domains, Electron J. Differ. Eq., 72 (2003). https://doi.org/10.1023/A:1022197004856 |
[29] |
L. T. Hoang, T. V. Nguyen, P. V. Tuoc, Gradient estimates and global existence of smooth solutions to a cross-diffusion system, SIAM J. Math. Anal., 47 (2015), 2122–2177. https://doi.org/10.1137/140981447 doi: 10.1137/140981447
![]() |
[30] |
P. V. Tuoc, Global existence of solutions to Shigesada-Kawasaki-Teramoto cross-diffusion systems on domains of arbitrary dimensions, Proc. Amer. Math. Soc., 135 (2007), 3933–3941. https://doi.org/10.1090/S0002-9939-07-08978-2 doi: 10.1090/S0002-9939-07-08978-2
![]() |
[31] |
D. Le, T. T. Nguyen, Global existence for a class of triangular parabolic systems on domains of arbitrary dimension, Proc. Amer. Math. Soc., 133 (2005), 1985–1992. https://doi.org/10.1090/S0002-9939-05-07867-6 doi: 10.1090/S0002-9939-05-07867-6
![]() |
[32] |
Y. Lou, W. M. Ni, Y. Wu, On the global existence of a cross-diffusion system, Discrete Cont. Dyn-A, 4 (1998), 193–203. https://doi.org/10.3934/dcds.1998.4.193 doi: 10.3934/dcds.1998.4.193
![]() |
[33] | Y. Yamada, Global solutions for the Shigesada-Kawasaki-Teramoto model with cross- diffusion, Recent progress on reaction-diffusion systems and viscosity solutions, World Scientific River Edge, NJ, 2009,282–299. https://doi.org/10.1142/9789812834744_0013 |
[34] |
Y. Tao, M. Winkler, Boundedness and stabilization in a population model with cross-diffusion for one species, Proc. London Math. Soc., 119 (2019), 1598–1632. https://doi.org/10.1112/plms.12276 doi: 10.1112/plms.12276
![]() |
[35] |
Y. Liu, Y. Tao, Dynamics in a parabolic elliptic two species population competition model with cross-diffusion for one species, J. Math. Anal. Appl., 456 (2017), 1–15. https://doi.org/10.1016/j.jmaa.2017.05.058 doi: 10.1016/j.jmaa.2017.05.058
![]() |
[36] |
A. A. Khan, R. Amin, S. Ullah, Numerical simulation of a Caputo fractional epidemic model for the novel coronavirus with the impact of environmental transmission, J. Alex. Eng., 61 (2022), 729–746. https://doi.org/10.1016/j.aej.2021.10.008 doi: 10.1016/j.aej.2021.10.008
![]() |
[37] | A. Alshabanat, M. Jleli, S. Kumar, B. Samet, Generalization of Caputo-Fabrizio fractional derivative and applications to electrical circuits, Front. Phy., 8 (2020). https://doi.org/10.3389/fphy.2020.00064 |
[38] | P. Veeresha, D. G. Prakasha, S. Kumar, A fractional model for propagation of classical optical solitons by using nonsingular derivative, Math. Meth. Appl. Sci., 2020, 1–15. https://doi.org/10.1002/mma.6335 |
[39] |
H. Mohammadi, S. Kuma, S. Rezapour, S. Etemad, A theoretical study of the Caputo-Fabrizio fractional modeling for hearing loss due to Mumps virus with optimal control, Chaos Soliton. Fract., 144 (2021), 110668. https://doi.org/10.1016/j.chaos.2021.110668 doi: 10.1016/j.chaos.2021.110668
![]() |
[40] |
E. F. D. Goufo, S. Kumar, S. B. Mugisha, Similarities in a fifth-order evolution equation with and with no singular kernel, Chaos Soliton. Fract., 130 (2020), 109467. https://doi.org/10.1016/j.chaos.2019.109467 doi: 10.1016/j.chaos.2019.109467
![]() |
[41] |
B. Ghanbari, S. Kumar, R. Kumar, A study of behaviour for immune and tumor cells in immunogenetic tumour model with non-singular fractional derivative, Chaos Soliton. Fract., 133 (2020), 109619. https://doi.org/10.1016/j.chaos.2020.109619 doi: 10.1016/j.chaos.2020.109619
![]() |
[42] |
S. Kumar, A. Kumar, B. Samet, J. F. Gomez-Aguilar, M. S. Osman, A chaos study of tumor and effector cells in fractional tumor-immune model for cancer treatment, Chaos Soliton. Fract., 141 (2020), 110321. https://doi.org/10.1016/j.chaos.2020.110321 doi: 10.1016/j.chaos.2020.110321
![]() |
[43] |
S. Kumar, R. Kumar, R. P. Agarwal, B. Samet, A study of fractional Lotka-Volterra population model using Haar wavelet and Adam-Bashforth-Moulton methods, Math. Meth. Appl. Sci., 43 (2020), 5564–5578. https://doi.org/10.1002/mma.6297 doi: 10.1002/mma.6297
![]() |