Constrained Shape Derivative for the Spectral Laplace-Dirichlet Problem Using the Minimax Method

Mame Gor Ngom¹, Bakary Kourouma¹, Ibrahima Faye²

¹Département de Mathematiques, Université Gamal Abdel Nasser de Conakry, Conakry, Guinea.
²Laboratoire d’Informatique, de Mathématiques et Applications (LIMA), U.F.R-SATIC, Université Alioune Diop, Bambey, Senegal.

DOI: 10.4236/jamp.2025.137136                

Abstract

In shape optimization, we often seek to optimize the shape of a domain in order to achieve a goal. For example, the sound of a drum depends on the shape of its membrane, via the eigenvalues of the Laplacian. We use the shape derivative to adjust this shape. The shape derivative for an eigenvalue problem is frequently used in shape optimisation, in particular when we want to minimise or maximise an eigenvalue, often the first of an elliptic operator, typically the Laplacian, under geometric constraints. This paper gives a simple method for calculating the shape derivative of a Dirichlet eigenvalue, but also the shape derivative of an objective function constrained to an eigenvalue problem.

Keywords

Shape Optimization, Shape Derivative, Eigenvalue Problem

1. Introduction

In modern applied mathematics, more and more problems involve not just finding a function that satisfies an equation, but determining the very shape of the domain in which that function lives. These shapes are no longer fixed: they become variables in their own right, malleable and deformable, which we seek to adjust to achieve a given objective. Imagine a musical instrument, a drum for example. Its timbre depends closely on the shape of its membrane. Altering its contour, even slightly, changes its sound. Behind this acoustic intuition lies a profound mathematical truth: the eigenvalues of the Laplacian, which govern the frequencies of vibration, depend on the geometry of the domain. And if we want to optimise this geometry to obtain a deeper, higher or purer sound, we need to know how these eigenvalues change when we modify the shape of the domain. This is where the shape derivative comes in. There are several works in the literature on the shape derivative of the first eigenvalue. These include [1]-[6]. Here, we use other techniques to calculate the shape derivative of the first eigenvalue for the Laplace-Dirichlet problem. Next, we calculate the shape derivative of an objective function, taking the spectral problem as a constraint. Readers interested in the approach developed can consult the work of M. C. Delfour [7] [8] or our own [9] [10]. The shape derivative is a mathematical tool used to study the infinitesimal variation of a functional $J\left(\Omega \right)$ when the domain is slightly deformed. This deformation is generally given by a vector field which defines a family of domains ${\Omega }_{ϵ}=\left(I+ϵV\right)\left(\Omega \right)$ , for $ϵ$ small or more generally ${\Omega }_{ϵ}=T_{ϵ}\left(\Omega \right)$ where the family of diffeomorphisms $\left(T_{ϵ}\right)$ is induced by the vector field $V$ . The shape derivative of $J\left(\Omega \right)$ is then defined by

$DJ\left(\Omega ,V\right):={\frac{\text{d}}{\text{d}ϵ}J\left({\Omega }_{ϵ}\right)|}_{ϵ=0}=\underset{ϵ\to 0}{\lim }\frac{J\left({\Omega }_{ϵ}\right)-J\left(\Omega \right)}{ϵ}.$

This is the analogue of the directional derivative, but in the space of shapes (domains). The aim of this work is to calculate, using the min-max method [7], the shape derivative of the objective function $J\left(\Omega \right)$ defined by

$J\left(\Omega \right)={\displaystyle {\int }_{\Omega }{\left|\nabla \eta -A\right|}^2\text{d}x}+{\displaystyle {\int }_{\Omega }{\left|\eta -{\eta }_d\right|}^2\text{d}x}$ (1.1)

where $\eta$ is the solution of the following eigenvalue problem

$\{\begin{array}{l}-\text{Δ}\eta -\lambda \eta =0\text{in}\Omega \\ \eta =0\text{on}\partial \Omega ,\end{array}$ (1.2)

where $\Omega$ is a regular open domain of ${ℝ}^N$ , $A$ a given vector of ${ℝ}^N$ and ${\eta }_d$ a given function. The existence and uniqueness of this problem are detailed in [11]. Note that the problem under study has no source term. In other words, $\lambda =\lambda \left(\Omega \right)$ is considered an eigenvalue and therefore depends on the domain $\Omega$ . The first step is to determine the derivative of this eigenvalue with respect to the domain. The rest of the manuscript is organized as follows: Section 2 presents a reminder of the methodology used to calculate the derivative. Section 3 is devoted to the study of the shape derivative of the volume function. Numerical simulations are proposed in Section 4, while the final section is devoted to results relating to the shape derivative in the context of the eigenvalue problem.

2. Methodology for Calculating Shape Derivatives

In this subsection, we describe how to calculate the topological derivative using the min-max approach, see e.g. [7] [9]. To begin with, we will look at the following definitions and notations.

Definition 2.1 A Lagrangian function is a function of the form

$\left(ϵ,x,y\right)\mapsto ℒ\left(ϵ,x,y\right):\left[0,\tau \right]\times X\times Y\to ℝ,\tau >0$

where X is a vector espace, Y a non empty subset of vector space and the function $y\mapsto ℒ\left(ϵ,x,y\right)$ is affine.

Associate with the parameter $ϵ$ the parametrized minimax

$ϵ\mapsto g\left(ϵ\right)=\underset{x\in X}{\text{inf}}\underset{y\in Y}{\text{sup}}ℒ\left(ϵ,x,y\right):\left[0,\tau \right]\to ℝ\text{and}dg\left(0\right)=\underset{ϵ\to 0^{+}}{\text{lim}}\frac{g\left(ϵ\right)-g\left(0\right)}{ϵ}.$

When the limits exist, we will use the following notations

$d_{ϵ}ℒ\left(0,x,y\right)=\underset{ϵ\to 0^{+}}{\text{lim}}\frac{ℒ\left(ϵ,x,y\right)-ℒ\left(0,x,y\right)}{ϵ}$

$\phi \in X,d_xℒ\left(ϵ,x,y;\phi \right)=\underset{\theta \to 0^{+}}{\text{lim}}\frac{ℒ\left(ϵ,x+\theta \phi ,y\right)-ℒ\left(ϵ,x,y\right)}{\theta }$

$\varphi \in Y,d_yℒ\left(ϵ,x,y;\varphi \right)=\underset{\theta \to 0^{+}}{\text{lim}}\frac{ℒ\left(ϵ,x,y+\theta \varphi \right)-ℒ\left(ϵ,x,y\right)}{\theta }.$

The state equation at $ϵ\ge 0$

$\text{Find}{x}^{ϵ}\in X\text{such}\text{that}\text{for}\text{all}\psi \in Y,d_yℒ\left(ϵ,x^{ϵ},0;\psi \right)=0.$ (2.1)

The set of states $x^{ϵ}$ at $ϵ\ge 0$ is denoted

$E\left(ϵ\right)=\left\{x^{ϵ}\in X,\forall \psi \in Y,d_yℒ\left(ϵ,x^{ϵ},0;\psi \right)=0\right\}.$ (2.2)

The adjoint equation at $ϵ\ge 0$ is

$\text{Find}{p}^{ϵ}\in Y\text{such}\text{that}\text{for}\text{all}\phi \in X,d_xℒ\left(ϵ,x^{ϵ},p^{ϵ};\phi \right)=0.$ (2.3)

The set of solutions $p^{ϵ}$ at $ϵ\ge 0$ is denoted

$Y\left(ϵ,x^{ϵ}\right)=\left\{p^{ϵ}\in Y,\forall \phi \in X,d_xℒ\left(ϵ,x^{ϵ},p^{ϵ};\phi \right)=0\right\}.$ (2.4)

Finally the set of minimisers for the minimax is given by

$X\left(ϵ\right)=\left\{x^{ϵ}\in X,g\left(ϵ\right)=\underset{x\in X}{\text{inf}}\underset{y\in Y}{\text{sup}}ℒ\left(ϵ,x,y\right)=\underset{y\in Y}{\text{sup}}ℒ\left(ϵ,x^{ϵ},y\right)\right\}.$ (2.5)

Lemma 2.1 (Constrained infimum and minimax)

We have the following assertions

1) $\underset{x\in X}{\text{inf}}\underset{y\in Y}{\text{sup}}ℒ\left(ϵ,x,y\right)=\underset{x\in E\left(ϵ\right)}{\text{inf}}ℒ\left(ϵ,x,0\right)$

2) The minimax $g\left(ϵ\right)=+\infty$ if and only if $E\left(ϵ\right)=\varnothing$ . And in this case we have $X\left(ϵ\right)=X$ .

3) If $E\left(ϵ\right)\ne \varnothing$ , then

$X\left(ϵ\right)=\left\{x^{ϵ}\in E\left(ϵ\right):ℒ\left(ϵ,x^{ϵ},0\right)=\underset{x\in E\left(ϵ\right)}{\text{inf}}ℒ\left(ϵ,x,0\right)\right\}\subset E\left(ϵ\right)$

and $g\left(ϵ\right)<+\infty$ .

Proof. See [7]. ■

First Results of the Proposed Approach

We need the following assumption for everything that follows:

Hypothesis (H0)

Let X be a vector space.

