Statistics of kinetic and thermal energy dissipation rates in two-dimensional turbulent Rayleigh–Bénard convection

2017 ◽  
Vol 814 ◽  
pp. 165-184 ◽  
Author(s):  
Yang Zhang ◽  
Quan Zhou ◽  
Chao Sun
Keyword(s):  
Boundary Layer ◽  
Energy Dissipation ◽  
Rayleigh Number ◽  
Small Scale ◽  
Two Dimensional ◽  
Number Range ◽  
Dissipation Rates ◽  
Log Normal ◽  
Rayleigh Benard ◽  
Gl Theory

We investigate the statistical properties of the kinetic $\unicode[STIX]{x1D700}_{u}$ and thermal $\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D703}}$ energy dissipation rates in two-dimensional (2-D) turbulent Rayleigh–Bénard (RB) convection. Direct numerical simulations were carried out in a box with unit aspect ratio in the Rayleigh number range $10^{6}\leqslant Ra\leqslant 10^{10}$ for Prandtl numbers $Pr=0.7$ and 5.3. The probability density functions (PDFs) of both dissipation rates are found to deviate significantly from a log-normal distribution. The PDF tails can be well described by a stretched exponential function, and become broader for higher Rayleigh number and lower Prandtl number, indicating an increasing degree of small-scale intermittency with increasing Reynolds number. Our results show that the ensemble averages $\langle \unicode[STIX]{x1D700}_{u}\rangle _{V,t}$ and $\langle \unicode[STIX]{x1D700}_{\unicode[STIX]{x1D703}}\rangle _{V,t}$ scale as $Ra^{-0.18\sim -0.20}$, which is in excellent agreement with the scaling estimated from the two global exact relations for the dissipation rates. By separating the bulk and boundary-layer contributions to the total dissipations, our results further reveal that $\langle \unicode[STIX]{x1D700}_{u}\rangle _{V,t}$ and $\langle \unicode[STIX]{x1D700}_{\unicode[STIX]{x1D703}}\rangle _{V,t}$ are both dominated by the boundary layers, corresponding to regimes $I_{l}$ and $I_{u}$ in the Grossmann–Lohse (GL) theory (J. Fluid Mech., vol. 407, 2000, pp. 27–56). To include the effects of thermal plumes, the plume–background partition is also considered and $\langle \unicode[STIX]{x1D700}_{\unicode[STIX]{x1D703}}\rangle _{V,t}$ is found to be plume dominated. Moreover, the boundary-layer/plume contributions scale as those predicted by the GL theory, while the deviations from the GL predictions are observed for the bulk/background contributions. The possible reasons for the deviations are discussed.

Observations of Small-Scale Turbulence and Energy Dissipation Rates in the Cloudy Boundary Layer

2006 ◽  
Vol 63 (5) ◽  
pp. 1451-1466 ◽  
Author(s):  
Holger Siebert ◽  
Katrin Lehmann ◽  
Manfred Wendisch
Keyword(s):  
Boundary Layer ◽  
Energy Dissipation ◽  
Wind Velocity ◽  
Turbulence Theory ◽  
Horizontal Wind ◽  
Inertial Subrange ◽  
Intermittent Flow ◽  
Small Scale ◽  
Dissipation Rates ◽  
Wide Range

Abstract Tethered balloon–borne measurements with a resolution in the order of 10 cm in a cloudy boundary layer are presented. Two examples sampled under different conditions concerning the clouds' stage of life are discussed. The hypothesis tested here is that basic ideas of classical turbulence theory in boundary layer clouds are valid even to the decimeter scale. Power spectral densities S( f ) of air temperature, liquid water content, and wind velocity components show an inertial subrange behavior down to ≈20 cm. The mean energy dissipation rates are ∼10−3 m2 s−3 for both datasets. Estimated Taylor Reynolds numbers (Reλ) are ∼104, which indicates the turbulence is fully developed. The ratios between longitudinal and transversal S( f ) converge to a value close to 4/3, which is predicted by classical turbulence theory for local isotropic conditions. Probability density functions (PDFs) of wind velocity increments Δu are derived. The PDFs show significant deviations from a Gaussian distribution with longer tails typical for an intermittent flow. Local energy dissipation rates ɛτ are derived from subsequences with a duration of τ = 1 s. With a mean horizontal wind velocity of 8 m s−1, τ corresponds to a spatial scale of 8 m. The PDFs of ɛτ can be well approximated with a lognormal distribution that agrees with classical theory. Maximum values of ɛτ ≈ 10−1 m2 s−3 are found in the analyzed clouds. The consequences of this wide range of ɛτ values for particle–turbulence interaction are discussed.

scholarly journals Prandtl–Blasius temperature and velocity boundary-layer profiles in turbulent Rayleigh–Bénard convection

2010 ◽  
Vol 664 ◽  
pp. 297-312 ◽  
Author(s):  
QUAN ZHOU ◽  
RICHARD J. A. M. STEVENS ◽  
KAZUYASU SUGIYAMA ◽  
SIEGFRIED GROSSMANN ◽  
DETLEF LOHSE ◽  
...  
Keyword(s):  
Boundary Layer ◽  
Rayleigh Number ◽  
Reference Frames ◽  
Temperature Profiles ◽  
Number Range ◽  
Velocity Boundary Layer ◽  
Bénard Convection ◽  
Benard Convection ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

The shapes of the velocity and temperature profiles near the horizontal conducting plates' centre regions in turbulent Rayleigh–Bénard convection are studied numerically and experimentally over the Rayleigh number range 108 ≲ Ra ≲ 3 × 1011 and the Prandtl number range 0.7 ≲ Pr ≲ 5.4. The results show that both the temperature and velocity profiles agree well with the classical Prandtl–Blasius (PB) laminar boundary-layer profiles, if they are re-sampled in the respective dynamical reference frames that fluctuate with the instantaneous thermal and velocity boundary-layer thicknesses. The study further shows that the PB boundary layer in turbulent thermal convection not only holds in a time-averaged sense, but is most of the time also valid in an instantaneous sense.

Torsional oscillations of the large-scale circulation in turbulent Rayleigh–Bénard convection

2008 ◽  
Vol 607 ◽  
pp. 119-139 ◽  
Author(s):  
DENIS FUNFSCHILLING ◽  
ERIC BROWN ◽  
GUENTER AHLERS
Keyword(s):  
Rayleigh Number ◽  
Aspect Ratio ◽  
Large Scale ◽  
Small Scale ◽  
Large Scale Circulation ◽  
Number Range ◽  
Bénard Convection ◽  
Benard Convection ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

Measurements over the Rayleigh-number range 108 ≲ R ≲ 1011 and Prandtl-number range 4.4≲σ≲29 that determine the torsional nature and amplitude of the oscillatory mode of the large-scale circulation (LSC) of turbulent Rayleigh–Bénard convection are presented. For cylindrical samples of aspect ratio Γ=1 the mode consists of an azimuthal twist of the near-vertical LSC circulation plane, with the top and bottom halves of the plane oscillating out of phase by half a cycle. The data for Γ=1 and σ=4.4 showed that the oscillation amplitude varied irregularly in time, yielding a Gaussian probability distribution centred at zero for the displacement angle. This result can be described well by the equation of motion of a stochastically driven damped harmonic oscillator. It suggests that the existence of the oscillations is a consequence of the stochastic driving by the small-scale turbulent background fluctuations of the system, rather than a consequence of a Hopf bifurcation of the deterministic system. The power spectrum of the LSC orientation had a peak at finite frequency with a quality factor Q≃5, nearly independent of R. For samples with Γ≥2 we did not find this mode, but there remained a characteristic periodic signal that was detectable in the area density ρp of the plumes above the bottom-plate centre. Measurements of ρp revealed a strong dependence on the Rayleigh number R, and on the aspect ratio Γ that could be represented by ρp ~ Γ2.7±0.3. Movies are available with the online version of the paper.

