Skip to content
BY 4.0 license Open Access Published by De Gruyter July 17, 2020

Positive Solutions for Resonant (p, q)-equations with convection

  • Zhenhai Liu EMAIL logo and Nikolaos S. Papageorgiou

Abstract

We consider a nonlinear parametric Dirichlet problem driven by the (p, q)-Laplacian (double phase problem) with a reaction exhibiting the competing effects of three different terms. A parametric one consisting of the sum of a singular term and of a drift term (convection) and of a nonparametric perturbation which is resonant. Using the frozen variable method and eventually a fixed point argument based on an iterative asymptotic process, we show that the problem has a positive smooth solution.

MSC 2010: 35J60; 35J91; 35J92; 35D30; 35D35

1 Introduction

Let Ω ⊆ ℝN be a bounded domain with a C2-boundary Ω. In this paper we study the following parametric singular double phase Dirichlet problem with gradient dependence (convection)

Δpu(z)Δqu(z)=λ[u(z)η+r(z)|Du(z)|τ1]+f(z,u(z))inΩ,u|Ω=0,u>0,λ>0,1<q,τ<p,0<η<1. (Eλ)

For every r ∈ (0, ∞) by Δr we denote the r-Laplace differential operator defined by

Δru=div(|Du|r2Du)uW01,r(Ω).

In problem (Eλ) we have the sum of two such operators with different exponents. Hence the differential operator of the problem is not homogeneous. In the reaction (right hand side of (Eλ)), we have the combined effects of three terms, each of different nature. There is a parametric contribution which is the sum of a singular term and of a gradient dependent term (a drift term). Both are multiplied with the parameter λ > 0. So, we have dependence on the gradient of u (convection). The third term is the perturbation f(z, x) which is a Carathéodory function (that is, for all x ∈ ℝ zf(z, x) is measurable and for a.a zΩ xf(z, x) is continuous). We assume that f(z, ⋅) exhibits (p-1)-linear growth as x → +∞ and it is resonant with respect to a nonprincipal variational eigenvalue of (Δp, W01,p (Ω)). We are looking for positive solutions of the problem.

The presence of the drift term makes the problem nonvariational and so ultimately our method of proof will be topological. Our approach uses the so called “frozen variable method”. More specifically, we fix (freeze) the gradient (drift) term and this way we have a variational problem. However, the presence of the singular term is a source of difficulties, since the energy (Euler) functional of this variational problem is not C1 and so the results and methods of critical point theory can not be used directly on this function. We need to find a way to bypass the singularity and work with a C1-functional. This is done by introducing an auxiliary problem which we solve and then use its solution and suitable truncation techniques to neutralize the singularity. We are able to show that for all small parameter values, the “frozen problem” has a positive solution. Next, we need to find a canonical way to choose a solution of the “frozen problem”. To this end, we show that each such problem has a minimal positive solution (a smallest positive solution). We choose this solution and we have the minimal positive solution map. We show that this map has a fixed point using the Leray-Schauder Alternative Principle (see Theorem 4). To show that the minimal solution map is compact (requirement in the Leray-Schauder theorem), we employ an iterative asymptotic process.

Singular problems with convection, were studied recently by Bai-Gasiński-Papageorgiou [1], Liu-Motreanu-Zeng [2] and Papageorgiou-Rădulescu-Repovš [3]. These works deal with equations driven by the p-Laplacian (hence the differential operator of the problem is homogeneous and this property is exploited in their proofs) and their hypotheses require that the perturbation term asymptotically as x → +∞ stays below λ͡1(p), the principal eigenvalue of (−Δp, W01,p (Ω)) and no interaction is allowed (uniform nonresonance). This way the “frozen problem” becomes coercive and the direct method of the calculus of variations can be used to produce a solution, In contrast here we have resonance with respect to a nonprincipal variational eigenvalue and so the frozen problem has an indefinite energy functional. Finally, we mention that resonant singular problems driven by the p-Laplacian only and with no convection, were investigated by Papageorgiou-Rădulescu-Repovš [4] and by Papageorgiou-Vetro-Vetro [5], while resonant (p, q)-equations were studied by Papageorgiou-Zhang [6]. Additional relevant results can be found in the works of Gasiński-Winkert [7]. Marano-Winkert [8], Papageorgiou-Rădulescu-Repovš [9, 10, 11, 12], Papageorgiou-Vetro-Vetro [13, 14], Ragusa-Tachikawa [15], Tang-Lin-Yu [16], Vetro [17, 18, 19].

2 Mathematical Background, Hypotheses

The main spaces in the study of problem (Eλ) are the Sobolev space W01,p (Ω) and the Banach space C01 (Ω) = {uC1(Ω) : uΩ = 0}. By ∥ ⋅ ∥, we denote the norm of the Sobolev space W01,p (Ω). From the Poincaré inequality, we have

u=Dupfor all uW01,p(Ω).

The Banach space C01 (Ω) is ordered with positive (order) cone C+ = {u C01 (Ω) : u(z) ≥ 0 for all zΩ}. This cone has a nonempty interior given by

intC+={uC+:u(z)>0for all zΩ,un|Ω<0}

with n(⋅) being the outward unit normal on Ω.

Let r ∈ (1, ∞). By r′ we denote the conjugate exponent (that is, 1r+1r=1 ). Also by r*, we denote the critical Sobolev exponent for r, that is,

r=NrNr,if r<N,+,if Nr.

By Ar : W01,r (Ω) → W−1,r(Ω) = W01,r (Ω)*, we denote the nonlinear map defined by

Ar(u),h=Ω|Du|r2(Du,Dh)RNdzfor all u,hW01,r(Ω).

The next proposition recalls the main properties of this map (see, for example, Gasiński-Papageorgiou [20], Problem 2.192, p.219)

Proposition 2.1

The operator Ar(⋅) is bounded (that is, maps bounded sets to bounded sets), continuous, strictly monotone (hence maximal monotone too) and of type (S)+, that is

unw uinW01,r(Ω),lim supnAr(un),unu0unuinW01,r(Ω).

For x ∈ ℝ. we set x± = max{±x, 0}. Then for u W01,p (Ω), we define u±(z) = u(z)± for all zΩ. We have

u±W01,p(Ω),u=u+u,|u|=u++u.

If u W01,p (Ω), then

[u)={hW01,p(Ω):u(z)h(z) for a.a. zΩ}.

If X is a Banach space and φC1(X, ℝ), then we set

Kφ={uX:φ(u)=0}(the critical set of φ).

We say that φ(⋅) satisfies the “C-condition”, if the following property holds:

“Every sequence {un} ⊆ X such that {φ(un)}n≥1 ⊆ ℝ is bounded and (1 + ∥unX)φ′(un) → 0 in X* as n → ∞, admits a strongly convergent subsequence.”

Next, we recall some basic facts about the spectrum of (−Δp, W01,p (Ω)). So, we consider the following nonlinear eigenvalue problem

Δpu(z)=λ^|u(z)|p2u(z) in Ω,u|Ω=0. (2.1)

We say that λ̂ ∈ ℝ is an “eigenvalue” of (−Δp, W01,p (Ω)), if problem (2.1) admits a nontrivial solution û, known as an “eigenfunction” corresponding to λ̂. By σ̂(p) we denote the set of eigenvalues (the “spectrum”) of (−Δp, W01,p (Ω)). There is a smallest eigenvalue denoted by λ̂1(p) which has the following properties:

  1. λ̂1(p) > 0 and it is isolated (that is, there exists ε > 0 such that σ̂(p)⋂(λ̂1(p), λ̂1(p) + ϵ) = ∅).

  2. λ̂1(p) is simple (that is, if û, W01,p (Ω) are eigenfunctions for λ̂1(p), then we have û = θv̂ for some θ ∈ ℝ∖{0}).

  3. λ^1(p)=infDuppupp:uW01,p(Ω),u0. (2.2)

The infimum in (2.2) is realized on the one dimensional eigenspace corresponding to λ̂1(p). By û1(p) we denote the Lp-normalized (that is, ∥û1(p)∥p = 1), positive eigenfunction corresponding to λ̂1(ρ). The nonlinear regularity and the nonlinear maximum principle (see, for example, Gasiński-Papageorgiou [21], pp 737-738), imply that û1 ∈ intC+. Since σ̂(p) ⊆ (0, +∞) is closed and λ̂1(p) is isolated, the second eigenvalue is well-defined by

λ^2(p)=infλ^σ^(p):λ^>λ^1(p).

Using the Ljusternik-Schnirelmann minimax scheme, we produce a whole strictly increasing sequence {λ̂k(p)}k∈ℕ of eigenvalues such that λ̂k(p) → +∞. These eigenvalues are known as “variational eigenvalues” of (−Δp, W01,p (Ω)). We do not know if the variational eigenvalues exhaust σ̂(p). This is the case if p = 2 (linear eigenvalue problem) or if N = 1 (ordinary differential equations). By the nonlinear regularity theory, we know that every eigenfunction û C01 (Ω). Moreover, if λ̂σ̂(p)∖{λ̂1(p)}, then every eigenfunction corresponding to λ̂ is nodal (that is, sign changing).

We will also use a weighted version of the eigenvalue problem (2.1). So, let mL(Ω), m(z) ≥ 0 for a.a. zΩ, m ≠ 0. We consider the following nonlinear eigenvalue problem

Δpu(z)=λ^m(z)|u(z)|p2u(z)in Ω,u|Ω=0. (2.3)

The same results hold for the eigenvalues λ̂(p, m) of problem (2.3). In this case the Rayleigh quotient in the variational characterization of λ̂1(p, m) > 0, is given by

R(u)=DuppΩm(z)|u|pdzfor all uW01,p(Ω).

Then using this fact, we can easily prove the following monotonicity property of the map mλ̂1(p, m).

Proposition 2.2

If m, L(Ω), 0 ≤ m(z) ≤ (z) for a.a. zΩ, m ≠ 0, m, then λ̂1(p, ) < λ̂1(p, m).

Remark 2.1

Evidently, if m(z) = 1 for a.a. zΩ, then λ̂1(p, m) = λ̂1(p).

Now, we introduce the hypotheses on the data of (Eλ).

H0 : rL(Ω), r(z) ≥ 0 for a.a. zΩ, r ≠ 0.

H1 : f : Ω × ℝ → ℝ is a Carathéodory function such that f(z, 0) = 0 for a.a. zΩ and

  1. 0 ≤ f(z, x) ≤ a(z)(1 + xp−1) for a.a. zΩ, all x ≥ 0, with aL(Ω);

  2. there exists m ∈ ℕ, m ≥ 2 such that

    limx+f(z,x)xp1=λ^m(p)uniformly for a.a. zΩ;
  3. if F(z, x) = 0x f(z, s)ds, then there exists θ ∈ (q, p) such that

    0<c0lim infx+pF(z,x)f(z,x)xxθuniformly for a.a. zΩ;
  4. limx0+f(z,x)xq1=0, uniformly for a.a. zΩ.

Remark 2.2

Since our goal is to find positive solutions and the above hypotheses concern the positive semi-axis ℝ+ = [0, +∞), without any loss of generality, we may assume that

f(z,x)=0for a.a. zΩ,all x0. (2.4)

Hypothesis H1(ii) implies that the problem is resonant with respect to a nonprincipal variational eigenvalue λ̂m(p).

As we already mentioned in the Introduction eventually our proof will be topological and to reach that point we will use the frozen variable method.

The topological tool which we will use is the so-called “Leray-Schauder Alternative Principle” (see Papageorgiou-Rădulescu-Repovš [22], Proposition 3.2.22, p.198). So, let X be a Banach space and ξ : XX a map. We say that ξ(⋅) is “compact”, if it is continuous and maps bounded sets to relatively compact sets. The Leray-Schauder Alternative Principle asserts the following:

Theorem 2.3

If X is a Banach space, ξ : XX is compact and

D(ξ)={uX:u=tξ(u)forsomet(0,1)},

then the following alternative holds

  1. D(ξ) is unbounded; or

  2. ξ(⋅) has a fixed point.

