Vertical structure of conventionally neutral atmospheric boundary layers.

SignificanceThe presented model describes the vertical structure of conventionally neutral atmospheric boundary layers. Due to the complicated interplay between buoyancy, shear, and Coriolis effects, analytical descriptions have been limited to the mean wind speed. We introduce an analytical approach based on the Ekman equations and the basis function of the universal potential temperature flux profile that allows one to describe the wind and turbulent shear stress profiles and hence capture features like the wind veer profile. The analytical profiles are validated against high-fidelity large-eddy simulations and atmospheric measurements. Our findings contribute to the scientific community's fundamental understanding of atmospheric turbulence with direct relevance for weather forecasting, climate modeling, and wind energy applications.

Aspect Ratio Dependence of Heat Transfer in a Cylindrical Rayleigh-Bénard Cell.

While the heat transfer and the flow dynamics in a cylindrical Rayleigh-Bénard (RB) cell are rather independent of the aspect ratio Γ (diameter/height) for large Γ, a small-Γ cell considerably stabilizes the flow and thus affects the heat transfer. Here, we first theoretically and numerically show that the critical Rayleigh number for the onset of convection at given Γ follows Ra_{c,Γ}∼Ra_{c,∞}(1+CΓ^{-2})^{2}, with C≲1.49 for Oberbeck-Boussinesq (OB) conditions. We then show that, in a broad aspect ratio range (1/32)≤Γ≤32, the rescaling Ra→Ra_{ℓ}≡Ra[Γ^{2}/(C+Γ^{2})]^{3/2} collapses various OB numerical and almost-OB experimental heat transport data Nu(Ra,Γ). Our findings predict the Γ dependence of the onset of the ultimate regime Ra_{u,Γ}∼[Γ^{2}/(C+Γ^{2})]^{-3/2} in the OB case. This prediction is consistent with almost-OB experimental results (which only exist for Γ=1, 1/2, and 1/3) for the transition in OB RB convection and explains why, in small-Γ cells, much larger Ra (namely, by a factor Γ^{-3}) must be achieved to observe the ultimate regime.

Universal Wind Profile for Conventionally Neutral Atmospheric Boundary Layers.

Conventionally neutral atmospheric boundary layers (CNBLs), which are characterized with zero surface potential temperature flux and capped by an inversion of potential temperature, are frequently encountered in nature. Therefore, predicting the wind speed profiles of CNBLs is relevant for weather forecasting, climate modeling, and wind energy applications. However, previous attempts to predict the velocity profiles in CNBLs have had limited success due to the complicated interplay between buoyancy, shear, and Coriolis effects. Here, we utilize ideas from the classical Monin-Obukhov similarity theory in combination with a local scaling hypothesis to derive an analytic expression for the stability correction function ψ=-c_{ψ}(z/L)^{1/2}, where c_{ψ}=4.2 is an empirical constant, z is the height above ground, and L is the local Obukhov length based on potential temperature flux at that height, for CNBLs. An analytic expression for this flux is also derived using dimensional analysis and a perturbation method approach. We find that the derived profile agrees excellently with the velocity profile in the entire boundary layer obtained from high-fidelity large eddy simulations of typical CNBLs.

Diffusion-Free Scaling in Rotating Spherical Rayleigh-Bénard Convection.

Direct numerical simulations are employed to reveal three distinctly different flow regions in rotating spherical Rayleigh-Bénard convection. In the high-latitude region I vertical (parallel to the axis of rotation) convective columns are generated between the hot inner and the cold outer sphere. The mid-latitude region I I is dominated by vertically aligned convective columns formed between the Northern and Southern hemispheres of the outer sphere. The diffusion-free scaling, which indicates bulk-dominated convection, originates from this mid-latitude region. In the equator region I I I , the vortices are affected by the outer spherical boundary and are much shorter than in region I I .

Transition to the Ultimate Regime in Two-Dimensional Rayleigh-Bénard Convection.

The possible transition to the so-called ultimate regime, wherein both the bulk and the boundary layers are turbulent, has been an outstanding issue in thermal convection, since the seminal work by Kraichnan [Phys. Fluids 5, 1374 (1962)PFLDAS0031-917110.1063/1.1706533]. Yet, when this transition takes place and how the local flow induces it is not fully understood. Here, by performing two-dimensional simulations of Rayleigh-Bénard turbulence covering six decades in Rayleigh number Ra up to 10^{14} for Prandtl number Pr=1, for the first time in numerical simulations we find the transition to the ultimate regime, namely, at Ra^{*}=10^{13}. We reveal how the emission of thermal plumes enhances the global heat transport, leading to a steeper increase of the Nusselt number than the classical Malkus scaling Nu∼Ra^{1/3} [Proc. R. Soc. A 225, 196 (1954)PRLAAZ1364-502110.1098/rspa.1954.0197]. Beyond the transition, the mean velocity profiles are logarithmic throughout, indicating turbulent boundary layers. In contrast, the temperature profiles are only locally logarithmic, namely, within the regions where plumes are emitted, and where the local Nusselt number has an effective scaling Nu∼Ra^{0.38}, corresponding to the effective scaling in the ultimate regime.

Roughness-Facilitated Local 1/2 Scaling Does Not Imply the Onset of the Ultimate Regime of Thermal Convection.

In thermal convection, roughness is often used as a means to enhance heat transport, expressed in Nusselt number. Yet there is no consensus on whether the Nusselt vs Rayleigh number scaling exponent (Nu∼Ra^{ß}) increases or remains unchanged. Here we numerically investigate turbulent Rayleigh-Bénard convection over rough plates in two dimensions, up to Ra≈10^{12}. Varying the height and wavelength of the roughness elements with over 200 combinations, we reveal the existence of two universal regimes. In the first regime, the local effective scaling exponent can reach up to 1/2. However, this cannot be explained as the attainment of the so-called ultimate regime as suggested in previous studies, because a further increase in Ra leads to the second regime, in which the scaling saturates back to a value close to the smooth wall case. Counterintuitively, the transition from the first to the second regime corresponds to the competition between bulk and boundary layer flow: from the bulk-dominated regime back to the classical boundary-layer-controlled regime. Our study demonstrates that the local 1/2 scaling does not necessarily signal the onset of ultimate turbulence.

Confined Rayleigh-Bénard, Rotating Rayleigh-Bénard, and Double Diffusive Convection: A Unifying View on Turbulent Transport Enhancement through Coherent Structure Manipulation.

Many natural and engineering systems are simultaneously subjected to a driving force and a stabilizing force. The interplay between the two forces, especially for highly nonlinear systems such as fluid flow, often results in surprising features. Here we reveal such features in three different types of Rayleigh-Bénard (RB) convection, i.e., buoyancy-driven flow with the fluid density being affected by a scalar field. In the three cases different stabilizing forces are considered, namely (i) horizontal confinement, (ii) rotation around a vertical axis, and (iii) a second stabilizing scalar field. Despite the very different nature of the stabilizing forces and the corresponding equations of motion, at moderate strength we counterintuitively but consistently observe an enhancement in the flux, even though the flow motion is weaker than the original RB flow. The flux enhancement occurs in an intermediate regime in which the stabilizing force is strong enough to alter the flow structures in the bulk to a more organized morphology, yet not too strong to severely suppress the flow motions. Near the optimal transport enhancements all three systems exhibit a transition from a state in which the thermal boundary layer (BL) is nested inside the momentum BL to the one with the thermal BL being thicker than the momentum BL. The observed optimal transport enhancement is explained through an optimal coupling between the suction of hot or fresh fluid and the corresponding scalar fluctuations.

Heat transport in bubbling turbulent convection.

Boiling is an extremely effective way to promote heat transfer from a hot surface to a liquid due to numerous mechanisms, many of which are not understood in quantitative detail. An important component of the overall process is that the buoyancy of the bubble compounds with that of the liquid to give rise to a much-enhanced natural convection. In this article, we focus specifically on this enhancement and present a numerical study of the resulting two-phase Rayleigh-Bénard convection process in a cylindrical cell with a diameter equal to its height. We make no attempt to model other aspects of the boiling process such as bubble nucleation and detachment. The cell base and top are held at temperatures above and below the boiling point of the liquid, respectively. By keeping this difference constant, we study the effect of the liquid superheat in a Rayleigh number range that, in the absence of boiling, would be between 2 × 10(6) and 5 × 10(9). We find a considerable enhancement of the heat transfer and study its dependence on the number of bubbles, the degree of superheat of the hot cell bottom, and the Rayleigh number. The increased buoyancy provided by the bubbles leads to more energetic hot plumes detaching from the cell bottom, and the strength of the circulation in the cell is significantly increased. Our results are in general agreement with recent experiments on boiling Rayleigh-Bénard convection.

Breakdown of the large-scale circulation in Γ=1/2 rotating Rayleigh-Bénard flow.

Experiments and simulations of rotating Rayleigh-Bénard convection in cylindrical samples have revealed an increase in heat transport with increasing rotation rate. This heat transport enhancement is intimately related to a transition in the turbulent flow structure from a regime dominated by a large-scale circulation (LSC), consisting of a single convection roll, at no or weak rotation to a regime dominated by vertically aligned vortices at strong rotation. For a sample with an aspect ratio Γ=D/L=1 (D is the sample diameter and L is its height) the transition between the two regimes is indicated by a strong decrease in the LSC strength. In contrast, for Γ=1/2, Weiss and Ahlers [J. Fluid Mech. 688, 461 (2011)] revealed the presence of a LSC-like sidewall temperature signature beyond the critical rotation rate. They suggested that this might be due to the formation of a two-vortex state, in which one vortex extends vertically from the bottom into the sample interior and brings up warm fluid while another vortex brings down cold fluid from the top; this flow field would yield a sidewall temperature signature similar to that of the LSC. Here we show by direct numerical simulations for Γ=1/2 and parameters that allow direct comparison with experiment that the spatial organization of the vertically aligned vortical structures in the convection cell do indeed yield (for the time average) a sinusoidal variation of the temperature near the sidewall, as found in the experiment. This is also the essential and nontrivial difference with the Γ=1 sample, where the vertically aligned vortices are distributed randomly.

Spatial distribution of heat flux and fluctuations in turbulent Rayleigh-Bénard convection.

We numerically investigate the radial dependence of the velocity and temperature fluctuations and of the time-averaged heat flux j ¯(r) in a cylindrical Rayleigh-Bénard cell with aspect ratio Γ=1 for Rayleigh numbers Ra between 2×10^{6} and 2×10^{9} at a fixed Prandtl number Pr=5.2. The numerical results reveal that the heat flux close to the sidewall is larger than in the center and that, just as the global heat transport, it has an effective power law dependence on the Rayleigh number, j ¯(r)∝Ra{γ{j}(r)}. The scaling exponent γ{j}(r) decreases monotonically from 0.43 near the axis (r≈0) to 0.29 close to the sidewalls (r≈D/2). The effective exponents near the axis and the sidewall agree well with the measurements of Shang et al. [Phys. Rev. Lett. 100, 244503 (2008)] and the predictions of Grossmann and Lohse [Phys. Fluids 16, 1070 (2004)]. Extrapolating our results to large Rayleigh number would imply a crossover at Ra≈10^{15}, where the heat flux near the axis would begin to dominate. In addition, we find that the local heat flux is more than twice as high at the location where warm or cold plumes go up or down than in plume depleted regions.

Logarithmic temperature profiles in turbulent Rayleigh-Bénard convection.

We report results for the temperature profiles of turbulent Rayleigh-Bénard convection (RBC) in the interior of a cylindrical sample of aspect ratio Γ≡D/L=0.50 (D and L are the diameter and height, respectively). Both in the classical and in the ultimate state of RBC we find that the temperature varies as A×ln(z/L)+B, where z is the distance from the bottom or top plate. In the classical state, the coefficient A decreases in the radial direction as the distance from the side wall increases. For the ultimate state, the radial dependence of A has not yet been determined. These findings are based on experimental measurements over the Rayleigh-number range 4×10(12)≲Ra≲10(15) for a Prandtl number Pr≃0.8 and on direct numerical simulation at Ra=2×10(12), 2×10(11), and 2×10(10), all for Pr=0.7.

Thermal boundary layer profiles in turbulent Rayleigh-Bénard convection in a cylindrical sample.

We numerically investigate the structures of the near-plate temperature profiles close to the bottom and top plates of turbulent Rayleigh-Bénard flow in a cylindrical sample at Rayleigh numbers Ra = 10(8) to Ra = 2 × 10(12) and Prandtl numbers Pr = 6.4 and Pr = 0.7 with the dynamical frame method [Zhou and Xia, Phys. Rev. Lett. 104, 104301 (2010)], thus extending previous results for quasi-two-dimensional systems to three-dimensional systems. The dynamical frame method shows that the measured temperature profiles in the spatially and temporally local frame are much closer to the temperature profile of a laminar, zero-pressure gradient boundary layer (BL) according to Pohlhausen than in the fixed reference frame. The deviation between the measured profiles in the dynamical reference frame and the laminar profiles increases with decreasing Pr, where the thermal BL is more exposed to the bulk fluctuations due to the thinner kinetic BL, and increasing Ra, where more plumes are passing the measurement location.

Connecting flow structures and heat flux in turbulent Rayleigh-Bénard convection.

The aspect ratio (Γ) dependence of the heat transfer (Nusselt number Nu in dimensionless form) in turbulent (two-dimensional) Rayleigh-Bénard convection is numerically studied in the regime 0.4≤Γ≤1.25 for Rayleigh numbers 10(7)≤Ra≤Ra(9) and Prandtl numbers Pr=0.7 (gas) and 4.3 (water). Nu(Γ) shows a very rich structure with sudden jumps and sharp transitions. We connect these structures to the way the flow organizes itself in the sample and explain why the aspect ratio dependence of Nu is more pronounced for small Pr. Even for fixed Γ different turbulent states (with different resulting Nu) can exist, between which the flow can or cannot switch. In the latter case the heat transfer thus depends on the initial conditions.

Effect of aspect ratio on vortex distribution and heat transfer in rotating Rayleigh-Bénard convection.

Numerical and experimental data for the heat transfer as a function of the Rossby number Ro in turbulent rotating Rayleigh-Bénard convection are presented for the Prandtl number Pr=4.38 and the Rayleigh number Ra=2.91×10(8) up to Ra=4.52×10(9). The aspect ratio Γ≡D/L, where L is the height and D the diameter of the cylindrical sample, is varied between Γ=0.5 and 2.0. Without rotation, where the aspect ratio influences the global large-scale circulation, we see a small-aspect-ratio dependence in the Nusselt number for Ra=2.91×10(8). However, for stronger rotation, i.e., 1/Ro>>1/Ro(c), the heat transport becomes independent of the aspect ratio. We interpret this finding as follows: In the rotating regime the heat is mainly transported by vertically aligned vortices. Since the vertically aligned vortices are local, the aspect ratio has a negligible effect on the heat transport in the rotating regime. Indeed, a detailed analysis of vortex statistics shows that the fraction of the horizontal area that is covered by vortices is independent of the aspect ratio when 1/Ro>>1/Ro(c). In agreement with the results of Weiss et al. [Phys. Rev. Lett. 105, 224501 (2010)], we find a vortex-depleted area close to the sidewall. Here we show that there is also an area with enhanced vortex concentration next to the vortex-depleted edge region and that the absolute widths of both regions are independent of the aspect ratio.

Flow reversals in thermally driven turbulence.

We analyze the reversals of the large-scale flow in Rayleigh-Bénard convection both through particle image velocimetry flow visualization and direct numerical simulations of the underlying Boussinesq equations in a (quasi-) two-dimensional, rectangular 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. The Rayleigh vs Prandtl number space is mapped out. The occurrence of reversals sensitively depends on these parameters.

Finite-size effects lead to supercritical bifurcations in turbulent rotating Rayleigh-Bénard convection.

In turbulent thermal convection in cylindrical samples with an aspect ratio Γ≡D/L (D is the diameter and L the height), the Nusselt number Nu is enhanced when the sample is rotated about its vertical axis because of the formation of Ekman vortices that extract additional fluid out of thermal boundary layers at the top and bottom. We show from experiments and direct numerical simulations that the enhancement occurs only above a bifurcation point at a critical inverse Rossby number 1/Ro(c), with 1/Ro(c)∝1/Γ. We present a Ginzburg-Landau-like model that explains the existence of a bifurcation at finite 1/Ro(c) as a finite-size effect. The model yields the proportionality between 1/Ro(c) and 1/Γ and is consistent with several other measured or computed system properties.

Transitions between turbulent states in rotating Rayleigh-Bénard convection.

Weakly rotating turbulent Rayleigh-Bénard convection was studied experimentally and numerically. With increasing rotation and large enough Rayleigh number a supercritical bifurcation from a turbulent state with nearly rotation-independent heat transport to another with enhanced heat transfer is observed at a critical inverse Rossby number 1/Roc approximately 0.4. The strength of the large-scale convection roll is either enhanced or essentially unmodified depending on parameters for 1/Ro<1/Roc, but the strength increasingly diminishes beyond 1/Roc where it competes with Ekman vortices that cause vertical fluid transport and thus heat-transfer enhancement.

Prandtl-, Rayleigh-, and Rossby-number dependence of heat transport in turbulent rotating Rayleigh-Bénard convection.

Experimental and numerical data for the heat transfer as a function of the Rayleigh, Prandtl, and Rossby numbers in turbulent rotating Rayleigh-Bénard convection are presented. For relatively small Ra approximately 10(8) and large Pr modest rotation can enhance the heat transfer by up to 30%. At larger Ra there is less heat-transfer enhancement, and at small Pr less, similar 0.7 there is no heat-transfer enhancement at all. We suggest that the small-Pr behavior is due to the breakdown of the heat-transfer-enhancing Ekman pumping because of larger thermal diffusion.