[论文解读] The Augmented Fast Marching Method for Level Set Reinitialization

David Salač|arXiv (Cornell University)|Nov 29, 2011

Advanced Numerical Methods in Computational Mathematics参考文献 11被引用 3

一句话总结

本文提出了一种新型的水平集重初始化技术——增强型快速推进法（AFMM），可高精度计算符号距离函数及其一阶与二阶导数。通过强制满足 |∇ϕ|² = 1，同时确保 ∇(|∇ϕ|²) = 0 和 ∇∇(|∇ϕ|²) = 0，该方法即使在粗网格上也能保证曲率场的平滑性，实现水平集与梯度场的二阶精度，曲率场的一阶精度。

ABSTRACT

Including derivative information in the modelling of moving interfaces has been proposed as one method to increase the accuracy of numerical schemes with minimal additional cost. Here a new level set reinitialization technique using the fast marching method is presented. This augmented fast marching method will calculate the signed distance function and up to the second-order derivatives of the signed distance function for arbitrary interfaces. In addition to enforcing the condition $| ablaϕ|^2=1$, where $ϕ$ is the level set function, the method ensures that $ abla(| ablaϕ|)^2=0$ and $ abla abla(| ablaϕ|)^2=0$ are also satisfied. Results indicate that for both two- and three-dimensional interfaces the resulting level set and curvature field are smooth even for coarse grids. Convergence results show that using first-order upwind derivatives and the augmented fast marching method result in a second-order accurate level set and gradient field and a first-order accurate curvature field.

研究动机与目标

• 解决移动界面问题中对精确且平滑水平集重初始化的需求，特别是在曲率及其导数至关重要的场景中。
• 克服标准快速推进法在重初始化过程中无法保持高阶导数信息的局限性。
• 确保水平集函数及其导数在满足狭义方程与几何约束的前提下保持平滑且一致。
• 即使在粗网格上，也实现水平集与梯度场的二阶精度，同时保持曲率场的一阶精度。
• 支持复杂界面（如囊泡）的数值模拟，其中弯曲力依赖于水平集函数的高阶导数。

提出的方法

• 通过求解 Eikonal 方程 |∇ϕ| = 1，同时对 |∇ϕ|² 的梯度与海森矩阵施加附加约束，将快速推进法（FMM）扩展以包含导数信息。
• 引入增强状态变量，以在快速推进传播过程中追踪符号距离函数的一阶与二阶导数。
• 使用迎风有限差分法近似 ϕ、∇ϕ 与 ∇∇ϕ 的空间导数，确保求解过程的因果性与稳定性。
• 在整个算法中保持三个节点集合——已接受、待评估与遥远节点——以确保正确的传播顺序与收敛性。
• 通过在每个网格点强制要求 ∇(|∇ϕ|²) = 0 与 ∇∇(|∇ϕ|²) = 0 来实现几何一致性，从而确保梯度模长的平滑性。
• 在二维与三维域中应用该方法，对 Eikonal 方程与导数约束采用一致的离散化方案。

实验结果

研究问题

• RQ1 能否将快速推进法扩展，以高精度计算符号距离函数及其一阶与二阶导数？
• RQ2 强制满足 ∇(|∇ϕ|²) = 0 与 ∇∇(|∇ϕ|²) = 0 是否能带来更平滑的曲率场并提升水平集模拟的数值稳定性？
• RQ3 在结构化笛卡尔网格上，该增强方法在水平集函数、其梯度与曲率上的收敛速率如何？
• RQ4 该方法是否能在粗网格上实现水平集与梯度场的二阶精度，同时保持曲率场的一阶精度？
• RQ5 在模拟囊泡等复杂界面时，该方法表现如何，其中曲率及其导数决定了物理作用力？

主要发现

• 增强型快速推进法成功实现了符号距离函数及其至二阶导数的高精度计算。
• 该方法确保 ∇(|∇ϕ|²) = 0 与 ∇∇(|∇ϕ|²) = 0 成立，即使在粗网格上也能生成平滑的曲率场。
• 收敛性研究显示，当使用一阶迎风差分时，水平集与梯度场达到二阶精度。
• 曲率场被证实为一阶精度，与该类方法的理论预期一致。
• 二维与三维的数值结果证实，该方法可生成平滑、稳定且几何一致的水平集表示。
• 通过保留高阶导数信息，该方法可实现由曲率相关力主导的界面（如囊泡膜）的精确模拟。