# Viscous Potential Flow Analysis of Electroaerodynamic Instability of a Liquid Sheet Sprayed with an Air Stream

Authors: Mukesh Kumar Awasthi; Vineet K. Srivastava; M. Tamsir
Journal: Modelling and Simulation in Engineering
 (2013)

---

## Abstract

The instability of a thin sheet of viscous and dielectric liquid moving in the same direction as an air stream in the presence of a uniform horizontal electric field has been carried out using viscous potential flow theory. It is observed that aerodynamic-enhanced instability occurs if the Weber number is much less than a critical value related to the ratio of the air and liquid stream velocities, viscosity ratio of two fluids, the electric field, and the dielectric constant values. Liquid viscosity has stabilizing effect in the stability analysis, while air viscosity has destabilizing effect.


---

## Body

## 1. Introduction

The past decade has witnessed a rapid advancement in the study of instability of a liquid sheet because of its importance in several scientific and technological processes. Application of plane liquid sheets is an interesting phenomenon, which can be seen in power generation and propulsion systems [1], chemical and pharmaceutical processes [2], surface curtain coatings, and in the adhesive industry [3]. A host of efforts have been devoted to study the behavior of a thin liquid sheet sprayed with an air stream. The instability and breakup process of a thin inviscid liquid sheet in a stationary gaseous medium have been investigated by Squire [4] and Hagerty and Shea [5]. Their results show that the surface tension resists the development of instability in a liquid sheet.Fraser [6] has defined four modes of disintegration of a liquid sheet, namely, rim, wavy sheet, perforated sheet, and air impact. In rim, the disintegration takes place due to the contraction of the liquid sheet edges under the effect of surface tension. In wavy sheet, the disintegration occurs due to any small protuberance on the sheet which is subjected to two opposing forces: surface tension force, which draws the liquid back to the original undisturbed shape, and aerodynamic force, which pulls the liquid outward. If the aerodynamic force exceeds to the surface tension force, then any small disturbance present in the sheet will grow rapidly, causing sheet instability. In perforated sheetdisintegration, disturbances on the sheet puncture it when the sheet becomes thin enough, and the resulting holes expand regularly by surface tension until they coalesce, forming threads. In air impactdisintegration, the disruption of the liquid is very near to that of a twin-fluid nozzle, where two streams of air and liquid are caused to impinge together. Squire [4] has studied the disintegration according to the wavy sheet mode. Dombrowski and Johns [7] extended the above analysis including the effect of liquid viscosity, and their results are valid for very large Weber number. Joostan et al. [8] have studied the effect of capillary waves on a free and viscous liquid sheet. Rashed et al. [9] have studied the disintegration of the liquid sheet according to perforated sheet mode and found that instability occurs if the Weber number is less than a critical value related to the ratio of the air and liquid stream velocities.Ibrahim and Jackson [10] have studied the spatial instability of an incompressible liquid sheet injected into a compressible gas and found that the growth rate of disturbances increases initially with gas compressibility, reaches a maximum near gas Mach number 1, and then decreases as the gas Mach number is increased further. Ibrahim and Akpan [11] have studied the linear stability of inviscid liquid sheet emanated into an inviscid gas medium. A power series mathematical solution for the problem of instability of an inviscid liquid sheet of parabolic velocity profile emanated from a nozzle into an inviscid gas has been established by Ibrahim [12]. Nonlinear instability and breakup of an annular liquid sheet have been modeled by Ibrahim and Jog [13].There is increasing interest in the electrohydrodynamic stability as the electric field plays an important role in many practical problems of biophysics, chemical engineering, and various other fields. The effect of an applied electric field on a liquid sheet has several applications in astrophysics, chemical engineering, and industry [14]. If an electric field is applied on a liquid sheet, the instability conditions will be considerably modified. El-Sayed [15] has studied the instability of a thin sheet of dielectric liquid moving in the same direction as an air stream in the presence of a uniform horizontal electric field and observed that the electric field has stabilizing effect, and there exists a critical Weber number above which instability is suppressed by the surface tension effect.Viscous potential flow theory [16] has played an important role in studying various stability problems. In viscous potential flow, viscous term in the Navier-Stokes equation is identically zero when the vorticity is zero but the viscous stresses are not zero. Tangential stresses are not considered in the viscous potential flow theory, and viscosity enters through normal stress balance. In this theory, no-slip condition at the boundary is not enforced so that two-dimensional solutions satisfy three-dimensional solutions. Joseph et al. [17] studied viscous potential flow of Rayleigh-Taylor instability. Funada and Joseph [18] have done the viscous potential flow analysis of Kelvin-Helmholtz instability in a channel and found that the stability criterion for viscous potential flow is given by the critical value of the relative velocity. From the above study Funada and Joseph [18] concluded that the critical value is maximum when viscosity ratio equals to the density ratio. Awasthi and Agrawal [19] have studied the viscous potential flow analysis of Kelvin-Helmholtz instability of a cylindrical interface and found that the viscosity of the fluids stabilizes the interface. The rupture of a 3D stationary free liquid film under the competing effects of surface tension and van der Waals forces has been studied by Ardekani and Joseph [20] as a linearized stability problem in a purely irrotational analysis utilizing the dissipation method.In the present work, the viscous potential flow theory has been used to study the instability of thin liquid sheet moving in the same direction as an air stream and in the presence of uniform horizontal electric field. Both the liquid and air treats as incompressible and viscous with different kinematic viscosities and having horizontal velocities. In the present investigation, the sheet has been taken as two-dimensional, irrotational, and of constant thickness, considering only antisymmetric oscillations (predominant for instability as mentioned by Squire [4]) of the moving sheet. Normal mode technique has been used to solve the governing equations. A dispersion relation is derived and stability is discussed theoretically as well as numerically. Various neutral curves have been drawn to show the effect of various physical parameters such as Reynolds number, Weber number, electric field intensity, and ratio of permittivity of two fluids on the stability of the system. Finally, a comparison has been made between the results of present study and the results obtained by El-Sayed [15] for inviscid potential analysis.


## 2. Problems Formulation

Consider a two-dimensional irrotational sheet of viscous and dielectric liquid of densityρ(1), viscosity  μ(1), dielectric constantε(1), surface tensionT, and thicknessh=2a, moving with velocityU1along an air stream (on both sides of the sheet), which is of densityρ(2), viscosityμ(2), dielectric constantε(2), and moving with velocityU2. Both fluids are assumed to be incompressible and irrotational. The origin is located at the mid-plane of the liquid sheet so in the equilibrium state, the equation of interface can be expressed as  y=a. On applying small disturbances to the system, the equation of upper surface can be written as
(1)F(x,y,t)=y-a-η(x,t)=0,
whereηrepresents the displacement of liquid sheet from its equilibrium position andF represents the equation of the disturbed surface.The unit outward normal to the interface is given by(2)n=∇F|∇F|={1+(∂η∂x)2}-1/2(-∂η∂xex+ey),
and the unit outward normal of the first order is given by
(3)n=(-∂η∂xex+ey).In each fluid layer velocity is expressed as the gradient of the potential functionϕ(x,y,t),and the potential functions satisfy Laplace’s equation, that is,
(4)∇2ϕ(j)=0,j=1,2,
where∇2=∂2/∂x2+∂2/∂y2.In initial state, the potential function can be written as(5)ϕ0(j)=Ujx(j=1,2).The considered system is influenced by a uniform horizontal electric fieldE0in the positive x-direction
(6)E=E0ex.We have assumed that the quasistatic approximation is valid for the problem, and therefore, the electric field can be derived from electric scalar potential functionψ(x,y,t)such that
(7)Ej=E0ex-∇ψj,(j=1,2).Gauss’s law requires that the electric potentials also satisfy Laplace’s equation, that is,(8)∇2ψj=0,(j=1,2).


## 3.        Boundary Conditions

(i)
Kinematic condition that every particle on the interface remains on the interface leads to the following boundary conditions:(9)∂η∂t+U1∂η∂x=-∂ϕ(1)∂yaty=-a,∂η∂t+U2∂η∂x=-∂ϕ(2)∂yaty=a.
(ii)
The tangential component of the electric field must be continuous across the interface, that is,n∧⟦E⟧=0or
(10)∂ψ(2)∂x=∂ψ(2)∂x.
(iii)
There is discontinuity in the normal current across the interface; charge accumulation within a material element is balanced by conduction from bulk fluid on either side of the surface. Hence, the boundary condition corresponding to normal component of the electric field at the interface is given byn·⟦εE⟧=0  or
(11)ε(2)(∂ψ(2)∂y+E0∂η∂x)=ε(1)(∂ψ(1)∂y+E0∂η∂x).
(iv)
Interfacial condition for conservation of momentum is(12)p2-p1-2μ(2)∂2ϕ(2)∂y2+2μ(1)∂2ϕ(1)∂y2-E0[ε(2)∂ψ(2)∂x-ε(1)∂ψ(1)∂x]=T∂2η∂x2,
wherepj(j=1,2)   are the pressure for lower and upper fluids, respectively, and this pressure can be obtained using Bernoulli’s equation.Using Bernoulli’s equation for the pressure in (12) and subsequently linearizing it, we get
(13)ρ(2)(∂ϕ(2)∂t+U2∂ϕ(2)∂x)-ρ(1)(∂ϕ(1)∂t+U1∂ϕ(1)∂x)-2μ(2)∂2ϕ(2)∂y2+2μ(1)∂2ϕ(1)∂y2-E0[ε(2)∂ψ(2)∂x-ε(1)∂ψ(1)∂x]=T∂2η∂x2.


