MALAYSIANMATHEMATICAL
SCIENCESSOCIETY http://math.usm.my/bulletin
Sharp Remainder Terms of the Rellich Inequality and its Application
1ALNARDETALLA,2TOSHIOHORIUCHI AND3HIROSHIANDO
1Department of Mathematics, Central Mindanao University, University Town, 8710 Musuan, Bukidnon, Philippines
2,3Department of Mathematical Sciences, Ibaraki University Mito, Ibaraki, 310, Japan
1al [email protected],2[email protected]
Abstract. In this article we shall study the improvement of the Rellich inequality by adding terms with a singular weight of the type(log(1/|x|))−2in the right hand side. We show that this weight function is optimal in the sense that the improved inequality fails for any other weight more singular than this one. As an application, we use our improved inequality to an- alyze the behaviour of the first eigenvalue of the weighted eigenvalue problem for the opera- torLµu=∆ |∆u|p−2∆u
− µ/|x|2p
|u|p−2uasµincreases to((n−2p)/p)p((np−n)/p)p for 1<p<n/2.
2010 Mathematics Subject Classification: Primary: 35J70; Secondary: 35J60 Keywords and phrases: Rellich inequality, eigenvalue.
1. Introduction
In this paper we shall study the Rellich inequality (1.1)
Z
Ω
|∆u|pdx≥
n−2p p
p np−n
p pZ
Ω
|u(x)|p
|x|2p dx
for anyu∈W02,p(Ω), whereΩis a bounded domain inRnwith 0∈Ω,n≥3, and 1<p<
n/2. Here the best constantΛn,p= ((n−2p)/p)p((np−n)/p)pis given by the infimum of I(u) =
R
Ω|∆u|pdx R
Ω
|u(x)|p
|x|2p dx .
Moreover there exists no extremal function inW02,p(Ω)which attains the infimum of these problem. Roughly speaking, the candidates of the extremals are singular at the origin which are not in the classW02,p(Ω). Hence it is natural to consider that there exist “missing terms” in the right hand side of (1.1). In view of this, we shall investigate the Rellich inequalities (1.1) and improve them by finding out missing terms.
Received:June 22, 2009;Revised:May 31, 2010.
For the case of gradient, such improved Hardy inequalities are known. For example, Adimurthi, Chaudhuri and Ramaswamy [1] have proved that there exists a constantC>0, depending onn≥2, 1<p<nandR>supΩ(|x|e2/p)such that foru∈W01,p(Ω)
(1.2) Z
Ω
|∇u|pdx≥ n−p
p pZ
Ω
|u(x)|p
|x|p dx+C Z
Ω
|u(x)|p
|x|p
log R
|x|
−γ
dx
Our aim in this article is to achieve an optimal improvement of the inequality (1.1) by adding a second term involving the singular weight(log(1/|x|))−2, in the sense that the improved inequality holds for this weight but fails if the weight is more singular than this one.
As an application, we use our improved inequality to determine exactly when the first eigenvalue of the weighted eigenvalue problem for the operator
(1.3) Lµu=∆ |∆u|p−2∆u
− µ
|x|2p|u|p−2u will tend to 0 asµincreases toΛn,p.
This paper is organized in the following way. In§2 we shall describe our main results on the improvement of Rellich inequality. In§3 we shall prepare lemmas which are needed in the proof of the main theorem stated in§2. In§4 we shall prove the main results (Theorem 2.1 and Corollary 2.1). In§5 we shall apply our results to study the weighted eigenvalue problem.
2. Main results
Theorem 2.1. Let n≥3,0∈ΩandΩis a bounded domain inRn. (1) Noncritical case(1<p<n/2)
Assumeγ≥2, then there exists K=K(n)>0and C=C(n)>0such that if R>
KsupΩ|x|then Z
Ω
|∆u|pdx≥
n−2p p
p np−n
p pZ
Ω
|u(x)|p
|x|2p dx +C
Z
Ω
|u(x)|p
|x|2p
log R
|x|
−γ
dx (2.1)
for any u∈W02,p(Ω).
(2) Critical case(p=n/2)
Assume thatγ≥n/2. Then there exists K∗=K∗(n)>0and C∗=C∗(n)>0such that if R>K∗supΩ|x|, then
Z
Ω
|∆u|n2dx≥ n−2
√n nZ
Ω
|u(x)|n2
|x|n
log R
|x|
−γ
dx
+C∗ Z
Ω
|u(x)|n2
|x|n
log R
|x|
−γ−1
dx (2.2)
for any u∈W02,n/2(Ω).
Remark 2.1. In (2.1),γ≥2 is sharp. In (2.2),γ≥n/2 is also sharp, and((n−2)/√ n)nis best constant.
Remark 2.2. The functionγ→(log(R/r))−γis monotonically decreasing on[2,∞) ([n/2,∞)) provided thatR>supΩ|x|. Hence it suffices to assumeγ=2(γ=n/2)in noncritical case (critical case).
Remark 2.3. In the proof of the noncritical case, we will use decreasing rearrangement argument. Hence, the functiong(r) =r−2p(log(R/r))−2should be monotone decreasing andR≥re1/p. ThenK=e1/p.
Remark 2.4. In the proof of the critical case, we will also use decreasing rearrangement argument. Hence, the function g∗(r) =r−n(log(R/r))−(n/2)−1 should also be monotone decreasing andR≥re(1/2)+(1/n). Moreover we need the condition to absorb the error terms in the right hand side of (2.2) withC∗>0. HenceK∗≥e(1/2)+(1/n).
Remark 2.5. CandC∗may depend onRandγin a weak sense. Sinceg(r)andg∗(r)tends to zero asγ→∞orR→∞, we can takeCandC∗to be bigger.
Corollary 2.1. Let 1<p<n/2, and let Fp=
(
f:Ω→R+
f ∈L∞loc(Ω\ {0})with lim sup
|x|→0
|x|2pf(x)
log 1
|x|
2
<∞ )
.
If f ∈Fp, then there existsλ(f)>0 such that foru∈W02,p(Ω) (2.3)
Z
Ω
|∆u|pdx≥Λn,p Z
Ω
|u(x)|p
|x|2p dx+λ(f) Z
Ω
|u(x)|pf(x)dx
Iff ∈/Fpand if|x|2pf(x) (log(1/|x|))2tends to∞as|x| →0, then no inequality of type (2.3) can hold.
3. Preliminary lemmas
In this section we shall prepare fundamental lemmas which are needed to prove our main results.
Lemma 3.1. For any R>1, q≤0,ν∈(0,1)satisfying2ν−1+q=0, (3.1)
Z 1 0
|h0(r)|2
logR r
q
r dr≥ν2 Z 1
0
|h(r)|2
logR r
q−2
dr r holds for any h∈C([0,1])∩C1(0,1), with h(0) =h(1) =0.
Proof. Leth(r) = (log(R/r))νw(r). Thenw(0) =w(1) =0.
|h0(r)|2=
logR r
(ν−1)2
−ν
rw(r) +w0(r)logR r
2
≥ν r
2
|h(r)|2
logR r
−2
−ν r
d
drw2(r) logR r
2ν−1
. Hence for 2ν−1+q=0, we have
Z 1 0
|h0(r)|2
logR r
q
r dr≥ν2 Z 1
0
|h(r)|2
logR r
q−2
dr r .
In the next lemma, we assume−∆u>0 andu>0 since we will use symmetric decreasing rearrangement in the proof of the main theorem. See also Lemma 3.4.
Lemma 3.2. Assume that f ∈C2(B1) and u∈C02(B1) are radial satisfying f(r)>0,
∆f(r)≤0, u(r)>0, and−∆u>0, where r=|x|. Set u(r) =f(r)v(r). Then Z
B1
|∆u|n2dx≥ωn Z 1
0
rn−1|v(r)∆f(r)|n2dr
+n(n−2) 4 ωn
Z 1 0
(v0(r))2vn−42 (r)rn−1|∆f(r)|n−22 f(r)dr
+ωn Z 1
0
vn2(r)∂r rn−1
|∆f(r)|n−22 f0(r)−∂r|∆f(r)|n−22 f(r) dr (3.2)
Proof. Note that∆u= f∆v+2v0f0+v∆f <0 andv(r)∆f(r)≤0. By the inequality(1+ x)p≥1+px(x>−1), we have
|∆u|n2 =|v(r)∆f(r)|n2
1+ f(r)∆v(r) +2f0(r)v0(r) v(r)∆f(r)
n 2
≥ |v(r)∆f(r)|n2
1+n 2
f(r)∆v(r) +2f0(r)v0(r) v(r)∆f(r)
and Z
B1
|∆u|n2dx≥ωn Z 1
0
rn−1|v∆f|n2dr
−n 2ωn
Z 1 0
rn−1|v(r)∆f(r)|n−22 f(r)∆v(r) +2f0(r)v0(r) dr.
Then by using the integration by parts, we get the desired result.
The next lemma is the critical case of Theorem 2.1, whereΩ=B1anduis radial.
Lemma 3.3. Consider positive radially nonincreasing function u∈C20(B1). Then there exist K∗=K∗(n)>0and C∗=C∗(n)>0such that if R>K∗, then
Z
B1
|∆u|n2dx≥ n−2
√n nZ
B1
|u(x)|n2
|x|n
log R
|x|
−n
2
dx
+C∗ Z
B1
|u(x)|n2
|x|n
log R
|x|
−n2−1
dx.
(3.3)
A sketch of the proof of Lemma 3.3. Foru(r) =f(r)v(r), we setf(r) = (logR/r)a, 0<a<
1. Then we use Lemma 3.1 and Lemma 3.2 to prove inequality (3.3).
We recall the rearrangement of domains and functions. For a domainΩ, we define a ballΩ∗such that|Ω∗|=|Ω|with center at the origin. We denote by|u|∗ the symmetric decreasing rearrangement of functionu. It is well-known that the symmetric rearrangement does not change theLp-norm and increases the integralRΩ(|u|p)/(|x|2p)dx. The following lemma is due to G. Talenti. As for the proof refer to [5] and [7].
Lemma 3.4. [7]LetΩbe a domain onRn, n≥3and f∈C∞0(Ω). If u is the weak solution of Dirichlet problem−∆u= f inΩ, u|∂Ω=0; v is the weak solution of Dirichlet problem
−∆v=|f|∗inΩ∗, v|∂Ω∗ =0; then v≥ |u|∗in pointwise.
From this lemma we see thatu∗≤vinΩ∗and Z
Ω
|∆u|pdx= Z
Ω
|f|2dx= Z
Ω∗
|f|∗2dx Z
Ω∗
|∆v|pdx.
Further we see that Z
Ω
|u|p
|x|2pdx≤ Z
Ω∗
|u∗|p
|x|2pdx≤ Z
Ω∗
|v|p
|x|2pdx.
Hence we have
R
Ω|∆u|pdx R
Ω
|u|p
|x|2pdx
≥ R
Ω∗|∆v|pdx R
Ω∗
|v|p
|x|2pdx .
Using this lemma we can assume thatuis radial andΩis a ball in the proof of the main results.
4. Proof of main results
We are now ready to give the proof of Theorem 2.1. We organize the proof in the following way: First for noncritical case (1<p<n/2), we prove inequality (2.1). Then we show the sharpness ofΛn,p, and then we show the optimality ofγ. Secondly for critical case (p=n2), we first prove Inequality (2.2). Then we show the sharpness of((n−2)/√
n)n, and then we show the optimality ofγ.
Proof of Theorem 2.1. Case 1 (1<p<n/2): Let 1<p<n/2 andγ≥2. It suffices to prove (2.1) for smooth positive radially nonincreasing functionu defined on a ballB1and from Remark 3.2 we can further assumeγ=2. Foru∈C0∞(B1),u>0, radially nonincreasing, we define
(4.1) v(r) =u(r)r(n/p)−2, r=|x|.
From Lemma 3.4, we may assume−∆u>0 and
∆(u(r)) =∆
v(r)r2−np
=r2−np
δβ(v(r))−αv(r) r2
where
δβ(v(r)) =∆β(v(r)) =v”(r) +β−1 r v0(r), β=n+4−2n
p and α=(n−2p)(np−n)
p2 .
Then by using(1+x)p≥1+px (x>−1) and Lemma 3.1 (ν=1/2,q=0), we get Inequality (2.1) whereC= ((Λn,p)/α)((p−1)/p).
We construct a family of functions inW02,p(B1)which we will use to show the sharpness ofΛn,pand optimality ofγ. Forε>0 sufficiently small, let us define
(4.2) uε=
0, 0<r<ε2
r1−np log1
ε
−1 log r
ε2, ε2<r<ε r1−np log1
ε
−1
log1r, ε<r<1
Letwε=Rr1uε(ρ)dρ. Then we get (4.3)
Z
B1
|∆wε|pdx= 2 p+1
n(p−1) p
p
ωnlog1 ε+O
log1
ε −1!
,
(4.4)
Z
B1
|wε|p
|x|2p dx≥ 2 p+1
p n−2p
p
ωnlog1 ε+O
log1
ε −1!
.
Also we get (4.5)
Z
B1
|wε|p
|x|2p
log R
|x|
−γ
≥O
log1 ε
1−γ! .
SHARPNESS OFΛn,p.The sharpness ofΛn,pwill follow if we can show that
(4.6) inf
u∈W02,p(B1)\{0}
I(u):=
R
B1|∆u|pdx R
B1
|u|p
|x|2pdx
=Λn,p
Using the family of functionswε∈W02,p(B1)and from (4.3) and (4.4), we get
ε→0limI(wε)≤Λn,p
Also by Rellich inequality, we get lim
ε→0I(wε)≥Λn,p, hence lim
ε→0I(wε) =Λn,p. Thus sharp- ness follow.
OPTIMALITY OFγ. Suppose that 0≤γ<2. SinceΛn,pis the best constant for inequality (1.1), inequality (2.1) follows for the caseγ=0. So we assume 0<γ<2. Optimality will follow if we can prove that
(4.7) inf
u∈W02,p(B1)\{0}
Iγ(u) = R
B1|∆u|pdx−Λn,pR B1
|u|p
|x|2pdx R
B1
|u|p
|x|2p
log|x|R−γ
dx
=0.
Since 0<γ<2, from (4.3), (4.4) and (4.5),Iγ(wε)→0 asε→0 and hence optimality follows.
Case p=n/2: Let us assume p=n/2 andγ ≥n/2. From Remark 2.2, we can further assumeγ=n/2. In this case, inequality (2.2) follows immediately from Lemma 3.3 together with the rearrangement of function and domain arguments.
SHARPNESS. To show the sharpness of((n−2)/√
n)n, we consider the test function zε=
log R
r+ε n−2n
−
log R 1+ε
n−2n .
Then it is easy to verify by similar calculation as in the previous case(1<p<n/2)that
ε→0lim R
B1|∆zε|n2dx R
B1
|zε|n2
|x|n
log|x|R−γ
dx
= n−2
√n n
.
Hence sharpness follow.
OPTIMALITY.To show the optimality, we use the same test functionuε defined in (4.2) withp=n/2. Then forwε=Rr1uε(ρ)dρ, it is easy to verify by similar calculation as in the previous case(1<p<n/2)that
(4.8)
Z
B1
|∆wε|n2dx= 4
n+2(n−2)n2ωnlog1 ε+O
log1
ε −1!
and (4.9)
Z
B1
|wε|n2
|x|n
log R
|x|
−γ
dx≥O
log1 ε
n2+1−γ!
Suppose that 0<γ<n/2. Optimality will follow if we can prove that inf
u∈W2,
n 0 2(B1)\{0}
Iγ(u) = R
B1|∆u|n2dx R
B1
|u|n2
|x|n
log|x|R−γ
dx
=0
Since 0<γ<n/2, from (4.8) and (4.9),Iγ(wε)→0 asε→0 and hence optimality follows.
Proof of Corollary 2.1. If f ∈Fp, then
ε→0limsup
x∈Bε
f(x)|x|2p
log 1
|x|
2
<∞
and hence for sufficiently smallε, inBε, f(x)<C|x|−2p
log 1
|x|
−2
.
OutsideBε, both are bounded functions and henceCcan be chosen so that this inequality holds inΩ. Then (2.3) will follow from (2.1).
If f ∈/Fpand if|x|2pf(x) (log(1/|x|))2tends to∞as|x| →0, then we can write f(x) = h(x)/
|x|2p
log|x|12
, whereh(x)tends to infinity asxtends to 0. Then from the calcu- lation of Theorem (2.1) forε>0 sufficiently small, we get
(4.10)
Z
B1
|wε|ph(x)
|x|2p
log|x|12dx≥O
log1 ε
−1! mε
wheremε =min{infB
ε2h(x),infBε\B
ε2h(x),infB1\Bεh(x)}. Sincemε tends to∞asε→0, we conclude thatIf(wε)→0 asε→0 and Inequality (2.3) cannot hold for such f∈/Fp. 5. Application
Consider the weighted eigenvalue problem with a singular weight
∆ |∆u|p−2∆u
− µ
|x|2p|u|p−2u=λ|u|p−2u f inΩ u=∆u=0 on∂Ω (5.1)
wheref ∈Fpand Fp:=
f :Ω→R+
lim
|x|→0|x|2pf(x) =0 with f∈L∞loc Ω¯\ {0}
,
1<p<n/2, 0≤µ<Λn,pandλ ∈R. We look for a weak solutionu∈W =W2,p(Ω)∩ W01,p(Ω)of this problem and study the asymptotic behaviour of the first eigenvalues for different singular weights as µincreases to Λn,p, after which the operatorLµ is no more bounded from below. Here we define weak solution in the following way.
Definition 5.1. u∈W is said to be a weak solution of (5.1) iff for anyφ∈C2(Ω)¯ withφ=0 on∂Ω
Z
Ω
|∆u|p−2∆u∆φ− µ
|x|2p|u|p−2uφ
dx=λ Z
Ω
|u|p−2u fφdx.
Lemma 5.1. Forµ∗<Λn,pand u∈W , there is a v∈W such that v>0and satisfies R
Ω|∆u|pdx−µ∗RΩ |u|
p
|x|2pdx R
Ω|u|pf dx ≥ R
Ω|∆v|pdx−µ∗RΩ |v|
p
|x|2pdx R
Ω|v|pf dx .
Remark 5.1. Sinceλ is first eigenvalue anduis the corresponding eigenfunction, by using Lemma 5.1, we can assumeu>0 inΩ. Then by the elliptic regularity theory,uis smooth near the boundary. From the definition of weak solution, one can derive the boundary condition of (5.1) by using integration by parts.
Theorem 5.1. The above problem admits a positive weak solution u∈W for all1<p<n/2 and0≤µ<Λn,p, corresponding to the first eigenvalueλ =λµ1(f)>0. Moreover, asµ increases toΛn,p,λµ1(f)→λ(f)≥0for all f ∈Fpand the limitλ(f)>0if f ∈Fp. If
f ∈/Fpand if|x|2pf(x) (log(1/|x|))2tends to∞as|x| →0, then the limitλ(f) =0.
In order to prove the theorem we need the following results; the first one is a standard result from measure theory (see [4], Chapter 1, section 4) and the second one is due to Boccardo and Murat [4].
Lemma 5.2. Let(gm)m∈N⊂Lp(Ω),1≤p<∞, be such that, as m→∞, (i) gm*g weakly in Lp(Ω)and (ii) gm(x)→g(x)a.e. inΩ. Then
m→∞lim kgmkpp− kgm−gkpp
=kgkpp.
Lemma 5.3(see Remark 2.7 in [2]).
|∆um|p−2∆um− |∆u|p−2∆u
∆(um−u)≥
(s|∆(um−u)|p if p≥2;
s|∆(um−u)|2
(|∆um|+|∆u|)2−p if p≤2, for some s>0.
Remark 5.2. In the proof of Theorem 5.1,uwill be characterized as a solution of variational problem defined in (5.2) and (5.1) becomes Euler-Lagrange equation of this variational problem.
Proof of Theorem 5.1. We define the functional
(5.2) Jµ(u) =
Z
Ω
|∆u|p−µ|u|p
|x|2p
dx
By the Hardy inequality, one can check thatJµ is continuous, Gateaux differentiable and coercive onW. We minimize this functional over the manifold
M=
u∈W Z
Ω
|u(x)|pf(x)dx=1
,
and letλµ1>0 be the infimum. By standard arguments, we can choose (um)m∈N⊂M a minimizing sequence ofJµ such thatJµ(um)→λµ1andJµ0(um)→0 strongly inW0. The coercivity ofJµ implies that(um)m∈Nis a bounded sequence and hence we have for a sub- sequence, ask→∞,umk*uweakly inW. Moreover, we can also assume the following:
∆umk*∆uweakly inLp(Ω), umk*uweakly inLp(Ω,|x|−2p), umk→ustrongly inLp(Ω,f), umk→ustrongly inLp(Ω),
∇umk→∇ustrongly inLp(Ω).
(5.3)
SinceW is compactly embedded inLp(Ω,f), it follows thatMis weakly closed and hence u∈M. Alsoumksatisfies inD0(Ω)that
∆ |∆umk|p−2∆umk
= µ
|x|2p|umk|p−2umk+λmk|umk|p−2umkf+fmk
where fmk→0 strongly inW0andλmk→λ asmk→∞. Letφ∈C∞(Ω)such that for any ballBρ,Bρ∩suppφ=/0 and
(5.4) φ=
(0, |x|<ρ; 1, |x|>2ρ.
Using the test functionvmk=φ(umk−u)∈W, we obtain Z
Ω
|∆umk|p−2∆umk− |∆u|p−2∆u
∆vmkdx
=− Z
Ω
|∆u|p−2∆u∆vmkdx+µ Z
Ω
|umk|p−2umk
|x|2p vmkdx +λmk
Z
Ω
|umk|p−2umkf vmkdx+ Z
Ω
fmkvmkdx
From (5.3) one deduces that vmk →0 weakly in W and Lp(Ω,|x|−2p), and strongly in Lp(Ω,f)andLp(Ω). Hence formk→∞and by expanding∆vmk, we get
(5.5) lim
mk→∞
Z
Ω
emkφ=0, whereemk= |∆umk|p−2∆umk− |∆u|p−2∆u
∆(umk−u).
Fixθwith 0<θ<1. Then Z
Ω
eθm
k=
Z
Ω
eθm
kφ+ Z
Ω
eθm
k(1−φ)
≤ Z
Ω
eθm
kφ+ Z
Ω
emk θZ
Ω
(1−φ)1−θ1 1−θ
. From (5.4) and (5.5), we get
mlimk→∞
Z
Ω
eθm
k≤Cθρn(1−θ). which by lettingρtends to 0 implies that for any subsequencee0m
kofemk,
(5.6) e0m
k→0 a.e. inΩ We apply Lemma 5.3 on (5.5) and get
(5.7) ∆um0
k(x)→∆u(x) a.e.x∈Ω.
We apply Lemma 5.2 onumk and also on∆umkto obtain λµ1=k∆(umk−u)kpp−µkumk−ukp
Lp(Ω,|x|−2p)+k∆ukpp−µkukp
Lp(Ω,|x|−2p)+o(1)
≥(Λn,p−µ)kumk−ukp
Lp(Ω,|x|−2p)+λµ1+o(1),
whereo(1)→0 asmk→∞. Asµ<Λn,p,we conclude thatkumk−ukp
Lp(Ω,|x|−2p)→0 as mk→∞and alsok∆(umk−u)kLpp(Ω)→0 asmk→∞and hence we haveJµ(u) =λµ1 and λ=λµ1. Sinceuis the first eigenfunction, it is not difficult to show by using Lemma 5.1 that u>0 inΩ. Then by the Euler-Lagrange equation,uis a distribution solution of (5.1) and sinceu∈W, it is a weak solution to the eigenvalue problem (5.1), corresponding toλ=λµ1. The remaining part of the proof follows from the corollary of the main theorem.
References
[1] Adimurthi, N. Chaudhuri and M. Ramaswamy, An improved Hardy-Sobolev inequality and its application, Proc. Amer. Math. Soc.130(2002), no. 2, 489–505 (electronic).
[2] A. Detalla, T. Horiuchi and H. Ando, Missing terms in Hardy-Sobolev inequalities,Proc. Japan Acad. Ser. A Math. Sci.80(2004), no. 8, 160–165.
[3] A. Detalla, T. Horiuchi and H. Ando, Missing terms in Hardy-Sobolev inequalities and its application,Far East J. Math. Sci. (FJMS)14(2004), no. 3, 333–359.
[4] L. Boccardo and F. Murat, Almost everywhere convergence of the gradients of solutions to elliptic and para- bolic equations,Nonlinear Anal.19(1992), no. 6, 581–597.
[5] T. Horiuchi, Missing terms in generalized Hardy’s inequalities and its applications,J. Math. Kyoto Univ.43 (2003), no. 2, 235–260.
[6] M. Struwe,Variational Methods, second edition, Ergebnisse der Mathematik und ihrer Grenzgebiete (3), 34, Springer, Berlin, 1996.
[7] G. Talenti, Nonlinear elliptic equations, rearrangements of functions and Orlicz spaces,Ann. Mat. Pura Appl.
(4)120(1979), 160–184.