1) For all $ϵ\in \left[0,\tau \right]$ , $x^0\in X\left(0\right)$ , $x^{ϵ}\in X\left(ϵ\right)$ and $y\in Y$ , the function $\theta \mapsto ℒ\left(ϵ,x^0+\theta \left(x^{ϵ}-x^0\right),y\right):\left[0,1\right]\to ℝ$ is absolutely continuous. This implies that for almost all $\theta$ the derivative exists and is equal to $d_xℒ\left(ℒ,x^0+\theta \left(x^{ϵ}-x^0\right),y;x^{ϵ}-x^0\right)$ and it is the integral of its derivative. In particular

$ℒ\left(ϵ,x^{ϵ},y\right)=ℒ\left(ϵ,x^0,y\right)+{\displaystyle {\int }_0^1{d}_xℒ\left(ϵ,x^0+\theta \left(x^{ϵ}-x^0\right),y;x^{ϵ}-x^0\right)\text{d}\theta }.$

2) For all $ϵ\in \left[0,\tau \right]$ , $x^0\in X\left(0\right)$ , $x^{ϵ}\in X\left(ϵ\right)$ and $y\in Y$ , $\varphi \in X$ and for almost all $\theta \in \left[0,1\right]$ , $d_xℒ\left(ϵ,x^0+\theta \left(x^{ϵ}-x^0\right),y;\varphi \right)$ exist et the functions $\theta \mapsto d_xℒ\left(ϵ,x^0+\theta \left(x^{ϵ}-x^0\right),y;\varphi \right)$ belong to $L^1\left[0,1\right]$

Definition 2.2 Given $x^0\in X\left(0\right)$ and $x^{ϵ}\in X\left(ϵ\right)$ , the averaged adjoint equation is:

$\text{Find}{y}^{ϵ}\in Y,\forall \varphi \in X,{\displaystyle {\int }_0^1{d}_xℒ\left(ϵ,x^0+\theta \left(x^{ϵ}-x^0\right),y;\varphi \right)\text{d}\theta }=0$

and the set of solutions is noted $Y\left(ϵ,x^0,x^{ϵ}\right)$ .

$Y\left(0,x^0,x^0\right)$ clearly reduces to the set of standard adjoint states $Y\left(0,x^0\right)$ at $ϵ=0$ .

Theorem 2.2 Consider the Lagrangian functional

$\left(ϵ,x,y\right)\mapsto ℒ\left(ϵ,x,y\right):\left[0,\tau \right]\times X\times Y\to ℝ,\tau >0$

where X and Y are vector spaces and the function $y\mapsto ℒ\left(ϵ,x,y\right)$ is affine. Assume that (H0) and the following hypotheses are satisfied

(H1) for all $ϵ\in \left[0,\tau \right]$ , $g\left(ϵ\right)$ is finite, $X\left(ϵ\right)=\left\{x^{ϵ}\right\}$ and $Y\left(0,x^0\right)=\left\{p^0\right\}$ are singletons,

(H2) $d_{ϵ}ℒ\left(0,x^0,y^0\right)$ exists,

(H3) The following limit exists

$ℛ\left(x^0,y^0\right)=\underset{ϵ\to 0^{+}}{\text{lim}}{\displaystyle {\int }_0^1{d}_xℒ}\left(ϵ,x^0+\theta \left(x^{ϵ}-x^0\right),p^0;\frac{x^{ϵ}-x^0}{ϵ}\right)\text{d}\theta .$

Then, $dg\left(0\right)$ exists and $dg\left(0\right)=d_{ϵ}ℒ\left(0,x^0,p^0\right)+ℛ\left(x^0,p^0\right)$ .

Proof. Recall that $g\left(ϵ\right)=ℒ\left(ϵ,x^{ϵ},y\right)$ and $g\left(0\right)=ℒ\left(0,x^0,y\right)$ for each $y\in Y$ , then for a standard adjoint state $p^0$ at $ϵ=0$

$g\left(ϵ\right)-g\left(0\right)=ℒ\left(ϵ,x^{ϵ},p^0\right)-ℒ\left(ϵ,x^0,p^0\right)+\left(ℒ\left(ϵ,x^0,p^0\right)-ℒ\left(0,x^0,p^0\right)\right).$

Dividing by $ϵ>0$ , we obtain

$\begin{array}{c}\frac{g\left(ϵ\right)-g\left(0\right)}{ϵ}=\frac{ℒ\left(ϵ,x^{ϵ},p^0\right)-ℒ\left(ϵ,x^0,p^0\right)}{ϵ}+\frac{ℒ\left(ϵ,x^0,p^0\right)-ℒ\left(0,x^0,p^0\right)}{ϵ}\\ ={\displaystyle {\int }_0^1{d}_xℒ}\left(ϵ,x^0+\theta \left(x^{ϵ}-x^0\right),p^0;\frac{x^{ϵ}-x^0}{ϵ}\right)\text{d}\theta \\ +\frac{ℒ\left(ϵ,x^0,p^0\right)-ℒ\left(0,x^0,p^0\right)}{ϵ}.\end{array}$

Going to the limit when $ϵ$ goes to zero, we obtain

$\begin{array}{c}dg\left(0\right)=\underset{ϵ\to 0^{+}}{\text{lim}}{\displaystyle {\int }_0^1{d}_xℒ}\left(ϵ,x^0+\theta \left(x^{ϵ}-x^0\right),p^0;\frac{x^{ϵ}-x^0}{ϵ}\right)\text{d}\theta +d_{ϵ}ℒ\left(0,x^0,p^0\right)\\ =d_{ϵ}ℒ\left(0,x^0,p^0\right)+ℛ\left(x^0,p^0\right).\end{array}$

■

Corollary 2.3 Consider the Lagrangian functional

$\left(ϵ,x,y\right)\mapsto ℒ\left(ϵ,x,y\right):\left[0,\tau \right]\times X\times Y\to ℝ,\tau >0$

where X and Y are vector spaces and the function $y\mapsto ℒ\left(ϵ,x,y\right)$ is affine. Assume that (H0) and the following assumptions are satisfied:

(H1a) for all $ϵ\in \left[0,\tau \right]$ , $X\left(ϵ\right)\ne \varnothing$ , $g\left(ϵ\right)$ is finite, and for each $x\in X\left(0\right)$ , $Y\left(0,x\right)\ne \varnothing$ ,

(H2a) for all $x\in X\left(0\right)$ and $p\in Y\left(0,x\right)$ $d_{ϵ}ℒ\left(0,x,p\right)$ exists,

(H3a) there exist $x^0\in X\left(0\right)$ and $p^0\in Y\left(0,x^0\right)$ such that the following limit exists

$ℛ\left(x^0,p^0\right)=\underset{ϵ\to 0^{+}}{\text{lim}}{\displaystyle {\int }_0^1{d}_xℒ}\left(ϵ,x^0+\theta \left(x^{ϵ}-x^0\right),p^0;\frac{x^{ϵ}-x^0}{ϵ}\right)\text{d}\theta .$

Then, $dg\left(0\right)$ exists and there exist $x^0\in X\left(0\right)$ and $p^0\in Y\left(0,x^0\right)$ such that $dg\left(0\right)=d_{ϵ}ℒ\left(0,x^0,p^0\right)+ℛ\left(x^0,p^0\right)$ .

Proposition 2.4 Let X and Y be two vector spaces and $ℒ:\left[0,\tau \right]\times X\times Y\to ℝ$ be a Lagrangian. If the function $y\to ℒ\left(ϵ,x,y\right)$ is affine, then for all $y,\psi \in Y$ , we have

$d_yℒ\left(ϵ,x,y;\psi \right)=ℒ\left(ϵ,x,\psi \right)-ℒ\left(ϵ,x,0\right)=d_yℒ\left(ϵ,x,0;\psi \right).$

Proof. The Lagrangian $ℒ$ is affine in $y$ means

$ℒ\left(ϵ,x,y\right)={ℒ}_1\left(ϵ,x\right)+{ℒ}_2\left(ϵ,x,y\right)$ , où ${ℒ}_2$ is linear in $y$ .

Since ${ℒ}_1$ does not depend on $y$ , we have

$d_yℒ\left(ϵ,x,y;\psi \right)=d_y{ℒ}_2\left(ϵ,x,y;\psi \right).$

By definition, the following applies

$\begin{array}{l}{d}_y{ℒ}_2\left(ϵ,x,y;\psi \right)=\underset{\theta \to 0}{\text{lim}}\frac{{ℒ}_2\left(ϵ,x,y+\theta \psi \right)-{ℒ}_2\left(ϵ,x,y\right)}{\theta }\\ =\underset{\theta \to 0}{\text{lim}}\frac{{ℒ}_2\left(ϵ,x,y\right)+\theta {ℒ}_2\left(ϵ,x,\psi \right)-{ℒ}_2\left(ϵ,x,y\right)}{\theta },\text{because}{ℒ}_2\text{is}\text{linear}\text{in}y\\ =\underset{\theta \to 0}{\text{lim}}\frac{\theta {ℒ}_2\left(ϵ,x,\psi \right)}{\theta }={ℒ}_2\left(ϵ,x,\psi \right).\end{array}$

So for all $y\in Y$

$d_yℒ\left(ϵ,x,y;\psi \right)=d_y{ℒ}_2\left(ϵ,x,y;\psi \right)={ℒ}_2\left(ϵ,x,\psi \right),$ (2.6)

and it follows that

$d_yℒ\left(ϵ,x,0;\psi \right)=d_y{ℒ}_2\left(ϵ,x,y;\psi \right)={ℒ}_2\left(ϵ,x,\psi \right).$ (2.7)

