Friction dependence of shallow granular flows from discrete particle simulations

A shallow-layer model for granular flows is completed with a closure relation for the macroscopic bed friction or basal roughness obtained from micro-scale discrete particle simulations of steady flows. We systematically vary the bed friction by changing the contact friction coefficient between basal and flowing particles, while the base remains geometrically rough. By simulating steady uniform flow over a wide parameter range, we obtain a friction law that is a function of both flow and bed variables. Surprisingly, we find that the macroscopic bed friction is only weakly dependent on the contact friction of bed particles and predominantly determined by the properties of the flowing particles

Micro-macro transition and simplified contact models for wet granular materials

Wet granular materials in a quasi-static steady state shear flow have been studied with discrete particle simulations. Macroscopic quantities, consistent with the conservation laws of continuum theory, are obtained by time averaging and spatial coarse-graining. Initial studies involve understanding the effect of liquid content and liquid properties like the surface tension on the macroscopic quantities. Two parameters of the liquid bridge contact model have been studied as the constitutive parameters that define the structure of this model (i) the rupture distance of the liquid bridge model, which is proportional to the liquid content, and (ii) the maximum adhesive force, as controlled by the surface tension of the liquid. Subsequently a correlation is developed between these micro parameters and the steady state cohesion in the limit of zero confining pressure. Furthermore, as second result, the macroscopic torque measured at the walls, which is an experimentally accessible parameter, is predicted from our simulation results as a dependence on the micro-parameters. Finally, the steady state cohesion of a realistic non-linear liquid bridge contact model scales well with the steady state cohesion for a simpler linearized irreversible contact model with the same maximum adhesive force and equal energy dissipated per contact

A study of the anisotropy of stress in a fluid confined in a nanochannel

We present molecular dynamics simulations of planar Poiseuille flow of a Lennard-Jones fluid at various temperatures and body forces. Local thermostatting is used close to the walls to reach steady-state up to a limit body force. Macroscopic fields are obtained from microscopic data by time- and space-averaging and smoothing the data with a self-consistent coarse-graining method based on kernel interpolation. Two phenomena make the system interesting: (i) strongly confined fluids show layering, i.e., strong oscillations in density near the walls, and (ii) the stress deviates from the Newtonian fluid assumption, not only in the layered regime, but also much further away from the walls. Various scalar, vectorial, and tensorial fields are analyzed and related to each other in order to understand better the effects of both the inhomogeneous density and the anisotropy on the flow behavior and rheology. The eigenvalues and eigendirections of the stress tensor are used to quantify the anisotropy in stress and form the basis of a newly proposed objective, inherently anisotropic constitutive model that allows for non-collinear stress and strain gradient by construction

Macroscopic bulk cohesion and torque for wet granular materials

Wet granular materials in steady-state in a quasi-static flow have been studied with discrete particle simulations. The total torque is an experimentally accessible macroscopic quantity that can be used to investigate the shear strength, bulk cohesion and other properties of the materials. We report in this paper how the macroscopic bulk cohesion and torque required to rotate the system change with the liquid content. Consequently, micro-macro correlations are obtained for the macro properties as a function of the microscopic liquid bridge volume which is one factor dominating the contact force

Surface flow profiles for dry and wet granular materials by Particle Tracking Velocimetry; the effect of wall roughness

Two-dimensional Particle Tracking Velocimetry (PTV) is a promising technique to study the behaviour of granular flows. The aim is to experimentally determine the free surface width and position of the shear band from the velocity profile to validate simulations in a split-bottom shear cell geometry. The position and velocities of scattered tracer particles are tracked as they move with the bulk flow by analyzing images. We then use a new technique to extract the continuum velocity field, applying coarse-graining with the postprocessing toolbox MercuryCG on the discrete experimental PTV data. For intermediate filling heights, the dependence of the shear (or angular) velocity on the radial coordinate at the free surface is well fitted by an error function. From the error function, we get the width and the centre position of the shear band. We investigate the dependence of these shear band properties on filling height and rotation frequencies of the shear cell for dry glass beads for rough and smooth wall surfaces. For rough surfaces, the data agrees with the existing experimental results and theoretical scaling predictions. For smooth surfaces, particle-wall slippage is significant and the data deviates from the predictions. We further study the effect of cohesion on the shear band properties by using small amount of silicon oil and glycerol as interstitial liquids with the glass beads. While silicon oil does not lead to big changes, glycerol changes the shear band properties considerably. The shear band gets wider and is situated further inward with increasing liquid saturation, due to the correspondingly increasing trend of particles to stick together

Segregation of large particles in dense granular flows: A granular Saffman effect?

We report on the scaling between the lift force and the velocity lag experienced by a single particle of different size in a monodisperse dense granular chute flow. The similarity of this scaling to the Saffman lift force in (micro) fluids, suggests an inertial origin for the lift force responsible for segregation of (isolated, large) intruders in dense granular flows. We also observe an anisotropic pressure/stress field surrounding the particle, which potentially lies at the origin of the velocity lag. These findings are relevant for modelling and theoretical predictions of particle-size segregation. At the same time, the suggested interplay between polydispersity and inertial effects in dense granular flows with stress- and strain-gradients, implies striking new parallels between fluids, suspensions and granular flows with wide application perspectives

Compaction of granular material inside confined geometries

In both nature and engineering, loosely packed granular materials are often compacted inside confined geometries. Here, we explore such behaviour in a quasi-two dimensional geometry, where parallel rigid walls provide the confinement. We use the discrete element method to investigate the stress distribution developed within the granular packing as a result of compaction due to the displacement of a rigid piston. We observe that the stress within the packing increases exponentially with the length of accumulated grains, and show an extension to current analytic models which fits the measured stress. The micromechanical behaviour is studied for a range of system parameters, and the limitations of existing analytic models are described. In particular, we show the smallest sized systems which can be treated using existing models. Additionally, the effects of increasing piston rate, and variations of the initial packing fraction, are described

Mercury-DPM: Fast particle simulations in complex geometries

Mercury-DPM is a code for performing discrete particle simulations. That is to say, it simulates the motion of particles, or atoms, by applying forces and torques that stem either from external body forces, (e.g. gravity, magnetic fields, etc…) or from particle interactions. For granular particles, these are typically contact forces (elastic, viscous, frictional, plastic, cohesive), while for molecular simulations, forces typically stem from interaction potentials (e.g. Lennard-Jones). Often the method used in these packages is referred to as the discrete element method (DEM), which was originally designed for geotechnical applications. However, as Mercury-DPM is designed for simulating particles with emphasis on contact models, optimized contact detection for highly different particle sizes, and in-code coarse graining (in contrast to post-processing), we prefer the more general name discrete particle simulation. The code was originally developed for granular chute flows, and has since been extended to many other granular applications, including the geophysical modeling of cinder cone creation. Despite its granular heritage it is designed in a flexible way so it can be adapted to include other features such as long-range interactions and non-spherical particles, etc