Understanding the orientation dynamics of anisotropic colloidal particles is important for suspension rheology and particle self-assembly. However, even for the simplest case of dilute suspensions in shear flow, the orientation dynamics of non-spherical Brownian particles are poorly understood. Here we analytically calculate the time-dependent orientation distributions for non-spherical axisymmetric particles confined to rotate in the flow–gradient plane, in the limit of small but non-zero Brownian diffusivity. For continuous shear, despite the complicated dynamics arising from the particle rotations, we find a coordinate change that maps the orientation dynamics to a diffusion equation with a remarkably simple ratio of the enhanced rotary diffusivity to the zero shear diffusion: $D_{eff}^{r}/D_{0}^{r}=(3/8)(p-1/p)^{2}+1$ , where $p$ is the particle aspect ratio. For oscillatory shear, the enhanced diffusion becomes orientation dependent and drastically alters the long-time orientation distributions. We describe a general method for solving the time-dependent oscillatory shear distributions and finding the effective diffusion constant. As an illustration, we use this method to solve for the diffusion and distributions in the case of triangle-wave oscillatory shear and find that they depend strongly on the strain amplitude and particle aspect ratio. These results provide new insight into the time-dependent rheology of suspensions of anisotropic particles. For continuous shear, we find two distinct diffusive time scales in the rheology that scale separately with aspect ratio $p$ , as $1/D_{0}^{r}p^{4}$ and as $1/D_{0}^{r}p^{2}$ for $p\gg 1$ . For oscillatory shear flows, the intrinsic viscosity oscillates with the strain amplitude. Finally, we show the relevance of our results to real suspensions in which particles can rotate freely. Collectively, the interplay between shear-induced rotations and diffusion has rich structure and strong effects: for a particle with aspect ratio 10, the oscillatory shear intrinsic viscosity varies by a factor of ${\approx}2$ and the rotational diffusion by a factor of ${\approx}40$ .
more »
« less
Shear rheology of a dilute suspension of thin rings
The rheology of suspensions of rings (tori) rotating in an unbounded low Reynolds number simple shear flow is calculated using numerical simulations at dilute particle number densities ( n ≪ 1 ). Suspensions of non-Brownian rings are studied by computing pair interactions that include hydrodynamic interactions modeled using slender body theory and particle collisions modeled using a short-range repulsive force. Particle contact and hydrodynamic interactions were found to have comparable influences on the steady-state Jeffery orbit distribution. The average tilt of the ring away from the flow-vorticity plane increased during pairwise interactions compared to the tilt associated with Jeffery rotation and the steady-state orbit distribution. Particle stresses associated with the increased tilt during the interaction were found to be comparable to the stresses induced directly by particle contact forces and the hydrodynamic velocity disturbances of other particles. The hydrodynamic diffusivity coefficients in the gradient and vorticity directions were also obtained and were found to be two orders of magnitude larger than the corresponding values in fiber suspensions at the same particle concentrations. Rotary Brownian dynamics simulations of isolated Brownian rings were used to understand the shear rate dependence of suspension rheology. The orbit distribution observed in the regime of weak Brownian motion, P e ≫ ϕ T − 3, was surprisingly similar to that obtained from pairwise interaction calculations of non-Brownian rings. Here, the Peclet number P e is the ratio of the shear rate and the rotary diffusivity of the particle and ϕ T is the effective inverse-aspect ratio of the particle (approximately equal to 2 π times the inverse of its non-dimensional Jeffery time period). Thus, the rheology results obtained from pairwise interactions should retain accuracy even for weakly Brownian rings ( n ≪ 1 and ϕ T − 3 ≪ P e ).
more »
« less
- Award ID(s):
- 2206851
- NSF-PAR ID:
- 10417821
- Date Published:
- Journal Name:
- Journal of Rheology
- Volume:
- 67
- Issue:
- 3
- ISSN:
- 0148-6055
- Page Range / eLocation ID:
- 723
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Three-phase suspensions, of liquid that suspends dispersed solid particles and gas bubbles, are common in both natural and industrial settings. Their rheology is poorly constrained, particularly for high total suspended fractions (≳0.5). We use a dam-break consistometer to characterize the rheology of suspensions of (Newtonian) corn syrup, plastic particles and CO 2 bubbles. The study is motivated by a desire to understand the rheology of magma and lava. Our experiments are scaled to the volcanic system: they are conducted in the non-Brownian, non-inertial regime; bubble capillary number is varied across unity; and bubble and particle fractions are 0 ≤ ϕ gas ≤ 0.82 and 0 ≤ ϕ solid ≤ 0.37, respectively. We measure flow-front velocity and invert for a Herschel–Bulkley rheology model as a function of ϕ gas , ϕ solid , and the capillary number. We find a stronger increase in relative viscosity with increasing ϕ gas in the low to intermediate capillary number regime than predicted by existing theory, and find both shear-thinning and shear-thickening effects, depending on the capillary number. We apply our model to the existing community code for lava flow emplacement, PyFLOWGO, and predict increased viscosity and decreased velocity compared with current rheological models, suggesting existing models may not adequately account for the role of bubbles in stiffening lavas.more » « less
-
Hydrodynamic interactions generate a diffusive motion in particulates in a shear flow, which plays seminal roles in overall particulate rheology and its microstructure. Here we investigate the shear induced diffusion in a red-blood cell (RBC) suspension using a numerical simulation resolving individual motion and deformation of RBCs. The non-spherical resting shape of RBCs gives rise to qualitatively different regimes of cell dynamics in a shear flow such as tank-treading, breathing, tumbling and swinging, depending on the cell flexibility determined by the elastic capillary number. We show that the transition from tumbling to tank-treading causes a reduction in the gradient diffusivity. The diffusivity is computed using a continuum approach from the evolution of a randomly packed cell-layer width with time as well as by the dynamic structure factor of the suspension. Both approaches, although operationally different, match and show that for intermediate capillary numbers RBCs cease tumbling accompanied by a drop in the coefficient of gradient diffusivity. A further increase of capillary number increases the diffusivity due to increased deformation. The effects of bending modulus and viscosity ratio variations are also briefly investigated. The computed shear induced diffusivity was compared with values in the literature. Apart from its effects in margination of cells in blood flow and use in medical diagnostics, the phenomenon broadly offers important insights into suspensions of deformable particles with non-spherical equilibrium shapes, which also could play a critical role in using particle flexibility for applications such as label free separation or material processing.more » « less
-
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 and debris flows.more » « less
-
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.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1122/8.0000628