Now let us calculate $ℒ\left(ϵ,x,\psi \right)-ℒ\left(ϵ,x,0\right)$ .

$\begin{array}{c}ℒ\left(ϵ,x,\psi \right)-ℒ\left(ϵ,x,0\right)={ℒ}_1\left(ϵ,x\right)+{ℒ}_2\left(ϵ,x,\psi \right)-{ℒ}_1\left(ϵ,x\right)-{ℒ}_2\left(ϵ,x,0\right)\\ ={ℒ}_2\left(ϵ,x,\psi \right)-{ℒ}_2\left(ϵ,x,0\right).\end{array}$

But ${ℒ}_2\left(ϵ,x,0\right)=0$ because ${ℒ}_2$ is linear in $y$ and it follows

$ℒ\left(ϵ,x,\psi \right)-ℒ\left(ϵ,x,0\right)={ℒ}_2\left(ϵ,x,\psi \right).$ (2.8)

So the relationships (2.6), (2.7) and (2.8) give the result. ■

3. Shape Derivative for the Volume Function via Min-Max Approach

The following example is a special case where we consider a shape functional independent of a PDE constraint. It is given in [5] by

$J\left(\Omega \right)=\left|\Omega \right|\lambda \left(\Omega \right)$ (3.1)

$\left|\Omega \right|$ is the volume of $\Omega$ et $\lambda \left(\Omega \right)=\lambda$ is the first eigenvalue of the following eigenvalue problem

$-\text{Δ}\eta =\lambda \eta \text{in}\Omega \text{and}\eta =0\text{on}\partial \Omega .$ (3.2)

This problem has a sequence of real positive eigenvalues:

$0<{\lambda }_1\left(\Omega \right)<{\lambda }_2\left(\Omega \right)\le \cdots \to \infty ,$

with associated eigenfunctions ${\eta }_i\in H_0^1\left(\Omega \right)$ , orthonormal in $L^2\left(\Omega \right)$ , i.e.

${\displaystyle {\int }_{\Omega }{\eta }_i{\eta }_j\text{d}x}={\delta }_{ij}.$

In this case, we need the volume expression of the derivative of the eigenvalue. For this we have the following proposition.

Proposition 3.1 Let $\Omega$ be a bounded open domain. Assume that the first eigenvalue $\lambda =\lambda \left(\Omega \right)$ is simple. Then $\lambda \left(\Omega \right)$ is shape differentiable and the Eulerian derivative is

$\begin{array}{l}D\lambda \left(\text{Ω},V\right)={\lambda }^{\prime }\left(0\right)\\ ={\displaystyle {\int }_{\Omega }\left(-2\nabla {\eta }^0\cdot DV\left(0\right){\left|\nabla {\eta }^0\right|}^2+divV\left(0\right)\left({\left|\nabla {\eta }^0\right|}^2-{\lambda }^0{\left|{\eta }^0\right|}^2\right)\right)\text{d}x}.\end{array}$ (3.3)

If, in addition, $\Omega$ is convex or of class $C^2$ , then

$D\lambda \left(\Omega ,V\right)={\lambda }^{\prime }\left(0\right)=-{\displaystyle {\int }_{\partial \Omega }{\left(\frac{\partial \eta }{\partial n}\right)}^{2}V\cdot n\text{d}\sigma }.$ (3.4)

Proof. A proof of this proposition can also be found in the work of A. Henrot and M. Pierre [4] or that of S. Zhu [5]. We only prove the derivative formula in volume form (3.3). The reader can also consult [6]. In this proof, we will therefore use a rigorous demonstration based on the min-max method introduced by M. C. Delfour in [7]. To do this, we consider the diffeomorphism family ${\left(T_{ϵ}\right)}_{ϵ\ge 0}$ and define the perturbed domain by ${\Omega }_{ϵ}=T_{ϵ}\left(\Omega \right)$ . The perturbed state ${\eta }_{ϵ}$ is the solution of the variational problem

${\displaystyle {\int }_{{\Omega }_{ϵ}}\nabla {\eta }_{ϵ}\nabla v^{\prime }\text{d}x}={\displaystyle {\int }_{{\Omega }_{ϵ}}{\lambda }_{ϵ}{\eta }_{ϵ}{v}^{\prime }\text{d}x},\forall v^{\prime }\in H_0^1\left({\Omega }_{ϵ}\right).$

A change of variables allows the problem to be transported to the reference domain $\Omega$ . Thus, noting by ${\eta }^{ϵ}={\eta }_{ϵ}\circ T_{ϵ}$ and ${\lambda }^{ϵ}={\lambda }_{ϵ}\circ T_{ϵ}$ , ${\eta }^{ϵ}$ is the solution of the following variational equation

${\displaystyle {\int }_{\Omega }B\left(ϵ\right)\nabla {\eta }^{ϵ}\nabla v\text{d}x}={\displaystyle {\int }_{\Omega }{\lambda }^{ϵ}{\eta }^{ϵ}vJ\left(ϵ\right)\text{d}x}.$ (3.5)

For $\phi ={\eta }^{ϵ}$

${\displaystyle {\int }_{\Omega }B\left(ϵ\right)\nabla \phi \nabla v\text{d}x}={\displaystyle {\int }_{\Omega }{\lambda }^{ϵ}\phi vJ\left(ϵ\right)\text{d}x}.$ (3.6)

Deriving with respect to $ϵ$ and taking $v=\phi$ , we obtain

${\displaystyle {\int }_{\Omega }{B}^{\prime }\left(ϵ\right){\left|\nabla \phi \right|}^2\text{d}x}={\displaystyle {\int }_{\Omega }{\left({\lambda }^{ϵ}\right)}^{\prime }{\left|\phi \right|}^{2}J\left(ϵ\right)\text{d}x}+{\displaystyle {\int }_{\Omega }{\left({\lambda }^{ϵ}\right)}^{\prime }{\left|\phi \right|}^{2}divV\left(ϵ\right)J\left(ϵ\right)\text{d}x}.$ (3.7)

So at $ϵ=0$ with $\phi ={\eta }^0$ .

${\displaystyle {\int }_{\Omega }{B}^{\prime }\left(0\right){\left|\nabla {\eta }^0\right|}^2\text{d}x}={\displaystyle {\int }_{\Omega }{\left({\lambda }^{ϵ}\right)}^{\prime }\left(0\right){\left|{\eta }^0\right|}^2\text{d}x}+{\displaystyle {\int }_{\Omega }{\lambda }^0}{\left|{\eta }^0\right|}^{2}divV\left(0\right)\text{d}x.$ (3.8)

Using the normalization condition, we have

${\left({\lambda }^{ϵ}\right)}^{\prime }\left(0\right)={\displaystyle {\int }_{\Omega }{B}^{\prime }\left(0\right){\left|\nabla {\eta }^0\right|}^2\text{d}x}-{\displaystyle {\int }_{\Omega }{\lambda }^0{\left|{\eta }^0\right|}^{2}divV\left(0\right)\text{d}x}$ (3.9)

with $B^{\prime }\left(0\right)=divV\left(0\right)-DV{\left(0\right)}^{\text{*}}-DV\left(0\right)$ . Using the fact that

$\left(DV{\left(0\right)}^{\text{*}}+DV\left(0\right)\right){\left|\nabla {\eta }^0\right|}^2=2\nabla {\eta }^0\cdot DV{\left|\nabla {\eta }^0\right|}^2$

${\left({\lambda }^{ϵ}\right)}^{\prime }\left(0\right)={\displaystyle {\int }_{\Omega }\left(-2\nabla {\eta }^0\cdot DV\left(0\right){\left|\nabla {\eta }^0\right|}^2+divV\left(0\right)\left({\left|\nabla {\eta }^0\right|}^2-{\lambda }^0{\left|{\eta }^0\right|}^2\right)\right)\text{d}x}.$ (3.10)

■

Proposition 3.2 Let $J$ be the objective function defined by $J\left(\Omega \right)=\left|\Omega \right|\lambda \left(\Omega \right)$ . Then the derivative of $J$ in $\Omega$ in the direction $V$ is given by

$\begin{array}{l}DJ\left(\Omega ,V\right)\\ ={\displaystyle {\int }_{\Omega }\left[\left|\Omega \right|\left(-2\nabla {\eta }^0\cdot DV{\left|\nabla {\eta }^0\right|}^2+divV\left({\left|\nabla {\eta }^0\right|}^2-{\lambda }^0{\left|{\eta }^0\right|}^2\right)\right)+{\lambda }^{0}divV\right]\text{d}x}.\end{array}$ (3.11)

Proof. If we consider the objective function

$J\left(\Omega \right)=\left|\Omega \right|\lambda \left(\Omega \right),$ (3.12)

first, we know that the derivative of the volume functional is given by

