A novel finite element approximation of anisotropic curve shortening flow

Klaus Deckelnick,Robert Nürnberg

DOI: https://doi.org/10.4171/IFB/500

2023-01-27

Abstract:We extend the DeTurck trick from the classical isotropic curve shortening flow to the anisotropic setting. Here the anisotropic energy density is allowed to depend on space, which allows an interpretation in the context of Finsler metrics, giving rise to e.g.\ geodesic curvature flow in Riemannian manifolds. Assuming that the density is strictly convex and smooth, we introduce a novel weak formulation for anisotropic curve shortening flow. We then derive an optimal $H^1$--error bound for a continuous-in-time semidiscrete finite element approximation that uses piecewise linear elements. In addition, we consider some fully practical fully discrete schemes and prove their unconditional stability. Finally, we present several numerical simulations, including some convergence experiments that confirm the derived error bound, as well as applications to crystalline curvature flow and geodesic curvature flow.

Numerical Analysis,Differential Geometry

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop a new finite - element approximation method for simulating the anisotropic curve shortening flow. Specifically, the author extends the classical DeTurck trick to adapt to the anisotropic situation and introduces a new weak - solution formula to describe this flow. The following is a detailed explanation of the problem:

1. **Anisotropic energy density**: The energy functional \(E(\Gamma)\) considered in the paper is defined as:
\[
E(\Gamma)=\int_{\Gamma} a(z)\gamma(z, \nu) \, dH^{1}(z)
\]
where \(\gamma(z, p)\) is the anisotropic function, \(a(z)\) is the weight function, and \(\nu\) is the unit normal vector of the curve \(\Gamma\).

2. **Natural gradient flow**: The natural gradient flow of this energy functional follows the following evolution law:
\[
V_{\gamma}=\kappa_{\gamma}
\]
where \(V_{\gamma}\) is the anisotropic normal velocity and \(\kappa_{\gamma}\) is the anisotropic curvature.

3. **Application of the DeTurck trick**: To deal with the degeneracy problem in the above - mentioned evolution law, the author introduces the DeTurck trick, which makes the system strictly parabolic in the tangential direction through re - parameterization. This makes the numerical analysis simpler and more effective.

4. **Error estimation and stability**: The author derives an optimal \(H^{1}\)-error bound and proves the unconditional stability of several fully - discrete finite - element schemes.

5. **Numerical simulation**: The paper also shows several numerical simulation results, including convergence experiments to verify the error bounds and examples applied to crystal curvature flow and geodesic curvature flow.

In summary, the main objective of this paper is to propose and analyze a new method based on the DeTurck trick for solving the curve - shortening - flow problem of general spatially non - uniform anisotropy and to provide theoretical and numerical support.