[論文レビュー] A theoretical one-dimensional model for variable-density Rayleigh-Taylor turbulence

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

Combustion and flame dynamics被引用数 0

ひとこと要約

要約: 本論文は、可変密度のRayleigh–Taylor乱流に対するBelen’kii–Fradkinの乱流拡散係数モデルを再検討し、完全形と簡略化された自己相似ODEを導出・解析し、相似性を物理空間へ写像し、DNS/実験と比較して、質量保存補正が適合性を改善することを確認する。

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と簡略化ODEを導出・解析する。
• 相似変数を物理空間量へ写像し、RTデータとの結びつきを確立する。
• 非ボース-クネス条件下での混合層高度のスケーリングとスパイク/バブルの成長を調査する。
• DNSおよび実験観測とモデルを検証し、Atwood数全体への適用性を議論する。

提案手法

• Prandtl混合長に基づく乱流拡散係数を採用し D_t^* ∼ v_t h_t として v_t を局所エネルギーの考慮から h_t に関連付け、D_t^* = K ( (1/ρ̄) ∂ρ̄/∂y )^{1/2} g^{1/2} h_t^2 を得る。
• 3D RANS RT方程式を高Pe領域においてρ̄の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 が小さい場合（R-1が小さい）には h_T ∝ ln R となる解析解を提供し、完全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乱流拡散モデルはスパイク/バブルの非対称性、軽容液に向かう速度プロファイルのシフト、 heights の ln R スケーリングといったRTの主要特徴を再現できるか。
• RQ3 物理空間への写像（高さ、密度、拡散係数）と相似変数の関係はどうなるか、質量保存が簡略化解に与える影響は。
• RQ4 ln R スケーリングが頑健に観測されるAtwood数の範囲はどれか、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 数の増加とともに質量が light fluidへ移動することを示し、DNSの傾向と一致。
• スパイク高さは Atwood 数が大きいほど A 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).