We said in the Introduction that in order to solve the frozen problem using the critical point theory, we will need to use the solution of an auxiliary double phase Dirichlet problem. This is the following parametric purely singular problem

Δpu(z)Δqu(z)=λu(z)ηin Ω,u|Ω=0,u>0,λ>0. (Qλ)

From Proposition 11 of Papageorgiou-Rădulescu-Repovš [4], we have the following result concerning problem (Qλ).

Proposition 2.3

For every λ > 0, problem (Qλ) has a unique positive solution uλint C+.

Consider the Banach space C0(Ω) = {uC(Ω) : u∂Ω = 0}. This is an ordered Banach space with positive (order) cone K+ = {uC0(Ω) : u(z) ≥ 0 for all zΩ}. This cone has a nonempty interior given by

intK+={uK+:cud^ufor some cu>0},

with (z) = d(z, Ω) for all zΩ̄. According to Lemma 14.16, p. 355, of Gilbarg-Trudinger [23], there exists δ0 > 0 such that C2(Ωδ0), where Ωδ0 = {zΩ̄ : d(z, Ω) < δ0}. It follows that ∈ intC+ and so from Proposition 4.1.22, p.274, of Papageorgiou-Rădulescu-Repovš [22], we can find 0 < c1 < c2 such that

c1d^u¯λc2d^(recall that u¯λintC+)u¯λintK+.

Let s > N and consider u^1(p)1sK+. Then using Proposition 4.1.22, p. 274, of Papageorgiou-Rădulescu-Repovš [22], we know that we can find c3 > 0 such that

u^1(p)1sc3u¯λ,u¯ληc4u^1(p)ηsfor some c4>0.

From the Lemma in Lazer-McKenna [24], we have u^1(p)ηsLs(Ω),s>N. So, we have

u¯ληLs(Ω),s>N. (2.5)

In the next section, this fact will help us bypass the singular term and use variational tools on the “frozen” problem.

3 The “Frozen” Problem

In this section, we develop the “frozen variable method” for problem (Eλ). So, we fix v C01 (Ω̄) and consider the Carathéodory function

gvλ(z,x)=λxη+r(z)|Dv(z)|τ1+f(z,x).

As we already mentioned in the Introduction, due to the singular term λ xη, the function gvλ (z, x) leads to an energy function which is not C1. For this reason, we use (2.5) and consider the following truncation of gvλ (z, ⋅)

g^vλ(z,x)=λ[u¯λ(z)η+r(z)|Dv(z)|τ1]+f(z,x)if x<u¯λ(z),λ[xη+r(z)|Dv(z)|τ1]+f(z,x)if u¯λ(z)x. (3.6)

Then we consider the following parametric double phase Dirichlet problem (“the frozen problem”)

Δpu(z)Δqu(z)=g^vλ(z,u(z))in Ω,u|Ω=0,u>0,λ>0. ( Evλ )

This problem is variational and its energy (Euler) function is C1. So, we can use the results of critical point theory of ( Evλ ).

Let Svλ denote the set of positive solutions of ( Evλ ). We set G^vλ(z,x)=0xg^vλ(z,s)ds and consider the functional ψ^vλ:W01,p(Ω)R defined by

ψ^vλ(u)=1pDupp+1qDuqqΩG^vλ(z,u(z))dzfor all uW01,p(Ω).

On account of (2.5) (see also Papageorgiou- Smyrlis [25], Proposition 3), we have that

ψ^vλC1(W01,p(Ω)).

Proposition 3.1

If hypotheses H0, H1 hold and λ > 0, then the function ψ^vλ (⋅) satisfies the C-condition.

Proof

We consider a sequence {un}n≥1 W01,p (Ω) such that

|ψ^vλ(un)|c5for some c5>0,all nN, (3.7)
(1+un)(ψ^vλ)(un)0 in W1,p(Ω)=W01,p(Ω) as n. (3.8)

From (3.8) we have

|Ap(un),h+Aq(un),hΩg^vλ(z,un)hdz|εnh1+un,for all hW01,p(Ω),with εn0+. (3.9)

In (3.9) we test with h=unW01,p(Ω) and obtain

Dunppc6for some c6>0,all nN (see (2.4), (2.5) and (3.6)){un}n1W01,p(Ω)is bounded. (3.10)

Next, we show that {un+}n1W01,p(Ω) is bounded. Suppose that this is not true. Passing to a subsequence if necessary, we may assume that

un+as n. (3.11)

We set yn=un+un+ for all n ∈ ℕ. Then ∥yn∥ = 1, yn ≥ 0 for all n ∈ ℕ. So, we may assume that

ynw yin W01,p(Ω)and ynyin Lp(Ω)as n. (3.12)

From (3.9) and (3.10), we have

|Ap(un+),h+Aq(un+),hΩg^vλ(z,un+)hdz|c7hfor some c7>0,all hW01,p(Ω),all nN,Ap(yn),h+1un+pqAq(yn),hΩg^vλ(z,un+)un+p1hdzc7hun+p1,for all hW01,p(Ω). (3.13)

From (2.5),(3.6) and hypothesis H1(i), we see that

g^vλ(,un+())un+p1nNLp(Ω)is bounded. (3.14)

In (3.13) we choose h=ynyW01,p(Ω), pass to the limit as n → ∞ and use (3.11) (recall q < p), (3.12),(3.14), we obtain

limnAp(yn),yny=0ynyin W01,p(Ω)and so y=1,y0 (see Proposition 2.1) (3.15)

From (3.14) and hypothesis H1(ii) and if we pass to a subsequence if necessary, we have

g^vλ(,un+())un+p1w λ^m(p)yp1in Lp(Ω)as n.(see Aizicovici-Papageorgiou-Staicu [26], proof of Proposition 16) (3.16)

We return to (3.13), pass to the limit as n → ∞ and use (3.11) (recall q < p), (3.15), and (3.16). We obtain

Ap(y),h=λ^m(p)Ωyp1hdzfor all hW01,p(Ω),Δpy(z)=λ^m(p)y(z)p1for a.a zΩ,y|Ω=0,ymust be nodal (since m2),

a contradiction to (3.15).

