A structure-preserving parametric finite element method with optimal energy stability condition for anisotropic surface diffusion

TL;DR

This paper introduces a structure-preserving finite element method for anisotropic surface diffusion that guarantees energy stability and conserves area, validated through numerical experiments.

Contribution

The paper develops a novel parametric finite element method that ensures unconditional energy stability and geometric structure preservation for anisotropic surface diffusion.

Findings

Unconditional energy stability under specific anisotropic conditions

Exact conservation of enclosed area during evolution

Numerical results confirm efficiency and accuracy

Abstract

We propose and analyze a structure-preserving parametric finite element method (SP-PFEM) for the evolution of closed curves under anisotropic surface diffusion with surface energy density $\hat{\gamma}(\theta)$. Our primary theoretical contribution establishes that the condition $3\hat{\gamma}(\theta)-\hat{\gamma}(\theta-\pi)\geq 0$ is both necessary and sufficient for unconditional energy stability within the framework of local energy estimates. The proposed method introduces a symmetric surface energy matrix $\hat{\boldsymbol{Z}}_k(\theta)$ with a stabilizing function $k(\theta)$, leading to a conservative weak formulation. Its fully discretization via SP-PFEM rigorously preserves the two geometric structures: enclosed area conservation and energy dissipation unconditionally under our energy stability condition. Numerical results are reported to demonstrate the efficiency and accuracy of the proposed method, along with its area conservation and energy dissipation properties.