Home Optimal rearrangement problem and normalized obstacle problem in the fractional setting
Article Open Access

Optimal rearrangement problem and normalized obstacle problem in the fractional setting

  • Julián Fernández Bonder , Zhiwei Cheng and Hayk Mikayelyan EMAIL logo
Published/Copyright: March 31, 2020

Abstract

We consider an optimal rearrangement minimization problem involving the fractional Laplace operator (–Δ)s, 0 < s < 1, and the Gagliardo seminorm |u|s. We prove the existence of the unique minimizer, analyze its properties as well as derive the non-local and highly non-linear PDE it satisfies

(Δ)sUχ{U0}min{(Δ)sU+;1}=χ{U>0},

which happens to be the fractional analogue of the normalized obstacle problem Δu = χ{u>0}.

MSC 2010: 35R11; 35J60

1 Introduction

One of the classical problems in rearrangement theory is the minimization of the functional

Φ(f)=D|uf|2dx, (1.1)

where uf is the unique solution of the Dirichlet boundary value problem in a bounded domain D

Δuf=fin D,uf=0on D, (1.2)

and f belongs to the set

R¯β=fL(D):0f1,Dfdx=β.

Recall that 𝓡̄β is the closure in the weak* topology of the rearrangement class

Rβ:=fL(D):f=χE,|E|=β.

This minimization problem is related to the stationary heat equation

tu=0Δu=f

in the bounded domain D, which is under the action of the external heat source modeled by the force function f. The Dirichlet boundary condition, u = 0 on ∂D, models the constant temperature on the boundary of D. Different force functions f result different heat distributions uf. The minimizer of the functional (1.1) is the force function from a certain rearrangement class 𝓡, which is resulting the most uniformly distributed heat u.

The problem and its variations, such as the p–harmonic case, has been studied by several authors (see [5, 6, 7, 13, 18]), and the results, for this particular setting, can be formulated in the following theorem.

Theorem 1.1

There exists a unique solution ∈ 𝓡β of the minimization problem (1.1). For the function û = u there exists a constant α > 0 such that

  • 0 < ûα in D,

  • = χ{û<α},

  • û = α in { = 0}.

Moreover, the function U = αû is the minimizer of the functional

J(w)=D|w|2+2w+dx,

among functions wH1(D) with boundary values α on ∂D, and solves the normalized obstacle problem equation

ΔU=χ{U>0}inD. (1.3)

We refrain from presenting here details about the normalized obstacle problem (1.3), which is one of the classical free boundary problems (see [9]).

In recent years there has been a great development of nonlocal diffusion problems, mainly due to some interesting new applications to different fields of the natural sciences such as some physical models [12, 14, 15, 19, 24, 30], finance [2, 20, 27], fluid dynamics [10], ecology [17, 23, 26] and image processing [16]. For a comprehensive reference for non-local diffusion problems and their applications see [4].

It is also worth mentioning that the link between non-local diffusion problems and optimal design problems has been considered in recently in [22, 29].

Among these models for nonlocal diffusion, probably the most important one is given by the fractional laplacian (–Δ)s, (0 < s < 1) that is given (for smooth functions) as

(Δ)su(x):=p.v. Rnu(x)u(y)|xy|n+2sdy=limε0RnBε(x)u(x)u(y)|xy|n+2sdy.

This operator is given as the gradient of the nonlocal Gagliardo energy

|u|s2:=R2n|u(x)u(y)|2|xy|n+2sdxdy, (1.4)

that is the nonlocal analog of the Dirichlet energy u22.

In view of the increasing interest in analyzing nonlocal diffusion models, it naturally comes into attention considering problem (1.2) where the Laplace operator is replace by its fractional counterpart.

Therefore, in this paper, similar to the way it has been done in [25], we will consider an optimal rearrangement problem and derive a related free boundary problem.

More precisely, we consider the minimization problem

Φs(f)min,

where Φs(f)=|uf|s2,uf is the unique solution to

(Δ)suf=f in D and u=0 in Dc

and f ∈ 𝓡̄β.

We show existence and uniqueness of a solution to the fractional rearrangement optimization problem and show that if is the solution and û = u, then 0 ≤ ûα for some α > 0 and, moreover, Û = αû is the unique solution to the normalized fractional obstacle problem

χ{U>0}(Δ)sUχ{U0} in D and U=α in Dc.

Also, we analyze the behavior of such solutions as the fractional parameter s goes to 1.

Finally, we show that the solution to the fractional normalized obstacle problem is also the solution to the (highly nonlinear) equation

(Δ)sUχ{U0}min{(Δ)sU+;1}=χ{U>0},

in D with U = α in Dc.

Organization of the paper

In Section 2 we give a brief introduction to fractional calculus, in Section 3 we analyze the optimal rearrangement problem in the fractional setting and show its relation with the normalized fractional obstacle problem. In Section 4, we study the behavior of the optimal fractional rearrangement problem as s → 1. Finally, in Section 5, we further analyze the normalized fractional obstacle problem and derive a (highly) nonlinear equation that the solution satisfies.

2 Preliminaries

2.1 A very short tour through the basics of the fractional Laplacian

All the results in this section are either well-known or easily proved, so we just recall them for further references without any attempt of giving proofs.

The fractional order Sobolev spaces Hs(ℝn) (for 0 < s < 1) is defined as

Hs(Rn)={vL2(Rn):|v|s2<},

where |⋅|s is the Gagliardo energy given by (1.4). This space is a Hilbert space with inner product given by

