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.     

Critical review of the building downwash algorithms in AERMOD

The only documentation on the building downwash algorithm in AERMOD (American Meteorological Society/U.S. Environmental Protection Agency Regulatory Model), referred to as PRIME (Plume Rise Model Enhancements), is found in the 2000 A&WMA journal article by Schulman, Strimaitis and Scire. Recent field and wind tunnel studies have shown that AERMOD can overpredict concentrations by factors of 2 to 8 for certain building configurations. While a wind tunnel equivalent building dimension study (EBD) can be conducted to approximately correct the overprediction bias, past field and wind tunnel studies indicate that there are notable flaws in the PRIME building downwash theory. A detailed review of the theory supported by CFD (Computational Fluid Dynamics) and wind tunnel simulations of flow over simple rectangular buildings revealed the following serious theoretical flaws: enhanced turbulence in the building wake starting at the wrong longitudinal location; constant enhanced turbulence extending up to the wake height; constant initial enhanced turbulence in the building wake (does not vary with roughness or stability); discontinuities in the streamline calculations; and no method to account for streamlined or porous structures.

Implications: This paper documents theoretical and other problems in PRIME along with CFD simulations and wind tunnel observations that support these findings. Although AERMOD/PRIME may provide accurate and unbiased estimates (within a factor of 2) for some building configurations, a major review and update is needed so that accurate estimates can be obtained for other building configurations where significant overpredictions or underpredictions are common due to downwash effects. This will ensure that regulatory evaluations subject to dispersion modeling requirements can be based on an accurate model. Thus, it is imperative that the downwash theory in PRIME is corrected to improve model performance and ensure that the model better represents reality.     

Maximum ground-level concentrations with downwash--analysis

Equations derived previously for critical downwind distance xc' wind speed uc' and plume rise zc' the values that produce maximum ground-level concentrations (MGLC) chi c under downwash conditions, have been solved. Tables of chi c' xc' uc' and zc' and graphs of the relationships among uc and zc, for a range of stack heights hs' and building heights hb' are presented. Results for two types of sources--a turbine and a reciprocating engine--are discussed. Some comparisons are made to the U.S. Environmental Protection Agency's (EPA) SCREEN3 model.     

Comparison of the industrial source complex and AERMOD dispersion models: case study for human health risk assessment

Air quality models are typically used to predict the fate and transport of air emissions from industrial sources to comply with federal and state regulatory requirements and environmental standards, as well as to determine pollution control requirements. For many years, the U.S. Environmental Protection Agency (EPA) widely used the Industrial Source Complex (ISC) model because of its broad applicability to multiple source types. Recently, EPA adopted a new rule that replaces ISC with AERMOD, a state-of-the-practice air dispersion model, in many air quality impact assessments. This study compared the two models as well as their enhanced versions that incorporate the Plume Rise Model Enhancements (PRIME) algorithm. PRIME takes into account the effects of building downwash on plume dispersion. The comparison used actual point, area, and volume sources located on two separate facilities in conjunction with site-specific terrain and meteorological data. The modeled maximum total period average ground-level air concentrations were used to calculate potential health effects for human receptors. The results show that the switch from ISC to AERMOD and the incorporation of the PRIME algorithm tend to generate lower concentration estimates at the point of maximum ground-level concentration. However, the magnitude of difference varies from insignificant to significant depending on the types of the sources and the site-specific conditions. The differences in human health effects, predicted using results from the two models, mirror the concentrations predicted by the models.     

Maximum Ground-Level Concentrations with Downwash

Equations are derived from the Gaussian plume mode! and prescribe the critical downwind distance, wind speed, and plume rise values that result in maximum ground-level concentrations (MGLC) under downwash conditions. The derivations apply to bent-over plumes and encompass the Schulman-Scire and Huber-Snyder building downwash treatments.     

Evaluation of Downwash Modifications to the Industrial Source Complex Model

The Industrial Source Complex (ISC) model has been modified to account for Improved understanding of plume rise and downwash around buildings. Two sets of observations are used to compare the modified and the original ISC models, showing that the new model better simulates the observed concentrations during high-wind conditions when downwash usually occurs.     

Panel Discussion

A stack design procedure is developed which accounts for the effect of plume interception by downwind buildings, and which provides information on effluent concentrations in a form useful to planning authorities. The information presented in this paper is directed to engineers carrying out stack designs for locations where downwind buildings are of comparable height to the stack. A wind tunnel investigation using tracer gas techniques indicates that, for a plume at building height, downwash on the upwind face of a building causes the high concentrations observed near the roof to be transported to ground level. The effect of a plume on elevated points is determined by the concept of the minimum descent height of the maximum allowable ambient concentration isopleth. This minimum descent height, computed using Gaussian plume dispersion theory, defines a building height below which pollutant concentrations will always lie within safe limits. A case study is presented for the use of the design procedure for a small thermal power plant in an urban area.     

Modeling Emissions from Elevated Cylindrical Bins

Abstract

This paper demonstrates how wind tunnel modeling data that accurately describe plume characteristics near an unconventional emission source can be used to improve the near-field downwind plume profiles predicted by conventional air pollution dispersion models. The study considers a vertical, cylindrical-shaped, elevated bin similar to large product storage bins that can be found at many industrial plant sites. Two dispersion models are considered: the U.S. Environmental Protection Agency's ISC2(ST) model and the Ontario Ministry of the Environment and Energy's GAS model. The wind tunnel study showed that plume behavior was contrary to what was predicted using conventional dispersion models such as ISC2(ST) and GAS and default values of input parameters. The wind tunnel data were used to develop a protocol for correcting the dispersion models inputs, resulting in a substantial improvement in the accuracy of the dispersion estimates.     

Maximum Ground-Level Concentrations with Downwash—Analysis

ABSTRACT

Equations derived previously for critical downwind distance x , wind speed u , and plume rise z , the values that produce maximum ground-level concentrations (MGLC) X_c under downwash conditions, have been solved. Tables of %_c, x_c, u_c, and z_c, and graphs of the relationships among u_c and z_c for a range of stack heights h_s, and building heights h_b, are presented. Results for two types of sources— a turbine and a reciprocating engine—are discussed. Some comparisons are made to the U.S. Environmental Protection Agency's (EPA) SCREEN3 model.     

Observations Less than the Analytical Limit of Detection: A New Approach

Abstract

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.     

Flue gas discharge from cooling towers. Wind tunnel investigation of building downwash effects on ground-level concentrations

German power plants are required to meet new emission standards which limit the maximum sulfur dioxide (SOs) concentration in flue gas discharges to 400 mg m⁻³. To achieve this level of reduction in SO₂ concentration, wet scrubbing is necessary for large plants using lignite or hard coal.Wet scrubbing results in a significant reduction in the flue gas temperature leading to low effective stack heights. Instead of using stack gas reheating to achieve the plume rise necessary to satisfy local environmental standards, it was proposed to discharge the scrubbed flue gas from the existing natural-draft cooling towers (NDCT). This method should be effective in reducing local ground-level concentrations since NDCT-plumes are typically very buoyant (densimetric Froude number below 1 ) and normally reach considerable heights of rise. Only under strong wind conditions does the situation reverse itself. For such strong winds, the NDCT-plume is subject to tower and building downwash with the possibility of unacceptably high ground-level concentrations.For a 2700 MW_e lignite-fired power plant near Cologne, a wind tunnel study was carried out to investigate the effects of tower and building downwash effects on the ground-level concentrations of SO₂ produced by discharging the scrubbed flue gas from the natural-draft cooling towers. Also, a comparison was made between the ground-level concentrations produced by the cooling tower discharge method and those produced by a traditional stack. It was found that for low and intermediate wind speeds, the groundlevel concentrations are lower for the case of the cooling tower discharge. Only for strong winds, which occur only very rarely at most German sites, did the conventional stack discharge appear to be superior.     

Pseudo-source parameters for flares: Derivation,implementation, and comparison

The United States Environmental Protection Agency (US EPA) flare pseudo-source parameters are over 30 years old and few dispersion modellers understand their basis and underlying assumptions. The calculation of plume rise from the user inputs of pseudo-stack diameter, temperature and velocity have the most influence on air dispersion model predictions of ground-level concentrations. Regulatory jurisdictions across Canada, the United States and around the world have adopted their own approach to pseudo-source parameters for flares; all relate buoyancy flux to the heat release rate, none consider momentum flux and flare tip downwash as adopted by the Alberta Energy Regulator (AER). This paper derives the plume buoyancy flux for flares burning a gas in terms of combustion variables readily known or calculated without simplifying assumptions. Dispersion model prediction sensitivity to flared gas composition, temperature and velocity, and ambient conditions are now correctly handled by the AER approach. The AER flare pseudo-source parameters are based on both the buoyancy and momentum flux, thus conserving energy and momentum. The AER approach to calculate the effective source height for flares during varying wind speeds is compared to the US EPA approach. Instead of a constant source for all meteorological conditions, multiple co-located sources with varying effective stack height and diameter are used. AERMOD is run with the no stack tip downwash option as flare stack tip downwash is accounted for in the effective stack height rather than the AERMOD model calculating the downwash incorrectly using the pseudo-source parameters. The modelling approaches are compared for an example flare. Maximum ground level predictions change, generally increasing near the source and decreasing further away, with the AER flare pseudo-source parameters. It's time to update how we model flares.

Implications: What are the implications of continuing to model flare source parameters using the overly simplified US EPA approach? First, the regulators perpetuate the myths that the flare source height, temperature, diameter and velocity are constant for all wind speeds and ambient temperatures. Second, that it is acceptable to make simplifying assumptions that violate the conservation of momentum and energy principles for the sake of convenience. Finally, regulatory decisions based on simplified source modelling result in predictions that are not conservative (or realistic). The AER regulatory approach for flare source parameters overcomes all of these shortcomings. AERflare is a publicly available spreadsheet that provides the “correct” inputs to AERMOD.     

Maximum Ground-Level Concentrations with Downwash: The Urban Stability Mode

Abstract

Gaussian model-based equations for critical downwind distance, wind speed, and plume height that result in maximum ground-level concentrations (MGLC) under downwash conditions for the rural stability mode were presented in a previous paper. This paper presents general equations for the critical downwind distance xc for the urban stability mode. Specific examples are presented for Schulman-Scire and Huber-Snyder downwash treatments for building-enhanced and regular sigmas.     

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.     

Validation of ADMS against wind tunnel data of dispersion from chemical warehouse fires

Comparisons are presented of the predictions of the atmospheric dispersion modelling system (ADMS) and wind tunnel data for plume dispersion from chemical warehouse fires. The focus of the comparisons is dispersion from structurally intact buildings with open roofs and dispersion of plumes flush with the ground without obstacles, however, dispersion from building shells and doors is also considered. Both buoyancy driven and momentum driven flows are treated, although emphasis is on buoyancy driven flows as these are generally more likely to occur in warehouse fires. The study shows that the ADMS building module is able to reproduce many of the features of dispersion observed in the wind tunnel. These include a recirculating region behind the building in which material may be trapped, a main wake which brings material down towards the surface, and appropriate sensitivity to the buoyancy and momentum of the emitted material, and the location of sources on the building roof. The comparisons suggest that the ADMS building model can be used to predict dispersion from the stages of fire development studied. The precise level of agreement depends (but not in a systematic way) on the buoyancy flux parameter F_B, the momentum flux parameter F_M and the number of roof lights. There are some significant differences between the wind tunnel boundary layer and the simulated atmospheric boundary layer in ADMS which have to be considered when making wind tunnel model comparisons. These relate mainly to the near surface where the wind tunnel underestimates turbulent velocities, the boundary layer height which in the wind tunnel corresponds to an atmospheric boundary layer depth of 82.5 m (atmospheric boundary layers are frequently an order of magnitude deeper), and the boundary layer top where the ADMS boundary layer is capped by an inversion and has low turbulence levels whereas the wind tunnel boundary layer has higher levels of turbulence and no capping inversion.     

Control of Asphaltic Concrete Plants in Los Angeles County

Characteristics of maximum short-term ground level concentrations from an elevated point source, namely, the effective plume height, the critical wind speed, the distance to the point of maximum concentration, and the maximum concentration, are derived from the gaussian plume model. Both phases of plume development—before and after it has reached its final height—are considered. The plume rise treatment includes both thermal buoyancy and momentum effects. Certain limitations on critical wind speed are discussed. The dispersion model whose basis is established in this paper should be especially useful in applications where on site meteorological data are unavailable.     

Monte Carlo simulation of nitrogen oxides dispersion from a vehicular exhaust plume and its sensitivity studies

The pollutant dispersion behavior from the vehicular exhaust plume has a direct impact on human health, particularly to the drivers, bicyclists, motorcyclists, pedestrians, people working nearby and vehicle passengers. A two-dimensional pollutant dispersion numerical model was developed based on the joint-scalar probability density function (PDF) approach coupled with a k–ε turbulence model to simulate the initial dispersion process of nitrogen oxides, temperature and flow velocity distributions from a vehicular exhaust plume. A Monte Carlo algorithm was used to solve the PDF transport equations in order to obtain the dispersion distribution of nitrogen oxides concentration. The model was then validated by a series of sensitivity experimental studies in order to assess the effects of vehicular exhaust tailpipe velocities, wind speeds and chemistry on the initial dispersion of NO and NO₂ mass concentrations from the vehicular exhaust plume. The results show that the mass concentrations of nitrogen oxides decrease along the centerline of the vehicular exhaust plume in the downstream distance. The dispersion process can be enhanced when the vehicular exhaust tailpipe velocity is much larger than the wind speed. The oxidation reaction of NO plays an important role when the wind speed is large and the vehicular exhaust exit velocity is small, which leads to chemical reduction of NO, and the formation and accumulation of NO₂ in the exhaust plume. It is also found that the effect of vehicular exhaust-induced turbulence in the vicinity of the exhaust tailpipe exit is more dominant than the effect of wind turbulence, while the wind turbulence gradually shows a significant role for the dispersion of nitrogen oxides along with the development of exhaust plume. The range of dispersion of nitrogen oxides in the radial direction is increased along with the development of vehicular exhaust plume.     

Improving Plume Rise Prediction Accuracy for Stable Atmospheres with Complex Vertical Structure

Accurately predicting the rise of a buoyant exhaust plume is difficult when there are large vertical variations in atmospheric stability or wind velocity. Such conditions are particularly common near shoreline power plants. Simple plume rise formulas, which employ only a mean temperature gradient and a mean wind speed, cannot be expected to adequately treat an atmosphere whose lapse rate and wind velocity vary markedly with height. This paper tests the accuracy of a plume rise model which is capable of treating complex atmospheric structure because it integrates along the plume trajectory. The model consists of a set of ordinary differential equations, derived from the fluid equations of motion, with an entralnment parameterization to specify the mixing of ambient air into the plume. Comparing model predictions of final plume rise to field observations yields a root mean square difference of 24 m, which is 9 % of the average plume rise of 267 m. These predictions are more accurate than predictions given by three simpler models which utilize variants of a standard plume rise formula, the most accurate of the simpler models having a 12% error.     

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.     

Unsteady characteristics of the dispersion process in the vicinity of a pig barn. Wind tunnel experiments and comparison with field data

Two complementary methods, field experiments and physical modelling in a wind tunnel, have been used to investigate the dispersion of tracer-gas released from the ventilation system of a pig barn, under near-neutral stability conditions. In both cases, concentration fluctuations were measured and the deduced statistical results were compared. The choice of data processing applied to the time series of concentration was motivated by special issues in the assessment of odour annoyances: “where, how often, how long and how strong does it smell?” These features were described by the mean concentration distribution, the intermittency factor, the persistence and the 90-percentile. The good agreement between field and wind tunnel data confirmed the ability to replicate in wind tunnel the unsteady properties of a dispersion process, if the unsteady turbulent behaviour of the atmospheric boundary layer was properly modelled.A parametrical study of the influence on the dispersion process of the ratio between the exhaust velocity from the stack and the wind speed was then performed in wind tunnel. The fundamental outcome was that the near-field dispersion process under neutral stability conditions, despite the strong influence of the building wake, was for the most part driven by the meandering behaviour of the plume and not so much by the diffusion process.This study was also focused on the influence of the averaging time on the statistical results. The scatter generated by using dimensionless averaging times 200<T_a^*<400 (used during field experiments) instead of T_a^*→∞ (averaging time to ensure reproducible statistic results) was quantified in the wind tunnel. A degree of representativity of the results obtained from short-term samples, compared to fully converged statistical results was therefore assessed.