Next Article in Journal
A Flexible, Microfluidic, Dispensing System for Screening Drug Combinations
Next Article in Special Issue
Characterization and Optimization of Isotachophoresis Parameters for Pacific Blue Succinimidyl Ester Dye on a PDMS Microfluidic Chip
Previous Article in Journal
Shrinkage-Considered Mold Design for Improvement of Micro/Nano-Structured Optical Element Performance
Previous Article in Special Issue
Large-Scale Flow in Micro Electrokinetic Turbulent Mixer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Electro-Hydrodynamics of Emulsion Droplets: Physical Insights to Applications

1
School of Mechanical Engineering, Sungkyunkwan University, Suwon 16419, Korea
2
Faculty of Mechanical Engineering, University of Engineering and Technology, Lahore 54890, Pakistan
3
Institute of Quantum Biophysics, Sungkyunkwan University, Suwon 16419, Korea
*
Author to whom correspondence should be addressed.
Micromachines 2020, 11(10), 942; https://doi.org/10.3390/mi11100942
Submission received: 29 September 2020 / Revised: 14 October 2020 / Accepted: 14 October 2020 / Published: 18 October 2020
(This article belongs to the Special Issue Electrokinetics in Micro-/nanofluidic Devices)

Abstract

:
The field of droplet electrohydrodynamics (EHD) emerged with a seminal work of G.I. Taylor in 1966, who presented the so-called leaky dielectric model (LDM) to predict the droplet shapes undergoing distortions under an electric field. Since then, the droplet EHD has evolved in many ways over the next 55 years with numerous intriguing phenomena reported, such as tip and equatorial streaming, Quincke rotation, double droplet breakup modes, particle assemblies at the emulsion interface, and many more. These phenomena have a potential of vast applications in different areas of science and technology. This paper presents a review of prominent droplet EHD studies pertaining to the essential physical insight of various EHD phenomena. Here, we discuss the dynamics of a single-phase emulsion droplet under weak and strong electric fields. Moreover, the effect of the presence of particles and surfactants at the emulsion interface is covered in detail. Furthermore, the EHD of multi-phase double emulsion droplet is included. We focus on features such as deformation, instabilities, and breakups under varying electrical and physical properties. At the end of the review, we also discuss the potential applications of droplet EHD and various challenges with their future perspectives.

Graphical Abstract

1. Introduction

In recent years, electrohydrodynamics (EHD) of emulsion droplets has evolved in many ways as the research gained pace in this field. The field emerged with the seminal works of Taylor in 1966, who presented his leaky dielectric model (LDM) and predicted the single-phase droplet steady deformation [1,2]. Subsequently, the validation of the model was reported, and more accurate models emerged [3,4,5,6,7,8,9,10,11,12]. Recent EHD studies have explored more complex EHD interface topologies related to multiple phase emulsion droplets, and also covered the broader aspects, such as the emulsion instabilities, breakups, and particles manipulation at the emulsion interface forming novel colloidal assemblies [13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34]. Though it has been almost six decades since the research in this area commenced, it is only until recently that a few studies pertaining to the droplet EHD with applicative prospects have been reported [35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52].
This review has emphasis on the recent advancements for the essential physical insights of the droplet EHD problems. Since droplet EHD is a very broad area of research and it is not possible to cover all the aspects in a single review, here we mainly focus on the behavior of a single-phase droplet with particle-free or particle-covered interface and multi-phase (compound or double) emulsion droplets within the perspective of their deformation, instabilities, and breakups under weak and strong external stimuli of electric fields, considering only the most prominent works (See Figure 1). Here, we do not discuss the EHD of droplet coalescence, as it merits a separate review. In the latter part of the article, we also discuss the applications of droplet EHD with their challenges and future directions.

2. Approaches Used to Study the Droplet Electrohydrodynamics (EHD)

The approaches to study the fundamentals of droplet EHD are twofold: (i) theory and/or numerical method, and (ii) experiments. The basic analytical model is the LDM [1,2,53]. This model assumes no net charge on the interface and electroneutrality of the bulk liquids within the Stokes regime, however, the derivation of the LDM from the electro kinetic equations considering the net charge (diffuse charge or electric double layers) is still an ongoing struggle [10,54,55,56]. For numerical simulation, the governing equations of fluid flow and electro-quasi statics are implemented, and the two-phase problem is solved using different methods, such as volume of fluid (VOF) [57,58], boundary element method (BEM) [26,59,60,61], finite volume technique [62,63], level set methods (LS) [14,64,65,66], and phase field methods (PF) [67,68]. Each method has its advantages and disadvantages. Recently, CLSVOF (coupled level set and volume of fluid) method has emerged as a new approach for the interface capturing [69,70]. VOF method has its drawbacks to calculate its spatial derivatives across the interface. Then, a coupled LS and VOF approach is employed to overcome this deficiency. With its spatial gradient calculated accurately, the level-set function is smooth and continuous. For the experiments, a container (cuvette or glass-walled rectangular box) with electrodes embedded on the two opposite sides is filled with a continuous liquid and a dispersed liquid is dispensed using a micro-pipette. The high voltage is supplied across the electrodes and dynamics are captured by using cameras. In these studies, some authors have also introduced surfactants at the liquid-liquid interface [27,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89]. A summary of the most prominent works related to droplet EHD are summarized in Table 1.

3. Background and Theory

It was a general perception before the seminal work by Taylor that the liquid droplets submerged in another immiscible liquid behave as a perfect dielectric and that their exposure to an electric field will always result in prolate deformation (elongation along the direction of the electric field) of the interface [96]. This concept was presented as the electro-hydrostatic (EHS) theory. This theory always predicted that the electric stresses act normal to the interface. However, later it was found by the experiments of Allan and Mason that all the dielectrics do not always exhibit elongation along the electric field direction and contrary to EHS theory, some show oblate deformation (elongation perpendicular to the direction of the applied electric field) [97]. Motivated by these experimental findings, Taylor in his pioneer work introduced the LDM, assuming that liquids with finite conductivities permit slight current flow through them and result in accumulation of free charge at the liquid-liquid interface [1,2]. The interaction of this free charge with the electric field results in the creation of net tangential electric stresses in addition to the normal stresses. Since hydrodynamic stresses should also exist to balance the tangential electric stresses, Taylor concluded the existence of flow circulations, which were later confirmed experimentally. Taylor’s theory predicts the small deformation of emulsion droplets under Stokes flow reasonably well. Later, Torza et al. [7] studied the response of droplets both experimentally and theoretically under DC and AC electric fields and found that Taylor’s theory underestimates the deformation for most of the cases, though the sense of the deformation observed in experiments agreed well with Taylor’s predictions. Ajayi [3] improved Taylor’s theory and calculated the droplet deformation up to a second-order approximation. His theory agreed with the experimental results up to a larger limit of droplet deformation compared to Taylor’s theory. In addition to the steady-state deformation, Esmaeeli and Sharifi [5] predicted the transient deformation of droplet.

3.1. Non-Dimensional Parameters

For a droplet of radius, a, under an unperturbed external electric field strength E o , the characteristic velocity U is defined as:
U = ε 2 E o 2 a / μ 2
The droplet EHD problem depends on the non-dimensional parameters based on the liquid properties, such as conductivity ratio ( R = σ 1 / σ 2 ), permittivity ratio ( S = ε 1 / ε 2 ), and viscosity ratio ( M = μ 1 / μ 2 ) between the inside and the outside liquids. Here, the subscripts ‘1′ and ‘2′ denote the droplet and outside liquid, respectively. Non-dimensionalization of the governing equations results in other various non-dimensional parameters, such as flow Reynolds number ( R e ), electric Reynolds number ( R e E ), and electric capillary number ( C a E ) . The R e is the ratio of inertial to viscous forces and is expressed mathematically as:
R e = ρ 2 a U / μ 2
where ρ 2 is the density of outside liquid. LDM is valid only for low flow Reynolds number R e O ( 1 ) and for small droplet size.
The ratio of the viscous to the inertial force and surface tension is called the Ohnesorge Number (Oh), and controls the transient modes of deformation, such as monotonic and oscillating,
O h = μ / ρ γ a
Electric capillary number is defined as the ratio between the magnitude of electrical stresses ( ε 2 E 0 2 ) and capillary stresses ( γ / a ), where γ is the interfacial tension between two liquids, and controls/delineates the droplet deformation regimes, such as stable and unstable,
C a E = μ 2 U / γ = ε 2 E o 2 a / γ .
Electric Reynolds number is defined as the ratio of the charge relaxation time scale ( t c = ε 2 / σ 2 ) to the time scale of charge convection by flow ( t h = μ 2 / ε 2 E 0 2 ). A small value implies that the charge convection effects are negligible. Based on the medium properties, it is given as:
R e E = t c / t h = ε 2 2 E o 2 / μ 2 σ 2
Another important non-dimensional number is the Saville number (Sa), which is the ratio of the electric capillary number and electric Reynolds number and decides the droplet breakup mode transitions,
S a = R e E / C a E
LDM assumes the interface is instantly charged, i.e., the surface charge convection effects represented by electric Reynolds number ( R e E ) are negligible.

3.2. Govening Equations of Leaky Dielectric Model (LDM)

Here, we summarize the core governing equation of LDM. Note that the equations reproduced here are for a single-phase emulsion droplet.

3.2.1. Electric Field Equations and Their Solutions

For a general dynamic system, electric field ( E ) and magnetic field ( B ) are coupled in the governing Maxwell equations [98]. The Faraday’s law is given as:
B t + × E = 0 .
Saville [53] showed that with small electric currents and in the absence of external magnetic field, Maxwell equations can be simplified significantly, and the Faraday’s law becomes:
× E = 0 .
Equation (8) shows that the electric field is irrotational and thus, E = V , where V is the electric potential and the droplet is assumed to be placed in an unperturbed external electric field strength of magnitude E o .
Gauss’s law is given as:
Q = · ( ε E )
where Q is the volumetric charge density and ε is the liquid permittivity.
Conservation of charge is another important governing equation given as:
Q t + · I ρ = 0 ,
where I ρ is the current density defined as:
I ρ = σ E + Q u
The above equation shows that the charge movement is owing to the Ohmic conduction and the convection, respectively. Here, u is the liquid velocity and for incompressible liquids with uniform properties in each of the phases; Equations (9)–(11) are combined to give
D Q D t + σ ε Q = 0 ,
where D / D t = / t + u · is the material derivative. The solution of the above equation yields:
Q = Q o exp ( t / t c ) .
The above equation shows that the volumetric charge decays by time exponentially and becomes zero in the bulk. Where, t c = ε / σ is the charge relaxation time scale or the rate of charge decay. Provided that if t c is much smaller than the time scale of any other process of interest, then the charge immediately accumulates at the interface and Equation (11) becomes:
I ρ = σ E , · ( σ E ) = 0 .
and for liquids with constant properties:
· E = 0 .
As the electric field is irrotational and E = V , this equation results in:
2 V = 0 ,
where 2 is the Laplacian operator? For a single droplet with azimuthal symmetry, the electric potential V is only a function of r and θ in the spherical coordinates. For a droplet of radius a under an unperturbed external electric field strength E o , the solution of Equation (16) with boundary conditions [1,7],
V 1 ( 0 , θ )   should   be   bounded , V 1 ( a , θ ) = V 2 ( a , θ ) , σ 1 V 1 / r = σ 2 V 2 / r , V 2 ( , θ ) = E o r cos θ ,
leads to the electric potential at the droplet inside V 1 and droplet outside V 2 given as:
V 1 = 3 R + 2 ( r a ) cos θ   a E o ,
V 2 = [ ( r a ) R 1 R + 2 ( a r ) 2 ] cos θ   a E o .
The electric field at the inside and the outside of the droplet can be determined by taking the gradient of the electric potential given in Equations (18) and (19), respectively.
The Maxwell stress tensor is given as [99]:
T M = ε E E 1 2 ε E · E I .
The divergence of Maxwell stress tensor determines the electric force per unit volume of the liquid ( f E ) that is substituted in Navier–Stokes equation. Conveniently, electric stresses can be split into normal and tangential components in a t n coordinate system and their jumps across the interface have a direct impact on the flow field and deformation of the droplet at the steady state, respectively. The jumps in the tangential T M , t 12 and normal components T M , n 12 of electric stresses can be found by the evaluation of the stress tensor at both sides of the droplet interface where E n E r and E t E θ , where 12 denotes the jump of the property, “ 2 1 ” across the interface.
T M , θ 12 = 9 2 S R ( R + 2 ) 2 ε 2 E o 2 sin 2 θ ,
T M , r 12 = 9 2 ( R 2 + 1 2 S ) cos 2 θ + S 1 ( R + 2 ) 2 ε 2 E o 2 .
From these equations, one can deduce the direction of shear stresses that are particularly important in determining the flow circulations. When R > S , the direction of net shear stress is from the equator to the poles and vice versa.

3.2.2. Fluid Flow Equations and Their Solutions

The governing equation for fluid motion with velocity ( u ) , viscosity ( μ ), and pressure ( p ) is the conservation of mass and momentum.
· u = 0 ,
μ 2 u + f E = p .
The divergence of the Maxwell stress tensor from Equation (20) gives:
f E = · T M = Q E 1 2 E 2 ε .
As discussed above, if the charge instantly accumulates at the interface, for example, for a leaky dielectric system, the volumetric charge density ( Q ) is zero in the bulk and for liquids with constant properties, ε = 0 . Thus, f E = 0 in the bulk and Equation (24) reads
μ 2 u = p .
This shows that the electric force affects the solution only through the momentum jump boundary condition at the interface. The solution of the resulting Stokes equation can be obtained by introducing a stream function at the two sides of the droplet interface [100].
u r = 1 r 2 sin θ ψ θ ,
u θ = 1 r sin θ ψ r ,
D 4 ψ = 0 .
D 4 = D 2 ( D 2 ) , where D 2 = 2 / r 2 + ( sin θ / r 2 ) ( / θ ) [ ( 1 / sin θ ) ( / θ ) ] is a well-known differential operator appearing in Stokes equation. The general solution of ψ is made specific by the following boundary conditions:
u r , 1   and   u θ , 1   should   be   bounded   at   r = 0 , u θ , 1 ( a , θ ) = u θ , 2 ( a , θ ) , u r , 1 ( a , θ ) = u r , 2 ( a , θ ) = 0 , T H , θ 12 + T M , θ 12 = 0 , u r , 2 ( , θ ) = u θ , 2 ( , θ ) = 0 .
The solution using the above boundary conditions at a steady state gives:
ψ 1 = [ ( r a ) 3 ( r a ) 5 ] U m a x a 2 sin 2 θ cos θ ,
ψ 2 = [ ( a r ) 2 1 ] U m a x a 2 sin 2 θ cos θ ,
U m a x = 9 ε 2 E o 2 a ( S R ) 10 μ 2 ( 1 + M ) ( R + 2 ) 2
The velocity field is obtained from Equations (27) and (28).
u r , 1 = [ ( r a ) 3 ( r a ) ] U m a x ( 1 3 cos 2 θ ) ,
u θ , 1 = 1 2 [ 5 ( r a ) 3 3 ( r a ) ] U max sin 2 θ ,
u r , 2 = [ ( a r ) 2 ( a r ) 4 ] U m a x ( 1 3 cos 2 θ ) ,
u θ , 2 = ( r a ) 4 U max sin 2 θ .
From these equations, it is evident that the flow direction is from the poles to the equator if U m a x is positive and vice versa. On the other hand, the sign of U m a x depends on the sign of quantity ( S R ) .

3.2.3. Taylor Deformation Parameter ( D )

The deformation or distortion of the droplet interface is caused by the electric and the hydrodynamic fields being discontinuous at the interface, resulting in the interfacial stress jumps countered by the interfacial capillary stresses.
[ σ H , r 12 + T M , r 12 ] = γ κ ,
where σ H , r are the normal hydrodynamic stresses ( σ H , r 12 = T H , r 12 p 12 ). T H , r 12 is found from the velocity field and p 12 is obtained by integrating Equation (26). The coupling between hydrodynamic and electric fields occur at the interface.
For small deformations:
κ = 2 a + 8 D 3 a ( 3 cos 2 θ 1 ) ,
where the Taylor’s deformation parameter D is defined as:
D = ( L W ) / ( L + W ) ,
where L and W represent the length of the axis of the deformed droplet parallel and perpendicular to the electric field direction, respectively. Additionally:
σ H , r 12 = ( 2 + 3 M ) ( 1 3 cos 2 θ ) μ 2 U m a x a .
Substitution from Equations (22), (39) and (41) in Equation (38) gives a deformation of droplet in a uniform electric field:
D D C =   9 16 C a E f d , f d = 1 ( R + 2 ) 2 [ ( R 2 + 1 ) 2 S + 3 ( R S ) ( 3 M + 2 5 M + 5 ) ]
where C a E is the electric capillary number and f d is the deformation characteristic function When   f d > 0 , the droplet shows prolate deformation and, when f d < 0 , it shows oblate deformation. A leaky dielectric droplet remains spherical if f d = 0 (See Figure 2A). The deformation-circulation map depends on R and S of the system. If S / R < 1 , the flow direction is from equator to pole and resulting shape of the droplet is always prolate (Region I). On the other hand, if S / R > 1 , the flow direction is from the pole to equator and droplet shape can be prolate (Region II) or oblate (Region III) as shown in Figure 2B.
Ajayi [3] used a domain perturbation method to improve the deformation predicted from Equation (42) by incorporating second-order function of C a E , and the formulations are given as [59],
D D C = k 1 C a E + k 2 C a E 2 + O ( C a E 3 ) ,
where
k 1 = ( 9 / 16 ) f d , k 2 = k 1 [ ( 9 5 R 1 R + 2 1 16 ) f d + R S ( R + 2 ) 2 β ] , β = 23 20 139 210 ( 1 M 1 + M ) 27 700 ( 1 M 1 + M ) 2 .
The transient deformation, D D C ( t )   of the droplet has an exponential form given as [5]:
D D C ( t ) = D [ 1 exp ( t / τ ) ] , τ = ( 19 M + 16 ) ( 2 M + 3 ) a μ 2 ( 40 M + 4 ) γ .
Here, D is the steady state deformation as calculated from Equations (42) or (43).

3.2.4. Effect of AC Electric Field

The theoretical foundation was established by Torza et al. [7] who considered a droplet under an AC electric field, E e x t ( t ) = E o cos ω t , within the Stokes flow regime, where ω is the frequency of applied electric field. The net normal T M , r 12 and tangential electric stresses T M , θ 12 , the velocity field   u , and the deformation D have a time-independent steady component ( s ) that is a function of frequency in addition to another important parameter and a time-dependent component ( t ) with frequency twice that of frequency of an applied electric field.
A C = s + t .
Here represents one of T M , r 12 , T M , θ 12 , u , and D . For the steady component of the stress, the deformation is given by an equation similar to (42), namely:
D s =   9 16 C a E f s , f s = 1 R ( 11 M + 14 ) + 15 ( M + 1 ) + S ( 19 M + 16 ) + 15 Π 2 ( M + 1 ) ( 2 S + 1 ) 5 ( M + 1 ) [ ( R + 2 ) 2 + Π 2 ( S + 2 ) 2 ] ,
where f s is the deformation characteristic function for an AC electric field, which determines the sense of drop deformation, and Π = ε 2 ω / σ 2 t c / t e , where t c = ε 2 / σ 2 and t e = 1 / ω are the charge relaxation and the electric field time scales. For ω 0 , f d = f s and the steady deformation ( D s )   are same as the deformation under a DC electric field, provided that
E r m s = lim T 1 2 T T T ( E o cos ω t ) 2 d t = E o 2   ,
in place of E o in the Taylor’s solution ( D s = ( 1 / 2 ) D D C ). For ω , the steady deformation ( D s ) is same as the solution provided by Allan and Mason [97] ( D D C _ P D = ( 9 / 16 ) C a E ( S 1 ) 2 / ( S + 2 ) 2 ) for a perfect dielectric in a DC electric field, provided that E r m s = E o / 2 is used ( D s = ( 1 / 2 ) D D C _ P D ). In this case, the charge relaxation time scale ( t c ) is much larger than the timescale of the electric field ( t e ) and hence the system behaves as a perfect dielectric. The transition between prolate and oblate ellipse takes place at a critical frequency defined as
ω c r = σ 2 ε 2 5 ( M + 1 ) ( R 2 + 1 ) + ( 9 M + 6 ) R + ( 16 + 19 M ) S   5 ( M + 1 ) ( S 1 ) 2   .
For ω > ω c r , D s > 0 (prolate deformation), and for ω < ω c r , D s > 0 , or D s < 0 (oblate deformation) (See Figure 2C). At ω = ω c r , the drop remains spherical ( D s = 0 ). For the time-dependent part,
D t = 9 32 C a E [ Φ cos ( 2 ω t + α ) ( 1 + k 2 λ 2 2 ) 1 / 2 ] ,
where
Φ { f s 2 + Π 2 ( 19 M + 16 ) ( R S ) 2 [ ( M + 1 ) ( 20 S 1 ) 3 ] 25 ( M + 1 ) 2 [ R 2 ( 1 + 2 R ) 2 + Π 2 ( S + 2 ) 2 ] 2 } 1 / 2 , k ω μ 2 a / γ , cos α h * + h * ¯ 2 I , sin α h * h * ¯ 2 I 𝕚 , h * F r * λ 1 F θ * 1 + 𝕚 k λ 2 , F r * A 2 ( 5 2 S ) 2 A + 1 , F θ * R ( 1 + 𝕚 Π ) [ R + 2 + 𝕚 Π ( S + 2 ) ] 2 , λ 1 3 ( 3 M + 2 ) 5 ( M + 1 ) , λ 2 ( 19 M + 16 ) ( 2 M + 3 ) 20 ( M + 1 ) , A 1 + 𝕚 Π 2 + R + 𝕚 Π ( S + 2 ) ,
where 𝕚 is imaginary unit, and h * ¯ represents the complex conjugate of h * .
To evaluate how much oscillatory deformation ( D t ) affects the overall deformation compared to the steady deformation ( D s ), we consider the following ratio:
D t D s = Φ cos ( 2 ω t + α ) 2 f s ( 1 + k 2 λ 2 2 ) 1 / 2 ,
Obtained from Equations (47) and (50). The dimensionless parameter k , defined by Equation (51) and which plays an important part in characterizing the ratio ( D t / D s ), may be considered as the ratio of the oscillatory hydrodynamic stress ( ω μ 2 ) and the capillary pressure ( γ / a ) at the droplet interface. When ω μ 2 and γ / a remains fixed, k and D t / D s , the droplet surface being unable to respond to the oscillating stress.

4. Single-phase Emulsion Droplets

4.1. Deformation

As discussed in Section 3, the droplet can undergo prolate or oblate steady deformation. The EHD of single-phase emulsion droplets can also be categorized with respect to electric capillary number ( C a E ) and steady or transient deformation (D). Taylor’s LDM predicts the deformation of droplets within low-deformation limit, i.e., when C a E 1 adequately well. In this case, the droplet adopts a spheroidal shape. If the electric field strength is increased, the droplet is susceptible to instabilities and breakups, and LDM no longer remains valid. The steady and transient deformation of droplets at a low electric field is vastly studied and is well understood. In literature, various analytical, numerical, and experimental studies addressing different issues, such as effect of inertia, charge convection, charge relaxation, and type of electric field (AC or DC, uniform or Non-uniform), are found [4,5,6,8,9,11,24,54,61,68,102,103,104,105,106,107,108,109,110,111].

4.1.1. Under DC Electric Field

The transient deformation of a droplet essentially depends on two dimensionless numbers, Ohnesorge number (Oh) and the Reynolds number (Re) [68]. When R e 1 or O h > 1 , the droplet always has a steady mode represented by monotonic deformation during its transient evolution. Whereas, if R e 1 and O h < 1 , the droplet has an oscillating deformation before achieving the steady configuration. The droplet has a quasi-steady mode of deformation at the borderline region that distinguishes the monotonic and oscillating deformation regions (See Figure 3A,B). Feng and Scott showed that the droplet shape may change from oblate to prolate due to the inertial effects associated with the increase in the strength of the electric field [6]. The main effect of charge convection is to reduce the interfacial velocity that results in more deformation of a prolate droplet, while an oblate droplet undergoes lesser deformation. Lanauza [103] and Das et al. [4] showed that the surface charge convection causes formation of “charge shocks” near an oblate droplet’s equatorial region. The experimental data for the deformation of a silicone oil droplet in castor oil is as shown in Figure 3C. All the theoretical or numerical models fail to capture the dynamics except the one that considers the charge relaxation and convection effects.

4.1.2. Under AC Electric Field

Although the sense of steady-state flow circulation pattern in an AC electric field is the same as that in a DC electric field, the sense of D s is not the same as that of D D C . Instead, it depends on the field frequency ω being below or above a critical frequency ω c , defined by Equation (49) [7,101,112]. Esmaeeli and Halim [101] compared the transient deformations of droplet in a DC electric field to those in an AC electric field using direct numerical simulation (DNS). For the range of the physical parameters used in their study, the evolution of the drop deformation D AC with time t showed that the drop settled to its quasi-steady state in a relaxation time t r = ( μ 1 + μ 2 ) a / γ (See Figure 3D). Additionally, they showed that the kinetic energy of the fluid system is minimized when the rate of the deformation ( D A C / t ) is zero, where there is a state of minimum or maximum deformation. On the other hand, it is maximized when in the state of steady deformation during the deformation half-cycle (See Figure 3E). Esmaeeli [112] showed flow field in and outside of drop in an AC electric field, and evolution of the total deformation D A C versus nondimensional time ω t for three different deformation regions.

4.2. Instabilities and Breakups

The dynamics of a droplet under a strong electric field when C a E 1 depends on the ratio of electrical conductivity to electrical permittivity ( R / S ).

4.2.1. Prolate Deforming Droplets

