[論文レビュー] A generalized MBO diffusion generated motion for orthogonal matrix-valued fields

Braxton Osting, Dong Wang|arXiv (Cornell University)|Nov 3, 2017

Characterization and Applications of Magnetic Nanoparticles参考文献 34被引用数 6

ひとこと要約

本稿では、多様体上の直交行列値場におけるディリクレエネルギーの局所的最小値を計算するための一般化されたMBO拡散生成運動アルゴリズムを導入する。拡散ステップとSVDを用いた点毎の直交射影を交互に繰り返すことで、リャプノフ関数を介して無条件安定性を保証し、空間離散化された設定において有限時間収束を達成する。この方法は、平坦なトーラスや球面上でアレン・カーン方程式やギンツブルグ＝ランダウ方程式と類似した振る舞いを示す。

ABSTRACT

We consider the problem of finding stationary points of the Dirichlet energy for orthogonal matrix-valued fields. Following the Ginzburg-Landau approach, this energy is relaxed by penalizing the matrix-valued field when it does not take orthogonal matrix values. A generalization of the MBO diffusion generated motion is introduced that effectively finds local minimizers of this energy by iterating two steps until convergence. In the first step, as in the original method, the current matrix-valued field is evolved by the diffusion equation. In the second step, the field is pointwise reassigned to the closest orthogonal matrix, which can be computed via the singular value decomposition. We extend the Lyapunov function of Esedoglu and Otto to show that the method is non-increasing on iterates and hence, unconditionally stable. We also prove that spatially discretized iterates converge to a stationary solution in a finite number of iterations. The algorithm is implemented using the closest point method and non-uniform fast Fourier transform. We conclude with several numerical experiments on flat tori and closed surfaces, which, unsurprisingly, exhibit classical behavior from the Allen-Cahn and complex Ginzburg Landau equations, but also new phenomena.

研究の動機と目的

• 閉じた多様体上の直交行列値場におけるディリクレエネルギーの局所的最小値を求める安定で反復的なアルゴリズムの開発。
• SVDを用いた直交射影と拡散を組み合わせることで、MBOフレームワークを行列値場に一般化すること。
• エゼドグルとオットーの結果を一般化する拡張されたリャプノフ関数を用いて、無条件安定性を証明すること。
• 空間離散化された反復ステップが定常解に有限時間で収束することを確立すること。
• 平坦なトーラスおよび閉じた曲面上での数値的実装と検証を実施し、体積保存型ダイナミクスを含む。

提案手法

• 行列場を熱方程式（拡散ステップ）で進めるのと、各行列成分をSVDを用いて最も近い直交行列に射影する操作を交互に繰り返す。
• 直交性からの逸脱をペナルティ項 $ \|A^tA - I_n\|_F^2 $ を含むギンツブルグ＝ランダウ型の緩和項を用いてエネルギー関数を緩和する。
• 緩和エネルギーの勾配流は $ \partial_t A = \Delta A - \varepsilon^{-2} A(A^tA - I_n) $ として導出され、$ n=1 $ の場合にアレン・カーン方程式に一般化される。
• 反復間でエネルギーが単調に減少することを保証するためのリャプノフ関数を構築し、無条件安定性を証明する。
• 空間離散化には最近接点法を用い、拡散ステップにおける畳み込みの効率的計算には高速フーリエ変換を採用する。
• 曲面上の体積保存型ダイナミクスのため、クイックソートに基づくアルゴリズムを用いて、集合 $ \Omega_\lambda^+ $ が所望の体積を持つように閾値 $ \lambda $ を計算する。

実験結果

リサーチクエスチョン

• RQ1 MBOフレームワークは、直交制約を満たす行列値場に一般化可能であり、安定性と収束性を保つことができるか？
• RQ2 提案手法は、直交行列値場に適用した場合に、エネルギー減少と無条件安定性を維持するか？
• RQ3 一般化されたMBO法のダイナミクスは、高次元における古典的アレン・カーン方程式や複素ギンツブルグ＝ランダウ方程式とどのように比較できるか？
• RQ4 球面などの曲がった表面における直交行列値場における線状欠陥や渦度の振る舞いはいかなるものか？
• RQ5 本手法は、表面問題において体積制約を満たすように適合可能か？

主な発見

• リャプノフ関数の構築により、反復間でエネルギーが単調に減少することを証明し、一般化されたMBO法が無条件安定であることを確認した。
• 空間離散化された反復ステップが有限回のステップで定常解に収束することを確認し、収束保証を裏付けた。
• 平坦なトーラスおよび球面上での数値実験から、線状欠陥が円形配置に進化することが観察され、古典的ギンツブルグ＝ランダウダイナミクスと整合的であった。
• 本手法は、古典的アレン・カーン挙動を超えた新たな現象、例えば位相欠陥や行列場における相分離を効果的に捉えた。
• 最近接点法と非一様FFTを用いた効率的実装により、複雑な幾何構造上でも大規模なシミュレーションが可能となった。
• 曲面上の体積保存型ダイナミクスでは、クイックソートに基づく閾値化アルゴリズムが、所望の体積を高精度で維持できた。