Aero-optics of subsonic turbulent boundary layers

2012 ◽  
Vol 696 ◽  
pp. 122-151 ◽  
Author(s):  
Kan Wang ◽  
Meng Wang
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Boundary Layers ◽  
Dominant Role ◽  
Elevation Angle ◽  
Turbulent Boundary Layers ◽  
Density Field ◽  
Small Scale ◽  
Number Range ◽  
Scale Turbulence

AbstractCompressible large-eddy simulations are carried out to study the aero-optical distortions caused by Mach 0.5 flat-plate turbulent boundary layers at Reynolds numbers of ${\mathit{Re}}_{\theta } = 875$, 1770 and 3550, based on momentum thickness. The fluctuations of refractive index are calculated from the density field, and wavefront distortions of an optical beam traversing the boundary layer are computed based on geometric optics. The effects of aperture size, small-scale turbulence, different flow regions and beam elevation angle are examined and the underlying flow physics is analysed. It is found that the level of optical distortion decreases with increasing Reynolds number within the Reynolds-number range considered. The contributions from the viscous sublayer and buffer layer are small, while the wake region plays a dominant role, followed by the logarithmic layer. By low-pass filtering the fluctuating density field, it is shown that small-scale turbulence is optically inactive. Consistent with previous experimental findings, the distortion magnitude is dependent on the propagation direction due to anisotropy of the boundary-layer vortical structures. Density correlations and length scales are analysed to understand the elevation-angle dependence and its relation to turbulence structures. The applicability of Sutton’s linking equation to boundary-layer flows is examined, and excellent agreement between linking equation predictions and directly integrated distortions is obtained when the density length scale is appropriately defined.

Plume formation and resonant bifurcations in porous-media convection

1994 ◽  
Vol 272 ◽  
pp. 67-90 ◽  
Author(s):  
Michael D. Graham ◽  
Paul H. Steen
Keyword(s):  
Boundary Layer ◽  
Porous Media ◽  
Rayleigh Number ◽  
Scaling Laws ◽  
Thermal Plumes ◽  
Scaling Behaviour ◽  
Parametric Resonances ◽  
Classical Boundary ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

The classical boundary-layer scaling laws proposed by Howard for Rayleigh–Bénard convection at high Rayleigh number extend to the analogous case of convection in saturated porous media. We computationally study two-dimensional porous-media convection near the onset of this scaling behaviour. The main result of the paper is the observation and study of instabilities that lead to deviations from the scaling relations.At Rayleigh numbers below the scaling regime, boundary-layer fluctuations born at a Hopf bifurcation strengthen and eventually develop into thermal plumes. The appearance of plumes corresponds to the onset of the boundary-layer scaling behaviour of the oscillation frequency and mean Nusselt number, in agreement with the classical theory. As the Rayleigh number increases further, the flow undergoes instabilities that lead to ‘bubbles’ in parameter space of quasi-periodic flow, and eventually to weakly chaotic flow. The instabilities disturb the plume formation process, effectively leading to a phase modulation of the process and to deviations from the scaling laws. We argue that these instabilities correspond to parametric resonances between the timescale for plume formation and the characteristic convection timescale of the flow.

scholarly journals Linear biglobal analysis of Rayleigh–Bénard instabilities in binary fluids with and without throughflow

2012 ◽  
Vol 713 ◽  
pp. 216-242 ◽  
Author(s):  
Jun Hu ◽  
Daniel Henry ◽  
Xie-Yuan Yin ◽  
Hamda BenHadid
Keyword(s):  
Reynolds Number ◽  
Rayleigh Number ◽  
Aspect Ratio ◽  
Critical Rayleigh Number ◽  
Two Dimensional ◽  
Binary Fluids ◽  
Transverse Modes ◽  
Aspect Ratios ◽  
Rayleigh Benard ◽  
Mode A

AbstractThree-dimensional Rayleigh–Bénard instabilities in binary fluids with Soret effect are studied by linear biglobal stability analysis. The fluid is confined transversally in a duct and a longitudinal throughflow may exist or not. A negative separation factor $\psi = \ensuremath{-} 0. 01$, giving rise to oscillatory transitions, has been considered. The numerical dispersion relation associated with this stability problem is obtained with a two-dimensional Chebyshev collocation method. Symmetry considerations are used in the analysis of the results, which allow the classification of the perturbation modes as ${S}_{l} $ modes (those which keep the left–right symmetry) or ${R}_{x} $ modes (those which keep the symmetry of rotation of $\lrm{\pi} $ about the longitudinal mid-axis). Without throughflow, four dominant pairs of travelling transverse modes with finite wavenumbers $k$ have been found. Each pair corresponds to two symmetry degenerate left and right travelling modes which have the same critical Rayleigh number ${\mathit{Ra}}_{c} $. With the increase of the duct aspect ratio $A$, the critical Rayleigh numbers for these four pairs of modes decrease and closely approach the critical value ${\mathit{Ra}}_{c} = 1743. 894$ obtained in a two-dimensional situation, one of the mode (a ${S}_{l} $ mode called mode A) always remaining the dominant mode. Oscillatory longitudinal instabilities ($k\approx 0$) corresponding to either ${S}_{l} $ or ${R}_{x} $ modes have also been found. Their critical curves, globally decreasing, present oscillatory variations when the duct aspect ratio $A$ is increased, associated with an increasing number of longitudinal rolls. When a throughflow is applied, the symmetry degeneracy of the pairs of travelling transverse modes is broken, giving distinct upstream and downstream modes. For small and moderate aspect ratios $A$, the overall critical Rayleigh number in the small Reynolds number range studied is only determined by the upstream transverse mode A. In contrast, for larger aspect ratios as $A= 7$, different modes are successively dominant as the Reynolds number is increased, involving both upstream and downstream transverse modes A and even the longitudinal mode.

scholarly journals Diffusional and accretional growth of water drops in a rising adiabatic parcel: effects of the turbulent collision kernel

2009 ◽  
Vol 9 (7) ◽  
pp. 2335-2353 ◽  
Author(s):  
W. W. Grabowski ◽  
L.-P. Wang
Keyword(s):  
Energy Dissipation ◽  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
Collision Kernel ◽  
Small Scale ◽  
Large Set ◽  
Initiation Time ◽  
Cloud Droplets ◽  
Dissipation Rates ◽  
Speedup Factor

Abstract. A large set of rising adiabatic parcel simulations is executed to investigate the combined diffusional and accretional growth of cloud droplets in maritime and continental conditions, and to assess the impact of enhanced droplet collisions due to small-scale cloud turbulence. The microphysical model applies the droplet number density function to represent spectral evolution of cloud and rain/drizzle drops, and various numbers of bins in the numerical implementation, ranging from 40 to 320. Simulations are performed applying two traditional gravitational collection kernels and two kernels representing collisions of cloud droplets in the turbulent environment, with turbulent kinetic energy dissipation rates of 100 and 400 cm2 s−3. The overall result is that the rain initiation time significantly depends on the number of bins used, with earlier initiation of rain when the number of bins is low. This is explained as a combination of the increase of the width of activated droplet spectrum and enhanced numerical spreading of the spectrum during diffusional and collisional growth when the number of model bins is low. Simulations applying around 300 bins seem to produce rain at times which no longer depend on the number of bins, but the activation spectra are unrealistically narrow. These results call for an improved representation of droplet activation in numerical models of the type used in this study. Despite the numerical effects that impact the rain initiation time in different simulations, the turbulent speedup factor, the ratio of the rain initiation time for the turbulent collection kernel and the corresponding time for the gravitational kernel, is approximately independent of aerosol characteristics, parcel vertical velocity, and the number of bins used in the numerical model. The turbulent speedup factor is in the range 0.75–0.85 and 0.60–0.75 for the turbulent kinetic energy dissipation rates of 100 and 400 cm2 s−3, respectively.

On the Structure of Turbulence in the Bottom Boundary Layer of the Coastal Ocean

2005 ◽  
Vol 35 (1) ◽  
pp. 72-93 ◽  
Author(s):  
W. A. M. Nimmo Smith ◽  
J. Katz ◽  
T. R. Osborn
Keyword(s):  
Boundary Layer ◽  
Energy Spectrum ◽  
Coastal Ocean ◽  
Bottom Boundary Layer ◽  
Small Scale ◽  
Shear Stresses ◽  
Weak Turbulence ◽  
Conditional Sampling ◽  
Bottom Boundary ◽  
Dissipation Rates

Abstract Six sets of particle image velocimetry (PIV) data from the bottom boundary layer of the coastal ocean are examined. The data represent periods when the mean currents are higher, of the same order, and much weaker than the wave-induced motions. The Reynolds numbers based on the Taylor microscale (Reλ) are 300–440 for the high, 68–83 for the moderate, and 14–37 for the weak mean currents. The moderate–weak turbulence levels are typical of the calm weather conditions at the LEO-15 site because of the low velocities and limited range of length scales. The energy spectra display substantial anisotropy at moderate to high wavenumbers and have large bumps at the transition from the inertial to the dissipation range. These bumps have been observed in previous laboratory and atmospheric studies and have been attributed to a bottleneck effect. Spatial bandpass-filtered vorticity distributions demonstrate that this anisotropy is associated with formation of small-scale, horizontal vortical layers. Methods for estimating the dissipation rates are compared, including direct estimates based on all of the gradients available from 2D data, estimates based on gradients of one velocity component, and those obtained from curve fitting to the energy spectrum. The estimates based on vertical gradients of horizontal velocity are higher and show better agreement with the direct results than do those based on horizontal gradients of vertical velocity. Because of the anisotropy and low turbulence levels, a −5/3 line-fit to the energy spectrum leads to mixed results and is especially inadequate at moderate to weak turbulence levels. The 2D velocity and vorticity distributions reveal that the flow in the boundary layer at moderate speeds consists of periods of “gusts” dominated by large vortical structures separated by periods of more quiescent flows. The frequency of these gusts increases with Reλ, and they disappear when the currents are weak. Conditional sampling of the data based on vorticity magnitude shows that the anisotropy at small scales persists regardless of vorticity and that most of the variability associated with the gusts occurs at the low-wave-number ends of the spectra. The dissipation rates, being associated with small-scale structures, do not vary substantially with vorticity magnitude. In stark contrast, almost all the contributions to the Reynolds shear stresses, estimated using structure functions, are made by the high- and intermediate-vorticity-magnitude events. During low vorticity periods the shear stresses are essentially zero. Thus, in times with weak mean flow but with wave orbital motion, the Reynolds stresses are very low. Conditional sampling based on phase in the wave orbital cycle does not show any significant trends.

Passive Manipulation of Separation-Bubble Transition Using Surface Modifications

2009 ◽  
Vol 131 (2) ◽  
Author(s):  
Brian R. McAuliffe ◽  
Metin I. Yaras
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Shear Layer ◽  
Separation Bubble ◽  
Surface Modifications ◽  
Small Scale ◽  
Two Dimensional ◽  
Vortex Pairing ◽  
Separated Shear Layer ◽  
Scale Turbulence

