Velocity profile of a sand cloud blowing over a gravel surface

Particle dynamic analyzer (PDA) measurement technology was used to study the turbulent characteristics and the variation with height of the mean horizontal (in the downwind direction) and vertical (in the upward direction) particle velocity of a sand cloud blowing over a gravel surface. The results show that the mean horizontal particle velocity of the cloud increases with height, while the mean vertical velocity decreases with height. The variation of the mean horizontal velocity with height is, to some extent, similar to the wind profile that increases logarithmically with height in the turbulent boundary layer. The variation of the mean vertical velocity with height is much more complex than that of the mean horizontal velocity. The increase of the resultant mean velocity with height can be expressed by a modified power function. Particle turbulence in the downwind direction decreases with height, while that in the vertical direction is complex. For fine sands (0.2–0.3 mm and 0.3–0.4 mm), there is a tendency for the particle turbulence to increase with height. In the very near-surface layer (<4 mm), the movement of blown sand particles is very complex due to the rebound of particles on the bed and the interparticle collisions in the air. Wind starts to accelerate particle movement about 4 mm from the surface. The initial rebound on the bed and the interparticle collisions in the air have a profound effect on particle movement below that height, where particle concentration is very high and wind velocity is very low.     

稳态风沙流中瞬态输沙特征

风沙流中沙粒运动在来流风速不变时也会表现出非稳态特征。在风洞内利用粒子图像测速系统（PIV）测量了风沙运动的时间序列，并基于PIV测量技术提出风沙流中沙粒平均直径、数密度、平均水平速度和输沙通量等参数在某一时刻的计算方法，其中输沙通量的计算考虑沙粒大小垂向分布的影响。结果表明：来流风速不变时，沙粒平均直径、数密度、平均水平速度和输沙通量随时间具有明显的波动性；沙粒平均直径和平均水平速度的标准偏差一般随高度增加而增加，沙粒数密度和输沙通量标准偏差随高度增大而减小；这些参数的相对标准偏差均随高度增加而增大。     

Grain size and transport characteristics of non-uniform sand in aeolian saltation

Size frequency distributions of sediment particles in a wind tunnel containing a bed of non-uniform sand are investigated by re-interpreting existing experimental data using particle-size analysis. Each particle sample is classified into one of eight groups according to its size grading. The analysis reveals that the modal shape of the particle-size frequency distributions of the saltating sand at different elevations or longitudinal distances is similar to that of the mixed sand in the bed once the boundary layer is fully developed. The standard deviation of the grain-size frequency distribution increases with increasing elevation above the bed then stays constant, whereas its skewness decreases. The mean grain size decays exponentially with elevation. The aeolian sand mass flux is determined for each size grading at different vertical and horizontal measurement locations. The vertical profile of aeolian horizontal mass flux depends on the size grading. The distribution of the sand transport rate according to the mean grain size in each grading fits the normal distribution. A parameter w_i is defined to reflect the likelihood of saltation for sand particles of the i-th size grading, and the mean sand size corresponding to the maximum value of w_i is found to be 0.2 mm. In addition, wind velocity strongly influences the magnitudes of the particle-size distribution and the sand mass flux distribution in both vertical and longitudinal directions.     

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.     

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

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

风沙二相流运动特点的分析

本文以平坦沙漠上定常，充分发展的风沙二相流为例分析了颗粒的跃移运动，在这种运动中颗粒吸收气流的机械能和水平动量向下输送，并在落地时传递给床面。本文还用拟流体观点分析了固相运动的特点：在床面上有较大的滑移速度，有较大的垂向速度脉动，例题垂向的平均速度却为零，等等。     

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.     

Tracking of saltating sand trajectories over a flat surface embedded in an atmospheric boundary layer

