Search for: All records

Dense non-Brownian suspensions with conservative repulsive forces between the particles are known to exhibit shear thickening, where viscosity increases with applied stress due to a change in the dominant stress mechanism. At low stress, repulsion maintains liquid films that lubricate particle interactions, while higher stress overcomes the repulsion to generate frictional contacts and leads to greater flow resistance. Here, shear-thickened suspensions are studied in stress-controlled simulations incorporating hydrodynamic, electrostatic double-layer repulsion, and frictional forces; two-dimensional monolayers are studied for monodisperse and bidisperse suspensions with size ratios δ=as/al from 1.0 to 4.0, where as and al are the small and large particle radii. Small-particle fractions ζ=ϕs/ϕ=0.25,0.50, and 0.75 are considered. Total area fractions of 0.71≤ϕ≤0.82 are studied, with the larger values at greater size ratios. Flow curves for mono- and bidisperse systems under varying stress are analyzed, along with detailed structural comparisons for different interparticle friction. We examine the approach to shear jamming, through the emergence of rigid local clusters generated by the reduction of degrees of freedom by frictional contacts. The variance of the fraction of particles in rigid clusters increases sharply near the jamming solid fraction, consistent with a second-order phase transition description of the phenomenon. The contact fabric tensor is determined to provide a measure of the structural anisotropy. 

This paper focuses on the origin and implications of particle pressure and discontinuous shear thickening in concentrated suspensions. These properties are both related to the tendency of a flowing suspension to exert normal forces on the confining boundaries, thus providing a conceptual relation of the two seemingly distinct issues through a consideration of the pressure-volume relation of a flowing suspension. An overview of basic elements of suspension mechanics related to these topics is presented, including microstructure and continuum formulations based on single-phase and two-phase perspectives. The historical development of understanding of particle pressure and its influence on particle migration and that of discontinuous shear thickening are described. The mechanistic basis for the particle pressure in terms of suspension microstructure and the role of frictional contact interactions in shear thickening are described. A few open questions related to these topics are presented in conclusion. 

The onset and growth of rigid clusters in a two-dimensional (2D) suspension in shear flow are studied by numerical simulations. The suspension exhibits the lubricated-to-frictional rheology transition, but the key results here are for stresses above the levels that cause extreme shear thickening. At large solid fraction, $\varphi$, but below the stress-dependent jamming fraction, we find a critical ${\varphi }_c(\sigma ,\mu )$where $\sigma$is a dimensionless shear stress and $\mu$is the interparticle friction coefficient. For $\varphi >{\varphi }_c$, the proportion of particles in rigid clusters grows sharply, as $f_{\mathrm{rig}}\sim {|\varphi -{\varphi }_c|}^{\beta }$with $\beta =1/8$. The fluctuations in the fraction of particles in rigid clusters yield a susceptibility measure ${\chi }_{\mathrm{rig}}\sim {|\varphi -{\varphi }_c|}^{-\gamma }$with $\gamma =7/4$. The system is thus found to exhibit criticality. The results are shown to depend on an effective field $h\left(\mu \right)$, which provides data collapse near ${\varphi }_c$for both $f_{\mathrm{rig}}$and ${\chi }_{\mathrm{rig}}$. This behavior occurs over a range of stresses, with ${\varphi }_c(\sigma ,\mu )$increasing as the stress decreases. 

Developments in the last century, and especially in the last 50 years, have advanced understanding of suspension rheology greatly. Here, a limited review of suspension work over this period is presented, emphasizing advances over the last three decades in understanding of the particle pressure and strong shear thickening, which were motivated by crucial experimental observations, computational advances, and a critical review, all from the 1980s. This review serves as a preview to some outstanding challenges in suspension mechanics. This article considers primarily dispersions of spherical particles, which serve not only as a model material for understanding the rheology of more complex fluids of practical relevance, but also as a basic system for the study of nonequilibrium statistical physics. 

