Convergence analysis for the Barrett--Garcke--Nurnberg method of transport type for evolving curves
Genming Bai, Harald Garcke, Shravan Veerapaneni
TL;DR
This work establishes the first convergence result for a fully discrete Barrett--Garcke--Nürnberg scheme solving curve evolution under general, transport-dominated flows. By leveraging a projection-error framework and Gauss–Lobatto mass lumping, the authors derive discrete transport-energy estimates that compensate for the lack of $H^1$ parabolicity. They prove a sub-optimal $L^2$ convergence rate $\| \hat{e}_h^{m} \|_{L^2} \le C(\tau + h^k)$ for degree $k\ge 3$, under a mild time-step constraint, and show stability through high-order tangential estimates and an intrinsic $H^2$ stability of the discrete Laplacian. The results advance the theory of parametric finite element methods for transport-type interface evolution and provide a robust discretization framework for moving interfaces driven by general flows, with potential extensions to higher dimensions and related interface problems.
Abstract
In this paper, we propose a Barrett-Garcke-Nurnberg (BGN) method for evolving geometries under general flows and present the corresponding convergence analysis. Specifically, we examine the scenario where a closed curve evolves according to a prescribed background velocity field. Unlike mean curvature flow and surface diffusion, where the evolution velocities inherently exhibit parabolicity, this case is dominated by transport which poses a significant difficulty in establishing convergence proofs. To address the challenges imposed by this transport-dominant nature, we derive several discrete energy estimates of the transport type on discretized polynomial surfaces within the framework of the projection error. The use of the projection error is indispensable as it provides crucial additional stability through its orthogonality structure. We prove that the proposed method converges sub-optimally in the L2 norm, and this is the first convergence proof for a fully discrete numerical method solving the evolution of curves driven by general flows.
