Turbulent Rayleigh-Bénard convection displays a large-scale order in the form of rolls and cells on lengths larger than the layer height once the fluctuations of temperature and velocity are removed. These turbulent superstructures are reminiscent of the patterns close to the onset of convection. Here we report numerical simulations of turbulent convection in fluids at different Prandtl number ranging from 0.005 to 70 and for Rayleigh numbers up to 107. We identify characteristic scales and times that separate the fast, small-scale turbulent fluctuations from the gradually changing large-scale superstructures. The characteristic scales of the large-scale patterns, which change with Prandtl and Rayleigh number, are also correlated with the boundary layer dynamics, and in particular the clustering of thermal plumes at the top and bottom plates. Our analysis suggests a scale separation and thus the existence of a simplified description of the turbulent superstructures in geo- and astrophysical settings.

We discuss two aspects of turbulent Rayleigh-Bénard convection (RBC) on the basis of highresolution direct numerical simulations in a unique setting; a closed cylindrical cell of aspect ratio of one. First, we present a comprehensive comparison of statistical quantities such as energy dissipation rates and boundary layer thickness scales. Data are used from three simulation run series at Prandtl numbers P r that cover two orders of magnitude. In contrast to most previous studies in RBC the focus of the present work is on convective turbulence at very low Prandtl numbers including P r = 0.021 for liquid mercury or gallium and P r = 0.005 for liquid sodium. In this parameter range of RBC, inertial effects cause a dominating turbulent momentum transport that is in line with highly intermittent fluid turbulence both in the bulk and in the boundary layers and thus should be able to trigger a transition to the fully turbulent boundary layers of the ultimate regime of convection for higher Rayleigh number. Secondly, we predict the ranges of Rayleigh numbers for which the viscous boundary layer will transition to turbulence and the flow as a whole will cross over into the ultimate regime. These transition ranges are obtained by extrapolation from our simulation data. The extrapolation methods are based on the large-scale properties of the velocity profile. Two of the three methods predict similar ranges for the transition to ultimate convection when their uncertainties are taken into account. All three extrapolation methods indicate that the range of critical Rayleigh numbers Rac is shifted to smaller magnitudes as the Prandtl number becomes smaller.

The global transport of heat and momentum in turbulent convection is constrained by thin thermal and viscous boundary layers at the heated and cooled boundaries of the system. This bottleneck is thought to be lifted once the boundary layers themselves become fully turbulent at very high values of the Rayleigh numberRa—the dimensionless parameter that describes the vigor of convective turbulence. Laboratory experiments in cylindrical cells forRa≳1012have reported different outcomes on the putative heat transport law. Here we show, by direct numerical simulations of three-dimensional turbulent Rayleigh–Bénard convection flows in a slender cylindrical cell of aspect ratio1/10, that the Nusselt number—the dimensionless measure of heat transport—follows the classical power law ofNu=(0.0525±0.006)×Ra0.331±0.002up toRa=1015. Intermittent fluctuations in the wall stress, a blueprint of turbulence in the vicinity of the boundaries, manifest at allRaconsidered here, increasing with increasingRa, and suggest that an abrupt transition of the boundary layer to turbulence does not take place.

We present high-resolution direct numerical simulation studies of turbulent Rayleigh-Bénard convection in a closed cylindrical cell with an aspect ratio of one. The focus of our analysis is on the finest scales of convective turbulence, in particular the statistics of the kinetic energy and thermal dissipation rates in the bulk and the whole cell. The fluctuations of the energy dissipation field can directly be translated into a fluctuating local dissipation scale which is found to develop ever finer fluctuations with increasing Rayleigh number. The range of these scales as well as the probability of high-amplitude dissipation events decreases with increasing Prandtl number. In addition, we examine the joint statistics of the two dissipation fields and the consequences of high-amplitude events. We have also investigated the convergence properties of our spectral element method and have found that both dissipation fields are very sensitive to insufficient resolution. We demonstrate that global transport properties, such as the Nusselt number, and the energy balances are partly insensitive to insufficient resolution and yield correct results even when the dissipation fields are under-resolved. Our present numerical framework is also compared with high-resolution simulations which use a finite difference method. For most of the compared quantities the agreement is found to be satisfactory.

Statistical properties of turbulent Rayleigh-Bénard convection at low Prandtl numbers P r, which are typical for liquid metals such as mercury or gallium (P r 0.021) or liquid sodium (P r 0.005), are investigated in high-resolution three-dimensional spectral element simulations in a closed cylindrical cell with an aspect ratio of one and are compared to previous turbulent convection simulations in air for P r = 0.7. We compare the scaling of global momentum and heat transfer. The scaling exponent β of the power law N u = αRa β is β = 0.265 ± 0.01 for P r = 0.005 and β = 0.26 ± 0.01 for P r = 0.021, which are smaller than that for convection in air (P r = 0.7, β = 0.29 ± 0.01). These exponents are in agreement with experiments. Mean profiles of the root-mean-square velocity as well as the thermal and kinetic energy dissipation rates have growing amplitudes with decreasing Prandtl number which underlies a more vigorous bulk turbulence in the low-P r regime. The skin-friction coefficient displays a Reynolds-number dependence that is close to that of an isothermal, intermittently turbulent velocity boundary layer. The thermal boundary layer thicknesses are larger as P r decreases and conversely the velocity boundary layer thicknesses become smaller. We investigate the scaling exponents and find a slight decrease in exponent magnitude for the thermal boundary layer thickness as P r decreases, but find the opposite case for the velocity boundary layer thickness scaling. A growing area fraction of turbulent patches close to the heating and cooling plates can be detected by exceeding a locally defined shear Reynolds number threshold. This area fraction is larger for lower P r at the same Ra, but the scaling exponent of its growth with Rayleigh number is reduced. Our analysis of the kurtosis of the locally defined shear Reynolds number demonstrates that the intermittency in the boundary layer is significantly increased for the lower Prandtl number and for sufficiently high Rayleigh number compared to convection in air. This complements our previous findings of enhanced bulk intermittency in low-Prandtl-number convection.