Therefore, {un+}n1W01,p(Ω) is bounded and so from (3.10), it follows that

{un}n1W01,p(Ω)is bounded.

We may assume that

unw uin W01,p(Ω)and unuin Lp(Ω)as n. (3.17)

In (3.9) we choose h = unu W01,p (Ω), pass to the limit as n → ∞ and use (3.17), we have

limnAp(un),unu+Aq(un),unu=0,lim supnAp(un),unu+Aq(u),unu0 (since Aq()is monotone),lim supnAp(un),unu0(see (3.17)),unuin W01,p(Ω)as n(see Proposition 2.1).

This proves that the functional ψ^vλ (⋅) satisfies the C-condition.□

The next proposition will help us satisfy the mountain pass geometry for the functional ψ^vλ (⋅).

Proposition 3.2

If Hypotheses H0, H1 hold and λ > 0, then ψ^vλ (tû1(p)) → −∞ as t → +∞.

Proof

We have

ddx(F(z,x)xp)=f(z,x)xppxp1F(z,x)x2p=f(z,x)xpF(z,x)xp+1for a.a. zΩ,all x>0. (3.18)

By hypothesis H1(iii), we can find M > 0 and c8 ∈ (0,c0) such that

f(z,x)xpF(z,x)c8xθfor a.a. zΩ,all xM. (3.19)

Using (3.19) in (3.18), we have

ddx(F(z,x)xp)c8xp+1θfor a.a. zΩ,all xM,F(z,t)tpF(z,x)xpc8pθ1tpθ1xpθfor a.a. zΩ,all txM. (3.20)

Hypothesis H1(ii) implies that

limt+pF(z,x)tp=λ^m(p)uniformly for a.a. zΩ. (3.21)

In (3.20), we let t → +∞ and use (3.21). We obtain

λ^m(p)pF(z,x)xpc8pθ1xpθλ^m(p)xppF(z,x)xθpc8pθfor a.a. zΩ,all xM,lim supx+λ^m(p)xppF(z,x)xpc9uniformly for a.a. zΩwith c9=pc8pθ. (3.22)

For t > 0, we have

ψ^vλ(tu^1(p))=tppλ^1(p)+tqqDu^1(p)qqΩG^vλ(z,tu^1(p))dz.

Recall that û1(p) ∈ int C+. So, we can find t ≥ 1 big such that ūλtû1(p) (see [22], p.274). Then from (3.6) and since m ≥ 2, we have

ψ^vλ(tu^1(p))tppλ^m(p)+tqqDu^1(p)qqΩG^vλ(z,tu^1(p))dz1pΩλ^m(p)(tu^1(p))ppF(z,tu^1(p))+1qtqDu^1(p)qq,ψ^vλ(tu^1(p))tθ1pΩλ^m(p)(tu^1(p))pF(z,tu^1(p))tθdz+1qtθqDu^1(p)qqlim supt+ψ^vλ(tu^1(p))tθc10,

for some c10 > 0 (see (3.22), using Fatou′s Lemma and recall q < θ),

ψ^vλ(tu^1(p))as t+.

Remark 3.1

The above proof reveals that the resonance at λ̂m(p) occurs from the right of the eigenvalue in the sense that

pF(z,x)λ^m(p)xp+uniformly for a.a. zΩ,as x+.

Proposition 3.3

If hypotheses H0, H1 hold and λ > 0, then Kψ^vλ ⊆ [ūλ) ⋂ int C+.

Proof

If Kψ^vλ = ∅, then the result is trivially true. So, suppose that Kψ^vλ ≠ ∅ and let u Kψ^vλ .We have

Ap(u),h+Aq(u),h=Ωg^vλ(z,u)hdzfor all hW01,p(Ω). (3.23)

In (3.23), we choose h = (ūλu)+ W01,p (Ω). We have

Ap(u),(u¯λu)++Aq(u),(u¯λu)+=Ωλu¯λη+λr(z)|Dv|τ1+f(z,u)(u¯λu)+dz(see (3.6))Ωλu¯λη(u¯λu)+dz(see H0,H1(i))=Ap(u¯λ),(u¯λu)++Aq(u¯λ),(u¯λu)+(see Proposition 2.3)u¯λu. (3.24)

From (3.24), (3.6), and (3.23), it follows that

Δpu(z)Δqu(z)=λu(z)η+r(z)|Dv(z)|τ1+f(z,u(z))for a.a. zΩ. (3.25)

Theorem 7.1, p.286, of Ladyzhenskaya-Uraltseva [29], implies that uL(Ω). Let

Kλ(z)=λu(z)η+r(z)|Dv(z)|τ1+f(z,u(z)).

On account of (3.24), (2.5) and hypotheses H0, H1(i), we have that

KλLs(Ω),s>N.

We consider the following linear Dirichlet problem

Δw(z)=Kλ(z)in Ω,w|Ω=0. (3.26)

Theorem 9.15, p.241, of Gilbarg-Trudinger [23] implies that problem (3.26) admits a unique solution w W02,s (Ω). The Sobolev embedding theorem says that

W2,s(Ω)C1,α(Ω¯)compactly with α=1Ns(0,1),wC1,α(Ω¯)C01(Ω¯).

Let σ(z) = Dw(z). Then σC0,α(Ω̄,ℝN). Using this function, we can rewrite (3.25) as

div|Du|p2Du+|Du|q2Duσ=0in Ω.

Then the nonlinear regularity theory of Lieberman [27] and (3.24) imply that u ∈ intC+. So, we conclude that

Kψ^vλ[u¯λ)intC+.

Next, we are ready to show the nonemptiness of Svλ when λ > 0 is small.

Proposition 3.4

If hypotheses H0, H1 hold and v C01 (Ω̄), then there exists λ* > 0 such that for all λ ∈ (0,λ*), we have

Svλ[u¯λ)intC+.

Proof

Let u W01,p (Ω) with ∥u∥ ≤ 1. We have

