城市冠层中湍流运动的统计特征

对１９９７ 年夏天和冬天北京湍流运动的各种统计特征量进行了初步的统计分析。结果表明,城市冠层中湍流运动的各种统计特征量与平坦下垫面条件下边界层湍流运动的相比, 有不同的地方也有相似的地方; 无论白天还是夜晚, 垂直方向的湍流强度和湍流脉动风速标准差均小于水平方向的, 水平方向的相应湍流特征量则总是接近相等; 城市冠层中湍流脉动强度和标准差几乎均大于平坦下垫面边界层的; 平均风速ｕ≥１ ｍ ／ｓ 时的湍流统计特征量与ｕ＜ １ ｍ ／ｓ 时的有所不同; 城市冠层的阻力系数较大, 可达００６２５,Ｐａｎｏｆｓｋｙ 等提出的公式σｗ ／ｕ＊ ＝ １３ （１－ ３ｚ／Ｌ）１／３在城市冠层中并不适用。     

Multifractal Characteristics of Intermittent Turbulence in the Urban Canopy Layer

The multifractality of energy and thermal dissipation of fully developed intermittent turbulence is investigated in the urban canopy layer under unstable conditions by the singularity spectrum for the fractal dimensions of sets of singularities characterizing multifractals. In order to obtain high-order moment properties of smallscale turbulent dissipation in the inertial range, an ultrasonic anemometer with a high sampling frequency of 100 Hz was used. The authors found that the turbulent signal could be singular everywhere. Moreover, the singular exponents of energy and thermal dissipation rates are most frequently encountered at around 0.2, which is significantly smaller than the singular exponents for a wind tunnel at a moderate Reynolds number. The evidence indicates a higher intermittency of turbulence in the urban canopy layer at a high Reynolds number, which is demonstrated by the data with high temporal resolution. Furthermore, the temperature field is more intermittent than the velocity field. In addition, a large amount of samples could be used for verification of the results.     

Canopy Architecture and Turbulence Structure in a Coniferous Forest

Synchronous sonic anemometric measurements at five heightswithin a mixed coniferous forest were used to test two different parameterisations ofcanopy architecture in the application of a second-order turbulence closure model. Inthe computation of the leaf drag area, the aerodynamic sheltering was replaced with anarchitectural sheltering, assumed to be analogous to the clumping index defined in radiativetransfer theory. Consequently, the ratio of leaf area density and sheltering factor was approximatedby the effective leaf area or the mean contact number, both obtained from the inversion of non-destructive optical measurements. The first parameter represents the equivalentrandomly dispersed leaf area in terms of shading, the second is the average number of leavesthat a straight line intercepts penetrating the canopy with a certain zenith angle. Theselection of this direction was determined by the analysis of the mean angle of the wind vectorduring sweep events. The drag coefficient values obtained from the inversion of themomentum flux equation, using the two proposed parameterisations, are in good agreement withvalues found in the literature. The predicted profiles of turbulence statistics reasonablymatch actual measurements, especially in the case of the mean contact numberparameterisation. The vertical profile of leaf drag area, obtained by forcing the turbulence modelto match the observed standard deviation of vertical velocity (_w), is intermediatebetween the two empirical ones. Finally, the proposed canopy parameterisations were appliedto a Lagrangian transport model to predict vertical profiles of air temperature, H₂O andCO₂ concentration.     

Turbulence Structure Within and Above a Canopy of Bluff Elements

Measurements of turbulence structure in a wind-tunnel model canopy of bluff elements show many of the features associated with vegetation canopies and roughness sublayers but also display features more characteristic of the inertial sublayer (ISL). Points of similarity include the existence of an inflexion point in the space-time averaged streamwise velocity at the canopy top, the variation with height of turbulent second moments and the departure of the turbulent kinetic energy budget from local equilibrium in and just above the canopy. Quadrant analysis shows characteristic dominance of sweep over ejection events within the canopy although sweeps are more frequent than usually seen in vegetation canopies. Points of difference are a u′, w′ correlation coefficient that is closer to the ISL value than to most canopy data, and a turbulent Prandtl number midway between canopy and ISL values. Within the canopy there is distinct spatial partitioning into two flow regimes, the wake and non-wake regions. Both time-mean and conditional statistics take different values in these different regions of the canopy flow. We explain many of these features by appealing to a modified version of the mixing-layer hypothesis that links the dominant turbulent eddies to the instability of the inflexion point at canopy top. However, it is evident that these eddies are perturbed by the quasi-coherent wakes of the bluff canopy elements. Based upon an equation for the instantaneous velocity perturbation, we propose a criterion for deciding when the eddies linked to the inflexion point will dominate flow structure and when that structure will be replaced by an array of superimposed element wakes. In particular, we show that the resemblance of some features of the flow to the ISL does not mean that ISL dynamics operate within bluff-body canopies in any sense.     

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.     

用连续子波变换提取城市冠层大气湍流的相干结构

切变湍流的相干结构是湍流研究中的重大发现,它表明湍流在表面上看来不规则运动中具有可检测的有序运动,这种相干结构在切变湍流的脉动生成和发展中起着主宰作用.因此识别和提取相干结构对于认识和研究湍流是非常重要的.用数字滤波法将包含相干结构的大尺度信号提取出来以后,再用子波分析,根据子波能量极大值的判别方法,分别确定出大气湍流三个方向上的速度脉动信号相干结构的频率或时间尺度,然后由确定尺度上的连续子波反演公式,提取出大气湍流三个方向上的速度脉动信号相干结构所对应的波形.     

Vertical Structure Of Turbulence In Offshore Flow During Rasex

The adjustment of the boundary layer immediately downstream froma coastline is examined based on two levels of eddy correlation data collected on a mast at the shore and six levels of eddy correlation data and profiles of mean variables collected from a mast 2 km offshore during the Risø Air-Sea Experiment. The characteristics of offshore flow are studied in terms of case studies and inter-variable relationships for the entire one-month data set. A turbulent kinetic energy budget is constructed for each case study.The buoyancy generation of turbulence is small compared to shear generation and dissipation. However, weakly stable and weakly unstable cases exhibit completely different vertical structure. With flow of warm air from land over cooler water, modest buoyancy destruction of turbulence and reduced shear generation of turbulence over the less rough sea surface cause the turbulence to rapidly weaken downstream from the coast. The reduction of downward mixing of momentum by the stratification leads to smaller roughness lengths compared to the unstable case. Shear generation at higher levels and advection of stronger turbulence from land often lead to an increase of stress and turbulence energy with height and downward transport of turbulence energy toward the surface.With flow of cool air over a warmer sea surface, a convective internal boundary layer develops downstream from the coast. An overlying relatively thick layer of downward buoyancy flux (virtual temperature flux) is sometimes maintained by shear generation in the accelerating offshore flow.     

山谷城市的近地层大气湍流谱特征

山谷城市地区的大气湍流结构研究至今报道很少.本文根据兰州市一次观测实验,对各风速分量谱、温度谱以及动量和感热通量的协谱进行了分析.结果表明,尽管仪器的安装高度较低(6.45m),湍谱特征仍与平坦均匀下垫面得到的一些典型结果基本一致,从而认为,在类似地区进行湍流通量等的直接观测时,对观测高度和地点的选择,可以不必象一些文献要求的那样高.     

Mean Winds Through an Inhomogeneous Urban Canopy

The mean flow within inhomogeneous urban areas is investigated using an urban canopy model. The urban canopy model provides a conceptual and computational tool for representing urban areas in a way suitable for parameterisation within numerical weather prediction and urban air quality models. Average aerodynamic properties of groups of buildings on a neighbourhood scale can be obtained in terms of the geometry and layout of the buildings. These canopy parameters then determine the spatially averaged mean wind speeds within the canopy as a whole. Using morphological data for real cities, computations are performed for representative sections of cities. Simulations are performed to study transitions between different urban neighbourhoods, such as residential areas and city centres. Such transitions are accompanied by changes in mean building density and building height. These are considered first in isolation, then in combination, and the generic effects of each type of change are identified. The simulation of winds through a selection of downtown Los Angeles is considered as an example. An increase in canopy density is usually associated with a decrease in the mean wind speed. The largest difference between mean winds in canopies of different densities occurs near ground level. Winds generally decrease upon encountering a taller canopy of the same density, but this effect may be reversed very near the ground, with possible speed-ups if the canopy is especially tall. In the vicinity of a transition there is an overshoot in the mean wind speed in the bottom part of the canopy. Mechanisms for these effects are discussed.     

The Effect of the Vertical Source Distribution on Scalar Statistics within and above a Forest Canopy

