Optimization of the shape of regions supporting boundary conditions

International audienceThis article deals with the optimization of the shape of the regions assigned to different types of boundary conditions in the definition of a 'physical' partial differential equation. At first, we analyze a model situation involving the solution u Ω to a Laplace equation in a domain Ω; the boundary ∂Ω is divided into three parts Γ D , Γ and Γ N , supporting respectively homogeneous Dirichlet, homogeneous Neumann and inhomogeneous Neumann boundary conditions. The shape derivative J (Ω)(θ) of a general objective function J(Ω) of the domain is calculated in the framework of Hadamard's method when the considered deformations θ are allowed to modify the geometry of Γ D , Γ and Γ N (i.e. θ does not vanish on the boundary of these regions). The structure of this shape derivative turns out to depend very much on the regularity of u Ω near the boundaries of the regions Γ D , Γ and Γ N. For this reason, in particular, J (Ω) is difficult to calculate and to evaluate numerically when the transition Γ D ∩ Γ between homogeneous Dirichlet and homogeneous Neumann boundary conditions is subject to optimization. To overcome this difficulty, an approximation method is proposed, in which the considered 'exact' Laplace equation with mixed boundary conditions is replaced with a 'smoothed' version, featuring Robin boundary conditions on the whole boundary ∂Ω with coefficients depending on a small parameter ε. We prove the consistency of this approach in our model context: the approximate objective function Jε(Ω) and its shape derivative converge to their exact counterparts as ε vanishes. Although it is rigorously justified only in a model problem, this approximation methodology may be adapted to many more complex situations, for example in three space dimensions, or in the context of the linearized elasticity system. Various numerical examples are eventually presented in order to appraise the efficiency of the proposed approximation process