[論文レビュー] A new symmetric linearly implicit exponential integrator preserving polynomial invariants or Lyapunov functions for conservative or dissipative systems

Lu Li|arXiv (Cornell University)|Apr 25, 2021

Numerical methods for differential equations参考文献 30被引用数 9

ひとこと要約

本稿では、保存的・散逸的剛性系に対して多項式不変量またはリャプノフ関数を保存する対称的で線形的陰解法の指数積分法を提案する。極性化離散勾配と指数積分法を組み合わせることで、完全陰解法と比較して著しく低い計算コストでエネルギー保存と優れた長時間安定性を達成する。

ABSTRACT

We present a new linearly implicit exponential integrator that preserves the polynomial first integrals or Lyapunov functions for the conservative and dissipative stiff equations, respectively. The method is tested by both oscillated ordinary differential equations and partial differential equations, e.g., an averaged system in wind-induced oscillation, the Fermi-Pasta-Ulam systems, and the polynomial pendulum oscillators. The numerical simulations confirm the conservative properties of the proposed method and demonstrate its good behavior in superior running speed when compared with fully implicit schemes for long-time simulations.

研究の動機と目的

• 多項式エネルギー関数を有する剛性保存的・散逸的系の構造保存型積分法の開発。
• 完全陰解法の限界を克服し、計算コストを低減する線形的陰解法を導入。
• 幾何積分法の技術を用いて対称性と長時間安定性を確保。
• 極性化離散勾配と指数時刻積分を用いてエネルギーまたはリャプノフ関数の正確な保存。
• 完全陰解法および陽解法と比較して、長時間シミュレーションにおいて優れた性能を示すこと。

提案手法

• 変数定数の公式を用いた指数時刻ステップを用いて、対称的線形的陰解法の指数積分法を構築。
• 非線形性を複数時刻にわたり分割するために極性化離散勾配を適用し、線形的陰解法を可能にする。
• 多項式ポテンシャルエネルギー関数（例：三次、四次）の二次的極性化を用いて、エネルギーまたはリャプノフ関数を保存する離散勾配を定義。
• yn+2 = exp(2hJM)yn + 2hφ(2hJM)J∇¯U(yn, yn+1, yn+2) として手法を定式化し、∇¯U は極性化離散勾配である。
• 多点極性化と離散勾配条件を用いて、高次多項式ポテンシャルへと手法を一般化。
• 一貫した平均化と勾配再構築を用いて2h時間ステップ上で手法を構築し、対称性を保証。

実験結果

リサーチクエスチョン

• RQ1 線形的陰解法の指数積分法は、剛性保存的・散逸的系において多項式不変量またはリャプノフ関数を保存できるか？
• RQ2 長時間シミュレーションにおいて、提案手法は完全陰解法および陽解法と比較して精度と計算コストで優れているか？
• RQ3 積分法における対称性は、線形誤差増大および不変量の近似保存といった、より優れた長時間挙動をもたらすか？
• RQ4 有界な自由エネルギーまたは追加のスカラー変数を必要とせずに、高次多項式ポテンシャルへと一般化可能か？
• RQ5 多項式エネルギー関数を有する系に対して、手法の精度の順序および収束挙動は何か？

主な発見

• 提案されたLIEEP手法は、保存的系（例：α-Fermi–Pasta–Ulam系、多項式振り子振動子）において極性化エネルギーを正確に保存する。
• 散逸的系では、リャプノフ関数の単調減少が物理的散逸挙動と整合的であることを確認した。
• 時間方向に2次精度であることが、h = 1/2^i（i = 1 から 5）の収束プロットにより確認された。
• 数値結果から、LIEEP手法は完全陰解法のEAVF手法と比較して著しく低い計算コストを達成しており、特に長時間シミュレーションで顕著である。
• 適切な極性化エネルギー定式化を用いることで、特定の系において超収束挙動を示した。
• h = 0.3 を用いた多項式振り子振動子の数値解は、[0, 1000] の区間で真の力学をよく再現し、円筒型位相空間構造を保存した。