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

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

Modelling Pollutant Dispersion in a Street Network

This study constitutes a further step in the analysis of the performances of a street network model to simulate atmospheric pollutant dispersion in urban areas. The model, named SIRANE, is based on the decomposition of the urban atmosphere into two sub-domains: the urban boundary layer, whose dynamics is assumed to be well established, and the urban canopy, represented as a series of interconnected boxes. Parametric laws govern the mass exchanges between the boxes under the assumption that the pollutant dispersion within the canopy can be fully simulated by modelling three main bulk transfer phenomena: channelling along street axes, transfers at street intersections, and vertical exchange between street canyons and the overlying atmosphere. Here, we aim to evaluate the reliability of the parametrizations adopted to simulate these phenomena, by focusing on their possible dependence on the external wind direction. To this end, we test the model against concentration measurements within an idealized urban district whose geometrical layout closely matches the street network represented in SIRANE. The analysis is performed for an urban array with a fixed geometry and a varying wind incidence angle. The results show that the model provides generally good results with the reference parametrizations adopted in SIRANE and that its performances are quite robust for a wide range of the model parameters. This proves the reliability of the street network approach in simulating pollutant dispersion in densely built city districts. The results also show that the model performances may be improved by considering a dependence of the wind fluctuations at street intersections and of the vertical exchange velocity on the direction of the incident wind. This opens the way for further investigations to clarify the dependence of these parameters on wind direction and street aspect ratios.     

Flow in a Street Canyon for any External Wind Direction

An analytical model has been developed for the flow along a street canyon (of height H and width W), generated by an external wind blowing at any angle relative to the axis of the street. Initially, we consider the special case of a wind blowing parallel to the street. The interior of the street is decomposed into three regions, and the flow within each region is assumed to depend only on the external wind and the distance to the closest solid boundary. This decomposition leads to two different flow regimes: one for narrow streets (H/W > 1/2) and one for wide streets (H/W < 1/2). The theoretical model agrees well with results obtained from numerical simulations using a Reynolds-Averaged Navier–Stokes model. We then generalize the model to the case of arbitrary wind direction. Numerical solutions show that the streamlines of the mean flow in the street have a spiral form, and for most angles of incidence, the mass flux along the street scales on the component of the external wind resolved parallel to the street. We use this result to generalize the model derived for wind blowing parallel to the street, and the results from this model agree well with the numerical simulations. The model that has been developed can be evaluated rapidly using only very modest computing power, so it is suitable for use as an operational tool.     

Effects of Time-Dependent Inflow Perturbations on Turbulent Flow in a Street Canyon

Urban flow and turbulence are driven by atmospheric flows with larger horizontal scales. Since building-resolving computational fluid dynamics models typically employ steady Dirichlet boundary conditions or forcing, the accuracy of numerical simulations may be limited by the neglect of perturbations. We investigate the sensitivity of flow within a unit-aspect-ratio street canyon to time-dependent perturbations near the inflow boundary. Using large-eddy simulation, time-periodic perturbations to the streamwise velocity component are incorporated via the nudging technique. Spatial averages of pointwise differences between unperturbed and perturbed velocity fields (i.e., the error kinetic energy) show a clear dependence on the perturbation period, though spatial structures are largely insensitive to the time-dependent forcing. The response of the error kinetic energy is maximized for perturbation periods comparable to the time scale of the mean canyon circulation. Frequency spectra indicate that this behaviour arises from a resonance between the inflow forcing and the mean motion around closed streamlines. The robustness of the results is confirmed using perturbations derived from measurements of roof-level wind speed.     

Numerical Study on the Impact of Ground Heating and Ambient Wind Speed on Flow Fields in Street Canyons

