Analysis and reformulation of the k –ω turbulence model for buoyancy-driven thermal convection

The representation of buoyancy-driven turbulence in Reynolds-averaged Navier–Stokes models remains unresolved, with no widely accepted standard formulation. A key difficulty is the lack of analytical guidance for incorporating buoyant effects, particularly under unstable stratification. This study derives an analytical solution of the standard k–ω model for Rayleigh–Bénard convection in an infinite layer, where turbulent kinetic energy is generated solely by buoyancy. The solution provides explicit scaling relations among the Rayleigh (Ra), Prandtl (Pr), and Nusselt (Nu) numbers that capture the simulation trends: Nu∼Ra1/3Pr1/3 for Pr≪1 and Nu∼Ra1/3Pr−0.415 for Pr≫1. This framework quantifies the discrepancies in the conventional buoyancy treatment and clarifies their origin. Informed by this analysis, the buoyancy-related modeling terms are reformulated to recover the measured trends: namely Nu∼Pr1/8 for Pr≪1 and Nu∼Pr0 for Pr≫1 at moderate Ra. Only two dimensionless algebraic functions are introduced, which vanish in the absence of buoyancy, ensuring full compatibility with the standard closure. The corrected model is validated across a range of buoyancy-driven flows, including two-dimensional Rayleigh–Bénard convection, internally heated convection in two configurations, unstably stratified Couette flow, and vertically heated natural convection with varying aspect ratios. Across all cases, the corrected model provides significantly improved predictions of mean temperature fields and turbulent heat flux distributions.