[논문 리뷰] Small Rarefaction, Large Consequences: Limits of Navier Stokes Turbulence Simulations

Songyan Tian, Lei Wu|arXiv (Cornell University)|2026. 02. 09.

Gas Dynamics and Kinetic Theory인용 수 0

한 줄 요약

논문은 약한 희박효과가 난류 제트 충돌에서 국소적으로 지배적이 될 수 있음을 보여주며, 이는 직접 볼츠만 시뮬레이션에서 surface shear stress와 heat flux에 대한 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 난류 모델링의 한계를 평가한다.
• 비평형 응력이 표면 하중을 지배하는 구성 재현성의 저하를 демон스트레이션한다.
• Jet-impingement 설정에서 연속체(NS) 예측과 운동학(볼츠만) 해를 비교한다.
• 고속 공기역학 및 행성 탐사 응용에 대한 통찰을 제공한다.

제안 방법

• modified Rykov collision model을 사용하는 이산 속도 법으로 질소 배기가스에 대해 Boltzmann 방정식을 푼다.
• Explicit domain decomposition 없이 다중 스케일 연속-희박 다이나믹스를 포착하기 위해 coarse-grained GSIS–SST 프레임워크를 결합한다.
• SST 폐쇄를 넘어 난류와 희박화 간 상호 작용을 검증하기 위해 시간 해상 Boltzmann 시뮬레이션을 수행한다.
• 총 응력을 NS( laminar+SST )와 Boltzmann 비평형 구성요소로 분해하여 constitutive 관계를 분석한다.
• NS 응력이 부호 반전으로 인해 0에 수렴하는 영역을 식별하기 위해 국소 Knudsen 수와 기울기 길이 Kn를 평가한다.
• 주요 영역에서 NS-SST 예측과 대조하기 위해 직접적인 Boltzmann 결과를 사용한다.

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를 기반으로 한 NS 예측은 지배 영역에서 운동학적 해보다 표면 전단응력을 약 25–30% 정도 과소 예측한다.
• NS-SST는 표면 열유속의 피크를 대략 50% 정도 과소 예측하고 볼츠만 해에 비해 하류로 갈수록 더 빠르게 감소한다.
• 비평형(볼츠만) 응력은 NS 전단응력이 부호 반전에 의해 0에 접근하는 국부적으로 지배적이 되며 구성 재현성의 저하를 드러낸다.
• 직접 볼츠만 시뮬레이션은 희박효과가 약해도 응력 균형을 지배하는 고차 운동학적 효과가 국소적으로 존재한다는 것을 보여준다.
• 퇴화 메커니즘은 약한 희박효과가 불균일하고 고속인 흐름에서 대류 난류 예측에 큰 오차를 초래할 수 있는 이유를 설명한다.

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}