On the prediction of the kurtosis of velocity derivatives in a turbulent field

The prediction of the values of non-dimensional fourth-order moment (kurtosis) of the velocity derivative in a turbulent field is made under the assumption that the values of kurtosis depend on both the turbulence Reynolds number and the intermittency factor. The method consists of modeling a suitable probability density of the variable in a given turbulence Reynolds number and the intermittency factor. A crude model of the probability density function is derived, and the numerical calculations based on the model show excellent agreement with many of the experimental data. The analysis shows that the values of kurtosis depend strongly on the intermittency factor, and that depending on the value of the intermittency factor, it is entirely possible to have values of kurtosis as low as five in a flow with a turbulence Reynolds number of 5000.     

Intermittency of turbulence within open canopies

Eddy covariance data have been analyzed to examine intermittency and clustering properties of turbulence within open canopies. Intermittency consists of two aspects: one is related to amplitude variation and the other to clustering. Using the telegraph approximation (TA), the clustering properties have been separated from amplitude effects. Intermittency of canopy turbulence has been explored via clustering exponent, probability density distribution of inter-pulse period of TA, intermittency exponent and structure kurtosis. Intermittency and clustering properties of turbulence within open canopies show similar features to those within dense canopy but some differences are also noted. Unlike within a dense canopy, temperature does not show larger clustering than velocity, which seems to be due to a different thermal structure of the sub-canopy and larger vertical scale of canopy eddy within open canopies. Within the crown region, the inter-pulse probability distribution of TA does not show the ‘double regime’ which was observed within the crown of a dense canopy, indicating less influence of near-field source on canopy turbulence within open canopies. For TA series of the flow variables, intermittency exponent is higher for temperature than for two velocity components within open canopies, which are opposite within a dense canopy. When comparing intermittency for flow variables and their TA series, it is shown that amplitude variation mitigates intermittency for both velocity components and temperature although amplitude variations play a much larger role in velocity intermittency than in temperature counterpart. Kurtosis analysis demonstrates that structure kurtosis is higher at large scales in stable conditions than in unstable conditions, indicating the existence of global intermittency due to stable stratification. The intermittency features of canopy turbulence within open canopies have been discussed in comparison with those within a dense canopy.     

Global Intermittency and Collapsing Turbulence in the Stratified Planetary Boundary Layer

Direct numerical simulation of the turbulent Ekman layer over a smooth wall is used to investigate bulk properties of a planetary boundary layer under stable stratification. Our simplified configuration depends on two non-dimensional parameters: a Richardson number characterizing the stratification and a Reynolds number characterizing the turbulence scale separation. This simplified configuration is sufficient to reproduce global intermittency, a turbulence collapse, and the decoupling of the surface from the outer region of the boundary layer. Global intermittency appears even in the absence of local perturbations at the surface; the only requirement is that large-scale structures several times wider than the boundary-layer height have enough space to develop. Analysis of the mean velocity, turbulence kinetic energy, and external intermittency is used to investigate the large-scale structures and corresponding differences between stably stratified Ekman flow and channel flow. Both configurations show a similar transition to the turbulence collapse, overshoot of turbulence kinetic energy, and spectral properties. Differences in the outer region resulting from the rotation of the system lead, however, to the generation of enstrophy in the non-turbulent patches of the Ekman flow. The coefficient of the stability correction function from Monin–Obukhov similarity theory is estimated as \(\beta \approx 5.7\) in agreement with atmospheric observations, theoretical considerations, and results from stably stratified channel flows. Our results demonstrate the applicability of this set-up to atmospheric problems despite the intermediate Reynolds number achieved in our simulations.     

On The Influence Of The Eulerian Velocity PDF Closure On The Eddy Diffusion Coefficient