(u,v)Hs(Rn)=Rnu(x)v(x)dx+R2n(u(x)u(y))(v(x)v(y))|xy|n+2sdxdy.

For a brief summary of the properties of fractional order Sobolev spaces Hs, we refer to the survey article [12].

Further we denote by Hs(ℝn) the topological dual space of Hs(ℝn) and for a domain D ⊂ ℝn, we denote

H0s(D)={vHs(Rn):v=0 a.e. in Dc}.

Recall that for Lipschitz domains D, the space H0s (D) coincides with the closure of test functions with compact support inside D. We will also denote by Hs(D) the topological dual space of H0s (D). Observe that we have

H0s(D)Hs(Rn)L2(Rn)Hs(Rn)Hs(D),

with continuous inclusions. Moreover, since 𝓓 ⊂ Hs(ℝn) with a dense inclusion, then Hs(ℝn) ⊂ 𝓓′ and, if D is Lipschitz, then Hs(D) ⊂ 𝓓′(D).

Recall that if D is bounded, the following Poincaré type inequality holds true

u2C|u|s for all uH0s(D). (2.1)

An easy fact is that the Gagliardo semi-norm ||s2 is Gâteaux - differentiable in Hs(ℝn) and

limε0ε1(|u+εv|s2|u|s2)=2R2n(u(x)u(y))(v(x)v(y))|xy|n+2sdxdy, (2.2)

for every u, vHs(ℝn).

Furthermore, for a function uHs(ℝn) we can also define the fractional Laplace operator as

(Δ)su(x)=p.v.Rnu(x)u(y)|xy|n+2sdy=limε0(Δ)εsu(x), (2.3)

where

(Δ)εsu(x)=RnBε(x)u(x)u(y)|xy|n+2sdy

and the limit is understood in Hs(ℝn).

Moreover, it holds that

(Δ)su,v=12R2n(u(x)u(y))(v(x)v(y))|xy|n+2sdxdy12|u|s|v|s,

for any u, vHs(ℝn).

The interested reader may consult with [1] for much more on the fractional laplacian and an analysis of differences and similarities with the local Laplace operator.

For any fHs(D) we say uf H0s (D) solves the fractional boundary value problem in D with homogeneous Dirichlet boundary condition

(Δ)suf=fin D,uf=0in Dc, (2.4)

if the equation is satisfied in the sense of distributions. Equivalently, if

12R2n(uf(x)uf(y))(v(x)v(y))|xy|n+2sdxdy=Dfvdx (2.5)

for any v H0s (D). It is easily seen from Riesz representation Theorem, using Poincaré inequality (2.1), that for any fHs(D) there exists a unique uf H0s (D) satisfying (2.5).

To finish these preliminaries we refer the reader to [28], and recall that for fL(D) the weak solution of (2.4), uf Cloc0,δ (D) for some δ > 0, if s 12 and uf Cloc1,δ (D) for some δ > 0, if s > 12 . Moreover, uf is a strong solution to (2.4), namely the limit in (2.3) exists pointwise a.e. and the equation (2.4) is also satisfied pointwise a.e.

3 The optimal fractional rearrangement problem

Let us now introduce the fractional analogue of the optimal rearrangement problem given in Theorem 1.1. Throughout this paper D ⊂ ℝn will always denote a bounded domain. Given fL2(D), let uf be the solution of (2.4) and let us define the functional

Φs(f)=|uf|s2. (3.1)

We are going to consider the minimization of the functional Φs on the closed, convex set 𝓡̄β, for 0 < β < |D|. The main result of this section is the following theorem.

Theorem 3.1

There exists a unique minimizer ∈ 𝓡̄β ∖ 𝓡β such that

Φs(f^)Φs(f)

for any f ∈ 𝓡̄β. Moreover, for some α > 0 the function û = u satisfies the following conditions

{f^<1}{u^=α},{u^<α}{f^=1},f^>0and0u^αinD.

Remark 3.2

Observe that this result shows a remarkable difference with the local optimal rearrangement problem, since the optimal configuration for the fractional case is not a characteristic function. See Theorem 1.1.

For the proof of Theorem 3.1 we need a couple of lemmas.

Lemma 3.3

The set 𝓡̄βL(D) is convex and

ext(R¯β)=Rβ,

where for a convex set C, ext(C) denotes the extreme points of C.

Proof

The proof is standard and is omitted. For a more general result see [8, Lemma 3].□

Lemma 3.4

Let Φs be the functional defined in (3.1). Then Φs : 𝓡̄β → ℝ is strictly convex and sequentially lower semi-continuous with respect to the weak* topology. Moreover, there exists a unique minimizer of the functional Φs in 𝓡̄β.

Proof

The strict convexity is a direct consequence of the linearity uf1+f2 = uf1 + uf2 and the strict convexity of tt2. Moreover, from (2.5), Hölder’s inequality and (2.1), we obtain

|uf|sCf2.

Therefore, fuf is strongly continuous from L2(D) into H0s (D) and hence, Φs is strongly continuous from L2(D) into ℝ. Since Φs is convex, it follows that is sequentially weakly lower semicontinuous.

Finally, observe that if fn ∈ 𝓡̄β is such that fn f weakly* in L(D), then fnf weakly in L2(D), and so

Φs(f)lim infΦs(fn).

To finish the proof just notice that the existence of a minimizer follows from Banach-Alaoglu’s theorem and the uniqueness of the minimizer from the strict convexity of Φs and the convexity of 𝓡̄β.□

Now we are ready to prove Theorem 3.1.

Proof of Theorem 3.1

