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.

[1]  J. Hadamard Sur les problemes aux derive espartielles et leur signification physique , 1902 .

[2]  P. Perzyna Fundamental Problems in Viscoplasticity , 1966 .

[3]  J. Lions,et al.  Les inéquations en mécanique et en physique , 1973 .

[4]  Wing Kam Liu,et al.  Lagrangian-Eulerian finite element formulation for incompressible viscous flows☆ , 1981 .

[5]  S. Savage,et al.  The Mechanics of Rapid Granular Flows , 1984 .

[6]  R. Glowinski,et al.  Augmented Lagrangian and Operator-Splitting Methods in Nonlinear Mechanics , 1987 .

[7]  T. Papanastasiou Flows of Materials with Yield , 1987 .

[8]  David G. Schaeffer,et al.  Instability in the evolution equations describing incompressible granular flow , 1987 .

[9]  E. Bruce Pitman,et al.  Ill‐posedness in three‐dimensional plastic flow , 1988 .

[10]  B. Maury Characteristics ALE Method for the Unsteady 3D Navier-Stokes Equations with a Free Surface , 1996 .

[11]  C. Bellone,et al.  ON A ROLE , 1996 .

[12]  P. Coussot,et al.  A theoretical framework for granular suspensions in a steady simple shear flow , 1999 .

[13]  G. Grest,et al.  Granular flow down an inclined plane: Bagnold scaling and rheology. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  J. Rappaz,et al.  Numerical simulation of three‐dimensional free surface flows , 2003 .

[15]  Wolfgang Losert,et al.  Superstable granular heap in a thin channel. , 2003, Physical review letters.

[16]  H. Huppert,et al.  Axisymmetric collapses of granular columns , 2004, Journal of Fluid Mechanics.

[17]  G. Midi,et al.  On dense granular flows , 2003, The European physical journal. E, Soft matter.

[18]  Srinivasan Natesan,et al.  Arbitrary Lagrangian–Eulerian method for Navier–Stokes equations with moving boundaries , 2004 .

[19]  Jean-Pierre Vilotte,et al.  Spreading of a granular mass on a horizontal plane , 2004 .

[20]  E. J. Hinch,et al.  Study of the collapse of granular columns using two-dimensional discrete-grain simulation , 2005, Journal of Fluid Mechanics.

[21]  Crucial role of sidewalls in granular surface flows: consequences for the rheology , 2005, Journal of Fluid Mechanics.

[22]  R. Kerswell,et al.  Dam break with Coulomb friction: a model for granular slumping , 2005 .

[23]  Jean-Pierre Vilotte,et al.  On the use of Saint Venant equations to simulate the spreading of a granular mass , 2005 .

[24]  Roberto Zenit,et al.  Computer simulations of the collapse of a granular column , 2005 .

[25]  O. Cazacu,et al.  Compressible rigid viscoplastic fluids , 2006 .

[26]  E. Hinch,et al.  Raining into shallow water as a description of the collapse of a column of grains , 2005, Journal of Fluid Mechanics.

[27]  Olivier Pouliquen,et al.  A constitutive law for dense granular flows , 2006, Nature.

[28]  Herbert E. Huppert,et al.  Static and flowing regions in granular collapses down channels: Insights from a sedimenting shallow water model , 2007 .

[29]  F. Bouchut,et al.  On new erosion models of Savage–Hutter type for avalanches , 2008 .

[30]  Rich R. Kerswell,et al.  Planar collapse of a granular column: Experiments and discrete element simulations , 2008 .

[31]  Giovanni B. Crosta,et al.  Numerical modeling of 2‐D granular step collapse on erodible and nonerodible surface , 2009 .

[32]  R. Kerswell,et al.  Axisymmetric granular collapse: a transient 3D flow test of viscoplasticity. , 2008, Physical review letters.

[33]  I. Ionescu Onset and dynamic shallow flow of a viscoplastic fluid on a plane slope , 2010 .

[34]  N. Mangold,et al.  Erosion and mobility in granular collapse over sloping beds , 2010 .

[35]  A. Tamburrino,et al.  On the role of the ambient fluid on gravitational granular flow dynamics , 2010, Journal of Fluid Mechanics.

[36]  Giovanni B. Crosta,et al.  Numerical modeling of landquakes , 2010 .

[37]  Yasamin Khazraii,et al.  Numerical modeling , 2019, Ground Penetrating Radar: Improving sensing and imaging through numerical modeling.

[38]  Pierre-Yves Lagrée,et al.  The granular column collapse as a continuum: validity of a two-dimensional Navier–Stokes model with a μ(I)-rheology , 2011, Journal of Fluid Mechanics.

[39]  Anthony Wachs,et al.  A three-dimensional discrete-grain model for the simulation of dam-break rectangular collapses: comparison between numerical results and experiments , 2012 .

[40]  Eleonore Stutzmann,et al.  Numerical modeling of the Mount Steller landslide flow history and of the generated long period seismic waves , 2012 .

[41]  Frédéric Hecht,et al.  New development in freefem++ , 2012, J. Num. Math..

[42]  I. Ionescu Viscoplastic shallow flow equations with topography , 2013 .

[43]  A. Mangeney,et al.  Fundamental changes of granular flow dynamics, deposition, and erosion processes at high slope angles: Insights from laboratory experiments , 2014 .

[44]  Ioan R. Ionescu,et al.  Augmented Lagrangian for shallow viscoplastic flow with topography , 2013, J. Comput. Phys..

[45]  K. Allstadt,et al.  Numerical modeling of the Mount Meager landslide constrained by its force history derived from seismic data , 2015 .

[46]  D. G. Schaeffer,et al.  Well-posed and ill-posed behaviour of the ${\it\mu}(I)$ -rheology for granular flow , 2015, Journal of Fluid Mechanics.

[47]  J. Monnier,et al.  Inverse rheometry and basal properties inference for pseudoplastic geophysical flows , 2015 .

[48]  F. Bouchut,et al.  Viscoplastic modeling of granular column collapse with pressure-dependent rheology , 2015 .

[49]  E. Calder,et al.  Friction weakening in granular flows deduced from seismic records at the Soufrière Hills Volcano, Montserrat , 2015 .

[50]  M. Yamada,et al.  Estimation of dynamic friction of the Akatani landslide from seismic waveform inversion and numerical simulation , 2016 .

[51]  F. Bouchut,et al.  An analytic approach for the evolution of the static/flowing interface in viscoplastic granular flows , 2016 .

[52]  P. Richard,et al.  Granular and particle-laden flows: from laboratory experiments to field observations , 2017 .