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.     

Controlled Enrichment System for Experimental Fumigation of Plants in the Field with Sulfur Dioxide

A prototype Free-Air Controlled Enrichment (FACE) system was developed at Brookhaven National "Laboratory (BNL) for the experimental control of gas concentrations in an open field, without any form of enclosure. This FACE system consists of a toroidal plenum chamber, 15 m in diameter, with a series of 32 vertical vent pipes from which the release of a test gas can be controlled. Data on gas concentration at the center of the array and wind velocity are fed to a microprocessor where a proportional, integrative, differential (PID) algorithm is used to regulate a mass-flow controller. Data on wind direction are used to control the opening and closing of the vertical vent pipes to maintain an upwind release. In 72 hours of continuous operation during which wind velocity varied between 0.3 to 8 m sec^-1, the FACE prototype maintained 1 -min averaged concentrations within twenty percent of the 40 ppb set point 94 percent of the time.     

Comparison of a two-dimensional numerical dust transport model with experimental dust emissions from soil surfaces in a wind tunnel

Near-surface wind-tunnel fugitive dust concentration profiles arising from soil surfaces beds were compared to a finite difference numerical dust transport model. Comparisons of the type shown in this study were previously non-existent in the literature due to the lack of experimental wind-tunnel data for near-surface concentrations over a soil bed. However, in a previous study by the authors, near-surface steady-state concentration profiles were measured in order to obtain fugitive dust emission rates, thus allowing the comparison to models shown in this paper. The novel aspects of the current study include: comparison of concentration profiles of dust obtained experimentally in the wind tunnel with those calculated numerically; comparison of the calculated numerical fetch effect on dust emissions with that obtained in the wind tunnel; and comparison of the emission rates calculated numerically with those obtained experimentally in the wind tunnel. Initial comparisons with the model indicate good agreement implying that the physical mechanism of advection–diffusion is reasonably modeled with the choice of equations for the simple “steady-state” process near the surface. Furthermore, the numerical solutions presented in this paper provide a means to systematically explore the relative impact of varied surface boundary conditions upon the emission process and provide a potential link between wind-tunnel simulations and field scale models.     

Utilizing vortex behaviour to minimize drift

Abstract

The ULV spray emitted from a TBM flying in a cross wind was mapped by a scanning lidar system. The fate of the spray cloud for 2 min after release from the aircraft was followed as the material was transported downwind of the flight line. Vertical scans at 6 s intervals with 1 m^‐3 resolution provided detailed insight into the entrainment of the spray into the wing‐tip vortices and ultimate release to drift or deposit. Relative concentration, dosage and deposit profiles are presented for this cross‐wind case. Vortex lifetimes were found to be significantly different for the up‐wind and downwind vortices. The majority of the near field deposit was associated with the up‐wind vortex while the drift was linked to the down‐wind vortex.     

Numerical simulation of flow and dispersion around an isolated cubical building: The effect of the atmospheric stratification

The aim of this work is to investigate atmospheric flow and dispersion of contaminants in the vicinity of single buildings under different stability conditions. The mathematical model used is based on the solution of equations of conservation of mass, linear momentum and energy with the use of a non-standard κ–? turbulence model. The modifications proposed in the κ–? model are the inclusion of the Kato and Launder correction in the production of turbulent kinetic energy and the use of a modified wall function. Results are presented of numerical simulations of dispersion around a cubical obstacle, under neutral, stable and unstable atmospheric conditions. Experimental data from wind tunnel and field trials obtained by previous authors are used to validate the numerical results. The numerical simulation results show a reasonable level of agreement with field and wind tunnel concentration data. The deviation between model results and field experimental data is of the same order as the deviation between field and wind tunnel data.     

Wind-tunnel study of concentration fields in street canyons

The paper presents results from a case study of gaseous pollutant dispersion in street canyons. Tracer-gas experiments were performed in a neutrally stratified wind tunnel. Vehicle emissions were simulated as line sources. Concentration profiles along building walls were measured. A two-dimensional street canyon was considered as the reference case. The influence of systematic parameter variations on the concentration field is studied and discussed. Building dimensions, upwind building configuration, wind direction and roof geometry were found to be important parameters. Data sets from the study may be used for evaluation of numerical models and for expert estimates of air quality in the urban environment     