$\left|\Omega \right|\text{'}\left(\Omega ,V\left(0\right)\right)={\displaystyle {\int }_{\Omega }divV\left(0\right)\text{d}x}.$

Indeed, by definition of the shape derivative

$\left|\Omega \right|\text{'}\left(\Omega ,V\left(0\right)\right)=\underset{ϵ\to 0}{\text{lim}}\frac{\left|{\Omega }_{ϵ}\right|-\left|\Omega \right|}{ϵ}=\underset{ϵ\to 0}{\text{lim}}\frac{1}{ϵ}\left({\displaystyle {\int }_{{\Omega }_{ϵ}}1\text{d}x}-{\displaystyle {\int }_{\Omega }1\text{d}x}\right).$

We use a change of variables to transform the domain ${\Omega }_{ϵ}$ into $\Omega$

$\begin{array}{c}\left|\Omega \right|\text{'}\left(\Omega ,V\left(0\right)\right)=\underset{ϵ\to 0}{\text{lim}}{\displaystyle {\int }_{\Omega }\frac{J\left(ϵ\right)-1}{ϵ}\text{d}x}=\underset{ϵ\to 0}{\text{lim}}{\displaystyle {\int }_{\Omega }\frac{J\left(ϵ\right)-J\left(0\right)}{ϵ}\text{d}x}\\ ={\displaystyle {\int }_{\Omega }{J}^{\prime }\left(0\right)\text{d}x}={\displaystyle {\int }_{\Omega }divV\left(0\right)\circ T\left(0\right)J\left(0\right)\text{d}x}.\end{array}$

Since $T\left(0\right)=I$ and $J\left(0\right)=1$ , we have well $\left|\Omega \right|\text{'}\left(\Omega ,V\left(0\right)\right)={\displaystyle {\int }_{\Omega }divV\left(0\right)\text{d}x}$ . And finally, the derivative of the objective function $J$ is given by

$DJ\left(\Omega ,V\right)={\displaystyle {\int }_{\Omega }\left[\left|\Omega \right|\left(-2\nabla {\eta }^0\cdot DV{\left|\nabla {\eta }^0\right|}^2+divV\left({\left|\nabla {\eta }^0\right|}^2-{\lambda }^0{\left|{\eta }^0\right|}^2\right)\right)+{\lambda }^{0}divV\right]\text{d}x}.$

■

4. Numerical Tests

In shape optimization, the aim is to find the most efficient shape for a certain criterion in an admissible set of domains. In this section, we present numerical tests for minimizing the objective function defined in (3.1) taking into account that $\lambda$ is the first eigenvalue of the operator $-\text{Δ}$ .

Shape Optimization Implementation Details

The interest of the shape derivative is to give a direction of descent to evolve the domain in order to minimize the eigenvalue. The idea is to find a vector field V that gives a decrease in the objective function, i.e. we want to find a vector field V such that $DJ\left(\Omega ,V\right)\le 0$ . This can be achieved by solving an auxiliary limit problem of the type

$\text{Find}V\in \mathscr{H};b\left(V,\varphi \right)=DJ\left(\Omega ,\varphi \right)\text{for}\text{all}\varphi \in \mathscr{H}$

for a suitable Herbert space $\mathscr{H}$ . $b\left(\cdot ,\cdot \right):\mathscr{H}\times \mathscr{H}\to ℝ$ is a positive bilinear form that we are free to choose. Then $-V$ is a descent direction because

$DJ\left(\Omega ,-V\right)=-DJ\left(\Omega ,V\right)=-b\left(V,V\right)<0.$

We present a shape gradient descent algorithm, based on a volume expression rather than a boundary expression. In this case, it will be essential to update the finite element mesh after each iteration. To this end, we follow the approaches developed by V. H. Schulz, M. Siebenborn, K. Welker, [12] using the linear elasticity equation given as follows

$div\left(\sigma \right)=f^{elas}\text{in}\Omega$

$u=0\text{on}\partial \Omega$

$\sigma :=\lambda Tr\left(ϵ\right)I+2\mu ϵ$

$ϵ:=\frac{1}{2}\left(\nabla u+\nabla u^{\text{T}}\right)$

where $\sigma$ and $ϵ$ are the deformation and continuity tensors, $I$ is the identity tensor, $Tr$ is the trace operator on a tensor and $u$ is the displacement vector field [12] [13]. $\beta$ and $\mu$ designate the Lamé parameters, which can be expressed in terms of Young's modulus $E$ and Poisson’s ratio $\nu$ as follows

$\beta =\frac{\nu E}{\left(1+\nu \right)\left(1-2\nu \right)},\mu =\frac{E}{2\left(1+\nu \right)}.$

The Young’s modulus $E$ and Poisson’s ratio $\nu$ have effects on mesh deformation. $E$ indicates the stiffness of the material, allowing us to control the step size for shape updating, while $\nu$ gives the ratio between mesh expansion in other coordinate directions and compression in a particular direction. These coefficients can be found, for instance, in [12] [13]. The bilinear form of the linear elasticity equation is given by

$a\left(u,v\right)={\displaystyle {\int }_{\Omega }\sigma }\left(u\right):ϵ\left(v\right)\text{d}x.$

Thanks to the Riesz representation theorem, we obtain the direction of descent or the shape gradient or the solution of the linear elasticity equation $U:\Omega \to {ℝ}^2$ , which can be solved as follows

${\displaystyle {\int }_{\Omega }\sigma }\left(U\right):ϵ\left(V\right)\text{d}x=DJ\left(\Omega ,V\right),\forall V\in H_0^1\left(\Omega ,{ℝ}^2\right).$

The shape functional is not subject to any PDE constraints. It is therefore sufficient to solve the state equation and the deformation equation at each iteration step using the following algorithm:

Algorithm 1

Initialization while $\Vert DJ\left({\Omega }_k,V\right)\Vert >{ϵ}_{TOL}$ do 1. Calculate the solution ${\eta }_k$ of the spectral problem. 2. Calculate the gradient $W_k$ [via $DJ\left(\Omega ,V\right)$ and Linear elasticity]. 3. Deforming the domain ${\Omega }_k$ in ${\Omega }_{k+1}$ [via ALE.move and $W_k$ ]. end while

First, it is worth noting that all numerical simulations are conducted in two dimensions within the FEniCS framework, see for instance [13] to find out more. The algorithm used here is a slight modification of the one proposed in [6], in which $\lambda$ is assigned a value. In our case, since $\lambda$ is an eigenvalue, it is calculated directly. In the absence of a partial differential equation constraint, each iteration consists of solving the equation of state and the system associated with the direction of descent, or simply the linear elasticity equation. As we shall see, obtaining the new domain shape relies solely on moving the points on the boundary. We implement the general problem of minimizing the functional (3.1) depending on $\lambda$ , eigenvalue for (3.2). For each iteration, we need to solve the state equation (3.2). The shape derivative we use is given by the formula (3.11). To define the Lame parameters, we fixed the values $E=0.1$ and $\nu =0.4$ throughout the numerical tests. We set the time step to 0.001. We follow the Algorithm 1 to obtain a decrease in the functional introduced in (3.1). Tests were first carried out for 100 iterations, then extended to 500 and finally 1000 iterations. The convergence of the algorithm depends on the tolerance chosen which is set here at 0.088. The initial mesh is the unit square, as shown in Figure 1.

After some numerical testing, the vector fields look like this: Figure 2 shows the vector fields for 100, 500 and 1000 iterations.

Figure 3, Figure 4 and Figure 5 show the optimized domains, the solution of the state equation and the evolution of the objective function for 100, 500 and 1000 iterations respectively.

Figure 1. Initial mesh.

Figure 2. Vectors fields after 100, 500, and 1000 iterations.

Figure 3. (1): Optimized domain, (2): State, (3): Objective function for 100 iterations.

Figure 4. (1): Optimized domain, (2): State, (3): Objective function for 500 iterations.

Figure 5. (1): Optimized domain, (2): State, (3): Objective function for 1000 iterations.

We also carried out tests with the unit disk. Given the slow rate of convergence, the number of iterations was limited to 500. The simulation results are presented in Figure 6 below.

Figure 6. (1): Optimized domain, (2): State, (3): Objective function for 500 iterations.

5. Shape Derivative of an Eigenvalue Problem Using a Min-Max Approach

In this section, we focus on calculating the shape derivative in the presence of a spectral problem as a constraint. We use the same computational techniques developed in Section 2. The difference is that here, the domain-dependent $\lambda =\lambda \left(\Omega \right)$ is considered as the first eigenvalue of the Dirichlet Laplacian.

Let $\Omega$ be a bounded lipschitzian open domain in ${ℝ}^N$ with a regular boundary $\partial \Omega$ . Let $\eta =\eta \left(\Omega \right)$ be the solution to the following spectral problem

$-\text{Δ}\eta -\lambda \eta =0\text{in}\Omega ,\eta =0\text{on}\partial \Omega .$ (5.1)

is considered here as the first eigenvalue of the Laplacian. We consider the functional $J$ defined by

$J\left(\Omega \right)={\displaystyle {\int }_{\Omega }{\left|\nabla \eta -A\right|}^2\text{d}x}+{\displaystyle {\int }_{\Omega }{\left|\eta -{\eta }_d\right|}^2\text{d}x}$ (5.2)

where $\eta$ is the solution of (5.1), A a given vector of ${ℝ}^N$ and ${\eta }_d$ a given function. As in the case of the HeLmholtz equation [9], we want to calculate the shape derivative of the function $J$ . Recall that the shape derivative of $J$ in $\Omega$ in the direction $V$ is defined by

$DJ\left(\Omega ,V\right)=\underset{ϵ\searrow 0}{\text{lim}}\frac{J\left({\Omega }_{ϵ}\right)-J\left(\Omega \right)}{ϵ}$

when the limit exists where ${\Omega }_{ϵ}=T_{ϵ}\left(\Omega \right)$ . We will also follow M. C. Delfour [7] [8] to calculate the shape derivative of the $J$ functional. The functional in the domain ${\Omega }_{ϵ}$ is then defined by

$J\left({\Omega }_{ϵ}\right)={\displaystyle {\int }_{{\Omega }_{ϵ}}{\left|\nabla {\eta }_{ϵ}-A\right|}^2\text{d}x}+{\displaystyle {\int }_{{\Omega }_{ϵ}}{\left|{\eta }_{ϵ}-{\eta }_d\right|}^2\text{d}x}$ (5.3)

where ${\eta }_{ϵ}$ is the solution to the problem

$-\text{Δ}{\eta }_{ϵ}-{\lambda }_{ϵ}{\eta }_{ϵ}=0\text{in}{\Omega }_{ϵ},{\eta }_{ϵ}=0\text{on}\partial {\Omega }_{ϵ}.$ (5.4)

$\partial {\Omega }_{ϵ}$ is the boundary of ${\Omega }_{ϵ}$ . The variational formulation of (5.4) is to find ${\eta }_{ϵ}\in H_0^1\left({\Omega }_{ϵ}\right)$ such that

${\displaystyle {\int }_{{\Omega }_{ϵ}}\nabla {\eta }_{ϵ}\nabla v^{\prime }-{\lambda }_{ϵ}{\eta }_{ϵ}{v}^{\prime }\text{d}x}=0\text{for}\text{all}{v}^{\prime }\in H_0^1\left({\Omega }_{ϵ}\right).$ (5.5)

Note that this problem is not trivial, as the operator acts on domains that change ${\Omega }_{ϵ}$ . The eigenfunctions ${\eta }_{ϵ}$ live in different spaces $H_0^1\left({\Omega }_{ϵ}\right)$ . The normalization ${\Vert {\eta }_{ϵ}\Vert }_{L^2\left({\Omega }_{ϵ}\right)}=1$ also depends on $ϵ$ . To solve this problem, we transport the domain onto a fixed $\Omega$ domain, via the change of variables $x=T_{ϵ}^{-1}\left(y\right)$ . This will enable us to manipulate objects in a fixed frame. To this end, we introduce the following parameterization

$H_0^1\left({\Omega }_{ϵ}\right)=\left\{\phi \circ T_{ϵ}^{-1},\phi \in H_0^1\left(\Omega \right)\right\}.$ (5.6)

and to work in $H_0^1\left(\Omega \right)$ , we also introduce the compositions ${\eta }^{ϵ}={\eta }_{ϵ}\circ T_{ϵ}$ , ${\lambda }^{ϵ}={\lambda }_{ϵ}\circ T_{ϵ}$ and $v=v^{\prime }\circ T_{ϵ}$ . So the variational formulation of (5.4) amounts to looking for ${\eta }^{ϵ}\in H_0^1\left(\Omega \right)$ such that

${\displaystyle {\int }_{{\Omega }_{ϵ}}\nabla }\left({\eta }^{ϵ}\circ T_{ϵ}^{-1}\right)\nabla \left(v\circ T_{ϵ}^{-1}\right)-\left({\lambda }^{ϵ}\circ T_{ϵ}^{-1}\right)\left({\eta }^{ϵ}\circ T_{ϵ}^{-1}\right)\left(v\circ T_{ϵ}^{-1}\right)\text{d}x=0,\forall v\in H_0^1\left(\Omega \right).$ (5.7)

Applying the change of variables formula, we obtain

${\displaystyle {\int }_{\Omega }\left\{\nabla \left({\eta }^{ϵ}\circ T_{ϵ}^{-1}\right)\nabla \left(v\circ T_{ϵ}^{-1}\right)-\left({\lambda }^{ϵ}\circ T_{ϵ}^{-1}\right)\left({\eta }^{ϵ}\circ T_{ϵ}^{-1}\right)\left(v\circ T_{ϵ}^{-1}\right)\right\}\circ T_{ϵ}J\left(ϵ\right)\text{d}x}=0.$

In addition, we have

$\begin{array}{l}{\displaystyle {\int }_{\Omega }\left\{\left({\lambda }^{ϵ}\circ T_{ϵ}^{-1}\right)\left({\eta }^{ϵ}\circ T_{ϵ}^{-1}\right)\left(v\circ T_{ϵ}^{-1}\right)\right\}\circ T_{ϵ}J\left(ϵ\right)\text{d}x}\\ ={\displaystyle {\int }_{\Omega }\left({\lambda }^{ϵ}\circ T_{ϵ}^{-1}\circ T_{ϵ}\right)\left({\eta }^{ϵ}\circ T_{ϵ}^{-1}\circ T_{ϵ}\right)\left(v\circ T_{ϵ}^{-1}\circ T_{ϵ}\right)J\left(ϵ\right)\text{d}x}\\ ={\displaystyle {\int }_{\Omega }{\lambda }^{ϵ}{\eta }^{ϵ}vJ\left(ϵ\right)\text{d}x}.\end{array}$

Let $D{T}_{ϵ}\left(X\right)$ be the Jacobian matrix of $T_{ϵ}$ evaluated in $X$ , ${}^{\text{*}}D{T}_{ϵ}\left(X\right)$ the transpose matrix of $D{T}_{ϵ}\left(X\right)$ , then we have the following proposition

Proposition 5.1 We have

1) $\left(\nabla \varphi \right)\circ T_{ϵ}= {}^{\ast }D{T}_{ϵ}^{-1}\cdot \nabla \left(\varphi \circ T_{ϵ}\right),\forall \varphi \in C^1\left({ℝ}^N\right),$

2) $D\left(T\circ S\right)=\left\{\left(DT\right)\circ S\right\}\cdot DS,\forall \left(T,S\right)\in C^1\left({ℝ}^N;{ℝ}^N\right).$

Proof. See [1]. ■

So using the proposition, we can write

$\begin{array}{l}{\displaystyle {\int }_{\Omega }\left\{\nabla \left({\eta }^{ϵ}\circ T_{ϵ}^{-1}\right)\nabla \left(v\circ T_{ϵ}^{-1}\right)\right\}\circ T_{ϵ}J\left(ϵ\right)\text{d}x}\\ ={\displaystyle {\int }_{\Omega }D{T}_{ϵ}^{-1}{\left(D{T}_{ϵ}^{-1}\right)}^{\text{*}}\nabla {\eta }^{ϵ}\nabla vJ\left(ϵ\right)\text{d}x}\\ ={\displaystyle {\int }_{\Omega }B\left(ϵ\right)\nabla {\eta }^{ϵ}\nabla v\text{d}x}.\end{array}$

Finally, the previous equation becomes

${\displaystyle {\int }_{\Omega }B\left(ϵ\right)\nabla {\eta }^{ϵ}\nabla v\text{d}x}-{\displaystyle {\int }_{\Omega }{\lambda }^{ϵ}{\eta }^{ϵ}vJ\left(ϵ\right)\text{d}x}=0$ (5.8)

with $B\left(ϵ\right)=J\left(ϵ\right)D{T}_{ϵ}^{-1}{\left(D{T}_{ϵ}^{-1}\right)}^{\text{*}}$ , $J\left(ϵ\right)=\det D{T}_{ϵ}$ , $D{T}_{ϵ}$ the Jacobian matrix of $T_{ϵ}$ .

To write the functional in the fixed domain $\Omega$ , we first assume that A is written $A=\nabla =\tilde{\eta }$ and note by ${\tilde{\eta }}^{ϵ}=\tilde{\eta }\circ T_{ϵ}$ .

$\begin{array}{c}J\left({\Omega }_{ϵ}\right)={\displaystyle {\int }_{{\Omega }_{ϵ}}{\left|\nabla {\eta }_{ϵ}-A\right|}^2\text{d}x}+{\displaystyle {\int }_{{\Omega }_{ϵ}}{\left|{\eta }_{ϵ}-{\eta }_d\right|}^2\text{d}x}\\ ={\displaystyle {\int }_{{\Omega }_{ϵ}}{\left|\nabla {\eta }_{ϵ}-\nabla \tilde{\eta }\right|}^2\text{d}x}+{\displaystyle {\int }_{{\Omega }_{ϵ}}{\left|{\eta }_{ϵ}-{\eta }_d\right|}^2\text{d}x}\\ ={\displaystyle {\int }_{{\Omega }_{ϵ}}\left(\nabla {\eta }_{ϵ}-\nabla \tilde{\eta }\right)\cdot \left(\nabla {\eta }_{ϵ}-\nabla \tilde{\eta }\right)\text{d}x}+{\displaystyle {\int }_{{\Omega }_{ϵ}}{\left|{\eta }_{ϵ}-{\eta }_d\right|}^2\text{d}x}\\ ={\displaystyle {\int }_{\Omega }\left\{\left(\nabla {\eta }_{ϵ}-\nabla \tilde{\eta }\right)\cdot \left(\nabla {\eta }_{ϵ}-\nabla \tilde{\eta }\right)\right\}\circ T_{ϵ}J\left(ϵ\right)\text{d}x}\\ +{\displaystyle {\int }_{\Omega }{\left|{\eta }_{ϵ}\circ T_{ϵ}-{\eta }_d\circ T_{ϵ}\right|}^{2}J\left(ϵ\right)\text{d}x}.\end{array}$

Applying Proposition 5.1 once more, we obtain

$\begin{array}{l}g\left(ϵ\right)=J\left({\Omega }_{ϵ}\right)={\displaystyle {\int }_{\Omega }\left\{B\left(ϵ\right)\left(\nabla {\eta }^{ϵ}-\nabla {\tilde{\eta }}^{ϵ}\right)\right\}\cdot \left(\nabla {\eta }^{ϵ}-\nabla {\tilde{\eta }}^{ϵ}\right)\text{d}x}\\ +{\displaystyle {\int }_{\Omega }{\left|{\eta }^{ϵ}-{\eta }_d\circ T_{ϵ}\right|}^{2}J\left(ϵ\right)\text{d}x}.\end{array}$(5.9)

Since the set of perturbed states admits a solution, then $g\left(ϵ\right)<+\infty$ . So the $ϵ$ -dependent Lagrangian can be defined as follows:

$\begin{array}{c}ℒ\left(ϵ,\phi ,\varphi \right)={\displaystyle {\int }_{\Omega }\left\{B\left(ϵ\right)\left(\nabla \phi -\nabla {\tilde{\eta }}^{ϵ}\right)\right\}\cdot \left(\nabla \phi -\nabla {\tilde{\eta }}^{ϵ}\right)+{\left|\phi -{\eta }_d\circ T_{ϵ}\right|}^{2}J\left(ϵ\right)\text{d}x}\\ +{\displaystyle {\int }_{\Omega }B\left(ϵ\right)\nabla \phi \cdot \nabla \varphi -{\lambda }^{ϵ}\phi \varphi J\left(ϵ\right)\text{d}x}\end{array}$

$g\left(ϵ\right)=\underset{\phi \in H_0^1\left(\Omega \right)}{\inf }\underset{\varphi \in H_0^1\left(\Omega \right)}{\sup }ℒ\left(ϵ,\phi ,\varphi \right);dg\left(0\right)=\underset{ϵ\searrow 0}{\text{lim}}\left(g\left(ϵ\right)-g\left(0\right)\right)/ϵ=DJ\left(\Omega ,V\left(0\right)\right).$

Let $R\left(ϵ\right)$ be the following function

$\begin{array}{c}R\left(ϵ\right)={\displaystyle {\int }_{\Omega }2}\left(\frac{B\left(ϵ\right)-I}{ϵ}\right)\left\{\nabla \left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}\right)-\nabla {\tilde{\eta }}^{ϵ}\right\}\nabla \left({\eta }^{ϵ}-{\eta }^0\right)\text{d}x\\ +ϵ{\displaystyle {\int }_{\Omega }{\left|\nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\right|}^2\text{d}x}+{\displaystyle {\int }_{\Omega }2\left(\frac{J\left(ϵ\right)-1}{ϵ}\right)\left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}-{\eta }_d\circ T_{ϵ}\right)\left({\eta }^{ϵ}-{\eta }^0\right)\text{d}x}\\ -{\displaystyle {\int }_{\Omega }2\left(\nabla {\tilde{\eta }}^{ϵ}-\nabla \tilde{\eta }\right)\nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}-{\displaystyle {\int }_{\Omega }2\left({\eta }_d\circ T_{ϵ}-{\eta }_d\right)\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}\\ +ϵ{\displaystyle {\int }_{\Omega }{\left|\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right|}^2\text{d}x}+{\displaystyle {\int }_{\Omega }\frac{B\left(ϵ\right)-I}{ϵ}\nabla p^0\nabla \left({\eta }^{ϵ}-{\eta }^0\right)}\\ -{\lambda }^0\left(\frac{J\left(ϵ\right)-1}{ϵ}\right)p^0\left({\eta }^{ϵ}-{\eta }^0\right)\text{d}x-{\displaystyle {\int }_{\Omega }J\left(ϵ\right)}\left(\frac{{\lambda }^{ϵ}-{\lambda }^0}{ϵ}\right)p^0\left({\eta }^{ϵ}-{\eta }^0\right)\text{d}x.\end{array}$