The impact of ground heating on flow fields in street canyons under different ambient wind speed conditions was studied based on numerical methods.A series of numerical tests were performed,and three factors including height-to-width(H/W) ratio,ambient wind speed and ground heating intensity were taken into account.Three types of street canyon with H/W ratios of 0.5,1.0 and 2.0,respectively,were used in the simulation and seven speed values ranging from 0.0 to 3.0 m s 1 were set for the ambient wind speed.The ground heating intensity,which was defined as the difference between the ground temperature and air temperature,ranged from 10 to 40 K with an increase of 10 K in the tests.The results showed that under calm conditions,ground heating could induce circulation with a wind speed of around 1.0 m s 1,which is enough to disperse pollutants in a street canyon.It was also found that an ambient wind speed threshold may exist for street canyons with a fixed H/W ratio.When ambient wind speed was lower than the threshold identified in this study,the impact of the thermal effect on the flow field was obvious,and there existed a multi-vortex flow pattern in the street canyon.When the ambient wind speed was higher than the threshold,the circulation pattern was basically determined by dynamic effects.The tests on the impact of heating intensity showed that a higher ground heating intensity could strengthen the vortical flow within the street canyon,which would help improve pollutant diffusion capability in street canyons.     

Micrometeorological Measurements in a Street Canyon during the Joint ATREUS-PICADA Experiment

In order to investigate the microclimatic conditions in a street canyon, a physical model was used to conduct the Joint ATREUS-PICADA Experiment (JAPEX) in situ experimental campaign. Four lines of buildings simulated by steel containers were installed to form three parallel street canyons at 1:5 scale, with width/height aspect ratio approximately 0.40. The reference wind and atmospheric conditions were measured, as well as the flow velocity and direction in the street. Preliminary results concern street canyon ventilation and thermal effects on in-canyon airflow, and show that vortical motions appear for reference wind directions perpendicular to the street axis. The presence of adjacent rows of buildings did not appear to significantly influence the flow character within the canyon for the case of a low aspect ratio corresponding to a skimming flow regime. The flow structure was not significantly affected by the thermal effects although some slight interference occurred in the lower part of the canyon. An analysis of horizontal temperature gradients indicated that a thin boundary layer develops near the heated facade. These facts imply that the thermal effects are considerable only very close to the wall.     

Scalar Fluxes from Urban Street Canyons Part II: Model

A practical model is developed for the vertical flux of a scalar, such as heat, from an urban street canyon that accounts for variations of the flow and turbulence with canyon geometry. The model gives the magnitude and geometric dependence of the flux from each facet of the urban street canyon, and is shown to agree well with wind-tunnel measurements described in Part I. The geometric dependence of the flux from an urban street canyon is shown to be determined by two physical processes. Firstly, as the height-to-width ratio of the street canyon increases, so does the roughness length and displacement height of the surface. This increase leads to a reduction in the wind speed in the inertial sublayer above the street canyons. Since the speed of the circulations in the street are proportional to this inertial sublayer wind speed, the flux then reduces with the inertial sublayer wind speed. This process is dominant at low height-to-width ratios. Secondly, the character of the circulations within the street canyon also varies as the height-to-width ratio increases. The flow in the street is partitioned into a recirculation region and a ventilated region. When the street canyon has high height-to-width ratios the recirculation region occupies the whole street canyon and the wind speeds within the street are low. This tendency decreases the flux at high height-to-width ratios. These processes tend to reduce the flux density from the individual facets of the street canyon, when compared to the flux density from a horizontal surface of the same material. But the street canyon has an increased total surface area, which means that the total flux from the street canyon is larger than from a horizontal surface. The variations in scalar flux from an urban street canyon with geometry is over a factor of two, which means that the physical mechanisms responsible should be incorporated into energy balance models for urban areas.     

Corrections to Inertial-Range Power Spectra Measured by Csat3 and Solent Sonic Anemometers,1. Path-Averaging Errors

Sonic anemometer path-averaging errors are determined for measurements of inertial-range velocity and temperature power spectra as a function of k ₁ p, where k ₁ is the streamwise wave number and p is the sonic path length. The attenuation of vertical velocity spectra is found to be quite similar for the CSAT3 and Solent anemometers and to be insensitive to wind direction. The attenuation of the horizontal wind component spectra is noticeably greater for the Solent sonics than for the CSAT3, but those for the CSAT3 have a greater dependence on wind direction. The attenuation of sonic temperature spectra is also found to be quite similar for the CSAT3 and Solent R3 sonics and to be insensitive to wind direction, while that of the Solent R2 is less than that of the other two sonics and has a minor dependence on wind direction.     

台风“利奇马”引导气流的计算方法研究

为了更准确计算台风的引导气流,基于相关系数计算方法,评估了各层环境风与台风移动速度的相关性。结果表明,剔除相关性差的低层环境风后,计算得到的引导气流与台风移动的一致性有所提高;引导气流分量的峰值谷值在时间上滞后于台风的峰谷值;低层风之所以与台风相关性差或与低空急流和地形的干扰有关。针对性调整引导气流计算方法,有助于提高引导气流与台风移动速度的相关性。     

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.     

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.     

基于单多普勒雷达的光流法改进技术

本文对基于径向速度的光流法在实际应用中存在的误差进行了分析,并提出了有效的改进方案。合成数据试验结果表明:基于径向速度的光流法在反演低仰角(小于2.5°)数据和包含非线性变化的实际风场时存在较大的误差。在低仰角情况下,由于径向速度在垂直方向的分量非常小,导致垂直速度的反演包含较大的计算误差并影响水平速度的反演精度。而当实际风场存在较大非线性变化时(如锋面、切变线等),算法中的局部约束条件被破坏,导致切变区域反演精度降低或出现奇异值。针对这两种情况,本文提出了有效的改进算法。在低仰角情况下引入了简化的约束方程,并加入了径向速度一致性的约束条件来克服实际风场非线性带来的误差。通过理想数据试验和实际个例的检验发现,改进后的算法可以有效地提高风场反演的精度和应用范围。     

The influence of street architecture on flow and dispersion in street canyons

Summary The paper presents an overview of the influence of street architecture on the wind and turbulence patterns in street canyons and discusses the effects on local air quality. The findings of recent experimental and numerical studies are summarized and wind-tunnel data sets are presented that illustrate the flow-field variability. It is shown that small-scale features of the street architecture play an important role. The formation of a vortex inside the street canyon is affected by the roof configuration. In shorter street canyons, the flow component along the street becomes important for pollutant transport. These results are of importance for urban air quality modeling in particular when dealing with pollution problems caused by road traffic. Furthermore, the findings should be taken into account in fast response models that are used to assess critical areas in the case of accidental or non-accidental releases of hazardous material in urban areas.     

Remarks on the Definition and Estimation of Friction Velocity

One of the mainscaling parameters in similarity theory of the atmospheric boundary layer is friction velocity. Unfortunately, several definitions of friction velocity exist in the literature. Some authors use the component of the horizontal Reynolds stress vector in the direction of the mean wind vector to define friction velocity. Others define the friction velocity by means of the absolute value of the horizontal Reynolds stress vector. The two definitions coincide only if the direction of the mean wind vector is parallel to the horizontal Reynolds stress vector. In general, the second definition gives larger values for the friction velocity. Over complex terrain the situation is further complicated by the fact that the terrain following flow is not necessarily horizontal. Thus, several authors have proposed to use terrain following coordinate systems for the definition of friction velocity. By means of a large dataset of fast-response wind measurements with an ultrasonic anemometer the friction velocities resulting from the different definitions are compared. Furthermore, it is shown that friction velocity can be well estimated from horizontal wind speed, and even better from simple horizontal or vertical turbulence parameters.     

Nocturnal boundary-layer height: Observations by acoustic sounders and predictions in terms of surface-layer parameters

