Well-posed and ill-posed behaviour of the μ(I) -rheology for granular flow

Barker, T., Schaeffer, D. G., Bohorquez, P. and Gray, J. M. N. T. 2015. Well-posed and ill-posed behaviour of the μ(I) -rheology for granular flow. Journal of Fluid Mechanics 779 , pp. 794-818. 10.1017/jfm.2015.412

PDF - Published Version Available under License Creative Commons Attribution. Download (1MB)

Abstract

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 μ . 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 μ(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 μ . In this paper it is shown that the μ(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.

Item Type:	Article
Date Type:	Publication
Status:	Published
Schools:	Mathematics
Additional Information:	This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Publisher:	Cambridge University Press
ISSN:	0022-1120
Date of First Compliant Deposit:	22 September 2021
Last Modified:	02 May 2023 17:41
URI:	https://orca.cardiff.ac.uk/id/eprint/144012

Citation Data

Cited 108 times in Scopus. View in Scopus. Powered By Scopus® Data

Actions (repository staff only)

Edit Item