风沙跃移运动发展过程及静电力影响的数值模拟

建立了描述风沙跃移运动发展过程以及风场一沙粒互馈机制的风沙跃移运动数值模型,模型中考虑了风场一沙粒之间酌耦合作用以及静电力的影响,同时风场是随时间变化的．模拟了风沙跃移运动从起始阶段至风沙流达到自平衡状态的整个发展过程所需的时间,起沙率与单宽输沙率随时间变化曲线,以及达到稳定后的单位面积输沙率沿高度分布以及贴地风速廓线的影响．计算结果表明,当沙粒平均带电量为60μC／kg时,计算所得的单宽输沙率以及输沙率沿高度分布与实验结果吻合得较好．     

风沙两相流跃移层中沙粒相的速度分布

从单个跃移沙粒在气流中的运动方程出发导出了风沙两相流中沙粒相速度分布函数的Boltzmann方程，对风沙流研究中几种不同的分布函数及其相应的统计平均值等基本概念给出了严密的数学定义，指出了不同分布函数之间的区别和联系，在略去铅垂方向空气阻力的情况下，给出了沙粒相速度分布函数沿铅垂方向的边缘分布，作为风沙流中跃移理论的主要基础之一。利用结果对前人在风沙流研究中发现的某些重要规律和现象进行了解释。     

可侵蚀地表上方含尘大气运动的数值模型

在稀相气固两相流的双流体模型框架下，导出含尘大气运动的基本方程，其中的源项考虑相间双向耦合作用；通过对基本方程进行无量纲化，求得控制气相和弥散相动力学行为的相似准则，包括弗洛德数、罗斯贝数、颗粒雷诺数、颗粒惯性参数、颗粒剪切横移参数、颗粒旋转横移参数和沙尘质量载荷率等. 作为一个模型问题，研究可侵蚀地表上方充分发展湍流大气边界层流动对土壤颗粒的气动卷扬作用，其中计及气动阻力、萨弗曼升力和重力等因素的影响. 为了克服由于沙尘跃移运动轨迹交叉导致流动参数多值性所造成的困难，引进拉格朗日坐标下弥散相连续方程. 对两种不同风速和4种不同粒径条件下的沙尘运动进行数值模拟，给出相应情况中沙尘运动特性和密度分布剖面并讨论风速和粒径等参数的影响，还细致分析了含尘大气边界层中沙尘与气流之间的能量传递过程. 所得的方程、准则和方法可以为土壤风蚀和沙尘暴等自然灾害的预报提供理论基础.     

关于跃移沙破坏生物结皮时的临界速度的理论分析

利用经典的弹性薄板与塑性理论探讨了生物结皮在单颗跃移沙粒碰撞作用下的破坏问 题,建立了一个联结跃移沙粒的临界撞击速度和结皮性质的解析关系式. 种类、数量、发展 阶段等生物效应主要通过结皮的杨氏模量、Possion比、屈服应力等宏观参量引入. 结果表 明,刚度大、强度低的生物结皮更容易被跃移沙粒所破坏. 给出的关系式将有助于进一 步研究风沙运动和结皮之间的相互作用.     

风沙运动的沙粒带电机理及其影响的研究进展

土壤风蚀以及沙尘暴是人类面临的重大自然环境灾害之一,它直接关系着人类的生存环境和生活质量.许多国家和政府以及科学家们对此一直十分关注并开展了大量研究.由于大气和地表条件的多样性以及地表沙粒的运动与风场作用的相互耦合,导致风沙运动问题研究的复杂性.除此之外,伴随着对风沙运动机理的研究和对沙尘暴的观测,人们发现在风沙跃移运动和沙尘暴中的沙粒带有电荷并形成风沙电场.沙粒的运动不仅导致沙粒带电和风沙电场形成,而且也直接受到沙粒带电和风沙电场的影响.为了弄清风沙电的起因以及揭示其影响规律,研究人员对此进行了野外观测、实验室测量和理论推测等.本文将就这一问题重点介绍有关风沙运动中沙粒带电的机理、风沙电场分布规律、沙粒带电对风沙运动的微宏观物理量以及沙尘暴中沙粒带电对电磁波散射的影响等方面的实验和理论研究的基本概况和进展以及作者等人在此方面的研究工作.主要包含内容有:沙粒带电及其电场的基本实验测量,粒带电对风沙运动的影响,沙尘暴对电磁波衰减的影响研究以及有关风沙起电的机理研究等.同时,对这一研究领域中目前关注与存在的主要问题给予评述.图6参81      

