Processing math: 84%
Research article

Existence of solutions to mixed local and nonlocal anisotropic quasilinear singular elliptic equations

  • Received: 14 July 2023 Revised: 14 August 2023 Accepted: 15 August 2023 Published: 24 August 2023
  • MSC : 35J67, 35R11

  • In this paper, we consider the existence of positive solutions to mixed local and nonlocal singular quasilinear singular elliptic equations

    {Δpu(x)+(Δ)spu(x)=f(x)u(x)δ,xΩ,u(x)>0,xΩ,u(x)=0,xRNΩ,

    where Ω is a bounded smooth domain of RN(N>2), Δpu is an anisotropic p-Laplace operator, p=(p1,p2,...,pN) with 2p1p2pN, (Δ)sp is the fractional p-Laplace operator. The major results shows the interplay between the summability of the datum f(x) and the power exponent δ in singular nonlinearities.

    Citation: Labudan Suonan, Yonglin Xu. Existence of solutions to mixed local and nonlocal anisotropic quasilinear singular elliptic equations[J]. AIMS Mathematics, 2023, 8(10): 24862-24887. doi: 10.3934/math.20231268

    Related Papers:

    [1] Xiangrui Li, Shuibo Huang, Meirong Wu, Canyun Huang . Existence of solutions to elliptic equation with mixed local and nonlocal operators. AIMS Mathematics, 2022, 7(7): 13313-13324. doi: 10.3934/math.2022735
    [2] CaiDan LaMao, Shuibo Huang, Qiaoyu Tian, Canyun Huang . Regularity results of solutions to elliptic equations involving mixed local and nonlocal operators. AIMS Mathematics, 2022, 7(3): 4199-4210. doi: 10.3934/math.2022233
    [3] Xicuo Zha, Shuibo Huang, Qiaoyu Tian . Uniform boundedness results of solutions to mixed local and nonlocal elliptic operator. AIMS Mathematics, 2023, 8(9): 20665-20678. doi: 10.3934/math.20231053
    [4] Chao Wang, Fajie Wang, Yanpeng Gong . Analysis of 2D heat conduction in nonlinear functionally graded materials using a local semi-analytical meshless method. AIMS Mathematics, 2021, 6(11): 12599-12618. doi: 10.3934/math.2021726
    [5] Tomas Godoy . Existence of positive weak solutions for a nonlocal singular elliptic system. AIMS Mathematics, 2019, 4(3): 792-804. doi: 10.3934/math.2019.3.792
    [6] Luis M. Pedruelo-González, Juan L. Fernández-Martínez . Generalization of Snell's Law for the propagation of acoustic waves in elliptically anisotropic media. AIMS Mathematics, 2024, 9(6): 14997-15007. doi: 10.3934/math.2024726
    [7] Jiaqing Hu, Xian Xu, Qiangqiang Yang . Bifurcation results of positive solutions for an elliptic equation with nonlocal terms. AIMS Mathematics, 2021, 6(9): 9547-9567. doi: 10.3934/math.2021555
    [8] Mahmoud Soliman, Hamdy M. Ahmed, Niveen Badra, M. Elsaid Ramadan, Islam Samir, Soliman Alkhatib . Influence of the $ \beta $-fractional derivative on optical soliton solutions of the pure-quartic nonlinear Schrödinger equation with weak nonlocality. AIMS Mathematics, 2025, 10(3): 7489-7508. doi: 10.3934/math.2025344
    [9] Mengyu Wang, Xinmin Qu, Huiqin Lu . Ground state sign-changing solutions for fractional Laplacian equations with critical nonlinearity. AIMS Mathematics, 2021, 6(5): 5028-5039. doi: 10.3934/math.2021297
    [10] Naveed Iqbal, Azmat Ullah Khan Niazi, Ikram Ullah Khan, Rasool Shah, Thongchai Botmart . Cauchy problem for non-autonomous fractional evolution equations with nonlocal conditions of order $ (1, 2) $. AIMS Mathematics, 2022, 7(5): 8891-8913. doi: 10.3934/math.2022496
  • In this paper, we consider the existence of positive solutions to mixed local and nonlocal singular quasilinear singular elliptic equations

    {Δpu(x)+(Δ)spu(x)=f(x)u(x)δ,xΩ,u(x)>0,xΩ,u(x)=0,xRNΩ,

    where Ω is a bounded smooth domain of RN(N>2), Δpu is an anisotropic p-Laplace operator, p=(p1,p2,...,pN) with 2p1p2pN, (Δ)sp is the fractional p-Laplace operator. The major results shows the interplay between the summability of the datum f(x) and the power exponent δ in singular nonlinearities.



    Our main purpose of this study is to investigate the existence of positive solutions to the following mixed local and nonlocal quasilinear singular elliptic equation

    {Δpu(x)+(Δ)spu(x)=f(x)u(x)δ,xΩ,u(x)>0,xΩ,u(x)=0,xRNΩ, (1.1)

    where Ω is a bounded smooth domain of RN(N>2), Δpu is an anisotropic version of the p-Laplace operator, which is sometimes referred as the pseudo p-Laplace operator,

    Δpu(x)=Ni=1i[|iu(x)|pi2iu(x)],

    where p=(p1,p2,...,pN), pi2 for all i=1,2,...,N. The fractional p-Laplace operator (Δ)sp, (s(0,1),p1) is defined by

    (Δ)spu(x)=P.V.RN|u(x)u(y)|p2(u(x)u(y))|xy|N+psdy,

    where P.V. denotes the Cauchy principal value.

    Recently, there has been increasing attention focused on the study of elliptic operators that involve mixed local and nonlocal operators.These equations often arise spontaneously in the study of plasma physics and population dynamics[1,2]. For some other related results of mixed local and nonlocal equation, see [3,4,5,6,7,8] and the references therein. In the nonlocal case (0<s<1), Barrios et al. [9] investigated the existence and uniqueness results of positive solutions to the following problem with p=2,

    {(Δ)spu(x)=f(x)u(x)δ,xΩ,u(x)>0,xΩ,u(x)=0,xRNΩ. (1.2)

    In the case δ>0, the existence of solutions to problem (1.2) obtained by the range of δ to the summability of f. In case 0<δ<1 and 1<δ, Youssfi and Mahmoud [10] studied the existence of solutions to problem (1.2) with p=2 under some suitable assumptions on the datum f. For further information, readers may refer to the related work [11,12,13] and references therein.

    In the local case, Boccardo and Orsina [14] used the method of approximation to prove the existence of solutions to following the problem with p=2,

    {Δpu(x)=f(x)u(x)δ,xΩ,u(x)>0,xΩ,u(x)=0,xRNΩ. (1.3)

    They also studied the summability of these solutions when δ(0,). Giacomoni and Schindler [15] employ variational methods proved the existence of solution to quasilinear elliptic problem for pi=p(1,) with δ(0,1). During the past few years, there has been a vast amount of literature devoted to studying the anisotropic operator, which has numerous applications in fluid dynamics and physical phenomena, (we refer readers to [16,17,18,19] and references therein). Miri [19] further extended some results of [14] to an anisotropic quasilinear singular elliptic problem with variable exponent δ(x), and obtained existence of a solution to this problem. Bal and Garain [20, Theorems 2.7 and Theorems 2.9] established existence and uniqueness of solutions to the following mixed singular problems

    {L1u(x)=f1(x)u(x)δ+g1(x)u(x)γ,xΩ,u(x)>0,xΩ,u(x)=0,xΩ, (1.4)

    and

    {L2u(x)=f2(x)u(x)δ+g2(x)u(x)γ,xΩ,u(x)>0,xΩ,u(x)=0,xΩ, (1.5)

    where Ω is a bounded smooth subset of RN, N>2, δ>0, γ>1, fj, gj (j=1, 2) are nonnegative integrable functions,

    L1u(x)=div[w(x)|u(x)|p2u(x)],L2u(x)=Ni=1i[|iu(x)|pi2iu(x)].

    When gj=0 (j=1,2), they obtained a solution to problems (1.4) and (1.5) associated with the following minimizing problems

    v1(Ω):=infuW1,p{0}(Ω,ω){Ω|u|pωdx:Ω|u|1δf1dx=1},

    and

    v2(Ω):=infuW1,p{0}(Ω){Ni=1Ω|iu|pidx:Ω|u|1δf2dx=1}.

    Garain and Ukhlov [21, Theorems 2.13] proved the existence of solution to the following problem

    {Δpu(x)+(Δ)spu(x)=f(x)u(x)δ,xΩ,u(x)>0,xΩ,u(x)=0,xRNΩ. (1.6)

    It has been shown that problem (1.6) has a weak solution uW1,p{0}(Ω) when δ(0,1] and fLm(Ω){0} with m=(p1δ), where p=NpNp, while if δ(1,) and fL1(Ω){0}, then problem (1.6) has a weak solution uW1,p{loc}(Ω) with up+δ1pW1,p{0}(Ω). Moreover, they proved that mixed Sobolev inequalities are both necessary and sufficient for the existence of weak solutions to such singular problems. For related results about mixed local and nonlocal elliptic operators see [22,23,24,25,26,27,28,29,30] and references therein.

    Motivated by the results of the above cited papers, especially [20,21], the our purpose of this study is to establish the existence of solutions to problem (1.1) according to the range of the power exponent δ and to the summability of datum f(x). The main results as follows:

    Theorem 1.1. Let 0<δ<1 and 1<ˉp<N. Suppose that f>0, fLm(Ω) with m1. Then there exists a weak solution u to problem (1.1) such that

    (i) uL(Ω) if m>NˉpNˉp(Nˉp)pN, where ˉp satisfies

    1ˉp=1NNi=11pi.

    (ii) uLt(Ω) if ˉm<m<NˉpNˉp(Nˉp)pN, where

    ˉm=NˉpNˉppi(Nˉp)(1δpi)(Nˉp),t=m(1δpi)NˉpNˉp(m1)mpi(Nˉp).

    (iii) uW1,q0(Ω) if 1m<ˉm, where

    q=pim(1δpi)Nˉpm(1δ)[Nˉp(Nˉp)pi]piNˉp.

    Remark 1.2. Notice that when pi=2, the range of corresponding m values is exactly the summability of solutions obtained in [14].

    When pi=p, then problem (1.1) reduces to problem (1.6). Therefore

    (ⅰ) If fLm(Ω) with NpNp(1δ)(Np)<m<Np, then the solutions u to problem (1.6) satisfies uLt(Ω) with

    t=mN(1δpi)mpN.

    (ⅱ) If fLm(Ω) with 1m<NpNp(1δ)(Np), then the solutions u to problem (1.6) satisfies uW1,q{0}(Ω) with q=mN(1δp)m(1δ)N.

    Theorem 1.3. Suppose that δ=1 and 1<ˉp<N, f>0, fLm(Ω) with m>1. Then there exists a weak solution u to problem (1.1) such that

    (i) uL(Ω) if m>NˉpNˉp(Nˉp)pN.

    (ii) uLt(Ω) if 1m<NˉpNˉp(Nˉp)pN, where

    t=mpiNˉpmpi(Nˉp)Nˉp(m1).

    Theorem 1.4. Let δ>1 and 1<ˉp<N. Suppose that f>0, fLm(Ω) with m>1. Then there exists a weak solution u to problem (1.1) such that

    (i) uL(Ω) if m>NˉpNˉp(Nˉp)pN.

    (ii) uLt(Ω) if 1m<NˉpNˉp(Nˉp)pN, where

    t=m(1δpi)NˉpNˉp(m1)mpi(Nˉp).

    The order of the article is organized as follows: In Section 2, we provide basic notations and algebraic inequalities needed in this paper, as well as some definitions and useful lemmas. In Section 3, we present the proofs of Theorem 1.1, Theorem 1.3 and Theorem 1.4.

    In this article, we will use the following notations:

    For any v, we denote by v+=max{v,0}, v=max{v,0}. For p>1, we denote by p=pp1 to mean the conjugate exponent of p.

    Definition 2.1. Let p>1, ΩRN with N>2. The fractional Sobolev space Ws,p(Ω) is defined by

    Ws,p(Ω)={uLp(Ω):|u(x)u(y)||xy|Np+sLp(Ω×Ω)},

    with

    uWs,p(Ω)=(Ω|u(x)|pdx+ΩΩ|u(x)u(y)|p|xy|N+psdxdy)1p.

    The space Ws,p(RN) and Ws,ploc(Ω) are defined analogously. The space Ws,p0(Ω) is defined as

    Ws,p0(Ω)={uWs,p(RN):u=0onRNΩ}.

    Both Ws,p(Ω) and Ws,p0(Ω) are reflexive Banach spaces [31].

    Recall that the Lebesgue space Lpi(E) is defined as the space of pi-integrable functions u:ER with the finite norm

    uLpi(E)=(E|u(x)|pidx)1pi,

    where pi(1,+) for all i=1,2,...,N.

    The anisotropic Sobolev space is defined as follows:

    W1,pi(Ω)={uW1,1(Ω):iuLpi(Ω)},

    and

    W1,pi0(Ω)={uW1,10(Ω):iuLpi(Ω)},

    endowed with the norm

    uW1,pi0(Ω)=Ni=1iuLpi(Ω). (2.1)

    Definition 2.2. A function J: W1,pi0R is defined to be weakly lower semi-continuous if

    J(u)lim infnJ(un),

    for any sequence un approaching uW1,pi0 in the weak topology on W1,pi0.

    The zero Dirichlet boundary condition in this paper is defined as follows:

    Definition 2.3. We say that u0 in RNΩ if u=0 in RNΩ and for any ϵ>0, we have

    (uϵ)+W1,p0(Ω).

    We say that u=0 on RNΩ, if u is nonnegative and u0 in RNΩ.

    The definition of weak solution in this paper is defined as

    Definition 2.4. A positive function uW1,piloc(Ω)Lpi1(RN) is a weak solution to problem (1.1) if

    u>0inΩ,u=0inRNΩ,f(x)uδL1loc(Ω),

    for every ϕC1c(Ω), we have that

    Ni=1Ω|iu(x)|pi2iu(x)iϕdx+D(Ω)Ku(x,y)(ϕ(x)ϕ(y))dμ=Ωf(x)u(x)δϕdx, (2.2)

    where

    D(Ω)=RN×RN(Ωc×Ωc),

    and

    Ku(x,y)=|u(x)u(y)|p2(u(x)u(y)),dμ=|xy|Npsdxdy.

    Lemma 2.5. [19, Theorem 1.2] There exists a positive constant C, such that for every uW1,pi(Ω), we have

    upNLˉp(Ω)CNi=1iupiLpi(Ω), (2.3)

    where

    1ˉp=1NNi=11pi,

    and

    ˉp=NˉpNˉp.

    Lemma 2.6. [32, Lemma 2.1] Let 1<pi<. Then for ξ, ηRN, there exists a constant C=C(pi)>0 such that

    |ξ|pi2ξ|η|pi2η,ξη{cpi|ξη|pi,if pi2,cpi|ξη|2(|ξ|+|η|)2pi,if 1<pi<2. (2.4)

    Lemma 2.7. [33, Lemma 2.1] Let f: R+R+ be a non-increasing function such that

    |ψ(h)|Mψ(k)α|hk|βforallh>k>0

    where M>0, α>1 and β>0. Then ψ(d)=0, where dβ=Cψ(0)α12αβ(α1).

    For nN, f(x)L1(Ω) and f(x)>0, let fn(x):=min{f(x),n} and we consider the following approximated problem

    {Δpun(x)+(Δ)spun(x)=fn(x)(u++1n)δ,xΩ,un(x)>0,xΩ,un(x)=0,xRNΩ. (2.5)

    First, we consider the following useful result.

    Lemma 2.8. Let g(x)L(Ω), g(x)0. Then the following elliptic problem

    {Δpu(x)+(Δ)spu(x)=g(x),xΩ,u(x)>0,xΩ,u(x)=0,xRNΩ, (2.6)

    has a unique positive weak solution uW1,pi0(Ω).

    Proof. Existence: Define the energy functional J:W1,pi0(Ω)R as

    J(u):=J1(u)+J2(u)J3(u),

    where

    J1(u)=Ni=11piΩ|iu(x)|pidx,
    J2(u)=1pD(Ω)|u(x)u(y)|p|xy|N+psdxdy,

    and

    J3(u)=Ωg(x)u(x)dx.

    (ⅰ) By the Sobolev embedding theorem and gL(Ω), we have

    J(v)1pivpiW1,pi0(Ω)|Ω|p1pgL(Ω)vLp(Ω)asvpiW1,pi0(Ω),

    which implies the J is coercive.

    (ⅱ) J(v) is weakly lower semi-continuous in W1,pi0(Ω).

    It easy to see that J1(u) is differentiable, according to [34, Lemma 3.4], we know that Ji(u),i=1,2 are weakly lower semi-continuous. On the other hand, we know that J3(u) is a bounded linear functional. Thus J3(u) is continuous. Therefore, J is weakly lower semi-continuous and

    lim infnJ(un)=lim infnJ1(un)+lim infnJ2(un)limnJ3(un)J1(u)+J2J3(u)=J(u).

    Combining the above properties of J, we know that there exists a minimizer uW1,pi0(Ω) and which is also a critical point of J, which also is the solution to Eq (2.6).

    Uniqueness: Let u1, u2W1,pi0(Ω) be two solutions to problem (2.6). Thus, for any φW1,pi0(Ω), we have

    Ni=1Ω|iu1|pi2iu1iφdx+D(Ω)Ku1(x,y)(φ(x)φ(y))dμ=Ωgφdx, (2.7)

    and

    Ni=1Ω|iu2|pi2iu2iφdx+D(Ω)Ku2(x,y)(φ(x)φ(y))dμ=Ωgφdx. (2.8)

    Choosing φ=u1u2 and then subtracting (2.7) and (2.8), we obtain

    Ni=1Ω(|iu1|pi2iu1|iu2|pi2iu2)(iu1iu2)dx+D(Ω)(Ku1(x,y)Ku2(x,y))[(u1u2)(x)(u1u2)(y)]dμ=0. (2.9)

    Using Lemma 2.6, we get the first term of the left hand side of (2.9) is nonnegative. On the other hand, by the monotonicity of the function f(t)=tp1(p>1), we have

    [Ku1(x,y)K(u2(x,y))][(u1u2)(x)(u1u2)(y)]=[|u1(x)u1(y)|p2(u1(x)u1(y))|u2(x)u2(y)|p2(u2(x)u2(y))][(u1(x)u1(x))(u2(y)u2(y))]0. (2.10)

    Consequently,

    Ni=1Ω(|iu1|pi2iu1|iu2|pi2iu2)(iu1iu2)dx=0. (2.11)

    Therefore, u1(x)u2(x)=C for all xRN. Note that u1u2=0 on RNΩ since ui(x)=0 for xRNΩ. Thus u1(x)u2(x), which implies that the solution of (2.6) is unique.

    Boundedness: For any k>1, decompose RN as RN=AkAck, where

    Ak={xΩ:u(x)k},Ack={xΩ:0<u(x)<k}.

    Taking Gk(u):=(uk)+=max{uk,0} as a test function in (2.6), we have

    Ni=1Ω|iu(x)|pi2iu(x)iGk(u(x))dx+D(Ω)Ku(x,y)[Gk(u(x))Gk(u(y))]dμ=Ωg(x)Gk(u(x))dx. (2.12)

    Obviously,

        Ku(x,y)[Gk(u(x))Gk(u(y))]=|u(x)u(y)|p2(u(x)u(y))[(u(x)k)+(u(y)k)+]={|u(x)u(y)|p,                         if      u(x)>k,u(y)>k,|u(y)u(x)|p1(u(y)k),  if      u(y)>ku(x),|u(x)u(y)|p1(u(x)k),  if      u(x)>ku(y),0,                                               if      u(x)k,u(y)k,0. (2.13)

    Therefore, combining (2.13) and (2.12) with Sobolev embedding theorem, we have

    Ni=1Ω|iGk(u)|pidx=Ni=1Ak|iGk(u)|pidx+Ni=1Ack|iGk(u)|pidx=Ni=1Ak|iGk(u)|pidxdxΩg(x)Gk(u)dxgL(Ω)(ΩGk(u)ˉpdx)1ˉp|A(k)|ˉp1ˉpCgL(Ω)(Ni=1Ω|iGk(u)|pidx)1pN|A(k)|ˉp1ˉp.

    Therefore

    Ni=1Ω|iGk(u)|pidxCgpNpN1L(Ω)|A(k)|pN(ˉp1)ˉp(pN1). (2.14)

    For every 1k<h we know that A(h)A(k) and u(x)k(hk) in A(h), we get

    (hk)pN|A(h)|pNp(A(h)Gk(u)ˉpdx)pNˉp(A(k)Gk(u)ˉpdx)pNˉpCNi=1Ω|iGk(u)|pidxCgpNpN1L(Ω)|A(k)|pN(ˉp1)ˉp(pN1).

    Hence, we have

    |A(h)|CgˉppN1L(Ω)(hk)ˉp|A(k)|ˉp1pN1.

    Obviously,

    ˉp>pN.

    Hence, using Lemma 2.7 we obtain

    uL(Ω)C.

    Positivity: First, taking u(x):=min{u(x),0} as a test function in (2.6) and using g0, we have

    Ni=1Ω|iu(x)|pidx+R2NKu(x,y)(u(x)u(y))dμ=Ωgudx0, (2.15)

    where Ku(x,y)=|u(x)u(y)|p2(u(x)u(y)). Rewrite

    RN×RN=4i=1Ai.

    Denote,

    A1={(x,y)RN×RN:u(x)0,u(y)0},A2={(x,y)RN×RN:u(x)0,u(y)<0},A3={(x,y)RN×RN:u(x)<0,u(y)0},A4={(x,y)RN×RN:u(x)<0,u(y)<0}.

    Therefore,

    Ku(x,y)(u(x)u(y))=|u(x)u(y)|p2(u(x)u(y))(u(x)u(y))={0,if(x,y)A1,|u(y)|p,if(x,y)A2,|u(x)|p,if(x,y)A3,|u(x)u(y)|p,if(x,y)A4,0.

    Obviously,

    |u(x)u(y)|p2(u(x)u(y))(u(x)u(y))0. (2.16)

    Using (2.16) in (2.15) we obtain

    Ni=1Ω|iu|pidx=0.

    Therefore, u=C for all xRN. Note that u=0 on RNΩ since u:=min{u,0}. Thus u0 in Ω.

    Second, assume that there exists a point x0Ω such that u(x0)=infxΩu(x)=0, thus

    Ni=1Ωi[|iu(x0)|pi2iu(x0)]+D(Ω)|u(x0)u(y)|p2[u(x0)u(y)]|xy|N+psdy=D(Ω)|u(y)|p2[u(y)]|xy|N+psdy=R2N|u(y)|p1|xy|N+psdy<0.

    This is a contradiction since g(x0)0. Hence, u>0 in Ω.

    Lemma 2.9. For any nN, there exists a unique positive solution unW1,pi0(Ω)L(Ω) to problem (2.5). Moreover, The sequence {un} is increasing with respect to n and

    un(x)CK>0forKΩ.

    Proof. Step1. (Existence) Let nN. By Lemma 2.8, for every uW1,pi0(Ω)L(Ω), there exists a unique vW1,pi0(Ω)L(Ω) such that

    {Δpv(x)+(Δ)spv(x)=fn(x)(u+1n)δ,xΩ,v(x)>0,xΩ,v(x)=0,xRNΩ. (2.17)

    Define the operator T:uv=T(u), where v is the unique solution to (2.18). Choosing v as a test function in (2.18), using Sobolev imbedding theorem, we obtain

    Ni=1Ω|iv(x)|pidxΩnδ+1v(x)dxCnδ+1|Ω|ˉp1ˉp(Ni=1Ω|iv(x)|pidx)1pN.

    Thus

    Ni=1Ω|iv(x)|pidxCnpN(δ+1)pN1|Ω|pN(ˉp1)ˉp(pN1):=R, (2.18)

    which implies that the ball with radius R in W1,pi0(Ω) remains unchanged under T.

    Now, we have to prove the continuity and compactness of T, which is an operator from W1,pi0(Ω) to W1,pi0(Ω).

    (ⅰ) Continuity of T: In order to do this, we have to show that limkvkvW1,pi0(Ω)=0 if limkukuW1,pi0(Ω)=0, where vk=T(uk) and v=T(u).

    Choosing ˉvk(x)=vk(x)v(x) as a test function of the equations of vk and v respectively, using (2.19), we get

    Ω|iˉvk(x)|pidxNi=1Ω|i(vk(x)v(x))|pidxNi=1Ω[|ivk(x)|pi2ivk(x)|iv(x)|pi2iv(x)][ivk(x)iv(x)]dx+D(Ω)[Kvk(x,y)Kv(x,y)][(vkv)(x)(vkv)(y)]dμ=Ω[fn(x)(vk+1n)δf(x)vδ][vk(x)v(x)]dx,ifpi2. (2.19)

    Using Hölder and Sobolev inequalities we infer that

    |Ω[fn(x)(uk+1n)δf(x)uδ][vk(x)v(x)]dx|[Ω|fn(x)(uk+1n)δf(x)uδ|pidx]1piˉvkLpi(Ω)C[Ω|fn(x)(uk+1n)δf(x)uδ|pidx]1piiˉvkLpi(Ω). (2.20)

    By (2.19) and (2.20), using [10, Lemma 2.2], we find

    iˉvkLpi(Ω)Cn1pi1[Ω|1(uk+1n)δ1uδ|pidx]1(pi1)piCn1pi1[Ω|uδ(uk+1n)δ(uk+1n)δuδ|pidx]1(pi1)piCn1pi1+δ[Ω|uuk|pidx]1(pi1)pi, (2.21)

    since the pointwise convergence of uku in W1,pi0(Ω). we get

    limk+vkvW1,pi0(Ω)=0.

    Therefore, in the case pi2, the operator T is continuous from W1,pi0(Ω) to W1,pi0(Ω).

    (ⅱ) Compactness of T: To achieve this, we have to show that, for some vW1,pi0(Ω), it holds

    limk+vkvW1,pi0(Ω)=0.

    Let uk be a bounded sequence in W1,pi0(Ω) and vk:=T(uk). Then we have

    ukuinW1,pi0(Ω),ukuinLt(Ω),1<t<ˉp.

    According to (2.18), we have T(uk)W1,pi0(Ω)C. Therefore there exists a subsequence, still denoted by {uk}, such that

    T(uk)vW1,pi0(Ω),T(uk)vLt(Ω),1<t<ˉp.

    For any φW1,pi0(Ω),

    Ni=1Ω|ivk(x)|pi2ivk(x)iφdx+D(Ω)|vk(x)vk(y)|p2(vk(x)vk(y))(φ(x)φ(y))|xy|N+psdxdy=Ωfn(x)(uk+1n)δφdx. (2.22)

    Now, we show that as k, (2.22) converges to

    Ni=1Ω|iv(x)|pi2iv(x)iφdx+D(Ω)|v(x)v(y)|p2(v(x)v(y))(φ(x)φ(y))|xy|N+psdxdy=Ωfn(x)(u+1n)δφdx. (2.23)

    By the dominated convergence theorem, we have

    limkΩfn(x)(uk+1n)δφdx=Ωfn(x)(u+1n)δφdx,

    and

    Ni=1ivkNi=1ivpointwise almost everywhere in Ω.

    Therefore, for every φC1c(Ω), we have

    limkNi=1Ω|ivk(x)|pi2ivk(x)iφdx=Ni=1Ω|iv(x)|pi2iv(x)iφdx.

    Since φC1c(Ω) and vk is uniformly bounded in W1,pi0(Ω),

    {|vk(x)vk(y)|p2(vk(x)vk(y))|xy|N+psp}nNLp(RN×RN),

    by the pointwise convergence of vk(x) to v(x)

    |vk(x)vk(y)|p2(vk(x)vk(y))|xy|N+psp|v(x)v(y)|p2(v(x)v(y))|xy|N+psp a.e.inR2N.

    Since

    φ(x)φ(y)|xy|N+pspLp(R2N),

    we get that the (2.22) converges to the (2.23). Similarly, combining (2.21) and (2.18), we have

    limk+T(uk)T(u)W1,pi0(Ω)=0.

    Therefore, the operator T is continuous from W1,pi0(Ω) to W1,pi0(Ω). Then, Schauder fixed point theorem implies the existence of a fixed points un such that un=T(un), which is a weak solution to approximated problem (2.5).

    Step2. (Monotonicity) Since un and un+1 are positive solutions to problem (2.8), for any φW1,pi0(Ω), we have

    Ni=1Ω|iun(x)|pi2iun(x)iφdx+D(Ω)Kun(x,y)(φ(x)φ(y))dμ=Ωfn(x)(un(x)+1n)δφdx, (2.24)

    and

    Ni=1Ω|iun+1(x)|pi2iun+1(x)iφdx+D(Ω)Kun+1(x,y)(φ(x)φ(y))dμ=Ωfn+1(x)(un+1(x)+1n+1)δφdx. (2.25)

    Taking φ=(un(x)un+1(x))+ as test function in (2.24) and (2.25), we get

    Ni=1Ω|iun(x)|pi2iun(x)i(un(x)un+1(x))+dx+D(Ω)Kun(x,y)[(un(x)un+1(x))+(un(y)un+1(y))+]dμ=Ωfn(x)(un(x)+1n)δ(un(x)un+1(x))+dx, (2.26)

    and

    Ni=1Ω|iun+1(x)|pi2iun+1(x)i(un(x)un+1(x))+dx+D(Ω)Kun+1(x,y)[(un(x)un+1(x))+(un(y)un+1(y))+]dμ=Ωfn+1(x)(un+1(x)+1n+1)δ(un(x)un+1(x))+dx. (2.27)

    Since fn(x)fn+1(x) for xΩ, we have

    Ω[fn(x)(un(x)+1n)δfn+1(x)(un+1(x)+1n+1)δ](un(x)un+1(x))+dxΩfn+1(x)[1(un(x)+1n)δ1(un+1(x)+1n+1)δ](un(x)un+1(x))+dx=Ωfn+1(x)[(un+1(x)+1n+1)δ(un(x)+1n)δ(un(x)+1n)δ(un+1(x)+1n+1)δ](un(x)un+1(x))+dx0. (2.28)

    Subtracting (2.26) with (2.27) and using the (2.28), we obtain

    Ni=1Ω[|iun(x)|pi2iun(x)|iun+1(x)|pi2iun+1(x)]i(un(x)un+1(x))+dx+D(Ω)(Kun(x,y)Kun+1(x,y))[(un(x)un+1(x))+(un(y)un+1(y))+]dμ0. (2.29)

    Following the argument in the proof of [35, Lemma 9], we obtain

    D(Ω)(Kun(x,y)Kun+1(x,y))[(un(x)un+1(x))+(un(y)un+1(y))+]dμ0. (2.30)

    Therefore, applying (2.30) in (2.29) we get

    Ni=1Ω[|iun(x)|pi2iun(x)|iun+1(x)|pi2iun+1(x)]i(un(x)un+1(x))+dx0.

    Using Lemma 2.6 we obtain

    (un(x)un+1(x))+=CforallxRN.

    Note that un(x)=un+1(x)=0 on RNΩ thus C=0, which implies that un+1(x)un(x) in Ω.

    Step3. (UniformPositivity) Let u1 solves (2.6). By Lemma 2.8, for every KΩ, there exists a constant CK>0 such that u1CK>0 in K. Again, since the monotonicity of un, we have unu1 in K. Therefore, for any KΩ,

    un(x)CK>0,forxK.

    In order to prove the existence of positive solution to (1.1), we use the sequence of solutions un of problem (2.5). Then we need a priori estimates on un.

    Lemma 3.1. Let 0<δ<1 and 1<ˉp<N. Suppose that f>0, fLm(Ω) with

    m>ˉm=NˉpNˉppi(Nˉp)(1δpi)(Nˉp).

    Then, the sequence solutions un to the approximate problem (2.5) such that

    (i) unL(Ω) if m>NˉpNˉppN(Nˉp).

    (ii) unLt(Ω), where t=m(1δpi)NˉpNˉp(m1)pim(Nˉp) if

    NˉpNˉppi(Nˉp)(1δpi)(Nˉp)<m<NˉpNˉppN(Nˉp).

    Proof. (i) Let Ak={xΩ:un(x)k}. Choosing Gk(u):=(uk)+W1,pi0(Ω) as a test function in (2.5), we get

    Ni=1Ω|iun(x)|pi2iun(x)iGk(un(x))dx+D(Ω)Kun(x,y)[Gk(un(x))Gk(un(y))]dμ=Ωfn(x)(un(x)+1n)δGk(un(x))dx. (3.1)

    Foe any k>1, by (2.13) we know that,

    Kun(x,y)[Gk(un(x))Gk(un(y))]0.

    By Hölder inequality, Sobolev embedding theorem, fn(x)f(x) and (3.1), we have

    Ni=1Ω|iGk(un)|pidx=Ni=1Ak|iGk(un)|pidxΩfn(x)(un(x)+1n)δGk(un)dxAkf(x)G1δk(un)dx(Akf(x)mdx)1m(ΩGk(un)ˉpdx)1δˉp|A(k)|11m1δˉpC(Akf(x)mdx)1m(Ni=1Ω|iGk(un)|pidx)1δpN|A(k)|11m1δˉp. (3.2)

    Hence

    Ni=1Ω|iGk(un)|pidxC(Akf(x)mdx)pNm(pN+δ1)|A(k)|(11m1δˉp)pNpN+δ1.

    Let h>k1, we know that AhAk and Gk(un)hk for in Ω, we have that

    |hk|pN|Ah|pNˉp(A(h)Gk(un)ˉpdx)pNˉp(A(k)Gk(un)ˉpdx)pNˉpCNi=1A(k)|iGk(un)|pidxCfpNpN+δ1Lm(Ω)|A(k)|(11m1δˉp)pNpN+δ1.

    Therefore

    |Ah|CfˉppN+δ1Lm(Ω)|Ak|(11m1δˉp)ˉppN+δ1|hk|ˉp.

    Note that

    (11m1δˉp)ˉppN+δ1>1,

    if m>NˉpNˉppN(Nˉp). Hence, apply Lemma 2.7 with

    M=CfˉppN+δ1Lm(Ω)>0,α=(11m1δˉp)ˉppN+δ1>1,β=ˉp>0andψ(k)=|Ak|,

    there exists k0 such that ψ(k)0 for all kk0. Thus,

    esssupΩuk0.

    (ⅱ) Choose upi(γ1)+1n (γ>1) as test function in (2.17), we have

    Ni=1Ω|iun|pi2iuniupi(γ1)+1ndx+D(Ω)Kun(x,y)[un(x)pi(γ1)+1un(y)pi(γ1)+1]dμ=Ωfn(x)(un(x)+1n)δun(x)pi(γ1)+1dx. (3.3)

    According to [10, Lemma 2.2], we have

    Kun(x,y)[un(x)pi(γ1)+1un(y)pi(γ1)+1]=|un(x)un(y)|p2(un(x)un(y))[un(x)pi(γ1)+1un(y)pi(γ1)+1]C[un(x)+un(y)]pi(γ1)|un(x)un(y)|p0. (3.4)

    Combining (3.4) and (3.3), and using Hölder inequality, we get

    Ni=1Ω|iun|pi2iuniupi(γ1)+1ndx=Ni=1[pi(γ1)+1]Ω|iun|piupi(γ1)ndxΩf(x)(un+1n)δupi(γ1)+1ndxΩf(x)upi(γ1)+1δndxfLm(Ω)(Ωu[pi(γ1)+1δ]mndx)1m.

    By Sobolev inequality,

    Ni=1Ω[pi(γ1)+1]|iun|piupi(γ1)ndx=Ni=1Ω[pi(γ1)+1](1γ)pi|iuγn|pidxC(Ωuγˉpn)pNˉp.

    Therefore,

    (Ωuγˉpn)pNˉpCfLm(Ω)(Ωun(x)[pi(γ1)+1δ]m)1m. (3.5)

    Now we choose γ such that

    γˉp=[pi(γ1)+1δ]m,

    that is

    γ=m(1δpi)(Nˉp)Nˉp(m1)pim(Nˉp).

    Since γ>1, we know

    NˉpNˉppi(Nˉp)(1δpi)(Nˉp)<m.

    Thus pNˉp>1m gives

    (Ωun(x)γˉp)pNˉp1mCfLm(Ω).

    Therefore, un is uniformly bounded in Lt(Ω) with t=γˉp.

    Lemma 3.2. Let 0<δ<1 and 1<ˉp<N. Suppose that f>0 and fLm(Ω) with

    1m<NˉpNˉppi(Nˉp)(1δpi)(Nˉp).

    Then, the sequence solutions {un} to the approximate problem (2.5) are uniformly bounded in W1,q{0}(Ω) with

    q=pim(1δpi)Nˉpm(1δ)[Nˉp(Nˉp)pi]piNˉp.

    Proof. Similar to above taking upi(γ1)+1n as test function in (2.5) with δ+pi1piγ<1. However, this option is not acceptable, since the gradient of such a test function will be singular where un(x)=0. Hence, for n fixed, choose (un+ε)pi(γ1)+1εpi(γ1)+1 (0<ε<1n) as test function in (2.5), we get

    Ni=1[pi(γ1)+1]Ω|iun|pi(un+ε)pi(γ1)dxΩf(x)[(un+ε)pi(γ1)+1εpi(γ1)+1](un+1n)δdx.

    By fn(x)f(x) and ε<1n, we have

    Ni=1[p1(γ1)+1]Ω|iun|pi(un+ε)pi(γ1)dxΩf(x)(un+ε)pi(γ1)+1δdx. (3.6)

    By Sobolev inequality,

    Ni=1[pi(γ1)+1]Ω|iun|pi(un+ε)pi(γ1)dx=Ni=1[pi(γ1)+1]Ω(1γ)pi|i[(un+ε)γεγ]|pidxCNi=1Ω|i[(un+ε)γεγ]|pidxC(Ω[(un+ε)γεγ]ˉpdx)pNˉp. (3.7)

    Hence, by (3.6) and (3.7), we get

    (Ω[(un+ε)γεγ]ˉpdx)pNˉpCΩf(x)(un+ε)pi(γ1)+1δdx.

    Let ε0, we have

    (Ωuγˉpndx)ˉpˉpCΩf(x)upi(γ1)+1δndx. (3.8)

    If m=1, we choose γ=pi+δ1pi in the (3.8), so that unL(δ+pi1)Nˉppi(Nˉp)(Ω).

    If m>1, from the proof of Lemma 3.1, we get that un(x)Lt(Ω) with

    t=m(1δpi)NˉpNˉp(m1)pim(Nˉp).

    Since γ<1, by (3.6), we have

    Ni=1Ω|iun|pi(un+ε)pipiγdx=Ni=1Ω|iun|pi(un+ε)piγpidxC.

    We can apply Hölder inequality (since q<pi),

    Ω|iun|qdx=Ω|iun|q(un+ε)(1γ)q(un+ε)(1γ)qdx=Ω[|iun|pi(un+ε)(1γ)pi]qpi(un+ε)(1γ)qdx(Ω|iun|pi(un+ε)(1γ)pidx)qpi(Ω(un+ε)(1γ)qpipiqdx)1q¯piC(Ω(un+ε)(1γ)piqpiqdx)1q¯pi.

    Choice \gamma and q such that

    \begin{align*} \frac{(1-\gamma)p_{i}q}{p_{i}-q} = t. \end{align*}

    Therefore, u_{n}\in \mathcal{W}_{\text{{0}}}^{1, q}(\Omega) with

    \begin{align*} q = \frac{p_{i}m(1-\delta-p_{i})N{\bar{p}}}{m(1-\delta)\left[N\bar{p}-(N-\bar{p}){p_{i}} \right]-{p_{i}}N\bar{p}}. \end{align*}

    Lemma 3.3. Suppose that \delta = 1 and 1 < \bar{p} < N , f > 0 , f\in L^{m}(\Omega) with m > 1 . Then there exists a weak solution u_{n} to problem (2.5) such that

    (i) If m > \frac{N\bar{p}}{N{\bar{p}}-{p_{N}}(N-\bar{p})} , Then u_{n}\in {{L}^{\infty}(\Omega)} ;

    (ii) If 1\leq m < \frac{N\bar{p}}{N{\bar{p}}-{p_{N}}(N-\bar{p})} , Then u_{n}\in{{L}^{t}(\Omega)} , where

    \begin{align*} t = \frac{mp_{i}N\bar{p}}{p_{i}m(N-\bar{p})-N\bar{p}(m-1)}. \end{align*}

    Proof. The proof of (i) is identical to that of Lemma 3.1, so we will omit it.

    As for (ii) , observe that if m = 1 , then t = \frac{N\bar{p}}{N-\bar{p}} = \bar{p}^* . If m > 1 , similar to Lemma 3.1, Choosing u_{n}^{p_{i}(\gamma-1)+1} as test function in (2.5), we know that there is

    \begin{align*} \left(\int_{\Omega}u_{n}^{{\gamma\bar{p}^*}}dx\right)^{\frac{\bar{p}}{\bar{p}^*}} \leq&C\|f\|_{L^{m}(\Omega)}\left(\int_{\Omega}u_{n}(x)^{{p_{i}(\gamma-1)}{m^{\prime}}}\right) ^{\frac{1}{m^{\prime}}}. \end{align*}

    Choose \gamma such that

    \begin{align*} \gamma\bar{p}^* = \left[{p_{i}(\gamma-1)}\right]m^{\prime}. \end{align*}

    Obviously

    \begin{align*} \gamma = \frac{mp_{i}(N-\bar{p})}{mp_{i}(N-\bar{p})-N\bar{p}(m-1)}. \end{align*}

    Since \gamma > 1 , we arrive at 1 < m . Thus \frac{p_{N}}{\bar{p}^*} > \frac{1}{m^{\prime}} being

    \begin{align*} m < \frac{N\bar{p}}{N{\bar{p}}-{p_{N}}(N-\bar{p})}, \end{align*}

    so that u_{n}\in {L}^{t}(\Omega) with t = \gamma\bar{p}^* .

    Simple modifications to the proof of Lemma 3.1 enable us to demonstrate Lemma 3.4.

    Lemma 3.4. Suppose that \delta > 1 and 1 < \bar{p} < N , f > 0 , f\in L^{m}(\Omega) with m > 1 . Then there exists a weak solution u_{n} to problem (2.5) such that

    (i) If m > \frac{N\bar{p}}{N{\bar{p}}-{p_{N}}(N-\bar{p})} , then u_{n}\in {{L}^{\infty}(\Omega)} .

    (ii) If 1\leq m < \frac{N\bar{p}}{N{\bar{p}}-{p_{N}}(N-\bar{p})} , then u_{n}\in{{L}^{t}(\Omega)} with

    \begin{align*} t = \frac{m(1-\delta-p_{i})N\bar{p}}{N\bar{p}(m-1)-p_{i}m(N-\bar{p})}. \end{align*}

    Proof. The proof of (i) is identical to that given in Lemma 3.1, so we omit it.

    For (ii) , by [21, Lemma 3.7], we known, if m = 1 , the sequence u_{n}^{\frac{\delta+p_{i}-1}{p_{i}}} is uniformly bounded in \mathcal{W}_{0}^{1, p_{i}}(\Omega) , This also gives u_{n} is bounded in \mathcal{W}_{\text {loc}}^{1, p_{i}}(\Omega) .

    If 1 < m < \frac{N\bar{p}}{N{\bar{p}}-{p_{N}}(N-\bar{p})} , similar to Lemma 3.1, taking u_{n}^{p_{i}(\gamma-1)+1} as test function in (2.5), this time with \gamma > 1 since \gamma > {\frac{\delta+p_{i}-1}{p_{i}}} , we have

    \begin{align*} \left(\int_{\Omega}u_{n}^{{\gamma\bar{p}^*}}dx\right)^{\frac{\bar{p}}{\bar{p}^*}} \leq&C\|f\|_{L^{m}(\Omega)}\left(\int_{\Omega}u_{n}(x)^{[{p_{i}(\gamma-1)+ 1-\delta}]m^{\prime}}\right) ^{\frac{1}{m^{\prime}}}. \end{align*}

    Choosing \gamma in such a way that

    \begin{align*} \gamma\bar{p}^* = \left[{p_{i}(\gamma-1)+1-\delta}\right]m^{\prime}, \end{align*}

    since \gamma > {\frac{\delta+p_{i}-1}{p_{i}}} gives m > 1 , and by \frac{p_{N}}{\bar{p}^*} > \frac{1}{m^{\prime}} being

    \begin{align*} m < \frac{N\bar{p}}{N{\bar{p}}-{p_{N}}(N-\bar{p})}. \end{align*}

    Therefore, u_{n} is uniformly bounded in {{L}^{t}(\Omega)} as well.

    In this section, we give the proof of Theorem 1.1 by the approximate method.

    Proof of Theorem 1.1. Let f\in L^{m}(\Omega) . By Lemma 3.2 and 3.1, we know that the solutions u_{n} to problem (2.5) are bounded in \mathcal{W}_{0}^{1, p_{i}}(\Omega) . Then, the pointwise limit u in \mathcal{W}_{0}^{1, p_{i}}(\Omega)\cap{L}^{p_{i}-1}(\Omega) . For any \varphi\in\mathcal{W}_{0}^{1, p_{i}}(\Omega) ,

    \begin{align} &\sum\limits_{i = 1}^{N}\int_{\Omega}\left|{\partial_{i} u_{n}(x)}\right|^{p_{i}-2}{\partial_{i} u_{n}(x)}{\partial_{i}\varphi} dx+\int\int_{{\mathcal{D}(\Omega)}} \mathcal{K}u_{n}(x, y)(\varphi(x)-\varphi(y)) d\mu\\ = &\int_{\Omega}\frac{f_{n}(x)}{\left(u_{n}(x)+ \frac{1}{n}\right)^{\delta}}\varphi dx. \end{align} (3.9)

    Then, for any \varphi\in C_{c}^{1}(\Omega) , we get

    \begin{align} \lim _{n\rightarrow \infty}\sum\limits_{i = 1}^{N} \int_{\Omega}\left|{\partial_{i} u_{n}(x)}\right|^{p_{i}-2}{\partial_{i} u_{n}(x)} {\partial_{i}\varphi}dx = \sum\limits_{i = 1}^{N}\int_{\Omega} \left|{\partial_{i} u(x)} \right|^{p_{i}-2}{\partial_{i} u(x)}{\partial_{i}\varphi}dx. \end{align} (3.10)

    Since \left\{u_{n}\right\} is uniformly bounded in \mathcal{W}_{0}^{1, p_{i}}(\Omega) ,

    \begin{align*} \frac{\left|u_{n}(x)-u_{n}(y)\right|^{p-2}\left(u_{n}(x)- u_{n}(y)\right)}{|x-y|^{\frac{N+ps} {p^{\prime}}}}\in L^{p^{\prime}}\left(\mathbb{R}^{N}\times\mathbb{R}^{N}\right). \end{align*}

    By point-wise convergence of u_{n}(x) to u(x)

    \frac{\left|u_{n}(x)-u_{n}(y)\right|^{p-2}\left(u_{n}(x)- u_{n}(y)\right)}{|x-y|^{\frac{N+ps} {p^{\prime}}}}\rightarrow \frac{\left|u(x)-u(y)\right|^{p-2}\left(u(x)-u(y)\right)} {|x-y|^{\frac{N+ps} {p^{\prime}}}}\; \text{ a.e.} \; \text{ in }\; \mathbb{R}^{2N}.

    Then, we have

    \begin{align} \quad\lim _{n\rightarrow \infty}\int\int_{{\mathcal{D}(\Omega)}}\mathcal{K} u_{n}(x, y)(\varphi(x)-\varphi(y)) d\mu = \int\int_{{\mathcal{D}(\Omega)}}\mathcal{K}u(x, y) (\varphi(x)-\varphi(y))d\mu. \end{align} (3.11)

    By Lemma 2.9, for any K\Subset\Omega , u_{n}(x)\geq C_{K} > 0 with \operatorname{supp}(\varphi) = C_{K} > 0 . Therefore, for any \varphi\in C_{c}^{1}(\Omega) such that

    \begin{align*} \left|\frac{f_{n}(x)}{\left(u_{n}(x)+ \frac{1}{n}\right)^{\delta}}\varphi\right|\leq \frac{\|\varphi\|_{L^{\infty}(\Omega)}}{C_{K}^{\delta}}|f|\; \text{ in }\; \Omega. \end{align*}

    We conclude that

    \begin{align} \lim _{n\rightarrow \infty}\int_{\Omega}\frac{f_{n}}{\left(u_{n}(x)+ \frac{1}{n}\right)^{\delta}}\varphi dx = \int_{\Omega} \frac{f(x)}{u(x)^{\delta}}\varphi dx. \end{align} (3.12)

    Finally, passing to the limit in (3.9), we conclude that

    \begin{align*} \sum\limits_{i = 1}^{N}\int_{\Omega}\left|{\partial_{i} u(x)}\right|^{p_{i}-2}{\partial_{i} u(x)}{\partial\varphi}dx+&\int\int_{{\mathcal{D}(\Omega)}} \mathcal{K}u(x, y)(\varphi(x)-\varphi(y))d\mu = \int_{\Omega} \frac{f(x)}{\left(u(x)+\frac{1}{n}\right)^{\delta}}\varphi dx, \end{align*}

    for all \varphi\in C_{c}^{1}(\Omega) , which shows that u is a solution to problem (1.1) and u\in \mathcal{W}_{\text{{0}}}^{1, q}(\Omega) .

    Proof of Theorem 1.3 and Theorem 1.4. The proof of Theorem 1.3 and 1.4 are similar, here we omit the details.

    The manuscript establishes the existence of solutions to mixed local and nonlocal anisotropic quasilinear singular elliptic eqautions. The interplay between the integrability and the singularity power is investigated. This results generalizes and complements the existing results.

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

    This works was partially supported by Fundamental Research Funds for the Central Universities (No. 31920220041) and Innovation Team Project of Northwest Minzu University (No. 1110130131).

    The authors declare no competing interests.



    [1] D. Blazevski, D. Negrete, Local and nonlocal anisotropic transport in reversed shear magnetic fields Shearless Cantori and nondiffusive transport, Phys. Rev. E, 87 (2013), 063106. https://doi.org/10.1103/PhysRevE.87.063106 doi: 10.1103/PhysRevE.87.063106
    [2] S. Guo, Patterns in a nonlocal time-delayed reaction-diffusion equation, Z. Angew. Math. Phys., 69 (2018), 10. https://doi.org/10.1007/s00033-017-0904-7 doi: 10.1007/s00033-017-0904-7
    [3] S. Biagi, S. Dipierro, E. Valdinoci, E. Vecchi, Mixed local and nonlocal elliptic operators: Regularity and maximum principles, Commun. Partial Differ. Equ., 47 (2022), 585–629. https://doi.org/10.1080/03605302.2021.1998908 doi: 10.1080/03605302.2021.1998908
    [4] S. Biagi, S. Dipierro, E. Valdinoci, E. Vecchi. A Hong-Krahn-Szegö inequality for mixed local and nonlocal operators, Math. Engrg., 5 (2023), 1–25. https://doi.org/10.48550/arXiv.2110.07129 doi: 10.48550/arXiv.2110.07129
    [5] S. Dipierro, E. Valdinoci, Description of an ecological niche for a mixed local/nonlocal dispersal: An evolution equation and a new Neumann condition arising from the superposition of Brownian and Lévy processes, Phys. A, 575 (2021), 126052. https://doi.org/10.1016/j.physa.2021.126052 doi: 10.1016/j.physa.2021.126052
    [6] S. Dipierro, E. P. Lippi, E. Valdinoci, (Non)local logistic equations with neumann conditions, Ann. Inst. H. Poincaré Anal. Non Linéaire, 2021. https://doi.org/10.48550/arXiv.2101.02315 doi: 10.48550/arXiv.2101.02315
    [7] P. Garain, On a class of mixed local and nonlocal semilinear elliptic equation with singular nonlinearity, J. Geom. Anal., 33 (2023), 212. https://doi.org/10.1007/s12220-023-01262-5 doi: 10.1007/s12220-023-01262-5
    [8] B. Hu, Y. Yang, A note on the combination between local and nonlocal p-Laplacian operators, Complex Var. Elliptic Equ., 65 (2020), 1763–1776. https://doi.org/10.1080/17476933.2019.1701450 doi: 10.1080/17476933.2019.1701450
    [9] B. Barrios, I. De Bonis, M. Medina, I. Peral, Semilinear problems for the fractional laplacian with a singular nonlinearity, Open Math. J., 13 (2015), 390–407. https://doi.org/10.1515/math-2015-0038 doi: 10.1515/math-2015-0038
    [10] A. Youssfi, G. Mahmoud, Nonlocal semilinear elliptic problems with singular nonlinearity, Calc. Var. Partial Differential Equations, 60 (2021), 153. https://doi.org/10.1007/s00526-021-02034-1 doi: 10.1007/s00526-021-02034-1
    [11] B. Abdellaoui, K. Biroud, A. Primo, Nonlinear fractional elliptic problem with singular term at the boundary, Complex Var. Elliptic Equ., 64 (2019), 909–932. https://doi.org/10.1080/17476933.2018.1487410 doi: 10.1080/17476933.2018.1487410
    [12] A. Ghanmi, K. Saoudi, A multiplicity results for a singular problem involving the fractional p-Laplacian operator, Complex Var. Elliptic Equ., 61 (2016), 1199–1216. https://doi.org/10.1080/17476933.2016.1154548 doi: 10.1080/17476933.2016.1154548
    [13] J. Giacomoni, T. Mukherjee, K. Sreenadh, Positive solutions of fractional elliptic equation with critical and singular nonlinearity, Adv. Nonlinear Anal., 6 (2017), 327–354. https://doi.org/10.1515/anona-2016-0113 doi: 10.1515/anona-2016-0113
    [14] L. Boccardo, L. Orsina, Semilinear elliptic equations with singular nonlinearities, Calc. Var. Partial Differential Equations, 37 (2010), 363–380. https://doi.org/10.1007/s00526-009-0266-x doi: 10.1007/s00526-009-0266-x
    [15] J. Giacomoni, I. Schindler, P. Takac, Sobolev versus Hölder local minimizers and existence of multiple solutions for a singular quasilinear equation, Ann. Sc. Norm. Super. Pisa Cl. Sci., 6 (2007), 5393–5423. https://www.ams.org/mathscinet-getitem?mr = 2341518
    [16] A. Di Castro, Existence and regularity results for anisotropic elliptic problems, Adv. Nonlinear Stud., 9 (2009), 367–393. https://doi.org/10.1515/ans-2009-0207 doi: 10.1515/ans-2009-0207
    [17] P. Garain, On a degenerate singular elliptic problem, Math. Nachr., 295 (2022), 1354–1377. https://doi.org/10.48550/arXiv.1803.02102 doi: 10.48550/arXiv.1803.02102
    [18] P. Garain, T. Mukherjee, On a class of weighted p-Laplace equation with singular nonlinearity, Mediterr. J. Math., 17 (2022), 110. https://doi.org/10.1007/s00009-020-01548-w doi: 10.1007/s00009-020-01548-w
    [19] S. Miri, On an anisotropic problem with singular nonlinearity having variable exponent, Ric. Mat., 66 (2017), 415–424. https://doi.org/10.1007/s11587-016-0309-5 doi: 10.1007/s11587-016-0309-5
    [20] K. Bal, P. Garain, Weighted and anisotropic sobolev inequality with extremal, Manuscripta Math., 168 (2022), 101–117. https://doi.org/10.1080/17476933.2018.1487410 doi: 10.1080/17476933.2018.1487410
    [21] P. Garain, A. Ukhlov, Mixed local and nonlocal Sobolev inequalities with extremal and associated singular elliptic problems, Nonlinear Anal., 223 (2022), 113022. https://doi.org/10.1016/j.na.2022.113022 doi: 10.1016/j.na.2022.113022
    [22] C. Filippis, G. Mingione, Gradient regularity in mixed local and nonlocal problems, Math. Ann., (2022), 1–68. https://doi.org/10.1007/s00208-022-02512-7 doi: 10.1007/s00208-022-02512-7
    [23] P. Garain, J. Kinnunen, On the regularity theory for mixed local and nonlocal quasilinear elliptic equations, Trans. Amer. Math. Soc., 375 (2022), 5393–5423. https://doi.org/10.1090/tran/8621 doi: 10.1090/tran/8621
    [24] P. Garain, E. Lindgren, Higher Hölder regularity for mixed local and nonlocal degenerate elliptic equations, Calc. Var. Partial Differential Equations, 62 (2023), 67. https://doi.org/10.1007/s00526-022-02401-6 doi: 10.1007/s00526-022-02401-6
    [25] S. Huang, H. Hajaiej, Lazer-McKenna type problem involving mixed local and nonlocal elliptic operators, Res. Gate, 2023. https://doi.org/10.13140/RG.2.2.13140.68481 doi: 10.13140/RG.2.2.13140.68481
    [26] C. LaMao, S. Huang, Q. Tian, C. Huang, Regularity results of solutions to elliptic equations involving mixed local and nonlocal operators, AIMS Mathematics, 7 (2022), 4199–4210. https://doi.org/10.3934/math.2022233 doi: 10.3934/math.2022233
    [27] X. Li, S. Huang, M. Wu, C. Huang, Existence of solutions to elliptic equation with mixed local and nonlocal operators, AIMS Mathematics, 7 (2022), 13313–13324. https://doi.org/10.3934/math.2022735 doi: 10.3934/math.2022735
    [28] A. Salort, E. Veccht, On the mixed local-nonlocal Hénon equation, Differ. Integral Equ., 35 (2022), 795–818. https://doi.org/10.57262/die035-1112-795 doi: 10.57262/die035-1112-795
    [29] X. Su, E. Valdinoci, Y. Wei, J. Zhang, Regularity results for solutions of mixed local and nonlocal elliptic equations, Math. Z., 302 (2022), 1855–1878. https://doi.org/10.1007/s00209-022-03132-2 doi: 10.1007/s00209-022-03132-2
    [30] X. Zha, S. Huang, Q. Tian, Uniform boundedness results of solutions to mixed local and nonlocal elliptic operator, AIMS Mathematics, 8 (2023), 20665–20678. https://doi.org/10.3934/math.20231053 doi: 10.3934/math.20231053
    [31] E. Di Nezza, G. Palatucci, E. Valdinoci, Hitchhiker's guide to the fractional Sobolev spaces, Bull. Sci. Math., 136 (2012), 521–573. https://doi.org/10.1016/j.bulsci.2011.12.004 doi: 10.1016/j.bulsci.2011.12.004
    [32] L. Damascelli, Comparison theorems for some quasilinear degenerate elliptic operators and applications to symmetry and monotonicity results, Ann. Inst. H. Poincaré Anal. Non Linéaire, 15 (1998), 493–516. https://doi.org/10.1016/S0294-1449(98)80032-2 doi: 10.1016/S0294-1449(98)80032-2
    [33] S. Byun, K. Song, Mixed local and nonlocal equations with measure data, Calc. Var. Partial Differential Equations, 62 (2023), 14. https://doi.org/10.1007/s00526-022-02349-7 doi: 10.1007/s00526-022-02349-7
    [34] A. Moameni, K. Wong, Existence of solutions for nonlocal supercritical elliptic problems, J. Geom. Anal., 31 (2021), 164–186. https://doi.org/10.1007/s12220-019-00254-8 doi: 10.1007/s12220-019-00254-8
    [35] E. Lindgren, P. Lindqvist, Fractional eigenvalues, Calc. Var. Partial Differential Equations, 49 (2014), 795–826. https://doi.org/10.1007/s00526-013-0600-1 doi: 10.1007/s00526-013-0600-1
  • This article has been cited by:

    1. Maoji Ri, Yongkun Li, Multiple high‐energy solutions for an elliptic system with critical Hardy–Sobolev nonlinearity, 2024, 47, 0170-4214, 7684, 10.1002/mma.9997
  • 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(1742) PDF downloads(82) Cited by(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog