Numerical shape optimization of the Canham-Helfrich-Evans bending energy
arXiv:2107.13794 · doi:10.1016/j.jcp.2023.112218
Abstract
In this paper we propose a novel numerical scheme for the Canham-Helfrich-Evans bending energy based on a three-field lifting procedure of the distributional shape operator to an auxiliary mean curvature field. Together with its energetic conjugate scalar stress field as Lagrange multiplier the resulting fourth order problem is circumvented and reduced to a mixed saddle point problem involving only second order differential operators. Further, we derive its analytical first variation (also called first shape derivative), which is valid for arbitrary polynomial order, and discuss how the arising shape derivatives can be computed automatically in the finite element software NGSolve. We finish the paper with several numerical simulations showing the pertinence of the proposed scheme and method.
References in corpus (9)
- Modelling fluid deformable surfaces with an emphasis on biological interfaces
- Efficient PDE constrained shape optimization based on Steklov-Poincaré type metrics
- Viscous regularization and r-adaptive remeshing for finite element analysis of lipid membrane mechanics
- A stabilized finite element formulation for liquid shells and its application to lipid bilayers
- Bending models of lipid bilayer membranes: spontaneous curvature and area-difference elasticity
- Shape optimisation with nearly conformal transformations
- The Hellan-Herrmann-Johnson Method for Nonlinear Shells
- Nonlinear Conjugate Gradient Methods for PDE Constrained Shape Optimization Based on Steklov-Poincaré-Type Metrics
- Uniform convergence of discrete curvatures from nets of curvature lines