Existence, comparison, and convergence results for a class of elliptic hemivariational inequalities

In this paper we study a class of elliptic boundary hemivariational inequalities which originates in the steady-state heat conduction problem with nonmonotone multivalued subdifferential boundary condition on a portion of the boundary described by the Clarke generalized gradient of a locally Lipschitz function. First, we prove a new existence result for the inequality employing the theory of pseudomonotone operators. Next, we give a result on comparison of solutions, and provide sufficient conditions that guarantee the asymptotic behavior of solution, when the heat transfer coefficient tends to infinity. Further, we show a result on the continuous dependence of solution on the internal energy and heat flux. Finally, some examples of convex and nonconvex potentials illustrate our hypotheses.


Introduction
We consider a bounded domain Ω in R d whose regular boundary Γ consists of the union of three disjoint portions Γ i , i = 1, 2, 3 with |Γ i | > 0, where |Γ i | denotes the (d − 1)-dimensional Hausdorff measure of the portion Γ i on Γ. The outward normal vector on the boundary is denoted by n. We formulate the following two steady-state heat conduction problems with mixed boundary conditions: −∆u = g in Ω, where u is the temperature in Ω, g is the internal energy in Ω, b is the temperature on Γ 3 for (1) and the temperature of the external neighborhood of Γ 3 for (2), q is the heat flux on Γ 2 and α > 0 is the heat transfer coefficient on Γ 3 , which satisfy the hypothesis: g ∈ L 2 (Ω), q ∈ L 2 (Γ 2 ) and b ∈ H 1 2 (Γ 3 ). Throughout the paper we use the following notation where γ : V → L 2 (Γ) denotes the trace operator on Γ. In what follows, we write u for the trace of a function u ∈ V on the boundary. In a standard way, we obtain the following variational formulations of (1) and (2), respectively: The standard norms on V and V 0 are denoted by It is well known by the Poincaré inequality, see [5,Proposition 2.94], that on V 0 the above two norms are equivalent. Note that the form a is bilinear, symmetric, continuous and coercive with constant m a > 0, i.e.
It is well known that the regularity of solution to the mixed elliptic problems (1) and (2) is problematic in the neighborhood of a part of the boundary, see for example the monograph [12]. A regularity results for elliptic problems with mixed boundary conditions can be found in [1,2,14]. Moreover, sufficient hypothesis on the data in order to have H 2 regularity for elliptic variational inequalities is given in [23]. We remark that, under additional hypotheses on the data g, q and b, problems (1) and (2) can be considered as steady-state two phase Stefan problems, see, for example, [10,26,28,30].
The problems (3) and (4) have been extensively studied in several papers such as [10,26,27,28,29]. Some properties of monotonicity and convergence, when the parameter α goes to infinity, obtained in the aforementioned works, are recalled in the following result.
Theorem 1. If the data satisfy b = const. > 0, g ∈ L 2 (Ω) and q ∈ L 2 (Γ 2 ) with the properties q ≥ 0 on Γ 2 and g ≤ 0 in Ω, then The main goal of this paper is to study a generalization of problem (2) for which we provide sufficient conditions that guarantee the comparison properties and asymptotic behavior, as α → ∞, stated in Theorem 1. Moreover, for a more general problem, we also show a result on the continuous dependence of solution on the data g and q.
The mixed nonlinear boundary value problem for the elliptic equation under consideration reads as follows.
Here α is a positive constant while the function j : Γ 3 ×R → R, called a superpotential (nonconvex potential), is such that j(x, ·) locally Lipschitz for a.e. x ∈ Γ 3 and not necessary differentiable. Since in general j(x, ·) is nonconvex, so the multivalued condition on Γ 3 in problem (6) is described by a nonmonotone relation expressed by the generalized gradient of Clarke. Such multivalued relation in problem (6) is met in certain types of steady-state heat conduction problems (the behavior of a semipermeable membrane of finite thickness, a temperature control problems, etc.). Further, problem (6) can be considered as a prototype of several boundary semipermeability models, see [15,19,20,32], which are motivated by problems arising in hydraulics, fluid flow problems through porous media, and electrostatics, where the solution represents the pressure and the electric potentials. Note that the analogous problems with maximal monotone multivalued boundary relations (that is the case when j(x, ·) is a convex function) were considered in [3,9], see also references therein.
Under the above notation, the weak formulation of the elliptic problem (6) becomes the following boundary hemivariational inequality: Here and in what follows we often omit the variable x and we simply write j(r) instead of j(x, r). Observe that if j(x, ·) is a convex function for a.e. x ∈ Γ 3 , then the problem (7) reduces to the variational inequality of second kind: Note that when j(r) = 1 2 (r − b) 2 , problem (8) reduces to a variational inequality corresponding to problem (2). Several other examples of convex potentials can be found in various diffusion problems. For instance, the following convex functions: with suitable constants β > 0 and c ∈ R, appear in models which describe a free boundary problem with Tresca condition, see [4], the Stefan-Boltzman heat radiation law, and the natural convection, respectively, see [3,13], and the references therein for further applications and extensions. On the other hand, the stationary heat conduction models with nonmonotone multivalued subdifferential interior and boundary semipermeability relations can not be described by convex potentials. They use locally Lipschitz potentials and their weak formulations lead to hemivariational inequalities, see [19,Chapter 5.5.3] and [20]. We mention that theory of hemivariational and variational inequalities has been proposed in the 1980s by Panagiotopoulos, see [19,21,22], as variational formulations of important classes of inequality problems in mechanics. In the last few years, new kinds of variational, hemivariational, and variational-hemivariational inequalities have been investigated, see recent monographs [5,18,25], and the theory has emerged today as a new and interesting branch of applied mathematics.
The rest of the paper is structured as follows. In Section 2 we provide a new existence result for problem (7). In Section 3 we establish two comparison properties for solutions to problem (7). The convergence result of solution of problem (7) to the solution of problem (3), when the parameter α goes to infinity, is provided in Section 4. In Section 5 we study the continuous dependence of solution to problem (7) on the internal energy g and the heat flux q. The proofs are based on arguments of compactness, lower semicontinuity, monotonicity, various estimates, the theory of elliptic hemivariational inequalities and nonsmooth analysis [6,7,8,9,11,18,22,24,25,31]. Finally, in Section 6 we deliver several examples of convex and nonconvex potentials which satisfy the hypotheses on function j required in this paper.