ψ^vλ(u)=1pDupp+1qDuqq{u<u¯λ}G^vλ(z,u)dz{uu¯λ}G^vλ(z,u)dzλ{uu¯λ}u¯λ1ηdz{uu¯λ}F(z,u(z))dz1pDupp+1qDuqqλ(2η)1ηΩ|u|1ηdzλΩr(z)|Dv(z)|τ1udzΩF(z,u)dz(see hypotheses H1(i),(iv)). (3.27)

On account of hypotheses H1(i), (iv), we see that given ε > 0, we can find c11 > 0 such that

f(z,x)xε|x|q+c11|x|pfor a.a. zΩ,all xR(see (2.4)). (3.28)

Returning to (3.27) and using (3.28) with ε > 0 small, we obtain from the Poincaré inequality and ∥u∥ ≤ 1

ψ^vλ(u)12qDuqqc12Duppλc13Dup1ηfor some c12,c13>0. (3.29)

Let m=infDvqq:vW01,p(Ω),v=Dvp=1. Obviously, m > 0. We choose ρ ∈ (0,1) such that for all v W01,p (Ω) with ∥v∥ = 1 and u = ρv

12qDuqqc12Dupp=12qρqDvqqc12ρpρq{m2qc12ρpq}η^0>0(recallq<p)

and then choose λ* > 0 small so that λ c13ρ1−η < η̂0 for all λ ∈ (0,λ*). From (3.29), it follows that

ψ^vλ(0)=0<inf[ψ^vλ(u):u=ρ]for all λ(0,λ). (3.30)

Then (3.30) together with Propositions 3.1 and 3.2, permit the use of the mountain pass theorem. So, we can find u^vλW01,p(Ω) such that

u^vλKψ^vλ[u¯λ)intC+(see Proposition (3.3)),u^vλSvλ(see (3.6)) and Svλ[u¯λ)intC+.

The next proposition suggests a canonical way to choose a solution from Svλ as v C01 (Ω̄) moves.

Proposition 3.5

If hypotheses H0, H1 hold, v C01 (Ω̄) and λ ∈ (0,λ*), then Svλ admits a smallest element u~vλ int C+, that is

u~vλuforalluSvλ.

Proof

From Proposition 19 of Papageorgiou-Rădulescu-Repovš [4], we know that the solution set Svλ is downward directed (that is given u1, u2 Svλ we can find u Svλ such that uu1, uu2). Then Lemma 3.10, p.178 of Hu-Papageorgiou [28] implies that we can find a decreasing sequence {un}n≥1 Svλ such that infn1un=infSvλ. We have

Ap(un),h+Aq(un),h=Ωg^vλ(z,un)hdzfor all hW01,p(Ω),all nN, (3.31)
u¯λunu1for all nN(see Proposition (3.3)). (3.32)

If we test (3.31) with h = un W01,p (Ω) and use (3.32) and hypothesis H1(i), we infer that {un}n≥1 W01,p (Ω) is bounded. So, we may assume that

unw u~vλin W01,p(Ω)and unu~vλin Lp(Ω)as n. (3.33)

In (3.31), we choose h=unu~vλW01,p(Ω), pass to the limit as n → ∞ and use (3.33). Reasoning as in the last part of the proof of Proposition 3.1, we obtain

unu~vλin W01,p(Ω). (3.34)

Pass to the limit as n → ∞ in (3.31) and using (3.34), we obtain

Ap(u~vλ),h+Aq(u~vλ),h=Ωg^vλ(z,u~vλ)hdzfor all hW01,p(Ω), (3.35)
u¯λu~vλ(see (3.32)). (3.36)

From (3.35) and (3.36), it follows that

u~vλSvλand u~vλ=infSvλ.

We can define the minimal solution map ξλ : C01 (Ω̄) → C01 (Ω̄) by

ξλ(u)=u~vλintC+(λ(0,λ)).

Clearly a fixed point of this map will be a positive smooth solution of (Eλ) λ ∈ (0,λ*). To produce a fixed point of ξλ(⋅), we will use Theorem 2.3 (Leray-Schauder Alternative Principle). This is done in the next section.

4 Positive Solution

According to Theorem 2.3, to produce a fixed point for the minimal solution map, we need to show that ξλ(⋅) is compact and that D(ξλ) is bounded (λ ∈ (0, λ*)).

First, we show that ξλ(⋅) is compact. To this end the following proposition will be helpful.

Proposition 4.1

If hypotheses H0, H1 hold, vnv in C01 (Ω̄), λ ∈ (0, λ*) and u Svλ , then for m ≥ 2 big we can find un Svnλ , n ∈ ℕ such that unu in C01 (Ω̄).

Proof

From Proposition 3.3, we know that u ∈ [ūλ) ⋂ intC+. We consider the following Dirichlet problem

Δpw(z)Δqw(z)=g^vnλ(z,u(z))in Ω,w|Ω=0,u>0,λ>0.

We have g^vnλ(,u())0,g^vnλ(,u())0 (see hypothesis H1(i) and recall ūλ ∈ int C+, u¯λη Ls(Ω)). In fact we have

g^vnλ(,u())Ls(Ω)with s>N.

Note that s<N=NN1<p. So, it follows that

Ls(Ω)W1,p(Ω)=(W01,p(Ω))compactly (see Gasinski-Papageorgiou [21], Lemma 2.2.27, p.141),g^vnλ(,u())W1,p(Ω){0}for all nN.

We consider the nonlinear map V : W01,p (Ω) → W−1,p(Ω) defined by

V(u)=Ap(u)+Aq(u)for all uW01,p(Ω).

According to Proposition 2.1, this map is continuous, strictly monotone (hence maximal monotone too) and coercive. It follows that V(⋅) is surjective (see Papageorgiou-Rădulescu-Repovš [22], Corollary 2.8.7, p.135). So, we can find a unique wn W01,p (Ω) (uniqueness is a consequence of the strict monotonicity of V(⋅)),wn ≥ 0,wn ≠ 0 such that

V(wn)=g^vnλ(,u())for all nN. (4.37)

Clearly, {g^vnλ(,u())}nNLs(Ω) is bounded. So, by Theorem 9.16, p.241 and Lemma 9.17, p.242, of Gilbarg-Trudinger [23], there exists c15 > 0 such that

