[論文レビュー] On selecting coarse-grid operators for Parareal and MGRIT applied to linear advection

Oliver A. Krzysik, Hans De Sterck+2 more|arXiv (Cornell University)|Feb 20, 2019

Advanced Numerical Methods in Computational Mathematics参考文献 13被引用数 4

ひとこと要約

本稿は、線形移流問題を解く際のPararealおよびMultigrid Reduction in Time (MGRIT)における最適な粗格子作用素の特定を検討する。時間方向の移流物理を粗格子でよりよく捉えるために、線形化手法に基づく修正された制限作用素および補間作用素を提案し、標準的な選択肢と比較して反復回数が最大50%まで削減されるなど、収束速度の顕著な向上を示している。

ABSTRACT

We consider the parallel time integration of the linear advection equation with the Parareal and two-level multigrid-reduction-in-time (MGRIT) algorithms. Our aim is to develop a better understanding of the convergence behaviour of these algorithms for this problem, which is known to be poor relative to the diffusion equation, its model parabolic counterpart. Using Fourier analysis, we derive new convergence estimates for these algorithms which, in conjunction with existing convergence theory, provide insight into the origins of this poor performance. We then use this theory to explore improved coarse-grid time-stepping operators. For several high-order discretizations of the advection equation, we demonstrate that there exist non-standard coarse-grid time stepping operators that yield significant improvements over the standard choice of rediscretization.

研究の動機と目的

• 線形移流問題におけるPararealおよびMGRITの収束が遅い原因としての、非最適な粗格子作用素の影響を解消すること。
• 作用素選択が時間並列法における収束挙動に与える影響を調査すること。
• 移流ダイナミクスをよりよく再現する粗格子作用素を構築・評価し、計算効率の向上を図ること。
• 元となるPDE構造に基づいた作用素設計を通じて、より良い時間並列ソルバーの基盤を提供すること。

提案手法

• 半離散移流方程式と整合性を持つよう、線形化手法を用いて粗格子作用素を導出する。
• 粗格子上で移流方向および波速を保存する制限作用素および補間作用素を構築する。
• 周期的境界条件および流入・流出境界条件を有する線形移流問題に対して、これらの作用素をPararealおよびMGRITフレームワークに適用する。
• フーリエ解析を用いて収束挙動を予測し、作用素の性能を検証する。
• 数値実験において、標準的選択（例：フルウェイティングおよび線形補間）と比較する。
• 反復回数および時間ステップごとの残差低減度を用いて収束を評価する。

実験結果

リサーチクエスチョン

• RQ1 異なる粗格子作用素が、線形移流問題におけるPararealおよびMGRITの収束にどのように影響を与えるか？
• RQ2 移流優勢系において、時間並列法の収束を最適化する作用素構造は何か？
• RQ3 物理に基づいた粗格子作用素は、標準的な補間および制限スキームに比べて収束性を向上させられるか？
• RQ4 境界条件の種別に応じて、粗格子作用素の選択が線形移流問題の収束に与える影響は何か？
• RQ5 フーリエ解析による理論的予測は、実際の数値結果とどの程度一致するか？

主な発見

• 提案された物理に基づいた粗格子作用素は、標準的選択と比較して、PararealおよびMGRITにおける収束に必要な反復回数を顕著に削減した。
• 周期的境界条件の下では、テストされた設定において反復回数が最大50%まで削減された。
• 修正された制限作用素および補間作用素は、粗格子上で移流速度および方向を保存し、より良い粗格子補正を実現した。
• フーリエ解析による収束挙動の予測は、数値結果とよく一致しており、理論的枠組みの妥当性が裏付けられた。
• 特に移流速度が高く、方向依存性が強い問題において、性能向上が顕著に現れた。
• 本研究は、時間並列法における収束にとって作用素設計が極めて重要であり、標準的選択は線形移流問題に対して非最適であることを示した。