The proof will be divided into a series of claims.

  1. Du^f^dxDu^fdx for any fR¯β.

    Let us take Ψ: L2(D) → ℝ̄ defined as

    Ψ(f)=Φs(f)+ξR¯β(f),

    where ξ𝓡̄β(f) is the indicator function, i.e.

    ξR¯β(f)=0,if fR¯β,,if fR¯β..

    Observe that Ψ is strictly convex there. Moreover, it is easy to see that minimizes Ψ in L2(D). Thus

    0Ψ(f^),

    where

    Ψ(f^)=gL2(D):Ψ(f)Ψ(f^)Dg(ff^)dx, for any fL2(D)

    is the sub-differential of Ψ at .

    From (2.5) and (2.2) we get that

    Φs(f^)={2u^}.

    Moreover

    ξR¯β(f^)=gL2(D):ξR¯β(f)ξR¯β(f^)Dg(ff^)dx, for any fL2(D)=gL2(D):0Dg(ff^)dx, for any fR¯β.

    Therefore the equation

    0Ψ(f^)=Φs(f^)+ξR¯β(f^),

    implies that

    u^ξR¯β(f^)

    and thus the claim.

  2. There exists a function = χE ∈ 𝓡β such that

    Du^f~dxDu^fdx

    for any f ∈ 𝓡̄β.

    This follows from Claim 1, Lemma 3.3, and the fact that the minimum of the linear functional L(f) = ∫D ûf dx on a bounded closed convex set 𝓡̄β is attained in an extreme point = χE ∈ 𝓡β.

  3. There exists α > 0 such that

    {u^<α}E{u^α}.

    The proof is an immediate consequence of the bathtub principle for . See [21, Theorem 1.14].

  4. f^=1 in {u^<α}.

    The proof is again an immediate consequence of the bathtub principle for .

  5. {u^>α}{f^=0}.

    Since , ∈ 𝓡̄β, we have that

    β=Df^dx={u^<α}f^dx+{u^=α}f^dx+{u^>α}f^dx=Df~dx={u^<α}f~dx+{u^=α}f~dx+{u^>α}f~dx

    Therefore, by Claims 3 and 4, we obtain that

    {u^=α}f^dx+{u^>α}f^dx={u^=α}f~dx. (3.2)

    On the other hand, by Claims 1 and 2, we get

    Du^f^dx=Du^f~dx

    that together with (3.2) give us

    α{u^=α}f~dx=α{u^=α}f^dx+α{u^>α}f^dx{u^=α}u^f^dx+{u^>α}u^f^dx={u^=α}u^f~dx=α{u^=α}f~dx,

    which implies

    α{u^>α}f^dx={u^>α}u^f^dx,

    and thus the claim.

  6. {u^>α}=.

    For β > α let us take ϕ(x) = (û(x) – β)+. Since supp ϕ = ω ⊂ {û > α}, claim 5 implies that

    0=2(Δ)su^,ϕ=R2n(u^(x)u^(y))(ϕ(x)ϕ(y))|xy|n+2sdxdy=ωω(u^(x)u^(y))2|xy|n+2sdxdy0+ωRnω(u^(x)u^(y))(u^(x)β)|xy|n+2sdydx0+Rnωω(u^(x)u^(y))(βu^(y))|xy|n+2sdydx0+RnωRnω(u^(x)u^(y))(00)|xy|n+2sdydx=0.

    Thus, |ω| = |{û > β}| = 0 for any β > α. Moreover, since û Cloc0,δ (D) for some δ > 0, it is easy to see that the claim follows.

  7. |{f^=0}|=0.

    Since (–Δ)s û = L(D) and ≥ 0 it is enough to check > 0 point-wise.

    Taken Claim 4 we need to check this only in the set {û = α}. But

    f^(x)=p.v.Rnu^(x)u^(y)|xy|n+2sdy=limϵ0RnBϵ(x)u^(x)u^(y)|xy|n+2sdy>RnDα|xy|n+2sdy>0.

    This proves the claim.

    The proof of Theorem 3.1 is complete.□

4 The behavior of the optimal rearrangement problem as s → 1

In this section we analyze the behavior of the optimal fractional rearrangement problem as the fractional parameter s goes to 1. For that purpose, we need to consider here the normalizing constant C(n, s) that is defined as

C(n,s)=Rn1cos(ζ1)|ζ|n+2sdζ1

and we need to modify the definitions of the fractional laplacian and of the Gagliardo seminorm accordingly, namely, we consider

|u|s2=C(n,s)R2n|u(x)u(y)|2|xy|n+2sdxdy

and

(Δ)su(x)=p.v. C(n,s)Rnu(x)u(y)|xy|n+2sdy.

It is a well known fact that this normalizing constant behaves like (1–s) for s close to 1. Moreover, the following result holds

Proposition 4.1

Let uL2(ℝn) be fixed. Then we have that

|u|s2u22and(Δ)suΔuass1.

where the first limit is understood as a limit if uH1(ℝn) and as lim inf |u|s2 = ∞ if u ⧸ ∈ H1(ℝn) and the second limit is in the sense of distributions.

For a proof, see for instance [12] and [3].

Moreover, it is shown in [3] the following stronger statement.

Proposition 4.2

Given a sequence sk → 1 and {uk}k∈ℕL2(ℝn) such that

supkNuk2<andsupkN|uk|sk<,

then there exists a function uH1(ℝn) such that (up to a subsequence),

ukustronglyinLloc2(Rn)andu22lim infk|uk|sk2.

