Effects of inertia and viscoelasticity on sedimenting anisotropic particles

Abstract

An axisymmetric particle sedimenting in an otherwise quiescent Newtonian fluid, in the Stokes regime, retains its initial orientation. For the special case of a spheroidal geometry, we examine analytically the effects of weak inertia and viscoelasticity in driving the particle towards an eventual steady orientation independent of initial conditions. The generalized reciprocal theorem, together with a novel vector spheroidal harmonics formalism, is used to find closed-form analytical expressions for the O(Re) inertial torque and the O(De) viscoelastic torque acting on a sedimenting spheroid of an arbitrary aspect ratio. Here, Re = UL/nu is the Reynolds number, with U being the sedimentation velocity, L the semi-major axis and nu the fluid kinematic viscosity, and is a measure of the inertial forces acting at the particle scale. The Deborah number, De = (lambda U)/L, is a dimensionless measure of the fluid viscoelasticity, with lambda being the intrinsic relaxation time of the underlying microstructure. The analysis is valid in the limit Re, De << 1, and the effects of viscoelasticity are therefore modelled using the constitutive equation of a second-order fluid. The inertial torque always acts to turn the spheroid broadside-on, while the final orientation due to the viscoelastic torque depends on the ratio of the magnitude of the first (N-1) to the second normal stress difference (N-2), and the sign (tensile or compressive) of N1. For the usual case of near-equilibrium complex fluids - a positive and dominant N-1 (N-1 > 0, N-2 < 0 and vertical bar N-1/N-2 vertical bar > 1) - both prolate and oblate spheroids adopt a longside-on orientation. The viscoelastic torque is found to be remarkably sensitive to variations in kappa in the slender-fibre limit (kappa >> 1), where kappa = L/b is the aspect ratio, b being the radius of the spheroid (semi-minor axis). The angular dependence of the inertial and viscoelastic torques turn out to be identical, and one may then characterize the long-time orientation of the sedimenting spheroid based solely on a critical value (El(c)) of the elasticity number, El = De/Re. For El < El(c) (> El(c)), inertia (viscoelasticity) prevails with the spheroid settling broadside-on (longside-on). The analysis shows that El(c) similar to O[(1/ln kappa)] for kappa >> 1, and the viscoelastic torque thus dominates for a slender rigid fibre. For a slender fibre alone, we also briefly analyse the effects of elasticity on fibre orientation outside the second-order fluid regime.