[論文レビュー] Stochastic Variational Partitioned Runge-Kutta Integrators for Constrained Systems

Nawaf Bou‐Rabee, Houman Owhadi|ArXiv.org|Sep 14, 2007

Numerical methods for differential equations参考文献 19被引用数 12

ひとこと要約

本稿では、ホロノミック制約を受けるホワイトノイズの作用を受ける力学系に対して、確率的変分分割ルンゲ・クッタ（VPRK）積分法を導入する。変分積分法を確率的設定に拡張し、制約を伴う系に対し、第一および第二次の強収束性を持つRATTLE型積分法を構築した。その収束次数は後退誤差解析とVanden-EijndenとCicotti（2006）の手法により証明された。

ABSTRACT

Stochastic variational integrators for constrained, stochastic mechanical systems are developed in this paper. The main results of the paper are twofold: an equivalence is established between a stochastic Hamilton-Pontryagin (HP) principle in generalized coordinates and constrained coordinates via Lagrange multipliers, and variational partitioned Runge-Kutta (VPRK) integrators are extended to this class of systems. Among these integrators are first and second-order strongly convergent RATTLE-type integrators. We prove order of accuracy of the methods provided. The paper also reviews the deterministic treatment of VPRK integrators from the HP viewpoint.

研究の動機と目的

• ホロノミック制約を伴う確率的力学系の構造保存型積分法を開発する。標準的手法では確率的摂動を正しく捉えられない。
• 一般化座標および制約付き座標の両方で確率的ハミルトニアン＝ポントリャーギン原理を定式化することで、変分積分法理論を確率的設定に拡張する。
• 確率的摂動下でもシンプレクティック性と運動量を保存する第一および第二次の強収束型RATTLE型積分法を構築する。
• Vanden-EijndenとCicotti（2006）の手法を用いて、提案された積分法の強収束次数を厳密に証明し、信頼性の高い長時間シミュレーションを可能にする。
• 制約付き確率的系におけるエネルギー挙動、自己相関、経験的分布などの統計的性質の正確な計算を可能にする。

提案手法

• 一般化座標における確率的ハミルトニアン＝ポントリャーギン原理を定式化し、ラグランジュ乗数を用いた制約付き定式化と等価であることを示す。
• 離散的変分原理を用いて確率的変分積分法を導出し、シンプレクティック性および保存量に対する離散的ノイターの定理を保証する。
• 速度と位置の更新を分割し、制約多様体上への射影を施すことで、RATTLE型積分法を構築する。
• 離散的変分積分法に2次精度の分割ルンゲ・クッタ法を適用し、ラグランジュ乗数と運動量更新の整合性を確保する。
• 後退誤差解析とVanden-EijndenとCicotti（2006）の手法を用いて強収束次数を証明し、O(h²)までの誤差展開を含む。
• 非制約更新と射影ステップの両方が、平均二乗意味でシンプレクティック形式を保存することを示すことで、離散的フローマップの平均二乗シンプレクティック性を検証する。

実験結果

リサーチクエスチョン

• RQ1 ホロノミック制約を伴う確率的力学系に対して、変分積分法を拡張し、シンプレクティック性と運動量を保存することは可能か？
• RQ2 制約付き系に対する確率的VPRK積分法の強収束次数は何か？ そして、それを厳密に証明する方法は？
• RQ3 提案された積分法は、標準的な射影ベースまたは高次精度法と比較して、確率的エネルギーおよび運動量摂動をどれほど正確に捉えられるか？
• RQ4 SVIの構造保存性（正確なエネルギー挙動および統計的精度）は、制約とノイズの存在下でも維持可能か？

主な発見

• 第一次の強収束性を持ち、シンプレクティック性を保つRATTLE型積分法が、制約付き確率的力学系に対して構築された。強収束次数はO(h²)誤差展開により証明された。
• 積分法の離散的フローマップは平均二乗シンプレクティックであるため、確率的摂動下でも長期間にわたりシンプレクティック構造が保存される。
• 一般座標と制約付き座標における確率的ハミルトニアン＝ポントリャーギン原理の等価性が厳密に確立され、変分原理の一貫性のある定式化が可能になった。
• 確率的摂動を受ける球形振り子に対する数値実験により、強収束次数が1.0、弱収束次数が1.5であることが確認され、理論的予測と整合的であった。
• 積分法はエネルギーおよび運動量における確率的摂動を正しく捉えており、大きな時間刻みや小さなノイズ振幅でない場合でも、人工的なドリフトを生じさせない。
• この手法は、確率的効果の力学的性質を保存する。すなわち、エネルギーおよび運動量の変化は、すべてノイズに起因し、数値的アーチファクトによるものではない。