We studied the evolution of capillary bridges between nominally flat plates undergoing multiple cycles of compression and stretching in experiments and simulations. We varied the distance between the plates in small increments to study the full evolution of the bridge shape. Experiments show that contact angle hysteresis determines the shape of the bridge. In sliding drops, hysteresis can be modeled using a contact angle-dependent resistive force F̃R applied at the contact line. We developed a model that accurately captures the evolution of the bridge shape by combining F̃R and constrained energy minimization. Unlike previous work, this allows for both complete and partial contact line pinning. We also explored the effect of using nonparallel plates. The asymmetry in the bridge shape causes the movement of the center of mass of the bridge and can be explained by contact angle hysteresis. We find that even a slight misalignment between the flat plates can have a measurable effect. 

Discrete-particle simulations of bidisperse shear thickening suspensions are reported. The work considers two packing parameters, the large-to-small particle radius ratio ranging from [Formula: see text] (nearly monodisperse) to [Formula: see text], and the large particle fraction of the total solid loading with values [Formula: see text], 0.5, and 0.85. Particle-scale simulations are performed over a broad range of shear stresses using a simulation model for spherical particles accounting for short-range lubrication forces, frictional interaction, and repulsion between particles. The variation of rheological properties and the maximum packing fraction [Formula: see text] with shear stress [Formula: see text] are reported. At a fixed volume fraction [Formula: see text], bidispersity decreases the suspension relative viscosity [Formula: see text], where [Formula: see text] is the suspension viscosity and [Formula: see text] is the suspending fluid viscosity, over the entire range of shear stresses studied. However, under low shear stress conditions, the suspension exhibits an unusual rheological behavior: the minimum viscosity does not occur as expected at [Formula: see text], but instead decreases with further increase of [Formula: see text] to [Formula: see text]. The second normal stress difference [Formula: see text] acts similarly. This behavior is caused by particles ordering into a layered structure, as is also reflected by the zero slope with respect to time of the mean-square displacement in the velocity gradient direction. The relative viscosity [Formula: see text] of bidisperse rate-dependent suspensions can be predicted by a power law linking it to [Formula: see text], [Formula: see text] in both low and high shear stress regimes. The agreement between the power law and experimental data from literature demonstrates that the model captures well the effect of particle size distribution, showing that viscosity roughly collapses onto a single master curve when plotted against the reduced volume fraction [Formula: see text]. 

A colloidal motor driven by surface tension forces is theoretically designed by encapsulating an active Janus particle in a liquid drop which is immiscible in the suspending medium. The Janus particle produces an asymmetric flux of a solute species which induces surface tension gradients along the liquid–liquid interface between the drop and the surrounding fluid. The resulting Marangoni forces at the interface propel the compound drop/Janus particle system. The propulsion speeds of the motor are evaluated for a range of relative sizes and positions of the drop and the particle and across a range of transport properties of the drop and the suspending medium. It is demonstrated that the proposed design can produce higher propulsion velocities than the traditional Janus-particle-based colloidal motors propelled by neutral diffusiophoresis. 

Dense suspensions of particles in viscous liquid often demonstrate the striking phenomenon of abrupt shear thickening, where their viscosity increases strongly with increase of the imposed stress or shear rate. In this work, discrete-particle simulations accounting for short-range hydrodynamic, repulsive, and contact forces are performed to simulate flow of shear thickening bidisperse suspensions, with the packing parameters of large-to-small particle radius ratio δ = 3 and large particle fraction ζ = 0.15, 0.50, and 0.85. The simulations are carried out for volume fractions 0.54 ≤ ϕ ≤ 0.60 and a wide range of shear stresses. The repulsive forces, of magnitude F R , model the effects of surface charge and electric double-layer overlap, and result in shear thinning at small stress, with shear thickening beginning at stresses σ ∼ F R a −2 . A crossover scaling analysis used to describe systems with more than one thermodynamic critical point has recently been shown to successfully describe the experimentally-observed shear thickening behavior in suspensions. The scaling theory is tested here on simulated shear thickening data of the bidisperse mixtures, and also on nearly monodisperse suspensions with δ = 1.4 and ζ = 0.50. Presenting the viscosity in terms of a universal crossover scaling function between the frictionless and frictional maximum packing fractions collapses the viscosity for most of the suspensions studied. Two scaling regimes having different exponents are observed. The scaling analysis shows that the second normal stress difference N 2 and the particle pressure Π also collapse on their respective curves, with the latter featuring a different exponent from the viscosity and normal stress difference. The influence of the fraction of frictional contacts, one of the parameters of the scaling analysis, and its dependence on the packing parameters are also presented.