Preliminaries
In this section first recall standard notation and preliminary concepts, and then provide a new result on existence of solution to the elliptic hemivariational inequality (7).
Let (X, · X ) be a reflexive Banach space, X * be its dual, and ·, · denote the duality between X * and X. For a real valued function defined on X, we have the following definitions [6, Section 2.1] and [7,18].
Definition 2. A function ϕ : X → R is said to be locally Lipschitz, if for every x ∈ X there exist U x a neighborhood of x and a constant L x > 0 such that For such a function the generalized (Clarke) directional derivative of j at the point x ∈ X in the direction v ∈ X is defined by The generalized gradient (subdifferential) of ϕ at x is a subset of the dual space X * given by We shall use the following properties of the generalized directional derivative and the generalized gradient, see [18,Proposition 3.23].
Proposition 3. Assume that ϕ : X → R is a locally Lipschitz function. Then the following hold: (i) for every x ∈ X, the function X ∋ v → ϕ 0 (x; v) ∈ R is positively homogeneous, and subadditive, i.e., respectively.
(iv) for every x ∈ X, the gradient ∂ϕ(x) is a nonempty, convex, and weakly compact subset of X * .
(v) the graph of the generalized gradient ∂ϕ is closed in X × (weak -X * )-topology. Now, we pass to a result on existence of solution to the elliptic hemivariational inequality: We admit the following standing hypothesis.
Note that the existence results for elliptic hemivariational inequalities can be found in several contributions, see [5,16,17,18,19]. In comparison to other works, the new hypothesis is H(j)(d). Under this condition we will show both existence of solution to problem (9) and a convergence result when α → ∞. We underline that, if the hypothesis H(j)(d) is replaced by the relaxed monotonicity condition (see Remark 10 for details) for all r, s ∈ R, a.e. x ∈ Γ 3 with m j ≥ 0, and the following smallness condition m a > α m j γ 2 is assumed, then problem (9) is uniquely solvable, see [17,Lemma 20] for the proof. However, this smallness condition is not suitable in the study of problem (9) since for a sufficiently large value of α, it is not satisfied. In the following result we apply a surjectivity result in [18,Proposition 3.61] and partially follow arguments of [17,Lemma 20]. For completeness we provide the proof. Theorem 4. If H(j) holds, f ∈ V * 0 and α > 0, then the hemivariational inequality (9) has a solution.
Proof. Let ·, · stand for the duality pairing between It is obvious that the operator A is linear, bounded and coercive, i.e., From H(j)(a)-(c), by [18,Corollary 4.15], we infer that the functional J enjoys the following properties: (p1) J is well defined and Lipschitz continuous on bounded subsets of L 2 (Γ 3 ), hence also locally Lipschitz, We introduce the operator B : where γ * : L 2 (Γ) → V * 0 denotes the adjoint to the trace γ. We show that B is pseudomonotone and bounded from V 0 to 2 V * 0 , see [18, Definition 3.57]. By Proposition 3 (iv), it follows that the values of ∂J are nonempty, convex and weakly compact subsets of L 2 (Γ 3 ). Hence, the set Bv is nonempty, closed and convex in V * 0 for all v ∈ V 0 . The operator B is bounded which is a consequence of the following estimate where γ denotes the norm of the trace operator. In order to establish pseudomonotonicity of the operator B, we take into account [18, Proposition 3.58(ii)], and prove that B is generalized pseudomonotone. Let n ∈ Bv n , it follows v * n = α γ * η n with η n ∈ ∂J(γv n ). By the estimate (p3), it is clear that {η n } remains in a bounded subset of L 2 (Γ 3 ). Thus, at least for a subsequence, denoted in the same way, we may suppose that η n → η weakly in L 2 (Γ 3 ) with η ∈ L 2 (Γ 3 ). Using the compactness of the trace operator, we have γv n → γv in L 2 (Γ 3 ) Now, we employ the strong-weak closedness of the graph of ∂J, see Proposition 3 (v), to obtain η ∈ ∂J(γv). On the other hand, by v * n = α γ * η n , it This completes the proof that B is generalized pseudomonotone. Hence, the operator B is also pseudomonotone.
Subsequently, we note that A : V 0 → V * 0 is pseudomonotone, see [18,Theorem 3.69], since it is linear, bounded and nonnegative. Therefore, A is pseudomonotone and bounded as a multivalued operator from V 0 to 2 V * 0 , see [18,Section 3.4]. Since the sum of multivalued pseudomonotone operators remains pseudomonotone, see [18, Proposition 3.59 (ii)], we infer that A + B is bounded and pseudomonotone.
Next, we prove that the operator A + B is coercive. In view of the coercivity of A, it is enough to show that Hence, by the definition of the generalized gradient and the property (p2), we obtain Using the latter and the equality which proves (10). In consequence, we have We conclude that the multivalued operator A + B is bounded, pseudomonotone, and coercive, hence surjective, see [18,Proposition 3.61]. We infer that there exists u ∈ V 0 such that (A + B)u ∋ f . In the final step of the proof, we observe that any solution u ∈ V 0 to the inclusion (A + B)u ∋ f is a solution to problem (9). Indeed, we have Combining the latter with the definition of the generalized gradient and the property (p2), we obtain for all v ∈ V 0 . This means that u ∈ V 0 solves problem (9). This completes the proof.

