[論文レビュー] Convergence Analysis of a Class of Massively Parallel Direction Splitting Algorithms for the Navier-Stokes Equations

Jean‐Luc Guermond, P Minev+1 more|arXiv (Cornell University)|Jan 19, 2011

Advanced Numerical Methods in Computational Mathematics参考文献 22被引用数 5

ひとこと要約

本稿では、非圧縮性ナビエ＝ストークス方程式に対して、運動量と圧力補正ステップの両方を一次元の解法に分解する、大規模並列で方向分割されたアルゴリズムを提示する。この手法により、線形計算量、無条件安定性、および非線形収束性を達成する。本手法の精度は2次元でO(τ)、3次元でO(τ^{3/2})であることが証明されており、数値結果により、高精度および最大1024プロセッサまで最適な弱スケーラビリティが確認されている。

ABSTRACT

We provide a convergence analysis for a new fractional time-stepping technique for the incompressible Navier-Stokes equations based on direction splitting. This new technique is of linear complexity, unconditionally stable and convergent, and suitable for massive parallelization.

研究の動機と目的

• 非圧縮性ナビエ＝ストークス方程式の時間積分法を、計算的に効率的かつ高度に並列化可能であるように開発すること。
• 従来の投影法における多次元ポアソン方程式の解法という計算ボトルネックを、方向分割によって一次元楕円型問題に置き換えることにより克服すること。
• 新アルゴリズムの厳密な収束解析を確立し、2次元および3次元における誤差推定を含めること。
• リーマン・ドライブン・キャビティおよびバックワード・フェーシング・ステップなどの複雑な流れの大型シミュレーションにおいて、安定性、精度、スケーラビリティを実証すること。
• 並列アーキテクチャを用いて、大規模グリッド上での一時的流れの効率的シミュレーションを可能にすること。

提案手法

• 運動量方程式および圧力補正ステップの両方に方向分割を適用し、それぞれを一次元楕円型問題の系列に簡略化する。
• 各時間ステップにおいて、速度のための1次元移流拡散問題と、圧力補正のための1次元ポアソン型問題の系列を解く。
• 時間離散化には分数ステップ法を用い、速度と圧力の進化を分離しながらも安定性を保持する。
• 空間離散化には、微分の中央差分を用いたMAC（マーカー・アンド・セル）有限差分法を採用する。
• MPIを用いた領域分割により並列実装を行い、各サブドメインで独立して1次元系を解く。
• サブドメイン間の解を結合するためにスチェュール補完（Schur complement）を用い、グローバル系の効率的並列解法を実現する。

実験結果

リサーチクエスチョン

• RQ1 運動量と圧力補正ステップの両方に方向分割を適用することで、線形計算量と大規模並列性を達成できるか？
• RQ2 得られたアルゴリズムは、非圧縮性ナビエ＝ストークス方程式に対して無条件安定かつ収束するか？
• RQ3 本手法の収束速度は2次元および3次元でそれぞれどのように異なり、空間次元に依存するか？
• RQ4 本手法は大規模3次元シミュレーションにおいても高精度とスケーラビリティを維持できるか？
• RQ5 ベンチマーク問題において、標準的な投影法と比較して性能と精度はどのように異なるか？

主な発見

• 誤差推定により、本手法は無条件安定かつ収束し、2次元ではO(τ)、3次元ではO(τ^{3/2})の収束率を示すことが証明されている。
• 数値実験により、実際の収束速度が2次元および3次元の両方でO(τ^{3/2})を上回ることが確認され、実際には高次の挙動を示していることが示された。
• 最大1024プロセッサまで最適な弱スケーラビリティを達成しており、最大2×10^9ノードのグリッドにおいても性能の著しい低下がないことが確認された。
• Re=1000およびRe=5000の3次元リーマン・ドライブン・キャビティのシミュレーションでは、t=4、8、12における速度場に4桁の精度が達成された。
• バックワード・フェーシング・ステップのシミュレーションでは、KimとMoin（2000）のベンチマークデータと非常に良好に一致し、誤差は1%未満であった。
• 分割されたスキームと分割しないスキームの速度プロファイルは、小数第4位までしか差がなく、高精度であることが確認された。