野外近地表风沙流脉动特征分析

基于风沙流实时测量系统对流沙地表的输沙强度以及不同高度风速进行实时测量所获得的数据,分析野外实际风场中的风沙流脉动特征。结果表明:沙粒的存在削弱平均流场,但对风速的高阶统计矩和分布形式影响不大,不同高度、不同来流平均风速条件下,风沙流中风速仍近似符合高斯分布;瞬时输沙强度服从指数分布,滞后风速约1 s,脉动强烈;30 min输沙强度平均值和脉动标准差均随平均风速的增大而呈幂函数增加,但随风速脉动强度的变化呈先增大后减小的趋势。这意味着脉动风速对输沙率具有重要影响,需要在输沙率预测公式中予以考虑。     

不同沙源供给条件下砾石床面的风沙流结构与蚀积量变化风洞实验研究

通过对不同沙源供给条件下各种砾石床面的风沙流结构、床面风蚀及堆积沙量变化的风洞实验,结果表明,风沙流结构是判断戈壁风沙流饱和与不饱和的一个重要途径,不同的戈壁风沙流结构对床面输、阻沙特性具有不同的指示意义。近地表0～6 cm高度内的风沙流结构决定了床面的输、阻性质,而6 cm以上的风沙流结构反映了风力对沙物质的输送状况。沙源供给的丰富与否,决定了风沙流的饱和程度,以及风沙流在砾石床面产生的蚀积状况。同等风速条件下,饱和风沙流的输沙率是非饱和风沙流输沙率的2～8倍。在饱和风沙流情形下,床面过程总体以积沙为主,且随风力的增强,床面积沙量急剧增加。在不饱和风沙流情形下,砾石床面总体以风蚀和输送沙物质过程为主,风沙流结构在0～2 cm高度内反映出砾石床面具有明显的阻沙功能,在2～5 cm高度上出现最大输沙值。     

风沙流中沙粒拖曳力系数研究

 拖曳力系数是计算风沙流中沙粒受空气阻力的重要参数。考虑实际风沙流中沙粒浓度及沙粒形状,利用FLUENT软件首先计算了不同风速下距地表不同高度处两沙粒的拖曳力系数,给出了影响沙粒拖曳力系数的间距范围,然后计算了真实风沙流中不同高度沙粒拖曳力系数。结果表明,给定风速下拖曳力系数随距地面高度的增加先减小后增加,并将沙粒拖曳力系数拟合成距床面高度的函数,该函数与风速有关。     

几种典型戈壁床面风沙流特性比较

通过不同材料覆盖的戈壁床面风沙流特性风洞模拟实验,发现对于棱角状砾石戈壁床面,地表动力学粗糙度随风速的增加而增加;而卵石床面,地表动力学粗糙度随风速的增加呈减小趋势。戈壁床面风速随高度的分布同样满足对数规律,棱角状砾石床面对风速的减弱程度相对于卵石床面更趋于显著。沙粒与戈壁床面棱角状砾石发生碰撞时其起跳高度增大,引起含沙量随高度分布不再满足流沙地表的指数衰减规律,而呈“象鼻效应”,出现拐点。戈壁地表输沙率与风速服从幂函数关系,但其幂指数远大于流沙地表。输沙率与风速之间幂函数关系中幂指数的取值主要受控于地表粒度组成。     

脉动风场下跃移沙粒对风的响应时间

风沙流是大气边界层的典型气固两相流,输沙率对风速的响应滞后时间,被称为风沙流响应时间。通过野外测量发现风沙流响应时间与风速测量高度成正比,通过对瞬时风速进行信号分解并计算近地表输沙率与不同时间尺度风速信号之间的相关性,发现响应时间随风速频率的减小而增加。为进一步得到定量的规律,采用数值模拟的方法分析了周期来流风场中,风沙流响应时间随各种参数(周期\,湍流强度\,来流摩阻风速和测量高度等)的变化规律。模拟结果显示:响应时间与风速变化周期、来流摩阻风速成正比,与湍流强度和边界层高度成反比,输沙强度对跃移层外风速变化的响应时间沿高度先增加后减小。     

砾石级沙粒胶结体抗风蚀效益的实验研究

塔克拉玛干沙漠腹地部分垄间地表发育了一种由众多沙粒胶结而成的大颗粒物质,称为沙粒胶结体(sand cemented bodies,缩写为SCB),其直径达到粗沙级、极粗沙级和砾石级。为了研究其对地表风沙活动的影响,本研究以野外采集的砾石级沙粒胶结体(gravel-size sand cemented bodies,缩写为GSSCB)为实验材料,在净风和挟沙风条件下进行了GSSCB覆盖沙面的抗风蚀模拟实验。结果表明:床面的蚀积状态与来流条件、GSSCB覆盖度和风速均有关,净风时所有覆盖度床面均呈风蚀状态,挟沙风时随覆盖度和风速而变化,床面可呈3种状态-风蚀、蚀积平衡、风积;在风蚀状态时,床面风蚀率随覆盖度增大以指数形式降低,随风速增大而以多种函数形式增加,抗风蚀效率随覆盖度增大而逐渐增加,但不同覆盖度范围增加率不同;挟沙风条件下呈蚀积平衡状态时的床面覆盖度临界C_b值与风速大小有关,随风速增加呈幂函数形式增加;在挟沙风条件下,覆盖度大于C_b值时床面呈风积状态,积沙率与风速的关系较为复杂,80%覆盖度床面积沙率随风速增大呈对数形式增加,但40%覆盖度床面积沙率则随风速增加呈指数形式降低。可见,由于与砾石的物理性质相近,GSSCB覆盖确实具有与砾石相类似的抗风蚀效益,并且在一定覆盖度条件下还能捕获风沙流挟沙颗粒。因此,塔克拉玛干沙漠腹地丘间地天然发育的GSSCB对于地表蚀积过程具有重要影响。GSSCB可作为一种新型固沙技术进行开发。     

风沙流和净风场中瞬时水平风速廓线特征比较

为进一步认识风沙流和净风场中瞬时水平风速廓线特征的异同,在风洞中分别对风沙流和净风场中瞬时水平风速廓线进行了测量,风速采集时间间隔缩短至0.01s,分析了风沙流和净风场中瞬时水平风速、瞬时摩阻风速和瞬时空气动力学粗糙度的变化特征。结果表明：相同来流条件下,风沙流中水平风速脉动强度高于净风场,风沙流中瞬时摩阻风速、瞬时空气动力学粗糙度以及它们的脉动幅度均大于净风场;风沙流和净风场中瞬时摩阻风速概率密度分布均可以表示为正态分布,但其正态分布的特征值却存在一定差别;净风场中瞬时空气动力学粗糙度的概率密度分布表现出单调递减分布,而风沙流中瞬时空气动力学粗糙度的概率密度分布呈现出单峰分布。因此,在相同主流风速下风沙流和净风场中瞬时水平风速廓线特征有明显差别。     

兰新高铁烟墩风区戈壁近地表风沙流跃移质垂直分布特性

跃移质作为风沙流的主体,其近地表垂直分布规律是风沙物理学的重要研究内容,对防沙工程具有重要的指导意义。受研究条件与观测仪器限制,戈壁特别是极端大风区近地表风沙流结构特性研究较为薄弱。利用多梯度风蚀传感器与阶梯式集沙仪对兰新高铁烟墩风区戈壁近地表风沙流跃移质的垂直分布特性进行了观测研究。结果表明：兰新高铁烟墩风区戈壁沙粒发生跃移运动的2 m高临界风速达12 m·s^-1;戈壁近地表风沙流具有明显的阵性特征,沙粒跃移发生的时间比例在50%以下,与平均风速成正相关关系,与风速脉动强度无显著相关关系;2 m高阵风7级风速下,戈壁跃移沙粒主要集中于地表50 cm范围内,近地表风沙流结构呈"象鼻效应",跃移质最大质量通量出现在地表2.5~5 cm高度处,沙粒最大跃移高度可达2 m,且沙粒跃移高度随2 m高风速的增加呈指数规律递增。因此,兰新高铁烟墩风区2 m高阻沙栅栏不足以完全阻截戈壁风沙流,是造成烟墩风区兰新高铁轨道积沙的重要原因之一。     

偏西风作用下敦煌月牙泉金字塔沙山顶部风沙流初步观测研究

