力学与沙尘暴

沙尘暴是近年来广受关注的环境问题之一,不仅是大气、地学等相关学科研究的重要内容之一,也因为其中有许多基础科学的共性问题,特别是力学问题,而得到其他学科的关注.本书的目的是通过对沙尘暴及其相关问题的介绍,揭示其与力学之间的联系,指出这些力学问题的解决对于力学学科本身发展所将起到的带动作用.全书共包括25个章节,分别涉及风沙流形成与发展的微观机理、沙尘暴的测量与预报、风沙电现象及其对电磁波的影响、风沙地貌、防沙治沙工程等多个方面.     

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

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

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

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

AN ARITHMETIC METHOD OF SEGMENTING MOVING AEOLIAN SAND PARTICLES＇ IMAGES FROM BACKGROUND INFORMATION OF DIGITAL HIGH-SPEED PHOTOGRAPHY IMAGES

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

裸露地表土壤风蚀及其电场数值模拟

土壤风蚀及电场机理研究是认识其灾害本质,进而实现有效防治的基础。采用颗粒流体动力学方法,建立裸露地表土壤风蚀的气-固耦合模型,并基于水的电离理论,建立了运动土壤颗粒碰撞带电模型。数值计算结果不仅从理论上再现了裸露地表土壤风蚀的发生、发展及达到动态稳定的过程,而且定量分析了土壤风蚀过程颗粒的带电特性及其产生的风蚀电场。在风-颗粒相互耦合作用下,裸露地表土壤颗粒的风蚀过程最终达到动态稳定,土壤风蚀量随高度呈指数规律递减,且贴近地表层风蚀量较大。土壤颗粒在运动中的碰撞使得颗粒发生带电现象,粒径较大颗粒带正电,粒径较小颗粒带负电;由带电土壤颗粒产生的风蚀电场强度随高度减小,且在风蚀过程中逐渐增强,在风蚀达到动态平衡时风蚀电场趋于稳定。本文综合考虑了风-土壤-电场的土壤风蚀过程,这将为建立更加符合实际、较为全面的土壤风蚀模型提供理论基础。     

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

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

沙粒跃移运动的数值模拟

沙粒在风场作用下形成的风沙运动,直接导致一系列环境问题.现有的数值模拟并没有揭示风沙流的二维性质.本文提出一种合理有效的模拟方法来揭示风沙流时空的变化.气相运动采用大涡模型,而固相(沙粒)处理为高散体系(DEM),并利用已有的沙粒与床面碰撞的激溅函数,以获得整个风沙系统运动特性.结果表明:其发展过程与以往的试验结果相比较符合:风沙流达到稳定后,不同截面输沙率沿高度分布与平均输沙率相比差别很大,并大约在9cm高处,平均输沙率出现极值,单宽输沙率沿风向分布同样具有非均匀性.     

沙丘背风侧不同铁路结构形式对风沙环境的适应性分析

当铁路穿越大风沙漠地区时,风沙灾害会对铁路工程及其正常运营产生严重威胁,而设计一种合理的铁路结构形式能够减小风沙沉积对铁路工程的危害.在本文中,以敦煌至格尔木铁路沙山沟段落为研究对象,采用多相流的方法对越过沙丘的风沙运动过程进行数值模拟,分别讨论了风沙运动对位于沙丘背风坡的铁路路基工程和桥梁工程的影响.主要的模拟结果显示:路基工程明显降低了风速并且将沙丘后的回流区分成了两部分,而桥梁工程的导流效应则压缩了沙丘背风坡的回流区;轨道间的道碴增大了铁路表面的粗糙度,在轨道间有少量沙粒沉积,而路基工程两侧则有大量积沙;铁路表面的积沙量与摩阻风速呈现出非线性关系,随着摩阻风速的增大,路基工程沙粒沉积的增加速度大于风蚀能力的增加速度,而桥梁工程则正好相反.在防止风沙危害铁路方面,设置桥梁工程明显优于路基工程.本研究为风沙运动对铁路工程的影响提供了理论支持,也为今后的铁路工程设计提供了新的思路与研究工具.     

沙丘背风侧不同铁路结构形式对风沙环境的适应性分析

当铁路穿越大风沙漠地区时,风沙灾害会对铁路工程及其正常运营产生严重威胁,而设计一种合理的铁路结构形式能够减小风沙沉积对铁路工程的危害. 在本文中,以敦煌至格尔木铁路沙山沟段落为研究对象,采用多相流的方法对越过沙丘的风沙运动过程进行数值模拟,分别讨论了风沙运动对位于沙丘背风坡的铁路路基工程和桥梁工程的影响. 主要的模拟结果显示:路基工程明显降低了风速并且将沙丘后的回流区分成了两部分,而桥梁工程的导流效应则压缩了沙丘背风坡的回流区;轨道间的道碴增大了铁路表面的粗糙度,在轨道间有少量沙粒沉积,而路基工程两侧则有大量积沙;铁路表面的积沙量与摩阻风速呈现出非线性关系,随着摩阻风速的增大,路基工程沙粒沉积的增加速度大于风蚀能力的增加速度,而桥梁工程则正好相反. 在防止风沙危害铁路方面,设置桥梁工程明显优于路基工程. 本研究为风沙运动对铁路工程的影响提供了理论支持,也为今后的铁路工程设计提供了新的思路与研究工具.      

NUMERICAL SIMULATION OF WIND-BLOWN SAND MOVEMENT BASED ON SPH

运用光滑粒子流体动力学(smoothed particles hydrodynamics, SPH)方法对沙粒和气流的相互耦合运动特性进行了分析,研究提出了风沙流的SPH数值方法并进行了数值模拟. 首先提出了风沙流的SPH建模方法和基本理论,建立了风沙流动的SPH数值模拟平台. 其次通过建立风沙流的SPH 模型并施加边界条件,对自然风作用下沙粒的运动情况进行了数值模拟,详细分析了沙粒运动轨迹及特性,最后通过与相关研究成果对比分析,验证了完善后的SPH方法有效性. 通过考虑气流场的可变性,在风沙流SPH计算模型中引入了加载(起风)和卸载(停风)方式,观察并对比分析了沙粒的运动轨迹和特性. 为进一步研究风沙流的实时动态非线性行为提供了SPH理论基础和数值分析方法.     

风沙运动研究中的若干关键力学问题

针对沙漠化（严重时导致沙尘暴发生）这一国家关注的重大环境的问题,在简要介绍风沙运动力学机理研究的主要发展阶段和部分最新研究进展的基础上,提出了风沙运动研究中值得关注的若干关键力学问题。     

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.     

SOURCE PROFILES AND FINGERPRINTS OF FINE AND COARSE SANDS RESUSPENDED FROM SOILS SAMPLED IN CENTRAL INNER MONGOLIA

Fingerprints and source profiles of fine and coarse sands that originate from Central Inner Mongolia during Asian continental sandstorms （ACS） can be used to identify the odgin of Asian sands and to trace them as they travel downwind. Soil samples collected at various land surfaces in Central Inner Mongolia were resuspended using a dry powder atomizer in an enclosure chamber. The resuspended sands were then sampled by two dichotomous samplers situated at the bottom of the enclosure chamber for fine （PM2 5） and coarse （PM2.5-10） sands, respectively. The chemical composition of sands, including water-soluble ionic species, metallic contents, and carbonaceous contents, were further analyzed. Results from resuspension tests indicated that the soils contained considerably more coarse particles than fine Moreover, Mg, K, Al, and Fe in coarse sand had strong correlations with each other. The ratio of Mg, K. Fe （or Al） to Al （or Fe） and OC/EC in the coarse sands can be used as the fingerprints of Asian sands originating from Central Inner Mongolia.     

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.     

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.     

Probability of rebound and eject of sand particles in wind-blown sand movement

When incident particles impact into a sand bed in wind-blown sand movement, rebound of the incident particles and eject of the sand particles by the incident particles affect directly the development of wind sand flux. In order to obtain rebound and eject lift-off probability of the sand particles, we apply the particle-bed stochastic collision model presented in our pervious works to derive analytic solutions of velocities of the incident and impacted particles in the post-collision bed. In order to describe randomness inherent in the real particle-bed collision, we take the incident angle, the impact position and the direction of resultant action of sand particles in sand bed on the impacted sand particle as random variables, and calculate the rebound and eject velocities, angles and coefficients (ratio of rebound and eject velocity to incident velocity). Numerical results are found in accordance with current experimental results. The rebound and eject lift-off probabilities versus the incident and creeping velocities are predicted. The project was supported by the National Natural Science Foundation of China (10532040, 10601022). The English text was polished by Yunming Chen.     

Radiation of electromagneticwaves by a spheroidal charged drop oscillating in a uniform electrostatic field

The dispersion relation is derived for a spheroidal charged drop oscillating in a uniform electrostatic field and radiating electromagnetic waves in the first and second orders with respect to the dimensionless oscillation amplitude and the steady-state deformation, respectively. For an individual drop and a model cloud the intensity of the electromagnetic radiation and the frequency bandwidth are estimated as functions of the drop dimensions and charge and the external electric field strength.     

Research in mechanical model of bionic foot intruding into sands with different physical characteristics

In this study, a bionic foot with sand fixation and fluidization limitation functions was designed. Also a rectangular foot with the same sizes, named the common foot, was designed for comparison. Three kinds of quartz sands were selected to study how particle size, shape and compactness affected the intrusion performances of mechanical feet. The intrusion resistive forces and pressures of the bionic foot on these three kinds of quartz sands were all smaller compared with the common foot. Discrete element simulations showed particle disturbance areas were smaller and particle motion trends were more consistent under the bionic foot versus the common foot. The intrusion resistive forces of these two kinds of mechanical feet firstly increased and then decreased with the increasing particle sizes of quartz sands. Moreover, the intrusion resistive force on spherical particles was less than that of irregular particles for both the bionic foot and the common foot. The corresponding resistive forces of mechanical feet were characterized based on quartz sand compactness. The classic pressure-sinkage model was modified based on the intrusion tests, and the relationships between intrusion resistive force and mechanical foot depth were obtained.