No Arabic abstract
As shown in earlier work (Ahlers et al., J. Fluid Mech. 569, p.409 (2006)), non-Oberbeck Boussinesq (NOB) corrections to the center temperature in turbulent Rayleigh-Benard convection in water and also in glycerol are governed by the temperature dependences of the kinematic viscosity and the thermal diffusion coefficient. If the working fluid is ethane close to the critical point the origin of non-Oberbeck-Boussinesq corrections is very different, as will be shown in the present paper. Namely, the main origin of NOB corrections then lies in the strong temperature dependence of the isobaric thermal expansion coefficient beta(T). More precisely, it is the nonlinear T-dependence of the density rho(T) in the buoyancy force which causes another type of NOB effect. We demonstrate that through a combination of experimental, numerical, and theoretical work, the latter in the framework of the extended Prandtl-Blasius boundary layer theory developed in Ahlers et al., J. Fluid Mech. 569, p.409 (2006). The latter comes to its limits, if the temperature dependence of the thermal expension coefficient beta(T) is significant.
We numerically analyze Non-Oberbeck-Boussinesq (NOB) effects in two-dimensional Rayleigh-Benard flow in glycerol, which shows a dramatic change in the viscosity with temperature. The results are presented both as functions of the Rayleigh number (Ra) up to $10^8$ (for fixed temperature difference between the top and bottom plates) and as functions of non-Oberbeck-Boussinesqness or NOBness ($Delta$) up to 50 K (for fixed Ra). For this large NOBness the center temperature $T_c$ is more than 5 K larger than the arithmetic mean temperature $T_m$ between top and bottom plate and only weakly depends on Ra. To physically account for the NOB deviations of the Nusselt numbers from its Oberbeck-Boussinesq values, we apply the decomposition of $Nu_{NOB}/Nu_{OB}$ into the product of two effects, namely first the change in the sum of the top and bottom thermal BL thicknesses, and second the shift of the center temperature $T_c$ as compared to $T_m$. While for water the origin of the $Nu$ deviation is totally dominated by the second effect (cf. Ahlers et al., J. Fluid Mech. 569, pp. 409 (2006)) for glycerol the first effect is dominating, in spite of the large increase of $T_c$ as compared to $T_m$.
Non-Oberbeck-Boussinesq (NOB) effects on the flow organization in two-dimensional Rayleigh-Benard turbulence are numerically analyzed. The working fluid is water. We focus on the temperature profiles, the center temperature, the Nusselt number, and on the analysis of the velocity field. Several velocity amplitudes (or Reynolds numbers) and several kinetic profiles are introduced and studied; these together describe the various features of the rather complex flow organization. The results are presented both as functions of the Rayleigh number Ra (with Ra up to 10^8) for fixed temperature difference (Delta) between top and bottom plates and as functions of Delta (non-Oberbeck-Boussinesqness) for fixed Ra with Delta up to 60 K. All results are consistent with the available experimental NOB data for the center temperature Tc and the Nusselt number ratio Nu_{NOB}/Nu_{OB} (the label OB meaning that the Oberbeck-Boussinesq conditions are valid). Beyond Ra ~ 10^6 the flow consists of a large diagonal center convection roll and two smaller rolls in the upper and lower corners. In the NOB case the center convection roll is still characterized by only one velocity scale.
We report an experimental study of the three-dimensional spatial structure of the low frequency temperature oscillations in a cylindrical Rayleigh-B{e}nard convection cell. It is found that thermal plumes are not emitted periodically, but randomly and continuously, from the top and bottom plates. We further found that the oscillation of the temperature field does not originate from the boundary layers, but rather is a result of the horizontal motion of the hot ascending and cold descending fluids being modulated by the twisting and sloshing motion of the bulk flow field.
We report an experimental study aiming to clarify the role of boundary conditions (BC) in high Rayleigh number $10^8 < {rm{Ra}} < 3 times 10^{12}$ turbulent thermal convection of cryogenic helium gas. We switch between BC closer to constant heat flux (CF) and constant temperature (CT) applied to the highly conducting bottom plate of the aspect ratio one cylindrical cell 30 cm in size, leading to dramatic changes in the temperature probability density function and in power spectral density of the temperature fluctuations measured at the bottom plate, while the dynamic thermal behaviour of the top plate and bulk convective flow remain unaffected. Within our experimental accuracy, we find no appreciable changes in Reynolds number Re(Ra) scaling, in the dimensionless heat transfer efficiency expressed via Nusselt number Nu(Ra) scaling, nor in the rate of direction reversals of large scale circulation.
Recently, in Zhang et al. (2020), it was found that in rapidly rotating turbulent Rayleigh-Benard convection (RBC) in slender cylindrical containers (with diameter-to-height aspect ratio $Gamma=1/2$) filled with a small-Prandtl-number fluid ($Pr approx0.8$), the Large Scale Circulation (LSC) is suppressed and a Boundary Zonal Flow (BZF) develops near the sidewall, characterized by a bimodal PDF of the temperature, cyclonic fluid motion, and anticyclonic drift of the flow pattern (with respect to the rotating frame). This BZF carries a disproportionate amount ($>60%$) of the total heat transport for $Pr < 1$ but decreases rather abruptly for larger $Pr$ to about $35%$. In this work, we show that the BZF is robust and appears in rapidly rotating turbulent RBC in containers of different $Gamma$ and in a broad range of $Pr$ and $Ra$. Direct numerical simulations for $0.1 leq Pr leq 12.3$, $10^7 leq Ra leq 5times10^{9}$, $10^{5} leq 1/Ek leq 10^{7}$ and $Gamma$ = 1/3, 1/2, 3/4, 1 and 2 show that the BZF width $delta_0$ scales with the Rayleigh number $Ra$ and Ekman number $Ek$ as $delta_0/H sim Gamma^{0} Pr^{{-1/4, 0}} Ra^{1/4} Ek^{2/3}$ (${Pr<1, Pr>1}$) and the drift frequency as $omega/Omega sim Gamma^{0} Pr^{-4/3} Ra Ek^{5/3}$, where $H$ is the cell height and $Omega$ the angular rotation rate. The mode number of the BZF is 1 for $Gamma lesssim 1$ and $2 Gamma$ for $Gamma$ = {1,2} independent of $Ra$ and $Pr$. The BZF is quite reminiscent of wall mode states in rotating convection.