非定常阵风场对沙粒跃移运动的影响浅析

野外阵风具有大尺度结构的非定常特征和小尺度结构的高频脉动特性,是造成风沙流输沙量波动和沙粒间歇性跃移的主要原因. 阵风的非定常性和脉动越显著,采用平均风速及沙粒粒径均一化分析风沙问题的误差越大,是野外观测、风洞实验及数值模拟的结果彼此偏差极大的主要因素. 野外阵风条件下的风沙跃移研究目前仍缺乏能够用于规范数学模型的有关阵风场和输沙率特征的可靠实时测量,有赖于风沙流整场高频测量技术和风洞模拟技术的发展.     

波浪场中物质输移离散系数的理论分析

海岸水域波浪引起的物质输移扩散存在着不同于水流引起的物质输移扩散的特征. 通过解析的方法对垂向扩散方程进行求解, 推导出波浪和潮流共同作用下的海岸水域物质输移离散系数的理论公式. 分析中将总的离散系数分为潮流、质量输移流、纯波浪波动作用以及潮流与质量输移流、潮流与波浪波动、质量输移流与波浪波动的相互作用6部分, 通过分析各组成部分特征, 得出海岸水域物质输移中波浪所产生的离散效应主要是质量输移流贡献, 波浪水质点运动及波流相互作用对时间平均离散系数的贡献较小的结论.仅波浪的作用(不包括潮流)的有关结果与他人数值模拟和实验结果的吻合良好, 证明了该理论推导的正确性. 给出了时间平均离散系数及离散系数波动幅值随波浪周期和波高的变化特征,同时将结果应用于垂向环流, 给出了垂向环流对应于质量输移流所产生的离散系数的特征.      

考虑系留点垂向运动的飞艇系缆响应特性分析

为研究考虑系留点垂向运动的飞艇系缆的面内运动响应特性,将运动的飞艇作为系缆的边界条件处理,对系缆下端的系留点受到垂向激励时的系留飞艇系缆的运动状态进行了数值模拟,研究了系缆的动张力、运动情况。分析了飞艇攻角、系缆线密度等参数对系统的影响。研究结果表明:当系留点竖向简谐运动时,系缆上各点都做近似的简谐运动,在不同的线密度和攻角参数下,系缆各点的法向振幅分布不同,系缆可能处于二阶的振动状态;且随着系缆密度和攻角的增大,切向振幅增大,法向振幅减小;系留点激励还引起整个系统偏离平衡位置,这会引起飞艇动态高度的降低,在系留飞艇的高度控制中会有一定的影响。以上研究结果可为系留在舰船等运动中平台上的飞行器控制和系缆设计提供参考。     

山地下击暴流风剖面特性试验研究

为研究山地下击暴流的近地风场特征,建立了三维余弦山体模型,通过物理试验研究了不同初始出流风速和径向距离下山地下击暴流的风场特性、风剖面特性和山顶风速加速效应。研究结果表明:在山体迎风侧,气流挤压使得风剖面中最大风速值所在高度更贴近地面,在山体背风侧,由于山顶对气流的抬升作用,风剖面中最大风速值所在高度显著提高;距离风暴中心径向距离和出流风速大小会明显影响山地下击暴流的风速值,对山顶最大风速加速比影响较小;山地下击暴流山顶最大风速加速比在数值上近似等于1.0加上对应的山体坡度值。     

一种改进的高浓度风沙图像的动态灰度阈值分割算法

