Development of a single-parameter spring-dashpot rolling friction model for coarse-grained DEM

TL;DR

This paper introduces a simplified single-parameter spring-dashpot rolling friction model for coarse-grained DEM, reducing calibration complexity while accurately capturing particle behavior in large-scale granular simulations.

Contribution

A novel single-parameter rolling friction model based on the critical rolling angle is developed, simplifying implementation and calibration in DEM simulations.

Findings

Model achieves stable equilibrium without oscillations.

Successfully reproduces macroscopic behavior in DEM-CFD simulations.

Reduces calibration effort by using only one physically meaningful parameter.

Abstract

Simulating granular materials composed of non-spherical particles remains a major challenge in discrete element method (DEM) simulations due to the complexity of contact detection and rotational dynamics, rendering large-scale simulations computationally prohibitive. To address this limitation, rolling friction is commonly introduced as an approximation to account for particle shape effects by applying a resistive torque to spherical particles. Among existing rolling friction formulations, the spring-dashpot (S-D) type model is widely recognized for its numerical stability and realistic representation of rolling resistance. However, conventional S-D models require multiple empirical parameters that must be calibrated in an interdependent manner, leading to increased experimental effort, parameter ambiguity, and uncertainty in practical applications. To overcome these issues, this study proposes a new S-D type rolling friction model that reduces the parameter set to a single physically meaningful quantity: the critical rolling angle. Derived from theoretical considerations, this parameter characterizes the transition from static to rolling motion at particle contacts. The use of a single parameter simplifies implementation and eliminates the need for extensive calibration. Stability analysis demonstrates that the proposed model allows particles to reach a physically consistent equilibrium state without spurious rotational oscillations. For large-scale applications, the model is further integrated into a coarse-grained DEM framework. Validation using DEM-CFD simulations of an incinerator system confirms that the proposed approach successfully reproduces the macroscopic behavior of the original particle system. Overall, this study enhances the applicability of DEM for industrial-scale simulations involving non-spherical particles.