The effect of wind speed and bed slope on sand transport

This paper reports on a wind tunnel study of the effects of bed slope and wind speed on aeolian mass transport. The use of a sloping wind tunnel has enabled estimation of the friction angle α to be about 40° for saltating particles in the range 170–540 μm. A formula relating dimensionless mass transport to friction speed and bed slope is proposed, and mass transport data for five uniform sand samples and one non-uniform sand sample are shown to fit the equation well. In particular, the relationship reveals an overshoot in mass transport slightly above threshold collisions, a feature also evident when previous experimental data is re-examined. As the number of mid-air collisions between the saltating particles increases greatly with wind speed, the overshoot may occur as a result of increasing energy losses resulting from the collisions. Finally, it is demonstrated that data for saltating snow shows a similar overshoot in the dimensionless transport rate.     

Dimensionless input parameters in discrete element modeling and assessment of scaling techniques

A set of dimensionless input parameters were defined for DEM using a characteristic time which is a function of density and elastic modulus of particles and an arbitrary characteristic length. Dimensionless strain rate and mass damping ratio are inversely proportional to the characteristic time, and stress is normalized by elastic modulus to give dimensionless stress. It was demonstrated that the response of a model in the dimensionless scale is invariant with the choice of density, elastic modulus and the characteristic length if dimensionless strain rate and mass damping ratio are kept constant. Small time step is a prohibitive aspect of DEM. Scaling techniques are widely employed to enlarge the time step. Using the dimensionless scheme, it was learned that density scaling is equivalent to the use of a higher strain rate, and stiffness scaling results in a higher strain rate and an elevated stress state in the dimensionless scale.     

黏土颗粒间液桥从形成至断裂时毛细力的演化规律

通过土颗粒间形如液桥的水分毛细力认识和解释非饱和土的持水特性日趋为土力学界所关注。为研究黏土颗粒间液桥毛细力的演化规律,将黏土颗粒其简化为一对平行的片状颗粒,不考虑黏土颗粒表面的粗糙度。以Young-Laplace方程为基础,推得表征黏土颗粒间液桥毛细力关于固-液接触角和液桥体积的无量纲表达式;再令该无量纲表达式为0,绘制出固-液接触角与无量纲液桥体积的关系曲线,将1/3＜vLB＜1范围划分为液桥的断裂区、形成区和消失区,这种分区方法对描述黏土样在毛细阶段和吸附阶段的持水特性至关重要。最后,分析了液桥形成时毛细力与液桥体积和固-液接触角的关系,结果表明：当固-液接触角一定时,液桥的毛细力随其体积增大呈递增趋势,而当液桥体积一定时,其毛细力随固-液接触角增大呈先增后减的趋势。     

粗颗粒间液桥毛细力演化规律的动态计算方法

水力特性是非饱和土力学理论与工程的重要课题之一,土样水力特性的变化过程本质上可反映为土颗粒间毛细力的演化规律。为此,以粗颗粒为研究对象,将其简化为一对不等径球体颗粒,而其间的水分形态可视为形如圆环的液桥,不考虑颗粒重力和浮力的影响。以Young-Laplace方程为基础,先推得计算液桥毛细力的控制方程组,再结合液桥的无量纲体积最大、最小值及其外曲率半径的割线迭代算法提出求解毛细力的动态计算方法,进而研究毛细力与颗粒间距、颗粒半径比和液桥体积的无量纲关系,结果表明：当颗粒间距一定时,液桥的毛细力随其体积和颗粒半径比增大均呈递增趋势;当液桥体积一定时,其抗拉刚度随颗粒间距增大呈递减趋势。最后,利用已有文献中液桥毛细力与颗粒间距的实测关系,验证了该动态计算方法在表征液桥从形成至断裂时毛细力演化规律的有效性。     

泥沙起动判别标准探讨

分析了现有起动标准的不足之处。在制定起动标准时,不但要考虑水流条件,还必须考虑泥沙颗粒在床面的相对位置。对同一粒径的泥沙颗粒,其起动条件不是一个常量,而是位于一个区间,它随颗粒在床面的相对暴露度而变化。根据力学原理,推导了泥沙颗粒起动时的临界无量纲切应力公式,并对其系数取值进行了理论计算。计算结果表明,对同一粒径、在同一起动状态下,泥沙颗粒的起动临界条件并非一个常量。并从理论上分析得到了临界值的范围,对个别及少量起动,其无量纲临界切应力为0.021～0.042及0.041～0.062。     

双峰型非均匀沙推移运动特性及输移规律

采用图像识别与推移质动态监测技术,开展基于双峰型非均匀推移质的系列水槽试验.通过引入反映床面粗糙度、粘性底层特性与颗粒非均匀度η(粗细比)的综合水流强度函数Ψ_b、特征弗劳德数Fr_b,系统研究了不同水流强度与床沙组成条件下的推移质输移特性以及颗粒非均匀度对输沙率的影响.通过对关键因子的辨识与量纲分析,提出了双峰型非均匀推移质输移模式,建立了基于近壁特征因子的水流强度Ψ_b与非均匀推移质输移强度Φ'的函数关系.对双峰型底沙输移机理的分析表明,非均匀沙的组成特征使得η成为影响Φ'的重要参量;正是细粒对粗粒的解怙作用对粗沙运动产生重要影响,使推移质输移率与颗粒非均匀度间呈现驼峰关系,峰值对应的粗细比η_c约为3∶7.     

Fine particle deposition to initially starved,stationary, planar beds

Reported here are results from new flume experiments examining deposition and entrainment of inert, silt‐sized particles (with spherical diameters in the range from 20 to 60 μm) to and from planar, impermeable and initially starved beds underlying channel flows. Bed surfaces comprised smooth or fixed sand‐size granular roughness and provided hydraulically smooth to transitionally rough boundaries. Results of these experiments were analysed with a simple model that describes the evolution of vertically averaged concentration of suspended sediment and accommodates the simultaneous delivery to and entrainment of grains from the bed. The rate of particle arrival to a bed diminishes linearly, and the rate of particle entrainment increases by the 5/2 power, as the value of the dimensionless Saffman parameter S = u_*³/g’ν approaches a threshold value of order unity, where u is the conventional friction velocity of the turbulent channel flow, g’ is the acceleration due to gravity adjusted for the submerged buoyancy of individual particles and ν is the kinematic viscosity of the transporting fluid. This transport behaviour is consistent with the notion that non‐cohesive, silt‐sized particles can neither reach nor remain on an impermeable bed under flow conditions where mean lift imposed on stationary particles in the viscous sublayer equals or exceeds the submerged weight of individual particles. Within the size range of particles used in these experiments, particle size and the characteristic size of granular roughness, up to that of medium sand, did not affect rates of dimensionless arrival or entrainment to a significant degree. Instead, a new but consistent picture of fine‐particle transport is emerging. Silt‐sized material, at least, is subject to potentially significant interaction with the bed during intermittent suspension transport at intermediate flow speeds greater than the value required for initiation of transport (ca 20 cm sec^?1) but less than the value (ca 50 cm sec^?1) required by the Saffman criterion ensuring transport in fully passive suspension or, equivalently, ‘wash‐load’.     

EPS颗粒改良土作为寒区路基填料的抗冻性能研究

为研究EPS颗粒改良土作为寒区路基填料的抗冻性能,以冻胀率及割线模量为评测指标对改良土的抗冻性能进行了测试分析。通过冻胀率试验,获得了不同EPS颗粒掺量下冻胀率与冻融循环次数的关系。通过不同围压下的室内三轴试验,获得了改良前后土体割线模量随冻融循环次数增加的衰减规律,据此定义了无量纲单位残余割线模量比并建立了其与冻胀率的关系。结果表明：EPS颗粒的掺入对改良土的抗冻性能有明显的提升,有效地降低了改良土的冻胀率,随着EPS掺量增加,冻胀率逐渐减小,二者大致呈指数关系。冻胀率在冻融的初期发展较快,后期趋于平缓。改良后割线模量的衰减有大幅度的降低,割线模量随着冻融循环次数增加而逐渐减小,残余割线模量比与冻胀率之间近似呈二次抛物线关系。利用灰色理论可以对改良土割线模量的衰减规律进行预测,模型预测结果精度较高,可以作为路基工程冻害防治的参考。     

Dynamics of settling-driven convection beneath a sediment-laden buoyant overflow: Implications for the length-scale of deposition in lakes and the coastal ocean

A series of laboratory experiments was conducted in order to determine how settling-driven convection influences the length-scale over which the majority of particles settle beneath a buoyant sediment-laden plume spreading over a denser saline layer. This system is analogous to sediment-laden river water spreading into a lake or the coastal ocean. The key dimensionless parameter that controls the settling dynamics of such flows is the density ratio, defined as the ratio of density differences due to the added salt and sediment. For a buoyant plume, this ratio has to be greater than unity, so that the experiments in the current study were performed for density ratios between one and five. When density ratio was close to one, settling-driven convection was vigorous and the length-scale over which sedimentation occurred was very small. A strong secondary turbidity current was generated in this case. On the other hand, for larger values of density ratio, the predicted length-scale over which a secondary plume was generated increased in proportion to the density ratio. A complete mathematical expression for this length-scale was developed using recent theory that described the timescale over which settling-driven convection evolved. The theoretically predicted propagation length-scale showed very good quantitative agreement with laboratory experiments. The use of the dimensionless density ratio allows the expression to predict which sediment-laden river plumes in lakes and the coastal ocean could quickly form secondary turbidity currents.     

外壁恒温条件下冻结管壁热流密度变化规律数值计算研究

应用相似理论将冻结温度场方程和参数无量纲化,然后建立无量纲的数值计算模型,研究了在外表面恒温情况下冻结管壁的无量纲热流密度与各无量纲影响参数间的关系.结果表明:无量纲热流密度随未冻土与冻土的导热系数之比线性增加,与冻土与未冻土的容积比热之比成线性关系,随无量纲结冰潜热线性增长,与无量纲盐水温度成线性关系;可表示为无量纲时间的负三分之一次方的线性函数.最后对167706次数值计算结果进行了回归分析,得到具有较高精度的回归计算公式.     

Attrition efficiency in the decontamination of stormwater sediments

The purpose of this research is to propose a laboratory method based on attrition and sieving for the treatment of runoff water sediments in the aim of developing a pilot unit. The attrition process serves to remove fine particles and pollutants from the surface of coarse stormwater sediment particles. In all cases, the efficiency of pollutant removal is dependent upon various parameters, including cutoff threshold, residence time, solid density, temperature and impeller speed. This article presents the optimization of these various parameters along with method efficiency; for this work, several sediments have been tested. The results indicate that an attrition scrubber may be effectively used to remediate contaminated sediment and that reuse is definitely possible. A model of the method will also be proposed to study the behavior of fine particles and pollutants.     

A simple model for the prediction of the deflation threshold shear velocity of dry loose particles

A very important parameter in aeolian equations is the deflation threshold shear velocity, which quantifies the instant of particle motion. In this paper, a simple model is presented for the prediction of the threshold shear velocity of dry loose particles. It has the same functional form as the widely used models of Bagnold (1941) and Greeley & Iversen (1985), but differs in its treatment of the so‐called threshold parameter. As the new expression was based on the moment balance equation used by Greeley & Iversen, it includes a function for the aerodynamic forces, including the drag force, the lift force and the aerodynamic moment force, and a function for the interparticle forces. The effect of gravitation is incorporated in both functions. However, rather than using an implicit function for the effect of the aerodynamic forces as in the Greeley & Iversen model, a constant aerodynamic coefficient was introduced. From consideration of the van der Waals' force between two particles, it was also shown that the function for the interparticle cohesion force is inversely proportional to the particle diameter squared. The model was calibrated on data reported by Iversen & White (1982). The new expression compared, at least for terrestrial conditions, very well with the Greeley & Iversen model, although it is much simpler. It was finally validated with data from wind‐tunnel experiments on different fractions of dune sand and sandy loam soil aggregates. The soil aggregates were treated as individual particles with a density equal to their bulk density. The good agreement between observations and predictions means that, when predicting mass transport of particles above a given soil, minimally dispersed particle‐size distributions should be considered rather than the granulometric composition of the soil.     

