Jump at the onset of saltation

Toward a large-scale particle-based parallel simulator of Aeolian sand transport, including a model for mobile sand availability

We develop a numerical tool for particle-based simulations of Aeolian sand transport. Our model combines a Discrete-Element-Method for the sand particles with an efficient hydrodynamic description of the average turbulent horizontal wind velocity field over the granular bed, which has been developed in previous work and accounts for the two-way coupling of the granular and fluid phases. However, here we implement our model within the open source library LAMMPS for granular massively parallel simulations and incorporate a new grid coarsening scheme for the wind model. We show that our model quantitatively reproduces observed values of the steady-state (saturated) sand flux under various flow conditions. Furthermore, we model different conditions of mobile sand availability and find a strong dependence of the sand flux on this availability.

A Numerical Study of Aeolian Sand Particle Flow Incorporating Granular Pseudofluid Optimization and Large Eddy Simulation

A numerical investigation of aeolian sand particle flow in atmospheric boundary layer is performed with a Eulerian–Eulerian granular pseudofluid model. In this model, the air turbulence is modelled with a large eddy simulation, and a kinetic–frictional constitutive model incorporating frictional stress and the kinetic theory of granular flow is applied to describe the interparticle movement. The simulated profiles of streamwise sand velocity and sand mass flux agree well with the reported experiments. The quantitative discrepancy between them occurs near the sand bed surface, which is due to the difference in sand sample, but also highlights the potential of the present model in addressing near-surface mass transport. The simulated profiles of turbulent root mean square (RMS) particle velocity suggest that the interparticle collision mainly account for the fluctuation of sand particle movement.

Numerical simulation of sand transfer in wind storm using the Eulerian-Lagrangian two-phase flow model

In this paper a two-dimensional gas-solid flow model is used to investigate the sand particles carrying velocity of the Iran eastern desert area around the railway track as a case study. Reynolds-averaged Navier-Stokes (RANS) equations and Discrete Phase Method (DPM) are used to simulate the characteristic movement of sand particles in wind flow. A random sample is gathered from the sand near the railway in Iran deserts. The sample is classified based on weight and diameter according to AASHTOO T27 and sand distribution is determined. Using simulations, the carrying velocity of sand in each category in wind storm is determined. Finally, the sand distribution of the sample is imported to the model by the Rosin-Rummler dissipation model. The behavior of sand particles in storm considering wind blowing scheme of desert is studied parametrically. The results can be used for estimating the sand mitigation of a special desert and land desertification control around railway tracks.

Onset and cessation of motion in hydrodynamically sheared granular beds

We performed molecular dynamics simulations of granular beds driven by a model hydrodynamic shear flow to elucidate general grain-scale mechanisms that determine the onset and cessation of sediment transport. By varying the Shields number (the nondimensional shear stress at the top of the bed) and particle Reynolds number (the ratio of particle inertia to viscous damping), we explore how variations of the fluid flow rate, particle inertia, and fluid viscosity affect the onset and cessation of bed motion. For low to moderate particle Reynolds numbers, a critical boundary separates mobile and static states. Transition times between these states diverge as this boundary is approached both from above and below. At high particle Reynolds number, inertial effects become dominant, and particle motion can be sustained well below flow rates at which mobilization of a static bed occurs. We also find that the onset of bed motion (for both low and high particle Reynolds numbers) is described by Weibullian weakest-link statistics and thus is crucially dependent on the packing structure of the granular bed, even deep beneath the surface.

Bursts in discontinuous Aeolian saltation

Close to the onset of Aeolian particle transport through saltation we find in wind tunnel experiments a regime of discontinuous flux characterized by bursts of activity. Scaling laws are observed in the time delay between each burst and in the measurements of the wind fluctuations at the fluid threshold Shields number θ_c. The time delay between each burst decreases on average with the increase of the Shields number until sand flux becomes continuous. A numerical model for saltation including the wind-entrainment from the turbulent fluctuations can reproduce these observations and gives insight about their origin. We present here also for the first time measurements showing that with feeding it becomes possible to sustain discontinuous flux even below the fluid threshold.

Numerical simulation of subaqueous chute flows of granular materials