wnW2,s(Ω)and wnW2,s(Ω)c15for all nN.

By the Sobolev embedding theorem, we have

W2,s(Ω)C1,α(Ω¯)α=1Ns(0,1)compactly.

So, we see that at least for a subsequence, we have

wnu~in C1,α(Ω¯)as n. (4.38)

Passing to the limit as n → ∞ in (4.37) and using (4.38), we obtain

V(u~)=g^vλ(,u()),Δpu~Δqu~=g^vλ(z,u~)in Ω,u~|Ω=0,u~=u(since the problem has a unique solution and uSvλ).

Therefore for the original sequence, we have

wnuin C01(Ω¯).

Next, we consider the following Dirichlet problem

ΔpuΔqu=g^vnλ(z,wn)in Ω,u~|Ω=0.

Reasoning as above, we see that this problem has a unique solution wn1W2,s(Ω),wn10,wn10 and

wn1uin C01(Ω¯) as n.

We continue this way and generate a sequence {wnk}n,kN such that

Ap(wnk)+Aq(wnk)=g^vnλ(,wnk1())for all n,kN(wn0=wn),wnk0, (4.39)
wnkuin C01(Ω¯)as nfor all kN. (4.40)

For each n ∈ ℕ, we consider the sequence {wnk}kN. We claim that this sequence is bounded in W01,p (Ω). Arguing by contradiction, we may assume that

wnkas kand {wnk}kNis increasing. (4.41)

We set yk=wnkwnk,kN. Then ∥yk∥ = 1,yk ≥ 0 for all k ∈ ℕ. So, we may assume that

ykw yin W01,p(Ω)and ykyin Lp(Ω)as k,y0. (4.42)

From (4.39) we have

Ap(yk)+1wnkpqAq(yk)=g^vnλ(,wnk1())wnkp1for all kN. (4.43)

On account of hypothesis H1(i) and of (4.38), (4.39), we see that

g^vnλ(,wnk1())wnkp1kNLp(Ω)is bounded. (4.44)

we act on (4.43) with yky W01,p (Ω) and then pass to the limit as k → ∞ and use (4.41), (4.42), and (4.44). We obtain

limkAp(yk),yky=0(recall q<p)ykyin W01,p(Ω)as kand y=1,y0(see Proposition 2.1). (4.45)

From (4.44) and hypothesis H1(ii), at least for a subsequence, we have

g^vnλ(,wnk1())wnkp1w λ^m(p)(ηy)p1 in Lp(Ω),η(0,1)as k(see Aizicovici-Papageorgiou-Staicu [26], proof of Proposition 16). (4.46)

In (4.43), we pass to the limit as k → ∞ and use (4.45), (4.41) (recall q < p) and (4.46), we obtain

Ap(y)=λ^m(p)(ηy)p1,Δpy(z)=λ^m(p)(ηy(z))p1in Ω,y|Ω=0.

If m ≥ 2 is big so that η>(λ^1(p)λ^m(p))1p1, then it follows that y = 0 or y is nodal, both possibilities leading to a contradiction (see (4.45)). This proves that for every n ∈ ℕ, {wnk}kN W01,p (Ω) is bounded. Then Theorem 7.1, p.286, of Ladyzhenskaya-Uraltseva [29] implies that {wnk}kN L(Ω) is bounded. The nonlinear regularity theory of Lieberman [27] implies that for each n ∈ ℕ, we can find α ∈ (0, 1) and c16 > 0 such that

wnkC01,α(Ω¯),wnkC01,α(Ω¯)c16for all kN.

As before, the compact embedding of C01,α(Ω¯) into C01 (Ω̄), implies that at least for a subsequence, we have

wnkunin C01(Ω¯)as k. (4.47)

Passing to the limit as k → ∞ in (4.39), we obtain

Ap(un)+Aq(un)=g^vnλ(,un())kN,Δpun(z)Δqun(z)=g^vnλ(z,un(z))in Ω,un|Ω=0. (4.48)

Moreover, as in the proof of Proposition 3.3, using (3.6) and the fact that f ≥ 0 (see hypothesis H1(i)), we have

u¯λwnkfor all n,kN,u¯λunfor all nN,unSvnλint C+for all nN.

As above, using (4.48) and a contradiction argument, we show that {un}n≥1 W01,p (Ω) is bounded, hence relatively compact in C01 (Ω̄) (nonlinear regularity). Then the double limit lemma (see Gasiński-Papageorgiou [30], Problem 1.175, p.61), implies that

unuin C01(Ω¯)as n,unSvnλfor all nN.

Using this proposition, we can show that the minimal solution map ξλ(⋅) is compact.

Proposition 4.2

If hypotheses H0, H1 hold, m ≥ 2 is big and λ ∈ (0, λ*), then the map ξλ : C01 (Ω̄) → C01 (Ω̄) is compact.

Proof

First, we show that ξλ(⋅) maps bounded sets in C01 (Ω̄) to relatively compact sets in C01 (Ω̄). So, let D C01 (Ω̄) be bounded. Then as in the proof of Proposition 4.1, with a contradiction argument and using the fact that m ≥ 2 is big, we show that ξλ(D) ⊆ W01,p (Ω) is bounded and from this by the nonlinear regularity theory (see [27, 29]), we obtain that ξλ(D) ⊆ C01 (Ω̄) relatively compact.

Next, we show that ξλ(⋅) is continuous. So, let vnv in C01 (Ω̄) and let u~vλ = ξλ(v). According to Proposition 4.1, we can find un Svnλ , n ∈ ℕ such that

unu~vλ. (4.49)

We know that

ξλ(vn)unfor all nN, (4.50)

and from the first part if the proof, we have that

ξλ(vn)nNC01(Ω¯)is relatively compact.

So, by passing to a subsequence if necessary, we have

ξλ(vn)uvλin C01(Ω¯)as n (4.51)

Clearly, uvλSvλ and from (4.49), (4.50),(4.51), we have

uvλu~vλuvλ=u~vλ=ξλ(v).

So, for the original sequence, we have

ξλ(vn)ξλ(v)in C01(Ω¯)as n,ξλ()is continuous.

Therefore, we have proved that the minimal map ξλ(⋅) is compact.

We introduce the following set

D(ξλ)={uC01(Ω¯):u=tξλ(u),0<t<1}.

Proposition 4.3

If hypotheses H0, H1 hold, m ≥ 2 is big and λ ∈ (0, λ*), then D(ξλ) ⊆ C01 (Ω̄) is bounded.

Proof

As before (see the proof of Proposition 4.1), it suffices to show that D(ξλ) ⊆ W01,p (Ω) is bounded. Then the nonlinear regularity theory (see [27, 29]) implies that D(ξλ) ⊆ C01 (Ω̄) is bounded (in fact, relatively compact).

Again, we argue indirectly. So, suppose we could find {un}n≥1D(ξλ) such that

unas n. (4.52)

we have

1tnun=ξλ(un)with 0<tn<1for all nN.

Therefore, we have

1tnp1Ap(un)+1tnq1Aq(un)=tnηunη+r(z)|Dun|τ1+f(,1tnun())for all nN,Ap(un)+tnpqAq(un)=tnp+η1unη+tnp1r(z)|Dun|τ1+tnp1f(,1tnun())for all nN, (4.53)

Let yn=unun,nN. Then we have ∥yn∥ = 1, yn ≥ 0 for all n ∈ ℕ and so we may assume that

ynw yin W01,p(Ω)and ynyin Lp(Ω)as n. (4.54)

From (4.53) we have

Ap(yn)+(tnun)pqAq(yn)=tnηunp1unη+tnp1unpτr(z)|Dyn|τ1+f(,1tnun())1tnunp1for all nN. (4.55)

Note that

0tnηunp1unη1unp1u~ληLs(Ω)s>N(see Proposition 3.3 and (2.5)). (4.56)

Also, from hypothesis H1(i), we have

f(,1tnun())1tnunp1n1Lp(Ω)is bounded. (4.57)

We test (4.55) with yny W01,p (Ω), pass to the limit as n → ∞ and use (4.52) (recall q, τ < p), (4.54), (4.56) and (4.57). We obtain

limnAp(yn),yny=0,ynyin W01,p(Ω)as nand so y=1,y0(see Proposition 2.1). (4.58)

From (4.57) and hypothesis H1(ii), we have

f(,1tnun())1tnunp1w λ^m(p)yp1in Lp(Ω)as n(see[26]). (4.59)

If in (4.55) we pass to the limit as n → ∞ and use (4.58), (4.59) and (4.52) (recall q, τ < p), we obtain

Ap(y)=λ^m(p)yp1,Δpy(z)=λ^m(p)y(z)p1 in Ω,y|Ω=0,y is zero or nodal (since m2),

both contradicting (4.58). This proves that D(ξλ) ⊆ W01,p (Ω) is bounded, from which it follows the boundedness in C01 (Ω̄) as explained in the beginning of the proof.□

On account of Propositions 4.2 and 4.3, we can apply Theorem 2.3 and produce a fixed point for ξλ(⋅), λ ∈ (0, λ*). So, we can state the following existence theorem for problem (Eλ).

Theorem 4.1

If hypotheses H0, H1 hold and m ≥ 2 is big, then there exists λ* > 0 such that for all λ ∈ (0, λ*) problem (Eλ) has a positive solution uλint C+.

Remark 4.2

There are three interesting open problems concerning problem (Eλ)

  1. What can be said if we have resonance with respect to the principal eigenvalue λ1^ (p) > 0 (that is, m = 1)?

  2. Is it possible to have τ = p in the drift term?

  3. Is it possible to have an indefinite drift term (that is, drop the hypothesis that r ≥ 0)? In the direction see Hu-Papageorgiou [31].

Acknowledgement

The authors wish to thank the two knowledgeable referees for their corrections and remarks.

The work was supported by NNSF of China Grant No. 11671101, NSF of Guangxi Grant No. 2018GXNSFDA138002.

References

[1] Y.R. Bai, L. Gasiński, N.S. Papageorgiou, Nonlinear Dirichlet problems with the combined effects of singular and convection terms. Electron. J. Differential Equations 2019, Paper No. 57, 13 pp.Search in Google Scholar

[2] Z.H. Liu, D. Motreanu, S.D. Zeng, Positive solutions for nonlinear singular elliptic equations of p -Laplacian type with dependence on the gradient. Calc. Var. Partial Differential Equations 58(2019), no.1, Art. 28, 22 pp.10.1007/s00526-018-1472-1Search in Google Scholar

[3] N.S. Papageorgiou, V.D. Rădulescu, D.D. Repovš, Positive solutions for nonlinear Neumann with sigular terms and convection, J. Math. Pures Appl. 136 (2020) 1-21.10.1016/j.matpur.2020.02.004Search in Google Scholar

[4] N.S. Papageorgiou, V.D. Rădulescu, D.D. Repovš, Nonlinear nonhomogeneous singular problems. Calc. Var. Partial Differential Equations 59 (2020), no. 1, Art. 9, 31 pp.10.1007/s00526-019-1667-0Search in Google Scholar

[5] N.S. Papageorgiou, C. Vetro, F. Vetro, Parametric nonlinear singular Dirichlet problems. Nonlinear Anal. Real World Appl. 45(2019), 239-254.10.1016/j.nonrwa.2018.07.006Search in Google Scholar

[6] N.S. Papageorgiou, C. Zhang, Noncoercive resonant (p, q)-equations with concave terms, Adv. Nonlinear Anal.,9(2020),228-249.10.1515/anona-2018-0175Search in Google Scholar

[7] L. Gasiński, P. Winkert, Existence and uniqueness results for double phase problems with convection term, J. Differential Equations 268 (2020), no. 8, 4183-4193.10.1016/j.jde.2019.10.022Search in Google Scholar

[8] S. A. Marano, P. Winkert, On a quasilinear elliptic problem with convection term and nonlinear boundary condition, Nonlinear Anal. 187 (2019), 159-169.10.1016/j.na.2019.04.008Search in Google Scholar

