Discrete variational calculus for double-bracket dissipation

Bloch らによる論文「Discrete variational calculus for double-bracket dissipation」の詳細な技術的概要を以下に示す。

1. 問題定義

本論文は、ダブルブラケット散逸を受ける機械系を数値シミュレーションする際の課題に取り組んでいる。これらはエネルギーが軌道に沿って減少する散逸系であるが、特定の幾何学的構造、具体的には基礎となるリー代数の余随伴軌道が不変でなければならない。

• 文脈: 標準的な幾何学的積分器（シンプレクティック積分器等）は、保存系に対して幾何学的性質（シンプレクティシティ、運動量写像）を保存する点で優れている。しかし、散逸を加えると通常、これらの不変量が破壊される。
• 具体的な課題: 減衰器を備えた人工衛星、地球物理学的流体、プラズマ物理学などの系では、ダイナミクスはリー代数の双対空間（$\mathfrak{g}^*$）上で進化します。エネルギーが散逸する一方で、系はカシミール不変量のレベルセットである特定の余随伴軌道上に厳密に留まらなければならない。
• 既存手法の限界: 汎用積分器（ルンゲ・クッタ法等）は、高次精度であっても、余随伴軌道を厳密に保存できず、カシミール不変量のドリフトや、長期的な挙動の質的な誤差を引き起こすことが多い。

2. 手法

著者らは、強制オイラー・ポアンカレ系およびダブルブラケット散逸を伴うリー・ポアンカレ系に特化した離散変分積分器を開発した。

A. 理論的基盤

1. ダブルブラケット散逸: 本論文は、Bloch ら (1996) によって確立された枠組みを利用しており、散逸力は写像 $\phi: \mathfrak{g} \to \mathfrak{g}$ を通じて定義される。連続的な運動方程式は以下の通りである：
   $$\dot{\mu} = \text{ad}^*_{\frac{\delta H}{\delta \mu}} \mu + \text{ad}^*_{\phi(\frac{\delta H}{\delta \mu})} \mu$$
   この定式化により、ハミルトニアン $H$ が減少する（$\dot{H} \leq 0$）一方で、カシミール関数 $C$ は保存され（$\dot{C} = 0$）、軌道が余随伴軌道上に留まることが保証される。

2. 離散変分原理:

   • 著者らは、リー群 $G$ 上で定義された離散ラグランジアン $L_d: G \to \mathbb{R}$ を定義する。
   • 散逸写像 $\phi_d$ を取り入れた特定のバリエーション $\delta w_k$ を用いた離散一般化変分原理を導入する。
   • 離散作用は、以下の形式の変分に対して最小化される：
     $$\delta w_k = -\eta_k w_k + w_k (\text{Ad}_{\phi_d(w_k)} \eta_{k+1})$$
     ここで、$\eta_k$ は境界で消えるリー代数における変分である。

3. 射影写像:
   実用的なアルゴリズムを構築するために、著者らは射影写像 $\tau: \mathfrak{g} \to G$（リー代数からリー群への解析的局所微分同相写像、例えば指数写像またはケーリー写像）を使用する。

   • 離散運動量 $M_k$ は以下のように更新される：
     $$M_{k+1} = \text{Ad}^*_{w_k \phi_d(w_k)} M_k$$
   • この更新則は、群の増分（$w_k$）をリー代数の速度（$\xi_k$）に関連付けるために射影写像を用いて連続変分原理を離散化することで導出される。

B. アルゴリズムの構築

提案された**離散ダブルブラケット（DDB）**法は、以下の手順で行われる：

1. 離散運動量 $M_0 = \mu(t_0)$ を初期化する。
2. 各時間ステップ $k$ について：
   • $M_k = (d\tau^{-1}_{h\xi_k})^* \frac{\partial l}{\partial \xi}(\xi_k)$ となるリー代数要素 $\xi_k$ を解く。
   • 群の増分 $w_k = \tau(h\xi_k)$ を計算する。
   • 運動量を更新する：$M_{k+1} = \text{Ad}^*_{w_k \phi_d(w_k)} M_k$。
3. 写像 $\phi_d$ は、射影を用いて連続的な $\phi$ から構成される：$\phi_d(w_k) = \tau(h\phi(\xi_k))$。

3. 主な貢献

1. 余随伴軌道の厳密な保存: 主な貢献は、ダブルブラケット散逸の存在下で余随伴軌道を厳密に保存する幾何学的積分器であることである。カシミール不変量がドリフトする標準的な手法とは異なり、この手法は数値解が初期条件で定義された多様体上に厳密に留まることを保証する。
2. 変分的導出: この手法は離散変分原理から導出されており、変分積分器の安定性と構造保存特性を継承している。
3. 射影写像による柔軟性: この枠組みは、異なる射影写像（指数写像、ケーリー写像等）の使用を可能にし、異なる計算コストと精度次数（例えば、2 次対称スキーム等）を持つ積分器の構築を可能にする。
4. 剛体への適用: 理論は散逸剛体に明示的に適用され、$SO(3)$ に対する具体的な行列形式が提供されている。

4. 結果と数値シミュレーション

著者らは、DDB 法（指数写像およびケーリー写像の両方を使用）を以下と比較して、散逸剛体をテストケースとして手法を検証した：

• 4 次ルンゲ・クッタ法（RK4）
• 3 段階ロバト IIIC
• DOP853（8 次、参照として使用）

主要な知見:

• カシミールの保存: すべてのシミュレーションにおいて、DDB 法はカシミール不変量（角運動量の大きさ $\|M\|$）を機械精度で保存した。これに対し、RK4 およびロバト IIIC はカシミールに顕著なドリフトを示し、軌道が余随伴球面から外れる原因となった。
• 漸近挙動:
  • シナリオ A（孤立平衡点）: 異なる慣性モーメントを持つ剛体の場合、系は最大慣性軸（最小エネルギー）周りの回転に収束する。DDB 法は正しい極限点に収束したが、非幾何学的な手法はカシミールの誤差によりドリフトした。
  • シナリオ B（大円上の平衡点）: 対称な物体（2 つの等しい慣性モーメント）の場合、系は大円上の平衡点の集合に収束する。DDB 法はこの円上の点に正しく収束したが、他の手法は多様体から外れてドリフトした。
• 効率性: DDB 法（特にケーリー射影を用いたもの）は計算効率が良く、高次のロバト IIIC 法よりも高速であることが多く、かつ質的な精度が優れていた。
• 収束性: 手法は参照解への収束を示し、長期的なシミュレーションにおいて、DDB 法は基礎となる幾何学の保存により、最終運動量の誤差が著しく低かった。

5. 意義

この研究は、幾何学的数値積分における重要なギャップを埋めるものである：

• 物理的忠実性: 内部減衰器を備えた人工衛星やプラズマ閉じ込めなどの系では、角運動量（または同等のカシミール）の保存は物理的に必須である。標準的な積分器はこれを違反し、長期的なシミュレーション時間において非物理的な結果をもたらす。
• 堅牢性: 提案された積分器は、人工的なドリフトなしに系の漸近安定性を正しく捉える、長期的な積分に対して堅牢である。
• 一般化可能性: この枠組みは剛体に限定されない。著者らは、このアプローチをリー代数束や理想流体に拡張できることに言及しており、オイラー方程式の無限次元カシミール不変量を保存しながら、散逸を伴う複雑な流体力学をシミュレートする道を開いている。

結論として、本論文は、ダブルブラケット散逸を効果的に処理する厳密な変分に基づく幾何学的積分器を提示しており、散逸機械系の数値シミュレーションが連続的な物理に対して幾何学的に忠実であることを保証している。