[論文レビュー] Relaxed Multirate Infinitesimal Step Methods

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

Numerical methods for differential equations被引用数 3

ひとこと要約

本稿では、一般化加法的ルンゲ・クッタ理論を活用して、再帰的フラックススプリットリング多時間スケール法フレームワークを拡張することで、常微分方程式系に対する4次精度の新しいクラスの緩和型多時間スケール無限小ステップ（RMIS）法を提案する。本手法は、線形安定性を損なわず、サブサイクリング、適応的時間スケール分離、埋め込み誤差推定を可能にする。多時間スケール問題における数値的収束性および効率性のテストにより検証されている。

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 version of the Multirate Infinitesimal Step 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.

研究の動機と目的

• 多時間スケールの常微分方程式系のための高次精度時間積分法の開発を目的とする。
• 線形安定性を保ちつつ、先進的機能を可能にする緩和型定式化を用いて、多時間スケール無限小ステップフレームワークを拡張することを目的とする。
• 積分中に高速/低速時間スケールの分離を動的に変更しても安定性を損なわないようにすることを目的とする。
• 安定で高次精度のフレームワーク内にサブサイクリングと時間誤差推定の埋め込みを組み込むことを目的とする。
• 理論的裏付けがあり、多様な多時間スケール問題に適用可能な柔軟な多時間スケール時間積分法を提供することを目的とする。

提案手法

• 本手法は、再帰的フラックススプリットリング多時間スケール無限小ステップ定式化に基づき、一般化加法的ルンゲ・クッタ理論を用いて緩和段を導入することで拡張されている。
• 任意の時間スケール分離においても線形安定性を維持しながら高次精度を達成できる緩和係数を導入している。
• フレームワークはテレスコピック拡張をサポートしており、異なるサブサイクリング比を有する複数の時間ステップでの効率的統合を可能にしている。
• 埋め込みを組み込むことで、基本手法に追加コストをかけずに時間誤差推定が可能になっている。
• 本手法は、ステップ間で高速/低速時間スケール分離を動的に適応させることができ、非一様時間ステッピング戦略をサポートしている。
• 理論的展開により、すべての向上が線形安定性を損なわず実現されていることが保証されており、これは多時間スケール手法における重要な制約である。

実験結果

リサーチクエスチョン

• RQ1 時間スケールの動的適応を可能にしつつ、線形安定性を維持する4次精度の多時間スケール法を構築できるか？
• RQ2 緩和型多時間スケールフレームワーク内にサブサイクリングと誤差推定を組み込む方法は何か？ その際、安定性が損なわれないか？
• RQ3 緩和型定式化は、時間スケール分離の変動に柔軟に対応できる多時間スケール手法の柔軟性をどの程度向上させるか？
• RQ4 提案手法は、標準的な多時間スケールテスト問題において、収束特性と効率性をどのように示すか？
• RQ5 テレスコピック拡張と一般化加法的ルンゲ・クッタ構造を保持しつつ、高次精度を維持できるか？

主な発見

• 提案された緩和型多時間スケール無限小ステップ法は、すべてのテスト設定において4次精度を達成するとともに、線形安定性を維持している。
• 本手法は、複数段階の時間ステッピングを有するサブサイクリングを効果的にサポートしており、剛性および非剛性成分の効率的統合を可能にしている。
• 時間誤差推定は、緩和フレームワーク内に埋め込まれたペアを用いて達成されており、適応的時間ステッピングが可能になっている。
• フレームワークは、ステップ間で高速/低速時間スケール分離を動的に変更しても安定性の低下を示さない。
• 数値的テストにより、多時間スケールテスト問題において期待される4次収束率が確認され、精度と効率の両方が裏付けられている。
• 本手法は、多様な多時間スケールシナリオにおいても頑健な性能を示しており、理論的柔軟性と実用的有用性が検証された。