Large-Eddy Simulation of Pollutant Removal from a Three-Dimensional Street Canyon

Large-eddy simulations were conducted to investigate the mechanism of pollutant removal from a three-dimensional street canyon. Five block configurations with aspect ratios (building height to length) of 1, 2, 4, 8 and $\infty $ were used to create an urban-like array. A pollutant was released from a ground-level line source at the centre of the target canyon floor. For smaller aspect ratios, the relative contribution of the turbulent mass flux to net mass flux at the roof level, which was spatially averaged along the roof-level ventilation area, was closer to unity, indicating that turbulent motions mainly affected pollutant removal from the top of the canyon. As aspect ratio increased, the relative contribution became smaller, owing to strong upwind motions. However, the relative contribution again reached near unity for the infinite aspect ratio (i.e. a two-dimensional street canyon) because of lowered lateral flow convergence. At least 75 % of total emissions from the three-dimensional street canyon were attributable to turbulent motions. Pollutant removal by turbulent motions was related to the coherent structures of low-momentum fluid above the canyons. Though the coherent structure size of the low-momentum fluid differed, the positions of low-momentum fluid largely corresponded to instantaneous high concentrations of pollutant above the target canyon, irrespective of canyon geometry.     

Large-Eddy Simulation for the Mechanism of Pollutant Removal from a Two-Dimensional Street Canyon

Large-eddy simulation (LES) is conducted to investigate the mechanism of pollutant removal from a two-dimensional street canyon with a building-height to street-width (aspect) ratio of 1. A pollutant is released as a ground-level line source at the centre of the canyon floor. The mean velocities, turbulent fluctuations, and mean pollutant concentration estimated by LES are in good agreement with those obtained by wind-tunnel experiments. Pollutant removal from the canyon is mainly determined by turbulent motions, except in the adjacent area to the windward wall. The turbulent motions are composed of small vortices and small-scale coherent structures of low-momentum fluid generated close to the plane of the roof. Although both small vortices and small-scale coherent structures affect pollutant removal, the pollutant is largely emitted from the canyon by ejection of low-momentum fluid when the small-scale coherent structures appear just above the canyon where the pollutant is retained. Large-scale coherent structures also develop above the canyon, but they do not always affect pollutant removal.     

On the Mechanism of Air Pollutant Removal in Two-Dimensional Idealized Street Canyons: A Large-Eddy Simulation Approach

Flow resistance, ventilation, and pollutant removal for idealized two-dimensional (2D) street canyons of different building-height to street-width (aspect) ratios $AR$ are examined using the friction factor $f$ , air exchange rate (ACH), and pollutant exchange rate (PCH), respectively, calculated by large-eddy simulation (LES). The flows are basically classified into three characteristic regimes, namely isolated roughness, wake interference, and skimming flow, as functions of the aspect ratios. The LES results are validated by various experimental and numerical datasets available in the literature. The friction factor increases with decreasing aspect ratio and reaches a peak at $AR = 0.1$ in the isolated roughness regime and decreases thereafter. As with the friction factor, the ACH increases with decreasing aspect ratio in the wake interference and skimming flow regimes, signifying the improved aged air removal for a wider street canyon. The PCH exhibits a behaviour different from its ACH counterpart in the range of aspect ratios tested. Pollutants are most effectively removed from the street canyon with $AR = 0.5$ . However, a minimum of PCH is found nearby at $AR = 0.3$ , at which the pollutant removal is sharply weakened. Besides, the ACH and PCH are partitioned into the mean and turbulent components to compare their relative contributions. In line with our earlier Reynolds-averaged Navier–Stokes calculations (Liu et al., Atmos Environ 45:4763–4769, 2011), the current LES shows that the turbulent components contribute more to both ACH and PCH, consistently demonstrating the importance of atmospheric turbulence in the ventilation and pollutant removal for urban areas.     

Large-Eddy Simulation of Flow and Pollutant Transport in Urban Street Canyons with Ground Heating

