Skip to main content Accessibility help

Well-posed and ill-posed behaviour of the ${\it\mu}(I)$ -rheology for granular flow

  • T. Barker (a1), D. G. Schaeffer (a2), P. Bohorquez (a3) and J. M. N. T. Gray (a1)


In light of the successes of the Navier–Stokes equations in the study of fluid flows, similar continuum treatment of granular materials is a long-standing ambition. This is due to their wide-ranging applications in the pharmaceutical and engineering industries as well as to geophysical phenomena such as avalanches and landslides. Historically this has been attempted through modification of the dissipation terms in the momentum balance equations, effectively introducing pressure and strain-rate dependence into the viscosity. Originally, a popular model for this granular viscosity, the Coulomb rheology, proposed rate-independent plastic behaviour scaled by a constant friction coefficient ${\it\mu}$ . Unfortunately, the resultant equations are always ill-posed. Mathematically ill-posed problems suffer from unbounded growth of short-wavelength perturbations, which necessarily leads to grid-dependent numerical results that do not converge as the spatial resolution is enhanced. This is unrealistic as all physical systems are subject to noise and do not blow up catastrophically. It is therefore vital to seek well-posed equations to make realistic predictions. The recent ${\it\mu}(I)$ -rheology is a major step forward, which allows granular flows in chutes and shear cells to be predicted. This is achieved by introducing a dependence on the non-dimensional inertial number $I$ in the friction coefficient ${\it\mu}$ . In this paper it is shown that the ${\it\mu}(I)$ -rheology is well-posed for intermediate values of $I$ , but that it is ill-posed for both high and low inertial numbers. This result is not obvious from casual inspection of the equations, and suggests that additional physics, such as enduring force chains and binary collisions, becomes important in these limits. The theoretical results are validated numerically using two implicit schemes for non-Newtonian flows. In particular, it is shown explicitly that at a given resolution a standard numerical scheme used to compute steady-uniform Bagnold flow is stable in the well-posed region of parameter space, but is unstable to small perturbations, which grow exponentially quickly, in the ill-posed domain.

  • View HTML
    • Send article to Kindle

      To send this article to your Kindle, first ensure is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

      Note you can select to send to either the or variations. ‘’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

      Find out more about the Kindle Personal Document Service.

      Well-posed and ill-posed behaviour of the ${\it\mu}(I)$ -rheology for granular flow
      Available formats

      Send article to Dropbox

      To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

      Well-posed and ill-posed behaviour of the ${\it\mu}(I)$ -rheology for granular flow
      Available formats

      Send article to Google Drive

      To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

      Well-posed and ill-posed behaviour of the ${\it\mu}(I)$ -rheology for granular flow
      Available formats


This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (, which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.

Corresponding author

Email address for correspondence:


Hide All
Ancey, C., Coussot, P. & Evesque, P. 1999 A theoretical framework for granular suspensions in a steady simple shear flow. J. Rheol. 43, 16731699.
Andreotti, B., Forterre, Y. & Pouliquen, O. 2013 Granular Media. Between Fluid and Solid. Cambridge University Press.
Bocquet, L., Colin, A. & Ajdari, A. 2009 Kinetic theory of plastic flow in soft glassy materials. Phys. Rev. Lett. 103, 036001.
Bohorquez, P., Sanmiguel-Rojas, E., Sevilla, A., Jimenez-Gonzalez, J. & Martínez-Bazán, C. 2011 Stability and dynamics of the laminar wake past a slender blunt-based axisymmetric body. J. Fluid Mech. 676, 110144.
Bouzid, M., Trulsson, M., Claudin, P., Clément, E. & Andreotti, B. 2013 Non-local rheology of granular flows across yield conditions. Phys. Rev. Lett. 111 (23), 238301.
Boyer, F., Guazzelli, E. & Pouliquen, O. 2011a Unifying suspension and granular rheology. Phys. Rev. Lett. 107, 188301.
Boyer, F., Pouliquen, O. & Guazzelli, E. 2011b Dense suspensions in rotating-rod flows: normal stresses and particle migration. J. Fluid Mech. 686, 525.
Campbell, C. S. 2002 Granular shear flows at the elastic limit. J. Fluid Mech. 465, 261291.
Cawthorn, C. J.2010 Several applications of a model for dense granular flows. PhD thesis, University of Cambridge.
Couturier, E., Boyer, F., Pouliquen, O. & Guazzelli, E. 2011 Suspensions in a tilted trough: second normal stress difference. J. Fluid Mech. 686, 2639.
da Cruz, F., Emam, S., Prochnow, M., Roux, J. & Chevoir, F. 2005 Rheophysics of dense granular materials: discrete simulation of plane shear flows. Phys. Rev. E 72, 021309.
Edwards, A. N. & Gray, J. M. N. T. 2015 Erosion–deposition waves in shallow granular free-surface flows. J. Fluid Mech. 762, 3567.
Ehlers, W. & Volk, W. 1998 On theoretical and numerical methods in the theory of porous media based on polar and non-polar elasto-plastic solid materials. Intl J. Solids Struct. 35, 45974617.
Favero, J. L., Secchi, A. R., Cardozo, N. S. M. & Jasak, H. 2010 Viscoelastic flow analysis using the software OpenFOAM and differential constitutive equations. J. Non-Newtonian Fluid Mech. 165 (23–24), 16251636.
Ferziger, J. H. & Perić, M. 2002 Computational Methods for Fluid Dynamics. Springer.
Forterre, Y. 2006 Kapiza waves as a test for three-dimensional granular flow rheology. J. Fluid Mech. 563, 123132.
Forterre, Y. & Pouliquen, O. 2003 Long-surface-wave instability dense granular flows. J. Fluid Mech. 486, 2150.
Fowler, A. C. 1997 Mathematical Models in the Applied Sciences. Cambridge University Press.
GDR-MiDi 2004 On dense granular flows. Eur. Phys. J. E 14 (4), 341365.
Goldhirsch, I. 2003 Rapid granular flows. Annu. Rev. Fluid Mech. 35 (1), 267293.
Gray, J. M. N. T. 1999 Loss of hyperbolicity and ill-posedness of the viscous-plastic sea ice rheology in uniaxial divergent flow. J. Phys. Oceanogr. 29 (11), 29202929.
Gray, J. M. N. T. & Edwards, A. N. 2014 A depth-averaged ${\it\mu}(I)$ -rheology for shallow granular free-surface flows. J. Fluid Mech. 755, 503534.
Henann, D. L. & Kamrin, K. 2013 A predictive, size-dependent continuum model for dense granular flows. Proc. Natl Acad. Sci. USA 110, 67306735.
Howell, D., Behringer, R. P. & Veje, C. 1999 Stress fluctuations in a 2D granular couette experiment: a continuous transition. Phys. Rev. Lett. 82 (26), 52415244.
Jasak, H.1996 Error analysis and estimation in the Finite Volume method with applications to fluid flows. PhD thesis, Imperial College, University of London.
Jenkins, J. T. 2006 Dense shearing flows of inelastic disks. Phys. Fluids 18, 103307.
Jenkins, J. T. & Savage, S. B. 1983 A theory for the rapid flow of identical, smooth, nearly elastic, spherical-particles. J. Fluid Mech. 130, 187202.
Jiang, Y. M. & Liu, M. 2003 Granular elasticity without the Coulomb condition. Phys. Rev. Lett. 91, 144301.
Jop, P., Forterre, Y. & Pouliquen, O. 2005 Crucial role of sidewalls in granular surface flows: consequences for the rheology. J. Fluid Mech. 541, 167192.
Jop, P., Forterre, Y. & Pouliquen, O. 2006 A constitutive relation for dense granular flows. Nature 44, 727730.
Joseph, D. D. & Saut, J. C. 1990 Short-wave instabilities and ill-posed initial-value problems. Theor. Comput. Fluid Dyn. 1, 191227.
Kamrin, K. 2010 Nonlinear elasto-plastic model for dense granular flow. Intl J. Plast. 26, 167188.
Kamrin, K. & Koval, G. 2012 Non-local constitutive relation for steady granular flow. Phys. Rev. Lett. 108 (17), 178301.
Lagrée, P.-Y., Staron, L. & Popinet, S. 2011 The granular column collapse as a continuum: validity of a two-dimensional Navier–Stokes model with a ${\it\mu}(I)$ -rheology. J. Fluid Mech. 686, 378408.
Leray, J. 1953 Hyperbolic Differential Equations. Princeton Institute for Advanced Study.
Pitman, E. B. & Schaeffer, D. G. 1987 Stability of time dependent compressible granular flow in two dimensions. Commun. Pure Appl. Maths 40 (4), 421447.
Pouliquen, O., Cassar, C., Jop, P., Forterre, Y. & Nicolas, N. 2006 Flow of dense granular material: towards simple constitutive laws. J. Stat. Mech. 2006, P07020.
Pouliquen, O. & Forterre, Y. 2002 Friction law for dense granular flows: application to the motion of a mass down a rough inclined plane. J. Fluid Mech. 453, 133151.
Pouliquen, O. & Forterre, Y. 2009 A non-local rheology for dense granular flows. Phil. Trans. R. Soc. Lond. A 367, 50915107.
Razis, D., Edwards, A. N., Gray, J. M. N. T. & van der Weele, K. 2014 Arrested coarsening of granular roll waves. Phys. Fluids 26, 123305.
Savage, S. & Sayed, M. 1984 Stresses developed by dry cohesionless granular materials sheared in an annular shear cell. J. Fluid Mech. 142, 391430.
Schaeffer, D. 1987 Instability in the evolution-equations describing incompressible granular flow. J. Differ. Equ. 66 (1), 1950.
Schaeffer, D. G. & Pitman, E. B. 1988 Ill-posedness in three-dimensional plastic flow. Commun. Pure Appl. Maths 41 (7), 879890.
Schmid, P. & Henningson, D. 2001 Stability and Transition in Shear Flows. Springer.
Staron, L., Lagrée, P.-Y. & Popinet, S. 2012 The granular silo as a continuum plastic flow: the hour-glass vs the clepsydra. Phys. Fluids 24, 103301.
Sun, J. & Sundaresan, S. 2011 A constitutive model with microstructure evolution for flow of rate-independent granular materials. J. Fluid Mech. 682, 590616.
Tejchman, J. & Gudehus, G. 2001 Shearing of a narrow granular layer with polar quantities. Intl J. Numer. Anal. Meth. Geomech. 25, 128.
Theofilis, V. 2003 Advances in global linear instability analysis of non-parallel and three-dimensional flows. Prog. Aerosp. Sci. 39, 249315.
Theofilis, V. 2011 Global linear instability. Annu. Rev. Fluid Mech. 43, 319352.
Toiya, M., Stambaugh, J. & Losert, W. 2004 Transient and oscillatory granular shear flow. Phys. Rev. E 93, 088001.
Vardoulakis, I. & Aifantis, E. C. 1991 A gradient flow theory of plasticity for granular materials. Acta Mechanica 87, 197217.
Vardoulakis, I., Goldscheider, M. & Gudehus, G. 1978 Formation of shear bands in sand bodies as a bifurcation problem. Intl J. Numer. Anal. Meth. Geomech. 2, 99128.
Weller, H. G., Tabor, G., Jasak, H. & Fureby, C. 1998 A tensorial approach to computational continuum mechanics using object-oriented techniques. Comput. Phys. 12 (6), 620631.
Woodhouse, M. J., Thornton, A. R., Johnson, C. G., Kokelaar, B. P. & Gray, J. M. N. T. 2012 Segregation-induced fingering instabilities in granular free-surface flows. J. Fluid Mech. 709, 543580.
Wu, W., Bauer, E. & Kolymbas, D. 1996 Hypoplastic constitutive model with critical state for granular materials. Mech. Mater. 23 (1), 4569.
MathJax is a JavaScript display engine for mathematics. For more information see

JFM classification

Type Description Title

Barker et al. supplementary movie
The temporal growth of the error between the base state and transient pressure for Bagnold flow down an inclined plane.

 Video (1.8 MB)
1.8 MB


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed