[論文レビュー] Relaxed Multirate Infinitesimal Step Methods for Initial-Value Problems

Jean Sexton, Daniel R. Reynolds|arXiv (Cornell University)|Aug 10, 2018

Numerical methods for differential equations参考文献 44被引用数 4

ひとこと要約

本稿では、一般化加法ルンゲ・クッタ（GARK）理論を用いて、再帰的フラックススプリット multirate（RFSMR）フレームワークを拡張することで、第四次精度を達成する新しい multirate 時間積分法—リラクスト multirate 無限小ステップ（RMIS）法—を提案する。本手法は、サブサイクリング、時間スケールの可変分離、埋め込み誤差推定を可能にするとともに、線形安定性を損なわず、明著な速い・遅いダイナミクスを有する multirate ODE システムにおいて高次精度を実現する。

ABSTRACT

This work focuses on the construction of a new class of fourth-order accurate methods for multirate time evolution of systems of ordinary differential equations. We base our work on the Recursive Flux Splitting Multirate (RFSMR) version of the Multirate Infinitesimal Step (MIS) methods and use recent theoretical developments for Generalized Additive Runge-Kutta methods to propose our higher-order Relaxed Multirate Infinitesimal Step extensions. The resulting framework supports a range of attractive properties for multirate methods, including telescopic extensions, subcycling, embeddings for temporal error estimation, and support for changes to the fast/slow time-scale separation between steps, without requiring any sacrifices in linear stability. In addition to providing rigorous theoretical developments for these new methods, we provide numerical tests demonstrating convergence and efficiency on a suite of multirate test problems.

研究の動機と目的

• 広範な時間スケールを有するマルチフィジックスシミュレーションにおいて、高次精度で安定的かつ効率的な時間積分法の必要性に対応すること。
• 時間スケールの分離が変化する状況において、安定性や精度を犠牲にする既存の multirate 法の限界を克服すること。
• 速い成分の明示的および陰的積分を両立させつつ第四次精度を維持するフレームワークの開発。
• 時間誤差推定、トーレスコピック拡張、統合中の時間スケール分離比の動的調整といった実用的機能の実装。
• ステップ間で速い・遅い時間スケール比が変化しても線形安定性が保たれることの保証。

提案手法

• 速い成分と遅い成分を明確に分離した加法的分割 ODE として multirate 問題を定式化：y' = f^f(y) + f^s(y)。
• GARK（一般化加法ルンゲ・クッタ）理論を用いて、multirate スキームの次数条件および安定性特性を導出する。
• 標準的な MIS 法の厳密な結合制約を緩和することで、より高い次数精度を達成する RMIS 法を構築する。
• 二段階構造を採用：遅い成分に対して外側のルンゲ・クッタ法を、速い成分に対して内側の方法を適用し、誤差推定用の埋め込み段階を設ける。
• 1回の遅いステップあたり m 回の速いステップを実行し、時間刻み比 m ≈ dt_s / dt_f を用いてサブサイクリングを実現する。
• 適切な係数設計により、元となるベース法の安定性特性を保ちながら線形安定性を確保する。

実験結果

リサーチクエスチョン

• RQ1 時間スケール分離が変化する状況下でも線形安定性を維持できる第四次精度 multirate 法を構築することは可能か？
• RQ2 安定性や精度を損なわず、サブサイクリングと時間誤差推定を multirate スキームに統合することは可能か？
• RQ3 GARK 理論は、緩い結合制約を許容する高次 multirate 法の実現にどのように寄与するか？
• RQ4 統合中に速い・遅い時間スケール比を動的に変更しても、収束性と安定性を保てるか？
• RQ5 標準的な multirate 法や単一時間刻み法と比較して、新規手法が multirate ODE 問題において実現する実用的効率向上はどの程度か？

主な発見

• 提案されたリラクスト multirate 無限小ステップ（RMIS）法は、明確に分離された速い・遅い成分を有する multirate ODE システムにおいて第四次精度を達成する。
• 本手法は、サブサイクリング、トーレスコピック拡張、埋め込み誤差推定を可能にし、線形安定性を損なわず、適応的時間刻みの実現を可能にする。
• フレームワークは、ステップ間での時間スケール分離比の動的調整を可能とし、複雑なシミュレーションにおける柔軟性を高める。
• 数値実験により、最適な収束率が確認され、単一時間刻み法と比較して multirate テスト問題において顕著な効率向上が示された。
• 理論的分析により、緩い結合仮定下でもすべての必要な次数条件および線形安定性制約を満たすことが確認された。
• 本手法は、速い成分の明示的および陰的積分の両方を対応可能とし、剛性および非剛性システムへの応用範囲を広げる。