Stable discretizations of elastic flow in Riemannian manifolds

John W. Barrett,Harald Garcke,Robert Nürnberg
DOI: https://doi.org/10.1137/18M1227111
2019-06-14
Abstract:The elastic flow, which is the $L^2$-gradient flow of the elastic energy, has several applications in geometry and elasticity theory. We present stable discretizations for the elastic flow in two-dimensional Riemannian manifolds that are conformally flat, i.e.\ conformally equivalent to the Euclidean space. Examples include the hyperbolic plane, the hyperbolic disk, the elliptic plane as well as any conformal parameterization of a two-dimensional manifold in ${\mathbb R}^d$, $d\geq 3$. Numerical results show the robustness of the method, as well as quadratic convergence with respect to the space discretization.
Numerical Analysis,Differential Geometry
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to achieve a stable discretization of elastic flow in two - dimensional Riemannian manifolds. Specifically, the authors focus on how to numerically simulate the L2 - gradient flow of elastic energy in conformally flat Riemannian manifolds and ensure that this numerical method is stable. ### Analysis of the Main Problem 1. **Definition and Importance of Elastic Flow**: - Elastic flow is the L2 - gradient flow of elastic energy and has wide applications in geometry and elasticity theory. - For a two - dimensional Riemannian manifold \((M, g)\), the elastic flow can be expressed as: \[ V_g=-(\kappa_g)_{s_g s_g}-\frac{1}{2}\kappa_g^3 - S_0\kappa_g \] where \(V_g\) is the normal velocity of the curve with respect to the metric \(g\), \(\kappa_g\) is the geodesic curvature, and \(S_0\) is the sectional curvature. 2. **Limitations of Existing Methods**: - Early methods would lead to errors in the normal direction when dealing with curves on hypersurfaces embedded in three - dimensional Euclidean space. - These methods cannot guarantee stability, especially when calculating elastic flow in non - embedded manifolds. 3. **Contributions of This Paper**: - New numerical methods are proposed, which can be directly calculated on the manifold without being embedded in three - dimensional Euclidean space. - By introducing semi - discrete continuous - time finite - element approximations, the stability of these methods is proved. - The effectiveness and robustness of these methods in calculating elastic flow in complex manifolds such as the hyperbolic plane are demonstrated. ### Summary of Mathematical Formulas - **Elastic Energy**: \[ W_{g,\lambda}(\vec{x})=\frac{1}{2}\int_I(\kappa_g^2 + 2\lambda)|\vec{x}_\rho|_g d\rho \] where \(\kappa_g\) is the geodesic curvature of the curve with respect to the metric \(g\). - **Strong Form of Elastic Flow**: \[ V_g = g^{1/2}(\vec{x})\vec{x}_t\cdot\vec{\nu}=-(\kappa_g)_{s_g s_g}-\frac{1}{2}\kappa_g^3 - S_0\kappa_g+\lambda\kappa_g \] - **Weak Form**: The paper proposes two weak forms (U) and (W), based on \(\kappa\) and \(\kappa_g\) as variables respectively. For example, the weak form (W) based on \(\kappa_g\) includes: \[ \int_I g(\vec{x})\vec{x}_t\cdot\vec{\nu}\chi|\vec{x}_\rho|_g d\rho=\int_I g^{-1/2}(\vec{x})[\kappa_g]_\rho\chi_\rho|\vec{x}_\rho|^{-1} d\rho-\frac{1}{2}\int_I g^{1/2}(\vec{x})\kappa_g^3\chi|\vec{x}_\rho| d\rho-\int_I S_0(\vec{x})g^{1/2}(\vec{x})\kappa_g\chi|\vec{x}_\rho| d\rho \] Through these methods, the authors have successfully solved the stability problem of numerical simulation of elastic flow in complex geometric structures.