Saltation is a major mechanism for the transport of soil particles. In the present study, we carried out wind tunnel tests to examine the saltating trajectories of two types of natural sand collected from a beach (diameter, d = 300–500 μm and 200–300 μm respectively) as well as sand from the Taklimakan desert (d = 100–125 μm) in an atmospheric boundary layer. Consecutive images of saltating particles were recorded using a high-speed digital camera at a rate of 2000 fps with a spatial resolution of 1024 × 1024 pixels. The high temporal resolution of the acquired images enabled us to study the particle motion very close to the surface. The saltating particle trajectories were reconstructed from consecutive images, and the physical quantities characterizing the initial and final stages of the particle flight in the windward direction at friction velocities of about 10%–25% above the threshold friction velocity (u / u_t = 1.11–1.26) were analyzed statistically. In addition, the transverse deviation of the saltating particles from the main streamwise direction was evaluated. The results shed new light on the complicated motions involved in sand saltation and should prove useful in the evaluation and formulation of theoretical models.     

基于起跳初速度分布的沙颗粒浓度廓线的数值模拟

跃移层内沙颗粒浓度分布是风沙两相流相互作用的结果,准确的沙颗粒浓度分布有助于弄清风沙互馈机制及沙颗粒间相互作用机制。由于沙颗粒浓度分布与沙颗粒起跳初速度分布以及气流运动密切相关,本文基于特定的沙颗粒起跳初速度分布函数,通过构建的沙颗粒在气流中运动的物理模型,并利用四阶精度的Adams-Bashforth-Moulton方法对所构建运动模型进行求解,统计分析稳定状态下两相流中沙颗粒运动轨迹的分布,分析其浓度廓线的垂向分布规律。计算结果表明跃移层内沙颗粒浓度分布廓线与高程呈负指数或伽马分布关系;高度一定时沙颗粒浓度廓线随摩阻风速的增大而减小,随颗粒直径的增大而增大。     

Numerical simulation of aeolian dusty sand transport in a marginal desert region at the early entrainment stage

Aeolian dusty sand transport in the marginal region of a desert is described numerically from first suspension motion at the early entrainment stage to the unsteady state within a moderate range (1000 m long × 500 m high). A two-dimensional model is built for wind-blown dust flow, and the calculations are carried out using Fluent software. The simulation results describe an integrated picture of aeolian dusty sand transport including uplift, suspension, diffusion, deposition and its space–time concentration. According to the features of a sand–dust storm near the surface, a volume concentration expression of ejection grains is developed as a boundary condition in the simulation approach. The model is verified by comparing results with both experimental data from a wind tunnel and an analytic solution. Uniform dust sizes and R–R distributions are used in the simulation. Gas–solid two-phase flow patterns with these grains are obtained in the downwind space, including the turbulence intensity, gas phase stream functions and solid volume concentration distributions. The influence of wind velocities and grain sizes is analyzed. From the simulation results, spatial distributions of dust volume concentration in the early entrainment stage are described clearly. Different from coarse sands, there is a clear band of uniformly saturated dust concentration in the region directly above the surface.     

Vertical profiles of aeolian sand mass flux

Vertical profiles of the horizontal mass flux of blown sand are investigated experimentally using a passive vertical array in a wind tunnel. Considering lower sampling efficiency of the sand trap in the near-bed region, this investigation is complemented by the measurements of the longitudinal profiles of mass flux made using a horizontal sand trap. The experiments were conducted with two test sands and five different stream velocities.In the upper part of the vertical profile, the measured data exhibit an exponential decay distribution with a positive deviation occurring in the near-bed region. The measured longitudinal profiles are similar to the measured vertical profiles. Linking both profiles and the modes of sand transport, it is possible that saltating sand grains give rise to the well-known exponential decay distribution of mass flux, and that creeping and reptating grains force a deviation from it. A simple equation applicable for both the vertical and the longitudinal sand mass flux variations is introduced and the parameters are estimated from experimental data.     

风沙流中跃移沙粒冲击速度和角度联合概率分布风洞实验

近地表跃移沙粒速度和角度概率分布是连接风沙运动宏观研究和微观研究的桥梁，对沙粒跃移轨迹和输沙率分布有很大影响。但是，目前跃移沙粒冲击速度概率分布和角度概率分布之间的关系还不明确。使用相位多普勒粒子分析仪（PDPA）测量了风洞沙床面1 mm高度处的跃移沙粒速度，通过分析各冲击速度对应的冲击角度的概率分布研究了沙粒冲击速度和角度的联合概率分布。结果表明：沙粒冲击速度概率分布可用对数正态函数描述，冲击角度概率分布可用指数函数描述；各冲击速度对应的冲击角度都符合指数函数分布，且其衰减速度随冲击速度增加而加快。这表明，冲击速度概率分布和冲击角度概率分布之间具有很强的相关性。最后分别给出了沙粒冲击速度和角度的独立假设概率分布和条件联合概率分布。     

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

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

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

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

