Optimal design problems for the first $p-$fractional eigenvalue with mixed boundary conditions
Julian Fernandez Bonder, Julio D. Rossi, Juan F. Spedaletti

TL;DR
This paper investigates the optimal shape design for the first eigenvalue of the fractional p-Laplacian with mixed boundary conditions, establishing existence results and asymptotic bounds as the fractional parameter approaches 1.
Contribution
It introduces a new optimal shape problem for fractional p-Laplacian eigenvalues with mixed boundary conditions, proving existence and analyzing asymptotic behavior.
Findings
Existence of an optimal design for the shape problem.
Asymptotic bounds for the eigenvalue as s approaches 1.
Asymptotic bounds are independent of the measure constraint .
Abstract
In this paper we study an optimal shape design problem for the first eigenvalue of the fractional laplacian with mixed boundary conditions. The optimization variable is the set where the Dirichlet condition is imposed (that is restricted to have measure equal than a prescribed quantity, ). We show existence of an optimal design and analyze the asymptotic behavior when the fractional parameter obtaining asymptotic bounds that are independent of .
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsNonlinear Partial Differential Equations · Advanced Mathematical Modeling in Engineering · Numerical methods in inverse problems
Optimal design problems for the first fractional eigenvalue
with mixed boundary conditions
Julián Fernández Bonder, Julio D. Rossi and Juan F. Spedaletti
Departamento de Matemática FCEN - Universidad de Buenos Aires and IMAS - CONICET. Ciudad Universitaria, Pabellón I (C1428EGA) Av. Cantilo 2160. Buenos Aires, Argentina.
[email protected] http://mate.dm.uba.ar/ jfbonder Departamento de Matemática FCEN - Universidad de Buenos Aires and IMAS - CONICET. Ciudad Universitaria, Pabellón I (C1428EGA) Av. Cantilo 2160. Buenos Aires, Argentina.
[email protected] http://mate.dm.uba.ar/ jrossi Departamento de Matemática, Universidad Nacional de San Luis and IMASL - CONICET. Ejército de los Andes 950 (D5700HHW), San Luis, Argentina.
Abstract.
In this paper we study an optimal shape design problem for the first eigenvalue of the fractional laplacian with mixed boundary conditions. The optimization variable is the set where the Dirichlet condition is imposed (that is restricted to have measure equal than a prescribed quantity, ). We show existence of an optimal design and analyze the asymptotic behavior when the fractional parameter obtaining asymptotic bounds that are independent of .
Key words and phrases:
Shape optimization, Fractional laplacian, Gamma convergence
2010 Mathematics Subject Classification:
35P30, 35J92, 49R05
1. Introduction
The purpose of this paper is to analyze some optimization problems related to the first eigenvalue of the fractional laplacian with mixed boundary conditions of Neumann and Dirichlet type, where the optimization variable is the region in which the Dirichlet condition is imposed.
Let us be more specific. Let be an open, bounded and Lipschitz domain. The complement of is then divided into two sets
[TABLE]
and , for , we consider the fractional laplacian operator of order with homogeneous Neumann condition on and homogeneous Dirichlet condition on . That is
[TABLE]
where is a normalization constant given by
[TABLE]
here is the Gamma function.
This operator is thought as acting in the fractional Sobolev space
[TABLE]
where, as usual, the fractional Sobolev space is defined as
[TABLE]
Hence, and
[TABLE]
With these definitions, the associated eigenvalue problem is
[TABLE]
The set does not enter in our formulation and is interpreted as a Neumann condition over it, since the system does not interact with .
This problem is analogous in the nonlocal fractional setting to the eigenvalue problem for the local Laplacian
[TABLE]
Here . When we have the usual Dirichlet problem (with a positive first eigenvalue, ); while for we have the Neumann problem (and here the first eigenvalue is zero, ).
Recall that the fractional Sobolev space is a Banach space when one considers the norm
[TABLE]
Moreover, for it is a reflexive, uniformly convex and separable Banach space. The term
[TABLE]
is called the Gagliardo seminorm of . We refer the interested reader to [6] for a throughout introduction to these spaces and operators.
It is straightforward to see that is an eigenvalue of (1.2) if and only if is a critical value of the functional
[TABLE]
restricted to the unit ball of . Moreover, eigenfunctions of (1.2) associated to are critical points of on the unit ball of with critical value .
Of particular importance is the first eigenvalue of (1.2), that is given by
[TABLE]
It is an easy consequence of the direct method of the Calculus of Variations (c.f. Section 2) that the following assertions hold:
- •
the above infimum is achieved, and we can assume that a function that realize the infimum (we will called such a function an extremal) is normalized in , that is, ;
- •
the number is in fact an eigenvalue of (1.2);
- •
is the first (smallest) eigenvalue, i.e. if is an eigenvalue of (1.2), then ;
- •
any eigenfunction associated to has constant sign;
- •
if .
The main problem that we address here is the optimization of the principal eigenvalue with respect to the region where the Dirichlet data is imposed. With that in mind, we fix a constant , define the class
[TABLE]
and consider the optimization problems
[TABLE]
As we will see, since pushing to infinity force that in the limit, so in order to recover a nontrivial constant for the minimization problem one has to restrict the sets to be uniformly bounded. Therefore, given large, we define
[TABLE]
and
[TABLE]
This value is strictly positive.
The main results in this paper are contained in the following theorem.
Theorem 1.1**.**
Let be bounded and open. Take and fix . Let , , be given by (1.5) and (1.6). Then the following hold:
- •
* while ;*
- •
*, where is the first Dirichlet eigenvalue of the local *laplacian in ;
- •
.
Moreover, for every there exists an optimal set for the constant .
Finally, given a sequence of quasi-optimal sets for , i.e. such that
[TABLE]
where as , we have that “surrounds” the boundary of in the sense that for any and any , there exists such that if , .
Notice that the first Neumann eigenvalue for the Laplacian in is (with a constant function as eigenfunction). Therefore, as we have that for every set of measure in
[TABLE]
we conclude that for close to we have that asymptotically lies between the first Dirichlet eigenvalue and the first Neumann eigenvalue for the local Laplacian in . In fact, for every there exists such that for it holds
[TABLE]
The surprising fact of these asymptotic bounds as is that they are independent of the size of the Dirichlet part in our nonlocal problem and they are also independent of the radius of the ball that bounds everything.
To emphasize that in the limit as for the lower bound for the eigenvalues we obtain a local problem with Neumann boundary conditions in we remark that with minor modifications of our arguments we can deal with the eigenvalue problem with a potential . In fact, let us consider a potential such that
[TABLE]
for some numbers and then we take
[TABLE]
with
[TABLE]
Associated with this functional we have the optimal constants
[TABLE]
and
[TABLE]
In this case, we obtain that
[TABLE]
as . Here and are the first eigenvalues for the local Laplacian with the potential with Dirichlet and Neumann boundary condition respectively, that are given by
[TABLE]
and
[TABLE]
Notice that in this case .
When we consider a potential one can also check that if we don’t constraint into a large ball we still have
[TABLE]
In fact, this limit can be deduced taking the limit as in the following two inequalities
[TABLE]
and
[TABLE]
We include some details in Section 5.
With these preliminaries, our second result is the following:
Theorem 1.2**.**
Let be bounded and open. Let be fixed and take . Let a potential such that (1.7) is satisfied.
Let , , be the constants defined in (1.8) and (1.9). Then the following hold
- •
, where is given by (1.10);
- •
, where is given by (1.11).
- •
.
A very brief comment on related bibliography is in order. Optimal configurations related to eigenvalue problems is by now a classical subject, just to mention a few references we quote [3, 4, 5, 8, 11, 12]. On the other hand, nonlocal problems are quite popular nowadays, we just refer to [6] and for references concerning eigenvalues for the nonlocal Laplacian to [2, 7] and references therein.
Organization of the paper
After this introduction, the rest of the paper is organized as follows: In Section 2 we revise some preliminary notions on fractional Sobolev spaces that are needed in the paper. In Section 3 we study the maximization problem and in Section 4 the minimization problem. Finally, in Section 5, we prove Theorem 1.2. Since the proof is similar to the one of Theorem 1.1 we only sketch it and stress the differences.
2. Preliminaries
In this section, we review some definitions on fractional Sobolev spaces and on the fractional Laplace operator. We believe that most of these results are known to experts and constitute part of the “folklore” on the subject, but we have chosen to include some proofs of the facts that are needed just for the reader’s convenience.
2.1. A probabilistic interpretation for the mixed boundary conditions
Let and take such that and .
Suppose that measures the density of some substance in space and time . Assume that the probability that a particle jumps from a position to a different position by time unit is given by a kernel . Assume moreover that the kernel is symmetric, i.e. .
Now, we further assume that the particles inside do not interact with the particles in and inside the density is zero (every particle that jumps into is automatically killed). Then, the conservation law for the mass gives rise to the equation
[TABLE]
In this paper we consider the case where the kernel behaves like a power of the distance, that is
[TABLE]
which means that the particle tends with high probability to stay close to the original position but can jump with positive probability to positions far away. This type of kernel give rise to the so-called fractional Laplace operator defined as in the introduction,
[TABLE]
where p.v. stands for principal value.
In some models, the kernel depends not only on the distance but also on the difference in concentration, making more likely to jump when the difference in concentration is large. i.e.
[TABLE]
Here we consider the case where
[TABLE]
and these type of kernels give rise to the so-called fractional laplace operator defined as
[TABLE]
One key factor to analyze the behavior of the solutions of (2.1) is the first eigenvalue of the associated operator, i.e. the smallest value of such that there exists a nontrivial solution to
[TABLE]
so the purpose of this paper is to analyze this problem and in particular the optimization of this first eigenvalue with respect to the region and the asymptotic behavior of this optimal eigenvalue when the fractional parameter tends to 1.
2.2. The fractional laplacian
In this subsection, we recall some basic facts about the fractional laplacian given in (2.2).
Lemma 2.1** ([7], Lemma 2.2).**
Let be fixed and let be open. For every , the fractional laplacian given by (2.2) defines a distribution . Moreover,
[TABLE]
for every , where is defined in (1.1).
Remark 2.2*.*
From Lemma 2.1 one observe that the fractional laplacian is a bounded operator between and its dual .
Remark 2.3*.*
The choice of the constant in the definition of the fractional laplacian is made in order for the operators to converge to the local laplacian as . In fact, it is easy to see, from the Gamma convergence results in [10] that for any , if is the weak solution of
[TABLE]
then as strongly in where is the weak solution to
[TABLE]
With all of these preliminaries, we establish the definition of weak solution for mixed boundary value problem for the fractional laplacian.
Definition 2.4**.**
Let be fixed and let be open and let . Given (or more generally, ), we say that is a weak solution of
[TABLE]
if the equality holds in the distributional sense. That is, if
[TABLE]
for every , where is given in (1.1).
The next Poincaré-type inequality, although its proof is elementary, is crucial in the remaining of the paper. This inequality is stablished in [7, Proposition 2.10]. We include here a proof for the reader’s convenience.
Proposition 2.5**.**
Let be a bounded, open set and let be a measurable bounded set. Then
[TABLE]
for every , where
[TABLE]
Proof.
The proof is rather simple. In fact, given we have
[TABLE]
Now, just observe that for every , one has
[TABLE]
and the proof is complete. ∎
We now want to remove the hypothesis that is bounded in Proposition 2.5.
Corollary 2.6**.**
Assume that is only measurable, then there is a positive constant such that
[TABLE]
Proof.
Take large so that has positive measure. Hence
[TABLE]
where is obtained from Proposition 2.5. The proof is complete. ∎
Remark 2.7*.*
if , one can take in the above proof such that . Hence one arrives at
[TABLE]
where . Therefore one can take
[TABLE]
for large.
From Corollary 2.6 it is not difficult to show the existence of an extremal for the constant . The main difficulty is that, even if is smooth, since we do not want to make any regularity assumptions on we cannot assume that the injection is compact.
Theorem 2.8**.**
Let be a bounded domain with Lipschitz boundary. Then, given measurable with positive measure, there exists such that
[TABLE]
where is given in (1.1). Moreover, the extremal can be taken to be normalized in , i.e. .
Proof.
The scheme of the proof is standard. Let be a normalized minimizing sequence for , i.e.
[TABLE]
Then, is bounded and hence, up to some subsequence, we can assume that weakly in .
Since in for every implies that in , so .
Moreover, weakly in and since is Lipschitz, this implies that the injection is compact (see [6, Theorem 7.1]) and so .
Therefore
[TABLE]
The proof is complete. ∎
Remark 2.9*.*
The Lipschitz regularity on is needed in order for the compactness of the embedding to hold. See [6, Theorem 7.1]. In fact, what is needed is that be a bounded extension domain. That is the existence of a bounded extension operator . Lipschitz boundary imply that is a bounded extension domain. See [6, Theorem 5.4].
3. The maximization problem
In this section we study the problem of maximization of . That is, given we define
[TABLE]
In particular, we are interested in the asymptotic behavior of the constant when . As we will see, the behavior of when is independent of the value of the constant .
We begin with a simple lemma.
Lemma 3.1**.**
Let be measurable sets with positive measure. Then .
Proof.
Let be such that . Then, and so
[TABLE]
taking infimum on the conclusion follows. ∎
As a consequence of this Lemma we have the following upper bound for
[TABLE]
Observe that is the first eigenvalue of the fractional laplacian with Dirichlet Boundary conditions. i.e. is the first eigenvalue of the problem
[TABLE]
Problem (3.1) was studied, among others, in [2] where, based on the famous results of Bourgain-Brezis-Mironescu [1], the asymptotic behavior of the eigenvalues when is obtained. Namely,
[TABLE]
where is the first eigenvalue of the local laplacian in with homogeneous Dirichlet boundary conditions, i.e.
[TABLE]
On the other hand, take where is the usual fattening of and is chosen in such a way as .
Then
[TABLE]
and so
[TABLE]
Now, Let be the eigenfunction associated to normalized in . Then, by [1, Theorem 4] it follows that there exist a sequence and a function , such that in and
[TABLE]
where is given in (1.1). Hence
[TABLE]
with . Therefore, is the normalized eigenfunction of the laplacian in and one can conclude that
[TABLE]
In order to finish the study of the maximization problem, we need to study the quasi-optimal Dirichlet sets . Notice that, since we are considering a maximization problem, it is not clear that there is an optimal set for . Hence, we deal with quasi-optimal sequences as , that is, we take such that
[TABLE]
where as . Our next result says that the quasi-optimal sets “covers” the boundary of as .
Theorem 3.2**.**
Given any quasi-optimal configuration for , any and , there exists such that for every it holds
[TABLE]
Proof.
Assume that the conclusion is false. Then, there exists a sequence , sets and such that
[TABLE]
Now, by (3.2), we have
[TABLE]
On the other hand, let and consider the following eigenvalue problem (that was already mentioned in the introduction)
[TABLE]
Let us denote by the first eigenvalue of (3.4), i.e.
[TABLE]
where the infimum is taken over all functions such that on .
It is straightforward to see that
[TABLE]
and that the above inequalities are strict when .
Now, denote by the eigenfunction of (3.4) associated to extended by 0 on normalized in . Then and so
[TABLE]
This contradicts (3.3) and the proof is complete. ∎
4. The minimization problem
Now we analize the minimization problem
[TABLE]
As we mentioned in the introduction, this problem is of little interest since taking a sequence of domains to infinity makes the eigenvalues go to zero.
Proposition 4.1**.**
Let be bounded and open and let be fixed. Then
[TABLE]
for every , where is the constant defined in (4.1).
Proof.
Take such that and define where so .
Now, let be defined as
[TABLE]
Observe that for large and .
So,
[TABLE]
with
[TABLE]
Now, just observe that for and and so there exists a constant independent of such that
[TABLE]
This completes the proof. ∎
We consider now the problem defined by
[TABLE]
In this case the admissible sets are forced to be inside of the ball . Note that in this case we have . Taking the limit as in we obtain that this quantity tends to zero. This is contained in the following proposition.
Proposition 4.2**.**
Let be bounded and open, let be fixed. Then
[TABLE]
for every , where is the constant defined in (4.2).
Proof.
Let be such that . Now we define the function
[TABLE]
so and .
We can then estimate as follows,
[TABLE]
On the other side
[TABLE]
then
[TABLE]
where the constant does not depend on . This implies the desired result. ∎
To end this section, we show that for any fixed there exists an optimal configuration for the constant . Let us observe that this does not follows by a direct application of the Direct Method of the Calculus of Variations since for a minimizing sequence of domains, we do not have enough compactness in a topology of domains and, what is more important, the associated eigenfunctions do not lie in the same functional space.
We begin with the following Lemma.
Lemma 4.3**.**
Given , there exists , such that is a solution to the following minimization problem
[TABLE]
Proof.
The proof is a direct consequence of the Bathtube Principle [9, Theorem 1.14]. In fact, let
[TABLE]
and observe that the level sets of , , has finite measure. So, the Bathtube principle says that the problem
[TABLE]
has a solution of the form , where , and for some .
This completes the proof of the Lemma. ∎
We can now prove the existence of the optimal configuration for .
Theorem 4.4**.**
Let be a bounded domain, let and let be such that . Then, there exists an optimal set , that is, such that and
[TABLE]
Proof.
Let be a minimizing sequence for , that is
[TABLE]
Let be the normalized, nonnegative eigenfunction associated to . Then, it is easy to see that the sequence is bounded in and so, there exists such that, up to a subsequence,
[TABLE]
Now, let be the optimal set given in Lemma 4.3. Let us check that this set is optimal for .
In fact, passing to a further subsequence, if necessary, we can assume that there exists , such that
[TABLE]
Moreover, this function satisfies
[TABLE]
Observe now, that given , the function belongs to . So, for any , we have that
[TABLE]
We can then apply Fatou’s Lemma to conclude that
[TABLE]
and so, from our choice of (c.f. Lemma 4.3), we obtain
[TABLE]
Also, by the weak semicontinuity of the Gagliardo seminorm, we have
[TABLE]
Observe now that if we extend to by zero, we have that
[TABLE]
and so (4.3) and (4.4) imply that and that
[TABLE]
The proof is complete. ∎
5. Proof of Theorem 1.2
The results for the maximization problem follows exactly as those in Section 3 with the obvious modifications.
It remains to check the results about the minimization problem.
For this, we first observe that, by Lemma 3.1 (with the obvious modifications to include the potential function ), it holds that
[TABLE]
for every , from where it follows that
[TABLE]
for any .
Recall that is the first eigenvalue of the fractional laplacian plus a potential with homogeneous Neumann boundary conditions and is given by
[TABLE]
It is not difficult to check (c.f. with Remark 2.3) that converges to the first eigenvalue of the local laplacian plus a potential with homogeneous Neumann boundary conditions. i.e.
[TABLE]
where
[TABLE]
Let us now give an estimate for . This estimate uses the ideas from Proposition 4.2.
Let be such that . So .
Denote by be the eigenfunction associated to normalized such that and extend to by zero. Since we have that and hence
[TABLE]
Now
[TABLE]
Using the results of [1] it follows that
[TABLE]
Finally
[TABLE]
(recall that , from (1.1)). Combining all this we arrive at
[TABLE]
Putting together (5.1), (5.2) and (5.3) we arrive at
[TABLE]
To finish this section just observe that the existence of an optimal configuration for the constant follows without change as in the proof of Theorem 4.4.
Acknowledgements
This paper was partially supported by Universidad de Buenos Aires under grant UBACyT 20020130100283BA, by ANPCyT under grant PICT 2012-0153 and by CONICET under grant PIP2015 11220150100032CO. J. Fernández Bonder and J. D. Rossi are members of CONICET.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Jean Bourgain, Haïm Brezis, and Petru Mironescu. Another look at Sobolev spaces, 2001. Original research article appeared at in Optimal Control and Partial Differential Equations IOS Press ISBN 1 58603 096 5.
- 2[2] Lorenzo Brasco, Enea Parini, and Marco Squassina. Stability of variational eigenvalues for the fractional p 𝑝 p -Laplacian. Discrete Contin. Dyn. Syst. , 36(4):1813–1845, 2016.
- 3[3] Dorin Bucur and Giuseppe Buttazzo. Variational methods in some shape optimization problems . Appunti dei Corsi Tenuti da Docenti della Scuola. [Notes of Courses Given by Teachers at the School]. Scuola Normale Superiore, Pisa, 2002.
- 4[4] Dorin Bucur and Paola Trebeschi. Shape optimisation problems governed by nonlinear state equations. Proc. Roy. Soc. Edinburgh Sect. A , 128(5):945–963, 1998.
- 5[5] S. Chanillo, D. Grieser, M. Imai, K. Kurata, and I. Ohnishi. Symmetry breaking and other phenomena in the optimization of eigenvalues for composite membranes. Comm. Math. Phys. , 214(2):315–337, 2000.
- 6[6] Eleonora Di Nezza, Giampiero Palatucci, and Enrico Valdinoci. Hitchhiker’s guide to the fractional Sobolev spaces. Bull. Sci. Math. , 136(5):521–573, 2012.
- 7[7] J. Fernandez Bonder and J. Spedaletti. Some nonlocal optimal design problems. Ar Xiv e-prints , January 2016.
- 8[8] Antoine Henrot and Edouard Oudet. Minimizing the second eigenvalue of the Laplace operator with Dirichlet boundary conditions. Arch. Ration. Mech. Anal. , 169(1):73–87, 2003.
