Open access peer-reviewed chapter

Electrified Pressure-Driven Instability in Thin Liquid Films

Written By

Hadi Nazaripoor, Adham Riad and Mohtada Sadrzadeh

Submitted: 25 August 2017 Reviewed: 21 November 2017 Published: 20 December 2017

DOI: 10.5772/intechopen.72618

From the Edited Volume

Electric Field

Edited by Mohsen Sheikholeslami Kandelousi

Chapter metrics overview

1,016 Chapter Downloads

View Full Metrics


The electrified pressure-driven instability of thin liquid films, also called electrohydrodynamic (EHD) lithography, is a pattern transfer method, which has gained much attention due to its ability in the fast and inexpensive creation of novel micro- and nano-sized features. In this chapter, the mathematical model describing the dynamics and spatiotemporal evolution of thin liquid film is presented. The governing hydrodynamic equations, intermolecular interactions, and electrostatic force applied to the film interface and assumptions used to derive the thin film equation are discussed. The electrostatic conjoining/disjoining pressure is derived based on the long-wave limit approximation since the film thickness is much smaller than the characteristic wavelength for the growth of instabilities. An electrostatic model, called an ionic liquid (IL) model, is developed which considers a finite diffuse electric layer with a comparable thickness to the film. This model overcomes the lack of assuming very large and small electrical diffuse layer, as essential elements in the perfect dielectric (PD) and the leaky dielectric (LD) models, respectively. The ion distribution within the IL film is considered using the Poisson-Nernst-Planck (PNP) model. The resulting patterns formed on the film for three cases of PD-PD, PD-IL, and IL-PD double layer system are presented and compared.


  • thin liquid films
  • electrohydrodynamic instabilities
  • electrokinetic
  • perfect dielectric
  • ionic liquids
  • micro- and nano-patterning

1. Introduction

For the past decades, researchers and scientists have been experimenting and exploring the use of electric fields in diverse range of applications: In health and biology like treating cancer [1] and cell sorting [2], in engineering and technological applications like enhancing the heat transfer [3, 4, 5, 6], colloidal hydrodynamics and stability [7, 8], and lithography [9]. The electric field is defined as a force field arising from the electric charges. Depending on the nature of material (ability to polarize) and the inherent or attained surface charges, the response in the electric field varies.

Surface instabilities can be either triggered using external mechanical, thermal, and electrical forces or via intermolecular interactions like van der Waals forces [9, 10, 11]. Development of these instabilities leads to the film disruption and formation of patterns, which is of interest in many applications. In the coating and cooling processes, determining the drainage time (i.e., the time when the film breakdown occurs) and in lithography providing insight regarding different morphological structures of the film interface are few examples in numerous applications. Proper implementation of the patterning process highly depends on the knowledge of the dynamics, instability, and morphological evolution of the interface or film. Interfacial tension and viscosity of liquid film are known as dampening factors for fluctuations on the free surfaces. In small-scale systems, intermolecular forces, which depend on material properties of the substrate, liquid film, and the bounding layer, also play a dominant role in the creation and amplifying the instabilities in thin liquid films. Electrically triggered instability of thin liquid films or often call electrohydrodynamic (EHD) patterning has gained extensive attention because of its ability in the creation of micro- and nano-sized structures ranging from single and bifocal microlens arrays [12], micro and nanochannels [13], and mushroom-shaped microfibers [14].

To investigate the electrically induced instabilities of thin films, it is necessary to have an electrostatic model to find the Maxwell stress acting on the interface as a primary driving force in the system. The molten polymer film is typically assumed to be a perfect dielectric (PD), with no free charge, or a leaky dielectric (LD) which has an infinitesimal amount of charges. This requires assuming the very large and very small electric diffuse layer compared to the film thickness as the main characteristics of the PD and LD models, respectively. In the nanofilms, the film thickness is comparable to the formed diffuse layer during the evolution process which violates the PD and LD assumptions. A general model is needed to bridge the gap between PD and LD models. In ionic liquid (IL) films, with a finite amount of free charges or ions, the diffuse layer thickness is comparable to the film thickness. In this chapter, an electrostatic model is presented for different cases of PD-PD, PD-IL, and IL-PD systems to find the net electrostatic force acting on the interface. Furthermore, the spatiotemporal evolution of the interface is investigated under each condition.


2. Deriving nonlinear equation for thin film dynamics

In this part, the nonlinear governing equation is derived for the two-dimensional (2D) thin film. Figure 1 shows the schematic of the 2D thin liquid film and bounding media sandwiched between two electrodes.

Figure 1.

Schematic of the thin film sandwiched between two electrodes.

The evolution of the thin film is described using mass and momentum balances for both thin film and bounding media. The boundary conditions are a no-slip condition on the walls, no penetration (two media are immiscible), and stress balance (normal and tangent) at the interface. It is assumed that fluid is Newtonian and incompressible. Detailed mathematical representations are as follows:

ρ i u i t + u i . u i = P i + . μ i u i + u i T + f e E1
. ρ i u i = 0 E2

where f e = ϕ accounts for external body force in which ϕ is the conjoining/disjoining pressure, μ is the dynamic viscosity, and i = 1 , 2 denotes the i th fluid phase. Here, we chose 1 for liquid film (molten polymer) and 2 for bounding media. Above equations are solved by the following boundary conditions:

u i = 0 at z = 0 & z = d E3

No slip condition and penetration at the interface:

u relative = 0 at z = h x y t E4

and kinematic boundary condition for the vertical component of velocity:

w = h t + u h x + v h y at z = h x y t E5

Normal and tangent stress balances at the interface:

n . σ 1 . n σ 2 . n = κγ + f e . n E6
t i . σ 1 . n σ 2 . n = f e . t i E7

where n and t i are normal and tangent vectors of the interface, respectively. γ is surface tension which is assumed constant here (Note that it can vary by location) and κ is the mean interfacial curvature of the film interface. Generally, in a 3D analysis, the film interface can be presented mathematically in the form of G x y z t = h x y t z . Using differential geometry, normal and tangent vectors and surface curvature can be found as bellow:

n = n x i x + n y i y + n z i z E8
n x = 1 C n h x n y = 1 C n h y n z = 1 C n E9

The scaling for the unit normal vector is C n = 1 + h x 2 + h y 2 and the tangent vector t i :

t i = t i x i x + t i y i y + t i z i z E10
t 1 x = 1 C t 1 h y t 1 y = 1 C t 1 h x t 1 z = 0 E11
t 2 x = 1 C t 2 h x t 2 y = 1 C t 2 h y t 2 z = 1 C t 2 h x 2 + h y 2 E12

The scaling for the tangent vectors are C t 1 = h x 2 + h y 2 and C t 2 = h x 2 + h y 2 + h x 2 + h y 2 2 . Finally, the surface curvature which is presented by two principle radii of curvature R 1 and R 2 is given by:

κ = 1 R 1 + 1 R 2 E13

Principle radii of curvature for the interface can be defined based on the film thickness as:

1 R 1 + 1 R 2 = 2 h x 2 1 + h y 2 + 2 h y 2 1 + h x 2 2 2 h x y h x h y 1 + h x 2 + h y 2 3 2 E14

In the case of 2D analysis, G x z t = h x t z , above relations are simplified as follows:

n x = 1 C n h x ; n y = 0 ; n z = 1 C n ; C n = 1 + h x 2 E15
t 1 x = 0 t 1 y = 1 C t 1 h x t 1 z = 0 E16

t 2 x = 1 C t 2 h x t 2 y = 0 t 2 z = 1 C t 2 h x 2 E17
C t 1 = | h x | & C t 2 = h x 1 + h x 2 ; κ = 2 h x 2 1 + h x 2 3 2 E18

Hence, normal and tangent components of stress balance Eqs. (6) and (7) for the film in the 2D case are given by

p 1 p 2 + 2 μ 1 h x 1 + h x 2 u 1 z + w 1 x + r u 2 z + w 2 x + 2 μ 1 1 h x 2 1 + h x 2 w 1 z r w 2 z
= h xx 1 + h x 2 3 2 γ E19
1 h x 2 u 1 z + w 1 x r u 2 z + w 2 x + 2 h x w 1 z u 1 x r w 2 z u 2 x = 0 E20

where r = μ 2 μ 1 1 is the bounding media to film dynamic viscosity ratio which is the small quantity in our case. For instance, water viscosity at 25 ° C is 0.894 [mPa.s or cP] and for castor oil is 985 [mPa.s or cP]. Therefore, the bounding media can be considered as an inactive media and there is no need to solve the evolution equations (Navier–Stokes and continuity) for both media. Also, as film thickness h 0 1 , inertial effects are negligible ( Re 1 , creeping flow). Here and thereafter, we just consider film evolution equations and do not use the subscript 1 for the film velocity components.

The subscripts in above equations x and z represent derivatives with respect to x and z , respectively. In what follows, scaling parameters are used to nondimensionalize governing equations and boundary conditions (capital letters are dimensionless):

x = λXy = λYz = h 0 Zu = u 0 Uv = u 0 V E21

where λ is the maximum wavelength of the surface instabilities, h 0 is the initial thickness of the film, u 0 is the maximum value of lateral velocity, and k is defined as k = h 0 λ . By using definitions in Eq. (21) and applying to mass continuity equation Eq. (2):

w = u 0 kW E22

Time is also re-scaled by employing long length scale λ and characteristic velocity u 0

t = h 0 k u 0 T E23

As a result of the long wave approximation theory k = h 0 λ 1 , flow within the film is locally parallel [15]. Pressure and conjoining/disjoining pressure are scaled as

p ϕ = k h 0 μ 1 u 0 P Φ E24

P Φ are dimensionless pressure and conjoining/disjoining pressure.

Stress balance in normal direction, Eq. (19), can be written in the scaled form as follows:

P 1 P ext + 2 H X 1 + k 2 H X 2 k 2 U Z + k W X + 2 1 k 2 H X 2 1 + k 2 H X 2 k 2 W Z = H XX 1 + k 2 H YY + H YY 1 + k 2 H XX 2 k 4 H XY H X H Y 1 + k 2 H X 2 + k 2 H Y 2 3 2 k 3 γ μ 1 u 0 E25

and in tangential direction, Eq. (20), is scaled as below:

1 k 2 H X 2 U Z + k 2 W X + 2 k 2 H X W Z U X = 0 E26

where P ext represents scaled form of p 2 in Eq. (19). P 1 P ext can be called gage pressure ( P ) without loss of generality. Rewriting Navier–Stokes and mass continuity equations, Eqs. (1) and (2), by using above definitions we have,

kRe U T + U U X + V U Y + W U Z = P + Φ X + k 2 U XX + U YY + U ZZ E27
kRe V T + U V X + V V Y + W V Z = P + Φ Y + k 2 V XX + V YY + V ZZ E28
kRe W T + U W X + V W Y + W W Z = P + Φ Z + k 4 W XX + W YY + k 2 W ZZ E29
U X + W Z = 0 E30

and boundary conditions, Eqs. (3) and (5),

