Bistability and charge-density blowup in the onset of drop Quincke rotation

Particles in a sufficiently strong electric field spontaneously rotate, provided that charge relaxation is slower in the particle than in the suspending fluid. It has long been known that drops also exhibit such “Quincke rotation,” with the electrohydrodynamic flow induced by electrical shear stresses at the interface leading to an increased critical field. However, the hysteretic onset of this instability, observed for sufficiently low-viscosity drops, has so far eluded theoretical understanding—including simulations that have struggled in this regime owing to charge-density steepening driven by surface convection. Here, we conduct a numerical study of the leaky-dielectric model in a simplified two-dimensional setting involving a circular drop, considering arbitrary viscosity ratios and field strengths. As the viscosity of the drop is decreased relative to the suspending fluid, the pitchfork bifurcation marking the onset of drop rotation is found to transition from supercritical to subcritical, giving rise to a field-strength interval of bistability. In this subcritical regime, the critical field is always large enough that, at the bifurcation, the symmetric base-state solution exhibits equatorial charge-density blowup singularities of the type recently described by Peng []. As the rotation speed increases along the initially unstable solution branch from the bifurcation, the singularities gradually shift from the equator and ultimately disperse once the rotational component of the flow is strong enough to eliminate the surface stagnation points.