[论文解读] Analysis of a space-time unfitted finite element method for PDEs on evolving surfaces

Arnold Reusken, Hauke Sass|arXiv (Cornell University)|Jan 1, 2024

Advanced Numerical Methods in Computational Mathematics被引用 1

一句话总结

本文提出了一种用于在演化曲面上求解PDE的时空非拟合有限元方法，基于水平集表示的高阶等参数映射。该方法在能量范数和L2范数下实现了最优收敛率，严格的误差分析证实了在标准假设下对光滑解的最优阶收敛，同时在拓扑变化情况下仍保持鲁棒性。

ABSTRACT

In this paper we analyze a space-time unfitted finite element method for the discretization of scalar surface partial differential equations on evolving surfaces. For higher order approximations of the evolving surface we use the technique of (iso)parametric mappings for which a level set representation of the evolving surface is essential. We derive basic results in which certain geometric characteristics of the exact space-time surface are related to corresponding ones of the numerical surface approximation. These results are used in an error analysis of a higher order space-time TraceFEM.

研究动机与目标

• 开发一种基于欧拉型、非拟合方法的高阶时空有限元方法，用于在演化曲面上求解PDE。
• 建立基于等参数映射和水平集表示的时空TraceFEM的严格误差分析。
• 解决非拟合方法在演化曲面上的几何逼近误差与积分一致性问题。
• 展示该方法在具有拓扑奇异性（如合并或分裂）问题中的鲁棒性。
• 为数值实验中观察到的最优收敛率提供理论依据。

提出的方法

• 该方法采用固定的、非拟合的体积分域时空网格，标准有限元空间被限制在演化曲面上。
• 水平集函数 φ(x,t) 定义演化曲面 S 为其零水平集，从而实现精确的几何表示。
• 通过基于近似水平集函数 φh ≈ φ 的等参数映射，构建高阶曲面近似 Sh。
• 时空TraceFEM在近似曲面 Sh 上使用边界有限元空间，通过体积法向导数项实现稳定化。
• 通过积分变换与分部积分公式，将几何量（如曲面测度和梯度）在精确曲面 S 与近似曲面 Sh 之间进行转换。
• 误差分析依赖于曲面测度变化、共法向跳跃和插值误差的估计，利用了近期几何分析的成果。

实验结果

研究问题

• RQ1 高阶时空非拟合有限元方法是否能在光滑几何与解的演化曲面上实现PDE的最优收敛？
• RQ2 基于水平集的等参数映射带来的几何逼近误差如何影响整体收敛率？
• RQ3 体积法向导数稳定化在控制条件数与实现高阶精度中起什么作用？
• RQ4 当存在拓扑奇异性（如合并或分裂）时，该方法是否仍能保持最优收敛？
• RQ5 尽管理论上预期，为何能量范数在时间节点上未表现出超收敛？

主要发现

• 该方法在能量范数下实现了最优收敛阶 ks，误差界 (7.35) 与数值实验结果一致，适用于 ks = 1, 2, 3, 4。
• L2 误差预计以阶数 ks+1 收敛，尽管目前尚未严格证明。
• 定理 7.8 中的误差界是该类高阶欧拉型时空非拟合方法的首个严格结果。
• 该方法在具有拓扑奇异性的问题中仍保持鲁棒性，但在这些情况下无法期望全局高阶收敛。
• 数值结果表明质量误差 m(t) ∼ h²（当 kg,s = ks = 1 且 ∆t ∼ h 时），表明质量守恒良好。
• 尽管理论上预期，但在时间节点上未观察到超收敛，该现象仍未得到解释。