Towards dense, realistic granular media in 2D

The development of an applicable theory for granular matter - with both qualitative and quantitative value - is a challenging prospect, given the multitude of states, phases and (industrial) situations it has to cover. Given the general balance equations for mass, momentum and energy, the limiting case of dilute and almost elastic granular gases, where kinetic theory works perfectly well, is the starting point.\ud \ud In most systems, low density co-exists with very high density, where the latter is an open problem for kinetic theory. Furthermore, many additional nonlinear phenomena and material properties are important in realistic granular media, involving, e.g.:\ud \ud (i) multi-particle interactions and elasticity\ud (ii) strong dissipation,\ud (iii) friction,\ud (iv) long-range forces and wet contacts,\ud (v) wide particle size distributions and\ud (vi) various particle shapes.\ud \ud \ud Note that, while some of these issues are more relevant for high density, others are important for both low and high densities; some of them can be dealt with by means of kinetic theory, some cannot.\ud \ud This paper is a review of recent progress towards more realistic models for dense granular media in 2D, even though most of the observations, conclusions and corrections given are qualitatively true also in 3D.\ud \ud Starting from an elastic, frictionless and monodisperse hard sphere gas, the (continuum) balance equations of mass, momentum and energy are given. The equation of state, the (Navier–Stokes level) transport coefficients and the energy-density dissipation rate are considered. Several corrections are applied to those constitutive material laws - one by one - in order to account for the realistic physical effects and properties listed above

Clustering Instabilities, Arching, and Anomalous Interaction Probabilities as Examples for Cooperative Phenomena in Dry Granular Media

In a freely cooling granular material fluctuations in density and temperature cause position dependent energy loss. Due to strong local dissipation, pressure and energy drop rapidly and material moves from `hot' to `cold' regions, leading to even stronger dissipation and thus causing the density instability. The assumption of `molecular chaos' is valid only in the homogeneous cooling regime. As soon as the density instability occurs, the impact parameter is not longer uniformly distributed. The pair-correlation and the structure functions show that the molecular chaos assumption --- together with reasonable excluded volume modeling --- is important for short distances and irrelevant on large length scales. In this study, the probability distribution of the collision frequency is examined for pipe flow and for freely cooling granular materials as well. Uncorrelated events lead to a Poisson distribution for the collision frequencies. In contrast, the fingerprint of the cooperative phenomena discussed here is a power-law decay of the probability for many collisions per unit time. Keywords: discrete element method, event driven simulations, clustering instability, arching, shock waves, power-law distribution, cooperative phenomena.Comment: 27 pages 14 figs (2 color

How to handle the inelastic collapse of a dissipative hard-sphere gas with the TC model

The inelastic hard sphere model of granular material is simple, easily accessible to theory and simulation, and captures much of the physics of granular media. It has three drawbacks, all related to the approximation that collisions are instantaneous: 1) The number of collisions per unit time can diverge, i.e. the ``inelastic collapse'' can occur. 2) All interactions are binary, multiparticle contacts cannot occur and 3) no static limit exists. We extend the inelastic hard sphere model by defining a duration of contact t_c such that dissipation is allowed only if the time between contacts is larger than t_c. We name this generalized model the ``TC model'' and discuss it using examples of dynamic and static systems. The contact duration used here does not change the instantaneous nature of the hard sphere contacts, but accounts for a reduced dissipation during ``multiparticle contacts''. Kinetic and elastic energies are defined as well as forces and stresses in the system. Finally, we present event-driven numerical simulations of situations far beyond the inelastic collapse, possible only with the TC model.Comment: 15 pages, Latex, 14 bw.ps figures + 2 col.ps figures, to be published in Granular Matter 1(3) 199

Cohesive, frictional powders: contact models for tension

The contacts between cohesive, frictional particles with sizes in the range 0.1–10 μm are the subject of this study. Discrete element model (DEM) simulations rely on realistic contact force models—however, too much details make both implementation and interpretation prohibitively difficult. A rather simple, objective contact model is presented, involving the physical properties of elastic–plastic repulsion, dissipation, adhesion, friction as well as rolling- and torsion-resistance. This contact model allows to model bulk properties like friction, cohesion and yield-surfaces. Very loose packings and even fractal agglomerates have been reported in earlier work. The same model also allows for pressure-sintering and tensile strength tests as presented in this study

Mean Field theory for a driven granular gas of frictional particles

We propose a mean field (MF) theory for a homogeneously driven granular gas of inelastic particles with Coulomb friction. The model contains three parameters, a normal restitution coefficient $r_n$, a maximum tangential restitution coefficient $r_t^m$, and a Coulomb friction coefficient $\mu$. The parameters can be tuned to explore a wide range of physical situations. In particular, the model contains the frequently used $\mu \to \infty$ limit as a special case. The MF theory is compared with the numerical simulations of a randomly driven monolayer of spheres for a wide range of parameter values. If the system is far away from the clustering instability ($r_n \approx 1$), we obtain a good agreement between mean field and simulations for $\mu=0.5$ and $r_t^m=0.4$, but for much smaller values of $r_n$ the agreement is less good. We discuss the reasons of this discrepancy and possible refinements of our computational scheme.Comment: 6 pages, 3 figures (10 *.eps files), elsart style (macro included), in Proceedings of the International Conference "Statistical Mechanics and Strongly Correlated Systems", University of Rome "La Sapienza" (Italy), 27-29 September 199

Macroscopic model with anisotropy based on micro-macro informations

Physical experiments can characterize the elastic response of granular materials in terms of macroscopic state-variables, namely volume (packing) fraction and stress, while the microstructure is not accessible and thus neglected. Here, by means of numerical simulations, we analyze dense, frictionless, granular assemblies with the final goal to relate the elastic moduli to the fabric state, i.e., to micro-structural averaged contact network features as contact number density and anisotropy. The particle samples are first isotropically compressed and later quasi-statically sheared under constant volume (undrained conditions). From various static, relaxed configurations at different shear strains, now infinitesimal strain steps are applied to "measure" the effective elastic response; we quantify the strain needed so that plasticity in the sample develops as soon as contact and structure rearrangements happen. Because of the anisotropy induced by shear, volumetric and deviatoric stresses and strains are cross-coupled via a single anisotropy modulus, which is proportional to the product of deviatoric fabric and bulk modulus (i.e. the isotropic fabric). Interestingly, the shear modulus of the material depends also on the actual stress state, along with the contact configuration anisotropy. Finally, a constitutive model based on incremental evolution equations for stress and fabric is introduced. By using the previously measured dependence of the stiffness tensor (elastic moduli) on the microstructure, the theory is able to predict with good agreement the evolution of pressure, shear stress and deviatoric fabric (anisotropy) for an independent undrained cyclic shear test, including the response to reversal of strain

Constitutive relations for the shear band evolution in granular matter under large strain

A so-called “split-bottom ring shear cell” leads to wide shear bands under slow, quasi-static deformation. Unlike normal cylindrical Couette shear cells or rheometers, the bottom plate is split such that the outer part of it can move with the outer wall, while the other part (inner disk) is immobile. From discrete element simulations (DEM), several continuum fields like the density, velocity, deformation gradient and stress are computed and evaluated with the goal to formulate objective constitutive relations for the powder flow behavior. From a single simulation, by applying time- and (local) space-averaging, a non-linear yield surface is obtained with peculiar stress dependence.\ud \ud The anisotropy is always smaller than the macroscopic friction coefficient. However, the lower bound of anisotropy increases with the strain rate, approaching the maximum according to a stretched exponential with a specific rate that is consistent with a shear path of about one particle diameter