Shock waves, dead zones and particle-free regions form when a thin surface avalanche of granular material flows around an obstacle or over a change in the bed topography. Understanding and modelling these flows is of considerable practical interest for industrial processes, as well as for the design of defences to protect buildings, structures and people from snow avalanches, debris flows and rockfalls. These flow phenomena also yield useful constitutive information that can be used to improve existing avalanche models. In this paper a simple hydraulic theory, first suggested in the Russian literature, is generalized to model quasi-two-dimensional flows around obstacles. Exact and numerical solutions are then compared with laboratory experiments. These indicate that the theory is adequate to quantitatively describe the formation of normal shocks, oblique shocks, dead zones and granular vacua. Such features are generated by the flow around a pyramidal obstacle, which is typical of some of the defensive structures in use today.

Shock formations are observed in granular avalanches when supercritical flow merges into a region of subcritical flow. In this paper we employ a shock-capturing numerical scheme for the one-dimensional Savage-Hutter theory of granular flow to describe this phenomenon. A Lagrangian moving mesh scheme applied to the nonconservative form of the equations reproduces smooth solutions of these free boundary problems very well, but fails when shocks are formed. A nonoscillatory central (NOC) difference scheme with TVD limiter or WENO cell reconstruction for the conservative equations is therefore introduced. For the avalanche free boundary problems it must be combined with a front-tracking method, developed here, to properly describe the margin evolution. It is found that this NOC scheme combined with the front-tracking module reproduces both the shock wave and the smooth solution accurately. A piecewise quadratic WENO reconstruction improves the smoothness of the solution near local extrema. The schemes are checked against exact solutions for (1) an upward moving shock wave, (2) the motion of a parabolic cap down an inclined plane, and (3) the motion of a parabolic cap down a curved slope ending in a flat run-out region, where a shock is formed as the avalanche comes to a halt. c 2002 Elsevier Science

No abstract

We use the non-Cartesian, topography-based equations of mass and momentum balance for gravity driven frictional flows of Luca et al. (Math. Mod. Meth. Appl. Sci. 19, 127-171 (2009)) to motivate a study on various approximations of avalanche models for single-phase granular materials. By introducing scaling approximations we develop a hierarchy of model equations which differ by degrees in shallowness, basal curvature, peculiarity of constitutive formulation (non-Newtonian viscous fluids, Savage-Hutter model) and velocity profile parametrization. An interesting result is that differences due to the constitutive behaviour are largely eliminated by scaling approximations. Emphasis is on avalanche flows; however, most equations presented here can be used in the dynamics of other thin films on arbitrary surfaces.

No abstract