Predicting plume meandering and averaging time effects on mean and fluctuating concentrations in atmospheric dispersion simulated in a water channel

Plume meandering and averaging time effects were measured directly using a high spatial resolution, high frequency, linescan laser-induced fluorescence (LIF) technique for measuring scalar concentrations in a plume dispersing in a water channel. Post-processing of the collected data removed time dependent background dye levels and corrected for attenuation across the laser beam to produce accurate measurements over long sample times in both a rough surface boundary-layer shear flow and shear free grid-generated turbulent flow. The data were used to verify the applicability of a meandering plume model for predicting the properties of mean and fluctuating concentrations. The centroid position of the crosswind concentration profile was found to have a Gaussian probability density function and the instantaneous plume spread about the centroid fluctuated log-normally. A modified travel-time power law model for averaging time adjustment was developed and compared to the widely used, but much less accurate, 0.2 power-law model.     

Lidar measurements of plume statistics

Surface-layer aerosol diffusion experiments have been conducted using artificial smoke plume releases at ground level over flat and homogeneously vegetated terrain at the Meppen proving grounds in the Federal Republic of Germany (1989). At fixed downwind locations in the range out to 800 m from the source, instantaneous crosswind plume profiles were detected repetitively at high spatial (1.5 m) and temporal (3 sec) intervals by use of a mini LIDAR system. The experiments were accompanied by measurement of the surface-layer mean wind and turbulence quantities by sonic anemometers. On the basis of measured crosswind concentration profiles, the following statistics were obtained: 1) Mean profile, 2) Root mean square profile, 3) Fluctuation intensities, and 4) Intermittency factors. Furthermore, some experimentally determined probability density functions (pdf's) of the fluctuations are presented. All the measured statistics are referred to a fixed and a moving frame of reference, the latter being defined as a frame of reference from which the (low frequency) plume meander is removed. Finally, the measured statistics are compared with statistics on concentration fluctuations obtained with a simple puff diffusion model (RIMPUFF) developed at Risø.     

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.     

Analysis of wind tunnel dispersion data over two-dimensional obstacles

An analysis of wind tunnel data of dispersion from elevated point sources over a flat floor and gently sloping, 2-D hills is performed. (The data were obtained by Khurshudyan et al., 1981.) Formulas for the mixing lengths and the Lagrangian time scales are tested, suitable for use in various dispersion models. Some expressions for the vertical first moments of the plume concentration distribution suitable for shear flow (Hunt, 1985) are also tested.Then, a normalization is suggested, based on the source mean flow and turbulence parameters, for the ground-level concentration maximum value and position along the plume centerline. Using this normalization, the maximum position is almost constant, regardless of source height variation and of whether the hill is present or not, at least for source positions upstream or over the top.The maximum values allow the determination of normalized terrain amplification factors TAF_N, which are shown to be in most cases closer to one than the corresponding TAF obtained without normalization.     

Characteristics of canopy turbulence over a deciduous forest on complex terrain

Canopy turbulence plays an important role in mass and energy exchanges at the canopy-atmosphere interface. Despite extensive studies on canopy turbulence over a flat terrain, less attention has been given to canopy turbulence in a complex terrain. The purpose of this study is to scrutinize characteristics of canopy turbulence in roughness sublayer over a hilly forest terrain. We investigated basic turbulence statistics, conditionally sampled statistics, and turbulence spectrum in terms of different atmospheric stabilities, wind direction and vertical structures of momentum fluxes. Similarly to canopy turbulence over a homogeneous terrain, turbulence statistics showed coherent structure. Both quadrant and spectrum analysis corroborated the role of intermittent and energetic eddies with length scale of the order of canopy height, regardless of wind direction except for shift of peak in vertical wind spectrum to relatively high frequency in the down-valley wind. However, the magnitude of the momentum correlation coefficient in a neutral condition was smaller than typical value over a flat terrain. Further scrutiny manifested that, in the up-valley flow, temperature skewness was larger and the contribution of ejection to both momentum and heat fluxes was larger compared to the downvalley flow, indicating that thermal instability and weaker wind shear in up-valley flow asymmetrically affect turbulent transport within the canopy.     

Flow over the summit of an isolated hill

Observations of the mean flow and turbulence statistics over the summit of an isolated, roughly circular hill, Nyland hill, are presented, Nyland hill rises 70 m above the surrounding terrain and has a base diameter of about 500 m. The summit of the hill is very smooth and allows representative measurements to be made close to the surface. The flow speed 8 m above the summit is increased by a factor of 2 over the upstream speed 8 m above level terrain, and flow separation occurs in the lee of the hill. The mean velocity profile over the summit shows an increase in velocity with height up to about 2 m and then a near constant velocity between 2 and 16 m. The flow perturbation relative to the upstream profile is thus a maximum at about 2 m. The measurements of turbulence structure show how the influence of the hill depends on the length scale of the turbulent eddies involved. Scales greater than the scale of the hill are modified through the flow speed-up whilst scales shorter than the hill suffer complex changes. The short-scale turbulence over the summit is only in local equilibrium in the lowest fraction of a metre. Above this equilibrium region, there is a complex adjustment towards the rapid distortion dynamics which appear to dominate at heights above about 8 m. The detailed results are compared with previous studies and available theories.     

Effect of Vertical Wind Shear on Concentration Fluctuation Statistics in a Point Source Plume

Measurements of concentration fluctuation intensity, intermittency factor, and integral time scale were made in a water channel for a plume dispersing in a well-developed, rough surface, neutrally stable, boundary layer, and in grid-generated turbulence with no mean velocity shear. The water-channel simulations apply to full-scale atmospheric plumes with very short averaging times, on the order of 1–4 min, because plume meandering was suppressed by the water-channel side walls. High spatial and temporal resolution vertical and crosswind profiles of fluctuations in the plume were obtained using a linescan camera laser-induced dye tracer fluorescence technique. A semi-empirical algebraic mean velocity shear history model was developed to predict these concentration statistics. This shear history concentration fluctuation model requires only a minimal set of parameters to be known: atmospheric stability, surface roughness, vertical velocity profile, and vertical and crosswind plume spreads. The universal shear history parameter used was the mean velocity shear normalized by surface friction velocity, plume travel time, and local mean wind speed. The reference height at which this non-dimensional shear history was calculated was important, because both the source and the receptor positions influence the history of particles passing through the receptor position.     

The vertical structure of concentration fluctuation statistics in plumes dispersing in the atmospheric surface layer

This paper describes a study of the vertical structure of concentration fluctuations in a neutrally buoyant plume from an elevated point source in slightly convective to moderately stable meteorological conditions at ranges of between 12.5 and 100 m for a range of source heights between 1 and 5 m. Observations were made of concentration fluctuations in a dispersing plume using a vertical array of sixteen very fast-response photoionization detectors placed at heights between 0.5 and 16 m. Vertical profiles of a number of concentration statistics were extracted, namely, mean concentration, fluctuation intensity, intermittency factor, peak-to-mean concentration ratio, mean dissipation rate of concentration variance, and various concentration time and length scales of dominant motions in the plume (e.g., integral macro-scale, in-plume mid-scale and Taylor micro-scale). The profiles revealed a similarity to corresponding crosswind profiles for a fully elevated plume, but showed greater and greater departure from the latter shapes once the plume had grown in the vertical so that its lower dege began to interact progressively more strongly with the ground. The evolution of the concentration probability density function at a fixed range, but with decreasing height from the ground, is similar to that obtained at a fixed height but with increasing distance from the source. Concentration power spectra obtained at different heights all had an extensive inertial-convective subrange spanning at least two decades in frequency, but spectra measured near the ground had a greater proportion of the total concentration variance in the lower frequencies (energetic subrange), with a correspondingly smaller proportion in the higher frequencies (inertial-convective subrange). It is believed that these effects result from the increased mean shear near the surface, and blocking by the surface. The effect of enhanced shear-induced molecular diffusion on concentration fluctuations is examined.     

A wind tunnel study of turbulent dispersion over two- and three-dimensional gentle hills from upwind point sources in neutral flow

A study of turbulent dispersion over hills for upstream, elevated sources was conducted, based on wind tunnel tracer gas (CO₂) experiments over a gentle 2-D ridge and a 3-D circular hill, both having a cosine-square cross-section. The concentration measurements were made with four different source locations for each hill case (2-D or 3-D), and the study focused on dispersion parameters under the influence of the presence of the hills in order to provide a better understanding of the mechanisms involved.The wind tunnel measurements show that, in the case of gentle hills, the topographic impact on turbulent dispersion from upstream sources is only moderate and is more pronounced for the 3-D than for the 2-D hill. The perturbation in mean flow introduced by the hills, including streamline divergence/convergence, is shown to dominate the changes in the dispersion due to the hills in this case. The plume spread, both in the lateral and the vertical, is enhanced over the upwind hill foot and reduced over the hill top in response to the mean flow slow-down and speed-up at these places, and is further enhanced or reduced due to streamline divergence/convergence in the vertical over the hills as well as in the horizontal over the 3-D hill. These results are also compared with cases of turbulent dispersion over more steep hills (Snyder and Britter, 1987).     

The effects of mesoscale turbulence on plume rise in a land/lake environment

An extensive data base of LIDAR and photographically measured plumes from the Nanticoke thermal generating station (located on the north shore of Lake Erie) has been accumulated over several years of full-scale atmospheric experiments. Analysis of the data shows that the empirical entrainment constant used in the Briggs simple plume rise formulation has a value of 0.8, which is significantly larger than 0.6 found by Briggs (1975). A hypothesis is proposed that mesoscale eddies present in the land/lake environment are responsible for vertical plume meander leading to the observed enhanced vertical plume spread and lower mean plume rise at this site. This is substantiated through theoretical considerations and experimental data. The results are significant in the computation of ground-level concentrations at such complex terrain sites.     

Comparison Between Large-Eddy Simulation and Reynolds-Averaged Navier–Stokes Computations for the MUST Field Experiment. Part II: Effects of Incident Wind Angle Deviation on the Mean Flow and Plume Dispersion

Large-eddy simulations (LES) and Reynolds-averaged Navier–Stokes (RANS) computations of pollutant dispersion are reported for the Mock Urban Setting Test (MUST) field experiment flow. In particular we address the effects of incident wind angle deviation on the mean velocity and on the mean concentration fields. Both computational fluid dynamical methods are assessed by comparing the simulation results with experimental field data. The comparative analysis proposes to relate the plume deflection with the flow channelling effects. The results show that the plume deflection angle varies with the altitude. As the ground is approached the plume is shown to be almost aligned with the street canyon direction and independent of the incident wind directions considered. At higher altitudes well above the obstacles, the plume direction is aligned with the mean wind direction as in dispersion over flat terrain. The near-ground plume deflection is the consequence of a strong channelling effect in the region near the ground. The mean concentration profiles predicted by LES and RANS are both in good qualitative agreement with experimental data but exhibit discrepancies that can be partly explained by the influence of small incident wind angle deviation effects. Compared to RANS, LES predicts a higher channelling and thus a higher deflection of the plume. Results on the fluctuating intensity of the concentration obtained from LES show a satisfactory agreement with experiments. This information is not available from RANS for which only the mean concentration modelling is considered.     

Atmospheric boundary-layer measurements of concentration statistics from isolated and multiple sources

Atmospheric tracer dispersion experiments have been carried out to measure the statistical characteristics (variance, frequency distribution, spectrum) of the concentration downwind of a pair of partly overlapping plumes. By releasing different tracer substances from each source, it was possible to identify the contributions of the two sources at a given measurement point, both separately and jointly, and thus to compare and interpret the joint statistics in terms of those from the individual sources.Statistics for the individual sources agree well with, and support, existing wind tunnel and theoretical results. Nondimensionalization of the data using the mean concentration and the lateral width of the plume as concentration and length scales successfully removes much of the variation due to changes in atmospheric and surface conditions.Measurements of the correlation between the concentration contributions from separated sources are consistent with recent wind tunnel measurements. Entirely new measurements of the frequency distribution of the combined concentration from a pair of sources show that in many situations, high concentrations relative to the mean occur much less frequently than for an isolated source. Generally the extent of the reduction in frequency of occurrence is inversely related to the degree of correlation between concentrations from the two sources.     

Multiscaling properties of concentration fluctuations in dispersing plumes revealed using an orthonormal wavelet decomposition

The dynamical characteristics of concentration fluctuations in a dispersing plume over the energetic and inertial-convective range of scales of turbulent motion are studied using a multiscale analysis technique that is based on an orthonormal wavelet representation. It is shown that the Haar wavelet concentration spectrum is similar to the Fourier concentration spectrum in that both spectra exhibit an extensive inertial-convective subrange spanning about two decades in frequency, with a scaling exponent of -5/3. Analysis of the statistical properties (e.g., fluctuation intensity, skewness, and kurtosis) of the concentration wavelet coefficients (i.e., the concentration discrete detailed signal) suggests that the small scales are always more intermittent than the large scales. The degree of intermittency increases monotonically with decreasing scale within the inertial-convective subrange, reaching a plateau at the very small scales associated with the beginning of the near-dissipation subrange. The probability density function (pdf) of the concentration discrete detailed signal displays stretched exponential tails with an intermittency exponent (tail slope) q that increases as ^a , where is the scale or dilation and a is a power-law exponent that is dependent on downwind distance, plume height, and stratification strength with typical values in the range from about 0.25 to 0.35. It is shown that the concentration variance cascade process requires a phase coherency of eddies between different scales at the small-scale end of the inertial-convective subrange.The variation of the concentration wavelet statistics with height above the ground is investigated. The increased mean shear near the ground smooths the fine-scale plume structure for scales within the inertial-convective subrange, producing a weaker spatiotemporal intermittency in the concentration field compared to that measured higher up in the plume. The pdf of the concentration detailed signal at a fixed scale possesses less elongated tails with decreasing height z. The intermittency exponent q is found to decrease roughly linearly with increasing z.Finally, the results of the wavelet decomposition are combined to provide a conceptual model of the turbulent transport, stirring, and mixing regimes in a dispersing plume. The implications of the results for contaminant texture in a plume are discussed.     

Dispersion of a Point-Source Release of a Passive Scalar Through an Urban-Like Array for Different Wind Directions

The dispersion of a point-source release of a passive scalar in a regular array of cubical, urban-like, obstacles is investigated by means of direct numerical simulations. The simulations are conducted under conditions of neutral stability and fully rough turbulent flow, at a roughness Reynolds number of Re _τ  = 500. The Navier–Stokes and scalar equations are integrated assuming a constant rate release from a point source close to the ground within the array. We focus on short-range dispersion, when most of the material is still within the building canopy. Mean and fluctuating concentrations are computed for three different pressure gradient directions (0°, 30°, 45°). The results agree well with available experimental data measured in a water channel for a flow angle of 0°. Profiles of mean concentration and the three-dimensional structure of the dispersion pattern are compared for the different forcing angles. A number of processes affecting the plume structure are identified and discussed, including: (i) advection or channelling of scalar down ‘streets’, (ii) lateral dispersion by turbulent fluctuations and topological dispersion induced by dividing streamlines around buildings, (iii) skewing of the plume due to flow turning with height, (iv) detrainment by turbulent dispersion or mean recirculation, (v) entrainment and release of scalar in building wakes, giving rise to ‘secondary sources’, (vi) plume meandering due to unsteady turbulent fluctuations. Finally, results on relative concentration fluctuations are presented and compared with the literature for point source dispersion over flat terrain and urban arrays.     

SOME MEASUREMENTS OF TURBULENCE CHARACTERISTICS OVER COMPLEX TERRAIN

Results are presented from a windturbulence measurements campaign at the summit of ahill in complex terrain in Andros island (Greece)where a wind turbine park has been constructed. Meanturbulence parameters and power spectra from three 30 mmasts located at the summit, upwind and downwind ofthe hill are analysed to ascertain the differencesfrom respective parameters in flat terrain and toassess the influence of the irregular topography. Thevariances of the horizontal (vertical) wind components are found tobe smaller (larger) than the onesfrom flat terrain. Of the available correlations forthe power spectra, the Von Karman one was found togive better results, although the low frequency partsshow pronounced lags. The asymmetry of the upwindterrain for different wind directions and especiallyits slope is clearly felt both by variances and byspectra. This influence, as expected, is found toincrease with increasing stability. These results havean important significance for wind turbine design andwind energy utilisation.     

Flow Over Partially Forested Ridges

Numerical simulations of flow over hills that are partially covered with a forest canopy are performed. This represents a much more realistic situation than previous studies that have generally concentrated on hills that are fully-forested. The results show that the flow over the hill is sensitive to where on the hill the forest is positioned. In particular, for low slopes flow separation is predominantly located within the forest on the lee slope. This has implications for the transport of scalars in the forest canopy. For large hills the results show more variability in scalar concentrations within the canopy compared to either a fully-forested hill or a patch of forest over flat terrain. These results are likely to have implications for a range of applications including the siting and interpretation of flux measurements over forests in complex terrain, predicting wind damage to trees and wind-farm developments. Calculation of the hill-induced pressure drag and canopy-plus-surface stress shows a strong sensitivity to the position of the forest relative to the hill. Depending on the position of the forest the individual drag terms may be strongly enhanced or reduced and may even change sign. The net impact is generally to reduce the total drag compared to an equivalent fully-forested hill, but the amount of the reduction depends strongly on the position of the forest canopy on the hill. In many cases with large, wide hills there is a clear separation of scales between the adjustment of the canopy to a forest edge (of order 6 ? 8L _c, where L _c is the canopy adjustment length scale) and the width of the hill. This separation means that the hill-induced pressure and flow fields and the forest-edge induced pressure and flow fields can in some sense be considered as acting separately. This provides a means of explaining the combined effects of partial forestation and terrain. It also offers a simple method for modelling the changes in drag over a hill due to partial forest cover by considering the impact of the hill and the partial canopy separately. Scaling arguments based on this idea successfully collapse the modelled drag over a range of different hill widths and heights and for different canopy parameters. This offers scope for a relatively simple parametrization of the effects of partial forest cover on the drag over a hill.     

Concentration fluctuation measurements in a plume dispersing in a stable surface layer

A set of tracer experiments studying concentration fluctuations in a pollutant plume dispersing near the surface in a stably stratified nocturnal boundary layer is described, and the results are compared with those obtained in near-neutral stability conditions by Mylne and Mason (1991). The results highlight the importance of slow meandering of the plume which is characteristic of stable conditions. This meandering makes it impossible to conduct experiments under near-stationary conditions, resulting in considerable statistical variability in the results, but is important in reducing time-averaged concentrations. Spectral characteristics of the plume and general fluctuation statistics are qualitatively similar to those in near-neutral stability, but there are significant quantitative differences. Fluctuation time scales are shown to be substantially longer under stable conditions. This difference cannot be fully explained by the reduced windspeed alone, indicating that the length scale of plume elements is also longer. Some of the differences observed in stable conditions, particularly the longer time scales, are shown to substantially increase the potential hazard due to fluctuations in practical applications. A conceptual model of plume dispersion is described, which explains the observed plume structure under different conditions by relating it to the turbulent velocity spectra.     

Comparison of Wind-tunnel and Water-channel Simulations of Plume Dispersion through a Large Array of Obstacles with a Scaled Field Experiment

We report on measurements of the near-field dispersion of contaminant plumes in a large array of building-like obstacles at three scales; namely, at full-scale in a field experiment, at 1:50 scale in a wind-tunnel simulation, and at 1:205 scale in a water-channel simulation. Plume concentration statistics extracted from the physical modelling in the wind-tunnel and water-channel simulations are compared to those obtained from a field experiment. The modification of the detailed structure of the plume as it interacts with the obstacles is investigated. To this purpose, measurements of the evolution of the mean concentration, concentration fluctuation intensity, concentration probability density function, and integral time scale of concentration fluctuations in the array plume obtained from the field experiment and the scaled wind-tunnel and water-channel experiments are reported and compared, as well as measurements of upwind and within-array velocity spectra. Generally, the wind-tunnel and water-channel results on the modification of the detailed plume structure by the obstacles were qualitatively similar to those observed in the field experiments. However, with the appropriate scaling, the water-channel simulations were able to reproduce quantitatively the results of the full-scale field experiments better than the wind-tunnel simulations.