Through experiments using two-dimensional particle-image velocimetry (PIV), this paper examines the nature of transition in a separation bubble and manipulations of the resultant breakdown to turbulence through passive means of control. An airfoil was used that provides minimal variation in the separation location over a wide operating range, with various two-dimensional modifications made to the surface for the purpose of manipulating the transition process. The study was conducted under low-freestream-turbulence conditions over a flow Reynolds number range of 28,000–101,000 based on airfoil chord. The spatial nature of the measurements has allowed identification of the dominant flow structures associated with transition in the separated shear layer and the manipulations introduced by the surface modifications. The Kelvin–Helmholtz (K-H) instability is identified as the dominant transition mechanism in the separated shear layer, leading to the roll-up of spanwise vorticity and subsequent breakdown into small-scale turbulence. Similarities with planar free-shear layers are noted, including the frequency of maximum amplification rate for the K-H instability and the vortex-pairing phenomenon initiated by a subharmonic instability. In some cases, secondary pairing events are observed and result in a laminar intervortex region consisting of freestream fluid entrained toward the surface due to the strong circulation of the large-scale vortices. Results of the surface-modification study show that different physical mechanisms can be manipulated to affect the separation, transition, and reattachment processes over the airfoil. These manipulations are also shown to affect the boundary-layer losses observed downstream of reattachment, with all surface-indentation configurations providing decreased losses at the three lowest Reynolds numbers and three of the five configurations providing decreased losses at the highest Reynolds number. The primary mechanisms that provide these manipulations include: suppression of the vortex-pairing phenomenon, which reduces both the shear-layer thickness and the levels of small-scale turbulence; the promotion of smaller-scale turbulence, resulting from the disturbances generated upstream of separation, which provides quicker transition and shorter separation bubbles; the elimination of the separation bubble with transition occurring in an attached boundary layer; and physical disturbance, downstream of separation, of the growing instability waves to manipulate the vortical structures and cause quicker reattachment.

The effect of distant side-walls on the evolution and stability of finite-amplitude Rayleigh-Bénard convection

1981 ◽  
Vol 378 (1775) ◽  
pp. 539-566 ◽  
Keyword(s):  
Rayleigh Number ◽  
Prandtl Number ◽  
Finite Amplitude ◽  
Two Dimensional ◽  
Bénard Convection ◽  
Benard Convection ◽  
Side Walls ◽  
Rayleigh Benard Convection ◽  
Rayleigh Numbers ◽  
Rayleigh Benard

A recent study by Cross et al . (1980) has described a class of finite-amplitude phase-winding solutions of the problem of two-dimensional Rayleigh-Bénard convection in a shallow fluid layer of aspect ratio 2 L (≫ 1) confined laterally by rigid side-walls. These solutions arise at Rayleigh numbers R = R 0 + O ( L -1 ) where R 0 is the critical Rayleigh number for the corresponding infinite layer. Nonlinear solutions of constant phase exist for Rayleigh numbers R = R 0 + O ( L -2 ) but of these only the two that bifurcate at the lowest value of R are stable to two-dimensional linearized disturbances in this range (Daniels 1978). In the present paper one set of the class of phase-winding solutions is found to be stable to two-dimensional disturbances. For certain values of the Prandtl number of the fluid and for stress-free horizontal boundaries the results predict that to preserve stability there must be a continual readjustment of the roll pattern as the Rayleigh number is raised, with a corresponding increase in wavelength proportional to R - R 0 . These solutions also exhibit hysteresis as the Rayleigh number is raised and lowered. For other values of the Prandtl number the number of rolls remains unchanged as the Rayleigh number is raised, and the wavelength remains close to its critical value. It is proposed that the complete evolution of the flow pattern from a static state must take place on a number of different time scales of which t = O(( R - R 0 ) -1 ) and t = O(( R - R 0 ) -2 ) are the most significant. When t = O(( R - R 0 ) -1 ) the amplitude of convection rises from zero to its steady-state value, but the final lateral positioning of the rolls is only completed on the much longer time scale t = O(( R - R 0 ) -2 ).

Sign in / Sign up

Export Citation Format

Share Document