Implementation of numerically exact Coulomb friction for discrete element simulations in three dimensions

United States

As a follow-up of an earlier work on the numerically exact Coulomb friction in two-dimensional simulations, we present here the relations and implementation for three-dimensional discrete element particles.

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

We investigate the avalanches of spherical and non-spherical granular particles inside half-filled rotating drums. The time series of the center of gravity of the particle assemblies are obtained via image analysis and their single-sided amplitude (SSA) spectra are analyzed. The spectra features of this new indicator turn out to be characteristic for the avalanches, in terms of the existence of peaks in the low-frequency range and the decay rate of high frequency components. The SSA spectrum has a peak for the packings of non-spherical particles but not for the spherical particles. The high frequency part is characterized by a power law decay { $1/f^a \ (a>0)$ }. A $1/f$-decay is found only for the spherical particles. For the packings of cornered particles, the exponents significantly deviate from $a=1$. As $1/f$ spectra are often associated with self-organized criticality and therefore a scale invariance of the dynamics, we may conclude that there is no scale-invariant structure for granular avalanches. Considering the small number of particles and the regularity of convex particle shapes being used, the spectral features revealed in this study could be utilized for validating particle simulations.

Experimental study of particle shape dependence of avalanches inside a rotating drum

The present contribution tackles the issue of incipient conditions for initiation of erosion by a fluid flow at the surface of cohesive materials. To this end, a typical assessment procedure consists in subjecting a soil sample to progressive hydrodynamic stresses induced by a submerged impinging jet flow whose injection velocity is gradually increased. This paper presents the results of an extensive use of this protocol both in experiments and numerical simulations, the latter being based on a coupled DEM and LBM approach. Here we consider the specific case of weakly cemented soils, either made experimentally of glass beads bonded by solid bridges or modelled numerically by a solid bond rheology with a parabolic yield condition involving the micromechanical traction, shearing and bending of the bonds. The results show that, as expected, the hydrodynamic stress for erosion onset substantially increases with solid cohesion as compared to cohesion-less cases but can, however, be satisfactorily predicted by a simple extension of the usual Shields criterion that only applies for cohesion-less granular sediments. This extension includes a cohesion number, the granular Bond number, with a simple definition based on tensile yield values.

Extending the Shields criterion to erosion of weakly cemented granular soils

The unknown particle scale mechanics of slow granular flows have hindered our understanding of its features, such as rate independence. The DEM simulations unveil the signatures of stick-slip pattern during flow. Furthermore, we analyzed the time evolution of contact networks and constructed a cascade failure mechanism - from the isostatic criterion - which demonstrates the emergence of stick-slip motion. This insight into micromechanics provides phenomenological reasoning for rate independence and large fluctuations in stress.

Micromechanics of slow granular flows

In this study, we investigate the flow rate of model granular material in a silo under the influence of mechanical vibrations. Experimental measurements and discrete element simulations (DEM) are performed in a quasi-2D silo. The influence on the flow rate of the opening size and the vibration applied on the entire silo is studied. Two distinct regimes are evidenced, governed by the Froude number Fr and the relative frequency. In the first regime, a decreased flow rate is observed when increasing the vibration intensity. This behavior is explained by the presence of reorganizations induced by the vibration, leading to a more homogeneous but also slower flow. In the second regime, an increased flow rate is evidenced when increasing the vibration intensity. We find this behavior comes from the intermittent nature of the flow, where the flow rate is directly controlled by the propagation of shock waves all along the silo.

Discharge of a granular silo under mechanical vibrations

Wave velocity and attenuation are two of the most important attributes of geomaterials under seismic loads. To understand how they are influenced by particle and microstructural characteristics, Discrete Elements Method simulations are carried out to simulate the propagation of linear and nonlinear waves in isotropic, frictional granular packings. The focus is on the effects of confining pressure and source characteristics (i.e., input strain and frequency) on wave velocity and attenuation. Three different approaches are used to quantify the seismic quality factor Q or damping in time and frequency domain. Our simulation results verify previous findings on the scaling of wave velocity with respect to confining pressure and coordination number. Quality factors measured both in time and frequency domains increase with confining pressure and show an interesting non-monotonic behavior with input frequency. As expected, wave velocity and quality factor for S-waves decrease with an increasing input shear strain amplitude, while for P-waves they first decrease and then increase at input strain amplitude of 6×10-4. Change of participation number shows that the inhomogeneity of the system increases with the increase of the input strain.

Elastic wave velocity and attenuation in granular material

