A convergent algorithm for mean curvature flow of surfaces with Dirichlet boundary conditions

We establish convergence results for a spatial semidiscretization of Mean Curvature Flow (MCF) for surfaces with fixed boundaries. Our analysis is based on Huisken's evolution equations for the mean curvature and the normal vector, enabling precise control of discretization errors and yielding optimal error estimates for discrete spaces with piecewise polynomials of degree $p \geq 2$. Building on techniques recently developed by Kov\'acs, Li, Lubich, and collaborators for closed surfaces, we extend these ideas to surfaces with boundaries by formulating appropriate boundary conditions for both the mean curvature and the normal vector. These boundary treatments are essential for proving convergence. The core of our analysis involves a classical error splitting strategy using auxiliary discrete functions that approximate the surface geometry, the mean curvature, and the normal vector. We estimate two types of errors for each variable to rigorously assess both stability and consistency. To effectively handle boundary conditions for the normal vector, we introduce a nonlinear Ritz projection into the analysis. As a result, we derive optimal $H^1$ error estimates for the surface position, velocity, mean curvature, and normal vector. Our theoretical findings are corroborated by numerical experiments.