Modified kinetic theory applied to the shear flows of granular materials

Abstract

Granular materials are characterized by large collections of discrete particles, where the particle-particle interactions are significantly more important than the particle-fluid interactions. The current kinetic theory captures fairly accurately the granular flow behavior in the dilute case, when only binary interactions are significant, but is not accurate at all in the dense flow regime, where multi-particle interactions and contacts must be modeled. To improve the kinetic theory results for granular flows in the dense flow regime, we propose a Modified Kinetic Theory (MKT) model that utilizes the contact duration or cutoff time to account for the complex particle-particle interactions in the dense regime. The contact duration model, also called TC model, was originally proposed by Luding and McNamara [“How to handle the inelastic collapse of a dissipative hard-sphere gas with the TC model,” Granular Matter 1, 113 (1998)] to solve the inelastic collapse issue existing in the inelastic hard sphere model. This model defines a cutoff time tc such that dissipation is not counted if the time between two consecutive contacts is less than tc. As shown in their study, the use of a cutoff time tc can also reduce the dissipation during multi-particle contacts. In this paper we relate the TC model with the Discrete Element Method (DEM) by choosing the cutoff time tc to be the duration of contact calculated from the linear-spring-dashpot soft-sphere model of the DEM. We examine two types of granular flows, simple shear flow and the plane shear flow, and compare the results of the classical kinetic theory model, the present MKT model, and the DEM model. We show that the MKT model entails a significant improvement over the kinetic theory model for simple shear flows at inertial regimes. With the MKT model the calculations are close to the DEM results at solid fractions as high as 0.57. Even for the plane shear flows, where shear rate and solid fraction are inhomogeneous, the results of the MKT model agree very well with the DEM results.