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,x∈RN∖Ω,
where Ω is a bounded smooth domain of RN(N>2), −Δ→pu is an anisotropic p-Laplace operator, →p=(p1,p2,...,pN) with 2≤p1≤p2≤⋅⋅⋅≤pN, (−Δ)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
[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,x∈RN∖Ω,
where Ω is a bounded smooth domain of RN(N>2), −Δ→pu is an anisotropic p-Laplace operator, →p=(p1,p2,...,pN) with 2≤p1≤p2≤⋅⋅⋅≤pN, (−Δ)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,x∈RN∖Ω, | (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)=N∑i=1∂i[|∂iu(x)|pi−2∂iu(x)], |
where →p=(p1,p2,...,pN), pi≥2 for all i=1,2,...,N. The fractional p-Laplace operator (−Δ)sp, (s∈(0,1),p≥1) is defined by
(−Δ)spu(x)=P.V.∫RN|u(x)−u(y)|p−2(u(x)−u(y))|x−y|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,x∈RN∖Ω. | (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,x∈RN∖Ω. | (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)|p−2∇u(x)],L2u(x)=N∑i=1∂i[|∂iu(x)|pi−2∂iu(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(Ω):=infu∈W1,p{0}(Ω,ω){∫Ω|∇u|pωdx:∫Ω|u|1−δf1dx=1}, |
and
v2(Ω):=infu∈W1,→p{0}(Ω){N∑i=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,x∈RN∖Ω. | (1.6) |
It has been shown that problem (1.6) has a weak solution u∈W1,p{0}(Ω) when δ∈(0,1] and f∈Lm(Ω)∖{0} with m=(p∗1−δ)′, where p∗=NpN−p, while if δ∈(1,∞) and f∈L1(Ω)∖{0}, then problem (1.6) has a weak solution u∈W1,p{loc}(Ω) with up+δ−1p∈W1,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, f∈Lm(Ω) with m≥1. Then there exists a weak solution u to problem (1.1) such that
(i) u∈L∞(Ω) if m>NˉpNˉp−(N−ˉp)pN, where ˉp satisfies
1ˉp=1NN∑i=11pi. |
(ii) u∈Lt(Ω) if ˉm<m<NˉpNˉp−(N−ˉp)pN, where
ˉm=NˉpNˉp−pi(N−ˉp)−(1−δ−pi)(N−ˉp),t=m(1−δ−pi)NˉpNˉp(m−1)−mpi(N−ˉp). |
(iii) u∈W1,q0(Ω) if 1≤m<ˉ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 f∈Lm(Ω) with NpNp−(1−δ)(N−p)<m<Np, then the solutions u to problem (1.6) satisfies u∈Lt(Ω) with
t=mN(1−δ−pi)mp−N. |
(ⅱ) If f∈Lm(Ω) with 1≤m<NpNp−(1−δ)(N−p), then the solutions u to problem (1.6) satisfies u∈W1,q{0}(Ω) with q=mN(1−δ−p)m(1−δ)−N.
Theorem 1.3. Suppose that δ=1 and 1<ˉp<N, f>0, f∈Lm(Ω) with m>1. Then there exists a weak solution u to problem (1.1) such that
(i) u∈L∞(Ω) if m>NˉpNˉp−(N−ˉp)pN.
(ii) u∈Lt(Ω) if 1≤m<NˉpNˉp−(N−ˉp)pN, where
t=mpiNˉpmpi(N−ˉp)−Nˉp(m−1). |
Theorem 1.4. Let δ>1 and 1<ˉp<N. Suppose that f>0, f∈Lm(Ω) with m>1. Then there exists a weak solution u to problem (1.1) such that
(i) u∈L∞(Ω) if m>NˉpNˉp−(N−ˉp)pN.
(ii) u∈Lt(Ω) if 1≤m<NˉpNˉp−(N−ˉp)pN, where
t=m(1−δ−pi)NˉpNˉp(m−1)−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′=pp−1 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(Ω)={u∈Lp(Ω):|u(x)−u(y)||x−y|Np+s∈Lp(Ω×Ω)}, |
with
‖u‖Ws,p(Ω)=(∫Ω|u(x)|pdx+∫Ω∫Ω|u(x)−u(y)|p|x−y|N+psdxdy)1p. |
The space Ws,p(RN) and Ws,ploc(Ω) are defined analogously. The space Ws,p0(Ω) is defined as
Ws,p0(Ω)={u∈Ws,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:E→R with the finite norm
‖u‖Lpi(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(Ω)={u∈W1,1(Ω):∂iu∈Lpi(Ω)}, |
and
W1,pi0(Ω)={u∈W1,10(Ω):∂iu∈Lpi(Ω)}, |
endowed with the norm
‖u‖W1,pi0(Ω)=N∑i=1‖∂iu‖Lpi(Ω). | (2.1) |
Definition 2.2. A function J: W1,pi0→R is defined to be weakly lower semi-continuous if
J(u)≤lim infn→∞J(un), |
for any sequence un approaching u∈W1,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 u≤0 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 u≤0 in RN∖Ω.
The definition of weak solution in this paper is defined as
Definition 2.4. A positive function u∈W1,piloc(Ω)∩Lpi−1(RN) is a weak solution to problem (1.1) if
u>0inΩ,u=0inRN∖Ω,f(x)uδ∈L1loc(Ω), |
for every ϕ∈C1c(Ω), we have that
N∑i=1∫Ω|∂iu(x)|pi−2∂iu(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)|p−2(u(x)−u(y)),dμ=|x−y|−N−psdxdy. |
Lemma 2.5. [19, Theorem 1.2] There exists a positive constant C, such that for every u∈W1,pi(Ω), we have
‖u‖pNLˉp∗(Ω)≤CN∑i=1‖∂iu‖piLpi(Ω), | (2.3) |
where
1ˉp=1NN∑i=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
⟨|ξ|pi−2ξ−|η|pi−2η,ξ−η⟩≥{cpi|ξ−η|pi,if pi≥2,cpi|ξ−η|2(|ξ|+|η|)2−pi,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)α|h−k|βforallh>k>0, |
where M>0, α>1 and β>0. Then ψ(d)=0, where dβ=Cψ(0)α−12αβ(α−1).
For n∈N, 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,x∈RN∖Ω. | (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,x∈RN∖Ω, | (2.6) |
has a unique positive weak solution u∈W1,pi0(Ω).
Proof. Existence: Define the energy functional J:W1,pi0(Ω)→R as
J(u):=J1(u)+J2(u)−J3(u), |
where
J1(u)=N∑i=11pi∫Ω|∂iu(x)|pidx, |
J2(u)=1p∫∫D(Ω)|u(x)−u(y)|p|x−y|N+psdxdy, |
and
J3(u)=∫Ωg(x)u(x)dx. |
(ⅰ) By the Sobolev embedding theorem and g∈L∞(Ω), we have
J(v)≥1pi‖v‖piW1,pi0(Ω)−|Ω|p−1p‖g‖L∞(Ω)‖v‖Lp(Ω)→∞as‖v‖piW1,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 infn→∞J(un)=lim infn→∞J1(un)+lim infn→∞J2(un)−limn→∞J3(un)≥J1(u)+J2−J3(u)=J(u). |
Combining the above properties of J, we know that there exists a minimizer u∈W1,pi0(Ω) and which is also a critical point of J, which also is the solution to Eq (2.6).
Uniqueness: Let u1, u2∈W1,pi0(Ω) be two solutions to problem (2.6). Thus, for any φ∈W1,pi0(Ω), we have
N∑i=1∫Ω|∂iu1|pi−2∂iu1∂iφdx+∫∫D(Ω)Ku1(x,y)(φ(x)−φ(y))dμ=∫Ωgφdx, | (2.7) |
and
N∑i=1∫Ω|∂iu2|pi−2∂iu2∂iφdx+∫∫D(Ω)Ku2(x,y)(φ(x)−φ(y))dμ=∫Ωgφdx. | (2.8) |
Choosing φ=u1−u2 and then subtracting (2.7) and (2.8), we obtain
N∑i=1∫Ω(|∂iu1|pi−2∂iu1−|∂iu2|pi−2∂iu2)(∂iu1−∂iu2)dx+∫∫D(Ω)(Ku1(x,y)−Ku2(x,y))[(u1−u2)(x)−(u1−u2)(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)=tp−1(p>1), we have
[Ku1(x,y)−K(u2(x,y))][(u1−u2)(x)−(u1−u2)(y)]=[|u1(x)−u1(y)|p−2(u1(x)−u1(y))−|u2(x)−u2(y)|p−2(u2(x)−u2(y))]⋅[(u1(x)−u1(x))−(u2(y)−u2(y))]≥0. | (2.10) |
Consequently,
N∑i=1∫Ω(|∂iu1|pi−2∂iu1−|∂iu2|pi−2∂iu2)(∂iu1−∂iu2)dx=0. | (2.11) |
Therefore, u1(x)−u2(x)=C for all x∈RN. Note that u1−u2=0 on RN∖Ω since ui(x)=0 for x∈RN∖Ω. Thus u1(x)≡u2(x), which implies that the solution of (2.6) is unique.
Boundedness: For any k>1, decompose RN as RN=Ak∪Ack, where
Ak={x∈Ω:u(x)≥k},Ack={x∈Ω:0<u(x)<k}. |
Taking Gk(u):=(u−k)+=max{u−k,0} as a test function in (2.6), we have
N∑i=1∫Ω|∂iu(x)|pi−2∂iu(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)|p−2(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)|p−1(u(y)−k), if u(y)>k≥u(x),|u(x)−u(y)|p−1(u(x)−k), if u(x)>k≥u(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
N∑i=1∫Ω|∂iGk(u)|pidx=N∑i=1∫Ak|∂iGk(u)|pidx+N∑i=1∫Ack|∂iGk(u)|pidx=N∑i=1∫Ak|∂iGk(u)|pidxdx≤∫Ωg(x)Gk(u)dx≤‖g‖L∞(Ω)(∫ΩGk(u)ˉp∗dx)1ˉp∗|A(k)|ˉp∗−1ˉp∗≤C‖g‖L∞(Ω)(N∑i=1∫Ω|∂iGk(u)|pidx)1pN|A(k)|ˉp∗−1ˉp∗. |
Therefore
N∑i=1∫Ω|∂iGk(u)|pidx≤C‖g‖pNpN−1L∞(Ω)|A(k)|pN(ˉp∗−1)ˉp∗(pN−1). | (2.14) |
For every 1≤k<h we know that A(h)⊂A(k) and u(x)−k≥(h−k) in A(h), we get
(h−k)pN|A(h)|pNp∗≤(∫A(h)Gk(u)ˉp∗dx)pNˉp∗≤(∫A(k)Gk(u)ˉp∗dx)pNˉp∗≤CN∑i=1∫Ω|∂iGk(u)|pidx≤C‖g‖pNpN−1L∞(Ω)|A(k)|pN(ˉp∗−1)ˉp∗(pN−1). |
Hence, we have
|A(h)|≤C‖g‖ˉp∗pN−1L∞(Ω)(h−k)ˉp∗|A(k)|ˉp∗−1pN−1. |
Obviously,
ˉp∗>pN. |
Hence, using Lemma 2.7 we obtain
‖u‖L∞(Ω)≤C. |
Positivity: First, taking u−(x):=min{u(x),0} as a test function in (2.6) and using g≥0, we have
N∑i=1∫Ω|∂iu−(x)|pidx+∫∫R2NKu(x,y)(u−(x)−u−(y))dμ=∫Ωgu−dx≤0, | (2.15) |
where Ku(x,y)=|u(x)−u(y)|p−2(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)|p−2(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)|p−2(u(x)−u(y))(u−(x)−u−(y))≥0. | (2.16) |
Using (2.16) in (2.15) we obtain
N∑i=1∫Ω|∂iu−|pidx=0. |
Therefore, u−=C for all x∈RN. Note that u−=0 on RN∖Ω since u−:=min{u,0}. Thus u≥0 in Ω.
Second, assume that there exists a point x0∈Ω such that u(x0)=infx∈Ωu(x)=0, thus
N∑i=1∫Ω∂i[|∂iu(x0)|pi−2∂iu(x0)]+∫∫D(Ω)|u(x0)−u(y)|p−2[u(x0)−u(y)]|x−y|N+psdy=∫∫D(Ω)|−u(y)|p−2[−u(y)]|x−y|N+psdy=−∫∫R2N|u(y)|p−1|x−y|N+psdy<0. |
This is a contradiction since g(x0)≥0. Hence, u>0 in Ω.
Lemma 2.9. For any n∈N, there exists a unique positive solution un∈W1,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 n∈N. By Lemma 2.8, for every u∈W1,pi0(Ω)∩L∞(Ω), there exists a unique v∈W1,pi0(Ω)∩L∞(Ω) such that
{−Δ→pv(x)+(−Δ)spv(x)=fn(x)(u+1n)δ,x∈Ω,v(x)>0,x∈Ω,v(x)=0,x∈RN∖Ω. | (2.17) |
Define the operator T:u↦v=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
N∑i=1∫Ω|∂iv(x)|pidx≤∫Ωnδ+1v(x)dx≤Cnδ+1|Ω|ˉp∗−1ˉp∗(N∑i=1∫Ω|∂iv(x)|pidx)1pN. |
Thus
N∑i=1∫Ω|∂iv(x)|pidx≤CnpN(δ+1)pN−1|Ω|pN(ˉp∗−1)ˉp∗(pN−1):=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 limk→∞‖vk−v‖W1,pi0(Ω)=0 if limk→∞‖uk−u‖W1,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)|pidx≤N∑i=1∫Ω|∂i(vk(x)−v(x))|pidx≤N∑i=1∫Ω[|∂ivk(x)|pi−2∂ivk(x)−|∂iv(x)|pi−2∂iv(x)][∂ivk(x)−∂iv(x)]dx+∫∫D(Ω)[Kvk(x,y)−Kv(x,y)][(vk−v)(x)−(vk−v)(y)]dμ=∫Ω[fn(x)(vk+1n)δ−f(x)vδ][vk(x)−v(x)]dx,ifpi≥2. | (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δ|p∗′idx]1p∗′i‖ˉvk‖Lp∗i(Ω)≤C[∫Ω|fn(x)(uk+1n)δ−f(x)uδ|p∗′idx]1p∗′i‖∂iˉvk‖Lpi(Ω). | (2.20) |
By (2.19) and (2.20), using [10, Lemma 2.2], we find
‖∂iˉvk‖Lpi(Ω)≤Cn1pi−1[∫Ω|1(uk+1n)δ−1uδ|p∗′idx]1(pi−1)p∗′i≤Cn1pi−1[∫Ω|uδ−(uk+1n)δ(uk+1n)δuδ|p∗′idx]1(pi−1)p∗′i≤Cn1pi−1+δ[∫Ω|u−uk|p∗′idx]1(pi−1)p∗′i, | (2.21) |
since the pointwise convergence of uk→u in W1,pi0(Ω). we get
limk→+∞‖vk−v‖W1,pi0(Ω)=0. |
Therefore, in the case pi≥2, the operator T is continuous from W1,pi0(Ω) to W1,pi0(Ω).
(ⅱ) Compactness of T: To achieve this, we have to show that, for some v∈W1,pi0(Ω), it holds
limk→+∞‖vk−v‖W1,pi0(Ω)=0. |
Let uk be a bounded sequence in W1,pi0(Ω) and vk:=T(uk). Then we have
uk⇀uinW1,pi0(Ω),uk→uinLt(Ω),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)⇀v∈W1,pi0(Ω),T(uk)→v∈Lt(Ω),1<t<ˉp∗. |
For any φ∈W1,pi0(Ω),
N∑i=1∫Ω|∂ivk(x)|pi−2∂ivk(x)∂iφdx+∫∫D(Ω)|vk(x)−vk(y)|p−2(vk(x)−vk(y))(φ(x)−φ(y))|x−y|N+psdxdy=∫Ωfn(x)(uk+1n)δφdx. | (2.22) |
Now, we show that as k→∞, (2.22) converges to
N∑i=1∫Ω|∂iv(x)|pi−2∂iv(x)∂iφdx+∫∫D(Ω)|v(x)−v(y)|p−2(v(x)−v(y))(φ(x)−φ(y))|x−y|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
N∑i=1∂ivk→N∑i=1∂ivpointwise almost everywhere in Ω. |
Therefore, for every φ∈C1c(Ω), we have
limk→∞N∑i=1∫Ω|∂ivk(x)|pi−2∂ivk(x)∂iφdx=N∑i=1∫Ω|∂iv(x)|pi−2∂iv(x)∂iφdx. |
Since φ∈C1c(Ω) and vk is uniformly bounded in W1,pi0(Ω),
{|vk(x)−vk(y)|p−2(vk(x)−vk(y))|x−y|N+psp′}n∈N∈Lp′(RN×RN), |
by the pointwise convergence of vk(x) to v(x)
|vk(x)−vk(y)|p−2(vk(x)−vk(y))|x−y|N+psp′→|v(x)−v(y)|p−2(v(x)−v(y))|x−y|N+psp′ a.e.inR2N. |
Since
φ(x)−φ(y)|x−y|N+psp∈Lp(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
N∑i=1∫Ω|∂iun(x)|pi−2∂iun(x)∂iφdx+∫∫D(Ω)Kun(x,y)(φ(x)−φ(y))dμ=∫Ωfn(x)(un(x)+1n)δφdx, | (2.24) |
and
N∑i=1∫Ω|∂iun+1(x)|pi−2∂iun+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
N∑i=1∫Ω|∂iun(x)|pi−2∂iun(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
N∑i=1∫Ω|∂iun+1(x)|pi−2∂iun+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))+dx≤0. | (2.28) |
Subtracting (2.26) with (2.27) and using the (2.28), we obtain
N∑i=1∫Ω[|∂iun(x)|pi−2∂iun(x)−|∂iun+1(x)|pi−2∂iun+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
N∑i=1∫Ω[|∂iun(x)|pi−2∂iun(x)−|∂iun+1(x)|pi−2∂iun+1(x)]∂i(un(x)−un+1(x))+dx≤0. |
Using Lemma 2.6 we obtain
(un(x)−un+1(x))+=Cforallx∈RN. |
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 u1≥CK>0 in K. Again, since the monotonicity of un, we have un≥u1 in K. Therefore, for any K⋐Ω,
un(x)≥CK>0,forx∈K. |
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, f∈Lm(Ω) with
m>ˉm=NˉpNˉp−pi(N−ˉp)−(1−δ−pi)(N−ˉp). |
Then, the sequence solutions un to the approximate problem (2.5) such that
(i) un∈L∞(Ω) if m>NˉpNˉp−pN(N−ˉp).
(ii) un∈Lt(Ω), where t=m(1−δ−pi)NˉpNˉp(m−1)−pim(N−ˉp) if
NˉpNˉp−pi(N−ˉp)−(1−δ−pi)(N−ˉp)<m<NˉpNˉp−pN(N−ˉp). |
Proof. (i) Let Ak={x∈Ω:un(x)≥k}. Choosing Gk(u):=(u−k)+∈W1,pi0(Ω) as a test function in (2.5), we get
N∑i=1∫Ω|∂iun(x)|pi−2∂iun(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
N∑i=1∫Ω|∂iGk(un)|pidx=N∑i=1∫Ak|∂iGk(un)|pidx≤∫Ωfn(x)(un(x)+1n)δGk(un)dx≤∫Akf(x)G1−δk(un)dx≤(∫Akf(x)mdx)1m(∫ΩGk(un)ˉp∗dx)1−δˉp∗|A(k)|1−1m−1−δˉp∗≤C(∫Akf(x)mdx)1m(N∑i=1∫Ω|∂iGk(un)|pidx)1−δpN|A(k)|1−1m−1−δˉp∗. | (3.2) |
Hence
N∑i=1∫Ω|∂iGk(un)|pidx≤C(∫Akf(x)mdx)pNm(pN+δ−1)|A(k)|(1−1m−1−δˉp∗)pNpN+δ−1. |
Let h>k≥1, we know that Ah⊂Ak and Gk(un)≥h−k for in Ω, we have that
|h−k|pN|Ah|pNˉp∗≤(∫A(h)Gk(un)ˉp∗dx)pNˉp∗≤(∫A(k)Gk(un)ˉp∗dx)pNˉp∗≤CN∑i=1∫A(k)|∂iGk(un)|pidx≤C‖f‖pNpN+δ−1Lm(Ω)|A(k)|(1−1m−1−δˉp∗)pNpN+δ−1. |
Therefore
|Ah|≤C‖f‖ˉp∗pN+δ−1Lm(Ω)|Ak|(1−1m−1−δˉp∗)ˉp∗pN+δ−1|h−k|ˉp∗. |
Note that
(1−1m−1−δˉp∗)ˉp∗pN+δ−1>1, |
if m>NˉpNˉp−pN(N−ˉp). Hence, apply Lemma 2.7 with
M=C‖f‖ˉp∗pN+δ−1Lm(Ω)>0,α=(1−1m−1−δˉp∗)ˉp∗pN+δ−1>1,β=ˉp∗>0andψ(k)=|Ak|, |
there exists k0 such that ψ(k)≡0 for all k≥k0. Thus,
esssupΩu≤k0. |
(ⅱ) Choose upi(γ−1)+1n (γ>1) as test function in (2.17), we have
N∑i=1∫Ω|∂iun|pi−2∂iun∂iupi(γ−1)+1ndx+∫∫D(Ω)Kun(x,y)[un(x)pi(γ−1)+1−un(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)+1−un(y)pi(γ−1)+1]=|un(x)−un(y)|p−2(un(x)−un(y))[un(x)pi(γ−1)+1−un(y)pi(γ−1)+1]≥C[un(x)+un(y)]pi(γ−1)|un(x)−un(y)|p≥0. | (3.4) |
Combining (3.4) and (3.3), and using Hölder inequality, we get
N∑i=1∫Ω|∂iun|pi−2∂iun∂iupi(γ−1)+1ndx=N∑i=1[pi(γ−1)+1]∫Ω|∂iun|piupi(γ−1)ndx≤∫Ωf(x)(un+1n)δupi(γ−1)+1ndx≤∫Ωf(x)upi(γ−1)+1−δndx≤‖f‖Lm(Ω)(∫Ωu[pi(γ−1)+1−δ]m′ndx)1m′. |
By Sobolev inequality,
N∑i=1∫Ω[pi(γ−1)+1]|∂iun|piupi(γ−1)ndx=N∑i=1∫Ω[pi(γ−1)+1](1γ)pi|∂iuγn|pidx≥C(∫Ωuγˉp∗n)pNˉp∗. |
Therefore,
(∫Ωuγˉp∗n)pNˉp∗≤C‖f‖Lm(Ω)(∫Ω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(m−1)−pim(N−ˉp). |
Since γ>1, we know
NˉpNˉp−pi(N−ˉp)−(1−δ−pi)(N−ˉp)<m. |
Thus pNˉp∗>1m′ gives
(∫Ωun(x)γˉp∗)pNˉp∗−1m′≤C‖f‖Lm(Ω). |
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 f∈Lm(Ω) with
1≤m<NˉpNˉp−pi(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 δ+pi−1pi≤γ<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
N∑i=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
N∑i=1[p1(γ−1)+1]∫Ω|∂iun|pi(un+ε)pi(γ−1)dx≤∫Ωf(x)(un+ε)pi(γ−1)+1−δdx. | (3.6) |
By Sobolev inequality,
N∑i=1[pi(γ−1)+1]∫Ω|∂iun|pi(un+ε)pi(γ−1)dx=N∑i=1[pi(γ−1)+1]∫Ω(1γ)pi|∂i[(un+ε)γ−εγ]|pidx≥CN∑i=1∫Ω|∂i[(un+ε)γ−εγ]|pidx≥C(∫Ω[(un+ε)γ−εγ]ˉp∗dx)pNˉp∗. | (3.7) |
Hence, by (3.6) and (3.7), we get
(∫Ω[(un+ε)γ−εγ]ˉp∗dx)pNˉp∗≤C∫Ωf(x)(un+ε)pi(γ−1)+1−δdx. |
Let ε→0, we have
(∫Ωuγˉp∗ndx)ˉpˉp∗≤C∫Ωf(x)upi(γ−1)+1−δndx. | (3.8) |
If m=1, we choose γ=pi+δ−1pi in the (3.8), so that un∈L(δ+pi−1)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(m−1)−pim(N−ˉp). |
Since γ<1, by (3.6), we have
N∑i=1∫Ω|∂iun|pi(un+ε)pi−piγdx=N∑i=1∫Ω|∂iun|pi(un+ε)piγ−pidx≤C. |
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−γ)qpipi−qdx)1−q¯pi≤C(∫Ω(un+ε)(1−γ)piqpi−qdx)1−q¯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
![]() |
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 |