Numerical Simulation of Drag Partition over Rough Surfaces

We present a numerical simulation of drag partition over rough surfaces. A computational fluid dynamics model is applied with high resolution to simulatingturbulent flows over arrays of roughness elements positioned on asmooth surface. The skin drag on the surface and the pressure drag on the roughnesselements are computed. The simulated drag partition compares well with wind-tunnelmeasurements and theoretical estimates for similar rough surfaces. This confirms that the computational approach offers an alternative to wind-tunnel and field experiments in studying drag and drag partition. The model is then applied to studying drag partition over rough surfaces with various roughness configurations. It is shown that drag partition depends not only on the magnitude of the roughness frontal area but also on the sizes and arrangement of roughness elements, because (1) the pressure drag coefficient is sensitive to roughness-element dimensions and (2) the arragement of roughness elements lead to different interferences of turbulent wakes. The impact ofthe latter factor is not insignificant.     

Splash-saltation trajectories of soil particles under wind-driven rain

It is usually recognized that relatively large amounts of soil particles cannot be transported by raindrop splashes under windless rain. However, the splash-saltation process can cause net transportation in the prevailing wind direction since variations in splash-saltation trajectory due to the wind are expected in wind-driven rain. Therefore, determining the combined effect of rain and wind on the process should enable improvement of the estimation of erosion for any given prediction technique. This paper presents experimental data on the effects of slope aspect, slope gradient, and horizontal wind velocity on the splash-saltation trajectories of soil particles under wind-driven rain. In a wind tunnel facility equipped with a rainfall simulator, the rains driven by horizontal wind velocities of 6, 10, and 14 m s⁻¹ were allowed to impact three agricultural soils packed into 20×55 cm soil pans placed at both windward and leeward slopes of 7%, 15%, and 20%. Splash-saltation trajectories were measured by trapping the splashed particles at distances downwind on a 7-m uniform slope segment in the upslope and downslope directions, respectively, for windward and leeward slopes. Exponential decay curves were fitted for the mass distribution of splash-saltation sediment as a function of travel distance, and the average splash-saltation trajectory was derived from the average value of the fitted functions. The results demonstrated that the average trajectory of a raindrop-induced and wind-driven soil particle was substantially affected by the wind shear velocity, and it had the greatest correlation (r=0.96 for all data) with the shear velocity; however, neither slope aspect nor slope gradient significantly predicted the splash-saltation trajectory. More significantly, a statistical analysis conducted with nonlinear regression model of C₁(u_*²/g) showed that average trajectory of splash saltation was approximately three times greater than that of typical saltating sand grain.     

A parameterisation for the threshold shear velocity to initiate deflation of dry and wet sediment

A new parameterisation for the threshold shear velocity to initiate deflation of dry and wet particles is presented. It is based on the balance of moments acting on particles at the instant of particle motion. The model hence includes a term for the aerodynamic forces, including the drag force, the lift force and the aerodynamic-moment force, and a term for the interparticle forces. The effect of gravitation is incorporated in both terms. Rather than using an implicit function for the effect of the aerodynamic forces as reported earlier in literature, a constant aerodynamic coefficient was introduced. From consideration of the van der Waals force between two particles, it was further shown that the effect of the interparticle cohesion force between two dry particles on the deflation threshold should be inversely proportional to the particle diameter squared. The interparticle force was further extended to include wet bonding forces. The latter were considered as the sum of capillary forces and adhesive forces. A model that expresses the capillary force as a function of particle diameter squared and the inverse of capillary potential was deduced from consideration of the well-known model of Fisher and the Young–Laplace equation. The adhesive force was assumed to be equal to tensile strength, and a function which is proportional to particle diameter squared and the inverse of the potential due to adhesive forces was derived. By combining the capillary-force model and the adhesive force model, the interparticle force due to wet bonding was simplified and written as a function of particle diameter squared and the inverse of matric potential. The latter was loglinearly related to the gravimetric moisture content, a relationship that is valid in the low-moisture content range that is important in the light of deflation of sediment by wind. By introducing a correction to force the relationship to converge to zero moisture content at oven dryness, the matric potential–moisture content relationship contained only one unknown model parameter, viz. moisture content at −1.5 MPa. Working out the model led to a rather simple parameterisation containing only three coefficients. Two parameters were incorporated in the term that applies to dry sediment and were determined by using experimental data as reported by Iversen and White [Sedimentology 29 (1982) 111]. The third parameter for the wet-sediment part of the model was determined from wind-tunnel experiments on prewetted sand and sandy loam aggregates. The model was validated using data from wind-tunnel experiments on the same but dry sediment, and on data obtained from simulations with the model of Chepil [Soil Sci. Soc. Am. Proc. 20 (1956) 288]. The experiments showed that soil aggregates should be treated as individual particles with a density equal to their bulk density. Furthermore, it was shown that the surface had to dry to a moisture content of about 75% of the moisture content at −1.5 MPa before deflation became sustained. The threshold shear velocities simulated with our model were found to be in good agreement with own observations and with simulations using Chepil's model.     

