Bond anisotropy and cohesion of wet granular materials

Scaling behavior of the tensile strength of viscocohesive granular aggregates

We numerically analyze the tensile strength of a single wet agglomerate modeled as a viscocohesive aggregate impacting a flat surface by using the discrete-element simulations. The viscocohesive agglomerate composed of primary spherical particles with the inclusion of the interstitial liquid in the form of the capillary bridges characterized by the cohesive and viscous forces between particles is extracted from a cuboidal sample of granular materials by applying a spherical probe. The tensile strength is measured from the impact test of a wet agglomerate by systematically varying different values of the surface tension of the interstitial liquid, the liquid viscosity, and the impact speed. We show that the tensile stress increases immediately when the collision occurs between the agglomerate and the flat surface. The peak of the tensile stress obtained after the collision, then decreases smoothly with increasing the particle movement. The maximum tensile stress is defined to be the tensile strength of such agglomerate. It is remarkable that the normalized tensile strength of such agglomerate can be well described as a function of a dimensionless impact number that incorporates the capillary number and Stokes number (calculated from the surface tension and the viscosity of the liquid and the impact rate of the agglomerate), thus providing the confirmation for the unified representation of the liquid properties and the impact rate of wet granular media.

Evolution of granular materials under isochoric cyclic simple shearing

By means of 3D particle dynamics simulations, we analyze the microstructure of granular materials subjected to isochoric (constant volume) cyclic shearing, which drives the system towards a liquefaction state characterized by loops of jamming-unjamming transition with periodic loss of strength and irreversible accumulation of shear strain. We first show that the macroscopic response obtained by these simulations agrees well with the most salient features of the well-known cyclic behavior of granular materials both before and after liquefaction. Then we investigate the evolution of particle connectivity, force transmission, and anisotropies of contact and force networks. The onset of liquefaction is marked by partial collapse of the force-bearing network with rapid drop of the coordination number and nonrattler fraction of particles, and significant broadening of the contact force probability density function, which begins in the preliquefaction period. We find that the jamming transition in each cycle occurs for a critical value of the coordination number that can be interpreted as the percolation threshold of the contact network and appears to be independent of the initial mean stress, void ratio, and cyclic shear amplitude. We show that upon unjamming in each cycle an isotropic loss of contacts occurs and is followed by the development of high contact anisotropy and a large proportion of particles with only two or three contacts. The higher mobility of the particles also involves a lower degree of frustration of particle rotations and thus lower friction mobilization and tangential force anisotropy. These findings are relevant to both undrained cyclic deformations of saturated soils and rheology of dense non-Brownian suspensions where volume change is coupled with pore liquid drainage conditions.

The role of inter-particle friction on rheology and texture of wet granular flows

In order to get insight into the rheology and texture of rough unsaturated granular flows, we study the effects of the inter-particle friction coefficient on the macroscopic attributes and the texture variables of steady-state shearing flow of wet granular materials by relying on three-dimensional (3D) particle dynamics simulations. The macroscopic attributes are characterized by the macroscopic friction coefficient, macroscopic cohesion, and the packing fraction. The microstructural variables are characterized by the fabric and force anisotropies, the coordination number, and the stress transmission ratio. We show that the macroscopic observables behave as a function of the inertial number as a dry case for different values of the inter-particle friction coefficient. In particular, the macroscopic friction coefficient increases and the packing fraction decreases rapidly for small values of the friction coefficient, then they almost reach plateaus for higher values of the friction coefficient. Interestingly, all the macroscopic observables nicely behave as a function of the small values of the friction coefficient. Similarly, we also observe these characteristics for the fabric and force anisotropies and the coordination number as well as the stress transmission ratio which reflects the intermediate relationship between the microstructure and the mechanical behavior of such flows.

Additive rheology of complex granular flows

Granular flows are omnipresent in nature and industrial processes, but their rheological properties such as apparent friction and packing fraction are still elusive when inertial, cohesive and viscous interactions occur between particles in addition to frictional and elastic forces. Here we report on extensive particle dynamics simulations of such complex flows for a model granular system composed of perfectly rigid particles. We show that, when the apparent friction and packing fraction are normalized by their cohesion-dependent quasistatic values, they are governed by a single dimensionless number that, by virtue of stress additivity, accounts for all interactions. We also find that this dimensionless parameter, as a generalized inertial number, describes the texture variables such as the bond network connectivity and anisotropy. Encompassing various stress sources, this unified framework considerably simplifies and extends the modeling scope for granular dynamics, with potential applications to powder technology and natural flows.