In most of the studies related to single-phase droplets, the droplet generally settles to an equilibrium shape under a relatively weak electric field. However, if the strength of the electric field is high and R / S > 1 , then the droplet becomes unstable and thereby disintegrates or breaks into daughter droplets. This topic is covered by different authors who have tried to determine the critical parameters to distinguish between the stable and unstable regimes [6,8,9,15,22,26,59,61,113]. The numerical models are developed assuming ellipsoidal or slender shapes. Sherwood [61] studied the droplets under a strong electric field using boundary integral method. Similar to Sherwood, Lac and Homsy [59] used an axisymmetric boundary integral method and studied droplet stability by changing different parameters. Bentenitis and Krause [8] theoretically studied the large deformation of the droplets, extending the conventional leaky dielectric model by analytical solutions of an electric field and flow field in a framework of spheroidal deformation. Figure 4A shows a transient evolution of a droplet that continuously stretches under the applied electric field and eventually becomes unstable, ejecting daughter droplets at the side apex. It was reported that as the droplet’s semi-major axis becomes greater than 1.5 times its undeformed radius, the droplet becomes unstable [15]. For the case of a castor oil droplet in silicone oil, this instability occurred at C a E ~ 0.25 (See Figure 4B). However, the critical electric capillary number for a conducting droplet, such as a water droplet in silicone oil, is well known to take place at C a E ~ 0.2 ± 0.2 [6,61].
Two important modes of breakups are reported, such as end pinching and tip streaming [114,115]. In the former mode, the droplet is converted into bulbous-shaped lobes that eventually disintegrate, whereas in the latter case, the droplet develops sharp cone-like tips ejecting jets. These jets eventually disintegrate into tiny droplets. The transition between these modes is rather complex and depends on various parameters, such as droplet conductivity ( R ), viscosity ratio ( M ) between the liquids, and applied electric field ( E 0 ). Figure 4C shows the effect of viscosity ratio and droplet conductivity. At M = 1 , the droplet forms small daughter droplets at the tip with small protrusion. As the M increases, the jet length and daughter droplet size increases. This shows the importance of viscous stresses to balance the electrical tangential stresses to form a jet. Droplet conductivity effect is represented in terms of S a ( S a 1 / σ 2 ). A high-conductivity droplet ejects thinner jets and smaller daughter droplets via tip streaming. For a low-conducting droplet, the jet vanishes, and the breakup is via end-pinching due to rapid domination of capillary stresses. It is also notable that with a decrease in droplet conductivity, the jet length first increases, reaches a maximum ( S a 1 ) , and then it decreases [114]. It is also reported that a change of mode can also take place due to surface charge convection effects represented by R e E , as shown in Figure 4D [115].

4.2.2. Oblate Deforming Droplets

The droplet undergoes electro-rotation under a strong electric field if R / S < 1 due to the surface charge convection effects (See Figure 5A) [18,90,116,117,118,119]. The phenomenon is similar to the rotation of solid particles in fluid known as Quincke rotation, where the induced dipole tries to align itself with the electric field, and the resulting torque causes a rotation perpendicular to the direction of the applied electric field. It takes place when the applied electric field is above the certain threshold value ( E q ). The value for a rigid sphere is given as [118,120]:
E q = 2 μ 2 σ 2 ( R + 2 ) 2 3 ε 1 ε 2 ( 1 R / S ) .
The equation is obtained by the balancing of electrical and viscous torques and is independent of the size of sphere. The equation is often used to estimate the threshold of electro-rotation for the liquid droplets and comparison with experimental results shows a relatively good approximation [117]. Results show that the droplet attains a steady tilt relative to the applied electric field and can undergo irregular rotational motions (See Figure 5B). The electro-rotation stabilizes the droplet against breaking by suppressing the droplet deformation. Recently, Brosseau et al. [18] observed very interesting phenomena, such as droplet dimpling and equatorial streaming (Figure 5C,D). In the dimpling mode, the droplet forms a torus shape; the torus subsequently breaks into relatively larger drops. This takes place for droplets with M 1 and C a E ~ O ( 1 ) . In the equatorial streaming mode, the droplet flattens with the sharp edges. Concentric waves in the form of rings emanate from the thin edges that break into tiny droplets. The critical conditions for this mode are R 1 , C a E 1 and M 0.1 .

5. Particle-Covered Droplets

5.1. Deformation

A colloidal particle or Pickering droplet is formed when the particle-laden droplet or suspension, for which R / S < 1 , is subjected to an electric field [21,91,92,121,122,123,124,125,126,127,128,129,130,131,132,133,134,135]. Particles can be trapped at the interface more rapidly under a weak electric field (~ 0.1 kV/mm) and once trapped, they remain there. Depending on the strength of the electric field, particle concentration, electric conductivity, and size, various particle assemblies can be formed at the interface. Mikkelsen et al. [134] used low to high particle coverage (~0.1 to ~0.8), and analyzed the conductive and non-conductive particle conditions, as shown in Figure 6A. For high particle concentration of non-conducting particles, such as polystyrene under a weak electric field, the oblate deformation of silicone oil increased as more charges accumulated at the droplet surface and the EHD flows were suppressed. On the other hand, as the electrical conductivity of particles entrapped at the silicone droplet interface increased, the oblate deformation of the silicone oil droplet changed to prolate deformation. In this case, the dipole moment becomes aligned with the electric field, unlike the case of the silicone droplet with the clean interface or covered with non-conducting particles. For low particle concentrations in a range from ~0.1 to ~0.5 under a weak electric field, the particles at the interface can produce different assemblies, such as “belts” (low-conductivity particles) or “chains” (high-conductivity particles). Belts may form a static or dynamic sinusoid with a change in particle properties and electric field strength [91,134]. High-conductivity particles with regular size, such as spheres, organize into regular chains aligned with the direction of an applied electric field, while random-size particles form random assemblies (See Figure 6B). The particle belts can open and close over the entire surface of the droplet by application and removal of the electric field, respectively [21] (See Figure 6C).

5.2. Instabilities and Breakups

For high particle coverage under a strong electric field, the droplet adopts peculiar drum-like shapes and may even start rotating (electro-rotation phenomenon as discussed in Section 4.1.1) or implode. The particle layer buckles, leading to various outcomes, such as ejection of particle clusters or formation of ephemeral wings [122], as shown in Figure 6D. For low particle concentrations under a strong electric field, droplets with chains (high-conductivity particles) eventually break at the side ends (tip streaming), as shown in Figure 6E. Whereas, those forming belts (low-conductivity particles), either redistribute to form counter-rotating vortices, lose their structure, or eject.

6. Surfactant-Laden Droplets

6.1. Deformation

The influence of interfacial properties on the droplet deformation and breakup under an electric field has also been a topic of profound interest for decades [27,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87]. The surface-active molecules prompt to change in the interracial and rheological properties, leading to interfacial perturbations and breakups [23,71,136,137,138]. It is important to note that under a weak electric field, the steady deformation for a surfactant-laden aqueous droplet is lower at the same C a E (calculated in terms of surfactant-free droplets), and this difference is more prominent for larger deformation [139,140] (See Figure 7A). This could be the result of the surfactant harassing the EHD flows, though the already weak flows for conducting droplets cannot fully justify it. The possible intriguing effect is the surfactant dilution with an increase in interfacial area. The larger the deformation, larger is the increased interfacial area, leading to an increased dilution and hindering the deformation significantly. The transient deformation showed an interesting result that the droplet-stretching velocity decreased with an increase in the non-dimensional surfactant concentration ( C * ), however, as the surfactant concentration increased beyond critical micelle concentration ( C * > 1 ), there was no further decrease [27] (See Figure 7B). The understanding of this dynamics is rather complex. The decrease in velocity for ( C * < 1 ) is attributed to the re-distribution of surfactant molecules with an increase in interfacial area and non-availability of surfactant molecules from the bulk causing local surface tension gradients and Marangoni flows. These flows oppose the droplet stretching and hence the droplet stretching velocity decrease. For ( C * > 1 ), the rate of exchange of the surfactant molecules from the bulk to the interface is still slow and saturation in terms of Marangoni effect occurs. It is noteworthy that the exchange of surfactant molecules in case of an electric field not only depends on viscosity ratio ( M ), but also on the permittivity ( S ) and conductivity ( R ) of the liquids, as the circulation inside and outside the droplet relies on all these parameters and can affect the droplet transient deformation.

6.2. Instabilities and Breakups

Unlike the surfactant-free droplet, where the breakup modes are controlled by ratios of viscosity, conductivity, and permittivity, dynamics of the interface manipulated by surfactant is another important criterion in the breakup of a surfactant-laden droplet. The breakup of a surfactant-laden aqueous droplet under a strong electric field can show different modes such as: conical shape-conical jetting, ellipsoidal shape-conical jetting, ellipsoidal shape-filamentous breakup, and cylindrical shape-filamentous breakup [83] (Figure 7C). The conversion criterion between these breakup modes with respect to electric capillary number and the increase in the surfactant concentration is specified in Figure 7D. It has been shown that this behavior is similar to the tensile-fracture process in the elastic–plastic mechanics and thus the stress-strain relationship is used to analyze the axial stretching.

7. Multi-Phase Emulsion Droplets

In the past few years, multi-phase droplets, such as double-emulsion droplets or compound droplets with three liquid components (core, shell, and ambient) and two interfaces (core/shell and shell/ambient), have gathered special attention due to their vast applications in various scientific areas, such as in targeted drug delivery [141], artificial supply of oxygen in blood [142], etc. The foundation of the theory of a double-emulsion droplet under an electric field is that of a single-emulsion droplet as discussed in Section 3. A more comprehensive coverage of the governing equations is not included in this review and can be found in literature [143,144]. An A/B/A type (core/shell/ambient) of a double-emulsion droplet has the same liquid in the core, and ambient, and is essentially a two-phase problem. On the other hand, an A/B/C type has three distinct liquid components. In this section, we first consider the dynamics of double-emulsion droplets or compound droplets within a weak electric field ( C a E 1 ) and, lately, under a strong electric field.

7.1. Deformation

The dynamics of a double-emulsion droplet within a weak electric field or small deformation limit is extensively reported and well understood [33,67,93,94,143,144,145,146,147,148,149,150,151]. A/B/C types of a double-emulsion droplets can show four distinct modes of deformation of core/shell-shell/ambient, i.e., prolate-prolate (PP), prolate-oblate (PO), oblate-prolate (OP), and oblate-oblate (OO), which entirely depends on the ratio ( R / S ) across each of the interface [67,93,143]. If R / S < 1 , the interface deforms oblate and if R / S > 1 , it deforms prolate. For an A/B/A type, only two modes of core-shell deformation, i.e., PO and OP, are possible. The use of an electric field in the centering of the core droplet can be found extensively in literature and is well understood [42,145,148,149,150]. Oguz and Sadhal [145] conducted the earliest study on the effects of a weak electric field on the dynamics of compound droplets. They showed that under an electric field, stable equilibrium configuration of a compound drop is possible, which may not be possible without using an electric field. T. Sukada et al. [147] showed that deformation and flow strength are positively related to the strength of an electric field and core-to-shell droplet volume ratio, using theoretical analysis and experimental observations.
Behjatian and Esmaeeli [143] studied the A/B/C type concentric emulsion analytically under a DC electric field using domain perturbation method and showed various free-charge distributions surrounding the interfaces depending on the relative magnitudes of R and S (See Figure 8A). They also studied the deformation-time history of the double droplets under weak electric fields [93]. Their analysis was based on the Laplace equation of charge distribution and considering the droplets to have no net charge. Same as in LDM, electro-neutrality was assumed for the bulk liquids, and electric and hydrodynamic tractions became discontinuous at the two interfaces. For R i j < S i j , the charge sense on the upper half of the core/shell and shell/ambient interface are the same with respect to each other and as that on the supply electrode (+), while the lower half has the negative charge (−). The oblate deformation of the interfaces and the OO configuration is expected. For R i j > S i j , the charges are reversed. Here, i j = 12 . for core/shell and i j = 23 for shell/ambient, respectively. The prolate deformation of the interfaces and the PP configuration is expected. For OP and PO configurations, the core/shell and shell/ambient interfaces have opposite charges with respect to each other. By changing the direction of electric field, the charge sense can be reversed but the direction of the tangential electric stress remains the same. Unlike the case for a single-phase emulsion droplet, it is possible to get two counter-rotating vertices in a multi-phase emulsion droplet, as shown in Figure 8B. The generation of toroidal vortices depends on the electrical properties of the constituting liquids; however, the precise criterion has still not yet been established. The extent of the deformation of the inner core droplet also depends on the electrical conductivity of the shell. Soni et al. [144] extended the analytical analysis to an AC electric field. They showed that in the limit of zero frequency, their results reduced to those derived by Behjatian and Esmaelli [143]. Abbasi et al. [14] showed that for shells comprised of fewer conductivity liquids and R 23 1 R 12 , the electric potential lines are densely concentrated in the shell region and reach the core/shell interface easily, causing deformation of both the core and the shell droplets. On the other hand, if shell liquid has large electric conductivity and R 23 1 R 12 , it behaves like a shield and electric potential lines cannot penetrate through it, resulting in a negligible core deformation (See Figure 8C). While Behjatian and Esmaeeli reported a steady droplet configuration under a weak electric field, it was later reported in an experimental study [13] that owing to electrophoretic effects and the presence of net charge arising from the existence of electric double/diffuse charge layers on the droplet, the steady configuration was not possible for some of the cases. As the charged core droplet experiences the electric field, it moves towards the droplet side apex and is eventually dispensed into the ambient.

7.2. Instabilities and Breakups

This feature of multi-phase droplets is relatively less explored. Only a few studies have explored the breakups within a certain range of electric and hydrodynamic properties [13,14,16,25,30,95,152]. In one of the earliest studies, Ha and Yang established the breakups of conducting emulsion droplets comprised of inhomogeneous inner phase encapsulated by an outer membrane phase suspended in silicone oil [25]. Castor oil and aqueous phase were used to form oil-in-water (o/w) or water-in-oil (w/o) emulsions, depending on their fraction of volume, comprising either the membrane or non-uniform inner droplet phase. They recognized that the existence of non-uniformity in the form of multiple inner droplets can significantly alter the breakup modes and stability as compared to homogenous Newtonian conducting droplets. Pinch off and tip streaming modes of emulsion-breaking were observed. Soni et al. [95] explored the breaking of a A/B/C type droplet comprised of NaCl aqueous/castor oil/silicone oil in an AC electric field exhibiting PP mode of deformation, as shown in Figure 9A. The breaking of the core occurred at a larger strength of electric field compared to a NaCl aqueous droplet in a bulk castor oil. This was due to less penetration of electric potential as the shell was acting as a shield. Recently, Abbasi et al. explored the dynamics of a double-emulsion droplet with a low electrical conductivity shell beyond the small deformation limit in two separate studies and explored various breakup modes of a core droplet under varying physical and electrical parameters, such as radius ratio, electric capillary number, and/or interfacial tension [13,16]. In the former, breakups of a A/B/A type of double-emulsion comprised of castor oil/silicone oil/castor oil (core/shell/ambient) for which PO mode of deformation occurs is thoroughly explored. Four distinct breakups were reported, such as uni-directional and three different bi-directional breakups. The charge convection effects caused a change in the breakup of the core from the bulbous-shaped lobes to the formation of conical ends. The breakup modes were then plotted on a core-to-shell droplet radius ratio vs. electric capillary number parametric space. In the later study, a surfactant was introduced at the aqueous core phase of an A/B/C type of double-emulsion to change the core/shell interfacial tension [16]. The core droplet instability lead to the formation of asymmetric Janus and snowman-like stable or unstable ternary droplets, as shown in Figure 9B. By theoretical analysis, the length of stretched semi-major axis of the core at which the instability occurred was established to be 1.5 times its initial undeformed radius, and breakup modes were delineated on a radius ratio vs. electric capillary number parametric space at different non-dimensional surfactant concentrations (See Figure 9C). It is pertinent to mention here that the dynamics of a double-emulsion droplet exhibiting OO mode of deformation under a strong electric field still remains unexplored.

8. Applications

8.1. Mass Production of Micro/Nano Droplets

EHD-based electric shear stress is used to produce tiny droplets by cutting off the liquid interface formed under the nozzle, which is called electrospray [153]. Also, in droplet-based microfluidic systems, electric fields are used to generate smaller droplets [154]. In addition to generating simple tiny single droplets or particles, these techniques can be used to make more complex droplets, such as double droplets [141,155], and thus many studies have been conducted in the fields of bioengineering [156], medical [157], and pharmaceutical [141]. Similarly, breakup at high voltage was observed in free-emulsion droplets [15] or emulsion droplets anchored on a nozzle [158,159] (Figure 10A,B). The droplet that is stretched under a strong electric field produces tiny daughter droplets at both ends, which is called tip streaming [15,25]. Using this phenomenon, it is possible to generate many small droplets ranging from sub-nano to micro size. Moreover, under special conditions, the droplets show equatorial streaming, which can be used to create a large number of micro droplets [18] (Figure 10C). In a more complex case, the condition in which a compound droplet breaks up was analyzed by Su et al., through a numerical analysis method [160].

8.2. Functional Droplet/Particle Synthesis

The electric field can be used to produce and manipulate various types of droplets that are not spherical. Brosseau et al. used cross-linking UV curable polymer (NOA 81, Norland Optical Adhesive, Norland Products Inc., East Windsor, NJ, USA) drop to make spheroidal shapes via EHD deformation [164]. On the other hand, in the case of droplets covered with high-concentration particles, the unique shapes deformed by the EHD flow can be fabricated due to particle jamming even if the electric field is removed [121,128,130] (Figure 10D). Manipulating particles anchored on the droplet interface is one of the promised techniques for designing functional droplet. Using active structuring of a particle-covered droplet under the electric field, Rozynek et al. demonstrated a millimeter-sized optical diaphragm [162] (Figure 10E). Li et al. redistributed aluminum nanoparticles on an o/w emulsion droplet via an electric field for generating a Janus droplet [163] (Figure 10F). Moreover, an AC electric field stabilized the double-emulsion droplet by centering the core droplet, which can be used to synthesize a polymer shell or hollow particle [42].

8.3. Micro-Rotor

When R / S < 1 , a droplet or particle shows a unique phenomenon, called Quincke rotation, under a strong electric field. Particles on the electrode show a rolling or hovering motion with Quincke rotation [165,166]. Many theoretical or experimental researches have been conducted for Quincke rotation of polymer particles to make a rotor [165,167], motor [168], pump [169], and so on. Similarly, a single droplet with/without particles also shows electro-rotation with ellipsoidal deformation (Figure 11A). Electrohydrodynamic propulsion of two Pickering droplets placed at a close proximately takes place at different angles relative to an applied electric field with the propulsion velocities approaching 0.1 mm/s [92]. The electric field is below the threshold for the Quincke rotation and the drops may co-rotate or counter-rotate perpendicular to the field (See Figure 11B).

8.4. Encapsulation/Target Delivery

The double emulsion type capsulation technique is an important technique in the process of drug delivery or producing functional particles, but the resulting double emulsion may be too unstable or too stable for these purposes. However, if an electric field is used, the core droplet can be moved to the center [94,148] or ejected outward [13,14]. Tucker–Schwartz et al. showed that in the process of encapsulation of a material having a density difference with a photocurable polymer, centering through an electric field has a great influence on the stability of the capsule after the curing process [42]. Jia et al. used an AC field to initiate the release process and control the release direction of actives encapsulated in single- and dual-core double-emulsion droplets [172]. Deng et al. demonstrated a two-step selective controlled release under an AC field based on the discrepancy of the shell thickness or core conductivity of the encapsulated droplet [173]. Moreover, some studies showed that the electric field can be used to manipulate the particles inside the droplet or on the droplet interface. Based on this techniques, Nudurupati et al. demonstrated that it was able to extract particles inside the droplet [170] (Figure 11C). Rozynek et al. suggested that the release rate of some areas of the droplet, named releasing holes, can be increased by rearranging the particles using an electric field [171] (Figure 11D).

9. Conclusions and Future Perspectives

In this review, we discussed the droplet EHD in detail. We discussed the features such as droplet deformation, instabilities, and breakups for single-phase and multi-phase droplets. The single-phase droplet deforms prolate or oblate and undergoes instabilities above a threshold value of the electric field. The low limit deformation is well-predicted by so-called Taylor’s leaky dielectric model, which is the most established theoretical model to date for the droplet electrohydrodynamics, neglects the effects of diffuse charge layers near the interfaces, and thus is incomplete to explain the electro-migration of droplets under a direct current electric field. Therefore, an effort to derive a complete model incorporating the electro-kinetic theory for liquid-liquid interface is required. Interesting phenomena, such as tip streaming or equatorial streaming, can be used to produce many tiny droplets. Particles at the emulsion interface can be manipulated to produce intriguing assemblies, such as belts, chains, and caps, by controlling the physical and electrical parameters. High particle concentration can produce colloids undergoing rotating flows under a strong electric field. Use of surfactant molecules alter the deformation and droplet breakup dynamics. Development of simulation models are required to understand the electrohydrodynamics of more complex interfaces, and those laden with surfactant or coated with particles merit further investigations. A double-emulsion droplet shows different breakups and can be controlled to produce Janus and ternary droplets. Simulating three-phase emulsion droplets under an electric field is likely to display further interesting dynamics.
The applications of droplet EHD can be broad, such as mass production of micro droplets, functional particle synthesis, development of micro-rotors, rigid ellipsoids, or complex shapes, and so on. Such a method could then be implemented on a microfluidic platform, if a large quantity of rigid, deformed complex shapes is to be produced. However, unfortunately the demonstration of droplet EHD on a mass scale to produce species or development of new products has still not been fully demonstrated. Still, there remain many issues and challenges that are either less explored or not explored at all and warrant further investigations.

Author Contributions

Conceptualization, M.S.A. and J.L.; investigation, M.S.A., R.S., and S.C.; writing—original draft preparation, M.S.A. and R.S.; writing—review and editing, M.S.A., R.S., S.C., and J.L.; supervision, J.L.; project administration, J.L.; funding acquisition, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2020R1A2C3010568) and the Korea Environment Industry & Technology Institute (KEITI) through its Ecological Imitation-based Environmental Pollution Management Technology Development Project, funded by the Korea Ministry of Environment (MOE) (2019002790003). The study was also supported by the New and Renewable Energy Core Technology Program of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) with a financial grant from the Ministry of Trade, Industry and Energy, Republic of Korea (No. 20193010014740).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Taylor, G. Studies in electrohydrodynamics. I. The circulation produced in a drop by electrical field. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1966, 291, 159–166. [Google Scholar]
  2. Melcher, J.R.; Taylor, G.I. Electrohydrodynamics: A Review of the Role of Interfacial Shear Stresses. Annu. Rev. Fluid Mech. 1969, 1, 111–146. [Google Scholar] [CrossRef]
  3. Ajayi, O. A note on Taylor’s electrohydrodynamic theory. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1978, 364, 499–507. [Google Scholar]
  4. Das, D.; Saintillan, D. A nonlinear small-deformation theory for transient droplet electrohydrodynamics. J. Fluid Mech. 2017, 810, 225–253. [Google Scholar] [CrossRef] [Green Version]
  5. Esmaeeli, A.; Sharifi, P. Transient electrohydrodynamics of a liquid drop. Phys. Rev. E 2011, 84, 036308. [Google Scholar] [CrossRef]
  6. Feng, J.Q.; Scott, T.C. A computational analysis of electrohydrodynamics of a leaky dielectric drop in an electric field. J. Fluid Mech. 1996, 311, 289–326. [Google Scholar] [CrossRef]
  7. Torza, S.; Cox, R.; Mason, S. Electrohydrodynamic deformation and bursts of liquid drops. Phil. Trans. R. Soc. Lond. A 1971, 269, 295–319. [Google Scholar]
  8. Bentenitis, N.; Krause, S. Droplet deformation in dc electric fields: The extended leaky dielectric model. Langmuir 2005, 21, 6194–6209. [Google Scholar] [CrossRef]
  9. Zhang, J.; Zahn, J.D.; Lin, H. Transient solution for droplet deformation under electric fields. Phys. Rev. E 2013, 87, 043008. [Google Scholar] [CrossRef] [Green Version]
  10. Zholkovskij, E.K.; Masliyah, J.H.; Czarnecki, J. An electrokinetic model of drop deformation in an electric field. J. Fluid Mech. 2002, 472, 1–27. [Google Scholar] [CrossRef]
  11. Moriya, S.; Adachi, K.; Kotaka, T. Deformation of droplets suspended in viscous media in an electric field. 1. Rate of deformation. Langmuir 1986, 2, 155–160. [Google Scholar] [CrossRef]
  12. Zabarankin, M. Small deformation theory for two leaky dielectric drops in a uniform electric field. Proc. R. Soc. A Math. Phys. Eng. Sci. 2020, 476, 20190517. [Google Scholar] [CrossRef] [Green Version]
  13. Abbasi, M.S.; Song, R.; Kim, H.; Lee, J. Multimodal breakup of a double emulsion droplet under an electric field. Soft Matter 2019, 15, 2292–2300. [Google Scholar] [CrossRef]
  14. Abbasi, M.S.; Song, R.; Kim, J.; Lee, J. Electro-hydrodynamic behavior and interface instability of double emulsion droplets under high electric field. J. Electrost. 2017, 85, 11–22. [Google Scholar] [CrossRef]
  15. Abbasi, M.S.; Song, R.; Kim, S.M.; Kim, H.; Lee, J. Mono-emulsion droplet stretching under direct current electric field. Soft Matter 2019, 15, 2328–2335. [Google Scholar] [CrossRef] [PubMed]
  16. Abbasi, M.S.; Song, R.; Lee, J. Breakups of an encapsulated surfactant-laden aqueous droplet under a DC electric field. Soft Matter 2019, 15, 8905–8911. [Google Scholar] [CrossRef] [PubMed]
  17. Bararnia, H.; Ganji, D. Breakup and deformation of a falling droplet under high voltage electric field. Adv. Powder Technol. 2013, 24, 992–998. [Google Scholar] [CrossRef]
  18. Brosseau, Q.; Vlahovska, P.M. Streaming from the Equator of a Drop in an External Electric Field. Phys. Rev. Lett. 2017, 119, 034501. [Google Scholar] [CrossRef]
  19. Collins, R.T.; Jones, J.J.; Harris, M.T.; Basaran, O.A. Electrohydrodynamic tip streaming and emission of charged drops from liquid cones. Nat. Phys. 2008, 4, 149. [Google Scholar] [CrossRef]
  20. Das, S.; Thaokar, R.M. Large-deformation electrohydrodynamics of an elastic capsule in a DC electric field. J. Fluid Mech. 2018, 841, 489–520. [Google Scholar] [CrossRef] [Green Version]
  21. Dommersnes, P.; Rozynek, Z.; Mikkelsen, A.; Castberg, R.; Kjerstad, K.; Hersvik, K.; Fossum, J.O. Active structuring of colloidal armour on liquid drops. Nat. Commun. 2013, 4, 2066. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Dubash, N.; Mestel, A.J. Breakup behavior of a conducting drop suspended in a viscous fluid subject to an electric field. Phys. Fluids 2007, 19, 072101. [Google Scholar] [CrossRef]
  23. Eggleton, C.D.; Tsai, T.-M.; Stebe, K.J. Tip streaming from a drop in the presence of surfactants. Phys. Rev. Lett. 2001, 87, 048302. [Google Scholar] [CrossRef] [Green Version]
  24. Ha, J.-W.; Yang, S.-M. Deformation and breakup of Newtonian and non-Newtonian conducting drops in an electric field. J. Fluid Mech. 2000, 405, 131–156. [Google Scholar] [CrossRef]
  25. Ha, J.-W.; Yang, S.-M. Breakup of a multiple emulsion drop in a uniform electric field. J. Colloid Interface Sci. 1999, 213, 92–100. [Google Scholar] [CrossRef] [Green Version]
  26. Karyappa, R.B.; Deshmukh, S.D.; Thaokar, R.M. Breakup of a conducting drop in a uniform electric field. J. Fluid Mech. 2014, 754, 550–589. [Google Scholar] [CrossRef]
  27. Luo, X.; Yan, H.; Huang, X.; Yang, D.; Wang, J.; He, L. Breakup characteristics of aqueous droplet with surfactant in oil under direct current electric field. J. Colloid Interface Sci. 2017, 505, 460–466. [Google Scholar] [CrossRef]
  28. Mikkelsen, A.; Dommersnes, P.; Rozynek, Z.; Gholamipour-Shirazi, A.; Carvalho, M.d.S.; Fossum, J.O.J.M. Mechanics of Pickering drops probed by electric Field–Induced Stress. Materials 2017, 10, 436. [Google Scholar] [CrossRef] [Green Version]
  29. Mikkelsen, A.; Wojciechowski, J.; Rajňák, M.; Kurimský, J.; Khobaib, K.; Kertmen, A.; Rozynek, Z.J.M. Electric field-driven assembly of sulfonated polystyrene microspheres. Materials 2017, 10, 329. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Spasic, A.M.; Jovanovic, J.M.; Manojlovic, V.; Jovanovic, M. Breaking of double emulsions based on electrohydrodynamics principles. J. Colloid Interface Sci. 2016, 479, 165–172. [Google Scholar] [CrossRef]
  31. Ahn, K.; Kerbage, C.; Hunt, T.P.; Westervelt, R.; Link, D.R.; Weitz, D. Dielectrophoretic manipulation of drops for high-speed microfluidic sorting devices. Appl. Phys. Lett. 2006, 88, 024104. [Google Scholar] [CrossRef] [Green Version]
  32. Fernández De La Mora, J. The fluid dynamics of Taylor cones. Annu. Rev. Fluid Mech. 2007, 39, 217–243. [Google Scholar] [CrossRef]
  33. Reznik, S.; Yarin, A.; Zussman, E.; Bercovici, L. Evolution of a compound droplet attached to a core-shell nozzle under the action of a strong electric field. Phys. Fluid 2006, 18, 062101. [Google Scholar] [CrossRef]
  34. Sankaran, S.; Saville, D. Experiments on the stability of a liquid bridge in an axial electric field. Phys. Fluids A Fluid Dyn. 1993, 5, 1081–1083. [Google Scholar] [CrossRef]
  35. Geng, H.; Feng, J.; Stabryla, L.M.; Cho, S.K. Dielectrowetting manipulation for digital microfluidics: Creating, transporting, splitting, and merging of droplets. Lab Chip 2017, 17, 1060–1068. [Google Scholar] [CrossRef]
  36. Jia, Y.; Ren, Y.; Hou, L.; Liu, W.; Deng, X.; Jiang, H. Sequential coalescence enabled two-step microreactions in triple-core double-emulsion droplets triggered by an electric field. Small 2017, 13, 1702188. [Google Scholar] [CrossRef]
  37. Liu, L.; Yang, J.-P.; Ju, X.-J.; Xie, R.; Liu, Y.-M.; Wang, W.; Zhang, J.-J.; Niu, C.H.; Chu, L.-Y. Monodisperse core-shell chitosan microcapsules for pH-responsive burst release of hydrophobic drugs. Soft Matter 2011, 7, 4821–4827. [Google Scholar] [CrossRef]
  38. Loscertales, I.G.; Barrero, A.; Guerrero, I.; Cortijo, R.; Marquez, M.; Ganan-Calvo, A. Micro/nano encapsulation via electrified coaxial liquid jets. Science 2002, 295, 1695–1698. [Google Scholar] [CrossRef]
  39. Mura, S.; Nicolas, J.; Couvreur, P. Stimuli-responsive nanocarriers for drug delivery. Nat. Mater. 2013, 12, 991. [Google Scholar] [CrossRef]
  40. Seifert, T.; Sowade, E.; Roscher, F.; Wiemer, M.; Gessner, T.; Baumann, R.R. Additive manufacturing technologies compared: Morphology of deposits of silver ink using inkjet and aerosol jet printing. Ind. Eng. Chem. Res. 2015, 54, 769–779. [Google Scholar] [CrossRef]
  41. Song, R.; Abbasi, M.S.; Lee, J. Fabrication of 3D printed modular microfluidic system for generating and manipulating complex emulsion droplets. Microfluid. Nanofluid. 2019, 23, 92. [Google Scholar] [CrossRef]
  42. Tucker-Schwartz, A.K.; Bei, Z.; Garrell, R.L.; Jones, T.B. Polymerization of Electric Field-Centered Double Emulsion Droplets to Create Polyacrylate Shells. Langmuir 2010, 26, 18606–18611. [Google Scholar] [CrossRef]
  43. Xie, J.; Jiang, J.; Davoodi, P.; Srinivasan, M.P.; Wang, C.-H. Electrohydrodynamic atomization: A two-decade effort to produce and process micro-/nanoparticulate materials. Chem. Eng. Sci. 2015, 125, 32–57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Bhaumik, S.K.; Roy, R.; Chakraborty, S.; DasGupta, S. Low-voltage electrohydrodynamic micropumping of emulsions. Sens. Actuators B Chem. 2014, 193, 288–293. [Google Scholar] [CrossRef]
  45. Wehking, J.D.; Kumar, R. Droplet actuation in an electrified microfluidic network. Lab Chip 2015, 15, 793–801. [Google Scholar] [CrossRef]
  46. Xi, H.-D.; Guo, W.; Leniart, M.; Chong, Z.Z.; Tan, S.H. AC electric field induced droplet deformation in a microfluidic T-junction. Lab Chip 2016, 16, 2982–2986. [Google Scholar] [CrossRef] [PubMed]
  47. Guan, X.; Hou, L.; Ren, Y.; Deng, X.; Lang, Q.; Jia, Y.; Hu, Q.; Tao, Y.; Liu, J.; Jiang, H. A dual-core double emulsion platform for osmolarity-controlled microreactor triggered by coalescence of encapsulated droplets. Biomicrofluidics 2016, 10, 034111. [Google Scholar] [CrossRef] [Green Version]
  48. Huo, M.; Guo, Y.J.P. Electric Field Enhances Shear Resistance of Polymer Melts via Orientational Polarization in Microstructures. Polymers 2020, 12, 335. [Google Scholar] [CrossRef] [Green Version]
  49. Im, D.J.; Noh, J.; Moon, D.; Kang, I.S. Electrophoresis of a charged droplet in a dielectric liquid for droplet actuation. Anal. Chem. 2011, 83, 5168–5174. [Google Scholar] [CrossRef]
  50. Lecuyer, S.; Ristenpart, W.; Vincent, O.; Stone, H. Electrohydrodynamic size stratification and flow separation of giant vesicles. Appl. Phys. Lett. 2008, 92, 104105. [Google Scholar] [CrossRef]
  51. Li, M.; Li, D. Vortices around Janus droplets under externally applied electrical field. Microfluid. Nanofluid. 2016, 20, 79. [Google Scholar] [CrossRef] [Green Version]
  52. Nguyen, V.D.; Byun, D. Mechanism of electrohydrodynamic printing based on ac voltage without a nozzle electrode. Appl. Phys. Lett. 2009, 94, 173509. [Google Scholar] [CrossRef]
  53. Saville, D.A. Electrohydrodynamics: The Taylor-Melcher Leaky Dielectric Model. Annu. Rev. Fluid Mech. 1997, 29, 27–64. [Google Scholar] [CrossRef]
  54. Schnitzer, O.; Yariv, E. The Taylor–Melcher leaky dielectric model as a macroscale electrokinetic description. J. Fluid Mech. 2015, 773, 1–33. [Google Scholar] [CrossRef] [Green Version]
  55. Bazant, M.Z. Electrokinetics meets electrohydrodynamics. J. Fluid Mech. 2015, 782, 1–4. [Google Scholar] [CrossRef] [Green Version]
  56. Mori, Y.; Young, Y.N. From electrodiffusion theory to the electrohydrodynamics of leaky dielectrics through the weak electrolyte limit. J. Fluid Mech. 2018, 855, 67–130. [Google Scholar] [CrossRef] [Green Version]
  57. Tomar, G.; Gerlach, D.; Biswas, G.; Alleborn, N.; Sharma, A.; Durst, F.; Welch, S.W.J.; Delgado, A. Two-phase electrohydrodynamic simulations using a volume-of-fluid approach. J. Comput. Phys. 2007, 227, 1267–1285. [Google Scholar] [CrossRef]
  58. Emdadi, M.; Pournaderi, P. Numerical simulation of conducting droplet impact on a surface under an electric field. Acta Mech. 2020, 231, 1083–1103. [Google Scholar] [CrossRef]
  59. Lac, E.; Homsy, G. Axisymmetric deformation and stability of a viscous drop in a steady electric field. J. Fluid Mech. 2007, 590, 239–264. [Google Scholar] [CrossRef]
  60. Lanauze, J.A.; Walker, L.M.; Khair, A.S. Nonlinear electrohydrodynamics of slightly deformed oblate drops. J. Fluid Mech. 2015, 774, 245–266. [Google Scholar] [CrossRef]
  61. Sherwood, J. Breakup of fluid droplets in electric and magnetic fields. J. Fluid Mech. 1988, 188, 133–146. [Google Scholar] [CrossRef] [Green Version]
  62. Fernández, A.; Tryggvason, G.; Che, J.; Ceccio, S.L. The effects of electrostatic forces on the distribution of drops in a channel flow: Two-dimensional oblate drops. Phys. Fluids 2005, 17, 093302. [Google Scholar] [CrossRef]
  63. Haywood, R.; Renksizbulut, M.; Raithby, G. Transient deformation of freely-suspended liquid droplets in electrostatic fields. AIChE J. 1991, 37, 1305–1317. [Google Scholar] [CrossRef]
  64. Olsson, E.; Kreiss, G. A conservative level set method for two phase flow. J. Comput. Phys. 2005, 210, 225–246. [Google Scholar] [CrossRef]
  65. Olsson, E.; Kreiss, G.; Zahedi, S. A conservative level set method for two phase flow II. J. Comput. Phys. 2007, 225, 785–807. [Google Scholar] [CrossRef]
  66. Paknemat, H.; Pishevar, A.; Pournaderi, P. Numerical simulation of drop deformations and breakup modes caused by direct current electric fields. Phys. Fluids 2012, 24, 102101. [Google Scholar] [CrossRef]
  67. Soni, P.; Juvekar, V.A.; Naik, V.M. Investigation on dynamics of double emulsion droplet in a uniform electric field. J. Electrost. 2013, 71, 471–477. [Google Scholar] [CrossRef]
  68. Jiang, Z.; Gan, Y.; Luo, Y. Effect of viscosity ratio on the dynamic response of droplet deformation under a steady electric field. Phys. Fluids 2020, 32, 053301. [Google Scholar]
  69. Huang, X.; He, L.; Luo, X.; Yin, H.; Yang, D. Deformation and coalescence of water droplets in viscous fluid under a direct current electric field. Int. J. Multiphas. Flow 2019, 118, 1–9. [Google Scholar] [CrossRef]
  70. Huang, X.; He, L.; Luo, X.; Yin, H. Droplet dynamic response in low-viscosity fluid subjected to a pulsed electric field and an alternating electric field. AIChE J. 2020, 66, e16869. [Google Scholar] [CrossRef]
  71. Ha, J.-W.; Yang, S.-M. Effect of nonionic surfactant on the deformation and breakup of a drop in an electric field. J. Colloid Interface Sci. 1998, 206, 195–204. [Google Scholar] [CrossRef]
  72. Ha, J.-W.; Yang, S.-M. Effects of surfactant on the deformation and stability of a drop in a viscous fluid in an electric field. J. Colloid Interface Sci. 1995, 175, 369–385. [Google Scholar] [CrossRef]
  73. Gu, Y.; Li, D. Electric charge on small silicone oil droplets dispersed in ionic surfactant solutions. J. Colloids Surf. A Physicochem. Eng. Asp. 1998, 139, 213–225. [Google Scholar] [CrossRef]
  74. Edwards, M.; Wu, X.; Wu, J.-S.; Huang, J.; H, K. Electric-field effects on a droplet microemulsion. Phys. Rev. E 1998, 57, 797. [Google Scholar] [CrossRef]
  75. Mousavichoubeh, M.; Shariaty-Niassar, M.; Ghadiri, M. The effect of interfacial tension on secondary drop formation in electro-coalescence of water droplets in oil. Chem. Eng. Sci. 2011, 66, 5330–5337. [Google Scholar] [CrossRef]
  76. Thiam, A.R.; Bremond, N.; Bibette, J. Adhesive emulsion bilayers under an electric field: From unzipping to fusion. Phys. Rev. Lett. 2011, 107, 068301. [Google Scholar] [CrossRef] [PubMed]
  77. Lee, S.-H.; Nguyen, X.H.; Ko, H.S. Study on droplet formation with surface tension for electrohydrodynamic inkjet nozzle. J. Mech. Sci. Technol. 2012, 26, 1403–1408. [Google Scholar] [CrossRef]
  78. Nganguia, H.; Young, Y.-N.; Vlahovska, P.M.; Blawzdziewicz, J.; Zhang, J.; Lin, H. Equilibrium electro-deformation of a surfactant-laden viscous drop. Phys. Fluids 2013, 25, 092106. [Google Scholar] [CrossRef] [Green Version]
  79. Wuzhang, J.; Song, Y.; Sun, R.; Pan, X.; Li, D. Electrophoretic mobility of oil droplets in electrolyte and surfactant solutions. Electrophoresis 2015, 36, 2489–2497. [Google Scholar] [CrossRef]
  80. He, L.; Yan, H.; Luo, X.; Wang, J.; Huang, X.; Cao, J.; Yang, D. An experimental investigation on the deformation of alkali-surfactant-polymer droplet under AC electric field. Colloid Polym. Sci. 2016, 294, 1651–1658. [Google Scholar] [CrossRef]
  81. Tuček, J.; Beránek, P.; Vobecká, L.; Slouka, Z.; Přibyl, M. Electric field driven addressing of oil-in-water droplets in the presence of gradients of ionic and nonionic surfactants. IEEE Trans. Ind. Appl. 2016, 52, 4337–4344. [Google Scholar] [CrossRef]
  82. Mandal, S.; Bandopadhyay, A.; Chakraborty, S. Dielectrophoresis of a surfactant-laden viscous drop. Phys. Fluids 2016, 28, 062006. [Google Scholar] [CrossRef]
  83. Luo, X.; Huang, X.; Yan, H.; Yang, D.; Wang, J.; He, L. Breakup modes and criterion of droplet with surfactant under direct current electric field. Chem. Eng. Res. Des. Stud. 2018, 132, 822–830. [Google Scholar] [CrossRef]
  84. Ervik, Å.; Penne, T.E.; Hellesø, S.M.; Munkejord, S.T.; Müller, B. Influence of surfactants on the electrohydrodynamic stretching of water drops in oil. Int. J. Multiphas. Flow 2018, 98, 96–109. [Google Scholar] [CrossRef] [Green Version]
  85. Poddar, A.; Mandal, S.; Bandopadhyay, A.; Chakrabort, S. Sedimentation of a surfactant-laden drop under the influence of an electric field. arXiv 2018, arXiv:1801.06317. [Google Scholar] [CrossRef] [Green Version]
  86. Sorgentone, C.; Tornberg, A.-K.; Vlahovska, P.M. A 3D boundary integral method for the electrohydrodynamics of surfactant-covered drops. J. Comput. Phys. 2019, 389, 111–127. [Google Scholar] [CrossRef] [Green Version]
  87. Poddar, A.; Mandal, S.; Bandopadhyay, A.; Chakraborty, S. Electrorheology of a dilute emulsion of surfactant-covered drops. J. Fluid Mech. 2019, 881, 524–550. [Google Scholar] [CrossRef] [Green Version]
  88. Ahmad, K.; Ho, C.; Fong, W.; Toji, D. Properties of palm oil-in-water emulsions stabilized by nonionic emulsifiers. J. Colloid Interface Sci. 1996, 181, 595–604. [Google Scholar] [CrossRef]
  89. Mandal, S.; Ghosh, U.; Chakraborty, S. Effect of surfactant on motion and deformation of compound droplets in arbitrary unbounded Stokes flows. J. Fluid Mech. 2016, 803, 200–249. [Google Scholar] [CrossRef]
  90. Das, D.; Saintillan, D. Electrohydrodynamics of viscous drops in strong electric fields: Numerical simulations. J. Fluid Mech. 2017, 829, 127–152. [Google Scholar] [CrossRef] [Green Version]
  91. Ouriemi, M.; Vlahovska, P.M. Electrohydrodynamics of particle-covered drops. J. Fluid Mech. 2014, 751, 106–120. [Google Scholar] [CrossRef]
  92. Dommersnes, P.; Mikkelsen, A.; Fossum, J.O. Electro-hydrodynamic propulsion of counter-rotating Pickering drops. Eur. Phys. J. Spec. Top. 2016, 225, 699–706. [Google Scholar] [CrossRef] [Green Version]
  93. Behjatian, A.; Esmaeeli, A. Transient electrohydrodynamics of compound drops. Acta Mech. 2015, 226, 2581. [Google Scholar] [CrossRef]
  94. Schoeler, A.M.; Josephides, D.N.; Chaurasia, A.S.; Sajjadi, S.; Mesquida, P. Electrophoretic manipulation of multiple-emulsion droplets. Appl. Phys. Lett. 2014, 104. [Google Scholar] [CrossRef] [Green Version]
  95. Soni, P.; Dixit, D.; Juvekar, V.A. Effect of conducting core on the dynamics of a compound drop in an AC electric field. Phys. Fluids 2017, 29, 112108. [Google Scholar] [CrossRef]
  96. O’Konski, C.T.; Thacher, H.C., Jr. The distortion of aerosol droplets by an electric field. J. Phys. Chem. 1953, 57, 955–958. [Google Scholar] [CrossRef]
  97. Allan, R.; Mason, S. Particle behaviour in shear and electric fields II. Rigid rods and spherical doublets. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1962, 267, 62–76. [Google Scholar]
  98. Jackson, J.D. Classical Electrodynamics; John Wiley & Sons: Hoboken, NJ, USA, 1998. [Google Scholar]
  99. Stratton, J.A. Electromagnetic Theory; John Wiley & Sons: Hoboken, NJ, USA, 2007; Volume 33. [Google Scholar]
  100. Lamb, H.J.P. Hydrodynamics; Cambridge University Press: Cambridge, UK, 1932; pp. 134–139. [Google Scholar]
  101. Esmaeeli, A.; Halim, M.A. Electrohydrodynamics of a liquid drop in AC electric fields. Acta Mech. 2018, 229, 3943–3962. [Google Scholar] [CrossRef]
  102. Dubash, N.; Mestel, A. Behaviour of a conducting drop in a highly viscous fluid subject to an electric field. J. Fluid Mech. 2007, 581, 469–493. [Google Scholar] [CrossRef]
  103. Lanauze, J.A.; Walker, L.M.; Khair, A.S. The influence of inertia and charge relaxation on electrohydrodynamic drop deformation. Phys. Fluids 2013, 25, 112101. [Google Scholar] [CrossRef]
  104. Nganguia, H.; Young, Y.-N.; Layton, A.; Lai, M.-C.; Hu, W.-F. Electrohydrodynamics of a viscous drop with inertia. Phys. Rev. E 2016, 93, 053114. [Google Scholar] [CrossRef] [Green Version]
  105. Vizika, O.; Saville, D.A. The electrohydrodynamic deformation of drops suspended in liquids in steady and oscillatory electric fields. J. Fluid Mech. 1992, 239, 1–21. [Google Scholar] [CrossRef]
  106. Jalaal, M.; Khorshidi, B.; Esmaeilzadeh, E. An experimental study on the motion, deformation and electrical charging of water drops falling in oil in the presence of high voltage D.C. electric field. Exp. Therm. Fluid Sci. 2010, 34, 1498–1506. [Google Scholar] [CrossRef]
  107. Mandal, S.; Chakrabarti, S.; Chakraborty, S. Effect of nonuniform electric field on the electrohydrodynamic motion of a drop in Poiseuille flow. Phys. Fluids 2017, 29, 052006. [Google Scholar] [CrossRef]
  108. Mandal, S.; Sinha, S.; Bandopadhyay, A.; Chakraborty, S. Drop deformation and emulsion rheology under the combined influence of uniform electric field and linear flow. J. Fluid Mech. 2018, 841, 408–433. [Google Scholar] [CrossRef]
  109. Esmaeeli, A.; Behjatian, A. Transient electrohydrodynamics of a liquid drop at finite Reynolds numbers. J. Fluid Mech. 2020, 893, 26. [Google Scholar] [CrossRef]
  110. Berg, G.; Lundgaard, L.E.; Abi-Chebel, N. Electrically stressed water drops in oil. Chem. Eng. Process. Process Intensif. 2010, 49, 1229–1240. [Google Scholar] [CrossRef]
  111. Gong, H.; Peng, Y.; Yang, Z.; Shang, H.; Zhang, X. Stable deformation of droplets surface subjected to a high-voltage electric field in oil. Colloids Surf. A Physicochem. Eng. Asp. 2015, 468, 315–321. [Google Scholar] [CrossRef]
  112. Esmaeeli, A. Transient electrohydrodynamics of a liquid drop in AC electric fields. Eur. Phys. J. E 2018, 41, 135. [Google Scholar] [CrossRef] [PubMed]
  113. Collins, R.T.; Sambath, K.; Harris, M.T.; Basaran, O.A. Universal scaling laws for the disintegration of electrified drops. Proc. Natl. Acad. Sci. USA 2013, 110, 4905. [Google Scholar] [CrossRef] [Green Version]
  114. Gawande, N.; Mayya, Y.; Thaokar, R. Jet and progeny formation in the Rayleigh breakup of a charged viscous drop. J. Fluid Mech. 2020. [Google Scholar] [CrossRef]
  115. Sengupta, R.; Walker, L.M.; Khair, A.S. The role of surface charge convection in the electrohydrodynamics and breakup of prolate drops. J. Fluid Mech. 2017, 833, 29–53. [Google Scholar] [CrossRef]
  116. Ha, J.-W.; Yang, S.-M. Electrohydrodynamics and electrorotation of a drop with fluid less conductive than that of the ambient fluid. Phys. Fluids 2000, 12, 764–772. [Google Scholar] [CrossRef] [Green Version]
  117. Salipante, P.F.; Vlahovska, P.M. Electrohydrodynamics of drops in strong uniform dc electric fields. Phys. Fluids 2010, 22, 112110. [Google Scholar] [CrossRef]
  118. Salipante, P.F.; Vlahovska, P.M. Electrohydrodynamic rotations of a viscous droplet. Phys. Rev. E 2013, 88, 043003. [Google Scholar] [CrossRef] [Green Version]
  119. Vlahovska, P.M. Electrohydrodynamic instabilities of viscous drops. Phys. Rev. Fluids 2016, 1, 060504. [Google Scholar] [CrossRef]
  120. Masliyah, J.H.; Bhattacharjee, S. Electrokinetic and Colloid Transport Phenomena; John Wiley & Sons: Hoboken, NJ, USA, 2006. [Google Scholar]
  121. Abbasi, M.S.; Farooq, H.; Ali, H.; Kazim, A.H.; Nazir, R.; Shabbir, A.; Cho, S.; Song, R.; Lee, J. Deformation of Emulsion Droplet with Clean and Particle-Covered Interface under an Electric Field. Materials 2020, 13, 2984. [Google Scholar] [CrossRef]
  122. Ouriemi, M.; Vlahovska, P.M. Electrohydrodynamic Deformation and Rotation of a Particle-Coated Drop. Langmuir 2015, 31, 6298–6305. [Google Scholar] [CrossRef] [Green Version]
  123. Ha, J.-W.; Yang, S.-M. Electrohydrodynamic effects on the deformation and orientation of a liquid capsule in a linear flow. Phys. Fluids 2000, 12, 1671–1684. [Google Scholar] [CrossRef] [Green Version]
  124. Nudurupati, S.; Janjua, M.; Aubry, N.; Singh, P.J.E. Concentrating particles on drop surfaces using external electric fields. Electrophoresis 2008, 29, 1164–1172. [Google Scholar] [CrossRef]
  125. Nudurupati, S.; Janjua, M.; Singh, P.; Aubry, N. Effect of parameters on redistribution and removal of particles from drop surfaces. Soft Matter 2010, 6, 1157–1169. [Google Scholar] [CrossRef]
  126. Hwang, K.; Singh, P.; Aubry, N. Destabilization of Pickering emulsions using external electric fields. Electrophoresis 2010, 31, 850–859. [Google Scholar] [CrossRef] [PubMed]
  127. Chen, G.; Tan, P.; Chen, S.; Huang, J.; Wen, W.; Xu, L. Coalescence of pickering emulsion droplets induced by an electric field. Phys. Rev. Lett. 2013, 110, 064502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  128. Cui, M.; Emrick, T.; Russell, T.P. Stabilizing liquid drops in nonequilibrium shapes by the interfacial jamming of nanoparticles. Science 2013, 342, 460–463. [Google Scholar] [CrossRef]
  129. Rozynek, Z.; Dommersnes, P.; Mikkelsen, A.; Michels, L.; Fossum, J. Electrohydrodynamic controlled assembly and fracturing of thin colloidal particle films confined at drop interfaces. Eur. Phys. J. Spec. Top. 2014, 223, 1859–1867. [Google Scholar] [CrossRef]
  130. Rozynek, Z.; Mikkelsen, A.; Dommersnes, P.; Fossum, J.O. Electroformation of Janus and patchy capsules. Nat. Commun. 2014, 5, 3945. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  131. Rozynek, Z.; Castberg, R.; Kalicka, A.; Jankowski, P.; Garstecki, P. Electric field manipulation of particles in leaky dielectric liquids. Arch. Mech. 2015, 67, 385–399. [Google Scholar]
  132. Dommersnes, P.; Fossum, J. Surface structuring of particle laden drops using electric fields. Eur. Phys. J. Spec. Top. 2016, 225, 715–728. [Google Scholar] [CrossRef]
  133. Mikkelsen, A.; Rozynek, Z.; Khobaib, K.; Dommersnes, P.; Fossum, J.O. Transient deformation dynamics of particle laden droplets in electric field. Colloids Surf. A Physicochem. Eng. Asp. 2017, 532, 252–256. [Google Scholar] [CrossRef]
  134. Mikkelsen, A.; Khobaib, K.; Eriksen, F.K.; Måløy, K.J.; Rozynek, Z. Particle-covered drops in electric fields: Drop deformation and surface particle organization. Soft Matter 2018, 14, 5442–5451. [Google Scholar] [CrossRef] [Green Version]
  135. Amah, E.; Shah, K.; Fischer, I.; Singh, P. Electrohydrodynamic manipulation of particles adsorbed on the surface of a drop. Soft Matter 2016, 12, 1663–1673. [Google Scholar] [CrossRef]
  136. Teigen, K.E.; Song, P.; Lowengrub, J.; Voigt, A. A diffuse-interface method for two-phase flows with soluble surfactants. J. Comput. Phys. 2011, 230, 375–393. [Google Scholar] [CrossRef] [Green Version]
  137. Yeh, C.-H.; Lee, M.-H.; Lin, Y.-C. Using an electro-spraying microfluidic chip to produce uniform emulsions under a direct-current electric field. Microfluid. Nanofluid. 2012, 12, 475–484. [Google Scholar] [CrossRef] [Green Version]
  138. Renardy, Y.Y.; Renardy, M.; Cristini, V. A new volume-of-fluid formulation for surfactants and simulations of drop deformation under shear at a low viscosity ratio. Eur. J. Mech. B Fluids 2002, 21, 49–59. [Google Scholar] [CrossRef] [Green Version]
  139. Ervik, Å.; Hellesø, S.M.; Munkejord, S.T.; Müller, B. Experimental and Computational Studies of Water Drops Falling through Model Oil with Surfactant and Subjected to an Electric Field. In Proceedings of the 18th IEEE International Conference on Dielectric Liquids (ICDL), Bled, Slovenia, 29 June–3 July 2014; pp. 1–6. [Google Scholar]
  140. Teigen, K.E.; Munkejord, S.T. Sharp-interface simulations of drop deformation in electric fields. IEEE Trans. Dielectr. Electr. Insul. 2009, 16, 475–482. [Google Scholar] [CrossRef]
  141. Enayati, M.; Ahmad, Z.; Stride, E.; Edirisinghe, M. One-step electrohydrodynamic production of drug-loaded micro- and nanoparticles. J. R. Soc. Interface 2010, 7, 667–675. [Google Scholar] [CrossRef] [Green Version]
  142. Enayati, M.; Chang, M.-W.; Bragman, F.; Edirisinghe, M.; Stride, E.J.C. Electrohydrodynamic preparation of particles, capsules and bubbles for biomedical engineering applications. Colloids Surf. A Physicochem. Eng. Asp. 2011, 382, 154–164. [Google Scholar] [CrossRef]
  143. Behjatian, A.; Esmaeeli, A. Electrohydrodynamics of a compound drop. Phys. Rev. E 2013, 88, 033012. [Google Scholar] [CrossRef]
  144. Soni, P.; Thaokar, R.M.; Juvekar, V.A. Electrohydrodynamics of a concentric compound drop in an AC electric field. Phys. Fluids 2018, 30, 032102. [Google Scholar] [CrossRef] [Green Version]
  145. Gouz, H.N.; Sadhal, S.S. Fluid dynamics and stability analysis of a compound droplet in an electric field. Q. J. Mech. Appl. Math. 1989, 42, 65–83. [Google Scholar] [CrossRef]
  146. Ha, J.-W.; Yang, S.-M. Fluid dynamics of a double emulsion droplet in an electric field. Phys. Fluids 1999, 11, 1029–1041. [Google Scholar] [CrossRef] [Green Version]
  147. Tsukada, T.; Mayama, J.; Sato, M.; Hozawa, M. Theoretical and Experimental Studies on the Behavior of a Compound Drop under a Uniform DC Electric Field. J. Chem. Eng. Jpn. 1997, 30, 215–222. [Google Scholar] [CrossRef] [Green Version]
  148. Bei, Z.; Jones, T.B.; Harding, D.R. Electric field centering of double-emulsion droplets suspended in a density gradient. Soft Matter 2010, 6, 2312–2320. [Google Scholar] [CrossRef]
  149. Bei, Z.-M.; Jones, T.; Tucker-Schwartz, A.; Harding, D. Electric field mediated droplet centering. Appl. Phys. Lett. 2008, 93, 184101. [Google Scholar] [CrossRef] [Green Version]
  150. Bei, Z.-M.; Jones, T.; Tucker-Schwartz, A. Forming concentric double-emulsion droplets using electric fields. J. Electrost. 2009, 67, 173–177. [Google Scholar] [CrossRef]
  151. Santra, S.; Das, S.; Chakraborty, S. Electrically modulated dynamics of a compound droplet in a confined microfluidic environment. J. Fluid Mech. 2020, 882. [Google Scholar] [CrossRef]
  152. Santra, S.; Das, S.; Chakraborty, S. Electric field-induced pinch-off of a compound droplet in Poiseuille flow. Phys. Fluids 2019, 31, 062004. [Google Scholar] [CrossRef] [Green Version]
  153. Pan, Z.H.; Men, Y.F.; Senapati, S.; Chang, H.C. Immersed AC electrospray (iACE) for monodispersed aqueous droplet generation. Biomicrofluidics 2018, 12, 044113. [Google Scholar] [CrossRef]
  154. Yin, S.; Huang, Y.; Wong, T.N.; Chong, W.H.; Ooi, K.T. Electric Scissors for Precise Generation of Organic Droplets in Microfluidics: A Universal Approach that Goes beyond Surface Wettability. J. Phys. Chem. C 2019, 123, 25643–25650. [Google Scholar] [CrossRef]
  155. Azizian, P.; Azarmanesh, M.; Dejam, M.; Mohammadi, M.; Shamsi, M.; Sanati-Nezhad, A.; Mohamad, A.A. Electrohydrodynamic formation of single and double emulsions for low interfacial tension multiphase systems within microfluidics. Chem. Eng. Sci. 2019, 195, 201–207. [Google Scholar] [CrossRef]
  156. Wang, H.; Zhao, Z.; Liu, Y.X.; Shao, C.M.; Bian, F.K.; Zhao, Y.J. Biomimetic enzyme cascade reaction system in microfluidic electrospray microcapsules. Sci. Adv. 2018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  157. Wu, D.; Yu, Y.R.; Zhao, C.; Shou, X.; Piao, Y.; Zhao, X.; Zhao, Y.J.; Wang, S.Q. NK-Cell-Encapsulated Porous Microspheres via Microfluidic Electrospray for Tumor Immunotherapy. ACS Appl. Mater. Interfaces 2019, 11, 33716–33724. [Google Scholar] [CrossRef] [PubMed]
  158. Wang, D.B.; Wang, J.F.; Yongphet, P.; Wang, X.Y.; Zuo, Z.W.; Li, B.; Zhang, W. Experimental study on electric-field-induced droplet generation and breakup in an immiscible medium. Exp. Fluids 2020, 61, 1–12. [Google Scholar] [CrossRef]
  159. Lee, K.H.; Yang, G.L.Z.; Wyslouzil, B.E.; Winter, J.O. Electrohydrodynamic Mixing-Mediated Nanoprecipitation for Polymer Nanoparticle Synthesis. ACS Appl. Polym. Mater. 2019, 1, 691–700. [Google Scholar] [CrossRef]
  160. Su, Y.; Yu, T.; Wang, G.; Zhang, C.; Liu, Z. Numerical simulation of electrohydrodynamics of a compound drop based on the ternary phase field method. Sci. Prog. 2019. [Google Scholar] [CrossRef]
  161. Lee, S.J.; Kang, J.Y.; Choi, W.; Kwak, R. Simultaneous electric production and sizing of emulsion droplets in microfluidics. Soft Matter 2020, 16, 614–622. [Google Scholar] [CrossRef]
  162. Rozynek, Z.; Khobaib, K.; Mikkelsen, A. Opening and Closing of Particle Shells on Droplets via Electric Fields and Its Applications. ACS Appl. Mater. Interfaces 2019, 11, 22840–22850. [Google Scholar] [CrossRef]
  163. Li, M.Q.; Li, D.Q. Redistribution of charged aluminum nanoparticles on oil droplets in water in response to applied electrical field. J. Nanoparticle Res. 2016. [Google Scholar] [CrossRef] [Green Version]
  164. Brosseau, Q.; Hickey, G.; Vlahovska, P.M. Electrohydrodynamic Quincke rotation of a prolate ellipsoid. Phys. Rev. Fluids 2017. [Google Scholar] [CrossRef]
  165. Pradillo, G.E.; Karani, H.; Vlahovska, P.M. Quincke rotor dynamics in confinement: Rolling and hovering. Soft Matter 2019, 15, 6564–6570. [Google Scholar] [CrossRef] [PubMed]
  166. Jákli, A.; Senyuk, B.; Liao, G.; Lavrentovich, O.D. Colloidal micromotor in smectic A liquid crystal driven by DC electric field. Soft Matter 2008, 4, 2471–2474. [Google Scholar] [CrossRef]
  167. Das, D.; Lauga, E. Active Particles Powered by Quincke Rotation in a Bulk Fluid. Phys. Rev. Lett. 2019, 122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  168. Zrinyi, M.; Nakano, M. Toward Colloidal Motors. Period. Polytech.-Chem. Eng. 2017, 61, 15–18. [Google Scholar] [CrossRef] [Green Version]
  169. Belovs, M.; Cebers, A. Quincke rotation driven flows. Phys. Rev. Fluids 2020. [Google Scholar] [CrossRef]
  170. Nudurupati, S.; Janjua, M.; Singh, P.; Aubry, N. Electrohydrodynamic removal of particles from drop surfaces. Phys. Rev. E 2009. [Google Scholar] [CrossRef]
  171. Rozynek, Z.; Jozefczak, A. Patchy colloidosomes—An emerging class of structures. Eur. Phys. J. Spec. Top. 2016, 225, 741–756. [Google Scholar] [CrossRef] [Green Version]
  172. Jia, Y.K.; Ren, Y.; Hou, L.K.; Liu, W.Y.; Jiang, T.Y.; Deng, X.K.; Tao, Y.; Jiang, H.Y. Electrically controlled rapid release of actives encapsulated in double-emulsion droplets. Lab Chip 2018, 18, 1121–1129. [Google Scholar] [CrossRef] [PubMed]
  173. Deng, X.K.; Ren, Y.K.; Hou, L.K.; Liu, W.Y.; Jiang, T.Y.; Jiang, H.Y. Compound-Droplet-Pairs-Filled Hydrogel Microfiber for Electric-Field-Induced Selective Release. Small 2019. [Google Scholar] [CrossRef]