Our study employed large-eddy simulation (LES) based on a one-equation subgrid-scale model to investigate the flow field and pollutant dispersion characteristics inside urban street canyons. Unstable thermal stratification was produced by heating the ground of the street canyon. Using the Boussinesq approximation, thermal buoyancy forces were taken into account in both the Navier–Stokes equations and the transport equation for subgrid-scale turbulent kinetic energy (TKE). The LESs were validated against experimental data obtained in wind-tunnel studies before the model was applied to study the detailed turbulence, temperature, and pollutant dispersion characteristics in the street canyon of aspect ratio 1. The effects of different Richardson numbers (Ri) were investigated. The ground heating significantly enhanced mean flow, turbulence, and pollutant flux inside the street canyon, but weakened the shear at the roof level. The mean flow was observed to be no longer isolated from the free stream and fresh air could be entrained into the street canyon at the roof-level leeward corner. Weighed against higher temperature, the ground heating facilitated pollutant removal from the street canyon.     

Numerical studies on flow fields around buildings in an Urban street canyon and cross-road

The questions on how vortices are constructed and on the relationship between the flow patterns and concentration distributions in real street canyons are the most pressing questions in pollution control studies. In this paper, the very large eddy simulation (VLES) and large eddy simulation (LES) are applied to calculate the flow and pollutant concentration fields in an urban street canyon and a cross-road respectively. It is found that the flow separations are not only related to the canyon aspect ratios, but also with the flow velocities and wall temperatures. And the turbulent dispersions are so strongly affected by the flow fields that the pollutant concentration distributions can be distinguished from the different aspect ratios, flow velocities and wall temperatures.     

城市街道峡谷汽车尾气污染的数值模拟

建立了一个简单的三维街道峡谷空气污染模式,并用实测资料进行了验证。利用建立的模式设计了7种试验方案,对街道峡谷内的污染状况进行了模拟。以CO为模拟对象的数值试验结果表明,街道峡谷上空的风速风向条件是决定街道峡谷内的污染状况的重要因素。峡谷上空风向与街道轴线的夹角越大、风速越小,则街道地面CO浓度越高。以现有的兰州典型车流量和排放因子,兰州街道地面CO浓度容易超标;若不控制车流量,到2008年,即使兰州上路汽车排放达标,但街道地面CO浓度仍然容易超标。     

Flow and Pollutant Transport in Urban Street Canyons of Different Aspect Ratios with Ground Heating: Large-Eddy Simulation

A validated large-eddy simulation model was employed to study the effect of the aspect ratio and ground heating on the flow and pollutant dispersion in urban street canyons. Three ground-heating intensities (neutral, weak and strong) were imposed in street canyons of aspect ratio 1, 2, and 0.5. The detailed patterns of flow, turbulence, temperature and pollutant transport were analyzed and compared. Significant changes of flow and scalar patterns were caused by ground heating in the street canyon of aspect ratio 2 and 0.5, while only the street canyon of aspect ratio 0.5 showed a change in flow regime (from wake interference flow to skimming flow). The street canyon of aspect ratio 1 does not show any significant change in the flow field. Ground heating generated strong mixing of heat and pollutant; the normalized temperature inside street canyons was approximately spatially uniform and somewhat insensitive to the aspect ratio and heating intensity. This study helps elucidate the combined effects of urban geometry and thermal stratification on the urban canyon flow and pollutant dispersion.     

Comprehensive Evaluation of Fast-Response,Reynolds-Averaged Navier–Stokes,and Large-Eddy Simulation Methods Against High-Spatial-Resolution Wind-Tunnel Data in Step-Down Street Canyons

Three computational fluid dynamics (CFD) methods with different levels of flow-physics modelling are comprehensively evaluated against high-spatial-resolution wind-tunnel velocity data from step-down street canyons (i.e., a short building downwind of a tall building). The first method is a semi-empirical fast-response approach using the Quick Urban Industrial Complex (QUIC-URB) model. The second method solves the Reynolds-averaged Navier–Stokes (RANS) equations, and the third one utilizes a fully-coupled fluid-structure interaction large-eddy simulation (LES) model with a grid-turbulence inflow generator. Unlike typical point-by-point evaluation comparisons, here the entire two-dimensional wind-tunnel dataset is used to evaluate the dynamics of dominant flow topological features in the street canyon. Each CFD method is scrutinized for several geometric configurations by varying the downwind-to-upwind building-height ratio (\(H_\mathrm{d}/H_\mathrm{u}\)) and street canyon-width to building-width aspect ratio (S / W) for inflow winds perpendicular to the upwind building front face. Disparities between the numerical results and experimental data are quantified in terms of their ability to capture flow topological features for different geometric configurations. Overall, all three methods qualitatively predict the primary flow topological features, including a saddle point and a primary vortex. However, the secondary flow topological features, namely an in-canyon separation point and secondary vortices, are only well represented by the LES method despite its failure for taller downwind building cases. Misrepresentation of flow-regime transitions, exaggeration of the coherence of recirculation zones and wake fields, and overestimation of downwards vertical velocity into the canyon are the main defects in QUIC-URB, RANS and LES results, respectively. All three methods underestimate the updrafts and, surprisingly, QUIC-URB outperforms RANS for the streamwise velocity component, while RANS is superior to QUIC-URB for the vertical velocity component in the street canyon.     

城市湍流边界层内汽车尾气扩散规律数值模拟研究

以纳维斯托克斯方程组、大气平流扩散方程、湍流动能及湍流动能耗散率方程组为基础．采用伪不定常方法,建立了一个数值模式．利用该模式列城市湍流边界层内流场结构及汽车排放污染物扩散规律进行了研究。结果表明：街谷内会形成一个涡旋型流场．汽车排放污染物浓度在地面及建筑物背风面产生堆积,且其沿高度方向的梯度变化在背风面大．迎风而小。随着街谷两侧建筑物屋顶风速的增大,峡谷内形成的涡旋流场的强度增大,污染物扩散速率增大：当屋顶来流与街道之间的夹角逐渐增大时．涡旋中心位置由街道中心偏向于背风面及更高层且污染物扩散速度加快。     

Turbulent Transfer Between Street Canyons and the Overlying Atmospheric Boundary Layer

The turbulent exchange of momentum between a two-dimensional cavity and the overlying boundary layer has been studied experimentally, using hot-wire anemometry and particle image velocimetry (PIV). Conditions within the boundary layer were varied by changing the width of the canyons upstream of the test canyon, whilst maintaining the square geometry of the test canyon. The results show that turbulent transfer is due to the coupling between the instabilities generated in the shear layer above the canyons and the turbulent structures in the oncoming boundary layer. As a result, there is no single, unique velocity scale that correctly characterizes all the processes involved in the turbulent exchange of momentum across the boundary layer. Similarly, there is no single velocity scale that can characterize the different properties of the turbulent flow within the canyon, which depends strongly on the way in which turbulence from the outer flow is entrained into the cavity and carried round by the mean flow. The results from this study will be useful in developing simple parametrizations for momentum exchange in the urban canopy, in situations where the street geometry consists principally of relatively long, uniform streets arranged in grid-like patterns; they are unlikely to be applicable to sparse geometries composed of isolated three-dimensional obstacles.     

Effects of Street-Bottom and Building-Roof Heating on Flow in Three-Dimensional Street Canyons

Using a computational fluid dynamics(CFD)model,the effects of street-bottom and building-roof heating on flow in three-dimensional street canyons are investigated.The building and street-canyon aspect ratios are one.In the presence of street-bottom heating,as the street-bottom heating intensity increases,the mean kinetic energy increases in the spanwise street canyon formed by the upwind and downwind buildings but decreases in the lower region of the streamwise street canyon.The increase in momentum due to buoyancy force intensifies mechanically induced flow in the spanwise street canyon.The vorticity in the spanwise street canyon strengthens.The temperature increase is not large because relatively cold above-roof-level air comes into the spanwise street canyon.In the presence of both street-bottom and building-roof heating,the mean kinetic energy rather decreases in the spanwise street canyon.This is caused by the decrease in horizontal flow speed at the roof level,which results in the weakening of the mean flow circulation in the spanwise street canyon.It is found that the vorticity in the spanwise street canyon weakens.The temperature increase is relatively large compared with that in the street-bottom heating case,because relatively warm above-roof-level air comes into the spanwise street canyon.     

Three-Dimensional Mapping of Air Flow at an Urban Canyon Intersection

In this experimental work both qualitative (flow visualisation) and quantitative (laser Doppler anemometry) methods were applied in a wind tunnel in order to describe the complex three-dimensional flow field in a real environment (a street canyon intersection). The main aim was an examination of the mean flow, turbulence and flow pathlines characterising a complex three-dimensional urban location. The experiments highlighted the complexity of the observed flows, particularly in the upwind region of the intersection. In this complex and realistic situation some details of the upwind flow, such as the presence of two tall towers, play an important role in defining the flow field within the intersection, particularly at roof level. This effect is likely to have a strong influence on the mass exchange mechanism between the canopy flow and the air aloft, and therefore the distribution of pollutants. This strong interaction between the flows inside and outside the urban canopy is currently neglected in most state-of-the-art local scale dispersion models.     

