Granular thermostat implementation within the soft-sphere Discrete Element Method (DEM) framework, considerations and limitations
Abstract
Dense granular flows are commonly investigated using the soft-sphere Discrete Element Method (DEM), whereas large-scale applications generally require continuum models. Granular temperature, defined as the variance of particle velocity fluctuations, is a relevant variable to collapse rheological scaling. However, experimental forcing approaches used to control temperature also affects other aspects of the granular state. An alternative is to use a thermostat algorithm to control the fluctuation energy directly. Although thermostat algorithms are well established in Molecular Dynamics (MD), their behaviour in dissipative DEM systems has received comparatively little attention. We show that conventional Langevin and Nosé-Hoover thermostats exhibit complementary limitations: Langevin control requires strong coupling to offset collisional dissipation, thereby damping particle dynamics, whereas Nosé-Hoover does not independently disrupt the correlations and segregation generated by repeated inelastic collisions, leading to non-ergodic and numerically unstable states. To address these limitations, we introduce two pairwise hybrid formulations combining deterministic temperature regulation with stochastic decorrelation. Both enforce the prescribed temperature, while variation of the stochastic decorrelation timescale modifies velocity statistics and associated mesoscopic properties such as diffusivity. Finally, application to pressure-controlled simple shear shows that increasing granular temperature at fixed inertial number reduces the apparent friction, consistent with previously reported trends. The framework therefore provides controlled reference states for comparing differently forced granular systems and identifying the variables required for temperature-dependent constitutive models.
Related
- Mechanistic transition between inter- and trans-granular creep cracking via a unified crystal plasticity-fracture framework
- The Sign Premise Is a State Clause: Pumped Baths, the de Sitter Thermostat, and the Limits of Scale Without Law
- Disentangling universality from shape in Platonic granular packings
- Rheology of dense vibrated granular flows: Nonmonotonic response controlled by granular temperature
- Dynamical and structural properties of an absorbing phase transition: a case study from granular systems
BibTeX
@misc{previtali2026granular,
title = {{Granular thermostat implementation within the soft-sphere Discrete Element Method (DEM) framework, considerations and limitations}},
author = {Marco Previtali and Herbert Eric Huppert and Sergio Andres Galindo Torres},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.21297},
eprint = {2609.21297},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.21297},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Granular thermostat implementation within the soft-sphere Discrete Element Method (DEM) framework, considerations and limitations
AU - Marco Previtali
AU - Herbert Eric Huppert
AU - Sergio Andres Galindo Torres
PY - 2026
DA - 2026/09/18
PB - arXiv
DO - 10.48550/arxiv.2609.21297
UR - https://arxiv.org/abs/2609.21297
AB - Dense granular flows are commonly investigated using the soft-sphere Discrete Element Method (DEM), whereas large-scale applications generally require continuum models. Granular temperature, defined as the variance of particle velocity fluctuations, is a relevant variable to collapse rheological scaling. However, experimental forcing approaches used to control temperature also affects other aspects of the granular state. An alternative is to use a thermostat algorithm to control the fluctuation energy directly. Although thermostat algorithms are well established in Molecular Dynamics (MD), their behaviour in dissipative DEM systems has received comparatively little attention. We show that conventional Langevin and Nosé-Hoover thermostats exhibit complementary limitations: Langevin control requires strong coupling to offset collisional dissipation, thereby damping particle dynamics, whereas Nosé-Hoover does not independently disrupt the correlations and segregation generated by repeated inelastic collisions, leading to non-ergodic and numerically unstable states. To address these limitations, we introduce two pairwise hybrid formulations combining deterministic temperature regulation with stochastic decorrelation. Both enforce the prescribed temperature, while variation of the stochastic decorrelation timescale modifies velocity statistics and associated mesoscopic properties such as diffusivity. Finally, application to pressure-controlled simple shear shows that increasing granular temperature at fixed inertial number reduces the apparent friction, consistent with previously reported trends. The framework therefore provides controlled reference states for comparing differently forced granular systems and identifying the variables required for temperature-dependent constitutive models.
ER -