Frictional hysteresis and particle deposition in granular free-surface flows

Abstract

Shallow granular avalanches on slopes close to repose exhibit hysteretic behaviour. For instance, when a steady-uniform granular flow is brought to rest it leaves a deposit of thickness $h_{stop}(\unicode[STIX]{x1D701})$ on a rough slope inclined at an angle $\unicode[STIX]{x1D701}$ to the horizontal. However, this layer will not spontaneously start to flow again until it is inclined to a higher angle $\unicode[STIX]{x1D701}_{start}$, or the thickness is increased to $h_{start}(\unicode[STIX]{x1D701})>h_{stop}(\unicode[STIX]{x1D701})$. This simple phenomenology leads to a rich variety of flows with co-existing regions of solid-like and fluid-like granular behaviour that evolve in space and time. In particular, frictional hysteresis is directly responsible for the spontaneous formation of self-channelized flows with static levees, retrogressive failures as well as erosion–deposition waves that travel through the material. This paper is motivated by the experimental observation that a travelling-wave develops, when a steady uniform flow of carborundum particles on a bed of larger glass beads, runs out to leave a deposit that is approximately equal to $h_{stop}$. Numerical simulations using the friction law originally proposed by Edwards et al. (J. Fluid Mech., vol. 823, 2017, pp. 278–315) and modified here, demonstrate that there are in fact two travelling waves. One that marks the trailing edge of the steady-uniform flow and another that rapidly deposits the particles, directly connecting the point of minimum dynamic friction (at thickness $h_{\ast }$) with the deposited layer. The first wave moves slightly faster than the second wave, and so there is a slowly expanding region between them in which the flow thins and the particles slow down. An exact inviscid solution for the second travelling wave is derived and it is shown that for a steady-uniform flow of thickness $h_{\ast }$ it produces a deposit close to $h_{stop}$ for all inclination angles. Numerical simulations show that the two-wave structure deposits layers that are approximately equal to $h_{stop}$ for all initial thicknesses. This insensitivity to the initial conditions implies that $h_{stop}$ is a universal quantity, at least for carborundum particles on a bed of larger glass beads. Numerical simulations are therefore able to capture the complete experimental staircase procedure, which is commonly used to determine the $h_{stop}$ and $h_{start}$ curves by progressively increasing the inclination of the chute. In general, however, the deposit thickness may depend on the depth of the flowing layer that generated it, so the most robust way to determine $h_{stop}$ is to measure the deposit thickness from a flow that was moving at the minimum steady-uniform velocity. Finally, some of the pathologies in earlier non-monotonic friction laws are discussed and it is explicitly shown that with these models either steadily travelling deposition waves do not form or they do not leave the correct deposit depth $h_{stop}$.