Rheology of a dilute viscoelastic suspension of spheroids in unconfined shear flow

Abstract

The rheology of a dilute viscoelastic suspension of spheroids subjected to unconfined shear flow is studied by numerical simulations. To highlight the effect of the suspending fluid rheology, two viscoelastic constitutive equations, i.e., the Giesekus and the Phan-Thien-Tanner models, have been selected. Simulations are performed for a spheroid with two aspect ratios (4 and 8). The spherical particle case is also investigated for comparison. The Deborah number D e is varied between 0 and 4. The particle contribution to the viscosity is weakly affected by the particle shape and orientation. In contrast, spheroids oriented with major axis out of the vorticity direction significantly reduce the particle contribution to the first and second normal stress difference. The maximum reduction is found in the flow-alignment regime. Simulations of the transient dynamics of the suspension show that the initial distribution of the particle orientation has a remarkable influence on the evolution of the stress coefficients. The time needed to obtain steady-state rheological properties is at least one order of magnitude higher than the fluid characteristic time

A 3D DLM/FD method for simulating the motion of spheres in a bounded shear flow of Oldroyd-B fluids

“…For Wi=0.75 and 1, those numerical results of the kayaking motion suggest that for a rigid long body rotating in a bounded shear flow of an Oldroyd-B fluid, its long axis migrates out of the shear plane. Those numerical results do not conflict with those of an ellipsoid in a bounded shear flow of a Giesekus fluid obtained by D'Avino et al in [46] and [47] since they did not include the effect of the particle and fluid inertia in their simulations.…”

On the rheology of particle suspensions in viscoelastic fluids

“…Particle interactions have a profound effect on these orientation distributions as measured experimentally and broad distributions extending all the way from those peaked at the vorticity direction (elasticity dominated) to those peaked along the flow direction (particle interaction dominated) have been measured, see Iso et al. The only existing publication heretofore, which attempts to relate these orientation distributions to rheology, is the simulation paper of Avino et al However, these authors apparently neglected the particle induced fluid stress in making their rheological calculations, and therefore this must be revisited. Note that all of these effects, as they relate to suspension rheology, are now accessible with the numerical methods and mathematical foundation discussed in the earlier section of this viewpoint and certainly will be examined in the coming years.…”

Rotational dynamics of a neutrally buoyant prolate spheroid in viscoelastic shear flows at finite Reynolds numbers

“…D'Avino et al. (2014) and D'Avino, Greco & Maffettone (2015) systemically analysed the detailed dynamics of prolate particles in a sheared Giesekus fluid by finite element simulations. Due to the fluid elasticity, the particle rotation could be categorized into four different modes: (i) in low Deborah number flow, the prolate particles align their symmetry axes along the vorticity direction and behave in a ‘log-rolling’ mode.…”