Processing math: 92%
Research article

An analysis of the isoparametric bilinear finite volume element method by applying the Simpson rule to quadrilateral meshes

  • In this work, we construct and study a special isoparametric bilinear finite volume element scheme for solving anisotropic diffusion problems on general convex quadrilateral meshes. The new scheme is obtained by employing the Simpson rule to approximate the line integrals in the classical isoparametric bilinear finite volume element method. By using the cell analysis approach, we suggest a sufficient condition to ensure the coercivity of the new scheme. The sufficient condition has an analytic expression, which only involves the anisotropic diffusion tensor and the geometry of quadrilateral mesh. This yields that for any diffusion tensor and quadrilateral mesh, we can directly judge whether this sufficient condition is satisfied. Specifically, this condition covers the traditional h1+γ-parallelogram and some trapezoidal meshes with any full anisotropic diffusion tensor. An optimal H1 error estimate of the proposed scheme is also obtained for a quasi-parallelogram mesh. The theoretical results are verified by some numerical experiments.

    Citation: Shengying Mu, Yanhui Zhou. An analysis of the isoparametric bilinear finite volume element method by applying the Simpson rule to quadrilateral meshes[J]. AIMS Mathematics, 2023, 8(10): 22507-22537. doi: 10.3934/math.20231147

    Related Papers:

    [1] Muhammad Farman, Ali Akgül, Kottakkaran Sooppy Nisar, Dilshad Ahmad, Aqeel Ahmad, Sarfaraz Kamangar, C Ahamed Saleel . Epidemiological analysis of fractional order COVID-19 model with Mittag-Leffler kernel. AIMS Mathematics, 2022, 7(1): 756-783. doi: 10.3934/math.2022046
    [2] Minghung Lin, Yiyou Hou, Maryam A. Al-Towailb, Hassan Saberi-Nik . The global attractive sets and synchronization of a fractional-order complex dynamical system. AIMS Mathematics, 2023, 8(2): 3523-3541. doi: 10.3934/math.2023179
    [3] Rahat Zarin, Amir Khan, Aurangzeb, Ali Akgül, Esra Karatas Akgül, Usa Wannasingha Humphries . Fractional modeling of COVID-19 pandemic model with real data from Pakistan under the ABC operator. AIMS Mathematics, 2022, 7(9): 15939-15964. doi: 10.3934/math.2022872
    [4] Pratap Anbalagan, Evren Hincal, Raja Ramachandran, Dumitru Baleanu, Jinde Cao, Michal Niezabitowski . A Razumikhin approach to stability and synchronization criteria for fractional order time delayed gene regulatory networks. AIMS Mathematics, 2021, 6(5): 4526-4555. doi: 10.3934/math.2021268
    [5] Muhammad Farman, Ali Akgül, Sameh Askar, Thongchai Botmart, Aqeel Ahmad, Hijaz Ahmad . Modeling and analysis of fractional order Zika model. AIMS Mathematics, 2022, 7(3): 3912-3938. doi: 10.3934/math.2022216
    [6] Jingfeng Wang, Chuanzhi Bai . Global Mittag-Leffler stability of Caputo fractional-order fuzzy inertial neural networks with delay. AIMS Mathematics, 2023, 8(10): 22538-22552. doi: 10.3934/math.20231148
    [7] Yuehong Zhang, Zhiying Li, Wangdong Jiang, Wei Liu . The stability of anti-periodic solutions for fractional-order inertial BAM neural networks with time-delays. AIMS Mathematics, 2023, 8(3): 6176-6190. doi: 10.3934/math.2023312
    [8] Weerawat Sudsutad, Jutarat Kongson, Chatthai Thaiprayoon, Nantapat Jarasthitikulchai, Marisa Kaewsuwan . A generalized Gronwall inequality via ψ-Hilfer proportional fractional operators and its applications to nonlocal Cauchy-type system. AIMS Mathematics, 2024, 9(9): 24443-24479. doi: 10.3934/math.20241191
    [9] Rana Safdar Ali, Saba Batool, Shahid Mubeen, Asad Ali, Gauhar Rahman, Muhammad Samraiz, Kottakkaran Sooppy Nisar, Roshan Noor Mohamed . On generalized fractional integral operator associated with generalized Bessel-Maitland function. AIMS Mathematics, 2022, 7(2): 3027-3046. doi: 10.3934/math.2022167
    [10] Yonghong Liu, Ghulam Farid, Dina Abuzaid, Hafsa Yasmeen . On boundedness of fractional integral operators via several kinds of convex functions. AIMS Mathematics, 2022, 7(10): 19167-19179. doi: 10.3934/math.20221052
  • In this work, we construct and study a special isoparametric bilinear finite volume element scheme for solving anisotropic diffusion problems on general convex quadrilateral meshes. The new scheme is obtained by employing the Simpson rule to approximate the line integrals in the classical isoparametric bilinear finite volume element method. By using the cell analysis approach, we suggest a sufficient condition to ensure the coercivity of the new scheme. The sufficient condition has an analytic expression, which only involves the anisotropic diffusion tensor and the geometry of quadrilateral mesh. This yields that for any diffusion tensor and quadrilateral mesh, we can directly judge whether this sufficient condition is satisfied. Specifically, this condition covers the traditional h1+γ-parallelogram and some trapezoidal meshes with any full anisotropic diffusion tensor. An optimal H1 error estimate of the proposed scheme is also obtained for a quasi-parallelogram mesh. The theoretical results are verified by some numerical experiments.



    As we all know, the study of variable exponent function space inspired by nonlinear elasticity theory and nonstandard growth differential equations is one of the key contents of harmonic analysis in the past three decades, attracting extensive attention from many scholars. In [19], the theory of function spaces with variable exponent was progressed since some elementary properties were established by Kováčik and Rákosník, and they studied many basic properties of variable exponent Lebesgue spaces and Sobolev spaces on Rn. Later, the Lebesgue spaces with variable exponent Lp()(Rn) were extensively investigated; see [7,8,22]. In [14], Izuki first introduced the Herz spaces with variable exponent ˙Kα,qp()(Rn), which are generalizations of the Herz spaces ˙Kα,qp(Rn), and considered the boundedness of commutators of fractional integrals on Herz spaces with variable exponent. In [13], Izuki introduced the Herz-Morrey spaces with variable exponent M˙Kα,λq,p()(Rn), which are generalizations of the Herz-Morrey spaces M˙Kα,λq,p(Rn), and studied the boundedness of vector valued sublinear operators on Herz-Morrey spaces with variable exponent M˙Kα,λq,p()(Rn). On the other hand, in the study of boundary value problems for the Laplace equation on Lipschitz domains, the classical theory of Muckenhoupt weights is a powerful tool in harmonic analysis; see [21]. Generalized Muckenhoupt weights with variable exponent have been intensively studied; see [4,5].

    In [11], Hardy defined the classical Hardy operators as:

    P(f)(x):=1xx0f(t)dt,x>0. (1.1)

    In [6], Christ and Grafakos defined the ndimensional Hardy operators as:

    H(f)(x):=1|x|n|t|<|x|f(t)dt,xRn{0}, (1.2)

    and established the boundedness of P(f)(x) in Lp(Rn), getting the best constants.

    In [9], under the condition of 0β<n and |x|=ni=1x2i, Fu et al. defined the ndimensional fractional Hardy operators and its adjoint operators as:

    Hβf(x):=1|x|nβ|t|<|x|f(t)dt,Hβf(x):=|t||x|f(t)|t|nβdt,xRn{0}, (1.3)

    and established the boundedness of their commutators in Lebesgue spaces and homogeneous Herz spaces.

    Let L1loc(Rn) be the collection of all locally integrable functions on Rn. Given a function bL1loc(Rn) and mN, Wang et al. [23] defined the mth order commutators of ndimensional fractional Hardy operators and adjoint operators as:

    Hmβ,bf(x):=1|x|nβ|t|<|x|(b(x)b(t))mf(t)dt (1.4)

    and

    Hmβ,bf(x):=|t||x|(b(x)b(t))mf(t)|t|nβdt,xRn{0}. (1.5)

    Obviously, when m=0, H0β,b=Hβ, H0β,b=Hβ, and when m=1, H1β,b=Hβ,b, H1β,b=Hβ,b. More important results with regard to these commutators, see [20,26,27].

    Due to the need of future calculation in this paper, let 0<β<n, and the fractional integral operator Iβ is defined as:

    Iβ(f)(x):=Rnf(y)|xy|nβdy,xRn. (1.6)

    Let 0β<n and fL1loc(Rn), and the fractional maximal operator Mβ is defined as:

    Mβf(x):=supxB1|B|1βnB|f(y)|dy,xRn, (1.7)

    where the supremum is taken over all balls BRn containing x. When β=0, we simply write M instead of M0, which is exactly the Hardy-Littlewood maximal function.

    Let fL1loc(Rn) and BMO(Rn) consist of all fL1loc(Rn) with BMO(Rn)<. b is a bounded mean oscillation function if bBMO<, and the bBMO is defined as follow:

    bBMO:=supBB|b(x)bB|dx, (1.8)

    where the supremum is taken all over the balls BRn and bB:=|B|1Bb(y)dy. For a comprehensive review of the bounded mean oscillation (BMO) space, please see the book [10].

    Recently, Muhammad Asim et al. established the estimates of fractional Hardy operators on weighted variable exponent Morrey-Herz spaces in [1]. Amjad Hussain et al. established the boundedness of the commutators of the Fractional Hardy operators on weighted variable Herz-Morrey spaces in [12]. Motivated by the mentioned work, in this paper, we will give the boundedness of the mth order commutators of ndimensional fractional Hardy operators Hmβ,b and its adjoint operators Hmβ,b on weighted variable exponent Morrey-Herz space M˙Kα,λq,p()(ω).

    The paper is organized as follows. In Section 2, we provide some preliminary knowledge. The main results and their proofs are given in Section 3. In Section 4, we provide the conclusion of this paper. Throughout this paper, we use the following symbols and notations:

    (1) For a constant R>0 and a point xRn, we write B(x,R):={yRn:|xy|<R}.

    (2) For any measurable set ERn, |E| denotes the Lebesgue measure, and χE means the characteristic function.

    (3) Given kZ, we write Bk:=¯B(0,2k)={xRn:|x|2k}.

    (4) We define a family {Ak}k= by Ak:=BkBk1={xRn:2k1<|x|2k}. Moreover χk denotes the characteristic function of Ak, namely, χk:=χAk.

    (5) For any index 1<p(x)<, p(x) denotes its conjugate index, namely, 1p(x)+1p(x)=1.

    (6) If there exists a positive constant C independent of the main parameters such that ACB, then we write AB. Additionally, AB means that both AB and BA hold.

    Definition2.1. ([7]) Let p():Rn[1,) be a measurable function.

    (ⅰ) The Lebesgue space with variable exponent Lp()(Rn) is defined by

    Lp()(Rn):={fismeasurablefunction:Rn(|f(x)|λ)p(x)dx<forsomeconstantλ>0}.

    (ⅱ) The spaces with variable exponent Lp()loc(Rn) are defined by

    Lp()loc(Rn):={fismeasurablefunction:fLp()(K)forallcompactsubsetsKRn}.

    The Lebesgue space Lp()(Rn) is a Banach space with the norm defined by

    fLp()(Rn):=inf{λ>0:Rn(|f(x)|λ)p(x)dx1}.

    Definition2.2. ([7]) (ⅰ) The set P(Rn) consists of all measurable functions p():Rn[1,) satisfying

    1<pp(x)p+<,

    where

    p:=essinf{p(x):xRn}>1,p+:=esssup{p(x):xRn}<.

    (ⅱ) The set B(Rn) consists of all measurable function p()P(Rn) satisfying that the Hardy-Littlewood maximal operator M is bounded on Lp()(Rn).

    Definition2.3. ([7]) Suppose that p() is a real-valued function on Rn. We say that

    (ⅰ) Clogloc(Rn) is the set of all local log-Hölder continuous functions p() satisfying

    |p(x)p(y)|Clog(|xy|),|xy|12,x,yRn. (2.1)

    (ⅱ) Clog0(Rn) is the set of all local log-Hölder continuous functions p() satisfying at origin

    |p(x)p0|Clog(e+1|x|),xRn. (2.2)

    (ⅲ) Clog(Rn) is the set of all local log-Hölder continuous functions satisfying at infinity

    |p(x)p|Clog(e+|x|),xRn. (2.3)

    (ⅳ) Clog(Rn)=Clog(Rn)Clogloc(Rn) denotes the set of all global log-Hölder continuous functions p().

    It was proved in [7] that if p()P(Rn)Clog(Rn), then the Hardy-Littlewood maximal operator M is bounded on Lp()(Rn).

    Definition2.4. ([21]) Given a non-negative, measure function ω, for 1<p<, ωAp if

    [ω]Ap:=supB(1|B|Bω(x)dx)(1|B|Bω(x)1pdx)p1<,

    where the supremum is taken over all balls BRn. Especially, we say ωA1 if

    [ω]A1:=supB1|B|Bω(x)dxessinf{ω(x):xB}<.

    These weights characterize the weighted norm inequalities for the Hardy-Littlewood maximal operator, that is, ωAp, 1<p<, if and only if M:Lp(ω)Lp(ω).

    Definition2.5. ([15]) Suppose that p()P(Rn). A weight ω is in the class Ap() if

    supB:ball|B|1ω1p()χBLp()ω1p()χBLp()<. (2.4)

    Obviously, if p()=p,1<p<, then the above definition reduces to the classical Muckenhoupt Ap class.

    From [15], if p(),q()P(Rn), and p()q(), then A1Ap()Aq().

    Definition2.6. ([15]) Let 0<β<n and p1(),p2()P(Rn) such that 1p2(x)=1p1(x)βn. A weight ω is said to be an A(p1(),p2()) weight if

    χBLp2()(ωp2())χB(Lp1()(ωp1())C|B|1βn. (2.5)

    Definition2.7. ([25]) Let p()P(Rn) and ωAp(). The weighted variable exponent Lebesgue space Lp()(ω) denotes the set of all complex-valued measurable functions f satisfying

    Lp()(ω):={f:fω1p()Lp()(Rn)}.

    This is a Banach space equipped with the norm

    fLp()(ω):=fω1p()Lp()(Rn).

    Definition2.8. ([1]) Let ω be a weight on Rn, 0λ<, 0<q<, p()P(Rn), and α():RnR with α()L(Rn). The weighted variable exponent Morrey-Herz space M˙Kα(),λq,p()(ω) is the set of all measurable functions f given by

    M˙Kα(),λq,p()(ω):={fLp()loc(Rn{0},ω):fM˙Kα(),λq,p()(ω)<},

    where

    fM˙Kα(),λq,p()(ω):=supk0Z2k0λ{k0k=2kα()qfχkqLp()(ω)}1q.

    It is noted that M˙Kα(),0q,p()(ω)=˙Kα()q,p()(ω) is the variable exponent weighted Herz space defined in [2].

    Definition2.9. ([15]) Let M be the set of all complex-valued measurable functions defined on Rn and X be a linear subspace of M.

    (1) The space X is said to be a Banach function space if there exists a function X:M[0,] satisfying the following properties: Let f,g,fjM(j=1,2,). Then

    (a) fX holds if and only if fX<.

    (b) Norm property:

    ⅰ. Positivity: fX0.

    ⅱ. Strict positivity: fX=0 holds if and only if f(x)=0 for almost every xRn.

    ⅲ. Homogeneity: λfX=|λ|fX holds for all λC.

    ⅳ. Triangle inequality: f+gXfX+gX.

    (c) Symmetry: fX=|f|X.

    (d) Lattice property: If 0g(x)f(x) for almost every xRn, then gXfX.

    (e) Fatou property: If 0fj(x)fj+1(x) for all j, and fj(x)f(x) as j for almost every xRn, then limjfjX=fX.

    (f) For every measurable set FRn such that |F|<, χFX is finite. Additionally, there exists a constant CF>0 depending only on F so that F|h(x)|dxCFhX holds for all hX.

    (2) Suppose that X is a Banach function space equipped with a norm X. The associated space X is defined by

    X:={fM:fX<},

    where

    fX:=supg{|Rnf(x)g(x)dx|:gX1}.

    Lemma2.1. ([3]) Let X be a Banach function space, and then we have the following:

    (ⅰ) The associated space X is also a Banach function space.

    (ⅱ) (X) and X are equivalent.

    (ⅲ) If gX and fX, then

    Rn|f(x)g(x)|dxfXgX (2.6)

    is the generalized Hölder inequality.

    Lemma2.2. ([15]) If X is a Banach function space, then we have, for all balls B,

    1|B|1χBXχBX. (2.7)

    Lemma2.3. ([16]) Let X be a Banach function space. Suppose that the Hardy-Littlewood maximal operator M is weakly bounded on X, that is,

    χ{Mf>λ}Xλ1fX

    is true for all fX and all λ>0. Then, we have

    supB:ball1|B|χBXχBX<. (2.8)

    Lemma2.4. ([15]) Given a function W such that 0<W(x)< for almost every xRn, WXloc(Rn) and W1(X)loc(Rn),

    (ⅰ) X(Rn,W) is Banach function space equipped with the norm

    fX(Rn,W):=fWX, (2.9)

    where

    X(Rn,W):={fM:fWX}. (2.10)

    (ⅱ) The associated space X(Rn,W1) of X(Rn,W) is also a Banach function space.

    Lemma2.5. ([15]) Let X be a Banach function space and M be bounded on X. Then, there exists a constant δ(0,1) for all BRn and EB,

    χEXχBX(|E||B|)δ. (2.11)

    The paper [19] shows that Lp()(Rn) is a Banach function space and the associated space Lp()(Rn) with equivalent norm.

    Remark2.6. ([1]) Let p()P(Rn), and by comparing the Lp()(ωp()) and Lp()(ωp()) with the definition of X(Rn,W), we have the following:

    (1) If we take W=ω and X=Lp()(Rn), then we get Lp()(Rn,ω)=Lp()(ωp()).

    (2) If we consider W=ω1 and X=Lp()(Rn), then we get Lp()(Rn,ω1)=Lp()(ωp()).

    By virtue of Lemma 2.4, we get

    (Lp()(Rn,ω))=(Lp()(ωp()))=Lp()(ωp())=Lp()(Rn,ω1).

    Lemma2.7. ([17]) Let p()P(Rn)Clog(Rn) be a log-Hölder continuous function both at infinity and at origin, if ωp2()Ap2() implies ωp2()Ap2(). Thus, the Hardy-Littlewood operator is bounded on Lp2()(ωp2()), and there exist constants δ1,δ2(0,1) such that

    χELp2()(ωp2())χBLp2()(ωp2())=χE(Lp2()(ωp2()))χB(Lp2()(ωp2()))C(|E||B|)δ1, (2.12)

    and

    χE(Lp2()(ωp2()))χB(Lp2()(ωp2()))C(|E||B|)δ2, (2.13)

    for all balls BRn and all measurable sets EB.

    Lemma2.8. ([15]) Let p1()P(Rn)Clog(Rn) and 0<β<np+1. Define p2() by 1p1(x)1p2()=βn. If ωA(p1(),p2()), then Iβ is bounded from Lp1()(ωp1()) to Lp2()(ωp2()).

    Lemma2.9. ([24, Corollary 3.11]) Let bBMO(Rn),mN, and k,jZ with k>j. Then, we have

    C1bmBMO(Rn)supB1χBLp()(ω)(bbB)mχBLp()(ω)CbmBMO(Rn). (2.14)
    (bbBj)mχBkLp()(ω)C(kj)mbmBMO(Rn)χBkLp()(ω). (2.15)

    Proposition3.1. ([12] Let q()P(Rn), 0<p<, and 0λ<. If α()L(Rn)Clog(Rn), then

    fpM˙Kα(),λp,q()(ωq())=supk0Z2k0λpk0j=2jα()pfχjpLq()(ωq())max{supk0Zk0<02k0λp(k0j=2jα(0)pfχjpLq()(ωq())),supk0Zk00(2k0λp(1j=2jα(0)pfχjpLq()(ωq()))+2k0λp(k0j=02jα()pfχjpLq()(ωq())))}.

    Theorem3.1. Let 0<q1q2<, p2()P(Rn)Clog(Rn) and p1() be such that 1p2()=1p1()βn. Also, let ωp2()A1, bBMO(Rn), λ>0 and α()L(Rn)Clog(Rn) be log-Hölder continuous at the origin, with α(0)α()<λ+nδ2β, where δ2(0,1) is the constant appearing in (2.13). Then,

    Hmβ,bfM˙Kα(),λq2,p2()(ωp2())bmBMOfM˙Kα(),λq1,p1()(ωp1()). (3.1)

    Proof. For arbitrary fM˙Kα(),λq1,p1()(ωp1()), let fj=fχj=fχAj for every jZ, and then

    f(x)=j=f(x)χj(x)=j=fj(x). (3.2)

    By the inequality of Cp, it is not difficult to see that

    |Hmβ,bf(x)χk(x)|1|x|nβ|t|<|x||b(x)b(t)|m|f(t)|dtχk(x)1|x|nβB(0,|x|)|b(x)b(t)|m|f(t)|dtχk(x)1|x|nβBk|b(x)b(t)|m|f(t)|dtχk(x)2k(nβ)kj=Aj|b(x)b(t)|m|f(t)|dtχk(x)2k(nβ)kj=Aj|b(x)bAj|m|f(t)|dtχk(x)+2k(nβ)kj=Aj|b(t)bAj|m|f(t)|dtχk(x)=E1+E2. (3.3)

    For E1, by the generalized Hölder inequality, we have

    E1=2k(nβ)kj=Aj|b(x)bAj|m|f(t)|dtχk(x)2k(nβ)kj=|b(x)bAj|mχk(x)fjLp1()(ωp1())χj(Lp1()(ωp1())). (3.4)

    By taking the (Lp2()(ωp2()))-norm on both sides of (3.4) and using (2.15) of Lemma 2.9, we get

    E1(Lp2()(ωp2()))2k(nβ)kj=|b(x)bAj|mχk(Lp2()(ωp2()))fjLp1()(ωp1())χj(Lp1()(ωp1()))2k(nβ)kj=(kj)mbmBMOχk(Lp2()(ωp2()))fjLp1()(ωp1())χj(Lp1()(ωp1())). (3.5)

    For E2, by the generalized Hölder inequality, we have

    E2=2k(nβ)kj=Aj|b(t)bAj|m|f(t)|dtχk(x)2k(nβ)kj=|b(t)bAj|mχj(x)(Lp1()(ωp1()))fjLp1()(ωp1())χk(x). (3.6)

    By taking the (Lp2()(ωp2()))-norm on both sides of (3.6) and using (2.14) of Lemma 2.9, we get

    E2(Lp2()(ωp2()))2k(nβ)kj=|b(x)bAj|mχj(Lp1()(ωp1()))fjLp1()(ωp1())χk(Lp2()(ωp2()))2k(nβ)kj=bmBMOχj(Lp1()(ωp1()))fjLp1()(ωp1())χk(Lp2()(ωp2())). (3.7)

    Hence, from inequalities (3.3), (3.5) and (3.7), we get

    Hmβ,bf(x)χk(Lp2()(ωp2()))2k(nβ)fjLp1()(ωp1()){kj=(kj)mbmBMOχk(Lp2()(ωp2()))χj(Lp1()(ωp1()))+kj=bmBMOχj(Lp1()(ωp1()))χk(Lp2()(ωp2()))}2k(nβ)bmBMOkj=(kj)mfjLp1()(ωp1())χj(Lp1()(ωp1()))χk(Lp2()(ωp2())). (3.8)

    By virtue of Lemma 2.5, we have

    χBkXχkX(|Bk||Ak|)δ=CχBkXCχkX. (3.9)

    Note that χjLp2()(ωp2())χBjLp2()(ωp2()) and χBj(x)2jβIβ(χBj) (see [18, p. 350]). By applying (2.8), (3.9) and Lemma 2.8, we obtain

    χjLp2()(ωp2())χBjLp2()(ωp2())2jβIβ(χBj)Lp2()(ωp2())2jβχBjLp1()(ωp1())2jβχjLp1()(ωp1())2j(nβ)χj1(Lp1()(ωp1())). (3.10)

    By virtue of (2.7) and (2.8), combining (2.13) and (3.10), we have

    2k(βn)χj(Lp1()(ωp1()))χkLp2()(ωp2())=2kβχj(Lp1()(ωp1()))2knχkLp2()(ωp2())2kβχj(Lp1()(ωp1()))χk1(Lp2()(ωp2()))=2kβχj(Lp1()(ωp1()))χj1(Lp2()(ωp2()))χj(Lp2()(ωp2()))χk(Lp2()(ωp2()))2kβ2nδ2(jk)χj(Lp1()(ωp1()))χj1(Lp2()(ωp2()))2kβ2nδ2(jk)2j(nβ)χj1Lp2()(ωp2())χj1(Lp2()(ωp2()))=2kβ2nδ2(jk)2jβ(2jnχjLp2()(ωp2())χj(Lp2()(ωp2())))12(βnδ2)(kj). (3.11)

    Hence by virtue of (3.8) and (3.11), we have

    Hmβ,bf(x)χk(Lp2()(ωp2()))bmBMOkj=(kj)m2(βnδ2)(kj)fjLp1()(ωp1()). (3.12)

    In order to estimate fjLp1()(ωp1()), we consider two cases as below.

    Case 1: For j<0, we get

    fjLp1()(ωp1())=2jα(0)(2jα(0)q1fjq1Lp1()(ωp1()))1q12jα(0)(ji=2iα(0)q1fiq1Lp1()(ωp1()))1q1=2j(λα(0)){2jλ(ji=2iα()q1fiq1Lp1()(ωp1()))1q1}2j(λα(0))fM˙Kα(),λq1,p1()(ωp1()). (3.13)

    Case 2: For j0, we get

    fjLp1()(ωp1())=2jα()(2jα()q1fjq1Lp1()(ωp1()))1q12jα()(ji=2iα()q1fiq1Lp1()(ωp1()))1q1=2j(λα()){2jλ(ji=2iα()q1fiq1Lp1()(ωp1()))1q1}2j(λα())fM˙Kα(),λq1,p1()(ωp1()). (3.14)

    Now, by virtue of the condition q1q2 and the definition of weighted variable exponent Morrey-Herz space along with the use of Proposition 3.1, we get

    Hmβ,bfq1M˙Kα(),λq2,p2()(ωp2())=supk0Z2k0λq1k0k=2kα()q1Hmβ,bfχkq1Lp2()(ωp2())max{supk0Zk0<02k0λq1k0k=2kα(0)q1Hmβ,bfχkq1Lp2()(ωp2()),supk0Zk002k0λq1(1k=2kα(0)q1Hmβ,bfχkq1Lp2()(ωp2())+k0k=02kα()q1Hmβ,bfχkq1Lp2()(ωp2()))}=max{J1,J2+J3}, (3.15)

    where

    J1=supk0Zk0<02k0λq1k0k=2kα(0)q1Hmβ,bfχkq1Lp2()(ωp2()),J2=supk0Zk002k0λq11k=2kα(0)q1Hmβ,bfχkq1Lp2()(ωp2()),J3=supk0Zk002k0λq1k0k=02kα()q1Hmβ,bfχkq1Lp2()(ωp2()).

    First, we estimate J1. Since α(0)α()<nδ2+λβ, combining (3.12) and (3.13), we get

    J1supk0Zk0<02k0λq1k0k=2kα(0)q1(kj=(kj)mbmBMO2(βnδ2)(kj)fLp1()(ωp1()))q1supk0Zk0<02k0λq1k0k=2kα(0)q1(kj=(kj)mbmBMO2(βnδ2)(kj)2j(λα(0))fM˙Kα(),λq1,p1()(ωp1()))q1bmq1BMOfq1M˙Kα(),λq1,p1()(ωp1())supk0Zk0<02k0λq1k0k=2kα(0)q1(kj=(kj)m2(βnδ2)(kj)2j(λα(0)))q1bmq1BMOfq1M˙Kα(),λq1,p1()(ωp1())supk0Zk0<02k0λq1k0k=2kλq1(kj=(kj)m2(jk)(nδ2+λβα(0)))q1bmq1BMOfq1M˙Kα(),λq1,p1()(ωp1()).

    The estimate of J2 is similar to that of J1.

    Lastly, we estimate J3. Since α(0)α()<nδ2+λβ, combining (3.12) and (3.14), we get

    J3supk0Zk002k0λq1k0k=02kα()q1(kj=(kj)mbmBMO2(βnδ2)(kj)fLp1()(ωp1()))q1supk0Zk002k0λq1k0k=02kα()q1(kj=(kj)mbmBMO2(βnδ2)(kj)2j(λα())fM˙Kα(),λq1,p1()(ωp1()))q1bmq1BMOfq1M˙Kα(),λq1,p1()(ωp1())supk0Zk002k0λq1k0k=02kα()q1(kj=(kj)m2(βnδ2)(kj)2j(λα()))q1bmq1BMOfq1M˙Kα(),λq1,p1()(ωp1())supk0Zk002k0λq1k0k=02kλq1(kj=(kj)m2(jk)(nδ2+λβα()))q1bmq1BMOfq1M˙Kα(),λq1,p1()(ωp1()).

    The desired result is obtained by combining the estimates of J1, J2 and J3.

    Theorem3.2. Let 0<q1q2<, p2()P(Rn)Clog(Rn) and p1() be such that 1p2()=1p1()βn. Also, let ωp2()A1, bBMO(Rn), λ>0 and α()L(Rn)Clog(Rn) be log-Hölder continuous at the origin, with λnδ1<α(0)α(), where δ1(0,1) is the constant appearing in (2.12). Then,

    Hmβ,bfM˙Kα(),λq2,p2()(ωp2())bmBMOfM˙Kα(),λq1,p1()(ωp1()). (3.16)

    Proof. From an application of the inequality of Cp, it is not difficult to see that

    |Hmβ,bf(x)χk(x)|RnBk|t|βn|b(x)b(t)|m|f(t)|dtχk(x)j=k+1Aj|t|βn|b(x)b(t)|m|f(t)|dtχk(x)j=k+1Aj|t|βn|b(x)bAj|m|f(t)|dtχk(x)+j=k+1Aj|t|βn|b(t)bAj|m|f(t)|dtχk(x)=F1+F2. (3.17)

    For F1, by the generalized Hölder inequality, we have

    F1j=k+12j(nβ)Aj|b(x)bAj|m|f(t)|dtχk(x)j=k+12j(nβ)|b(x)bAj|mχk(x)fjLp1()(ωp1())χj(Lp1()(ωp1())). (3.18)

    By taking the (Lp2()(ωp2()))-norm on both sides of (3.18) and using (2.15) of Lemma 2.9, we get

    F1(Lp2()(ωp2()))j=k+12j(nβ)|b(x)bAj|mχk(Lp2()(ωp2()))fjLp1()(ωp1())χj(Lp1()(ωp1()))j=k+12j(nβ)(jk)mbmBMOχk(Lp2()(ωp2()))fjLp1()(ωp1())χj(Lp1()(ωp1())). (3.19)

    For F2, by the generalized Hölder inequality, we have

    F2j=k+12j(nβ)Aj|b(t)bAj|m|f(t)|dtχk(x)j=k+12j(nβ)|b(t)bAj|mχj(x)(Lp1()(ωp1()))fjLp1()(ωp1())χk(x). (3.20)

    By taking the (Lp2()(ωp2()))-norm on both sides of (3.20) and using (2.15) of Lemma 2.9, we get

    F2(Lp2()(ωp2()))j=k+12j(nβ)|b(t)bAj|mχj(Lp1()(ωp1()))fjLp1()(ωp1())χk(Lp2()(ωp2()))j=k+12j(nβ)bmBMOχj(Lp1()(ωp1()))fjLp1()(ωp1())χk(Lp2()(ωp2())). (3.21)

    Hence, from inequalities (3.17), (3.19) and (3.21), we get

    Hmβ,bf(x)χk(Lp2()(ωp2()))fjLp1()(ωp1()){j=k+12j(nβ)(jk)mbmBMOχk(Lp2()(ωp2()))χj(Lp1()(ωp1()))+j=k+12j(nβ)bmBMOχj(Lp1()(ωp1()))χk(Lp2()(ωp2()))}bmBMOj=k+12j(nβ)(jk)mfjLp1()(ωp1())χj(Lp1()(ωp1()))χk(Lp2()(ωp2())). (3.22)

    On the other hand, by (2.7) and (2.8), combining (2.12) and (3.10), we have

    2j(nβ)χkLp2()(ωp2())χj(Lp1()(ωp1()))=2jβχkLp2()(ωp2())2jnχj(Lp1()(ωp1()))2jβχkLp2()(ωp2())χj1Lp1()(ωp1())=2jβχj1Lp1()(ωp1())χjLp2()(ωp2())χkLp2()(ωp2())χjLp2()(ωp2())2jβ2nδ1(kj)χj1Lp1()(ωp1())χjLp2()(ωp2())2jβ2nδ1(kj)2j(nβ)χj1Lp1()(ωp1())χj1(Lp1()(ωp1()))=2jβ2nδ1(kj)2jβ(2jnχjLp1()(ωp1())χj(Lp1()(ωp1())))12nδ1(kj). (3.23)

    Hence combining (3.22) and (3.23), we obtain

    Hmβ,bf(x)χk(Lp2()(ωp2()))bmBMOj=k+1(jk)m2nδ1(kj)fjLp1()(ωp1()). (3.24)

    Next, by virtue of the condition q1q2 and the definition of weighted variable exponent Morrey-Herz space along with the use of Proposition 3.1, we get

    Hmβ,bfq1M˙Kα(),λq2,p2()(ωp2())=supk0Z2k0λq1k0k=2kα()q1Hmβ,bfχkq1Lp2()(ωp2())max{supk0Zk0<02k0λq1k0k=2kα(0)q1Hmβ,bfχkq1Lp2()(ωp2()),supk0Zk002k0λq1(1k=2kα(0)q1Hmβ,bfχkq1Lp2()(ωp2())+k0k=02kα()q1Hmβ,bfχkq1Lp2()(ωp2()))}=max{Y1,Y2+Y3}, (3.25)

    where

    Y1=supk0Zk0<02k0λq1k0k=2kα(0)q1Hmβ,bfχkq1Lp2()(ωp2()),Y2=supk0Zk002k0λq11k=2kα(0)q1Hmβ,bfχkq1Lp2()(ωp2()),Y3=supk0Zk002k0λq1k0k=02kα()q1Hmβ,bfχkq1Lp2()(ωp2()).

    First, we estimate Y1. Since λnδ1<α(0)α(), combining (3.24) and (3.13), we get

    Y1supk0Zk0<02k0λq1k0k=2kα(0)q1(j=k+1(jk)mbmBMO2nδ1(kj)fLp1()(ωp1()))q1supk0Zk0<02k0λq1k0k=2kα(0)q1(j=k+1(jk)mbmBMO2nδ1(kj)2j(λα(0))fM˙Kα(),λq1,p1()(ωp1()))q1bmq1BMOfq1M˙Kα(),λq1,p1()(ωp1())supk0Zk0<02k0λq1k0k=2kα(0)q1(j=k+1(jk)m2nδ1(kj)2j(λα(0)))q1bmq1BMOfq1M˙Kα(),λq1,p1()(ωp1())supk0Zk0<02k0λq1k0k=2kλq1(j=k+1(jk)m2(jk)(λnδ1α(0)))q1bmq1BMOfq1M˙Kα(),λq1,p1()(ωp1()).

    The estimate of Y_{2} is similar to that of Y_{1} .

    Lastly, we estimate Y_{3} . Since \lambda-n\delta_{1} < \alpha(0)\leq \alpha(\infty) , combining (3.24) and (3.14), we get

    \begin{align*} Y_{3}&\lesssim \sup\limits_{\substack{k_{0}\in\mathbb{Z}\\k_{0}\geq0}}2^{-k_{0}\lambda q_{1}}\sum\limits_{k = 0}^{k_{0}}2^{k\alpha(\infty) q_{1}}\Big(\sum\limits_{j = k+1}^{\infty}(j-k)^{m}\|b\|_{\mathrm{BMO}}^{m}2^{n\delta_{1}(k-j)}\|f\|_{L^{p_{1}(\cdot)}(\omega^{p_{1}(\cdot)})}\Big)^{q_{1}}\\ &\lesssim \sup\limits_{\substack{k_{0}\in\mathbb{Z}\\k_{0}\geq0}}2^{-k_{0}\lambda q_{1}}\sum\limits_{k = 0}^{k_{0}}2^{k\alpha(\infty) q_{1}}\Big(\sum\limits_{j = k+1}^{\infty}(j-k)^{m}\|b\|_{\mathrm{BMO}}^{m}2^{n\delta_{1}(k-j)}2^{j(\lambda-\alpha(\infty))}\|f\|_{\mathrm{M\dot{K}}_{q_{1},p_{1}(\cdot)}^{\alpha(\cdot), \lambda}(\omega^{p_{1}(\cdot)})}\Big)^{q_{1}}\\ &\lesssim \|b\|_{\mathrm{BMO}}^{mq_{1}}\|f\|^{q_{1}}_{\mathrm{M\dot{K}}_{q_{1},p_{1}(\cdot)}^{\alpha(\cdot), \lambda}(\omega^{p_{1}(\cdot)})}\sup\limits_{\substack{k_{0}\in\mathbb{Z}\\k_{0}\geq0}}2^{-k_{0}\lambda q_{1}}\sum\limits_{k = 0}^{k_{0}}2^{k\alpha(\infty) q_{1}}\Big(\sum\limits_{j = k+1}^{\infty}(j-k)^{m}2^{n\delta_{1}(k-j)}2^{j(\lambda-\alpha(\infty))}\Big)^{q_{1}}\\ &\lesssim \|b\|_{\mathrm{BMO}}^{mq_{1}}\|f\|^{q_{1}}_{\mathrm{M\dot{K}}_{q_{1},p_{1}(\cdot)}^{\alpha(\cdot), \lambda}(\omega^{p_{1}(\cdot)})}\sup\limits_{\substack{k_{0}\in\mathbb{Z}\\k_{0}\geq0}}2^{-k_{0}\lambda q_{1}}\sum\limits_{k = 0}^{k_{0}}2^{k\lambda q_{1}}\Big(\sum\limits_{j = k+1}^{\infty}(j-k)^{m}2^{(j-k)(\lambda-n\delta_{1}-\alpha(\infty))}\Big)^{q_{1}}\\ &\lesssim \|b\|_{\mathrm{BMO}}^{mq_{1}}\|f\|^{q_{1}}_{\mathrm{M\dot{K}}_{q_{1},p_{1}(\cdot)}^{\alpha(\cdot), \lambda}(\omega^{p_{1}(\cdot)})}. \end{align*}

    The desired result is obtained by combining the estimates of Y_{1} , Y_{2} and Y_{3} .

    This paper considers the boundedness for m th order commutators of n- dimensional fractional Hardy operators \mathcal{H}^{m}_{_{\beta, b}} and adjoint operators \mathcal{H}_{\beta, b}^{\ast m} on weighted variable exponent Morrey-Herz spaces \mathrm{M\dot{K}}_{q, p(\cdot)}^{\alpha(\cdot), \lambda}(\omega) . When m = 0 , our main result holds on weighted variable exponent Morrey-Herz space for fractional Hardy operators and generalizes the result of Asim et al. in [1, Theorems 4.2 and 4.3]. When m = 1 , our main result holds on weighted variable exponent Morrey-Herz space for commutators of the fractional Hardy operators and generalizes the result of Hussain et al. in [12, Theorems 18 and 19].

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This paper is supported by the National Natural Science Foundation of China (Grant No. 12161071), Qinghai Minzu University campus level project (Nos. 23GH29, 23GCC10).

    All authors declare that they have no conflicts of interest.



    [1] P. Zhu, R. Li, Generalized difference methods for second order elliptic partial differential equations. Ⅱ. Quadrilateral subdivision, Numer. Math. J. Chin. Univ., 4 (1982), 360–375.
    [2] R. E. Bank, D. J. Rose, Some error estimates for the box method, SIAM J. Numer. Anal., 24 (1987), 777–787. https://doi.org/10.1137/0724050 doi: 10.1137/0724050
    [3] S. Chou, Q. Li, Error estimates in L^2, H^1 and L^{\infty} in covolume methods for elliptic and parabolic problems: A unified approach, Math. Comput., 69 (2000), 103–120. https://doi.org/10.1090/S0025-5718-99-01192-8 doi: 10.1090/S0025-5718-99-01192-8
    [4] Z. Cai, On the finite volume element method, Numer. Math., 58 (1990), 713–735. https://doi.org/10.1007/BF01385651 doi: 10.1007/BF01385651
    [5] I. Mishev, Finite volume element methods for non-definite problems, Numer. Math., 83 (1999), 161–175. https://doi.org/10.1007/s002110050443 doi: 10.1007/s002110050443
    [6] P. Chatzipantelidis, R. Lazarov, Error estimates for a finite volume element method for elliptic PDEs in nonconvex polygonal domains, SIAM J. Numer. Anal., 42 (2005), 1932–1958. https://doi.org/10.1137/S0036142903427639 doi: 10.1137/S0036142903427639
    [7] S. Chou, X. Ye, Unified analysis of finite volume methods for second order elliptic problems, SIAM J. Numer. Anal., 45 (2007), 1639–1653. https://doi.org/10.1137/050643994 doi: 10.1137/050643994
    [8] R. Li, Z. Chen, W. Wu, Generalized difference methods for differential equations: Numerical analysis of finite volume methods, New York: Marcel Dekker, 2000.
    [9] Y. Lin, J. Liu, M. Yang, Finite volume element methods: An overview on recent developments, Int. J. Numer. Anal. Mod. B, 4 (2013), 14–34.
    [10] Z. Zhang, Q. Zou, Some recent advances on vertex centered finite volume element methods for elliptic equations, Sci. China Math., 56 (2013), 2507–2522. https://doi.org/10.1007/s11425-013-4740-8 doi: 10.1007/s11425-013-4740-8
    [11] J. Xu, Q. Zou, Analysis of linear and quadratic simplicial finite volume methods for elliptic equations, Numer. Math., 111 (2009), 469–492. https://doi.org/10.1007/s00211-008-0189-z doi: 10.1007/s00211-008-0189-z
    [12] Z. Chen, R. Li, A. Zhou, A note on the optimal L^{2}-estimate of the finite volume element method, Adv. Comput. Math., 16 (2002), 291–303. https://doi.org/10.1023/A:1014577215948 doi: 10.1023/A:1014577215948
    [13] R. E. Ewing, T. Lin, Y. Lin, On the accuracy of the finite volume element method based on piecewise linear polynomials, SIAM J. Numer. Anal., 39 (2002), 1865–1888. https://doi.org/10.1137/S0036142900368873 doi: 10.1137/S0036142900368873
    [14] C. Erath, D. Praetorius, Adaptive vertex-centered finite volume methods for general second-order linear elliptic partial differential equations, IMA J. Numer. Anal., 39 (2019), 983–1008. https://doi.org/10.1093/imanum/dry006 doi: 10.1093/imanum/dry006
    [15] Y. Li, R. Li, Generalized difference methods on arbitrary quadrilateral networks, J. Comput. Math., 17 (1999), 653–672.
    [16] Z. Zhang, Q. Zou, Vertex-centered finite volume schemes of any order over quadrilateral meshes for elliptic boundary value problems, Numer. Math., 130 (2015), 363–393. https://doi.org/10.1007/s00211-014-0664-7 doi: 10.1007/s00211-014-0664-7
    [17] T. Schmidt, Box schemes on quadrilateral meshes, Computing, 51 (1993), 271–292. https://doi.org/10.1007/BF02238536 doi: 10.1007/BF02238536
    [18] Q. Hong, J. Wu, A Q_1-finite volume element scheme for anisotropic diffusion problems on general convex quadrilateral mesh, J. Comput. Appl. Math., 372 (2020), 112732. https://doi.org/10.1016/j.cam.2020.112732 doi: 10.1016/j.cam.2020.112732
    [19] J. Lv, Y. Li, L^2 error estimates and superconvergence of the finite volume element methods on quadrilateral meshes, Adv. Comput. Math., 37 (2012), 393–416. https://doi.org/10.1007/s10444-011-9215-2 doi: 10.1007/s10444-011-9215-2
    [20] Y. Lin, M. Yang, Q. Zou, L^2 error estimates for a class of any order finite volume schemes over quadrilateral meshes, SIAM J. Numer. Anal., 53 (2015), 2030–2050. https://doi.org/10.1137/140963121 doi: 10.1137/140963121
    [21] C. Nie, S. Shu, H. Yu, W. Xia, Superconvergence and asymptotic expansions for bilinear finite volume element approximation on non-uniform grids, J. Comput. Appl. Math., 321 (2017), 323–335. https://doi.org/10.1016/j.cam.2016.12.024 doi: 10.1016/j.cam.2016.12.024
    [22] W. He, Z. Zhang, Q. Zou, Maximum-norms error estimates for high-order finite volume schemes over quadrilateral meshes, Numer. Math., 138 (2018), 473–500. https://doi.org/10.1007/s00211-017-0912-8 doi: 10.1007/s00211-017-0912-8
    [23] Z. Chen, J. Wu, Y. Xu, Higher-order finite volume methods for elliptic boundary value problems, Adv. Comput. Math., 37 (2012), 191–253. https://doi.org/10.1007/s10444-011-9201-8 doi: 10.1007/s10444-011-9201-8
    [24] X. Wang, Y. Li, L^2 error estimates for high order finite volume methods on triangular meshes, SIAM J. Numer. Anal., 54 (2016), 2729–2749. https://doi.org/10.1137/140988486 doi: 10.1137/140988486
    [25] Y. Zhou, J. Wu, A unified analysis of a class of quadratic finite volume element schemes on triangular meshes, Adv. Comput. Math., 46 (2020), 71. https://doi.org/10.1007/s10444-020-09809-8 doi: 10.1007/s10444-020-09809-8
    [26] X. Wen, Y. Zhou, A coercivity result of quadratic finite volume element schemes over triangular meshes, Adv. Appl. Math. Mech., 15 (2023), 901–931. https://doi.org/10.4208/aamm.OA-2021-0311 doi: 10.4208/aamm.OA-2021-0311
    [27] M. Yang, A second-order finite volume element method on quadrilateral meshes for elliptic equations, ESAIM: M2AN, 40 (2006), 1053–1067. https://doi.org/10.1051/m2an:2007002 doi: 10.1051/m2an:2007002
    [28] J. Lv, Y. Li, Optimal biquadratic finite volume element methods on quadrilateral meshes, SIAM J. Numer. Anal., 50 (2012), 2379–2399. https://doi.org/10.1137/100805881 doi: 10.1137/100805881
    [29] Y. Zhou, Y. Zhang, J. Wu, A polygonal finite volume element method for anisotropic diffusion problems, Comput. Math. Appl., 140 (2023), 225–236. https://doi.org/10.1016/j.camwa.2023.04.025 doi: 10.1016/j.camwa.2023.04.025
    [30] Y. Zhang, X. Wang, Unified construction and L^2 analysis for the finite volume element method over tensorial meshes, Adv. Comput. Math., 49 (2023), 2. https://doi.org/10.1007/s10444-022-10004-0 doi: 10.1007/s10444-022-10004-0
    [31] Y. Zhou, Y. Jiang, Q. Zou, Three dimensional high order finite volume element schemes for elliptic equations, Numer. Methods Partial Differ. Eq., 39 (2023), 1672–1705. https://doi.org/10.1002/num.22950 doi: 10.1002/num.22950
    [32] Y. Zhou, J. Wu, A new high order finite volume element solution on arbitrary triangular and quadrilateral meshes, Appl. Math. Lett., 134 (2022), 108354. https://doi.org/10.1016/j.aml.2022.108354 doi: 10.1016/j.aml.2022.108354
    [33] S. Shu, H. Yu, Y. Huang, C. Nie, A symmetric finite volume element scheme on quadrilateral grids and superconvergence, Int. J. Numer. Anal. Mod., 3 (2006), 348–360.
    [34] Q. Hong, J. Wu, Coercivity results of a modified Q_{1}-finite volume element scheme for anisotropic diffusion problems, Adv. Comput. Math., 44 (2018), 897–922. https://doi.org/10.1007/s10444-017-9567-3 doi: 10.1007/s10444-017-9567-3
    [35] F. Fang, Q. Hong, J. Wu, Analysis of a special Q_{1}-finite volume element scheme for anisotropic diffusion problems, Numer. Math. Theor. Meth. Appl., 12 (2019), 1141–1167. https://doi.org/10.4208/nmtma.OA-2018-0080 doi: 10.4208/nmtma.OA-2018-0080
    [36] S. Chou, S. He, On the regularity and uniformness conditions on quadrilateral grids, Comput. Methods Appl. Mech. Eng., 191 (2002), 5149–5158. https://doi.org/10.1016/S0045-7825(02)00357-2 doi: 10.1016/S0045-7825(02)00357-2
    [37] P. Ciarlet, The finite element method for elliptic problems, Amsterdam: North-Holland, 1978.
    [38] D. Kershaw, Differencing of the diffusion equation in Lagrangian hydrodynamic codes, J. Comput. Phys., 39 (1981), 375–395. https://doi.org/10.1016/0021-9991(81)90158-3 doi: 10.1016/0021-9991(81)90158-3
    [39] G. Yuan, Z. Sheng, Monotone finite volume schemes for diffusion equations on polygonal meshes, J. Comput. Phys., 227 (2008), 6288–6312. https://doi.org/10.1016/j.jcp.2008.03.007 doi: 10.1016/j.jcp.2008.03.007
  • This article has been cited by:

    1. Rishi Kumar Pandey, Kottakkaran Sooppy Nisar, Enhanced numerical techniques for solving generalized rotavirus mathematical model via iterative method and ρ-Laplace transform, 2024, 12, 26668181, 100963, 10.1016/j.padiff.2024.100963
    2. Huda Alsaud, Muhammad Owais Kulachi, Aqeel Ahmad, Mustafa Inc, Muhammad Taimoor, Investigation of SEIR model with vaccinated effects using sustainable fractional approach for low immune individuals, 2024, 9, 2473-6988, 10208, 10.3934/math.2024499
    3. Cicik Alfiniyah, Tutik Utami, Nashrul Millah, Reuben Iortyer Gweryina, Optimal control and stability analysis of an alcoholism model with treatment centers, 2025, 22150161, 103311, 10.1016/j.mex.2025.103311
  • Reader Comments
  • © 2023 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(1495) PDF downloads(116) Cited by(1)

Figures and Tables

Figures(5)  /  Tables(12)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog