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.     

CFD model simulation of dispersion from chlorine railcar releases in industrial and urban areas

To assist in emergency response decisions and planning in case of releases of pressurized liquefied chlorine from railroad tank cars in industrial sites and cities, the FLACS Computational Fluid Dynamics (CFD) model has been used to simulate the transport and dispersion of the dense chlorine cloud. Two accident locations are studied: an actual railcar accident at an industrial site in Festus, MO, and a hypothetical railcar accident at a rail junction in the Chicago urban area. The results show that transport of a large dense gas release at ground level in an industrial site or large city could initially extend a hundred meters or more in the upwind and crosswind directions. The dense cloud may follow terrain drainage, such as river channels. Near the source, the obstacles tend to slow down the dense gas cloud and may constrain it and cause increased concentrations. Farther downwind, the obstacles may cause enhanced mixing and dilution once the cloud has grown larger. In some cases, significant amounts of cloud mass may become “trapped” in obstacle wakes for many minutes after the main cloud has passed. Although the CFD model can account for the details of the flow and dispersion much better than standard widely-used simple dense gas models, many similarities are found among the various models in their simulated variations with downwind distance of the maximum cloud centerline concentration.     

Overview of Petroleum Environmental Research Forum (PERF) dense gas dispersion modeling project

This paper provides a background for and an overview of the results of a comprehensive study of transport and dispersion of dense gas plumes over rough surfaces typical of industrial sites. The Petroleum Environmental Research Forum (PERF) 93-16 project involved model development and evaluations using observations from three wind tunnels and from the Kit Fox field experiment. Detailed discussions of the results of the research are given in the other papers in this special issue. The wind tunnel experiments produced data showing that the resulting best-fit vertical entrainment formula was close to (i.e., within about 30%) the vertical entrainment formulas already in use by current models, which were derived primarily from observations over smooth surfaces. Observations from the Kit Fox field experiment demonstrated the validity of the entrainment curves derived from the wind tunnel data. The Kit Fox data were also used to evaluate algorithms for along-wind dispersion and cloud advection speeds for short-duration releases typical of an industrial site, and to evaluate the HEGADAS dense gas dispersion model.     

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.     

Plume rise in a shear layer with neutral stability

Several solutions have been published to predict the rise of buoyant plumes in a shear layer with a power-law velocity profile. Each of these solutions is either a special case or is based on oversimplifying assumptions. In this paper, solutions to the plume-rise equations are given for buoyant and nonbuoyant plumes with initial vertical momentum. Solutions are given for both point sources and sources with a finite initial size under neutral stability. For a constant wind speed, these solutions simplify to the conventional plume-rise equations.     

Use of the Kit Fox field data to analyze dense gas dispersion modeling issues

The 1995 Kit Fox dense gas field data set consists of 52 trials where short-duration CO₂ gas releases were made at ground level over a rough surface during neutral to stable conditions. The experiments were intended to demonstrate the effects on dense gas clouds of relatively large roughnesses typical of industrial process plants. Fast response concentration observations were made by 80 samplers located on four downwind lines (25, 50, 100, and 225 m), including profile observations on three towers on each of the closest three arcs. Detailed meteorological measurements were made on several towers within and outside of the roughness arrays. The data analysis emphasized the variation of maximum concentration with surface roughness, the dependence of cloud advection speed on cloud depth, the variation of the three components of dispersion with ambient turbulence, and the dependence of vertical entrainment rate on ambient friction velocity and cloud Richardson number. The Kit Fox data were used to evaluate a specific dense gas dispersion model (HEGADAS 3+), with emphasis on whether it would be able to account for the increased roughness. The model was able to satisfactorily simulate the observed concentrations, with a mean bias of about 5% and with about 90% of the predictions within a factor of two of the observations.     

A controlled field experiment on groundwater contamination by a multicomponent DNAPL: creation of the emplaced-source and overview of dissolved plume development

A unique field experiment has been undertaken at the CFB Borden research site to investigate the development of dissolved chlorinated solvent plumes from a residual dense non-aqueous phase liquid (DNAPL) source. The "emplaced-source" tracer test methodology involved a controlled emplacement of a block-shaped source of sand containing chlorinated solvents below the water table. The gradual dissolution of this residual DNAPL solvent source under natural aquifer conditions caused dissolved solvent plumes of trichloromethane (TCM), trichloroethene (TCE) and perchloroethene (PCE) to continuously develop down gradient. Source dissolution and 3-D plume development were successfully monitored via 173 multilevel samplers over a 475-day tracer test period prior to site remediation research being initiated. Detailed groundwater level and hydraulic conductivity data were collected. Development of plumes with concentrations spanning 1-700,000 micrograms/1 is described and key processes controlling their migration identified. Plumes were observed to be narrow due to the weakness of transverse dispersion processes and long due to advection and significant longitudinal dispersion, very limited sorptive retardation and negligible, if any, attenuation due to biodegradation or abiotic reaction. TCM was shown to be essentially conservative, TCE very nearly conservative and PCE, consistent with its greater hydrophobicity, more retarded yet having a greater mobility than observed in previous Borden field tests. The absence of biodegradation was ascribed to the prevailing aerobic conditions and lack of any additional biodegradable carbon substrates. The transient groundwater flow regime caused significant transverse lateral plume movement, plume asymmetry and was likely responsible for most of the, albeit limited, transverse horizontal plume spreading. In agreement with the widespread incidence of extensive TCE and PCE plumes throughout the industrialized world, the experiment indicates such solvent plumes are likely to be highly mobile and persistent, at least in aquifers that are aerobic and have low sorption potential (low foc content).     

Atmospheric dispersion in the presence of a three-dimensional cubical obstacle: Modelling of mean concentration and concentration fluctuations

This paper presents computational simulations of atmospheric dispersion experiments conducted around isolated obstacles in the field. The computational tool used for the simulations was the code ADREA-HF, which was especially developed for the simulation of the dispersion of positively or negatively buoyant gases in complicated geometries. The field experiments simulated involve a single cubic obstacle normal to the mean wind direction and two upwind sources of ammonia and propane, with the ammonia source located at different lateral positions [Mavroidis, I., Griffiths, R.F., Hall, D.J., 2003. Field and wind tunnel investigations of plume dispersion around single surface obstacles. Atmospheric Environment 37, 2903–2918]. Concentrations and concentration fluctuations for both gases were calculated by the model and compared with the experimental results. Certain modelling aspects were studied, such as the effect of using different turbulence closure schemes in the computations. Furthermore, specific characteristics of dispersion were investigated using the computational tool, such as the effect of the lateral displacement of a source on the concentration fluctuations intensity, the effects of natural variability and the sensitivity of concentrations to wind direction fluctuations. The results showed a good level of agreement between calculated and measured concentrations and concentration fluctuations when ensemble averaged data were available from the field experiments. Differences observed between measured and predicted concentrations and concentration fluctuations, in the case of laterally displaced sources, were mainly attributed to the specificities of the experimental cases, such as the interaction of a laterally displaced plume with an obstacle, and to the variability observed in the field. The effect of this variability is indicated by the difference between the predicted-to-observed ratios of ensemble-averaged centreline values for propane and the respective ratios from the single ammonia experiment with co-located gas sources, the latter being higher by upto 30% for concentrations and 70% for concentration fluctuations. Using the computational tool it was shown that, for a laterally displaced source, a change of 5° in the mean wind direction can lead up to a 100% variation in the measured concentrations.     

Remote monitoring of air pollutant emissions from point sources by a mobile lidar/sodar system

This paper describes remote monitoring of air pollutant emissions by a mobile lidar (light detection and ranging)/ sodar (sound detection and ranging) system. First, measurements are carried out in the flue gas plume of a public power plant. The investigations focus mainly on quantifying SO2 emissions, but the uncertainties of such measurements are also emphasized. Furthermore, an example providing valuable data sets for the development and validation of plume dispersion models is outlined with measurements of the dilution of SO2 along the plume axis. Series of repeated determinations of SO2 emissions show a large variation in the obtained flux values, with moderate margins of error. Incomplete recording of the plume within the individual lidar scans, induced by strong looping movements of the flue gas plume, predominantly causes the variations of flux values. Therefore, the highest flux values determined are considered to be the most exact. This is verified by a comparison of measured fluxes with in situ measurements made by the plant operators. The results further indicate that lidar measurements illustrate the location and dimension of aerosol plumes better than the location and dimension of the plumes of gaseous compounds. The wind direction affecting the plume at any moment can be determined faster by lidar than by sodar because the latter requires much longer time intervals of signal averaging. Measurements show higher concentrations of SO2 compared with results from a Gaussian plume model for periods of less than 5 min after dispersion. The findings emphasize the suitability of remote sensing for detecting emissions and for investigating the propagation and dilution of air pollutant plumes.     

Experiments on vertical transverse mixing in a large-scale heterogeneous model aquifer

Vertical transverse mixing is known to be a controlling factor in natural attenuation of extended biodegradable plumes originating from continuously emitting sources. We perform conservative and reactive tracer tests in a quasi two-dimensional 14 m long sand box in order to quantify vertical mixing in heterogeneous media. The filling mimics natural sediments including a distribution of different hydro-facies, made of different sand mixtures, and micro-structures within the sand lenses. We quantify the concentration distribution of the conservative tracer by the analysis of digital images taken at steady state during the tracer-dye experiment. Heterogeneity causes plume meandering, leading to distorted concentration profiles. Without knowledge about the velocity distribution, it is not possible to determine meaningful vertical dispersion coefficients from the concentration profiles. Using the stream-line pattern resulting from an inverse model of previous experiments in the sand box, we can correct for the plume meandering. The resulting vertical dispersion coefficient is approximately approximately 4 x 10(-)(9) m(2)/s. We observe no distinct increase in the vertical dispersion coefficient with increasing travel distance, indicating that heterogeneity has hardly any impact on vertical transverse mixing. In the reactive tracer test, we continuously inject an alkaline solution over a certain height into the domain that is occupied otherwise by an acidic solution. The outline of the alkaline plume is visualized by adding a pH indicator into both solutions. From the height and length of the reactive plume, we estimate a transverse dispersion coefficient of approximately 3 x 10(-)(9) m(2)/s. Overall, the vertical transverse dispersion coefficients are less than an order of magnitude larger than pore diffusion coefficients and hardly increase due to heterogeneity. Thus, we conclude for the assessment of natural attenuation that reactive plumes might become very large if they are controlled by vertical dispersive mixing.     

A simple model for the dispersion of pollutants from a road tunnel portal

The dispersion of pollutants from a roadway tunnel portal is mainly determined by the interaction between the ambient wind and the jet stream from the tunnel portal. In principal, Eulerian microscale models by solving the conservation equations for mass, momentum, and energy, are thus able to simulate effects such as buoyancy etc. properly. However, for engineering applications such models need too much CPU time, and are not easy to handle by non-scientific personnel. Only a few dispersion models, applicable for regulatory purposes, have so far appeared in the literature. These models are either empirical models not always applicable for different sites, or they do not capture important physical effects like buoyancy phenomena. Here, a rather simple model is presented, which takes into account most of the important processes considered to govern the dispersion of a jet stream from portals. These are the exit velocity, the buoyancy, the influence of ambient wind direction fluctuations on the position of the jet stream, and traffic induced turbulence. Although the model contains some heuristic elements, it was successfully tested against tracer experiments taken near a motorway tunnel portal in Austria. The model requires relatively little CPU time. Current limitations of the model include the neglect of terrain, building, and vehicle effects on the dispersion, and the neglect of the horizontal dispersion arising from entrainment of ambient air in the jet stream. The latter could lead to an underestimation of plume spreads for higher wind speeds. The validation of the model will be the focus of future research. The experimental data set is also available for the scientific community.     

The absorption of Hydrogen chloride by aqueous aerosols

A general gas/droplet interaction model has been used to examine the absorption of gaseous hydrogen chloride in aqueous aerosols formed in the presence of typical ambient concentrations of soluble and insoluble particulate over the relative humidity range 60–99 % and temperature range 273–283 K. Absorption of hydrogen chloride in aqueous aerosols formed on inert fugitive ash from a point source has also been simulated using simple puff and constant angle dispersion models. For ambient atmospheres it is predicted that appreciable uptake of hydrogen chloride occurs only for relative humidities approaching 99%. Equilibration times increase with particle size and gas concentration and lie in the range < 1–600 s. For diluting point source emissions, it is predicted that conditions close to equilibrium exist throughout most of the dispersion period. Hydrogen chloride absorption increases with initial plume water content and ambient relative humidity and decreases with ambient temperature. Substantial temporary removal of gaseous hydrogen chloride is predicted under most conditions simulated but the period of uptake is significant only for slowly diluting plumes at ambient relative humidities > 95%. Initial plume particulate diameters > 5 μm are required for enhanced deposition and substantial permanent removal of gaseous HCl. Absorption in aqueous ambient and plume aerosols is thus unlikely to constitute a general highly efficient process for the removal of hydrogen chloride from the atmosphere.     

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.     

Water flume study of the enhancement of buoyant rise in pairs of merging plumes

When multiple stacks are grouped or ganged together at a site, the effluent plumes are often observed to merge downwind, forming a single buoyant plume whose rate of rise is enhanced relative to the rise of the plumes individually. The magnitude of this rise enhancement depends on many factors, and the few available models for rise enhancement do not always agree with one another. In the present study the rise behaviour of pairs of merging, buoyant plumes was studied by physical modelling in a water flume at 1:500 scale. The experiments were conducted at several stack separation distances and various exit velocity ratios for stack pairs aligned with, or perpendicular to, the ambient flow. Limited experiments were also done with the stacks aligned at other angles to the flow. The stack releases were made buoyant by heating the source water, and the resulting plumes were measured with an array of sensitive temperature probes. From these measurements it was possible to determine the plume structure and rise rates. For small stack separations when the stacks are aligned with the ambient flow, the experimental results show that the enhanced rise is close to, and sometimes above, the maximum theoretical rise enhancement factor of 2^1/3. For the perpendicular orientation there is little or no rise enhancement. The rise enhancement for other stack orientations is somewhere between these two extremes. A plausible physical explanation for the observed behaviour is given, based on initial momentum shielding and line vortex dynamics in the merging plumes.     

A new Lagrangian stochastic model for the gravity-slumping/spreading motion of a dense gas

A new dispersion model for dense gas which is released into the atmosphere on the flat terrain is constructed within the Lagrangian framework. Using the hydrostatic assumption for pressure distribution within cloud due to density variation, slumping motion is successfully incorporated into the Lagrangian model with entrainment effect naturally considered. Turbulence suppression due to stable stratification within cloud is also taken into consideration in the model formulation. Various results including time variant and maximum concentration predictions by the proposed model are compared with the available measured data in the experiment conducted in Thorney Island in 1984 with good agreement.     

工业点源大气污染扩散空间信息系统

开发了一个基于高斯扩散的大气污染扩散空间信息系统,用于模拟工业点源污染对区域大气质量的影响。该工业点源污染模型包括工业点源数据库、扩散参数、气象条件和大气质量评价4个主要数据库。用该模型计算上海市主要工业区的SO2排放,结果表明,该模型为模拟SO2污染扩散提供了一个有效便捷的方法。     

Measurements and simulations of the visual effects of particulate plumes

In this paper we present field measurements of the visual effects of particulate plumes from two power plants and a copper smelter. The measurements were conducted at downwind distances ranging from 7 to 34 km and for sun-observer angles ranging from 40 to 160°. The visual effects of the power plant plumes were relatively small due to atmospheric dispersion (Kincaid power plant plume in February 1981) or hazy background (Labadie power plant plume in August 1981). The plume from the San Manuel smelter was more visible because of the clean environment. The measurements of plume contrasts range from − 0.15 to + 0.15. Further development of the EPA plume visibility model to improve the treatment of multiple scattering of light and incorporate light absorption by carbonaceous aerosols is described. Teleradiometer measurements and model simulations are in reasonable agreement for cases in which experimental uncertainties are small. The model appears to underpredict forward scattering of light by plume particles.     

Measurement of atmospheric pollutants associated with oil and natural gas exploration and production activity in Pennsylvania’s Allegheny National Forest

Oil and natural gas exploration and production (E&P) activities generate emissions from diesel engines, compressor stations, condensate tanks, leaks and venting of natural gas, construction of well pads, and well access roads that can negatively impact air quality on both local and regional scales. A mobile, autonomous air quality monitoring laboratory was constructed to collect measurements of ambient concentrations of pollutants associated with oil and natural gas E&P activities. This air-monitoring laboratory was deployed to the Allegheny National Forest (ANF) in northwestern Pennsylvania for a campaign that resulted in the collection of approximately 7 months of data split between three monitoring locations between July 2010 and June 2011. The three monitoring locations were the Kane Experimental Forest (KEF) area in Elk County, which is downwind of the Sackett oilfield; the Bradford Ranger Station (BRS) in McKean County, which is downwind of a large area of historic oil and gas productivity; and the U.S. Forest Service Hearts Content campground (HC) in Warren County, which is in an area relatively unimpacted by oil and gas development and which therefore yielded background pollutant concentrations in the ANF. Concentrations of criteria pollutants ozone and NO₂ did not vary significantly from site to site; averages were below National Ambient Air Quality Standards. Concentrations of volatile organic compounds (VOCs) associated with oil and natural gas (ethane, propane, butane, pentane) were highly correlated. Applying the conditional probability function (CPF) to the ethane data yielded most probable directions of the sources that were coincident with known location of existing wells and activity. Differences between the two impacted and one background site were difficult to discern, suggesting the that the monitoring laboratory was a great enough distance downwind of active areas to allow for sufficient dispersion with background air such that the localized plumes were not detected.

ImplicationsMonitoring of pollutants associated with oil and natural gas exploration and production activity at three sites within the Allegheny National Forest (ANF) showed only slight site-to-site differences even with one site far removed from these activities. However, the impact was evident not in detection of localized plumes but in regional elevated ethane concentrations, as ethane can be considered a tracer species for oil and natural gas activity. The data presented serve as baseline conditions for evaluation of impacts from future development of Marcellus or Utica shale gas reserves.     

Effects of sorption on transverse mixing in transient flows

Transverse mixing has been identified as a potentially limiting factor for natural attenuation of plumes originating from continuously emitting sources. Under steady-state flow conditions, dispersion is the only process leading to lateral mixing. This process is very slow and cannot explain the lateral spread of plumes observed in the field. When the flow direction fluctuates with time, transverse dispersion is slightly enhanced, but not very pronounced. Under these flow conditions, however, sorption can contribute to mixing into the mean transverse direction. If the reacting compounds differ in their strength of sorption, chromatographic mixing and separation alternate in time-periodic flows. For instantaneous sorption, the plumes may overlap within a stripe of fixed width. In contrast to sorption in local equilibrium, kinetic sorption contributes to mixing also for compounds with identical sorption strength. I derive an analytical expression for the equivalent transverse dispersion coefficient of a kinetically sorbing compound in a spatially uniform flow field undergoing sinusoidal fluctuations in time. This expression may be used for reactive transport calculations in an equivalent domain with constant flow. The effects are the strongest for compounds with a dimensionless partitioning coefficient of about unity, slow sorption kinetics, and slowly fluctuating velocities. For realistic parameters, kinetic sorption contributes to transverse mixing in the same range as heterogeneity.