On non-Oberbeck-Boussinesq effects in three-dimensional Rayleigh-Benard convection in glycerol

Rayleigh-Benard convection in glycerol (Prandtl number Pr = 2547.9) in a cylindrical cell with an aspect ratio of Gamma = 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(5) <= Ra <= 10(9). For the highest temperature differences, Delta = 80 K, the viscosity at the top is similar to 360% times higher than at the bottom, while the differences of the other material properties are less than 15%. 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-Benard convection. Under NOB conditions, the temperature in the centre of the cell T-c increases with increasing Delta 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 T-c, while models based on the temperature or velocity scales predict T-c 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 R e and the Nusselt number Nu were investigated: R e 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 R a in the NOB case is Nu proportional to Ra-0.298 and is in excellent agreement with the experimental data.