Small deformation theory for two leaky dielectric drops in a uniform electric field

Abstract

A small deformation theory for two non-identical spherical drops freely suspended in an ambient fluid and subjected to a uniform electric field is presented. The three phases are assumed to be leaky dielectric (slightly conducting) viscous incompressible fluids and the nonlinear effects of inertia and surface charge convection are neglected. The deformed shapes of the drops are linearized with respect to the electric capillary number that characterizes the balance between the electric stress and the surface tension. When the two drops are sufficiently far apart, their deformed shapes are predicted by Taylor’s small deformation theory—depending on Taylor’s discriminating function, the drops may become prolate, oblate or remain spherical. When the two drops get closer to each other, in addition to becoming prolate/oblate, they start translating and developing an egg shape. (Since there is no net charge, the centre of mass of the two drops remains stationary.) The extent of each of these ‘modes’ of deformation depends on the distance between the drops’ centres and on drop-to-ambient fluid ratios of electric conductivities, dielectric constants and viscosities. The predictions of the small deformation theory for two drops perfectly agree with the existing results of two-drop dynamics simulation based on a boundary-integral equation approach. Moreover, while previous works observed only three types of behaviour for two identical drops—the drops may either become prolate or oblate and move towards each other or become prolate and move away from each other—the small deformation theory predicts that non-identical drops may, in fact, become oblate and move away from each other when the drop-to-ambient fluid conductivity ratios are reciprocal and the drop-to-ambient fluid viscosity ratios are sufficiently large. The presented theory also readily yields an analytical insight into the interplay among different modes of drop deformation and can be used to guide the selection of the phases’ electromechanical properties for two-drop dynamics simulations.