Numerical models of granular materials are useful in tribology, and can be used to predict wear and friction in contacts. DEM-like simulations are used to model particles of third-body, which are partly wear debris from rubbing bodies. It has been shown that the third-body particles can have different flow regimes, depending on their mechanical properties. Among the different characteristics of flow regimes, agglomerate size seems to be crucial. A method based on vortex analysis used in fluid mechanics allows characterizing this cluster size. The results show that different vortex sizes can be observed during the simulation. In particular, it is observed that some vortexes of a characteristic size persist over time, and could be representative of agglomerates. These results pave the way for a better characterisation of the different flow regimes.

Detecting vorticity in cohesive deformable granular material

We studied experimentally the effect of grain size and maximum load on the compaction and subsequent relaxation of a granular column when it is subjected to vertical uniaxial compression. The experiments were performed using two different types of grains: 1) solid glass beads, and 2) porous beads that consist of agglomerates of glass powder. We found that the compression force increases non-linearly with time, with sudden drops for the case of glass beads and periodic undulations for dust particles. Whereas the grain size effect is small in the average force load, the fluctuations become larger as the grain size increases. On the other hand, the relaxation process is well described by the Maxwell model with three different relaxation time scales.

Grain size effect on the compression and relaxation of a granular column: solid particles vs dust agglomerates

The dissipation of energy during the compaction of granular materials was studied by performing
confined drop tower experiments on Ottawa sand. Energy dissipated due to breakage was quantified by evaluating the creation of new surfaces at varying drop heights. Post-compaction particle size distributions (PSD)
were measured and the amount of breakage was quantified by the position of the current PSD relative to the
pre-compaction and ultimate PSD. Our observations revealed that the percentage of input energy dissipated
due to breakage accounted for less than 0.5% of the total energy budget and was a constant proportion regardless of the total energy applied to the system. We also evaluated the effects of die wall friction by measuring
post-compaction PSD in various positions within the sample.

Energy dissipation due to breakage during confined compaction of granular materials

This presentation shows the results of five triaxial compression tests performed on approximately 9000 oblate spheroids (lentils) with different initial orientations with respect to the loading axis (0º, 30°, 45º, 60º, and 90°). Typical stress-strain information is complemented with x-ray tomography scans every 1 % strain. Starting from an initial labelling of particles, a new technique is used to obtain unprecedented detail of tracking of all the particles through time. This rich dataset is analysed from the perspective of inter-particle contacts (building on previous metrological work on this subject) revealing that the mean contact orientation in both samples rotates towards the direction of σ1 at a rate depending on the initial orientation. Different trends are observed for the alignment of contact orientation. 


Contact evolution in granular materials with inherently anisotropic fabric

Several granular materials are hydro-sensitive, i.e., contact with water severely affects their morphology and mechanical behaviour. To this broad class belong a series of materials of great relevance for the pharmaceutical and food industry. Former studies have been conducted on the effect that moisture or humidity have on granular flow or on the individual grain mechanical response, but the processes occurring at the micro-level and their influence on the overall granular packing behaviour is yet to be fully understood. This study presents an experimental investigation of the response of a water-sensitive material (couscous) exposed to high relative humidity (97%). 4D (3D plus time) x-ray tomographies were acquired in operando. A data treatment approach based on Dicrete Digital Image Correlation (dDIC) is developed and detailed here. This allows for following each individual grain throughout the 4+ days of the humidification test. This, in turn, allows the study of the inter- and intra-granular strain, which is particularly pronounced in water sensitive materials. The expansion of the individual grains also severely affects the contact network, whose evolving properties are here analysed and correlated to the macroscopic (sample-scale) deformation.

Implementation of numerically exact Coulomb friction for discrete element simulations in three dimensions

Downloads

Next from Powders & Grains 2021

Experimental study of particle shape dependence of avalanches inside a rotating drum

Similar lecture

Investigating the effect of porosity on the soil water retention curve using the multiphase Lattice Boltzmann Method

Stay up to date with the latest Underline news!

PRESENTATIONS

CONFERENCES

COMPANY

RESOURCES

Implementation of numerically exact Coulomb friction for discrete element simulations in three dimensions

.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

Experimental study of particle shape dependence of avalanches inside a rotating drum

Similar lecture

Investigating the effect of porosity on the soil water retention curve using the multiphase Lattice Boltzmann Method

Stay up to date with the latest Underline news!

PRESENTATIONS

CONFERENCES

COMPANY

RESOURCES

Downloads