[論文レビュー] High Performance Computing with a Conservative Spectral Boltzmann Solver

Irene M. Gamba, J. Haack|arXiv (Cornell University)|Nov 1, 2012

Gas Dynamics and Kinetic Theory参考文献 8被引用数 8

ひとこと要約

本稿では、速度空間におけるフーリエ基底の畳み込みを用いた、第二階層の精度と保存性を備えたスペクトル法を、ボルツマン方程式に適用する。この手法により、非平衡状態の気体流れの高精度なシミュレーションが可能となり、スパコン上で近似的な線形スケーリングを達成し、古典的流体力学では記述できない急激な壁加熱に起因する衝撃波形成問題を解消できる。

ABSTRACT

We present new results building on the conservative deterministic spectral method for the space inhomogeneous Boltzmann equation developed by Gamba and Tharkabhushaman. This approach is a two-step process that acts on the weak form of the Boltzmann equation, and uses the machinery of the Fourier transform to reformulate the collisional integral into a weighted convolution in Fourier space. A constrained optimization problem is solved to preserve the mass, momentum, and energy of the resulting distribution. We extend this method to second order accuracy in space and time, and explore how to leverage the structure of the collisional formulation for high performance computing environments. The locality in space of the collisional term provides a straightforward memory decomposition, and we perform some initial scaling tests on high performance computing resources. We also use the improved computational power of this method to investigate a boundary-layer generated shock problem that cannot be described by classical hydrodynamics.

研究の動機と目的

• 質量、運動量、エネルギーの保存を保ちながら、高精度に非平衡状態の希薄ガス流れをシミュレートする課題に対処すること。
• ボルツマン方程式における五次元の衝突積分の計算ボトルネックを、スペクトル法によって克服すること。
• 非一様グリッドを用いて、空間的・時間的第二階層の精度を持つ保存的スペクトル法を拡張すること。
• スパコン上での大規模シミュレーションを可能にするスケーラブルな高性能コンピューティング（HPC）実装を実現すること。
• 古典的流体力学の適用範囲外にある境界層衝撃問題を調査すること。

提案手法

• 方程式の弱形式を用いて、ボルツマン衝突積分をフーリエ空間における重み付き畳み込みに再定式化する。
• L2最適化問題を用いて、質量・運動量・エネルギーの保存を、不変多様体への衝突出力の射影によって強制する。
• 畳み込みを効率的に計算するために高速フーリエ変換（FFT）を用い、構造化された重み処理により、O(N^{2d})の全操作を回避する。
• 空間的精度を向上させるために、非一様グリッドを用いた第二階層の有限体積法を実装する。
• 空間的データ交換のため、ゴーストセル法を用いて計算領域を複数プロセッサに分割する。
• メモリ使用量が制限される状況では、六次元重み配列 Ĝ(ζ,ξ) を事前に計算・保存するか、必要に応じてオンザフライで計算することにより、メモリアクセスを最適化する。

実験結果

リサーチクエスチョン

• RQ1 ボルツマン方程式に対する保存的スペクトル法が、空間的・時間的第二階層の精度を保ちつつ保存則を満たすことは可能か？
• RQ2 大規模な非平衡希薄ガス力学的問題に適用した場合、高性能コンピューティングシステムにおけるスケーリング特性はいかほどか？
• RQ3 古典的流体力学の有効範囲外にある急激な壁温度変化に起因する衝撃波形成を、本手法は解像できるか？
• RQ4 メモリ制限のあるHPC環境において、畳み込み重みを事前に計算するか、オンザフライで計算するかの性能的トレードオフは何か？
• RQ5 不連続な分布関数が境界付近に存在する場合、スペクトル法はギブス現象を回避できるか？

主な発見

• 本手法は、Lonestarスパコン上で近似的な線形スケーリングを達成し、1ステップの壁時刻を1ノードで203秒から64ノードで7.7秒に短縮した。
• ノード数を4から64に倍増するごとに計算時間が半減しており、強力な並列効率と最小限の通信オーバーヘッドを示している。
• 壁温度の急激な上昇に起因する衝撃波形成を、バルク領域では平均自由行程あたり2グリッドポイント、壁付近では8グリッドポイントで正確に捉えた。
• 壁付近の分布関数に不連続性が存在しても、ギブス現象は観測されず、畳み込み重みによる安定化効果が示唆された。
• 極めて非平衡的な流れに対しても、安定的かつ高精度に保たれ、モンテカルロ法に比べて急勾配や不連続性を優れて解像できる。
• 大規模なNにおいて、六次元重み配列 Ĝ(ζ,ξ) の保存はメモリ制限に陥りやすくなるが、現代のHPCシステムでは、オンザフライ計算の方が保存よりも効率的である可能性がある。