Granular materials are abundant in nature, but we haven’t fully understood their rheological properties as complex interactions between particles are involved. Here, Vo et al. show that granular flows can be described by a generalized dimensionless number based on stress additivity.

Evolution of wet agglomerates inside inertial shear flow of dry granular materials

We use particle dynamics simulations to investigate the evolution of a wet agglomerate inside homogeneous shear flows of dry particles. The agglomerate is modeled by introducing approximate analytical expressions of capillary and viscous forces between particles in addition to frictional contacts. During shear flow, the agglomerate may elongate, break, or be eroded by loss of its capillary bonds and primary particles. By systematically varying the shear rate and surface tension of the binding liquid, we characterize the rates of these dispersion modes. All the rates increase with increasing inertial number of the flow and decreasing cohesion index of the agglomerate. We show that the data points for each mode collapse on a master curve for a dimensionless scaling parameter that combines the inertial number and the cohesion index. The erosion rate vanishes below a cutoff value of the scaling parameter. This leads to a power-law borderline between the vanishing erosion states and erosion states in the phase space defined by the inertial number and the cohesion index.

Agglomeration of wet particles in dense granular flows

In order to get insight into the wet agglomeration process, we numerically investigate the growth of a single granule inside a dense flow of an initially homogeneous distribution of wet and dry particles. The simulations are performed by means of the discrete element method and the binding liquid is assumed to be transported by the wet particles, which interact via capillary and viscous force laws. The granule size is found to be an exponential function of time, reflecting the conservation of the amount of liquid and the decrease of the number of available wet particles inside the flow during agglomeration. We analyze this behavior in terms of the accretion and erosion rates of wet particles for a range of different values of material parameters such as mean particle size, size polydispersity, friction coefficient and liquid viscosity. In particular, we propose a phase diagram of the granule growth as a function of the mean primary particle diameter and particle size span, which separates the parametric domain in which the granule grows from the domain in which the granule does not survive.

Root growth and force chains in a granular soil

Roots provide basic functions to plants such as water and nutrient uptake and anchoring in soil. The growth and development of root systems contribute to colonizing the surrounding soil and optimizing the access to resources. It is generally known that the variability of plant root architecture results from the combination of genetic, physiological, and environmental factors, in particular soil mechanical resistance. However, this last factor has never been investigated at the soil grain scale for roots. In this paper, we are interested in the effect of the disordered texture of granular soils on the evolution of forces experienced by the root cap during its growth. We introduce a numerical model in which the root is modeled as a flexible self-elongating tube that probes a soil composed of solid particles. By means of extensive simulations, we show that the forces exerted on the root cap reflect interparticle force chains. Our simulations also show that the mean force declines exponentially with root flexibility, the highest force corresponding to the soil hardness. Furthermore, we find that this functional dependence is characterized by a single dimensionless parameter that combines granular structure and root bending stiffness. This finding will be useful to further address the biological issues of mechanosensing and thigmomorphogenesis in plant roots.

Shear strength of wet granular materials: Macroscopic cohesion and effective stress : Discrete numerical simulations, confronted to experimental measurements

Rheometric measurements on assemblies of wet polystyrene beads, in steady uniform quasistatic shear flow, for varying liquid content within the small saturation (pendular) range of isolated liquid bridges, are supplemented with a systematic study by discrete numerical simulations. The numerical results agree quantitatively with the experimental ones provided that the intergranular friction coefficient is set to the value [Formula: see text], identified from the behaviour of the dry material. Shear resistance and solid fraction [Formula: see text] are recorded as functions of the reduced pressure [Formula: see text], which, defined as [Formula: see text], compares stress [Formula: see text], applied in the velocity gradient direction, to the tensile strength [Formula: see text] of the capillary bridges between grains of diameter a, and characterizes cohesion effects. The simplest Mohr-Coulomb relation with [Formula: see text]-independent cohesion c applies as a good approximation for large enough [Formula: see text] (typically [Formula: see text]. Numerical simulations extend to different values of μ and, compared to experiments, to a wider range of [Formula: see text]. The assumption that capillary stresses act similarly to externally applied ones onto the dry granular contact network (effective stresses) leads to very good (although not exact) predictions of the shear strength, throughout the numerically investigated range [Formula: see text] and [Formula: see text]. Thus, the internal friction coefficient [Formula: see text] of the dry material still relates the contact force contribution to stresses, [Formula: see text], while the capillary force contribution to stresses, [Formula: see text], defines a generalized Mohr-Coulomb cohesion c, depending on [Formula: see text] in general. c relates to [Formula: see text] , coordination numbers and capillary force network anisotropy. c increases with liquid content through the pendular regime interval, to a larger extent, the smaller the friction coefficient. The simple approximation ignoring capillary shear stress [Formula: see text] (referred to as the Rumpf formula) leads to correct approximations for the larger saturation range within the pendular regime, but fails to capture the decrease of cohesion for smaller liquid contents.

Rheology of granular materials composed of crushable particles

We investigate sheared granular materials composed of crushable particles by means of contact dynamics simulations and the bonded-cell model for particle breakage. Each particle is paved by irregular cells interacting via cohesive forces. In each simulation, the ratio of the internal cohesion of particles to the confining pressure, the relative cohesion, is kept constant and the packing is subjected to biaxial shearing. The particles can break into two or more fragments when the internal cohesive forces are overcome by the action of compressive force chains between particles. The particle size distribution evolves during shear as the particles continue to break. We find that the breakage process is highly inhomogeneous both in the fragment sizes and their locations inside the packing. In particular, a number of large particles never break whereas a large number of particles are fully shattered. As a result, the packing keeps the memory of its initial particle size distribution, whereas a power-law distribution is observed for particles of intermediate size due to consecutive fragmentation events whereby the memory of the initial state is lost. Due to growing polydispersity, dense shear bands are formed inside the packings and the usual dilatant behavior is reduced or cancelled. Hence, the stress-strain curve no longer passes through a peak stress, and a progressive monotonic evolution towards a pseudo-steady state is observed instead. We find that the crushing rate is controlled by the confining pressure. We also show that the shear strength of the packing is well expressed in terms of contact anisotropies and force anisotropies. The force anisotropy increases while the contact orientation anisotropy declines for increasing internal cohesion of the particles. These two effects compensate each other so that the shear strength is nearly independent of the internal cohesion of particles.

Viscoinertial regime of immersed granular flows

By means of extensive coupled molecular dynamics-lattice Boltzmann simulations, accounting for grain dynamics and subparticle resolution of the fluid phase, we analyze steady inertial granular flows sheared by a viscous fluid. We show that, for a broad range of system parameters (shear rate, confining stress, fluid viscosity, and relative fluid-grain density), the frictional strength and packing fraction can be described by a modified inertial number incorporating the fluid effect. In a dual viscous description, the effective viscosity diverges as the inverse square of the difference between the packing fraction and its jamming value, as observed in experiments. We also find that the fabric and force anisotropies extracted from the contact network are well described by the modified inertial number, thus providing clear evidence for the role of these key structural parameters in dense suspensions.

Shear strength and microstructure of polydisperse packings: The effect of size span and shape of particle size distribution

By means of extensive contact dynamics simulations, we analyzed the effect of particle size distribution (PSD) on the strength and microstructure of sheared granular materials composed of frictional disks. The PSDs are built by means of a normalized β function, which allows the systematic investigation of the effects of both, the size span (from almost monodisperse to highly polydisperse) and the shape of the PSD (from linear to pronouncedly curved). We show that the shear strength is independent of the size span, which substantiates previous results obtained for uniform distributions by packing fraction. Notably, the shear strength is also independent of the shape of the PSD, as shown previously for systems composed of frictionless disks. In contrast, the packing fraction increases with the size span, but decreases with more pronounced PSD curvature. At the microscale, we analyzed the connectivity and anisotropies of the contacts and forces networks. We show that the invariance of the shear strength with the PSD is due to a compensation mechanism which involves both geometrical sources of anisotropy. In particular, contact orientation anisotropy decreases with the size span and increases with PSD curvature, while the branch length anisotropy behaves inversely.

Unexpected Effects of Activator Molecules' Polarity on the Electroreological Activity of Titanium Dioxide Nanopowders

Titanium dioxide nanoparticles, obtained using the sol-gel method and modified with organic solvents, such as acetone, acetonitrile, benzene, diethyl ether, dimethyl sulfoxide, toluene, and chloroform, were used as the filler of polydimethylsiloxane-based electrorheological fluids. The effect of electric field strength on the shear stress and yield stress of electrorheological fluids was investigated, as well as the spectra of their dielectric relaxation in the frequency range from 25 to 10⁶ Hz. Modification of titanium dioxide by polar molecules was found to enhance the electrorheological effect, as compared with unmodified TiO₂, in accordance with the widely accepted concept of polar molecule dominated electrorheological effect (PM-ER). The most unexpected result of this study was an increase in the electrorheological effect during the application of nonpolar solvents with zero or near-zero dipole moments as the modifiers. It is suggested that nonpolar solvents, besides providing additional polarization effects at the filler particles interface, alter the internal pressure in the gaps between the particles. As a result, the filler particles are attracted to one another, leading to an increase in their aggregation and the formation of a network of bonds between the particles through liquid bridge contacts. Such changes in the electrorheological fluid structure result in a significant increase in the mechanical strength of the structures that arise when an electric field is applied, and an increase in the observed electrorheological effect in comparison with the unmodified titanium dioxide.

Stochastic disks that roll

We study a model of rolling particles subject to stochastic fluctuations, which may be relevant in systems of nano- or microscale particles where rolling is an approximation for strong static friction. We consider the simplest possible nontrivial system: a linear polymer of three disks constrained to remain in contact and immersed in an equilibrium heat bath so the internal angle of the polymer changes due to stochastic fluctuations. We compare two cases: one where the disks can slide relative to each other and the other where they are constrained to roll, like gears. Starting from the Langevin equations with arbitrary linear velocity constraints, we use formal homogenization theory to derive the overdamped equations that describe the process in configuration space only. The resulting dynamics have the formal structure of a Brownian motion on a Riemannian or sub-Riemannian manifold, depending on if the velocity constraints are holonomic or nonholonomic. We use this to compute the trimer's equilibrium distribution with and without the rolling constraints. Surprisingly, the two distributions are different. We suggest two possible interpretations of this result: either (i) dry friction (or other dissipative, nonequilibrium forces) changes basic thermodynamic quantities like the free energy of a system, a statement that could be tested experimentally, or (ii) as a lesson in modeling rolling or friction more generally as a velocity constraint when stochastic fluctuations are present. In the latter case, we speculate there could be a "roughness" entropy whose inclusion as an effective force could compensate the constraint and preserve classical Boltzmann statistics. Regardless of the interpretation, our calculation shows the word "rolling" must be used with care when stochastic fluctuations are present.

Flow of wet granular materials: A numerical study

We simulate dense assemblies of frictional spherical grains in steady shear flow under controlled normal stress P in the presence of a small amount of an interstitial liquid, which gives rise to capillary menisci, assumed isolated (pendular regime), and attractive forces, which are hysteretic: Menisci form at contact, but do not break until grains are separated by a finite rupture distance. The system behavior depends on two dimensionless control parameters, inertial number I and reduced pressure P*=aP/(πΓ), comparing confining forces ∼a2P to meniscus tensile strength F0=πΓa, for grains of diameter a joined by menisci with surface tension Γ. We pay special attention to the quasistatic limit of slow flow and observe systematic, enduring strain localization in some of the cohesion-dominated (P*∼0.1) systems. Homogeneous steady flows are characterized by the dependence of internal friction coefficient μ* and solid fraction Φ on I and P*. We also record normal stress differences, fairly small but not negligible and increasing for decreasing P*. The system rheology is moderately sensitive to saturation within the pendular regime, but would be different in the absence of capillary hysteresis. Capillary forces have a significant effect on the macroscopic behavior of the system, up to P* values of several units, especially for longer force ranges associated with larger menisci. The concept of effective pressure may be used to predict an order of magnitude for the strong increase of μ* as P* decreases but such a crude approach is unable to account for the complex structural changes induced by capillary cohesion, with a significant decrease of Φ and different agglomeration states and anisotropic fabric. Likewise, the Mohr-Coulomb criterion for pressure-dependent critical states is, at best, an approximation valid within a restricted range of pressures, with P*≥1. At small enough P*, large clusters of interacting grains form in slow flows, in which liquid bonds survive shear strains of several units. This affects the anisotropies associated with different interactions and the shape of function μ*(I), which departs more slowly from its quasistatic limit than in cohesionless systems (possibly explaining the shear banding tendency).

Centrifugation-assisted Assembly of Colloidal Silica into Crack-Free and Transferrable Films with Tunable Crystalline Structures

Self-assembly of colloidal particles into colloidal films has many actual and potential applications. While various strategies have been developed to direct the assembly of colloidal particles, fabrication of crack-free and transferrable colloidal film with controllable crystal structures still remains a major challenge. Here we show a centrifugation-assisted assembly of colloidal silica spheres into free-standing colloidal film by using the liquid/liquid interfaces of three immiscible phases. Through independent control of centrifugal force and interparticle electrostatic repulsion, polycrystalline, single-crystalline and quasi-amorphous structures can be readily obtained. More importantly, by dehydration of silica particles during centrifugation, the spontaneous formation of capillary water bridges between particles enables the binding and pre-shrinkage of the assembled array at the fluid interface. Thus the assembled colloidal films are not only crack-free, but also robust and flexible enough to be easily transferred on various planar and curved substrates.

Initiation of immersed granular avalanches

By means of coupled molecular dynamics-computational fluid dynamics simulations, we analyze the initiation of avalanches in a granular bed of spherical particles immersed in a viscous fluid and inclined above its angle of repose. In quantitative agreement with experiments, we find that the bed is unstable for a packing fraction below 0.59 but is stabilized above this packing fraction by negative excess pore pressure induced by the effect of dilatancy. From detailed numerical data, we explore the time evolution of shear strain, packing fraction, excess pore pressures, and granular microstructure in this creeplike pressure redistribution regime, and we show that they scale excellently with a characteristic time extracted from a model based on the balance of granular stresses in the presence of a negative excess pressure and its interplay with dilatancy. The cumulative shear strain at failure is found to be ≃ 0.2, in close agreement with the experiments, irrespective of the initial packing fraction and inclination angle. Remarkably, the avalanche is triggered when dilatancy vanishes instantly as a result of fluctuations while the average dilatancy is still positive (expanding bed) with a packing fraction that declines with the initial packing fraction. Another nontrivial feature of this creeplike regime is that, in contrast to dry granular materials, the internal friction angle of the bed at failure is independent of dilatancy but depends on the inclination angle, leading therefore to a nonlinear dependence of the excess pore pressure on the inclination angle. We show that this behavior may be described in terms of the contact network anisotropy, which increases with a nearly constant connectivity and levels off at a value (critical state) that increases with the inclination angle. These features suggest that the behavior of immersed granular materials is controlled not only directly by hydrodynamic forces acting on the particles but also by the influence of the fluid on the granular microstructure.

Polydirectional stability of granular matter

We investigate jammed granular matter in a slowly rotating drum partially filled with granular material and find a state of polydirectional stability. In this state, the material responds elastically to small stresses in a wide angular interval while it responds by plastic deformation when subjected to small stresses outside this interval of directions. We describe the evolution of the granulate by means of a rate equation and find quantitative agreement with the experiment. The state of polydirectional stability complements the fragile state, where the material responds elastically to small applied stresses only in a certain direction but even very small stresses in any other direction would lead to plastic deformations. Similar to fragile matter, polydirectionally stable matter is created in a dynamic process by self-organization.

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.

Stability, deformation, and variability of granular fills composed of polyhedral particles

By means of extensive contact dynamics simulations, we investigate the mechanical equilibrium and deformation of a granular material composed of irregular polyhedral particles confined between two horizontal frictional planes. We show that, as a consequence of mobilized wall-particle friction forces at the top and bottom boundaries, the transient deformation induced by a constant vertical load is controlled by the aspect ratio (thickness over width) of the packing as well as the stress ratio. The transient deformation declines considerably for increasingly smaller aspect ratios and grows with the stress ratio. From the simulation data for a large number of independent configurations, we find that sample-to-sample fluctuations of the deformation have a broad distribution and they scale with the average deformation. We also analyze the evolution of particle connectivity during settlement and with the applied force. The face-face and edge-face contacts between polyhedral particles concentrate strong force chains with a growing proportion as a function of the applied force.

Application of colloid probe atomic force microscopy to the adhesion of thin films of viscous and viscoelastic silicone fluids

The adhesive characteristics of thin films (0.2-2 μm) of linear poly(dimethylsiloxane) (PDMS) liquids with a wide range of molecular weights have been measured using an atomic force microscope with a colloid probe (diameters 5 and 12 μm) for different separation velocities. The data were consistent with a residual film in the contact region having a thickness of ∼6 nm following an extended dwell time before separation of the probe. It was possible to estimate the maximum adhesive force as a function of the capillary number, Ca, by applying existing theoretical models based on capillary interactions and viscous flow except at large values of Ca in the case of viscoelastic fluids, for which it was necessary to develop a nonlinear viscoelastic model. The compliance of the atomic force microscope colloid beam was an important factor in governing the retraction velocity of the probe and therefore the value of the adhesive force, but the inertia of the beam and viscoelastic stress overshoot effects were not significant in the range of separation velocities investigated.