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

Consequences of using different pair-correlation functions on the stability properties of the Homogeneous Cooling State for a monodisperse system of near-elastic disks

We show the differences in the stability properties of the Homogeneous Cooling State (HCS) of a two-dimensional monodisperse collection of rigid and near-elastic disks, obtained by using different formulae for the pair-correlation function.For an equation of state that takes into account the crystallization and ordering of the particles (and the respective pressure drop), the critical wavelength of the heat conduction mode is considerably modified in the transition zone, involving a bifurcation and an additional mode of instability. The theoretical predictions, using the improved equation of state are confirmed by numerical simulations. Nevertheless, some open questions remain

Evolution of swelling pressure of cohesive-frictional, rough and elasto-plastic granulates

The subject of this study is the modeling of the evolution of the swell-ing pressure of granulates with cohesive-frictional, rough and elasto-plastic “mi-croscopic” contact properties. The spherical particles are randomly arranged in a periodic cubic space with a fixed volume so that an increase of the particle size – i.e. swelling that can be caused by intake of some fluid – is accompanied by a de-crease of the void space. An analytical function is proposed that properly de-scribes the (macroscopic) void ratio as function of pressure for different micro-scopic contact properties

The critical-state yield stress (termination locus) of adhesive powders from a single numerical experiment

Dry granular materials in a split-bottom ring shear cell geometry show wide shear bands under slow, quasi-static, large deformation. This system is studied in the presence of contact adhesion, using the discrete element method (DEM). Several continuum fields like the density, the deformation gradient and the stress tensor are computed locally and are analyzed with the goal to formulate objective constitutive relations for the flow behavior of cohesive powders. From a single simulation only, by applying time- and (local) space-averaging, and focusing on the regions of the system that experienced considerable deformations, the critical-state yield stress (termination locus) can be obtained. It is close to linear, for non-cohesive granular materials, and nonlinear with peculiar pressure dependence, for adhesive powders—due to the nonlinear dependence of the contact adhesion on the confining forces. The contact model is simplified and possibly will need refinements and additional effects in order to resemble realistic powders. However, the promising method of how to obtain a critical-state yield stress from a single numerical test of one material is generally applicable and waits for calibration and validation

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

A master equation for force distributions in soft particle packings - Irreversible mechanical responses to isotropic compression and decompression

Mechanical responses of soft particle packings to quasi-static deformations are determined by the microscopic restructuring of force-chain networks, where complex non-affine displacements of constituent particles cause the irreversible macroscopic behavior. Recently, we have proposed a master equation for the probability distribution functions of contact forces and interparticle gaps [K. Saitoh et al., Soft Matter 11, 1253 (2015)], where mutual exchanges of contacts and interparticle gaps, i.e. opening and closing contacts, are also involved in the stochastic description with the aid of Delaunay triangulations. We describe full details of the master equation and numerically investigate irreversible mechanical responses of soft particle packings to cyclic loading. The irreversibility observed in molecular dynamics simulations is well reproduced by the master equation if the system undergoes quasi-static deformations. We also confirm that the degree of irreversible responses is a decreasing function of the area fraction and the number of cycles.Comment: 17 pages, 21 figures (6 figures are not displayed