[論文レビュー] Efficient supersonic flows through high-order guided equilibrium with lattice Boltzmann

Jonas Lätt, Christophe Coreixas+2 more|arXiv (Cornell University)|Oct 29, 2019

Lattice Boltzmann Simulation Studies参考文献 54被引用数 5

ひとこと要約

本稿では、3次元で39個の離散速度のみを用いて、高次元のラティスボルツマン法（LBM）を提案する。この手法は13モーメント平衡分布関数を用い、超音速・圧縮性流れを高精度にシミュレート可能である。Knudsen数に依存する緩和時間の解析的計算と、局所的粘性係数を調整するキネティックセンサーを組み合わせることで、安定した高レイノルズ数の超音速流れのシミュレーションが実現された。NACA0012の翼型をMa=1.5、Re=10⁴でシミュレートした結果、GPUハードウェアを用いることでCPU実装と比較して20倍の高速化を達成した。

ABSTRACT

A double-distribution-function based lattice Boltzmann method (DDF-LBM) is proposed for the simulation of polyatomic gases in the supersonic regime. The model relies on an extended equilibrium state that is constructed to reproduce the first 13 moments of the Maxwell-Boltzmann distribution exactly. This extends the validity of the standard 5-constraint (mass, momentum and energy) approach and leads to the correct simulation of thermal, compressible flows with only 39 discrete velocities in 3D. The stability of this BGK-LBM is reinforced by relying on Knudsen-number-dependent relaxation times that are computed analytically. Hence, high-Reynolds number, supersonic flows can be simulated in an efficient and elegant manner. While the 1D Riemann problem shows the ability of the proposed approach to handle discontinuities in the zero-viscosity limit, the simulation of the flow past a NACA0012 airfoil (Mach number $\mathrm{Ma}=1.5$, Reynolds number $\mathrm{Re=10^4}$) confirms the excellent behavior of this model in a low-viscosity and supersonic regime. The proposed model is substantially more efficient than the previous 5-moment D3Q343 DDF-LBM and opens up a whole new world of compressible flow applications that can be realistically tackled with a purely LB approach.

研究の動機と目的

• 標準的な5モーメントモデルの制限を超えて、超音速・圧縮性流れを安定かつ効率的にシミュレートできるラティスボルツマン法の開発。
• 熱的および圧縮性効果を正確に捉えることで、LBMの有効範囲を非圧縮および弱圧縮領域から拡大すること。
• 3次元で343個の離散速度を必要とする従来の高速度DDF-LBM方式と比較して、計算コストを低減すること。
• 純粋なラティスボルツマン法を用いて、実用的で産業スケールの高マッハ数流れシミュレーションを可能にすること。

提案手法

• 本手法は、Maxwell-Boltzmann分布の最初の13モーメントを正確に再現する拡張された平衡関数を備えた二重分布関数ラティスボルツマンモデル（DDF-LBM）を採用している。
• 3次元で39速度格子を用い、従来の343速度モデルと比較して自由度を顕著に低減しながらも、圧縮性・超音速領域での精度を維持している。
• Knudsen数に依存する緩和時間を解析的に計算することで、高レイノルズ数流れにおける数値的安定性を向上させている。
• 平衡からの逸脱に基づいて局所的粘性係数を増加させるキネティックセンサーを導入し、衝撃波や不連続領域における安定性を向上させている。
• アルゴリズムはCPUおよびGPUの両方で実装されており、GPU版ではマルチコアCPU実行と比較して20倍の高速化を達成している。
• 1次元Riemann問題と、Ma=1.5、Re=10⁴における2次元NACA0012翼型シミュレーションを通じて、モデルの妥当性が検証された。

実験結果

リサーチクエスチョン

• RQ1 13モーメント平衡関数に基づくラティスボルツマン法は、従来のモデルと比較してより少ない離散速度で、超音速・圧縮性流れを正確にシミュレート可能か？
• RQ2 Knudsen数に依存する緩和時間を用いることで、高レイノルズ数の超音速流れにおける安定性はどのように向上するか？
• RQ3 平衡からの逸脱に基づくキネティックセンサーは、衝撃波捕捉シナリオにおける数値的安定性をどの程度向上させるか？
• RQ4 提案されたLBMは、産業スケールの圧縮性流れシミュレーションに実用的な計算性能を提供できるか？
• RQ5 GPUによる加速を経て、13モーメントアプローチは実世界の応用に十分にスケーラブルかつ効率的か？

主な発見

• 提案された13モーメントLBMは、3次元で39個の離散速度のみを用いても、超音速流れの正確なシミュレーションを実現した。これは、従来のDDF-LBM方式が要請する343個の速度と比較して顕著な削減である。
• 零粘性極限における1次元Riemann問題のシミュレーションでは、不連続性を正確に扱えることが確認され、強力な衝撃波の取り扱い能力が裏付けられた。
• Ma=1.5、Re=10⁴におけるNACA0012翼型シミュレーションは、低粘性・超音速領域でも優れた安定性と精度を示した。
• GPU実装により、マルチコアCPU実行と比較して20倍の性能向上が達成され、NACA0012の全シミュレーションが4時間未満で実行可能となった。
• 標準的な多項式平衡関数LBMと比較して、計算コストはわずか5倍のオーバーヘッドに抑えられつつ、圧縮性流れに対する物理的精度ははるかに優れていた。
• 平衡からの逸脱に基づくキネティックセンサーは、精度を損なうことなく、動的かつ局所的な安定化機構を効果的に提供した。