粗糙元几何参数的交互作用对床面空气动力学粗糙度的影响

 为了全面地揭示粗糙元的所有几何参数的交互作用对空气动力学粗糙度的影响,利用风洞实验研究了粗糙元高度、密度、高度与间距比、孔隙度和方向比率等几何参数交互作用对空气动力学粗糙度的影响。结果表明,密实和孔隙粗糙元的无量纲空气动力学粗糙度（空气动力学粗糙度度/粗糙元高度）均可表示为粗糙元密度/等效密度的正比例函数,而比例系数反映了粗糙元几何参数交互作用。据此,该研究发展了一个全面反映粗糙元高度h、密度/等效密度λ、高度与行间距比Sp和方向比率AR等几何参数交互作用的空气动力学粗糙度模式:Z0h=-0.0028+0.5403S0.32p·AR-0.07·λ。该模式改进了模拟的精度,扩大了适用范围。     

The Effect of Stratification on the Aerodynamic Roughness Length and Displacement Height

The roughness length, z _0u , and displacement height, d _0u , characterise the resistance exerted by the roughness elements on turbulent flows and provide a conventional boundary condition for a wide range of turbulent-flow problems. Classical laboratory experiments and theories treat z _0u and d _0u as geometric parameters independent of the characteristics of the flow. In this paper, we demonstrate essential stability dependences—stronger for the roughness length (especially in stable stratification) and weaker but still pronounced for the displacement height. We develop a scaling-analysis model for these dependences and verify it against experimental data.     

Temporal and Spatial Variations of the Aerodynamic Roughness Length in the Ablation Zone of the Greenland Ice Sheet

To understand the response of the Greenland ice sheet to climate change the so-called ablation zone is of particular importance, since it accommodates the yearly net surface ice loss. In numerical models and for data analysis, the bulk aerodynamic method is often used to calculate the turbulent surface fluxes, for which the aerodynamic roughness length (z ₀) is a key parameter. We present, for the first time, spatial and temporal variations of z ₀ in the ablation area of the Greenland ice sheet using year-round data from three automatic weather stations and one eddy-correlation mast. The temporal variation of z ₀ is found to be very high in the lower ablation area (factor 500) with, at the end of the summer melt, a maximum in spatial variation for the whole ablation area of a factor 1000. The variation in time matches the onset of the accumulation and ablation season as recovered by sonic height rangers. During winter, snow accumulation and redistribution by snow drift lead to a uniform value of z ₀≈ 10⁻⁴ m throughout the ablation area. At the beginning of summer, snow melt uncovers ice hummocks and z ₀ quickly increases well above 10⁻² m in the lower ablation area. At the end of summer melt, hummocky ice dominates the surface with z ₀ > 5  ×  10⁻³ m up to 60 km from the ice edge. At the same time, the area close to the equilibrium line (about 90 km from the ice edge) remains very smooth with z ₀ = 10⁻⁵ m. At the beginning of winter, we observed that single snow events have the potential to lower z ₀ for a very rough ice surface by a factor of 20 to 50. The total surface drag of the abundant small-scale ice hummocks apparently dominates over the less frequent large domes and deep gullies. The latter results are verified by studying the individual drag contributions of hummocks and domes with a drag partition model.     

The Parameterisation of Scalar Transfer over Rough Ice

We test a surface renewal model that is widely used over snow and ice surfaces to calculate the scalar roughness length (z _s ), one of the key parameters in the bulk aerodynamic method. For the first time, the model is tested against observations that cover a wide range of aerodynamic roughness lengths (z ₀). During the experiments, performed in the ablation areas of the Greenland ice sheet and the Vatnajökull ice cap in Iceland, the surface varied from smooth snow to very rough hummocky ice. Over relatively smooth snow and ice with z ₀ below a threshold value of approximately 10^?3 m, the model performs well and in accord with earlier studies. However, with growing hummock size, z ₀ increases well above the threshold and the bulk aerodynamic flux becomes significantly smaller than the eddy-correlation flux (e.g. for z ₀ = 0.01 m, the bulk aerodynamic flux is about 50% smaller). Apparently, the model severely underpredicts z _s over hummocky ice. We argue that the surface renewal model does not account for the deep inhomogeneous roughness sublayer (RSL) that is generated by the hummocks. As a consequence, the homogeneous substrate ice grain cover becomes more efficiently ‘ventilated’. Calculations with an alternative model that includes the RSL and was adapted for use over hummocky ice, qualitatively confirms our observations. We suggest that, whenever exceedance of the threshold occurs (z ₀  >  10^?3 m, i.e., an ice surface covered with at least 0.3-m high hummocks), the following relation should be used to calculate scalar roughness lengths, ln (z _s /z ₀)  =  1.5  ? 0.2 ln (Re _*)  ? 0.11(ln (Re _*))².     

A straightforward approach for using single time domain simulations to assess characteristic extreme responses

Short-term wave design approach of marine structures, using nonlinear time domain simulations, is a design procedure that is recognized by various modern standard codes. One of the most challenging points of this approach is the evaluation of the characteristic extreme values for response parameters used in the design check equations. The most straightforward and recommended way to evaluate a response characteristic value is by fitting an extreme value probability distribution to the N-sample of extreme values extracted from N independent time domain simulations with duration equal to the short-term period indicated by the code, which is usually taken as 3 h. However, this procedure would not be practical for some types of marine structures, such as risers and mooring lines, under numerous design load cases and demanding huge finite element models. A more feasible approach would be to assess the response extreme value distribution using only a single short-term time domain simulation with duration shorter than 3 h. But reduced time simulations always introduce some additional statistical uncertainty into the extreme values estimates. This paper discusses a workable way of properly taking into account the statistical uncertainty associated with the simulation length in the assessment of a characteristic short-term extreme response value based on a single time series.     

独立坐标系模型对长度变形影响的研究

针对海拔较高的城市建立或改造独立坐标系,存在着如何选择模型及各模型对长度变形的影响等问题,通过算例对常见椭球变换和比例缩放模型进行较为全面系统的分析,归纳出常见模型在不同投影带区域、不同投影面高度对长度变形的影响差异,并得出有益的结论。     

城市下垫面空气动力学参数的确定

分析近10年北京325m气象塔常年观测资料,研究城市化发展对城市大气边界层动力学结构和特征的影响.结果表明,空气动力学粗糙度和零平面位移总体是逐年递增的.气象塔西南方向的高层建筑群是20世纪90年代初期逐渐建成,在SW方向,空气动力学粗糙度和零平面位移在90年代增加明显,增加的幅度1997年和1999年之间较大,这一点与气象塔西南方向城市化的加快相符合.而在NE、SE和NW方向,从1987年到1994年空气动力学粗糙度和零平面位移增加明显,1994年以后变化较小.     

Influence of length effect on embankment slope reliability in 3D

Embankment slopes composed of spatially variable soils have a variety of different failure modes that are affected by the correlation distances of the material properties and the geometry and total length of the slope. This paper examines the reliability of soil slopes for embankments of different length and uses parallel computing to analyse very long embankments (up to 100 times the embankment height) for a clay soil characterised by a spatially varying undrained shear strength. Based on a series of analyses using the 3D random finite element method (RFEM), it is first shown that the reliability of slopes of various length can be efficiently computed by combining simple probability theory with a detailed 3D RFEM analysis of a representative shorter slope of length 10 times the slope height. RFEM predictions of reliability indices for longer slopes are then compared with results obtained using Vanmarcke's (1977a) simplified 3D method and Calle's (1985) extended 2D approach. It is shown that these methods can give significantly different results, depending on the horizontal scale of fluctuation relative to the slope length, with RFEM predicting a lower slope reliability than the Vanmarcke and Calle solutions in all cases. The differences in the solutions are evaluated and attributed to differences in the assumed and computed failure surface geometries.     

城市坐标系的建立与申报

为了满足城市建设对测绘基准的要求,需建立相对独立的平面坐标系,我国测绘法等相关法律法规对独立平面坐标系的建立及管理有着明确的要求,本文结合新疆某县城相对独立平面坐标系的建立技术方案、申报材料准备、申报组织以及发布等方面进行了详细论述,希望能为测绘同行和测绘管理单位建立城市相对独立平面坐标系起到一定的参考作用。     

Estimation of hourly and daily evapotranspiration and soil moisture using downscaled LST over various urban surfaces

Surface moisture is important to link land surface temperature (LST) to people’s thermal comfort. In urban areas, the surface roughness from buildings and urban trees impacts wind speed, and consequently surface moisture. To find the role of surface roughness in surface moisture estimation, we developed methods to estimate daily and hourly evapotranspiration (ET) and soil moisture, based on a case study of Indianapolis, Indiana, USA. In order to capture the spatial and temporal variations of LST, hourly and daily LST was produced by downscaling techniques. Given the heterogeneity in urban areas, fractions of vegetation, soil, and impervious surfaces were calculated. To describe the urban morphology, surface roughness parameters were calculated from digital elevation model (DEM), digital surface model (DSM), and Terrestrial Light Detection and Ranging (LiDAR). Two source energy balance (TSEB) model was employed to generate ET, and the temperature vegetation index (TVX) method was used to calculate soil moisture. Stable hourly soil moisture fluctuated from 15% to 20%, and daily soil moisture increased due to precipitation and decreased due to seasonal temperature change. ET over soil, vegetation, and impervious surface in the urban areas yielded different patterns in response to precipitation. The surface roughness from high-rise has bigger influence on ET in central urban areas.