Throughout this section, we will denote by s the optimal load for Φs, ûs = us the solution to (2.4). Also, denote Φ(f) as

Φ(f)=D|uf|2dx,

where in this section, uf will denote the solution to

Δuf=fin Du=0on D.

Finally, denote by ∈ 𝓡̄β the solution to the minimization problem

Φ(f^)=inffR¯βΦ(f).

So the main result in this section is the following:

Theorem 4.3

Under the above notations, s weakly* in L as s → 1. Moreover we also obtain that

Φs(f^s)Φ(f^)andu^su^stronglyinL2(D),

as s → 1.

For the proof of Theorem 4.3 we need the concept of Γ–convergece. This concept was introduced by De Giorgi in the 60s and is now a well understood tool to deal with the convergence of minimum problems. For a throughout introduction to the subject, we cite [11]. Let us recall now the definition of Γ-$convergence and some of its properties.

Definition 4.4

Let X be a metric space and Fn, F : X → ℝ̄. We say that Fn Γ–converges to F, and is denoted by Fn Γ F, is the following two inequalities hold true

  • (lim inf–inequality) For any xX and any sequence {xn}n∈ℕX such that xnx in X, it holds that

    F(x)lim infnFn(xn).
  • (lim sup–inequality) For any xX, there exists a sequence {yn}n∈ℕX such that ynx in X and

    F(x)lim supnFn(yn).

The main feature of the Γ–convergence is that it implies the convergence of minima. In fact we have the following:

Theorem 4.5

Let X be a metric space and Fn, F: X be functions such that Fn Γ F. Moreover, assume that for each n ∈ ℕ, there exists xnX such that

Fn(xn)=infXFn

and that {xn}n∈ℕX is precompact. Then

limninfXFn=infXF

and every accumulation point of the sequence {xn}n∈ℕ is a minimum point of F.

The proof of Theorem 4.5 is easy and can be found in [11].

The following result is key in the proof of Theorem 4.3.

Theorem 4.6

Given 0 < s < 1, let fsL2(Ω) be such that fsf weakly in L2(Ω) and let us H0s (Ω) and u H01 (Ω) be the solutions to

(Δ)sus=fsinΩ,us=0inRnΩ

and

Δu=finΩ,u=0onΩ

respectively.

Then usu strongly in L2(Ω). Moreover,

|us|su2.

Proof

Let Fs, F: L2(ℝn) → ℝ̄ given by

Fs(v):=12|v|s2Ωfsvdxif vH0s(Ω),+otherwise

and

F(v):=12v22Ωfvdxif vH01(Ω),+otherwise.

Since ∫Ω fs vs dx → ∫Ω fv dx if vsv strongly in L2(ℝn), from Propositions 4.1 and 4.2 we can conclude that Fs Γ F as s → 1.

Now, observe that

Fs(us)=infvL2(Rn)Fs(v)andF(u)=infvL2(Rn)F(v).

The trivial estimate |us|s ≤ ∥fs2 imply that, for any sk → 1, the sequence {usk}k∈ℕL2(ℝn) is precompact. Then from Theorem 4.5 we obtain that usu strongly in L2(Ω).

Finally,

lim|us|s2=limΩfsusdx=Ωfudx=u22.

This completes the proof.□

Now we are ready to prove the main result of the section.

Proof of Theorem 4.3

Let s ∈ 𝓡̄β be the optimal load for Φs. Observe that, for a subsequence, s f weakly* in L(Ω) for some f ∈ 𝓡̄β. Moreover, this convergence also holds weakly in L2(Ω).

From Theorem 4.6 we have that ûsuf strongly in L2(Ω) and using Theorem 4.2 we get

infR¯βΦΦ(f)=ufs2lim inf|u^s|s2=lim infinfR¯βΦs.

On the other hand, let ∈ 𝓡̄β be the optimal load for Φ. Then, using the final part of Theorem 4.6, we obtain

lim supinfR¯βΦslimΦs(f^)=Φ(f^)=infR¯βΦ.

The proof is complete.□

5 The normalized fractional obstacle problem

This section is devoted to the study of the connection between the solutions to the optimal fractional rearrangement problem consider in Section 3 with solutions of the normalized fractional obstacle problem.

The fractional analogue of the classical obstacle problem has been well known in the literature, however its so called normalized version, i.e., the equation

Δu=χ{u>0}, (5.1)

has not been considered. Here we find the corresponding fractional analog of (5.1) and prove that the solution of the fractional rearrangement problem is a solution of the fractional normalized obstacle problem.

Our first result is the following.

Theorem 5.1

Let ∈ 𝓡̄β be the solution to the optimal fractional rearrangement problem and û : = u H0s (D) be given by (2.4). Let α > 0 be the constant given in Theorem 3.1. Then the function Û : = αû minimizes the functional

J(v)=|v|s2+Dv+dx

over the set Hα = {v Hlocs (ℝn): vα H0s (D)}. Moreover, Û verifies the inequalities

χ{U>0}(Δ)sUχ{U0}inD (5.2)

in the sense of distributions.

Finally, the minimizer of J in Hα is unique and is the unique solution to the inequality (5.2).

Proof

Let

I(v)=|v|s2+Df^vdx

and observe that, since 0 ≤ ≤ 1, for any vHα it follows that J(v) ≥ I(v).

Next, observe that I(Û) = J(Û) and so the set of inequalities

J(v)I(v)I(U^)=J(U^),for any vHα

imply the desired result.

Next, observe that the inequalities

χ{U^>0}(Δ)sU^χ{U^0}.

