Skip to content
BY 4.0 license Open Access Published by De Gruyter November 26, 2020

Convergence analysis for double phase obstacle problems with multivalued convection term

  • Shengda Zeng EMAIL logo , Yunru Bai , Leszek Gasiński and Patrick Winkert

Abstract

In the present paper, we introduce a family of the approximating problems corresponding to an elliptic obstacle problem with a double phase phenomena and a multivalued reaction convection term. Denoting by 𝓢 the solution set of the obstacle problem and by 𝓢n the solution sets of approximating problems, we prove the following convergence relation

∅≠w-lim supn→∞Sn=s-lim supn→∞Sn⊂S,

where w-lim supn→∞ 𝓢n and s-lim supn→∞ 𝓢n denote the weak and the strong Kuratowski upper limit of 𝓢n, respectively.

MSC 2010: 35J20; 35J25; 35J60

1 Introduction

Recently, based on a surjectivity result for pseudomonotone operators obtained by Le [25], the authors [44] have studied the nonemptyness, boundedness and closedness of the set of weak solutions to the following double phase problem with a multivalued convection term and obstacle effect

−div|∇u|p−2∇u+μ(x)|∇u|q−2∇u∈f(x,u,∇u)in Ω,u(x)≤Φ(x)in Ω,u=0on ∂Ω, (1.1)

where Ω ⊆ ℝN is a bounded domain with Lipschitz boundary ∂Ω, 1 < p < q < N, μ: Ω → [0, ∞) is Lipschitz continuous, f: Ω × ℝ × ℝN → 2ℝ is a multivalued function depending on the gradient of the solution and Φ: Ω → ℝ+ is a given function, see Section 3 for the precise assumptions.

As the obstacle effect leads to various difficulties in obtaining the exact and numerical solutions, it is reasonable to consider some appropriate approximating methods to overcome/avoid the obstacle effect. In the present paper, we are going to propose a family of approximating problems corresponding to (1.1) and deliver an important convergence theorem which indicates that the solution set of the obstacle problem can be approximated by the solutions of perturbation problems. More precisely, let {ρn} be a sequence of positive numbers such that ρn → 0 as n → ∞ and for each n ∈ ℕ, we consider the following problem

−div|∇u|p−2∇u+μ(x)|∇u|q−2∇u+1ρnu(x)−Φ(x)+∈f(x,u,∇u)in Ω,u=0on ∂Ω. (1.2)

Denoting by 𝓢 and 𝓢n the sets of solutions to problems (1.1) and (1.2), respectively, we shall establish the relations between the sets 𝓢, w-lim supn→∞ 𝓢n (being the weak Kuratowski upper limit of 𝓢n) and s-lim supn→∞ 𝓢n (being the strong Kuratowski upper limit of 𝓢n), see Definition 2.2.

The introduction of so-called double phase operators goes back to Zhikov [46] who described models of strongly anisotropic materials by studying the functional

u↦∫|∇u|p+μ(x)|∇u|qdx. (1.3)

The integral functional (1.3) is characterized by the fact that the energy density changes its ellipticity and growth properties according to the point in the domain. More precisely, its behavior depends on the values of the weight function μ(⋅). Indeed, on the set {x ∈ Ω : μ(x) = 0} it will be controlled by the gradient of order p and in the case {x ∈ Ω : μ(x) ≠ 0} it is the gradient of order q. This is the reason why it is called double phase.

Functionals of the expression (1.3) have been studied more intensively in the last five years. Concerning regularity results, we refer, for example, to the works of Baroni-Colombo-Mingione [4, 5, 6], Baroni-Kuusi-Mingione [7], Cupini-Marcellini-Mascolo [15], Colombo-Mingione [13], [14], Marcellini [28, 29] and the references therein.

Double phase differential operators and corresponding energy functionals appear in several physical applications. For example, in the elasticity theory, the modulating coefficient μ(⋅) dictates the geometry of composites made of two different materials with distinct power hardening exponents q and p, see Zhikov [47]. We also refer to other applications which can be found in the works of Bahrouni-Rădulescu-Repovš [1] on transonic flows, Benci-D’Avenia-Fortunato-Pisani [8] on quantum physics and Cherfils-Il’yasov [9] on reaction diffusion systems.

Existence and uniqueness results have been recently obtained by several authors. In the case of single-valued equations with or without convection term, we refer to Colasuonno-Squassina [12], Gasiński-Papageorgiou [16, 17], Gasiński-Winkert [19, 20, 21], Liu-Dai [27], Perera-Squassina [39], Papageorgiou-Vetro-Vetro [34, 35] and the references therein.

Finally, papers or monographs dealing with certain types of double phase problems or multivalued problems can be found in Bahrouni-Rădulescu-Repovš [1], Bahrouni-Rădulescu-Winkert [2], [3], Carl-Le-Motreanu [10], Cencelj-Rădulescu-Repovš [11], Clarke [22], Gasiński-Papageorgiou [18], Marino-Winkert [30], Papageorgiou-Rădulescu-Repovš [32, 33], Papageorgiou-Vetro-Vetro [37], Rădulescu [40], Vetro [41], Vetro-Vetro [42], Zhang-Rădulescu [45], Zeng-Bai-Gasiński-Winkert [43] and the references therein.

The paper is organized as follows. In Section 2 we recall the definition of the Musielak-Orlicz spaces L𝓗(Ω) and its corresponding Sobolev spaces W1,𝓗(Ω) and we recall the definition of the Kuratowski lower and upper limit, respectively. In Section 3 we present the full assumptions on the data of problem (1.2), give the definition of weak solutions for (1.1) as well as (1.2) and state and prove our main result, see Theorem 3.4.

2 Preliminaries

Let Ω be a bounded domain in ℝN and let 1 ≤ r ≤ ∞. In what follows, we denote by Lr(Ω) : = Lr(Ω;ℝ) and Lr(Ω;ℝN) the usual Lebesgue spaces endowed with the norm ∥⋅∥r. Moreover, W1,r(Ω) and W01,r (Ω) stand for the Sobolev spaces endowed with the norms ∥⋅∥1,r and ∥⋅∥1,r,0, respectively. For any 1 < r < ∞ we denote by r′ the conjugate of r, that is, 1r+1r′=1.

For the weight function Îź and powers p, q we will assume that:

H(μ): μ : Ω → ℝ+ := [0, ∞) is Lipschitz continuous and 1 < p < q < N are chosen such that

qp<1+1N.

We consider the function 𝓗: Ω × ℝ+ → ℝ+ defined by

H(x,t)=tp+μ(x)tqfor all (x,t)∈Ω×R+.

Based on the definition of 𝓗 we are able to introduce the Musielak-Orlicz space L𝓗(Ω) given by

LH(Ω)=u | u:Ω→R is measurable and ρH(u):=∫ΩH(x,|u|)dx<+∞,

endowed with the Luxemburg norm

∥u∥H=infτ>0|ρHuτ≤1.

We know that L𝓗(Ω) is uniformly convex and so a reflexive Banach space. In addition, we introduce the seminormed function space

Lμq(Ω)=u | u:Ω→R is measurable and ∫Ωμ(x)|u|qdx<+∞,

which is equipped with the seminorm ∥⋅∥q,μ given by

∼u∼q,Ο=∍ΊΟ(x)|u|qdx1q.

It is known that the embeddings

Lq(Ω)↪LH(Ω)↪Lp(Ω)∩Lμq(Ω)

are continuous, see Colasuonno-Squassina [12, Proposition 2.15 (i), (iv) and (v)]. Taking into account these embeddings we have the inequalities

min∥u∥Hp,∥u∥Hq≤∥u∥pp+∥u∥q,μq≤max∥u∥Hp,∥u∥Hq (2.1)

for all u ∈ L𝓗(Ω).

By W1,𝓗(Ω) we denote the corresponding Sobolev space which is defined by

W1,H(Ω)=u∈LH(Ω):|∇u|∈LH(Ω)

equipped with the norm

∥u∥1,H=∥∇u∥H+∥u∥H,

where ∥∇u∥𝓗 = ∥|∇u|∥𝓗.

By W01,H (Ω) we denote the completion of C0∞ (Ω) in W1,𝓗(Ω), that is,

W01,H(Ω)=C0∞(Ω)¯W1,H(Ω).

Besides, from condition H(Îź) and Colasuonno-Squassina [12, Proposition 2.18] we can see that

∥u∥1,H,0=∥∇u∥Hfor all u∈W01,H(Ω)

is an equivalent norm on W01,H (Ί). Now we are able to adapt (2.1) in terms of W01,H (Ί)-norm as follows

min∥u∥1,H,0p,∥u∥1,H,0q≤∥∇u∥pp+∥∇u∥q,μq≤max∥u∥1,H,0p,∥u∥1,H,0q (2.2)

for all u ∈ W01,H (Ω). Since both spaces W1,𝓗(Ω) and W01,H (Ω) are uniformly convex, we know that they are reflexive Banach spaces.

Furthermore, we have the following compact embedding

W01,H(Ω)↪Lr(Ω) (2.3)

for each 1 < r < p*, where p* is the critical exponent to p given by

p∗:=NpN−p, (2.4)

see Colasuonno-Squassina [12, Proposition 2.15].

Let us now consider the eigenvalue problem for the negative r-Laplacian with homogeneous Dirichlet boundary condition and 1 < r < ∞ which is defined by

−Δru=λ|u|r−2uin Ω,u=0on ∂Ω. (2.5)

From Lê [26] we know that the set σr being the set of all eigenvalues of (–Δr, W01,r (Ω)) has a smallest element λ1,r which is positive, isolated, simple and it can be variationally characterized through

λ1,r=inf∥∇u∥rr∥u∥rr:u∈W01,r(Ω),u≠0.

Now, let A: W01,H (Ω) → W01,H (Ω)* be the operator defined by

〈A(u),v〉H:=∫Ω|∇u|p−2∇u+μ(x)|∇u|q−2∇u⋅∇vdx, (2.6)

for u, v ∈ W01,H (Ω), where 〈 ⋅, ⋅〉𝓗 is the duality pairing between W01,H (Ω) and its dual space W01,H (Ω)*.

The properties of the operator A: W01,H (Ω) → W01,H (Ω)* can be summarized as follows, see Liu-Dai [27].

Proposition 2.1

The operator A defined by (2.6) is bounded, continuous, monotone (hence maximal monotone) and of type (S+).

Throughout the paper the symbols “⇀” and “→” stand for the weak and the strong convergence, respectively. Let (V, ∥⋅∥V) be a Banach space with its dual V* and denote by 〈⋅, ⋅〉 the duality pairing between V* and V. We end this section by recalling the following definition, see, for example, Papageorgiou-Winkert [38, Definition 6.7.4].

Definition 2.2

Let (X, τ) be a Hausdorff topological space and let {An} ⊂ 2X be a sequence of sets. We define the τ-Kuratowski lower limit of the sets An by

τ-lim infn→∞An:=x∈X∣x=τ-limn→∞xn,xn∈Anfor all n≥1,

and the τ-Kuratowski upper limit of the sets An

τ-lim supn→∞An:=x∈X∣x=τ-limk→∞xnk,xnk∈Ank,n1<n2<…<nk<….

If

A=τ-lim infn→∞An=τ-lim supn→∞An,

then A is called τ-Kuratowski limit of the sets An.

3 Main results

We assume the following hypotheses on the data of problem (1.2).

  1. The multivalued convection mapping f: Ω × ℝ × ℝN → 2ℝ has nonempty, compact and convex values such that

    1. the multivalued mapping x ↦ f(x, s, ξ) has a measurable selection for all (s, ξ) ∈ ℝ × ℝN;

    2. the multivalued mapping (s, ξ) ↦ f(x, s, ξ) is upper semicontinuous for almost all (a. a.) x ∈ Ω;

    3. there exists α ∈ Lq1q1−1 (Ω) and a1, a2 ≥ 0 such that

      |η|≤a1|ξ|pq1−1q1+a2|s|q1−1+α(x)

      for all η ∈ f(x, s, ξ), for a. a. x ∈ Ω, all s ∈ ℝ and all ξ ∈ ℝN, where 1 < q1 < p* with the critical exponent p* given in (2.4);

    4. there exist w ∈ L+1 (Ω) and b1, b2 ≥ 0 such that

      b1+b2λ1,p−1<1,

      and

      ηs≤b1|ξ|p+b2|s|p+w(x)

      for all η ∈ f(x, s, ξ), for a. a. x ∈ Ω, all s ∈ ℝ and all ξ ∈ ℝN, where λ1,p is the first eigenvalue of the Dirichlet eigenvalue problem for the p-Laplacian, see (2.5).

  2. Φ: Ω → [0, ∞) is such that Φ ∈ Lq1′ (Ω).

  3. {ρn} is a sequence with ρn > 0 for each n ∈ ℕ such that ρn → 0 as n → ∞.

Let K be a subset of W01,H (Ί) defined by

K:=u∈W01,H(Ω)|u(x)≤Φ(x)  for a.a. x∈Ω. (3.1)

Remark 3.1

  1. The set K is a nonempty, closed and convex subset of W01,H (Ί).

  2. From assumption H(Ό) we see that 0 ∈ K.

The weak solutions for problems (1.1) and (1.2) are understood in the following way.

Definition 3.2

  1. We say that u ∈ K is a weak solution of problem (1.1) if there exists η ∈ Lq1q1−1 (Ω) such that η(x) ∈ f(x, u(x), ∇ u(x)) for a. a. x ∈ Ω and

    ∫Ω|∇u|p−2∇u⋅∇(v−u)+μ(x)|∇u|q−2∇u⋅∇(v−u)dx=∫Ωη(x)(v−u)dx

    for all v ∈ K, where K is given by (3.1).

  2. We say that u ∈ W01,H (Ω) is a weak solution of problem (1.2) if there exists η ∈ Lq1q1−1 (Ω) such that η(x) ∈ f(x, u(x), ∇ u(x)) for a. a. x ∈ Ω and

    ∫Ω|∇u|p−2∇u+μ(x)|∇u|q−2∇u⋅∇vdx+1ρn∫Ωu(x)−Φ(x)+v(x)dx=∫Ωη(x)v(x)dx

    for all v ∈ W01,H (Ί).

It is straightforward, to prove the following lemma.

Lemma 3.3

If hypothesis H(Φ) holds, then the function B: Lq1(Ω) → Lq1′ (Ω) given by

〈Bu,v〉q1=∫Ωu(x)−Φ(x)+v(x)dx for all u,v∈Lq1(Ω), (3.2)

is bounded, demicontinuous and monotone, where 〈⋅, ⋅, 〉q1 denotes the duality pairing between Lq1(Ω) and its dual space Lq1′ (Ω).

Now, we can state the main result of this paper.

Theorem 3.4

If hypotheses H(Îź), H(f), H(ÎŚ), and H(0) hold, then

  1. for each n ∈ ℕ, the set 𝓢n of solutions to problem (1.2) is nonempty, bounded and closed;

  2. it holds

    ∅≠w-lim supn→∞Sn=s-lim supn→∞Sn⊂S;
  3. for each u ∈ s- lim supn→∞ 𝓢n and any sequence {u͠n} with

    u~n∈T(Sn,u) for each n∈N,

    there exists a subsequence of {u͠n} converging strongly to u in W01,H (Ω), where the set 𝓣(𝓢n, u) is defined by

    T(Sn,u):={u~∈Sn∣∥u−u~∥1,H,0≤∥u−v∥1,H,0 for all v∈Sn}.

Proof

  1. Let i: W01,H (Ω) → Lq1(Ω) be the embedding operator from W01,H (Ω) to Lq1(Ω) with its adjoint operator i*: Lq1′ (Ω) → W01,H (Ω)*. Since 1 < q1 < p* the embedding operator i is compact and so i* as well. From hypotheses H(f)(i) and (iii), we see that the Nemytskij operator N͠f: W01,H (Ω) ⊂ Lq1(Ω) → 2Lq1′(Ω) associated to the multivalued mapping f given by

    N~f(u):=η∈Lq1′(Ω)|η(x)∈f(x,u(x),∇u(x))  for a.a. x∈Ω

    for all u ∈ W01,H (Ί) is well-defined (see the proof of Proposition 3 in Papageorgiou-Vetro-Vetro [36]). The convexity and closedness of the values of f ensure that N͠f has closed and convex values as well. Moreover, by hypothesis H(f)(iv) we have

    ∥η∥q1′q1′=∫Ω|η(x)|q1′dx≤∫Ωa1|∇u(x)|pq1′+a2|u(x)|q1−1+α(x)q1′dx≤M0∫Ω|∇u(x)|p+|u(x)|q1+α(x)q1′dx=M0∥∇u∥pp+∥u∥q1q1+∥α∥q1′q1′. (3.3)

    Notice that the embeddings W01,H (Ω) ⊂ W01,p (Ω) ⊂ Lq1(Ω) are both continuous, so, N͠f(u) is bounded in Lq1′ (Ω) for each u ∈ W01,H (Ω).

    It is easy to see that u ∈ W01,H (Ί) is a weak solution of problem (1.2) (see Definition 3.2(b)), if and only if u solves the following inclusion:

    Find u ∈ W01,H (Ί) and Ρ ∈ N͠f(u) such that

    A(u)+1ρni∗B(u)−i∗N~f(u)∋0,

    where A: W01,H (Ω) → W01,H (Ω)* and B: Lq1(Ω) → Lq1′ (Ω) are given by (2.6) and (3.2), respectively.

    Then, using the same arguments as in the proof of Zeng-Gasiński-Winkert-Bai [44, Theorem 3.3], we can conclude that for each n ∈ ℕ, the set 𝓢n of solutions to problem (1.2) is nonempty, bounded and closed.

  2. First, we prove that the set w- lim supn→∞ 𝓢n is nonempty. Indeed, we have the following claims.

    1. The set ⋃n∈N 𝓢n is uniformly bounded in W01,H (Ω).

      Arguing by contradiction, suppose that ⋃n∈N 𝓢n is unbounded. Without any loss of generality (passing to a subsequence if necessary), we may assume that there exists a sequence {un} ⊂ W01,H (Ω) with un ∈ 𝓢n for each n ∈ ℕ such that

      ∥un∥1,H,0→∞ as n→∞.

      Hence, for each n ∈ ℕ, we are able to find ηn ∈ N͠f(un) such that

      ∫Ω|∇un|p−2∇un+μ(x)|∇un|q−2∇un⋅∇vdx+1ρn∫Ωun(x)−Φ(x)+v(x)dx=∫Ωηn(x)v(x)dx

      for all v ∈ W01,H (Ί). Inserting v = un into the inequality above, we get

      ∫Ω|∇un|p−2∇un+μ(x)|∇un|q−2∇un⋅∇undx−∫Ωηn(x)un(x)dx=−1ρn∫Ωun(x)−Φ(x)+un(x)dx.

      By the nonnegativity of Φ and the monotonicity of the function s ↦ s+, we have

      ∫Ω|∇un|p−2∇un+μ(x)|∇un|q−2∇un⋅∇undx−∫Ωηn(x)un(x)dx=−1ρn∫Ωun(x)−Φ(x)+−0−Φ(x)+un(x)dx≤0,

      thus

      ∥∇un∥pp+∥∇un∥q,μq−∫Ωηn(x)un(x)dx≤0. (3.4)

      However, by hypothesis H(f)(iv), we have

      ∫Ωηn(x)un(x)dx≤b1∥∇un∥pp+b2∥un∥pp+∥w∥1. (3.5)

      Applying (3.5) in (3.4), using the continuity of the embedding W01,H (Ω) ⊆ W01,p (Ω) as well as the estimate

      ∥u∥pp≤λ1,p−1∥∇u∥ppfor all u∈W01,p(Ω),

      we get

      0≥∥∇un∥pp+∥∇un∥q,μq−∫Ωηn(x)un(x)dx≥∥∇un∥pp+∥∇un∥q,μq−b1∥∇un∥pp−b2∥un∥pp−∥w∥1≥1−b1−b2λ1,p−1∥∇un∥pp+∥∇un∥q,μq−∥w∥1≥1−b1−b2λ1,p−1∥∇un∥pp+∥∇un∥q,μq−∥w∥1≥1−b1−b2λ1,p−1min∥un∥1,H,0p,∥un∥1,H,0q−∥w∥1,

      where the last inequality is obtained by (2.2). Since 1 < p < q < N and b1 + b2λ1,p−1 < 1, we can take R0 > 0 large enough such that for all R ≥ R0 it holds

      1−b1−b2λ1,p−1minRp,Rq−∥w∥1>0.

      Therefore, we are able to find N0 > 0 large enough such that ∥un∥1,𝓗,0 > R0 for all n ≥ N0 and

      0≥1−b1−b2λ1,p−1min∥un∥1,H,0p,∥un∥1,H,0q−∥w∥1>0

      for all n ≥ N0. This gives a contradiction, so Claim 1 is proved.

      Let {un} ⊂ W01,H (Ω) with un ∈ 𝓢n for each n ∈ ℕ be an arbitrary sequence. Claim 1 indicates that {un} is bounded in W01,H (Ω). Then, we may assume that along a relabeled subsequence we have

      un⇀u as n→∞ (3.6)

      for some u ∈ W01,H (Ω). This guarantees that the set w- lim supn→∞ 𝓢n is nonempty.

      Next, we are going to demonstrate that w- lim supn→∞ 𝓢n is a subset of 𝓢. Let u ∈ w- lim supn→∞ 𝓢n be arbitrary. Without loss of generality, we may suppose that there exists a subsequence {un} ⊂ W01,H (Ω) with un ∈ 𝓢n for all n ∈ ℕ, satisfying (3.6). Our goal is to prove that u ∈ 𝓢.

    2. u(x) ≤ Φ(x) for a.a. x ∈ Ω.

      For every n ∈ ℕ, we have

      1ρn∫Ωun(x)−Φ(x)+v(x)dx=〈Aun,−v〉H+∫Ωηn(x)v(x)dx. (3.7)

      It follows from Hölder’s inequality and (3.3) that

      ∫Ωηn(x)v(x)dx≤M01q1′∥∇un∥pp+∥un∥q1q1+∥α∥q1′q1′1q1′∥v∥q1. (3.8)

      Putting (3.8) into (3.7), employing the boundedness of A (see Proposition 2.1), the convergence (3.6), and the embedding (2.3), we have

      1ρn∫Ωun(x)−Φ(x)+v(x)dx≤∥Aun∥1,H,0∥v∥1,H,0+M01q1′∥∇un∥pp+∥un∥q1q1+∥α∥q1′q1′1q1′∥v∥q1≤M1∥v∥1,H,0

      for some M1 > 0, where M1 > 0 is independent of n, that is

      ∫Ωun(x)−Φ(x)+v(x)dx≤ρnM1∥v∥1,H,0

      for all v ∈ W01,H (Ί). Passing to the limit in the inequality above, using convergence (3.6), the compact embedding (2.3), and the Lebesgue Dominated Convergence Theorem, we conclude that

      ∫Ωu(x)−Φ(x)+v(x)dx=∫Ωlimn→∞un(x)−Φ(x)+v(x)dx=limn→∞∫Ωun(x)−Φ(x)+v(x)dx≤limn→∞ρnM1∥v∥1,H,0=0

      for all v ∈ W01,H (Ω). Therefore, we have (u(x) – Φ(x))+ = 0 for a.a. x ∈ Ω, thus, u(x) ≤ Φ(x) for a.a. x ∈ Ω.

    3. u ∈ 𝓢.

      For each n ∈ ℕ, we have

      〈Aun,un−v〉H=1ρn∫Ωun(x)−Φ(x)+(v(x)−un(x))dx+∫Ωηn(x)(un(x)−v(x))dx

      for all v ∈ W01,H (Ω). The latter combined with the monotonicity of s ↦ s+ gives

      〈Aun,un−v〉H≤1ρn∫Ωv(x)−Φ(x)+(v(x)−un(x))dx+∫Ωηn(x)(un(x)−v(x))dx

      for all v ∈ W01,H (Ί). Hence,

      〈Aun,un−v〉H−∫Ωηn(x)(un(x)−v(x))dx≤0 (3.9)

      for all v ∈ K, where K is defined in (3.1).

      Claim 2 indicates that u ∈ K, so, we put v = u in (3.9) to obtain

      〈Aun,un−u〉H−∫Ωηn(x)(un(x)−u(x))dx≤0,

      that is,

      lim supn→∞〈Aun−i∗ηn,un−u〉H≤0.

      It follows from the proof of Theorem 3.3 in Zeng-Gasiński-Winkert-Bai [44] that the multivalued mapping 𝓐 = A – i*N͠f is pseudomonotone. So, for each v ∈ K, there exists u* ∈ 𝓐u such that

      lim infn→∞〈Aun−i∗ηn,un−v〉H≥〈u∗(v),un−v〉.

      This means that for each v ∈ K, there is an element Ρ(v) ∈ N͠f(u) satisfying

      u∗(v)=Au−i∗η(v).

      For each v ∈ K, passing to the lower limit as n → ∞ in inequality (3.9), we are able to find an element η(v) ∈ N͠f(u) such that

      〈Au,v−u〉H−∫Ωη(v)(x)(v(x)−u(x))dx≥0. (3.10)

      We shall prove that u ∈ K is a weak solution to problem (1.1), namely, there exists an element Ρ* ∈ N͠f(u), which is independent of v, such that

      〈Au,v−u〉H−∫Ωη∗(x)(v(x)−u(x))dx≥0 (3.11)

      for all v ∈ K. Arguing by contradiction, suppose that for each Ρ ∈ N͠f(u), there is v ∈ K such that

      〈Au,v−u〉H−∫Ωη(x)(v(x)−u(x))dx<0.

      For any v ∈ K, let us consider the set Rv ⊂ N͠f(u) defined by

      Rv:=η∈N~f(u)∣〈Au,v−u〉H−∫Ωη(x)(v(x)−u(x))dx<0

      for all v ∈ K. We now assert that for each v ∈ K, the set Rv is weakly open. Let {ηn} ⊂ Rvc be such that ηn ⇀ η for some η ∈ Lq1′ (Ω) as n → ∞, where Rvc denotes the complement of Rv. Hence,

      〈Au,v−u〉H−∫Ωηn(x)(v(x)−u(x))dx≥0

      for all n ∈ ℕ. Passing to the limit in the inequality above, we obtain that η ∈ Rvc . Therefore, for every v ∈ K, the set Rv is weakly open in Lq1′ (Ω). Besides, we observe that {Rv}v∈K is an open covering of N͠f(u). The latter coupled with the facts that Lq1(Ω) is reflexive and N͠f(u) is weakly compact and convex in Lq1′ (Ω), ensures that {Rv}v∈K has a finite sub-covering of N͠f(u), let us say {Rv1, Rv2, …, Rvn} for some points {v1, v2, …, vn} ⊆ K. Let κ1, κ2, …, κn be a partition of unity for N͠f(u), where for each i = 1, 2, …, n, κi : N͠f(u) → [0, 1] is a weakly continuous function such that ∑i=1n κi(η) = 1 for all η ∈ N͠f(u), see, for example, Granas-Dugundji [23, Lemma 7.3].

      Also, we introduce a function 𝓜: N͠f(u) → W01,H (Ω) defined by

      M(η)=∑i=1nκi(η)vi for all η∈N~f(u).

      Obviously, the function 𝓜 is also weakly continuous due to the weak continuity of κi for i = 1, 2, …, n. For any η ∈ N͠f(u), we have

      〈Au−i∗η,M(η)−u〉H=〈Au−i∗η,∑i=1nκi(η)vi−u〉H=∑i=1nκi(η)〈Au−i∗η,vi−u〉H<0 (3.12)

      for all Ρ ∈ N͠f(u), where the last inequality is obtained by the use of Lemma 7.3(ii) of Granas-Dugundji [23].

      Let us define two multivalued functions Λ: K → 2N͠f(u) and Ψ: N͠f(u) → 2N͠f(u) by

      Λ(v):=η∈N~f(u)∣〈Au,v−u〉H−∫Ωη(x)(v(x)−u(x))dx≥0

      for all v ∈ K, and

      Ψ(η):=Λ(M(η)) for all η∈N~f(u).

      Then, Ψ has nonempty, weakly compact and convex values (by (3.10) and because N͠f(u) is bounded closed and convex in Lq1′ (Ω)) and Λ is upper semicontinuous from the normal topology of K to weak topology of Lq1′ (Ω). From Migórski-Ochal-Sofonea [31, Proposition 3.8], it is enough to verify that for each weakly closed set D in Lq1′ (Ω), the set

      Λ−(D):=v∈K∣Λ(v)∩D≠∅

      is closed in W01,H (Ω). Let {vn} ⊂ Λ–(D) be a sequence such that vn → v as n → ∞. Then, for each n ∈ ℕ, we are able to find ηn ∈ N͠f(u) satisfying

      〈Au,vn−u〉H−∫Ωηn(x)(vn(x)−u(x))dx≥0. (3.13)

      From the weak compactness of N͠f(u), without any loss of generality, we may suppose that ηn ⇀ η in Lq1′ (Ω), as n → ∞, for some η ∈ N͠f(u). Passing to the upper limit as n → ∞ for (3.13), we have

      〈Au,v−u〉H−∫Ωη(x)(v(x)−u(x))dx≥0,

      that is, η ∈ Λ(v). But, the weak closedness of D implies that η ∈ D. Therefore, η ∈ Λ(v) ∩ D and so v ∈ Λ–(D). Applying Migórski-Ochal-Sofonea [31, Proposition 3.8] derives that Λ is strongly-weakly upper semicontinuous. On the other hand, the continuity of 𝓜 and Theorem 1.2.8 of Kamenskii-Obukhovskii-Zecca [24] imply that Ψ is also strongly-weakly upper semicontinuous.

      We are now in a position to employ Tychonov fixed point principle, (see, for example, Granas-Dugundji [23, Theorem 8.6]) for function Ψ, to conclude that there exists Ρ ∈ N͠f(u) such that

      〈Au,M(η)−u〉H−∫Ωη(x)(M(η)(x)−u(x))dx≥0.

      This leads to a contraction with (3.12). Consequently, we infer that u ∈ K solves problem (1.1) as well, that means, there exists Ρ ∈ N͠f(u), which is independent of v, such that (3.11) holds.

      Consequently, we conclude that ∅ ≠ w- lim supn→∞ 𝓢n ⊂ 𝓢.

    4. It holds w- lim supn→∞ 𝓢n = s- lim supn→∞ 𝓢n.

      Since s- lim supn→∞ 𝓢n ⊂ w- lim supn→∞ 𝓢n, it is enough to verify the condition w- lim supn→∞ 𝓢n ⊂ s- lim supn→∞ 𝓢n. Let u ∈ w- lim supn→∞ 𝓢n be arbitrary. Without any loss of generality, there exists a sequence, still denoted by {un} with un ∈ 𝓢n such that un ⇀ u as n → ∞. We claim that un → u as n → ∞. For each n ∈ ℕ, it holds

      〈Aun,un−v〉H=−∫Ω(un(x)−Φ(x))+(un(x)−v(x))dx+∫Ωηn(x)(un(x)−v(x))dx

      for some ηn ∈ N͠f(un) and for all v ∈ W01,H (Ω). Inserting v = u into the above inequality and passing to the upper limit as n → ∞ for the resulting inequality, we can use the compact embedding (2.3) to get

      lim supn→∞〈Aun,un−u〉H≤0.

      The latter combined with the convergence un ⇀ u as n → ∞ and the fact that A is of type (S+) (see Proposition 2.1) implies that un → u as n → ∞. This means that u ∈ s- lim supn→∞ 𝓢n. Therefore s- lim supn→∞ 𝓢n = w- lim supn→∞ 𝓢n.

  3. Let u ∈ s- lim supn→∞ 𝓢n be arbitrary. Since 𝓢n is nonempty, bounded and closed, so, the set 𝓣(𝓢n, u) is nonempty. Let {u͠n} be any sequence such that

    u~n∈T(Sn,u) for each n∈N.

    It follows from Claim 1 that the sequence {uÍ n} is bounded. So, passing to a subsequence, we may assume, that

    u~n⇀u~ as n→∞

    for some u͠ ∈ W01,H (Ω). Thus, using the same argument as the proof of Claim 2, we get that u͠ ∈ K. Then, for each n ∈ ℕ, we have

    〈Au~n,u~n−v〉H=1ρn∫Ωu~n(x)−Φ(x)+(v(x)−u~n(x))dx+∫Ωηn(x)(u~n(x)−v(x))dx

    for all v ∈ W01,H (Ω). Proceeding in the same way as in the proof of Claim 3, we conclude that u͠ is a solution to problem (1.1) as well. Consequently, the desired conclusion is proved.□

Acknowledgment

The authors wish to thank the three knowledgeable referees for their useful remarks in order to improve the paper.

Project supported by the NNSF of China Grant No. 12001478, H2020-MSCA-RISE-2018 Research and Innovation Staff Exchange Scheme Fellowship within the Project No. 823731 CONMECH, and National Science Center of Poland under Preludium Project No. 2017/25/N/ST1/00611. It is also supported by the Startup Project of Doctor Scientific Research of Yulin Normal University No. G2020ZK07, and International Project co-financed by the Ministry of Science and Higher Education of Republic of Poland under Grant No. 3792/GGPJ/H2020/2017/0.

References

[1] A. Bahrouni, V. D. Rădulescu, D. D. Repovš, Double phase transonic flow problems with variable growth: nonlinear patterns and stationary waves, Nonlinearity 32 (2019), no. 7, 2481–2495.10.1088/1361-6544/ab0b03Search in Google Scholar

[2] A. Bahrouni, V. D. Rădulescu, P. Winkert, A critical point theorem for perturbed functionals and low perturbations of differential and nonlocal systems, Adv. Nonlinear Stud. 20 (2020), no. 3, 663–674.10.1515/ans-2020-2095Search in Google Scholar

[3] A. Bahrouni, V. D. Rădulescu, P. Winkert, Double phase problems with variable growth and convection for the Baouendi-Grushin operator, Z. Angew. Math. Phys., 71 (2020), no. 6, 183.10.1007/s00033-020-01412-7Search in Google Scholar

[4] P. Baroni, M. Colombo, G. Mingione, Harnack inequalities for double phase functionals, Nonlinear Anal. 121 (2015), 206–222.10.1016/j.na.2014.11.001Search in Google Scholar

[5] P. Baroni, M. Colombo, G. Mingione, Non-autonomous functionals, borderline cases and related function classes, St. Petersburg Math. J. 27 (2016), 347–379.10.1090/spmj/1392Search in Google Scholar

[6] P. Baroni, M. Colombo, G. Mingione, Regularity for general functionals with double phase, Calc. Var. Partial Differential Equations 57 (2018), no. 2, Art. 62, 48 pp.10.1007/s00526-018-1332-zSearch in Google Scholar

[7] P. Baroni, T. Kuusi, G. Mingione, Borderline gradient continuity of minima, J. Fixed Point Theory Appl. 15 (2014), no. 2, 537–575.10.1007/s11784-014-0188-xSearch in Google Scholar

[8] V. Benci, P. D’Avenia, D. Fortunato, L. Pisani, Solitons in several space dimensions: Derrick’s problem and infinitely many solutions, Arch. Ration. Mech. Anal. 154 (2000), no. 4, 297–324.10.1007/s002050000101Search in Google Scholar

[9] L. Cherfils, Y. Il′ yasov, On the stationary solutions of generalized reaction diffusion equations with p&q-Laplacian, Commun. Pure Appl. Anal. 4 (2005), no. 1, 9–22.10.3934/cpaa.2005.4.9Search in Google Scholar

[10] S. Carl, V. K. Le, D. Motreanu, “Nonsmooth Variational Problems and Their Inequalities”, Springer, New York, 2007.10.1007/978-0-387-46252-3Search in Google Scholar

[11] M. Cencelj, V. D. Rădulescu, D. D. Repovš, Double phase problems with variable growth, Nonlinear Anal. 177 (2018), part A, 270–287.10.1016/j.na.2018.03.016Search in Google Scholar

[12] F. Colasuonno, M. Squassina, Eigenvalues for double phase variational integrals, Ann. Mat. Pura Appl. (4) 195 (2016), no. 6, 1917–1959.10.1007/s10231-015-0542-7Search in Google Scholar

[13] M. Colombo, G. Mingione, Bounded minimisers of double phase variational integrals, Arch. Ration. Mech. Anal. 218 (2015), no. 1, 219–273.10.1007/s00205-015-0859-9Search in Google Scholar

[14] M. Colombo, G. Mingione, Regularity for double phase variational problems, Arch. Ration. Mech. Anal. 215 (2015), no. 2, 443–496.10.1007/s00205-014-0785-2Search in Google Scholar

[15] G. Cupini, P. Marcellini, E. Mascolo, Local boundedness of minimizers with limit growth conditions, J. Optim. Theory Appl. 166 (2015), no. 1, 1–22.10.1007/s10957-015-0722-zSearch in Google Scholar

[16] L. Gasiński, N.S. Papageorgiou, Constant sign and nodal solutions for superlinear double phase problems, Adv. Calc. Var., 10.1515/acv-2019-0040.Search in Google Scholar

[17] L. Gasiński, N.S. Papageorgiou, Positive solutions for nonlinear elliptic problems with dependence on the gradient, J. Differential Equations 263 (2017), 1451–1476.10.1016/j.jde.2017.03.021Search in Google Scholar

[18] L. Gasiński, N.S. Papageorgiou, “Nonsmooth critical point theory and nonlinear boundary value problems”, Chapman & Hall/CRC, Boca Raton, FL, 2005.10.1201/9780367801632Search in Google Scholar

[19] L. Gasiński, P. Winkert, Constant sign solutions for double phase problems with superlinear nonlinearity, Nonlinear Anal. 195 (2020), 111739.10.1016/j.na.2019.111739Search in Google Scholar

[20] 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

[21] L. Gasiński, P. Winkert, Sign changing solution for a double phase problem with nonlinear boundary condition via the Nehari manifold, https://arxiv.org/abs/2003.13241Search in Google Scholar

[22] F. H. Clarke, “Optimization and Nonsmooth Analysis”, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 1990.Search in Google Scholar

[23] A. Granas, J. Dugundji, “Fixed Point Theory”, Springer-Verlag, New York, 2003.10.1007/978-0-387-21593-8Search in Google Scholar