Figure 1. Hydrodynamics of emulsion droplets under an electric field.
Figure 1. Hydrodynamics of emulsion droplets under an electric field.
Micromachines 11 00942 g001
Figure 2. Steady deformation-circulation maps of single-phase emulsion droplet under weak DC/AC electric field. (A) Deformation regimes depending on conductivity ratio (R), permittivity ratio (S), and viscosity ratio (M) as predicted by leaky dielectric theory (LDM). (B) Flow patterns surrounding the droplet in DC electric field for different regimes as mentioned in (A) [68]. (C) Flow patterns surrounding the droplet in an AC electric field for different regimes [101].
Figure 2. Steady deformation-circulation maps of single-phase emulsion droplet under weak DC/AC electric field. (A) Deformation regimes depending on conductivity ratio (R), permittivity ratio (S), and viscosity ratio (M) as predicted by leaky dielectric theory (LDM). (B) Flow patterns surrounding the droplet in DC electric field for different regimes as mentioned in (A) [68]. (C) Flow patterns surrounding the droplet in an AC electric field for different regimes [101].
Micromachines 11 00942 g002
Figure 3. Transient deformation of single-phase emulsion droplet under a weak DC/AC electric field. (A) Three different modes of deformation encountered in transient deformation of a droplet. (B) The mapping of deformation modes as mentioned in (A) on Ohnesorge number (Oh) vs. Reynolds Number (Re) parametric space [68]. (C) The effect of charge relaxation and convection on transient deformation of a silicone oil droplet in castor oil [60]. (D) Evolution of the deformation with time. (E) Evolution of the kinetic energy of the fluid system with time [101].
Figure 3. Transient deformation of single-phase emulsion droplet under a weak DC/AC electric field. (A) Three different modes of deformation encountered in transient deformation of a droplet. (B) The mapping of deformation modes as mentioned in (A) on Ohnesorge number (Oh) vs. Reynolds Number (Re) parametric space [68]. (C) The effect of charge relaxation and convection on transient deformation of a silicone oil droplet in castor oil [60]. (D) Evolution of the deformation with time. (E) Evolution of the kinetic energy of the fluid system with time [101].
Micromachines 11 00942 g003
Figure 4. Breakup of a prolate droplet under a strong DC electric field. (A) Transient evolution of a castor droplet in silicone oil and breakup at the side ends. (B) Non-dimensional stretching vs. electric capillary number leading to instability (reproduced from [15] with permission from The Royal Society of Chemistry). (C) Effects of viscosity ratio and liquid conductivity (in terms of Saville number) on the droplet shapes [114]. (D) Droplet breakup mode transition depending on electric Reynolds Number [115].
Figure 4. Breakup of a prolate droplet under a strong DC electric field. (A) Transient evolution of a castor droplet in silicone oil and breakup at the side ends. (B) Non-dimensional stretching vs. electric capillary number leading to instability (reproduced from [15] with permission from The Royal Society of Chemistry). (C) Effects of viscosity ratio and liquid conductivity (in terms of Saville number) on the droplet shapes [114]. (D) Droplet breakup mode transition depending on electric Reynolds Number [115].
Micromachines 11 00942 g004
Figure 5. Instabilities of an oblate droplet under a strong electric field. (A) Electro-rotation of a droplet above the threshold electric field of Quincke effect [90]. (B) Irregular rotational motions [118]. (C) Equatorial streaming. (D) Mapping of different instabilities, such as electro-rotation, dimpling, and equatorial streaming [18].
Figure 5. Instabilities of an oblate droplet under a strong electric field. (A) Electro-rotation of a droplet above the threshold electric field of Quincke effect [90]. (B) Irregular rotational motions [118]. (C) Equatorial streaming. (D) Mapping of different instabilities, such as electro-rotation, dimpling, and equatorial streaming [18].
Micromachines 11 00942 g005
Figure 6. Self-assembly of particles on the droplet interface by electro-hydrodynamic flows. Under a weak electric field: (A) effect of particle electric conductivity and concentration on the particle structure (reproduced from [134] with permission from The Royal Society of Chemistry). (B) Effect of particle electric conductivity and size on the particle structure [91]. (C) Opening and closing of belt [21]. Under a high electric field: (D,E) effect of particle size, concentration, and electric conductivity on the particle structure [91].
Figure 6. Self-assembly of particles on the droplet interface by electro-hydrodynamic flows. Under a weak electric field: (A) effect of particle electric conductivity and concentration on the particle structure (reproduced from [134] with permission from The Royal Society of Chemistry). (B) Effect of particle electric conductivity and size on the particle structure [91]. (C) Opening and closing of belt [21]. Under a high electric field: (D,E) effect of particle size, concentration, and electric conductivity on the particle structure [91].
Micromachines 11 00942 g006
Figure 7. Surfactant-laden droplets under a DC electric field. (A) Steady deformation represented by degree of stretching (l) vs. electric capillary number ( C a E ) at different surfactant concentrations ( C * ). (B) Stretching dynamics [27]. (C) Breakup modes. (D) Conversion between breakup modes [83].
Figure 7. Surfactant-laden droplets under a DC electric field. (A) Steady deformation represented by degree of stretching (l) vs. electric capillary number ( C a E ) at different surfactant concentrations ( C * ). (B) Stretching dynamics [27]. (C) Breakup modes. (D) Conversion between breakup modes [83].
Micromachines 11 00942 g007
Figure 8. Dynamics of a double-emulsion droplet depending on relative magnitudes of permittivity ratio and conductivity ratio under a weak electric field. (A) Free electric charge distributions for four distinct cases. Solid and dashed arrows show the direction of tangential electric stresses at the interfaces. (B) Scenarios of different flow patterns [143]. (C) Contours of electric potential lines [14].
Figure 8. Dynamics of a double-emulsion droplet depending on relative magnitudes of permittivity ratio and conductivity ratio under a weak electric field. (A) Free electric charge distributions for four distinct cases. Solid and dashed arrows show the direction of tangential electric stresses at the interfaces. (B) Scenarios of different flow patterns [143]. (C) Contours of electric potential lines [14].
Micromachines 11 00942 g008
Figure 9. Breakups of a double-emulsion droplet under a strong electric field. (A) Emulsion with conducting symmetric core [95]. (B) Emulsion with surfactant-laden aqueous core. (C) Criterion for core instability and critical electric capillary number for breakups (reproduced from [16] with permission from The Royal Society of Chemistry).
Figure 9. Breakups of a double-emulsion droplet under a strong electric field. (A) Emulsion with conducting symmetric core [95]. (B) Emulsion with surfactant-laden aqueous core. (C) Criterion for core instability and critical electric capillary number for breakups (reproduced from [16] with permission from The Royal Society of Chemistry).
Micromachines 11 00942 g009
Figure 10. Applications of electro-hydrodynamics (EHD) of emulsion droplets. (A) Electric-field-induced breakup of a liquid column from a nozzle, and small droplet formation behavior [158]. (B) Electric production of emulsion droplets in a microfluidic system [161]. (C) Droplet generation via streaming from the equator of a drop in an electric field [18]. (D) Stabilizing the shape of deformed colloidal particles by the interfacial jamming of nanoparticles [121]. (E) Optical diaphragm based on the particle self-assembly [162]. (F) Generation of an aluminum Janus particle using an electric field [163].
Figure 10. Applications of electro-hydrodynamics (EHD) of emulsion droplets. (A) Electric-field-induced breakup of a liquid column from a nozzle, and small droplet formation behavior [158]. (B) Electric production of emulsion droplets in a microfluidic system [161]. (C) Droplet generation via streaming from the equator of a drop in an electric field [18]. (D) Stabilizing the shape of deformed colloidal particles by the interfacial jamming of nanoparticles [121]. (E) Optical diaphragm based on the particle self-assembly [162]. (F) Generation of an aluminum Janus particle using an electric field [163].
Micromachines 11 00942 g010
Figure 11. Applications of EHD of emulsion droplets. (A) Rotation of glass-sphere-covered air bubbles under a DC electric field [166]. (B) Pickering droplet propulsion [92]. (C) Removal or collection of particles from a w/o emulsion using an electric field [170]. (D) The releasing holes for material transportation from/to the capsule [171].
Figure 11. Applications of EHD of emulsion droplets. (A) Rotation of glass-sphere-covered air bubbles under a DC electric field [166]. (B) Pickering droplet propulsion [92]. (C) Removal or collection of particles from a w/o emulsion using an electric field [170]. (D) The releasing holes for material transportation from/to the capsule [171].
Micromachines 11 00942 g011
Table 1. Summary of the most prominent works related to droplet electrohydrodynamics.
Table 1. Summary of the most prominent works related to droplet electrohydrodynamics.
Ref.Approach aMethod b,c Electric   Capillary   Number   ( C a E ) Deformation TypeCharacteristicsRemarks
Single-phase droplet (clean)
[1]E and ALeaky-dielectric model, o/o emulsion (S/C and C/S), C a E < 1 BothDeformationDC field
[7]E and ALeaky-dielectric model, o/o emulsion (S/C, C/S, etc.)
w/o emulsion (W/C, W/S, etc.)
C a E < 1 BothDeformation,
pinch-off,
dimpling
DC and AC field
[18]E and ALeaky-dielectric model,
o/o emulsion (S/C),electrolyte (AOT, TBAB), linear stability analysis
C a E 4 Oblate
(R/S > 1)
Streaming from equator, Quincke rotation, dimplingDC field, micro-droplet generation
[24]ELeaky-dielectric model, w/o emulsion (W/S, G/S, etc.),
o/o emulsion (C/S), non-Newtonian emulsion (PAM and XAN solution)
C a E < 1 Prolate
(R/S < 1)
Deformation, pinch-off, tip streamingDC field, effect of viscosity
[90]A and NLeaky-dielectric model, boundary element method, interfacial charge convection C a E < 3 Oblate
(R/S > 1)
Deformation, Quincke rotationDC field, 3D simulation
Single-phase droplet (surfactant-laden)
[72]E and ALeaky-dielectric model,
linear stability analysis,
o/o emulsion (PVDF/PS, S/C, C/S),
w/o emulsion (W/S), surfactant (PS-b-PMMA, Tween 60, etc.),
non-Newtonian fluid
(XAN solution)
C a E < 1 BothDeformation, bulbous end, tip streaming DC field, effect of viscosity, non-ionic surfactant
[84]E and NLevel-set method,
ghost-fluid method,
w/o emulsion (W/Marcol 52), surfactant (Span 80)
C a E < 0.3Prolate
(R/S < 1)
Deformation, breakupDC field, non-ionic surfactant
Single-phase droplet (particle-covered)
[91]E and ALeaky-dielectric model,
o/o emulsion (S/C), particle (Al, glass, PMMA, PE, etc.)
C a E < 5 Oblate
(R/S > 1)
Deformation, Quincke rotationDC field, particle self-assembly
[92]E and ALeaky-dielectric model,
o/o emulsion (S/C),
particle (PE)
1 < C a E < 8Oblate
(R/S > 1)
Quincke rotationDC field, micro-motor
Multi-phase droplet (clean)
[13]E, A, and NLeaky dielectric model,
level-set method,
o/o/o emulsion (C/S/C), electrolyte (Red-O dye)
C a E 1 Prolate/OblateDeformation, pinch-off, bulbous end, tip streamingDC field, effect of viscosity, and R/S
[93]ALeaky-dielectric model,
closed form analytical solution,
o/o/o emulsion (C/S/C, C/Corn/S, S/C/S, etc.)
C a E < 0.1Both/BothDeformationDC field
[94]Eo/w/o emulsion (S/W/S),
w/o/w/o emulsion (W/S/W/S),
C a E < 0.1Oblate/Prolate, P/O/PEccentricity stability, electro-phoresisDC field
[95]E and NLeaky dielectric model,
phase field model,
w/o/o emulsion (W/C/S)
C a E < 1 Prolate/ProlateDeformation, pinch-off, tip streaming AC field
Multi-phase droplet (surfactant-laden)
[16]E and ALeaky dielectric model,
w/o/o emulsion (W/S/C),
surfactant (Tween 80)
C a E 1 Prolate/OblateDeformation, bulbous endDC field
[25]E and ALeaky dielectric model,
o/w/o emulsion (HC/WE/S), w/o/o emulsion (WE/C/S), surfactant (Span85, Brij58)
C a E < 1 Both/OblatePinch-off, tip streamingDC field, multi-core double emulsion
a E: Experimental, A: Analytical, N: Numerical. b S: Silicone oil, C: Castor oil, W: DI water, G: Glycerol solution, HC: Hexane and chlorobenzene, WE: Water-ethyl alcohol mixture. c PAM: Polyacrylamide, XAN: Xanthan gum, AOT: Dioctyl sulfosuccinate sodium salt, TBAB: Tetrabutylammonium bromide.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abbasi, M.S.; Song, R.; Cho, S.; Lee, J. Electro-Hydrodynamics of Emulsion Droplets: Physical Insights to Applications. Micromachines 2020, 11, 942. https://doi.org/10.3390/mi11100942

AMA Style

Abbasi MS, Song R, Cho S, Lee J. Electro-Hydrodynamics of Emulsion Droplets: Physical Insights to Applications. Micromachines. 2020; 11(10):942. https://doi.org/10.3390/mi11100942

Chicago/Turabian Style

Abbasi, Muhammad Salman, Ryungeun Song, Seongsu Cho, and Jinkee Lee. 2020. "Electro-Hydrodynamics of Emulsion Droplets: Physical Insights to Applications" Micromachines 11, no. 10: 942. https://doi.org/10.3390/mi11100942

APA Style

Abbasi, M. S., Song, R., Cho, S., & Lee, J. (2020). Electro-Hydrodynamics of Emulsion Droplets: Physical Insights to Applications. Micromachines, 11(10), 942. https://doi.org/10.3390/mi11100942

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop