Direct Numerical Simulation of the Differentially Heated Cavity and Comparison with the κ-ε Model for High Rayleigh Numbers

This study presents a numerical comparison between Direct numerical simulation (DNS) and the standard κ-ε turbulence model to evaluate natural convection in a two-dimensional, differentially heated, air-filled cavity over the Rayleigh number range 10³ to 10¹⁰. The objective is to assess the predictive capabilities of both methods across laminar and turbulent regimes, with a particular emphasis on the quantitative comparison of thermal characteristics under high Rayleigh number conditions. The Navier–Stokes and energy equations were solved using the finite element method with Boussinesq approximation, employing refined meshes near the hot and cold walls to resolve thermal and velocity boundary layers. The results indicate that for Ra ≤ 10⁶, the κ-ε model significantly underestimates temperature gradients, maximum velocities, and average Nusselt numbers, with errors up to 19.39%, due to isotropic assumptions and empirical formulation. DNS, in contrast, achieves global energy balance errors of less than 0.0018% across the entire range. As Ra increases, the κ-ε model predictions converge to DNS, with Nusselt number deviations dropping below 1.2% at Ra = 10¹⁰. Streamlines, temperature profiles, and velocity distributions confirm that DNS captures flow dynamics more accurately, particularly near the wall vortices. These findings validate DNS as a reference solution for high-Ra natural convection and establish benchmark data for assessing turbulence models in confined geometries