[24] M. Kamenskii, V Obukhovskii, P. Zecca, “Condensing Multivalued Maps and Semilinear Differential Inclusions in Banach Spaces”, de Gruyter, Berlin, 2001.10.1515/9783110870893Search in Google Scholar

[25] V. K. Le, A range and existence theorem for pseudomonotone perturbations of maximal monotone operators, Proc. Amer. Math. Soc. 139 (2011), no. 5, 1645–1658.10.1090/S0002-9939-2010-10594-4Search in Google Scholar

[26] A. Lê, Eigenvalue problems for the p-Laplacian, Nonlinear Anal. 64 (2006), no. 5, 1057–1099.10.1016/j.na.2005.05.056Search in Google Scholar

[27] W. Liu, G. Dai, Existence and multiplicity results for double phase problem, J. Differential Equations 265 (2018), no. 9, 4311–4334.10.1016/j.jde.2018.06.006Search in Google Scholar

[28] P. Marcellini, The stored-energy for some discontinuous deformations in nonlinear elasticity, in “Partial differential equations and the calculus of variations, Vol. II”, vol. 2, 767–786, Birkhäuser Boston, Boston, 1989.10.1007/978-1-4684-9196-8_32Search in Google Scholar

[29] P. Marcellini, Regularity and existence of solutions of elliptic equations with p, q-growth conditions, J. Differential Equations 90 (1991), no. 1, 1–30.10.1016/0022-0396(91)90158-6Search in Google Scholar

[30] G. Marino, P. Winkert, Existence and uniqueness of elliptic systems with double phase operators and convection terms, J. Math. Anal. Appl. 492 (2020), 124423, 13 pp.10.1016/j.jmaa.2020.124423Search in Google Scholar

[31] S. Migórski, A. Ochal, M. Sofonea, “Nonlinear Inclusions and Hemivariational Inequalities”, Springer, New York, 2013.10.1007/978-1-4614-4232-5Search in Google Scholar

[32] 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

[33] 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, Art. 108, 21 pp.10.1007/s00033-018-1001-2Search in Google Scholar

[34] N. S. Papageorgiou, C. Vetro, F. Vetro, Continuous spectrum for a two phase eigenvalue problem with an indefinite and unbounded potential, J. Differential Equations 268 (2020), no. 8, 4102–4118.10.1016/j.jde.2019.10.026Search in Google Scholar

[35] N. S. Papageorgiou, C. Vetro, F. Vetro, Multiple solutions for parametric double phase Dirichlet problems, Commun. Contemp. Math., 10.1142/S0219199720500066.Search in Google Scholar

[36] N. S. Papageorgiou, C. Vetro, F. Vetro, Nonlinear multivalued Duffing systems, J. Math. Anal. Appl. 468 (2018), no. 1, 376–390.10.1016/j.jmaa.2018.08.024Search in Google Scholar

[37] N. S. Papageorgiou, C. Vetro, F. Vetro, Relaxation for a Class of Control Systems with Unilateral Constraints, Acta Appl. Math. 167 (2020), no. 1, 99–115.10.1007/s10440-019-00270-4Search in Google Scholar

[38] N. S. Papageorgiou, P. Winkert, “Applied Nonlinear Functional Analysis. An Introduction”, De Gruyter, Berlin, 2018.10.1515/9783110532982Search in Google Scholar

[39] K. Perera, M. Squassina, Existence results for double-phase problems via Morse theory, Commun. Contemp. Math. 20 (2018), no. 2, 1750023, 14 pp.10.1142/S0219199717500237Search in Google Scholar

[40] V. D. Rădulescu, Isotropic and anistropic double-phase problems: old and new, Opuscula Math. 39 (2019), no. 2, 259–279.10.7494/OpMath.2019.39.2.259Search in Google Scholar

[41] C. Vetro, Parametric and nonparametric A-Laplace problems: Existence of solutions and asymptotic analysis, Asymptot. Anal., 10.3233/ASY-201612.Search in Google Scholar

[42] C. Vetro, F. Vetro, On problems driven by the (p(⋅), q(⋅))-Laplace operator, Mediterr. J. Math. 17 (2020), no. 1, 1–11.Search in Google Scholar

[43] S. D. Zeng, Y. R. Bai, L. Gasiński, P. Winkert, Existence results for double phase implicit obstacle problems involving multivalued operators, Calc. Var. Partial Differential Equations, 59:5 (2020), pages 18.10.1007/s00526-020-01841-2Search in Google Scholar

[44] S. D. Zeng, L. Gasiński, P. Winkert, Y. R. Bai, Existence of solutions for double phase obstacle problems with multivalued convection term, J. Math. Anal. Appl., 10.1016/j.jmaa.2020.123997.Search in Google Scholar

[45] Q. Zhang, V. D. Rădulescu, Double phase anisotropic variational problems and combined effects of reaction and absorption terms, J. Math. Pures Appl. (9) 118 (2018), 159–203.10.1016/j.matpur.2018.06.015Search in Google Scholar

[46] V. V. Zhikov, Averaging of functionals of the calculus of variations and elasticity theory, Izv. Akad. Nauk SSSR Ser. Mat. 50 (1986), no. 4, 675–710.10.1070/IM1987v029n01ABEH000958Search in Google Scholar

[47] V. V. Zhikov, On variational problems and nonlinear elliptic equations with nonstandard growth conditions, J. Math. Sci. 173 (2011), no. 5, 463–570.10.1007/s10958-011-0260-7Search in Google Scholar

Received: 2019-06-10
Accepted: 2020-09-27
Published Online: 2020-11-26

Š 2021 S. Zeng et al., published by De Gruyter

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

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