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

United States

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

Powders &amp; Grains 2021

[EPJ Web of Conferences, Volume 249 (2021)](https://www.epj-conferences.org/articles/epjconf/abs/2021/03/contents/contents.html)

Follow us on [Twitter](https://twitter.com/g2021_p)

**Meet the Organizers**

<center>
<figure class="video_container">
	  <iframe src="https://screencast-o-matic.com/watch/cr130qV10B3?v=6&ff=1&title=0&controls=1" width=826 height=526  frameborder="0" allowfullscreen="true"> </iframe>
</figure>
	</center>


Welcome to the 2021 Powders and Grains conference. Please review the event's **Code of Conduct** [here](https://drive.google.com/file/d/1P_bOrjI2jSE55twRPVzWl5VkfARR-qBL/view?usp=sharing). 


**Powders & Grains 2021 Walkthrough Screencast**

This video will demonstrate how to navigate the conference and virtual platform -- and network with other attendees and speakers.


<center>
<figure class="video_container">
	  <iframe src="https://screencast-o-matic.com/watch/cr1vnKV1pgn?v=6&ff=1&title=0&controls=1" width=826 height=526  frameborder="0" allowfullscreen="true"> </iframe>
</figure>
	</center>

To register for the event, please follow the link below. 

Powders & Grains 2021

Powders & Grains 2021 will be held in Buenos Aires (Argentina, VIRTUAL) July 5, 13, 21, 29 and Aug 6, 2021. The aim of the conference is to give an up-to-date picture of the broad research field of granular media, powders and grains. Contributions from experts around the world will focus on research topics that cover a wide range of open problems in the field.

technical paper

For homogeneous systems like classical fluid dynamics and structural mechanics, finite element method (FEM) grid generation has reached a mature state. On the other hand, for multi-physics-problems like fluids with a high density of immersed particles, many researchers may not even be aware of the types of instabilities which may be triggered by unsuitable meshes. We review common types of grid generation, point out previously unrecognised types of instabilities for particles in fluids as well as remedies to obtain particle-fluid simulations with higher stability and fewer redundant degrees of freedom.

Grid-algorithm improvements for dense suspensions of discrete element particles in finite element fluid simulations

Discrete element modelling (DEM) is widely used to simulate granular systems, nowadays routinely on graphical processing units. Graphics processing units (GPUs) are inherently designed for parallel computation, and recent advances in the architecture, compiler design and language development are allowing general-purpose computation to be computed on multiple GPUs. Application of DEM to bonded particle systems are much less common, with a number of open research questions remaining. This study outlines a Bonded-Particle Research DEM Framework, PyBONDEM-GPU, written in Python. This framework leverages the parallel nature of GPUs for computational speed-up and the rapid prototype flexibility of Python. Python is faster and easier to learn than classical compiled languages, making computational simulation development accessible to undergraduate and graduate engineers. PyBONDEM-GPU leverages the Numba-CUDA module to compile Python syntax for execution on GPUs. The framework enables research of fibre pull-out from fibre-matrix embeddings. Bonds are simulated between all interacting particles. The performance of PyBONDEM-GPU is compared against Python CPU implementations of PyBONDEM using the Numpy and Numba-CPU Python modules. PyBONDEM-GPU was found to be 1000 times faster than the Numpy implementation and 4 times faster than the Numba-CPU implementation to resolve forces and to integrate the equations of motion.

PyBONDEM-GPU: a discrete element bonded particle Python research framework – Development and examples

The present paper describes laboratory tests on steady granular flows down a smooth chute using X-ray radiography. By using an original granular chute set-up with a gate at its end, standing discontinuites in height, velocity and density, namely jumps, were produced during the granular flows. The X-ray radiography was successful to get consistent measurements of the free-surface and density profiles of the flows along the chute under different flow states: the (nearly) uniform flows before the jumps, the gradually-varied flows after the jumps, and the highly non-uniform flows across the jumps. This demonstrates the efficiency of the X-ray radiography to measure granular flows.

Measuring jumps during granular chute-flows using X-ray radiography

The transport of a colouring solute, driven by the buoyant displacement of microscopic suspended particles, and in the absence of net flow, is studied experimentally in a Hele Shaw cell. Initially, a sharp interface between a transparent fluid without particles and an underlying coloured suspension is obtained. From this situation, the suspended particles rise, carrying the solute in the form of a fingering instability across the interface, where a light transmission technique is used to measure the local solute concentration. This one attains an asymptotic value that increases with the solid fraction ϕ of suspended particles, and decreases with the distance to the interface. The solute mass discharge also increases with ϕ, always being relatively small (< 3%). The onset and development of the instability as the mechanism driving the transport of the solute is discussed.

Solute transport by suspended buoyant particles

The ability to predict thermal-induced stresses in granular materials is of practical importance across a range of disciplines ranging from process engineering to geotechnical engineering. This study presents an analytical formula to predict thermal-induced stress increments in mono-disperse granular materials subject to an initial isotropic stress state. A complementary series of DEM simulations were carried out to explore the applicability of the proposed analytical formula. The comparative analysis showed that the proposed expression can accurately predict stress changes in packings where there are negligible particle displacements as a consequence of the thermal loading (e.g. regular packings and medium/dense random packings); however large errors were observed in loose samples with a random packing.

Analytical and DEM studies of thermal stress in granular materials

The numerical modelling of geotechnical problems often poses major challenges when large displacements and strain localization are involved. Conventional continuum mechanical approaches like the finite element method (FEM) or the finite difference method (FDM) suffer from mesh distortion and numerical inaccuracy when large deformations are involved. In addition, they require the use of appropriate constitutive models to simulate the soil behaviour. The distinct element method is a promising alternative for large deformation analyses. It does not have the limitations resulting from the numerical discretization of the continuum- and not need a constitutive model since the macroscopic response results from the individual particle interaction. However, the maximum number of particles and therefore the domain of the simulation is nowadays limited by the available computational capacity. To overcome this limitation, a coupled DEM-FDM approach is proposed used to optimize the number of particles for a combined numerical domain consisting of areas of large and small displacements. The performance of the coupled DEM-FDM approach is investigated by simulating cone penetration tests in coarse grained soils.

Coupled DEM-FDM simulation of cone penetration tests in coarse grained soils

Using a discrete-element approach and a bonding interaction law, we model and test crushable inherently anisotropic structures reminiscent of the layering found in sedimentary and metamorphic rocks. By systematically modifying the level of inherent anisotropy, we characterize the evolution of the failure strength of circular rock samples discretized using a modified Voronoi tesselation under diametral point loading at different orientations relative to the sample’s layers. We characterize the failure strength, which can dramatically increase as the loading becomes orthogonal to the rock layers. We also describe the evolution of the macroscopic failure modes as a function of the loading orientation and the energy consumption at fissuring. Our simulation strategy lets us conclude that the length of bonds between Voronoi cells controls the energy being consumed in fissuring the rock sample, although failure modes and strength are considerably changing. We end up this work showing that the microstructure is largely affected by the level of inherent anisotropy and loading orientation.

Strength and energy consumption of inherently anisotropic rocks at failure

In this presentation a general concept of frictional heat and the associated rise in contact temperature is discussed, the theoretical/analytical models for predicting the rise in contact temperature are presented. Similarly, a finite element method model is discussed. Results for a maximum contact temperature for a friction slide between two spherical particles show that the TK model predictions are in good agreement with the Finite element analysis.

Friction-induced heat generation between two particles

The objective of this work was to measure the coefficient of restitution of fine granite sand dust (passing sieve ASTM N°30) and to study its influence on the prediction of classification parameters in an elutriator device. Three types of restitution coefficients were evaluated in the computational fluid dynamic model: experimental, elastic, and half the experimental value. The influence of process variables, such as the inclination angle of the equipment and the air inlet velocity, on the efficiency of the obtained fine and coarse fractions was also studied. The results indicate that the coefficient of restitution significantly influences the efficiency parameters calculated for the inclined tests, given the greater contact that exists between the particles and the equipment wall. The experimental value of the coefficient achieved almost total congruence with the coarse fraction collected for the tested velocities. For the vertical tests, the variation of the coefficient did not show significant differences for both collected fractions. The results obtained in this work allowed to gain a deeper insight into the studied gas-solid systems, and to improve the numerical predictions on fine particles classification processes widely used in the ceramic, mining and metallurgical industries.

Influence of the coefficient of restitution on the classification of granite sand using an elutriator

The use of the Discrete Element Method to model engineering structures implementing granular materials has proven to be an ecient method to response under various behaviour conditions. However, the computational cost of the simulations increases rapidly, as the number of particles and particle shape complexity increases. An aordable solution to render problems computationally tractable is to use graphical processing units (GPU) for computing. Modern GPUs oer up 10496 compute cores, which allows for a greater parallelisation relative to 32-cores oered by high-end Central Processing Unit (CPU) compute. This study outlines the application of BlazeDEM-GPU, using an RTX 2080Ti GPU (4352 cores), to investigate the influence of the modelling of particle shape on the lateral pull behaviour of granular ballast systems used in railway applications.
The idea is to validate the model and show the benefits of simulating non-spherical shapes in future large-scale tests. The algorithm, created to generate the shape of the ballast based on real grain scans, and using polyhedral shape approximations of varying degrees of complexity is shown. The particle size is modelled to scale. A preliminary investigation of the effect of the grain shape is conducted, where a sleeper lateral pull test is carried out in a spherical grains sample, and various polyhedral grains samples.
Results show that elementary polyhedral shape representations recreate some of the characteristic responses in the lateral pull test,
such as stick/slip phenomena and force chain distributions, which looks promising for future works on railway simulations.
These responses cannot be recreated with simple spherical grains, even when heuristics are added, which requires additional calibration and approximations.
The significant reduction in time when using non-spherical grains also implies that larger granular systems can be investigated.

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

Downloads

Next from Powders & Grains 2021

Grid-algorithm improvements for dense suspensions of discrete element particles in finite element fluid simulations

Similar lecture

An average model for disordered sphere packings

Stay up to date with the latest Underline news!

PRESENTATIONS

CONFERENCES

COMPANY

RESOURCES

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

.css-70qvj9{display:-webkit-box;display:-webkit-flex;display:-ms-flexbox;display:flex;-webkit-align-items:center;-webkit-box-align:center;-ms-flex-align:center;align-items:center;}Downloads

Next from Powders & Grains 2021

Grid-algorithm improvements for dense suspensions of discrete element particles in finite element fluid simulations

Similar lecture

An average model for disordered sphere packings

Stay up to date with the latest Underline news!

PRESENTATIONS

CONFERENCES

COMPANY

RESOURCES

Downloads