In this paper we report on numerical studies of unsteady, gravity-driven flow of a subaqueous erodible granular bed on an inclined plane. According to our simulations, the evolution of the flow can be partitioned in three phases. In the first phase, due to the onset of an interfacial instability, the material interface deforms into a series of long waves. In the second phase, these waves are transformed to skewed vortex ripples that grow in time and eventually coalesce. The computed wavelengths of these ripples are in good agreement with previously reported experimental measurements. In the third phase of the flow evolution, the high fluid velocities wash out the vortex ripples and a layer of rapidly moving particles is formed at the material interface. The predicted granular velocities comprise two segments: a concave one at the vicinity of the material interface, where the maximum is attained, followed by a slightly convex one, where they decrease monotonically to zero. The same trend has been reported in experimental results for the corresponding steady flows. Finally, we investigate via a parametric study the effect of the configuration stresses, which represent contact forces between grains. As it turns out, such stresses have a stabilizing effect, in the sense that increasing their magnitude inhibits the formation of vortex ripples.

Analysis of wind-blown sand movement over transverse dunes

Wind-blown sand movement often occurs in a very complicated desert environment where sand dunes and ripples are the basic forms. However, most current studies on the theoretic and numerical models of wind-blown sand movement only consider ideal conditions such as steady wind velocity, flat sand surface, etc. In fact, the windward slope gradient plays a great role in the lift-off and sand particle saltation. In this paper, we propose a numerical model for the coupling effect between wind flow and saltating sand particles to simulate wind-blown sand movement over the slope surface and use the SIMPLE algorithm to calculate wind flow and simulate sands transport by tracking sand particle trajectories. We furthermore compare the result of numerical simulation with wind tunnel experiments. These results prove that sand particles have obvious effect on wind flow, especially that over the leeward slope. This study is a preliminary study on windblown sand movement in a complex terrain, and is of significance in the control of dust storms and land desertification.

Intermittency of aeolian saltation

Saltation motion of sand grains in a steady wind was measured using a high-speed camera at very high frequency in a wind tunnel. A Heaviside-type function was defined to quantificationally describe an inherent property of saltation, i.e. intermittency. Kurtosis and periodicity of state function are statistical manifestations of intermittency. In addition, the strong autocorrelation of time series of volume concentration clearly confirms that saltation is not a completely random process at the timescale of subsecond. Formation mechanism, especially turbulent structures responsible for intermittent saltation, remains to be revealed from the viewpoint of classical mechanics.

Direct numerical simulations of aeolian sand ripples

Aeolian sand beds exhibit regular patterns of ripples resulting from the interaction between topography and sediment transport. Their characteristics have been so far related to reptation transport caused by the impacts on the ground of grains entrained by the wind into saltation. By means of direct numerical simulations of grains interacting with a wind flow, we show that the instability turns out to be driven by resonant grain trajectories, whose length is close to a ripple wavelength and whose splash leads to a mass displacement toward the ripple crests. The pattern selection results from a compromise between this destabilizing mechanism and a diffusive downslope transport which stabilizes small wavelengths. The initial wavelength is set by the ratio of the sediment flux and the erosion/deposition rate, a ratio which increases linearly with the wind velocity. We show that this scaling law, in agreement with experiments, originates from an interfacial layer separating the saltation zone from the static sand bed, where momentum transfers are dominated by midair collisions. Finally, we provide quantitative support for the use of the propagation of these ripples as a proxy for remote measurements of sediment transport.

Numerical modeling of wind-blown sand on Mars

Recent observation results show that sand ripples and dunes are movable like those on Earth under current Martian climate. And the aeolian process on Mars therefore is re-attracting the eyes of scientific researchers in different fields. In this paper, the spatial and temporal evolution of wind-blown sand on Mars is simulated by the large-eddy simulation method. The simulations are conducted under the conditions of both friction wind speed higher and lower than the "fluid threshold", respectively. The fluid entrainment of the sand particles, the processes among saltation sand particles and sand bed, and the negative feedback of sand movement to flow field are considered. Our results show that the "overshoot" phenomenon also exists in the evolution of wind-blown sand on Mars both temporally and spatially; impact entrainment affects the sand transport rate on Mars when the wind speed is smaller or larger than the fluid threshold; and both the average saltation length and height are one order of magnitudes larger than those on Earth. Eventually, the formulas describing the sand transport rate, average saltation length and height on Mars are given, respectively.

Analytical model for flux saturation in sediment transport

