[論文レビュー] Convergence Analysis of A Second-order Semi-implicit Projection Method for Landau-Lifshitz Equation

Jingrun Chen, Cheng Wang+1 more|arXiv (Cornell University)|Feb 26, 2019

Electromagnetic Simulation and Numerical Methods参考文献 35被引用数 4

ひとこと要約

本稿では、2次精度の後退差分公式（BDF2）と片側補外を用い、磁化長の保存を保証するための射影ステップを組み合わせた、Landau-Lifshitz方程式の2次精度半陽的投影法を提案する。本手法は、時間・空間両方で2次収束を示し、時間刻みや空間刻みの制限なしに一意解の存在を保証するという、微磁気モデルにおいて初の結果である無条件安定性を確立する。

ABSTRACT

The numerical approximation for the Landau-Lifshitz equation, the dynamics of magnetization in a ferromagnetic material, is taken into consideration. This highly nonlinear equation, with a non-convex constraint, has several equivalent forms, and involves solving an auxiliary problem in the infinite domain. All these features have posed interesting challenges in developing numerical methods. In this paper, we first present a fully discrete semi-implicit method for solving the Landau-Lifshitz equation based on the second-order backward differentiation formula and the one-sided extrapolation (using previous time-step numerical values). A projection step is further used to preserve the length of the magnetization. Subsequently, we provide a rigorous convergence analysis for the fully discrete numerical solution, with second-order accuracy in both time and space, provided that the spatial step-size is the same order as the temporal step-size. And also, the unique solvability of the numerical solution is theoretically justified, which turns out to be the first such result for the micromagnetics model. All these theoretical properties are verified by numerical examples in both one- and three- dimensional spaces.

研究の動機と目的

• 非凸制約と無限領域における補助問題を有する、極めて非線形な系であるLandau-Lifshitz方程式のための堅牢な数値スキームの開発。
• 時間刻みが $k = \mathcal{O}(h^2)$ のような制限を伴う陰解法の限界を克服すること。
• 時間と空間の両方で2次精度を達成するとともに、無条件安定性および数値解の一意解の存在を保証すること。
• 誤差関数の線形化された安定性解析と、2組の補助的近似解を用いた厳密な収束解析を実施すること。
• 1次元および3次元領域における数値実験を通じて、ドメインウォールの動的挙動を含め、本手法の有効性を検証すること。

提案手法

• 時間離散化に2次精度の後退差分公式（BDF2）を用いた完全離散化された半陽的スキームを構築する。
• 右辺の非線形係数を、直前の時間ステップの値を用いた片側補外により線形化する。
• 各時間ステップで $|\mathbf{m}| = 1$ の制約を満たすために、射影ステップを適用し、磁化の大きさを保存する。
• 空間離散化には、時間と空間の両方で均等メッシュを用いた有限差分法を採用する。
• 誤差を評価するために、2組の補助的近似解を用いた厳密な収束解析を実施し、特に射影ステップにおける $W_h^{1,\infty}$ 評価に注目する。
• 線形化された系の単調性解析により、時間刻みや空間刻みの制限なしに一意解の存在を証明する。

実験結果

リサーチクエスチョン

• RQ1 制限的な時間刻み条件を回避できる、Landau-Lifshitz方程式の2次精度半陽的スキームを設計できるか？
• RQ2 本問題において、時間と空間の両方で無条件安定性と2次収束を確立できるか？
• RQ3 特に微磁気モデルにおいて、時間刻みや空間刻みの制限なしに数値スキームの一意解の存在を証明できるか？
• RQ4 本手法は、ドメインウォール移動のような複雑な磁化ダイナミクスのシミュレーションにおいて、どのように性能を発揮するか？
• RQ5 射影ステップにおける $W_h^{1,\infty}$ 誤差推定値が収束解析において果たす役割は何か？

主な発見

• 数値実験により、$L^\infty$ ノルムにおいて約1.92の収束次数を示し、時間と空間の両方で2次収束が達成されたことが確認された。
• 無条件安定性が数値的にも確認された：表6に示すように、時間刻み $k$ を空間刻み $h$ より著しく大きくしても、誤差は有界に保たれた。
• 時間刻みや空間刻みの制限なしに、数値解の一意解の存在が厳密に証明された。これは、微磁気モデルにおいて初の結果である。
• 本手法は、磁性ナノストリップにおける横方向ドメインウォールのダイナミクスを的確に捉え、外部磁場50 Oeの下で255 m/sの安定した速度で一定のプロファイルを維持して移動した。
• 図4の可視化結果から、射影ステップが磁化長を効果的に保存しており、$k=1/256$ および $h=1/32$ の条件下で正確な解と数値解の間に良好な一致が得られた。
• 収束解析は、射影ステップにおける誤差の $W_h^{1,\infty}$ 推定値に強く依存しており、これが安定性解析における最良の誤差バインディングを可能にした。