are the Euler-Lagrange equation for the functional J based on the variation uε(x) = u(x) + εϕ(x), with ϕ Cc (D).

Now, the uniqueness of minimizer for J is an immediate consequence of the strict convexity of J.

Assume that the function U satisfies the inequalities (5.2), but the unique minimizer of the convex functional J is the function VU.

Since J is strictly convex and J(V) < J(U) by taking Uε = U + ε(VU) we will obtain

limε0+J(Uε)J(U)ε<0.

Thus for ψ = VU

0>limε0+J(Uε)J(U)ε=R2n(u(x)u(y))(ψ(x)ψ(y))|xy|n+2sdxdy+Dχ{U>0}ψ+χ{U=0}ψ+dx=R2n(u(x)u(y))(ψ+(x)ψ+(y))|xy|n+2sdxdy+Dχ{U0}ψ+dx0(R2n(u(x)u(y))(ψ(x)ψ(y))|xy|n+2sdxdy+Dχ{U>0}ψdx0)0,

where the last inequality follows from (5.2). This is a contradiction and the result follows.□

Remark 5.2

This result again shows an interesting difference between the classical obstacle problem and the fractional normalized version. Observe that in the positivity set, we still have −(−Δ)s Û = 1, but in the zero set the function Û is not sharmonic (even if it is identically zero!). The free boundary condition on {Û > 0} is given by the fact that (−Δ)s Û is a function bounded by 0 and 1 across the free boundary.

The results in Theorem 5.1 are not completely satisfactory, since we do not obtain an equation satisfied by Û but only the inequalities (5.2).

Our last result shows that in fact Û is the solution to a fully nonlinear equation.

Theorem 5.3

Let Û be solution of the normalized fractional obstacle problem given by Theorem 5.1. Then Û is a solution to

(Δ)sUχ{U0}min{(Δ)sU+;1}=χ{U>0},inD,U=αinDc. (5.3)

Moreover, problem (5.3) is equivalent to (5.2). Finally, U verifies (5.3) if and only if it is a minimizer of J in Hα, where J and Hα are given in Theorem 5.1.

Before we start with the proof, let us observe that for uHs(ℝn),

|u±|s|u|s

and hence (−Δ)s u±Hs(ℝn). On the other hand (−Δ)s u+ is a distribution and the expression

min{(Δ)su+;1)=max{(Δ)su+,1}=1((Δ)su++1)+

makes in general no sense, unless (−Δ)s u+ is a signed measure in D. Let us further observe that since

χ{u0}(Δ)su+0,

we need to search for solutions of (5.3) only among functions u, such that (−Δ)s u ≤ 0 in D. This leads us to the introduction of fractional subharmonic functions in D, which form a convex subset of Hs(D)

Hsubs(D)=uHlocs(Rn):(Δ)su0 in D.

Here the inequality (−Δ)s u ≤ 0 should be understood in the sense of distributions.

The following lemma is essential for the equation (5.3) to make sense.

Lemma 5.4

Let u Hsubs (D). Then u+ Hsubs (D).

Proof

If u is smooth, then the fractional laplacian has pointwise values. In this case, we simply compute:

  • For x ∈ {u ≤ 0},

    (Δ)su+(x)=p.v.Rnu+(y)|xy|n+2sdy0.
  • For x ∈ {u > 0},

    (Δ)su+(x)=p.v.Rnu(x)u+(y)|xy|n+2sdy=p.v.Rnu(x)u(y)|xy|n+2sdyp.v.Rnu(y)|xy|n+2sdy0.

For a general uHs(ℝn), we take {ρε}ε>0 a smooth family of approximations of the identity such that ρε(z) = ρε(−z) and define uε = uρε.

The result of the lemma will follow from the identity

(Δ)suε,ϕ=(Δ)su,ϕε, (5.4)

for every ϕ Cc (ℝn).

Indeed, assuming (5.4), if (−Δ)s u ≤ 0, then (−Δ)s uε ≤ 0 for every ε > 0. Hence, from the smooth case we conclude that (−Δ)s uε+ ≤ 0 and since uε+ u+ in Hs the result is proved.

It remains to prove (5.4). For that purpose, it is useful to introduce the notation

Dsu(x,y)=u(x)u(y)|xy|n2+s,

the Hölder quotient of order s of u. Then, using that ρε(−z) = ρε(z), we observe that

(Δ)2uε,ϕ=R2nDsuε(x,y)Dsϕ(x.y)dxdy=R2nRnDsu(xz,yz)ρε(z)Dsϕ(x,y)dzdxdy=R2nRnDsu(x,y)ρε(z)Dsϕ(x+z,y+z)dzdxdy=R2nRnDsu(x,y)ρε(z)Dsϕ(xz,yz)dzdxdy=R2nDsu(x,y)Dsϕε(x,y)dxdy=(Δ)2u,ϕε.

The proof is now complete.□

Corollary 5.5

If u Hsubs (D) then min{−(−Δ)s u+;1} ∈ L(D).

Corollary 5.5 allows us to formulate the following normalized fractional obstacle problem:

For α > 0 solve

(Δ)sUχ{U0}min{(Δ)sU+;1}=χ{U>0}, (5.5)

among functions UHs(ℝn), such that U = α in Dc and U Hsubs (D).

The weak formulation of the equation (5.5) is

12R2n(U(x)U(y))(ϕ(x)ϕ(y))|xy|n+2sdxdy=D[χ{U0}(x)min{(Δ)sU+;1}+χ{U>0}(x)]ϕ(x)dx,

for any ϕ H0s (D).

Now we are ready to prove Theorem 5.3.

Proof

Proof of Theorem 5.3. We only need to show that problems (5.2) and (5.5) are equivalent

For convenience let us break down the proof into several claims.

  1. Assume that U Hlocs (ℝn) is a solution of (5.5). Then $U ≥ 0.

    Observe first the following general fact:

    (U(x)U(y))(U(x)U(y))=(U+(x)U+(y))(U(x)U(y))(U(x)U(y))2.

    This simple identity implies that

    (Δ)sU,U=(Δ)sU+,U|U|s2. (5.6)

    Now, let us take ϕ = U as a test function in the weak formulation of (5.5). Then we obtain

    (Δ)sU,U=D[χ{U0}min((Δ)sU+;1)+χ{U>0}]Udx=Dmin{(Δ)sU+;1}Udx=Dmax{(Δ)2U+;1}Udx (5.7)

    But

    Dmax{(Δ)sU+;1}Udx(Δ)sU+,U. (5.8)

    Therefore, combining (5.6), (5.7) and (5.8), we arrive at

    |U|s20,

    and so the claim is proved.

  2. (5.5) implies (5.2).

    It is immediate from Claim 1.

  3. (5.2) implies U ≥ 0.

    The argument is similar to the one of Claim 6 in Theorem 3.1. Let U Hlocs (ℝn) be a solution to (5.2). Take β < 0 and ϕ = (Uβ), so ω = supp ϕ ⊂ {U < 0}. Then

    0=2(Δ)sU,ϕ=R2n(U(x)U(y))(ϕ(x)ϕ(y))|xy|n+2sdxdy=ωω(U(x)U(y))(U(y)U(x))|xy|n+2sdxdy0+ωRω(U(x)U(y))(βU(x))|xy|n+2sdxdy0+Rωω(U(x)U(y))(U(y)β)|xy|n+2sdxdy0+RωRω(U(x)U(y))(00)|xy|n+2sdxdy=00.

    Thus, |ω| = |{U < β}| = 0 for any β < 0.

  4. (5.2) implies (5.5).

    Can be verified directly.□

Acknowledgement

The research of Zhiwei Cheng and Hayk Mikayelyan has been partly supported by the National Science Foundation of China (grant no.1161101064). Julián F. Bonder is supported by by grants UBACyT 20020130100283BA, CONICET PIP 11220150100032CO and ANPCyT PICT 2012-0153.

This research was done while J. F. Bonder was a visiting Professor at University of Nottingham at Ningbo, China (UNNC). He wants to thank the institution for the support, the atmosphere and the hospitality that make the stay so enjoyable.

References

[1] Nicola Abatangelo and Enrico Valdinoci, Getting acquainted with the fractional Laplacian, Contemporary research in elliptic PDEs and related topics, volume 33 of Springer INdAM Ser., pages 1–105. Springer, Cham, 2019.10.1007/978-3-030-18921-1_1Search in Google Scholar

[2] Vedat Akgiray and G. Geoffrey Booth, The siable-law model of stock returns, Journal of Business & Economic Statistics, 6(1):51–57, 1988.10.1080/07350015.1988.10509636Search in Google Scholar

[3] Jean Bourgain, Haim Brezis, and Petru Mironescu, Another look at Sobolev spaces, Optimal control and partial differential equations, pages 439–455. IOS, Amsterdam, 2001.Search in Google Scholar

[4] Claudia Bucur and Enrico Valdinoci, Nonlocal diffusion and applications, volume 20 of Lecture Notes of the Unione Matematica Italiana, Springer, [Cham]; Unione Matematica Italiana, Bologna, 2016.10.1007/978-3-319-28739-3Search in Google Scholar

[5] G. R. Burton, Rearrangements of functions, maximization of convex functionals, and vortex rings, Math. Ann., 276(2):225–253, 1987.10.1007/BF01450739Search in Google Scholar

[6] G. R. Burton, Variational problems on classes of rearrangements and multiple configurations for steady vortices, Ann. Inst. H. Poincaré Anal. Non Linéaire, 6(4):295–319, 1989.10.1016/s0294-1449(16)30320-1Search in Google Scholar

[7] G. R. Burton and J. B. McLeod, Maximisation and minimisation on classes of rearrangements, Proc. Roy. Soc. Edinburgh Sect. A, 119(3-4):287–300, 1991.10.1017/S0308210500014840Search in Google Scholar

[8] G. R. Burton and E. P. Ryan, On reachable sets and extremal rearrangements of control functions, SIAM J. Control Optim., 26(6):1481–1489, 1988.10.1137/0326086Search in Google Scholar

[9] L. A. Caffarelli, The obstacle problem revisited, J. Fourier Anal. Appl., 4(4-5):383–402, 1998.10.1007/BF02498216Search in Google Scholar

[10] Peter Constantin, Euler equations, Navier-Stokes equations and turbulence, Mathematical foundation of turbulent viscous flows, volume 1871 of Lecture Notes in Math., pages 1–43. Springer, Berlin, 2006.10.1007/11545989_1Search in Google Scholar

[11] Gianni Dal Maso, An introduction to Γ-convergence, volume 8 of Progress in Nonlinear Differential Equations and their Applications, Birkhäuser Boston, Inc., Boston, MA, 1993.10.1007/978-1-4612-0327-8Search in Google Scholar

