As a specific type of shape gradient descent algorithm, shape gradient flow is widely used for shape optimization problems constrained by partial differential equations. In this approach, the constraint partial differential equations could be solved by finite element methods on a domain with a solution-driven evolving boundary. Rigorous analysis for the stability and convergence of such finite element approximations is still missing from the literature due to the complex nonlinear dependence of the boundary evolution on the solution. In this article, rigorous analysis of numerical approximations to the evolution of the boundary in a prototypical shape gradient flow is addressed. First-order convergence in time and kth order convergence in space for finite elements of degree k2 are proved for a linearly semi-implicit evolving finite element algorithm up to a given time. The theoretical analysis is consistent with the numerical experiments, which also illustrate the effectiveness of the proposed method in simulating two- and three-dimensional boundary evolution under shape gradient flow. The extension of the formulation, algorithm and analysis to more general shape density functions and constraint partial differential equations is also discussed.  

　　Publication:  

　　IMA Journal of Numerical Analysis, drad080 (27 October 2023)  

　　https://doi.org/10.1093/imanum/drad080  

　　Author:  

　　Wei Gong  

　　The State Key Laboratory of Scientific and Engineering Computing, Institute of Computational Mathematics and National Center for Mathematics and Interdisciplinary Sciences, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, 100190 Beijing, China  

　　Email: wgong@lsec.cc.ac.cn  

　　Buyang Li  

　　Department of Applied Mathematics, The Hong Kong Polytechnic University, Hung Hom, Hong Kong  

　　Qiqi Rao  

　　Department of Applied Mathematics, The Hong Kong Polytechnic University, Hung Hom, Hong Kong