[論文レビュー] Error analysis for a finite difference scheme for axisymmetric mean curvature flow of genus-0 surfaces

Klaus Deckelnick, Robert Nürnberg|arXiv (Cornell University)|Oct 19, 2020

Geometric Analysis and Curvature Flows参考文献 16被引用数 4

ひとこと要約

本稿では、 genus-0 表面の軸対称平均曲率流れに対して、回転軸における退化を丁寧に取り扱うことで、有限差分スキームを提示する。離散 L2 および H1 範囲における最適誤差境界を確立し、数値的収束実験および自己収縮解や非埋め込み表面のシミュレーションによって検証されている。

ABSTRACT

We consider a finite difference approximation of mean curvature flow for axisymmetric surfaces of genus zero. A careful treatment of the degeneracy at the axis of rotation for the one dimensional partial differential equation for a parameterization of the generating curve allows us to prove error bounds with respect to discrete $L^2$- and $H^1$-norms for a fully discrete approximation. The theoretical results are confirmed with the help of numerical convergence experiments. We also present numerical simulations for some genus-0 surfaces, including for a non-embedded self-shrinker for mean curvature flow.

研究の動機と目的

• 回転軸に対して直角に交わる開曲線である genus-0 表面の軸対称平均曲率流れに対して、安定的かつ高精度な有限差分スキームを開発すること。
• 回転軸において半径座標が消えるため、支配方程式に特異性が生じる問題に対処すること。
• 境界での特異的挙動にもかかわらず、完全離散有限差分スキームについて、離散 L2 および H1 範囲における厳密な誤差推定を確立すること。
• 滑らかでない初期データおよび特異的初期データに対する数値的収束実験を通じて、理論的誤差境界を検証すること。
• 非埋め込み表面、特に円錐特異性を有する非埋め込み自己収縮解に対しても、スキームの頑健性を示すこと。

提案手法

• DeTurckのトリックを用いて、生成曲線を極座標と軸座標で記述する1次元PDEに定式化し、厳密な放物型性を達成する。
• PDEを発散型に変換することで、自然な変分法的定式化と安定な有限差分離散化を可能にする。
• 時間方向に半implicitな後退オイラー法を適用し、空間方向に区分的線形有限差分を用いることで、安定性を確保する。
• 回転軸における特異性に対処するため、L’Hôpitalの定理を用いて速度および曲率項の整合性のある境界条件を導出する。
• 座標変換と軸付近のメッシュ細分化を用いることで、精度を維持し、メッシュのねじれを防ぐ。
• 均一な空間的・時間的グリッドを用いた完全離散スキームを実装し、数値実験による収束性を検証する。

実験結果

リサーチクエスチョン

• RQ1 回転軸における退化を考慮した場合、genus-0 表面の軸対称平均曲率流れに対して、有限差分スキームの最適誤差境界を確立できるか？
• RQ2 回転軸におけるPDEの特異性を数値的に取り扱うにはどうすればよいか？ これにより精度と安定性を維持できるか？
• RQ3 非埋め込み表面や、90°接触角条件を満たさない初期データに対しても、有限差分スキームは収束性と安定性を保つのか？
• RQ4 自己相似収縮解、特に円錐特異性を有する非埋め込み自己収縮解を、正確にシミュレートできるか？
• RQ5 離散L2およびH1範囲におけるスキームの収束速度は何か？ また、理論的予測と一致するか？

主な発見

• 回転軸における退化にもかかわらず、完全離散有限差分スキームについて、離散L2およびH1範囲における最適誤差境界が証明された。
• 数値実験により、空間方向で2次収束、時間方向で1次収束が確認された（J = 512、Δt = 10−4 で実施）。
• 表面積は時間とともに線形に減少し、消滅時刻は約1.0となり、自己収縮解の理論的予測と整合的である。
• 円錐特異性を有する非埋め込み初期データに対しても、スキームは正常に処理できる：外向きの円錐は速やかに滑らかになり、内向きの円錐は収縮半円に近づく。
• 初期に接触角条件が満たされなくても、数値解は極限において90°接触角条件を満たしており、頑健性と自己修復性を示している。
• L’Hôpitalの定理を用いて境界条件を一貫的に導出することで、従来の手法が失敗した領域においても誤差制御が可能となり、genus-0 表面への理論的拡張が達成された。