Direct numerical simulations of Rayleigh–Bénard convection in water with non-Oberbeck–Boussinesq effects

Abstract

Rayleigh–Bénard convection in water is studied by means of direct numerical simulations, taking into account the variation of properties. The simulations considered a three-dimensional (3-D) cavity with a square cross-section and its two-dimensional (2-D) equivalent, covering a Rayleigh number range of 106⩽Ra⩽10910^{6}\leqslant Ra\leqslant 10^{9} and using temperature differences up to 60 K. The main objectives of this study are (i) to investigate and report differences obtained by 2-D and 3-D simulations and (ii) to provide a first appreciation of the non-Oberbeck–Boussinesq (NOB) effects on the near-wall time-averaged and root-mean-squared (r.m.s.) temperature fields. The Nusselt number and the thermal boundary layer thickness exhibit the most pronounced differences when calculated in two dimensions and three dimensions, even though the RaRa scaling exponents are similar. These differences are closely related to the modification of the large-scale circulation pattern and become less pronounced when the NOB values are normalised with the respective Oberbeck–Boussinesq (OB) values. It is also demonstrated that NOB effects modify the near-wall temperature statistics, promoting the breaking of the top–bottom symmetry which characterises the OB approximation. The most prominent NOB effect in the near-wall region is the modification of the maximum r.m.s. values of temperature, which are found to increase at the top and decrease at the bottom of the cavity.