Little is known about in-canopy processes that may alter forest–atmosphere exchanges of trace gases and aerosols. To improve our understanding of in-canopy mixing, we use large-eddy simulation to study the effect of scalar source/sink distributions on scalar concentration moments, fluxes, and correlation coefficients within and above an ideal forest canopy. Scalars are emitted from: (1) the ground, (2) the canopy, and (3) both the ground and the canopy; a scalar is also deposited onto the canopy. All scalar concentration moments, fluxes, and correlation coefficients are affected by the source location/distribution, as is the scalar segregation intensity. We conclude that vertical source/sink distribution has a profound impact on scalar concentration profiles, fluxes, correlation coefficient, and scalar segregation.     

Temperature Profiles Within and Above an Urban Canopy

The seasonal trend of vertical temperature profiles within and above an urban canopy has been investigated. We measured air temperatures and wind velocities along a 29-m tower in a residential area of Tokyo, Japan continually for 14 months. The height of the daily maximum temperature Z_Tmax varied with the season; Z_Tmax was at the roof level in winter but near the ground in summer. The seasonal change of Z_Tmax is likely due to the change of height at which solar energy is absorbed. At the time of the maximum temperature, the atmosphere above the canopy is always unstable whereas the air within the canopy is unstable in summer but stable in winter.     

Different Turbulent Regimes and Vertical Turbulence Structures of the Urban Nocturnal Stable Boundary Layer

Turbulence in the nocturnal boundary layer(NBL) is still not well characterized, especially over complex underlying surfaces. Herein, gradient tower data and eddy covariance data collected by the Beijing 325-m tower were used to better understand the differentiating characteristics of turbulence regimes and vertical turbulence structure of urban the NBL. As for heights above the urban canopy layer(UCL), the relationship between turbulence velocity scale(V_TKE) and wind speed(V) was con...     

A k - *epsiv Turbulence Closure Model For The Atmospheric Boundary Layer Including Urban Canopy

A numerical model for the computation of the wind field,air temperature and humidity in the atmospheric boundary layer (ABL) including the urbancanopy was developed for urban climate simulation. The governing equations of the modelare derived by applying ensemble and spatial averages to the Navier–Stokes equation, continuityequation and equations for heat and water vapour transfer in the air. With the spatial averagingprocedure, effects of buildings and other urban structures in the urban canopy can be accounted for byintroducing an effective volume function, defined as the ratio between the volume of air in acomputational mesh over the total volume of the mesh. The improved k - model accounts for the anisotropyof the turbulence field under density stratification. In the improved k - model, the transportof momentum and heat in the vertical direction under density stratification is evaluated based onthe assumption of a near-equilibrium shear flow where transport effects on the stresses andheat fluxes are negligible. The heating processes at surfaces of buildings and ground are alsomodelled. The comparison of the computational results obtained with the present modeland existing observational data and numerical models shows that the present model is capableof predicting the structure of turbulence in the urban canopy layer under density stratification.Numerical experiments with the new model show that the flow behaviour of the air in the urbancanopy layer is strongly affected by the existence of buildings and density stratification.     

南京夏季城郊湍流统计特性及湍流通量对比分析

利用2010年南京夏季城市热岛三维结构试验中的湍流观测资料,对南京市委党校和南京信息工程大学观测场2个观测点的湍流统计量和湍流通量进行了对比分析。结果表明：城市平均风速大于郊区;城郊2地的湍流强度都呈现出Ix〉Iy〉Iz的趋势,水平方向上的湍强城郊相差不大,而垂直方向上,郊区的Iz要明显低于城市;城、郊的风速归一化标准差σu/u＊、σv/u＊、σw/u＊都符合MoninObukhov相似理论的＂1/3＂定律,城市的垂直风速归一化标准差大于郊区,而水平方向上的风速归一化标准差则明显小于其他下垫面;郊区风速归一化标准差与常熟农田的拟合结果较为相似,而城市风速归一化标准差与长白山森林的拟合结果更为接近;夏季城市以感热通量为主,而郊区湍流能量的输送以潜热通量为主,且城市的湍流热通量受太阳辐射的影响更大。     

Turbulence in an almond orchard: Vertical variations in turbulent statistics