The flux-gradient model, often used to describe turbulent dispersion, implicitly defines an eddy diffusion coefficient K that is known to be related to the Eulerian probability density function (pdf) of the turbulent velocity field. In the strict limit of applicability of Fick's law, the relationship between K and the pdf is used to investigate the influence of non-Gaussianity on dispersion in homogeneous turbulence. A bi-Gaussian pdf is used as a closure model that allows for separate studies of skewness and kurtosis variations. The choice of model parameters can have a significant influence on K, especially when the pdf is bimodal. Both arbitrariness of the closure and bimodality are then reduced using the maximum entropy criterion for the selection of the free parameter of the closure scheme, together with the assumption that the model is valid only for those values of the parameters for which a unimodal pdf is possible. The variations of K are found to be sensitive to both skewness and kurtosis showing a more complex behaviour than that found in literature.     

On Reynolds Averaging of Turbulence Time Series

We show that validity of Reynolds averaging for estimating the (ensemble) mean of a turbulence time series requires that the series values be both stationary and uncorrelated. In strict statistical terminology, these two conditions are jointly designated as independent identically distributed (i.i.d.). Moreover, we show that when the series values are correlated, knowledge of the correlation between the values is needed to obtain a reliable estimate of the mean. Last, we contend that a viable averaging algorithm must be Reynolds number (Re) dependent, requiring one version for low Re (Gaussian) turbulence and another for high Re (non-Gaussian) turbulence. Alternatively the median (as opposed to the mean) is recommended as a measure of the central tendency of the turbulence probability density function.     

On the lognormality of the small-scale structure of turbulence

Higher-order moments of turbulent velocity gradients and their behavior with Reynolds number were measured in the nearly isotropic turbulent field generated by a square-mesh grid and in a turbulent boundary layer along a flat plate with zero pressure gradient. Hot-wire anemometry and instrumentation combining analog and digital methods were used to measure moments up to the fourteenth order. Measurements of such high-order moments required that particular attention be given to their validity. Involved herein was the evaluation of such effects as nonlinearity, averaging intervals, and the adequacy of the statistics for the tails of the probability density distributions. The results obtained are compared with those of other investigators for a variety of flow configurations in the laboratory as well as in the atmosphere. The concept of the intermittency of the small-scale structure and the theoretical approach involving lognormality of the probability density distribution of the dissipation rate are evaluated.     

A two-dimensional trajectory-simulation model for non-Gaussian,inhomogeneous turbulence within plant canopies

We report a two-dimensional (alongwind u, vertical w) trajectory-simulation model, consistent with Thomson's (1987) well-mixed criteria, that allows for the non-Gaussian turbulence typical of flow within a plant canopy. The effect of non-Gaussian turbulence was examined by formulating a non-Gaussian u, w joint probability density function (PDF) as the sum of two Gaussian joint-PDFs. The resultant PDF reproduced the desired means, variances, skewnesses, and kurtoses, and the correct covariance. In prediction of the location of maximum concentration downwind of a line source in homogeneous, slightly non-Gaussian turbulence, it proved advantageous to incorporate skewness and kurtosis. However, in the case of inhomogeneous, highly non-Gaussian turbulence, the addition of skewness and kurtosis in the model resulted in substantially worse agreement with measurements than the results of the model using Gaussian PDFs. This may be due to inaccuracy in our PDF formulation. Dispersion predictions from the model with Gaussian PDFs were generally not statistically different from measurements. These results indicate that a two-dimensional Gaussian trajectory-simulation approach is adequate to predict mean concentrations and fluxes resulting from sources within plant canopies.     

Multifractal Characteristics of Intermittent Turbulence in the Urban Canopy Layer

The multifractality of energy and thermal dissipation of fully developed intermittent turbulence is investigated in the urban canopy layer under unstable conditions by the singularity spectrum for the fractal dimensions of sets of singularities characterizing multifractals. In order to obtain high-order moment properties of smallscale turbulent dissipation in the inertial range, an ultrasonic anemometer with a high sampling frequency of 100 Hz was used. The authors found that the turbulent signal could be singular everywhere. Moreover, the singular exponents of energy and thermal dissipation rates are most frequently encountered at around 0.2, which is significantly smaller than the singular exponents for a wind tunnel at a moderate Reynolds number. The evidence indicates a higher intermittency of turbulence in the urban canopy layer at a high Reynolds number, which is demonstrated by the data with high temporal resolution. Furthermore, the temperature field is more intermittent than the velocity field. In addition, a large amount of samples could be used for verification of the results.     

Surface-layer intermittency investigated with conditional sampling

A conditional sampling technique is used to provide statistics of surface-layer plume properties. A selection criterion based on the high-frequency variance of the horizontal wind component enables an accurate division of plume and nonplume states. The intermittency factor derived with this technique closely matches values obtained using other techniques at various heights in the atmospheric boundary layer. The intermittency factor in addition to other plume statistics are found to be stability dependent. Conditional averages are used to produce scatter diagrams from which the interrelationships between properties of both the plume and nonplume states can be examined. Several provocative relationships discovered in this way are discussed.An extensive investigation into the bimodal nature of the fine structure of turbulence is described. These results provide the most compelling support for the division of surface-layer turbulence into separate states. Length scales derived from the second moments of distributions fitted to conditionally sampled data are found to correlate with external parameters of the flow.Department of Atmospheric Sciences contribution number 514.     

Boundedness of Turbulent Temperature Probability Distributions,and their Relation to the Vertical Profile in the Convective Boundary Layer

Higher-order moments, minima and maxima of turbulent temperature and water vapour mixing ratio probability density functions measured with an eddy-covariance system near the ground were related to each other and to vertical boundary-layer profiles of the same scalars obtained through airborne soundings. The dependence of kurtosis on squared skewness showed a kurtosis intercept below the Gaussian expectation, suggesting a compression of the probability density function by the presence of natural boundaries. This hypothesis was corroborated by comparing actual minima and maxima of turbulent fluctuations to estimates obtained from the first four sample moments by fitting a four-parameter beta distribution. The most sharply defined boundaries were found for the minima of temperature datasets during the day, indicating that negative temperature fluctuations at the sensor are limited by the availability of lower temperatures in the boundary layer. By comparison to vertical profiles, it could be verified that the turbulent minimum of temperature near the ground is close to the minimum of potential temperature in the boundary layer. The turbulent minimum of water vapour mixing ratio was found to be equal to the mixing ratio at a height above the minimum of the temperature profile. This height roughly agrees with the top of the non-local unstable domain according to bulk Richardson number profiles. We conclude that turbulence statistics measured near the surface cannot be solely explained by local effects, but contain information about the whole boundary layer including the entrainment zone.     

A model for sensible heat flux probability density function for near-neutral and slightly-stable atmospheric flows

The probability density function for sensible heat flux was measured above a uniform dry lakebed (Owens lake) in Owens Valley, California. It was found that for moderately stable to near neutral atmospheric stability conditions, the probability density function exhibits well defined exponential tails. These exponential tails are consistent with many laboratory boundarylayer measurements and numerical simulations. A model for the sensible heat flux probability density function was developed and tested. A key assumption in the model derivation was the near Gaussian statistics of the vertical velocity and temperature fluctuations. This assumption was verified from time series measurements of temperature and vertical velocity. The parameters for the sensible heat flux probability density function model were also derived from mean meteorological and surface conditions using surface-layer similarity theory. It was found that the best agreement between modeled and measured sensible heat flux probability density function was at the tails. Finally, a relation between the intermittency parameter, the probability density function, and the mean meteorological conditions was derived. This relation rigorously links the intermittency parameter to mean meteorological conditions.     

大气边界层湍流标量场的概率分布及其特征分析

利用2004年11月在白洋淀地区和2005年1月在中国科学院大气物理研究所北京325 m气象塔的47 m高度由超声风温仪和水汽二氧化碳分析仪观测的湍流脉动资料,分析了大气边界层不同下垫面湍流标量场(温度、水汽和二氧化碳)的概率分布及其特征。标量场的概率分布通常不同于高斯分布,而且还会产生偏斜,可以用指数分布描述。因此,标量的偏斜度通常不为0,陡峭度也往往比3大。非0偏斜度的出现可能是由湍流时间序列中的相干结构和间歇性部分造成的。其中,相干结构的存在使概率分布偏斜,但是,它们对偏斜度的贡献相对较小,而与概率分布的长尾现象有关的间歇性则会使偏斜度大大增加。温度、水汽和二氧化碳的平均偏斜度和陡峭度反应了标量场与稳定度、下垫面、天气条件、源汇等因素之间的关系。在不同下垫面,温度和感热通量的偏斜度随稳定度变化比较一致;水汽通量的偏斜度在稳定和不稳定条件下都为正,而水汽本身在不稳定条件下可能出现负的偏斜度;二氧化碳和二氧化碳通量的偏斜度受下垫面影响很大,在不同下垫面,偏斜度与稳定度之间的关系并不一致;而3个标量的陡峭度随稳定度的变化不显著,它们与相应的偏斜度之间存在平方关系。     

Statistical characteristics of concentration fluctuations in dispersing plumes in the atmospheric surface layer

Measurements have been made of concentration fluctuations in a dispersing plume from an elevated point source in the atmospheric surface layer using a recently developed fast-response photoionization detector. This detector, which has a frequency response (–6 dB point) of about 100 Hz, is shown to be capable of resolving the fluctuation variance contributed by the energetic subrange and most of the inertial-convective subrange, with a reduction in the fluctuation variance due to instrument smoothing of the finest scales present in the plume of at most 4%.Concentration time series have been analyzed to obtain the statistical characteristics of both the amplitude and temporal structure of the dispersing plume. We present alongwind and crosswind concentration fluctuation profiles of statistics of amplitude structure such as total and conditional fluctuation intensity, skewness and kurtosis, and of temporal structure such as intermittency factor, burst frequency, and mean burst persistence time. Comparisons of empirical concentration probability distributions with a number of model distributions show that our near-neutral data are best represented by the lognormal distribution at shorter ranges, where both plume meandering and fine-scale in-plume mixing are equally important (turbulent-convective regime), and by the gamma distribution at longer ranges, where internal structure or spottiness is becoming dominant (turbulent-diffusive regime). The gamma distribution provides the best model of the concentration pdf over all downwind fetches for data measured under stable stratification. A physical model is developed to explain the mechanism-induced probabilistic schemes in the alongwind development of a dispersing plume, that lead to the observed probability distributions of concentration. Probability distributions of concentration burst length and burst return period have been extracted and are shown to be modelled well with a powerlaw distribution. Power spectra of concentration fluctuations are presented. These spectra exhibit a significant inertial-convective subrange, with the frequency at the spectral peak decreasing with increasing downwind fetch. The Kolmogorov constant for the inertial-convective subrange has been determined from the measured spectra to be 0.17±0.03.     

The Effects of Thermal Stratification on Clustering Properties of Canopy Turbulence

The intermittent structure of turbulence within the canopy sublayer (CSL) is sensitive to the presence of foliage and to the atmospheric stability regime. How much of this intermittency originates from amplitude variability or clustering properties remains a vexing research problem for CSL flows. Using a five-level set of measurements collected within a dense hardwood canopy, the clustering properties of CSL turbulence and their dependence on atmospheric stability are explored using the telegraphic approximation (TA). The binary structure of the TA removes any amplitude variability from turbulent excursions but retains their zero-crossing behaviour, and thereby isolating the role of clustering in intermittency. A relationship between the spectral exponents of the actual and the TA series is derived across a wide range of atmospheric stability regimes and for several flow variables. This relationship is shown to be consistent with a relationship derived for long-memory and monofractal processes such as fractional Brownian motion (fBm). Moreover, it is demonstrated that for the longitudinal and vertical velocity components, the vegetation does not appreciably alter fine-scale clustering but atmospheric stability does. Stable atmospheric stability conditions is characterized by more fine scale clustering when compared to other atmospheric stability regimes. For scalars, fine-scale clustering above the canopy is similar to its velocity counterpart but is significantly increased inside the canopy, especially under stable stratification. Using simplified scaling analysis, it is demonstrated that clustering is much more connected to space than to time within the CSL. When comparing intermittency for flow variables and their TA series, it is shown that for velocity, amplitude variations modulate intermittency for all stability regimes. However, amplitude variations play only a minor role in scalar intermittency. Within the crown region of the canopy, a ‘double regime’ emerges in the inter-pulse duration probability distributions not observed in classical turbulence studies away from boundaries. The double regime is characterized by a power-law distribution for shorter inter-pulse periods and a log-normal distribution for large inter-pulse periods. The co-existence of these two regimes is shown to be consistent with near-field/far-field scaling arguments. In the near-field, short inter-pulse periods are controlled by the source strength, while in the far-field long inter-pulse periods are less affected by the precise source strength details and more affected by the transport properties of the background turbulence.     

Some similarity states of stably stratified homogeneous turbulence

The decay of statistically homogeneous velocity and density fluctuations in a stably stratified fluid is considered. Over decay times long compared with the turbulence time scale but short compared with the period of internal gravity waves, three distinct high Reynolds number similarity states may develop. These. similarity states are a consequence of the invariance of the low wavenumber coefficients of the three-dimensional kinetic or potential energy spectrum, and their preferential development depends on the relative magnitudes of the initial kinetic and potential energy per unit mass of the fluid. When the turbulence has decayed over a time comparable with the period of the gravity waves, the three similarity states mentioned above are disrupted. Evidence will be presented of a new similarity state which then develops asymptotically. In this similarity state, the time decay exponent of the total energy per unit mass of the turbulence is reduced by a factor of two from its value for decaying isotropic turbulence, and the associated vertical integral scale approaches a constant independent of time.     

Intermittency of Turbulence in the Atmospheric Boundary Layer: Scaling Exponents and Stratification Influence

We show the relationship between the intermittency of turbulence and the type of stratification for different atmospheric situations during the SABLES98 field campaign. With this objective, we first demonstrate the scaling behaviour of the velocity structure functions corresponding to these situations; next, we analyze the curvature of the scaling exponents of the velocity structure functions versus the order of these functions (ζ _p vs. p), where ζ _p are the exponents of the power relation for the velocity structure function with respect to the scale. It can be proved that this curve must be concave, under the assumption that the incompressible approximation does not break down at high Reynolds numbers. The physical significance of this kind of curvature is that the energy dissipation rate increases as the scale of the turbulent eddies diminishes (intermittency in the usual sense). However, the constraints imposed by stability, preventing full development of the turbulence, allow the function ζ _p versus p to show any type of curvature. In this case, waves of high frequency trapped by the stability, or bursts of turbulence caused by the breaking up of internal waves, may produce a redistribution of energy throughout the scaling range. Due to this redistribution, the variation with the scale of the energy dissipation rate may be smaller (decreasing the intermittency) and, even in more stable situations, this rate may diminish (instead of increasing) as the scale diminishes (convex form of the curve ζ _p vs. p).     

The effects of source size on concentration fluctuations in plumes

