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.     

Evaluation of results of a numerical simulation of dispersion in an idealised urban area for emergency response modelling

WinMISKAM is evaluated from an emergency response perspective. Comparisons are made between ground level concentrations observed during selected Mock Urban Setting Test (MUST) field trials and predictions generated by the model. The model was driven by 5 min averaged on-site meteorological data, and used minimum grid spacing of 0.5 m in both the horizontal and vertical. The code was found to perform well, with 46% of all predictions (paired in time and space) and 83% of arc maxima predictions within a factor of two of observed concentrations. The model was found to perform better for neutral cases than stable cases with 27% of stable case predictions and 57% of neutral case predictions within a factor of two when compared in time and space.     

Three-dimensional modeling of air flow and pollutant dispersion in an urban street canyon with thermal effects

Effects of excess ground and building temperatures on airflow and dispersion of pollutants in an urban street canyon with an aspect ratio of 0.8 and a length-to-width ratio of 3 were investigated numerically. Three-dimensional governing equations of mass, momentum, energy, and species were modeled using the RNG k-epsilon turbulence model and Boussinesq approximation, which were solved using the finite volume method. Vehicle emissions were estimated from the measured traffic flow rates and modeled as banded line sources, with a street length and bandwidths equal to typical vehicle widths. Both measurements and simulations reveal that pollutant concentrations typically follow the traffic flow rate; they decline as the height increases and are higher on the leeward side than on the windward side. Three-dimensional simulations reveal that the vortex line, joining the centers of cross-sectional vortexes of the street canyon, meanders between street buildings and shifts toward the windward side when heating strength is increased. Thermal boundary layers are very thin. Entrainment of outside air increases, and pollutant concentration decreases with increasing heating condition. Also, traffic-produced turbulence enhances the turbulent kinetic energy and the mixing of temperature and admixtures in the canyon. Factors affecting the inaccuracy of the simulations are addressed.     

温度层结对城市街区流场及污染物扩散影响的数值模拟研究

采用标准k-ε湍流模型研究了温度层结对三维街区流场和污染物扩散的影响.结果表明,温度层结对街区流场和污染物均有一定影响.随着不稳定性的增加,气流涡旋中心向地面靠近.中性温度层结下,污染物随着街区内的涡旋先向背风侧迁移,然后主要随气流向下游迁移,很少向上游街区迁移.而不稳定温度层结下,上游街区污染物浓度也随之增加.根据污...     

Flow and dispersion in an urban cubical cavity

Flow and dispersion in an urban cubical cavity are numerically investigated using a Reynolds-averaged Navier–Stokes equations (RANS) model with the renormalization group (RNG) k–? turbulence closure model. The urban cubical cavity is surrounded by flank walls that are parallel to the streamwise direction, called end-walls, as well as upstream and downstream walls. A primary vortex and secondary vortices including end-wall vortices are formed in the cavity. Because of the end-wall drag effect, the averaged mean-flow kinetic energy in the cavity is smaller than that in an urban street canyon that is open in the along-canyon direction. A trajectory analysis shows that the end-wall vortices cause fluid particles to move in the spanwise direction, indicating that flow in the cavity is essentially three-dimensional. The iso-surfaces of the Okubo–Weiss criterion capture cavity vortices well. The pollutant concentration is high near the bottom of the upstream side in the case of continuous pollutant emission, whereas it is high near the center of the primary vortex in the case of instantaneous pollutant emission. To get some insight into the degree of pollutant escape from the cavity according to various meteorological factors, extensive numerical experiments with different ambient wind speeds and directions, inflow turbulence intensities, and cavity-bottom heating intensities are performed. For each experiment, we calculate the time constant, which is defined as the time taken for the pollutant concentration to decrease to e^?1 of its initial value. The time constant decreases substantially with increasing ambient wind speed, and tends to decrease with increasing inflow turbulence intensity and cavity-bottom heating intensity. The time constant increases as the ambient wind direction becomes oblique. High ambient wind speed is found to be the most crucial factor for ventilating the cavity, thus improving air quality in an urban cubical cavity.     

孤立与非孤立城市街道峡谷内污染物扩散

通过求解二维不可压N-S方程、k-ε方程及污染物对流扩散方程,模拟了孤立街道峡谷与非孤立街道峡谷内的流场及交通污染物浓度场.计算结果与风洞试验结果总体趋势一致.非孤立街道峡谷内污染物壁面浓度要大于孤立街道峡谷内的壁面浓度.通过计算街道峡谷建筑屋顶高度处的垂直方向污染物通量,说明了湍流扩散是污染物扩散出街道峡谷的主要原因,其污染物通量总为正,而平均流通量可以为负.非孤立街道峡谷由于平均流流动和湍流流动的总扩散通量减少,造成污染物在街道峡谷内集聚,从而理论上解释了非孤立街道峡谷与孤立街道峡谷污染扩散的差别.     

Characteristics of puff dispersion in an urban environment

Instantaneous releases of sulfur hexafluoride tracer were carried out as part of the Joint Urban 2003 field campaign in Oklahoma City. Data from 10 fast-response tracer samplers were used to examine the crosswind and along-wind spread of the tracer, the decay of tracer concentrations, and the retention of tracer within approximately 1 km of the release locations. The time variation of the median values of the tracer concentrations, normalized by the peak value observed at a given sampler, could be approximately described by an exponential decay with characteristic decay times on the order of 1–2 min. The longer times were found for early morning releases and the shorter times were associated with later morning or afternoon releases, suggesting that atmospheric stability or the depth of the mixed layer may affect puff dispersion even in urban environments. The median retention times required for 99% of the exposure to be realized at a given location were found to be correlated reasonably well with the median decay times. These characteristic time scales should be regarded as lower limits for concentration decay because the analysis excluded a number of anomalous cases in which the decaying concentrations exhibited an extended tail that indicated a very slow ventilation rate. The median values of the along-wind dispersion parameter σ_x grouped into downwind distance ranges can be described by a linear variation with distance with an initial “hold up” contribution due to building effects of about 30–45 m, but there are considerable variations about this relationship. Downwind 0.5–1 km from the release point the lateral puff dispersion (σ_y) was roughly 70% of the along-wind dispersion.     

Estimating micrometeorological inputs for modeling dispersion in urban areas during stable conditions

We examine the performance of three methods to estimate the surface friction velocity and the Monin–Obukhov (MO) length in stable conditions. Estimates from these methods are compared with measurements made at two urban sites: the Wilmington site located in the middle of an urban area, and the VTMX site located on a sloping, smooth area in Salt Lake City. The first method uses the mean wind at a single height (Single U or SU), the second uses the wind speed at a single level and the temperature difference between two levels (U delta T or UDT), and the third method uses two levels of wind speed and temperature (delta U delta T or DUDT). The performance of the SU and UDT methods in estimating u_* are comparable. The SU method yields better estimates of the MO length than the UDT method does. The DUDT method performs poorly in estimating both u_* and L. The major conclusions of this study are that (1) measurements of mean winds and temperatures at one or two levels at an urban location can provide adequate estimates of micrometeorological variables required in modeling dispersion in the stable boundary layer, and (2) methods based on using differences in temperatures and velocities between two levels can provide unreliable estimates of these variables because these differences can be overwhelmed by inevitable uncertainties in the measurement of mean variables.     

Air pollution dispersion within urban street canyons

A semi-empirical mathematical model, Urban Street Model (USM), is proposed to efficiently estimate the dispersion of vehicular air pollution in cities. This model describes urban building arrangements by combining building density, building heights and the permeability of building arrangements relative to wind flow. To estimate the level of air pollution in the city of Krasnoyarsk (in Eastern Siberia), the spatial distribution of pollutant concentrations off roadways is calculated using Markov's processes in USM. The USM-predicted numerical results were compared with field measurements and with results obtained from other frequently used models, CALINE-4 and OSPM. USM consistently yielded the best results. OSPM usually overestimated pollutant concentration values. CALINE-4 consistently underestimated these values. For OSPM, the maximum differences were 160% and for CALINE-4 about 400%. Permeability and building density are necessary parameters for accurately modeling urban air pollution and influencing regulatory requirements for building planning.     

Modeling reactive pollutant dispersion in an urban street canyon

Reactive pollutant dispersion in an urban street canyon with a street aspect ratio of one is numerically investigated using a computational fluid dynamics (CFD) model. The CFD model developed is a Reynolds-averaged Navier–Stokes equations (RANS) model with the renormalization group (RNG) k–ε turbulence model and includes transport equations for NO, NO₂, and O₃ with simple photochemistry. An area emission source of NO and NO₂ is considered in the presence of background O₃ and street bottom heating (ΔT=5 °C) with an ambient wind perpendicular to the along-canyon direction. A primary vortex is formed in the street canyon and the line connecting the centers of cross-sectional vortices meanders over time and in the canyon space. The cross-canyon-averaged temperature and reactive pollutant concentrations oscillate with a period of about 15 min. The averaged temperature is found to be in phase with NO and NO₂ concentrations but out of phase with O₃ concentration. The photostationary state defect is small in the street canyon except for near the roof level and the upper downwind region of the canyon and its local minimum is observed near the center of the primary vortex. The budget analysis of NO (NO₂) concentration shows that the magnitude of the advection or turbulent diffusion term is much larger (larger) than that of the chemical reaction term and that the advection term is largely balanced by the turbulent diffusion term. On the other hand, the budget analysis of O₃ concentration shows that the magnitude of the chemical reaction term is comparable to that of the advection or turbulent diffusion term. The inhomogeneous temperature distribution itself affects O₃ concentration to some extent due to the temperature-dependent photolysis rate and reaction rate constant.     

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.     

Pollutant dispersion in idealized city models with different urban morphologies

The mechanism of pollutant dispersion in idealized city models is investigated numerically by the introduction of a uniformly distributed pollutant source at street pedestrian level. We first study three short city forms with a single main street or two crossing streets, characterized by street length/street height ratios of L/H = 6 or 7 and a street height/street width ratio of H/W = 1, including a sharp-edged round city model, a smooth-edged round city model, and a sharp-edged square city model. For short city models with a single street and a parallel approaching wind, pollutant dilution mainly depends on the horizontal flow rate which decreases along the street. This decreasing rate is smallest for the smooth-edged round city model, which results in the lowest street concentrations. For city models with two crossing streets and the approaching wind parallel to the main street, the differences in overall city form result in different dispersion processes. For a sharp-edged round city model with two crossing streets, an approaching wind slightly non-parallel to the main street generates a lower pollutant concentration in the entire street volume. We also studied a sharp-edged round city model with one narrow street (L/H = 6; H/W = 6.7), finding that the uniformly distributed pollutants are transported from two street entries to the city centre, and are then removed out across the street roof. In contrast to the short city models we studied a single-street sharp-edged long rectangular city model (L/H = 21.7; H/W = 1) in which the horizontal flow rate remained nearly constant in a region far from the two entries. Within this region the turbulence across the street roof contributed more to the pollutant removal than vertical mean flows.     

Influence of trees on the dispersion of pollutants in an urban street canyon—Experimental investigation of the flow and concentration field

Flow field and concentration measurements have been performed in an idealized model of an urban street canyon with one row of trees arranged along the center axis. The model was set up in an atmospheric boundary layer wind tunnel and the approach flow was directed perpendicular to the street axis. A line source embedded in the bottom of the street was used to release tracer gas for the simulation of traffic exhaust emissions. Trees with spherical crowns were modeled and positioned inside the street canyon, varying crown diameter, crown permeability, trunk height and tree spacing. Traffic-induced turbulence was simulated by rotating belts with thin plates. Concentrations were measured at the facades of the street canyon. For small tree crowns, only little changes in concentration were measured, however, increasing crown diameters led to increasing concentrations at the leeward street canyon wall associated with a reduction of local concentrations at the windward wall. For some cases, a variation of trunk height led to a modification of the concentration pattern on the walls. Increasing the tree spacing resulted in a noticeable concentration decrease. When compared to the situation with standing (but emitting) traffic, the traffic-induced turbulence by two-way car movements always contributed to a more homogenous concentration field inside the street canyon yielding to reduced mean concentration levels.     

A comprehensive experimental databank for the models verification of urban car emission dispersion

A summary presentation is made of representative samples from a comprehensive experimental databank on car exhaust dispersion in urban street canyons. Physical modelling, under neutral stratification conditions, was used to provide visualisation, pollutant concentration and velocity measurements above and inside test canyons amidst surrounding urban roughness. The study extended to two different canyon aspects ratios, in combination with different roof configurations on the surrounding buildings. To serve as a reliable basis for validation and testing of urban pollution dispersion codes, special emphasis was placed in this work on data quality assurance.     

街道峡谷机动车尾气污染扩散模拟与影响因子分析研究

为预测和分析街道峡谷污染物浓度,研究了街道峡谷污染物浓度影响因子.利用重庆市交通干线街道峡谷两侧NOx浓度的监测数据,验证了街道峡谷机动车尾气污染扩散模型--OSPM模型.风速转换系数修正后的OSPM模型的模拟值与实测值的R达0.862 58;风场因子验证了风速转换系数修正后的OSPM模型能较好地模拟重庆市街道峡谷的污染物浓度,一定程度上能满足环境空气质量评价要求.同时,通过分析OSPM模型的影响因子,提出了控制街道峡谷机动车尾气污染状况的建议.     

Assessment of traffic-related air pollution in the urban streets before and during the 2008 Beijing Olympic Games traffic control period

In order to investigate the air quality and the abatement of traffic-related pollution during the 2008 Olympic Games, we select 12 avenues in the urban area of Beijing to calculate the concentrations of PM₁₀, CO, NO₂ and O₃ before and during the Olympic traffic controlling days, with the OSPM model.Through comparing the modeled results with the measurement results on a representative street, the OSPM model is validated as sufficient to predict the average concentrations of these pollutants at street level, and also reflects their daily variations well, i.e. CO presents the similar double peaks as the traffic flow, PM₁₀ concentration is influenced by other sources. Meanwhile, the model predicts O₃ to stay less during the daytime and ascend in the night, just opposite to NO₂, which reveals the impact of photochemical reactions. In addition, the predicted concentrations on the windward side often exceed the leeward side, indicating the impact of the special street shape, as well as the wind.The comparison between the predicted street concentrations before and during the Olympic traffic control period shows that the overall on-road air quality was improved effectively, due to the 32.3% traffic flow reduction. The concentrations of PM₁₀, CO and NO₂ have reduced from 142.6 μg m⁻³, 3.02 mg m⁻³ and 118.7 μg m⁻³ to 102.0 μg m⁻³, 2.43 mg m⁻³ and 104.1 μg m⁻³. However, the different pollutants show diverse changes after the traffic control. PM₁₀ decreases most, and the reduction effect focusing on the first half-day even clears the morning peak, whereas CO and NO₂ have even reductions to minify the daily fluctuations on the whole. Opposite to the other pollutants, ozone shows an increase of concentration. The average reduction rate of PM₁₀, CO, NO₂ and O₃ are respectively 28%, 19.3%, 12.3% and −25.2%. Furthermore, the streets in east, west, south and north areas present different air quality improvements, probably induced by the varied background pollution in different regions around Beijing, along with the impact of wind force. This finding suggests the pollution control in the surrounding regions, not only in the urban area.     

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.     

Coupled Monte Carlo simulation and Copula theory for uncertainty analysis of multiphase flow simulation models

Environmental Science and Pollution Research - Simulation-optimization techniques are effective in identifying an optimal remediation strategy. Simulation models with uncertainty, primarily in the...     

Multipoint moss passive samplers assessment of urban airborne polycyclic aromatic hydrocarbons: concentrations profile and distribution along Warsaw main streets

Polycyclic aromatic hydrocarbons (PAHs) distribution along 28 km of Warsaw main street have been surveyed in July 2000 using moss passive samplers as a simple and economic surrogate of direct air sampling. Altogether 74 samplers at 39 crossroads with traffic lights were placed on the lamp post approximately/=3.5 m above ground. PAHs levels determined in samplers are in range from 828 to 3573 ng/g moss dry weight. The spatial spread of pollution within this range is statistically close to normal distribution with mean value of 2332 ng/g. Variability within and between study areas are rationalized in terms of urban environmental factors. PAHs concentrations profiles across the town have appeared uniform. The dominant compounds are phenantrene, fluoranthene and pyrene. Their contribution is 49-68% of total PAHs burden.