Shear strength properties of wet granular materials

Effect of cohesion on shear banding in quasistatic granular materials

Cohesive powders have widely different bulk behavior due to their peculiar interactions. We use discrete element simulations to investigate the effect of contact cohesion on the steady state flow of dense powders in a slowly sheared split-bottom Couette cell, which imposes a wide stable shear band. The intensity of cohesive forces can be quantified by the granular Bond number (Bo), namely the ratio between maximum attractive force and average force due to external compression. We find that the shear banding phenomenon is almost independent of cohesion for Bond numbers Bo<1, however for Bo≥1 cohesive forces start to play an important role, as both width and center position of the band increase. Inside the shear band, the mean normal contact force is independent of cohesion and depends only on the confining stress. In contrast, when the behavior is analyzed focusing on the eigendirections of the local strain rate tensor, a dependence on cohesion shows up. Forces carried by contacts along the compressive and tensile directions are symmetric about the mean force (larger and smaller respectively), while the force along the third, neutral direction follows the mean force. This anisotropy of the force network increases with cohesion, just like the heterogeneity in all (compressive, tensile and neutral) directions.

Aggregation and separation in ternary particle/oil/water systems with fully wettable particles

We report that a variety of ternary particle/liquid/liquid mixtures heavily aggregate or separate completely if (1) the particles are fully or almost fully wetted by one fluid, and (2) if the wetting fluid volume fraction is comparable to the particle volume fraction. Aggregation and separation do not happen if the particles are partially wetted by both fluids, in which case Pickering emulsions appear at all compositions. Rheological and geometric criteria for aggregation are proposed and compared with a state diagram of a ternary system composed of oil, water, and hydrophilic glass particles. Analogies are drawn to wet granulation and spherical agglomeration, two particle processing operations in which wetting phenomena are important.

Synthesis and mechanical response of disordered colloidal micropillars

A method for synthesizing and uniaxially compressing free-standing colloidal micropillars is presented. The mechanical response of the micropillars is strongly dependent upon their initial defect population and water content.

Erosion dynamics of a wet granular medium

Liquid may give strong cohesion properties to a granular medium, and confer a solidlike behavior. We study the erosion of a fixed circular aggregate of wet granular matter subjected to a flow of dry grains inside a half-filled rotating drum. During the rotation, the dry grains flow around the fixed obstacle. We show that its diameter decreases linearly with time for low liquid content, as wet grains are pulled out of the aggregate. This erosion phenomenon is governed by the properties of the liquids. The erosion rate decreases exponentially with the surface tension while it depends on the viscosity to the power -1. We propose a model based on the force fluctuations arising inside the flow, explaining both dependencies: The capillary force acts as a threshold and the viscosity controls the erosion time scale. We also provide experiments using different flowing grains, confirming our model.

Cohesive granular materials composed of nonconvex particles

The macroscopic cohesion of granular materials made up of sticky particles depends on the particle shapes. We address this issue by performing contact dynamics simulations of 2D packings of nonconvex aggregates. We find that the macroscopic cohesion is strongly dependent on the strain and stress inhomogeneities developing inside the material. The largest cohesion is obtained for nearly homogeneous deformation at the beginning of unconfined axial compression and it evolves linearly with nonconvexity. Interestingly, the aggregates in a sheared packing tend to form more contacts with fewer neighboring aggregates as the degree of nonconvexity increases. We also find that shearing leads either to an isotropic distribution of tensile contacts or to the same privileged direction as that of compressive contacts.

Sculpting sandcastles grain by grain: self-assembled sand towers

We study the spontaneous formation of granular towers produced when dry sand is poured on a wet sand bed. When the liquid content of the bed exceeds a threshold value W*, the impacting grains have a nonzero probability to stick on the wet grains due to instantaneous liquid bridges created during the impact. The trapped grains become wet by the capillary ascension of water and the process continues, giving rise to stable narrow towers. The growth velocity is determined by the surface liquid content which decreases exponentially as the tower height augments. This self-assembly mechanism (only observed in the funicular and capillary regimes) could theoretically last while the capillary rise of water is possible; however, the structure collapses before reaching this limit. The collapse occurs when the weight of the tower surpasses the cohesive stress at its base. The cohesive stress increases as the liquid content of the bed is reduced. Consequently, the highest towers are found just above W*.

Entangled granular media

