[論文レビュー] Symmetric parareal algorithms for Hamiltonian systems

Xing Dai, Bris, C. Le+2 more|arXiv (Cornell University)|Nov 29, 2010

Numerical methods for differential equations参考文献 32被引用数 10

ひとこと要約

この論文では、ハミルトニアン系に対する対称パラレルアルゴリズムを導入し、エネルギー保存則と軌道精度を長時間にわたり維持する。対称時間刻み法とエネルギー多様体への粗解像の幾何的射影を統合することで、1000プロセッサを用いて逐次解法にほぼ匹敵する精度を達成しながら66倍の高速化を実現。標準パラレルアルゴリズムに比べ、ハミルトニアン力学における安定性と幾何的忠実性が著しく向上する。

ABSTRACT

The parareal in time algorithm allows to efficiently use parallel computing for the simulation of time-dependent problems. It is based on a decomposition of the time interval into subintervals, and on a predictor-corrector strategy, where the propagations over each subinterval for the corrector stage are concurrently performed on the processors. In this article, we are concerned with the long time integration of Hamiltonian systems. Geometric, structure-preserving integrators are preferably employed for such systems because they show interesting numerical properties, in particular excellent preservation of the total energy of the system. Using a symmetrization procedure and/or a (possibly also symmetric) projection step, we introduce here several variants of the original plain parareal in time algorithm [Lions, Maday and Turinici 2001, Baffico, Bernard, Maday, Turinici and Zerah 2002, Bal and Maday 2002] that are better adapted to the Hamiltonian context. These variants are compatible with the geometric structure of the exact dynamics, and are easy to implement. Numerical tests on several model systems illustrate the remarkable properties of the proposed parareal integrators over long integration times. Some formal elements of understanding are also provided.

研究の動機と目的

• ハミルトニアン系に適用した際の標準パラレルアルゴリズムの長時間安定性とエネルギー保存則の欠如を是正すること。
• ハミルトニアン力学の幾何的構造、特にシンプレクティック性と時間反転対称性を尊重するパラレル変種を開発すること。
• 対称化と対称的射影ステップの組み合わせにより、長時間シミュレーションにおける収束性と不変量保存の向上を図ること。
• 提案されたアルゴリズムが、逐次的高精度ソルバーと同等の精度を維持しながら、ハミルトニアン問題において顕著な並列高速化を達成できることを示すこと。

提案手法

• 予測補正ステップの対称化により時間反転対称性と幾何的構造の改善を確保する、標準パラレルアルゴリズムの適合。
• 粗解像において各時間サブインターバルの終了時に、解を一定エネルギー多様体上に射影する対称的射影ステップの導入。
• 良好な幾何的性質（例：シンプレクティック積分法）を有する細かいソルバと、予測フェーズ用の粗く単純なソルバを用いる。
• 計算コストを低減するため、射影処理を粗解像でのみ実行し、サブインターバルあたり1〜2回の反復で十分である。
• 細かい補正をサブインターバルごとに並列に実行することで、アルゴリズムが並列性を維持し、対称性と幾何的整合性を保つ。
• 対称化と射影の組み合わせにより、標準パラレルアルゴリズムが抱える共鳴問題を効果的に解消し、エネルギー以外の不変量の安定化を実現する。

実験結果

リサーチクエスチョン

• RQ1 パラレルアルゴリズムの対称化は、ハミルトニアン系における幾何的構造保存に寄与するか？
• RQ2 粗解像においてエネルギー多様体上に対称的射影を適用することで、長時間安定性と精度が向上するか？
• RQ3 対称化と射影の組み合わせは、長時間シミュレーションにおける収束性と不変量保存にどのように影響するか？
• RQ4 新しい対称パラレルアルゴリズムでは、逐次的高精度ソルバーと同等の精度を維持しながら、どの程度の高速化が達成できるか？

主な発見

• 幾何的射影を施した対称パラレルアルゴリズムは、長時間積分において逐次的高精度ソルバーと同等の軌道精度を達成する。
• 1000プロセッサを用いた場合、約66倍の高速化が達成され、15回のパラレル反復で収束に至る。
• 非線形的対称射影に要する平均反復回数はサブインターバルあたりたった1.12回であり、計算コストの増加は顕著でない。
• アルゴリズムは全エネルギーの安定化に加え、他の不変量の安定化も実現し、長期シミュレーションの信頼性が向上する。
• 対称化と射影の組み合わせにより、ハミルトニアン系における標準パラレルアルゴリズムが抱える共鳴問題が効果的に解消される。
• モデル系に対する数値実験により、本手法のロバスト性と長期積分における優れた幾何的性質が確認された。