Constitutive relations for compressible granular flow in the inertial regime

Schaeffer, D. G., Barker, T., Tsuji, D., Gremaud, P., Shearer, M. and Gray, J. M. N. T. 2019. Constitutive relations for compressible granular flow in the inertial regime. Journal of Fluid Mechanics 874 , pp. 926-951. 10.1017/jfm.2019.476

PDF - Published Version Available under License Creative Commons Attribution. Download (688kB)

Abstract

Granular flows occur in a wide range of situations of practical interest to industry, in our natural environment and in our everyday lives. This paper focuses on granular flow in the so-called inertial regime, when the rheology is independent of the very large particle stiffness. Such flows have been modelled with the μ(I),Φ(I)-rheology, which postulates that the bulk friction coefficient μ (i.e. the ratio of the shear stress to the pressure) and the solids volume fraction ϕ are functions of the inertial number I only. Although the μ(I),Φ(I)-rheology has been validated in steady state against both experiments and discrete particle simulations in several different geometries, it has recently been shown that this theory is mathematically ill-posed in time-dependent problems. As a direct result, computations using this rheology may blow up exponentially, with a growth rate that tends to infinity as the discretization length tends to zero, as explicitly demonstrated in this paper for the first time. Such catastrophic instability due to ill-posedness is a common issue when developing new mathematical models and implies that either some important physics is missing or the model has not been properly formulated. In this paper an alternative to the μ(I),Φ(I)-rheology that does not suffer from such defects is proposed. In the framework of compressible I-dependent rheology (CIDR), new constitutive laws for the inertial regime are introduced; these match the well-established μ(I) and Φ(I) relations in the steady-state limit and at the same time are well-posed for all deformations and all packing densities. Time-dependent numerical solutions of the resultant equations are performed to demonstrate that the new inertial CIDR model leads to numerical convergence towards physically realistic solutions that are supported by discrete element method simulations.

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:	17 May 2023 23:12
URI:	https://orca.cardiff.ac.uk/id/eprint/144014

Citation Data

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

Actions (repository staff only)

Edit Item