Processing math: 100%
Research article

Global behavior of solutions to an SI epidemic model with nonlinear diffusion in heterogeneous environment

  • Received: 28 November 2021 Revised: 07 January 2022 Accepted: 18 January 2022 Published: 26 January 2022
  • MSC : 35J60, 35B32, 92D25

  • 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

    Related Papers:

    [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(n22)2Rn|ϕ|2|x|2dx,

    for n3 and ϕC0(Rn{0}). The constant (n22)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(n4)4)2Rn|ϕ|2|x|4dx,

    where n5 and ϕC0(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)24Mρα2ϕ2dvg,

    where αR, C+α1>0, ϕC0(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)pMραpϕpdvg,

    where 1p<, C+1+αp>0, and ϕC0(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)pMρα|ϕ|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 n2, 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,qR with q>1+C and any compactly supported smooth function ϕC0(Mρ1{0}), the following inequality holds:

    (qC1)M|ϕ|pρqdvg|p|M|ϕ|p1ρq1|ϕ|dvg+MH|ϕ|pρq1dvg.

    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μ=1Vefdvg, 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=ΔfVVg1V2V+Ricmf,

    where

    Δf:=efdiv(ef)=Δf, (1.1)

    is the f-Laplacian (see [18,19,20,21]),

    Ricmf=Ric+2f1mnff,

    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[efV2(φ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μ=1Vefdvg. That is, for any smooth functions u,vC0(M), we have

    MuLf,Vvdμ=MV3uVvVdμ=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:

    mm1Mφ2efdvgMef(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

    mm1MVφ2efdvgMefV(^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 (n2) 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 WC1(M). Then, for any compactly supported smooth functions uVC10(M) and φVC20(M), the following identity holds:

    MWV3|(uV)|2dμ=MWV3|(uV)(uV)(φV)1(φV)|2dμMefV(uV)2(φV)1div(efV2W(φV))dμ, (2.1)

    where dμ=1Vefdvg 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μ=MefV2W|(uV)|2dvg2MefV2W(uV),(uV)(φV)1(φV)dvg+MefV2W(uV)2(φV)2|(φV)|2dvg=MefV2W|(uV)|2dvgMefV2W(φV)1(uV)2,(φV)dvg+MefV2W(uV)2(φV)2|(φV)|2dvg=MefV2W|(uV)|2dvg+M(uV)2div(efV2W(φV)1(φV))dvg+MefV2W(uV)2(φV)2|(φV)|2dvg=MefV2W|(uV)|2dvg+M(φV)1(uV)2div(efV2W(φV))dvgM(φV)2(uV)2(φV),efV2W(φV)dvg+MefV2W(uV)2(φV)2|(φV)|2dvg=MWV3|(uV)|2dμ+MefV(uV)2(φV)1div(efV2W(φV))dμ.

    Then, we can get

    MWV3|(uV)|2dμ=MWV3|(uV)(uV)(φV)1(φV)|2dμMefV(uV)2(φV)1div(efV2W(φ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 (n2) 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 uVC10(M) and φVC20(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μ=1Vefdvg 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 (n2) 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 uVC10(M) and φVC20(M), the following inequality holds:

    MV3|(uV)|2dμMV(uV)2(φV)1(Lf,Vφ)dμ, (2.3)

    where dμ=1Vefdvg and dvg is the Riemannian volume element related to g.

    As the special case that φV=(ρV)1C2 in the Corollary 2.3, we have the following.

    Corollary 2.4. Let (Mn,g,dμ) be an n-dimensional (n2) 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 uVC10(M(ρV)1{0}), the following inequality holds:

    MV3|(uV)|2dμ(C1)24MV3(uV)2(ρV)2dμ,

    where dμ=1Vefdvg and dvg is the Riemannian volume element related to g.

    Proof. Taking φV=(ρV)1C2, we have

    (φV)1(Lf,Vφ)=ef(ρV)C12div(efV2(ρV)1C2)=(ρV)C12efdiv((1C2)efV2(ρV)C12(ρV))=(1C2)(ρV)C12ef(ρV)C12,efV2(ρV)(1C2)(ρV)C12(ρV)C12efdiv(efV2(ρV))=1C24V2(ρV)2|(ρV)|2+C12(ρV)1Lf,Vρ.

    Then, we have

    (φV)1(Lf,Vφ)=1C24V2(ρV)2|(ρV)|2+C12(ρV)1Lf,Vρ1C24V2(ρV)2+C2C2V2(ρV)2=C22C+14V2(ρV)2=(C1)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φ)(C1)24V2(ρV)2. (2.4)

    Substituting (2.4) into (2.3), we can get

    MV3|(uV)|2dμ(C1)24MV3(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 (n2) 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,qR,q>1+C, and any compactly supported smooth functions ϕVC0(M(ρV)1{0}), the following inequality holds:

    (qC1)MV3|ϕV|p(ρV)qdμ|p|MV3|ϕV|p1(ρV)1q|(ϕV)|dμ+MHV3|ϕV|p(ρV)1qdμ. (3.1)

    Proof. For a vector field X on M, we denote by div X the divergence of X. Note that

    div(efV2|ϕV|p(ρV)1q(ρV))=pefV2|ϕV|p1|ϕV|,(ρV)1q(ρV)+|ϕV|pdiv(efV2(ρV)1q(ρV))=pefV2|ϕV|p1|ϕV|,(ρV)1q(ρV)+|ϕV|p(ρV)1qdiv(efV2(ρV))+(1q)efV2|ϕV|p(ρV)q(ρV),(ρV)=pefV2|ϕV|p1|ϕV|,(ρV)1q(ρV)+ef|ϕV|p(ρV)1q(Lf,Vρ)+(1q)efV2|ϕV|p(ρV)q|(ρV)|2.

    By direct computation, we have

    div(efV2|ϕV|p(ρV)1q(ρV))=pefV2|ϕV|p1|ϕV|,(ρV)1q(ρV)+ef|ϕV|p(ρV)1q(Lf,Vρ)+(1q)efV2|ϕV|p(ρV)q|(ρV)|2pefV2|ϕV|p1|ϕV|,(ρV)1q(ρV)+efV2|ϕV|p(ρV)1q(H+C(ρV)1)+(1q)efV2|ϕV|p(ρV)q=pefV2|ϕV|p1|ϕV|,(ρV)1q(ρV)+(1+Cq)efV2|ϕV|p(ρV)q+efV2H|ϕV|p(ρV)1q, (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

    (qC1)MefV2|ϕV|p(ρV)qdvgpMefV2|ϕV|p1|ϕV|,(ρV)1q(ρV)dvg+MefV2H|ϕV|p(ρV)1qdvg|pMefV2|ϕV|p1|ϕV|,(ρV)1q(ρV)dvg|+MefV2H|ϕV|p(ρV)1qdvg|p|MefV2|ϕV|p1(ρV)1q|(ϕV)||(ρV)|dvg+MefV2H|ϕV|p(ρV)1qdvg=|p|MefV2|ϕV|p1(ρV)1q|(ϕV)|dvg+MefV2H|ϕV|p(ρV)1qdvg.

    Then, we have

    (qC1)MV3|ϕV|p(ρV)qdμ|p|MV3|ϕV|p1(ρV)1q|(ϕV)|dμ+MHV3|ϕV|p(ρV)1qdμ. (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 (n2) 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,qR,1<p<+,0qp, and any compactly supported smooth function ϕVC0(M(ρV)1{0}), we have

    (i) When Lf,VρC(ρV)1V2 in the sense of distributions, where C<q1 is a constant, the following inequality holds:

    MV3|ϕV|p(ρV)qdμ(pqC1)q(MV3|ϕV|pdμ)pqp(MV3|(ϕV)|pdμ)qp. (3.4)

    (ii) When Lf,VρC(ρV)1V2 in the sense of distributions, where C>q1 is a constant, the following inequality holds:

    MV3|ϕV|p(ρV)qdμ(pC+1q)q(MV3|ϕV|pdμ)pqp(MV3|(ϕV)|pdμ)qp. (3.5)

    Proof. (ⅰ) In the special case that H=0 in (3.1), we can get

    (qC1)MV3|ϕV|p(ρV)qdμ|p|MV3|ϕV|p1(ρV)1q|(ϕV)|dμ. (3.6)

    It follows from the Hölder inequality that

    MefV2|ϕV|p1(ρV)1q|(ϕV)|dvg=M(efV2)1p(efV2)p1p|ϕV|p1(ρV)1q|(ϕV)|dvg(MefV2|(ϕV)|pdvg)1p(MefV2|ϕV|p(ρV)p(1q)p1dvg)p1p(MefV2|(ϕV)|pdvg)1p(MefV2|ϕV|p(ρV)qdvg)q1q(MefV2|ϕV|pdvg)(1p(q1)q(p1))p1p. (3.7)

    Substituting (3.7) into (3.6), we can get

    (qC1)MefV2|ϕV|p(ρV)qdvg|p|(MefV2|(ϕV)|pdvg)1p(MefV2|ϕV|p(ρV)qdvg)q1q(MefV2|ϕV|pdvg)(1p(q1)q(p1))p1p|p|(MefV2|(ϕV)|pdvg)1p(MefV2|ϕV|p(ρV)qdvg)q1q(MefV2|ϕV|pdvg)pqpq.

    Then,

    MV3|ϕV|p(ρV)qdμ(pqC1)q(MV3|ϕV|pdμ)pqp(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 (n2) 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, φC0(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μ+Ap1αMV3G(ρV)α+1|φV|pdμ. (3.8)

    Proof. It follows from |(ρV)|=1 and Lf,Vρ(C(ρV)1+G)V2 that

    efdiv(efV2(ρV)(ρV))=(ρV),V2(ρV)+ef(ρV)div(efV2(ρ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)MefV2(ρV)α|φV|pdvg+MefGV2(ρV)α+1|φV|pdvgM(ρV)α|φV|pdiv(efV2(ρV)(ρV))dvg=M(ρV)(ρV),efV2((ρV)α|φV|p)dvg=αMefV2(ρV)α1|φV|p(ρV)(ρV),(ρV)dvgpMefV2|φV|p2(φV)(ρV)α(ρV)(φV),(ρV)dvg=αMefV2(ρV)α|φV|pdvgpMefV2(ρV)α+1|φV|p2(φV)(ρV),(φV)dvg, (3.10)

    which implies

    (1+C+α)MefV2(ρV)α|φV|pdvg+MefV2G(ρV)α+1|φV|pdvgpMefV2(ρV)α+1|φV|p2(φV)(ρV),(φV)dvg. (3.11)

    It follows from the Hölder inequality that

    pMefV2(ρV)α+1|φV|p2(φV)(ρV),(φV)dvg|pMefV2(ρV)α+1|φV|p2(φV)(ρV),(φV)dvg|pMefV2(ρV)α+1|φV|p1|(ρV),(φV)|dvg=pM(efV2(ρV)α|φV|p)p1p(efV2(ρV)α+p|(ρV),(φV)|p)1pdvgp(MefV2(ρV)α|φV|pdvg)p1p(MefV2(ρV)α+p|(ρV),(φV)|pdvg)1p.

    Thus,

    pMefV2(ρV)α+1|φV|p2(φV)(ρV),(φV)dvgp(MefV2(ρV)α|φV|pdvg)p1p(MefV2(ρV)α+p|(ρV),(φV)|pdvg)1p. (3.12)

    Then, we use the Young inequality in (3.12), which is described as follows: Denoting

    Φ=:(MefV2(ρV)α+p|(ρV),(φV)|pdvg)1pandΨ=:(MefV2(ρV)α|φV|pdvg)p1p,

    then for any ϵ>0,

    ΦΨ=ϵΦΨϵ1p(ϵΦ)p+1q(Ψϵ)qwithq=pp1relating top,

    we can get

    (MefV2(ρV)α|φV|pdvg)p1p(MefV2(ρV)α+p|(ρV),(φV)|pdvg)1p(p1)pϵpp1MefV2(ρV)α|φV|pdvg+1pϵpMefV2(ρV)α+p|(ρV),(φV)|pdvg,

    namely,

    pMefV2(ρV)α+1|φV|p2(φV)(ρV),(φV)dvg(p1)ϵpp1MefV2(ρV)α|φV|pdvg+ϵpMefV2(ρV)α+p|(ρV),(φV)|pdvg. (3.13)

    Hence, putting (3.11) and (3.13) together, we can get

    MefV2(ρV)α+p|(ρV),(φV)|pdvgϵp(1+C+αp1ϵpp1)MefV2(ρV)α|φV|pdvg+ϵpMefV2G(ρV)α+1|φV|pdvg. (3.14)

    We take Y(ϵ) to be the function Y(ϵ)=ϵp(1+C+αp1ϵpp1). By direct computation, we can conclude that Y(ϵ) reaches its maximum value when ϵ=(p1+C+α)p1p(for details, see [25, Theorem 3]).

    Finally, the required inequality can be determined by substituting ϵ=(p1+C+α)p1p into (3.14) as follows:

    MV3(ρV)α+p|(ρV),(φV)|pdμApαMV3(ρV)α|φV|pdμ+Ap1α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,qR,1<p<+,0qp, and any compactly supported smooth function ϕVC0(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)pMV3(ρ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)pMV3(ρV)α|ϕV|pdμ. (3.16)

    Proof. (ⅰ) It follows from |(ρV)|=1 and Lf,VρC(ρV)1V2 that

    div(efV2(ρV)(ρV))=efV2|(ρV)|2+(ρV)div(efV2(ρV))(1+C)efV2.

    Thus,

    div(efV2(ρV)(ρV))(1+C)efV2. (3.17)

    Multiplying (3.17) by (ρV)α|ϕV|pand integrating over M yields

    (1+C)MefV2(ρV)α|ϕV|pdvgM(ρV)α|ϕV|pdiv(efV2(ρV)(ρV))dvg=M(ρV)(ρV),efV2((ρV)α|ϕV|p)dvg=αMefV2(ρV)α|ϕV|pdvgpMefV2|ϕV|p2(ϕ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|MefV2(ρV)α|ϕV|pdvgpMefV2|ϕV|p2(ϕV)(ρV)α+1(ρV),(ϕV)dvg|pMefV2|ϕV|p2(ϕV)(ρV)α+1(ρV),(ϕV)dvg|pMefV2|ϕV|p1(ρV)α+1|(ρV),(ϕV)|dvgpM(efV2(ρV)α|ϕV|p)p1p(efV2(ρV)α+p|(ρV),(ϕV)|p)1pdvgp(MefV2(ρV)α|ϕV|pdvg)p1p(MefV2(ρV)α+p|(ρV),(ϕV)|pdvg)1p. (3.19)

    It follows from the Young inequality that

    |C+α+1|MefV2(ρV)α|ϕV|pdvgp(MefV2(ρV)α|ϕV|pdvg)p1p(MefV2(ρV)α+p|(ρV),(ϕV)|pdvg)1p(p1)[ϵ1(MefV2(ρV)α|ϕV|pdvg)p1p]pp1+[ϵ(MefV2(ρV)α+p|(ρV),(ϕV)|pdvg)1p]p=(p1)ϵpp1MefV2(ρV)α|ϕV|pdvg+ϵpMefV2(ρV)α+p|(ρV),(ϕV)|pdvg. (3.20)

    Thus, for any ϵ>0, from (3.20), we have

    MefV2(ρV)α+p|(ρV),(ϕV)|pdvgϵp(|C+α+1|(p1)ϵpp1)MefV2(ρV)α|ϕV|pdvg.

    Taking

    ϵ=(p|C+α+1|)p1p,

    in the above inequality, we can get

    MV3(ρV)α+p|(ρV),(ϕV)|pdμ(|C+α+1|p)pMV3(ρ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 (n2) 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 ϕVC0(Mρ1{0}):

    MV3(ρV)α(sinhβ(ρV))|(ρV)|2p|(ϕV)|pdμ(C+1+α+βpp)pMV3(ρV)αp(sinhβ(ρV))|ϕV|p|(ρV)|2dμ, (3.21)

    where p,α,βR,1p<, and C+1+α+βp>0.

    Proof. By direct computation, we have

    MefV2(ρV)αpsinhβ(ρV)|ϕV|p|(ρV)|2dvg1CMef(ρV)αp+1sinhβ(ρV)|ϕV|p(Lf,Vρ)dvg=1CM(ρV)αp+1|ϕV|psinhβ(ρV)div(efV2(ρV))dvg=1CM(sinhβ(ρV)|ϕV|p(ρV)αp+1),efV2(ρV)dvg=M(αp+1+β(ρV)coth(ρV))CefV2(ρV)αpsinhβ(ρV)|ϕV|p|(ρV)|2dvgpCMefV2(ρV)αp+1sinhβ(ρV)|ϕV|p2(ϕ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

    MefV2(ρV)αpsinhβ(ρV)|ϕV|p|(ρV)|2dvg(αp+1+β)CMefV2(ρV)αpsinhβ(ρV)|ϕV|p|(ρV)|2dvgpCMefV2(ρV)αp+1sinhβ(ρV)|ϕV|p2(ϕV)(ρV),(ϕV)dvg.

    Then,

    (C+αp+1+β)MefV2(ρV)αpsinhβ(ρV)|ϕV|p|(ρV)|2dvgpMefV2(ρV)αp+1sinhβ(ρV)|ϕV|p1(ρV),(ϕV)dvg. (3.23)

    It follows from the Hölder inequality that

    (C+αp+1+β)MefV2(ρV)αpsinhβ(ρV)|ϕV|p|(ρV)|2dvgpMefV2(ρV)αp+1sinhβ(ρV)|ϕV|p1(ρV),(ϕV)dvg|pMefV2(ρV)αp+1sinhβ(ρV)|ϕV|p1(ρV),(ϕV)dvg|pM(efV2(ρV)αpsinhβ(ρV)|ϕV|p|(ρV)|2)p1p(efV2(ρV)αsinhβ(ρV)|(ρV)|2p|(ϕV)|p)1pdvgp(MefV2(ρV)αpsinhβ(ρV)|ϕV|p|(ρV)|2dvg)p1p(MefV2(ρV)αsinhβ(ρV)|(ρV)|2p|(ϕV)|pdvg)1p. (3.24)

    namely,

    MV3(ρV)αsinhβ(ρV)|(ρV)|2p|(ϕV)|pdμ(C+1+α+βpp)pMV3(ρ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 (n2) 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{3C,(7C)/3}α2. Then, the following inequality holds for ϕVC0(M(ρV)1{0}):

    MV3(ρV)α|Lf,Vϕ|2dμ(C+1α2)2MV3(ρV)α2|(ϕV)|2dμ, (4.1)

    where dμ=1Vefdvg 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:

    MefV2(ρV)α2|(ϕV)|2dvg=MefV2(ρV)α2(ϕV),(ϕV)dvg=M(ϕV),efV2(ρV)α2(ϕV)dvg=M(ϕV)div(efV2(ρV)α2(ϕV))dvg=M(ϕV)(efV2(ρV)α2,(ϕV)+(ρV)α2div(efV2(ϕV)))dvg=MefV2(ϕV)(ρV)α2,(ϕV)dvgMef(ρV)α2(ϕV)(Lf,Vϕ)dvg=(α2)MefV2(ϕV)(ρV)α3(ρV),(ϕV)dvgMef(ρV)α2(ϕV)(Lf,Vϕ)dvg=(α2)2MefV2(ρV)α3(ρV),(ϕV)2dvgMef(ρV)α2(ϕV)(Lf,Vϕ)dvg=(α2)2M(ρV)α3efV2(ρV),(ϕV)2dvgMef(ρV)α2(ϕV)(Lf,Vϕ)dvg=(α2)2M(ϕV)2div(efV2(ρV)α3(ρV))dvgMef(ρV)α2(ϕV)(Lf,Vϕ)dvg. (4.2)

    It follows from |(ρV)|=1 and (ρV)Lf,VρCV2 that

    MefV2(ρV)α2|(ϕV)|2dvg=(α2)2M(ϕV)2div(efV2(ρV)α3(ρV))dvgMef(ρV)α2(ϕV)(Lf,Vϕ)dvg=(2α)2M((α3)efV2(ρV)α4(ϕV)2|(ρV)|2+(ρV)α3(ϕV)2div(efV2(ρV)))dvgMef(ρV)α2(ϕV)(Lf,Vϕ)dvg=(2α)2M((α3)efV2(ρV)α4(ϕV)2+ef(ρV)α3(ϕV)2(Lf,Vρ))dvgMef(ρV)α2(ϕV)(Lf,Vϕ)dvg(2α)(C+α3)2MefV2(ρV)α4(ϕV)2dvgMef(ρV)α2(ϕV)(Lf,Vϕ)dvg.

    Thus,

    MefV2(ρV)α2|(ϕV)|2dvg(2α)(C+α3)2MefV2(ρV)α4(ϕV)2dvgMef(ρV)α2(ϕV)(Lf,Vϕ)dvg. (4.3)

    Applying the Hölder inequality to the above, we have

    Mef(ρV)α2(ϕV)(Lf,Vϕ)dvg|Mef(ρV)α2(ϕV)(Lf,Vϕ)dvg|(MefV2(ρV)α|Lf,Vϕ|2dvg)12(MefV2(ϕV)2(ρV)α4dvg)12.

    Thus,

    Mef(ρV)α2(ϕV)(Lf,Vϕ)dvg(MefV2(ρV)α|Lf,Vϕ|2dvg)12(MefV2(ϕV)2(ρV)α4dvg)12. (4.4)

    Then, substituting (4.4) into (4.3),

    (MefV2(ρV)α|Lf,Vϕ|2dvg)12(MefV2(ϕV)2(ρV)α4dvg)12MefV2(ρV)α2|(ϕV)|2dvg+(2α)(C+α3)2MefV2(ρV)α4(ϕV)2dvg.

    We denote

    A=MefV2(ρV)α2|(ϕV)|2dvg,B=MefV2(ρV)α4(ϕV)2dvg,D=MefV2(ρV)α|Lf,Vϕ|2dvg,

    and then can get

    DA(AB+(2α)(C+α3)2BA).

    Denoting further S=AB,

    DA(S+(2α)(C+α3)21S):=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=MefV2(ρV)α2|(ϕV)|2dvg(C+α32)2MefV2(ρ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 α(7C)/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)2MV3(ρ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
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2403) PDF downloads(109) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog