[论文解读] On the self-similar solutions of the crystalline mean curvature flow in three dimensions

Norbert Požár|arXiv (Cornell University)|Jun 7, 2018

Geometric Analysis and Curvature Flows参考文献 43被引用 5

一句话总结

本文提出三维晶体平均曲率流的正亏格自相似收缩解，利用它们验证基于最小化运动格式的有限元法（FEM）水平集算法。FEM离散化在三维中改善了边缘处理，无需正则化即可准确捕捉复杂动力学，如晶面破裂、弯曲及拓扑变化。

ABSTRACT

We present two types of self-similar shrinking solutions of positive genus for the crystalline mean curvature flow in three dimensions analogous to the solutions known for the standard mean curvature flow. We use them to test a numerical implementation of a level set algorithm for the crystalline mean curvature flow in three dimensions based on the minimizing movements scheme of A.~Chambolle, extit{Interfaces Free Bound.~6}~(2004). We implement a finite element method discretization that seems to improve the handling of edges in three dimensions compared to the standard finite difference method and illustrate its behavior on a few examples.

研究动机与目标

• 构建三维晶体平均曲率流的正亏格自相似收缩解。
• 通过最小化运动格式测试并验证水平集方法在晶体平均曲率流中的数值实现。
• 比较有限元法（FEM）与有限差分法（FDM）在三维中处理边缘和奇点的性能。
• 展示该方法在无需正则化的情况下捕捉晶面破裂、弯曲及拓扑变化（如颈部缩颈）的能力。

提出的方法

• 采用Chambolle（2004）提出的最小化运动格式实现水平集方法，用于晶体平均曲率流。
• 在总变差能量中采用有限元法（FEM）离散化，替代标准有限差分法（FDM），以改善三维中的边缘处理。
• 使用Bregman迭代格式求解最小化运动框架中产生的子问题。
• 采用避免重距离化伪影的方案对距离函数进行重初始化。
• 使用自相似解——六边形Wulff形状、环形结构及海绵状结构——作为基准测试案例。
• 通过收敛率和计算成本评估数值性能，指标包括每时间步的Bregman迭代次数和每次迭代耗时。

实验结果

研究问题

• RQ1 能否在三维晶体平均曲率流中构造出正亏格的自相似收缩解？
• RQ2 FEM离散化在三维晶体平均曲率流中与FDM相比，能否更有效地保持锐利边缘和顶点？
• RQ3 该数值格式能否在无需正则化的情况下准确模拟晶面破裂、弯曲及拓扑变化（如缩颈）？
• RQ4 FEM-based水平集方法在3D晶体流中的收敛行为和计算效率如何？

主要发现

• 对于非立方各向异性，FEM离散化在视觉质量上优于FDM，尤其在保持锐利边缘和顶点方面，尽管计算成本更高。
• 收敛性研究证实，该方法在空间和时间上均达到一阶精度。
• 对于三维六边形Wulff形状，最大误差随网格加密而减小，收敛率趋近于$ M^{-1} $。
• 该方法成功模拟了复杂动力学：晶面破裂（图14）、晶面弯曲（图15）及拓扑变化（如缩颈，图16、18），且无需正则化。
• 重距离化伪影极少，重初始化方案在所有测试案例中均表现稳健。
• 性能基准显示，FEM每次Bregman迭代耗时高于FDM，但几何保真度更优。