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.     

A Comparison Of The Detailed Structure In Dispersing Tracer Plumes Measured In Grid-Generated Turbulence With A Meandering Plume Model Incorporating Internal Fluctuations

A meandering plume model that explicitly incorporatesinternal fluctuations has been developed and used to model the evolutionof concentration fluctuations in point-source plumes in grid turbulenceobtained from a detailed water-channel simulation. This fluctuating plumemodel includes three physical parameters: the mean plume spread in fixedcoordinates, which represents the outer plume length scale; the meaninstantaneous plume spread in coordinates attached to the instantaneousplume centroid, which represents the inner plume length scale; and, theconcentration fluctuation intensity in the meandering reference frame,which represents the in-plume fluctuation scale. These parameters arespecified in terms of a set of coupled dynamical equations that modeltheir development with downstream distance from the source. Explicitexpressions for the concentration moments of arbitrary integral orderand the concentration probability density function have been obtainedfrom the fluctuating plume model. Detailed comparisons of model predictionsagainst water-channel measurements for the first four concentrationmoments and the concentration probability distributions generally showvery good overall quantitative agreement. Exact quantitative conditions,expressed in terms of the physical parameters of the fluctuating plumemodel, have been derived for the emergence of off-centreline peaks inthe concentration variance profile. These quantitative conditions havebeen illustrated in terms of a diagram of states of the dispersing plume,and the qualitatively different regimes of plume concentration variancebehaviour on this state diagram have been identified and characterized.     

Extension of a fluctuating plume model of tracer dispersion to a sheared boundary layer and to a large array of obstacles

Fluctuating plume models provide a useful conceptual paradigm in the understanding of plume dispersion in a turbulent flow. In particular, these models have enabled analytical predictions of higher-order concentration moments, and the form of the one-point concentration probability density function (PDF). In this paper, we extend the traditional formalism of these models, grounded in the theory of homogeneous and isotropic turbulent flow, to two cases: namely, a simple sheared boundary layer and a large array of regular obstacles. Some very high-resolution measurements of plume dispersion in a water channel, obtained using laser-induced fluorescence (LIF) line-scan techniques are utilised. These data enable us to extract time series of plume centroid position (plume meander) and dispersion in the relative frame of reference in unprecedented detail. Consequently, experimentally extracted PDFs are able to be directly compared with various theoretical forms proposed in the literature. This includes the PDF of plume centroid motion, the PDF of concentration in the relative frame, and a variety of concentration moments in the absolute and relative frames of reference. The analysis confirms the accuracy of some previously proposed functional forms of model components used in fluctuating plume models, as well as suggesting some new forms necessary to deal with the complex boundary conditions in the spatial domain.     

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.     

A Concentration pdf for the Relative Dispersion of a Contaminant Plume in the Atmosphere

Observations of the dispersion of a contaminant plume in theatmospheric boundary layer, obtained using a Lidar, are analysedin a coordinate frame relative to the instantaneous centre of massof the plume. To improve the estimates of relative dispersionstatistics, maximum entropy inversion is used to remove noise fromthe Lidar concentration profiles before carrying out the analysis.A parametric form is proposed for the probability density function(pdf) of concentration, consisting of a mixture of a betadistribution and of a generalised Pareto distribution (GPD). Thispdf allows for the possibility of a unimodal or bimodaldistribution, and is shown to give a satisfactory fit toobservations from a range of positions relative to the source. Thevariation of the fitted parameters with crossplume location isanalysed, and the maximum possible concentration is found todecrease away from the plume centre.     

A Study of Concentration Fluctuations in Instantaneous Clouds Dispersing in the Atmospheric Surface Layer for Relative Turbulent Diffusion: Basic Descriptive Statistics

A series of tracer experiments studying the statistical properties of concentration fluctuations in clouds dispersing in the atmospheric surface layer is described and analyzed. Experiments were conducted at downwind fetches between about 200 and 1200 m, under a wide range of atmospheric conditions ranging from very unstable to moderately stable stratification. The present experiments have addressed basic requirements not met by past field experiments involving instantaneously released clouds; namely, the experiments provided repeat realizations of instantaneously released clouds measured with high-resolution concentration detectors, accompanied by the contemporaneous acquisition of high-quality meteorological and turbulence measurements.Extensive analyses are performed on the cloud concentration data in the framework of relative diffusion. Ensembles of cloud concentration realizations have been constructed. From these ensembles, crosswind and time profiles of the ensemble-mean concentration, concentration variance, ensemble-mean dosage, and dosage variance are obtained. The behaviour of the time profiles of the integral time scale of cloud concentration fluctuations is studied. The use of surface-layer similarity theory for the analysis of the downwind variation of a number of cloud quantities (e.g., cloud size and duration, cloud centre ensemble-mean concentration and dosage, cloud centre concentration and dosage variance, cloud centre integral time scale) is shown to be an effective basis for ordering these quantities. Furthermore, a number of approximate universal relationships describing the behavior of these cloud quantities has been derived. Finally, it is shown that the scaled crosswind and time profiles of ensemble-mean concentration and concentration variance as well as the scaled time profiles of the concentration fluctuation integral time scale exhibit self-similar forms that are independent of atmospheric stratification and downwind fetch.     

Simulation of Time Series of Concentration Fluctuations in Atmospheric Dispersion Using a Correlation-distortion Technique

Short-duration fluctuations in the concentration of airborne substances can be important in a variety of atmospheric dispersion problems, especially when assessing the risks posed by harmful materials. This paper discusses a simulation technique for generating individual realisations of fluctuating concentration time series in dispersing plumes based on target probability distributions and spectral functions. The scheme uses a correlation-distortion approach to simulate these time series. Gaussian processes with modified spectral characteristics are generated and then transformed to yield non-Gaussian processes with the desired spectral characteristics. The simulation scheme is initially developed for a single receptor, and is then extended to model pairs of correlated time series at two receptors. In fact, the modelling technique can be generalised to an arbitrary number of receptors and this provides, in principal, an approach that is applicable to a wide class of similar problems (such as the modelling of instantaneous puff releases or the response of line-of-sight detection systems). The simulation technique is illustrated using observations made during recent field experiments, conducted both in the United Kingdom and in the U.S.A., investigating the short-range dispersion of a passive tracer.     

The Interference of Higher-Order Statistics of the Concentration Field Produced by Two Point Sources According to a Generalized Fluctuating Plume Model

The higher-order correlation functions for the concentrationfluctuations arising from a two-point-source configuration have beencalculated analytically within the context of the phenomenology of afluctuating plume model (viz., a meandering plume model that explicitlyincorporates internal fluctuations). Explicit expressions for thesecond-, third-, and fourth-order correlationfunctions between the concentrationfluctuations produced by two point sources are given in terms of the sourceseparation d and the five physically based parameters that define thegeneralized fluctuating plume model: namely, the absolute plume dispersion,_a, which determines the outer plume length scale; the relative plume dispersion, _r, which determines the inner plume length scale; the fluctuation intensity, i_r, in relative coordinates, which determines the internal concentration fluctuation level; the correlation coefficient, r,between the positions of the centroids of the two interfering plumes; and,the correlation coefficient, r^*, between the concentration fluctuationsof the two plumes in relative coordinates, which determines the degree ofinternal mixing of the two scalars. Furthermore, the form of the totalconcentration probability density function arising from the interferenceproduced by two point sources is presented. Predictions for the second-ordercorrelation function, , and for the total concentration probabilitydensity function have been compared with some new experimental data fora two-point-source configuration in grid turbulence generated in awater-channel simulation. These results are in good agreement with the dataand suggest that the analytical model for the second-order correlationfunction and the total concentration probability density function canreproduce many qualitative trends in the interaction of plumes from twosources.     

A Stochastic Time Series Model for Threshold Crossing Statistics of Concentration Fluctuations in non-Intermittent Plumes

A numerical stochastic model is developed for the upcrossing rate across a specified threshold concentration. The model assumes that the concentration time series at a given spatial point within a dispersing plume can be approximated as a first-order Markovian process designed to be consistent with a given time-invariant concentration probability density function (pdf). The model requires only the specification of a concentration pdf with a given mean and variance and a concentration fluctuation integral time scale. Predicted upcrossing rates are compared with atmospheric plume concentration data obtained from a point source near the ground. For this data set, a log-normal pdf is found to give better estimates of the threshold crossing rate than a gamma pdf.     

Spatial Resolution in the Measurement of Concentration Fluctuations

In a turbulent flow, a miscible contaminant is confined to sheets and strands of the very thin Batchelor conduction cut-off length. This fact has been surmised for some time and has recently been observed directly through high precision measurements. This fine-scaled texture of the contaminant concentration field makes it an extraordinary challenge to achieve (experimentally) adequate continuum scale resolution, particularly in important environmental flows such as the atmospheric boundary layer. In this paper, an extrapolation scheme is proposed whereby the systematic measurement (with known sample volumes) of the lower-order moments of the concentration fluctuations are used to approximate the true, perfectly resolved values, and hence to approximate the true probability density function. Such a scheme relies on empirical data, and so the need for more experiments designed to investigate the effects of spatial resolution cannot be over-emphasized.     

Computational Fluid Dynamical Modelling of Concentration Fluctuations in an Idealized Urban Area

Turbulent mixing induces variability in concentration that is important in many applications, such as reactive plumes, risk assessments or odour impact analyses (when the effects can have time scales on the order of a second). In urban canopies, the variability may be modified by the presence of buildings. Our purpose is to study concentration fluctuation variance in built-up areas using an Eulerian approach. We performed numerical simulations with the computational fluid dynamics model Mercure_Saturne, which is a three-dimensional model adapted to atmospheric flow and pollutant dispersion. We use a k − ϵ turbulence closure and predict the concentration variance with a transport equation model. The model performance is evaluated with the near-full scale experiment MUST (Mock Urban Setting Test), a field experiment conducted in Utah’s West Desert Test Center. The modelled root-mean-square of the concentration fluctuations is compared to measurements for 20 of the MUST trials. The model shows good agreement with the measurements, with the fraction of predictions within a factor of two of observations of 60.1%, with better results for horizontal lines of detectors than for the detectors on vertical masts (with fractions of predictions within a factor of two of observations of respectively 66.4% and 52.6%). The influence of different parameters on the fluctuation variance is also studied and we show the importance of taking into account the stability of the stratification when modelling the turbulent kinetic energy.     

A Simple Model of Concentration Fluctuations in Neutrally Buoyant Clouds

This paper presents a new model of concentration fluctuations for neutrally buoyant gas clouds dispersing in a wind tunnel. It is derived from a series of exact results, which apply in the hypothetical case when there is no molecular diffusion, coupled with a probability density function model previously used to describe steady releases of contaminant. A simple self-similar relationship between the evolution of the concentration intensity and mean is established. As a first step the time independent variant of the model, applicable to a continuous plume, is tested against some previously published experimental data for steady wind-tunnel releases. Comparisons of experimental results and model predictions at different downwind positions, heights and source geometry are presented. Then, the results for the time dependent model, applicable to instantaneous releases, are discussed. The experimental evidence presented here supports the self-similar relationship established earlier. The implications for modelling higher moments of concentration and the fixed point probability density function are investigated.     

Probabilistic Model for Concentration Fluctuations in Compact-Source Plumes in an Urban Environment

A comprehensive model for the prediction of concentration fluctuations in plumes dispersing in the complex and highly disturbed wind flows in an urban environment is formulated. The mean flow and turbulence fields in the urban area are obtained using a Reynolds-averaged Navier-Stokes (RANS) flow model, while the standard k-ϵ turbulence model (k is the turbulence kinetic energy and ϵ is the viscous dissipation rate) is used to close the model. The RANS model provides a specification of the velocity statistics of the highly disturbed wind flow in the urban area, required for the solution of the transport equations for the mean concentration and concentration variance (both of which are formulated in the Eulerian framework). A physically-based formulation for the scalar dissipation time scale t _d , required for the closure of the transport equation for , is presented. This formulation relates t _d to an inner time scale corresponding to “internal” concentration fluctuation associated with relative dispersion, rather than an outer time scale associated with the entire portion of the fluctuation spectrum. The two lowest-order moments of concentration ( and ) are used to determine the parameters of a pre-chosen functional form for the concentration probability density function (clipped-gamma distribution). Results of detailed comparisons between a water-channel experiment of flow and dispersion in an idealized obstacle array and the model predictions for mean flow, turbulence kinetic energy, mean concentration, concentration variance, and concentration probability density function are presented.     

Micrometeorological Measurements in a Street Canyon during the Joint ATREUS-PICADA Experiment

In order to investigate the microclimatic conditions in a street canyon, a physical model was used to conduct the Joint ATREUS-PICADA Experiment (JAPEX) in situ experimental campaign. Four lines of buildings simulated by steel containers were installed to form three parallel street canyons at 1:5 scale, with width/height aspect ratio approximately 0.40. The reference wind and atmospheric conditions were measured, as well as the flow velocity and direction in the street. Preliminary results concern street canyon ventilation and thermal effects on in-canyon airflow, and show that vortical motions appear for reference wind directions perpendicular to the street axis. The presence of adjacent rows of buildings did not appear to significantly influence the flow character within the canyon for the case of a low aspect ratio corresponding to a skimming flow regime. The flow structure was not significantly affected by the thermal effects although some slight interference occurred in the lower part of the canyon. An analysis of horizontal temperature gradients indicated that a thin boundary layer develops near the heated facade. These facts imply that the thermal effects are considerable only very close to the wall.     

Measurements of Damping of Temperature Fluctuations in a Tube

Using field measurements we evaluate a previously deduced semi-empirical model for the dampening of temperature fluctuations in a circular tube. The measurements show that the model is a reasonably good approximation, and the previous recommendation of using a tube length of the order of 1,000 times the inner diameter of the tube to remove most of the temperature fluctuations seems to hold. Part of the difference between the empirical results and the model is possibly due to deviations from the idealized conditions assumed in the model such as no inclusion of the influence from the ambient atmospheric stability. Limitations of the measurements are also likely to affect the evaluation.     

Wind-tunnel Modelling of Dispersion from a Scalar Area Source in Urban-Like Roughness

A wind-tunnel study was conducted to investigate ventilation of scalars from urban-like geometries at neighbourhood scale by exploring two different geometries a uniform height roughness and a non-uniform height roughness, both with an equal plan and frontal density of λ _p = λ _f = 25%. In both configurations a sub-unit of the idealized urban surface was coated with a thin layer of naphthalene to represent area sources. The naphthalene sublimation method was used to measure directly total area-averaged transport of scalars out of the complex geometries. At the same time, naphthalene vapour concentrations controlled by the turbulent fluxes were detected using a fast Flame Ionisation Detection (FID) technique. This paper describes the novel use of a naphthalene coated surface as an area source in dispersion studies. Particular emphasis was also given to testing whether the concentration measurements were independent of Reynolds number. For low wind speeds, transfer from the naphthalene surface is determined by a combination of forced and natural convection. Compared with a propane point source release, a 25% higher free stream velocity was needed for the naphthalene area source to yield Reynolds-number-independent concentration fields. Ventilation transfer coefficients w _T /U derived from the naphthalene sublimation method showed that, whilst there was enhanced vertical momentum exchange due to obstacle height variability, advection was reduced and dispersion from the source area was not enhanced. Thus, the height variability of a canopy is an important parameter when generalising urban dispersion. Fine resolution concentration measurements in the canopy showed the effect of height variability on dispersion at street scale. Rapid vertical transport in the wake of individual high-rise obstacles was found to generate elevated point-like sources. A Gaussian plume model was used to analyse differences in the downstream plumes. Intensified lateral and vertical plume spread and plume dilution with height was found for the non-uniform height roughness.