Interactive discovery of sequential patterns in time series of wind data

Wind speed and direction vary over space and time due to the interactions between different pressures and temperature gradients within the atmospheric layers. Near the earth’s surface, these interactions are modulated by topography and artificial structures. Hence, characterizing wind behaviour over large areas and long periods is a complex but essential task for various energy-related applications. In this study, we present a novel approach to discover wind patterns by integrating sequential pattern mining and interactive visualization techniques. The approach relies on the use of the Linear time Closed pattern Miner sequence algorithm in conjunction with a time sliding window that allows the discovery of all sequential patterns present in the data. These patterns are then visualized using integrated 2D and 3D coordinated multiple views and visually explored to gain insight into the characteristics of the wind from a spatial, temporal and attribute (type of wind pattern) point of view. This proposed approach is used to analyse 10 years of hourly wind speed and direction data for 29 weather stations in the Netherlands. The results show that there are 15 main sequential patterns in the data. The spatial task shows that weather stations located in the same region do not necessarily experience similar wind pattern. For within the selected time interval, similar wind patterns can be observed in different stations and in the same station at different times of occurrence. The attribute task discovered that the repetitive occurrences of chosen pattern indicate as regular wind behaviour at different weather stations that persisted continuously over time. The results of these tasks show that the proposed interactive discovery facilitates the understanding of wind dynamics in space and time.     

中国中部典型高山气象站风速的均一性检验及其变化特点

采用二相回归方法并结合台站历史沿革信息,在对中国中部典型高山站南岳和庐山1960-2017年平均风速资料进行均一性检验和订正的基础上,分析其变化特征及其与周边低海拔台站的差异,并利用NCEP/NCAR再分析风速资料对其差异进行验证。结果表明：南岳站平均风速序列存在一个由测风仪器变更而导致的非均一点,而庐山站不存在非均一点;南岳和庐山年及四季平均风速均显著高于周边台站,且高山站以春季和夏季风速最大,而低海拔台站各季节风速差异较小;近58 a高山站及周边低海拔台站的年及四季平均风速均呈显著的减小趋势,但高山站的减小速率显著高于低海拔台站;同区域NCEP/NCAR的1000 hPa和850 hPa平均风速变化的差异与高山站和低海拔台站的差异基本一致,说明中低空和地面风速的这种差异在中国中部地区具有一定的普遍性。     

冷空气强风在大型城市中的精细结构和形成机制

地面强风可直接造成人员伤亡和经济损失,严重影响出行安全、工农业生产等社会秩序。强风的发生与天气系统和复杂下垫面的共同作用密切相关,在城市区域尤为明显。受数值模拟技术和计算资源的限制,对实际天气条件下大范围城区的强风现象进行建筑物分辨率的大规模数值模拟研究仍是一个挑战。本研究利用中尺度气象模式嵌套流体计算动力模式的超高分辨率局地气象预报系统,对强冷空气过程造成广州市区的一次强风事件进行数值模拟,深入探讨强风的精细结构和形成机制。结果表明,伴随着强冷空气入侵,广州市区的平均风速和风场高频扰动均明显增强,且在城市冠层顶尤为明显,呈现区域不均匀的三维结构,数值模拟与地面观测相一致。较大范围的强风速和阵风主要出现在建筑物较为低矮的老城区上空,并持续影响下游河道等开阔区域。在高层建筑密集的新城区,虽然整体风速明显减弱,但能在平行风向的街道狭管和下游区域形成局地强风。特别是,超高层建筑群引起显著的垂直环流,导致强风扰动向下传播,造成最大风速达8 m s?1的地面局地强风,阵风指数接近2。上游建筑群引起的风场扰动呈现大尺度湍流结构,能沿着平均气流传播影响数公里之远的下游地区。当风场扰动经过广州塔等单体超高层建筑时,可在其两侧绕流区再次加强,形成局地强风。局地强风和阵风还出现在垂直于风向排列的沿江高层建筑群的侧边,与建筑屏风的阻挡效应和缺口溢出有关。研究结果促进认识城市强风的时空特征和物理机制,有助于提升城市气象的精细化预报水平。     

延庆-张家口地区复杂地形冬季山谷风特征分析

