[논문 리뷰] A theoretical one-dimensional model for variable-density Rayleigh-Taylor turbulence

Chian Yeh Goh, Guillaume Blanquart|arXiv (Cornell University)|2026. 02. 10.

Combustion and flame dynamics인용 수 0

한 줄 요약

논문은 가변 밀도 Rayleigh–Taylor 난류에 대한 Belen’kii–Fradkin 난류 확산율 모델을 재검토하고, 완전 및 간략화된 자체유사 ODE를 개발하며, 유사성을 물리 공간으로 매핑하고 DNS/실험과 비교하며, 질량 보존 보정이 합의를 개선하는 h ~ (ln R) g t^2 스케일링을 확인한다.

ABSTRACT

In an early theoretical work published in 1965, Belen'kii & Fradkin proposed a turbulent diffusivity model for Rayleigh--Taylor (RT) mixing. We review its derivation and present alternative arguments leading to the same final similarity equation. The original work then introduced an approximation that led to a simplified ordinary differential equation (ODE), which was used primarily to derive the important scaling result, $h \sim (\ln R)gt^2$. Here, we extend the analysis by examining the solutions to both the full similarity ODE and the simplified ODE in detail. It is shown that the full similarity equation captures many now well-known features of non-Boussinesq RT flows, including asymmetric spike and bubble growth and a systematic shift of velocity statistics toward the light-fluid side. Comparisons of the theoretical model with numerical and experimental studies show reasonable agreement in both spatial profiles and growth trends of mixing layer heights. We further show that a global mass correction applied to the simplified solution closely approximates the full solution, highlighting that, to leading order, RT mixing is governed by the competing dynamics between diffusion of $\ln \barρ$ and mass conservation.

연구 동기 및 목표

• 정적 배경을 가진 표준 RT에 대해 Belen’kii–Fradkin 난류 확산율 프레임워크를 명확화하고 간소화한다.
• 1D VD RT 난류에서 평균 밀도에 대해 완전한 및 단순화된 자체유사 ODE를 도출하고 분석한다.
• 이론과 RT 데이터를 연결하기 위해 유사성 변수와 물리 공간 양을 매핑한다.
• 비보우시네스 조건에서의 혼합층 높이 스케일링 및 스파이크/버블 성장 연구.
• DNS 및 실험 관찰과의 모델 검증 및 Atwood 수 전반에 걸친 적용 가능성 논의.

제안 방법

• 프란틸 혼합 길이 기반 난류 확산율을 채택하고 D_t^* ∼ v_t h_t 이며, v_t를 h_t를 통해 국부 에너지 고려로 연결하여 D_t^* = K ( (1/ρ̄) ∂ρ̄/∂y )^{1/2} g^{1/2} h_t^2 로 얻는다.
• 높은 Pe에서 3D RANS RT 형식을 ρ̄에 대한 1D 확산 문제로 축소하여 ∂ρ̄/∂t = ∂/∂y ( D_t ∂ρ̄/∂y ) 를 얻는다.
• 자체유사 유사방정식: ∂φ/∂τ = ∂/∂ξ [ φ^{-1/2} (∂φ/∂ξ)^{3/2} ], η = ξ/(a τ^{2/5})를 만족하는 φ(η)를 구한다.
• 물리적 공간으로 변환하여 무차원 높이 λ_i를 실제 높이 h_i와 연결: h_i(t) = K^2 λ_i λ_t^4 g t^2를 통해 h ~ g t^2 스케일링과 밀도비 비율 요인을 회복한다.
• Δρ/ρ_L이 작을 경우(Small R-1) h_T ∝ ln R에 대한 해를Analytical하게 제공하고, 완전한 ODE 해와 단순화된 ODE 해를 비교하며 질량 보정을 통해 전역 질량 보존을 강제하는 방법을 포함한다.

Figure 1: Solutions to full and simplified ODEs. (a,b) Full ODE: dimensionless diffusivity $x=(\varphi^{\prime}/\varphi)^{1/2}$ and density $\varphi$ ; (c,d) simplified ODE: dimensionless diffusivity $\hat{x}$ and density $\hat{\varphi}$ . Legend: $A=$ 0.01 (black solid), 0.2 (blue dashed), 0.5 (gre

실험 결과

연구 질문

• RQ1 완전한 자체유사 ODE와 단순화된 ODE가 Atwood 수 전체에서 혼합층 높이 및 밀도 프로파일에 대해 일관된 예측을 내는가?
• RQ2 1D 난류 확산율 모델이 스파이크/버블 비대칭, 가벼운 유체 쪽으로의 속도 프로파일 이동, 그리고 h의 ln R 스케일링과 같은 RT의 핵심 특징을 재현하는가?
• RQ3 물리 공간 매핑(높이, 밀도, 확산율)이 유사성 변수와 어떤 관계를 가지며 질량 보존이 간소화된 해에 어떤 영향을 주는가?
• RQ4 ln R 스케일링이 견고하게 관찰되는 Atwood 수(R) 범위는 어디이며 DNS/실험과의 예측 차이는 어떤가?
• RQ5 질량 보정된 간소화된 해가 전체 ODE 해 결과를 어느 정도 재현하는가?

주요 결과

• 완전한 ODE는 비보우시네스 RT 특징을 포착한다: 스파이크/버블 비대칭, 가벼운 유체 쪽으로의 속도 왜곡, 그리고 밀도 프로파일의 적절한 표현에서의 수렴.
• 높이는 h_i ∝ λ_i λ_t^4 g t^2로 스케일링되며, Atwood-number 효과를 비차원 요인과 연결한다.
• 작은-R에서의 해석적 결과는 h_T ∝ ln R로 나와 이전의 Belen’kii–Fradkin 발견 및 최근 RT 데이터와 일치한다.
• 단순화된 ODE에 질량 보정을 추가하면 전체 ODE 및 DNS와의 합치를 회복하며 질량 보존의 핵심 역할을 강조한다.
• 정규화된 확산율 x와 밀도 φ 프로파일은 Atwood 수 증가에 따라 질량이 가벼운 유체로 이동하는 경향을 드러내며 DNS 경향과 일치한다.
• 스파이크 높이는 Atwood 수가 증가할수록 버블 높이보다 g t^2 스케일링으로부터 더 많이 벗어나며 RT 관찰과 일치한다.
• 모델의 스파이크/버블 성장률 및 비율 α_s/α_b에 대한 예측은 실험 및 DNS에서 보고된 값의 범위와 일치한다.

Figure 2: Normalized solution to the full ODE: (a) normalized diffusivity $x/2\lambda_{T}^{2}$ , and (b) mole fraction $X$ . Legend: $A=$ 0.01 (black solid), 0.2 (blue dashed), 0.5 (green dash-dotted), and 0.8 (red dotted).