Optimization of the shape of regions supporting boundary conditions

• Published in: Numerische Mathematik Vol. 146; no. 1; pp. 51 - 104
• Main Authors: Dapogny, Charles, Lebbe, Nicolas, Oudet, Edouard
• Format: Journal Article
• Language: English
• Published: Berlin/Heidelberg Springer Berlin Heidelberg 01.09.2020; Springer Nature B.V; Springer Verlag
• Subjects: Approximation; Boundary conditions; Context; Dirichlet problem; Domains; Laplace equation; Mathematical and Computational Engineering; Mathematical and Computational Physics; Mathematical Methods in Physics; Mathematical models; Mathematics; Mathematics and Statistics; Numerical Analysis; Numerical and Computational Physics; Optimization; Optimization and Control; Partial differential equations; Simulation; Theoretical; 35B65; 49M41; 65K05; 49Q10
• ISSN: 0029-599X; 0945-3245
• DOI: 10.1007/s00211-020-01140-0

This 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.