A new Lagrangian particle model for the simulation of dense gas dispersion

A Lagrangian stochastic model (MicroSpray), able to simulate the airborne dispersion in complex terrain and in presence of obstacles, was modified to simulate the dispersion of dense gas clouds. This is accomplished by taking into account the following processes: negative buoyancy, gravity spreading and the particle's reflection at the bottom computational boundary. Elevated and ground level sources, continuous and instantaneous emissions, time varying sources, plumes with initial momentum (horizontal, vertical or oblique in any direction), plumes without initial momentum are considered. MicroSpray is part of the model system MSS, which also includes the diagnostic MicroSwift model for the reconstruction of the 3-D wind field in presence of obstacles and orography. To evaluate the MSS ability to simulate the dispersion of heavy gases, its simulation performances are compared in detail to two field experiments (Thorney Island and Kit Fox) and to a chlorine railway accident (Macdona). Then, a comprehensive analysis considering several experiments of the Modelers Data Archive is presented. The statistical analysis on the overall available data reveals that the performance of the new MicroSpray version for dense gas releases is generally reliable. For instance, the agreement between concentration predictions and observations is within a factor of two in the 72% up to 99% of the occurrences for the case studies considered. The values of other performance measures, such as correlation coefficient, geometric mean bias and geometric variance, mostly set in the ranges indicated as good-model performances in the specialized literature.     

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.     

Application of a model system for the study of transport and diffusion in complex terrain to the TRACT experiment

The transport and diffusion processes of a tracer gas released near the ground in the Rhine valley region, in Central Europe, during the 1992 TRACT field experiment, are simulated by a computational model system for complex terrain. This system (RMS) is composed of the prognostic mesoscale model RAMS, the Lagrangian stochastic dispersion model SPRAY and the interface code MIRS, which links RAMS to SPRAY. Three flow simulations were performed, with different initialisations and the one showing the best agreement with the measured flow was selected for the simulation of the TRACT tracer experiment. Tracer concentrations measured by an array of samplers at ground level and by an airplane aloft, are used to evaluate the 3-D concentration field simulated by the model system. The analysis of the simulation results generated by RMS shows that our model system very well reproduces the general behaviour of the contaminant plume, the temporal and spatial distribution of the concentration and the location of the concentration maxima.     

Analytical modelling of stable isotope fractionation of volatile organic compounds in the unsaturated zone

Analytical models were developed that simulate stable isotope ratios of volatile organic compounds (VOCs) near a point source contamination in the unsaturated zone. The models describe diffusive transport of VOCs, biodegradation and source ageing. The mass transport is governed by Fick's law for diffusion. The equation for reactive transport of VOCs in the soil gas phase was solved for different source geometries and for different boundary conditions. Model results were compared to experimental data from a one-dimensional laboratory column and a radial-symmetric field experiment. The comparison yielded a satisfying agreement. The model results clearly illustrate the significant isotope fractionation by gas phase diffusion under transient state conditions. This leads to an initial depletion of heavy isotopes with increasing distance from the source. The isotope evolution of the source is governed by the combined effects of isotope fractionation due to vaporisation, diffusion and biodegradation. The net effect can lead to an enrichment or depletion of the heavy isotope in the remaining organic phase, depending on the compound and element considered. Finally, the isotope evolution of molecules migrating away from the source and undergoing degradation is governed by a combined degradation and diffusion isotope effect. This suggests that, in the unsaturated zone, the interpretation of biodegradation of VOC based on isotopic data must always be based on a model combining gas phase diffusion and degradation.     

Modification of an Ambient Air Quality Model for Assessment of U. S. Naval Aviation Emittants

Previous studies have indicated that the EPA-HIWAY model significantly overestimates the pollutant concentrations for stable atmospheric conditions, especially under parallel wind-road orientation angles with low wind speed. This overestimation is due to the fact that the model's dispersion parameters do not properly account for the traffic-induced turbulence near roadways. In this paper, the Pasquill- Gifford dispersion curves used by the model are modified based on the recent studies that have quantified the nature of the trafficinduced turbulence and its influence on the pollutant dispersion in the near-field. The results show that the model performance is significantly improved when these new dispersion curves in conjunction with an aerodynamic drag factor, which in a rough way accounts for the change in the mean wind field due to the moving vehicles, are used in the HIWAY model.     

Experimental analysis of particle concentration heterogeneity in a ventilated scale chamber

The mixing processes of the aerosol particles from an outdoor environment in a ventilated scale chamber were experimentally studied. The particles were classified into five groups by size: 0.3–0.5 μm, 0.5–1.0 μm, 1.0–3.0 μm, 3.0–5.0 μm and 5.0–10.0 μm. The developing process for the concentration of each particle group was measured in different kinds of flow fields.The results show that the flow field configuration can effectively influence the dispersion time rate of the particles at certain positions. The increase in particle diameter can decrease the dispersion time rate. When the gas flow velocity is high, the particle dispersion time rate is independent of particle size; but when the gas flow velocity is low, particle size can significantly affect the particle dispersion time rate because the turbulent diffusion becomes important in the air and particle transport. The uniformity of the particle concentration for certain positions in steady state tends to be controlled by the inflow velocity, flow field configuration and the particle diameters.     

Turbulence and dispersion modeling near highways

Traffic-induced turbulence plays a dominant role in the dispersion of pollutants near highways. The formulations for velocity deficit and turbulence in vehicle wakes, developed from theoretical and physical modeling studies of Eskridge and his colleagues at US EPA about 20 years ago, are discussed. The vehicle wake parameterizations incorporated in ROADWAY-2, a near-highway pollutant dispersion model, and its evaluation results are described. The first field measurements of velocities and turbulence in the vehicle wake, using a towed array of 3-D sonic anemometers, are analyzed, and the results are presented and discussed. Specific recommendations are made for additional work in field measurements, laboratory studies, and mathematical model development and evaluation.     

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.     

Application of dispersion modeling to indoor gas release scenarios

Many complex models are available to study the dispersion of contaminants or ventilation effectiveness in indoor spaces. Because of the computationally complex numerical schemes employed, most of these models require mainframe computers or workstations. However, simple design tools or guidelines are needed, in addition to complicated models. A dispersion model based on the basic governing equations was developed and uses an analytical solution. Because the concentration is expressed by an analytical solution, the grid size and time steps are user definable. A computer program was used to obtain numerical results and to obtain release history from a thermodynamic source model. The model can be used to estimate three-dimensional spatial and temporal variations in concentrations resulting from transient gas releases in an enclosure. The model was used to study a gas release scenario from a pressurized cylinder into a large ventilated building, in this case, a transit parking and fueling facility.     

An atmospheric dispersion model for the environmental impact assessment of thermal power plants in Japan--a method for evaluating topographical effects

An atmospheric dispersion model was developed for the environmental impact assessment of thermal power plants in Japan, and a method for evaluating topographical effects using this model was proposed. The atmospheric dispersion model consists of an airflow model with a turbulence closure model based on the algebraic Reynolds stress model and a Lagrangian particle dispersion model (LPDM). The evaluation of the maximum concentration of air pollutants such as SO2, NOx, and suspended particulate matter is usually considered of primary importance for environmental impact assessment. Three indices were therefore estimated by the atmospheric dispersion model: the ratios (alpha and beta, respectively) of the maximum concentration and the distance of the point of the maximum concentration from the source over topography to the respective values over a flat plane, and the relative concentration distribution [gamma(x)] along the ground surface projection of the plume axis normalized by the maximum concentration over a flat plane. The atmospheric dispersion model was applied to the topography around a power plant with a maximum elevation of more than 1,000 m. The values of alpha and beta evaluated by the atmospheric dispersion model varied between 1 and 3 and between 1 and 0.4, respectively, depending on the topographical features. These results and the calculated distributions of y(x) were highly similar to the results of the wind tunnel experiment. Therefore, when the slope of a hill or mountain is similar to the topography considered in this study, it is possible to evaluate topographical effects on exhaust gas dispersion with reasonable accuracy using the atmospheric dispersion model as well as wind tunnel experiments.     

A numerical experiment of roadside diffusion under traffic-produced flow and turbulence

Roadside air pollution due to heavy traffic is one of the unsettled issues in the atmospheric environment in urban areas. As a practical application of a Computational Fluid Dynamics (CFD) model, a coupled mesoscale-CFD model was applied to the Ikegamicho area of Kawasaki City, Japan. For this study, the effects of traffic-produced flow and turbulence (TPFT) on the dispersion of the pollutants near the heavy traffic road were mainly investigated in an actual urban area. First, a series of preliminary CFD calculations was conducted for a road tunnel field experiment to obtain a fitting parameter for the traffic-produced flow. The calculation was then performed for 24 h in December 2005 around Ikegamicho, and the results were compared with the data at a roadside monitoring post in the area, located 10 m from the boundary of the ground road. In general, the effect of traffic-produced flow and turbulence was limited at the downstream side of the roads. The maximum concentration of NO_x was reduced and smoothed out along the traffic flow by the traffic-produced flow and turbulence on the road. The effects of traffic-produced turbulence on the dispersion of pollutants were greater than those of traffic-produced flow; however, the effects of traffic-produced flow were not negligible. The concentration of pollutants was not particularly dependent on the turbulent Schmidt number because most of the emission sources were introduced as volume sources in the present calculations, and the effect caused by differences in the material diffusion coefficient was not particularly significant at the outside of the road.     

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.     

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.     

Modeling the impacts of traffic emissions on air toxics concentrations near roadways

The dispersion formulation incorporated in the U.S. Environmental Protection Agency's AERMOD regulatory dispersion model is used to estimate the contribution of traffic-generated emissions of select VOCs – benzene, 1,3-butadiene, toluene – to ambient air concentrations at downwind receptors ranging from 10-m to 100-m from the edge of a major highway in Raleigh, North Carolina. The contributions are computed using the following steps: 1) Evaluate dispersion model estimates with 10-min averaged NO data measured at 7 m and 17 m from the edge of the road during a field study conducted in August, 2006; this step determines the uncertainty in model estimates. 2) Use dispersion model estimates and their uncertainties, determined in step 1, to construct pseudo-observations. 3) Fit pseudo-observations to actual observations of VOC concentrations measured during five periods of the field study. This provides estimates of the contributions of traffic emissions to the VOC concentrations at the receptors located from 10 m to 100 m from the road. In addition, it provides estimates of emission factors and background concentrations of the VOCs, which are supported by independent estimates from motor vehicle emissions models and regional air quality measurements. The results presented in the paper demonstrate the suitability of the formulation in AERMOD for estimating concentrations associated with mobile source emissions near roadways. This paper also presents an evaluation of the key emissions and dispersion modeling inputs necessary for conducting assessments of local-scale impacts from traffic emissions.     

城市生活垃圾焚烧处理过程中重金属迁移规律研究

通过采集厦门2个垃圾焚烧发电厂进厂垃圾、渗滤液、飞灰、炉渣和烟气样品,分析垃圾组成成分及各组分的重金属含量,结果表明,垃圾中以厨余、橡塑类和纸类为主,共占到垃圾干基的78．08％,重金属含量大小次序为Zn〉Cu〉Pb〉Cr〉Ni〉Cd〉Hg;垃圾渗滤液中除Zn外,其他金属的含量都较低,一厂渗滤液中重金属Ni、Zn迁移量较大,分别达到24．46％和8．52％。二厂渗滤液中重金属的迁移有相同的趋势,但含量相对较高。垃圾焚烧后其重金属主要分布在飞灰和废渣中,烟气中的含量非常少,不同金属的含量均有差别,这与金属的性质有很大关系;通过浸出毒性分析,飞灰中重金属酸溶态含量多,容易浸出,属于危险废物。同时,不同烟气处理工艺产生的飞灰的重金属浸出量有很大差别。     

An Atmospheric Dispersion Model for the Environmental Impact Assessment of Thermal Power Plants in Japan—A Method for Evaluating Topographical Effects

Abstract

An atmospheric dispersion model was developed for the environmental impact assessment of thermal power plants in Japan, and a method for evaluating topographical effects using this model was proposed. The atmospheric dispersion model consists of an airflow model with a turbulence closure model based on the algebraic Reynolds stress model and a Lagrangian particle dispersion model (LPDM). The evaluation of the maximum concentration of air pollutants such as SO₂, NO_x, and suspended particulate matter is usually considered of primary importance for environmental impact assessment. Three indices were therefore estimated by the atmospheric dispersion model: the ratios (α and β, respectively) of the maximum concentration and the distance of the point of the maximum concentration from the source over topography to the respective values over a flat plane, and the relative concentration distribution [γ(x)] along the ground surface projection of the plume axis normalized by the maximum concentration over a flat plane. The atmospheric dispersion model was applied to the topography around a power plant with a maximum elevation of more than 1000 m. The values of α and β evaluated by the atmospheric dispersion model varied between 1 and 3 and between 1 and 0.4, respectively, depending on the topographical features. These results and the calculated distributions of γ(x) were highly similar to the results of the wind tunnel experiment. Therefore, when the slope of a hill or mountain is similar to the topography considered in this study, it is possible to evaluate topographical effects on exhaust gas dispersion with reasonable accuracy using the atmospheric dispersion model as well as wind tunnel experiments.     

Derivation of exemption formulas for air quality regulatory applications

The regulatory agencies and the industries have the responsibility for assessing the environmental impact from the release of air pollutants, and for protecting environment and public health. The simple exemption formula is often used as a criterion for the purpose of screening air pollutants. That is, the exemption formula is used for air quality review and to determine whether a facility applying for and described in a new, modified, or revised air quality plan is exempted from further air quality review. The Bureau of Ocean Energy Management’s (BOEM) air quality regulations are used to regulate air emissions and air pollutants released from the oil and gas facilities in the Gulf of Mexico. If a facility is not exempt after completing the air quality review, a refined air quality modeling will be required to regulate the air pollutants. However, at present, the scientific basis for BOEM’s exemption formula is not available to the author. Therefore, the purpose of this paper is to provide the theoretical framework and justification for the use of BOEM’s exemption formula. In this paper, several exemption formulas have been derived from the Gaussian and non-Gaussian dispersion models; the Gaussian dispersion model is a special case of non-Gaussian dispersion model. The dispersion parameters obtained from the tracer experiments in the Gulf of Mexico are used in the dispersion models. In this paper, the dispersion parameters used in the dispersion models are also derived from the Monin-Obukhov similarity theory. In particular, it has been shown that the total amount of emissions from the facility for each air pollutant calculated using BOEM’s exemption formula is conservative.

Implications:?The operation of offshore oil and gas facilities under BOEM’s jurisdiction is required to comply with the BOEM’s regulations. BOEM’s air quality regulations are used to regulate air emissions and air pollutants released from the oil and gas facilities in the Gulf of Mexico. The exemption formulas have been used by BOEM and other regulatory agencies as a screening tool to regulate air emissions emitted from the oil and gas and other industries. Because of the BOEM’s regulatory responsibility, it is important to establish the scientific basis and provide the justification for the exemption formulas. The methodology developed here could also be adopted and used by other regulatory agencies.     

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.     

Numerical modeling on air quality in an urban environment with changes of the aspect ratio and wind direction

Due to heavy traffic emissions within an urban environment, air quality during the last decade becomes worse year by year and hazard to public health. In the present work, numerical modeling of flow and dispersion of gaseous emissions from vehicle exhaust in a street canyon were investigated under changes of the aspect ratio and wind direction. The three-dimensional flow and dispersion of gaseous pollutants were modeled using a computational fluid dynamics (CFD) model which was numerically solved using Reynolds-averaged Navier–Stokes (RANS) equations. The diffusion flow field in the atmospheric boundary layer within the street canyon was studied for different aspect ratios (W/H?=?1/2, 3/4, and 1) and wind directions (θ?=?90°, 112.5°, 135°, and 157.5°). The numerical models were validated against wind tunnel results to optimize the turbulence model. The numerical results agreed well with the wind tunnel results. The simulation demonstrated that the minimum concentration at the human respiration height within the street canyon was on the windward side for aspect ratios W/H?=?1/2 and 1 and wind directions θ?=?112.5°, 135°, and 157.5°. The pollutant concentration level decreases as the wind direction and aspect ratio increase. The wind velocity and turbulence intensity increase as the aspect ratio and wind direction increase.