Search

We study the dynamics of thermal and momentum boundary regions in three-dimensional direct numerical simulations of Rayleigh-B\'enard convection for the Rayleigh number range $10^5 \le Ra \le 10^{11}$ and $Pr=0.7$. Using a Cartesian slab with horizontal periodic boundary conditions and an aspect ratio of 4, we obtain statistical homogeneity in the horizontal $x$- and $y$-directions, thus approximating best an extended convection layer relevant for most geo- and astrophysical flow applications. We observe upon canonical use of combined long-time and area averages, with averaging periods of at least 100 free-fall times, that a global coherent mean flow is practically absent and that the magnitude of the velocity fluctuations is larger than the mean by up to 2 orders of magnitude. The velocity field close to the wall is a collection of differently oriented local shear-dominated flow patches interspersed by extensive shear-free incoherent regions which can be as large as the whole cross section, unlike for a closed cylindrical convection cell of aspect ratio of the order 1. The incoherent regions occupy a $60\%$ area fraction for all Rayleigh numbers investigated here. Rather than resulting in a pronounced mean with small fluctuations about such a mean, as found in small-aspect-ratio convection, the velocity field is dominated by strong fluctuations of all three components around a non-existent or weak mean. We discuss the consequences of these observations for convection layers with larger aspect ratios, including boundary layer instabilities and the resulting turbulent heat transport., Comment: 18 pages, 10 figures

Inverse cascades of kinetic energy and thermal variance in the subset of vertically homogeneous modes in spectral space are found to cause a slow aggregation to a pair of convective supergranules that eventually fill the whole horizontally extended, three-dimensional, turbulent Rayleigh-B\'{e}nard convection layer when a heat flux is prescribed at the top and bottom. An additional weak rotation of the layer around the vertical axis stops this aggregation at a scale that is smaller than the lateral domain extension and ceases the inverse cascade for the thermal variance. The inverse cascade for the kinetic energy remains intact, even for times at which the root mean square values of temperature and velocity have reached the statistically steady state. This kinetic energy inverse cascade sustains the horizontally extended convection patterns which are best visible in the temperature field. The resulting characteristic length of the aggregated convection patterns depends on the thermal driving and linearly on the strength of rotation. Our study demonstrates the importance of inverse energy cascades beyond the two-dimensional turbulence case in a three-dimensional convection flow that is subject to a multi-scale energy injection by thermal plumes and driven by boundary heat fluxes as typically present in natural geo- and astrophysical systems, such as solar convection., Comment: 15 pages, 10 figures

Turbulent convection processes in nature are often found to be organized in a hierarchy of plume structures and flow patterns. The gradual aggregation of convection cells or granules to a supergranule which eventually fills the whole horizontal layer is reported and analysed in spectral element direct numerical simulations of three-dimensional turbulent Rayleigh-B\'{e}nard convection at an aspect ratio of $60$. The formation proceeds over a time span of more than $10^4$ convective time units for the largest accessible Rayleigh number and occurs only when the turbulence is driven by a constant heat flux which is imposed at the bottom and top planes enclosing the convection layer. The resulting gradual inverse cascade process is observed for both temperature variance and turbulent kinetic energy. An additional analysis of the leading Lyapunov vector field for the full turbulent flow trajectory in its high-dimensional phase space demonstrates that turbulent flow modes at a certain scale continue to give rise locally to modes with longer wavelength in the turbulent case. As a consequence successively larger convection patterns grow until the horizontal extension of the layer is reached. This instability mechanism, which is known to exist near the onset of constant heat flux-driven convection, is shown here to persist into the fully developed turbulent flow regime thus connecting weakly nonlinear pattern formation with the one in fully developed turbulence. We discuss possible implications of our study for observed, but not yet consistently numerically reproducible, solar supergranulation which could lead to improved simulation models of surface convection in the Sun., Comment: 15 pages, 11 figures

Turbulent Rayleigh-B\'enard 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. They are analyzed by numerical simulations of turbulent convection in fluids at different Prandtl number ranging from 0.005 to 70 and for Rayleigh numbers up to $10^7$. For each case, 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 found to be 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., Comment: 16 pages (incl. Supplementary Material), 12 figures (all with downsized figure size)

We discuss two aspects of turbulent Rayleigh-B\'{e}nard convection (RBC) on the basis of high-resolution 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 $Pr$ 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 $Pr=0.021$ for liquid mercury or gallium and $Pr=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 $Ra_c$ is shifted to smaller magnitudes as the Prandtl number becomes smaller.

The boundary layer structure of the velocity and temperature fields in turbulent Rayleigh-Benard flows in closed cylindrical cells of unit aspect ratio is revisited from a transitional and turbulent viscous boundary layer perspective. When the Rayleigh number is large enough, the dynamics at the bottom and top plates can be separated into an impact region of downwelling plumes, an ejection region of upwelling plumes and an interior region away from the side walls. The latter is dominated by the shear of the large-scale circulation (LSC) roll which fills the whole cell and continuously varies its orientation. The working fluid is liquid mercury or gallium at a Prandtl number Pr=0.021 for Rayleigh numbers between Ra=3e+5 and 4e+8. The generated turbulent momentum transfer corresponds to macroscopic flow Reynolds numbers with values between 1800 and 46000. It is shown that the viscous boundary layers for the largest Rayleigh numbers are highly transitional and obey properties that are directly comparable to transitional channel flows at friction Reynolds numbers Re_tau slightly below 100. The transitional character of the viscous boundary layer is also underlined by the strong enhancement of the fluctuations of the wall stress components with increasing Rayleigh number. An extrapolation of our analysis data suggests that the friction Reynolds number Re_tau in the velocity boundary layer can reach values of 200 for Ra beyond 1e+11. Thus the viscous boundary layer in a liquid metal flow would become turbulent at a much lower Rayleigh number than for turbulent convection in gases and gas mixtures., Comment: 16 pages, 15 figures in reduced quality, accepted for Phys. Rev. Fluids

An extreme dissipation event in the bulk of a closed three-dimensional turbulent convection cell is found to be correlated with a strong reduction of the large-scale circulation flow in the system that happens at the same time as a plume emission event from the bottom plate. The reduction in the large-scale circulation opens the possibility for a nearly frontal collision of down- and upwelling plumes and the generation of a high-amplitude thermal dissipation layer in the bulk. This collision is locally connected to a subsequent high-amplitude energy dissipation event in the form of a strong shear layer. Our analysis illustrates the impact of transitions in the large-scale structures on extreme events at the smallest scales of the turbulence, a direct link that is observed in a flow with boundary layers. We also show that detection of extreme dissipation events which determine the far-tail statistics of the dissipation fields in the bulk requires long-time integrations of the equations of motion over at least hundred convective time units., Comment: accepted for Phys. Rev. E, all figures are downsized

Statistical properties of turbulent Rayleigh-Benard convection at low Prandtl numbers (Pr), which are typical for liquid metals such as mercury, gallium or liquid sodium, 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. We compare the scaling of global momentum and heat transfer. The scaling exponents are found to be 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-Pr 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 Pr 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 Pr 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 Pr at the same Ra. 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.

We compute fully local boundary layer scales in three-dimensional turbulent Rayleigh-Benard convection. These scales are directly connected to the highly intermittent fluctuations of the fluxes of momentum and heat at the isothermal top and bottom walls and are statistically distributed around the corresponding mean thickness scales. The local boundary layer scales also reflect the strong spatial inhomogeneities of both boundary layers due to the large-scale, but complex and intermittent, circulation that builds up in closed convection cells. Similar to turbulent boundary layers, we define inner scales based on local shear stress which can be consistently extended to the classical viscous scales in bulk turbulence, e.g. the Kolmogorov scale, and outer scales based on slopes at the wall. We discuss the consequences of our generalization, in particular the scaling of our inner and outer boundary layer thicknesses and the resulting shear Reynolds number with respect to Rayleigh number. The mean outer thickness scale for the temperature field is close to the standard definition of a thermal boundary layer thickness. In the case of the velocity field, under certain conditions the outer scale follows a similar scaling as the Prandtl-Blasius type definition with respect to Rayleigh number, but differs quantitatively. The friction coefficient c_epsilon scaling is found to fall right between the laminar and turbulent limits which indicates that the boundary layer exhibits transitional behavior. Additionally, we conduct an analysis of the recently suggested dissipation layer thickness scales versus Rayleigh number and find a transition in the scaling. We also performed one study of aspect ratio equal to three in the case of Ra=1e+8., Comment: 29 pages, 16 figures, some figures in reduced quality, abstract is shortened

We present high-resolution direct numerical simulation studies of turbulent Rayleigh-Benard 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 also have 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., Comment: 33 pages, 24 figures