## 4. Normal Mode Analysis and Dispersion Relation

Now the normal mode technique has been used to find the solution of the governing equations. Let the interface elevation be represented by(14)η=Cexp[i(nx-ωt)]+c.c,
whereCis constant,kis the real wave number,ωis the growth rate, and c.c. refers to the complex conjugate of the preceding term.On solving (3) and (8) with the help of boundary conditions, we get
(15)ϕ(1)=i(ωn-U1)Csinhnycoshnaexp[i(nx-ωt)]+c.c,ϕ(2)=-i(ωn-U2)Cexp{-n(y-a)}×exp[i(nx-ωt)]+c.c,ψ(1)=iE0(ε(2)-ε(1))(ε(2)+ε(1)cothna)Csinhnysinhna×exp[i(nx-ωt)]+c.c,ψ(2)=iE0(ε(2)-ε(1))(ε(2)+ε(1)cothna)Cexp{-n(y-a)}×exp[i(nx-ωt)]+c.c.Substituting the values ofη, ϕ(1), ϕ(2), ψ(1), and ψ(2) in (13), we get the dispersion relation
(16)ρ(1)(ωn-U1)2tanhna+ρ(2)(ωn-U2)2-2iμ(1)(ωn-U1)ntanhna-2iμ(2)×(ωn-U2)n-E02(ε(2)-ε(1))2ε(2)+ε(1)cothna=Tn.Dividing (16) byρ(1)U12  and callingρ(2)/ρ(1)=ρ, μ(2)/μ(1)=μ, ε(2)/ε(1)=ε, σ=ω/nU1, k=na, andE2=E02ε(1)/ρ(1)U12, we can get
(17)(σ-1)2tanhk+ρ(σ-U)2-i2kRe×[(σ-1)tanhk+μ(σ-U)]=[E2(ε-1)2ε+cothk+Wk],
whereW=T/ρ(1)U12arepresents the Weber number and Re=ρ(1)U1a/μ(1) denotes Reynolds number.Equation (17) can be written as
(18)D(σ,k)=a0σ2+(a1+ib1)σ+(a2+ib2)=0,
where a0=tanhk+ρ, a1=-2(tanhk+ρU), b1=-(2k/Re)(tanhk+μ), a2=(tanhk+ρU2)-E2(ε-1)2/(ε+cothk)-Wk, and b2=(2k/Re)(tanhk+μU).Letσ=σR+iσΙ, separating the real and imaginary parts of (18), we have
(19)a0(σR2-σI2)+(a1σR-b1σI)+a2=0,σR=-a1σI+b22a0σI+b1.Eliminating the value ofσR from the previous equations, we get a quartic equation inσΙ  as
(20)A4σI4+A3σI3+A2σI2+A1σI+A0=0,
where A4=-4a03, A3=-8a02b1, A2=4a02a2-5a0b12-a0a12, A1=4a0a2b1-b13-a12b1, and A0=a0b22-a1b1b2+a2b12.From (20), we can get the value of maximum growth rateσIm  and corresponding wave number km. The neutral curves can be obtained from (20) by putting σI=0, so (20) reduces to
(21)A0=0,thatis,a0b22-a0b1b2-a2b12=0.Putting the values ofa0, a1, a2, b1, and b2 in the above equation, we get
(22)(1-U)2{μ2tanhk+ρtanh2k}=[E2(ε-1)2ε+cothk+Wk]{tanhk+μ}2.In practical applications,kis very small [4], then taking the approximationtanhk≈kin (22), we have
(23)(1-U)2{μ2+ρk}=[E2(ε-1)2εk+1+W]{k+μ}2.Hence, for instability it must be(24)(1-U)2W>[E2(ε-1)2W(εk+1)+1]{k+μ}2{μ2+ρk}.From the expression (24) the following conclusions can be easily drawn.(1)
If fluids are inviscid, that is,μ=0, it is reduced to dispersion relation as obtained by El-Sayed [15].(2)
In the absence of electric field and viscosity, (24) will represent the same dispersion relation as obtained by Rashed et al. [9].(3)
In the limit whenE, U, andμ tend to zero, (24) is reduced to dispersion relation as obtained by Squire [4].The condition of instability of the liquid sheet can also be written as(25)W<[(1-U)2{μ2+ρk}{k+μ}2-E2(ε-1)2(εk+1)].Hence, the instability occurs if Weber number is less than the critical value of Weber numberWc; otherwise system is stable whereWc  is given by
(26)Wc=[(1-U)2{μ2+ρk}{k+μ}2-E2(ε-1)2(εk+1)].


