Contact Dynamics modeling of viscoelastic granular materials using irregular polyhedral particles

United States

Viscoelastic granular materials are present in several disciplines. One example is asphalt mixture employed in road construction. In the last three decades, discrete element modeling has been positioned as a valid tool for the analysis of this multiphase material at the grain-scale. All this despite the simplification of the
shape of the particles used in these studies. In this work, it is proposed a simplified procedure for the generation of viscoelastic granular samples composed of irregular polyhedra. The numerical aggregates were generated by a Poisson-Voronoi tessellation based on the particle size distribution (PSD) and statistic data of aggregates, without using complex imaging technics. This procedure set the porosity of the packing, while controlling the PSD. Using this procedure implies a significant computational-time reduction by skipping several preparation stages for polyhedral samples, such as deposition by gravity and compaction. This approach can be used for the study granular materials as inclusion in a solid matrix as concrete or asphalt mixtures, particle breaking, and fatigue damage of viscoelastic materials.

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.

poster

Copper mining produces significant amounts of associate waste, leading to the construction of huge mining waste dumps in the highly seismic zone of the Andes Mountains. To analyze the seismic stability of these deposits, it is necessary to characterize the mechanical properties of these materials in laboratory, which is quite difficult due to the large size of particles composing these structures. In this study, the global stability of mine waste dumps is analyzed by means of discrete element simulations. Three analogues upscaled particle size distributions were generated from the original PSD. The particles composing the dump were modeled using regular polygons. The modeled dumps were subject to a seismic loading by applying a velocity signal on the boundary walls, which was integrated directly from accelerometers database. The results show that failure in these structures is produced by infinite slope, for samples scaled by a factor until ten. This methodology is encouraging for the analysis of the global stability of mine waste dumps and for forecasting the runout distance after failure.

Discrete modeling of waste rock dumps stability under seismic loading

To understand the typically heterogeneous flowing behavior of granular materials, it is important to combine flow tests with three-dimensional imaging. To probe the flow behavior of granular materials over a wide range of flow rates, it is imperative to be able to impose such flow rates in a well controlled manner while performing imaging tests that are compatible with all imposed flow rates. Achieving both flow control and bulk imaging capacity is challenging for a number of reasons. Here, we describe the design of a setup in which we are able to do imaging while imposing a constant overall shear rate on a granular material. We characterize the setup in which flow tests will be performed, which consists of a bottom-driven cone-plate or double-cone design. We show that the setup can be integrated in x-ray microtomography devices to aid particle tracking based flow measurements. The design is also compatible with typical rheometer setups. We also perform high
speed imaging of a granular flow in an ultra-fast x-ray scanner, for which we provide proof-of-principle data in a simplified shear setup. The designed flow geometry is also compatible with said high speed imaging facility, where particle image velocimetry can be employed to extract quantitative flow field data.

Improved evaluation of granular media flows using an X-ray scanning compatible cone-plate setup

Granular materials can segregate spontaneously due to differences in particle properties when subjected to vibrations, shear strain or because of the equipment geometries. Although the difference in particle size is the most critical factor that drives segregation, the effects of large density difference may also be detrimental for a lot of industries. In this work, we experimentally investigate density-driven segregation in bi-disperse mixtures of particles having the same size but different density when subjected to non-uniform shear rates. We found that the features of the segregation process are related to the density ratio as well as to the dimensionless loaded mass. The experimental outcomes are then compared with the solution of a simple density-driven segregation model. The model can successfully capture the main features of segregation driven by density for a range of density ratios.

Shear-driven density segregation: an experimental study

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

Some years ago, Harth et al. experimentally explored the steady state dynamics of a heated granular gas of rod-like particles in microgravity [K. Harth et al. Phys. Rev. Lett. 110, 144102 (2013)]. Here, we report numerical results that quantitatively reproduce their experimental findings and provide additional insight into the process. A system of sphero-cylinders is heated by the vibration of three flat side walls, resulting in one symmetrically heated direction, one non-symmetrically heated direction, and one non-heated direction. In the non-heated direction, the speed distribution follows a stretched exponential distribution p(v) ∝ exp(−(|v|/C)^1.5 ).
In the symmetrically heated direction, the velocity statistics at low speeds is similar but it develops pronounced exponential tails at high speeds. In the non-symmetrically heated direction (not accessed experimentally), the distribution also follows p(v) ∝ exp(−(|v|/C)^1.5 ), but the velocity statistics of rods moving toward the vibrating wall resembles the indirectly excited direction, whereas the velocity statistics of those moving away from the wall resembles the direct excited direction.

Continuously heated granular gas of elongated particles

Poroelasticity theory predicts wave velocities in a saturated porous medium through a coupling between the bulk deformation of the solid skeleton and porous fluid flow. The challenge emerges below the characteristic wavelengths at which hydrodynamic interactions between grains and pore fluid become important. We investigate the pressure and volume fraction dependence of compressional- and shear-wave velocities in fluid-saturated, random, isotropic, frictional granular packings. The lattice Boltzmann method (LBM) and discrete element method (DEM) are two-way coupled to capture the particle-pore fluid interactions; an acoustic source is implemented to insert a traveling wave from the fluid reservoir to the saturated medium. We extract wave velocities from the acoustic branches in the wavenumber-frequency space, for a range of confining pressures and volume fractions. For random isotropic granular media, the pressure-wave velocity data collapse on a single curve when scaled properly by the volume fraction.

Direct numerical simulation of wave propagation in saturated random granular packings using coupled LBM-DEM

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

Tsunami-like waves can arise from large-scale geological events such as the collapse of a mountain flank or a volcano into a lake or an ocean. We here mimic this situation at the laboratory scale. The landslide is modeled by the sudden release of a rectangular granular column, which impacts a still water layer and generates a wave. Through experiments varying both the column dimensions and the initial water depth, three regimes of nonlinear waves of different shapes are reported. The amplitude and width of the generated waves are governed by a local Froude number based on the velocity of the granular front at the water surface.

On water waves generated by gravity driven granular collapse

One striking difference between aggregates of flexible frictional fibres and other granular materials like rigid spheres is the effective cohesion of their assembly. While glue or capillary bridges are needed to shape aggregates of spherical particles and build sandcastles, for fibres, no need for glue to build a nest. Here we study an assembly of mono disperse flexible fibres. We first use X-ray microtomography to characterise the geometry of the initial assembly, the number of contact points and mean curvatures of the fibres. Using forcedisplacement measurements, we characterise the macroscopic cohesive strength of the aggregate by varying the geometry of the fibres, the fibres mechanicals properties and the packing of the preparation. Finally, we relate the macroscopic mechanical behaviour of the assembly with the filament reorganisation at the microscopic scale.

Cohesion of bird nests

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

Contact Dynamics modeling of viscoelastic granular materials using irregular polyhedral particles

Downloads

Next from Powders & Grains 2021

Discrete modeling of waste rock dumps stability under seismic loading

Similar lecture

Revisiting Hardin’s parameters for the quantification of particle breakage – A statistical entropy approach