The transport of sediment by a fluid along the surface is responsible for dune formation, dust entrainment, and a rich diversity of patterns on the bottom of oceans, rivers, and planetary surfaces. Most previous models of sediment transport have focused on the equilibrium (or saturated) particle flux. However, the morphodynamics of sediment landscapes emerging due to surface transport of sediment is controlled by situations out of equilibrium. In particular, it is controlled by the saturation length characterizing the distance it takes for the particle flux to reach a new equilibrium after a change in flow conditions. The saturation of mass density of particles entrained into transport and the relaxation of particle and fluid velocities constitute the main relevant relaxation mechanisms leading to saturation of the sediment flux. Here we present a theoretical model for sediment transport which, for the first time, accounts for both these relaxation mechanisms and for the different types of sediment entrainment prevailing under different environmental conditions. Our analytical treatment allows us to derive a closed expression for the saturation length of sediment flux, which is general and thus can be applied under different physical conditions.

Numerical simulation of wind sand movement in straw checkerboard barriers

Straw checkerboard barrier (SCB) is the most representative antidesertification measure and plays a significant role in antidesertification projects. Large-eddy simulation and discrete-particle tracing were used to numerically simulate the wind sand movement inside the straw checkerboard barrier (SCB), study the movement characteristics of sand particles, find the transverse velocities of sand particles and flow field, and obtain the contour of the transverse velocity of coupled wind field within the SCB. The results showed that 1) compared with that at the inlet of the SCB, the sand transport rate inside the SCB greatly decreases and the speed of sand grain movement also evidently drops, indicating that the SCB has very good sand movement preventing and fixing function; 2) within the SCB there exists a series of unevenly distributed eddies of wind sand flow, their strength decreases gradually with increasing the transverse distance; 3) affected by eddies or reflux, sand particles carried by the wind sand flow have to drop forward and backward the two interior walls inside the SCB, respectively, forming a v-shaped sand trough; 4) the sand transport rate gradually decreases with increasing number of SCBs, which reveals that the capacity of the wind field to transport sand particles decreases. This research is of significance in sandstorm and land desertification control.

Midair collisions enhance saltation

Here we address the old question in aeolian particle transport about the role of midair collisions. We find that, surprisingly, these collisions do enhance the overall flux substantially. The effect depends strongly on restitution coefficient and wind speed. We can explain this observation as a consequence of a soft bed of grains which floats above the ground and reflects the highest flying particles. We make the unexpected observation that the flux is maximized at an intermediate restitution coefficient of about 0.7, which is comparable to values experimentally measured for collisions between sand grains.

Incident angle of saltating particles in wind-blown sand

Incident angle of saltating particles plays a very important role in aeolian events. In this paper, the incident angles of sand particles near the sand bed were measured in wind tunnel. It reveals that the incident angles range widely from 0° to 180° and thereby the means of angles are larger than published data. Surprisingly, it is found the proportion that angles of 5°–15° occupy is far below previous reports. The measuring height is probably the most important reason for the measurement differences between this study and previous investigations.

Influence of sand grain diameter and wind velocity on lift-off velocities of sand particles

In this paper, the velocities of sand particles near the sand bed in the saltation cloud were measured in a wind tunnel through an improved experimental scheme of the Particle Image Velocimetry (PIV) system. The influences of the diameter of sand particles in the saltation cloud and wind velocity on the probability distribution function (PDF) of lift-off velocities of sand particles were investigated. Results demonstrate that for the sand particles saltating above the sand bed with the mean grain diameter (d m = 0.3 mm), smaller and larger ones have the same velocity distribution, and wind velocity has no obvious influence on the distribution shape of the lift-off velocities, i.e., the PDFs of the horizontal and vertical lift-off velocities both follow a lognormal distribution, but the diameter of sand particles in the saltation cloud and wind velocity have an influence on the parameters of the PDF of horizontal and vertical lift-off velocities. Eventually, we present formulas to describe the PDF of lift-off velocities of sand particles with regard to the influence of wind velocity and the diameter of sand particles in the saltation cloud above the sand bed with d m = 0.3 mm.

The physics of wind-blown sand and dust

The transport of sand and dust by wind is a potent erosional force, creates sand dunes and ripples, and loads the atmosphere with suspended dust aerosols. This paper presents an extensive review of the physics of wind-blown sand and dust on Earth and Mars. Specifically, we review the physics of aeolian saltation, the formation and development of sand dunes and ripples, the physics of dust aerosol emission, the weather phenomena that trigger dust storms, and the lifting of dust by dust devils and other small-scale vortices. We also discuss the physics of wind-blown sand and dune formation on Venus and Titan.