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.     

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

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

CFD modelling of small particle dispersion: The influence of the turbulence kinetic energy in the atmospheric boundary layer

When considering the modelling of small particle dispersion in the lower part of the Atmospheric Boundary Layer (ABL) using Reynolds Averaged Navier Stokes simulations, the particle paths depend on the velocity profile and on the turbulence kinetic energy, from which the fluctuating velocity components are derived to predict turbulent dispersion. It is therefore important to correctly reproduce the ABL, both for the velocity profile and the turbulence kinetic energy profile.For RANS simulations with the standard k–ε model, Richards and Hoxey (1993. Appropriate boundary conditions for computational wind engineering models using the k–ε turbulence model. Journal of Wind Engineering and Industrial Aerodynamics 46–47, 145–153.) proposed a set of boundary conditions which result in horizontally homogeneous profiles. The drawback of this method is that it assumes a constant profile of turbulence kinetic energy, which is not always consistent with field or wind tunnel measurements. Therefore, a method was developed which allows the modelling of a horizontally homogeneous turbulence kinetic energy profile that is varying with height.By comparing simulations performed with the proposed method to simulations performed with the boundary conditions described by Richards and Hoxey (1993. Appropriate boundary conditions for computational wind engineering models using the k–ε turbulence model. Journal of Wind Engineering and Industrial Aerodynamics 46–47, 145–153.), the influence of the turbulence kinetic energy on the dispersion of small particles over flat terrain is quantified.     

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.     

Dilution of exhaust from a rooftop stack on a cubical building in an urban environment

The dispersion of a tracer gas emitted from a short stack on a cubical building was investigated using field and wind tunnel experiments. Air samples were obtained on the roof and the leeward side of the building. Dilution data were compared with estimates of minimum dilution (D_min) obtained with design formulas of Wilson/Chui/Lamb and Halitsky. The Halitsky model produced conservative predictions of D_min. Estimates of D_min obtained with the Wilson/Chui/Lamb model were in reasonable agreement with the field data. Wind tunnel dilution values were usually within a factor of two of the field data. The accuracy of the wind tunnel data improved as distance from the source increased. Wind tunnel dilution was found to be strongly dependent on the ratio of exhaust speed to wind speed, M, particularly for M values in the range of 2–4. This dependence is believed to be associated with the wind tunnel modelling of the stack exhaust, and is thus probably not a feature of the full-scale situation.     

Numerical simulation of dispersion around an isolated cubic building: Comparison of various types of k–ɛ models

Prediction accuracy of flow and dispersion around a cubic building with a flush vent located on its roof was examined using various k–? models, and numerical results were compared with wind-tunnel data. Four types of turbulence models, i.e., the standard k–? model, the RNG k–? model, the k–? model with Launder and Kato modification and the Realizable k–? model were compared in this study. The standard k–? model provided inadequate results for the concentration field, because it could not reproduce the basic flow structure, such as the reverse flow on the roof. However, revised k–? models provided concentrations in better agreement with the experimental data. The effect of an oblique wind angle and vent locations on the prediction accuracy was also investigated. It was confirmed that the prediction accuracy of the velocity field strongly affected that of the concentration field. The RNG model showed general agreement with the experiment, and was the best of the turbulence models tested. However, it becomes clear that the results for all CFD models show poor prediction accuracy of concentration distribution at the side and leeward surfaces of the building since they all underestimate the concentration diffusion on these regions. The concentrations predicted by all CFD models were less diffusive than those of the experiment.     

Evaluation of a Gaussian and a Lagrangian model against a roadside data set,with emphasis on low wind speed conditions

The evaluation of the high percentiles of concentration distributions is required by most national air quality guidelines, as well as the EU directives. However, it is problematic to compute such high percentiles in stable, low wind speed or calm conditions. This study utilizes the results of a previous measurement campaign near a major road at Elimäki in southern Finland in 1995, a campaign specifically designed for model evaluation purposes. In this study, numerical simulations were performed with a Gaussian finite line source dispersion model CAR-FMI and a Lagrangian dispersion model GRAL, and model predictions were compared with the field measurements. In comparison with corresponding results presented previously in the literature, the agreement of measured and predicted data sets was good for both models considered, as measured using various statistical parameters. For instance, considering all NO_x data (N=587), the so-called index of agreement values varied from 0.76 to 0.87 and from 0.81 to 1.00 for the CAR-FMI and GRAL models, respectively. The CAR-FMI model tends to slightly overestimate the NO_x concentrations (fractional bias FB=+14%), while the GRAL model has a tendency to underestimate NO_x concentrations (FB=−16%). The GRAL model provides special treatment to account for enhanced horizontal dispersion in low wind speed conditions; while such adjustments have not been included in the CAR-FMI model. This type of Lagrangian model therefore predicts lower concentrations, in conditions of low wind speeds and stable stratification, in comparison with a standard Lagrangian model. In low wind speed conditions the meandering of the flow can be quite significant, leading to enhanced horizontal dispersion. We also analyzed the difference between the model predictions and measured data in terms of the wind speed and direction. The performance of the CAR-FMI model deteriorated as the wind direction approached a direction parallel to the road, and for the lowest wind speeds. However, the performance of the GRAL model varied less with wind speed and direction; the model simulated better the cases of low wind speed and those with the wind nearly parallel to the road.     

Numerical simulation of scalar dispersion downstream of a square obstacle using gradient-transport type models

We present a numerical study of scalar transport released from a line source downstream of a square obstacle to investigate the capabilities and limitations of gradient-transport modeling in predicting atmospheric dispersion. The standard k–? and k–ω models and a Reynolds Stress Transport closure are employed and compared to predict the time-averaged turbulent flow field, while a standard gradient–diffusion model is initially adopted to relate the scalar flux to mean gradients of the concentration field. The analysis of two algebraic closures for turbulent scalar fluxes based on the generalized-gradient–diffusion hypothesis and its quadratic extension is also presented. In spite of the rather simple flow setup, where both the flow and the scalar fields can be assumed homogeneous in the spanwise direction, the analysis clarifies several critical issues concerning gradient-transport type models. We established the dominant role of predicted turbulent kinetic energy on scalar dispersion when a scalar diffusivity is employed, irrespectively of the Reynolds stress closure adopted for the averaged momentum equation. Moreover, the standard gradient–diffusion hypothesis failed to predict the streamwise component of the scalar flux, which is characterized by a counter-gradient-transport mechanism. Although the resulting contribution in the averaged scalar transport equation is small in the present flow configuration, this limitation can become severe for strongly inhomogeneous flows in the presence of point sources, where the spread of the scalar plume is essentially three-dimensional. The predictive capabilities of gradient-transport type modeling are found clearly improved using algebraic closures, which appear to represent a promising tool for predicting atmospheric dispersion in complex flows when unsteady transport mechanisms are not dominant.     

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.     

Modelling dispersion of traffic pollution in a deep street canyon: Application of CFD and operational models

In this study, numerical modelling of the flow and concentration fields has been undertaken for a deep street canyon in Naples (Italy), having aspect ratio (i.e. ratio of the building height H to the street width W) H/W = 5.7. Two different modelling techniques have been employed: computational fluid dynamics (CFD) and operational dispersion modelling. The CFD simulations have been carried out by using the RNG k–? turbulence model included in the commercial suite FLUENT, while operational modelling has been conducted by means of the WinOSPM model. Concentration fields obtained from model simulations have been compared with experimental data of CO concentrations measured at two vertical locations within the canyon. The CFD results are in good agreement with the experimental data, while poor agreement is observed for the WinOSPM results. This is because WinOSPM was originally developed and tested for street canyons with aspect ratio H/W ≌ 1. Large discrepancies in wind profiles simulated within the canyon are observed between CFD and OSPM models. Therefore, a modification of the wind profile within the canyon is introduced in WinOSPM for extending its applicability to deeper canyons, leading to an improved agreement between modelled and experimental data. Further development of the operational dispersion model is required in order to reproduce the distinct air circulation patterns within deep street canyons.     

A numerical study of interacting buoyant cooling-tower plumes

The compact design of mechanical cooling towers necessitates that the plumes are issued into the cross-wind in close proximity. An improved understanding of the interaction of adjacent plumes is therefore required for better design of such cooling towers, which may lead to a reduction in their environmental impact. This paper presents the results of a numerical investigation into the interaction of two adjacent plumes in a cross-flow. The numerical model simulates small-scale wind tunnel experiments of a cooling tower arrangement. The computations are performed for three-dimensional, turbulent, buoyant and interacting plumes, and for a single plume for comparison. Two double-source arrangements, namely, tandem and side-by-side, with respect to the oncoming atmospheric boundary layer are considered. A low Reynolds number k–ε turbulence model is used with two discretisation schemes, hybrid and QUICK, and the results are compared. Comparisons are also made with the experimental results. The results show that the interaction of side-by-side plumes is dominated by the interaction of the rotating vortex pairs within the plumes. A tandem source arrangement leads to early merging and efficient rise enhancement. Comparisons of the predicted results with experimental data show good agreement for the plume rise.     

A hybrid plume model for local-scale atmospheric dispersion

This paper presents a local-scale dispersion model, based on atmospheric boundary layer scaling theory. In the vicinity of the source, Gaussian equations are used in both the horizontal and vertical directions. After a specified transition distance, gradient transfer theory is applied in the vertical direction, while the horizontal dispersion is still assumed to be Gaussian. The dispersion parameters and eddy diffusivity are modelled in a form, which facilitates the use of a meteorological pre-processor. We present a novel model of the vertical eddy diffusivity (K_z), which is a continuous function of height in various atmospheric scaling regions. The model also includes a treatment of the dry deposition of gases and particulate matter. The accuracy of the numerical model was analysed by comparing the model predictions with two analytical solutions; the numerical deviations from these solutions were less than 2% for the computational regime. The model has been tested against the Kincaid experimental field data. The agreement of the predictions and the data is good on the average, although the internal variation of the predictions versus data scatter plot is substantial.     

Comparison of emission rate values for odour and odorous chemicals derived from two sampling devices

Field and laboratory measurements identified a complex relationship between odour emission rates provided by the US EPA dynamic emission chamber and the University of New South Wales wind tunnel. Using a range of model compounds in an aqueous odour source, we demonstrate that emission rates derived from the wind tunnel and flux chamber are a function of the solubility of the materials being emitted, the concentrations of the materials within the liquid; and the aerodynamic conditions within the device – either velocity in the wind tunnel, or flushing rate for the flux chamber. The ratio of wind tunnel to flux chamber odour emission rates (OU m^?2 s) ranged from about 60:1 to 112:1. The emission rates of the model odorants varied from about 40:1 to over 600:1.These results may provide, for the first time, a basis for the development of a model allowing an odour emission rate derived from either device to be used for odour dispersion modelling.     

Air ventilation impacts of the “wall effect” resulting from the alignment of high-rise buildings

The objective of this study is to investigate the air ventilation impacts of the so called “wall effect” caused by the alignment of high-rise buildings in complex building clusters. The research method employs the numerical algorithm of computational fluid dynamics (CFD – FLUENT) to simulate the steady-state wind field in a typical Hong Kong urban setting and investigate pollutant dispersion inside the street canyon utilizing a pollutant transport model. The model settings of validation study were accomplished by comparing the simulation wind field around a single building block to wind tunnel data. The results revealed that our model simulation is fairly close to the wind tunnel measurements. In this paper, a typical dense building distribution in Hong Kong with 2 incident wind directions (0° and 22.5°) is studied. Two performance indicators are used to quantify the air ventilation impacts, namely the velocity ratio (VR) and the retention time (T_r) of pollutants at the street level. The results indicated that the velocity ratio at 2 m above ground was reduced 40% and retention time of pollutants increased 80% inside the street canyon when high-rise buildings with 4 times height of the street canyon were aligned as a “wall” upstream. While this reduction of air ventilation was anticipated, the magnitude is significant and this result clearly has important implications for building and urban planning.     

Wind tunnel experiment of tracer gas diffusion within unstable boundary layer over coastal region

A wind tunnel experiment was carried out to simulate stack gas diffusion within an unstable atmospheric boundary layer over a coastal region. The wind tunnel floor, 4 m leeward of the entrance of the test section, was heated to 90°C over a length of 6 m in the streamwise direction, and wind tunnel experiments were performed under the flat plate condition with a prototype-to-model length scale ratio of 1200. Three similarity criteria of flow fields in the wind tunnel and in atmosphere, viz., bulk Richardson number, surface Reynolds number and the ratio of the Peclet number to the Richardson number, were considered in the wind tunnel experiment. Tracer gas was released along the coastline at a height of 10 cm, which corresponded to 120 m in height in atmosphere. The obtained wind tunnel experimental results of ground level concentration were compared with 30-min average values of the field experiments, viz., the data from the Tokai 82 field experiment. The maximum ground level concentration and its location were accurately simulated when there was close similarity between the wind tunnel and atmospheric flow conditions. The maximum concentration increased and occurred closer to the source when the level of convection was relatively stronger in atmosphere.     

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.     

An ensemble Kalman filter for atmospheric data assimilation: Application to wind tunnel data

In the previous work (Zheng et al., 2007, Zheng et al., 2009), a data assimilation method, based on ensemble Kalman filter, has been applied to a Monte Carlo Dispersion Model (MCDM). The results were encouraging when the method was tested by the twin experiment and a short-range field experiment. In this technical note, the measured data collected in a wind tunnel experiment have been assimilated into the Monte Carlo dispersion model. The uncertain parameters in the dispersion model, including source term, release height, turbulence intensity and wind direction have been considered. The 3D parameters, i.e. the turbulence intensity and wind direction, have been perturbed by 3D random fields. In order to find the factors which may influence the assimilation results, eight tests with different specifications were carried out. Two strategies of constructing the 3D perturbation field of wind direction were proposed, and the result shows that the two level strategy performs better than the one level strategy. It is also found that proper standard deviation and the correlation radius of the perturbation field play an important role for the data assimilation results.     

Computational Fluid Dynamics (CFD) Simulations on the Effect of Rough Surface on Atmospheric Turbulence Flow Above Hilly Terrain Shapes

The behavioral distribution of the atmospheric turbulence flow over the terrain with changes in a rough surface has become one of the most important topics of air pollution research, among such other topics as transportation and dispersion pollutants. In this study, a computational model on atmospheric turbulence flow over a terrain hill shaped with rough surface was investigated under neutral atmospheric conditions. The flow was assumed to be 2D and modeled using computational fluid dynamics (CFD) models, which were numerically solved using Reynolds-averaged Navier-Stokes equations. Rough surface conditions were modeled using a number of windbreak fences regularly spaced on the hill. The mean velocity and turbulent structures such as turbulence intensity and turbulent kinetic energy were investigated in the upwind and downwind regions over the hill, and the numerical models were validated against the wind-tunnel results to optimize the turbulence model. The computational results agreed well with the results obtained from the wind tunnel experiments. The computational results indicate that the mean velocity was observed to increase dramatically around the crest of the upwind slope of the hill. A thick internal boundary layer was observed with a fence on the crest and downwind region of the hill. The reversed flow and recirculation zone were formed in the wake region behind the hill. It was thus determined that turbulent kinetic energy decreases as the mean velocity increases.     

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.     

Experimental investigation of the residence of contaminants in the wake of an obstacle under different stability conditions

The detrainment behaviour of contaminants in the wake of an isolated building was investigated in the field under atmospheric stability conditions ranging from very stable to very unstable. The model building used was a 2 m cube and two orientations were investigated, with the cube either normal or at 45° to the wind. Tracer gas was first entrained into the wake from a source located a short distance upwind of the cube, the gas being released continuously for a limited period in order to fill the wake. Thereafter, the source was switched off, and the concentration (measured using several fast-response gas detectors located in the wake) was observed to decay in an exponential manner. This procedure was repeated in a total of 118 experiments to provide confidence in statistics. The residence time (T_d), which is defined as the time it takes for the concentration to decay to 1/e of its original value, was measured. The decay duration (t), which is the time it takes for the gas to become fully detrained from the wake, was found to be greater in stable atmospheric conditions, mainly due to the lower wind speeds and higher concentrations observed under these conditions. However, the non-dimensional residence time (τ) was found to be independent of atmospheric stability. The values of τ for a cube normal (τ=6.2) or at 45° to the flow (τ=9.5) are in very good agreement with values calculated using empirical formulae derived from wind tunnel experiments.