Characterization of near-road pollutant gradients using path-integrated optical remote sensing

Understanding motor vehicle emissions, near-roadway pollutant dispersion, and their potential impact to near-roadway populations is an area of growing environmental interest. As part of ongoing U.S. Environmental Protection Agency research in this area, a field study was conducted near Interstate 440 (I-440) in Raleigh, NC, in July and August of 2006. This paper presents a subset of measurements from the study focusing on nitric oxide (NO) concentrations near the roadway. Measurements of NO in this study were facilitated by the use of a novel path-integrated optical remote sensing technique called deep ultraviolet differential optical absorption spectroscopy (DUV-DOAS). This paper reviews the development and application of this measurement system. Time-resolved near-road NO concentrations are analyzed in conjunction with wind and traffic data to provide a picture of emissions and near-road dispersion for the study. Results show peak NO concentrations in the 150 ppb range during weekday morning rush hours with winds from the road accompanied by significantly lower afternoon and weekend concentrations. Traffic volume and wind direction are shown to be primary determinants of NO concentrations with turbulent diffusion and meandering accounting for significant near-road concentrations in off-wind conditions. The enhanced source capture performance of the open-path configuration allowed for robust comparisons of measured concentrations with a composite variable of traffic intensity coupled with wind transport (R2 = 0.84) as well as investigations on the influence of wind direction on NO dilution near the roadway. The benefits of path-integrated measurements for assessing line source impacts and evaluating models is presented. The advantages of NO as a tracer compound, compared with nitrogen dioxide, for investigations of mobile source emissions and initial dispersion under crosswind conditions are also discussed.     

Impact of using prognostic and objective wind fields on the photochemical modeling of Athens, Greece

A three-dimensional Eulerian photochemical model is used to follow the dynamics of ozone, NO_x, and CO over the Athens area, for 25 May 1990, the day considered in the APSIS project. A unique aspect of this work lies in the study of the impacts of the wind field preparation methods on the concentrations predicted by the model. Three sets of wind fields are developed. The first one used is derived from a prognostic meteorological model. The second one is calculated from available wind observations using objective: methods. For these two cases, a previous day is simulated, using the same conditions, to develop preconditioned initial conditions for the following day. For the third simulation, again two days are simulated, this time using the observed winds for each of the two days modeled. The predictions using the prognostically derived and the objective analysis wind fields are significantly different, particularly for the primary pollutants. Comparing predictions to the observations did not favor any particular method of wind field preparation. In this case, when using the prognostically derived field, the simulations are very sensitive to boundary conditions. In contrast, when using the wind fields constructed by objective methods, the simulations became most sensitive to emissions and initial conditions. This comes directly from the different residence times in the domain, which are governed by the wind speed.     

Estimating fugitive dust emission rates using an environmental boundary layer wind tunnel

Emissions from fugitive dust due to erosion of “natural” wind-blown surfaces are an increasingly important part of PM₁₀ (particulate matter with sizes of 10 μm aerodynamic diameter) emission inventories. These inventories are particularly important to State Implementation Plans (SIP), the plan required for each state to file with the Federal government indicating how they will comply with the Federal Clean Air Act (FCAA). However, techniques for determining the fugitive dust contribution to over all PM₁₀ emissions are still in their developmental stages. In the past, the methods have included field monitoring stations, specialized field studies and field wind-tunnel studies. The measurements made in this paper allow for systematic determination of PM₁₀ emission rates through the use of an environmental boundary layer wind tunnel in the laboratory. Near surface steady-state concentration profiles and velocity profiles are obtained in order to use a control volume approach to estimate emission rates. This methodology is applied to soils retrieved from the nation's single largest PM₁₀ source, Owens (dry) Lake in California, to estimate emission rates during active storm periods. The estimated emission rates are comparable to those obtained from field studies and lend to the validity of this method for determining fugitive dust emission rates.     

Assimilation of conventional and satellite wind observations in a mesoscale atmospheric model for studying atmospheric dispersion