Numerical and experimental studies on flow and pollutant dispersion in urban street canyons

In this study numerical simulations and water tank experiments were used to investigate the flow and pollutant dispersion in an urban street canyon. Two types of canyon geometry were tested. The studies indicate that in a step-up notch canyon （higher buildings on the downstream side of the canyon）, the height and shape of the upstream lower buildings plays an important role in flow pattern and pollutant dispersion, while in a step-down notch canyon （lower buildings on the downstream side）, the downstream lower buildings have little influence. The studies also show that the substitution of tall towers for parailelepiped buildings on one side of the canyon may enhance the street ventilation and decrease the pollutant concentration emitted by motor vehicles.     

Large-Eddy Simulation of Flow and Pollutant Dispersion in High-Aspect-Ratio Urban Street Canyons with Wall Model

A large-eddy simulation (LES) with a one-equation subgrid-scale (SGS) model was developed to investigate the flow field and pollutant dispersion inside street canyons of high aspect ratio (AR). A 1/7th power-law wall model was implemented near rigid walls to mitigate the demanding near-wall resolution requirements in LES. This LES model had been extensively validated against experimental results for street canyons of AR = 1 and 2 before it was applied to the cases of AR = 3 and 5. A ground-level passive pollutant line source, located in the middle of the street, was used to simulate vehicular emissions. Three and five vertically aligned primary recirculations were developed in the street canyons of AR 3 and 5, respectively. The ground-level mean wind speed was less than 0.5% of the free stream value, which makes it difficult for the pollutant to be transported upward for removal. High pollutant concentration and variance were found near the buildings where the air flow is upwards. It was found that the velocity fluctuation, pollutant concentration and variance were all closely related to the interactions between the primary recirculations and/or the free surface layer. Several quantities, which are non-linear functions of AR, were introduced to quantify the air quality in street canyons of different configurations.     

Large-Eddy Simulation of Turbulent Organized Structures within and above Explicitly Resolved Cube Arrays

Large-eddy simulations have been performed for fully developed turbulent flow within and above explicitly resolved simple cube arrays. The results from our model, hereafter LES-CITY, are shown to agree with laboratory experiments. We investigated the systematic influence of cube density on turbulent flow characteristics by performing numerical experiments for cube areal densities from 0 to 44%. The following results were obtained: (1) The dispersive momentum flux was quite large within the canopy layer due to a mean stream re-circulation, whereas it was smaller above the canopy. The spatial variation of temporally averaged momentum in the roughness sub-layer was 20% or less of the total kinematic surface drag. (2) The temporally and spatially-averaged flow structure confirmed the existence of conventionally described canyon flow regimes; isolated, interfacial, and wake. However, the intermittency of the canyon flow for all cube densities was quite large and the stream patterns were never persistent. (3) Turbulent organized structures (TOS) similar to those observed in turbulent surface-layer flows were simulated, which are characterized by longitudinally-elongated low speed streaks and the corresponding shorter streamwise vortices. The streaks in sparse and dense canopy flows were likely to be aligned to the street line and to the roof lines, respectively. Such heterogeneity of TOS partially accounts for the large spatial variation of momentum flux. (4) In contrast to the mixing layer analogy of vegetation flows, the TOS and the resulting turbulent statistics of urban flow above the canopy resembled those in surface layers. The recirculation within the canopy significantly influenced the turbulent statistical properties.     

A Laboratory Model for the Flow in Urban Street Canyons Induced by Bottom Heating

Water tank experiments are carried out to investigate the convection flow induced by bottom heating and the effects of the ambient wind on the flow in non-symmetrical urban street canyons based on the PIV (Particle Image Visualization) technique. Fluid experiments show that with calm ambient wind,the flows in the street canyon are completely driven by thermal force, and the convection can reach the upper atmosphere of the street canyon. Horizontal and vertical motions also appear above the roofs of the buildings. These are the conditions which favor the exchange of momentum and air mass between the street canyon and its environment. More than two vortices are induced by the convection, and the complex circulation pattern will vary with time in a wider street canyon. However， in a narrow street canyon， just one vortex appears. With a light ambient wind, the bottom heating and the associated convection result in just one main vortex. As the ambient wind speed increases, the vortex becomes more organized and its center shifts closer to the leeward building.     

Experimental study of the transfer velocity for urban surfaces with a water evaporation method

A major problem in urban climate modelling is determining how the heat fluxes from various canyon surfaces are affected by canyon flow. To address this problem, we developed a water evaporation method involving filter paper to study the distribution of the convective transfer velocity in urban street canyons. In this method, filter paper is pasted onto a building model and the evaporation rate from the paper is measured with an electric balance. The method was tested on 2D (two-dimensional) street canyon models and 3D model arrangements. Moreover, in this technique, it is easy to restrict the flux within an arbitrary surface in question. That is, the evaporation distribution on a surface can be studied by using several small pieces of filter paper. In the 2D case, the wall transfer velocity was strongly dependent on the canyon aspect ratio for perpendicular wind directions and it varied widely with height within both windward and leeward wall surfaces. For 3D cubic arrays, the relation to canyon aspect ratio was largely different from that of the 2D canyon. And, as a case study, the variation of wind direction was investigated for a city-like setting. The area-averaged transfer velocity was insensitive to wind direction but its local deviation was significant. Finally, we measured the transfer velocity for a clustered block array surrounded by relatively wide streets. The effect of spatial heterogeneity on the transfer velocity was significant. Moreover, for a fixed total building volume, the transfer velocity was considerably larger when the building height varied than when it was uniform. Therefore, the water evaporation method with filter paper is expected to be useful for studying the transfer velocity and ventilation rates in urban areas with various canyon shapes.     

二维街谷地面加热引起的流场特征的水槽实验研究

利用拖曳式水槽,采用激光粒子成像速度场测量系统(PIV),模拟了街谷存在地面加热时流场特征;讨论了环境风场对其的影响。我们发现在静风条件下,街谷中环流完全由热力驱动,对流活动可伸展至街谷上方;在建筑物层顶以上,也可发现水平和垂直方向的运动。这些对流活动有助于基本风场为零时,街谷内外动量和物质的交换。当街谷较宽时,对流形成的涡旋可能为两个以上,形态较为复杂并随时间变化,当街谷变窄时,涡旋蜕化成只有一个。当有弱环境风场存在时,街谷中的对流呈现为一个主涡旋,随着风速增加,涡旋形状更加规则,其中心并向下风向移动。     

An Improved Three-Dimensional Simulation of the Diurnally Varying Street-Canyon Flow

The impact of diurnal variations of the heat fluxes from building and ground surfaces on the fluid flow and air temperature distribution in street canyons is numerically investigated using the PArallelized Large-eddy Simulation Model (PALM). Simulations are performed for a 3 by 5 array of buildings with canyon aspect ratio of one for two clear summer days that differ in atmospheric instability. A detailed building energy model with a three-dimensional raster-type geometry—Temperature of Urban Facets Indoor-Outdoor Building Energy Simulator (TUF-IOBES)—provides urban surface heat fluxes as thermal boundary conditions for PALM. In vertical cross-sections at the centre of the spanwise canyon the mechanical forcing and the horizontal streamwise thermal forcing at roof level outweigh the thermal forces from the heated surfaces inside the canyon in defining the general flow pattern throughout the day. This results in a dominant canyon vortex with a persistent speed, centered at a constant height. Compared to neutral simulations, non-uniform heating of the urban canyon surfaces significantly modifies the pressure field and turbulence statistics in street canyons. Strong horizontal pressure gradients were detected in streamwise and spanwise canyons throughout the day, and which motivate larger turbulent velocity fluctuations in the horizontal directions rather than in the vertical direction. Canyon-averaged turbulent kinetic energy in all non-neutral simulations exhibits a diurnal cycle following the insolation on the ground in both spanwise and streamwise canyons, and it is larger when the canopy bottom surface is paved with darker materials and the ground surface temperature is higher as a result. Compared to uniformly distributed thermal forcing on urban surfaces, the present analysis shows that realistic non-uniform thermal forcing can result in complex local airflow patterns, as evident, for example, from the location of the vortices in horizontal planes in the spanwise canyon. This study shows the importance of three-dimensional simulations with detailed thermal boundary conditions to explore the heat and mass transport in an urban area.