Next Article in Journal
On Singular Distributions With Statistical Structure
Previous Article in Journal
Stochastic Aspects of Proportional Vitalities Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Constrained Variational-Hemivariational Inequalities on Nonconvex Star-Shaped Sets

by
Stanisław Migórski
1,2,* and
Long Fengzhen
3
1
College of Sciences, Beibu Gulf University, Qinzhou 535000, China
2
Faculty of Mathematics and Computer Science, Jagiellonian University in Krakow, ul. Lojasiewicza 6, 30348 Krakow, Poland
3
School of Science, Nanjing University of Science and Technology, Nanjing 210094, China
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(10), 1824; https://doi.org/10.3390/math8101824
Submission received: 24 September 2020 / Revised: 13 October 2020 / Accepted: 14 October 2020 / Published: 17 October 2020
(This article belongs to the Section Difference and Differential Equations)

Abstract

:
In this paper, we study a class of constrained variational-hemivariational inequality problems with nonconvex sets which are star-shaped with respect to a certain ball in a reflexive Banach space. The inequality is a fully nonconvex counterpart of the variational-hemivariational inequality of elliptic type since it contains both, a convex potential and a locally Lipschitz one. Two new results on the existence of a solution are proved by a penalty method applied to a variational-hemivariational inequality penalized by the generalized directional derivative of the distance function of the constraint set. In the first existence theorem, the strong monotonicity of the governing operator and a relaxed monotonicity condition of the Clarke subgradient are assumed. In the second existence result, these two hypotheses are relaxed and a suitable hypothesis on the upper semicontinuity of the operator is adopted. In both results, the penalized problems are solved by using the Knaster, Kuratowski, and Mazurkiewicz (KKM) lemma. For a suffciently small penalty parameter, the solution to the penalized problem solves also the original one. Finally, we work out an example on the interior and boundary semipermeability problem that ilustrate the applicability of our results.

1. Introduction

In this paper, we are initially motivated by the investigation of the class of variational-hemivariational inequalities considered in [1]. Let V be a reflexive Banach space. Consider an operator A : V V * , functions φ : V R , j : V R and a set K V . The constrained variational-hemivariational inequality studied in the paper reads as follows: find an element u K such that
A u , v u + φ ( v ) φ ( u ) + j 0 ( u ; v u ) f , v u for all v K .
Particular forms of problem (1) contain various formulations investigated in the literature: the elliptic variational inequalities of the first and second kind, see [2,3,4,5], the elliptic hemivariational inequalities, see [6,7,8,9], and the elliptic equations, see [6,10,11]. Moreover, the quasi-variational inequalities corresponding to problem (1) and its variants can be treated by a fixed point technique, see [1,5]. In all aforementioned papers, the usual hypotheses for existence (and uniqueness) of a solution involve the function φ , which is supposed to be convex, the function j which is locally Lipschitz and in general nonconvex, the operator A pseudomonotone and strongly monotone, and a nonempty, closed and convex set of constraints K.
In the current paper, we treat the counterpart of problem (1) where the set K represents a set of admissible constraints which is star-shaped with respect to a certain ball in V. Note that for this class of nonconvex sets, some particular versions have been studied earlier in Section 7.4 of [8] if j = 0 , in Section 7.3 of [8,12] if φ = 0 , and Section 7.2 of [8] when φ = j = 0 . The first novelty of our contribution is to study the class of fully variational-hemivariational inequalities involving nonconvex constraints. In contrast to contributions in [8,12,13] which are based on surjectivity methods for the multivalued pseudomonotone operators for the existence proof of the penalized problem, in this paper we employ another method based on the Knaster, Kuratowski, and Mazurkiewicz (KKM) theory for the penalized problem. The second novelty is to study the constrained variational-hemivariational inequality on star-shaped sets without hypothesis on the strong monotonicity of the operator and without the relaxed monotonicity condition of the generalized subgradient. To the best of our knowledge, it is a new approach in the examination of this class of variational-hemivariational inequalities. We will use the celebrated lemma by Fan which is a milestone in the KKM theory and it is sufficient for our purpose.
For related results on variational-hemivariational equalities on nonconvex star-shaped constraint sets, we refer to [12,13] for stationary problems, and to [14,15] for evolution problems. Numerous applications of variational-hemivariational inequalities to problems of nonsmooth contact mechanics, economics, etc. can be found in classical monographs [8,9,16], and in two recent books [17,18], and the references therein. A unified method, based on the hemivariational inequality formulation, to study contact problems of viscoelasticity is given in [19], the abstract elliptic variational-hemivariational inequalities in reflexive Banach spaces with applications can be found in [1], and the variational-hemivariational inqualities which model fluid flow in mechanics were treated in [20,21] and very recently in [22,23]. Other recent developments on variational methods in the study of existence and multiplicity of solutions, see [24,25].
The paper is structured as follows. Notations, basic definitions and preliminaries are provided in Section 2. Section 3 contains problem formulation with the first existence result, Theorem 1, whose proof is provided in Section 4. Finally, the second existence theorem is proved in Section 5.

2. Basic Material

