The Modelling of Turbulence from Traffic in Urban Dispersion Models — Part I: Theoretical Considerations

The modelling of pollutant dispersion at the street scale in an urban environment requires the knowledge of turbulence generated by the traffic motion in streets. In this paper, a theoretical framework to estimate mechanical turbulence induced by traffic in street canyons at low wind speed conditions is established. The standard deviation of the velocity fluctuations is adopted as a measure of traffic-produced turbulence (TPT). Based on the balance between turbulent kinetic energy production and dissipation, three different parameterisations for TPT suitable for different traffic flow conditions are derived and discussed. These formulae rely on the calculations of constants that need to be estimated on the basis of experimental data. One such estimate has been made with the help of a wind tunnel data set corresponding to intermediate traffic densities, which is the most common regime, with interacting vehicle wakes.     

Evaluation of Reynolds stress, k-ε and RNG k-ε turbulence models in street canyon flows using various experimental datasets

In this study the Reynolds-averaged Navier-Stokes computational fluid dynamics methodology is used, which has proved to be a powerful tool for the simulations of the airflow and pollutant dispersion in the atmospheric environment. The interest is focused on the urban areas and more specifically on the street canyons, several types of which are examined in order to evaluate the performance of various turbulence models, including a Reynolds-stress model and variations of the k-ε model. The results of the two-dimensional simulations are compared with measurements from a diversity of independent street canyon experimental datasets, covering a wide range of aspect ratios, free stream velocities and roughnesses. This way more general and reliable conclusions can be reached about the applicability, accuracy and ease of use of each turbulence model. In this work, the renormalization group k-ε presented better results in most cases examined, while the Reynolds-stress model did not stand up for the expectations and also exhibited convergence problems.     

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.     

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.     

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.     

A Computational Fluid Dynamics Approach for Urban Area Transport and Dispersion Modeling

A simulation tool has been developed to model the wind fields, turbulence fields, and the dispersion of Chemical, Biological, Radiological and Nuclear (CBRN) substances in urban areas on the building to city blocks scale. A Computational Fluid Dynamics (CFD) approach has been taken that naturally accounts for critical flow and dispersion processes in urban areas, such as channeling, lofting, vertical mixing and turbulence, by solving the steady-state, Reynolds-Averaged Navier–Stokes (RANS) equations. Rapid generation of high quality cityscape volume meshes is attained by a unique voxel-based model generator that directly interfaces with common Geographic Information Systems (GIS) file formats. The flow and turbulence fields are obtained by solving the steady-state RANS equations using a collocated, pressure-based approach formulated for unstructured and polyhedral mesh elements. Turbulence modeling is based upon the Renormalization Group variant of the k–ε model (k–ε RNG). Neutrally buoyant simulations are made by prescribing velocity boundary condition profiles found by a power–law relationship, while turbulence quantities boundary conditions are defined by a prescribed mixing length in conjunction with the assumption of turbulence equilibrium. Dispersion fields are computed by solving an unsteady transport equation of a dilute gas, formulated in a Eulerian framework, using the velocity and turbulence fields found from the steady-state RANS solution. In this paper the model is explained and detailed comparisons of predicted to experimentally obtained velocity, turbulence and dispersion fields are made to neutrally stable wind tunnel and hydraulic flume experiments.     

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.

     

A CFD-based wind solver for an urban fast response transport and dispersion model

In many cities, ambient air quality is deteriorating leading to concerns about the health of city inhabitants. In urban areas with narrow streets surrounded by clusters of tall buildings, called street canyons, air pollution from traffic emissions and other sources may accumulate resulting in high pollutant concentrations. For various situations, including the evacuation of populated areas in the event of an accidental or deliberate release of chemical, biological and radiological agents, it is important that models should be developed that produce urban flow fields quickly. Various computational techniques have been used to calculate these flow fields, but these techniques are often computationally intensive. Most fast response models currently in use are at a disadvantage in these cases as they are unable to correlate highly heterogeneous urban structures with the diagnostic parameterizations on which they are based. In this paper, a novel variant of the popular projection method for solving the Navier–Stokes equations has been developed and applied to produce fast and reasonably accurate parallel computational fluid dynamics (CFD) solutions for flow in complex urban areas. This model, called QUIC-CFD represents an intermediate balance between fast (on the order of minutes for a several block problem) and reasonably accurate solutions. This paper details the solution procedure and validates this model for various simple and complex urban geometries.     

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.     

A building-based data capture and data mining technique for air quality assessment

Recently, a building-based air quality model system which can predict air quality in front of individual buildings along both sides of a road has been developed. Using the Macau Peninsula as a case study, this paper shows the advantages of building-based model system in data capture and data mining. Compared with the traditional grid-based model systems with input/output spatial resolutions of 1–2 km, the building-based approach can extract the street configuration and traffic data building by building and therefore, can capture the complex spatial variation of traffic emission, urban geometry, and air pollution. The non-homogeneous distribution of air pollution in the Macau Peninsula was modeled in a high-spatial resolution of 319 receptors·km^-2. The spatial relationship among air quality, traffic flow, and urban geometry in the historic urban area is investigated. The study shows that the building-based approach may open an innovative methodology in data mining of urban spatial data for environmental assessment. The results are particularly useful to urban planners when they need to consider the influences of urban form on street environment.     

Three-Dimensional Modelling of Concentration Fluctuations in Complicated Geometry

The strong fluctuating component in the measured concentration time series of a dispersing gaseous pollutant in the atmospheric boundary layer, and the hazard level associated to short-term concentration levels, demonstrate the necessity of calculating the magnitude of turbulent fluctuations of concentration using computational simulation models. Moreover the computation of concentration fluctuations in cases of dispersion in realistic situations, such as built-up areas or street canyons, is of special practical interest for hazard assessment purposes. In this paper, the formulation and evaluation of a model for concentration fluctuations, based on a transport equation, are presented. The model is applicable in cases of complex geometry. It is included in the framework of a computational code, developed for simulating the dispersion of buoyant pollutants over complex geometries. The experimental data used for the model evaluation concerned the dispersion of a passive gas in a street canyon between 4 identical rectangular buildings performed in a wind tunnel. The experimental concentration fluctuations data have been derived from measured high frequency concentrations. The concentration fluctuations model is evaluated by comparing the model's predictions with the observations in the form of scatter plots, quantile-quantile plots, contour plots and statistical indices as the fractional bias, the geometrical mean variance and the factor-of-two percentage. From the above comparisons it is concluded that the overall model performance in the present complex geometry case is satisfactory. The discrepancies between model predictions and observations are attributed to inaccuracies in prescribing the actual wind tunnel boundary conditions to the computational code.     

Multiscale Plume Transport from the Collapse of the World Trade Center on September 11, 2001

The collapse of the world trade center (WTC) produced enhanced levels of airborne contaminants in New York City and nearby areas on September 11, 2001 through December, 2001. This catastrophic event revealed the vulnerability of the urban environment, and the inability of many existing air monitoring systems to operate efficiently in a crisis. The contaminants released circulated within the street canyons, but were also lifted above the urban canopy and transported over large distances, reflecting the fact that pollutant transport affects multiple scales, from single buildings through city blocks to mesoscales. In this study, ground-and space-based observations were combined with numerical weather forecast fields to initialize fine-scale numerical simulations. The effort is aimed at reconstructing pollutant dispersion from the WTC in New York City to surrounding areas, to provide means for eventually evaluating its effect on population and environment. Atmospheric dynamics were calculated with the multi-grid Regional Atmospheric Modeling System (RAMS), covering scales from 250 m to 300 km and contaminant transport was studied using the Hybrid Particle and Concentration Transport (HYPACT) model that accepts RAMS meteorological output. The RAMS/HYPACT results were tested against PM2.5 observations from the roofs of public schools in New York City (NYC), Landsat images, and Multi-angle Imaging SpectroRadiometer (MISR) retrievals. Calculations accurately reproduced locations and timing of PM2.5 peak aerosol concentrations, as well as plume directionality. By comparing calculated and observed concentrations, the effective magnitude of the aerosol source was estimated. The simulated pollutant distributions are being used to characterize levels of human exposure and associated environmental health impacts.     

Comparison of atmospheric modelling systems simulating the flow,turbulence and dispersion at the microscale within obstacles

Three different modelling techniques to simulate the pollutant dispersion in the atmosphere at the microscale and in presence of obstacles are evaluated and compared. The Eulerian and Lagrangian approaches are discussed, using RAMS6.0 and MicroSpray models respectively. Both prognostic and diagnostic modelling systems are considered for the meteorology as input to the Lagrangian model, their differences and performances are investigated. An experiment from the Mock Urban Setting Test field campaign observed dataset, measured within an idealized urban roughness, is used as reference for the comparison. A case in neutral conditions was chosen among the available ones. The predicted mean flow, turbulence and concentration fields are analysed on the basis of the observed data. The performances of the different modelling approaches are compared and their specific characteristics are addressed. Given the same flow and turbulence input fields, the quality of the Lagrangian particle model is found to be overall comparable to the full-Eulerian approach. The diagnostic approach for the meteorology shows a worse agreement with observations than the prognostic approach but still providing, in a much shorter simulation time, fields that are suitable and reliable for driving the dispersion model.     

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

针对不同的城市街道峡谷结构形状,通过求解二维不可压缩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.     

An approach to simulate wind fields around an urban environment for wind energy application

This work proposes an approach to simulate wind flow fields around an urban environment with the aim of evaluating the potential impact of buildings on the general wind patterns and power production using the current generation of commercial wind turbines. The simulation process was performed with the aid of accessible computational tools that can potentially render the proposed procedure applicable in other cases of interest. The roughness of the urban environment was defined as the association of roughness map, topography, and an alternative process for obtaining the volumetry of buildings. A case study was conducted in a region located at the district of Boa Viagem (Recife-PE) for assessing the applicability of the approach. Scenarios were designed in order to simulate wind flow patterns and pre-identify sites that have suitable wind energy potential for electric power production by investigating the combination of wind speed magnitude and turbulence intensity. From the results obtained, it was possible to identify zones of potential wind sources that are not detected in classical wind atlas probably due to the influence of the built environment on local wind flow patterns.     

春季不同天气城市街头绿地内PM2.5质量浓度分布特征研究

于2012年春季,采用水平分布监测法监测了3种天气下距交通污染源不同距离的街头绿地中的PM2.5质量浓度,研究了不同天气下PM2.5的日变化、水平分布规律及绿地对PM2.5的净化效应,为城市街头绿地、城市公园建设及市民合理选择休闲锻炼时间和地点提供理论依据。结果表明：1）晴天和多云时,PM2.5质量浓度上午高于下午,7：00浓度最高,15：00最低;雨后阴天基本保持上升趋势。2）PM2.5日均质量浓度为晴天（61.67μg·m-3）〈多云（187.98μg·m-3）〈雨后阴天（291.48μg·m-3）。晴天除5：00和7：00外,其他时刻均达到国家二级标准,且13：00—15：00达到国家一类功能区空气质量要求;多云和雨后阴天PM2.5质量浓度分别超过国家二级浓度限值150.6%和288.6%。3）观测时段内,无论哪种天气,5：00、7：00、11：00和15：00绿地的净化功能较强,19：00净化功能均最差,所有监测点无一例外的表现为负效应。4）3种天气下,PM2.5质量浓度在距离道路10~25 m最高,绿地的净化效应最差,55 m外基本可以形成稳定的森林内环境。5）在南方高湿环境下,空气相对湿度是影响PM2.5质量浓度的主要因素,晴天和多云天气PM2.5浓度与相对湿度呈显著正相关,而雨后阴天二者呈负相关关系。6）在一定阈值内,街头绿地能够缓解PM2.5污染,为居民提供良好的休闲环境。从游憩时间来看,市民可以选择晴天进入街头绿地休闲,多云和雨后阴天尽量减少外出;从活动最佳地段来看,距离污染源55 m以上适宜休闲锻炼;从街头绿地的规划面积来看,半径以不小于55 m为宜。     

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.     

Turbulence Perturbations in the Neutrally Stratified Surface Layer due to the Interaction of a Katabatic Flow with a Steep Ridge

Turbulence measurements were performed in Antarctica, on the Nansen Ice Sheet, dominated by westerly katabatic winds. The measurements were taken at two sites aligned with the katabatic wind fall-line. The measuring stations were located in the middle of a wide, flat iced area at a distance of 14 km from the base of a sloping surface and at the top of a steep ridge (Inexpressible Island). The aim was to investigate the perturbation of turbulence close to the ground generated by the interaction of the flow with the ridge. We present an analysis comparing the data measured at the upstream unperturbed station with those at the top of the obstacle. Moments and spectra of velocity components have been calculated for almost steady periods. The topography and roughness change produce a combined effect on the flow acceleration (of the opposite sign) and on the turbulent stresses (of the same sign). Spectra of velocity components measured at the top of the ridge and scaled by unperturbed quantities evidence an increment of energy in the high frequency subrange with respect to the up-stream flow. Moreover, the horizontal velocity components display a shift in turbulence maximum towards higher frequencies. The vertical velocity spectrum exhibits an energy increment at low frequencies with respect to the upstream spectrum.     

A Numerical Study for Turbulent Flow and Thermal Influence over Inhomogenous Canopy of Roughness Elements

A large-eddy simulation with transitional structure function(TSF) subgrid model we previously proposed was performed to investigate the turbulent flow with thermal influence over an inhomogeneous canopy, which was represented as alternative large and small roughness elements. The aerodynamic and thermodynamic effects of the presence of a layer of large roughness elements were modelled by adding a drag term to the three-dimensional Navier–Stokes equations and a heat source/sink term to the scalar equation, respectively. The layer of small roughness elements was simply treated using the method as described in paper (Moeng 1984, J. Atmos Sci. 41, 2052–2062) for homogeneous rough surface. The horizontally averaged statistics such as mean vertical profiles of wind velocity, air temperature, et al., are in reasonable agreement with Gao et al.(1989, Boundary layer meteorol. 47, 349–377) field observation (homogeneous canopy). Not surprisingly, the calculated instantaneous velocity and temperature fields show that the roughness elements considerably changed the turbulent structure within the canopy. The adjustment of the mean vertical profiles of velocity and temperature was studied, which was found qualitatively comparable with Belcher et al. (2003, J Fluid Mech. 488, 369–398)’s theoretical results. The urban heat island(UHI) was investigated imposing heat source in the region of large roughness elements. An elevated inversion layer, a phenomenon often observed in the urban area (Sang et al., J Wind Eng. Ind. Aesodyn. 87, 243–258)’s was successfully simulated above the canopy. The cool island(CI) was also investigated imposing heat sink to simply model the evaporation of plant canopy. An inversion layer was found very stable and robust within the canopy.