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

Abstract

The representation of buoyancy-driven turbulence in Reynolds-averaged Navier--Stokes (RANS) 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\'enard 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/3 Pr1/3 for Pr 1 and Nu Ra1/3 Pr-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\'enard 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.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…