A mesoscale atmospheric model PSU/NCAR MM5 is used to provide operational weather forecasts for a nuclear emergency response decision support system on the southeast coast of India. In this study the performance of the MM5 model with assimilation of conventional surface and upper-air observations along with satellite derived 2-d surface wind data from QuickSCAT sources is examined. Two numerical experiments with MM5 are conducted: one with static initialization using NCEP FNL data and second with dynamic initialization by assimilation of observations using four dimensional data assimilation (FDDA) analysis nudging for a pre-forecast period of 12 h. Dispersion simulations are conducted for a hypothetical source at Kalpakkam location with the HYSPLIT Lagrangian particle model using simulated wind field from the above experiments. The present paper brings out the differences in the atmospheric model predictions and the differences in dispersion model results from control and assimilation runs. An improvement is noted in the atmospheric fields from the assimilation experiment which has led to significant alteration in the trajectory positions, plume orientation and its distribution pattern. Sensitivity tests using different PBL and surface parameterizations indicated the simple first order closure schemes (Blackadar, MRF) coupled with the simple soil model have given better results for various atmospheric fields. The study illustrates the impact of the assimilation of the scatterometer wind and automated weather stations (AWS) observations on the meteorological model predictions and the dispersion results.     

Development of the ClearSky smoke dispersion forecast system for agricultural field burning in the Pacific Northwest

The post-harvest burning of agricultural fields is commonly used to dispose of crop residue and provide other desired services such as pest control. Despite careful regulation of burning, smoke plumes from field burning in the Pacific Northwest commonly degrade air quality, particularly for rural populations. In this paper, ClearSky, a numerical smoke dispersion forecast system for agricultural field burning that was developed to support smoke management in the Inland Pacific Northwest, is described. ClearSky began operation during the summer through fall burn season of 2002 and continues to the present. ClearSky utilizes Mesoscale Meteorological Model version 5 (MM5v3) forecasts from the University of Washington, data on agricultural fields, a web-based user interface for defining burn scenarios, the Lagrangian CALPUFF dispersion model and web-served animations of plume forecasts. The ClearSky system employs a unique hybrid source configuration, which treats the flaming portion of a field as a buoyant line source and the smoldering portion of the field as a buoyant area source. Limited field observations show that this hybrid approach yields reasonable plume rise estimates using source parameters derived from recent field burning emission field studies. The performance of this modeling system was evaluated for 2003 by comparing forecast meteorology against meteorological observations, and comparing model-predicted hourly averaged PM_2.5 concentrations against observations. Examples from this evaluation illustrate that while the ClearSky system can accurately predict PM_2.5 surface concentrations due to field burning, the overall model performance depends strongly on meteorological forecast error. Statistical evaluation of the meteorological forecast at seven surface stations indicates a strong relationship between topographical complexity near the station and absolute wind direction error with wind direction errors increasing from approximately 20° for sites in open areas to 70° or more for sites in very complex terrain. The analysis also showed some days with good forecast meteorology with absolute mean error in wind direction less than 30° when ClearSky correctly predicted PM_2.5 surface concentrations at receptors affected by field burns. On several other days with similar levels of wind direction error the model did not predict apparent plume impacts. In most of these cases, there were no reported burns in the vicinity of the monitor and, thus, it appeared that other, non-reported burns were responsible for the apparent plume impact at the monitoring site. These cases do not provide information on the performance of the model, but rather indicate that further work is needed to identify all burns and to improve burn reports in an accurate and timely manner. There were also a number of days with wind direction errors exceeding 70° when the forecast system did not correctly predict plume behavior.     

Evolution of stratospheric aerosols in the post-Pinatubo period measured by solar occultation

This paper presents particle size distributions of stratospheric aerosols derived from solar occultation data measured by the instrument ORA (an acronym for Occultation RAdiometer) during the period August 1992–May 1993. Starting from the UV-visible wavelength dependence of extinction coefficient profiles, an algorithm is developed that allows to retrieve the three parameters of an equivalent log-normal distribution and that makes use of vertical regularization. Comparison of retrieved mode radius and particle number density with existing data is found to be satisfactory. The evolution of the stratospheric aerosols is clearly influenced by sedimentation and coagulation as expected but also by vertical circulation. In a simple 1-D model, we derive the vertical wind profile and we interpret the temporal evolution of the particle mode radius.     

Intercomparison of nocturnal mixing height estimate methods for urban air pollution modelling

One of the most important meteorological input parameters for three-dimensional photochemical air pollution models is the mixing height h, which has a strong influence on the shape and intensity of the vertical diffusivity K_z and, as a consequence, on ground-level air concentrations of primary and secondary pollutants. A number of indirect algorithms for the estimate of h in nocturnal, stable conditions, when the mixing is dominated by mechanical turbulence, are reviewed and compared with mixing heights derived from wind (SODAR) and temperature (RASS) profiles measured in the Milan urban area during spring and summer 1996. Mixing heights derived from temperature soundings correlate positively with those derived from wind soundings only when a stable layer is superimposed to a quasi-adiabatic layer, while the correlation is very weak in the presence of a ground-based inversion. In general, indirect algorithms perform very poorly if compared with RASS-based estimates, and reasonably well if compared with SODAR-based estimates. Among the others, Benkley and Schulman (1979, Journal of Applied Meteorology 18, 772–780) method, which makes use of wind speed observed at 10 m height, and Nieuwstadt (1984, Boundary-Layer Meteorology 30, 31–55), which makes use of friction velocity and Monin–Obukhov length, give the best results.     

Dispersion of Traffic and Space Heating Emissions In Urban Settings

A long-term dispersion model is presented for traffic and space heating emissions in urban areas, allowing fast assessment of the spatial-averaged and center-maximum pollutant concentrations.

The assumption of study areas with circular shape and normal emissions density profiles is made for the purpose of streamlining model inputs with the inventory data normally available. In addition, the rather typical assumptions of Gaussian dispersion, narrow plume, flat or gently rolling terrain, homogeneous wind field and nonreactive pollutants are made. Values of σz from Briggs correlation are used with an initial value of 30 to account for building effects.

Meterological data inputs are reduced to six parameters, inventory data inputs to two, while computations are simplified to a degree that use of a digital computer is not required.

The model is well suited to yield separate assessments for individual types of sources and control measures, as well as to reveal sensitivities from parameters such as city size, or emission density levels and distribution patterns. Its predictions are virtually identical to those of the CDM-2 UNAPMAR model for study areas with circular shape and normal emissions density profiles, and as results do not appear overly sensitive to shape and distribution patterns, the model is believed to be valid for most urban areas.     

Characterization of Near-Road Pollutant Gradients Using Path-Integrated Optical Remote Sensing

Abstract

Understanding motor vehicle emissions, near-roadway pollutant dispersion, and their potential impact to near-roadway populations is an area of growing environmental interest. As part of ongoing U.S. Environmental Protection Agency research in this area, a field study was conducted near Interstate 440 (I-440) in Raleigh, NC, in July and August of 2006. This paper presents a subset of measurements from the study focusing on nitric oxide (NO) concentrations near the roadway. Measurements of NO in this study were facilitated by the use of a novel path-integrated optical remote sensing technique called deep ultraviolet differential optical absorption spectroscopy (DUV-DOAS). This paper reviews the development and application of this measurement system. Time-resolved near-road NO concentrations are analyzed in conjunction with wind and traffic data to provide a picture of emissions and near-road dispersion for the study. Results show peak NO concentrations in the 150 ppb range during weekday morning rush hours with winds from the road accompanied by significantly lower afternoon and weekend concentrations. Traffic volume and wind direction are shown to be primary determinants of NO concentrations with turbulent diffusion and meandering accounting for significant near-road concentrations in off-wind conditions. The enhanced source capture performance of the open-path configuration allowed for robust comparisons of measured concentrations with a composite variable of traffic intensity coupled with wind transport (R² = 0.84) as well as investigations on the influence of wind direction on NO dilution near the roadway. The benefits of path-integrated measurements for assessing line source impacts and evaluating models is presented. The advantages of NO as a tracer compound, compared with nitrogen dioxide, for investigations of mobile source emissions and initial dispersion under crosswind conditions are also discussed.     

Sowing simulation tests of a pneumatic drill equipped with systems aimed at reducing the emission of abrasion dust from maize dressed seed

The utilization of dressed seed for spring sowing is a widespread practice to control some pests with reduced doses of chemical products. However some insecticides employed in maize seed dressing, namely belonging to the neonicotinoid family and fipronil, have been claimed to play a role in the decline of honeybees (Apis mellifera L.). Pneumatic drills used in maize sowing are charged with contributing to the dispersion of the abrasion dust produced by dressed seeds, favoring the contamination of the honeybee habitat. Different devices similar to air deflectors have been introduced on pneumatic drills in order to reduce dust drift. During previous field tests carried out by the authors during recent years reduction of dust concentration both in the air and at soil surface has been shown as a consequence of their application. As field tests are affected by the variability of environmental parameters (namely wind speed and direction) the results are not always reliable, comparable and of a general validity. This paper refers to a sowing simulation test system in which pneumatic drills can be tested at a fixed point under controlled conditions of the main environmental parameters. In the test area, protected by external influences, artificial wind conditions are created by means of a fan. The drill, suitably placed in the test area, operates the seed distribution "sur place" by means of an electric engine connected to the drill's driving wheel. A 22.5 m long sampling area, leeward with respect to the drill position, has been identified. Along the sampling area a series of Petri dishes has been placed, with the aim of capturing the depositing dust and providing the concentration of the active ingredients (a.i.) at ground level. At the same time, three air samplers with PTFE diskette filters have been used for the detection of the a.i. The test system has been used for the test of a pneumatic drill, equipped with and without air deflectors, using maize seed dressed with four a. i. (imidacloprid, clothianidin, thiametoxam, fipronil). The results showed regularly decreasing of the concentrations as distance increased, both in the air and at ground level. Moreover, the difference determined by the adoption of the drift reducing device (air deflectors) resulted clear and it can be quantified at around 50 % of the a.i. amounts observed without deflectors. Finally, the paper proposes a data processing method that, from the values observed at fixed point, provides the theoretical a.i. concentration behavior that would occur in field, under the same conditions of wind speed and direction and working speed. The obtained results are coherent with previous field test.     

Wind tunnel study of air entrainment into two-dimensional dense gas plumes at the Chemical Hazards Research Center

Experiments in a neutrally stable wind tunnel boundary layer were made for two-dimensional (quasi-line) sources of carbon dioxide dispersing over two types of uniformly spaced (billboard) surface roughness elements. Velocity and concentration measurements were made with each surface roughness over a wide range of source Richardson number by varying carbon dioxide release rate and wind speed. Concentration measurements were made with a FID gas analyzer using an ethane tracer in the source gas, and velocity measurements were made with independent LDV and HWA systems. For each surface roughness, this paper describes the wind tunnel boundary layer and presents alongwind and vertical concentration profiles in the gas plume. Vertical velocity and concentration profiles were measured at selected downwind distances, and the profiles were integrated to confirm the consistency of the measurements with the mass of carbon dioxide released. The data are intended for development of improved vertical turbulent entrainment correlations for use in dense gas dispersion models applied to hazardous chemical consequence analyses.     

Turbulent Plume in a Turbulent Cross Flow: Comparison of Wind Tunnel Tests with Field Observations

Plume rise downwind of a large stationary gas turbine was measured in the field and the conditions were then scaled in the laboratory. For the laboratory, the plume exit conditions, wind velocity and temperature profiles, and wind direction were matched. It was found that for high temperature exhaust, the buoyancy is best matched by calculating a dimensionless density difference. With properly calculated buoyancy length scales, the plume trajectories were compared and were found to agree quite well. The probability distributions of the entrainment constant and the average values of the entrapment constant with downwind distance were compared. The field data showed about 15% greater plume rise. The median entrainment constant was about 10% greater for the lab test and the shape of the probability distribution matched very closely.     

KAPPA-G,A Non-Gaussian Plume Dispersion Model: Description and Evaluation Against Tracer Measurements

This paper presents a mixed methodology for the simulation of atmospheric disperson phenomena in which vertical diffusion is computed using an analytical solution of the K-theory equation, while horizontal diffusion is simulated by the Gaussian formula. This new formulation, while maintaining a simple analytical form for the concentration field, incorporates the effects of power-law vertical profiles of both wind speed and eddy diffusivity. The performance of this approach, which has been implemented into a full computer package (KAPPA-G), is evaluated by comparison with data from SF₆ tracer experiments.