Modelling and Simulation in Engineering

Volume 2013 (2013), Article ID 528723, 6 pages

http://dx.doi.org/10.1155/2013/528723

Research Article

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

Mukesh Kumar Awasthi,¹ Vineet K. Srivastava,² and M. Tamsir³

¹Department of Mathematics, University of Petroleum and Energy Studies, Dehradun 248007, India

²ISRO Telemetry, Tracking and Command Network (ISTRAC), Bangalore 560058, India

³Department of Mathematics, Graphic Era University, Dehradun 248002, India

Received 19 May 2013; Accepted 31 July 2013

Academic Editor: Abdelali El Aroudi

Copyright © 2013 Mukesh Kumar Awasthi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.

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 sheet disintegration, 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 impact disintegration, 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  , viscosity  , dielectric constant, surface tension, and thickness, moving with velocityalong an air stream (on both sides of the sheet), which is of density, viscosity, dielectric constant, and moving with velocity. 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  . On applying small disturbances to the system, the equation of upper surface can be written as

whererepresents the displacement of liquid sheet from its equilibrium position and represents the equation of the disturbed surface.

The unit outward normal to the interface is given by

and the unit outward normal of the first order is given by

In each fluid layer velocity is expressed as the gradient of the potential functionand the potential functions satisfy Laplace’s equation, that is,

where.

In initial state, the potential function can be written as

The considered system is influenced by a uniform horizontal electric fieldin the positive -direction

We have assumed that the quasistatic approximation is valid for the problem, and therefore, the electric field can be derived from electric scalar potential functionsuch that

Gauss’s law requires that the electric potentials also satisfy Laplace’s equation, that is,

3. Boundary Conditions

(i)Kinematic condition that every particle on the interface remains on the interface leads to the following boundary conditions:

(ii)The tangential component of the electric field must be continuous across the interface, that is, or

(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 by   or

(iv)Interfacial condition for conservation of momentum is

where   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

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

whereis constant,is the real wave number,is the growth rate, and . refers to the complex conjugate of the preceding term.

On solving (3) and (8) with the help of boundary conditions, we get

Substituting the values of, , , , and in (13), we get the dispersion relation

Dividing (16) by  and calling, , , , , and, we can get

whererepresents the Weber number and denotes Reynolds number.

Equation (17) can be written as

where , , , , and .

Let, separating the real and imaginary parts of (18), we have

Eliminating the value of from the previous equations, we get a quartic equation in  as

where , , , , and .

From (20), we can get the value of maximum growth rate  and corresponding wave number . The neutral curves can be obtained from (20) by putting , so (20) reduces to

Putting the values of , , , , and in the above equation, we get

In practical applications,is very small [4], then taking the approximationin (22), we have

Hence, for instability it must be

From the expression (24) the following conclusions can be easily drawn.

(1)If fluids are inviscid, that is, , 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 when , , 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

Hence, the instability occurs if Weber number is less than the critical value of Weber number ; otherwise system is stable where  is given by

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.

, , , , , , and .

The neutral curves of Weber number for the different values of velocity ratio of two fluids have been drawn in Figure 1. 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 when, , .

The effect of Weber number on the neutral curves of velocity ratio of two fluids has been shown in Figure 2. 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 when, , .

In Figure 3, 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 when, , .

The effect of electric field intensityon the neutral curves of velocity ratio has been shown in Figure 4 when Weber number  and viscosity ratio . 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 when, , .

In Figure 5, the variation of critical value of velocity ratiowith the permittivity ratio of two fluidsfor different values of electric field intensity, , and   has been shown for Weber number. 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 thatplays 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 when, , .

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.

References

A. H. Lefebvre, Atomization and Sprays, Hemisphere, New York, NY, USA, 1989.

K. Masters, Spray Drying Handbook, Wiley, New York, NY, USA, 4th edition, 1985.

S. P. Lin and W. C. Liu, “Instability of film coating of wires and tubes,” AIChE Journal, vol. 21, no. 4, pp. 775–782, 1975.

H. B. Squire, “Investigation of the instability of a moving liquid film,” British Journal of Applied Physics, vol. 4, no. 6, article no. 302, pp. 167–169, 1953.

W. W. Hagerty and J. F. Shea, “A study of the stability of plane fluid sheets,” Journal of Applied Mechanics, vol. 22, p. 509, 1955.

R. P. Fraser, “The fluid kinetics of applications of pesticidal chemicals,” Advances in Pest Control Research, vol. 11, pp. 1–106, 1958.

N. Dombrowski and W. R. Johns, “The aerodynamic instability and disintegration of viscous liquid sheets,” Chemical Engineering Science, vol. 18, no. 3, pp. 203–214, 1963.

J. G. H. Joosten, A. Vrij, and H. M. Fijnaut, in Physicochemical Hydrodynamics, Proceedings of the International Conference on Physical Chemistry and Hydrodynamics, D. B. Spalding, Ed., vol. II, p. 639, Advance Publications, Guernsey, UK, 1978.

M. I. I. Rashed, M. A. Ghazi, and H. Elbanna, “Aerodynamic instability of a liquid sheet sprayed with an air stream,” Journal of Physics D, vol. 12, no. 10, article 009, pp. 1679–1684, 1979.

E. A. Ibrahim and S. L. Jackson, “Spatial instability of a liquid sheet in a compressible gas,” Journal of Colloid and Interface Science, vol. 180, no. 2, pp. 629–631, 1996.

E. A. Ibrahim and E. T. Akpan, “Liquid sheet instability,” Acta Mechanica, vol. 131, no. 3-4, pp. 153–167, 1998.

E. A. Ibrahim, “Instability of a liquid sheet of parabolic velocity profile,” Physics of Fluids, vol. 10, no. 4, pp. 1034–1036, 1998.

A. A. Ibrahim and M. A. Jog, “Nonlinear instability of an annular liquid sheet exposed to gas flow,” International Journal of Multiphase Flow, vol. 34, no. 7, pp. 647–664, 2008.

J. R. Melcher, Continuum Electromechanics, MIT Press, Cambridge, Mass, USA, 1981.

M. F. El-Sayed, “Electro-aerodynamic instability of a thin dielectric liquid sheet sprayed with an air stream,” Physical Review E, vol. 60, no. 6B, pp. 7588–7591, 1999.

D. D. Joseph and T. Y. Liao, “Potential flows of viscous and viscoelastic fluids,” Journal of Fluid Mechanics, vol. 265, pp. 1–23, 1994.

D. D. Joseph, J. Belanger, and G. S. Beavers, “Breakup of a liquid drop suddenly exposed to a high-speed airstream,” International Journal of Multiphase Flow, vol. 25, no. 6-7, pp. 1263–1303, 1999.

T. Funada and D. D. Joseph, “Viscous potential flow analysis of Kelvin-Helmholtz instability in a channel,” Journal of Fluid Mechanics, vol. 445, pp. 261–283, 2001.

M. K. Awasthi and G. S. Agrawal, “Viscous Potential flow analysis of Kelvin-Helmholtz instability of cylindrical interface,” International Journal of Applied Mathematics and Computation, vol. 3, no. 2, pp. 131–138, 2011.

A. M. Ardekani and D. D. Joseph, “Instability of stationary liquid sheets,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 13, pp. 4992–4996, 2009.