skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Average stress in a dilute suspension of rigid spheroids in a second-order fluid in a linear flow
The microhydrodynamics of particle suspensions in polymeric fluids has a wide range of applications in industry and biology. To discern the dynamics of particles in such systems, it is important to analyze the stress response of the suspension to applied flow fields. While such investigations have been theoretically done for suspensions of rigid spheres in weakly viscoelastic fluids, the effect of nonsphericity of particles on the stress remains relatively unexplored. The interplay between the response of the polymeric fluid and the particle orientation yields rich physics. The viscoelastic torques make the particle inhabit a preferred orientation in a given flow, resulting in time-dependent stresses. In this paper, we determine the average extra stress in a dilute suspension of rigid, non-Brownian spheroids in a second-order fluid subject to shear and extensional flows. We perform this task by examining the flow around a single spheroid in the limit of small Weissenberg number (Wi≪1) and perform an ensemble average of the stress tensor over all particle configurations. There are two contributions to the extra stress: one from the force dipole on the particles (stresslet) and another from the fluctuations in the velocity in the bulk fluid (fluid-induced particle stress), the latter of which does not arise in a zero Reynolds number Newtonian fluid. We present results for the 𝑂⁡(𝜙⁢Wi) corrections to the long-time effective shear viscosity, normal stress coefficients, and extensional viscosities in the suspension in shear, uniaxial extensional, and planar extensional flows, where 𝜙 is the particle volume fraction. To elucidate the effect of particle shape on the effective viscosity, we repeat this analysis for different aspect ratios (𝐴𝑅) for prolate (needlelike) and oblate (disklike) spheroids.  more » « less
Award ID(s):
2341154
PAR ID:
10647837
Author(s) / Creator(s):
; ;
Publisher / Repository:
American Physical Society
Date Published:
Journal Name:
Physical Review Fluids
Volume:
10
Issue:
5
ISSN:
2469-990X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This dataset is for the following study: Anomalous crystalline ordering of particles in a viscoelastic fluid under high shear by Sijie Sun, Nan Xue, Stefano Aime, Hyoungsoo Kim, Jizhou Tang, Gareth H. McKinley, Howard A. Stone, and David A. Weitz. In this research, the authors investigate the high-shear-rate behavior of particle suspensions in viscoelastic fluids using a fully immersed parallel plate geometry. They discover an unexpected particle separation within the suspension, leading to the formation of a crystalline layer in the center of the cell. This solid layer disrupts the flow instability and introduces a new single-frequency component to the torque fluctuations, reflecting a dominant velocity pattern in the flow. The findings reveal the complex interplay between particles and the suspending viscoelastic fluid under extreme shear conditions. The SI videos, raw data and code of postprocessing are included 
    more » « less
  2. 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
  3. This paper evaluates the behavior of a single rigid ellipsoidal particle suspended in homogeneous viscous flow with a power-law generalized Newtonian fluid rheology using a custom-built finite element analysis (FEA) simulation. The combined effects of the shear-thinning fluid rheology, the particle aspect ratio, the initial particle orientation, and the shear-extensional rate factor in various homogeneous flow regimes on the particles dynamics and surface pressure evolution are investigated. The shear-thinning fluid behavior was found to modify the particle’s trajectory and alter the particle’s kinematic response. Moreover, the pressure distribution over the particle’s surface is significantly reduced by the shear-thinning fluid rheology. The FEA model is validated by comparing results of the Newtonian case with results obtained from the well-known Jeffery’s analytical model. Furthermore, Jeffery’s model is extended to define the particle’s trajectory in a special class of homogeneous Newtonian flows with combined extension and shear rate components typically found in axisymmetric nozzle flow contractions. The findings provide an improved understanding of key transport phenomenon related to physical processes involving fluid–structure interaction such as that which occurs within the flow field developed during material extrusion–deposition additive manufacturing of fiber reinforced polymeric composites. These results provide insight into important microstructural formations within the print beads. 
    more » « less
  4. The stability of a rigid particle in yield stress fluids, comprised of soft particle glasses (SPGs), is investigated in shear flow under an applied external force, such as weight, using particle dynamics simulations. Results provide the critical force threshold, in terms of the dynamic yield stress and the flow strength, required to initiate sedimentation of the rigid particle over a wide range of shear rates and volume fractions. The streamlines of the SPGs show local disturbances when the rigid particle settles. The form of these disturbances is consistent with the microdynamics and microstructure response of the neighboring soft particles of the sedimenting rigid particle. Sedimenting particle induces non-affine displacement to the suspensions at low shear rates and high applied forces, while these dynamical events are localized and suppressed at high shear rates. Stability diagrams, which provide the conditions of the sedimentation of the rigid particle, are presented in terms of the applied force and the shear rate. These individual stability diagrams at each volume fraction map onto a universal stability diagram when the external force is scaled by the dynamic yield stress and shear rate with a ratio of the solvent viscosity to the low-frequency modulus of the SPGs. 
    more » « less
  5. A ring with a cross-section that has a blunt inner and sharper outer edge can attain an equilibrium orientation in a Newtonian fluid subject to a low Reynolds number simple shear flow. This may be contrasted with the continuous rotation exhibited by most rigid bodies. Such rings align along an orientation when the rotation due to fluid vorticity balances the counter-rotation due to the extensional component of the simple shear flow. While the viscous stress on the particle tries to rotate it, the pressure can generate a counter-vorticity torque that aligns the particle. Using boundary integral computations, we demonstrate ways to effectively control this pressure by altering the geometry of the ring cross-section, thus leading to alignment at moderate particle aspect ratios. Aligning rings that lack fore–aft symmetry can migrate indefinitely along the gradient direction. This differs from the periodic spatial trajectories of fore–aft asymmetric axisymmetric particles that rotate in periodic orbits. The mechanism for migration of aligned rings along the gradient direction is elucidated in this work. The migration speed can be controlled by varying the cross-sectional shape and size of the ring. Our results provide new insights into controlling motion of individual particles and thereby open new pathways towards manipulating macroscopic properties of a suspension. 
    more » « less