U i = 0 at Z = 0 & Z = d h 0 E31
W = H T + U H X atZ = H X T E32

In the EHD pattern evolution process, the electric field destabilizes the interface of the fluids, causing the interface deformation with the height of h = f x y t . The characteristic wavelength of the growing instabilities, λ , is much larger than the initial film thickness ( k 1 ), so a “long-wave approximation” [16] is used to simplify the governing equations. To avoid losing the effects of interfacial tension, this assumption is made: k 3 γ = O 1 .

P + Φ X = U ZZ E33
P + Φ Z = 0 E34
U X + W Z = 0 E35
P = H XX at Z = H X T E36
U Z = 0 at Z = H X T E37

P + Φ is called modified pressure P ¯ . Eq. (34) shows that, the modified pressure does not change across the film thickness. From Eq. (31), (33), and (37), velocity component U is found as follows:

U = P ¯ X Z 1 2 Z H E38

Using Eqs. (38) and (35) leads to velocity component W,

W = 1 2 Z 2 1 3 Z P ¯ XX + P ¯ X H X E39

Substituting velocity component values U and W that are evaluated at Z = H X T into Eq. (32) gives the film evolution equation:

H T H 3 3 P ¯ X X = 0 E40

Replacing the modified pressure with pressure and conjoining/disjoining pressure, then using Eq. (36) for pressure results in the following thin film equation,

H T + X H 3 ∂Ψ X = 0 Ψ = 1 3 H XX Φ E41

and in a two-dimensional form and considering y-direction,