金字塔沙山丘体高大,形态复杂,坡面风沙动力过程实地观测困难,至今尚无坡面风沙流野外观测数据。通过新研制自动集沙仪,对偏西风作用下敦煌月牙泉北侧金字塔沙山顶部风沙流进行了实时观测,共观测8组近地表51.4 cm高度内输沙数据,各观测时段中2 m高度平均风速7.46～14.15 m·s^-1,各组观测时长18~95 min。结果表明：偏西风纵向气流作用下金字塔沙山风沙流结构与平坦地表一致,即输沙量随高度增加呈单一的指数规律递减,且风沙流主要在近地表30 cm高度内输移。风沙流结构特征值λ以大于1为主,且随风速增大而增大,表明金字塔沙山坡面过程主要表现为风蚀过程。偏西风作用下金字塔沙山风沙流中沙粒以细沙为主,平均粒径随高度增加而减小;沙山迎风坡沙粒从坡脚到沙脊线逐渐粗化,风蚀作用增强,细沙从而随风沙流搬运至沙山顶,对金字塔沙山多以细沙为主的格局形成起重要作用。     

基于高速摄影技术的风沙流输沙通量剖面的风洞实验测量

风沙研究者非常重视对输沙通量随高度变化特征的研究,并为寻找可靠的测量手段付出了不懈的努力。基于高速摄影技术获得的沙粒平均水平速度与沙粒数的垂直剖面,推导了较低风速下环境风洞内输沙通量的垂直剖面。结果表明：沙粒平均水平速度随高度呈幂函数增加,颗粒浓度随高度的算数平方根呈指数衰减。由颗粒平均水平速度剖面与浓度剖面的乘积可获得输沙通量剖面。所获得的输沙通量随高度变化曲线在距床面1~3 mm处均有一个明显的拐点,拐点上方输沙通量随高度呈指数衰减。在床面与拐点之间输沙通量没有明显的变化趋势,这可能是由于气流中颗粒间的碰撞以及颗粒与床面碰撞的影响。平均跃移高度和相对衰减系数是描述输沙通量随高度变化的两个重要参数,两者有着很好的相关性,表明了随着风速增加和沙粒粒径减小跃移颗粒可以达到更大的高度,随着风速减小与粒径增大,输沙通量迅速衰减。     

Experimental investigation of the concentration profile of a blowing sand cloud

Detailed wind tunnel tests were carried out to establish the mean downwind velocity and transport rate of different-sized loose dry sand at different free-stream wind velocities and heights, as well as to investigate the vertical variation in the concentration of blowing sand in a cloud. Particle dynamic analyzer (PDA) technology was used to measure the vertical variation in mean downwind velocity of a sand cloud in a wind tunnel. The results reveal that within the near-surface layer, the decay of blown sand flux with height can be expressed using an exponential function. In general, the mean downwind velocity increases with height and free-stream wind velocity, but decreases with grain size. The vertical variation in mean downwind velocity can be expressed by a power function. The concentration profile of sand within the saltation layer, calculated according to its flux profile and mean downwind profile, can be expressed using the exponential function: c_z=ae^−bz, where c_z is the blown sand concentration at height z, and a and bare parameters changing regularly with wind velocity and sand size. The concentration profiles are converted to rays of straight lines by plotting logarithmic concentration values against height. The slope of the straight lines, representing the relative decay rate of concentration with height, decreases with an increase in free-stream wind velocity and grain size, implying that more blown sand is transported to greater heights as grain size and wind speed increase.     

风沙流中不同粒径组沙粒的输沙量垂向分布实验研究

在非均匀沙床面上, 风沙流中不同粒径组沙粒的输沙量垂向分布, 是非均匀风沙运动研究的重点。研究首先通过风洞实验, 收集了风洞中垂线垂向输沙量分布沙样, 然后对集沙沙样进行了沙粒粒度分析实验, 实验分析结果得出了不同粒径组沙粒的输沙量垂向分布规律, 基于稳定平衡风沙跃移运动模型和本文实验结果, 最后数值模拟研究了不同粒径组沙粒输沙量垂向分布, 与沙粒起跳速度和角度之间的关系。本文研究结果得出, 在非均匀风沙流中, 粗粒径组沙粒垂向输沙量上部符合指数递减分布但近床面区偏离指数分布, 呈现为偏大型分布, 粗粒径组对应的沙粒起跳速度和角度分布均为指数函数; 细粒径组沙粒垂向输沙量在整 个高度上均符合指数递减规律, 细粒径组沙粒对应的起跳速度分布为指数函数, 起跳角度分布为高斯函数。沙粒的平均起跳速度, 在0.4u*~2.2u* 之间变化, 随着气流风速(u*) 和沙粒粒径的增加而减小。     