Acoustic sounder measurements of the stable boundary-layer height taken during the EPRI Plume Model Validation and Development Project experiment are examined. Comparison of simultaneous measurements by two sodars located 15 km apart shows good agreement. Several widely used diagnostic formulas for estimation of the boundary-layer height, based on wind speed and surface-layer parameters, such as friction velocity and Monin-Obukhov length, are tested against the sodar data. Of these, best performance is found using a simple linear relationship with friction velocity or, alternatively, wind speed at 10 m height. No evidence is found to support the more often used Zilitinkevich (1972) formula. Tests using selected data from the Cabauw site in the Netherlands confirm the results found on the basis of EPRI data.     

Large-Scale Structures over a Single Street Canyon Immersed in an Urban-Type Boundary Layer

An analysis of the dynamics of the flow over a street canyon immersed in an atmospheric boundary layer is presented, using particle image velocimetry measurements in a wind tunnel. Care was taken to generate a 1:200 model scale urban type boundary layer that is correctly scaled to the size of the canyon buildings. Using proper orthogonal decomposition (POD) of the velocity field and conditional averaging techniques, it is first shown that the flow above the opening of the canyon consists of a shear layer separating from the upstream obstacle, animated by a coherent flapping motion and generating large-scale vortical structures. These structures are alternately injected into the canyon or shed off the obstacle into the outer flow. It is shown that unsteady fluid exchanges between the canyon and the outer flow are mainly driven by the shear layer. Finally, using POD, the non-linear interaction between the large-scale structures of the oncoming atmospheric boundary layer and the flow over the canyon is demonstrated.     

街谷环流和热力结构的数值模拟

文章介绍一种用于模拟街谷流场和温度场的动力学模式和热力学模式.应用动力模式模拟了方柱体塔楼和圆柱体塔楼形成的流场,应用动力和热力模式模拟了街谷中流场和温度场的日变化过程.计算实例表明,上述模式可用于城市街谷和建筑群风环境和热力环境研究以及街谷中空气污染物传输和扩散的计算.     

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

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

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.     

Mean Flow and Turbulence Characteristics in an Urban Roughness Sublayer

In this study, a detailed model of an urban landscape has been re-constructed inthe wind tunnel and the flow structure inside and above the urban canopy has beeninvestigated. Vertical profiles of all three velocity components have been measuredwith a Laser-Doppler velocimeter, and an extensive analysis of the measured meanflow and turbulence profiles carried out. With respect to the flow structure inside thecanopy, two types of velocity profiles can be distinguished. Within street canyons,the mean wind velocities are almost zero or negative below roof level, while closeto intersections or open squares, significantly higher mean velocities are observed.In the latter case, the turbulent velocities inside the canopy also tend to be higherthan at street-canyon locations. For both types, turbulence kinetic energy and shearstress profiles show pronounced maxima in the flow region immediately above rooflevel.Based on the experimental data, a shear-stress parameterization is proposed, inwhich the velocity scale, u_s, and length scale, z_s, are based on the level and magnitude of the shear stress peak value. In order to account for a flow region inside the canopy with negligible momentum transport, a shear stress displacement height, d_s, is introduced. The proposed scaling and parameterization perform well for the measured profiles and shear-stress data published in the literature.The length scales derived from the shear-stress parameterization also allowdetermination of appropriate scales for the mean wind profile. The roughnesslength, z₀, and displacement height, d₀, can both be described as fractions of the distance, z_s - d_s, between the level of the shear-stress peak and the shear-stress displacement height. This result can be interpreted in such a way that the flow only feels the zone of depth z_s - d_s as the roughness layer. With respect to the lower part of the canopy (z < d_s) the flow behaves as a skimming flow. Correlations between the length scales z_s and d_s and morphometric parameters are discussed.The mean wind profiles above the urban structure follow a logarithmic windlaw. A combination of morphometric estimation methods for d₀ and z₀ with wind velocity measurements at a reference height, which allow calculation of the shear-stress velocity, u_*, appears to be the most reliable and easiest procedure to determine mean wind profile parameters. Inside the roughnesssublayer, a local scaling approach results in good agreement between measuredand predicted mean wind profiles.