μ h t + x h 3 ψ x + y h 3 ψ y = 0 ψ = 1 3 ( γ h xx + h yy ϕ E42

2.1. Interaction potentials

Conjoining pressure, ϕ , (force acting on the film interface per unit area) is defined as the gradient of excess intermolecular interactions, ϕ = ΔG . In EHD, the electric field induces a pressure at the film interface, which is added to the natural interactions to generate excess intermolecular interactions. The conjoining pressure, ϕ , is a summation of these interaction potentials: van der Waals, electrostatic, and Born repulsive interaction potentials and is given as:

ϕ = ϕ vdW + ϕ Br + ϕ EL E43

The van der Waals interaction is the summation of Keeson, Debye, and London dispersion forces [17, 18]. This interaction is defined as ϕ vd W L = A L / 6 π h 3 and ϕ vd W U = A U / 6 π d h 3 for the lower and upper electrodes, respectively. A L and A U are effective Hamaker constant lower and upper electrodes, which depend on the materials used for electrodes and fluid layers. For instance, for three-layered systems, it is defined: A 213 = A 33 A 11 A 22 A 11 in which 1, 2, and 3 denote substrate, liquid film, and bounding fluid, respectively.

The van der Waals interaction becomes singular as h 0 and h d . To avoid nonphysical penetration of liquid to solid phase, in case of film rupture and touching, a cutoff distance, l 0 , is defined for which a short-range repulsive force, called Born repulsion, acts on the film interface [19]. This is used to maintain a minimum equilibrium liquid thickness on both electrodes and is defined as ϕ B r L = 8 B L / h 9 and ϕ B r U = 8 B U / d h 9 for the lower and upper electrodes, respectively. Coefficients B L and B U are found by setting the net conjoining pressure equal to zero at h = l 0 and h = d l 0 for lower and upper surfaces.

The electrostatic conjoining/disjoining pressure depends on the electrostatic property of liquid thin film and bounding layer, which is discussed in the following sections. Throughout this study, it is assumed that electric breakdown does not occur during the EHD patterning process.

2.2. Numerical scheme

Different numerical methods are used to track the free interfaces [20, 21], and particularly, in the EHD patterning process [22, 23, 24, 25]. The numerical methods are applied as versatile tools to monitor and visualize the transient evolution of liquid film subject to an electric field. Here, the thin film equation, Eq. (42) is solved numerically to obtain the transient behavior using finite difference method and adaptive time step solver. More details about the numerical scheme are available in [26].


3. Perfect dielectric film and bounding media

In this part, we consider the film and bounding media as two perfect dielectric media and solve Laplace equation in 1D. Governing equations and applied boundary conditions [27, 28, 29] are as follows:

2 ψ i = 0 E44
ψ 1 = ψ low at z = 0 E45
ψ 1 = ψ 2 at z = h x y t E46
n . ε D 1 n . ε D 2 = 0 at z = h x y t E47
ψ 2 = 0 at z = d E48

Solving Eq. (44) with mentioned boundary conditions give rise to electric potential distribution across the domain as follows:

ψ 1 = ε 2 ε 1 ψ low h ε 2 ε 1 1 + d z + ψ low E49
ψ 2 = ψ low h ε 2 ε 1 1 + d z d E50

The Maxwell stress tensor is defined as follows:

T = ε E x 2 1 2 E 2 E x E y E x E z E y E x E y 2 1 2 E 2 E y E z E z E x E z E y E z 2 1 2 E 2

where E is defined as E 2 = E x 2 + E y 2 + E z 2 and electric field components can be found as follows:

E x = ψ x E y = ψ y E z = ψ z E51

Traction force vector acting on the interface because of Maxwell stress is found as follows:

F = S n . T ¯ dS E52

where S is the interface area and n is the unit normal vector of the interface.

n = n x i x + n y i y + n z i z

where n is the unit normal vector into the medium. In the long-wave limit approximation, the traction forces for each layer are given as:

F 1 = 1 2 ε 1 E 1 z 2 S i z E53
F 2 = 1 2 ε 2 E 2 z 2 S i z E54

Therefore, the net force per unit area, conjoining pressure is

ϕ ES = 1 2 ε 1 ε 1 ε 2 1 ψ low ε 1 ε 2 d h ε 1 ε 2 1 2 E55

The net electrostatic force acting on the interface depends on the electric permittivity of each layer, applied potential, electrodes separation distance, and interface height.

3D snapshots were plotted to investigate the effect that a transverse electric field has on the nondimensional structural height variations over time of a PD film, with an initial thickness of 30 nm. Figure 2(a) shows the formation of small random disturbances in the liquid film that are conical in shape. These evolve with time as liquid flows from regions of lower thickness to regions of higher thickness causing them to increase in size and length and become more pillar-shaped as shown in Figure 2(b). However, once the pillars reach the upper electrode, their height ceases to increase, and the pillars begin to increase their cross-sectional area, as shown in Figure 2(c).

Figure 2.

PD film with h 0 = 30 nm , 3D spatiotemporal evolution of a PD liquid PD film (images a–c). Nondimensional times, T are: (a) 3 × 10 5 (b) 3.5 × 10 5 (c) 4.5 × 10 5 .

A detailed spatiotemporal evolution for the liquid–liquid interface instabilities in a 2D domain highlighting the different patterns formation using different initial film thickness is presented in Figure 3. The relative electric permittivity ratio of bilayer system is 0.6 and the initial lower layer film thickness is increased from 20 nm to 85 nm. As a result, four patterns of pillars, bicontinuous, holes, and roll-like features are formed. At lower thickness, 20 nm, pillars are formed (image a(ii)) but merging of neighbor pillars results in coarse final structure (image a(iii)). Bicontinuous structures are formed when the film thickness is increased to 50 nm (image b(i–iv)) which is not desired. The bicontinuous structures behave similarly to an air-in-liquid dispersion that also takes place in air-polymer systems with high filling ratios. The further increase in the film thickness leads to a columnar holes formation in the film (images c(i–iv)). In very thick films, 75 nm, the role-like features form that are spaced with micrometer distance. This type of features can be used as nanochannels in practical applications. The roll-like structures, in d(i), seem to be an organized version of the bicontinuous structures generated by the same phase inversion mechanism shown in b(i).

Figure 3.

Base case PD-PD bilayer, (a–d) images from i to iii show the 2D spatiotemporal evolution for liquid–liquid interface instabilities in a 2 domain when ϵ r = 0.6 and ϵ l = 2.5. Initial mean film thicknesses ( h 0 )are: A(i–iv) 20 nm, b(i–iv) 50 nm, c(i–iv) 70 nm, and d(i–iv) 85 nm. Initial electric field intensities ( E 0 ) are 294, 250, 227, and 212 MV/m, respectively. Reprinted with permission from (HADI NAZARIPOOR et al. LANGMUIR 2014, 30, 14734–14744). Copyright (2014) American Chemical Society.

In the EHL process, both the initial layer thickness and electric permittivity ratio play an important role in the process and can develop different patterns on the film. Figure 4 shows the different types of patterns formed as a function of the relative electric permittivity ratios of layers ‘ ϵ r ’ and nondimensional initial mean film thickness of h 0 / d , in a liquid–liquid interface. As shown in Figure 4, changing the electric permittivity ratios may lead to having two to four kinds of structures depending on the initial layer thickness. In electric permittivity ratios greater than one, two main shapes of structures form, pillars in filling ratios ( h 0 / d ) less than 0.5 and bicontinuous structures in relatively thicker films. However, as shown in Figure 4, in electric permittivity ratios less than one, four different shapes of structures can form. The generated map shown in Figure 4 provides details about the threshold values of the filling ratio in which the transition between structures happens at different electric permittivity ratios. From thermodynamics, the shape that leads to the lowest free energy in the system determines the final pattern. Figure 4 also shows that there is a critical initial film thickness for each value of ϵ r , where below it only pillars are formed.

Figure 4.

A parametric map that shows the different types of patterns formed, in 2D, as a function of the relative electric permittivity ratios of layers ‘ ϵ r ’ and nondimensional initial mean film thickness of h 0 / d , in a liquid–liquid interface. Reprinted with permission from (HADI NAZARIPOOR et al. LANGMUIR 2014, 30, 14734–14744). Copyright (2014) American Chemical Society.


4. Perfect dielectric film and ionic conductive bounding media

In this section, an ionic conductor is chosen as a bounding media and a dielectric media for thin film, for instance, salt water–oil system. Saltwater behaves like conductors due to having free ions, so Poisson equation is solved to find the electric potential distribution over the domain, and for the oil part similar to the previous section, Laplace equation is considered. Detailed mathematical procedure is

ε i 2 ψ i = ρ f i E56
ψ 1 = 0 at z = 0 E57
ψ 1 = ψ 2 at z = h x y t E58
n . ε D 1 n . ε D 2 = 0 at z = h x y t E59
ψ 2 = ψ up at z = d E60

The space charge density of the mobile ions ρ f is zero for dielectric media as they do not have any ions, but for the IL, ρ f 2 is given by [18],

ρ f 2 = i = 1 N z i e n i exp z i e ψ ψ k B T E61

Eq. (61) is called the Boltzmann distribution. Here, z i is the valence of species i , e is the magnitude of electron charge, 1.602 × 1 0 19 C , n i is the bulk ionic number concentration (say, in m 3 ), N is the number of species in the electrolyte, k B is the Boltzmann constant, 1.378 × 1 0 23 J / K , and T is the temperature in ( K ).

Ionic number concentration, n , is given by

n = 1000 N A M

with Avogadro number, N A = 6.022 × 1 0 23 mo l 1 and M being the electrolyte molar concentration (mol/L).


Generally, ρ f is defined based on ionic number concentration n i of species i th as follows:

ρ f = i = 1 N z i e n i E62

Ionic number concentration, n i , near a charged surface can be found by considering ion conservation Eq. (63) for the steady-state and equilibrium condition which leads to zero-ion flux, j i = 0 . Here, R i represents for reactions that produce/consume species i in the electrolyte and set to zero.

n i t = . j i + R i E63

Ion flux, j i , for species i , is defined by the Nernst-Plank relationship Eq. (64)

j i = n i u D i n i z i e n i D i k B T ψ E64

At the equilibrium, there is no fluid velocity, u = 0 , so for the 1D case (z-direction) it is simplified as follows:

d n i dz + z i e n i k B T dz = 0 E65

Eq. (65) has an analytical solution with the following boundary conditions:

n i = n i and ψ = ψ at z 0 E66
n i = n i and ψ = ψ at z = 0 E67

where z shows distance from charged surface, ψ and n i are bulk (electro neutral) electric potential and ion number concentration, respectively. ψ ∞i is also called reference potential, ψ ref . From now we use ψ ref instead of ψ . By using Eqs. (66) and (67), the solution is:

n i = n i exp z i e ψ ψ ref k B T E68

Substituting Eq. (68) into Eq. (62) results in Eq. (61).

Proof. end

Using Boltzmann distribution (relation Eq. (61)) for the free space charge density, ρ f , in the Poisson Eq. (56) give rise to well-known Poisson-Boltzmann (PB) equation

ε 2 2 ψ 2 = i = 1 N z i e n i exp z i e ψ 2 ψ ref k B T E69

and for the 1D case

ε 2 d 2 ψ 2 d z 2 = i = 1 N z i e n i exp z i e ψ 2 ψ ref k B T E70

PB equation is simplified for the monovalent, N = 2 , symmetric, z : z , electrolyte solution which is known as the Gouy-Chapman theory

ε 2 d 2 ψ 2 d z 2 = 2 ze n sinh ze ψ 2 ψ ref k B T E71

Governing equations in the long-wave limit condition [29] is simplified as follows:

d 2 ψ 1 d z 2 = 0 E72
ε 2 d 2 ψ 2 d z 2 = 2 ze n sinh ze ψ ψ ref k B T E73

and boundary conditions

ψ 1 = 0 at z = 0 E74
ψ 1 = ψ 2 at z = h x y t E75
ε 1 d ψ 1 dz = ε 2 d ψ 2 dz at z = h x y t E76
ψ 2 = ψ up at z = d E77

The following scaling parameters are defined to nondimensionalize the governing equations and the boundary conditions:

Ψ = zeψ k B T , Z = z d E78

Also, a dimensionless parameter is defined, κ 1 = ε 2 k B T 2 e 2 z 2 n , called Debye length. By using these definitions, we can rewrite equations and boundary conditions in scaled form

d 2 Ψ 1 d Z 2 = 0 E79
d 2 Ψ 2 d Z 2 = κd 2 sinh Ψ Ψ ref E80

and boundary conditions

Ψ 1 = 0 at Z = 0 E81
Ψ 1 = Ψ 2 at Z = h x y t d = Z E82
ε 1 d Ψ 1 dZ = ε 2 d Ψ 2 dZ at Z = h x y t d = Z E83

Solution is as follows,

Ψ 1 = Ψ s Z Z at 0 Z Z E84
Ψ 2 = Ψ up + 2 ln 1 + exp κd Z Z tanh Ψ s Ψ up 4 1 exp κd Z Z tanh Ψ s Ψ up 4 at Z Z 1 E85

where Ψ s results from solution of

Ψ s + ε 2 ε 1 2 κd sinh Ψ s Ψ up 2 Z = 0 E86

In dimensional form

ψ 1 = ψ s h z at 0 z h E87
ψ 2 = ψ up + 2 k B T ze ln 1 + exp κ z h tanh k B T ze ψ s ψ up 4 1 exp κ z h tanh k B T ze ψ s ψ up 4 at h z d E88

After finding the electric potential distribution, we can calculate the net electrostatic force acting on the interface. By using definitions in Eqs. (53) and (54), the net force is given by,

1 S i z . F net = 1 2 ε 1 ε 1 ε 2 1 ψ s h 2 E89

So conjoining pressure becomes [29]:

ϕ ES = 1 2 ε 1 ε 1 ε 2 1 ψ s h 2 E90

Figure 5(a) shows the effects of changing the electrolytes molarity on the resultant electrostatic pressure acting on the IL-PD interface. The negative values’ curves have an electric permittivity ratio ( ϵ r = 2.5 ) and are associated with upward disjoining pressure, while the positive values’ curves have an electric permittivity ratio ( ϵ r = 0.6 ) and are associated with downward joining pressure. However, both curves experience a decay in the value of the electrostatic pressure as the interface height increases. Figure 5(b) outlines the effect of changing the applied voltage on the electrostatic pressure for both IL-PD bilayer and PD-PD bilayer. The electrolytic molarity and electric permittivity ratio ( ϵ r = 0.6 ) were kept constant for all applied voltages. The graph shows that the IL-PD bilayer experiences a much larger pressure than that of the PD bilayer for the same applied voltage (Notice at V = 10 V).

Figure 5.

Electrostatic pressure distribution versus interface height. (a) Effects of electrolyte molarity in IL-PD bilayers and (b) comparison between PD-PD and PD-IL bilayers. (a) Molarity M = 0.001 , 0.0001 , and 0.00001 mol / L and applied potential of 0.25 V and ϵ 1 = 2.5 . (b) M = 0.0001 mol / L for IL-PD, ϵ 1 = 2.5 and ϵ 1 = 4.17 . Reprinted with permission from (HADI NAZARIPOOR et al. LANGMUIR 2014, 30, 14734–14744). Copyright (2014) American Chemical Society.

Figure 6 shows a 2D spatiotemporal evolution for uniform electric field liquid–liquid interface instabilities, where (b) is the PD-PD bilayer and (c) is the IL-PD bilayer. When visually compared to one another, the number of pillars in the b(i) 3D representation of the PD-PD bilayer appears to be in the same order of magnitude of that of the c(i) 3D representation of the PD-IL bilayer. However, the difference between them is much higher than that since the physical domain size for PD-IL is actually 10 times smaller than that of the PD-PD. Moreover, Figure 6 b(ii) and c(ii) confirms this size difference and shows that the average center distance between the pillars for PD-IL bilayer is 210 nm while that of the PD-PD bilayer has an average center distance of 1336 nm. The formed pillars also disordered and dispersed randomly in the PD-IL case.

Figure 6.

(a) Schematic of the IL-PD bilayer. 3D and 2D snapshots of interface morphology of (b) PD-PD, (c) IL-PD bilayer, (i) in a 3D domain, (ii) in physical domain, when ϵ r = 0.6 and M = 100 mmol / L for IL-PD bilayer. Initial mean film thickness is: h 0 = 20 nm and nondimensional time for the plots are T = (b) 3.4 × 10 5 (c) 13.5 . ψ up = 20 V and electric field intensities E 0 are (a) 294 and (b) 995.


5. Ionic conductive film and perfect dielectric bounding media

In the previous section, the complicated case of nonlinear Poisson-Boltzmann equation Eq. (59) was considered for the IL layer. In case of low applied potential (less than 25 mV), the Debye Hückel approximation [10, 30, 31] is used to linearize the Poisson-Boltzmann equation (Eq. (80)).

d 2 ψ 1 d z 2 = κ 2 ψ 1 ψ ref E91

and ions conservation within the layer is satisfied as follows:

0 h ρ f dz = 0 h ψ 1 ψ ref dz = 0 E92

and boundary conditions

ψ 1 = ψ l at z = 0 E93
Ψ 1 = Ψ 2 at z = h x y t E94
ε 1 d ψ 1 dZ = ε 2 d ψ 2 dZ at z = h x y t E95
ψ 2 = 0 at z = 1 E96

where κ = 2000 e 2 N A M ε 1 ε 0 k B T 1 / 2 is the inverse of Debye length. Solving Eqs. (91)(96) results in the electric potential distribution within PD and IL layers as follows:

ψ 2 = z d h d ψ s E97
ψ 1 = ψ s ψ l cosh κh + 1 2 sinh κh sinh κz ψ s ψ l 2 cosh κz + ψ s + ψ l 2 E98

where ψ s is the interface potential and is determined by the following equation:

ψ s = ψ l 1 + cosh κh 1 + cosh κh + 2 ε 2 ε 1 κ d h sinh κh E99

The electrostatic pressure, in this case, is given by,

ϕ EL = 1 2 ε 2 ε 0 ε 2 ε 1 1 ψ s d h 2 E100

A 50-nm thick IL film surrounded by a 50 nm PD media, such as air, was considered to plot the effects of changing the electrolytic molarity have on the potential distribution across the film thickness. The plots obtained for the electric potential distribution are linear as shown in Figure 7(a). For molarity (M) of 0.01 mol / L , the reduction in electric potential across the film thickness is almost zero and, thus the electric field is zero. On the other hand, the molarity of 0.00001 mol / L causes the electric potential distribution to behave similarly to a PD.

Figure 7.

(a) Effect of molarity in electric potential distribution. Molarity, M = 0.01 , 0.001 , 0.0001 , and 0.00001 mol / L and corresponding Debye lengths, κ 1 are: 0.5 , 1.7 , 5.4 , and 17 nm . (b) Electrostatic pressure, left axis, and spinodal parameter, right axis, distributions versus film thickness for three molarity values of M = 0.01 , 0.001 , and 0.00001 mol / L . ϕ l = 20 V , ε 1 = 2.5 and ε 2 = 1 [31].

Figure 7(b) shows the changes in the electrostatic pressure and the spinodal parameter across the IL film for different molarities. The graph shows that as the molarity increases the concentration of ions increases and consequently the electrostatic pressure and its corresponding force increase. This relationship is even more significant for higher IL thicknesses. Additionally, the negative values for the pressure indicate that the forces are pushing the interface toward the upper electrode.

Figure 8(a) shows the schematic view of the IL film bounded with a PD film while Figure 8(b) and (c) compare the structural variations between a PD film and IL film. Additionally, Figure 8(c) emphasizes the influence of altering the film thickness on the morphology of the IL film. It is found that a faster growth of instabilities is developed by increasing the initial thickness of IL film from h 0 = 30 to 50 nm, which mostly leads to more rapid patterning. On the other hand, the formed patterns are less stable as they tend to coalesce at more advanced stages.

Figure 8.

(a) A schematic view of the IL (film)–PD (bounding layer) bilayer and electric potential distribution. 3D snapshot of interface morphology of (b) PD film (c) IL film. Initial film thickness, h 0 = (b) 30 nm c(i) 30 nm c(ii) 50 nm . ψ L = 20 V and c M = 0.0001 mol / L [31].

As a way of understanding the IL films ability to make smaller sized patterns in the EHD patterning process, molarity is increased. Figure 9 portrays the influence of molarity and the initial film thickness on the number of pillars in the EHD patterning process. In the IL film interface, the presence of large numbers of pillars indicates the creation of smaller sized pillars, given a fixed area. Changing the molarity affects the conductivity of the IL films, which in turn affect the numbers and speed of formation of the pillars. Figure 9(a) shows a plot of the number of pillars with respect to nondimensional time, for a PD film of the same initial thickness, h 0 = 30 nm. The end point of each curve represents the final number of pillars formed when the quasi-stable conditions are established. As shown in Figure 9(a), the final number of pillars formed for the PD film is 16, while that of the IL films is 20, 38, 41, and 43 for the following molarities, respectively: M = 0.00001, 0.0001, 0.001, and 0.01 mol/L.

Figure 9.

Number of formed pillars with (a) changes in molarity with a constant initial film thickness of h 0 = 30 nm and (b) changes in initial film thickness with a constant molarity of M = 0.0001 mol/L [31].

Therefore, one can deduce that when the molarity is initially increased by 10 times, the number of formed pillars is almost doubled. However, any further increase in molarity causes a slight increase in the number of formed pillars, finally reaching a plateau. Additionally, it can be noted that increasing the molarity decreases the time gap between the initial and final pillar formed. Figure 9(b) shows how the initial film thickness affects the number of pillars formed. In order to cancel the effect of changing the film thickness, nondimensional time was normalized with the initial film thickness and was plotted against the number of pillars formed, keeping the molarity constant. The results showed that the higher initial film thicknesses tend to have a faster time evolution when compared to smaller film thicknesses.


6. Conclusion

In the confined liquid layer system, the electrically induced instabilities are enlarged leading to interface deformation and pattern generation. This technique has been employed as a fast and inexpensive method for patterning of molten polymer films. In this chapter, the electrified pressure-driven instability of thin liquid films and bilayers is discussed under the long-wave approximation limit. It is shown that the difference between electrical properties of film layers results in a net electrostatic force acting normal to the interface due to Maxell stress. The thin film equation governs the dynamics and instability of film layers, and its derivation is discussed. An extensive numerical study is performed to generate a map, in bilayer systems, based on the filling ratio and the electric permittivity ratio of layers. This map provides a baseline for IL-PD and PD-IL bilayers and is used as a predictive model for the formation of various structures in bilayer systems. To create different morphologies with lower pattern size, the net electrostatic force is increased by introducing an ionic conductive property of liquid layers. Thus, an electrostatic model is developed to find the net electrostatic force ionic liquid films and bilayers. The developed model is integrated to the thin film equation, and the spatiotemporal evolution of the interface is presented to show the compact and smaller sized pattern formation compared to the base cases of perfect dielectric films.



A Hamaker constant
d i Electrodes distance, i = 1; 2
D Diffusion coefficient
e Electron charge magnitude
f e → External body force
F E Electrostatic force
F os Osmotic force
h Interface height
h 0 Initial mean film thickness
I Identity tensor
k B Boltzmann constant
l 0 Born repulsion cut off distance
L Domain length
L s Scaling factor for length
m the mobility of charges/ions
M Molarity
n → Normal vector to the interface
n ± Number concentration of ions/charges, positive or negative
n ∞ Bulk number concentration of ions/charges
N A Avogadro number
P Pressure
Pos Osmotic pressure
Species production rate in chemical reactions
t → i Tangential vector to the interface , i = 1 and 2
T Temperature
Ts Scaling factor for time
u → i Velocity vector for film and bounding layer, i = 1 and 2
u relative Relative velocity
x x direction in Cartesian coordinate
y y direction in Cartesian coordinate
z z direction in Cartesian coordinate
s Amplitude for growth of instabilities
εi Electric permittivity of film and bounding layer, i =1 and 2
εr Electric permittivity ratio of layers
ε 0 Free space electric permittivity
φ Conjoining/disjoining pressure
φBr Born repulsion pressure
φEL Electrostatic pressure
Φs Scaling factor for conjoining pressure
φT Thermocapillary pressure
φvdW van der Waals pressure
γ Relative interfacial tension between film and bounding layer
κ Inverse of Debye length
κv Wave number for growth of instabilities
κ ∗ Mean interfacial curvature of the film interface
λ max Maximum wavelength for growth of instabilities
λ Wavelength for growth of instabilities
μ Dynamic viscosity
ν Kinematic viscosity
ρ Density
ρf Free charge density
ψi Electric potential within the layers, i =1 and 2
ψr Reference electric potential
ψl Electric potential of the lower electrode
ψs Interface electric potential
ψup Electric potential of the upper electrode
Ψ Nondimensional electric potential
∆ t Time step
σ Conductivity
τc Charge relaxation time
τp Process time


  1. 1. Pavesi A, Adriani G, Tay A, Warkiani ME, Yeap WH, Wong SC, Kamm RD. Engineering a 3D microfluidic culture platform for tumor-treating field application. Scientific Reports. 2016
  2. 2. Shields CW, Reyes CD, López GP. Microfluidic cell sorting: A review of the advances in the separation of cells from debulking to rare cell isolation. Lab on a Chip. 2015;15(5):1230-1249
  3. 3. Safarnia H, Sheikholeslami M, Ganji DD. Electrohydrodynamic nanofluid flow and forced convective heat transfer in a channel. The European Physical Journal Plus. 2016;131:96
  4. 4. Sheikholeslami M, Chamkha AJ. Electrohydrodynamic free convection heat transfer of a nanofluid in a semi-annulus enclosure with a sinusoidal wall. Numerical Heat Transfer, Part A: Applications. 2016;69:781-793
  5. 5. Sheikholeslami M, Hayat T, Alsaedi A, Abelman S. Numerical analysis of EHD nanofluid force convective heat transfer considering electric field dependent viscosity. International Journal of Heat and Mass Transfer. 2017;108:2558-2565
  6. 6. Sheikholeslami M. Influence of coulomb forces on Fe3O4-H2O nanofluid thermal improvement. International Journal of Hydrogen Energy. 2017;42(2):821-829
  7. 7. Teubner M. The motion of charged colloidal particles in electric fields. The Journal of Chemical Physics. 1982;76(11):5564-5573
  8. 8. de Ven, Van and T. G., M.. Colloidal Hydrodynamics. San Diego, CA: Academic Press; 1989
  9. 9. Ning W, Russel WB. Micro- and Nano-patterns created via Electrohydrodynamic instabilities. Nano Today. 2009;4:180-192
  10. 10. Gambaryan-Roisman T. Modulation of Marangoni convection in liquid films. Advances in Colloid and Interface Science. 2015;222:319-331
  11. 11. Mukherjee R, Sharma A. Instability, self-organization and pattern formation in thin soft films. Soft Matter. 2015;11:8717--8740
  12. 12. Hong H, Tian H, Shao J, Ding Y, Jiang C, Liu H. Fabrication of bifocal microlens arrays based on controlled electrohydrodynamic reflowing of pre-patterned polymer. Journal of Micromechanics and Microengineering. 2014;24:095027
  13. 13. Li H, Yu W, Wang Y, Bu H, Liu Z, Abraham E, Desmulliez MPY. Simulation of the electrohydrodynamic instability process used in the fabrication of hierarchic and hollow micro/nanostructures. RSC Advances. 2014;4:13774-13781
  14. 14. Hu H, Tian H, Li X, Shao J, Ding Y, Liu H, An N. Biomimetic mushroom-shaped microfibers for dry adhesives by electrically induced polymer deformation. ACS Applied Materials & Interfaces. 2014;6:14167–14173
  15. 15. Atta A, Crawford DG, Koch CR, Bhattacharjee S. Influence of electrostatic and chemical heterogeneity on the electric-field-induced destabilization of thin liquid films. Langmuir. 2011;27(20):2472-2485
  16. 16. Oron A, Davis SH, Bankoff SG. Long-scale evolution of thin liquid films. Review of Modern Physics. 1997;69:931-980
  17. 17. Israelachvili JN. Intermolecular and Surface Forces. Burlington, MA: Academic Press; 2011
  18. 18. Masliyah JH, Bhattacharjee S. Electrokinetic and Colloid Transport Phenomena. Hoboken, NJ: Wiley-Interscience; 2006
  19. 19. Vladimir SM. Dewetting of solid surface: Analogy with Spinodal decomposition. Journal of Colloid and Interface Science. 1993;156:491-497
  20. 20. Harlow FH, Eddie Welch J. Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface. Physics of Fluids. 1965;8:2182-2189
  21. 21. Hirt CW, Nichols BD. Volume of fluid (VOF) method for the dynamics of free boundaries. Journal of Computational Physics. 1981;39:201–225
  22. 22. Craster RV, Matar OK. Electrically induced pattern formation in thin leaky dielectric films. Physics of Fluids (1994-present). 2005;17:1131–1198
  23. 23. Verma R, Sharma A, Kargupta K, Bhaumik J. Electric field induced instability and pattern formation in thin liquid films. Langmuir. 2005;21(8):3710-3721
  24. 24. Tian H, Shao J, Ding Y, Li X, Liu H. Numerical characterization of electrohydrodynamic micro- or nanopatterning processes based on a phase-field formulation of liquid dielectrophoresis. Langmuir. 2013;29(15):4703-4714
  25. 25. Yang Q, Li BQ, Ding Y. A numerical study of nanoscale electrohydrodynamic patterning in a liquid film. Soft Matter. 2013;9:3412-3423
  26. 26. Hadi N. Electrohydrodynamic and Thermocapillary Instability of Thin Liquid Films [thesis]. University of Alberta: University of Alberta; 2016
  27. 27. Nazaripoor H, Koch CR, Sadrzadeh M, Bhattacharjee S. Thermo-Electrohydrodynamic patterning in Nanofilms. Langmuir. 2016;32(23). DOI: 10.1021/acs.langmuir.6b01810
  28. 28. Nazaripoor H, Koch CR, Sadrzadeh M. Enhanced electrically induced micropatterning of confined thin liquid films: Thermocapillary role and its limitations. Industrial & Engineering Chemistry Research. 2017;56. DOI: 10.1021/acs.iecr.7b02814
  29. 29. Nazaripoor H, Koch CR, Bhattacharjee S. Electrical perturbations of ultrathin bilayers: Role of ionic conductive layer. Langmuir. 2014;30(49):14734-14744. DOI: 10.1021/la503839x
  30. 30. Nazaripoor H, Koch CR, Sadrzadeh M, Bhattacharjee S. Compact micro/nano electrohydrodynamic patterning: Using a thin conductive film and a patterned template. Soft Matter. 2016;12(4)
  31. 31. Nazaripoor H, Koch CR, Sadrzadeh M, Bhattacharjee S. Electrohydrodynamic patterning of ultra-thin ionic liquid films. Soft Matter. 2015;11(11):2193-2202

Written By

Hadi Nazaripoor, Adham Riad and Mohtada Sadrzadeh

Submitted: 25 August 2017 Reviewed: 21 November 2017 Published: 20 December 2017