Intermittent saltation

During a typical wind erosion event, large variations in wind strength produce temporal variations in saltation activity. The focus of this paper is on a special type of unsteady behaviour - intermittent saltation - a process characterized by bursts of blowing soil interspersed with periods of inactivity. We report here measurements from a field study designed to measure intermittent saltation during three separate 1-h periods. Our measurements show that natural wind erosion events consist of intermittent bursts of blowing soil often occupying a small fraction of the total time. We have managed to describe the level of intermittency by a simple and universal mathematical expression. We find that the level of intermittency is governed by whether typical wind fluctuations span the gap between the mean wind speed and threshold wind speed. We propose a nondimensional number which expresses the ratio of these velocity scales, called the relative wind strength, and find that the level of intermittency can be described by a simple distribution function of the relative wind strength.     

潜孔锤反循环钻进中心通道内岩屑运移数学模型

针对潜孔锤反循环钻进中上返岩屑回落堵塞中心通道问题,将气力输送理论应用于反循环携屑。采用气力输送理论对潜孔锤反循环钻进中心通道内岩屑颗粒进行受力分析;从动力学角度建立潜孔锤反循环钻井中心通道内岩屑运移数学模型;通过编程计算求解岩屑在不同初始条件下速度随时间、空间变化特性。结果表明:初速度大于沉降末速度的岩屑颗粒作减速运动,初速度小于沉降末速度的岩屑颗粒作加速运动,最终岩屑颗粒沉降末速度趋于同一个值,与初速度无关;气固速度比与岩屑尺寸有关,随着岩屑颗粒直径增大,气固比降低。该数学模型为改善反循环钻进携屑能力提供了理论依据。      

考虑起始坡降双层地基的大应变非线性固结

天然软土成层分布特性及土中渗流存在起始水力坡降的现象已被人们熟知。但变荷载下能同时考虑黏土中起始水力坡降、软土非线性压缩渗透特性及大应变特性的双层地基固结理论还鲜见报道。在拉格朗日坐标系中建立以超静孔压为变量的双层软土地基大应变非线性固结模型并给出其有限差分解。通过与考虑起始水力坡降的单层地基大应变非线性固结数值计算结果对比,验证了差分解的可靠性。着重分析了上、下土层起始坡降无量纲参数R1、R2对双层地基固结性状的影响,分析在大应变与小应变假定下双层地基超静孔压消散及固结沉降变形的异同。结果表明：上层土无量纲参数R1对双层地基固结性状的影响程度较下层土无量纲参数R2显著;大应变假定下双层地基渗流前锋的下移速度要快于小应变假定下的移动速度;大应变假定下考虑起始水力坡降的双层软土地基超静孔压消散速率要比小应变假定下快,且大应变假定下考虑起始水力坡降的双层地基最终沉降量要比小应变假定下大。     

Settling velocity and drag coefficient of platy shell fragments

Settling velocity of bioclastic particles in coastal and shallow marine environments is essential for interpreting depositional facies and processes. There is, however, a paucity of accurate formulae for predicting the settling velocities and drag coefficients of platy biogenic particles in particular. This study provides experimental settling data based on 320 platy shell fragments from a sediment core recovered in Li'an Lagoon, south-eastern Hainan Island, China. The results indicate that the settling velocities of platy shell fragments are strongly correlated with nominal diameters and Corey shape factors (ranging from 0.02 to 0.20 in this study). On this basis, a practical equation of acceptable accuracy was established for platy particles, relating dimensionless settling velocities to dimensionless diameters and Corey shape factors. Similarly, another empirical formula for quickly calculating the equivalent diameter of platy shell fragments in practice was proposed as well. Regarding the strong dependence of the drag coefficients using equivalent spherical diameters to Corey shape factors, the drag coefficient based on the diameter of the equivalent maximum projected area remains almost constant and is hence physically well-suited for the definition of grain drag coefficients. The settling data of this study has extended the lower Corey shape factors limit of bioclastic particles, and the equations presented here can be used for quantitative interpretations of sedimentary records, modelling of depositional processes and investigations of other platy particles.     

Prediction of velocity-dip-position over entire cross section of open channel flows using entropy theory

An analytical model to predict the velocity-dip-position is presented in this study. Unlike the previous studies where empirical or semiempirical models were suggested, in this study the model is derived by using entropy theory. Using the principle of maximum entropy, the model for dip-position is derived by maximizing the Shannon entropy function after assuming dimensionless dip-position as a random variable. No estimation of empirical parameter is required for calculating dip-position from the proposed model. The model is capable of predicting the velocity-dip-position at any section from sidewall region (along lateral direction) of an open channel with any aspect ratio. The ratio of mean to maximum value of dip-position is analyzed from data, and it is found that ratio is almost constant for narrow open channels and it increases with aspect ratio for wide open channels. A relation is also proposed to predict this ratio in case of wide open channels. The developed model of velocity-dip-position is tested with existing experimental data for a wide range of aspect ratio and is also compared with other empirical models. The present model shows good agreement with the observed data and is comparable with the existing models.     

Threshold for particle entrainment into suspension

ABSTRACT Laboratory observations regarding the limit conditions for particle entrainment into suspension are presented. A high‐speed video system was used to investigate conditions for the entrainment of sediment particles and glass beads lying over a smooth boundary as well as over a rough bed. The results extend experimental conditions of previous studies towards finer particle sizes. A criterion for the limit of entrainment into suspension is proposed which is a function of the ratio between the flow shear velocity and particle settling velocity. Observations indicate that particles totally immersed within the viscous sublayer can be entrained into suspension by the flow, which contradicts the conclusions of previous researchers. A theoretical analysis of the entrainment process within the viscous sublayer, based on force–balance considerations, is used to show that this phenomenon is related to turbulent flow events of high instantaneous values of the Reynolds stress, in agreement with previous observations. In the case of experiments with a rough bed, a hiding effect was observed, which tends to preclude the entrainment of particles finer than the roughness elements. This implies that, as the ratio between particle and roughness element sizes becomes smaller, progressively higher bed shear stresses are required to entrain particles into suspension. On the other hand, an overexposure effect was also observed, which indicates that a particle moving on a smooth bed is more prone to be entrained than the same particle moving on a bed formed by identical particles.     

Saltation threshold on Earth, Mars and Venus

New formulations valid for wide ranges of particle diameter and density and gas density are presented for prediction of saltation threshold speed for small particles. A low-air-density wind tunnel was used to extend the range of previous investigations and to separate the effects of Reynolds number and interparticle forces of cohesion. The new formulations are used to predict saltation threshold for atmospheric conditions on the surface of the Earth, Mars, and Venus.     

Settling velocities of particulate systems 15: Velocities in turbulent Newtonian flows

The success in previous papers of this series in obtaining closed form solutions to the sedimentation of individual particles and to the suspensions of particles in fluids motivated the application of the same procedure to the turbulent flow of Newtonian fluids in conduits.The purpose of this paper is to present explicit equations for the Fanning friction factor and of the average velocity of a Newtonian fluid in a cylindrical tube, for any value of the Reynolds number and dimensionless roughness.Previously published equations for the friction coefficient are either implicit or explicit in the friction coefficient but dependent on the Reynolds number. Since the Reynolds number depends on the pipe diameter and fluid velocity, in both cases iterative method must be used to calculate the flowrate.By defining a dimensionless pipe diameter and a dimensionless average fluid velocity, and using conduit-flow material properties to define characteristic size and velocity parameters, explicit equations were obtained for the friction coefficient and for the average flow velocity valid for different roughness and flow regimes.