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.     

Application of the k–ε turbulence model to the high Reynolds number skimming flow field of an urban street canyon

A two-dimensional, steady, k–ε turbulence model was used to investigate the high Reynolds number skimming flow field of an urban street canyon. We describe the critical canyon width-to-height ratios that distinguish a cascade of vortex patterns that form in an urban street canyon. Details of the flow field are reported that includes the structure of the mean flow field, turbulent kinetic energy, turbulent length scale, turbulent eddy viscosity, and Reynolds stress for three typical different aspect ratios, W/H, of a street canyon. The consequences of vortex layering on vertical transport are explored.     

Impact of building facades and ground heating on wind flow and pollutant transport in street canyons

This paper investigates the impacts of building facades and ground heating on the wind flow and pollutant transport in street canyons using the computational fluid dynamic (CFD) technique. Street canyons of H/W (H representing the building height and W the street width) varied from 0.1 to 2, which covered the basic flow regimes of skimming flow (H/W=1 or 2), wake interference flow (H/W=0.5), and isolated roughness flow (H/W=0.1), were examined in a series of sensitivity tests. Heating that occurred on different surfaces, including ground surface and building façades, posed considerable effects on the street canyon wind flow and pollutant transport compared with those under isothermal conditions. The CFD results showed that the mechanically induced wind flow and pollutant transport were complicated by the buoyancy under temperature stratification. Individual street canyons of different H/W and surface-heating scenarios exhibited their unique wind flow structure and pollutant transport behaviors. Two counter-rotating vortices were calculated in the street canyons of H/W=1, in which the zone of higher pollutant concentration under isothermal conditions was switched from the leeward side to the windward side. In the street canyon of H/W=2, the recirculating wind pattern was perturbed by surface heating that led to the development of either one primary vortex or three closely coupled vortices. Because of the complicated wind structure, the zones of higher pollutant concentration located either on the leeward or windward ground level were subjected to the surface-heating scenarios. Only two vortices were developed inside the street canyon of H/W=0.5. The large primary vortex, centered inside the street canyon, extended above the roof level of the street canyon. Meanwhile, a small secondary vortex was found at the ground-level windward corner whose size results as a function of surface-heating configurations. Finally, in the street canyon of H/W=0.1, an isolated clockwise-rotating vortex was developed beside the leeward building while the wind in the windward side blew in the prevailing wind direction. As a result, air pollutant emitted at the street centerline was unlikely to be carried into the leeward vortex. Instead, it was dispersed rapidly on the windward side before being removed from the street canyon.     

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.     

Strategic guidelines for street canyon geometry to achieve sustainable street air quality

This paper is concerned with the motion of air within the urban street canyon and is directed towards a deeper understanding of pollutant dispersion with respect to various simple canyon geometries and source positions. Taking into account the present days typical urban configurations, three principal flow regimes “isolated roughness flow”, “skimming flow” and “wake interference flow” (Boundary Layer Climates, 2nd edition, Methuen, London) and their corresponding pollutant dispersion characteristics are studied for various canopies aspect ratios, namely relative height (h₂/h₁), canyon height to width ratio (h/w) and canyon length to height ratio (l/h). A field-size canyon has been analyzed through numerical simulations using the standard k-ε turbulence closure model. It is found that the pollutant transport and diffusion is strongly dependent upon the type of flow regime inside the canyon and exchange between canyon and the above roof air. Some rules of thumbs have been established to get urban canyon geometries for efficient dispersion of pollutants.     

On the escape of pollutants from urban street canyons

Pollutant transport from urban street canyons is numerically investigated using a two-dimensional flow and dispersion model. The ambient wind blows perpendicular to the street and passive pollutants are released at the street level. Results from the control experiment with a street aspect ratio of 1 show that at the roof level of the street canyon, the vertical turbulent flux of pollutants is upward everywhere and the vertical flux of pollutants by mean flow is upward or downward. The horizontally integrated vertical flux of pollutants by mean flow at the roof level of the street canyon is downward and its magnitude is much smaller than that by turbulent process. These results indicate that pollutants escape from the street canyon mainly by turbulent process and that the net effect of mean flow is to make some escaped pollutants reenter the street canyon. Further experiments with different inflow turbulence intensities, inflow wind speeds, and street aspect ratio confirm the findings from the control experiment. In the case of two isolated buildings, the horizontally integrated vertical flux of pollutants by mean flow is upward due to flow separation but the other main results are the same as those from the control experiment.     

Impacts of Upstream Building Width and Upwind Building Arrangements on Airflow and Pollutant Dispersion in a Street Canyon

A two-dimensional numerical model for evaluating the wind flow and pollutant dispersion within a street canyon was first developed using the FLUENT code, which was then validated against a wind tunnel experiment. Then, the effects of the upstream building width and upwind building arrangement on the airflow and pollutant dispersion inside an isolated street canyon were investigated numerically. The numerical results revealed that: (1) the in-canyon vortex center shifts downwards as the upstream building width increases; (2) the recirculation zone covers the entire upstream building roof for the cases when W/H = 0.5, 1.0, 1.5, and 2.0 (W is the upstream building width and H is the building height), whereas the flow reattaches the upstream building roof for the cases when W/H = 2.5 and 3.0; (3) when the upstream building width is shorter than the critical width W_C (= 2H), an increase in the upstream building width leads to an increase in the pollution level on the leeward wall of the canyon and a decrease in the roof-level concentrations at the upstream building; (4) when the upstream building width is longer than the critical width, the roof-level concentrations at the upstream building are negligibly small and the pollution level on the leeward wall of the canyon is almost unaffected by a further increase in the upstream building width; (5) when the buildings are placed upwind of the canyon, the flow attaches the upstream building roof and, therefore, almost none of the pollutants are distributed on the upstream building roof; and (6) the pollution levels inside the canyon and on the downstream building roof increase significantly with the number of upwind buildings.     

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.     

A measurement campaign in a street canyon in Helsinki and comparison of results with predictions of the OSPM model

In 1997, a measuring campaign was conducted in a street canyon (Runeberg St.) in Helsinki. Hourly mean concentrations of CO, NO_x, NO₂ and O₃ were measured at street and roof levels, the latter in order to determine the urban background concentrations. The relevant hourly meteorological parameters were measured at roof level; these included wind speed and direction, temperature and solar radiation. Hourly street level measurements and on-site electronic traffic counts were conducted throughout the whole of 1997; roof level measurements were conducted for approximately two months, from 3 March to 30 April in 1997. CO and NO_x emissions from traffic were computed using measured hourly traffic volumes and evaluated emission factors. The Operational Street Pollution Model (OSPM) was used to calculate the street concentrations and the results were compared with the measurements. The overall agreement between measured and predicted concentrations was good for CO and NO_x (fractional bias were −4.2 and +4.5%, respectively), but the model overpredicted the measured NO₂ concentrations (fractional bias was +22%). The agreement between the measured and predicted values was also analysed in terms of its dependence on wind speed and direction; the latter analysis was performed separately for two categories of wind velocity. The model qualitatively reproduces the observed behaviour very well. The database, which contains all measured and predicted data, is available for further testing of other street canyon dispersion models. The dataset contains a larger proportion of low wind speed cases, compared with other available street canyon measurement datasets.     

Strategic guidelines for street canyon geometry to achieve sustainable street air quality

This paper is concerned with the motion of air within the urban street canyon and is directed towards a deeper understanding of pollutant dispersion with respect to various simple canyon geometries and source positions. Taking into account the present days typical urban configurations, three principal flow regimes “isolated roughness flow”, “skimming flow” and “wake interference flow” (Boundary Layer Climates, 2nd edition, Methuen, London) and their corresponding pollutant dispersion characteristics are studied for various canopies aspect ratios, namely relative height (h₂/h₁), canyon height to width ratio (h/w) and canyon length to height ratio (l/h). A field-size canyon has been analyzed through numerical simulations using the standard k-ε turbulence closure model. It is found that the pollutant transport and diffusion is strongly dependent upon the type of flow regime inside the canyon and exchange between canyon and the above roof air. Some rules of thumbs have been established to get urban canyon geometries for efficient dispersion of pollutants.     

A field study of factors influencing the concentrations of a traffic-related pollutant in the vicinity of a complex urban junction

The paper describes a field study focused on the dispersion of a traffic-related pollutant within an area close to a busy intersection between two street canyons in Central London. Simultaneous measurements of airflow, traffic flow and carbon monoxide concentrations ([CO]) are used to explore the causes of spatial variability in [CO] over a full range of background wind directions. Depending on the roof-top wind direction, evidence of both flow channelling and recirculation regimes were identified from data collected within the main canyon and the intersection. However, at the intersection, the merging of channelled flows from the canyons increased the flow complexity and turbulence intensity. These features, coupled with the close proximity of nearby queuing traffic in several directions, led to the highest overall time-average measured [CO] occurring at the intersection. Within the main street canyon, the data supported the presence of a helical flow regime for oblique roof-top flows, leading to increased [CO] on the canyon leeward side. Predominant wind directions led to some locations having significantly higher diurnal average [CO] due to being mostly on the canyon leeward side during the study period. For all locations, small changes in the background wind direction could cause large changes in the in-street mean wind angle and local turbulence intensity, implying that dispersion mechanisms would be highly sensitive to small changes in above roof flows. During peak traffic flow periods, concentrations within parallel side streets were approximately four times lower than within the main canyon and intersection which has implications for controlling personal exposure. Overall, the results illustrate that pollutant concentrations can be highly spatially variable over even short distances within complex urban geometries, and that synoptic wind patterns, traffic queue location and building topologies all play a role in determining where pollutant hot spots occur.     

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.     

Analysis of air quality within a street canyon using statistical and dispersion modelling techniques

The dispersion model, ADMS-Urban, alongside the statistical modelling technique, generalized additive modelling, have been used to predict hourly NO_x and nitrogen dioxide (NO₂) concentrations at a busy street canyon location and the results compared with measurements. Generalized additive models (GAMs) were constructed for NO₂ and NO_x concentrations using input data required to run ADMS-Urban. Bivariate polar plots have been produced from the wind flow (speed and direction) and pollution data (measured and predicted concentrations) to provide further information regarding the complex wind-pollutant interactions in an urban street canyon. The predictions made with the GAMs show excellent agreement with measured concentrations at this location, reproducing both the magnitude of NO_x and NO₂ concentrations and also the wind speed-wind direction dependence of pollutant sources within the canyon. However, the predictions made with ADMS-Urban under-estimated the measured NO_x by 11% and NO₂ by 21% and there are clear differences in the bivariate polar plots. Several sensitivity tests were carried out with ADMS-Urban in an attempt to produce predictions in closer agreement to those measured at Gillygate. Increasing the primary NO₂ fraction in ADMS-Urban (from 10% to 20%) had a considerable effect on the predictions made with this model, increasing NO₂ predictions by ∼20%. However, the bivariate plots still showed major differences to those of the measurements. This work illustrates that generalized additive modelling is a useful tool for investigating complex wind-pollutant interactions within a street canyon.     

On the correlation of air and pollutant exchange for street canyons in combined wind-buoyancy-driven flow

The ventilation and pollutant transport in a two-dimensional (2D) street canyon of building-height-to-street-width (aspect) ratio h/b = 1 under different unstable stratifications were examined. To characterize the combined wind-buoyancy-driven flow and pollutant transport at different Richardson number Ri, a computational fluid dynamics (CFD) model based on the Reynolds-averaged Navier–Stokes (RANS) equations with the Renormalization Group (RNG) k ? ε turbulence model was adopted. Unlike the isothermal condition, a secondary recirculation is initiated at the ground-level windward corner of the street canyon once the unstable stratification is switched on (Ri < 0). It traps the ground-level pollutant leading to elevated pollutant concentration there. As Ri further decreases, the enlarging secondary recirculation enables direct pollutant removal from its core to the shear layer that offsets the ground-level pollutant accumulation. The ventilation and pollutant removal performance under different unstable stratifications are compared by the air (ACH) and pollutant (PCH) exchange rates, and pollutant retention time (τ). Both the mean and turbulent components of ACH are found to increase with decreasing Ri, suggesting that unstable stratification promotes ventilation in street canyons. Moreover, the CFD results agree well with our theoretical model that ACH² varies linearly with Ri. Turbulent transport originally dominates the pollutant removal under isothermal condition. However, progressive domination of pollutant removal by mean wind can be observed with decreasing stability (decreasing Ri from 0 to ?10.6). The critical value is estimated to be Ri = ?8, below which mean wind is the major pollutant removal carrier. Reduction in τ is also observed with decreasing Ri. Hence, in unstable stratification, pollutant resides shorter time in the street canyon compared with its isothermal counterpart, and the ventilation and pollutant removal are more favorable.     

A multi-approach monitoring of particulate matter, metals and PAHs in an urban street canyon

For the first time until now, the results from a prediction model (Atmospheric Dispersion Modelling System (ADMS)-Road) of pollutant dispersion in a street canyon were compared to the results obtained from biomonitors. In particular, the instrumental monitoring of particulate matter (PM10) and the biomonitoring of 14 polycyclic aromatic hydrocarbons (PAHs) and 11 metals by Quercus ilex leaves and Hypnum cupressiforme moss bags, acting as long- and short-term accumulators, respectively, were carried out. For both PAHs and metals, similar bioaccumulation trends were observed, with higher concentrations in biomonitors exposed at the leeward canyon side, affected by primary air vortex. The major pollutant accumulation at the leeward side was also predicted by the ADMS-Road model, on the basis of the prevailing wind direction that determines different exposure of the street canyon sides to pollutants emitted by vehicular traffic. A clear vertical (3, 6 and 9 m) distribution gradient of pollutants was not observed, so that both the model and biomonitoring results suggested that local air turbulences in the street canyon could contribute to uniform pollutant distribution at different heights.     

Removal of lead and zinc ions from aqueous solution using Amasya zeolites from Turkey

A 3-D Eulerian-Lagrangian approach to moving vehicles is presented that takes into account the traffic-induced flow rate and turbulence. The method is applied to pollutant dispersion in an individual street canyon and a system of two street canyons forming a perpendicular intersection. The approach is based on computational fluid dynamics (CFD) calculations using a Eulerian approach for continuous phase and a Lagrangian approach for moving vehicles. The wind speed was assigned values of 4, 7 and 12 m/s. One-way and two-way traffic with different traffic rates per lane is considered. In the case of the intersection, a longitudinal wind direction was assumed. Predictions show differences in the pollutant dispersion in the case of one-way and two-way traffic.     

Street canyon ventilation and atmospheric turbulence

Operational models for pollutant dispersion in urban areas require an estimate of the turbulent transfer between the street canyons and the overlying atmospheric flow. To date, the mechanisms that govern this process remain poorly understood. We have studied the mass exchange between a street canyon and the atmospheric flow above it by means of wind tunnel experiments. Fluid velocities were measured with a Particle Image Velocimetry system and passive scalar concentrations were measured using a Flame Ionisation Detector. The mass-transfer velocity between the canyon and the external flow has been estimated by measuring the cavity wash-out time. A two-box model, used to estimate the transfer velocity for varying dynamical conditions of the external flow, has been used to interpret the experimental data. This study sheds new light on the mechanisms which drive the ventilation of a street canyon and illustrates the influence of the external turbulence on the transfer process.     

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 Lagrangian stochastic model for predicting concentration fluctuations in urban areas

A combined Lagrangian stochastic model with a micromixing sub-model is used to estimate the fluctuating concentrations observed in two wind tunnel experiments. The Lagrangian stochastic model allows fluid trajectories to be simulated in the inhomogeneous flow, while the mixing model accounts for the dissipation of fluctuations using the interaction by exchange with the mean (IEM) mechanism. The model is first tested against the open terrain, wind tunnel data of Fackrell, J.E. and Robins, A.E. [1982. Concentration fluctuations and fluxes in plumes from point sources in a turbulent boundary layer. Journal of Fluid Mechanics 117, 1–26] and shows good agreement with the observed mean concentrations and fluctuation intensities. The model is then compared with the wind tunnel simulation of a two-dimensional street canyon by Pavageau, M. and Schatzmann, M. [1999. Wind tunnel measurements of concentration fluctuations in an urban street canyon. Atmospheric Environment 33, 3961–3971]. Despite the limitations of the k–ε turbulence scheme and the IEM mixing mechanism, the model reproduces the fluctuation intensity pattern within the canyon well. Overall, the comparison with both sets of wind tunnel experiments are encouraging, and the simplicity of the model means that predictions in a complex, three-dimensional geometry can be produced in a practicable amount of time.     

Operational air pollution modelling in the UK—Street canyon applications and challenges

Local air quality management requires the use of screening and advanced modelling tools that are able to predict roadside pollution levels under a variety of meteorological and traffic conditions. So far, more than 200 air pollution hotspots have been identified by local authorities in the UK, many of them associated with NO₂ and/or PM₁₀ exceedences in heavily trafficked urban streets that may be classified as street canyons or canyon intersections. This is due to the increased traffic-related emissions and reduced natural ventilation in such streets. Specialised dispersion models and empirical adjustment factors have been commonly used to account for the entrapment of pollutants in street canyons. However, most of the available operational tools have been validated using experimental datasets from relatively deep canyons (H/W⩾1) from continental Europe. The particular characteristics of low-rise street canyons (H/W<1), which are a typical feature of urban/sub-urban areas in the UK, have been rarely taken into account.The main objective of this study is to review current practice and evaluate three widely used regulatory dispersion models, WinOSPM, ADMS-Urban 2.0 and AEOLIUS Full. The model evaluation relied on two comprehensive datasets, which included CO, PM₁₀ and NO_x measurements, traffic information and relevant meteorological data from two busy street canyons in Birmingham and London for a 1-year period. The performance of the selected models was tested for different times of the day/days of the week and varying wind conditions. Furthermore, the ability of the models to reproduce roadside NO₂/NO_x concentration ratios using simplified chemistry schemes was evaluated for one of the sites. Finally, advantages and limitations of the current regulatory street canyon modelling practice in the UK, as well as needs for future research, have been identified and discussed.