skip to main content


Title: Order and density fluctuations near the boundary in sheared dense suspensions
We introduce a novel approach to reveal ordering fluctuations in sheared dense suspensions, using line scanning in a combined rheometer and laser scanning confocal microscope. We validate the technique with a moderately dense suspension, observing modest shear-induced ordering and a nearly linear flow profile. At high concentration ( ϕ = 0.55) and applied stress just below shear thickening, we report ordering fluctuations with high temporal resolution, and directly measure a decrease in order with distance from the suspension’s bottom boundary as well as a direct correlation between order and particle concentration. Higher applied stress produces shear thickening with large fluctuations in boundary stress which we find are accompanied by dramatic fluctuations in suspension flow speeds. The peak flow rates are independent of distance from the suspension boundary, indicating that they likely arise from transient jamming that creates solid-like aggregates of particles moving together, but only briefly because the high speed fluctuations are interspersed with regions flowing much more slowly, suggesting that shear thickening suspensions possess complex internal structural dynamics, even in relatively simple geometries.  more » « less
Award ID(s):
1809890
NSF-PAR ID:
10445505
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Frontiers in Physics
Volume:
10
ISSN:
2296-424X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We report direct measurements of spatially resolved stress at the boundary of a shear-thickening cornstarch suspension revealing persistent regions of high local stress propagating in the flow direction at the speed of the top boundary. The persistence of these propagating fronts enables precise measurements of their structure, including the profile of boundary stress measured by boundary stress microscopy (BSM) and the nonaffine velocity of particles at the bottom boundary of the suspension measured by particle image velocimetry (PIV). In addition, we directly measure the relative flow between the particle phase and the suspending fluid (fluid migration) and find the migration is highly localized to the fronts and changes direction across the front, indicating that the fronts are composed of a localized region of high dilatant pressure and low particle concentration. The magnitude of the flow indicates that the pore pressure difference driving the fluid migration is comparable to the critical shear stress for the onset of shear thickening. The propagating fronts fully account for the increase in viscosity with applied stress reported by the rheometer and are consistent with the existence of a stable jammed region in contact with one boundary of the system that generates a propagating network of percolated frictional contacts spanning the gap between the rheometer plates and producing strong localized dilatant pressure. 
    more » « less
  2. The rheological behaviour of dense suspensions of ideally conductive particles in the presence of both electric field and shear flow is studied using large-scale numerical simulations. Under the action of an electric field, these particles are known to undergo dipolophoresis (DIP), which is the combination of two nonlinear electrokinetic phenomena: induced-charge electrophoresis (ICEP) and dielectrophoresis (DEP). For ideally conductive particles, ICEP is predominant over DEP, resulting in transient pairing dynamics. The shear viscosity and first and second normal stress differences$N_1$and$N_2$of such suspensions are examined over a range of volume fractions$15\,\% \leq \phi \leq 50\,\%$as a function of Mason number$Mn$, which measures the relative importance of viscous shear stress over electrokinetic-driven stress. For$Mn < 1$or low shear rates, the DIP is shown to dominate the dynamics, resulting in a relatively low-viscosity state. The positive$N_1$and negative$N_2$are observed at$\phi < 30\,\%$, which is similar to Brownian suspensions, while their signs are reversed at$\phi \ge 30\,\%$. For$Mn \ge 1$, the shear thickening starts to arise at$\phi \ge 30\,\%$, and an almost five-fold increase in viscosity occurs at$\phi = 50\,\%$. Both$N_1$and$N_2$are negative for$Mn \gg 1$at all volume fractions considered. We illuminate the transition in rheological behaviours from DIP to shear dominance around$Mn = 1$in connection to suspension microstructure and dynamics. Lastly, our findings reveal the potential use of nonlinear electrokinetics as a means of active rheology control for such suspensions.

     
    more » « less
  3. 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]. 
    more » « less
  4. Monodisperse suspensions of Brownian colloidal spheres crystallize at high densities, and ordering under shear has been observed at densities below the crystallization threshold. We perform large-scale simulations of a model suspension containing over [Formula: see text] particles to quantitatively study the ordering under shear and to investigate its link to the rheological properties of the suspension. We find that at high rates, for [Formula: see text], the shear flow induces an ordering transition that significantly decreases the measured viscosity. This ordering is analyzed in terms of the development of layering and planar order, and we determine that particles are packed into hexagonal crystal layers (with numerous defects) that slide past each other. By computing local [Formula: see text] and [Formula: see text] order parameters, we determine that the defects correspond to chains of particles in a squarelike lattice. We compute the individual particle contributions to the stress tensor and discover that the largest contributors to the shear stress are primarily located in these lower density, defect regions. The defect structure enables the formation of compressed chains of particles to resist the shear, but these chains are transient and short-lived. The inclusion of a contact friction force allows the stress-bearing structures to grow into a system-spanning network, thereby disrupting the order and drastically increasing the suspension viscosity. 
    more » « less
  5. null (Ed.)
    The geometric organization and force networks of 3D dense suspensions that exhibit both shear thinning and thickening have been examined as a function of varying strength of interparticle attractive interactions using lubrication flow discrete element simulations. Significant rearrangement of the geometric topology does not occur at either the local or global scale as these systems transition across the shear thinning and shear thickening regimes. In contrast, massive rearrangements in the balance of attractive, lubrication, and contact forces are observed with interesting behavior of network growth and competition. In agreement with prior work, in shear thinning regions the attractive force is dominant, however as the shear thickening region is approached there is growth of lubrication forces. Lubrication forces oppose the attraction forces, but as viscosity continues to increase under increasing shear stress, the lubrication forces are dominated by contact forces that also resist attraction. Contact forces are the dominant interactions during shear thickening and are an order of magnitude higher than their values in the shear-thinning regime. At high attractive interaction strength, contact networks can form even under shear thinning conditions, however high shear stress is still required before contact networks become the driving mechanism of shear thickening. Analysis of the contact force network during shear thickening generally indicates a uniformly spreading network that rapidly forms across empty domains; however the growth patterns exhibit structure that is significantly dependent upon the strength of interparticle interactions, indicating subtle variations in the mechanism of shear thickening. 
    more » « less