## 5. Results and Discussion

In this section the numerical computation has been carried out using the expressions (26) and (23) presented in the previous section. Water and air have been taken as working fluids identified with phase 1 and phase 2, respectively. Following parametric values have been taken.ρ
(
1
)
=
1.0
gm
/
c
m
3, ρ(2)=0.001gm/cm3, μ(1)=0.01poise, μ(2)=0.00018poise, ε(1)=80.37F/cm, ε(2)=1.0F/cm, and T=72.3dyne/cm.The neutral curves of Weber number for the different values of velocity ratio of two fluids have been drawn in Figure1. The region above the curves denotes the stable region, while below region is the region of instability. As velocity ratio increases, stable region decreases, and therefore, the velocity ratio of two fluids has destabilizing effect. In other words, one can say that the air velocity has destabilizing effect while liquid velocity plays stabilizing role.Figure 1
The neutral curves of Weber number for the different values of velocity ratio whenE=1, μ=0.018, ρ=0.0012.The effect of Weber number on the neutral curves of velocity ratio of two fluids has been shown in Figure2. Neutral curves for velocity ratio divide the plane into the stable region (below the curve) and unstable region (above the curve). As Weber number increases, the region of stability also increases for the same value of other parameters. The Weber number is directly proportional to the surface tension of the liquid. Therefore, as surface tension increases, the disintegration of the liquid sheet decreases. Therefore, it is concluded that the effect of surface tension resists the occurrence and development of instability.Figure 2
The neutral curves of velocity ratio for the different values of Weber number whenE=0, μ=0.018, ρ=0.0012.In Figure3, the effect of ratio of viscosity of two fluids has been studied. As viscosity ratio increases, the stable region decreases, and this shows that viscosity ratio of two fluids has destabilizing effect on the aerodynamic instability of the liquid sheet. The viscosity ratio is directly proportional to the air viscosity and inversely proportional to the liquid viscosity. Hence, air viscosity has destabilizing effect on the aerodynamic stability analysis of a liquid sheet, while liquid viscosity plays a stabilizing role. As the air viscosity increases, the air flow is impeded and instability occurs.Figure 3
The neutral curves of velocity ratio for the different values of viscosity ratio of two fluids whenE=5, W=100, ρ=0.0012.The effect of electric field intensityEon the neutral curves of velocity ratio has been shown in Figure 4 when Weber numberW=100  and viscosity ratio μ=0.018. It has been noticed that the increase in electric field intensity enlarges the stable region. Hence electric field has stabilizing effect on the instability of the liquid sheet. If electric field is present in the analysis, the term contributed from the applied electric field is added in the right hand side of (23) and so that critical value of velocity ratio increases. This shows that electric field has stabilizing influence.Figure 4
The neutral curves of velocity ratio for the different values of electric field intensity whenW=100, μ=0.018, ρ=0.0012.In Figure5, the variation of critical value of velocity ratioUwith the permittivity ratio of two fluidsεfor different values of electric field intensityE=10, 12, and 15  has been shown for Weber numberW=100. The figure shows that as the ratio of the permittivities of the two fluids increases, the critical value of relative velocity first decreases and then increases. It shows thatεplays dual role, that is, destabilizing as well as stabilizing effect in the stability analysis.Figure 5
The neutral curves of velocity ratio versus ratio of dielectric constants of two fluids for the different values of electric field intensity whenW=100, μ=0.018, ρ=0.0012.


## 6. Conclusions

The instability of a thin sheet of viscous and dielectric liquid moving in the same direction as an air stream in the presence of a uniform horizontal electric field has been investigated using viscous potential flow theory. Both fluids are taken as incompressible and viscous with different kinematic viscosities. A dispersion relation has been derived, and stability is discussed theoretically as well as numerically. Stability criterion is given in the terms of velocity ratio of two fluids as well as critical value of applied electric field. The system is unstable when the Weber number is less than the critical value of Weber number; otherwise, it is stable. It has been found that aerodynamic-enhanced instability occurs if the velocity ratio is greater than a critical value related to the ratio of the air and liquid viscosities, Weber number, the electric field, and the dielectric constant values. The electric field is found to have a stabilizing effect while dielectric constant plays dual role in the stability criterion. Surface tension has stability effect on the aerodynamic instability of a liquid sheet. Liquid viscosity has stabilizing effect in the stability analysis, while air viscosity has destabilizing effect.



---

*Source: 528723-2013-08-29.xml*