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

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.