On Runge-Kutta convolution quadrature based fractional variational integrators

技術要約：ルンゲ・クッタ型畳み込み量子化に基づく分数階変分積分器

問題提起
分数階微分方程式は、粘弾性材料や異常拡散などの、長期記憶を持つシステムのモデリングに不可欠である。分数階微積分は、変分原理を散逸系（分数階減衰）へと拡張するが、これらのシステムに対して高次数値積分器を構築することは依然として困難である。先行研究[23]では、保存的なラグランジアンのガラーキン変分法と、後退微分畳み込み量子化（BDFCQ）を組み合わせた。しかし、この手法は「飽和」現象に直面する。すなわち、BDFCQ自体は高次の精度を持つポテンシャルがあるものの、結果として得られる分数階変分積分器（FVI）は2次の精度に制限されてしまう。この制限は、BDFCQの多ステップ構造が主要な節点値のみに作用し、保存的な作用のガラーキン離散化に内在する内部ステージ情報を活用できないことに起因する。したがって、分数階部分の離散化が、保存的な離散化の内部ステージ構造と構造的に互換性を持つような、分数階変分積分器が必要とされている。

手法
著者らは、ルンゲ・クッタ型畳み込み量子化（RKCQ）と高次ガラーキン変分法を組み合わせた、新しいクラスの分数階変分積分器（FVI）を提案している。その手法は、以下の3つの柱に基づいている。

1. 拡張変分フレームワーク： 制限付きハミルトンの原理[22]に従い、分数階減衰を組み込むために状態変数（$x$ および $y$）を二重化してシステムをモデル化する。作用汎関数は、保存的な部分（$L_{con}$）と分数的な部分（$L_{fra}$）に分割される。
2. 離散化戦略：
   • 保存的な部分： 次数 $s$ の多項式補間とロバット（Lobatto）求積法を用いた高次ガラーキンの変分法によって近似される。
   • 分数的な部分： ロバットIIIC法に基づくRKCQを用いて近似される。BDFCQとは異なり、RKCQは内部のルンゲ・クッタ・ステージ（$t_k + c_i h$）において畳み込み積分を評価する。これにより、分数微分近似が、ガラーキンによる保存的作用の離散化で使用される内部節点と完全に一致することが可能となる。
3. 方程式の導出： 離散制限付きハミルトンの原理を適用することで、分数階離散オイラー・ラグランジュ（FDEL）方程式を導出する。この導出において鍵となるのは、RKCQの設定において確立された新しい離散恒等式である。
   • 畳み込み重みの離散半群特性（補題 4.1）。
   • 非局所的な分数微分の変分フレームワーク内での扱いを可能にするための、非対称な離散部分積分公式（補題 4.2 および 4.3）。

主な貢献
本論文の主な貢献は以下の4点である。

1. 互換性のある積分器の構築： RKCQ（具体的にはロバットIIIC）のステージ構造が、ガラーキン離散化の内部節点と一致するように、FVIを構築した。この構造的互換性により、従来のBDFCQベースのアプローチで見られた精度の飽和を回避できる。
2. 理論的基礎： RKCQの設定内でFDEL方程式を導出するために必要な、必須の離散恒等式（半群特性および部分積分）を確立した。
3. 収束解析：
   • 理論的： $\beta = 2\alpha \in (0, 1]$ である二次形式のラグランジアンに対し、2ステージ・ロバットIIICスキームが2次の収束性を有することを厳密に証明した。
   • 数値的： 数値実験により、滑らかな問題に適用した場合、3ステージおよび4ステージのロバットIIICスキームが、それぞれ4次および6次の高次精度を達成することを実証した。
4. 中点スキーム： 中点法に基づく分数階積分器（MIDCQ）の簡潔な議論と実装についても述べており、これも2次の精度を達成する。

結果
提案されたスキームは、2次元強制調和振動子と分数階バグリー・トーヴィック（Bagley-Torvik）方程式の2つのモデルを用いてテストされた。

• 精度： 2ステージ・ロバットIIICスキームにおいて、数値結果は理論通りの2次収束を確認した。3ステージおよび4ステージのスキームについては、観測された収束次数はそれぞれ4次および6次であり、これは基礎となるガラーキン法の理論的期待値と一致している。
• エネルギー減衰： 提案された積分器は、分数階減衰システムの期待されるエネルギー減衰挙動を正確に再現した。長時間のシミュレーションにおいて、提案されたFVIは、古典的なRK4法と比較して優れたエネルギー減衰特性を維持した。
• 正則性への依存性： 本研究は、高次収束が解の正則性に依存することを強調した。不連続な強制力（例：ステップ関数）を持つ問題では、収束次数は1次に飽和したが、これは基礎となるRKCQ理論の正則性要件と一致している。一方、滑らかな強制力に対しては、フルな高次精度が実現された。
• BDFCQとの比較： RKCQのアプローチが、内部ステージ情報を完全に活用することで、BDFCQのアプローチにおける次数制限を克服することを結果が裏付けている。

意義と主張
本研究の主要な意義は、分数階微分の近似と保存的な離散化との間の構造的互換性にある。ロバットIIICに基づくRKCQを利用することで、手法はガラーキン離散ラグランジアンで使用されるものと同じステージノードを自然に取り込むことができる。これにより、保存的な項と分数的な項の両方を共通の格子上で評価することが可能となり、BDFCQのような多ステップ型のアプローチを悩ませる精度の飽和を回避できる。

著者らは、高次スキームにおける4次および6次の収束を数値的証拠が支持しているものの、これらの高次に対する厳密な理論的証明は依然として未解決の問題であると謙虚に述べている。また、現在の実装の計算コストは、畳み込み重みの履歴依存性により $O(N^2)$ であるが、将来的な作業として高速畳み込み技術を用いることで $O(N \log^2 N)$ に削減できる可能性があることを認めている。結論として、本アプローチは、分数階減衰ラグランジアン系に対する高次で構造保存的な積分器を構築するための堅牢なフレームワークを提供するものである。