In this part of the paper we recall the standard notations and definitions from [10,17,26,27].
Let ( X , · X ) be a Banach space. By ( X * , · X * ) we denote its dual space while the symbol · , · X * × X stands for the duality pairing between X * and X.
Nonlinear operators and the KKM lemma. Let X be a reflexive Banach space and T : X 2 X * be a multivalued mapping. A mapping T : X 2 X * is called bounded if the image of each bounded set in X remains in a bounded subset of X * . A mapping T : X 2 X * is called pseudomonotone, provided the following conditions are satisfied
(i)
T has nonempty, bounded, closed and convex values.
(ii)
T is upper semicontinuous (u.s.c.) from each finite dimensional subspace of X into X * equipped with its weak topology.
(iii)
if u n X , u n u weakly in X, u n * T u n and lim sup u n * , u n u X * × X 0 , then to each y X , there exists u * ( y ) T u such that u * ( y ) , u y X * × X lim inf u n * , u n y X * × X .
A mapping T : X 2 X * is called generalized pseudomonotone, if for each sequences { u n } X , u n u weakly in X, { u n * } X * , u n * u * weakly in X * , u n * T u n and lim sup u n * , u n u X * × X 0 , we have u * T u and u n * , u n X * × X u * , u X * × X .
The following relations concern the classes of pseudomonotone and generalized pseudomonotone mappings, see Propositions 1.3.65 and 1.3.68 in [10]. If T : X 2 X * is a pseudomonotone mapping, then it is also generalized pseudomonotone. If T : X 2 X * is a bounded, generalized pseudomonotone mapping with nonempty, closed and convex values, then T is pseudomonotone.
A single-valued mapping A : X X * is said to be pseudomonotone, provided it is bounded and if u n u weakly in X with lim sup A u n , u n u X * × X 0 imply A u , u v X * × X lim inf A u n , u n v X * × X for all v X . Equivalently, see Proposition 3.66 in [17], a single-valued mapping A is pseudomonotone, if and only if it is bounded and u n u weakly in X with lim sup A u n , u n u X * × X 0 imply lim A u n , u n u X * × X = 0 and A u n A u weakly in X * . A mapping A : X X * is demicontinuous, if it is continuous as a map from X to X * furnished with the weak topology.
We recall below the KKM lemma in a version stated by Fan in Lemma 1 of [28]. For various extensions to the Fan-KKM theorem, we refer to [29] and the references therein.
Lemma 1.
Let X be a subset of a Haussdorf topological vector space Y. For any x X , let a set F ( x ) in Y be given such that:
(a)
for every x X , F ( x ) is a closed set in Y,
(b)
convex hull of any finite set { x 1 , , x r } of X is contained in i = 1 r F ( x i ) ,
(c)
F ( x ) is a compact set at least for one x X .
Then x X F ( x ) .
The Clarke generalized subgradient and tangent cones. Let X be a Banach space, h : X R be a locally Lipschitz function, and x, v X . The Clarke generalized directional derivative of h at x in the direction v is given by
h 0 ( x ; v ) = lim sup y x , λ 0 h ( y + λ v ) h ( y ) λ .
The Clarke generalized subgradient of h at x is defined by
h ( x ) = { ζ X * h 0 ( x ; v ) ζ , v X * × X for all v X } .
It is well known that
h 0 ( x ; v ) = max { ζ , v ζ h ( x ) } .
Let B ( u 0 , ϱ ) be the closed ball in a normed space E with centre u 0 E and radius ϱ > 0 . A nonempty set K E is called star-shaped with respect to a ball B ( u 0 , ϱ ) , if t v + ( 1 t ) w K for all v K , w B ( u 0 , ϱ ) , t [ 0 , 1 ] . When a set K is star-shaped with respect to a closed ball, we just say that K is star-shaped. Next, we denote by d : E R the distance function of K defined by
d ( u ) = inf v K v u E for all u E .
The Clarke generalized directional derivative of the function d is well defined since d is a Lipschitz continuous. Recall that for a star-shaped set, the Clarke directional derivative of the function d enjoys the following discontinuity property, see Lemma 7.2, p. 224 in [8].
Lemma 2.
Let E be a reflexive Banach space and K E be a closed set which is star-shaped with respect to a ball B ( u 0 , ϱ ) with some u 0 K and ϱ > 0 . Then
d 0 ( u , u 0 u ) d ( u ) ϱ for all u K , d 0 ( u , u 0 u ) = 0 for all u K .
Finally, we shortly recall a material on tangent cones needed in what follows. Let K E be a set of a Banach space E and u K . The Bouligand (contingent) cone to K at the point u is defined by
K K ( u ) = { v E lim inf t 0 d ( u + t v ) t = 0 }
and the (Clarke) tangent cone to K at u is given by
T K ( u ) = { v E d 0 ( u ; v ) = 0 } = { v E d 0 ( u ; v ) 0 } .
It is well-known that T K ( u ) K K ( u ) . The set K is said to be regular at u K when T K ( u ) = K K ( u ) . We also know that if K is closed, convex and u K , then K is regular at u, see Theorem 10.39 in [30]. Further, it follows from Proposition 2.9 in [30] that if K is a convex set in E and u K , then K K ( u ) is a convex cone and
K u + K K ( u ) .
Equivalent definitions and properties of these and other cones can be found in [26,30] and Section 5.7 of [27].

3. Formulation of the Problem

In this section we consider the constrained problem in which the set of admissible elements is nonconvex. The main goal is to prove the existence of solution.
Let ( V , · V ) be a reflexive Banach space which is continuously and compactly embedded in a Hilbert space ( H , · H ) . The duality pairing between V * and V is denoted by · , · , and · , · H stands for the inner product in H. Let K 0 be nonempty, closed, star-shaped with respect to a closed ball B ( u 0 , ϱ ) in H, where u 0 V and ϱ > 0 . Let K and T K ( u ) for u K denote the realization of K 0 and T K 0 ( u ) in V, i.e.,
K : = K 0 V , T K ( u ) : = T K 0 ( u ) V ,
where T K 0 ( u ) stands for the Clarke tangent cone of K 0 at u.
Problem 1.
Find an element u K such that
A u f , v u + j 0 ( u ; v u ) + φ ( v ) φ ( u ) 0 for all v u + T K ( u ) .
The hypotheses on the data of Problem 1 are as follows.
H(A):
A : V V * is a mapping such that
(a)
A is pseudomonotone,
(b)
A is strongly monotone with constant m A > 0 , i.e.,
A v 1 A v 2 , v 1 v 2 m A v 1 v 2 V 2 for all v 1 , v 2 V .
H(j):
j : V R is a mapping such that
(a)
j is locally Lipschitz,
(b)
j ( v ) V * c 0 + c 1 v V for   all v V with c 0 , c 1 0 ,
(c)
there   exist α j 0 such   that
j 0 ( v 1 ; v 2 v 1 ) + j 0 ( v 2 ; v 1 v 2 ) α j v 1 v 2 V 2 for   all v 1 , v 2 V .
H(φ):
φ : V R is a convex and lower semicontinuous function.
H(K):
K is a nonempty, closed and star-shaped in V.
H(f):
f V * .
In hypothesis H ( j ) , the notation j 0 and j stand for the Clarke generalized directional derivative and the Clarke generalized subgradient, respectively, of the function j. We write j ( v ) V * = sup { v * V * v * j ( v ) } . Condition H ( j ) –(c) is known in the literature as a relaxed monotonicity condition, it holds with α j = 0 when j is a convex function, see [1,7,17,18] and the references therein. Examples of functions j that satisfy H ( j ) –(c) with single-valued and multivalued generalized subgradient can be found in Section 7.4 of [17] and Examples 16 and 17 in [1].
A sufficient condition for hypothesis H ( K ) reads as follows, see Theorem 7.4 of [8] and [12,13]. Let K i , i = 1 , , k , be nonempty, closed, convex subsets of V such that there is u 0 i = 1 k int K i . Then the set K : = i = 1 k is star-shaped with respect to a certain ball with center at u 0 .
We state below our first existence result.
Theorem 1.
Let the hypotheses H ( A ) , H ( j ) , H ( φ ) , H ( K ) , and H ( f ) hold, and assume the following smallness condition
m A > α j .
Then, Problem 1 admits a solution u K .
The proof of this theorem will be given in the next section. The motivation to study Problem 1 is given in the remark below.
Remark 1.
If K is a nonempty, closed and convex set, then any solution to Problem 1 is a solution to the classical variational-hemivariational inequality:
f i n d   a n   e l e m e n t u K s u c h   t h a t A u f , v u + j 0 ( u ; v u ) + φ ( v ) φ ( u ) 0 for all v K .
In fact, let u K solve Problem 1 and v K . The set K is regular, being convex, and by (3), we have
K u + K K ( u ) = u + T K ( u ) .
Therefore, we deduce that v u + T K ( u ) . Then, u K is a solution to problem (5).

4. Proof of Theorem 1

In this section, we assume the hypotheses of Theorem 1. The proof is based on the penalty method, where the penalty parameter is taken to be small and does not necessary tend to zero. Let d d K 0 : H R be the distance function of the set K 0 defined by d K 0 ( v ) = inf { v w H w K 0 } . Let λ > 0 represent a penalty parameter. Consider the penalized problem corresponding to Problem 1:
Problem 2.
Find an element u λ V such that
A u λ f , v u λ + j 0 ( u λ ; v u λ ) + φ ( v ) φ ( u λ ) + 1 λ d 0 ( u λ ; v u λ ) 0 for all v V .
We formulate a result on the generalized pseudomonotonicity property of a multivalued mapping.
Lemma 3.
Under the assumptions of Theorem 1, the multivalued mapping A + j + 1 λ d : V 2 V * is generalized pseudomonotone.
Proof. 
Let { w n } V , { w n * } V * with w n * A w n + j ( w n ) + 1 λ d ( w n ) , w n w weakly in V, w n * w * weakly in V * and
lim sup w n * , w n w 0 .
We need to show that w * A w + j ( w ) + 1 λ d ( w ) and w n * , w n w * , w .
We have w n * = A w n + ξ n + 1 λ η n with ξ n j ( w n ) and η n d ( w n ) . From hypothesis H ( j ) –(b) and the estimate d ( v ) H 1 for all v V , we infer that
ξ n ξ weakly   in V * and η n η weakly   in H ,
where ξ V * and η H . Since the embedding V H is compact, we get
w n w in H .
Moreover, by the closedness of the graph of d : H 2 H in H × ( w H ) -topology, see Proposition 3.23(v) in [17], we deduce that
η d ( w ) .
Subsequently, using the relation
A w n + j ( w n ) ( A w + j ( w ) ) , w n w + 1 λ d ( w n ) d ( w ) , w n w H = A w n + j ( w n ) + 1 λ d ( w n ) ( A w + j ( w ) + 1 λ d ( w ) ) , w n w ,
we have
( m A α j ) w n w V 2 + 1 λ η n η , w n w H w n * , w n w A w + j ( w ) + 1 λ d ( w ) , w n w .
Take the upper limit of both sides,
( m A α j ) lim sup n w n w V 2 + 1 λ lim n η n η , w n w H lim sup n w n * , w n w lim n A w + j ( w ) + 1 λ d ( w ) , w n w .
Exploiting (6)–(8), the boundeness of the mapping A + j + 1 λ d , and the smallness condition (4), we deduce
w n w in V .
Since every single-valued pseudomonotone operator is demicontinuous, see Theorem 3.69(ii) in [17], we infer A w n A w weakly in V * . From the convergences (7) and (9), we have ξ j ( w ) . Hence, taking the limit in the equality w n * = A w n + ξ n + 1 λ η n , we get
w * = A w + ξ + 1 λ η A w + j ( w ) + 1 λ d ( w ) .
Further, by (9), it is obvious that w n * , w n w * , w , which completes the proof. □
We continue the proof with three main steps.
Step 1. We show the existence of solution to Problem 2, for every λ > 0 fixed. For simplicity, we skip λ in this part of the proof. We define the multivalued mapping F : V 2 V by
F ( v ) = { u V M v ( u ) 0 } for v V
and
M v ( u ) : = A u f , v u + j 0 ( u ; v u ) + φ ( v ) φ ( u ) + 1 λ d 0 ( u ; v u )
for v V . We note that
u V solves Problem   2 u V satisfies M v ( u ) 0 for   all v V u F ( u ) for   all v V u v V F ( v ) .
Now, we prove that v V F ( v ) . We will apply the KKM lemma, see Lemma 1, with the space V endowed with the weak topology. We shall verify that the mapping F defined by (10) enjoys the properties:
(a)
for every v V , the set F ( v ) is closed in V,
(b)
for any finite set { v 1 , , v r } V , we have co { v 1 , , v r } i = 1 r F ( v i ) ,
(c)
there is v 0 V such that F ( v 0 ) is compact in V.
We show that the set F ( v ) in bounded in V for all v V . Let v V and u F ( v ) . Thus
A u f , v u + j 0 ( u ; v u ) + φ ( v ) φ ( u ) + 1 λ d 0 ( u ; v u ) 0 .
We show that u stays in a bounded subset of V. First, by hypothesis H ( j ) –(b), (c), we have
j 0 ( u ; v u ) = j 0 ( u ; v u ) + j 0 ( v ; u v ) j 0 ( v , u v ) α j u v V 2 + | j ( v ) , u v | α j u v V 2 + ( c 0 + c 1 v V ) u v V .
Next, using H ( φ ) and Proposition 5.2.25 in [27], it follows that φ admits an affine minorant, that is, we can find l V * and b R such that φ ( v ) l , v + b for all v V . Hence
φ ( u ) l , u b .
We employ the strong monotonicity of the mapping A, (11)–(13) to deduce
m A u v V 2 A u A v , u v = A u , u v A v , u v j 0 ( u ; v u ) + φ ( v ) φ ( u ) + f , u v + 1 λ d 0 ( u ; v u ) A v , u v α j u v V 2 + ( c 0 + c 1 v V ) u v V + | φ ( v ) | + l V * u V + | b | + f A v V * u v V + 1 λ d 0 ( u ; v u ) .
We exploit the global Lipschitz property of the function d, see Lemma 2.1 in [14], to infer that d ( v ) H 1 for all v V . Hence
d 0 ( u ; v u ) u v H c e u v V for   all u , v V ,
where c e > 0 denotes the embedding constant V H . Using (14) in the estimate, we obtain
( m A α j ) u v V 2 ( c 0 + c 1 v V ) u v V + | φ ( v ) | + l V * u V + | b | + f A v V * u v V + c e λ u v V .
By the smallness condition (4), since λ > 0 is fixed, we know that u v V is bounded, and thus u V is also bounded. This completes the proof that F ( v ) is a bounded subset of V.
First, we establish the property (a). Let λ > 0 and v V be fixed. It is enough to prove that the set F ( v ) is sequentially weakly closed in V. Let { u k } F ( v ) and u k u weakly in V, as k . We prove that u F ( v ) . We have M v ( u k ) 0 for all v V , that is,
0 A u k f , v u k + j 0 ( u k ; v u k ) + φ ( v ) φ ( u k ) + 1 λ d 0 ( u k ; v u k ) .
From Proposition 3.23(ii) in [17], there are ξ k V , η k H such that ξ k j ( u k ) , η k d ( u k ) with
j 0 ( u k ; v u k ) = ξ k , v u k a n d d 0 ( u k ; v u k ) = η k , v u k H = η k , v u k .
Subsequently, by (16), it follows
A u k + ξ k + 1 λ η k , u k v f , u k v + φ ( v ) φ ( u k ) for   all v V
and
lim sup k A u k + ξ k + 1 λ η k , u k v f , u v + φ ( v ) φ ( u ) for   all v V .
Here, we have used the fact that φ is weakly sequentially l.s.c. (being convex and sequentially l.s.c. by H ( φ ) ). Take v = u in (17) to get
lim sup k A u k + ξ k + 1 λ η k , u k u 0 .
Let u k * = A u k + ξ k + 1 λ η k A u k + j ( u k ) + 1 λ d ( u k ) . By the boundedness of the mapping A + j + 1 λ d , we may assume that
u k * u * weakly in V *
for some u * V * . Summing up, we obtain
{ u k } V , { u k * } V * , u k * = A u k + ξ k + 1 λ η k A u k + j ( u k ) + 1 λ d ( u k ) , u k u weakly   in V , u k * u * weakly   in V * and lim sup k u k * , u k u 0 .
Hence, by the generalized pseudomonotonicity of the mapping A + j + 1 λ d , see Lemma 3, we have
u * A u + j ( u ) + 1 λ d ( u ) ,
u k * , u k u * , u .
By (19), we know that u * = A u + ξ 0 + 1 λ η 0 with some ξ 0 j ( u ) , η 0 d ( u ) and
ξ 0 , v u j 0 ( u ; v u ) , η 0 , v u d 0 ( u ; v u ) for   all v V .
From (20), it follows
u k * , u k v u * , u v .
Using the latter, we can pass to the limit in inequality (17) to get
A u + ξ 0 + 1 λ η 0 , u v f , u v + φ ( v ) φ ( u ) for   all v V .
Hence and by (21) implies that u V satisfies
0 A u f , v u + j 0 ( u ; v u ) + φ ( v ) φ ( u ) + 1 λ d 0 ( u ; v u ) for   all v V .
This means that 0 M v ( u ) for all v V , that is, u F ( v ) for all v V . Thus, the set F ( v ) is sequentially weakly closed in V.
Second, we show the property (b), that is, F : V 2 V is a KKM map. Let r N , { z 1 , , z r } V be an arbitrary finite set. Let z = i = 1 r λ i z i with λ i ( 0 , 1 ) and i = 1 r λ i = 1 . We suppose by contradiction that
z i = 1 r F ( z i ) .
Then, for all i = 1 , , r , we have z F ( z i ) . So, for all i = 1 , , r , it holds M z i ( z ) < 0 . By this inequality and the convexity of the function z M z ( u ) for all u V , we have
0 = M z ( z ) = M i = 1 r λ i z i ( z ) i = 1 r λ i M z i ( z ) < 0
which is a contradiction. This proves (b).
We show property (c): the set F ( v ) is weakly compact in V for all v V . Let v V and { w n } V , w n F ( v ) . Since the sequence { w n V } is bounded in V by a constant independent of n N , by the reflexivity of V, it is clear that there is a subsequence { w n k } { w n } such that w n k w 0 weakly in V, with w 0 V . Because F ( v ) is sequentially weakly closed in V, see property (a), we have w 0 F ( v ) . Therefore, F ( v ) is weakly compact in V for all v V .
Having verified properties (a)–(c), for any λ > 0 fixed, by the KKM lemma, we get that v V F ( v ) . This means that, for any λ > 0 fixed, there exists u λ V solution to problem (2). This finishes Step 1.
Step 2. We show that the solution u λ V to Problem 2 obtained in Step 1 satisfies u λ K for λ > 0 sufficiently small. We claim that there is a constant M > 0 such that u λ V M for all λ > 0 . In fact, to obtain the uniform estimate of { u λ } λ > 0 in V, we choose v : = u 0 in Problem 2, where u 0 V is the center of the ball B ( u 0 , ϱ ) to get
A u λ f , u 0 u λ + j 0 ( u λ ; u 0 u λ ) + φ ( u 0 ) φ ( u λ ) + 1 λ d 0 ( u λ ; u 0 u λ ) 0 .
Analogously as in estimate (15), we have
( m A α j ) u λ u 0 V 2 ( c 0 + c 1 u 0 V ) u λ u 0 V + | φ ( u 0 ) | + l V * u λ V + | b | + f A u 0 V * u λ u 0 V + 1 λ d 0 ( u λ ; u 0 u λ ) .
Hence
( m A α j ) u λ V 2 d 1 u λ V + d 2 + 2 λ d 0 ( u λ ; u 0 u λ ) ,
where
d 1 = 2 ( c 0 + c 1 u 0 V ) + l V * + f A u 0 V * , d 2 = 2 ( ( c 0 + c 1 u 0 V ) u 0 V + | b | + | φ ( u 0 ) | + f A u 0 V * u 0 V + ( m A α j ) u 0 V ) .
By Lemma 2, it follows that d 0 ( u λ ; u 0 u λ ) 0 , so we can skip the last term in estimate (23). We deduce that u λ V M with M > 0 independent of λ > 0 .
Let λ 0 : = 2 ρ d 1 M + d 2 . We claim that for all λ ( 0 , λ 0 ) , it holds u λ K . We continue by contradiction and assume that
λ ( 0 , λ 0 ) such   that u λ K .
Since u λ K , by Lemma 2, we have
d 0 ( u λ ; u 0 u λ ) d ( u λ ) ρ ρ .
Again, by (23), we infer that
0 d 1 M + d 2 2 λ ρ ,
which implies
λ 2 ρ d 1 M + d 2 ,
a contradiction with the choice of λ . Hence
u λ K for all λ ( 0 , λ 0 ) ,
which proves the claim.
Step 3. Fix λ ( 0 , λ 0 ) with λ 0 defined in Step 2. We will show that u : = u λ solves Problem 1. From Step 1 we know that u λ V is a solution to Problem 2, and by Step 2, it is obvious that u λ K . Thus u = u λ K satisfies
A u f , v u + j 0 ( u ; v u ) + φ ( v ) φ ( u ) + 1 λ d 0 ( u ; v u ) 0 for all v V .
We choose v u + T K ( u ) as the test function in the latter. Since v u T K ( u ) , we have d 0 ( u ; v u ) = 0 . Hence
A u f , v u + j 0 ( u ; v u ) + φ ( v ) φ ( u ) 0 for all v u + T K ( u ) .
Finally, u K is a solution to Problem 1. This finishes the proof. □

5. Second Existence Result

In this section we deliver an existence result for Problem 1 under hypotheses different than the ones of Section 3. We do not assume the strong monotonicity of the mapping A, the relaxed monotonicity of the generalized subgradient j , and consider the nonconvex star-shaped admissible set of constraints.
We admit the following assumptions.
H(A)1:
A : V V * is a mapping such that
(a)
for all v V , V u A u , v u R is weakly upper semicontinuous,
(b)
for any v V , there exists m v > 0 such that A u , u v α A u 2 for
all u m v , where α A > 0 .
H(j)1:
j : V R is a mapping such that
(a)
j is locally Lipschitz,
(b)
j ( v ) V * c 0 + c 1 v V for   all v V with c 0 , c 1 0 ,
(c)
lim sup j 0 ( u n ; v u n ) j 0 ( u ; v u ) for all v V and u n u weakly in V.
The following example provides a sufficient condition for H ( j ) 1 –(c). Let Y be a reflexive Banach space, ψ : Y R be a locally Lipschitz function such that ψ or ψ is regular, and M : V Y be given by M v = L v + v 0 , where L : V Y represents a linear, compact operator and v 0 Y is fixed. In this situation, the function j : V R defined by j ( v ) = ψ ( M v ) for v V satisfies H ( j ) 1 –(c). In applications, M is a compact trace operator or a compact embedding operator, see, for instance, Section 6.
We obtain the second existence result by employing the approach used in Theorem 1.
Theorem 2.
Under hypotheses H ( A ) 1 , H ( j ) 1 , H ( φ ) , H ( K ) , H ( f ) and the smallness condition
α A > c 1 ,
Problem 1 has a solution u K .
Proof. 
We use the notation of Section 4 and treat Problem 1 by the penalized inequality stated in Problem 2. We follow the three steps in the proof of Theorem 1. We only indicate below the new ingredients of the proof.
(i) We prove that the set F ( v ) in bounded in V for all v V , where the multivalued mapping F : V 2 V is given by (10). Let v V and u F ( v ) . Hence
A u , u v f , u v + j 0 ( u ; v u ) + φ ( v ) φ ( u ) + 1 λ d 0 ( u ; v u ) .
From this inequality, we use H ( J ) 1 –(b), (13) and (14) to obtain
A u , u v f V * u V + f V * v V + ( c 0 + c 1 u V ) u V + ( c 0 + c 1 u V ) v V + | φ ( v ) | + l V * u V + | b | + c e λ u V + c e λ v V .
By the condition H ( A ) 1 –(a) and the smallness condition (24), we deduce that u V is bounded by a constant which depends on v but is independent of u. This completes the proof that F ( v ) is a bounded set in V.
(ii) We prove that for λ > 0 and v V fixed, the set F ( v ) is sequentially weakly closed in V. Let { u k } F ( v ) and u k u weakly in V, as k . This means that M v ( u k ) 0 for all v V , and
0 A u k f , v u k + j 0 ( u k ; v u k ) + φ ( v ) φ ( u k ) + 1 λ d 0 ( u k ; v u k ) .
We show that M v ( · ) : V R is sequentially weakly upper semicontinuous, that is,
lim sup M v ( u k ) M v ( u ) .
We take upper limit in (27) to get
0 lim sup A u k , v u k lim f , v u k + lim sup j 0 ( u k ; v u k ) + φ ( v ) + lim sup ( φ ( u k ) ) + 1 λ lim sup d 0 ( u k ; v u k ) .
Again, by the compactness of the embedding V H , we apply Proposition 3.23(ii) in [17] to deduce
lim sup d 0 ( u k ; v u k ) d 0 ( u ; v u ) .
We use H ( A ) 1 –(a), H ( j ) 1 –(c), H ( φ ) and (29) in (28) to obtain
0 A u f , v u + j 0 ( u ; v u ) + φ ( v ) φ ( u ) + 1 λ d 0 ( u ; v u ) .
Hence 0 M v ( u ) and u F ( v ) . This completes this proof.
(iii) We show that the solution u λ V to Problem 2 satisfies u λ V M 1 for all λ > 0 with M 1 > 0 independent of λ . We take v : = u 0 in the penalized Problem 2, where, recall, u 0 V is the center of the ball B ( u 0 , ϱ ) . We have
0 A u λ f , u 0 u λ + j 0 ( u λ ; u 0 u λ ) + φ ( u 0 ) φ ( u λ ) + 1 λ d 0 ( u λ ; u 0 u λ ) .
Simarily as in (26), we infer that
A u λ , u λ v f V * u λ V + f V * v V + ( c 0 + c 1 u λ V ) u λ V + ( c 0 + c 1 u λ V ) v V + | φ ( v ) | + l V * u λ V + | b | + c e λ u λ V + c e λ v V .
Again, by the condition H ( A ) 1 –(a) and the smallness condition (24), we deduce that u λ V is bounded by a constant M 1 > 0 which is independent of λ .
The remaining parts of the proof follow the lines in the proof of Theorem 1. □
From Remark 1 and Theorem 2, the deduce the following result on the existence of solution to the variational-hemivariational inequality with the convex constraint set.
Corollary 1.
Assume hypotheses H ( A ) 1 , H ( j ) 1 , H ( φ ) , H ( f ) and the smallness condition (24). If K is a nonempty, closed and convex set, then Problem 5 has a solution u K .
We conclude the section with comments on the assumptions H ( A ) –(a) and H ( A ) 1 –(a).
(1) If A : V V * is monotone, bounded and ( w V ) × ( w V * ) continuous, then H ( A ) –(a) holds. This follows by an observation that every ( w V ) × ( w V * ) continuous operator is demicontinuous. Then, the notions of demicontinuity and hemicontinuity coincide for monotone operators, see Exercise VI.9 in [10]. Finally, by Theorem 3.69(i) in [17], a bounded, monotone and hemicontinuous mapping is pseudomonotone, that is, H ( A ) –(a) holds. Note that linear and bounded operators, and most quasilinear differential operators are weakly-weakly continuous, see e.g., [11].
(2) If A : V V * is monotone and ( w V ) × ( w V * ) continuous, then H ( A ) 1 –(a) holds. To prove this, first, we observe that if A is monotone, then for all u n u weakly in V, we have lim sup A u n , u u n 0 . Indeed, we proceed by contradiction. Suppose that there are { w n } V , w V such that w n w weakly in V and lim sup A w n , w w n > 0 . The latter is equivalent to lim inf A w n , w n w < 0 . On the other hand by the monotonicity of A, we have
A w , w n w A w n , w n w .
We take lower limit in this inequality, and deduce
0 = lim inf A w , w n w lim inf A w n , w n w < 0 ,
which is a contradiction. Second, to prove H ( A ) 1 –(a), let v V , u n u weakly in V. Then, A u n A u weakly in V * and
lim sup A u n , v u n lim sup ( A u n , v u + A u n , u u n ) lim sup A u n , v u + lim sup A u n , u u n A u , v u ,
where we have used the inequality lim sup A u n , u u n 0 . Hence, H ( A ) 1 –(a) is verified.
(3) It is clear that if the mapping A is strongly monotone with constant m A > 0 , see H ( A ) –(b), then A is coercive in the sense that
A u , u m A u V 2 A 0 V * u V for   all u V ,
and
A u , u v m A u v V 2 A v V * u v V for   all u , v V ,
compare with H ( A ) 1 –(b).

6. Semipermeability Model

In this section we provide an illustrative model which weak formulation leads to a constrained variational-hemivariational inequality. Based on this model, we justify the nature of operator M which may appear in applications.
Consider the following semipermeability model for the stationary heat conduction problem. Let Ω R d be a bounded domain with the regular boundary Γ which consists of three mutually disjoint and relatively open subsets Γ 1 , Γ 2 and Γ 3 with Γ = Γ ¯ 1 Γ ¯ 2 Γ ¯ 3 and m ( Γ 1 ) > 0 . Consider the following boundary value problem.
Problem 3.
Find a temperature u : Ω R such that u U and
A u + j 1 ( u ) + g 1 ( u ) f 1 in Ω u = 0 on Γ 1 u ν A + j 2 ( u ) f 2 on Γ 2 u ν A + g 2 ( u ) f 3 on Γ 3 ,
where A : V V * is a given linear mapping, ν denotes the outward normal on the boundary, and u ν A is the conormal derivative with respect to A and represents the heat flux through a part of the boundary. Further, U V is a nonempty, closed set of constraints which can be convex or nonconvex, and V = { v H 1 ( Ω ) v = 0 o n Γ 1 } . We denote by i : V L 2 ( Ω ) the embedding operator and by γ : V L 2 ( Γ ) the trace operator. It is well known that both operators are linear and compact.
Problem 3 is motivated by several kinds of semipermeability relations which arise in several situations in hydraulics, flow problems through porous media and electrostatics, the solution can represent temperature, pressure and the electric potential Chapter I in [31] (where the monotone semipermeability relations were considered with convex potentials), and Chapter 5.5.3 of [8,16] (where nonmonotone subdifferential conditions were treated). Mappings g 1 and j 1 describe the interior semipermeability phenomena while g 2 and j 2 provide the boundary semipermeability relations in the subdifferential form. Additional constraints for the temperature (or the pressure of the fluid in a fluid flow model) are represented by the condition u U . The set U can be employed to introduce a bilateral obstacle which means that we look for the temperature within prescribed bounds in the domain Ω . The function f 1 corresponds to the density of heat sources in the domain. The multivalued subdifferential boundary conditions on Γ 2 (and Γ 3 ) describe the nonmonotone (and monotone, respectively) behavior of a semipermeable membrane (a wall) of finite thickness, and appear in a temperature control problem, see [32].
We need the following hypotheses on the data.
H(A)2:
A : V V * is a mapping such that A = i , j = 1 d D i a i j ( x ) D j , and
(i)
a i j L ( Ω ) for i, j = 1 , , d .
(ii)
i , j = 1 d a i j ( x ) ξ i ξ j α 0 ξ 2 for all ξ R d , a.e. x Ω with α 0 > 0 .
H(j1):
j 1 : R R is such that
(i)
j 1 is locally Lipschitz.
(ii)
| j 1 ( r ) | c 0 j + c 1 j | r | for all r R with c 0 j , c 1 j 0 .
(iii)
( j 1 ( r 1 ) j 1 ( r 2 ) ) ( r 1 r 2 ) β 1 j | r 1 r 2 | 2 all r 1 , r 2 R with β 1 j 0 .
H(j2):
j 2 : R R is such that
(i)
j 2 is locally Lipschitz.
(ii)
| j 2 ( r ) | c 2 j + c 3 j | r | for all r R with c 2 j , c 3 j 0 .
(iii)
( j 2 ( r 1 ) j 2 ( r 2 ) ) ( r 1 r 2 ) β 2 j | r 1 r 2 | 2 all r 1 , r 2 R with β 2 j 0 .
H(g1):
g 1 : R R is such that
(i)
g 1 is convex and l.s.c.
(ii)
| g 1 ( r ) | c 0 g + c 1 g | r | for all r R with c 0 g , c 1 g 0 .
H(g2):
g 2 : × R R is such that
(i)
g 2 is convex and l.s.c.
(ii)
| g 2 ( r ) | c 2 g + c 3 g | r | for all r R with c 2 g , c 3 g 0 .
H(f):
f 1 L 2 ( Ω ) , f 2 L 2 ( Γ 2 ) , f b L 2 ( Γ 3 ) .
(H0):
α 0 > β 1 j i 2 + β 2 j γ 2 .
Assume that U is a convex set. Under the hypotheses above, by a standard procedure, we obtain the following weak formulation of Problem 3.
Problem 4.
Find u U such that
A u f , v u + Ω j 1 0 ( i u ; i v i u ) d x + Γ 2 j 2 0 ( γ u ; γ v γ u ) d Γ + Ω ( g 1 ( i v ) g 1 ( i u ) ) d x + Γ 3 ( g 2 ( γ v ) g 2 ( γ u ) ) d Γ 0 for all v U ,
where f V * is given by
f , v = Ω f 1 i v d x + Γ 2 f 2 γ v d Γ + Γ 3 f 3 γ v d Γ for v V .
We introduce the functionals G 1 , G 2 , J 1 , J 2 : V R defined by
G 1 ( v ) = Ω g 1 ( i v ) d x , G 2 ( v ) = Γ 3 g 2 ( γ v ) d Γ , J 1 ( v ) = Ω j 1 ( i v ) d x , J 2 ( v ) = Γ 2 j 2 ( γ v ) d Γ
for v V . Let φ = G 1 + G 2 and j = J 1 + J 2 . We observe that from Propositions 3.37(i) and 3.46(iv) in [17], we get
j 0 ( u ; v u ) J 1 0 ( u ; v u ) + J 2 0 ( u ; v u ) Ω j 1 0 ( i u ; i v i v ) d x + Γ 2 j 2 0 ( γ u ; γ v γ u ) d Γ for   all u , v V .
Using the last inequality, definitions of the convex subdifferential and the Clarke subgradient, we easily deduce that if u U solves the problem: find an element u U such that
A u f , v u + φ ( v ) φ ( u ) + j 0 ( u ; v u ) 0 for all v U ,
then u is a solution to Problem 4.
If the set of constraints U is a nonempty, closed and nonconvex subset of V, we can derive a weak formulation as in Problem 4. Then, Theorems 1 and 2 can be applied in this situation.

7. Conclusions

In this paper, we have given some sufficient conditions for the existence of solution to a class of variational-hemivariational inequality problems involving the Clarke tangent cone of the constraint set in a reflexive Banach space. The main feature of this class is the nonconvexity of the constraint set, and the presence of two potential which are convex and locally Lipschitz, respectively. In the proofs, we have employed the well-known KKM lemma combined with the penalty method without making the small parameter tend to zero.
It is an intersting open problem to establish existence results without the smallness hypotheses (4) and (24). The results will find important applications to model the semipermeable media, contact problems in solid and fluid mechanics, etc. Moreover, it would be desirable to extend the results with nonconvex constraints sets to second order evolution problems motivated by dynamic contact models in viscoelasticity, thermoviscoelasticty, see [17,19], and nonstationary fluid models, see [20,21].

Author Contributions

Conceptualization, S.M. and L.F.; methodology, S.M. and L.F.; formal analysis, S.M. and L.F.; investigation, S.M. and L.F.; resources, S.M. and L.F.; writing—original draft preparation, S.M. and L.F.; writing—review and editing, S.M. and L.F.; supervision, S.M. All authors have read and agreed to the published version of the manuscript.

Funding

The project is supported by the European Union’s Horizon 2020 Research and Innovation Programme under the Marie Skłodowska-Curie grant agreement No. 823731 CONMECH, NSF of Guangxi Grant No. 2018GXNSFAA281353, the Beibu Gulf University Project No. 2018KYQD03, and the projects financed by the Ministry of Science and Higher Education of Republic of Poland under Grants Nos. 4004/GGPJII/H2020/2018/0 and 440328/PnH2/2019.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Migórski, S.; Ochal, A.; Sofonea, M. A class of variational-hemivariational inequalities in reflexive Banach spaces. J. Elast. 2017, 127, 151–178. [Google Scholar] [CrossRef] [Green Version]
  2. Baiocchi, C.; Capelo, A. Variational and Quasivariational Inequalities: Applications to Free-Boundary Problems; John Wiley: Chichester, UK, 1984. [Google Scholar]
  3. Han, W.; Sofonea, M. Quasistatic Contact Problems in Viscoelasticity and Viscoplasticity; Studies in Advanced Mathematics 30; American Mathematical Society: Providence, RI, USA; International Press: Somerville, MA, USA, 2002. [Google Scholar]
  4. Lions, J.L.; Stampacchia, G. Variational inequalities. Comm. Pure Appl. Math. 1967, 20, 493–519. [Google Scholar] [CrossRef] [Green Version]
  5. Sofonea, M.; Matei, A. Mathematical Models in Contact Mechanics; London Mathematical Society Lecture Note Series 398; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  6. Carl, S.; Le, V.K.; Motreanu, D. Nonsmooth Variational Problems and their Inequalities; Springer: New York, NY, USA, 2007. [Google Scholar]
  7. Han, W.; Migórski, S.; Sofonea, M. (Eds.) Advances in Variational and Hemivariational Inequalities with Applications. Theory, Numerical Analysis, and Applications. In Advances in Mechanics and Mathematics; Springer: Berlin/Heidelberg, Germany, 2015; Volume 33. [Google Scholar]
  8. Naniewicz, Z.; Panagiotopoulos, P.D. Mathematical Theory of Hemivariational Inequalities and Applications; Dekker: New York, NY, USA, 1995. [Google Scholar]
  9. Panagiotopoulos, P.D. Inequality Problems in Mechanics and Applications; Birkhäuser: Boston, MA, USA, 1985. [Google Scholar]
  10. Denkowski, Z.; Migórski, S.; Papageorgiou, N.S. An Introduction to Nonlinear Analysis: Applications; Kluwer Academic Publishers: Boston, MA, USA, 2003. [Google Scholar]
  11. Francǔ, J. Weakly continuous operators, Applications to differential equations. Appl. Math. 1994, 39, 45–56. [Google Scholar] [CrossRef]
  12. Naniewicz, Z. Hemivariational inequality approach to constrained problems for star-shaped admissible sets. J. Optim. Theory Appl. 1994, 83, 97–112. [Google Scholar] [CrossRef]
  13. Goeleven, D. On the hemivariational inequality approach to nonconvex constrained problems in the theory of von Kárman plates. Z. Angew. Math. Mech. (ZAMM) 1995, 75, 861–866. [Google Scholar] [CrossRef]
  14. Gasiński, L.; Liu, Z.; Migórski, S.; Ochal, A.; Peng, Z. Hemivariational inequality approach to evolutionary constrained problems on star-shaped sets. J. Optimiz. Theory Appl. 2015, 164, 514–533. [Google Scholar] [CrossRef]
  15. Peng, Z.; Gasinski, L.; Migórski, S.; Ochal, A. A class of evolution variational inequalities with nonconvex constraints. Optimization 2019, 68, 1881–1895. [Google Scholar] [CrossRef]
  16. Panagiotopoulos, P.D. Hemivariational Inequalities, Applications in Mechanics and Engineering; Springer: Berlin, Germany, 1993. [Google Scholar]
  17. Migórski, S.; Ochal, A.; Sofonea, M. Nonlinear Inclusions and Hemivariational Inequalities. Models and Analysis of Contact Problems; Advances in Mechanics and Mathematics 26; Springer: New York, NY, USA, 2013. [Google Scholar]
  18. Sofonea, M.; Migórski, S. Variational-hemivariational Inequalities with Applications; Monographs and Research Notes in Mathematics; Chapman & Hall/CRC: Boca Raton, FL, USA, 2017. [Google Scholar]
  19. Migórski, S.; Ochal, A. A unified approach to dynamic contact problems in viscoelasticity. J. Elast. 2006, 83, 247–275. [Google Scholar] [CrossRef]
  20. Dudek, S.; Migórski, S. Evolutionary Oseen model for generalized Newtonian fluid with multivalued nonmonotone friction law. J. Math. Fluid Mech. 2018, 20, 1317–1333. [Google Scholar]
  21. Migórski, S.; Dudek, S. A new class of variational-hemivariational inequalities for steady Oseen flow with unilateral and frictional type boundary conditions. Z. Angew. Math. Mech. 2020, 100, e201900112. [Google Scholar] [CrossRef]
  22. Migórski, S.; Pa̧czka, D. On steady flow of non-Newtonian fluids with frictional boundary conditions in reflexive Orlicz spaces. Nonlinear Anal. Ser. Real World Appl. 2018, 39, 337–361. [Google Scholar] [CrossRef]
  23. Migórski, S.; Pa̧czka, D. Frictional contact problems for steady flow of incompressible fluids in Orlicz spaces. In Current Trends in Mathematical Analysis and Its Interdisciplinary Applications; Birkhäuser: Basel, Switzerland, 2019; Chapter 1; pp. 1–53. [Google Scholar]
  24. Anh, P.N.; Thach, H.T.C.; Kim, J.K. Proximal-like subgradient methods for solving multi-valued variational inequalities. Nonlinear Funct. Anal. Appl. 2020, 25, 437–451. [Google Scholar]
  25. Papageorgiou, N.S.; Scapellato, A. Nonlinear Robin problems with general potential and crossing reaction. Rend. Lincei Mat. Appl. 2019, 30, 1–29. [Google Scholar] [CrossRef]
  26. Clarke, F.H. Optimization and Nonsmooth Analysis; Wiley Interscience: New York, NY, USA, 1983. [Google Scholar]
  27. Denkowski, Z.; Migórski, S.; Papageorgiou, N.S. An Introduction to Nonlinear Analysis: Theory; Kluwer Academic Publishers: Boston, MA, USA, 2003. [Google Scholar]
  28. Fan, K. A generalization of Tychonoff’s fixed point theorem. Math. Ann. 1961, 142, 305–310. [Google Scholar] [CrossRef]
  29. Park, S. Evolution of the 1984 KKM theorem of Ky Fan. Fixed Point Theory Appl. 2012, 2012, 146. [Google Scholar] [CrossRef] [Green Version]
  30. Clarke, F.H. Functional Analysis, Calculus of Variations and Optimal Control; Springer: London, UK, 2013. [Google Scholar]
  31. Duvaut, G.; Lions, J.-L. Inequalities in Mechanics and Physics; Springer: Berlin, Germany, 1976. [Google Scholar]
  32. Migórski, S.; Ochal, A. Boundary hemivariational inequality of parabolic type. Nonlinear Anal. Theory Methods Appl. 2004, 57, 579–596. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Migórski, S.; Fengzhen, L. Constrained Variational-Hemivariational Inequalities on Nonconvex Star-Shaped Sets. Mathematics 2020, 8, 1824. https://doi.org/10.3390/math8101824

AMA Style

Migórski S, Fengzhen L. Constrained Variational-Hemivariational Inequalities on Nonconvex Star-Shaped Sets. Mathematics. 2020; 8(10):1824. https://doi.org/10.3390/math8101824

Chicago/Turabian Style

Migórski, Stanisław, and Long Fengzhen. 2020. "Constrained Variational-Hemivariational Inequalities on Nonconvex Star-Shaped Sets" Mathematics 8, no. 10: 1824. https://doi.org/10.3390/math8101824

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop