Numerical simulations of the effect of building configurations and wind direction on fine particulate matters dispersion in a street canyon

To explore the effect of traffic emissions on air quality within street canyon, the wind flow and pollutant dispersion distribution in urban street canyons of different H/W, building gap and wind direction are studied and discussed by 3D computational fluid dynamics simulations. The largest PM_2.5 concentrations are 46.4, 37.5, 28.4 µg/m³ when x = ? 88, ? 19.3, ? 19.3 m in 1.5 m above the ground level and the ratio of H/W is 1:1, 1:2 and 2:1, respectively. The flow around the top of the building and clearance flow between the buildings in street canyon influence by different H/W, which affected the diffusion of fine particulate matters. The largest PM_2.5 concentrations are 88.1, 31.6 and 33.7 µg/m³ when x = 148.0, ? 92.3 and ? 186.7 m above the ground level of 1.5 m height and the building gap of 0, 20 and 40%, respectively. The air flows are cut by the clearance in the street canyons, and present the segmental characteristics. The largest PM_2.5 concentrations are 10.6, 11.2 and 16.0 µg/m³ when x = 165.3 m, x = 58.0 and 1.5 m above the ground level of 1.5 m height and wind direction of the parallel to the street, perpendicular to the street and southwest, respectively. Modelled PM_2.5 concentrations are basic agreement with measured PM_2.5 concentrations for southwest wind direction. These results can help analyze the difussion of PM_2.5 concentration in street canyons and urban planning.

     

On thermally forced flows in urban street canyons

During sunny days with periods of low synoptic wind, buoyancy forces can play a critical role on the air flow, and thus on the dispersion of pollutants in the built urban environments. Earlier studies provide evidence that when a surface inside an urban street canyon is at a higher temperature than that of local ambient air, buoyancy forces can modify the mechanically-induced circulation within the canyons (i.e., gaps between buildings). The aspect ratio of the urban canyon is a critical factor in the manifestation of the buoyancy parameter. In this paper, computational fluid dynamics simulations are performed on urban street canyons with six different aspect ratios, focusing on the special case where the leeward wall is at a greater temperature than local ambient air. A non-dimensional measure of the influence of buoyancy is used to predict demarcations between the flow regimes. Simulations are performed under a range of buoyancy conditions, including beyond those of previous studies. Observations from a field experiment and a wind tunnel experiment are used to validate the results.     

街道峡谷内汽车排放污染物浓度分布的观测与数值模拟

在上海某典型街道峡谷内按一定的空间布点,在一定时段内同时对各布点进行采样并做一氧化碳浓度分析,同时记录车辆种类、车流量、气象条件等,分析街道峡谷内污染物浓度的分布.运用风向频率加权(WDFW)方法,结合大气流动和污染物扩散的CFD模型进行数值模拟计算.结果表明,数值模拟结果和现场观测结果较吻合,建筑物低的一侧污染物浓度远高于建筑物高的一侧污染物浓度,两侧的污染物浓度随着高度的增加而降低.     

Large-eddy simulation of turbulent flow and dispersion over a complex urban street canyon

Turbulent flow and dispersion characteristics over a complex urban street canyon are investigated by large-eddy simulation using a modified version of the Fire Dynamics Simulator. Two kinds of subgrid scale (SGS) models, the constant coefficient Smagorinsky model and the Vreman model, are assessed. Turbulent statistics, particularly turbulent stresses and wake patterns, are compared between the two SGS models for three different wind directions. We found that while the role of the SGS model is small on average, the local or instantaneous contribution to total stress near the surface or edge of the buildings is not negligible. By yielding a smaller eddy viscosity near solid surfaces, the Vreman model appears to be more appropriate for the simulation of a flow in a complex urban street canyon. Depending on wind direction, wind fields, turbulence statistics, and dispersion patterns show very different characteristics. Particularly, tall buildings near the street canyon predominantly generate turbulence, leading to homogenization of the mean flow inside the street canyon. Furthermore, the release position of pollutants sensitively determines subsequent dispersion characteristics.     

Application of CFD modelling to air quality in Kuwait City

High-resolution computational fluid dynamics (CFD) simulations have been performed to assess the dispersion of air pollutants (CO₂) emanating from traffic in a busy street and in the vicinity of a complex configuration of buildings located in Salmiya, Kuwait City. New buildings are planned for this area, and the work here includes predictions for the dispersion of pollutants after the buildings’ completion. The CFD simulations are based on calculated CO₂ concentration levels for traffic counts taken on location in Salmiya with the existing configuration of buildings. As the computer code used in this work has been evaluated previously, it will be applied here to predict with confidence any potential air pollution problem areas on the addition of the new buildings. It was found for very light wind, that the proposed new buildings help reduce pollution in the vicinity of residential buildings within the configuration of buildings, but as the wind becomes moderate to strong, there was a tendency for the pollutant to get trapped in the residential area. Results are given for both exceptionally high ambient temperatures and very light wind, which are not often reported in studies found in the literature.     

城市街道内机动车排放污染物NO对流扩散的特征

采用现场观测和数值模拟的方法研究了城市街道内机动车排放污染物中的NO扩散特征。结果表明:城市街道中机动车排放污染物的对流扩散取决于屋顶风向和风速,随着建筑物顶部气流速度的增大,街道内同样位置的污染物浓度相对减小;当风向垂直于街道轴线时,街道内同样位置的污染物浓度最大;同时街道内机动车排放的污染物浓度与车流量成正比关系,即机动车流量越大污染物浓度越高。     

不同结构形状的街道峡谷内污染物扩散

针对不同的城市街道峡谷结构形状,通过求解二维不可压缩N-S方程和K-ε湍流模型方程及污染物对流扩散方程,数值模拟了街道峡谷内的流场及机动车排放污染物浓度场,从而说明了街道峡谷的结构是影响街道峡谷内污染气体扩散的主要因素之一。     

On the heating environment in street canyon

This study investigates the impact of building aspect ratio (building-height-to-street-canyon-width-ratio), wind speed and surface and air-temperature difference (Δθ_s−a) on the heating environment within street canyon. The Reynolds-averaged Navier-Stokes (RANS) and energy transport equations were solved with Renormalization group (RNG) theory version of k-e{\varepsilon} turbulence model. The validation process demonstrated that the model could be trusted for simulating air-temperature and velocity trends. The temperature and velocity patterns were discussed in idealized street canyons of different aspect ratios (0.5–2.0) with varying ambient wind speeds (0.5–1.5 m/s) and Δθ_s−a (2–8 K). Results show that air-temperatures are directly proportional to bulk Richardson number (R _b ) for all but ground heating situation. Conversely, air-temperatures increase significantly across the street canyon with a decrease in ambient wind speed; however, the impact of Δθ_s−a was negligible. Clearly, ambient wind speed decreases significantly as it passes over higher AR street canyons. Notably, air-temperatures were the highest when the windward wall was heated and the least during ground heating. Conversely, air-temperatures were lower along the windward side but higher within the street canyon when the windward wall was heated.     

The role of materials selection in the urban heat island effect in dry mid-latitude climates

This work investigates the role of materials selected for different urban surfaces (e.g. on building walls, roofs and pavements) in the intensity of the urban heat island (UHI) phenomenon. Three archetypal street-canyon geometries are considered, reflecting two-dimensional canyon arrays with frontal packing densities (λ_f) of 0.5, 0.25 and 0.125 under direct solar radiation and ground heating. The impact of radiative heat transfer in the urban environment is examined for each of the different built packing densities. A number of extreme heat scenarios were modelled in order to mimic conditions often found at low- to mid-latitudes dry climates. The investigation involved a suite of different computational fluid dynamics (CFD) simulations using the Reynolds-Averaged Navier–Stokes equations for mass and momentum coupled with the energy equation as well as using the standard k-ε turbulence model. Results indicate that a higher rate of ventilation within the street canyon is observed in areas with sparser built packing density. However, such higher ventilation rates were not necessarily found to be linked with lower temperatures within the canyon; this is because such sparser geometries are associated with higher heat transfer from the wider surfaces of road material under the condition of direct solar radiation and ground heating. Sparser canyon arrays corresponding to wider asphalt street roads in particular, have been found to yield substantially higher air temperatures. Additional simulations indicated that replacing asphalt road surfaces in streets with concrete roads (of different albedo or emissivity characteristics) can lead up to a ~5 °C reduction in the canyon air temperature in dry climates. It is finally concluded that an optimized selection of materials in the urban infrastructure design can lead to a more effective mitigation of the UHI phenomenon than the optimisation of the built packing density.     

Flows across high aspect ratio street canyons: Reynolds number independence revisited

The Reynolds number for flow in a street canyon, Re?=?U_refH/ν (where U_ref is a reference velocity, H the street canyon height, and ν the kinematic viscosity), cannot be matched between reduced-scale experiments and full-scale field measurements. This mismatch is often circumvented by satisfying the Re independence criterion, which states that above a critical Re (Re_c), the flow field remains invariant with Re. Re_c?=?11,000 is often adopted in reduced-scale experiments. In deep street canyons with height-to-width aspect ratio ≥?1.5, reduced-scale experiments have shown two recirculation vortices induced by the mean flows, but full-scale field measurements have observed only one vortex. We investigated this discrepancy by conducting water channel experiments with Re between 10⁴ and 10⁵ at three aspect ratios. The canyons with aspect ratio 1.0 have Re_c?=?11,000, the canyons with aspect ratio 1.5 have Re_c between 31,000 and 58,000, while the canyons with aspect ratio 2.0 have Re_c between 57,000 and 87,000. Therefore, the widely adopted Re_c?=?11,000 is not applicable for canyons with aspect ratio greater than 1.5. Our results also confirm that there is only one vortex in deep canyons at high Re. This single-vortex flow regime could change our fundamental understanding of deep canyons, which are often assumed to exhibit multiple-vortex flow regimes. Applications such as numerical model validation based on the multiple-vortex regime should be revisited. Our experimental data with Re up to 10⁵ could be used to validate numerical models at high Re.     

Evaluation of the QUIC-URB fast response urban wind model for a cubical building array and wide building street canyon

This paper describes the QUIC-URB fast response urban wind modeling tool and evaluates it against wind tunnel data for a 7 × 11 cubical building array and wide building street canyon. QUIC-URB is based on the Röckle diagnostic wind modeling strategy that rapidly produces spatially resolved wind fields in urban areas and can be used to drive urban dispersion models. Röckle-type models do not solve transport equations for momentum or energy; rather, they rely heavily on empirical parameterizations and mass conservation. In the model-experiment comparisons, we test two empirical building flow parameterizations within the QUIC-URB model: our implementation of the standard Röckle (SR) algorithms and a set of modified Röckle (MR) algorithms. The MR model attempts to build on the strengths of the SR model and introduces additional physically based, but simple parameterizations that significantly improve the results in most regions of the flow for both test cases. The MR model produces vortices in front of buildings, on rooftops and within street canyons that have velocities that compare much more favorably to the experimental results. We expect that these improvements in the wind field will result in improved dispersion calculations in built environments.     

Scalar mixing in an urban canyon

The scalar dynamics within a unit-aspect-ratio street canyon are studied using large-eddy simulation. The key processes of ventilation and mixing are analysed with the canyon-averaged concentration, mean tracer age and variance. The results are sensitive to the source location and can be classified according to the streamline geometry. The canyon-averaged concentrations for the corner vortices, vortex sea and central vortex do not converge to the same value at large times, though the mean decay rates do. The variance measured with respect to the canyon average shows two distinct decay regimes: the early regime reflects large-scale straining and enhanced diffusion across streamlines, while the late regime is associated with escape from the canyon, i.e., ventilation. Analytical predictions for the variance-decay or mixing time scales are verified for the early regime. It is argued that the presence of an open boundary at the roof level suppresses rapid mixing of the scalar field and is responsible for differences with respect to scalar dynamics within closed domains.

     

Experimental Study of the Criteria of Flow Laminarization in Two-Dimensional Dense Gas Plumes

Laminarization of flow in a two-dimensional dense gas plume was experimentally investigated in this study. The plume was created by releasing CO₂ through a ground-level line source into a simulated turbulent boundary layer over an aerodynamically rough surface in a meteorological wind tunnel. The bulk Richardson number (Ri_*), based on negative plume buoyancy, plume thickness, and friction velocity, was varied over a wide range so that the effects of stable stratification on plume laminarization could be observed. A variety of ambient wind speeds as well as three different sizes of roughness arrays were used so that possible effects of roughness Reynolds number (Re_*) on plume laminarization could also be identified. Both flow visualization methods and quantitative measurements of velocity and intermittency of turbulence were used to provide quantitative assessments of plume laminarization.Flow visualization provided an overall picture of how the plume was affected by the negative buoyancy. With increasing Ri_*, both the plume depth and the vertical mixing were significantly suppressed, while upstream propagation of the plume from the source was enhanced. The most important feature of the flow revealed by visualization was the laminarization of flow in the lower part of the plume, which appeared to be closely related to both Ri_* and Re_*.Measurements within the simulated dense gas plumes revealed the influence of the stable stratification on mean velocity and turbulence intensity profiles. Both the mean velocity and turbulence intensity were significantly reduced near the surface; and these reductions systematically depended on Ri_*. The roughness Reynolds number also had considerable influence on the mean flow and turbulence structure of the dense gas plumes.An intermittency analysis technique was developed and applied to the digitized instantaneous velocity signals. It not only confirmed the general flow picture within the dense plume indicated by the flow visualization, but also clearly demonstrated the changes of flow regime with variations in Ri_* and Re_*. Most importantly, based on this intermittency analysis, simple criteria for characterizing different flow regimes are formulated; these may be useful in predicting when plume laminarization might occur.     

Wind-induced pressure at a tunnel portal

In order to properly size the mechanical ventilation system of a tunnel, it is essential to estimate the wind-driven pressure difference that might rise between its two portals. In this respect, we explore here the pressure distribution over a tunnel portal under the influence of an incident atmospheric boundary layer and, in particular, its dependency on wind direction and on tunnel geometry. Reduced scale models of generic configurations of a tunnel portal are studied in an atmospheric wind tunnel. Pressure distributions over the front section of different open cavities are measured with surface taps, which allows us to infer the influence of the tunnel aspect ratio and wind direction on a pressure coefficient \(C_{P}\), defined as a spatially and time averaged non-dimensional pressure. Experiments reveal that the magnitude of the coefficient \(C_{P}\), as a function of the wind direction, is significantly influenced by the portal height-to-width ratio and almost insensitive to its length. The experimental data set is completed by hot-wire anemometry measurements providing vertical distribution of velocity statistics. The same configurations are simulated by numerically solving the Reynolds-averaged Navier–Stokes equations, adopting the standard \(k - \varepsilon\) turbulence model. Despite some discrepancies between numerical and experimental estimates of some flow parameters (namely the turbulent kinetic energy field), the numerical estimates of the pressure coefficients \(C_{P}\) show very good agreement with experimental data. The latter is also compared to the predictions of an analytical model, based on the estimate of a spatially averaged velocity within an infinitely long street canyon. The results of the model, which takes into account varying canyon aspect ratios, are in reasonable agreement with experimental data for all cases studied. Notably, its predictions are significantly better than those provided by the simple analytical relations usually adopted as a reference in tunnel ventilation studies.     

Upstream Turbulence Effect on Pollution Dispersion

Experiments have been carried out to investigate turbulence at and above roof-level in an urban environment, and to predict the behaviour of street pollution from experiments using dye dispersion, for different roughness conditions and bed geometries. The flow in the boundary layer above an idealised urban environment has been simulated in a laboratory water flume. Comparisons have been made for the same model street canyon with and without the presence of upstream roughness. In the tests reported here, model street canyons were aligned perpendicular to the flow direction, and velocity measurements made within and above the model street canyons using a laser Doppler velocimeter (LDV). Flow visualisation techniques have also been used to confirm the gross flow features from streak images. Turbulence generated from the upstream roughness has a significant effect on the turbulence production and dispersion behaviour of the dye simulating pollution in street canyons.     

Evaluation of CFD Simulation using RANS Turbulence Models for Building Effects on Pollutant Dispersion

CFD evaluations were performed to examine the applicability of the RANS methods in simulating pollutant dispersion near, within and over three typical building configurations: (1) an isolated building, (2) a building array and (3) an urban intersection. The CFD results are compared with values obtained from wind tunnel tests. In some situations major differences between the wind tunnel tests and the CFD results were observed. The main source of difference between the CFD and wind tunnel results was inadequate modelling of local flow patterns using the RANS turbulence models. Also inappropriate evaluation of high intermittent turbulent mixing in the RANS approach may lead to either over-prediction or under-prediction of the concentration level, by up to a factor of 10, depending on the case investigated.     

An assessment of the air pollution data from two monitoring stations in Kuwait

This article assesses the air pollution data from two monitoring stations in Kuwait. The measurements cover major pollutants, i.e., CO, CO₂, methanated and non-methanated hydrocarbons, NO _x , SO₂, O₃, and particulate matter (PM10). The data also includes meteorological parameters, i.e., solar intensity, temperature, wind speed, and wind direction, and has been collected over a period 4 years, from 2001 to 2004. Data analysis includes the assessment of annual hourly averages and 1-h maxima. Typical pollutant concentration trends, similar to those previously reported for Kuwait and for other locations around the world, are observed except for particulate matter measurements, which have higher values because of proximity to the desert. Emissions of nitrogen oxides show a consistent increase over the years. This is caused by the increase in the number of motor vehicles and the expansion in power generation and industrial activities. The data collected is a subset of the air quality criteria, as defined by the US EPA (United States Environmental Protection Agency).     

Effects of SO<Subscript>2</Subscript> contamination on rising CO<Subscript>2</Subscript> drops under high pressure

Although the flow dynamics of pure liquid drops in other liquids has been well researched, little attention has been paid to the impacts of impurities. Hence, most of research is not directly applicable to the real world. To address this gap, we conducted numerical experiments simulating the rise of pure and contaminated drops. It was selected to study liquid CO₂ drops contaminated with SO₂ under high pressure because such mixtures mimic potential scenarios in which drops may leak from carbon capture and storage (CCS) facilities or pipelines. First, numerical simulation experiments were performed to validate our method by comparing our results with previous research on pure drops. Second, the validated numerical approach was applied to simulations of contaminated drops to investigate how contaminants affect rising drops. The results show that the SO₂ contamination caused changes in deformation, breakup phenomena, rising velocities, surrounding flow fields and drag coefficients. Most importantly, the contamination resulted in the formation of smaller “child drops”; such breakup is not observed in pure CO₂ drops. The formation of child drops in turn affects the streamlines, patterns and areas of wakes behind the contaminated drops. The addition of contaminants also enhances the dissolution rate, which is affected by the contaminant concentration and by the flow dynamics of the rising drop. Our results would improve understanding the rise of impure CO₂ drops, such as drops potentially leaked by future CCS operations.     

Large-eddy simulation of low-level jet-like flow in a canopy

The computational method of Large-Eddy Simulations has been used to study the weak, neutrally stable drainage flow within tree canopies. The computational results show that a secondary velocity maximum that resembles a jet is formed within the canopy under the nocturnal flow conditions. This jet-like flow is important in the analysis and measurements of the net ecosystem-atmosphere exchange (NEE) for carbon dioxide (CO₂). A uniformly distributed, plane source was placed within the canopy in order to simulate the nocturnal production of CO₂. The NEE is calculated as the sum of the integration of the rate of change of the concentration of CO₂ over the computational domain, the vertical turbulent flux measured directly by eddy-covariance (EC) method, and the advection terms, which are not taken into account in the EC method. Numerical results of the velocity and concentration fields, within and above the canopy, are presented and their impact on the CO₂ transport is investigated in detail. The computational results show that 15–20% of NEE is drained out by the advection process under the canopy. The results also show that the turbulent fluctuations in the lateral direction are also significant and may result in 2–5% CO₂ transport.     

Transport and Diffusion of Ozone in the Nocturnal and Morning Planetary Boundary Layer of the Phoenix Valley

The evolution of ozone (O ₃) in the nocturnal and morning-transitional planetary boundary layer (PBL) of the Phoenix valley was measured as a part of the `Phoenix Sunrise Experiment 2001' of the U.S. Department of Energy conducted in June 2001. The goal of the field program was to study the transport, distribution and storage of ozone and its precursors in the urban boundary layer over a diurnal cycle. The ground level O <sub>3</sub> as well as mean meteorological variables and turbulence were measured over the entire period, and vertical profiling (using a tethered balloon) was made during the morning transition period. Approximately half of the observational days showed the usual diurnal cycle of high O <sub>3</sub> during the day and low O <sub>3</sub> at night, with nitrogen oxides (NO <sub>x</sub> = NO <sub>2</sub> + NO) showing an out of phase relationship with O <sub>3</sub>. The rest of the days were signified by an anomalous increase of O <sub>3</sub> in the late evening ( 2200 LST), concomitant with a sudden drop of temperature, an enhancement of wind speed and Reynolds stresses, a positive heat flux and a change of wind direction. NO <sub>x</sub> measurements indicated the simultaneous arrival of an `aged' air mass, which was corroborated by the wind predictions of a mesoscale numerical model. In all, the results indicate that the recirculation of O ₃ rich air masses is responsible for the said high-O ₃ events. Such air masses are produced during the transport of O ₃ precursors by the upslope flow toward mountainous suburbs during the day, and they return back to the city at night via downslope winds (i.e. mountain breeze). The corresponding flow patterns, and hence the high-O ₃ events, are determined by background meteorological conditions. The vertical profiling of O ₃ and flow variables during the morning transition points to a myriad of transport, mixing and chemical processes that determine the fate of tropospheric O ₃. How well such processes are incorporated and resolved in predictive O ₃ models should determine the accuracy of their predictions.