Flow in a differentially rotated cylindrical drop at moderate Reynolds number

Galerkin finite-element approximations are combined with computer-implemented perturbation methods for tracking families of solutions to calculate the steady axisymmetric flows in a differentially rotated cylindrical drop as a function of Reynolds number Re, drop aspect ratio and the rotation ratio between the two end disks. The flows for Reynolds numbers below 100 are primarily viscous and reasonably described by an asymptotic analysis. When the disks are exactly counter-rotated, multiple steady flows are calculated that bifurcate to higher values of Re from the expected solution with two identical secondary cells stacked symmetrically about the axial midplane. The new flows have two cells of different size and are stable beyond the critical value Rec. The slope of the locus of Rec for drops with aspect ratio up to 3 disagrees with the result for two disks of infinite radius computed assuming the similarity form of the velocity field. Changing the rotation ratio from exact counter-rotation ruptures the junction of the multiple flow fields into two separated flow families.

Flow in a differentially rotated cylindrical drop at low Reynolds number

A liquid drop held captive between parallel disks that are differentially rotated is a model for the swirling flows induced by crystal rotation in the floating-zone process for growing semiconductor materials. An asymptotic analysis for a cylindrical drop is presented that elucidates the structure of the axisymmetric cellular motions caused by disk rotation at low Reynolds number. Variations of meniscus shape induced by these flows are described in the limit of small capillary number. Most cellular flow fields break the bifurcation point that corresponds to the Plateau–Rayleigh limit for the length of a static drop into two disjoint shape families and lower the maximum stable drop length. This effect is studied by a singular bifurcation analysis.