We investigate the structures of the near-plate velocity and temperature profiles at different horizontal positions along the conducting bottom (and top) plate of a Rayleigh-B{e}nard convection cell, using two-dimensional (2D) numerical data obtained
at the Rayleigh number Ra=10^8 and the Prandtl number Pr=4.4 of an Oberbeck-Boussinesq flow with constant material parameters. The results show that most of the time, and for both velocity and temperature, the instantaneous profiles scaled by the dynamical frame method [Q. Zhou and K.-Q. Xia, Phys. Rev. Lett. 104, 104301 (2010) agree well with the classical Prandtl-Blasius laminar boundary layer (BL) profiles. Therefore, when averaging in the dynamical reference frames, which fluctuate with the respective instantaneous kinematic and thermal BL thicknesses, the obtained mean velocity and temperature profiles are also of Prandtl-Blasius type for nearly all horizontal positions. We further show that in certain situations the traditional definitions based on the time-averaged profiles can lead to unphysical BL thicknesses, while the dynamical method also in such cases can provide a well-defined BL thickness for both the kinematic and the thermal BLs.
Results on the Prandtl-Blasius type kinetic and thermal boundary layer thicknesses in turbulent Rayleigh-Benard convection in a broad range of Prandtl numbers are presented. By solving the laminar Prandtl-Blasius boundary layer equations, we calculat
e the ratio of the thermal and kinetic boundary layer thicknesses, which depends on the Prandtl number Pr only. It is approximated as $0.588Pr^{-1/2}$ for $Prll Pr^*$ and as $0.982 Pr^{-1/3}$ for $Pr^*llPr$, with $Pr^*= 0.046$. Comparison of the Prandtl--Blasius velocity boundary layer thickness with that evaluated in the direct numerical simulations by Stevens, Verzicco, and Lohse (J. Fluid Mech. 643, 495 (2010)) gives very good agreement. Based on the Prandtl--Blasius type considerations, we derive a lower-bound estimate for the minimum number of the computational mesh nodes, required to conduct accurate numerical simulations of moderately high (boundary layer dominated) turbulent Rayleigh-Benard convection, in the thermal and kinetic boundary layers close to bottom and top plates. It is shown that the number of required nodes within each boundary layer depends on Nu and Pr and grows with the Rayleigh number Ra not slower than $simRa^{0.15}$. This estimate agrees excellently with empirical results, which were based on the convergence of the Nusselt number in numerical simulations.
We analyze the reversals of the large scale flow in Rayleigh-Benard convection both through particle image velocimetry flow visualization and direct numerical simulations (DNS) of the underlying Boussinesq equations in a (quasi) two-dimensional, rect
angular geometry of aspect ratio 1. For medium Prandtl number there is a diagonal large scale convection roll and two smaller secondary rolls in the two remaining corners diagonally opposing each other. These corner flow rolls play a crucial role for the large scale wind reversal: They grow in kinetic energy and thus also in size thanks to plume detachments from the boundary layers up to the time that they take over the main, large scale diagonal flow, thus leading to reversal. Based on this mechanism we identify a typical time scale for the reversals. We map out the Rayleigh number vs Prandtl number phase space and find that the occurrence of reversals very sensitively depends on these parameters.
When the classical Rayleigh-Benard (RB) system is rotated about its vertical axis roughly three regimes can be identified. In regime I (weak rotation) the large scale circulation (LSC) is the dominant feature of the flow. In regime II (moderate rotat
ion) the LSC is replaced by vertically aligned vortices. Regime III (strong rotation) is characterized by suppression of the vertical velocity fluctuations. Using results from experiments and direct numerical simulations of RB convection for a cell with a diameter-to-height aspect ratio equal to one at $Ra sim 10^8-10^9$ ($Pr=4-6$) and $0 lesssim 1/Ro lesssim 25$ we identified the characteristics of the azimuthal temperature profiles at the sidewall in the different regimes. In regime I the azimuthal wall temperature profile shows a cosine shape and a vertical temperature gradient due to plumes that travel with the LSC close to the sidewall. In regime II and III this cosine profile disappears, but the vertical wall temperature gradient is still observed. It turns out that the vertical wall temperature gradient in regimes II and III has a different origin than that observed in regime I. It is caused by boundary layer dynamics characteristic for rotating flows, which drives a secondary flow that transports hot fluid up the sidewall in the lower part of the container and cold fluid downwards along the sidewall in the top part.
We studied the properties of the large-scale circulation (LSC) in turbulent Rayleigh-Benard (RB) convection by using results from direct numerical simulations in which we placed a large number of numerical probes close to the sidewall. The LSC orient
ation is determined by either a cosine or a polynomial fit to the azimuthal temperature or azimuthal vertical velocity profile measured with the probes. We study the LSC in Gamma=D/L=1/2 and Gamma=1 samples, where D is the diameter and L the height. For Pr=6.4 in an aspect ratio Gamma=1 sample at $Ra=1times10^8$ and $5times10^8$ the obtained LSC orientation is the same, irrespective of whether the data of only 8 or all 64 probes per horizontal plane are considered. In a Gamma=1/2 sample with $Pr=0.7$ at $Ra=1times10^8$ the influence of plumes on the azimuthal temperature and azimuthal vertical velocity profiles is stronger. Due to passing plumes and/or the corner flow the apparent LSC orientation obtained using a cosine fit can result in a misinterpretation of the character of the large-scale flow. We introduce the relative LSC strength, which we define as the ratio between the energy in the first Fourier mode and the energy in all modes that can be determined from the azimuthal temperature and azimuthal vertical velocity profiles, to further quantify the large-scale flow. For $Ra=1times10^8$ we find that this relative LSC strength is significantly lower in a Gamma=1/2 sample than in a Gamma=1 sample, reflecting that the LSC is much more pronounced in a Gamma=1 sample than in a Gamma=1/2 sample. The determination of the relative LSC strength can be applied directly to available experimental data to study high Rayleigh number thermal convection and rotating RB convection.
Results from direct numerical simulation for three-dimensional Rayleigh-Benard convection in samples of aspect ratio $Gamma=0.23$ and $Gamma=0.5$ up to Rayleigh number $Ra=2times10^{12}$ are presented. The broad range of Prandtl numbers $0.5<Pr<10$ i
s considered. In contrast to some experiments, we do not see any increase in $Nu/Ra^{1/3}$, neither due to $Pr$ number effects, nor due to a constant heat flux boundary condition at the bottom plate instead of constant temperature boundary conditions. Even at these very high $Ra$, both the thermal and kinetic boundary layer thicknesses obey Prandtl-Blasius scaling.
The effect of rotation on the boundary layers (BLs) in a Rayleigh-Benard (RB) system at a relatively low Rayleigh number, i.e. $Ra = 4times10^7$, is studied for different Pr by direct numerical simulations and the results are compared with laminar BL
theory. In this regime we find a smooth onset of the heat transfer enhancement as function of increasing rotation rate. We study this regime in detail and introduce a model based on the Grossmann-Lohse theory to describe the heat transfer enhancement as function of the rotation rate for this relatively low Ra number regime and weak background rotation $Rogtrsim 1$. The smooth onset of heat transfer enhancement observed here is in contrast to the sharp onset observed at larger $Ra gtrsim 10^8$ by Stevens {it{et al.}} [Phys. Rev. Lett. {bf{103}}, 024503, 2009], although only a small shift in the Ra-Ro-Pr phase space is involved.
