An algorithm is proposed for generalized mean curvature flow of closed two-dimensional surfaces, which include inverse mean curvature flow, powers of mean and inverse mean curvature flow, etc. Error estimates are proven for semi- and full discretisations for the generalized flow. The algorithm proposed and studied here combines evolving surface finite elements, whose nodes determine the discrete surface, and linearly implicit backward difference formulae for time integration. The numerical method is based on a system coupling the surface evolution to non-linear second-order parabolic evolution equations for the normal velocity and normal vector. Convergence proofs are presented in the case of finite elements of polynomial degree at least two and backward difference formulae of orders two to five. The error analysis combines stability estimates and consistency estimates to yield optimal-order $H^1$-norm error bounds for the computed surface position, velocity, normal vector, normal velocity, and therefore for the mean curvature. The stability analysis is performed in the matrix-vector formulation, and is independent of geometric arguments, which only enter the consistency analysis. Numerical experiments are presented to illustrate the convergence results, and also to report on monotone quantities, e.g.~Hawking mass for inverse mean curvature flow. Complemented by experiments for non-convex surfaces.
翻译:对于封闭二维表面的普通平均曲线流,包括反平均曲线流、平均和反平均曲线流的功率等等,提出了一种算法。对于一般流动的半和完全分解,可以证明错误估计。这里所建议和研究的算法结合了变化中的表面有限元素,这些元素的节点决定离散表面,以及时间整合的线性隐含后向偏差公式。数字方法基于将表面进化与非线性二级的正常速度和正常矢量的二次分向演进方程相结合的系统。在单向水平的有限元素中至少提供了两个和后向偏差公式。错误分析结合了稳定性估计和一致性估计,以生成计算表层位置、速度、正常矢量、正常速度和平均曲度的最佳偏差界限。稳定性分析是在正常速度和正常矢量的组合中进行的,并且独立于只进入单向多向水平水平的多向值的组合度公式的参数参数的对比证据。 数值分析还以正向正向正态的轨道试验展示结果。