基于2016年12月—2017年2月和2017年12月—2018年2月两年冬季的近地面自动气象站逐时观测数据以及张家口探空数据分析延庆-张家口一带（包括张家口崇礼、赤城、海坨、小五台山区,延怀、怀涿、洋河、蔚县盆地以及北京延庆、昌平、怀柔部分平原地区）复杂地形的风场精细化时、空分布特征,揭示不同复杂地形下局地风场的时、空变化规律,加深对复杂地形动力、热力作用对近地面风场影响的认识,为冬季山区风场预报以及复杂地形数值模式改进提供参考。结果表明:晴朗小风天风持续性作为矢量平均风速和标量平均风速的比值,可以作为研究风场变化规律的重要参数。根据风持续性的日变化特征,可以将研究区域内所有站点分为10种类型,分别代表不同局地地形特征的影响,风持续与风向变化的相关也很强。研究区域主要有3种类型的地形风:斜坡风、峡谷风以及较大尺度的山区平原风。不同地形特征下的风场、风持续性存在明显不同的日变化特征,山风和谷风相互转化的时间也不同,山区最早,盆地次之,平原区最晚;山风时段持续时间较谷风时段长,风速小;晴朗小风天实测风反映了实际风场的特征,而排除环境背景风场,弱化地形动力作用后整个冬季的局地风作为理论山谷风,更能反映热力作用下的山谷风特征。      

双雷达风场反演拼图在登陆台风“莫兰蒂”(1614)强降水精细预报中的同化应用试验

利用福建龙岩、漳州、泉州新一代多普勒天气雷达和厦门海沧双偏振雷达探测资料,采用动态地球坐标系下双雷达三维风场反演与拼图技术,基于天气研究和预报模式（Weather Research and Forecasting,WRF）及其资料同化系统,对登陆台风“莫兰蒂”（1614）引起的2016年9月14—15日福建强降水过程进行了双雷达风场反演拼图资料检验及其三维变分同化对强降水精细预报影响的数值试验,结果发现:（1）动态地球坐标系下双雷达反演风场能合理反映实际风场分布状况,其误差相对较小。相较厦门翔安风廓线雷达及厦门探空秒级测风数据,反演风风向（风速）平均绝对误差分别为7.8°（2.6 m/s）及3.4°（1.1 m/s）;（2）反演风场水平方向稀疏化对同化及预报结果极为重要,过密的反演风场资料会给同化及预报结果带来负效果。文中采用18、6、2 km 3重嵌套,在3重嵌套区域均进行同化以及仅在2 km区域进行同化两种情况下,均表现为当反演风场资料水平分辨率提高到0.1°时,同化分析及预报的台风环流开始受到负影响;且当反演风场资料水平分辨率越高时,负效果越明显。敏感性试验结果显示,分辨率取0.2°时数值预报效果最好;（3）以美国国家环境预报中心全球预报系统（National Centers for Environmental Prediction/Global Forecast System,NCEP/GFS）0.5°×0.5°分析场为初值,基于3个不同起报时刻（2016年9月14日14时、20时及15日02时）（北京时,下同）模拟的福建省境内台风内核雨带和螺旋雨带逐时演变、台风路径与强度、逐时降水TS评分和空间相关差异显著,其中14日14时起报试验效果最好;而14日20时起报试验效果最差,这与该试验初始台风大风轴风速明显偏大有关;（4）在上述3个不同起报时刻试验基础上,分别增加双雷达反演风场资料的三维变分同化后,福建境内地面风场和台风内核雨带、螺旋雨带逐时分布、逐时降水TS评分和空间相关、台风环流结构以及U、V风垂直廓线分布均有明显改善,最大正影响时效可达24 h;但仅对1—6 h时效内台风路径有改善。      

河西走廊风速变化及风能资源研究

利用河西走廊地区25年风速气候资料和风塔周年资料,研究了该区域近地面风速及风能的演变和分布。结果表明:河西走廊绿洲内风速下降十分明显,而其它高山站风速比较稳定,没有明显的减少趋势;风速变化具有周期振荡特点;垂直风速差由近地面向上减小,高层风速极大(小)值滞后于低层;风速随高度按自然对数规律增大,风能距地面8 m层内随高度变化迅速;该地区4~12 m/s风居多,是风能的主要贡献者;2~4月风速最大,1,5月最小;该地区风能丰富,10~70 m层内年风能储量在2200 khW/m2以上。     

风浪频谱中的特征量

在三参量风浪频谱的基础上对谱参量进行深入地研究,给出了谱参量与风场要素、波场要素的关系,提供了依据风场要素、波场要素及实测波浪资料计算谱参量的方法。从而可以依据上述因素直接计算出三参量风浪频谱。此外还根据谱宽度的变化,描述了风浪频谱的成长方式,解释了传统的波陡、波龄关系中经验常数的不同选取所代表的物理背景。