A Lagrangian shape and topology optimization framework based on semi-discrete optimal transport

Charles Dapogny,Bruno Levy,Edouard Oudet
2024-09-12
Abstract:This article revolves around shape and topology optimization, in the applicative context where the objective and constraint functionals depend on the solution to a physical boundary value problem posed on the optimized domain. We introduce a novel framework based on modern concepts from computational geometry, optimal transport and numerical analysis. Its pivotal feature is a representation of the optimized shape by the cells of an adapted version of a Laguerre diagram. Although such objects are originally described by a collection of seed points and weights, recent results from optimal transport theory suggest a more intuitive parametrization in terms of the seed points and measures of the associated cells. The polygonal mesh of the shape induced by this diagram serves as support for the deployment of the Virtual Element Method for the numerical solution of the physical boundary value problem at play and the calculation of the objective and constraint functionals. The sensitivities of the latter are derived next; at first, we calculate their derivatives with respect to the positions of the vertices of the Laguerre diagram by shape calculus techniques; a suitable adjoint methodology is then developed to express them in terms of the seed points and cell measures of the diagram. The evolution of the shape is realized by first updating the design variables according to these sensitivities and then reconstructing the diagram with efficient algorithms from computational geometry. Our shape optimization strategy is versatile: it can be applied to a wide gammut of physical situations. It is Lagrangian by essence, and it thereby benefits from all the assets of a consistently meshed representation of the shape. Yet, it naturally handles dramatic motions, including topological changes, in a very robust fashion. These features, among others, are illustrated by a series of 2d numerical examples.
Optimization and Control
What problem does this paper attempt to address?
### Problems Addressed by the Paper This paper focuses on the issues of shape and topology optimization, particularly in the context where the objective and constraint functions depend on the solutions of physical boundary value problems defined over the optimization domain. The authors introduce a new numerical framework based on modern concepts from computational geometry, optimal transport, and numerical analysis. The core feature of this framework is the representation of optimized shapes through cells of an adapted version of the Laguerre diagram. Specifically, the framework addresses the following key issues: 1. **Shape Representation and Update**: - Using the Laguerre diagram to represent optimized shapes, which is defined by a set of seed points and weights. - Proposing a more intuitive parameterization method using the latest results from optimal transport theory, representing the Laguerre diagram with seed points and cell metrics. - Achieving shape evolution by reconstructing the updated Laguerre diagram through efficient algorithms from computational geometry. 2. **Solving Physical Boundary Value Problems**: - Utilizing the discrete polygonal mesh induced by the Laguerre diagram, employing the Virtual Element Method (VEM) to numerically solve boundary value problems describing physical phenomena. - This method can handle various physical situations, such as thermodynamics and structural mechanics problems. 3. **Sensitivity Analysis of Optimization Criteria**: - Calculating the sensitivity of optimization criteria and their derivatives with respect to the positions of the Laguerre diagram vertices using shape calculus techniques. - Developing appropriate adjoint methods to express these sensitivities as functions of seed points and cell metrics. 4. **Shape Evolution in the Optimization Process**: - Updating design variables based on sensitivity information, then regenerating the Laguerre diagram to achieve iterative shape optimization. - This method is inherently Lagrangian, allowing precise tracking of shape changes at each stage while naturally handling severe deformations and topological changes. 5. **Application Scope**: - The framework has broad applicability and can be applied to various physical situations (such as heat conduction and structural mechanics) and optimization settings (including single-phase and multi-phase problems). ### Summary This paper proposes a new framework for shape and topology optimization, combining the latest advances in computational geometry, optimal transport, and numerical analysis. By using the Laguerre diagram and the Virtual Element Method, this framework can accurately represent shapes, efficiently handle physical boundary value problems, and achieve robust shape optimization. This approach is not only applicable to various physical problems but also naturally handles complex shape changes and topological variations.