We study the geometrically induced cohesion of ensembles of granular "u particles" that mechanically entangle through particle interpenetration. We vary the length-to-width ratio l/w of the u particles and form them into freestanding vertical columns. In a laboratory experiment, we monitor the response of the columns to sinusoidal vibration (with peak acceleration Γ). Column collapse occurs in a characteristic time τ which follows the relation τ∝exp(Γ/Δ). Δ resembles an activation energy and is maximal at intermediate l/w. A simulation reveals that optimal strength results from competition between packing and entanglement.

Simulation of cemented granular materials. II. Micromechanical description and strength mobilization at the onset of macroscopic yielding

This is the second of two papers investigating the mechanical response of cemented granular materials by means of contact dynamics simulations. In this paper, a two-dimensional polydisperse sample with high void ratio is sheared in a load-controlled simple shear numerical device until the stress state of the sample reaches the yield stress. We first study the stress transmission properties of the granular material in terms of the fabric of different subsets of contacts characterized by the magnitude of their normal forces. This analysis highlights the existence of a peculiar force carrying structure in the cemented material, which is reminiscent of the bimodal stress transmission reported for cohesionless granular media. Then, the evolution of contact forces and torques is investigated trying to identify the micromechanical conditions that trigger macroscopic yielding. It is shown that global failure can be associated to the apparition of a group of particles whose contacts fulfill at least one of the local rupture conditions. In particular, these particles form a large region that percolates through the sample at the moment of failure, evidencing the relationship between macroscopic yielding and the emergence of large-scale correlations in the system.

Simulation of cemented granular materials. I. Macroscopic stress-strain response and strain localization

This is the first of two papers investigating the mechanical response of cemented granular materials by means of contact dynamics simulations. In this paper, a two-dimensional polydisperse sample with high-void ratio is constructed and then sheared in a simple shear numerical device at different confinement levels. We study the macroscopic response of the material in terms of mean and deviatoric stresses and strains. We show that the introduction of a local force scale, i.e., the tensile strength of the cemented bonds, causes the material to behave in a rigid-plastic fashion, so that a yield surface can be easily determined. This yield surface has a concave-down shape in the mean:deviatoric stress plane and it approaches a straight line, i.e., a Coulomb strength envelope, in the limit of a very dense granular material. Beyond yielding, the cemented structure gradually degrades until the material eventually behaves as a cohesionless granular material. Strain localization is also investigated, showing that the strains concentrate in a shear band whose thickness increases with the confining stress. The void ratio inside the shear band at the steady state is shown to be a material property that depends only on contact parameters.

Stress-strain behavior and geometrical properties of packings of elongated particles

We present a numerical analysis of the effect of particle elongation on the quasistatic behavior of sheared granular media by means of the contact dynamics method. The particle shapes are rounded-cap rectangles characterized by their elongation. The macroscopic and microstructural properties of several packings subjected to biaxial compression are analyzed as a function of particle elongation. We find that the shear strength is an increasing linear function of elongation. Performing an additive decomposition of the stress tensor based on a harmonic approximation of the angular dependence of branch vectors, contact normals, and forces, we show that the increasing mobilization of friction force and the associated anisotropy are key effects of particle elongation. These effects are correlated with partial nematic ordering of the particles which tend to be oriented perpendicular to the major principal stress direction and form side-to-side contacts. However, the force transmission is found to be mainly guided by cap-to-side contacts, which represent the largest fraction of contacts for the most elongated particles. Another interesting finding is that, in contrast to shear strength, the solid fraction first increases with particle elongation but declines as the particles become more elongated. It is also remarkable that the coordination number does not follow this trend so that the packings of more elongated particles are looser but more strongly connected.

Bond anisotropy and cohesion of wet granular materials

We analyse the Coulomb cohesion of wet granular materials and its relationship with the distribution of capillary bonds between particles. We show that, within a harmonic representation of the bond and force states, the shear strength can be expressed as a state equation in terms of internal anisotropy parameters. This formulation involves a dependence of the shear strength on loading direction and leads to the fragile behaviour of granular materials. The contact dynamics simulations of a wet material, in which a capillary force law is prescribed, are in excellent agreement with the predictions of this model. We find that the fragile character decreases as the local adhesion is increased or the mean stress is decreased.

Process modeling in the pharmaceutical industry using the discrete element method

