[論文レビュー] 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を強制することで、粗いグリッド上でも滑らかな曲率場を確保し、レベルセットおよび勾配場で2次精度、曲率で1次精度を達成する。

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方程式および幾何的制約と整合的で滑らかに保たれることを保証すること。
• 粗いグリッド上でも、レベルセットおよび勾配場で2次精度、曲率で1次精度を達成する。
• 曲率の高次導関数が物理的力に寄与する複雑な界面、例えばビシクルの数値的シミュレーションを支援すること。

提案手法

• Eikonal方程式 |∇ϕ| = 1 を解く際、|∇ϕ|² の勾配およびヘッセ行列に関する追加制約を課すことで、フットマーチング法（FMM）に導関数情報を拡張する。
• フットマーチング伝搬中に、符号付き距離関数の一次および二次の偏導関数を追跡するための拡張状態変数を導入する。
• 空間的微分を上流有限差分法で近似し、解法プロセスにおける因果性および安定性を確保する。
• 正しい伝搬順序と収束を保証するため、常に3つのノード集合（受理済み、試行中、遠方）を維持する。
• 各グリッド点で ∇(|∇ϕ|²) = 0 および ∇∇(|∇ϕ|²) = 0 を要求することで、勾配の大きさの滑らかさを強制し、幾何的整合性を確保する。
• 2次元および3次元領域の両方で手法を適用し、Eikonal方程式および導関数制約に対して一貫した離散化スキームを用いる。

実験結果

リサーチクエスチョン

• RQ1 フットマーチング法を、符号付き距離関数に加え、その一次および二次の偏導関数を高い精度で計算できるように拡張できるか？
• RQ2 ∇(|∇ϕ|²) = 0 および ∇∇(|∇ϕ|²) = 0 を強制することで、レベルセットシミュレーションにおける曲率場の滑らかさと数値的安定性が向上するか？
• RQ3 構造的直交グリッド上での、拡張手法のレベルセット関数、勾配、曲率の収束速度は何か？
• RQ4 粗いメッシュ上でも、レベルセットおよび勾配場で2次精度、曲率で1次精度を達成できるか？
• RQ5 曲率およびその導関数が物理的力を支配する複雑な界面、例えばビシクルのシミュレーションにおいて、本手法はどのように性能を発揮するか？

主な発見

• 拡張型フットマーチング法は、符号付き距離関数およびレベルセット関数の2次までの偏導関数を高い精度で計算できることを確認した。
• 本手法により、∇(|∇ϕ|²) = 0 および ∇∇(|∇ϕ|²) = 0 が満たされ、粗いグリッド上でも滑らかな曲率場が得られることを示した。
• 収束解析の結果、1次上流差分を用いた場合、レベルセットおよび勾配場で2次精度が達成された。
• 曲率場は理論的期待と整合的に1次精度であることが判明した。
• 2次元および3次元両方の数値的結果から、本手法が滑らかで安定的かつ幾何的に整合性のあるレベルセット表現を生成することが確認された。
• 高次導関数情報を保持できることから、本手法は曲率依存力に支配される界面、例えばビシクル膜の正確なシミュレーションを可能にした。