Abstract
In this paper we analyze iterations of the obstacle problem for two different operators. We solve iteratively the obstacle problem from above or below for two different differential operators with obstacles given by the previous functions in the iterative process. When we start the iterations with a super or a subsolution of one of the operators this procedure generates two monotone sequences of functions that we show that converge to a solution to the two membranes problem for the two different operators. We perform our analysis in both the variational and the viscosity settings.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
The main goal in this paper is to find solutions to the two membranes problem as limits of sequences obtained by iterating the obstacle problem. Let us first describe the obstacle problem (from above or below) and the then two membranes problem (as described in [8]).
1.1 The obstacle problem
The obstacle problem is one of the main problems in the mathematical study of variational inequalities and free boundary problems. The problem is to find the equilibrium position of an elastic membrane whose boundary is fixed, and which is constrained to lie above a given obstacle. In mathematical terms, given an operator L (notice that here we can consider fully nonlinear problems of the form \(L u = F(D^2u, Du, u, x)\)) that describes the elastic configuration of the membrane, a bounded Lipschitz domain \(\Omega \) and a boundary datum f, the obstacle problem from below (here solutions are assumed to be above the obstacle) reads as
or equivalently
We will denote by
the solution to (1.1).
When the operator L is in divergence form and is associated to an energy functional E(u) the problem becomes a variational problem and a solution can be obtained minimizing E in the set of functions in the appropriate Sobolev space that are above the obstacle and take the boundary datum. We refer to Sect. 2 for details.
The obstacle problem can be also stated as follows: we look for the smallest supersolution of L (with boundary datum f) that is above the obstacle. This formulation is quite convenient when dealing with fully nonlinear problems using viscosity solutions. We will assume here that the problem (1.1) has a unique continuous viscosity solution (for general theory of viscosity solutions we refer to [9] and [12]). This is guaranteed if L has a comparison principle and one can construct barriers close to the boundary so that the boundary datum f is taken continuously. See Sect. 3.
We can also consider the obstacle problem from above (here solutions are assumed to be below the obstacle)
or equivalently
In this case, the obstacle problem can be viewed as follows: we look for the largest subsolution of L (with boundary datum g) that is below the obstacle. We will denote by
the solution to (1.2).
For general references on the obstacle problem (including regularity of solutions, that in some cases are proved to be optimally \(C^{1,1}\)) we just mention [5, 6, 13, 24], the survey [21] and references therein.
1.2 The two membranes problem
Closely related to the obstacle problem, one of the systems that attracted the attention of the PDE community is the two membranes problem. This problem models the behaviour of two elastic membranes that are clamped at the boundary of a prescribed domain, they are assumed to be ordered, one membrane is above the other, and they are subject to different external forces (the membrane that is on top is pushed down and the one that is below is pushed up). The main assumption here is that the two membranes do not penetrate each other (they are assumed to be ordered in the whole domain). This situation can be modeled by two obstacle problems; the lower membrane acts as an obstacle from below for the free elastic equation that describes the location of the upper membrane, while, conversely, the upper membrane is an obstacle from above for the equation for the lower membrane. When the equations that obey the two membranes have a variational structure this problem can be tackled using calculus of variations (one aims to minimize the sum of the two energies subject to the constraint that the functions that describe the position of the membranes are always ordered inside the domain, one is bigger or equal than the other), see [27]. On the other hand, when the involved equations are not variational the analysis relies on monotonicity arguments (using the maximum principle). Once existence of a solution (in an appropriate sense) is obtained a lot of interesting questions arise, like uniqueness, regularity of the involved functions, a description of the contact set, the regularity of the contact set, etc. See [7, 8, 26], the dissertation [28] and references therein. More concretely, given two differential operators \(L_1 (D^2u, D u, u, x)\) and \( L_2 (D^2v, D v, v,x)\) the mathematical formulation the two membranes problem (with Dirichlet boundary conditions) is the following:
Definition 1.1
A pair (u, v) is called a solution to the two membranes problem if it solves
With our previous notation, this system can be written as
Remark that, in general, the two membranes problem as stated in Definition 1.1 does not have uniqueness. To see this, just take a solution to \(L_1(D^2 u)=0\) with \(u|_{\partial \Omega } = f\) and v the solution to the obstacle problem for \(L_2(D^2 v)\) with u as obstacle from above, that is, \(v=\overline{\text {O}}(L_2,u, g)\). One can easily check that the pair (u, v) is a solution to the general formulation of the two membranes problem stated Definition 1.1. Analogously, one can take a solution to \(L_2(D^2 v)=0\) with \(v|_{\partial \Omega } = g\) and u the solution to the corresponding obstacle problem, \(u=\underline{\text {O}}(L_1,v, f)\), to obtain another solution to the two membranes problem. In general, these two pairs do not coincide. For example, consider \(\Omega \) as the interval (0, 1), \(f\equiv 1\), \(g\equiv 0\) on \(\partial \Omega \) and the operators \(L_{1}(u'')=u''-10\) and \(L_{2}(v'')=v''+2.\) The solution to \(L_{1}(u'')=0\) with \(u(0)=u(1)=1\) is \(\widehat{u}(x)=- 5x(1-x)+1\) and the solution to \(L_{2}(v'')=0\) with \(v(0)=v(1)=0\) is \(\tilde{v}(x)=x(1-x)\) for \(x\in (0,1)\). Let \(\widehat{v}=\overline{\text {O}}(L_2,\widehat{u}, g)\) and \(\tilde{u}=\underline{\text {O}}(L_1,\tilde{v}, f)\). Both pairs \((\widehat{u},\widehat{v})\) and \((\tilde{u},\tilde{v})\) are solutions to the two membranes problem but they are different. Since \(\widehat{u}\) is a solution and \(\tilde{u}\) is a supersolution to \(L_{1}\) with the same boundary datum, by the comparison principle, we get \(\widehat{u}\le \tilde{u}\). On the other hand, since \(\tilde{u}\) is the solution to the obstacle problem for \(L_1(u'')\) with \(\tilde{v}\) as obstacle from below, \(\tilde{u}\ge \tilde{v}\). Since \(\widehat{u}\) and \(\tilde{v}\) are not ordered, we obtain that \(\widehat{u}\not \equiv \tilde{u}\).
The two membranes problem for the Laplacian with a right hand side, that is, for \(L_1(D^2u)=-\Delta u +h_1\) and \(L_2(D^2v)=-\Delta v-h_2\), was considered in [27] using variational arguments. Latter, in [7] the authors solve the two membranes problem for two different fractional Laplacians of different order (two linear non-local operators defined by two different kernels). Notice that in this case the problem is still variational. In these cases an extra condition appears, namely, the sum of the two operators vanishes, \( L_1 (D^2 u) + L_2 (D^2 v) =0\), inside \(\Omega \). Moreover, this extra condition together with the variational structure is used to prove a \(C^{1,\gamma }\) regularity result for the solution.
The two membranes problem for a fully nonlinear operator was studied in [7, 8, 26]. In particular, in [8] the authors consider a version of the two membranes problem for two different fully nonlinear operators, \(L_1(D^2 u)\) and \(L_2(D^2 v)\). Assuming that \(L_1\) is convex and that \( L_2(X) = -L_1(-X),\) they prove that solutions are \(C^{1,1}\) smooth.
We also mention that a more general version of the two membranes problem involving more than two membranes was considered by several authors (see for example [1, 10, 11]).
1.3 Description of the main results
As we mentioned at the beginning of the introduction, our main goal is to obtain solutions to the two membranes problems as limits of iterations of the obstacle problem. For iterations of the obstacle problem in a different context we refer to [4].
Let us consider two different operators \(L_1\) and \(L_2\) (with boundary data f and g respectively, we assume that \(f>g\)) and generate two sequences iterating the obstacle problems from above and below. Given an initial function \(v_{0}\), take \(u_0\) as the solution to the obstacle problem from below for \(L_1\) with boundary datum f and obstacle \(v_0\), that is,
Now, we use this \(u_0\) as the obstacle from above for \(L_2\) with datum g and obtain
We can iterate this procedure (solving the obstacle problem for \(L_1\) or \(L_2\) with boundary data f or g and obstacle the previous u or v) to obtain two sequences \(\{u_n\}_n\), \(\{v_n\}_n\), given by
Our main goal here is to show that, when the initial function \(v_0\) is a subsolution for \(L_2\) with boundary datum g, then both sequences of functions \(\{u_{n}\}_{n}\), \(\{v_{n}\}_{n}\) are nondecreasing sequences that converge to a limit pair that gives a solution of the two membranes problem, i.e., there exists a pair of functions \((u_{\infty },v_{\infty })\) such that
and the limit functions satisfy
An analogous result can be obtained when we start the iteration with \(u_0\) a supersolution to \(L_1\) with boundary datum f and consider \(\{u_n\}_n\), \(\{v_n\}_n\), given by
In this case the sequences are non-increasing sequences and also converge to a solution to the two membranes problem.
We will study this iterative scheme using two very different frameworks for the solutions. First, we deal with variational methods and understand solutions in the weak sense (this imposes that \(L_1\) and \(L_2\) must be in divergence form) and next we deal with viscosity solutions (allowing \(L_1\) and \(L_2\) to be general elliptic operators).
Variational operators. To simplify the notation and the arguments, when we deal with the problem using variational methods, we will concentrate in the particular choice of \(L_1\) and \(L_2\) as two different \(p-\)Laplacians; that is, we let
where \(\Delta _{p}w = \hbox {div} (|\nabla w|^{p-2} \nabla w)\), \(\Delta _{q}w = \hbox {div} (|\nabla w|^{q-2} \nabla w)\) are two different \(p-\)Laplacians and \(h_{p}\), \(h_{q}\) are two given functions. These operators are associated to the energies
These energies are naturally well-posed in the Sobolev spaces \(W^{1,p} (\Omega )\) and \(W^{1,q} (\Omega )\), respectively.
In this variational framework the obstacle problem can be solved minimizing the energy in the appropriate set of functions. In fact, the solution to the obstacle problem from below, \(u=\underline{\text {O}}(\mathcal {L}_{p},\varphi ,f )\), can be obtained as the minimizer of the energy \(E_{p}\) among functions are constrained to be above the obstacle, that is,
with
The minimizer in (1.3) exists and is unique provided the set \(\underline{\Lambda }^{\,p}_{f,\varphi }\) is not empty. We refer to Sect. 2 for details and precise definitions.
Due to the fact that the problem satisfied by the limit of the sequences constructed iterating the obstacle problem from above and below does not have uniqueness (see the comments after Definition 1.1) we need to rely on monotonicity of the sequences to obtain convergence (rather than use compactness arguments that provide only convergence along subsequences).
We notice at this point that, in this variational setting, one can solve the two membranes problem just minimizing the total energy \( E(u,v) = E_p (u) + E_q(v) \) in the subset of \(W^{1,p}(\Omega )\times W^{1,q}(\Omega )\) given by \(\{(u,v): u-f\in W^{1,p}_{0}(\Omega ), v-g\in W^{1,q}_{0}(\Omega ), u\ge v \text { a.e }\Omega \}\). The fact that there is a unique pair that minimizes E(u, v) follows from the strict convexity of the functional using the direct method of calculus of variations. One can check that the minimizing pair is in fact a solution to the two membranes problem given in Definition 1.1. However, in general, there are other solutions to the two membranes problem in the sense of Definition 1.1 that are not minimizers of E(u, v). Since solutions according to Definition 1.1 are, in general, not unique, we observe that the limit that we prove to exist for the sequences that we construct iterating the obstacle problem from above and below does not necessarily converge to the unique minimizer to E(u, v).
Fully nonlinear operators. On the other hand, when we deal with viscosity solutions to the obstacle problems we use the general framework described in [12] and we understand sub and supersolutions applying the operator to a smooth test function that touches the graph of the sub/supersolution from above or below at some point in the domain. Remark that here we can consider fully nonlinear elliptic operators that need not be in divergence form.
When we develop the viscosity theory of the iterations in this general viscosity setting it is more difficult to obtain estimates in order to pass to the limit in the sequences \(\{u_n\}_n\), \(\{v_n\}_n\) and hence we have to rely again on monotonicity (that is obtained using comparison arguments).
In the viscosity framework there is no gain in considering two particular operators. Therefore, we will state and prove the results for two general second order elliptic operators that satisfy the comparison principle. However, if one looks for a model problem one can think on two normalized \(p-\)Laplacians, that is, let
Here \(\Delta w = trace (D^2w)\) (the usual Laplacian), \(\Delta _\infty w = \langle D^2 w \frac{\nabla w}{|\nabla w|}, \frac{\nabla w}{|\nabla w|} \rangle \) (the normalized infinity Laplacian), \(\alpha _i, \beta _i\) are nonnegative coefficients and \(h_p\), \(h_q\) are continuous functions. This operator, the normalized \(p-\)Laplacian, appears naturally when one considers game theory to solve nonlinear PDEs, we refer to [19, 22, 23, 27] and the recent books [3, 18].
Let us finish the introduction with a short paragraph concerning the equivalence between weak and viscosity solutions for divergence form operators. For the Dirichlet problem \(\Delta u =0\) the notions of weak and viscosity solutions coincide (and in fact the Dirichlet problem has a unique classical solution), see [16] and [25]. Moreover, the equivalence between weak and viscosity solutions include quasi-linear equations, [17, 20], and some non-local equations, [2, 14].
The paper is organized in two sections; in Sect. 2 we deal with that variational setting of the problem and in Sect. 3 we analyze the problem using viscosity theory.
2 Variational solutions
In this section we will obtain two sequences of solutions to variational obstacle problems whose limits constitute a pair of functions that is a variational solution to the two membranes problem.
Let us start by introducing some notations and definitions. We follow [15]. For short, we write
for the operators and their associated energies. In this variational context, we now introduce the definition of weak super and subsolutions.
Definition 2.1
Given \(p\in (1,\infty )\) and \(h_{p}\in W^{-1,p}(\Omega )\) we say that \(u\in W^{1,p}(\Omega )\) is a weak supersolution (resp. subsolution) for \(\mathcal {L}_{p}\) if
for all non-negative \(w\in W^{1,p}_0(\Omega )\) and in this case we written \(\mathcal {L}_{p}u\ge 0\) (resp. \(\mathcal {L}_{p}u\le 0\)) weakly in \(\Omega \). We say that \(u\in W^{1,p}(\Omega )\) is a weak solution for \(\mathcal {L}_{p}\) if
for all \(w\in W^{1,p}_0(\Omega )\) and we write \(\mathcal {L}_{p}u= 0\) weakly in \(\Omega \).
Definition 2.2
Given \(p,\,q\in (1,\infty )\), let \(f\in W^{1,p}(\Omega )\), \(\varphi \in L^{1}(\Omega )\), \(h_{p}\in W^{-1,p}\) and
When \(\underline{\Lambda }^{\,p}_{f,\varphi } \ne \emptyset \), we say \(u\in \underline{\Lambda }^{\,p}_{f,\varphi }\) is the solution of the \(p\,\)-Laplacian upper obstacle problem in a variational sense with obstacle \(\varphi \) and boundary datum f if
i.e, u minimizes the energy associated with the operator \(\mathcal {L}_{p}:=-\Delta _{p}+h_{p}\) in the set of functions that are above the obstacle and take the boundary datum, \(\underline{\Lambda }^{\,p}_{f,\varphi }\). In this case, we denote \(u=\underline{\text {O}}(\mathcal {L}_{p},\varphi , f)\). In Proposition 1 we show that this solution exists and is unique provided the set \(\underline{\Lambda }^{\,p}_{f,\varphi }\) is not empty.
Analogously, for \(g\in W^{1,q}(\Omega )\), \(\varphi \in L^{1}(\Omega )\), \(h_{q}\in W^{-1,q}\), we define
When \(\overline{\Lambda }^{\,q}_{g,\varphi } \ne \emptyset \), we say that \(v\in \overline{\Lambda }^{\,q}_{g,\varphi }\) is the solution of the \(q\,\)-Laplacian lower obstacle problem in a varational sense with obstacle \(\varphi \) and boundary g data if
For short, we denote v as \(v=\overline{\text {O}}(\mathcal {L}_{q},\varphi ,g)\).
Proposition 1
Let be \(h_{p}\in W^{-1,p}(\Omega )\) and \(f\in W^{1,p}(\Omega )\). If \( \underline{\Lambda }^{\,p}_{f,\varphi }\ne \emptyset \) there exists a unique \(u\in \underline{\Lambda }^{\,p}_{f,\varphi }\) such that \(E_{p}(u)\le E_{p}(w)\) for all \(w\in \underline{\Lambda }^{\,p}_{f,\varphi }\), i.e, there exists a unique \(u\in \underline{\Lambda }^{\,p}_{f,\varphi }\) that minimizes the energy \(E_p\) in \(\underline{\Lambda }^{\,p}_{f,\varphi }\). That is, \(u=\underline{\text {O}}(\mathcal {L}_{p},\varphi , f)\). This function satisfies
Analogously, let \(h_{q}\in W^{-1,q}(\Omega )\) and \(g\in W^{1,q}(\Omega )\), if \(\overline{\Lambda }^{\,q}_{g,\varphi }\ne \emptyset \), there exists a unique \(v\in \overline{\Lambda }^{\,q}_{g,\varphi }\) such that \(E_{q}(v)\le E_{q}(w)\) for all \(w\in \overline{\Lambda }^{\,q}_{g,\varphi }\). That is, \(v=\overline{\text {O}}(\mathcal {L}_{q},\varphi ,g)\). This function satisfies
Proof
Let us prove existence and uniqueness of the solution to the obstacle problem from below. The proof is contained in Theorem 3.21 in [15] but we reproduce some details here for completeness. Suppose that \(\underline{\Lambda }^{\,p}_{f,\varphi }\ne \emptyset \). Since \(E_p\) is a coercive and strictly convex functional in \(W^{1,p}(\Omega )\), we obtain that all minimizing sequences \(u_j\in \underline{\Lambda }^{\,p}_{f,\varphi }\) converge weakly to the same limit \(u\in W^{1,p}(\Omega )\) (that we want to show that is the unique minimizer of \(E_p\) in \( \underline{\Lambda }^{\,p}_{f,\varphi }\)). Let us show that \(u\in \underline{\Lambda }^{\,p}_{f,\varphi }\). Rellich–Kondrachov’s Theorem implies that, taking a subsequence, \(u_j\rightarrow u\) strongly in \(L^p(\Omega )\). Then, taking a subsequence, we have that \(u_j\rightarrow u\) a.e. in \(\Omega \). Thus, since \(u_j\ge \varphi \) a.e. \(\Omega \) and \(u_j-f\in W^{1,p}_0(\Omega )\) for all \(j\ge 1\), we get \(u\ge \varphi \) a.e. \(\Omega \) and \(u-f\in W^{1,p}_0(\Omega )\).
Let us prove that u verifies the properties stated in the proposition. Let us start with the fact that \(\mathcal {L}_p u\ge 0\) weakly in \(\Omega \). Since \(\underline{\Lambda }^{\,p}_{f,\varphi }\) is a convex set, \(u+t(v-u)\in \underline{\Lambda }^{\,p}_{f,\varphi }\) for all \( v\in \underline{\Lambda }^{\,p}_{f,\varphi }\) and \(t\in [0,1]\). Thanks to the fact that \(E_{p}\) reaches its minimum in the set \(\underline{\Lambda }^{\,p}_{f,\varphi }\) at u, we have that
satisfies \(i(t)\ge i(0)\) and therefore \(i'(0)\ge 0\). Let us compute
Taking \(t=0\), we get
If we consider \(w\in W^{1,p}_0(\Omega )\) \(w\ge 0\) a.e. \(\Omega \), and \(v=u+tw\in \underline{\Lambda }^{\,p}_{f,\varphi }\) for \(t\ge 0\) small enough, if we come back to (2.4), we obtain
with \(t\ge 0\). This implies that
weakly in \(\Omega \).
Finally, we show that
weakly in \(\{u(x) >\varphi (x)\}\) when \(\varphi \) is continuous. From Section 3.26 in [15] we have that the solution to the obstacle problem with continuous obstacle becomes also continuous (after a redefinition in a set of measure zero). Therefore, for a continuous obstacle the set \(\{x:u (x) > \varphi (x) \}\) is an open set. Take a ball
and consider a nonnegative \(w\in W^{1,p}_0(B_r) \cap C_0(\overline{B_r})\). Using that \(\varphi < u\) in \(\partial B_r\), we get \(v=u+tw\in \underline{\Lambda }^{\,p}_{f,\varphi }\) for |t| small enough. Then, let us come back to (2.4) and obtain
In this case t can be positive or negative. Thus, we get
which implies
weakly in \(B_r\) and therefore in \(\{u>\varphi \}\).
The other case (an obstacle from above) is analogous. \(\square \)
This minimizer of the energy is in fact the infimum of weak supersolutions are above/below the obstacle.
Proposition 2
Let \(u=\underline{\text {O}}(\mathcal {L}_{p},\varphi , f)\). This function satisfies
Analogously, let \(v=\overline{\text {O}}(\mathcal {L}_{q},\varphi ,g)\), then
Proof
Let us consider
We will prove that \(\overline{u}=u\) in \(\Omega \).
Let us start proving that \(u\ge \overline{u}\). By Proposition 1, we have
and then we get
Now, let us prove that \(u\le \overline{u}\). Given \(w\in \underline{\Lambda }^{\,p}_{f,\varphi }\) such that \(\mathcal {L}_{p}w\ge 0\) weakly in \(\Omega \), we have
and then, using (2.4), we get
Let us consider \(z=\min \{w,u\}\in \underline{\Lambda }^{\,p}_{f,\varphi }\), then \(u-z\in W^{1,p}_0(\Omega )\), \(u-z\ge 0\). Then
and, using (2.7), we obtain
Substracting these inequalities we conclude that
which implies
using that \((|b|^{p-2}b-|a|^{p-2}a)(a-b)\le 0\) (with a strict inequality for \(a\ne b\)) we obtain \(|\{u>z\}|=0\). Then, \(u=z\le w\). Thus, we have obtained that \(u\le \overline{u}\).
The other case is analogous. \(\square \)
Remark 1
Notice that Proposition 2 implies that the solution to the obstacle problem is monotone with respect to the obstacle in the sense that, for \(\varphi _1\ge \varphi _2\), we have \(u_1=\underline{\text {O}}(\mathcal {L}_{p},\varphi _1, f) \ge u_2=\underline{\text {O}}(\mathcal {L}_{p},\varphi _2, f)\).
Now, we are ready to introduce the definition of a weak solution to the two membranes problem.
Definition 2.3
Let \(f\in W^{1,p}(\Omega )\) and \(g\in W^{1,q}(\Omega )\) be two functions such that \(f\ge g\) in \(\partial \Omega \) in the sense of traces. Take \(h_p\in W^{-1,p}(\Omega )\) and \(h_{q}\in W^{-1,q}(\Omega )\). We say that the pair (u, v) with \(u\in \underline{\Lambda }^{\,p}_{f,v} \), \(v\in \overline{\Lambda }^{\,q}_{g,u}\) is a solution of the two membranes problem if
i.e.,
Remark 2
In general, given a pair (f, g), the solution to the two membranes problem is not unique.
2.1 Iterative method
The main result of this section reads as follows.
Theorem 2.4
For \(p\ge q \) let \(f\in W^{1,p}(\Omega )\) and \(g\in W^{1,q}(\Omega )\) be two functions such that \(f\ge g\) in \(\partial \Omega \) in the sense of traces and take \(h_p\in W^{-1,p}(\Omega )\) and \(h_{q}\in W^{-1,q}(\Omega )\). Let us consider \(v_{0}\) a weak subsolution of \( \mathcal {L}_{q}v=0\) in \(\Omega \) such that \(v_{0}-g\in W^{1,q}_{0}(\Omega )\) and then define inductively the sequences
Both sequences of functions \(\{u_{n}\}_{n=0}^{\infty }\subset W^{1,p}(\Omega )\), \(\{v_{n}\}_{n=0}^{\infty }\subset W^{1,q}(\Omega )\) converge strongly in \(W^{1,p}(\Omega )\) and \(W^{1,q}(\Omega )\), respectively. Moreover, the limits of the sequences are a solution of the two membranes problem. That is, there exists a pair of functions \(u_{\infty }\in \underline{\Lambda }^{\,p}_{f,v_{\infty }}\) and \(v_{\infty }\in \overline{\Lambda }^{\,q}_{g,u_{\infty }}\) such that
and, in addition, the limit pair \((u_{\infty },v_{\infty })\) satisfies
Proof
Notice that from the definition of the sequences as solutions to the obstacle problem we have
a.e. \(\Omega \) for each \(n\in \mathbb {N}\). Since we assumed that the boundary data are ordered, \(f\ge g\), and functions in \(W^{1,p}(\Omega )\) and in \(W^{1,q}(\Omega )\) have a trace on the boundary, we can also say that \(u_n\ge v_n\) a.e. on the boundary (with respect to the \((N-1)\)-dimensional measure). From the construction of the sequences we obtain that the sets \(\underline{\Lambda }^{\,p}_{f,v_n}\) and \( \overline{\Lambda }^{\,q}_{g,u_n}\) are not empty, in fact we have that \(u_{n-1} \in \underline{\Lambda }^{\,p}_{f,v_n}\) (since \(v_n\) solves an obstacle problem from above with obstacle \(u_{n-1}\), we have \(v_n \le u_{n-1}\)) and \(v_n \in \overline{\Lambda }^{\,q}_{g,u_n}\) (since \(u_n\) solves an obstacle problem from below with obstacle \(v_{n}\), we have \(v_n\le u_n\)). Hence the sequences \(\{u_{n}\}\) and \(\{v_{n}\}\) are well defined.
We will divide our arguments in several steps.
First step. First, we prove a monotonicity result. The two sequences are non-decreasing.
Let us see that \(v_{0}\le v_{1}\) a.e. \(\Omega \). This is due to the fact that \(v_{0}\) is a weak subsolution of \(\mathcal {L}_{q}v=0\) and \(v_{0}\in \overline{\Lambda }^{q}_{g,u_{0}}\) and by (2.6) we have
Then, we conclude that \(v_0\le v_1\).
Also we can see that \(u_{0}\le u_{1}\) a.e. \(\Omega \). In fact, from (2.5), we have
and then, using that \(v_{1}\ge v_{0}\) a.e. \(\Omega \), we obtain that \(u_{1}\in \underline{\Lambda }^{q}_{f,v_{0}}\) and \(u_{1}\) is a weak supersolution for \(\mathcal {L}_{p}u=0\). Then, \(u_0\le u_1\) a.e. in \(\Omega \).
Now, to obtain the general case, we just use an inductive argument. Suppose \(v_{n-1}\le v_{n}\) and \(u_{n-1}\le u_{n}\) a.e. \(\Omega \). Since \(v_{n}\le u_{n}\), we have \(v_{n}\in \overline{\Lambda }^{q}_{g,u_{n}}\). From (2.6) we have
which implies \(v_{n}\le v_{n+1}\) a.e. \(\Omega \). On the other hand, the last inequality and the fact that \(v_{n+1}\le u_{n+1}\) give us \(u_{n+1}\in \underline{\Lambda }^{p}_{f,v_{n}}\). Hence, (2.5) implies that
and we conclude that \(u_{n}\le u_{n+1}\) a.e. \(\Omega \).
Second step. Our next step is to show that the sequences \(\{u_{n}\}_{n=1}^{\infty }\) and \(\{v_{n}\}_{n=1}^{\infty }\) are bounded in \(W^{1,p}(\Omega )\) and \(W^{1,q}(\Omega )\), respectively. To prove this fact let us take \(w\in W^{1,q}(\Omega )\) the solution to \(\mathcal {L}_{q}w=0\) with boundary datum g. That is, w satisfies
By (2.6), we have that \(\mathcal {L}_{q}v_{n}\le 0\) weakly in \(\Omega \), and thanks to \(v_{n}-w\in W^{1,q}_{0}(\Omega )\) and the comparison principle for \(\mathcal {L}_{q}\), we obtain that
a.e. \(\Omega \). Let us consider \(z=\underline{\text {O}}(\mathcal {L}_{p}, w,f)\). Since \(v_{n}\le w\) we have \(z\in \underline{\Lambda }^{p}_{f,v_{n}}\) for all \(n\in \mathbb {N}\). Now, since \(u_n\) is the solution to the obstacle problem, \(u_{n}\) minimizes the energy among functions in the set \(\underline{\Lambda }^{p}_{f,v_{n}}\). Then, as z is in the former set, we get
for all \(n\in \mathbb {N}\), that is, we have
We can rewrite this inequality to obtain
Using that \(u_{n}-f\in W^{1,p}_{0}(\Omega )\), from Poincare’s inequality, we obtain
Then, by the triangle inequality we get
If we come back to (2.8), we obtain
Now, if we use Young’s inequality \(ab\le \varepsilon a^p+C(\varepsilon )b^{p'}\) with \(\varepsilon =\frac{1}{2}\) in the last term, we conclude
Thus, we get the desired bound for \(u_n\), there exists a constant C independent of n such that
Now, we prove that \(\{v_{n}\}\) is bounded in \(W^{1,q}(\Omega )\). Since \(v_{0}\le v_{n}\le u_{n-1}\) we have that \(v_{0}\in \underline{\Lambda }^{q}_{g,u_{n-1}}\) for all \(n\in \mathbb {N}\). This implies
That is,
Now, we just proceed as before to obtain
Third step. We show strong convergence to a solution of the two membranes problem.
Since \(\{u_{n}\}\) and \(\{v_{n}\}\) are bounded in \(W^{1,p}(\Omega )\) and \(W^{1,q}(\Omega )\) respectively, there exist subsequences \(\{u_{n_{j}}\}\subset \{u_{n}\}\) and \(\{v_{n_{j}}\}\subset \{v_{n}\}\) such that
as \( n_{j}\rightarrow \infty \).
By the Rellich–Kondrachov theorem, there exists a subsequence \(\{u_{n_j}\}\) that converges strongly in \(L^p(\Omega )\) and a subsequence \(\{v_{n_j}\}\) that converges strongly in \(L^q(\Omega )\). Extracting again a subsequence if needed, we get
as \( n_{j}\rightarrow \infty \). Using that \(\{u_n\}\) and \(\{v_n\}\) are increasing we obtain that the entire sequences converge pointwise (and weakly) to the unique limits (that are given by the supremums of the sequences). This implies that \(u_{\infty }\ge v_{\infty }\) a.e. \(\Omega \), \(u_{\infty }=f\) and \(v_{\infty }=g\) a.e. \(\partial \Omega \) in the sense of traces. This gives us \(u_{\infty }\in \underline{\Lambda }^{p}_{f,v_{\infty }}\) and \(v_{\infty }\in \overline{\Lambda }^{q}_{g,u_{\infty }}\).
Let us define \(u=\underline{\text {O}}(\mathcal {L}_{p},v_{\infty },f)\) and \(v=\overline{\text {O}}(\mathcal {L}_{q},u_{\infty },g)\). Our next step is to prove that \(u=u_{\infty }\) and \(v=v_{\infty }\) a.e. \(\Omega \) and to get that the convergence is strong in the corresponding Sobolev spaces.
First, we begin with the case of the upper membranes \(u_{n}\). Since \(u=\underline{\text {O}}(\mathcal {L}_{p},v_{\infty },f)\) and \(u_{\infty }\in \underline{\Lambda }^{p}_{f,v_{\infty }}\) we have that \(E_{p}(u)\le E_{p}(u_{\infty })\). On the other hand, since \(\{v_{n}\}\) is an increasing sequence we have \(u_{\infty }\ge v_{\infty }\ge v_{n}\) for each \(n\in \mathbb {N}\) wich implies \(u\in \underline{\Lambda }^{p}_{f,v_{n}}\). Then, we have that \(E_{p}(u_{n})\le E_{p}(u)\).
On the other hand, by the semicontinuity of the norm, \(u_{n}=u_{\infty }\) in \(\partial \Omega \) and by weak convergence, we obtain
Then, we have
Thus, \(E_{p}(u_{\infty })=E_{p}(u)\). And since the energy minimizer is unique, \(u_{\infty }=u\) a.e. \(\Omega \). Moreover, by (2.9) and the weak convergence, we have obtained that \(\{u_{n}\}\) is a minimizing sequence of the energy and the gradient of \( u_{n}\) converges strongly to the gradient of \(u_{\infty }\) as n goes to infinity. Then, we conclude that
Once we have seen that \(u_{\infty }=\underline{\text {O}}(\mathcal {L}_{p},v_{\infty },f)\), let us prove that \(v_{\infty }=\overline{\text {O}}(\mathcal {L}_{q},u_{\infty },g)\). Again, since \(v_{\infty }\in \overline{\Lambda }^{q}_{g,u_{\infty }}\), we have that \(E_{p}(v)\le E_{p}(v_{\infty })\).
Note that we cannot repeat the previous argument step by step because we do not know if v is in \(\overline{\Lambda }^{q}_{g,u_{n-1}}\). Though, it is enough to change a little our strategy (we use at this point that \(p\ge q\)). We construct a subsequence \(\{\tilde{v}_{n}\}\) such that \(\tilde{v}_{n}\in \overline{\Lambda }^{q}_{g,u_{n-1}}\) and \(E_{p}(\tilde{v}_{n})\) converges to \(E_{p}(v)\) as n goes to infinity. We define \(\{\tilde{v}_{n}\}\) as
Since \(p\ge q\) we have that \(\tilde{v}_{n}\in W^{1,q}(\Omega )\), \(\tilde{v}_{n}-g\in W^{1,q}_0(\Omega )\) and \(E_{q}(\tilde{v}_{n})\) goes to \(E_{q}(v)\) as \(n\rightarrow \infty \) due to the Rellick–Kondrachov Compactness Theorem and the fact that \(u_{n}\) converges to \(u_{\infty }\) in \(W^{1,p}(\Omega )\). Also we have the inequality \(\tilde{v}_{n}\le u_{n}\) because \(v\le u_{\infty }\). Then, \(\tilde{v}_{n}\in \overline{\Lambda }^{q}_{g,u_{n}}\) which implies \(E_{q}( v_{n+1})\le E_{q}(\tilde{v}_{n})\). Besides, by the semicontinuity of the norm, the fact that \(v_{n}=v_{\infty }\) on \(\partial \Omega \) and the weak convergence, we obtain
As a consequence,
Then, \(E_{q}(v_{\infty })=E_{q}(v)\). By the uniqueness of the minimizer, \(v_{\infty }=v\) a.e. \(\Omega \). Furthermore, from the same reasons as in the case of the upper membranes, we have that the entire sequence converges strongly,
This ends the proof. \(\square \)
If \(p<q\) we can also construct a pair of sequences that converges in \(W^{1,p}(\Omega )\) and \(W^{1,q}(\Omega )\) to a solution of the two membranes problem. In this case we just have to start the iterations of the obstacle problems with \(u_{0}\) a weak supersolution to \( \mathcal {L}_{p}u=0\).
Theorem 2.5
For \(p<q\) let \(f\in W^{1,p}(\Omega )\) and \(g\in W^{1,q}(\Omega )\) be two functions such that \(f\ge g\) in \(\partial \Omega \) in the sense of traces and take \(h_p\in W^{-1,p}(\Omega )\) and \(h_{q}\in W^{-1,q}(\Omega )\). Take \(u_{0}\) a weak supersolution to \( \mathcal {L}_{p}u=0\) in \(\Omega \) such that \(u_{0}-f\in W^{1,p}_{0}(\Omega )\) and then let
Both sequences \(\{u_{n}\}_{n=0}^{\infty }\subset W^{1,p}(\Omega )\) and \(\{v_{n}\}_{n=0}^{\infty }\subset W^{1,q}(\Omega )\) converge strongly in \(W^{1,p}(\Omega )\) and \(W^{1,q}(\Omega )\), respectively. Moreover, the limits of the sequences are a solution of the two membranes problem. That is, there exists a pair of functions \(u_{\infty }\in \underline{\Lambda }^{\,p}_{f,v_{\infty }}\) and \(v_{\infty }\in \overline{\Lambda }^{\,q}_{g,u_{\infty }}\) such that
and, in addition, the limit pair \((u_{\infty },v_{\infty })\) satisfies
Proof
The difference between these sequences and the sequences defined in Theorem 2.4 is that, here, \(\{u_{n}\}\) and \(\{v_{n}\}\) are decreasing sequences.
Due to the monotonicity of the sequences, we can prove that \(u_{n}\) converges weakly to some \(u_{\infty }\) in \(W^{1,p}(\Omega )\) and \(v_{n}\) converges to some \(v_{\infty }\) weakly in \(W^{1,q}(\Omega )\) as n goes to infinity.
The strong convergence of \(v_{n}\) to \(v_{\infty }\) is given in the same way that we got strong convergence for \(u_{n}\) to \(u_{\infty }\) in Proposition 2.4 thanks to \(\{v_{n}\}\) is a decreasing sequence. On the other hand, for \(\{u_{n}\}\) we can reproduce the proof for \(\{v_{n}\}\) in proposition (2.4) because \(W^{1,q}(\Omega )\hookrightarrow W^{1,p}(\Omega )\) continuously. \(\square \)
Remark 3
An alternative idea to deal with the case \(p<q\) runs as follows: It can be easily proved that
Then, if \(p<q\) we consider
For the problem with the new set of parameters, \(p'\), \(q'\), \(h_{p'}\), \(h_{q'}\), \(f'\) and \(g'\), we can apply the iterative method described in Theorem 2.4 and get a pair \((u_{\infty }',v_{\infty }')\) such that \(u_{\infty }\in \underline{\Lambda }^{p'}_{f',v_{\infty }'}\) and \(v_{\infty }'\in \overline{\Lambda }^{q'}_{g',u_{\infty }'}\) such that \(u_{\infty }'=\underline{\text {O}}(-\Delta _{p'}+h_{p'},v_{\infty }',f')\) and \(v_{\infty }'=\overline{\text {O}}(-\Delta _{q'}+h_{q'},u_{\infty }',g')\). Thus, \(u_{\infty }=-v_{\infty }'\) and \(v_{\infty }=-u_{\infty }'\) are in \(\underline{\Lambda }^{p}_{f,v_{\infty }}\) and \(\overline{\Lambda }^{q}_{g,u_{\infty }}\) respectively and
i.e, \((u_{\infty },v_{\infty })\) is a solution for the original two membranes problem.
Remark 4
We remark that the limit depends strongly on the initial function from where we start the iterations. We would like to highlight that \(\{u_{n}\}\) and \(\{v_{n}\}\) converge to \(u_{\infty }\) and \(v_{\infty }\) strongly in \(W^{1,p}(\Omega )\) and \(W^{1,q}(\Omega )\) respectively due to the fact that both sequences are monotone. The monotonicity, in turn, depends on the initial function for the iteration. Specifically, when \(p\ge q\), the monotonicity arises because the initial datum \(v_{0}\) is a weak subsolution of \(\mathcal {L}_{q}v=0\). On the other hand, when \(p\le q\) it stems from the fact that \(u_{0}\) is a weak supersolution of \(\mathcal {L}_{p}u=0\).
If one considers the pair (u, v) obtained as follows: let u be the solution to \( \mathcal {L}_{p}u=0\) in \(\Omega \) with \(u=f\) on \(\partial \Omega \) and v the solution to the corresponding obstacle problem \(v=\overline{\text {O}}( \mathcal {L}_{q},u,g)\), then, as we have mentioned in the introduction this pair (u, v) is a solution to the two membranes problem. Now, if one starts the iteration procedure with \(u_0 =u\) we obtain that the sequences converge after only one step. In fact we get \(v_1 =\overline{\text {O}}( \mathcal {L}_{q},u,g)\) and next the iteration gives again, that is, we have \(u =\underline{\text {O}}( \mathcal {L}_{p},v_1,f)\). This is due to the fact that u is a solution to \( \mathcal {L}_{p}u=0\) in the whole \(\Omega \) with \(u\ge v_1\) and therefore it is the solution to the obstacle problem (with \(v_1\) as obstacle from below). Analogously, if one starts with v the solution to \( \mathcal {L}_{q}v = 0\) in \(\Omega \) with \(v=g\) on \(\partial \Omega \) we obtain u as the solution to the corresponding obstacle problem \(u=\underline{\text {O}}( \mathcal {L}_{p},v,f)\). This pair (u, v) is also a solution to the two membranes problem, but, in general the two pairs are different.
2.2 Extension to the two membranes problem for nonlocal operators
If we consider
with
the fractional \(p-\)Laplacian. Notice that \(\mathcal {L}_{p}\) and \(\mathcal {L}_{q}\) are associated with the energies
and
Our iterative method of the Theorem 2.4 gives us a pair \((u_{\infty },v_{\infty })\) that is solution of the “new” two membranes problem whenever
To be able to define a variational solution of the obstacle problem for these new \(\mathcal {L}_{p}\) and \(\mathcal {L}_{q}\) is needed that \(\underline{\Lambda }^{p}_{f,\varphi }\) and \(\overline{\Lambda }^{q}_{g,\varphi }\) are closed sets under the weak topology. This is given by Mazur’s Theorem if \(W^{s,p}(\Omega )\) and \(W^{t,q}(\Omega )\) are reflexive Banach space and \(\underline{\Lambda }^{p}_{f,\varphi }\) and \(\overline{\Lambda }^{q}_{g,\varphi }\) are closed convex sets (with the norm) in \(W^{s,p}(\Omega )\) and \(W^{t,q}(\Omega )\) respectively.
3 Viscosity solutions
In this section we aim to prove similar convergence results for iterations of the obstacle problem to a pair that solves the two membranes problems when the involved operators are not variational. Here we will consider two nonlinear elliptic operators and we understand solutions in the viscosity sense. To define the sequences as solutions to the corresponding obstacle problems we need to assume that the involved operators verify the following set of conditions:
Hypothesis The operators \(\mathcal {L}_{1}\) and \(\mathcal {L}_{2}\) verify
-
\(\mathcal {L}w= F(D^{2}w,\nabla w)-h(x)\) with F continuous in both coordinates, such that the comparison principle holds.
-
If f is continuous, there exists a unique \(w\in C(\overline{\Omega })\) solution to \(\mathcal {L}w=0\) in \(\Omega \), \(w=f\) in \(\partial \Omega \).
-
If \(\varphi \), \(\psi \) and f, g are continuous, then the corresponding solutions to the obstacle problems, \(u=\underline{\mathcal {O}}(\mathcal {L}_{1},\varphi ,f)\) and \(v=\overline{\mathcal {O}}(\mathcal {L}_{2},\psi ,g)\) are continuous. See Definition 3.4 below.
To be precise, working in the viscosity sense, we need to introduce the definition of semicontinuous functions.
Definition 3.1
\(f:\Omega \longrightarrow \mathbb {R}\) is a lower semicontinuous function, l.s.c, at \(x\in \Omega \) if for each \(\varepsilon >0\) there exists a \(\delta >0\) such that
If a function f is l.s.c. at every \(x\in \Omega \), we say that f is l.s.c. in \(\Omega .\) The lower semicontinuous envelope of f is
On the other hand, \(g:\Omega \longrightarrow \mathbb {R}\) is a upper semicontinuous function, u.s.c, at \(x\in \Omega \) if for each \(\varepsilon >0\) there exists a \(\delta >0\) such that
If a functions g is u.s.c. at every \(x\in \Omega \), we say that it is u.s.c. in \(\Omega .\) The upper semicontinuous envelope of g is
Now, we give the precise definition of sub and supersolutions in the viscosity sense.
Definition 3.2
The function \(u:\Omega \longrightarrow \mathbb {R}\) is called a \(\mathcal {L}\)-viscosity supersolution if its lower semicontinuous envelope function, \(u_{*}\), satisfies the following: for every \(\phi \in C^{2}(\overline{\Omega })\) such that \(\phi \) touches \(u_{*}\) at \(x\in \Omega \) strictly from below, that is, \(u_{*}-\phi \) has a strict minimum at x with \(u_{*}(x)=\phi (x)\), we have
In this case, we write \(\mathcal {L}u\ge 0\) in the viscosity sense.
Conversely, u is called a \(\mathcal {L}\)-viscosity subsolution if its upper semicontinuous envelope function, \(u^{*}\), satisfies that for every \(\phi \in C^{2}(\overline{\Omega })\) such that \(\phi \) touches \(u^{*}\) at \(x\in \Omega \) strictly from above, that is, \(\phi -u^{*}\) has a strict minimum at x with \(u^{*}(x)=\phi (x)\), we have
and we write \(\mathcal {L}u\le 0\) in the viscosity sense.
Finally, u is a \(\mathcal {L}\)-viscosity solution if it is both a \(\mathcal {L}\)-viscosity supersolution and a \(\mathcal {L}\)-viscosity subsolution and we denote \(\mathcal {L}u=0\) in the viscosity sense.
Now, let us introduce three results concerning the limit of non-decreasing sequences of viscosity sub or supersolutions. Although these results are not difficult to prove, they will be crucial for the proof of our main Theorem 3.9 in this section.
Proposition 3
Let be \(\{v_{n}\}_{n=0}^{\infty }\subset C(\overline{\Omega })\) a non-decreasing sequence of functions, continuous up to the boundary such that \(v_{n}\) is a \(\mathcal {L}\)-viscosity subsolution for each \(n\in \mathbb {N}\) and \(v_{n}\longrightarrow v_{\infty }\) pointwise. Then, \(v_{\infty }\) is a \(\mathcal {L}\)-viscosity subsolution.
Proof
Let be \(\phi \in C^{2}(\overline{\Omega })\) such that \(\phi -v_{\infty }^{*}\) has a strict minimum at \(x_{0}\in \Omega \) and \(\phi (x_{0})=v_{\infty }^{*}(x_{0})\).
Since \(v_{\infty }^{*}\) is upper semicontinuous, \(\phi -v_{\infty }^{*}\) is lower semicontinuous. Then, fixed \(r>0\) small enough, \(\phi -v_{\infty }^{*}\) reaches a minimum in \(D_{r}:=\overline{\Omega }\setminus B_{r}(x_{0})\). Say in \(z_{\infty }\in D_{r}\). Thanks to \(- v_{\infty }^{*}=(-v_{\infty })_{*}\) and that the sequence \(\{v_{n}\}\) is non-decreasing, we have
for each \(z\in D_{r}\) and for each \(n\in \mathbb {N}\).
On the other hand, by the definition of pointwise value of the lower semicontinuous envelope,
Then, given \(\varepsilon >0\) there exists a sequence \(\{y_{k}\}\) within \(B_{r/2}(x_{0})\) such that
Thus, there is \(k_{0}\in \mathbb {N}\) such that \(|(\phi -v_{\infty })(y_{k_{0}})|\le 2\varepsilon /3\). Since \(v_{n}\) converges to \(v_{\infty }\) pointwise, there exists \(n_{0}\in \mathbb {N}\) such that \(|(\phi -v_{\infty })(y_{k_{0}})-(\phi -v_{n})(y_{k_{0}})|\le \varepsilon /3\) for all \(n\ge n_{0}\).
In short, for each \(\varepsilon >0\) there exists a \(y\in B_{r/2}(x_{0})\) and \(n_{0}\in \mathbb {N}\) such that
Since, \(\phi -v_{n}\) is a continuous function, given \(\varepsilon <(\phi -v_{\infty }^{*})(z_{\infty })\) one can see that the minimum of \(\phi -v_{n}\) in \(\overline{\Omega }\) is reached at some \(x_{n}^{r}\in B_{r/2}(x_{0})\) for all \(n\ge n_{0}\) thanks to the above inequality and (3.10). Note that \(n_{0}\) depends on r, so we write \(n_{0}=n_{0}(r)\). We have \(n_{0}(r)\rightarrow \infty \) and \(x_{n}^{r}\rightarrow x_{0}\) as \(r\rightarrow 0^{+}\). Due to the fact that \(v_{n}\) is a \(\mathcal {L}\)- viscosity subsolution, we have
Since \(\mathcal {L}\) is a continuous operator we have
Therefore \(v_{\infty }\) is a \(\mathcal {L}\)-viscosity subsolution. \(\square \)
Lemma 3.3
Let \(\{u_{n}\}\subset C(\overline{\Omega })\) a non-decreasing sequence of continuous up to the boundary functions. If \(u_{n}\) converge pointwise to some function \(u_{\infty }\), then \(u_{\infty }\) is a lower semicontinuous function in \(\overline{\Omega }\).
Proof
Let be \(x_{0}\in \Omega \). Since \(u_{n}\longrightarrow u_{\infty }\) pointwise, we have that for each \(\varepsilon >0\) there exists a \(n_{0}\in \mathbb {N}\) such that \(u_{\infty }(x_{0})\le u_{n_{0}}(x_{0})+\varepsilon .\) On the other hand, due to the fact that \(u_{n_{0}}\) is continuous there exist a \(\delta \) such that \(u_{n_{0}}(x_{0})\le u_{n_{0}}(y)+\varepsilon \) for all \(y\in B_{\delta }(x_{0})\). Moreover, thanks to the fact that the sequence \(\{u_{n}\}\) is non-decreasing we get that, for every \(\varepsilon >0\) there exists \(\delta >0\) such that
The proof is finished. \(\square \)
Proposition 4
Let \(\{u_{n}\}_{n=0}^{\infty }\subset C(\overline{\Omega })\) be a non-decreasing sequence of continuous up to the boundary functions such that \(u_{n}\) is a \(\mathcal {L}\)-viscosity supersolution for each \(n\in \mathbb {N}\) and \(u_{n}\longrightarrow u_{\infty }\) pointwise. Then, \(u_{\infty }\) is a \(\mathcal {L}\)-viscosity supersolution.
Proof
By Lemma 3.3 we know that \(u_{\infty }\) is a lower semicontinuous function. To prove that \(u_{\infty }\) is a \(\mathcal {L}\)-viscosity supersolution let us consider \(\phi \in C^{2}(\overline{\Omega })\) such that \(u_{\infty }-\phi \) has a strict minimum at \(x_{0}\in \Omega \) such that \((u_{\infty }-\phi )(x_{0})=0\). We want to obtain that \((\mathcal {L}\phi )(x_{0})\ge 0\).
Let \(r>0\) be a fixed small radius. Since \(u_{\infty }-\phi \) is lower semicontinuous, these functions reach a positive minimum in \(\overline{\Omega }\setminus B_{r}(x_{0})\), say in \(z_{\infty }\in \overline{\Omega }{\setminus } B_{r}(x_{0})\), then, we have
Also, since \(u_{n}-\phi \) is continuous, there exists \(z_{n}\in \overline{\Omega }\setminus B_{r}(x_{0})\) where \(u_{n}-\phi \) reaches a minimum in \(\overline{\Omega }\setminus B_{r}(x_{0})\). We claim that \((u_{n}-\phi )(z_{n})\) converges to \((u_{\infty }-\phi )(z_{\infty })\) as n to infinity. Let us continue the proof assuming this fact and prove it after finishing our argument.
The claim implies that there exists \(n_{0}=n_{0}(r)\) such that
for all \(x\in \overline{\Omega }{\setminus } B_{r}(x_{0})\) and \(n\ge n_{0}\).
On the other hand, since \(u_{n}\) converges to \(u_{\infty }\) pointwise and \((u_{\infty }-\phi )(x_{0})=0\), there exists \(n_{1}=n_{1}(r)\) such that
Thus, the above and (3.11) imply \(u_{n}-\phi \) reaches its minimum in \(\overline{\Omega }\) in the interior of the ball \(B_{r}(x_{0})\) for all \(n\ge n_{3}=\max \{n_{0},n_{1}\}\). Call this point \(x_{n}^{r}\). Then, since \(u_{n}\) is a \(\mathcal {L}\)-viscosity supersolution, we get
Due to the fact that \(x_{n}^{r}\rightarrow x_{0}\) as \(r\rightarrow 0^{+}\), we have
Therefore, \(u_{\infty }\) is a \(\mathcal {L}\)-viscosity supersolution in \(\Omega \).
Now, we will prove the claim. Let us define \(w_{\infty }=u_{\infty }-\phi \), \(w_{n}=u_{n}-\phi \) and \(a_{n}=w_{n}(z_{n})\). Thanks to \(\{u_{n}\}\) is a pointwise increasing sequence, \(\{a_{n}\}\) is also an increasing sequence,
Moreover, the sequence is bounded above by \(w_{\infty }(z_{\infty })\) due to the fact that \(u_{n}\) converges pointwise to \(u_{\infty }\). Hence,
Then, there exists a number \(a_{\infty }\le w_{\infty }(z_{\infty })\) such that \(a_{n}\) converges to \(a_{\infty }\) as n tends to infinity.
Suppose that \(a_{\infty }<w_{\infty }(z_{\infty })\). Then, because \(\{z_{n}\}\) is a collection of points in a compact set, there exists a subsequence \(\{z_{n_{j}}\}\subset \{z_{n}\}\) and a point \(\widetilde{z}\in \overline{\Omega }\setminus B_{r}(x_{0})\) such that
Since \(w_{n}\) is a lower semicontinuous function, for each \(\varepsilon >0\), there exists a \(\delta =\delta (\varepsilon ,\widetilde{z})\) such that
Due to the fact that \(z_{n_{j}}\) converges to \(\widetilde{z}\), we can suppose \(z_{n_{j}}\in B_{\delta }(\widetilde{z})\) for all \(n_{j}\). Then
For each \(n_{j}\) there exists a \(m_{j}\in \mathbb {N}\) such that
Observe that we can construct a sequence \(\{(n_{j},k_{j})\}\) where \(k_{j}\ge m_{j}\) for all \(j\in \mathbb {N}\) and \((n_{j},k_{j})\longrightarrow (\infty ,\infty )\) as \(j\rightarrow \infty \) such that
Since \(a_{\infty }<w_{\infty }(z_{\infty })\) there exists \(\varepsilon _{0}>0\) such that
For that \(\varepsilon _{0}\) the above argument implies that there exists a sequence \(\{(n_{j},k_{j})\}\) that goes to \((\infty ,\infty )\) as \(j\rightarrow \infty \) such that
and therefore,
On the other hand, since \(\{a_{n}\}\) is a increasing sequence and \(a_{k_{j}}=w_{k_{j}}(z_{k_{j}})\) we have that
Thanks to the fact that for each \(j\in \mathbb {N}\) the functions \(w_{k_{j}}\) are continuous, \(\overline{\Omega }\setminus B_{r}(x_{0})\) is connected and \(z_{n_{j}}\) and \(z_{k_{j}}\) are in \(\overline{\Omega }\setminus B_{r}(x_{0})\), there exists \(y_{j}\in \overline{\Omega }\setminus B_{r}(x_{0})\) such that
Thus, there is a subsequence \(\{y_{j_{i}}\}\subset \{y_{j}\}\) and a point \(y\in \overline{\Omega }\setminus B_{r}(x_{0})\) such that \(y_{j_{i}}\rightarrow y\) as \(i\rightarrow \infty \). Due to the fact that \(w_{\infty }\) is a lower semicontinuous function we have
However, this is a contradiction because \(y\in \overline{\Omega }\setminus B_{r}(x_{0})\) and that implies
Therefore, we have the desired equality \(a_{\infty }=w_{\infty }(z_{\infty })\). \(\square \)
Now, we are ready to introduce the definition of a solution to the obstacle problem in the viscosity sense.
Definition 3.4
Given \(\varphi :\overline{\Omega }\rightarrow \mathbb {R}\) and \(f\in C(\partial \Omega )\) such that \(f\ge \varphi ^{*}\) in \(\partial \Omega \). We say that \(u:\Omega \longrightarrow \mathbb {R}\) is a viscosity solution of the \(\mathcal {L}\)-lower obstacle problem with obstacle \(\varphi \) and boundary datum f and we denote \(u=\underline{\mathcal {O}}(\mathcal {L},\varphi ,f)\) if u satisfies:
Whereas, given \(\psi \) and g such that \(f\le \psi _{*}\) on \(\partial \Omega \), we say that v is a viscosity solution of the \(\mathcal {L}\)-upper obstacle problem with obstacle \(\psi \) and boundary datum g, and we denote \(v=\overline{\mathcal {O}}(\mathcal {L},\psi ,g)\), if v satisfies:
Remark 5
Note that, when u is a viscosity solution of the \(\mathcal {L}\)-lower obstacle problem with obstacle \(\varphi \) and boundary datum f, then, its lower and upper semicontinuous envelopes, \(u_{*}\) and \(u^{*}\), are also viscosity solutions for the same problem.
The same applies to the upper obstacle, if \(v=\overline{\mathcal {O}}(\mathcal {L},\psi ,g)\), then, \(v_{*}=\overline{\mathcal {O}}(\mathcal {L},\psi ,g)\) and \(v^{*}=\overline{\mathcal {O}}(\mathcal {L},\psi ,g)\).
However, we have uniqueness of solutions to the obstacle problem up to semicontinuous envelopes.
Lemma 3.5
Given \(\varphi \) and f defined as before
-
(a)
if \(f\ge \varphi ^{*}\) in \(\partial \Omega \), there exists at most one lower semicontinuous function u such that \(u=\underline{\mathcal {O}}(\mathcal {L},\varphi ,f)\).
-
(b)
if \(f\le \varphi ^{*}\) in \(\partial \Omega \), there exists at most one upper semicontinuous function u such that \(u=\overline{\mathcal {O}}(\mathcal {L},\varphi ,f)\).
Proof
We will only prove item (a), the other item is analogous. (a) Suppose that there exists two lower semicontinuous functions solutions \(u_1\) and \(u_2\), let us prove that \(u_2\ge u_1\).
In the set \(\{u_1=\varphi ^{*}\}\) we have \(u_2\ge \varphi ^{*}=u_1\). In the open set \(\{u_1>\varphi ^{*}\}\) we have that \(\mathcal {L}u_2\ge 0\) and \(\mathcal {L}u_1=0\) and \(u_2\ge u_1\) in \(\partial \{u_1>\varphi ^{*}\}\). Then, by the Comparison Principle we obtain that \(u_2\ge u_1\) in \(\{u_1>\varphi ^{*}\}\). Thus, we conclude that
The reverse inequality follows interchanging the roles of \(u_1\) and \(u_2\). \(\square \)
In view of the previous lemma, from now we will suppose that \(u=\underline{\mathcal {O}}(\mathcal {L},\varphi ,f)\) is lower semicontinuous and \(v=\overline{\mathcal {O}}(\mathcal {L},\varphi ,f)\) is upper semicontinuous.
Next, we show that when the obstacles and the boundary data are ordered then the solutions to the obstacle problems are also ordered.
Lemma 3.6
-
(a)
Given \(\varphi _1\) and \(\varphi _2\) functions such that \(\varphi _1\le \varphi _2\) (a.e.) and \(f_1, f_2\in C(\partial \Omega )\) such that \(f_1\le f_2\), if \(f_1\ge \varphi _1^{*}\) and \(f_2\ge \varphi _2^{*}\) on \(\partial \Omega \), let us consider \(u_1=\underline{\mathcal {O}}(\mathcal {L},\varphi _1,f)_1\) and \(u_2=\underline{\mathcal {O}}(\mathcal {L},\varphi _2,f_2)\), then
$$\begin{aligned} u_1\le u_2\;\;\text {in}\;\;\Omega . \end{aligned}$$ -
(b)
Given \(\psi _1\) and \(\psi _2\) functions such that \(\psi _1\le \psi _2\) (a.e.) and \(g_1, g_2\in C(\partial \Omega )\) such that \(g_1\le g_2\), if \(g_1\le \psi _{1{*}}\) and \(g_2\le \psi _{2{*}}\) on \(\partial \Omega \), let us consider \(v_1=\overline{\mathcal {O}}(\mathcal {L},\psi _1,g_1)\) and \(v_2=\overline{\mathcal {O}}(\mathcal {L},\psi _2,g_2)\), then
$$\begin{aligned} v_1\le v_2 \;\;\text {in}\;\;\Omega . \end{aligned}$$
Proof
We will prove item (a), the other case is analogous.
(a) We will consider two cases:
In the set \(\{u_1=\varphi _1^{*}\}\) we have \(u_2\ge \varphi _2^{*}\ge \varphi _1^{*}=u_1\).
In the open set \(\{u_1>\varphi _1^{*}\}\) we have that \(\mathcal {L}u_2\ge 0\) and \(\mathcal {L}u_1=0\) and \(u_2\ge u_1\) in \(\partial \{u_1>\varphi _1^{*}\}\). Then, by the Comparison Principle we obtain that \(u_2\ge u_1\) in \(\{u_1>\varphi ^{*}\}\).
Thus \(u_2\ge u_1\) in \(\Omega \). \(\square \)
Lemma 3.7
Given \(\varphi \) and f as before, let us consider \(u=\underline{\mathcal {O}}(\mathcal {L},\varphi ,f)\), then
In the same way, if we take \(v=\overline{\mathcal {O}}(\mathcal {L},\varphi ,f)\), then
Proof
Let us prove the first claim. Consider
Since \(u\ge \varphi ^{*}\) in \(\Omega \), \(u=f\) on \(\partial \Omega \) and \(\mathcal {L}u\ge 0\) in \(\Omega \) in viscosity sense, we get
Now, let us consider
In the set \(\{u=\varphi ^{*}\}\) we have \(w\ge \varphi ^{*}=u\).
In the open set \(\{u>\varphi ^{*}\}\) we have that \(\mathcal {L}w\ge 0\) and \(\mathcal {L}u=0\) in the viscosity sense and \(w\ge u\) in \(\partial \{u>\varphi ^{*}\}\). Then, by the Comparison Principle we obtain that \(w\ge u\) in \(\{u_1>\varphi ^{*}\}\).
Then, \(u\le w\) in \(\Omega \). Taking minimum we get
The other case is analogous. \(\square \)
Now, we introduce the definition of a solution to the two membranes problem in the viscosity sense (this is just Definition 1.1 with solutions understood in the viscosity sense).
Definition 3.8
Given \(f,g\in C(\partial \Omega )\) with \(f\ge g\), let \(\mathcal {L}_{1}\) and \(\mathcal {L}_{2}\) by two operators that satisfy the hypothesis at the begining of the section. We say a pair of functions (u, v) is a solution of the two membranes problem with boundary data (f, g) if
Our main result is the following.
Theorem 3.9
Given f, g and \(\mathcal {L}_{1},\mathcal {L}_{2}\), let us consider \(v_{0}\in C(\overline{\Omega })\) a \(\mathcal {L}_{2}\)-viscosity subsolution such that \(v_{0}\le g\) in \(\partial \Omega \) and then define inductively the sequences
Both sequences of functions \(\{u_{n}\}_{n=0}^{\infty }\subset C(\overline{\Omega })\), \(\{v_{n}\}_{n=0}^{\infty }\subset C(\overline{\Omega })\) converge to some limit functions \(u_{\infty }\) and \(v_{\infty }\), respectively. Moreover, theses limits are a solution of the two membranes problem with boundary data (f, g), that is,
In addition, the functions \(u_{\infty }\) and \(v_{\infty }\) are lower semicontinuous in \(\bar{\Omega }\) and continuous in the interior of the set where \(u_{\infty }\) touches \(v_{\infty }\), i.e, in the interior of \(\{ u_{\infty }=v_{\infty }\}\).
Proof
As before, we divide the proof in several steps.
First step. Let us start proving that \(\{u_{n}\}_{n=0}^{\infty }\) and \(\{v_{n}\}_{n=0}^{\infty }\) are non-decreasing.
Let us see that \(v_0\le v_1\). Recall that \(v_1=\overline{\mathcal {O}}(\mathcal {L}_{2},u_0,g)\) and \(u_0=\underline{\mathcal {O}}(\mathcal {L}_{1},v_0,f)\), then \(v_1\le u_0\) and \(v_0\le u_0\).
In the set \(\{v_1=u_0\}\) we have \(v_0\le u_0=v_1\).
In the open set \(\{v_1<u_0\}\) we have that \(\mathcal {L}_2 v_1= 0\) and \(\mathcal {L}_2 v_0 \le 0\) and \(v_1\ge v_0\) in \(\partial \{v_1>u_0\}\). Then, by the Comparison Principle we obtain that \(v_1\ge v_0\) in \(\{v_1>u_0\}\).
Thus, \(v_0\le v_1\) in \(\Omega \).
Now, using Lemma 3.6 we have \(u_0\le u_1\). By induction we can continue and obtain \(u_n\le u_{n+1}\) and \(v_n\le v_{n+1}\) for all \(n\ge 1\).
Second step. Let us prove that \(\{u_{n}\}_{n=0}^{\infty }\) and \(\{v_{n}\}_{n=0}^{\infty }\) are bounded.
Let \(w\in C(\overline{\Omega })\) be a \( \mathcal {L}_2\)-solution with \(w=g\) in \(\partial \Omega \). Then, since \(\mathcal {L}_2 w= 0\), we get \(\mathcal {L}_2 v_n \le 0\). Moreover, we have \(w\ge v_n\) in \(\partial \Omega \). Hence, applying the Comparison Principle, we get \(w\ge v_{n}\) in \(\bar{\Omega }\) for all \(n\ge 0\). This implies, together with the fact that \(v_{n}\ge v_{0}\), that \(\{v_n\}\) is bounded.
Now, if we consider \(z=\underline{\mathcal {O}}(\mathcal {L}_{1},w,f )\), using that \(v_n\le w\) by Lemma 3.6\(u_n\le z\) and \(z\in C(\overline{\Omega })\). This implies that \(\{u_n\}\) is bounded.
Third step. Let us define
for all \(x\in \overline{\Omega }\). Using that \(\{u_n\}\) are continuous and the sequence is increasing we obtain that \(u_{\infty }\) is lower semicontinuous. We also have \(u_{\infty }=f\) and \(v_{\infty }=g\) in \(\partial \Omega \). Let us consider
This function v is well defined and upper semicontinuous. Our goal is to prove that \(v=v_{\infty }\).
Using that \(u_n\le u_{\infty }\), we get \(v_n\le v\) for all \(n\ge 0\) by Lemma 3.6. Then, we have that \(v_{\infty }\le v\) and \(v_{\infty }^*\le v^*=v\). This implies that \(v_{\infty }^* \le u_{\infty }\). On the other hand, by Proposition 3, we have that \(\mathcal {L}_2 v_{\infty }^*\le 0\) in \(\Omega \). We also have \(v_{\infty }^*\ge g\) in \(\partial \Omega \). Now, we consider \(w\in C(\Omega )\) a solution to
By the Comparison Principle we have that \(v_n\le w\) in \(\overline{\Omega }\). Then \(v_{\infty }\le w\) in \(\overline{\Omega }\) and \(v_{\infty }^*\le w^*=w\) in \(\overline{\Omega }\). This implies \(v^{*}_{\infty }\le g\) on \(\partial \Omega \) and, therefore, \(v^{*}_{\infty }= g\) on \(\partial \Omega \). Finally, we will prove that \(\mathcal {L}_2v_{\infty }^*=0\) in the (open) set \(\{v_{\infty }^*<u_{\infty }\}\). Let us consider \(x_0\in \{v_{\infty }^*<u_{\infty }\}\), and call \(\delta =u_{\infty }(x_0)-v_{\infty }^*(x_0)\). Let us consider \(n_0\ge 0\) such that \(u_{\infty }(x_0)-u_{n_0}(x_0)<\frac{\delta }{4}\). Using that \(u_{n_0}\) is continuous, there exists \(\eta _1=\eta _1(n_0)>0\) such that \(|u_{n_0}(x_0)-u_{n_0}(y)|<\frac{\delta }{4}\) for all \(y\in B_{\eta _1}(x_0)\). On the other hand, there exists \(\eta _2>0\) such that \(v_{\infty }^*(y)<v_{\infty }^*(x_0)+\frac{\delta }{4}\) for all \(y\in B_{\eta _2}(x_0)\). Gathering the previous estimates we obtain
for all \(y\in B_{\eta }(x_0)\) with \(\eta =\min \{\eta _1,\eta _2\}\). Now, we have that \(u_n(y)\ge u_{n_0}(y)\) and \(v_n(y)\le v_{\infty }^*(y)\) for all \(n\ge n_0\). Thus
for all \(n\ge n_0\) and for all \(y\in B_{\eta }(x_0)\). Then \(B_{\eta }(x_0)\subseteq \{u_n>v_n\}\) for all \(n\ge n_0\). Using that \(\mathcal {L}_2 v_n =0\) in \(B_{\eta }(x_0)\) and taking the limit we obtain \(\mathcal {L}_2 v_{\infty }=0\) in \(B_{\eta }(x_0)\) which is the same that \(\mathcal {L}_2 v_{\infty }^*=0\) in \(B_{\eta }(x_0)\) because \(v_{\infty }\) is lower semicontinuous. As a consequence, \(\mathcal {L}_2 v_{n}\ge 0\) in \(B_{\eta }(x_{0})\). In particular, \(\mathcal {L}_2 v_{n}\ge 0\) in \(B_{\eta }(x_{0})\) for all \(n\ge n_{0}\). Then, by Proposition 4, \(\mathcal {L}_{2} v_{\infty }\ge 0\) in \(B_{\eta }(x_{0})\) in the viscosity sense. Since previously we just proved that \(v_{\infty }\) is a \(\mathcal {L}_{2}\)-viscosity subsolution in \(\Omega \), we obtain that \(v_{\infty }\) is a \(\mathcal {L}_{2}\)-viscosity solution in \(B_{\eta }(x_{0})\). Therefore, \(v_{\infty }\) is a \(\mathcal {L}_{2}\)-viscosity solution in the set \(\{v_{\infty }^{*}<u_{\infty }\}\).
Putting all together, since \(u_{\infty }\) is lower semicontinuous, we obtain
By uniqueness of the obstacle problem we get \(v_{\infty }=v=\overline{\mathcal {O}}(\mathcal {L}_{2},u_{\infty },g)\).
Now, let us define \(u=\underline{\mathcal {O}}(\mathcal {L}_{1},v_{\infty },f)\). Using \(v_n\le v_{\infty }^*\), we have \(u_{n-1}\le u\). Then, taking the limit, \(u_{\infty }\le u\). On the other hand we have that \(\mathcal {L}_{1}u_{\infty }\ge 0\) in \(\Omega \) in viscosity sense due to Proposition 4 and \(u_{\infty }=f\) in \(\partial \Omega \). Moreover, since \(u_{\infty }\) is lower semicontinuous, taking \(\limsup \) in the inequality \(u_{n}\ge v_{n}\), we obtain \(u_{\infty }\ge v_{\infty }^*\) in \(\overline{\Omega }\). Then, using the Lemma 3.7, we get \(u\le u_{\infty }\). Thus, we conclude that \(u_{\infty }=u\).
Finally, we have that \(u_{\infty }\) and \(v_{\infty }\) are lower semicontinuous functions in \(\bar{\Omega }\) by Lemma 3.3. Moreover, \(u_{\infty }\) is continuous in the interior of \(\{u_{\infty }=v_{\infty }^*\}\) because \(u_{\infty }\) is a lower semicontinuous and, by definition, \(v_{\infty }^*\) is an upper semicontinuous function. This also implies that \(v_{\infty }^{*}\) is a continuous function in the interior of \(\{u_{\infty }=v_{\infty }^*\}\). Then, \(v_{\infty }=v_{\infty }^{*}\) there and therefore \(v_{\infty }\) is a continuous function in the interior of the set \(\{u_{\infty }=v_{\infty }\}\). \(\square \)
Remark 6
If \(u_{\infty }\) or \(v_{\infty }\) are continuous at \(\partial \{u_{\infty }=v_{\infty }^{*}\}\), then, both functions \(u_{\infty }\) and \(v_{\infty }\) are continuous in \(\bar{\Omega }\) and therefore the sequences defined in the previous theorem converge uniformly in the whole \(\overline{\Omega }\).
In fact, without loss of generality, suppose that \(u_{\infty }\) is continuous in \(\partial \{u_{\infty }=v_{\infty }^*\}\). We will prove the above using the Comparison Principle for \(\mathcal {L}_{1}\). Since \(\mathcal {L}_{1}u_{\infty }=0\) in \(\{u_{\infty }>v_{\infty }^{*}\}\) in the viscosity sense, we have in particular that \(u_{\infty }^{*}\) is a \(\mathcal {L}_{1}\)-subsolution and \(u_{\infty *}\) a \(\mathcal {L}_{1}\)-supersolution in \(\{u_{\infty }>v_{\infty }^{*}\}\) in the viscosity sense. Moreover, since \(f>g\) in \(\partial \Omega \), the boundary \(\partial \{u_{\infty }>v_{\infty }^{*}\}\) is the disjoint union of \(\partial \Omega \) and \(\partial \{u_{\infty }=v_{\infty }^{*}\}\). In \(\partial \Omega \), we have \(u_{\infty }^{*}=u_{\infty *}\) by construction of \(u_{n}\). And in \(\partial \{u_{\infty }=v_{\infty }^{*}\}\) we have also that \(u_{\infty }^{*}=u_{\infty *}\) because we have supposed \(u_{\infty }\) continuous across that boundary. Then, \(u_{\infty }^{*}=u_{\infty *}\) in \(\partial \{u_{\infty }>v_{\infty }^{*}\}\). Therefore, by the Comparison Principle of \(\mathcal {L}_{1}\), \(u_{\infty }^{*}\le u_{\infty *}\) in \(\overline{\{u_{\infty }>v_{\infty }^{*}\}}\) and therefore \(u_{\infty }\) is continuous in \(\overline{\{u_{\infty }>v_{\infty }^{*}\}}\). Moreover, as we have seen in the proof of the above theorem, \(u_{\infty }\) is also continuous in \(\{u_{\infty }>v_{\infty }^{*}\}\). Then \(u_{\infty }\) is continuous in \(\overline{\Omega }.\) As a consequence, using the latest hypothesis concerning the operator \(\mathcal {L}_{2}\), \(v_{\infty }\) is continuous in the whole \(\overline{\Omega }\), because \(v_{\infty }\) is the solution of the upper obstacle problem with continuous boundary datum and continuous obstacle.
Finally, since the continuous functions \(u_{\infty }\) and \(v_{\infty }\) are the limit of the sequences of continuous functions \(\{u_{n}\}\) and \(\{v_{n}\}\) in a compact set \(\bar{\Omega }\), the convergence is uniform.
Remark 7
As happens in the variational setting, we also have here that the limit depends strongly on the initial function from where we start the iterations. Moreover, the convergence of \((u_{n},v_{n})\) to a solution of the two membranes problem relies on the monotonicity of the sequences \(\{u_{n}\}\) and \(\{v_{n}\}\). This property comes from the fact that the initial function \(v_{0}\) is a \(\mathcal {L}_{2}\)-viscosity subsolution. For this discussion in the variational context, we refer to the arguments given in Remark 4.
Data Availability
There is no data associated with this manuscript.
References
Azevedo, A., Rodrigues, J.F., Santos, L.: The N-membranes problem for quasilinear degenerate systems. Interfaces Free Bound. 7(3), 319–337 (2005)
Barrios, B., Medina, M.: Equivalence of weak and viscosity solutions in fractional non-homogeneous problems. Math. Ann. 381(3–4), 1979–2012 (2021)
Blanc, P., Rossi, J.D.: Game Theory and Partial Differential Equations. De Gruyter Series in Nonlinear Analysis and Applications, vol. 31 (2019). ISBN 978-3-11-061925-6. ISBN 978-3-11-062179-2 (eBook)
Blanc, P., Rossi, J.D.: Obstacle problems and maximal operators. Adv. Nonlinear Stud. 16(2), 355–362 (2016)
Caffarelli, L.A.: The regularity of free boundaries in higher dimensions. Acta Math. 139(3–4), 155–184 (1977)
Caffarelli, L.A.: The obstacle problem revisited. J. Fourier Anal. Appl. 4(4–5), 383–402 (1998)
Caffarelli, L., De Silva, D., Savin, O.: The two membranes problem for different operators. Ann. l’Institut Henri Poincare C Anal. Non Lineaire 34(4), 899–932 (2017)
Caffarelli, L., Duque, L., Vivas, H.: The two membranes problem for fully nonlinear operators. Discrete Contin. Dyn. Syst. 38(12), 6015–6027 (2018)
Caffarelli, L.A., Cabre, X.: Fully Nonlinear Elliptic Equations, vol. 43. Amer. Math. Soc. Colloquium Publications, Providence (1995)
Carillo, S., Chipot, M., Vergara-Caffarelli, G.: The N-membrane problem with nonlocal constraints. J. Math. Anal. Appl. 308(1), 129–139 (2005)
Chipot, M., Vergara-Caffarelli, G.: The N-membranes problem. Appl. Math. Optim. 13(3), 231–249 (1985)
Crandall, M.G., Ishii, H., Lions, P.L.: User guide to viscosity solutions of second order partial differential equations. Bull. Am. Math. Soc. 27, 1–67 (1992)
Daskalopoulos, P., Feehan., P.M.N.: \(C^{1,1}\) regularity for degenerate elliptic obstacle problems in mathematical finance, preprint
de Borbon, M., Del Pezzo, L.M., Ochoa, P.: Weak and viscosity solutions for non-homogeneous fractional equations in Orlicz spaces. Adv. Differ. Equ. 27(11–12), 735–780 (2022)
Heinonen, J., Kilpeläinen, T., Martio, O.: Nonlinear Potential Theory of Degenerate Elliptic Equations. Dover, Illinois (2006)
Ishii, H.: On the equivalence of two notions of weak solutions, viscosity solutions and distribution solutions. Funkcial Ekvac. Ser. Int. 38(1), 101–120 (1995)
Juutinen, P., Lindqvist, P., Manfredi, J.J.: On the equivalence of viscosity solutions and weak solutions for a quasi-linear equation. SIAM J. Math. Anal. 33(3), 699–717 (2001)
Lewicka, M.: A Course on Tug-of-War Games with Random Noise. Introduction and Basic Constructions. Universitext book series, Springer (2020)
Manfredi, J.J., Parviainen, M., Rossi, J.D.: On the definition and properties of p-harmonious functions. Ann. Scuola Norm. Super. Pisa 11, 215–241 (2012)
Medina, M., Ochoa, P.: On viscosity and weak solutions for non-homogeneous p-Laplace equations. Adv. Nonlinear Anal. 8, 468–481 (2019)
Monneau, R.: A Brief Overview on the Obstacle Problem. In: Casacuberta, C., Miró-Roig, R.M., Verdera, J., Xambó-Descamps, S. (eds.) European Congress of Mathematics Progress in Mathematics, vol. 202. Birkhäuser, Basel (2001)
Peres, Y., Schramm, O., Sheffield, S., Wilson, D.: Tug-of-war and the infinity Laplacian. J. Am. Math. Soc. 22, 167–210 (2009)
Peres, Y., Sheffield, S.: Tug-of-war with noise: a game theoretic view of the \(p\)-Laplacian. Duke Math. J. 145(1), 91–120 (2008)
Petrosyan, A., Shagholian, H., Uraltseva, N.: Regularity of Free Boundaries in Obstacle Type Problems. American Mathematical Society, Providence (2012)
Ren, J., Wu, J., Zheng, M.: On the equivalence of viscosity and distribution solutions of second-order PDEs with Neumann boundary conditions. Stoch. Process. Appl. 130(2), 656–676 (2020)
Silvestre, L.: The two membranes problem. Commun. Partial Differ. Equ. 30(1–3), 245–257 (2005)
Vergara-Caffarelli, G.: Regolarita di un problema di disequazioni variazionali relativo a due membrane. Atti Accad. Naz. Lincei Rend. Cl. Sci. Fis. Mat. Natur. (8) 50, 659–662 (1971) (Italian, with English summary)
Vivas, H.: The two membranes problem for fully nonlinear local and nonlocal operators. PhD Thesis dissertation. UT Austin (2018). https://repositories.lib.utexas.edu/handle/2152/74361
Acknowledgements
I. Gonzálvez was supported by the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No. 777822, and by grants CEX2019-000904-S, PID2019-110712GB-I00, PID2020-116949GB-I00, and RED2022-134784-T, all four funded by MCIN/AEI/10.13039/ 501100011033 (Spain). A. Miranda and J. D. Rossi were partially supported by CONICET PIP GI No 11220150100036CO (Argentina), PICT-2018-03183 (Argentina) and UBACyT 20020160100155BA (Argentina).
Funding
Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature. I. Gonzálvez was supported by the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No. 777822, and by grants CEX2019-000904-S, PID2019-110712GBI00, PID2020-116949GB-I00, and RED2022-134784-T, all four funded by MCIN/AEI/10.13039/ 501100011033 (Spain). A. Miranda and J. D. Rossi were partially supported by CONICET PIP GI No 11220150100036CO (Argentina), PICT-2018-03183 (Argentina) and UBACyT 20020160100155BA (Argentina).
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Gonzálvez, I., Miranda, A. & Rossi, J.D. Monotone iterations of two obstacle problems with different operators. J Elliptic Parabol Equ 10, 679–709 (2024). https://doi.org/10.1007/s41808-024-00268-6
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s41808-024-00268-6