skip to main content


Title: Ordered domains in sheared dense suspensions: The link to viscosity and the disruptive effect of friction
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
Award ID(s):
2026842
PAR ID:
10376796
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Rheology
Volume:
66
Issue:
5
ISSN:
0148-6055
Page Range / eLocation ID:
1055 to 1065
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. We investigate the rheological behavior of athermal particle suspensions using experiments and theory. A generalized version of the homogenization estimates of Ponte Castañeda and Willis [J. Mech. Phys. Solids, 43(12), 1919–1951 (1995)] is presented for the effective viscosity of athermal suspensions accounting for additional microstructural features (e.g., polydispersity) via an empirical parameter, [Formula: see text]. For the case of identically sized spheres dispersed with statistical isotropy in a Newtonian fluid, the parameter [Formula: see text] is estimated from the results of Batchelor and Green [J. Fluid Mech. 56(2), 375–400 (1972)] for the Huggins coefficient. Predictions for the macroscopic viscosity are found to be in good agreement with measurements for monodisperse polymethyl methacrylate (PMMA) spheres in glycerol, as well as for the empirical Krieger–Dougherty equation for the shear viscosity. The proposed estimates have the added benefit that they can also be used to get information on the statistics of the stress and strain-rate fields in the fluid and particle phases. In addition, results for the effective shear viscosity are used in combination with the linear comparison method of Ponte Castañeda [J. Mech. Phys. Solids 39(1), 45–71 (1991)] to generate the corresponding estimates for the effective macroscopic behavior and field statistics of particle suspensions in (viscoplastic) yield stress fluids. Good agreement is also found between the theoretical estimates and experimental results for the effective yield and flow stress of suspensions with monodisperse PMMA spheres in Carbopol. Finally, it is argued that the results for the phase averages and fluctuations of the stress and strain-rate fields can be used to provide a physical interpretation for the parameter [Formula: see text] in terms of the polydispersity of the suspension and its implications for the percolation threshold.

     
    more » « less
  4. We have determined the susceptibility of T4 DNA (166 kilobase pairs, kbp) to fragmentation under steady shear in a cone-and-plate rheometer. After shearing for at least 30 min at a shear rate of [Formula: see text], corresponding to a Reynolds number of [Formula: see text] and a Weissenberg number of [Formula: see text], [Formula: see text]% of the sample is broken into a polydisperse mixture with a number-averaged molecular weight of [Formula: see text] kbp and a polydispersity index of [Formula: see text], as measured by pulsed-field gel electrophoresis (with a 95% confidence interval). The molecular weight distributions observed here from a shear flow are similar to those produced by a (dominantly extensional) sink flow of DNA and are qualitatively different than the midpoint scission observed in simple extensional flow. Given the inability of shear flow to produce a sharp coil–stretch transition, the data presented here support a model where polymers can be fragmented in flow without complete extension. These results further indicate that DNA fragmentation by shear is unlikely to be a significant issue in microfluidic devices, and anomalous molecular weight observations in experiments are due to DNA processing prior to observation in the device. 
    more » « less
  5. Fungi such as Candida albicans exist in biofilm phenotypes, which present as viscoelastic materials; however, a method to measure linear viscoelastic moduli, yield stress, and yield strain is lacking. Characterization methods for fungal materials have been limited to techniques specific to particular industries. Here, we present a method to measure the shear stress, strain amplitude, and creep of C. albicans BWP17 biofilms. Our method includes features tailored to the analysis of fungi including an in vitro growth protocol attuned to the slow growth rates of C. albicans biofilms and a resultant cultured biofilm that has sufficient integrity to be transferred to the rheometer tooling without disrupting its structure. The method's performance is demonstrated by showing that results are insensitive to gap, evaporative sealant, length of experiment, and specimen radius. Multiscale imaging of the fungal biofilm showed complex entanglement networks at the hundred-micrometer scale. For a wild-type strain cultivated for 14 days, using small-amplitude oscillatory rheology, we found that the elastic (G′) and viscous (G″) moduli were nearly independent of frequency over the range 0.1–10 s −1 , with magnitudes of [Formula: see text] and [Formula: see text], respectively. The yield stress was approximately [Formula: see text]. We modeled the linear creep response of the fungal biofilm and found that C. albicans has a characteristic relaxation time of [Formula: see text] and a viscosity of [Formula: see text]. We applied this method to probe the effects of altered chitin deposition in the C. albicans cell wall. Differences between the biofilm's phenotypic cell shape and rheological properties in mutants with altered chitin synthase activity were resolved. Discovering how genotypic, phenotypic, and environmental factors impact the material properties of these microbial communities can have implications for understanding fungal biofilm growth and aid in the development of remediation strategies. 
    more » « less