A multilayer shallow model for dry granular ows with the (I)-rheology: application to granular collapse on erodible beds

Abstract

In this work we present a multilayer shallow model to approximate the Navier–Stokes equations with the μ(I)μ(I)-rheology through an asymptotic analysis. The main advantages of this approximation are (i) the low cost associated with the numerical treatment of the free surface of the modelled flows, (ii) the exact conservation of mass and (iii) the ability to compute two-dimensional profiles of the velocities in the directions along and normal to the slope. The derivation of the model follows Fernández-Nieto et al. (J. Comput. Phys., vol. 60, 2014, pp. 408–437) and introduces a dimensional analysis based on the shallow flow hypothesis. The proposed first-order multilayer model fully satisfies a dissipative energy equation. A comparison with steady uniform Bagnold flow – with and without the sidewall friction effect – and laboratory experiments with a non-constant normal profile of the downslope velocity demonstrates the accuracy of the numerical model. Finally, by comparing the numerical results with experimental data on granular collapses, we show that the proposed multilayer model with the μ(I)μ(I)-rheology qualitatively reproduces the effect of the erodible bed on granular flow dynamics and deposits, such as the increase of runout distance with increasing thickness of the erodible bed. We show that the use of a constant friction coefficient in the multilayer model leads to the opposite behaviour. This multilayer model captures the strong change in shape of the velocity profile (from S-shaped to Bagnold-like) observed during the different phases of the highly transient flow, including the presence of static and flowing zones within the granular column.