Computational study of erosion and deposition in transient granular flows on an erodible heap
Phys. Rev. E 112, 055408 – Published 10 November, 2025
DOI: https://doi.org/10.1103/rrcf-3j5f
Abstract
We study, using Discrete Element Method (DEM) simulations, the transient response of an initially steady surface granular flow on a heap of the same material, when the heap is tilted by a small, positive or negative, angle. The system does not have sidewalls, unlike many previous studies. Upon application of the tilt, the height of the free surface, the height of the interface between the flowing and static zones, and the angle of the interface vary exponentially with time while approaching a steady state. A parametric study is conducted considering different mass flow rates, interparticle friction coefficients, particle diameters, and tilt angles. Local rates of erosion and deposition (), calculated at the interface, follow predictions of the models of Bouchaud et al. [J. de Phys. I 4, 1383 (1994)] and Boutreux et al. [Phys. Rev. E 58, 4692 (1998)] with where is the erosion velocity, is the angle of the interface, and is the neutral angle at which there is no erosion or deposition. The erosion velocity is significantly larger than the deposition velocity for all the cases studied, with . The erosion velocities are nearly independent of the mass flow rate in the layer, the interparticle friction, and the tilt angle. The velocities increase with particle diameter and flowing layer thickness as .