Three-dimensional wind velocity components were measured above and within a uniform almond orchard. Turbulent statistics associated with the turbulent flow inside the canopy are examined in detail. Turbulence in an almond orchard is characterized by relatively high turbulent intensities and large skewness and kurtosis values. These results indicate that the frequency distribution of wind velocity components is non-Gaussian. Conditional sampling of the turbulent measurements show that large, infrequent sweeps provide the predominant mechanism for tangential momentum stress in the canopy crown. Deep inside the canopy, a secondary wind maximum and small, but positive, tangential momentum stresses are observed.     

The Spatial Variability of Turbulence above a Forest

Summary The spatial variability of turbulence above a forest has been examined. Two measurement towers were erected 800 m apart within a heterogeneous mixed forest located in the north east of the Netherlands. The measurements of u ^* /u were analysed and subsequently used to test a surface layer model. The model simulated the magnitude of the measurements reasonably well, but measured trends were not always reproduced by the model. The variable (du/dz)/u did not adapt as quickly to the new surface as u ^* /u. This is in agreement with Schmid (1994), and can be explained by a local decrease in mixing length. It is recommended to adapt the mixing length near a surface transition to improve the accuracy of surface layer models of heterogeneous landscapes. Received September 8, 1997 Revised April 15, 1998     

Turbulent Transport of Momentum and Scalars Above an Urban Canopy

Turbulent transport of momentum and scalars over an urban canopy is investigated using the quadrant analysis technique. High-frequency measurements are available at three levels above the urban canopy (47, 140 and 280 m). The characteristics of coherent ejection–sweep motions (flux contributions and time fractions) at the three levels are analyzed, particularly focusing on the difference between ejections and sweeps, the dissimilarity between momentum and scalars, and the dissimilarity between the different scalars (i.e., temperature, water vapour and $\hbox {CO}_{2})$ . It is found that ejections dominate momentum and scalar transfer at all three levels under unstable conditions, while sweeps are the dominant eddy motions for transporting momentum and scalars in the urban roughness sublayer under neutral and stable conditions. The flux contributions and time fractions of ejections and sweeps can be adequately captured by assuming a Gaussian joint probability density function for flow variables. However, the inequality of flux contributions from ejections and sweeps is more accurately reproduced by the third-order cumulant expansion method (CEM). The incomplete cumulant expansion method (ICEM) also works well except for $\hbox {CO}_{2}$ at 47 m where the skewness of $\hbox {CO}_{2}$ fluctuations is significantly larger than that for vertical velocity. The dissimilarity between momentum and scalar transfers is linked to the dissimilarity in the characteristics of ejection–sweep motions and is further quantified by measures of transport efficiencies. Atmospheric stability is the controlling factor for the transport efficiencies of momentum and heat, and fitted functions from the literature describe their behaviour fairly accurately. However, transport efficiencies of water vapour and $\hbox {CO}_{2}$ are less affected by the atmospheric stability. The dissimilarity among the three scalars examined in this study is linked to the active role of temperature and to the surface heterogeneity effect.     

The Effect of Vegetation Density on Canopy Sub-Layer Turbulence

The canonical form of atmospheric flows near theland surface, in the absence of a canopy, resembles a rough-wallboundary layer. However, in the presence of an extensive and densecanopy, the flow within and just above the foliage behaves as aperturbed mixing layer. To date, no analogous formulation existsfor intermediate canopy densities. Using detailed laser Dopplervelocity measurements conducted in an open channel over a widerange of canopy densities, a phenomenological model that describesthe structure of turbulence within the canopy sublayer (CSL) isdeveloped. The model decomposes the space within the CSL intothree distinct zones: the deep zone in which the flow field isshown to be dominated by vortices connected with vonKármán vortex streets, butperiodically interrupted by strong sweep events whose features areinfluenced by canopy density. The second zone, which is near thecanopy top, is a superposition of attached eddies andKelvin–Helmholtz waves produced by inflectional instability in themean longitudinal velocity profile. Here, the relative importanceof the mixing layer and attached eddies are shown to vary withcanopy density through a coefficient . We show that therelative enhancement of turbulent diffusivity over its surface-layer value near the canopy top depends on the magnitude of. In the uppermost zone, the flow follows the classicalsurface-layer similarity theory. Finally, we demonstrate that thecombination of this newly proposed length scale and first-orderclosure models can accurately reproduce measured mean velocity andReynolds stresses for a wide range of roughness densities. Withrecent advancement in remote sensing of canopy morphology, thismodel offers a promising physically based approach to connect theland surface and the atmosphere without resorting to empiricalmomentum roughness lengths.