紊流风场中起跳沙粒的轨迹特征

运用颗粒流体二相流的随机行走扩散模型(DRW)研究了紊流风场下起跳沙粒的运动轨迹特性。得出了轨迹特性(长度、高度、撞击角)随沙粒的大小和风速的变化规律,发现随沙粒粒径的减小紊流脉动对轨迹的影响显著增大。粒径在200μm以下时,紊流风场下计算得的轨迹形状与平均风场下得到的轨迹有明显差异。     

广西跨海大桥桥位区风场特征实测分析

为了全面准确地掌握广西跨海大桥的风环境和风特性参数,2012 年在桥位区设立了一座高度为 55 m 的测风塔,安装了梯度风和超声风观测仪器对其近地风场环境进行了连续 1 a 的现场观测,分析大桥桥位区实测的平均风场特征和湍流脉动风场特征,并推算了桥位区设计风速,结果表明：1）桥位区平均风速垂直变化曲线符合幂指数规律,大风状况时的风切变指数为 0.062,平均风速上半年＞下半年,累年平均和秋、冬、春季主导风向较稳定且均为偏北风;2）在台风“启德”过程观测到的风攻角基本不超过规范用于抗风性能检验规定的±3°范围,脉动风湍流表现为均匀平稳的湍流变化过程,三维方向的湍流强度变化趋势基本一致、湍流积分尺度随时间的变化不大,湍流功率谱在惯性副区内基本符合 Kolmogrove 谱-5/3 律;3）大桥桥面高度处的设计风速 50 a 一遇为 39.7 m/s,100 a 一遇为 42.7 m/s。     

新月形沙丘表面100cm高度内风沙流输沙量垂直分布函数分段拟合

为研究新月形沙丘表面不同层位风沙流输沙量的垂直分布函数,实测了塔克拉玛干沙漠腹地典型新月形沙丘表面100 cm高度内(以1 cm分隔)的输沙量。分段拟合分析表明:新月形沙丘迎风坡脚输沙量垂直分布规律不完全服从指数函数,出现与戈壁风沙流结构特征相似的"象鼻效应",在0~3 cm区间内输沙量逐渐增大,3cm以上输沙量随高度呈指数函数衰减;沙丘顶部0~10 cm区间输沙量随高度呈指数函数衰减,10 cm以上呈二次函数衰减;沙丘左翼端输沙量随高度呈幂函数分布,沙丘右翼端0~20 cm内以指数函数衰减,20 cm以上呈三次函数衰减;沙丘背风坡脚风沙流输沙量在0~60 cm和60 cm以上分别呈不同形式的三次函数分布。     

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.     

风沙流中风速纵向脉动的实验研究

在4种风速条件下对风沙流中风速纵向脉动特性进行了风洞实验研究。研究发现，风速纵向脉动很大．不容忽视；脉动绝对大小随进口风速提高有增强趋势，随着高度的增加则变化不大；垂线脉动强度平均值随主流速度基本成正向线性关系；以脉动强度与当地风速之比表示的相对脉动强度随高度增加而降低，并且不同风速条件下相对脉动强度在床面附近相差较大，远离床面则趋向一致。另外，采用幂函数形式可较好地描述风沙流中平均风速的垂线分布。     

稳态风沙流中沙粒体积浓度的模型分析

本文建立了一个简单有效的稳态风沙流中沙粒体积浓度的数学模型,包括3个部分:稳态风沙流的风速廓线描述、跃移沙粒轨迹的计算、床面沙粒起跳速度分布的描述。利用已知实验的风速廓线作为气流速度场输入参数,跃移沙粒轨迹的计算主要考虑重力和拖曳力,基于体积观点的床面起跳沙粒的水平速度和垂直速度概率密度分布函数分别采用正态分布和指数分布函数来描述。根据稳态风沙流中运动沙粒的动态平衡特征可推导计算沙粒体积浓度。计算的沙粒体积浓度与风洞实验结果基本一致,说明本模型有较好的预测能力。