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

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.

2010 ◽  
Vol 27 (10) ◽  
pp. 1652-1664 ◽  
Author(s):  
Ewan J. O’Connor ◽  
Anthony J. Illingworth ◽  
Ian M. Brooks ◽  
Christopher D. Westbrook ◽  
Robin J. Hogan ◽  
...  

Abstract A method of estimating dissipation rates from a vertically pointing Doppler lidar with high temporal and spatial resolution has been evaluated by comparison with independent measurements derived from a balloon-borne sonic anemometer. This method utilizes the variance of the mean Doppler velocity from a number of sequential samples and requires an estimate of the horizontal wind speed. The noise contribution to the variance can be estimated from the observed signal-to-noise ratio and removed where appropriate. The relative size of the noise variance to the observed variance provides a measure of the confidence in the retrieval. Comparison with in situ dissipation rates derived from the balloon-borne sonic anemometer reveal that this particular Doppler lidar is capable of retrieving dissipation rates over a range of at least three orders of magnitude. This method is most suitable for retrieval of dissipation rates within the convective well-mixed boundary layer where the scales of motion that the Doppler lidar probes remain well within the inertial subrange. Caution must be applied when estimating dissipation rates in more quiescent conditions. For the particular Doppler lidar described here, the selection of suitably short integration times will permit this method to be applicable in such situations but at the expense of accuracy in the Doppler velocity estimates. The two case studies presented here suggest that, with profiles every 4 s, reliable estimates of ε can be derived to within at least an order of magnitude throughout almost all of the lowest 2 km and, in the convective boundary layer, to within 50%. Increasing the integration time for individual profiles to 30 s can improve the accuracy substantially but potentially confines retrievals to within the convective boundary layer. Therefore, optimization of certain instrument parameters may be required for specific implementations.


2017 ◽  
Vol 814 ◽  
pp. 165-184 ◽  
Author(s):  
Yang Zhang ◽  
Quan Zhou ◽  
Chao Sun

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.


2006 ◽  
Vol 63 (4) ◽  
pp. 1291-1307 ◽  
Author(s):  
B. B. Balsley ◽  
R. G. Frehlich ◽  
M. L. Jensen ◽  
Y. Meillier

Abstract Some 50 separate high-resolution profiles of small-scale turbulence defined by the energy dissipation rate (ɛ), horizontal wind speed, and temperature from near the surface, through the nighttime stable boundary layer (SBL), and well into the residual layer are used to compare the various definitions of SBL height during nighttime stable conditions. These profiles were obtained during postmidnight periods on three separate nights using the Tethered Lifting System (TLS) during the Cooperative Atmosphere–Surface Exchange Study (CASES-99) campaign in east-central Kansas, October 1999. Although the number of profiles is insufficient to make any definitive conclusions, the results suggest that, under most conditions, the boundary layer top can be reasonably estimated in terms of a very significant decrease in the energy dissipation rate (i.e., the mixing height) with height. In the majority of instances this height lies slightly below the height of a pronounced minimum in wind shear and slightly above a maximum in N 2, where N is the Brunt–Väisälä frequency. When combined with flux measurements and vertical velocity variance data obtained from the nearby 55-m tower, the results provide additional insights into SBL processes, even when the boundary layer, by any definition, extends to heights well above the top of the tower. Both the TLS profiles and tower data are then used for preliminary high-resolution studies into various categories of SBL structure, including the so-called upside-down boundary layer.


2018 ◽  
Vol 856 ◽  
pp. 135-168 ◽  
Author(s):  
S. T. Salesky ◽  
W. Anderson

A number of recent studies have demonstrated the existence of so-called large- and very-large-scale motions (LSM, VLSM) that occur in the logarithmic region of inertia-dominated wall-bounded turbulent flows. These regions exhibit significant streamwise coherence, and have been shown to modulate the amplitude and frequency of small-scale inner-layer fluctuations in smooth-wall turbulent boundary layers. In contrast, the extent to which analogous modulation occurs in inertia-dominated flows subjected to convective thermal stratification (low Richardson number) and Coriolis forcing (low Rossby number), has not been considered. And yet, these parameter values encompass a wide range of important environmental flows. In this article, we present evidence of amplitude modulation (AM) phenomena in the unstably stratified (i.e. convective) atmospheric boundary layer, and link changes in AM to changes in the topology of coherent structures with increasing instability. We perform a suite of large eddy simulations spanning weakly ($-z_{i}/L=3.1$) to highly convective ($-z_{i}/L=1082$) conditions (where$-z_{i}/L$is the bulk stability parameter formed from the boundary-layer depth$z_{i}$and the Obukhov length $L$) to investigate how AM is affected by buoyancy. Results demonstrate that as unstable stratification increases, the inclination angle of surface layer structures (as determined from the two-point correlation of streamwise velocity) increases from$\unicode[STIX]{x1D6FE}\approx 15^{\circ }$for weakly convective conditions to nearly vertical for highly convective conditions. As$-z_{i}/L$increases, LSMs in the streamwise velocity field transition from long, linear updrafts (or horizontal convective rolls) to open cellular patterns, analogous to turbulent Rayleigh–Bénard convection. These changes in the instantaneous velocity field are accompanied by a shift in the outer peak in the streamwise and vertical velocity spectra to smaller dimensionless wavelengths until the energy is concentrated at a single peak. The decoupling procedure proposed by Mathiset al.(J. Fluid Mech., vol. 628, 2009a, pp. 311–337) is used to investigate the extent to which amplitude modulation of small-scale turbulence occurs due to large-scale streamwise and vertical velocity fluctuations. As the spatial attributes of flow structures change from streamwise to vertically dominated, modulation by the large-scale streamwise velocity decreases monotonically. However, the modulating influence of the large-scale vertical velocity remains significant across the stability range considered. We report, finally, that amplitude modulation correlations are insensitive to the computational mesh resolution for flows forced by shear, buoyancy and Coriolis accelerations.