PIV技术及其在风沙边界层研究中的应用

为了考察粒子图像测速度技术在风沙环境风洞中的测量精度及在风沙边界层研究中的应用潜力,通过筛选适当的示踪颗粒,借助PIV测量系统重新测量了风沙环境风洞中的风廓线,并获得了风沙边界层内跃移沙粒的速度和浓度分布规律。实验结果表明：PIV测得的风速廓线与标准风速廓线仪所测结果相当吻合(R2≥0.99);沙粒跃移的平均水平速度和相对浓度（灰度）沿距离沙面高度分别呈幂函数和负指数分布,沙粒速度随高度和自由风速的增加而增大,相对浓度随高度的增加而快速衰减,风速越小衰减越快,风速越大衰减越慢,这一结果与前人的结论一致。PIV系统为将来能够进行更加精确的风沙运动微观机理研究提供了技术保证。     

Wind tunnel experimental investigation of sand velocity in aeolian sand transport

Sand velocity in aeolian sand transport was measured using the laser Doppler technique of PDPA (Phase Doppler Particle Analyzer) in a wind tunnel. The sand velocity profile, probability distribution of particle velocity, particle velocity fluctuation and particle turbulence were analyzed in detail. The experimental results verified that the sand horizontal velocity profile can be expressed by a logarithmic function above 0.01 m, while a deviation occurs below 0.01 m. The mean vertical velocity of grains generally ranges from − 0.2 m/s to 0.2 m/s, and is downward at the lower height, upward at the higher height. The probability distributions of the horizontal velocity of ascending and descending particles have a typical peak and are right-skewed at a height of 4 mm in the lower part of saltation layer. The vertical profile of the horizontal RMS velocity fluctuation of particles shows a single peak. The horizontal RMS velocity fluctuation of sand particles is generally larger than the vertical RMS velocity fluctuation. The RMS velocity fluctuations of grains in both horizontal and vertical directions increase with wind velocity. The particle turbulence intensity decreases with height. The present investigation is helpful in understanding the sand movement mechanism in windblown sand transport and also provides a reference for the study of blowing sand velocity.     

Dynamic processes of a simple linear dune—a study in the Taklimakan Sand Sea, China

Samples of dune sands, surveys of the morphology and field measurements of wind velocity and direction of a simple linear dune in Taklimakan Sand Sea show that the airflow and sand flux vary with the change of wind direction on the dune surface. Decrease of the airflow stress on the lee flank does not result in much decrease of the sand flux because of the low threshold shear velocities and the airflow conditions. There are no significant relations between the sand flux on the lee flank and the angle of incidence of the airflow. The low threshold shear velocities and the maintenance of the sand flux at the lee flank are the main mechanisms keeping the linear shape of the dunes. Measurements of the sand flux shows that it reaches a maximum on the crest of the dune. The grain size of the transported sands has some differences compared to that of the dune surface. The sands transported are finer than that on the dune surface, but better sorted under the influence of the medium to low wind activity. The field experiment results exhibit that it is possible for the dunes to be shaped as linear dunes during the processes of accumulation and elongation.     

The flux profile of a blowing sand cloud: a wind tunnel investigation

The flux profile of a blowing sand cloud, or the variation of blown sand flux with height, is the reflection of blown sand particles that move in different trajectories, and also the basis for checking drifting sand. Here we report the wind tunnel results of systematic tests of the flux profiles of different sized sands at different free-stream wind velocities. The results reveal that within the 60-cm near-surface layer, the decay of blown sand flux with height can be expressed by an exponential function: q_h=aexp(−h/b), where, q_h is the blown sand transport rate at height h, a and b are parameters that vary with wind velocity and sand size. The significance of coefficient a and b in the function is defined: a represents the transport rate in true creep and b implies the relative decay rate with height of the blown sand transport rate. The true creep fraction, the ratio of the sand transported on the surface (h=0) to the total transport varies widely, decreasing with both sand size and wind speed. The flux profiles are converted to straight lines by plotting sand transport rate, q_h, on a log-scale. The slope of the straight lines that represents the relative decay rate with height of sand transport rate decreases with an increase in free-stream wind velocity and sand grain size, implying that relatively more of the blown sand is transported to greater heights as grain size and wind speed increase. The average saltating height represented by the height where 50% of the cumulative flux percentage occurs increases with both wind speed and grain size, implying that saltation becomes more intense as grain size and/or wind velocity increase.     

