The orientation dynamics of a massive ellipsoid in simple shear flow

The orientation dynamics of a massive rigid ellipsoid in simple shear flow of a Newtonian fluid is investigated in detail. The term `massive' refers to dominant particle inertia, as characterized by $St \gg 1$, $St = \dot{\gamma}l^2\rho_p/\nu \rho_f$ being the Stokes number; here, $\dot{\gamma}$ is the shear rate, $\nu$ is the kinematic viscosity, $l$ is a characteristic ellipsoidal dimension (taken to be the longest semi axis), and $\rho_p$ and $\rho_f$ are the particle and fluid densities, respectively. Fluid inertial effects are neglected, so the particle Reynolds number ($Re$) is zero. The equations of motion of the ellipsoid in this limit reduce to those governing an Euler top, supplemented by a weak viscous torque. The dynamics consists of a fast conservative motion on time scales of $O(\dot{\gamma}^{-1})$ that, for an ellipsoid, involves a combination of spin, precession and nutation, and a slower component driven by the viscous torque; the latter modulates the angular momentum and rotational kinetic energy on asymptotically longer time scales of $O(St{\dot{\gamma}}^{-1})$. The separation of time scales for large $St$ allows for use of the method of averaging for a triaxial ellipsoid, to derive an autonomous system of ODE's that govern motion on a four-dimensional slow manifold consisting of the three angular momenta and the rotational kinetic energy. There are three fixed points on the slow manifold, with rotation about the shortest axis being stable. The fixed point corresponding to rotation about the longest axis is a saddle point, while intermediate-axis-aligned rotation corresponds to a singular fixed point. For spheroids, the singular fixed point merges with one of the other two, leading to only two fixed points, one a stable node and the other a saddle; the former corresponds to rotation about the shorter axis.