On the Stability of Rotating Drops

We consider the equilibrium and stability of rotating axisymmetric fluid drops by appealing to a variational principle that characterizes the equilibria as stationary states of a functional containing surface energy and rotational energy contributions, augmented by a volume constraint. The linear stability of a drop is determined by solving the eigenvalue problem associated with the second variation of the energy functional. We compute equilibria corresponding to both oblate and prolate shapes, as well as toroidal shapes, and track their evolution with rotation rate. The stability results are obtained for two cases: (i) a prescribed rotational rate of the system ("driven drops"), or (ii) a prescribed angular momentum ("isolated drops"). For families of axisymmetric drops instabilities may occur for either axisymmetric or non-axisymmetric perturbations; the latter correspond to bifurcation points where non-axisymmetric shapes are possible. We employ an angle-arc length formulation of the problem which allows the computation of equilibrium shapes that are not single-valued in spherical coordinates. We are able to illustrate the transition from spheroidal drops with a strong indentation on the rotation axis to toroidal drops that do not extend to the rotation axis. Toroidal drops with a large aspect ratio (major radius to minor radius) are subject to azimuthal instabilities with higher mode numbers that are analogous to the Rayleigh instability of a cylindrical interface. Prolate spheroidal shapes occur if a drop of lower density rotates within a denser medium; these drops appear to be linearly stable. This work is motivated by recent investigations of toroidal tissue clusters that are observed to climb conical obstacles after self-assembly [Nurse et al., Journal of Applied Mechanics 79 (2012) 051013].

The shape and stability of rotating liquid drops

Equilibrium shapes and stability of rotating drops held together by surface tension are found by computer-aided analysis that uses expansions in finiteelement basis functions. Shapes are calculated as extrema of appropriate energies. Stability and relative stability are determined from curvatures of the energy surface in the neighbourhood of the extremum. Families of axisymmetric, two-, three- and four-lobed drop shapes are traced systematically. Bifurcation and turning points are located and the principle of exchange of stabilities is tested. The axisymmetric shapes are stable at low rotation rates but lose stability at the bifurcation to twolobed shapes. Two-lobed drops isolated with constant angular momentum are stable. The results bear on experiments designed to further those of Plateau (1863).