On non-Oberbeck–Boussinesq effects in three-dimensional Rayleigh–Bénard convection in glycerol

Rayleigh–Bénard convection in glycerol (Prandtl number Pr= 2547. 9 ) in a cylindrical cell with an aspect ratio of Γ = 1 was studied by means of three-dimensional direct numerical simulations (DNS). For that purpose, we implemented temperature-dependent material properties into our DNS code, by prescribing polynomial functions up to seventh order for the viscosity, the heat conductivity and the density. We performed simulations with the common Oberbeck–Boussinesq (OB) approximation and with non-Oberbeck–Boussinesq (NOB) effects within a range of Rayleigh numbers of 10⁵ ≤ Ra≤ 10⁹ . For the highest temperature differences, Δ = 80~K , the viscosity at the top is ~ 3600.167em % times higher than at the bottom, while the differences of the other material properties are less than 150.167em % . We analysed the temperature and velocity profiles and the thermal and viscous boundary-layer thicknesses. NOB effects generally lead to a breakdown of the top–bottom symmetry, typical for OB Rayleigh–Bénard convection. Under NOB conditions, the temperature in the centre of the cell Tc increases with increasing Δ and can be up to 15~K higher than under OB conditions. The comparison of our findings with several theoretical and empirical models showed that two-dimensional boundary-layer models overestimate the actual Tc , while models based on the temperature or velocity scales predict Tc very well with a standard deviation of 0. 4~K . Furthermore, the obtained temperature profiles bend closer towards the cold top plate and further away from the hot bottom plate. The situation for the velocity profiles is reversed: they bend farther away from the top plate and closer towards to the bottom plate. The top boundary layers are always thicker than the bottom ones. Their ratio is up to 2.5 for the thermal and up to 4.5 for the viscous boundary layers. In addition, the Reynolds number Re and the Nusselt number Nu were investigated: Re is higher and Nu is lower under NOB conditions. The Nusselt number Nu is influenced in a nonlinear way by NOB effects, stronger than was suggested by the two-dimensional simulations. The actual scaling of Nu with Ra in the NOB case is Nu∝ Ra0. 298 and is in excellent agreement with the experimental data.