Stable solutions of a scalar conservation law for particle-size segregation in dense granular avalanches

Abstract

Dense, dry granular avalanches are very efficient at sorting the larger particles towards the free surface of the flow, and finer grains towards the base, through the combined processes of kinetic sieving and squeeze expulsion. This generates an inversely graded particle-size distribution, which is fundamental to a variety of pattern formation mechanisms, as well as subtle size-mobility feedback effects, leading to the formation of coarse-grained lateral levees that create channels in geophysical flows, enhancing their run-out. In this paper we investigate some of the properties of a recent model [Gray, J. M. N. T. & Thornton, A. R. (2005) A theory for particle size segregation in shallow granular free-surface flows. Proc. R. Soc. 461, 1447–1473]; [Thornton, A. R., Gray, J. M. N. T. & Hogg, A. J. (2006) A three-phase mixture theory for particle size segregation in shallow granular free-surface flows. J. Fluid. Mech. 550, 1–25] for the segregation of particles of two sizes but the same density in a shear flow typical of shallow avalanches. The model is a scalar conservation law in space and time, for the volume fraction of smaller particles, with non-constant coefficients depending on depth within the avalanche. It is proved that for steady flow from an inlet, complete segregation occurs beyond a certain finite distance down the slope, no matter what the mixture at the inlet. In time-dependent flow, dynamic shock waves can develop; they are interfaces separating different mixes of particles. Shock waves are shown to be stable if and only if there is a greater concentration of large particles above the interface than below. Constructions with shocks and rarefaction waves are demonstrated on a pair of physically relevant initial boundary value problems, in which a region of all small particles is penetrated from the inlet by either a uniform mixture of particles or by a layer of small particles over a layer of large particles. In both cases, and under a linear shear flow, solutions are constructed for all time and shown to have similar structure for all choices of parameters.