[論文レビュー] An explicit and practically invariants-preserving method for conservative systems

Wenjun Cai, Yuezheng Gong+1 more|arXiv (Cornell University)|Sep 15, 2020

Numerical methods for differential equations参考文献 43被引用数 4

ひとこと要約

本稿では、高次精度の陽的ルンゲ・クッタ法と簡素化された射影技術を組み合わせることで、保存系における不変量を保存する明示的不変量保存法（EIP）を提案する。この手法は、元のルンゲ・クッタ法と同一の精度順序を維持し、不変量の誤差が $\mathcal{O}(h^{2(p+1)})$ のオーダーに抑えられ、十分に高次の方法では機械精度に達する。

ABSTRACT

An explicit numerical strategy that practically preserves invariants is derived for conservative systems by combining an explicit high-order Runge-Kutta (RK) scheme with a simple modification of the standard projection approach, which is named the explicit invariants-preserving (EIP) method. The proposed approach is shown to have the same order as the underlying RK method, while the error of invariants is analyzed in the order of $\mathcal{O}\left(h^{2(p+1)} ight),$ where $h$ is the time step and $p$ represents the order of the method. When $p$ is appropriately large, the EIP method is practically invariants-conserving because the error of invariants can reach the machine accuracy. The method is illustrated for the cases of single and multiple invariants, with regard to both ODEs and high-dimensional PDEs. Extensive numerical experiments are presented to verify our theoretical results and demonstrate the superior behaviors of the proposed method in a long time numerical simulation. Numerical results suggest that the fourth-order EIP method preserves much better the qualitative properties of the flow than the standard fourth-order RK method and it is more efficient in practice than the fully implicit integrators.

研究の動機と目的

• 非線形方程式の反復解法を必要とせずに、保存系における不変量を保存する明示的数値解法を開発すること。
• 標準的な陽的ルンゲ・クッタ法が線形および二次形式を超える一般の多項式不変量を保存できないという制限を克服すること。
• 簡素化された射影アプローチにより、保存量の誤差を機械精度まで低減することで、実用的な不変量保存を達成すること。
• 不変量保存法の適用範囲を高次元PDEおよび複数の不変量を有する系へ拡張すること。
• 完全に陰的解法および標準的な陽的解法と比較して、長時間にわたる数値的安定性と効率性に優れていることを示すこと。

提案手法

• 明示的高次ルンゲ・クッタ（RK）法と、非線形方程式の解法を回避する修正された射影技術を組み合わせる。
• 射影ステップでは、不変量の離散勾配に基づく簡素化されたラグランジュ乗数法を用いて、解を調整することで不変量を強制的に満たす。
• 標準的な射影法の実用的簡略化として導出され、元のRK法の精度順序を保持する。
• 単一および複数の不変量（ハミルトニアン系におけるエネルギーと質量を含む）を扱えるように一般化される。
• PDEに対しては、空間方向の半離散化後に適用され、グロス＝ピタエフケィ方程式のような高次元系でも不変量を保存する。
• 一般の不変量を二次形式に再定式化するための2次化戦略を用いて、射影による明示的処理を可能にする。

実験結果

リサーチクエスチョン

• RQ1 明示的数値解法は、基本積分法の精度順序を維持したまま、保存系における一般の多項式不変量を保存できるか？
• RQ2 非線形方程式の反復解法を回避しつつも、明示的であることを保ったまま、射影に基づく不変量保存アプローチをどのように簡略化できるか？
• RQ3 提案手法における不変量誤差の収束速度は何か？ また、高次精度スキームでは機械精度に達するか？
• RQ4 EIP法は、標準的な陽的解法および完全に陰的解法と比較して、長時間シミュレーションにおける性能と保存性に優れているか？
• RQ5 ボーズ＝アインシュタイン凝縮のような、複数の保存量を有する高次元PDEへ、この手法を効果的に拡張できるか？

主な発見

• EIP法は、元の明示的ルンゲ・クッタ法と同一の精度順序を維持し、不変量誤差が $\mathcal{O}(h^{2(p+1)})$ のオーダーに抑えられる。
• 十分に高次の方法（例：4次）では、不変量誤差が機械精度に達し、実用的に不変量を保存する。
• ボーズ＝アインシュタイン凝縮におけるリング状および3次元渦運動のODEおよびPDEに対する数値実験により、質量およびエネルギーの長期間にわたる機械精度の保存が確認された。
• 4次EIP法は、長時間シミュレーションにおいて、標準的な4次ルンゲ・クッタ法よりも質的流れの性質をよりよく保存する。
• 明示的構造と線形方程式の解法要求により、完全に陰的解法よりも実用的に効率的である。
• 2次元および3次元の動的系において、複数の不変量および高次元問題に対し、優れた保存性を維持し、一貫した誤差トレンドが観察された。