[論文レビュー] 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

ひとこと要約

本稿では、三次元における結晶的平均曲率流れの自己相似収縮解（正の genus を有するもの）を提示し、最小化の動きスキームに基づく有限要素法（FEM）を用いたレベルセット法アルゴリズムの妥当性を検証する。FEM 離散化は、有限差分法（FDM）に比べてエッジの取り扱いを改善し、正則化なしに複雑なダイナミクス（フェースの破壊、曲げ、トポロジー変化）を正確に捉えることができる。

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.

研究の動機と目的

• 三次元における結晶的平均曲率流れの正の genus を有する自己相似収縮解を構築すること。
• 最小化の動きスキームに基づくレベルセット法の数値実装を検証すること。
• 3次元におけるエッジおよび特異点の取り扱いにおいて、有限要素法（FEM）と有限差分法（FDM）の離散化を比較すること。
• 正則化なしにフェースの破壊、曲げ、トポロジー変化（例：首の細まり）を捉える能力を示すこと。

提案手法

• 結晶的平均曲率流れに対して、Chambolle（2004）の最小化の動きスキームに基づくレベルセット法を実装する。
• 全 Variation エネルギーの離散化に有限要素法（FEM）を用い、標準的な有限差分法（FDM）に代えて3次元におけるエッジの取り扱いを改善する。
• 最小化の動きフレームワークで生じる最小化サブプロブレムを解くために、Bregman反復スキームを用いる。
• 赤再距離化アーチファクトを回避するスキームを用いて距離関数を再初期化する。
• 自己相似解（六角形のWulff形状、ドーナツ型、スポンジ状構造）をベンチマークテストケースとして用いる。
• 収束率と計算コストを評価するための数値的性能評価を実施し、Bregman反復回数／時間ステップ、1反復あたりの時間などの指標を用いる。

実験結果

リサーチクエスチョン

• RQ1 三次元における結晶的平均曲率流れに対して、正の genus を有する自己相似収縮解を構築できるか？
• RQ2 FEMベースの離散化は、3次元結晶的平均曲率流れにおいて、鋭いエッジや頂点をFDMと比較してどれほどよく保持できるか？
• RQ3 正則化なしに、フェースの破壊、曲げ、首の細まりのようなトポロジー変化を正確にシミュレートできるか？
• RQ4 FEMベースのレベルセット法は、3次元結晶流れにおいて収束特性と計算効率の両面でどのように振る舞うか？

主な発見

• FEM 離散化は、立方体でない異方性に対して、視覚的品質が FDM を上回るが、計算コストは高くなる。特に鋭いエッジや頂点の保持が優れている。
• 自己相似解を用いた収束スタディにより、空間および時間において1次収束性を達成していることが確認された。
• 3次元における六角形Wulff形状では、メッシュの細分化に伴い最大誤差が減少し、$ M^{-1} $ の収束率に近づく。
• 本手法は、複雑なダイナミクスを成功裏にシミュレートした：フェースの破壊（図14）、フェースの曲げ（図15）、首の細まり（図16, 18）などのトポロジー変化を正則化なしに再現した。
• 赤再距離化アーチファクトは最小限に抑えられ、再初期化スキームはすべてのテストケースで安定に動作した。
• 性能ベンチマークでは、FEMはFDMに比べて1回のBregman反復に要する時間が長いが、幾何学的忠実性に優れた結果を得た。