scholarly journals Azimuthal diffusion of the large-scale-circulation plane, and absence of significant non-Boussinesq effects, in turbulent convection near the ultimate-state transition

2016 ◽  
Vol 791 ◽  
Author(s):  
Xiaozhou He ◽  
Eberhard Bodenschatz ◽  
Guenter Ahlers

We present measurements of the orientation ${\it\theta}_{0}$ and temperature amplitude ${\it\delta}$ of the large-scale circulation in a cylindrical sample of turbulent Rayleigh–Bénard convection (RBC) with aspect ratio ${\it\Gamma}\equiv D/L=1.00$ ($D$ and $L$ are the diameter and height respectively) and for the Prandtl number $Pr\simeq 0.8$. The results for ${\it\theta}_{0}$ revealed a preferred orientation with up-flow in the west, consistent with a broken azimuthal invariance due to the Earth’s Coriolis force (see Brown & Ahlers (Phys. Fluids, vol. 18, 2006, 125108)). They yielded the azimuthal diffusivity $D_{{\it\theta}}$ and a corresponding Reynolds number $Re_{{\it\theta}}$ for Rayleigh numbers over the range $2\times 10^{12}\lesssim Ra\lesssim 1.5\times 10^{14}$. In the classical state ($Ra\lesssim 2\times 10^{13}$) the results were consistent with the measurements by Brown & Ahlers (J. Fluid Mech., vol. 568, 2006, pp. 351–386) for $Ra\lesssim 10^{11}$ and $Pr=4.38$, which gave $Re_{{\it\theta}}\propto Ra^{0.28}$, and with the Prandtl-number dependence $Re_{{\it\theta}}\propto Pr^{-1.2}$ as found previously also for the velocity-fluctuation Reynolds number $Re_{V}$ (He et al., New J. Phys., vol. 17, 2015, 063028). At larger $Ra$ the data for $Re_{{\it\theta}}(Ra)$ revealed a transition to a new state, known as the ‘ultimate’ state, which was first seen in the Nusselt number $Nu(Ra)$ and in $Re_{V}(Ra)$ at $Ra_{1}^{\ast }\simeq 2\times 10^{13}$ and $Ra_{2}^{\ast }\simeq 8\times 10^{13}$. In the ultimate state we found $Re_{{\it\theta}}\propto Ra^{0.40\pm 0.03}$. Recently, Skrbek & Urban (J. Fluid Mech., vol. 785, 2015, pp. 270–282) claimed that non-Oberbeck–Boussinesq effects on the Nusselt and Reynolds numbers of turbulent RBC may have been interpreted erroneously as a transition to a new state. We demonstrate that their reasoning is incorrect and that the transition observed in the Göttingen experiments and discussed in the present paper is indeed to a new state of RBC referred to as ‘ultimate’.

2016 ◽  
Vol 790 ◽  
Author(s):  
Olga Shishkina ◽  
Susanne Horn

By means of direct numerical simulations (DNS) we investigate the effect of a tilt angle ${\it\beta}$, $0\leqslant {\it\beta}\leqslant {\rm\pi}/2$, of a Rayleigh–Bénard convection (RBC) cell of aspect ratio 1, on the Nusselt number $\mathit{Nu}$ and Reynolds number $\mathit{Re}$. The considered Rayleigh numbers $\mathit{Ra}$ range from $10^{6}$ to $10^{8}$, the Prandtl numbers range from 0.1 to 100 and the total number of the studied cases is 108. We show that the $\mathit{Nu}\,({\it\beta})/\mathit{Nu}(0)$ dependence is not universal and is strongly influenced by a combination of $\mathit{Ra}$ and $\mathit{Pr}$. Thus, with a small inclination ${\it\beta}$ of the RBC cell, the Nusselt number can decrease or increase, compared to that in the RBC case, for large and small $\mathit{Pr}$, respectively. A slight cell tilt may not only stabilize the plane of the large-scale circulation (LSC) but can also enforce an LSC for cases when the preferred state in the perfect RBC case is not an LSC but a more complicated multiple-roll state. Close to ${\it\beta}={\rm\pi}/2$, $\mathit{Nu}$ and $\mathit{Re}$ decrease with increasing ${\it\beta}$ in all considered cases. Generally, the $\mathit{Nu}({\it\beta})/\mathit{Nu}(0)$ dependence is a complicated, non-monotonic function of ${\it\beta}$.


2012 ◽  
Vol 706 ◽  
pp. 5-33 ◽  
Author(s):  
Nan Shi ◽  
Mohammad S. Emran ◽  
Jörg Schumacher

AbstractThe structure of the boundary layers in turbulent Rayleigh–Bénard convection is studied by means of three-dimensional direct numerical simulations. We consider convection in a cylindrical cell at aspect ratio one for Rayleigh numbers of $\mathit{Ra}= 3\ensuremath{\times} 1{0}^{9} $ and $3\ensuremath{\times} 1{0}^{10} $ at fixed Prandtl number $\mathit{Pr}= 0. 7$. Similar to the experimental results in the same setup and for the same Prandtl number, the structure of the laminar boundary layers of the velocity and temperature fields is found to deviate from the prediction of Prandtl–Blasius–Pohlhausen theory. Deviations decrease when a dynamical rescaling of the data with an instantaneously defined boundary layer thickness is performed and the analysis plane is aligned with the instantaneous direction of the large-scale circulation in the closed cell. Our numerical results demonstrate that important assumptions of existing classical laminar boundary layer theories for forced and natural convection are violated, such as the strict two-dimensionality of the dynamics or the steadiness of the fluid motion. The boundary layer dynamics consists of two essential local dynamical building blocks, a plume detachment and a post-plume phase. The former is associated with larger variations of the instantaneous thickness of velocity and temperature boundary layer and a fully three-dimensional local flow. The post-plume dynamics is connected with the large-scale circulation in the cell that penetrates the boundary region from above. The mean turbulence profiles taken in localized sections of the boundary layer for each dynamical phase are also compared with solutions of perturbation expansions of the boundary layer equations of forced or natural convection towards mixed convection. Our analysis of both boundary layers shows that the near-wall dynamics combines elements of forced Blasius-type and natural convection.


2010 ◽  
Vol 662 ◽  
pp. 409-446 ◽  
Author(s):  
G. SILANO ◽  
K. R. SREENIVASAN ◽  
R. VERZICCO

We summarize the results of an extensive campaign of direct numerical simulations of Rayleigh–Bénard convection at moderate and high Prandtl numbers (10−1 ≤ Pr ≤ 104) and moderate Rayleigh numbers (105 ≤ Ra ≤ 109). The computational domain is a cylindrical cell of aspect ratio Γ = 1/2, with the no-slip condition imposed on all boundaries. By scaling the numerical results, we find that the free-fall velocity should be multiplied by $1/\sqrt{{\it Pr}}$ in order to obtain a more appropriate representation of the large-scale velocity at high Pr. We investigate the Nusselt and the Reynolds number dependences on Ra and Pr, comparing the outcome with previous numerical and experimental results. Depending on Pr, we obtain different power laws of the Nusselt number with respect to Ra, ranging from Ra2/7 for Pr = 1 up to Ra0.31 for Pr = 103. The Nusselt number is independent of Pr. The Reynolds number scales as ${\it Re}\,{\sim}\,\sqrt{{\it Ra}}/{\it Pr}$, neglecting logarithmic corrections. We analyse the global and local features of viscous and thermal boundary layers and their scaling behaviours with respect to Ra and Pr, and with respect to the Reynolds and Péclet numbers. We find that the flow approaches a saturation state when Reynolds number decreases below the critical value, Res ≃ 40. The thermal-boundary-layer thickness increases slightly (instead of decreasing) when the Péclet number increases, because of the moderating influence of the viscous boundary layer. The simulated ranges of Ra and Pr contain steady, periodic and turbulent solutions. A rough estimate of the transition from the steady to the unsteady state is obtained by monitoring the time evolution of the system until it reaches stationary solutions. We find multiple solutions as long-term phenomena at Ra = 108 and Pr = 103, which, however, do not result in significantly different Nusselt numbers. One of these multiple solutions, even if stable over a long time interval, shows a break in the mid-plane symmetry of the temperature profile. We analyse the flow structures through the transitional phases by direct visualizations of the temperature and velocity fields. A wide variety of large-scale circulation and plume structures has been found. The single-roll circulation is characteristic only of the steady and periodic solutions. For other regimes at lower Pr, the mean flow generally consists of two opposite toroidal structures; at higher Pr, the flow is organized in the form of multi-jet structures, extending mostly in the vertical direction. At high Pr, plumes mainly detach from sheet-like structures. The signatures of different large-scale structures are generally well reflected in the data trends with respect to Ra, less in those with respect to Pr.


2009 ◽  
Vol 638 ◽  
pp. 383-400 ◽  
Author(s):  
ERIC BROWN ◽  
GUENTER AHLERS

In agreement with a recent experimental discovery by Xi et al. (Phys. Rev. Lett., vol. 102, 2009, paper no. 044503), we also find a sloshing mode in experiments on the large-scale circulation (LSC) of turbulent Rayleigh–Bénard convection in a cylindrical sample of aspect ratio one. The sloshing mode has the same frequency as the torsional oscillation discovered by Funfschilling & Ahlers (Phys. Rev. Lett., vol. 92, 2004, paper no. 1945022004). We show that both modes can be described by an extension of a model developed previously Brown & Ahlers (Phys. Fluids, vol. 20, 2008, pp. 105105-1–105105-15; Phys. Fluids, vol. 20, 2008, pp. 075101-1–075101-16). The extension consists of permitting a lateral displacement of the LSC circulation plane away from the vertical centreline of the sample as well as a variation of the displacement with height (such displacements had been excluded in the original model). Pressure gradients produced by the sidewall of the container on average centre the plane of the LSC so that it prefers to reach its longest diameter. If the LSC is displaced away from this diameter, the walls provide a restoring force. Turbulent fluctuations drive the LSC away from the central alignment, and combined with the restoring force they lead to oscillations. These oscillations are advected along with the LSC. This model yields the correct wavenumber and phase of the oscillations, as well as estimates of the frequency, amplitude and probability distributions of the displacements.


2017 ◽  
Vol 832 ◽  
pp. 721-744 ◽  
Author(s):  
Věra Musilová ◽  
Tomáš Králík ◽  
Marco La Mantia ◽  
Michal Macek ◽  
Pavel Urban ◽  
...  

We perform an experimental study of turbulent Rayleigh–Bénard convection up to very high Rayleigh number, $10^{8}<Ra<10^{14}$, in a cylindrical aspect ratio one cell, 30 cm in height, filled with cryogenic helium gas. We monitor temperature fluctuations in the convective flow with four small (0.2 mm) sensors positioned in pairs 1.5 cm from the sidewalls and 2.5 cm vertically apart and symmetrically around the mid-height of the cell. Based on one-point and two-point correlations of the temperature fluctuations, we determine different types of Reynolds numbers, $\mathit{Re}$, associated with the large-scale circulation (LSC). We observe a transition between two types of $\mathit{Re}(\mathit{Ra})$ scaling around $\mathit{Ra}=10^{10}{-}10^{11}$, which is accompanied by a scaling change of the skewness of the probability distribution functions (PDFs) of the temperature fluctuations. The $\mathit{Re}(\mathit{Ra})$ dependencies measured near the sidewall at Prandtl number $\mathit{Pr}\sim 1$ are consistent with the $\mathit{Ra}^{4/9}\mathit{Pr}^{-2/3}$ scaling above the transition, while for $\mathit{Ra}<10^{10}$, the $\mathit{Re}(\mathit{Ra})$ dependencies are steeper. It seems likely that this change in $\mathit{Re}(\mathit{Ra})$ scaling is linked to the previously reported change in the Nusselt number $\mathit{Nu}(\mathit{Ra})$ scaling. This behaviour is in agreement with independent cryogenic laboratory experiments with $\mathit{Pr}\sim 1$, but markedly different from the $\mathit{Re}$ scaling obtained in water experiments ($\mathit{Pr}\sim 3.3{-}5.6$). We discuss the results in comparison with different versions of the Grossmann–Lohse theory.


2015 ◽  
Vol 778 ◽  
Author(s):  
Jin-Qiang Zhong ◽  
Sebastian Sterl ◽  
Hui-Min Li

We present measurements of the azimuthal rotation velocity $\dot{{\it\theta}}(t)$ and thermal amplitude ${\it\delta}(t)$ of the large-scale circulation in turbulent Rayleigh–Bénard convection with modulated rotation. Both $\dot{{\it\theta}}(t)$ and ${\it\delta}(t)$ exhibit clear oscillations at the modulation frequency ${\it\omega}$. Fluid acceleration driven by oscillating Coriolis force causes an increasing phase lag in $\dot{{\it\theta}}(t)$ when ${\it\omega}$ increases. The applied modulation produces oscillatory boundary layers and the resulting time-varying viscous drag modifies ${\it\delta}(t)$ periodically. Oscillation of $\dot{{\it\theta}}(t)$ with maximum amplitude occurs at a finite modulation frequency ${\it\omega}^{\ast }$. Such a resonance-like phenomenon is interpreted as a result of optimal coupling of ${\it\delta}(t)$ to the modulated rotation velocity. We show that an extended large-scale circulation model with a relaxation time for ${\it\delta}(t)$ in response to the modulated rotation provides predictions in close agreement with the experimental results.


Sign in / Sign up

Export Citation Format

Share Document