[12] Eleonora Di Nezza, Giampiero Palatucci, and Enrico Valdinoci, Hitchhiker’s guide to the fractional Sobolev spaces, Bull. Sci. Math., 136(5):521–573, 2012.10.1016/j.bulsci.2011.12.004Search in Google Scholar

[13] Behrouz Emamizadeh and Yichen Liu, Constrained and unconstrained rearrangement minimization problems related to the p-Laplace operator, Israel J. Math., 206(1):281–298, 2015.10.1007/s11856-014-1141-9Search in Google Scholar

[14] A. Cemal Eringen, Nonlocal continuum field theories, Springer-Verlag, New York, 2002.Search in Google Scholar

[15] Giambattista Giacomin and Joel L. Lebowitz, Phase segregation dynamics in particle systems with long range interactions. I. Macroscopic limits, J. Statist. Phys., 87(1-2):37–61, 1997.10.1007/BF02181479Search in Google Scholar

[16] Guy Gilboa and Stanley Osher, Nonlocal operators with applications to image processing, Multiscale Model. Simul., 7(3):1005–1028, 2008.10.1137/070698592Search in Google Scholar

[17] Nicolas et al. Humphries, Environmental context explains lévy and brownian movement patterns of marine predators, Nature, 465:1066–1069, 2010.10.1038/nature09116Search in Google Scholar

[18] Bernhard Kawohl, Rearrangements and convexity of level sets in PDE, volume 1150 of Lecture Notes in Mathematics, Springer-Verlag, Berlin, 1985.10.1007/BFb0075060Search in Google Scholar

[19] Nikolai Laskin, Fractional quantum mechanics and Lévy path integrals, Phys. Lett. A, 268(4-6):298–305, 2000.10.1016/S0375-9601(00)00201-2Search in Google Scholar

[20] Sergei Levendorski, Pricing of the american put under lévy processes, International Journal of Theoretical and Applied Finance, 7(03):303–335, 2004.10.1142/S0219024904002463Search in Google Scholar

[21] Elliott H. Lieb and Michael Loss, Analysis, volume 14 of Graduate Studies in Mathematics, American Mathematical Society, Providence, RI, second edition, 2001.Search in Google Scholar

[22] Diego Marcon and Rafayel Teymurazyan, A free boundary problem with fractional diffusion. arXiv e-prints, page arXiv:2001.01145, Jan 2020.Search in Google Scholar

[23] A. Massaccesi and E. Valdinoci, Is a nonlocal diffusion strategy convenient for biological populations in competition? ArXiv e-prints, March 2015.10.1007/s00285-016-1019-zSearch in Google Scholar

[24] Ralf Metzler and Joseph Klafter, The random walk’s guide to anomalous diffusion: a fractional dynamics approach, Phys. Rep., 339(1):77, 2000.10.1016/S0370-1573(00)00070-3Search in Google Scholar

[25] Hayk Mikayelyan, Cylindrical optimal rearrangement problem leading to a new type obstacle problem, ESAIM Control Optim. Calc. Var., 24(2):859–872, 2018.10.1051/cocv/2017047Search in Google Scholar

[26] A. M. Reynolds and C. J. Rhodes, The Lévy flight paradigm: random search patterns and mechanisms. Ecology, 90(4):877–887, 2009.10.1890/08-0153.1Search in Google Scholar PubMed

[27] Wim Schoutens, Lévy Processes in Finance: Pricing Financial Derivatives, Willey Series in Probability and Statistics, Willey, New York, 2003.10.1002/0470870230Search in Google Scholar

[28] Luis Silvestre, Regularity of the obstacle problem for a fractional power of the laplace operator, Communications on Pure and Applied Mathematics, 60(1):67–112, 2007.10.1002/cpa.20153Search in Google Scholar

[29] Eduardo V. Teixeira and Rafayel Teymurazyan, Optimal design problems with fractional diffusions, J. Lond. Math. Soc. (2), 92(2):338–352, 2015.10.1112/jlms/jdv034Search in Google Scholar

[30] Kun Zhou and Qiang Du, Mathematical and numerical analysis of linear peridynamic models with nonlocal boundary conditions, SIAM J. Numer. Anal., 48(5):1759–1780, 2010.10.1137/090781267Search in Google Scholar

Received: 2019-05-21
Accepted: 2020-02-20
Published Online: 2020-03-31

© 2020 J.F. Bonder et al., published by De Gruyter

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

