Wake instability of a fixed spherical droplet with a high drop-to-fluid viscosity ratio

Abstract

Direct numerical simulations of a uniform flow past a fixed spherical droplet are performed to investigate the parameter range within which the axisymmetric flow becomes unstable due to an external flow bifurcation. The hydrodynamics is governed by three dimensionless parameters: the viscosity ratio μ*, the external Reynolds numbers Re^e, and internal Reynolds numbers Reⁱ, respectively. The drop-to-fluid density ratio is related to these parameters as ρ* = μ*Reⁱ/Re^e. This study focuses on highly viscous droplets with μ* ≥ 5, where wake instability is driven by the vorticity flux transferred from the droplet surface into the surrounding fluid. By analysing the wake structure, we confirm that the onset of the external bifurcation is linked to the tilting of the azimuthal vorticity ω_ϕ, in the wake and that the bifurcation occurs once the isocontours of ω_ϕ align nearly perpendicular to the symmetry axis. We propose an empirical criterion for predicting the onset of the external bifurcation, formulated in terms of the maximum vorticity on the external side of the droplet surface. This criterion is applicable for sufficiently high Reⁱ and holds over a wide range of μ* and Re^e. Additionally, we examine the bifurcation sequence for two specific external Reynolds numbers, Re^e = 300 and Re^e = 500, and show that, beyond a critical viscosity ratio, the axisymmetric wake first transitions to a steady planar-symmetric state before undergoing a secondary Hopf bifurcation. Finally, we highlight the influence of Reⁱ on external bifurcation and show that, at moderate Reⁱ, wake instability may set in at a lower vorticity threshold than predicted by our criterion. These findings provide new insights into the external flow bifurcation of viscous droplets.