We assume that the function $R\left(ϵ\right)$ admits a finite limit noted $R\left({\eta }^0,p^0\right)$ , then we have the following proposition.

Proposition 5.2 Let $\Omega$ be the solution to the problem ${\min }_{\Omega \in {\mathcal{O}}_{ad}}J\left(\Omega \right)$ , then the shape derivative of the functional $J\left(\Omega \right)$ defined in (5.2) is given by

$\begin{array}{c}DJ\left(\text{Ω},V\right)=ℛ\left({\eta }^0,p^0\right)+{\displaystyle {\int }_{\Omega }\left\{B^{\prime }\left(0\right)\left(\nabla {\eta }^0-\nabla \tilde{\eta }\right)\right\}\cdot \left(\nabla {\eta }^0-\nabla \tilde{\eta }\right)\text{d}x}\\ -{\displaystyle {\int }_{\Omega }2\left(\nabla {\eta }^0-\nabla \tilde{\eta }\right)\cdot \nabla \left(\left[\nabla \tilde{\eta }\cdot V\left(0\right)\right]\right)\text{d}x}\\ +{\displaystyle {\int }_{\Omega }{\left|{\eta }^0-{\eta }_d\right|}^{2}divV\left(0\right)-2\left({\eta }^0-{\eta }_d\right)\nabla {\eta }_d\cdot V\left(0\right)\text{d}x}\\ +{\displaystyle {\int }_{\Omega }{B}^{\prime }\left(0\right)\nabla {\eta }^0\cdot \nabla p^0}-{\lambda }^0{\eta }^0{p}^{0}divV\left(0\right)-{\left({\lambda }^{ϵ}\right)}^{\prime }\left(0\right){\eta }^0{p}^0\text{d}x.\end{array}$

