PERTURBATION OF THE FREE BOUNDARY IN ELLIPTIC PROBLEM WITH DISCONTINUITIES

Transcription

1 Electronic Journal of Differential Equations, Vol. 16 (16), No. 13, pp IN: URL: or PERTURBATION OF THE FREE BOUNDARY IN ELLIPTIC PROBLEM WITH DICONTINUITIE ABRI BENID Abstract. We study the discontinuous elliptic problem u = H(u µ) in Ω, u = h on Ω, where Ω is a regular bounded domain of R n, H is the Heaviside function,, µ are a positive real parameters and h is a given function. We prove the existence of solutions, and characterize the free boundaries {x Ω : u(x) = µ} using the perturbation of the boundary condition and smooth boundary of the domain. 1. Introduction Partial differential equations with discontinuous nonlinearities arise in models from many concrete problems in mathematical physics like those of combustion theory, porous media, plasma physics. In this article, we study the existence of solutions for the problem u = H(u µ) in Ω, u = h on Ω, where Ω is a smooth bounded domain of R n, H is the Heaviside function { 1 if t H(t) = if t <, (1.1) h is a given function, and, µ are a positive real parameters. This problem can be reformulated as an equivalent free boundary problem: Find u C (Ω \ w) C 1 (Ω) such that u = in w, u = in Ω \ w, u = h on Ω, (1.) where w = {x Ω : u(x) > µ} and w is the free boundary to be determined. The set {x Ω : u(x) = µ}, dividing the domain Ω into two (or more) regions where u = or u = is satisfied in the classical sense. 1 Mathematics ubject Classification. 34R35, 35J5. Key words and phrases. Perturbation; implicit function theorem; free boundary problem. c 16 Texas tate University. ubmitted April 14, 16. Published June 7, 16. 1

2 . BENID EJDE-16/13 To the best of our knowledge, no investigation has been devoted to establishing the existence of solutions to such problem when Ω is a general domain using perturbation methods. Here, we give a positive answer to this study. The problem (1.1) was investigated in the variational context when h = by Ambrosetti and Badial [3] and the extension to the results for p-laplacian operator (p > 1) has been studied by Arcoya and Calahorrano in [5]. Further works with the Neumann boundary conditions can be found in [11, 14]. When h does not vanish identically, and the domain Ω is the unit ball of R n, the existence of positive solutions and the behavior of the corresponding free boundaries with respect to h has been discussed by Alexander in [1]. In the case n =, Alexander and Fleishman [] studied the problem (1.1) when Ω is the unit square. Recently, the authors in [6] and [7] studied the more general following problem when Ω is the unit ball of R n, u = f(u)h(u µ) in Ω, u = h on Ω, where f is a given function. The methods used hinge on the parametrization of the free boundary which is the unknown of our problem. This technique reduces the study to solve a nonlinear integral equation and allows us to obtain positive result regarding the solvability of the equation of the free boundary under perturbation of the boundary without the additional requirement of the regularity of the boundary of the initial domain. This will be the approach of this paper. In general, the natural way to deal with the discontinuous elliptic problems are the variational methods, see for instance [9, 1]. But one soon realizes that the study entails serious difficulties mainly to characterize the variation of the free boundary. Our goal is different, having a solution of the given problem, we study the effect on the solutions under perturbations of the boundary conditions and a smooth boundary of the domain Ω. We point out also that in this work, we can not use the symmetry of solutions as in [1, 6, 7]. Hence, we have uses the perturbation techniques to overcome the encountered difficulties. For more works on the discontinuous elliptic problems, we invite the reader to consult [1, 3, 5, 6, 7, 15] and the references given here. Because the nonlinearity has a discontinuity at u = µ, so, a suitable concept of solution is needed. We say that a function u W,p (Ω), (p > 1) is a solution of problem (1.1) if u = H(u µ) a.e in Ω and the trace of u on Ω is equal to h. We need the following assumptions in this work. (H1) Assuming that the boundary Ω of the domain Ω can be parameterized as R + β(θ) where β C (), θ and R > where is the unit sphere. (H) For n =, µ > 4e R. and for n 3, µ > M n(r) := 1 n(n ) R ( n ) n ( n ) n n (H3) The function h is small enough with h < µ where h = max x Ω h(x).

3 EJDE-16/13 PERTURBATION OF THE FREE BOUNDARY 3 This article is organized as follows: ection collects some known results for the problem in a ball giving only slight information on the different methods of proof. ection 3 contains the statement of the essential result. ection 4 provides an approach for studying the problem (1.1). In ection 5, we treat the regularity of free boundary and finally, we give an appendix which contains some useful results.. Existence results in a ball In this section, we consider the problem u = H(u µ) u = h in B(, R), on B(, R), We assume that the function h satisfies (H3). Theorem.1 ([7]). uppose that there exists µ > such that µ > 4e, R for n =, µ > M n(r), for n 3, (.1) then the problem (.1) has at least two positive solutions and the free boundaries are analytic hypersurfaces. Remark.. In [7], we have treat the case n 3 and the more general problem u = f(u)h(u µ) in B(, 1), u = h on B(, 1), There the assumptions are (H4) The function f is k-lipstchitzian, non-decreasing, positive and there exist two strictly positive constants k, β > such that f(s) ks + β with k < min{ 1, 1}, where 1 is the first eigenvalue of under homogeneous Dirichlet boundary conditions. (H5) The function f is differentiable and constant on the interval of the form [, c] where c > β n k and the function h is small enough, h < µ where h = max x B(,1) h(x). (H6) There exists µ > such that f(µ) µ > M n(n ) n = ( n ) n ( n ), for n 3. n n In this paper, we put f(u) = to clarify the obtained the results. For n =, the treatment is similar with the adequate modification. For the convenience of the reader, we give some calculus related to the case Ω = B(, R) in the appendix. The proof of Theorem.1 is based on the transformation of our problem into an equivalent nonlinear integral equation which is solved with respect to the unknown free boundary {x Ω : u(x) = µ} by the application of the implicit function theorem. We can find the complete proof in [7]. Here, we give some ideas of demonstration that will be used throughout this paper. We will denoted the solution of (.1) by u and we consider the problem called the reduced problem u = H(u µ) u = on B(, R), in B(, R), (.)

4 4. BENID EJDE-16/13 We remark that since the Heaviside step function is monotone and not decreasing, the result of Gidas, Ni and Nirenberg [1] shows that all positive solutions of (.) are radial. Hence, we look for the free boundary in the form {(r, θ), θ } for some r (, R) and obtain all radial solutions of (.) by finding u, and r so that the differential equation r 1 n r (rn 1 u r ) = < r < r r 1 n u r (rn 1 r ) = r < r < R, (.3) u u (R) =, r () = is satisfied with the following transmission conditions on the free boundary u u (r ) = µ and r (r ) = u r (r + ), where u r (r ) denotes the left derivative of u and u r (r + ) denotes the right derivative at the value r = r. Using assumption (H), the resolution of (.3) with the previous conditions gives two r (, R) says r 1, r. Hence, a simple calculus shows that the free boundaries are spheres with radii r 1, r (, R). Remark.3. Let r denote one of the values r 1 or r. Using hypothesis (H), r Re 1/ for n = and ( ) 1 n r nr n for n 3. ee the Appendix. Now, to study the problem u = H(u µ) u = h on B(, R), in B(, R), (.4) we use the effect of perturbation on the solution in the boundary values. More precisely, when h, we look for the free boundary in the form r + b(θ), θ, where b(θ) is the perturbation caused by h. Let w = {(r, θ) (, R), r < r + b(θ), θ } Now, we denote by χ w the characteristic function of w. In the following result, we formulate a nonlinear equation for the unknown function b and we prove that by solving it, we can solve the problem (.4). Proposition.4 ([7]). Under assumption (H), the problem u = χ w (r, θ) in B(, R), u = h on B(, R) (.5) has a unique solution u C 1,α (Ω, R) with α = 1 n p. Moreover if u (r +b(θ), θ) = µ with h < µ, then u is a solution of (.4). As in [7], we conclude with problem (.4) by showing the existence of the function b via the implicit function theorem. We can represent the solution of (.5) in integral form, and we check the hypothesis of previous theorem. For more details, we refer the reader to [7]. In this article, we denote r + b(θ) := b (θ), θ.

5 EJDE-16/13 PERTURBATION OF THE FREE BOUNDARY 5 3. Perturbation of the free boundary In this section, we are concerned with problem (1.1). We introduce some sets necessary for the study using the perturbation of boundary. We assume that Ω is diffeomorphic to a ball. Then, we can construct a curvilinear parametrisation (r, θ) in a neighborhood of any set in Ω, r (, R] and θ. Let Ω β be a perturbation of B(, R) in the sense that the boundary Ω β of every smooth open set Ω β close to the ball B(, R) can be described by R + β(θ), θ and β C (). Hence, from now, Ω = Ω β, and Ω = B(, R). We denote also by Ω β the admissible perturbation of the ball B(, R) and we define the set of admissible surfaces in Ω β by β = {f C() : (f(θ), θ) Ω β for θ }. Now, for a function ψ β, we define the set Ω β,ψ = {(r, θ) Ω β, r < ψ(θ)}. We seek a solution in W,p (Ω β ), p > 1, then the boundary value function h which is a trace of W,p (Ω β ) function will be taken in the set H = {h W 1 p,p ( Ω β, R), p > n}. We denote by χ Ωβ,ψ the characteristic function of Ω β,ψ, then, we have the following result. Proposition 3.1. Assume (H1), (H) and that Then the problem µ > 4e R for n =, µ > M n(r), for n 3. u = χ Ωβ,ψ in Ω β, u = h on Ω β (3.1) has a unique solution u C 1,α (Ω β, R) with α = 1 n p. Moreover, if u(ψ(θ), θ) = µ with h < µ then u is a solution of (1.1). Proof. First, we see that χ Ωβ,ψ L p (Ω β ), p > 1. From [13, Theorem 9.15], there exists a unique solution of (3.1) in W,p (Ω β ). For p > n, W,p (Ω β ) C 1,α (Ω β, R) with α = 1 n p. Now, we remark that u satisfies u = in Ω β,ψ, u = in Ω β \ Ω β,ψ, u = h on Ω β. If we prove the existence of a function ψ such that u(ψ(θ), θ) = µ, then u will be a solution of u = in Ω β,ψ, u = µ on Ω β,ψ u = in Ω β \ Ω β,ψ, u = h on Ω β.

6 6. BENID EJDE-16/13 In Ω β,b, the function u satisfies The maximum principle implies u = in Ω β,ψ, u = µ on Ω β,ψ min Ω β,ψ u = min Ω β,ψ u = µ. Hence, u > µ in Ω β,ψ. In Ω β \ Ω β,b, we have As h < µ, then u = in Ω β \ Ω β,ψ, u = µ on Ω β,ψ, u = h on Ω β. max u = max u = µ Ω β \Ω β,ψ Ω β,ψ and consequently, u < µ in Ω β \ Ω β,ψ. Therefore, the function u satisfies u = H(u µ) in Ω β, u = h on Ω β. To conclude with the existence of solutions of problem (1.1), we need only to show the existence of the function ψ satisfying the equation u(ψ(θ), θ) = µ. Now, since the solution u depend on the domain Ω β, then we can not use the local methods directly to prove the existence of ψ. The variation of the domain Ω β suggests to use an adequate transformation which maps the changing domain into a fixed domain and solves the governing equations in the mapped domain. To exclude this difficulties, we proceed as the following. To each admissible perturbation (Ω, Ω β ) of the domain Ω correspond a transformation T β of the domain Ω β onto the initial domain Ω. T β : Ω β Ω, (r, θ) (r, θ) = (r + r β R, θ) where (r, θ) is the coordinates in Ω β and (r, θ) the coordinates in Ω. For a small β, the transformation T β is a diffeomorphism of class C of the domain Ω β into Ω. The mapping T β maps the class of admissible surfaces β into. Hence, Using the relation T β (θ, ψ(θ)) = (θ, f(θ)) where f(θ). u(t β (r, θ)) = u(r, θ), problem (3.1) with u(ψ(θ), θ) = µ is equivalent to the problem with the equation L β u = χ Ωβ,f in Ω, u = h on Ω (3.) u(f(θ), θ) = µ, where L β is a linear operator with continuous coefficients depending on β. The following lemma gives the exact expression of L β.

7 EJDE-16/13 PERTURBATION OF THE FREE BOUNDARY 7 Lemma 3.. The linear operator L β is given by L β = + δ β where δ β has the form δ β = β R ( + β R ) (n 1) β + r r R r + 1 [ r β ( 1 β [a r ij (θ) + R θ j r θ i R(1 + β R ) θ i + r β R θ j θ i r + r β ] + b i (θ) [ r R θ i θ j r R and = (n 1) + r r r + 1 r θ when θ is the Beltrami-Laplace operator. r + r β ) R θ i r ] ] r Remark 3.3. This formulas are obtained from the Laplacian in polar coordiantes. More precisely, from = (n 1) + r r r + 1 r θ (3.3) where θ = a ij (θ) θ i θ j + b i (θ) θ i. Proof of Lemma 3.. First, we have β θ i u(r, θ) = u(r, θ) where r = r + r β R, u r = u r r r + u θ θ r = u r + u β r R := Q, u r = Q r r r + Q θ θ r = u( β ) u β ( β ) r R r R R = u r + β ( β ) u + R R r. u = u r + u θ i = r β u θ i r θ i θ i θ i R θ i r + u := Z θ i u = Z = Z r + Z θ i θ j θ i θ j r θ j θ j θ j = r β ( u + r β u ) + Z R θ j r θ i R θ i r θ i = r β ( u 1 β u + R θ j r θ i R(1 + β R ) θ i r + r β u ) R θ i r + u + r β u θ i θ j R θ j θ i r + r β u R θ i θ j r. If we replace the previous terms in the expression (3.3), we find the formulas of L β. Thus, to solve problem (1.1), it is sufficient to prove the existence of function f such that u(f(θ), θ) = µ.

8 8. BENID EJDE-16/13 In fact, using the implicit function theorem, we prove that the equation u(f(θ), θ) µ = is uniquely solvable in a given small neighborhood. following section. (ee Theorem 4. below). This is the subject of the 4. olvability in the neighborhood of smooth free boundary To each function f, we associate the solution u of problem (3.). We have u δ β u = χ Ωβ,f in Ω, u = h on Ω (4.1) The solution u corresponding to (4.1) has an integral representation which is well defined [6, Theorem 4.1], given by u(x) = P (x, y)h(y)d G(x, y)χ Ωβ,f dy G(x, y)δ β (u)(y)dy. Ω Ω Now, we consider polar coordinates and we define the operator J : H D C(, R) by J(h, f, β)(θ) = u(f(θ), θ) µ, where D is the neighborhood of zero in C(); i.e., J(h, f, β)(θ) = P (f(θ), θ, θ )h(θ )dθ G(f(θ), θ, r, θ )χ Ωβ,f (r, θ )dr dθ Ω Ω δ β (u)(r, θ )G(f(θ), θ, r, θ )dr dθ µ. Ω Then we have the following result. Lemma 4.1. The operator J is continuously differentiable with respect to the second variable. Proof. Let D j J denote the Frechet derivative of J, with respect to the variable of order j (j = 1,, 3). Let ϕ(θ) be a small perturbation of f(θ), then J(h, f + ϕ, β)(θ) J(h, f, β)(θ) D Jϕ(θ) = P (f(θ) + ϕ(θ), θ, θ )h(θ )dθ Ω G(f(θ) + ϕ(θ), θ, r, θ )χ Ωf+ϕ (r, θ )dr dθ Ω l β (u)(r, θ )G(f(θ) + ϕ(θ), θ, r, θ )dr dθ P (f(θ), θ, θ )h(θ )dθ Ω Ω + G(f(θ), θ, r, θ )χ Ωf (r, θ )dr dθ + l β (u)(r, θ )G(f(θ), θ, r, θ )dr dθ Ω Ω + G(f(θ) + ϕ(θ), θ, r, θ )χ Ωf dr dθ Ω G(f(θ) + ϕ(θ), θ, r, θ )χ Ωf dr dθ Ω P Ω r (f(θ ), θ, θ )h(θ )ϕ(θ)dθ

9 EJDE-16/13 PERTURBATION OF THE FREE BOUNDARY Hence, f (r n 1 G ) r (f(θ), θ, r, θ )ϕ(θ)dr dθ f n 1 (θ)g(f(θ), θ, f(θ ), θ )ϕ(θ )dθ Ω l β (u)(r, θ ) G r ((fθ), θ, r, θ )ϕ(θ)dθ where J(h, f + ϕ, β)(θ) J(h, f, β)(θ) D Jϕ(θ) = I 1 + I + I 3 + I 4 I 1 = P (f(θ) + ϕ(θ), θ, θ )h(θ )dθ P (f(θ), θ, θ )h(θ )dθ Ω Ω P Ω r (f(θ ), θ, θ )h(θ )ϕ(θ)dθ I = G(f(θ) + ϕ(θ), θ, r, θ )χ Ωf dr dθ Ω + G(f(θ), θ, r, θ )χ Ωf (r, θ )dr dθ Ω f + (r n 1 G ) r (f(θ), θ, r, θ )ϕ(θ)dr dθ I 3 = l β (u)(r, θ )G(f(θ) + ϕ(θ), θ, r, θ )dr dθ Ω + l β (u)(r, θ )G(f(θ), θ, r, θ )dr dθ Ω + l β (u)(r, θ ) G Ω r ((fθ), θ, r, θ )ϕ(θ)dθ I 4 = G(f(θ) + ϕ(θ), θ, r, θ )[χ Ωf χ Ωf+ϕ ]dr dθ Ω + f n 1 (θ)g(f(θ), θ, f(θ ), θ )ϕ(θ )dθ = dθ [ f+ϕ (r ) n 1 G(f(θ) + ϕ(θ), θ, r, θ )dr f ] f n 1 (θ)g(f(θ), θ, f(θ ), θ )ϕ(θ ) Using Taylor s theorem, we have I 1, I, I 3, I 4 = o( ϕ ) when ϕ +. Hence, D J(h, f, β)ϕ(θ) = u (f(θ), θ)ϕ(θ) r f n 1 (θ )G(f(θ), θ, f(θ ), θ )ϕ(θ )dθ. Now, to solve the equation J(h, f, β)(θ) =, for θ, in the neighborhood of (h, b, ). we need to know the invertibility of J. We recall the reader that we consider here a general domain and that causes some difficulties to build clear conditions so that the operator D J(h, b, ) either invertible of into C() or not. We refer to [6] to see explicit conditions such that the proper operator is

10 1. BENID EJDE-16/13 invertible when the domain is a ball. Let the operator D J(h, b, )ϕ(θ) = u r (b (θ), θ)ϕ(θ) First, note that u satisfies the equation u = H(u µ) u = h b n 1 (θ )G(b (θ), θ, b (θ ), θ )ϕ(θ )dθ. in B(, R), on B(, R). Then, when h is small, the solution u is close in C 1,α to the solution of u = H(u µ) u = on B(, R). in B(, R), ince, u r (r, θ) < for r (, R), θ, then u r (b (θ), θ) < for b caused by the small perturbation h. For the other part, let Kϕ(θ) = b n 1 (θ )G(b (θ), θ, b (θ ), θ )ϕ(θ )dθ. Let L b () the space of functions belonging to the space L () with the inner product u, v = b n 1 u(θ)v(θ)dθ. (4.) We remark that the operator K is negative definite in the space L b (), (The Green function G is negative). ince, the function b is bounded in, the inner product (4.) is equivalent to the standard product in the space L (). In [6], we have proved that the operator D J(, r, ) is invertible (note that b = r when h = on Ω). In this case, we use the explicit eigenvalue of the operator K to conclude the nondegeneracy of the operator D J(, r, ). For more details, see [6]. Finally, using the fact that D J(, r, ) is invertible in L r (), we can be sure that D J(h, b, ) is invertible in L b () with the previous inner product. The preceding argument shows that D J(h, b, ) is invertible in L (). Hence, using the implicit function theorem, we have the following result. Theorem 4.. Under assumptions (H1) (H3), there exist a neighborhood V of (h, ) in H C () and a continuous mapping B : V C() such that (i) B(h, ) = b, (ii) J(h, B(h, β), β) =. We recall that h is a given function satisfying (H3) and b (θ) := r +b(θ), where r (, 1), b C 1,α (), α = 1 n p and θ. Proof of Theorem 4.. Because the operator J is invertible in the neighborhood of (h, b, ) and J(h, b, )(θ) = for θ, then the implicit function theorem implies the existence of function B depending on h and β such that B(h, β) satisfies J(h, B(h, β), β) =.

11 EJDE-16/13 PERTURBATION OF THE FREE BOUNDARY Regularity of the free boundary Theorem 5.1. If h and β are small enough, then the free boundary {x Ω/u(x) = µ} is an analytic hypersurface. Proof. First, let u be a solution of (4.1) and let Γ : = {x Ω/u(x) = µ} = {(r, θ) (, R), u(r, θ) = µ} = {(f(θ), θ), for θ }. when β is sufficiently close to in C () and f is close to b in C(), then the solution u is close to u in C 1,α (Ω), α (, 1). ince u r (b (θ), θ), then u (f(θ), θ) for θ. r The implicit function theorem gives that f C 1,α (). Hence, the free boundary Γ is an hypersurface of class C 1,α. Now, we conclude that the free boundary Γ is analytic by the application of Hodograph transformation. This method was using by the author in [7]. Final remarks. (1) o far, the problem (1.1) has been studied for a few class of domain, in particular the ball for uses the notion of symmetry. Hence, the main results of this paper remain true when Ω is a ring shaped domain. () The advantage of the perturbation procedure described in this paper is the the explicit formula of behavior of the free boundary with respect to h. For example, in dimension, taking Ω := {(r, θ)/r < + sin(θ), θ } and h(x) = cos(x), we can give an explicit formulation of the solution of problem (1.1) and the shape of free boundary. (3) When n = 1, the problem (1.1) becomes the following second order differential equation. For example, we can see easily that under a suitable conditions, the problem u = H(u µ) for x a, u(±a) =, (a R + (5.1) ). admits multiple solutions. ee [15]. In the other part, the problem u = H(µ u) for x a, u(±a) =, (a R + ). (5.) has a unique positive solution when the set {x Ω, u(x) = µ} is a given segment included in x a. In fact, the form of the discontinuity in the second member has a surprising effects on the set of solutions and their free boundaries. In higher dimensions, the result stays true also,we can see [8] for more details. (4) The regularity of the free boundary is preserved after perturbations. One of the important questions is when will the free boundary develop singularities. (5) It is also interesting to studied the bifurcation phenomenon. In the case when Ω is the unit ball, the study is in [7].

12 1. BENID EJDE-16/13 6. Appendix In this section, we give the proof that the condition in assumption (H6) in Remark. becomes µ > 4e, for n =, R µ > M n(r) := 1 n(n ) R ( n ), for n 3, n ( n n ) n as in Theorem.1 when Ω = B(, R) and f(u) =. Now, define a function u that satisfies u = H(u µ) in Ω, (6.1) u = on Ω. ince H(u µ) is monotone and not decreasing, the result of [1] shows that all positive solutions are radial. Hence, by the maximum principle, a positive radial solution u take on the value µ at only one value say r. Then we obtain all radial solutions of (.) by finding, u and r so that the two problems and r 1 n u (rn 1 r r ) = for < r < r, u () =, u(r ) = µ r 1 n u (rn 1 r r ) = for r < r < R, u(r) =, u(r ) = µ are satisfied. When n = and from the problems (6.), (6.3), we have { 4 u(r) = (r r ) + µ < r r, µ ln(r/r) ln(r/r ) r r < R. ince u C 1,α (Ω), we have (6.) (6.3) u r (r ) = u r (r+ ). We obtain r µ = r (ln(r/r )) which implies µ = r ln(r/r ) := g(r ). (6.4) Now, by considering the function g(ρ) in (, R), we have that f has a maximum value R e at ρ = Re 1. Hence, the equation (6.4) has two roots when µ > 4e R. imilarly, we treat the case n 3. The solution u is given by n (r r ) + µ < r r, u(r) = µr n r n R The transmission conditions imply that µr n n r n R n r r < R. u r (r+ ) = u r (r ),

13 EJDE-16/13 PERTURBATION OF THE FREE BOUNDARY 13 so Hence, ( n)µr 1 n r n R = r n n. n(n ) = µ r := g(r rn R n It follows that the function g has a minimum value reached at the point ( nr ) 1 n n roots when µ > M n(r). M n (R) = 1 n(n ) R ( n ) n ( n ) n n (6.5) which implies that the equation (6.5) has two Acknowledgments. The author would like to thank the professor B. Abdellaoui for the discussion and the suggestions about the problem in this paper. I want to thank also the anonymous referee for valuable suggestions and comments. References [1] R. Alexander; A discontinuous nonlinear eigenvalue/free boundary problem, Math. Meth in the Appl. ci., 4, (198), [] R. Alexander, B. Fleishman; Perturbation and bifurcation in a free boundary problem. J. Differential Equations, 45, no. 1, (198), [3] A. Ambrosetti, M. Badiale; The dual variational principle and elliptic problems with discontinuous nonlinearities, J. Math. Anal. Appl, 14, No., (1989), [4] A. Ambrosetti, P.H. Rabinowitz; Dual variational methods in critical point theory and applications, J. Funct. Anal, 14, (1973), [5] D. Arcoya, M. Calahorrano; ome discontinuous problems with quasilinear operator, J. Math. Anal. Appl. 187, (1994), [6]. Bensid,. M. Bouguima; On a free boundary problem, Nonlinear Anal. T.M.A, 68, (8), [7]. Bensid,. M. Bouguima; Existence and multiplicity of solutions to elliptic problems with discontinuities and free boundary conditions, Electron. J. Differential Equations., 56 (1), [8]. Bensid,. M. Bouguima; A note on discontinuous problem with a free boundary, J. Egyptian Math. oc., 19 (11), [9] G. Bonanno, P. Candito; Non-differentiable functionals and applications to elliptic problems with discontinuous nonlinearities. J. Differential Equations., 44, (8), no. 1, [1] G. Bonanno, P. Candito, D. Motreanu; A coincidence point theorem for sequentially continuous mappings. J. Math. Anal. Appl., 435, (16), no. 1, [11] P. Candito; Infinitely many solutions to the Neumann problem for elliptic equations involving the p-laplacian and with discontinuous nonlinearities, Proc. Edinb. Math. oc., II. er., 45 () (), [1] B. Gidas, W. N. Ni, L. Nirenberg; ymmetry and related properties via maximum principle. Commun. Math. Phys 68, (1979), [13] D. Gilbarg, N. Trudinger; Elliptic Partial Differential Equations of econd Order. pringer- Verlag, New York, [14]. A. Marano, D. Motreanu; Infinitely many critical points of non-differentiable functions and applications to a Neumann-type problem involving the p-laplacian, J. Differential Equations 18, 1 (), [15] C. A. tuart; Differential equations with discontinuous nonlinearities, Arch. Ration. Mech. Anal., 63, 1 (1976),

14 14. BENID EJDE-16/13 abri Bensid Department of Mathematics, Faculty of ciences, University of Tlemcen, B.P. 119, Tlemcen 13, Algeria address: edp