2010 ◽  
Vol 40 (11) ◽  
pp. 2401-2417 ◽  
Author(s):  
Pascale Bouruet-Aubertot ◽  
Hans van Haren ◽  
M. Pascale Lelong

Abstract Deep-ocean high-resolution moored temperature data are analyzed with a focus on superbuoyant frequencies. A local Taylor hypothesis based on the horizontal velocity averaged over 2 h is used to infer horizontal wavenumber spectra of temperature variance. The inertial subrange extends over fairly low horizontal wavenumbers, typically within 2 × 10−3 and 2 × 10−1 cycles per minute (cpm). It is therefore interpreted as a stratified inertial subrange for most of this wavenumber interval, whereas in some cases the convective inertial subrange is resolved as well. Kinetic energy dissipation rate ε is inferred using theoretical expressions for the stratified inertial subrange. A wide range of values within 10−9 and 4 × 10−7 m2 s−3 is obtained for time periods either dominated by semidiurnal tides or by significant subinertial variability. A scaling for ε that depends on the potential energy within the inertio-gravity waves (IGW) frequency band PEIGW and the buoyancy frequency N is proposed for these two cases. When semidiurnal tides dominate, ε ≃ (PEIGWN)3/2, whereas ε ≃ PEIGWN in the presence of significant subinertial variability. This result is obtained for energy levels ranging from 1 to 30 times the Garrett–Munk energy level and is in contrast with classical finescale parameterization in which ε ∼ (PEIGW)2 that applies far from energy sources. The specificities of the stratified bottom boundary layer, namely a weak stratification, may account for this difference.


2009 ◽  
Vol 9 (7) ◽  
pp. 2335-2353 ◽  
Author(s):  
W. W. Grabowski ◽  
L.-P. Wang

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.


2005 ◽  
Vol 35 (1) ◽  
pp. 72-93 ◽  
Author(s):  
W. A. M. Nimmo Smith ◽  
J. Katz ◽  
T. R. Osborn

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.


1970 ◽  
Vol 41 (1) ◽  
pp. 141-152 ◽  
Author(s):  
R. W. Stewart ◽  
J. R. Wilson ◽  
R. W. Burling

Derivatives of velocity signals obtained in a turbulent boundary layer are examined for correspondence to the lognormal distribution. It is found that there is rough agreement but that unlikely events at high values are much less common in the observed fields than would be inferred from the lognormal distribution. The actual distributions correspond more to those obtained from a random walk with a limited number of steps, so the difference between these distributions and the lognormal may be related to the fact that the Reynolds number is finite.The third-order structure function is examined, and found to be roughly consistent with the existence of an inertial subrange of a Kolmogoroff equilibrium reacute;gime over a range of scale which is a priori reasonable but which is far less extensive than the $-\frac{5}{3}$ region of the spectrum.


1970 ◽  
Vol 44 (1) ◽  
pp. 145-159 ◽  
Author(s):  
C. W. Van Atta ◽  
W. Y. Chen

Structure functions of turbulent velocity fluctuations up to fourth order have been measured at several heights in the atmospheric boundary layer over the open ocean, and the results are compared with theoretical predictions for separations in the inertial subrange. The behaviour of second- and third-order quantities shows substantial agreement with the predictions of Kolmogorov's original theory over a wide range of separations, but the results of a recent modification of the theory, attempting to account for intermittency in the local dissipation rate, are also consistent with the data over somewhat shorter separation intervals. The behaviour of the measured fourth-order structure function disagrees with that predicted from Kolmogorov's original work, but good agreement is found with the results of the modified theory.


2008 ◽  
Vol 8 (4) ◽  
pp. 14717-14763 ◽  
Author(s):  
W. W. Grabowski ◽  
L.-P. Wang

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.


Sign in / Sign up

Export Citation Format

Share Document