where ${\eta }^0$ solution of (5.1) and $p^0$ , the adjoint state solution of

${\displaystyle {\int }_{\Omega }2\left(\nabla {\eta }^0-\nabla \tilde{\eta }\right)\cdot \nabla {\phi }^{\prime }+2\left({\eta }^0-{\eta }_d\right){\phi }^{\prime }+\nabla p^0\cdot \nabla {\phi }^{\prime }-{\lambda }^0{p}^0{\phi }^{\prime }\text{d}x}=0.$

Proof. The derivative of the $ϵ$ -dependent Lagrangian with respect to $ϵ$ .

$\begin{array}{c}{d}_{ϵ}ℒ\left(ϵ,\phi ,\varphi \right)={\displaystyle {\int }_{\Omega }\left\{B^{\prime }\left(ϵ\right)\left(\nabla \phi -\nabla {\tilde{\eta }}^{ϵ}\right)\right\}\cdot \left(\nabla \phi -\nabla {\tilde{\eta }}^{ϵ}\right)\text{d}x}\\ -{\displaystyle {\int }_{\Omega }2B\left(ϵ\right)\left(\nabla \phi -\nabla {\tilde{\eta }}^{ϵ}\right)\cdot \nabla \left(\left[\nabla \tilde{\eta }\cdot V\left(ϵ\right)\right]\circ T_{ϵ}\right)\text{d}x}\\ +{\displaystyle {\int }_{\Omega }{\left|\phi -{\eta }_d\circ T_{ϵ}\right|}^{2}div{V}_{ϵ}J\left(ϵ\right)-2\left(\phi -{\eta }_d\circ T_{ϵ}\right)\nabla {\eta }_d\cdot V_{ϵ}J\left(ϵ\right)\text{d}x}\\ +{\displaystyle {\int }_{\Omega }B\left(ϵ\right)\nabla \phi \cdot \nabla \varphi -{\lambda }^{ϵ}\phi \varphi div{V}_{ϵ}J\left(ϵ\right)-{\left({\lambda }^{ϵ}\right)}^{\prime }\phi \varphi J\left(ϵ\right)\text{d}x}.\end{array}$

Taking $\phi ={\eta }^0$ and $\varphi =p^0$ , we obtain

$\begin{array}{c}{d}_{ϵ}ℒ\left(ϵ,{\eta }^0,p^0\right)={\displaystyle {\int }_{\Omega }\left\{B^{\prime }\left(ϵ\right)\left(\nabla {\eta }^0-\nabla {\tilde{\eta }}^{ϵ}\right)\right\}\cdot \left(\nabla {\eta }^0-\nabla {\tilde{\eta }}^{ϵ}\right)\text{d}x}\\ -{\displaystyle {\int }_{\Omega }2B}\left(ϵ\right)\left(\nabla {\eta }^0-\nabla {\tilde{\eta }}^{ϵ}\right)\cdot \nabla \left(\left[\nabla \tilde{\eta }\cdot V\left(ϵ\right)\right]\circ T_{ϵ}\right)\text{d}x\\ +{\displaystyle {\int }_{\Omega }{\left|{\eta }^0-{\eta }_d\circ T_{ϵ}\right|}^{2}div{V}_{ϵ}J\left(ϵ\right)-2\left({\eta }^0-{\eta }_d\circ T_{ϵ}\right)\nabla {\eta }_d\cdot V_{ϵ}J\left(ϵ\right)\text{d}x}\\ +{\displaystyle {\int }_{\Omega }{B}^{\prime }}\left(ϵ\right)\nabla {\eta }^0\cdot \nabla p^0-{\lambda }^0{\eta }^0{p}^{0}div{V}_{ϵ}J\left(ϵ\right)-{\left({\lambda }^{ϵ}\right)}^{\prime }\left(0\right){\eta }^0{p}^0{J}_{ϵ}\text{d}x.\end{array}$