[9] N.S. Papageorgiou, V.D. Rădulescu, D.D. Repovš, Double-phase problems with reaction of arbitrary growth. Z. Angew. Math. Phys. 69 (2018), no. 4, Paper No. 108, 21 pp.10.1007/s00033-018-1001-2Search in Google Scholar

[10] N.S. Papageorgiou, V.D. Rădulescu, D.D. Repovš, Double-phase problems and a discontinuity property of the spectrum. Proc. Amer. Math. Soc. 147 (2019), no. 7, 2899-2910.10.1090/proc/14466Search in Google Scholar

[11] N.S. Papageorgiou, V.D. Rădulescu, D.D. Repovš, Positive solutions for nonlinear parametric singular Dirichlet problems. Bull. Math. Sci. 9 (2019), no. 3, 1950011, 21 pp.10.1142/S1664360719500115Search in Google Scholar

[12] N.S. Papageorgiou, V.D. Rădulescu, D.D. Repovš, Ground state and nodal solutions for a class of double phase problems. Z. Angew. Math. Phys. 71 (2020), no. 1, Paper No. 15, 15 pp.10.1007/s00033-019-1239-3Search in Google Scholar

[13] N.S. Papageorgiou, C. Vetro, F. Vetro, Multiple solutions with sign information for semilinear Neumann problems with convection, Rev. Mat. Complut. 33 (2020), 19-38.10.1007/s13163-019-00312-3Search in Google Scholar

[14] N.S. Papageorgiou, C. Vetro, F. Vetro, Nonlinear Robin problems with unilateral constraints and dependence on the gradient, Electron. J. Differential Equations 2018 (2018), no. 182, 1-14.Search in Google Scholar

[15] M.A. Ragusa, A. Tachikawa, Regularity for minimizers for functionals of double phase with variable exponents. Adv. Nonlinear Anal. 9 (2020), no. 1, 710-728.10.1515/anona-2020-0022Search in Google Scholar

[16] X.H. Tang, X.Y. Lin, J.S. Yu, Existence of a bound state solution for quasilinear Schrödinger equations, J. Dyn. Differ. Equ. 31 (2019), 369-383.10.1007/s10884-018-9662-2Search in Google Scholar

[17] C. Vetro, Pairs of nontrivial smooth solutions for nonlinear Neumann problems, Appl. Math. Lett. 103:106171 (2020), 1-7.10.1016/j.aml.2019.106171Search in Google Scholar

[18] C. Vetro, Semilinear Robin problems driven by the Laplacian plus an indefinite potential, Complex Var. Elliptic Equ. 65 (2020), 573-587.10.1080/17476933.2019.1597066Search in Google Scholar

[19] C. Vetro, F. Vetro, On problems driven by the (p(⋅), q(⋅))-Laplace operator, Mediterr. J. Math. 17:24 (2020), 1-11.10.1007/s00009-019-1448-1Search in Google Scholar

[20] L. Gasiński, N.S. Papageorgiou, Exercises in analysis. Part 2. Nonlinear analysis. Problem Books in Mathematics. Springer, Cham, 2016. viii+1062 pp.10.1007/978-3-319-27817-9Search in Google Scholar

[21] L. Gasiński, N.S. Papageorgiou, Nonlinear analysis. Series in Mathematical Analysis and Applications, 9. Chapman & Hall/CRC, Boca Raton, FL, 2006. xii+971 ppSearch in Google Scholar

[22] N.S. Papageorgiou, V.D. Rădulescu, D.D. Repovš, Nonlinear Analysis: Theory and Methods. Springer Monographs in Mathematics. Springer, Cham, 2019. xi+577 pp.10.1007/978-3-030-03430-6Search in Google Scholar

[23] D. Gilbarg D, N.S. Trudinger, Elliptic partial differential equations of second order. Classics in Mathematics. Springer-Verlag, Berlin, 1998.Search in Google Scholar

[24] A.C. Lazer, P.J. McKenna, On a singular nonlinear elliptic boundary-value problem. Proc. Amer. Math. Soc. 111 (1991), no. 3, 721-730.10.1090/S0002-9939-1991-1037213-9Search in Google Scholar

[25] N.S. Papageorgiou, G. Smyrlis, A bifurcation-type theorem for singular nonlinear elliptic equations. Methods Appl. Anal. 22 (2015), no. 2, 147-170.10.4310/MAA.2015.v22.n2.a2Search in Google Scholar

[26] S. Aizicovici, N.S. Papageorgiou, V. Staicu, Degree Theory for Operators of Monotone Type and Nonlinear Elliptic Equations with Inequality Constraints, Memoirs of AMS, Vol.196, no.905, 2008.10.1090/memo/0915Search in Google Scholar

[27] G.M. Lieberman, The natural generalization of the natural conditions of Ladyzhenskaya and Uraľtseva for elliptic equations. Comm. Partial Differential Equations 16 (1991), no. 2-3, 311-361.10.1080/03605309108820761Search in Google Scholar

[28] S.C. Hu, N.S. Papageorgiou, Handbook of multivalued Analysis: Volume I Theory, Kluwer Academic Publishers, Dordrecht Boston, London, 1997.10.1007/978-1-4615-6359-4Search in Google Scholar

[29] O.A. Ladyzhenskaya, N.N. Uraltseva, Linear and quasilinear elliptic equations. Translated from the Russian by Scripta Technica, Inc. Translation editor: Leon Ehrenpreis Academic Press, New York-London 1968 xviii+495 pp.Search in Google Scholar

[30] L. Gasiński, N.S. Papageorgiou, Exercises in analysis. Part 1. Problem Books in Mathematics. Springer, Cham, 2014. x+1037 pp.10.1007/978-3-319-06176-4_1Search in Google Scholar

[31] S.C. Hu, N.S. Papageorgiou, Positive solution for nolinear Dirichlet problems with convection, Appl. Math. Optim., 10.1007/s00245-018-9534-5Search in Google Scholar

Received: 2020-04-07
Accepted: 2020-05-31
Published Online: 2020-07-17

© 2021 Zhenhai Liu and Nikolaos S. Papageorgiou, published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 29.5.2024 from https://www.degruyter.com/document/doi/10.1515/anona-2020-0108/html
Scroll to top button