Validation of a Lagrangian model plume rise scheme using the Kincaid data set

Correct prediction of the initial rise of a plume due to momentum and buoyancy effects is an important factor in dispersion modelling. A new plume rise scheme, based upon conservation equations of mass, momentum and heat, for the Lagrangian model, NAME, is described. The conservation equations are consistent with the well-known analytical plume rise formulae for both momentum- and buoyancy-dominated plumes. The performance of the new scheme is assessed against data from the Kincaid field experiment. Results show that the new scheme adds value to the model and significantly outperforms the previous plume rise scheme. Using data from assessments of atmospheric dispersion models using the Kincaid data set, it is shown that NAME is comparable to other models over short ranges.     

Stable plume rise in a shear layer

Solutions are given for plume rise assuming a power-law wind speed profile in a stably stratified layer for point and finite sources with initial vertical momentum and buoyancy. For a constant wind speed, these solutions simplify to the conventional plume rise equations in a stable atmosphere. In a shear layer, the point of maximum rise occurs further downwind and is slightly lower compared with the plume rise with a constant wind speed equal to the wind speed at the top of the stack. If the predictions with shear are compared with predictions for an equivalent average wind speed over the depth of the plume, the plume rise with shear is higher than plume rise with an equivalent average wind speed.     

A three-dimensional plume rise model for dry and wet plumes

We present a plume rise model which can be applied to situations with arbitrary wind fields and source exit directions and to both dry and wet plumes. The model is an integral model which considers plume properties averaged over the plume cross section. It is validated by means of water tank, wind tunnel, and field experiments (stacks and cooling towers).     

Cooling tower plume rise analyses by airborne lidar

As part of the Atmospheric Studies in Complex Terrain (ASCOT) program, five cooling-tower plume experiments were conducted at The Geysers geothermal area in California during August 1981. The experiments were designed to investigate plume behavior during conditions of valley nocturnal drainage flow and daytime vertical mixing on a mountain ridge. The Airborne Lidar Plume and Haze Analyzer (ALPHA-1) system was used during the experiments to observe plume geometry and aerosol structure of the boundary layer. Subvisible plume rise derived from the backscatter signatures is related to visible plume-rise results recently published. This analysis indicates that the lidar-observed plumes are about 10 times higher than visually observed plumes. The lidar-observed plume rise seemed to correspond with heights of elevated aerosol layers that typically indicate presence of temperature inversions. Pictorial displays are presented which provide information on atmospheric behavior over mountainous terrain.     

The effects of parametric uncertainties in simulations of a reactive plume using a Lagrangian stochastic model

A combined Lagrangian stochastic model with micro-mixing and chemical sub-models is used to investigate a reactive plume of nitrogen oxides (NO_x) released into a turbulent grid flow doped with ozone (O₃). Sensitivities to the model input parameters are explored for different source NO_x scenarios. The wind tunnel experiments of Brown and Bilger (1996) provide the simulation conditions for the first case study where photolysis reactions are not included and the main uncertainties occur in parameters defining the turbulence scales, source size and reaction rate of NO with O₃. Using nominal values of the parameters from previous studies, the model gives a good representation of the radial profile of the conserved mean scalar ${\overline{\Gamma }}_{{\text{NO}}_{\text{x}}}$ although slightly over predicts peak mean NO₂ concentrations ${\overline{\Gamma }}_{{\text{NO}}_{\text{2}}}$ compared to the experiments. The high dimensional model representation (HDMR) method is used to investigate the effects of uncertainties in model inputs on the simulation of chemical species concentrations. For this scenario, the Lagrangian velocity structure function coefficient has the largest impact on simulated ${\overline{\Gamma }}_{{\text{NO}}_{\text{x}}}$ profiles. Photolysis reactions are then included in a chemical scheme consisting of eight reactions between species NO, O, O₃ and NO₂. Independent and interactive effects of 22 input parameters are studied for two source NO_x scenarios using HDMR, including turbulence parameters, temperature dependant rate parameters, photolysis rates, temperature, fraction of NO in total NO_x at the source and background ozone concentration [O₃]. For this reactive case, the variance in the predicted mean plume centre ${\overline{\Gamma }}_{{\text{O}}_{\text{3}}}$ is caused by parameters describing both physical (mixing time-scale coefficient) and chemical processes (activation energy for the reaction O₃+NO). The variance in predicted plume centre ${\overline{\Gamma }}_{{\text{NO}}_{\text{2}}}$ and root mean square NO₂ concentration ${\gamma }_{{\text{NO}}_{\text{2}}}^{\prime }$, is strongly influenced by the fraction of NO in the source NO_x, and to a lesser extent the mixing time-scale coefficient. Adjusting the latter gives improved agreement with the Brown and Bilger experiment. Some weak parameter interactions are observed.     

Development and evaluation of the PRIME plume rise and building downwash model

A new Gaussian dispersion model, the Plume Rise Model Enhancements (PRIME), has been developed for plume rise and building downwash. PRIME considers the position of the stack relative to the building, streamline deflection near the building, and vertical wind speed shear and velocity deficit effects on plume rise. Within the wake created by a sharp-edged, rectangular building, PRIME explicitly calculates fields of turbulence intensity, wind speed, and streamline slope, which gradually decay to ambient values downwind of the building. The plume trajectory within these modified fields is estimated using a numerical plume rise model. A probability density function and an eddy diffusivity scheme are used for dispersion in the wake. A cavity module calculates the fraction of plume mass captured by and recirculated within the near wake. The captured plume is re-emitted to the far wake as a volume source and added to the uncaptured primary plume contribution to obtain the far wake concentrations. The modeling procedures currently recommended by the U.S. Environmental Protection Agency (EPA), using SCREEN and the Industrial Source Complex model (ISC), do not include these features. PRIME also avoids the discontinuities resulting from the different downwash modules within the current models and the reported overpredictions during light-wind speed, stable conditions. PRIME is intended for use in regulatory models. It was evaluated using data from a power plant measurement program, a tracer field study for a combustion turbine, and several wind-tunnel studies. PRIME performed as well as or better than ISC/SCREEN for nearly all of the comparisons.     

A puff pollutant dispersion model with wind shear and dynamic plume rise

A puff diffusion model, which includes wind shear and dynamic plume rise, is developed for numerical prediction of pollutant concentrations under unsteady and non-uniform flow conditions. The plume from a continuous source is treated as a series of puffs emitted successively from the source. Each puff is represented by a set of six tracer particles, which define the size, shape and location of the puff. Initially these particles are located at the surface of the source, on arbitrarily chosen orthogonal axes. The location of the particles is computed at each time step by taking into account advection, eddy diffusion, wind shear and entrainment of ambient air during plume rise. The concentration distribution of each puff is determined by fitting an ellipsoid to the cluster of the six particles and assuming a three-dimensional Gaussian distribution, with standard deviation equal to the half-lengths of the principal axes of the ellipsoid. The concentration at a point of interest is obtained by summing the contributions from nearby puffs. The effect of wind shear on the pollutant concentration is investigated by use of a typical wind shear encountered in the atmosphere. The results show that, at 600 m downstream from the source, the present model gives concentrations a factor of 2 higher and lower at one standard deviation below and above the plume center, respectively, than that of conventional models in which no wind shear is considered. The plume-rise formulation is calibrated against the observations compiled by Briggs and the model is used to predict the trajectory of a plume observed by Slawson and Csanady. Excellent agreement between the prediction and the observation can be achieved if an appropriate eddy diffusivity is chosen.     

Cooling tower visible plume rise analyses by time integrated photographs

Cooling tower plume rise measurements are required for evaluation and validation of cooling tower plume rise models. In August 1981, five cooling tower plume defining experiments were conducted at the Geysers Geothermal Area in California during both valley nocturnal drainage flow conditions, and diurnal limited vertical mixing conditions on a mountain ridge. Sequential 300-s time integrated photographs of the plumes were made for several hours during each experiment. Time lapse movies were also made during the daytime experiments. Meteorological conditions were measured during all experimental periods by tethersonde and rawinsonde. The time integrated photographs have been analyzed to provide values of plume rise heights vs time for each of the experiments. Predictions by one version of the Argonne National Laboratory plume rise model are within a factor of 2 of the plume rise calculated from the photographs.     

Numerical and approximate solutions for plume rise

Numerical and approximate analytical solutions are compared for turbulent plume rise in a crosswind. The numerical solutions were calculated using the plume rise model of Hoult, Fay and Forney (1969, J. Air Pollut. Control Ass.19, 585–590), over a wide range of pertinent parameters. Some wind shear and elevated inversion effects are included. The numerical solutions are seen to agree with the approximate solutions over a fairly wide range of the parameters. For the conditions considered in the study, wind shear effects are seen to be quite small. A limited study was made of the penetration of elevated inversions by plumes. The results indicate the adequacy of a simple criterion proposed by Briggs (1969, AEC Critical Review Series, USAEC Division of Technical Information extension, Oak Ridge, Tennesse).     

Modelling plume rise and dispersion from pool fires

This paper describes an investigation into the behaviour of smoke plumes from pool fires, and the subsequent generation of empirical models to predict plume rise and dispersion from such a combustion source. Synchronous video records of plumes were taken from a series of small-scale (0.06–0.25m²) outdoor methanol/toluene pool fire experiments, and used to produce sets of images from which plume dimensions could be derived. Three models were used as a basis for the multiple regression analysis of the data set, in order to produce new equations for improved prediction. Actual plume observations from a large (20.7 m×14.2 m) aviation fuel pool fire were also used to test the predictions. The two theoretically based models were found to give a better representation of plume rise and dispersion than the empirical model based on measurements of small-scale fires. It is concluded that theoretical models tested on small-scale fires (heat output ≈70 kW) can be used to predict plume behaviour from much larger combustion sources (heat output ≈70 MW) under near neutral atmospheric conditions.     

Effective stratification for plume rise and boundary layer growth computations

A method is proposed for evaluating the effective stratification compatible with well-known plume rise and entrainment growth formulations when the vertical distribution of density differs from idealized conditions. In addition consideration is given to reducing a measured non-constant profile of wind speed to one compatible with theory. A nomogram method based on the commonalities of the developments is offered for small volume data reduction.     

Vertical emission profiles for Europe based on plume rise calculations

The vertical allocation of emissions has a major impact on results of Chemistry Transport Models. However, in Europe it is still common to use fixed vertical profiles based on rough estimates to determine the emission height of point sources. This publication introduces a set of new vertical profiles for the use in chemistry transport modeling that were created from hourly gridded emissions calculated by the SMOKE for Europe emission model. SMOKE uses plume rise calculations to determine effective emission heights. Out of more than 40 000 different vertical emission profiles 73 have been chosen by means of hierarchical cluster analysis. These profiles show large differences to those currently used in many emission models. Emissions from combustion processes are released in much lower altitudes while those from production processes are allocated to higher altitudes. The profiles have a high temporal and spatial variability which is not represented by currently used profiles.     

A model for the maximum credible hourly impact on any ground receptor from point sources with thermal plume rise

A pollutant dispersion model is developed, allowing fast evaluation of the maximum credible 1-h average concentration on any given ground-level receptor, along with the corresponding critical meteorological conditions (wind speed and stability class) for stacks with buoyant plumes in urban or rural areas. Site-specific meteorological data are not required, as the computed concentrations are maximized against all credible combinations of wind speed, stability class, and mixing height. The analysis is based on the dispersion relations of Pasquill-Gifford and Briggs for rural and urban settings, respectively, the buoyancy induced dispersion correlation of Pasquill, the wind profile exponent values suggested by Irwin, the buoyant plume rise relations of Briggs, as well as the Benkley and Schulman's model for the minimum mixing heights. The model is particularly suited for air pollution management studies, as it allows fast screening of the maximum impact on any selected receptor and evaluation of the ways to have this impact reduced. It is also suited for regulatory purposes, as it can be used to define the minimum stack size requirements for a given source as a function of the exit gas volume and temperature, the pollutant emission rates and their hourly concentration standards, as well as the source location relative to sensitive receptors.     

Water tank measurements of buoyant plume rise and structure in neutral crossflows

In this paper, an experimental study of the rise and development of a single buoyant plume and a pair of in-line buoyant plumes is presented. The investigations were carried out at small scale in a water filled towing tank using both quantitative flow visualisation and local concentration measurements. The measured plume trajectories for a single plume were compared with the Briggs plume rise equation and predictions from a numerical integral model. Plume trajectories were studied for twin in-line plumes, with particular attention to changes in the plume trajectory, especially any additional rise that resulted from the interaction between the two plumes. Concentration field distributions in cross-sections through both single and interacting twin plumes were obtained from the local concentration measurement system. These showed how the interaction affected the plume structure, notably the double vortex system that occurs in a fully developed plume.