风沙跃移是塑造干旱区地貌的最主要的动力,其伴生的粉尘的释放、输送和沉降不但严重地影响了大气环境质量,也引起了全球气候 系统和海洋系统的变化. 单个沙粒轨迹形成是理解粒--床碰撞过程和风沙两相流耦合过程的重要纽带. 高浓度风沙图像处理算法是深入理解沙粒轨迹形成机制的关键技术. 为了实现对高浓度风沙图像精细处理, 本文提出了不依赖于经验参数的动态灰度阈值分割算法,它包括背景模板去噪、绿光通道灰度化处理、图像微分、灰度方差阈值目标检测和 最大类间方差灰度阈值分割等,其中背景模板去噪和灰度方差阈值目标检测等是新算法的主要亮点. 高浓度风沙图像分割实验显示, 扣减背景模板去除了条纹状和斑点状的稳定噪声;图像微分和灰度方差阈值目标检测显著地提高了暗沙粒识别的数量并有效地去除了随机噪声; 改进算法的沙粒有效识别个数(人 机判读坐标一致的沙粒数目)、查全率(计算机提取沙粒数目与人工判读的实际沙粒数比值)和查准率(有效沙粒数目与实际沙粒数目的比值) 分别为461, 71%和86%,显著地高于传统算法对应的85, 13%和82%,这表明新算法对高浓度风沙图像的分割效果良好. 避免单粒子分割和表观重叠问题是进一步完善高浓度风沙图像的分割算法的可能途径.      

Experimental analysis of sand particles＇ lift-off and incident velocities in wind-blown sand flux

The probability distributions of sand particles' lift-off and incident velocities in a wind-blown sand flux play very important roles in the simulation of the wind-blown sand movement. In this paper, the vertical and the horizontal speeds of sand particles located at 1.0 mm above a sand-bed in a wind-blown sand flux are observed with the aid of Phase Doppler Anemometry (PDA) in a wind tunnel. Based on the experimental data, the probability distributions of not only the vertical lift-off speed but also the lift-off velocity as well as its horizontal component and the incident velocity as well as its vertical and horizontal components can be obtained by the equal distance histogram method. It is found, according to the results of the χ ²-test for these probability distributions, that the probability density functions (pdf's) of the sand particles' lift-off and incident velocities as well as their vertical components are described by the Gamma density function with different peak values and shapes and the downwind incident and lift-off horizontal speeds, respectively, can be described by the lognormal and the Gamma density functions. These pdf's depend on not only the sand particle diameter but also the wind speed.     

Experimental analysis of sand particles'' lift-off and incident velocities in wind-blown sand flux

The probability distributions of sand particles' lift-off and incident velocities in a wind-blown sand flux play very important roles in the simulation of the wind-blown sand movement. In this paper, the vertical and the horizontal speeds of sand particles located at 1.0 mm above a sand-bed in a wind-blown sand flux are observed with the aid of Phase Doppler Anemometry (PDA) in a wind tunnel. Based on the experimental data, the probability distributions of not only the vertical lift-off speed but also the lift-off velocity as well as its horizontal component and the incident velocity as well as its vertical and horizontal components can be obtained by the equal distance histogram method. It is found, according to the results of the χ²-test for these probability distributions, that the probability density functions (pdf's) of the sand particles' lift-off and incident velocities as well as their vertical components are described by the Gamma density function with different peak values and shapes and the downwind incident and lift-off horizontal speeds, respectively, can be described by the lognormal and the Gamma density functions. These pdf's depend on not only the sand particle diameter but also the wind speed. The project supported by the National Natural Science Foundation of China (10532040) and the Hundred Talents Project, the Knowledge Innovation Project of Chinese Academy of Sciences (KZCX2-304). The English text was polished by Keren Wang.     

Analysis of sand particles’ lift-off and incident velocities in wind-blown sand flux

In the research of windblown sand movement,the lift-off and incident velocities of saltating sand particles play a significant role in bridging the spatial and temporal scales from single sand particle’s motion to windblown sand flux.In this paper,we achieved wind tunnel measurements of the movement of sand particles near sand bed through improving the wind tunnel experimental scheme of paticle image velocimetry(PIV) and data processing method.And then the influence of observation height on the probability distributions of lift-off and incident velocities of sand particles was analyzed.The results demonstrate that the observation height has no obvious influence on the distribution pattern of the lift-off and incident velocities of sand particles,i.e.,the probability distribution of horizontal and vertical velocities of lift-off and incident sand particles follow a Gaussian distribution and a negative exponential distribution,respectively.However,it influences the center of the Gaussian distribution,the decay constant and the amplitude of the negative exponential distribution.     

Effect of thermal diffusion and electrostatic force on evolution of wind-blown sand flow

A theoretical model is suggested to mathematically describe the effect of thermal diffusion from a sand-bed on evolution of a wind-blown sand flow.An upward wind field is engendered by the thermal diffusion and the coupling interaction among the horizontal and upward wind flow,saltating grains,and a kind of electrostatic force exerted on the grains are considered in this theoretical model.The numerical results show that the effect of the thermal diffusion on the evolution process of wind-blown grain flow is quite obvious and very similar to the effect of the electrostatic force on the evolution.Not only the time for the entire system to reach a steady state(called the duration time),the transport rate of grains,the mass-flux profiles and the trajectory of saltating grains are affected by the thermal diffusion and the electrostatic force exerted on saltating grains, but also the wind profiles and the temperature profiles at the steady state are affected by the wind-blown sand flow.     

风沙颗粒运动的数字高速摄影图像的分割算法

为了深入理解多因素驱动下风沙颗粒起动的动态演化规律, 需要准确地获得沙质床面附近沙粒群起跳的方式、速度和运动轨迹. 以连续强激光源为照明的数字高速摄影技术是研究这类问题的有效手段, 但由于风沙运动的高速摄影图像具有运动沙粒和静止床面的对比度小、相邻两帧图像相似性小等特点, 原始图像叠加算法难以有效实现目标与背景的分割. 该文提出了基于相邻的风沙运动图像灰度差值变化原理的图像分割算法. 实例显示,只要选择合适的相邻图像灰度差值阈值和自适应二值化处理方法就能实现图像中运动沙粒与床面分割. 当起跳沙粒浓度较低情况下, 基于MATLAB平台的最小距离匹配的粒子追踪算法(PTV算法)能较为准确地恢复床面附近沙粒的运动轨迹.      

Simultaneous PIV and PTV measurements of wind and sand particle velocities

Wind-blown sand is a typical example of two-phase particle-laden flows. Owing to lack of simultaneous measured data of the wind and wind-blown sand, interactions between them have not yet been fully understood. In this study, natural sand of 100–125 μm taken from Taklimakan Desert was tested at the freestream wind speed of 8.3 m/s in an atmospheric boundary layer wind tunnel. The captured flow images containing both saltating sand and small wind tracer particles, were separated by using a digital phase mask technique. The 2-D PIV (particle imaging velocimetry) and PTV (particle tracking velocimetry) techniques were employed to extract simultaneously the wind velocity field and the velocity field of dispersed sand particles, respectively. Comparison of the mean streamwise wind velocity profile and the turbulence statistics with and without sand transportation reveal a significant influence of sand movement on the wind field, especially in the dense saltating sand layer (y/δ < 0.1). The ensemble-averaged streamwise velocity profile of sand particles was also evaluated to investigate the velocity lag between the sand and the wind. This study would be helpful in improving the understanding of interactions between the wind and the wind-blown sand.     

Discrete particle simulation of mixed sand transport

An Eulerian/Lagrangian numerical simulation is performed on mixed sand transport. Volume averaged Navier–Stokes equations are solved to calculate gas motion, and particle motion is calculated using Newton's equation, involving a hard sphere model to describe particle-to-particle and particle-to-wall collisions. The influence of wall characteristics, size distribution of sand particles and boundary layer depth on vertical distribution of sand mass flux and particle mean horizontal velocity is analyzed, suggesting that all these three factors affect sand transport at different levels. In all cases, for small size groups, sand mass flux first increases with height and then decreases while for large size groups, it decreases exponentially with height and for middle size groups the behavior is in-between. The mean horizontal velocity for all size groups well fits experimental data, that is, increasing logarithmically with height in the middle height region. Wall characteristics greatly affects particle to wall collision and makes the flat bed similar to a Gobi surface and the rough bed similar to a sandy surface. Particle size distribution largely affects the sand mass flux and the highest heights they can reach especially for larger particles.