Continuum viscoplastic simulation of a granular column collapse on large slopes: μ(I) rheology and lateral wall effects

We simulate here dry granular flows resulting from the collapse of granular columns on an inclined channel (up to 22°) and compare precisely the results with laboratory experiments. Incompressibility is assumed despite the dilatancy observed in the experiments (up to 10%). The 2-D model is based on the so-called μ(I) rheology that induces a Drucker-Prager yield stress and a variable viscosity. A nonlinear Coulomb friction term, representing the friction on the lateral walls of the channel, is added to the model. We demonstrate that this term is crucial to accurately reproduce granular collapses on slopes ≳10°, whereas it remains of little effect on the horizontal slope. Quantitative comparison between the experimental and numerical changes with time of the thickness profiles and front velocity makes it possible to strongly constrain the rheology. In particular, we show that the use of a variable or a constant viscosity does not change significantly the results provided that these viscosities are of the same order. However, only a fine tuning of the constant viscosity (η=1 Pa s) makes it possible to predict the slow propagation phase observed experimentally at large slopes. Finally, we observed that small-scale instabilities develop when refining the mesh (also called ill-posed behavior, characterized in the work of Barker et al. [“Well-posed and ill-posed behaviour of the μ(I)-rheology for granular flow,” J. Fluid Mech. 779, 794–818 (2015)] and in the present work) associated with the mechanical model. The velocity field becomes stratified and the bands of high velocity gradient appear. These model instabilities are not avoided by using variable viscosity models such as the μ(I) rheology. However we show that the velocity range, the static-flowing transition, and the thickness profiles are almost not affected by them.