Articles in the same Issue

  1. Frontmatter
  2. On the moving plane method for boundary blow-up solutions to semilinear elliptic equations
  3. Regularity of solutions of the parabolic normalized p-Laplace equation
  4. Cahn–Hilliard equation on the boundary with bulk condition of Allen–Cahn type
  5. Blow-up solutions for fully nonlinear equations: Existence, asymptotic estimates and uniqueness
  6. Radon measure-valued solutions of first order scalar conservation laws
  7. Ground state solutions for a semilinear elliptic problem with critical-subcritical growth
  8. Generalized solutions of variational problems and applications
  9. Existence and non-existence results for Kirchhoff-type problems with convolution nonlinearity
  10. Nonlinear Sherman-type inequalities
  11. Global regularity for systems with p-structure depending on the symmetric gradient
  12. Homogenization of a net of periodic critically scaled boundary obstacles related to reverse osmosis “nano-composite” membranes
  13. Noncoercive resonant (p,2)-equations with concave terms
  14. Evolutionary quasi-variational and variational inequalities with constraints on the derivatives
  15. Sharp estimates on the first Dirichlet eigenvalue of nonlinear elliptic operators via maximum principle
  16. Localization and multiplicity in the homogenization of nonlinear problems
  17. Remarks on a nonlinear nonlocal operator in Orlicz spaces
  18. A Picone identity for variable exponent operators and applications
  19. On the weakly degenerate Allen-Cahn equation
  20. Continuity results for parametric nonlinear singular Dirichlet problems
  21. Construction of type I blowup solutions for a higher order semilinear parabolic equation
  22. Singularly perturbed Choquard equations with nonlinearity satisfying Berestycki-Lions assumptions
  23. Comparison results for nonlinear divergence structure elliptic PDE’s
  24. Constant sign and nodal solutions for parametric (p, 2)-equations
  25. Monotonicity formulas for coupled elliptic gradient systems with applications
  26. Berestycki-Lions conditions on ground state solutions for a Nonlinear Schrödinger equation with variable potentials
  27. A class of semipositone p-Laplacian problems with a critical growth reaction term
  28. The role of superlinear damping in the construction of solutions to drift-diffusion problems with initial data in L1
  29. Reconstruction of Tesla micro-valve using topological sensitivity analysis
  30. Lewy-Stampacchia’s inequality for a pseudomonotone parabolic problem
  31. Global well-posedness of nonlinear wave equation with weak and strong damping terms and logarithmic source term
  32. Regularity Criteria for Navier-Stokes Equations with Slip Boundary Conditions on Non-flat Boundaries via Two Velocity Components
  33. Homoclinics for singular strong force Lagrangian systems
  34. A constructive method for convex solutions of a class of nonlinear Black-Scholes equations
  35. On a class of nonlocal nonlinear Schrödinger equations with potential well
  36. Superlinear Schrödinger–Kirchhoff type problems involving the fractional p–Laplacian and critical exponent
  37. Regularity for minimizers for functionals of double phase with variable exponents
  38. Boundary blow-up solutions to the Monge-Ampère equation: Sharp conditions and asymptotic behavior
  39. Homogenisation with error estimates of attractors for damped semi-linear anisotropic wave equations
  40. A-priori bounds for quasilinear problems in critical dimension
  41. Critical growth elliptic problems involving Hardy-Littlewood-Sobolev critical exponent in non-contractible domains
  42. On the Sobolev space of functions with derivative of logarithmic order
  43. On a logarithmic Hartree equation
  44. Critical elliptic systems involving multiple strongly–coupled Hardy–type terms
  45. Sharp conditions of global existence for nonlinear Schrödinger equation with a harmonic potential
  46. Existence for (p, q) critical systems in the Heisenberg group
  47. Periodic traveling fronts for partially degenerate reaction-diffusion systems with bistable and time-periodic nonlinearity
  48. Some hemivariational inequalities in the Euclidean space
  49. Existence of standing waves for quasi-linear Schrödinger equations on Tn
  50. Periodic solutions for second order differential equations with indefinite singularities
  51. On the Hölder continuity for a class of vectorial problems
  52. Bifurcations of nontrivial solutions of a cubic Helmholtz system
  53. On the exact multiplicity of stable ground states of non-Lipschitz semilinear elliptic equations for some classes of starshaped sets
  54. Sign-changing multi-bump solutions for the Chern-Simons-Schrödinger equations in ℝ2
  55. Positive solutions for diffusive Logistic equation with refuge
  56. Null controllability for a degenerate population model in divergence form via Carleman estimates
  57. Eigenvalues for a class of singular problems involving p(x)-Biharmonic operator and q(x)-Hardy potential
  58. On the convergence analysis of a time dependent elliptic equation with discontinuous coefficients
  59. Multiplicity and concentration results for magnetic relativistic Schrödinger equations
  60. Solvability of an infinite system of nonlinear integral equations of Volterra-Hammerstein type
  61. The superposition operator in the space of functions continuous and converging at infinity on the real half-axis
  62. Estimates by gap potentials of free homotopy decompositions of critical Sobolev maps
  63. Pseudo almost periodic solutions for a class of differential equation with delays depending on state
  64. Normalized multi-bump solutions for saturable Schrödinger equations
  65. Some inequalities and superposition operator in the space of regulated functions
  66. Area Integral Characterization of Hardy space H1L related to Degenerate Schrödinger Operators
  67. Bifurcation of time-periodic solutions for the incompressible flow of nematic liquid crystals in three dimension
  68. Morrey estimates for a class of elliptic equations with drift term
  69. A singularity as a break point for the multiplicity of solutions to quasilinear elliptic problems
  70. Global and non global solutions for a class of coupled parabolic systems
  71. On the analysis of a geometrically selective turbulence model
  72. Multiplicity of positive solutions for quasilinear elliptic equations involving critical nonlinearity
  73. Lack of smoothing for bounded solutions of a semilinear parabolic equation
  74. Gradient estimates for the fundamental solution of Lévy type operator
  75. π/4-tangentiality of solutions for one-dimensional Minkowski-curvature problems
  76. On the existence and multiplicity of solutions to fractional Lane-Emden elliptic systems involving measures
  77. Anisotropic problems with unbalanced growth
  78. On a fractional thin film equation
  79. Minimum action solutions of nonhomogeneous Schrödinger equations
  80. Global existence and blow-up of weak solutions for a class of fractional p-Laplacian evolution equations
  81. Optimal rearrangement problem and normalized obstacle problem in the fractional setting
  82. A few problems connected with invariant measures of Markov maps - verification of some claims and opinions that circulate in the literature
Downloaded on 31.10.2025 from https://www.degruyterbrill.com/document/doi/10.1515/anona-2020-0067/html
Scroll to top button