The discrete element method (DEM) is widely used to model a range of processes across many industries. This paper reviews current DEM models for several common pharmaceutical processes including material transport and storage, blending, granulation, milling, compression, and film coating. The studies described in this review yielded interesting results that provided insight into the effects of various material properties and operating conditions on pharmaceutical processes. Additionally, some basic elements common to most DEM models are overviewed. A discussion of some common model extensions such as nonspherical particle shapes, noncontact forces, and interstitial fluids is also presented. While these more complex systems have been the focus of many recent studies, considerable work must still be completed to gain a better understanding of how they can affect the processing behavior of bulk solids.

Computer simulation of model cohesive powders: plastic consolidation, structural changes, and elasticity under isotropic loads

The quasistatic behavior of a simple two-dimensional model of a cohesive powder under isotropic loads is investigated by discrete element simulations. We ignore contact plasticity and focus on the effect of geometry and collective rearrangements on the material behavior. The loose packing states, as assembled and characterized in a previous numerical study [Gilabert, Roux, and Castellanos, Phys. Rev. E 75, 011303 (2007)], are observed, under growing confining pressure P , to undergo important structural changes, while solid fraction Phi irreversibly increases (typically, from 0.4-0.5 to 0.75-0.8). The system state goes through three stages, with different forms of the plastic consolidation curve, i.e., Phi as a function of the growing reduced pressure P;{*}=PaF_{0} , defined with adhesion force F0 and grain diameter a . In the low-confinement regime (I), the system undergoes negligible plastic compaction, and its structure is influenced by the assembling process. In regime II the material state is independent of initial conditions, and the void ratio varies linearly with lnP [i.e., Delta(1Phi)=lambdaDelta(lnP;{*}) ], as described in the engineering literature. Plasticity index lambda is reduced in the presence of a small rolling resistance (RR). In the last stage of compaction (III), Phi approaches an asymptotic, maximum solid fraction Phi_{max} , as a power law Phi_{max}-Phi proportional, variant(P;{*});{-alpha} , with alpha approximately 1 , and properties of cohesionless granular packs are gradually retrieved. Under consolidation, while the range xi of fractal density correlations decreases, force patterns reorganize from self-balanced clusters to force chains, with correlative evolutions of force distributions, and elastic moduli increase by a large amount. Plastic deformation events correspond to very small changes in the network topology, while the denser regions tend to move like rigid bodies. Elastic properties are dominated by the bending of thin junctions in loose systems. For growing RR those tend to form particle chains, the folding of which, rather than tensile ruptures, controls plastic compaction.

Morphological clues to wet granular pile stability

When a granular material such as sand is mixed with a certain amount of liquid, the surface tension of the latter bestows considerable stiffness to the material, which enables, for example, sand castles to be sculpted. The geometry of the liquid interface within the granular pile is of extraordinary complexity and strongly varies with the liquid content. Surprisingly, the mechanical properties of the pile are largely independent of the amount of liquid over a wide range. We resolve this puzzle with the help of X-ray microtomography, showing that the remarkable insensitivity of the mechanical properties to the liquid content is due to the particular organization of the liquid in the pile into open structures. For spherical grains, a simple geometric rule is established, which relates the macroscopic properties to the internal liquid morphologies. We present evidence that this concept is also valid for systems with non-spherical grains. Hence, our results provide new insight towards understanding the complex physics of a large variety of wet granular systems including land slides, as well as mixing and agglomeration problems.

Jamming transition and new percolation universality classes in particulate systems with attraction

We numerically study the jamming transition in particulate systems with attraction by investigating their mechanical response at zero temperature (T=0). We find three regimes of mechanical behavior separated by two critical transitions--connectivity and rigidity percolation. The transitions belong to different universality classes than their lattice counterparts, due to force balance constraints. We also find that these transitions are unchanged at low temperatures and resemble gelation transitions in experiments on colloidal and silica gels.

Strength and failure of cemented granular matter

Cemented granular materials (CGMs) consist of densely packed solid particles and a pore-filling solid matrix sticking to the particles. We use a sub-particle lattice discretization method to investigate the particle-scale origins of strength and failure properties of CGMs. We show that jamming of the particles leads to highly inhomogeneous stress fields. The stress probability density functions are increasingly wider for a decreasing matrix volume fraction, the stresses being more and more concentrated in the interparticle contact zones with an exponential distribution as in cohesionless granular media. Under uniaxial loading, pronounced asymmetry can occur between tension and compression both in strength and in the initial stiffness as a result of the presence of bare contacts (with no matrix interposed) between the particles. Damage growth is analyzed by considering the evolution of stiffness degradation and the number of broken bonds in the particle phase. A brutal degradation appears in tension as a consequence of brittle fracture in contrast to the more progressive nature of damage growth in compression. We also carry out a detailed parametric study in order to assess the combined influence of the matrix volume fraction and particle-matrix adherence. Three regimes of crack propagation can be distinguished corresponding to no particle damage, particle abrasion and particle fragmentation, respectively. We find that particle damage scales well with the relative toughness of the particle-matrix interface with respect to the particle toughness. This relative toughness is a function of both matrix volume fraction and particle-matrix adherence and it appears therefore to be the unique control parameter governing transition from soft to hard behavior.

Chaoticity of the wet granular gas

In this work we derive an analytic expression for the Kolmogorov-Sinai entropy of dilute wet granular matter, valid for any spatial dimension. The grains are modeled as hard spheres and the influence of the wetting liquid is described according to the capillary model, in which dissipation is due to the hysteretic cohesion force of capillary bridges. The Kolmogorov-Sinai entropy is expanded in a series with respect to density. We find a rapid increase of the leading term when liquid is added. This demonstrates the sensitivity of the granular dynamics to humidity, and shows that the liquid significantly increases the chaoticity of the granular gas.

Computer simulation of model cohesive powders: influence of assembling procedure and contact laws on low consolidation states

Molecular dynamics simulations are used to investigate the structure and mechanical properties of a simple two-dimensional model of a cohesive granular material. Intergranular forces involve elasticity, Coulomb friction, and a short-range attraction akin to the van der Waals force in powders. The effects of rolling resistance (RR) at intergranular contacts are also studied. The microstructure of the cohesive packing under low pressure is shown to depend sensitively on the assembling procedure which is applied to the initially isolated particles of a granular gas. While a direct compression produces a final equilibrated configuration with a similar density to that of cohesionless systems, the formation of large aggregates prior to the application of an external pressure results in much looser stable packings. A crucial state variable is the ratio P;{*}=PaF_{0} of applied pressure P , acting on grains of diameter a , to maximum tensile contact force F0 . At low P;{*} the force-carrying structure and force distribution are sensitive to the level of velocity fluctuations in the early stages of cluster aggregation. The coordination number of packings with RR approaches 2 in the limit of low initial velocities or large rolling friction. In general the force network is composed of hyperstatic clusters, typically comprising four to a few tens of grains, in which forces reach values of the order of F0 , joined by barely rigid arms, where contact forces are very small. Under growing P;{*} , it quickly rearranges into force chainlike patterns that are more familiar in dense systems. Density correlations are interpreted in terms of a fractal structure, up to a characteristic correlation length xi of the order of ten particle diameters for the studied solid fractions. The fractal dimension in systems with RR coincides, within measurement uncertainties, with the ballistic aggregation result, in spite of a possibly different connectivity, but is apparently higher without RR. Possible effects of micromechanical and assembling process parameters on mechanical strength of packings are evoked.

Stress transmission in wet granular materials

We analyze stress transmission in wet granular media in the pendular state by means of three-dimensional molecular-dynamics simulations. We show that the tensile action of capillary bonds induces a self-stressed particle network organized in two percolating "phases" of positive and negative particle pressures. Various statistical descriptors of the microstructure and bond force network are used to characterize this partition. Two basic properties emerge: 1) the highest particle pressure is located in the bulk of each phase; 2) the lowest pressure level occurs at the interface between the two phases, involving also the largest connectivity of the particles via tensile and compressive bonds. When a confining pressure is applied, the number of tensile bonds falls off and the negative phase breaks into aggregates and isolated sites.

Additive decomposition of shear strength in cohesive granular media from grain-scale interactions

We study cemented granular media by introducing cohesive bonding (sliding or rolling friction and tensile strength) between grains in the framework of the contact dynamics method. We find that, for a wide range of bond parameters, the macroscopic angle of friction at the peak state can be split into three distinct terms of collisional, frictional and dilational origins. Remarkably, the macroscopic tensile strength depends only on the bond tensile strength, and the friction angle at the peak state is proportional to the dilatancy angle which varies linearly with sliding friction.