Comparison results
In this section we study the following two problems under the standing hypothesis H(j) on the superpotential. For every α > 0, we consider the hemivariational inequality of the form and the weak form of the elliptic equation Recall that It follows from Theorem 4 that for each α > 0, problem (11) has a solution u α ∈ V 0 while [5, Corollary 2.102] entails that problem (12) has a unique solution u ∞ ∈ K. Moreover, it is easy to observe that problem (12) can be equivalently formulated as follows In what follows we need the hypothesis on the data.
where u α ∈ V 0 is a solution to problem (11) and u ∞ ∈ K is the unique solution to problem (12).
By the linearity of the form a, we easily obtain a(u α , −w + ) = −a(w + , w + ), Subsequently, H(j)(d) entails Using the coercivity condition (5) of the form a and (14), we deduce m a w + 2 V ≤ 0. Hence w + = 0 in Ω, and finally u α ≤ b in Ω.
(b) We denote w = u α − u ∞ . It is enough to show that w + = 0 in Ω. We observe that w Γ 1 = 0. This allows to choose v = −w + ∈ V 0 in problem (11) to obtain Next, part (a) of the proof shows that and w + Γ 3 = 0, and consequently w + ∈ K 0 . Since u ∞ ∈ K solves (12), taking v = w + ∈ K 0 in equality (12), and using the result in (15), it follows that Since u ∞ = b on Γ 3 , by (14), we get Again, by the coercivity of the form a, we have m a w + 2 V ≤ 0. Therefore, w + = 0 in Ω, and finally u α ≤ u ∞ in Ω. This completes the proof.

Note that properties (a) and (b) of Theorem 5 obtained for the hemivariational inequality (11) have been provided for linear elliptic problem (3) in properties (ii) and (iii) of Theorem 1.
In what follows, we comment on the monotonicity property analogous to condition (iv) stated for problem (3) in Theorem 1. for all c ≥ 1, all r, s ∈ R, a.e. x ∈ Γ 3 . Let u α i ∈ V 0 denote the unique solution to the inequality (11) corresponding to α i > 0, i = 1, 2. Then the following monotonicity property holds: Proof. Let 0 < α 1 ≤ α 2 and w = u α 1 − u α 2 in Ω. It is sufficient to prove that w + = 0 in Ω. Since w Γ 1 = 0, we have w + ∈ V 0 . We choose v = −w + ∈ V 0 in problem (11) for α 1 , and v = w + ∈ V 0 in problem (11) for α 2 to get By adding the last two inequalities, we have −a(w, w + ) + α 1 Using the coercivity of the form a, we deduce that w + = 0, which completes the proof.
In particular, taking c = 1 we obtain the condition equivalent to the relaxed monotonicity condition with m j = 0, which means that j(x, ·) is convex (see Remark 10).
We conclude that the monotonicity property of Proposition 6 holds for convex potentials, i.e., for variational inequalities. The proof of the monotonicity property for hemivariational inequalities remains an open problem.

Asymptotic behavior of solutions
In this section we investigate the asymptotic behavior of solutions to problem (11) when α → ∞. To this end, we need the following additional hypothesis on the superpotential j.
Proof. First, we prove the estimate on the sequence {u α } in V . We choose v = u ∞ − u α ∈ V 0 as a test function in problem (11) to obtain We observe that j 0 (x, u α ; v) = j 0 (x, u α ; b − u α ) on Γ 3 , and by H(j)(d), we have By the boundedness and coercivity of a, we infer with M > 0, and subsequently where C > 0 is independent of α. Hence, since a(v, v) ≥ 0, from (17), we have It follows from (18) that {u α } remains in a bounded subset of V . Thus, there exists u * ∈ V such that, by passing to a subsequence if necessary, we have Next, we will show that u * = u ∞ . We observe that u * ∈ V 0 because {u α } ⊂ V 0 and V 0 is sequentially weakly closed in V . Let w ∈ K and v = w − u α ∈ V 0 . From (11), we have Next, we use the weak lower semicontinuity of the functional V ∋ v → a(v, v) ∈ R and from (21), we deduce Subsequently, we will show that u * ∈ K. In fact, from (20), by the compactness of the trace operator, we have Passing to a subsequence if necessary, we may suppose that u α (x) → u * (x) for a.e. x ∈ Γ 3 and there exists h ∈ L 2 (Γ 3 ) such that |u α (x)| ≤ h(x) a.e. x ∈ Γ 3 . Using the upper semicontinuity of the function R × R ∋ (r, s) → j 0 (x, r; s) ∈ R for a.e. x ∈ Γ 3 , see Proposition 3(iii), we get Next, taking into account the estimate Consequently, from H(j)(d) and (19), we have Again by H(j)(d), we get j 0 (x, u * ; b − u * ) = 0 a.e. x ∈ Γ 3 . Using (H 1 ), we have u * (x) = b for a.e. x ∈ Γ 3 , which together with (22) implies Next, we prove that u * = u ∞ . To this end, let v := w − u * ∈ K 0 with arbitrary i.e., u * ∈ K is a solution to problem (12). By the uniqueness of solution to problem (12), we have u * = u ∞ and hence u α → u ∞ weakly in V , as α → ∞. From the uniqueness of solution to (12), we also infer that the whole sequence {u α } converges weakly in V to u ∞ . Finally, we prove the strong convergence u α → u ∞ in V , as α → ∞. Choosing v = u ∞ − u α ∈ V 0 in problem (11), we obtain Since u ∞ = b on Γ 3 , by H(j)(d) and the coercivity of the form a, we have Employing the weak continuity of both a(u ∞ , ·) and L, we conclude that u α → u ∞ in V , as α → ∞. This completes the proof.

Continuous dependence result
In this section we provide the result on continuous dependence of solution to problem (11) on the internal energy g and the heat flux q for fixed α > 0.
First, from the compactness of the embedding V into L 2 (Ω) and of the trace operator from V into L 2 (Γ), we obtain the following convergence result.
If g n → g weakly in L 2 (Ω), q n → q weakly in L 2 (Γ 2 ), and v n ∈ V , v n → v weakly in V , then L n (v n ) → L(v), as n → ∞, and there exists a constant C > 0 independent of n such that L n V * ≤ C for all n ∈ N.
The continuous dependence result reads as follows.
Theorem 9. Assume that α > 0 is fixed, L, L n ∈ V * , n ∈ N and H(j) holds. Let u n ∈ V 0 , n ∈ N, be a solution to problem (11) corresponding to L n , and lim L n (z n ) = L(z) for any z n → z weakly in V, as n → ∞.
Then, there exists a subsequence of {u n } which converges weakly in V to a solution of problem (11) corresponding to L. If, in addition, the following hypotheses hold j 0 (x, r; s − r) + j 0 (x, s; r − s) ≤ m j |r − s| 2 for all r, s ∈ R, a.e. x ∈ Γ 3 , (24) where m j ≥ 0, then problem (11) has the unique solution u and u n ∈ V 0 corresponding to L and L n , respectively, and the whole sequence {u n } converges to u in V , as n → ∞.
Proof. Let u n ∈ V 0 be a solution to problem (11) corresponding to L n , and u ∞ ∈ K be the solution to problem (12). We have Hence From hypothesis H(j)(d), since the form a is bounded and coercive, we get and subsequently for all n ∈ N with k 1 , k 2 > 0 independent of n. Hence, {u n } is uniformly bounded in V and also in V 0 . From the reflexivity of V 0 , there exist ξ ∈ V 0 and a subequence of {u n }, denoted in the same way, such that u n → ξ weakly in V 0 , as n → ∞.
We will show that ξ ∈ V 0 satisfies (11). We know that u n ∈ V 0 and a(u n , v) + α Taking the upper limit, we use the weak continuity of a(·, v) and (23) to get By the compactness of the trace operator from V into L 2 (Γ 3 ), we have u n Γ 3 → ξ Γ 3 in L 2 (Γ 3 ), as n → ∞, and at least for a subsequence, u n (x) → ξ(x) for a.e. x ∈ Γ 3 and |u n (x)| ≤ η(x) a.e. x ∈ Γ 3 , where η ∈ L 2 (Γ 3 ). Since the function R × R ∋ (r, s) → j 0 (x, r; s) ∈ R a.e. on Γ 3 is upper semicontinuous, see Proposition 3(iii), we obtain Recalling the estimate Using the latter in (26) entails which means that ξ ∈ V 0 is a solution to problem (11), and completes the first part of the proof.
Next, in addition, we assume (24) and (25). The existence of solution to (11) follows from the first part of the theorem. To prove uniqueness, let u 1 , u 2 ∈ V 0 solve (11). Then taking as test functions u 2 − u 1 ∈ V 0 for u 1 and u 1 − u 2 ∈ V 0 for u 2 , and adding corresponding inequalities, we obtain From the coercivity of the form a and (24), we have Hence, (m a − α m j γ 2 ) u 1 − u 2 2 V ≤ 0, and by the smallness condition (25), we get Hence, we deduce that solutions u, u n ∈ V 0 to (11) are unique, and by (27), we immediately have ξ = u.
Finally, we will show the strong convergence of {u n } to u in V . We choose suitable test functions from V 0 in (11) and (27) to obtain a(u n , u − u n ) + α Adding the two inequalities, we have Using the coercivity of the form a and (24), we get . From hypotheses (23) and (25), we deduce that u n − u V → 0, as n → ∞. Since u ∈ V 0 is unique, we infer that the whole sequence {u n } converges weakly in V to u. This proof is complete.
Remark 10. It is known that for a locally Lipschitz function j : R → R, the condition (24) is equivalent to the so-called relaxed monotonicity condition of the subdifferential for all r i ∈ R, η i ∈ ∂j(r i ), i = 1, 2. The latter was extensively used in the literature, see [18] and the references therein. Condition (24) can be verified by proving that the function R ∋ r → j(r) + m j 2 |r| 2 ∈ R is nondecreasing. An example of a nonconvex function which satisfies the condition (24) is given in Example 11. Note that if j : R → R is convex, then (24) and (28) hold with m j = 0. In fact, by convexity, for all r, s ∈ R which imply j 0 (r; s − r) + j 0 (s; r − s) ≤ 0. Therefore, for a convex function j : R → R, condition (24) or, equivalently, (28) reduces to monotonicity of the (convex) subdifferential, i.e., m j = 0.

Examples
The following examples provide nonconvex and convex functions which satisfies the hypotheses H(j), (H 1 ) and (24).
Example 11. Let j : R → R be the function defined by for r ∈ R with a constant b ∈ R. This function is nonconvex, locally Lipschitz and its subdifferential is given by for all r ∈ R. Hence, we have |∂j(r)| ≤ 1 + 2|b| + 2|r| for all r ∈ R. Moreover, using Proposition 3(ii), one has for all r ∈ R. Thus H(j) is satisfied. By the above formula, we also infer that (H 1 ) is satisfied. Further, we show that condition (24) holds with m j = 1. The condition (24) is equivalent to the relaxed monotonicity of the subdifferential (∂j(r) − ∂j(s))(r − s) ≥ −|r − s| 2 for all r, s ∈ R.
The latter means that i.e., the subdifferential ∂ψ of the function ψ : R → R defined by ψ(r) = j(r) + 1 2 r 2 is monotone. Now, the monotonicity of ∂ψ can be verified using the formula for all r ∈ R. We conclude that H(j), (H 1 ) and (24) are satisfied.
Example 12. (see [15,Example 3]) Let the function j : R → R be given by j(r) = min{j 1 (r), j 2 (r)} for r ∈ R, where j i : R → R are convex, quadratic and such that j ′ i (b) = 0, i = 1, 2. It is known, see [6, Theorem 2.5.1], that ∂j(r) ⊂ conv{j ′ 1 (r), j ′ 2 (r)} for all r ∈ R, so, the subgradient of j has at most a linear growth. Using the monotonicity of the subgradient of convex function, we get for all r ∈ R, i = 1, 2, and, by Proposition 3(iii), we have Hence, we deduce that condition H(j) is satisfied. Similarly, if j 0 (r; b − r) = 0 for all r ∈ R, then λj ′ 1 (r)(b − r) = 0 and (1 − λ)j ′ 2 (r)(b − r) = 0 for all r ∈ R with λ ∈ [0, 1], which is possible when r = b. So, j satisfies also (H 1 ). Further, it is easy to observe that in the case when the graphs of functions j 1 and j 2 have two common points, then the function j is nonconvex. Example 14. Let m 1 , m 2 , r 0 ∈ R be constants such that m 1 ≤ −r 0 < 0 and m 2 ≥ r 0 > 0. Consider the function j : R → R defined by The function j is convex, its subdifferential is given by for all r ∈ R, and its generalized directional derivative has the form for all r, s ∈ R. Hence, we obtain that j 0 (r; b − r) ≤ 0 for all r ∈ R. Similarly, if j 0 (r, b − r) = 0 for all r ∈ R, then r = b. We conclude that j satisfies H(j) and (H 1 ). Moreover, the function j, being convex, satisfies (24) with m j = 0, see Remark 10.
Example 15. We define j : R → R by for r ∈ R with a constant b ∈ R. Then, we have for all r ∈ R, and for all r ∈ R. Thus, j 0 (r; b−r) ≤ 0 for all r ∈ R. Also, we observe that if j 0 (r; b−r) = 0 for all r ∈ R, then r = b. In consequence, the properties H(j) and (H 1 ) are verified. Further, since j is convex, it satisfies (24) with m j = 0, see Remark 10.

Conclusions
We have studied the nonlinear elliptic problem with mixed boundary conditions involving a nonmonotone multivalued subdifferential boundary condition on a part of the boundary. Based on the notion of the Clarke generalized gradient, the variational form of the problem leads to an elliptic boundary hemivariational inequality. We have provided results on existence, comparison of solutions and continuous dependence on the data. Sufficient conditions have been found which guarantee the asymptotic behavior of solution, when the heat transfer coefficient tends to infinity, to a problem with the Dirichlet boundary condition. Under our hypotheses, the proof of the monotonicity property of Theorem 1(iv) for the elliptic hemivariational inequality (7) remains an interesting open problem. We have also given some examples of locally Lipschitz (nondifferentiable and nonconvex) functions to which our results can be applied.