The effects of source size on plume behaviour have been examined in a 1.2 m wind tunnel boundary layer for isokinetic sources with diameters from 3 to 35 mm at source heights of 230 mm and at ground level. Experimental measurements of mean concentration and the variance, intermittency and probability density functions of the concentration fluctuations were obtained. In addition, a fluctuating Gaussian plume model is presented which reproduces many of the observed features of the elevated emission. The mean plume width becomes independent of source size much more rapidly than the instantaneous plume width. Since it is the meandering of the instantaneous plume which generates most of the concentration fluctuations near the source, these are also dependent on source size. The flux of variance in the plume reaches a maximum, whose value is greatest for the smallest source size, close to the source and thereafter is monotonically decreasing. The intermittency factor reaches a minimum, whose value is lowest for the smallest source, and increases back towards one. Concentration fluctuations for the ground-level source are much less dependent on source size due to the effects of the surface.     

High-Order Statistics of Temperature Fluctuations in an Unstable Atmospheric Surface Layer over Grassland

Skewness(S) and kurtosis(K) of temperature in the surface layer over a grassland are investigated under unstable thermal stratifications. We find that both skewness and kurtosis generally obey Monin–Obukhov similarity theory and tend to be constant values(1.5 and 5.3, respectively) when the stability parameter z/L -2. Quantitative formulas of the similarity functions are proposed. The temperature probability density function(PDF) is close to Gaussian in near neutral stratification and non-Gaussian in unstable stratification. The influence of coherent motions on the PDF behavior is analyzed using the quadrant analysis technique. It shows that PDF behaviors are controlled by ejections and sweeps. The results also indicate that the PDF type of the ejections always follows a Gaussian distribution, while the PDF of the sweeps changes with stability.     

A Statistical Model for the Prediction of Wind-Speed Probabilities in the Atmospheric Surface Layer

Wind fields in the atmospheric surface layer (ASL) are highly three-dimensional and characterized by strong spatial and temporal variability. For various applications such as wind-comfort assessments and structural design, an understanding of potentially hazardous wind extremes is important. Statistical models are designed to facilitate conclusions about the occurrence probability of wind speeds based on the knowledge of low-order flow statistics. Being particularly interested in the upper tail regions we show that the statistical behaviour of near-surface wind speeds is adequately represented by the Beta distribution. By using the properties of the Beta probability density function in combination with a model for estimating extreme values based on readily available turbulence statistics, it is demonstrated that this novel modelling approach reliably predicts the upper margins of encountered wind speeds. The model’s basic parameter is derived from three substantially different calibrating datasets of flow in the ASL originating from boundary-layer wind-tunnel measurements and direct numerical simulation. Evaluating the model based on independent field observations of near-surface wind speeds shows a high level of agreement between the statistically modelled horizontal wind speeds and measurements. The results show that, based on knowledge of only a few simple flow statistics (mean wind speed, wind-speed fluctuations and integral time scales), the occurrence probability of velocity magnitudes at arbitrary flow locations in the ASL can be estimated with a high degree of confidence.     

Turbulence structure in the planetary boundary layer

This review of the last three years of progress in the understanding of wind profiles and the structure of turbulence in the planetary boundary layer is divided into three parts. The first part, by N. E. Busch, deals with the atmospheric surface layer below 30 m. It is shown that the Monin-Oboukhov similarity hypotheses fail at low frequencies and large wave-lengths, probably due to mesoscale influences. Also, it is suggested that the neutral surface layer is a poor reference state in some respects, because the structure of turbulence in unstable conditions is quite different from that in stable stratification. The second part, by H. Tennekes, is concerned with the intermittency of the dissipative structure of turbulence and its effects on the velocity and temperature structure functions. It is shown that the modified Kolmogorov-Oboukhov theory, which attempts to explain the consequences of the dissipative intermittency, is unable to predict the shape of the temperature structure functions. The third part of this review, by H. A. Panofsky, deals with wind profiles and turbulence structure above 30 m. It is shown that between 30 and 150 m, surface-layer formulas can be used, if such mesoscale effects as changes of terrain roughness are taken into account where needed. Experimental data on turbulence above 150 m are quite sparse; some of the current scaling laws that can be used in this region are described.