skip to main content

Title: Contact criterion for suspensions of smooth and rough colloids
We report a procedure to obtain the search distance used to determine particle contact in dense suspensions of smooth and rough colloids. This method works by summing physically relevant length scales in an uncertainty analysis and does not require detailed quantification of the surface roughness. We suspend sterically stabilized, fluorescent poly(methyl methacrylate) colloids in a refractive index-matched solvent, squalene, in order to ensure hard sphere-like behavior. High speed centrifugation is used to pack smooth and rough colloids to their respective jamming points, ϕ J . The jammed suspensions are subsequently diluted with known volumes of solvent to ϕ < ϕ J . Structural parameters obtained from confocal laser scanning micrographs of the diluted colloidal suspensions are extrapolated to ϕ J to determine the mean contact number at jamming, 〈 z 〉 J . Contact below jamming refers to nearest neighbors at a length scale below which the effects of hydrodynamic or geometric friction come into play. Sensitivity analyses show that a deviation of the search distance by 1% of the particle diameter results in 〈 z 〉 changing by up to 10%, with the error in contact number distribution being magnified in dense suspensions ( ϕ > 0.50) due to more » an increased number of nearest neighbors in the first coordination shell. « less
Award ID(s):
Publication Date:
Journal Name:
Soft Matter
Page Range or eLocation-ID:
4980 to 4989
Sponsoring Org:
National Science Foundation
More Like this
  1. We report the linear rheology for dense suspensions of sterically stabilized smooth and mesoscopically rough colloids interacting as hard particles. Small amplitude oscillatory measurements reveal that rough colloids at high volume fractions exhibit storage and loss moduli that are orders of magnitude greater than smooth colloids. Frequency-concentration superposition is used to collapse the viscoelasticity data onto a master curve, where shift factors suggest a more elastic microstructure and reduced cage volume for rough particles. A combination of the mode-coupling theory, hydrodynamic modeling, and the activated hopping theory shows that these rough particles with significantly reduced localization lengths tend to become trapped in their glassy cages for extended periods of time. High-frequency data show that rough colloids, but not smooth colloids, display a transition from a free-draining to a fully lubricated state above the crossover volume fraction and, furthermore, exhibit solidlike behavior. Scaling analyses support the idea that lubrication forces between interlocking asperities are enhanced, leading to rotational constraints and stress-bearing structures that significantly elevate the viscoelasticity of dense suspensions. The results provide a framework for how particle surface topology affects the linear rheology in applications such as coatings, cement, consumer products, and shock-absorbing materials.
  2. Debris flows are dense and fast-moving complex suspensions of soil and water that threaten lives and infrastructure. Assessing the hazard potential of debris flows requires predicting yield and flow behavior. Reported measurements of rheology for debris flow slurries are highly variable and sometimes contradictory due to heterogeneity in particle composition and volume fraction ( ϕ ) and also inconsistent measurement methods. Here we examine the composition and flow behavior of source materials that formed the postwildfire debris flows in Montecito, CA, in 2018, for a wide range of ϕ that encapsulates debris flow formation by overland flow. We find that shear viscosity and yield stress are controlled by the distance from jamming, Δ ϕ = ϕ m − ϕ , where the jamming fraction ϕ m is a material parameter that depends on grain size polydispersity and friction. By rescaling shear and viscous stresses to account for these effects, the data collapse onto a simple nondimensional flow curve indicative of a Bingham plastic (viscoplastic) fluid. Given the highly nonlinear dependence of rheology on Δ ϕ , our findings suggest that determining the jamming fraction for natural materials will significantly improve flow models for geophysical suspensions such as hyperconcentrated flows andmore »debris flows.« less
  3. The remarkable increase in the flow resistance of dense suspensions can hinder 3D-printing processes on account of flow cessation in the extruder, and filament fragility/rupture following deposition. Understanding the nature of rheological changes that occur is critical to manipulate flow conditions or to dose flow modifiers for 3D-printing. Therefore, this paper elucidates the influences of clay particulates on controlling flow cessation and the shape stability of dense cementing suspensions that typically feature poor printability. A rope coiling method was implemented with varying stand-off distances to probe the buckling stability and tendency to fracture of dense suspensions that undergo stretching and bending during deposition. The contributions of flocculation and short-term percolation due to the kinetics of structure formation to deformation rate were deconvoluted using a stepped isostress method. It is shown that the shear stress indicates a divergence with a power-law scaling when the particle volume fraction approaches the jamming limit; ϕ → ϕ j ≈ ϕ max . Such a power-law divergence of the shear stress decreases by a factor of 10 with increasing clay dosage. Such behavior in clay-containing suspensions arises from a decrease in the relative packing fraction ( ϕ / ϕ max ) and the formation ofmore »fractally-architected aggregates with stronger interparticle interactions, whose uniform arrangement controls flow cessation in the extruder and suspension homogeneity, thereby imparting greater buckling stability. The outcomes offer new insights for assessing/improving the extrudability and printability behavior during slurry-based 3D-printing process.« less
  4. This study explores thermal convection in suspensions of neutrally buoyant, non-colloidal suspensions confined between horizontal plates. A constitutive diffusion equation is used to model the dynamics of the particles suspended in a viscous fluid and it is coupled with the flow equations. We employ a simple model that was proposed by Metzger, Rahli & Yin ( J. Fluid Mech. , vol. 724, 2013, pp. 527–552) for the effective thermal diffusivity of suspensions. This model considers the effect of shear-induced diffusion and gives the thermal diffusivity increasing linearly with the thermal Péclet number ( Pe ) and the particle volume fraction ( ϕ ). Both linear stability analysis and numerical simulation based on the mathematical models are performed for various bulk particle volume fractions $({\phi _b})$ ranging from 0 to 0.3. The critical Rayleigh number $(R{a_c})$ grows gradually by increasing ${\phi _b}$ from the critical value $(R{a_c} = 1708)$ for a pure Newtonian fluid, while the critical wavenumber $({k_c})$ remains constant at 3.12. The transition from the conduction state of suspensions is subcritical, whereas it is supercritical for the convection in a pure Newtonian fluid $({\phi _b} = 0)$ . The heat transfer in moderately dense suspensions $({\phi _b} = 0.2\text{--}0.3)$more »is significantly enhanced by convection rolls for small Rayleigh number ( Ra ) close to $R{a_c}$ . We also found a power-law increase of the Nusselt number ( Nu ) with Ra , namely, $Nu\sim R{a^b}$ for relatively large values of Ra where the scaling exponent b decreases with ${\phi _b}$ . Finally, it turns out that the shear-induced migration of particles can modify the heat transfer.« less
  5. We perform computational studies of jammed particle packings in two dimensions undergoing isotropic compression using the well-characterized soft particle (SP) model and deformable particle (DP) model that we developed for bubbles and emulsions. In the SP model, circular particles are allowed to overlap, generating purely repulsive forces. In the DP model, particles minimize their perimeter, while deforming at the fixed area to avoid overlap during compression. We compare the structural and mechanical properties of jammed packings generated using the SP and DP models as a function of the packing fraction ρ, instead of the reduced number density φ. We show that near jamming onset the excess contact number Δz=z-z J and shear modulus G scale as Δρ 0.5 in the large system limit for both models, where Δρ=ρ-ρ J and z J ≈4 and ρ J ≈0.842 are the values at jamming onset. Δz and G for the SP and DP models begin to differ for ρ≥0.88. In this regime, Δz∼G can be described by a sum of two power-laws in Δρ, i.e. Δz∼G∼C 0 Δρ 0.5 +C 1 Δρ 1.0 to lowest order. We show that the ratio C 1 /C 0 is much larger for the DP modelmore »compared to that for the SP model. We also characterize the void space in jammed packings as a function of ρ. We find that the DP model can describe the formation of Plateau borders as ρ→1.0. We further show that the results for z and the shape factor A versus ρ for the DP model agree with recent experimental studies of foams and emulsions.« less