Fluidization of wet granulates under shear

Sticky Matters: Jamming and Rigid Cluster Statistics with Attractive Particle Interactions

While the large majority of theoretical and numerical studies of the jamming transition consider athermal packings of purely repulsive spheres, real complex fluids and soft solids generically display attraction between particles. By studying the statistics of rigid clusters in simulations of soft particles with an attractive shell, we present evidence for two distinct jamming scenarios. Strongly attractive systems undergo a continuous transition in which rigid clusters grow and ultimately diverge in size at a critical packing fraction. Purely repulsive and weakly attractive systems jam via a first-order transition, with no growing cluster size. We further show that the weakly attractive scenario is a finite size effect, so that for any nonzero attraction strength, a sufficiently large system will fall in the strongly attractive universality class. We therefore expect attractive jamming to be generic in the laboratory and in nature.

Rheology of dilute cohesive granular gases

Rheology of a dilute cohesive granular gas is theoretically and numerically studied. The flow curve between the shear viscosity and the shear rate is derived from the inelastic Boltzmann equation for particles having square-well potentials in a simple shear flow. It is found that (i) the stable uniformly sheared state only exists above a critical shear rate and (ii) the viscosity in the uniformly sheared flow is almost identical to that for uniformly sheared flow of hard core granular particles. Below the critical shear rate, clusters grow with time, in which the viscosity can be approximated by that for the hard-core fluids if we replace the diameter of the particle by the mean diameter of clusters.

Athermal rheology of weakly attractive soft particles

We study the rheology of a soft particulate system where the interparticle interactions are weakly attractive. Using extensive molecular dynamics simulations, we scan across a wide range of packing fractions (ϕ), attraction strengths (u), and imposed shear rates (γ[over ̇]). In striking contrast to repulsive systems, we find that at small shear rates generically a fragile isostatic solid is formed even if we go to ϕ≪ϕ_{J}. Further, with increasing shear rates, even at these low ϕ, nonmonotonic flow curves occur which lead to the formation of persistent shear bands in large enough systems. By tuning the damping parameter, we also show that inertia plays an important role in this process. Furthermore, we observe enhanced particle dynamics in the attraction-dominated regime as well as a pronounced anisotropy of velocity and diffusion constant, which we take as precursors to the formation of shear bands. At low enough ϕ, we also observe structural changes via the interplay of low shear rates and attraction with the formation of microclusters and voids. Finally, we characterize the properties of the emergent shear bands, and thereby, we find surprisingly small mobility of these bands, leading to prohibitively long time scales and extensive history effects in ramping experiments.

Arrest stress of uniformly sheared wet granular matter

We conduct extensive independent numerical experiments considering frictionless disks without internal degrees of freedom (rotation, etc.) in two dimensions. We report here that for a large range of the packing fractions below random-close packing, all components of the stress tensor of wet granular materials remain finite in the limit of zero shear rate. This is direct evidence for a fluid-to-solid arrest transition. The offset value of the shear stress characterizes plastic deformation of the arrested state which corresponds to dynamic yield stress of the system. Based on an analytical line of argument, we propose that the mean number of capillary bridges per particle, ν, follows a nontrivial dependence on the packing fraction, ϕ, and the capillary energy, ɛ. Most noticeably, we show that ν is a generic and universal quantity which does not depend on the driving protocol. Using this universal quantity, we calculate the arrest stress, σ(a), analytically based on a balance of the energy injection rate due to the external force driving the flow and the dissipation rate accounting for the rupture of capillary bridges. The resulting prediction of σ(a) is a nonlinear function of the packing fraction, ϕ, and the capillary energy, ɛ. This formula provides an excellent, parameter-free prediction of the numerical data. Corrections to the theory for small and large packing fractions are connected to the emergence of shear bands and of contributions to the stress from repulsive particle interactions, respectively.

Impact of attractive interactions on the rheology of dense athermal particles

Using numerical simulations, the rheological response of an athermal assembly of soft particles with tunable attractive interactions is studied in the vicinity of jamming. At small attractions, a fragile solid develops and a finite yield stress is measured. Moreover, the measured flow curves have unstable regimes, which lead to persistent shear banding. These features are rationalized by establishing a link between the rheology and the interparticle connectivity, which also provides a minimal model to describe the flow curves.

Universal features of the jamming phase diagram of wet granular materials

We investigate the influence of the shape of a particle on the structure of the jamming phase diagram of wet granular materials. We compute the jamming phase diagram of wet dimers (two fused disks) and compare it with that of the wet disks. Amplitude of the external force at solidification, i.e., the jamming force F(s), is computed as a function of the packing fraction ϕ, the capillary bridge energy ɛ, and the aspect ratio of dimers α. Based on data collapse, an equation for amplitude of the external force at solidification F(s)(ϕ,ɛ,α) is derived. F(s) has scaling and logarithmic relations with ϕ and ɛ, respectively, exactly the same type reported for wet disks earlier. Interestingly, F(s) does not depend on the aspect ratio of dimers α. The only difference is that wet dimers are found to be more stiffer than wet disks. However, the similarities of the equations describing F(s)(ϕ,ɛ,α) of wet dimers and disks imply that there exists, yet unknown, universal aspects of mechanical response of wet granular materials to the external forces, independent from the particle shape. In addition, we study local orientation of particles and its statistical properties.

From concave to convex: capillary bridges in slit pore geometry

We investigate the morphological evolution of nonaxisymmetric capillary bridges in slit-pore geometry as the height of the pore and aspect ratio of the bridge are varied. The liquid bridges are formed between two hydrophobic surfaces patterned with hydrophilic strips. The aspect ratio of the capillary bridges (length/width) is varied from 2.5 to 120 by changing the separation between the surfaces, the width of the strips, or the fluid volume. As the bridge height is increased, the aspect ratio decreases and we observe a large increase in the mean curvature of the bridge. More specifically, the following counterintuitive result is observed: the mean curvature of the bridges changes sign and goes from negative (concave bridge) to positive (convex bridge) when the height is increased at constant volume. These experimental observations are in quantitative agreement with Surface Evolver simulations. Scaling shows a collapse of the data indicating that this transition in the sign of the Laplace pressure is universal for capillary bridges with high aspect ratios. Finally, we show that the morphology diagrams obtained from our 3D analysis are considerably different from those expected from a 2D analysis.

Flow-induced shear instabilities of cohesive granulates

In this work we use a multiscale framework to calculate the fluidization threshold of three-dimensional cohesive granulates under shear forces exerted by a creeping flow. A continuum model of flow through porous media provides an analytical expression for the average drag force on a single grain. The balance equation for the forces and a force propagation model are then used to investigate the effects of porosity and packing structure on the stability of the pile. We obtain a closed-form expression for the instability threshold of a regular packing of monodisperse frictionless cohesive spherical grains in a planar fracture. Our result quantifies the compound effect of structural (packing orientation and porosity) and dynamical properties of the system on its stability.