skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Chrono::GPU: An Open-Source Simulation Package for Granular Dynamics Using the Discrete Element Method
We report on an open-source, publicly available C++ software module called Chrono::GPU, which uses the Discrete Element Method (DEM) to simulate large granular systems on Graphics Processing Unit (GPU) cards. The solver supports the integration of granular material with geometries defined by triangle meshes, as well as co-simulation with the multi-physics simulation engine Chrono. Chrono::GPU adopts a smooth contact formulation and implements various common contact force models, such as the Hertzian model for normal force and the Mindlin friction force model, which takes into account the history of tangential displacement, rolling frictional torques, and cohesion. We report on the code structure and highlight its use of mixed data types for reducing the memory footprint and increasing simulation speed. We discuss several validation tests (wave propagation, rotating drum, direct shear test, crater test) that compare the simulation results against experimental data or results reported in the literature. In another benchmark test, we demonstrate linear scaling with a problem size up to the GPU memory capacity; specifically, for systems with 130 million DEM elements. The simulation infrastructure is demonstrated in conjunction with simulations of the NASA Curiosity rover, which is currently active on Mars.  more » « less
Award ID(s):
1835674
PAR ID:
10349475
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Processes
Volume:
9
Issue:
10
ISSN:
2227-9717
Page Range / eLocation ID:
1813
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Hambleton, J. P. (Ed.)
    Discrete element method (DEM) has been widely applied to simulate granular soil behavior. However, traditional DEM uses sphere clusters to approximate realistic particles, which is computationally demanding when simulating many particles. This study explores the use of physics engine, a platform developed for simulating physical processes in video games, to simulate realistic particles. This paper compares realistic particle simulation methodologies using physics engine and discrete element method, including contact models, parameter settings, computational speeds, and simulation results. The results show that the physics engine and DEM achieve similar simulation outputs, while the physics engine runs significantly faster than DEM, because PhysX uses both CPUs (central processing units) and GPUs (graphics processing units) of computers, triangular face tessellations to represent realistic particles, and a simplified contact model to accelerate simulations. This study provides geo-mechanicians and DEM modelers with one more option for them to consider when they simulate realistic particles. 
    more » « less
  2. Abstract The soil contact model (SCM) is widely used in practice for off‐road wheeled vehicle mobility studies when simulation speed is important and highly accurate results are not a main concern. In practice, the SCM parameters are obtained via a bevameter test, which requires a complex apparatus and experimental procedure. Here, we advance the idea of running a virtual bevameter test using a high‐fidelity terramechanics simulation. The latter employs the “continuous representation model” (CRM), which regards the deformable terrain as an elasto‐plastic continuum that is spatially discretized using the smoothed particle hydrodynamics method. The approach embraced is as follows: a virtual bevameter test is run in simulation using CRM terrain to generate “ground truth” data; in a Bayesian framework, this data is subsequently used to calibrate the SCM terrain. We show that (i) the resulting SCM terrain, while leading to fast terramechanics simulations, serves as a good proxy for the more complex CRM terrain; and (ii) the SCM‐over‐CRM simulation speedup is roughly one order of magnitude. These conclusions are reached in conjunction with two tests: a single wheel test, and a full rover simulation. The SCM and CRM simulations are run in the open‐source software Chrono. The calibration is performed using PyMC, which is a Python package that interactively communicates with Chrono to calibrate SCM. The models and scripts used in this contribution are available as open source for unfettered use and distribution in a public repository. 
    more » « less
  3. Abstract Frictional sliding along grain boundaries in brittle shear zones can result in the fragmentation of individual grains, which ultimately can impact slip dynamics. During deformation at small scales, stick–slip motion can occur between grains when existing force chains break due to grain rearrangement or failure, resulting in frictional sliding of granular material. The rearrangement of the grains leads to dilation of the granular package, reducing the shear stress and subsequently leading to slip. Here, we conduct physical experiments employing HydroOrbs, an elasto-plastic material, to investigate grain comminution in granular media under simple shear conditions. Our findings demonstrate that the degree of grain comminution is dependent on both the normal force and the size of the grains. Using the experimental setup, we benchmark Discrete Element Method (DEM) numerical models, which are capable of simulating the movement, rotation, and fracturing of elasto-plastic grains subjected to simple shear. The DEM models successfully replicate both grain comminution patterns and horizontal force fluctuations observed in our physical experiments. They show that increasing normal forces correlate with higher horizontal forces and more fractured grains. The ability of our DEM models to accurately reproduce experimental results opens up new avenues for investigating various parameter spaces that may not be accessible through traditional laboratory experiments, for example, in assessing how internal friction or cohesion affect deformation in granular systems. 
    more » « less
  4. null (Ed.)
    Flowing granular materials segregate due to differences in particle size (driven by percolation) and density (driven by buoyancy). Modelling the segregation of mixtures of large/heavy particles and small/light particles is challenging due to the opposing effects of the two segregation mechanisms. Using discrete element method (DEM) simulations of combined size and density segregation we show that the segregation velocity is well described by a model that depends linearly on the local shear rate and quadratically on the species concentration for free surface flows. Concentration profiles predicted by incorporating this segregation velocity model into a continuum advection–diffusion–segregation transport model match DEM simulation results well for a wide range of particle size and density ratios. Most surprisingly, the DEM simulations and the segregation velocity model both show that the segregation direction for a range of size and density ratios depends on the local species concentration. This leads to a methodology to determine the combination of particle size ratio, density ratio and particle concentration for which a bidisperse mixture will not segregate. 
    more » « less
  5. In this study, deformation and failure mechanisms of mineralized tissue (bone) were investigated both experimentally and computationally by performing diametral compression tests on millimetric disk specimens and conducting finite element analysis in which a granular micromechanics-based nonlinear user-defined material model is implemented. The force–displacement relationship obtained in the simulation agreed well with the experimental results. The simulation was also able to capture location of the failure initiation observed in the experiment, which is inside out from the hole along the loading axis. Furthermore, propagation of micro-sized cracks into failure was observed both in the experiment using simultaneous slow-motion microscopy imaging and in the simulation analyzing the local distortion and local volume change within the specimen. The anisotropy evolution was found to be significant around the hole along the loading axis by evaluating the anisotropy index computed using finite element results. In conclusion, this work revealed that the prediction capability of granular micromechanics-based user-defined nonlinear material model (UMAT) is promising considering the match between the results and observations from the physical experiment and finite element analysis such as force–displacement relationship and failure initiation/pattern. This work has also shown that the tensile damage and failure of mineralized tissues can be characterized using diametral compression (split tension) test. 
    more » « less