[論文レビュー] Accuracy and Performance Evaluation of Low Density Internal and External Flow Predictions using CFD and DSMC

Surya Kiran Peravali, Vahid Jafari+5 more|arXiv (Cornell University)|Jan 1, 2024

Aerosol Filtration and Electrostatic Precipitation被引用数 1

ひとこと要約

本研究では、OpenFOAM と SPARTA を用いて、Knudsen 数が変化する低密度の内部および外部流れにおける CFD および DSMC シミュレーションの精度と計算性能を評価する。ハイブリッド CFD-DSMC 法は遷移領域において実用的で効率的な代替手法を提供し、最適なパrameter選定により統計誤差とエネルギー消費を顕著に低減しながら、実験データと高い一貫性を維持することができる。

ABSTRACT

The Direct Simulation Monte Carlo (DSMC) method was widely used to simulate low density gas flows with large Knudsen numbers. However, DSMC encounters limitations in the regime of lower Knudsen numbers (Kn<0.1). In such cases, approaches from classical computational fluid dynamics (CFD) relying on the continuum assumption are preferred, offering accurate solutions at acceptable computational costs. In experiments aimed at imaging aerosolized nanoparticles in vacuo a wide range of Knudsen numbers occur, which motivated the present study on the analysis of the advantages and drawbacks of DSMC and CFD simulations of rarefied flows in terms of accuracy and computational effort. Furthermore, the potential of hybrid methods is evaluated. For this purpose, DSMC and CFD simulations of the flow inside a convergent-divergent nozzle (internal expanding flow) and the flow around a conical body (external shock generating flow) were carried out. CFD simulations utilize the software OpenFOAM and the DSMC solution is obtained using the software SPARTA. The results of these simulation techniques are evaluated by comparing them with experimental data (1), evaluating the time-to-solution (2) and the energy consumption (3), and assessing the feasibility of hybrid CFD-DSMC approaches (4). Keywords: DSMC; SPARTA; Continuum assumption; Transition regime; Rarefied flow; high-performance computing

研究の動機と目的

• ナノ粒子イメージング実験に関連する Knudsen 数（Kn ≈ 0.01–10）の広い範囲において、CFD および DSMC 法の精度と計算効率を評価すること。
• メッシュサイズ、粒子数、時間ステップ、衝突モデル、境界条件といった主要なシミュレーションパrameter が解の精度とパフォーマンスに与える影響を調査すること。
• 中間 Knudsen 数領域における精度と計算コストのバランスを最適化するハイブリッド CFD-DSMC アプローチの最適構成を特定すること。
• CFD および DSMC シミュレーションのソリューション到達時間とエネルギー消費を定量的に評価し、ハイパフォーマンスコンピューティングシステムにおける持続可能性とスケーラビリティを検証すること。
• 物理的妥当性を保証するため、内部（ノズル）および外部（円錐体）の流れ設定において、シミュレーション結果を実験データと照合すること。

提案手法

• CFD には OpenFOAM、DSMC には SPARTA を用いて、収束・発散ノズル内の内部流れおよび円錐体周囲の外部流れの 2D および 3D シミュレーションを実施。
• CFD シミュレーションでは、スリップ境界条件を伴うナビエ–ストークス方程式と、不完全熱吸収を伴う拡散反射の Cercignani-Lampis-Lovelace（CLL）モデルを適用。
• 希薄領域では、ボルツマン方程式を解くために、可変硬球（VHS）衝突モデルと確率的衝突サンプリングを用いた直接シミュレーション・モンテカルロ（DSMC）法を適用。
• グリッド独立性の検証を実施し、CFD には構造的でボディに適合した非構造格子を、DSMC にはセルベースの粒子分布と適応的時間ステップを用いた。
• 統計誤差を最小限に抑えるために、セルあたりの粒子数（Nc）と時間ステップ数（NT）に対する感度分析を実施。
• HPC クラスタ HSUper および DESY の Maxwell を用いて、強スケーリングおよび弱スケーリングのパフォーマンスとエネルギー消費を評価。SPARTA の動的ロードバランシングを活用して最適なパフォーマンスを達成。

実験結果

リサーチクエスチョン

• RQ1 Knudsen 数が変化する低密度の内部および外部流れにおいて、CFD と DSMC 法の精度はどのように比較されるか？
• RQ2 メッシュ解像度、セルあたりの粒子数、時間ステップ、衝突モデルといった主要パrameter が、解の精度と計算コストに与える影響は何か？
• RQ3 CFD と DSMC シミュレーションのソリューション到達時間とエネルギー消費はどのように異なり、希薄流れ問題においてどちらの手法がより高い計算効率を示すか？
• RQ4 ハイブリッド CFD-DSMC アプローチは、遷移領域（Kn ≈ 0.01–10）において、計算効率を損なわずに精度とパフォーマンスをどのように向上させ得るか？
• RQ5 ノズルおよび円錐体構成におけるマッハ数プロファイルと流れの物理的性質について、シミュレーション結果は実験データとどの程度一致するか？

主な発見

• DSMC シミュレーションは Kn > 0.1 の高 Knudsen 数領域で高い精度を達成した一方、CFD 法は Kn < 0.01 において信頼性の高い結果を示し、流体力学的極限における連続体仮定の妥当性を裏付けた。
• すべてのケースにおいて、定常緩和を伴う VHS 衝突モデルが、他のモデルと比較して優れた計算速度を示し、精度への影響は最小限であった。
• 統計誤差解析から、特に並進温度の予測において、時間ステップ数（NT）を増やすのではなく、セルあたりの粒子数（Nc）を増やすことで統計的フラクチュエーションをより効果的に低減できることが明らかになった。
• ハイブリッド CFD-DSMC アプローチは、DESY におけるアerosol化ナノ粒子を含む広範な Knudsen 数を有する分子ビーム実験のシミュレーションにおいて、強く実現可能であることが示された。
• パフォーマンスベンチマークから、SPARTA の動的ロードバランシングが強スケーリングおよび弱スケーリングの両面で最適化を実現し、DSMC シミュレーションにおけるソリューション到達時間の短縮とエネルギー消費の低減に顕著な効果を示した。
• エネルギー消費と計算コストは、DSMC において Kn−4 に比例して増加することが判明し、Kn < 0.05 の領域ではシミュレーションが極めて高コストとなることが示された。これにより、中間領域におけるハイブリッド手法の必要性が強化された。