The distribution of velocity and energy of saltating sand grains in a wind tunnel

Sand transport by wind is a special case of two-phase flow of gas and solids, with saltating grains accounting for about 75% of the transport rate. This form of flow is not only the main external agent moulding aeolian landforms but also the motive force responsible for transport, sorting and deposition of aeolian sediments. High-speed multiflash photography is an effective method of studying the distribution of velocity and energy of saltating grains within the boundary layer of wind tunnel. The experimental wind shear velocities were set at 0.63, 0.64, 0.74 and 0.81 ms⁻¹. The statistical study of the results showed that there is a power function relation between mean velocity and height of saltating grains. As the height is divided into 0.5-cm intervals, the sand grain velocities at various levels are consistent with the Pearson VII distribution pattern. The variations in kinetic energy and total energy of sand grains with height accord with the pulse peak modified with power term (Pulsepow) law; the maximum values occur at heights of 6 cm or so and tend to shift upward with increasing wind velocity.     

Height profile of the mean velocity of an aeolian saltating cloud: Wind tunnel measurements by Particle Image Velocimetry

The velocity of saltating particles is an important parameter in studying the aeolian sand movement. We used Particle Image Velocimetry to measure the variation with height of the mean particle velocity of a saltating cloud over a loose sand surface in a wind tunnel. The results suggest that both the horizontal and vertical particle velocities fit the Gaussian distribution well, and that the mean particle velocity of a saltating cloud varies with wind velocity, particle size and the height above bed. The mean horizontal velocity is mainly the result of acceleration by the wind and increases with an increase in friction wind velocity but decreases with an increase in grain size because greater wind velocity causes more acceleration and finer particles are more easily accelerated at a given wind velocity. It also increases with an increase in height by a power function, in agreement with previous results obtained by other methods such as the high-speed multi-flash photographic method and Particle Dynamics Analyzer (PDA), reflecting, first, the increase in wind velocity with height through the boundary layer, and second, the longer trajectory-particle path length increases with height and affords a longer time for acceleration by the wind. An empirical model relating the mean horizontal particle velocity and height, friction wind velocity as well as particle size is developed. The ratio of the mean horizontal particle velocity to the clean wind velocity at the same height increases with height but decreases with grain size. The magnitude of mean vertical velocity is much less (one or two orders less) compared with the mean horizontal velocity. The average movement in the vertical direction of a saltating cloud is upward (the mean vertical velocity is positive). Although the upward velocity of a saltating particle should decrease with height due to gravity the mean vertical (upward) velocity (the average of both ascending and descending particles) generally shows a tendency to increase with height. It seems that at higher elevations the data are more and more dominated by the ‘high-flyers’. The underlying mechanism for the mean vertical velocity distribution patterns needs to be clarified by further study.     

A Wind Tunnel Investigation of the Shear Stress with A Blowing Sand Cloud

In a blowing sand system,the wind provides the driving forces for the particle movement while the moving particles exert the opposite forces to the wind by extracting its momentum.The wind-sand interaction that can be characterized by shear stress and force exerted on the wind by moving particles results in the modification of wind profiles.Detailed wind pro-files re-adapted to blown sand movement are measured in a wind tunnel for different grain size populations and at differ-ent free-stream wind velocities.The shear stress with a blowing sand cloud and force exerted on the wind by moving par-ticles are calculated from the measured wind velocity profiles.The results suggest that the wind profiles with presence of blowing sand cloud assume convex-upward curves on the u(z)-ln(z) plot compared with the straight lines characterizing the velocity profiles of clean wind,and they can be better fitted by power function than log-linear function.The exponent of the power function ranging from 0.1 to 0.17 tends to increase with an increase in wind velocity but decrease with an increase in particle size.The force per unit volume exerted on the wind by blown sand drift that is calculated based on the empirical power functions for the wind velocity profiles is found to decrease with height.The particle-induced force makes the total shear stress with blowing sand cloud partitioned into air-borne stress that results from the wind velocity gradient and grain-borne stress that results from the upward or downward movement of particles.The air-borne stress in-creases with an increase in height,while the grain-borne stress decreases with an increase in height.The air-borne shear stress at the top of sand cloud layer increases with both wind velocity and grain size,implying that it increases with sand transport rate for a given grain size.The shear stress with a blowing sand cloud is also closely related to the sand transport rate.Both the total shear stress and grain-borne stress on the grain top is directly proportional to the squ