Taking $ϵ$ equal to zero, we get

$\begin{array}{c}{d}_{ϵ}ℒ\left(0,{\eta }^0,p^0\right)={\displaystyle {\int }_{\Omega }\left\{B^{\prime }\left(0\right)\left(\nabla {\eta }^0-\nabla \tilde{\eta }\right)\right\}\cdot \left(\nabla {\eta }^0-\nabla \tilde{\eta }\right)\text{d}x}\\ -{\displaystyle {\int }_{\Omega }2}\left(\nabla {\eta }^0-\nabla \tilde{\eta }\right)\cdot \nabla \left(\left[\nabla \tilde{\eta }\cdot V\left(0\right)\right]\right)\text{d}x\\ +{\displaystyle {\int }_{\Omega }{\left|{\eta }^0-{\eta }_d\right|}^{2}divV\left(0\right)-2\left({\eta }^0-{\eta }_d\right)\nabla {\eta }_d\cdot V\left(0\right)\text{d}x}\\ +{\displaystyle {\int }_{\Omega }{B}^{\prime }}\left(0\right)\nabla {\eta }^0\cdot \nabla p^0-{\lambda }^0{\eta }^0{p}^{0}divV\left(0\right)-{\left({\lambda }^{ϵ}\right)}^{\prime }\left(0\right){\eta }^0{p}^0\text{d}x.\end{array}$

The derivative of the $ϵ$ -dependent Lagrangian with respect to $\phi$ in the direction ${\phi }^{\prime }$ :

$\begin{array}{c}{d}_{\phi }ℒ\left(ϵ,\phi ,\varphi ;{\phi }^{\prime }\right)={\displaystyle {\int }_{\Omega }2B}\left(ϵ\right)\left(\nabla \phi -\nabla {\tilde{\eta }}^{ϵ}\right)\cdot \nabla {\phi }^{\prime }+2\left(\phi -{\eta }_d\circ T_{ϵ}\right){\phi }^{\prime }J\left(ϵ\right)\text{d}x\\ +{\displaystyle {\int }_{\Omega }B}\left(ϵ\right)\nabla {\phi }^{\prime }\cdot \nabla \varphi -{\lambda }^{ϵ}{\phi }^{\prime }\varphi J\left(ϵ\right)\text{d}x.\end{array}$

The derivative of the $ϵ$ -dependent Lagrangian with respect to $\varphi$ in the direction ${\varphi }^{\prime }$ :

$d_{\varphi }ℒ\left(ϵ,\phi ,\varphi ;{\varphi }^{\prime }\right)={\displaystyle {\int }_{\Omega }B}\left(ϵ\right)\nabla \phi \cdot \nabla {\varphi }^{\prime }-{\lambda }^{ϵ}\phi {\varphi }^{\prime }J\left(ϵ\right)\text{d}x.$ (5.10)

The state equation at $ϵ\ge 0$ and the adjoint state equation at $ϵ=0$ are respectively given by

${\eta }^{ϵ}\in H_0^1\left(\Omega \right),\forall {\varphi }^{\prime }\in H_0^1\left(\Omega \right),{\displaystyle {\int }_{\Omega }B\left(ϵ\right)\nabla {\eta }^{ϵ}\cdot \nabla {\varphi }^{\prime }-{\lambda }^{ϵ}{\eta }^{ϵ}{\varphi }^{\prime }J\left(ϵ\right)\text{d}x}=0$ (5.11)

$p^0\in H_0^1\left(\Omega \right),\forall {\phi }^{\prime }\in H_0^1\left(\Omega \right),$

${\displaystyle {\int }_{\Omega }2}\left(\nabla {\eta }^0-\nabla \tilde{\eta }\right)\cdot \nabla {\phi }^{\prime }+2\left({\eta }^0-{\eta }_d\right){\phi }^{\prime }+\nabla p^0\cdot \nabla {\phi }^{\prime }-{\lambda }^0{p}^0{\phi }^{\prime }\text{d}x=0.$(5.12)

$\begin{array}{c}ℛ\left(ϵ\right)={\displaystyle {\int }_0^1{d}_{\phi }ℒ}\left(ϵ,{\eta }^0+\theta \left({\eta }^{ϵ}-{\eta }^0\right),p^0;\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}\theta \\ ={\displaystyle {\int }_{\Omega }2B}\left(ϵ\right)\left\{\nabla \left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}\right)-\nabla {\tilde{\eta }}^{ϵ}\right\}\cdot \nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x\\ +{\displaystyle {\int }_{\Omega }2}\left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}-{\eta }_d\circ T_{ϵ}\right)\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)J\left(ϵ\right)\text{d}x\\ +{\displaystyle {\int }_{\Omega }B}\left(ϵ\right)\nabla p^0\cdot \nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)-{\lambda }^{ϵ}{p}^0\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)J\left(ϵ\right)\text{d}x.\end{array}$

Taking ${\phi }^{\prime }=\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}$ in adjoint Equation (5.12), we have

$\begin{array}{l}{\displaystyle {\int }_{\Omega }2\nabla }\left({\eta }^0-\tilde{\eta }\right)\nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)+2\left({\eta }^0-{\eta }_d\right)\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x\\ +{\displaystyle {\int }_{\Omega }\nabla p^0\cdot \nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)-{\lambda }^0{p}^0\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}=0.\end{array}$

Consequently, we get

$\begin{array}{l}{\displaystyle {\int }_{\Omega }2\left\{\nabla \left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}\right)-\nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{2}\right)-\nabla {\tilde{\eta }}^{ϵ}+\nabla {\tilde{\eta }}^{ϵ}-\nabla \tilde{\eta }\right\}\cdot \nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}\\ +{\displaystyle {\int }_{\Omega }2\left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}-\frac{{\eta }^{ϵ}-{\eta }^0}{2}-{\eta }_d\circ T_{ϵ}+{\eta }_d\circ T_{ϵ}-{\eta }_d\right)\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}\\ +{\displaystyle {\int }_{\Omega }\nabla p^0\cdot \nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)-{\lambda }^0{p}^0\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}=0.\end{array}$

In other words,

$\begin{array}{l}{\displaystyle {\int }_{\Omega }2\left\{\nabla \left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}\right)-\nabla {\tilde{\eta }}^{ϵ}\right\}\cdot \nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}-ϵ{\displaystyle {\int }_{\Omega }{\left|\nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\right|}^2\text{d}x}\\ +{\displaystyle {\int }_{\Omega }2\left({\eta }_d\circ T_{ϵ}-{\eta }_d\right)\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}+{\displaystyle {\int }_{\Omega }2\left(\nabla {\tilde{\eta }}^{ϵ}-\nabla \tilde{\eta }\right)\cdot \nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}\\ +{\displaystyle {\int }_{\Omega }2\left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}-{\eta }_d\circ T_{ϵ}\right)\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}-ϵ{\displaystyle {\int }_{\Omega }{\left|\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right|}^2\text{d}x}\\ +{\displaystyle {\int }_{\Omega }\nabla p^0\cdot \nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)-{\lambda }^0{p}^0\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}=0.\end{array}$

Accordingly, the function $R\left(ϵ\right)$ can be rewritten as follows:

$\begin{array}{c}ℛ\left(ϵ\right)={\displaystyle {\int }_{\Omega }2\left(\frac{B\left(ϵ\right)-I}{ϵ}\right)\left\{\nabla \left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}\right)-\nabla {\tilde{\eta }}^{ϵ}\right\}\nabla \left({\eta }^{ϵ}-{\eta }^0\right)\text{d}x}\\ +ϵ{\displaystyle {\int }_{\Omega }{\left|\nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\right|}^2\text{d}x}+{\displaystyle {\int }_{\Omega }2}\left(\frac{J\left(ϵ\right)-1}{ϵ}\right)\left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}-{\eta }_d\circ T_{ϵ}\right)\left({\eta }^{ϵ}-{\eta }^0\right)\text{d}x\\ -{\displaystyle {\int }_{\Omega }2\left(\nabla {\tilde{\eta }}^{ϵ}-\nabla \tilde{\eta }\right)\nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}-{\displaystyle {\int }_{\Omega }2\left({\eta }_d\circ T_{ϵ}-{\eta }_d\right)\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x}\\ +ϵ{\displaystyle {\int }_{\Omega }{\left|\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right|}^2\text{d}x}+{\displaystyle {\int }_{\Omega }\frac{B\left(ϵ\right)-I}{ϵ}\nabla p^0\nabla \left({\eta }^{ϵ}-{\eta }^0\right)}\\ -\left(\frac{{\lambda }^{ϵ}J\left(ϵ\right)-{\lambda }^0}{ϵ}\right)p^0\left({\eta }^{ϵ}-{\eta }^0\right)\text{d}x.\end{array}$

$\begin{array}{c}ℛ\left(ϵ\right)={\displaystyle {\int }_{\Omega }2\left(\frac{B\left(ϵ\right)-I}{ϵ}\right)\left\{\nabla \left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}\right)-\nabla {\tilde{\eta }}^{ϵ}\right\}\nabla \left({\eta }^{ϵ}-{\eta }^0\right)\text{d}x}\\ +ϵ{\displaystyle {\int }_{\Omega }{\left|\nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\right|}^2\text{d}x}+{\displaystyle {\int }_{\Omega }2}\left(\frac{J\left(ϵ\right)-1}{ϵ}\right)\left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}-{\eta }_d\circ T_{ϵ}\right)\left({\eta }^{ϵ}-{\eta }^0\right)\text{d}x\\ -{\displaystyle {\int }_{\Omega }2}\left(\nabla {\tilde{\eta }}^{ϵ}-\nabla \tilde{\eta }\right)\nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x-{\displaystyle {\int }_{\Omega }2}\left({\eta }_d\circ T_{ϵ}-{\eta }_d\right)\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\text{d}x\\ +ϵ{\displaystyle {\int }_{\Omega }{\left|\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right|}^2\text{d}x}+{\displaystyle {\int }_{\Omega }\frac{B\left(ϵ\right)-I}{ϵ}\nabla p^0\nabla \left({\eta }^{ϵ}-{\eta }^0\right)}\\ -{\lambda }^0\left(\frac{J\left(ϵ\right)-1}{ϵ}\right)p^0\left({\eta }^{ϵ}-{\eta }^0\right)\text{d}x-{\displaystyle {\int }_{\Omega }J\left(ϵ\right)}\left(\frac{{\lambda }^{ϵ}-{\lambda }^0}{ϵ}\right)p^0\left({\eta }^{ϵ}-{\eta }^0\right)\text{d}x.\end{array}$

As in the case of the Helmholtz equation studied in [9], we have the following inequality

$\begin{array}{c}ℛ\left(ϵ\right)\le 2\Vert \frac{B\left(ϵ\right)-I}{ϵ}\Vert \Vert \nabla \left(\frac{{\eta }^{ϵ}+{\eta }^0}{2}\right)-\nabla {\tilde{\eta }}^{ϵ}\Vert \Vert \nabla \left({\eta }^{ϵ}-{\eta }^0\right)\Vert +ϵ\Vert \nabla {\left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)}^2\Vert \\ +2\Vert \frac{J\left(ϵ\right)-1}{ϵ}\Vert \Vert \frac{{\eta }^{ϵ}+{\eta }^0}{2}-{\eta }_d\circ T_{ϵ}\Vert \Vert {\eta }^{ϵ}-{\eta }^0\Vert +2\Vert \nabla {\tilde{\eta }}^{ϵ}-\nabla \tilde{\eta }\Vert \Vert \nabla \left(\frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\right)\Vert \\ +2\Vert \frac{{\eta }_d\circ T_{ϵ}-{\eta }_d}{ϵ}\Vert \Vert {\eta }^{ϵ}-{\eta }^0\Vert +ϵ{\Vert \frac{{\eta }^{ϵ}-{\eta }^0}{ϵ}\Vert }^2+\Vert \frac{B\left(ϵ\right)-I}{ϵ}\Vert \Vert \nabla p^0\Vert \Vert \nabla \left({\eta }^{ϵ}-{\eta }^0\right)\Vert \\ +{\lambda }^0\Vert \frac{J\left(ϵ\right)-1}{ϵ}\Vert \Vert p^0\Vert \Vert {\eta }^{ϵ}-{\eta }^0\Vert +\Vert J\left(ϵ\right)\Vert \Vert \frac{{\lambda }^{ϵ}-{\lambda }^0}{ϵ}\Vert \Vert p^0\Vert \Vert {\eta }^{ϵ}-{\eta }^0\Vert .\end{array}$

Now we need to show that the limit of $ℛ\left(ϵ\right)$ exists. First, we know the terms $\left(B\left(ϵ\right)-I\right)/ϵ$ and $\left(J\left(ϵ\right)-1\right)/ϵ$ are uniformly bounded. It will therefore suffice to show that ${\eta }^{ϵ}\to {\eta }^0$ in $H_0^1\left(\Omega \right)$ -strong and that the $H_0^1\left(\Omega \right)$ -norm of $\left({\eta }^{ϵ}-{\eta }^0\right)/ϵ$ is bounded. But the difficulty lies in the fact that the eigenvalue ${\lambda }^{ϵ}$ depends on the variable $ϵ$ and we have no information on whether $\left({\lambda }^{ϵ}-{\lambda }^0\right)/ϵ$ is bounded. So we just assume that $ℛ\left(ϵ\right)$ admits a finite limit denoted $R\left({\eta }^0,p^0\right)$ and in this case the shape derivative of the functional $J\left(\Omega \right)$ is given by:

$\begin{array}{c}DJ\left(\Omega ,V\right)=R\left({\eta }^0,p^0\right)+{\displaystyle {\int }_{\Omega }\left\{B^{\prime }\left(0\right)\left(\nabla {\eta }^0-\nabla \tilde{\eta }\right)\right\}\cdot \left(\nabla {\eta }^0-\nabla \tilde{\eta }\right)\text{d}x}\\ -{\displaystyle {\int }_{\Omega }2}\left(\nabla {\eta }^0-\nabla \tilde{\eta }\right)\cdot \nabla \left(\left[\nabla \tilde{\eta }\cdot V\left(0\right)\right]\right)\text{d}x\\ +{\displaystyle {\int }_{\Omega }{\left|{\eta }^0-{\eta }_d\right|}^{2}divV\left(0\right)}-2\left({\eta }^0-{\eta }_d\right)\nabla {\eta }_d\cdot V\left(0\right)\text{d}x\\ +{\displaystyle {\int }_{\Omega }{B}^{\prime }\left(0\right)\nabla {\eta }^0\cdot \nabla p^0}-{\lambda }^0{\eta }^0{p}^{0}divV\left(0\right)-{\left({\lambda }^{ϵ}\right)}^{\prime }\left(0\right){\eta }^0{p}^0\text{d}x.\end{array}$

where ${\eta }^0$ solution of (5.1) and $p^0$ , the adjoint state solution of

${\displaystyle {\int }_{\Omega }2}\left(\nabla {\eta }^0-\nabla \tilde{\eta }\right)\cdot \nabla {\phi }^{\prime }+2\left({\eta }^0-{\eta }_d\right){\phi }^{\prime }+\nabla p^0\cdot \nabla {\phi }^{\prime }-{\lambda }^0{p}^0{\phi }^{\prime }\text{d}x=0.$

■

6. Conclusion

This work provides a simple method for calculating the shape derivative in the context of Dirichlet eigenvalue problems. We have used alternative techniques to calculate the shape derivative of the first eigenvalue of the Laplace-Dirichlet problem, as well as that of an objective function constrained by a spectral problem. A fundamental tool in shape optimization, the shape derivative enables us to analyze how a quantity of interest varies under the effect of small deformations of the domain. This approach, inspired in particular by the work of Delfour, opens the way to a better understanding and control of geometry-dependent spectral problems. Some numerical tests were also carried out in two dimensions. The next step is to calculate the topological derivative using the same approach for an eigenvalue problem.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Delfour, M. and Zolésio, J.P. (2001) Shapes and Geometries: Metrics, Analysis, Differential Calculus, and Optimization (Advances in Design and Control). SIAM.
[2]	Henrot, A. (2006) Extremum Problems for Eigenvalues of Elliptic Operators. Birkha üser Verlag. https://doi.org/10.1007/3-7643-7706-2
[3]	Henrot, A. (2017) Shape Optimization and Spectral Theory. De Gruyter, 464.
[4]	Henrot, A. and Pierre, M. (2006) Variation et optimisation de formes: Une analyse géométrique, volume 48. Springer Science & Business Media. https://doi.org/10.1007/3-540-37689-5
[5]	Zhu, S. (2017) Effective Shape Optimization of Laplace Eigenvalue Problems Using Domain Expressions of Eulerian Derivatives. Journal of Optimization Theory and Applications, 176, 17-34. https://doi.org/10.1007/s10957-017-1198-9
[6]	Ngom, M.G., Faye, I. and Seck, D. (2024) On Shape and Topological Optimization Problems with Constraints Helmholtz Equation and Spectral Problems. Results in Applied Mathematics, 23, Article ID: 100474. https://doi.org/10.1016/j.rinam.2024.100474
[7]	Delfour, M.C. (2018) Control, Shape, and Topological Derivatives via Minimax Differentiability of Lagrangians. In: Falcone, M., Ferretti, R., Grüne, L. and McEneaney, W., Eds., Numerical Methods for Optimal Control Problems, Springer, 137-164. https://doi.org/10.1007/978-3-030-01959-4_7
[8]	Delfour, M.C. (2021) Topological Derivatives via One-Sided Derivative of Parametrized Minima and Minimax. Engineering Computations, 39, 34-59. https://doi.org/10.1108/ec-06-2021-0318
[9]	Ngom, M.G., Faye, I. and Seck, D. (2023) A Minimax Method in Shape and Topological Derivatives and Homogenization: The Case of Helmholtz Equation. Nonlinear Studies, 30, page.
[10]	Ngom, M.G., Faye, I. and Seck, D. (2025) Topological Derivative for Shallow Water Equations. Nonlinear Studies, 32, page.
[11]	Evans, L.C. (1998) Partial Differential Equations, Graduate Studies in Mathematics. American Mathematical Society, 19.
[12]	Schulz, V.H., Siebenborn, M. and Welker, K. (2016) Efficient PDE Constrained Shape Optimization Based on Steklov—Poincaré-Type Metrics. SIAM Journal on Optimization, 26, 2800-2819. https://doi.org/10.1137/15m1029369
[13]	Langtangen, H.P. and Logg, A. (2017) Solving PDEs in Python: The FEniCS Tutorial I. Springer.