Scaling of the velocity power spectra in turbulent thermal convection

Applicability of Taylor's hypothesis in thermally driven turbulence

In this paper, we show that, in the presence of large-scale circulation (LSC), Taylor's hypothesis can be invoked to deduce the energy spectrum in thermal convection using real-space probes, a popular experimental tool. We perform numerical simulation of turbulent convection in a cube and observe that the velocity field follows Kolmogorov's spectrum (k^-5/3). We also record the velocity time series using real-space probes near the lateral walls. The corresponding frequency spectrum exhibits Kolmogorov's spectrum (f^-5/3), thus validating Taylor's hypothesis with the steady LSC playing the role of a mean velocity field. The aforementioned findings based on real-space probes provide valuable inputs for experimental measurements used for studying the spectrum of convective turbulence.

Effect of polymer additives on heat transport and large-scale circulation in turbulent Rayleigh-Bénard convection

The present paper presents direct numerical simulations of Rayleigh-Bénard convection (RBC) in an enclosed cell filled with the polymer solution in order to investigate the viscoelastic effect on the characteristics of heat transport and large-scale circulation (LSC) of RBC. To overcome the difficulties in numerically solving a high Weissenberg number (Wi) problem of viscoelastic fluid flow with strong elastic effect, the log-conformation reformulation method was implemented. Numerical results showed that the addition of polymers reduced the heat flux and the amount of heat transfer reduction (HTR) behaves nonmonotonically, which firstly increases but then decreases with Wi. The maximum HTR reaches around 8.7% at the critical Wi. The nonmonotonic behavior of HTR as a function of Wi was then corroborated with the modifications of the period of LSC and turbulent energy as well as viscous boundary layer thickness. Finally, a standard turbulent kinetic energy (TKE) budget analysis was done for the whole domain, the boundary layer region, and the bulk region. It showed that the role change of elastic stress contributions to TKE is mainly responsible for this nonmonotonic behavior of HTR.

Probing the energy cascade of convective turbulence

The existence of a buoyancy-dominated scaling range in convective turbulence is a longstanding open question. We investigate this issue by considering the scale-by-scale energy budget in direct numerical simulations of Rayleigh-Bénard convection. We try to minimize the so-called Bolgiano length scale, the length scale at which buoyancy becomes dominant for scaling. Therefore, we deliberately choose modest Rayleigh numbers Ra=2.5×10(6) and 2.5×10(7). The budget reveals that buoyant forcing, turbulent energy transfer, and dissipation are contributing significantly over a wide range of scales. Thereby neither Kolmogorov-like (balance of turbulent transfer and dissipation) nor Bolgiano-Obukhov-like scaling (balance of turbulent transfer and buoyancy) is expected in the structure functions, which indeed reveal inconclusive scaling behavior. Furthermore, we consider the calculation of the Bolgiano length scale. To account for correlations between the dissipation rates of kinetic energy and thermal variance we propose to average the Bolgiano length scale directly. This gives an estimate, which is one order of magnitude larger than the previous estimate, and actually larger than the domain itself. Rather than studying the scaling of structure functions, we propose that the use of scale-by-scale energy budgets resolving anisotropic contributions is appropriate to consider the energy cascade mechanisms in turbulent convection.

Energy spectrum of buoyancy-driven turbulence

Using high-resolution direct numerical simulation and arguments based on the kinetic energy flux Π(u), we demonstrate that, for stably stratified flows, the kinetic energy spectrum E(u)(k)∼k(-11/5), the potential energy spectrum E(θ)(k)∼k(-7/5), and Π(u)(k)∼k(-4/5) are consistent with the Bolgiano-Obukhov scaling. This scaling arises due to the conversion of kinetic energy to the potential energy by buoyancy. For weaker buoyancy, this conversion is weak, hence E(u)(k) follows Kolmogorov's spectrum with a constant energy flux. For Rayleigh-Bénard convection, we show that the energy supply rate by buoyancy is positive, which leads to an increasing Π(u)(k) with k, thus ruling out Bolgiano-Obukhov scaling for the convective turbulence. Our numerical results show that convective turbulence for unit Prandt number exhibits a constant Π(u)(k) and E(u)(k)∼k(-5/3) for a narrow band of wave numbers.

Entropy and energy spectra in low-Prandtl-number convection with rotation

We present results for entropy and kinetic energy spectra computed from direct numerical simulations for low-Prandtl-number (Pr < 1) turbulent flow in Rayleigh-Bénard convection with uniform rotation about a vertical axis. The simulations are performed in a three-dimensional periodic box for a range of the Taylor number (0 ≤ Ta ≤ 10(8)) and reduced Rayleigh number r = Ra/Ra(∘)(Ta,Pr) (1.0 × 10(2) ≤ r ≤ 5.0 × 10(3)). The Rossby number Ro varies in the range 1.34 ≤ Ro ≤ 73. The entropy spectrum E(θ)(k) shows bisplitting into two branches for lower values of wave number k. The entropy in the lower branch scales with k as k(-1.4 ± 0.1) for r>10(3) for the rotation rates considered here. The entropy in the upper branch also shows scaling behavior with k, but the scaling exponent decreases with increasing Ta for all r. The energy spectrum E(v)(k) is also found to scale with the wave number k as k(-1.4 ± 0.1) for r>10(3). The scaling exponent for the energy spectrum and the lower branch of the entropy spectrum vary between -1.7 and -2.4 for lower values of r (<10(3)). We also provide some simple arguments based on the variation of the Kolmogorov picture to support the results of simulations.

Enhanced and reduced heat transport in turbulent thermal convection with polymer additives

We present an experimental study of turbulent Rayleigh-Bénard convection with polymer additives made in two convection cells, one with a smooth top and bottom plates and the other with a rough top and bottom plates. For the cell with smooth plates, a reduction of the measured Nusselt number (Nu) was observed. Furthermore, the amount of Nu reduction increases with increasing polymer concentration (c), reaching ~12% for c = 120 ppm and an apparent leveling off thereafter. For the cell with rough plates, however, an enhancement (~4%) of Nu was observed when the polymer concentration is greater than 120 ppm. This increase in Nu is corroborated by an increased large-scale circulation (LSC) velocity in the same cell when polymers are added. In contrast, the LSC velocity in the smooth cell is found to be essentially the same with and without polymers. It is further found that in the smooth cell the rms values of the global Nu, σ(Nu), and that of the local temperature, σ(T), both exhibit similar dependence on c as Nu itself. In contrast, σ(Nu) and σ(T) in the rough cell are found to be essentially independent of c.

Sonic anemometry to measure natural ventilation in greenhouses

The present work has developed a methodology for studying natural ventilation in Mediterranean greenhouses by means of sonic anemometry. In addition, specific calculation programmes have been designed to enable processing and analysis of the data recorded during the experiments. Sonic anemometry allows us to study the direction of the airflow at all the greenhouse vents. Knowing through which vents the air enters and leaves the greenhouse enables us to establish the airflow pattern of the greenhouse under natural ventilation conditions. In the greenhouse analysed in this work for Poniente wind (from the southwest), a roof vent designed to open towards the North (leeward) could allow a positive interaction between the wind and stack effects, improving the ventilation capacity of the greenhouse. The cooling effect produced by the mass of turbulent air oscillating between inside and outside the greenhouse at the side vents was limited to 2% (for high wind speed, u_o ≥ 4 m s⁻¹) reaching 36.3% when wind speed was lower (u_o = 2 m s⁻¹).

Energy spectra and fluxes for Rayleigh-Bénard convection

We compute the spectra and fluxes of the velocity and temperature fields in Rayleigh-Bénard convection in turbulent regime for a wide range of Prandtl numbers using pseudospectral simulations on 512(3) grids. Our spectral and flux results support the Kolmogorov-Obukhov (KO) scaling for zero Prandtl number and low Prandtl number (P=0.02) convection. The KO scaling for the velocity field in zero-Prandtl number and low-Prandtl number convection is because of the weak buoyancy in the inertial range (buoyancy is active only at the very low wave numbers). We also observe that for intermediate Prandtl numbers (P=0.2) the KO scaling fits better with the numerical results than the Bolgiano-Obukhov (BO) scaling. For large Prandtl number (P=6.8) , the spectra and flux results are somewhat inconclusive on the validity of the KO or BO scaling, yet the BO scaling is preferred over the KO scaling for these cases. The numerical results for P=1 is rather inconclusive.

Origin of the temperature oscillation in turbulent thermal convection

We report an experimental study of the three-dimensional spatial structure of the low-frequency temperature oscillations in a cylindrical Rayleigh-Bénard convection cell. Through simultaneous multipoint temperature measurements it is found that, contrary to the popular scenario, thermal plumes are emitted neither periodically nor alternately, but randomly and continuously, from the top and bottom plates. We further identify a new flow mode-the sloshing mode of the large-scale circulation (LSC). This sloshing mode, together with the torsional mode of the LSC, are found to be the origin of the oscillation of the temperature field.

Azimuthal motion, reorientation, cessation, and reversal of the large-scale circulation in turbulent thermal convection: a comparative study in aspect ratio one and one-half geometries

We report a systematic experimental study of the orientation and the flow strength of the large-scale circulation (LSC) in water-filled cylindrical Rayleigh-Bénard convection cells with aspect ratios 2.3, 1, and 0.5 by both direct velocity measurement and the indirect multithermal-probe measurement. Unlike its weak effect in the system's global heat transport, the aspect ratio Gamma is found to play an important role in the dynamics of the azimuthal motion of the LSC. It is found that in larger Gamma geometries the azimuthal motion of the LSC's vertical plane is confined in smaller azimuthal region than that in smaller Gamma geometries. The twisting motion between top and bottom parts of the LSC observed in the Gamma=1 geometry is found to be absent in the Gamma=1/2 case. It is found that in the Gamma=1/2 geometry the orientational change mid R:Deltavarphimid R: through a reorientation has an exponential distribution, in contrast to the power-law distribution for the Gamma=1 case. Despite the difference in orientational change, the occurrence of the reorientations is a Poisson process in both geometries. Using the conditional average of the time interval between adjacent cessations or reversals on the rebound flow strength, we demonstrate the possibility to empirically predict when the next cessation or reversal will most likely occur if the rebound flow strength of the preceding cessation or reversal is given.

Comparative experimental study of local mixing of active and passive scalars in turbulent thermal convection

We investigate experimentally the statistical properties of active and passive scalar fields in turbulent Rayleigh-Bénard convection in water, at Ra approximately 10;{10} . Both the local concentration of fluorescence dye and the local temperature are measured near the sidewall of a rectangular cell. It is found that, although they are advected by the same turbulent flow, the two scalars distribute differently. This difference is twofold, i.e., both the quantities themselves and their small-scale increments have different distributions. Our results show that there is a certain buoyant scale based on time domain, i.e., the Bolgiano time scale t_{B} , above which buoyancy effects are significant. Above t_{B} , temperature is active and is found to be more intermittent than concentration, which is passive. This suggests that the active scalar possesses a higher level of intermittency in turbulent thermal convection. It is further found that the mixing of both scalar fields are isotropic for scales larger than t_{B} even though buoyancy acts on the fluid in the vertical direction. Below t_{B} , temperature is passive and is found to be more anisotropic than concentration. But this higher degree of anisotropy is attributed to the higher diffusivity of temperature over that of concentration. From the simultaneous measurements of temperature and concentration, it is shown that two scalars have similar autocorrelation functions and there is a strong and positive correlation between them.

Numerical and experimental investigation of structure-function scaling in turbulent Rayleigh-Bénard convection

Direct numerical simulation and stereoscopic particle image velocimetry of turbulent convection are used to gather spatial data for the calculation of structure functions. We wish to add to the ongoing discussion in the literature whether temperature acts as an active or passive scalar in turbulent convection, with consequences for structure-function scaling. The simulation results show direct confirmation of the scalings derived by Bolgiano and Obukhov for turbulence with an active scalar for both velocity and temperature statistics. The active-scalar range shifts to larger scales when the forcing parameter (Rayleigh number) is increased. Furthermore, a close inspection of local turbulent length scales (Kolmogorov and Bolgiano lengths) confirms conjectures from earlier studies that the oft-used global averages are not suited for the interpretation of structure functions. In the experiment, a characterization of the domain-filling large-scale circulation of confined convection is carried out for comparison with other studies. The measured velocity fields are also used to calculate velocity structure functions, further confirming the Bolgiano-Obukhov scalings when interpreted with the local turbulent length scales found in the simulations. An extended self-similarity analysis shows that the relative scalings are different for the Kolmogorov and Bolgiano-Obukhov regimes.

Measured oscillations of the velocity and temperature fields in turbulent Rayleigh-Bénard convection in a rectangular cell

Temperature and velocity oscillations have been found in a rectangular Rayleigh-Bénard convection cell, in which one large-scale convection roll exists. At Rayleigh number Ra=8.9x10(11) and Prandtl number Pr=4, temperature oscillations can be observed in most parts of the system and the oscillation period remains almost constant, tT=74+/-2 s. Velocity oscillation can only be found in its horizontal component vy (perpendicular to the large-scale circulation plane) near the cell sidewall, its oscillation period is also constant, tv=65+/-2 s, at these positions. Temperature and velocity oscillations have different Ra dependences, which are, respectively, indicated by the Péclect number PeT=0.55Ra0.47 and Pev=0.28Ra0.50. In comparison to the case of a cylindrical cell, we find that velocity oscillations are affected by the system geometry.

Scaling laws in the central region of confined turbulent thermal convection

In confined turbulent thermal convection, the velocity is separated into two parts: one that is correlated with some function of the temperature fluctuations, and thus associated with the plume velocity, and the other part, the background velocity, which is uncorrelated with any function of the temperature fluctuations. As a result, one should focus on the plume velocity, and not the whole velocity, and the temperature when studying the scaling behavior. In this paper, a phenomenological theory for the scaling behavior in the central region of confined turbulent thermal convection is presented. The spatial (temporal) plume velocity structure functions are found to have the same scaling behavior as the spatial (temporal) temperature structure functions. For tau> or = taub, where the buoyant scale taub is determined in terms of measurable quantities, the scaling exponents of the temporal temperature structure functions and hence those of the temporal plume velocity structure functions are obtained. These results are checked against experimental measurements, and good agreement is found.

Rayleigh-Taylor turbulence is nothing like Kolmogorov turbulence in the self-similar regime

An increasing number of numerical simulations and experiments describing the turbulent spectrum of Rayleigh-Taylor (RT) mixing layers came to light over the past few years. Results reported in recent studies allow to rule out a Kolmogorov-like turbulence as a mechanism acting on a self-similar RT turbulent mixing layer. A different mechanism is presented, which complies with both numerical and experimental results and relates RT flow to other buoyant flows.

Cascades of velocity and temperature fluctuations in buoyancy-driven thermal turbulence

Direct multipoint measurements of the velocity and temperature fields have been made in a turbulent Rayleigh-Bénard convection cell. In the central region of the cell it is found that both velocity and temperature exhibit the same scaling behavior that one would find for the velocity and for a passive scalar in homogeneous and isotropic Navier-Stokes turbulence. This is despite the fact that energy is pumped into the system vertically via buoyancy. Near the cell's sidewall where thermal plumes abound, vertical velocity and temperature exhibit different scalings. A model that takes into account both buoyancy and energy dissipation is proposed and its predictions agree well with the sidewall experimental results.

Azimuthal motion of the mean wind in turbulent thermal convection

We present an experimental study of the azimuthal motion of the mean wind in turbulent thermal convection. The experiments were conducted with cylindrical convection cells of unity aspect ratio and over the range of the Rayleigh number from 1 x 10(9) to 1 x 10(10). The azimuthal angle of the circulation plane of the mean wind was measured using both the particle image velocimetry and flow-visualization techniques. It is found that the azimuthal motion consists of erratic fluctuations and a time-periodic oscillation. The orientation of the wind is found to be "locked," i.e., it fluctuates about a preferred direction most of the time with all other orientations appearing as "transient states," and large excursions of the azimuthal angle often result in a net rotation which takes the wind back to the preferred orientation. The rate of erratic rotation of the circulation plane is found to have a strong dependence on Ra. Our result suggests that the oscillatory motion of the wind in its vertically oriented circulation plane and the orientational oscillation of the circulation plane itself have the same dynamic origin.

Scaling of the Reynolds number in turbulent thermal convection

A riddle in turbulent thermal convection is the apparent dispersion from 0.42 to 0.5 in the value of the scaling exponent of experimentally measured Reynolds number Re approximately Ragamma, where Ra is the Rayleigh number. The measured Re may be divided into two groups: one based on the circulation frequency of the mean wind and the other based on a directly measured velocity. With new experimental results we show that in frequency measurements the dispersion in gamma is a result of the evolution in the circulation path of the wind, and that in the velocity measurements it is caused by the inclusion of a counterflow in the mean velocity. When these factors are properly accounted for both groups give gamma=0.5, which may imply that a single mechanism is driving the flow for both low and high values of Ra.

Three-dimensional flow structures and dynamics of turbulent thermal convection in a cylindrical cell

The technique of particle image velocimetry is used to study the velocity field of turbulent Rayleigh-Bénard convection in an aspect-ratio-1 cylindrical cell filled with water. By measuring the two-dimensional (2D) velocity vector map in different vertical cross sections of the cell, we investigate the 3D structures and dynamics of turbulent thermal convection. The experiment reveals how thermal plumes synchronize their emissions and organize their motions spatially between the top and bottom plates, leading to an oscillatory motion in the bulk region of the fluid with a period equal to twice the plume's cell-crossing time. From the measured instantaneous velocity vector map, we find the phase relationship between the velocity components along different directions and at different positions in a 2D plane. These phase relations illustrate how the convecting fluid in different regions of the cell interact with each other and generate a synchronized and coherent motion in a closed system.

23/9 dimensional anisotropic scaling of passive admixtures using lidar data of aerosols

In buoyancy-driven flows, another dimensional quantity appears in addition to the energy flux. Classically, this leads to the prediction that at large scales, isotropic Bolgiano-Obukhov (BO) scaling can dominate isotropic Kolmogorov scaling. We investigate this in the atmosphere by using state-of-the-art high-powered lidar data. We examine simultaneous horizontal and vertical sections of passive scalar surrogates over the ranges 100 m to 120 km and 3 m to 4.5 km , respectively. Overall, this spans the crucial "mesoscale" and involves nearly 1000 times more data than the largest relevant experiments to date. Rather than a transition from one isotropic regime to another, we find that the two regimes always coexist in an anisotropic Corrsin-Obukhov law with the Kolmogorov holding in the horizontal, and the BO holding in the vertical. The stratification is quantified by an elliptical dimension D(el) found to be equal to 2.55+/-0.02 . This anisotropic scaling is very close to that predicted by the 23/9 dimensional unified scaling model of the atmosphere and is consistent with observations of the horizontal wind.

Instantaneous measurement of velocity fields in developed thermal turbulence in mercury

Using ultrasonic velocimetry we measured the vertical profile of the velocity fluctuation in high-Rayleigh-number thermal convection in a cell with aspect ratio of 0.5, filled with a low-Prandtl-number fluid, mercury. The intriguing fluctuating dynamics of the mean flow and universal nature of the kinetic energy cascade are elucidated utilizing spectral decomposition and reconstruction. The scaling properties of the structure functions and the energy spectrum are directly calculated without the use of Taylor's frozen-flow hypothesis. Despite the complex nature of the mean flow, it is found that the energy cascade process exhibits universal laws in thermal turbulence.

Temperature oscillations in turbulent Rayleigh-Bénard convection

A systematic study of temperature oscillations in turbulent thermal convection was carried out in two aspect-ratio-one convection cells filled with water. Temperature correlation functions and local velocity fluctuations were measured over varying Rayleigh numbers and spatial positions across the entire cell. These measurements fully characterize the spatial structure of the temperature oscillation and reveal the mixing and emission dynamics of the thermal plumes near the conducting surface. A sharp transition from a random chaotic state to a correlated turbulent state of finite coherence time is found when the Rayleigh number becomes larger than a critical value Ra(c) approximately equal 5 x 10(7). Above Ra(c) the measured temperature correlation functions show a well-defined oscillation with a finite coherence time. The oscillation period is found to be twice as large as the cell crossing time. The experiment demonstrates how the thermal plumes in a closed cell organize themselves both in space and time and generate coherent oscillations in a turbulent environment.

Prandtl number dependence of the viscous boundary layer and the Reynolds numbers in Rayleigh-Bénard convection

We report results from high Prandtl number turbulent thermal convection experiments. The viscous boundary layer and the Reynolds number are measured in four different fluids over wide ranges of the Prandtl number Pr and the Rayleigh number Ra, all in a single convection cell of unity aspect ratio. We find that the normalized viscous layer thickness may be represented as delta(v)/L=0.65Pr(0.24)Ra(-0.16). The Reynolds number based on the oscillation frequency of the large-scale flow is found as Re(o)(Ra,Pr)=1.1Ra(0.43)Pr(-0.76) and that based on the rms velocity Re(rms)(Ra,Pr)=0.84Ra(0.40)Pr-0.86. Both the Ra and the Pr exponents of Re(V(m))(Ra,Pr) based on the maximum velocity of the circulating wind appear to vary across the range of Pr covered, changing from 0.5 to 0.68 and -0.88 to -0.95, respectively, as Pr is increased from 6 to 1027.