[論文レビュー] Small Rarefaction, Large Consequences: Limits of Navier Stokes Turbulence Simulations

Songyan Tian, Lei Wu|arXiv (Cornell University)|Feb 9, 2026

Gas Dynamics and Kinetic Theory被引用数 0

ひとこと要約

弱い稀人人化効果が乱流ジェット衝突で局所的に支配的となり、直接ボルツマン計算によって表面せん断応力と熱流束のNavier–Stokes予測に大きな誤差を生じさせることを示す。

ABSTRACT

We conduct numerical simulations of rocket plume impingement on a lunar landing surface using two complementary frameworks: the Boltzmann equation, which naturally captures rarefied gas dynamics, and the Navier Stokes (NS) equations, the conventional workhorse for turbulent flow simulations. We show that subtle rarefaction effects, long considered negligible in turbulent regimes, can become locally dominant within shear layers where viscous stresses predicted by the NS constitutive relation undergo sign reversals. This phenomenon, which we term constitutive degeneracy, produces order-one relative errors in predicted surface shear stress and heat flux. Our results demonstrate that turbulence can expose hidden limits of NS equations with broad implications for high-speed aerodynamics and planetary exploration.

研究の動機と目的

• 航空宇宙シナリオに関連する乱流・不均一流における稀人人化効果の調査を動機づける。
• 局所的な稀人人化の存在下でのNavier–Stokes乱流モデル化の限界を評価する。
• 平衡を超える非平衡応力が表面荷重を支配する構成定数の退化を示す。
• ジェット衝突の設定で連続体（NS）予測と運動論（Boltzmann）解を比較する。
• 高速 aerodynamics および惑星探査応用への洞察を提供する。

提案手法

• 窒素排気に対して修正Rykov衝突モデルを用いた離散速度法でボルツマン方程式を解く。
• Explicitな領域分解なしに多スケール連続-稀薄動力学を捉える粗視化GSIS–SSTフレームワークを組み合わせる。
• SST閉塞を超える乱流と稀 People稀有人の相互作用を検証する時間分解ボルツマンシミュレーションを実施する。
• 全応力をNS（層流+SST）と Boltzmann非平衡成分に分解して組成関係を分析する。
• 局所Knudsen数と勾配長Knを評価し、NS応力が消滅または符号を反転する縮退領域を特定する。
• ボルツマンの直接的結果をNS-SST予測と主要領域で対比する。

Figure 1: (a) The density gradient and streamlines. (b) Shear stress at the ground surface. (c) The turbulent-to-laminar viscosity ratio $\mu_{r}$ and the local Knudsen number $\text{Kn}_{gll}$ , in the left and right half-domains, respectively. (d) Heat flux at the ground surface. (e) Contour of th

実験結果

リサーチクエスチョン

• RQ1 弱い稀人人化効果は局所的せん断層における重回に関する平均的荷重に影響を与えるか。
• RQ2 流れの転向や循環領域でNS構成関係が退化し、表面応力と熱流束に大きな誤差を生むか。
• RQ3 ジェット衝突で運動論からの高次非平衡応力はNS予測とどう異なるか。
• RQ4 運動論的シミュレーションは、連続体乱流モデルで捉えられない工学的に重要な表面量の顕著な差異を明らかにできるか。

主な発見

• NS-SSTに基づく予測は、 dominant regions において表面せん断応力を運動論解に対して約25–30%低く見積もる。
• NS-SSTは表面熱流束のピークを約50%過小評価し、ボルツマン解と比較して下流での減衰が速い。
• 非平衡応力（ボルツマン）は、NSせん断応力が符号反転によりほぼ0になる局所で局所的に支配的となり、構成定数退化を露呈する。
• 直接的ボルツ曼シミュレーションは、弱い稀人人化にもかかわらず、特定領域で高次運動論効果が応力平衡を支配することを示す。
• 退化メカニズムは、局所的な弱い稀人人化が不均一・高速度の flows における大きな誤差を生む理由を説明する。

Figure 2: Direct numerical simulation of the Boltzmann equation using the transient GSIS solver in turbulent-model-free mode [ Zeng2023GSIS ] . (a,b) The NS and Boltzmann shear stress, and their relative strength. (c) The Reynolds shear stress $R_{12}$ and the turbulence production term $\text{Prod}