The Instability of an Electrohydrodynamic Viscous Liquid Micro-Cylinder Buried in a Porous Medium: Effect of Thermosolutal Marangoni Convection

The electrohydrodynamic (EHD) thermosolutal Marangoni convection of viscous liquid, in the presence of an axial electric field through a micro cylindrical porous flow, is considered. It is assumed that the surface tension varies linearly with both temperature and concentration. The instability of the interface is investigated for the free surface of the fluid. The expression of the free surface function is derived taking into account the independence of the surface tension of the heat and mass transfer. The transcendental dispersion relation is obtained considering the dependence of the surface tension on the heat and mass transfer. Numerical estimations for the roots of the transcendental dispersion relation are obtained indicating the relation between the disturbance growth rate and the variation of the wave number. It is found that increasing both the temperature and concentration at the axial microcylinder has a destabilizing effect on the interface, according to the reduction of the surface tension. The existence of the porous structure restricts the flow and hence has a stabilizing effect. Also, the axial electric field has a stabilizing effect. Some of previous analytical and experimental results are recovered upon appropriate data choices.


Introduction
The fluid flow in microjet has various applications in medical, biomedical, computer chips and chemical separations.Advent of microelectrical mechanical systems is one of the major advances in industrial technologies.Also, microjets are fundamental to connecting different devices, utilized in biochemical reaction chamber, in physical particle separation, in inkjet print heads, in infrared detectors, in diode lasers, in miniature gas chromatographs, and in heat exchanger for cooling computer chips.Understanding the flow characteristics through a microjet is very important in determining the free surface instability of the fluid flow.The augmentation of heat transfer from a solid cylinder wrapped with a porous layer was considered by Bhattacharyya and Singh [1].An experimental analysis of unsteady heat and moisture transfer around a heated cylinder buried into a porous medium was performed by Moya et al. [2].Furlani [3] studied the temporal instability of an infinite Newtonian cylindrical microjet (as a cylindrical fluid surface) that is subjected to a sinusoidal variation of surface tension along its length.Furlani and Hanchak [4] developed the same problem to study numerically the nonlinear analysis of the deformation and breakup of viscous microjets.The instability of non-Newtonian liquid jets (viscoelastic and power law) was investigated by Gao [5] and Gao and Ng [6].The dispersion relation between the growth rate and the wave number for a non-Newtonian cylindrical liquid jets was derived and the instability behavior of viscoelastic jets was investigated.
The interaction of electric fields with moving fluids makes EHD a very complicated phenomenon.The process of EHD is dependent on many parameters and properties of the liquid and environment of the flow.On the other hand, EHD can be considered as the branch of fluid mechanics concerned with the electrical force effects or as the part that is involved in the influence of moving media on electric fields.It is thus concerned with the interaction between electrical fields and free or polarized charges in fluids.In the EHD stability, the electric field plays an important role in many practical problems of chemical engineering, biophysics and many other fields.Mohamed and Nayyar [7] have investigated the stability of a cylindrical jet of incompressible inviscid liquid in the presence of an axial electric field.They showed that the uniform axial electric field has a strong stabilizing influence on the cylindrical interface for short and long wavelengths in all symmetric and asymmetric modes of perturbation.
Marangoni boundary layers are dissipative layers which may occur along liquid-liquid or liquid-gas interfaces.When a free liquid surface is present, the surface tension variation is resulting from the temperature gradients along the surface.It can also induce a motion within the fluid which is called thermocapillary flow (thermal Marangoni convection) [8].The surface tension gradients that are responsible for Marangoni convection depend on both temperature and/or concentration gradients.Earliest work in this field was first investigated by Napolitano [9,10].Surfactants or surface contaminants are used to control the stability of the free surface by lowering surface tension.Variations in surface tension that are introduced by surfactant can substantially alter interfacial evolution and flow.The effect of surfactants on liquid has been studied by many authors [11][12][13][14][15]. McGough and Basaran [16] studied the breakup mechanisms of the fluid threads containing surfactant by solving the Navier-Stokes equations.The effect of the Prandtl number on Marangoni convection over a flat plate was studied by Christopher and Wang [17].Physically, the heat and mass transfer across the interface play an essential role in determining the flow field in some cases.For instance, when the fluid is boiling, whether it is a film or bubble, the motion of the film and bubbles is mainly dependent on the effect of mass and heat transfer.In chemical industry, the interfacial stability problems are important in the presence of mass and heat transfer.The mechanism of heat and mass transfer across an interface is of a great importance in numerous industrial and environmental problem processes.These include the design of many types of contacting equipment, for example, boilers, condensers, evaporators, gas absorbers, pipelines, chemical reactors, and nuclear reactors.The effect of heat and mass transfer in the absence of gravity on the surface tension was studied by Straub [18].
The flow through porous media is usually described by the Darcy's law that relates the movement of the fluid to the pressure gradient acting on a parcel of the fluid.Meanwhile, the Darcy's equations are not the fundamental equations for the flow through porous media.It is an approximation for the balance of linear momentum for the fluid flow through a porous solid within the context of mixtures.For example, the interactive force between the fluid and the porous medium, the frictional effects due to viscosity, and the inertial nonlinearities are not included in the Darcy's model [19].These effects can be included in more generality model like Brinkman-Darcy's model.The basic assumptions that lead to the Brinkman-Darcy equation were illustrated by Rajagopal [19].Elcoot and Moatimid [20] studied the instability of finitely conducting cylindrical flows through porous media, under the influence of an axial electrostatic field.They found that, under certain conditions, the field may have a stabilizing or destabilizing effect.An experimental study of the steam injection into a porous media and the stability of the interface between the steam and the water was carried out by Catton and Chung [21].
The aim of this work is to discuss the effect of axial electric field with the existence of the porous structure on the stability of the microcylindrical flow of a viscous liquid past a microsolid cylindrical surface.Also, the effect of the Marangoni thermosolutal convection (according to the dependence of the surface tension on the heat and mass transfer) is discussed.The main construction of the problem is as follows.Firstly, we studied the problem according to spatial surface tension.So, the expression of the free surface function is obtained.Secondly, the problem is studied in the presence of the heat and mass transfer rates with thermosolutal surface tension.The dispersion relation, in the latter case, is obtained as a transcendental function of the growth rate and the wave number.Finally, we studied the effect of the porous medium and the axial electric field on the stability of viscous liquid interface and tabulated results are obtained.
To clarify the problem, in Section 2, the physical description of the problem including the basic equations that govern the motion is presented.Section 3 is devoted to introducing the interfacial conditions for both the fluid and the electric field together with the boundary conditions.The method of solution according to the normal modes technique is presented in Section 4. In Section 5, the instability analysis in the absence of the heat and mass transfer is studied and the expression of the free surface function is derived.The transcendental dispersion relation in the presence of the heat and mass transfer is obtained in Section 6.Throughout Section 7, we introduced numerical estimations and gave a discussion for the stability picture according to the relation between the growth rate and the wave number.Finally, in Section 8, we give concluding remarks for this study based on the obtained results of the stability analysis.

Formulation of the Problem
Consider a microcylindrical layer of viscous fluid with density , viscosity , and radius  0 flowing past an axial solid cylinder of infinitely small radius.The flow is surrounded by a vacuum (which has zero density and zero viscosity) with hypothetical sufficiently large radius  2 .The axial microsolid cylinder is elongated along the -axis.It is earthed and has a radius  1 (where  1 ≪  2 ).The axial microcylindrical solid rod, which induces temperature and concentration gradients, is considered to be at a constant temperature  1 and constant concentration  1 .The subscripts (1) and (2) refer to the liquid and the surrounding vacuum, respectively.The fluid and the vacuum are influenced by an external uniform electric field of intensity  0 which acts along the positive -axis.The dielectric constants for the liquid and the vacuum are  1 and  2 , respectively.The interface between the fluid and the surrounding vacuum is considered to be free of charges.The flow jet is considered as a microcylindrical porous structure with Darcy's coefficient ], while its porosity is assumed to be unity (this will be explained physically later).Also, we neglect the effects of gravity.Considering the previous hypotheses, the model may be illustrated graphically as in Figure 1.After making a small disturbance, the interface function becomes and the unit outward normal vector n and the unit tangential one t are given by where r and ẑ are the unit vectors along the radial and axial cylindrical coordinate directions and ℎ  = ℎ/.According to the assumption of the viscosity of the fluid, considering the existence of the porous structure and the heat and mass transfer, the basic equations of the problem may take the following form [3,19].
The continuity equation (incompressibility condition) requires The balance of linear momentum (Navier-Stokes) gives The balance of energy is And finally the concentration equation becomes In the previous equations ( 3)-( 6),  is the porosity of the porous medium,  eff is the effective viscosity of the fluid ( eff = / and if  = 1 then  eff =  as we mention in the next paragraph), and (/) = (/)+(V⋅∇) V. Also,  is the thermal diffusivity,  is the mass diffusivity,  is the temperature,  is the concentration,  is the pressure, V is the velocity vector, ] is the Darcy's coefficient, and  is the density.
The basic assumptions that lead to the Brinkman-Darcy equation and the form of the basic assumptions that lead to the Brinkman-Darcy equation were illustrated by Rajagopal [19], and can be summarized in the following points.
(1) The porous medium is solid and thus the balance of linear momentum of the porous medium can be ignored.
(2) The interactive force between the fluid and the porous medium is due to the frictional forces only and this force is proportional to the flow velocity which is represented by the term ]V, where ] = /,  is the fluid viscosity and  is the permeability of the porous medium.
(3) The frictional effects due to viscosity were taken into account by the term  eff ∇ 2 V, where  eff = /.
(4) The flow is unsteady and sufficiently fast, so that the inertial nonlinearities cannot be ignored; thus the term (1/ 2 )(V ⋅ ∇)V needs to be retained.
According to the previous assumptions, the balance of linear momentum can be written as in (4).Also, we want to confirm the following points.
(1) The velocity of the flow through the porous medium is the average of the fluid velocity.This quantity has been given various names by different authors, such as seepage velocity, filtration velocity, superficial velocity, or Darcy velocity [22].The seepage velocity V is related by the velocity of the fluid  (where  is the flow velocity in the absence of the porous medium) by V = .
(2) The porosity  for natural media does not normally exceed 0.6.For beds of solid spheres of uniform diameter  can vary between the limits 0.2595 (rhombohedral packing) and 0.4764 (cubic packing).Nonuniformity of grain size tends to lead to smaller porosities than those for uniform grains, because smaller grains fill the pores formed by larger grains.For man-made materials such as metallic foams  can approach the value 1 (see the tabulated values for the porosity and permeability in [22, page 5]).
(3) The term (1/ 2 )(V ⋅ ∇)V is inappropriate if the flow is sufficiently slow.Then the inertial nonlinearities can be neglected [19].However, this term needs to be retained in the case of highly porous media where the flow is sufficiently fast and at least the irrotational part of this term needs to be retained [22] and for incompressible fluids this term reduces to and for highly porous media we can assume  = 1 (as we assume in the present problem).( 4) Brinkman sets the viscosity  and the effective viscosity  eff equal to each other but in general that is not true [22].The difference between  and  eff is due to the momentum dispersion and  eff is dependent on the type of the porous medium as well as the strength of flow, and it is common practice to take  eff equal to  for high porosity media [23].
(5) Experimental checks of Brinkman's theory have been indirect and few in number.The Brinkman's equation reduces to Darcy's equation when the spatial length scale is much greater than ( eff /) 1/2 .Therefore the term ∇ 2 V is negligible in comparison with the term proportional to V [22].But, the Brinkman's model holds only for particles whose size is of order  3 , where  ≪ 1 is the distance between neighboring particles, and for large particles the fluid filtration is governed by Darcy's law [22].
Because of the quasi-static approximation and in poorly conducting fluid [24,25], the induced magnetic field is negligible.Also, we deal with very small electric conductivity and there is no applied magnetic field.So, the electric field is conservative.Since, the electric field is curl free vector the Maxwell's equations are ∇× = 0 and ∇⋅ = 0. So, the electric field has an electric scalar potential  such that  =  0 ẑ−∇.Therefore the electric potential   satisfies the Laplace's equation The boundary and interfacial conditions for the problem may be displayed in the following section.

Boundary Conditions
The boundary conditions adopted here must be classified into three categories.The first relates the hydrodynamic part at the boundaries.The second is concerned with the electric part at the boundaries.Meanwhile, the third relates the combination condition of the electric and hydrodynamic balance of the stress tensor components at the interface between the liquid and the surrounding vacuum.
(1) For the hydrodynamic part, because  (in (1)) is a scalar function that is always equal to zero at any point on the fluid interface, its time derivative following any material point on the interface is obviously equal to zero, which means that there is no phase transformation occurring [26, Chapter 2, page 75].So, the continuity of the normal velocity at the interface requires where  is the radial distance and ℎ is the radius for the disturbed fluid surface.
According to the microscale of the cylindrical fluid radius, the standard no-slip boundary condition in the classical fluid mechanics does not apply, so that the shear stress at the boundary vanishes (i.e., there is no considerable surface fraction between the fluid and the axial bar) and hence V  / = 0 at the axial bar.Meanwhile, the normal velocity must be vanishing at the axial bar.So the boundary conditions for the velocity at the surface of the microcylindrical axis may be written as follows: where  1 is the microradius of the axial microcylindrical bar such that  1 is sufficiently small ( 1 ≪ 1) and V  , V  are the radial and axial velocity components, respectively.
(2) For the electric part, because there is no surface charges accumulated at the interface, the normal electric displacement must be continuous across the interface.So that  ⋅ ‖‖ = 0 and the interfacial condition for the normal electric field displacement, in the linear form, yields Since the axial microsolid bar and the vacuum at infinity have no electric field, the electric potential at the boundaries satisfies (3) The balance of the stress components at the interface may be presented as follows.
At the interface between the liquid and vacuum, the fluid and the electrical stresses must be balanced.The components of these stresses consist of the electric hydrodynamic stresses together with the surface tension force [27].The combination of the electric stress  electro  and hydrodynamic stress  hydro  parts may be presented in the total stress as The electric force density   can be written in terms of the electric field as   = (∇ ⋅ ).Since ∇ ×  = 0, the electric force can be written as   = (∇ ⋅ ) + (∇ × ) ×  and the Maxwell stress tensor  electro  , in the index notation, may be written as [28]  electro meanwhile, the hydrodynamic stress tensor may be described by the constitutive relation where Δ  is the usual kronecker delta and V  , V  ,   ,   are a general notations for the velocity components and the coordinates, respectively.According to the simple interface description, which involves only interfacial tension, the forces acting on any segment of an interface are of two kinds.First, there are the bulk pressure and stresses that act on the faces of the interface element and produce a net effect that is proportional to the surface area; second, there is a tensile force that is due to surface or interfacial tension that acts in the plane of the interface at the edges of the surface element and is specified by means of the magnitude of the surface or interfacial tension as a force per unit length [26, Chapter 2, page 76].So, the stress balance can be written as [3] where  is the curvature of the interface, ∇  is the surface gradient at the interface, and () is the surface tension.
The previous condition can be decomposed into the normal and tangential stress tensor components as follows.
(i) The normal stress component requires (ii) The tangential stress components yield where the surface curvature and the gradient at the interface may be given as follows: Finally, the normal and tangential stress tensor conditions as given in ( 17) and (18) (at  = ℎ) may be written as follows: We study the solution of the equations of motion (4) according to the boundary and interfacial conditions in two cases.The first case is obtained in the absence of the heat and mass transfer where the spatial surface tension is independent of the rate of heat and mass transfer.In this case, we follow the same analysis of Furlani [3] to gain the function of the free surface of the microcylindrical fluid flow.In the second one, taking into account the dependence of the surface tension on the heat and mass transfer, the surface tension becomes a function of the time and the axial distance.So, in this case, the dispersion relation may be obtained as a transcendental function and its roots can be computed numerically hence the stability picture may be obtained graphically.Before studying the previous cases, we obtain the solution of the electric field according to the boundary and interfacial conditions.

Solution of the Eclectic Field
Our analysis will be based on the normal modes technique.Therefore, the electric potential may take the following form: where  is the frequency of the surface wave and  is its wave number which is assumed to be real and positive.Using (7) with the boundary conditions ( 10)-( 12), we may obtain the solution for the electric potential functions as follows: where  0 and  0 are the modified Bessel's functions of first and second kinds of order zero, respectively.The dashes denote the differentiation with respect to  and  = √ −1.Also, the functions  1 (),  2 (),  1 (), and  2 () are given in the appendix.

Instability in the Absence of Heat and Mass Transfer
This section is devoted to solving the equations of motion and the continuity equation (see ( 3) and ( 4)) according to the interfacial and boundary conditions (see (8) and ( 9)) with employing the normal and tangential stress tensor interfacial conditions.We follow the same procedure given by Furlani [3] to get the function that describes the free surface.So, we expand the velocity components and the pressure as a power of  (depending on the microradius of the flow jet), as follows: From the continuity equation (3), one gets where V 0 (, ), V 2 (, ) are perturbed velocity functions and  0 (, ),  2 (, ) are the perturbed pressure functions.It is important to notice here that the expansions in ( 24)-( 25) are compatible with the boundary conditions in (9) and with the continuity equation (3).Where  1 and  are a microradii of the axial microcylinder and the microporous cylinder that contains the fluid, respectively.So, all terms containing () and ( 1 ) or higher orders can be neglected according to their tendency to zero.
According to the previous expansions, the equation of motion (4) for V  can be neglected because it is identical to the lowest order of .Therefore the equation of motion for V  can be written as The second order term of the velocity V 2 can be determined from the tangential stress condition as given in (21); on neglecting the terms of () and all higher orders, it can be written as follows: The zero order term of the pressure  0 can be determined from the normal stress condition, as given in (20), in the form Now, substituting from ( 27) and ( 28) into ( 26), we obtain The kinematic condition equation (8) gives Also, the axial velocity V 0 (, ), the radial position ℎ(, ), and the surface tension () can be represented in a perturbed and unperturbed parts as follows: ℎ (, ) =  0 (1 +  (, )) , where  0 ,  0 , and  0 are the constant unperturbed axial streaming velocity, radius, and surface tension, respectively.The analysis may be simplified by using the transformation By using ( 31)-(34) into ( 29) and (30), after linearizing the resulting equations, we can obtain Equations ( 35) can be combined into one equation that determines the free surface function (, ) as follows: (36) Initially, the disturbed fluid radius (, ) is independent of the time (also its partial derivative with respect to the time).So, the disturbed fluid radius (, ) has the following initial conditions: Assuming that the surface tension  1 () has a spatially periodic form, which is independent of the heat and mass transfer.Then, the surface tension function may be written in the following form: where Δ is a constant variation of the surface tension.
The initial value problem given by (36) and the initial conditions equation (37) can be solved by decomposing the solution into a homogenous and particular parts,  ℎ (, ) and   (, ), respectively.Based on the form of the surface tension  1 (), which is given in (38), the homogenous solution of the wave function  ℎ (, ) may be written in the form The final form of the general solution of the free surface function, (, ) =  ℎ (, ) +   (, ), can be written as follows: where, the roots  ± and  0 are defined in the appendix.
It is worthwhile to mention here that the final solution in (40) reduces to the same solution obtained earlier by Furlani [3] for the flow of cylindrical microjet surface in the absence of the electric field ( 0 = 0).It is worthwhile to mention here that the final solution in (40) reduces to the same solution that was obtained earlier by Furlani [3] for the flow of cylindrical microjet surface in the absence of the electric field ( 0 = 0) with neglecting the porous structure (] = 0) and with taking the real part only of the free surface function.

Instability of Thermosolutal Marangoni Convection
If the two-fluid system is a multicomponent system, it is often the case that there may be a preferential concentration of one or more of the components at the interface (e.g., if we consider a system of pure  and pure , which are immiscible, with a third solute component  that is soluble in  and/or  but that is preferentially attracted to the interface), and then the interfacial tension will also be a function of the (surfaceexcess) concentration of these solute components.Both the temperature and the concentrations of adsorbed species can be functions of position on the interface, thus leading to spatial gradients of the surface tension [26].So, to study the effect of the heat and mass transfer on the interfacial instability of the microcylindrical surface, the surface tension is assumed to be a function of the temperature and the solute concentration and the surface tension will depend on temperature and concentration linearly [29], as follows: where   and   denote the temperature and concentration coefficients of the surface tension.Also,  0 and  0 are the reference temperature and concentration, respectively.Note that (41) is the same as (33) with  1 () being  1 (, ) that depends on the heat and mass transfer as in (41). 0 is the unperturbed constant value of the surface tension.
To determine the form of the surface tension, we need firstly to obtain the solution of the heat and concentration equations.Again, the analysis will be based on the normal modes technique as given by  (, , ) =   ()  + ,  (, , ) =   ()  + .
If we consider the diffusion of the thermal energy and the concentration in an infinite isolated column of stationary flow of radius  0 then there are no temperature and concentration gradients at  =  0 .So, the temperature and concentration boundary conditions are as follows [3]: Mathematical Problems in Engineering Equations ( 43) and (44) represent the standard Bessel's differential equation.Therefore, the general solution of the temperature and concentration function may be written as follows: (, , ) On combining the following equations: (1) the general solution of the heat and mass transfer (see (46)), (2) the surface tension expression (see (41)), (3) the expression of the free surface function, which may be given based on the normal modes technique as (, ) =  + , into the free surface deferential equation given by (36).
We finally get the transcendental dispersion relation of the instability of the interfacial free surface as where the functions  1 () and  1 (, ) are defined in the appendix.
It is convenient to write the stability criterion in an appropriate dimensionless form.This can be done in a number of ways depending primarily on the choice of the characteristic length, time, temperature, and mass.Consider the following dimensionless forms depending on the characteristic length =  0 , the characteristic time = 1/ ω, the characteristic temperature =  1 , and the characteristic mass =  0 / ω2 , where ω is a characteristic value of .The other dimensionless quantities are given by The superscript asterisks refer to the dimensionless quantities.From now on, these will be omitted for simplicity.The dispersion relation in (47) and its coefficients  1 () and  1 (, ) may be written as follows: where the functions () and (, ) are defined in the appendix.
The dimensionless transcendental dispersion relation in (49) describes the instability of the interfacial free surface of the microfluid jet.We will discuss the instability criteria numerically and graphically in the next section.

Numerical Estimations
In the case of the thermosolutal Marangoni convection, the surface tension depends on the heat and mass transfer and the dispersion relation is the transcendental equation (49).The stability criteria of the system may be discussed by getting the values of the growth rate  from the roots of the transcendental equation (49).The positive values of  indicate a disturbance that grows with time (instability) and the negative values indicate decays with time (stability).So, to study the instability of the microcylindrical surface, we indicate the relation between the growth rate  and the wave number  graphically.In all coming curves, the positive values region of  indicates the instability region, denoted by letter .Meanwhile, the negative values region of  indicates the stability region, denoted by letter .
It is well known that increasing the wave number enhances the stability of the system.As seen in Figure 2, the stability occurs at  ≅ 0.7, when the intensity of the electric field  0 = 1000.Meanwhile, at  0 = 1 the stability occurs at  ≅ 1.This means that the stability occurs at less value of  and at the higher value of  0 , which means that the system is more stable at highly electric field intensity and the electric field has a stabilizing effect.This is an early result obtained by many researchers [7,20,24,25].Also, we can notice that the maximum instability in all figures, except Figure 2, according to the peak of the curve occurs at  ≅ 0.6.But at high values of the electric field, the peak of the curve approaches the stability region and the maximum instability occurs at  ≅ 0.4.This implies that existence of the electric field causes a high stability of the interface.
Figure 3 illustrates the electric field and the growth rate relation.It is clear from this figure that the maximum value of the electric field  0 is congruent to the zero value of the growth rate .The zero value the growth rate  corresponds to the steady state of the interface.This result confirms the same result of Figure 2; that is, the electric field has a stabilizing effect.Figure 4 illustrates the effect of the dielectric constants difference ε =  1 −  2 on the stability picture.The dielectric difference ε has a stabilizing effect, where the higher value of ε is identical to the minimum peak value of the curve.
Figures 5-7 illustrate the effect of Darcy's coefficient, the viscosity, and the density of fluid on the stability behavior.The of Darcy's coefficient, the viscosity, and the density of fluid have stabilizing effect for  ≤ 1 according to reducing the peak of the curve and decreasing the number of the points in the instability region.Figure 5 illustrates the effect of the existence of the porous structure, according to the different values of the Darcy's coefficient, on the stability picture.It is clear that the microcylindrical surface is more stable, in case of the flow through a porous medium for majority values of the wave number, due to an increase in the values of the Darcy's coefficients in comparison with the case of the pure flow (] = 0).This is due to the fact that the greatest peak value of the curve in the instability region is at ] = 0.
Hence the Darcy's coefficient ] is equal to /, where  is the permeability of the porous medium, so an increase in the values of the Darcy's coefficients corresponds to a decrease in the permeability of the porous medium.This in turn restricts the streaming velocity of the fluid flow.It is known that when the streaming velocity decreases, the stability of the system increases.In other words, when the permeability of the medium increases ( increases and ] decreases) the holes of the porous medium are very large and the resistance of the medium may be neglected so that the streaming velocity  increases and causes instability of the system.The same result was illustrated experimentally by Catton and Chung [21], where the authors studied the interfacial stability between the water and the steam.Two porous media were studied experimentally: a bead of glass with diameter 0.6, 2, and 6 and with permeability 0.14, 1.57, and 18.0 (×10 −9 m 2 ), respectively, and a bead of steel with diameter 4 and with permeability 9.02 (×10 −9 m 2 ).It was found that in the more permeable cylinder the steam bubble is chaotic just after it has been formed.This means that the interface becomes unstable and the chaotic turbulence becomes more and more profound as the water flow rate increases (according to increasing the permeability of the porous medium).The viscosity  and the density  have a stabilizing effect, for larger values of the wave number, as shown in Figures 6  and 7.This influence occurs due to an increase in the viscosity damping and in the inertia.So, the interface of the microcylindrical surface becomes stable with an increase in the viscosity and the density.
The effect of heat transfer on the free surface is illustrated in Figure 8.To study the instability influence, according to the heat transfer, we notice that the temperature coefficient is   = ±(1/ 0 )(/) [30], where   has a positive or negative sign according to the increasing or decreasing of the surface tension.When  <  1 , the axial microsolid cylinder, which induces the heat and mass gradients, is heated and then the temperate transfers from it to the free surface by convection currents.This in turn tends to evaporate the fluid easily and provokes higher instabilities at the interface.So, in this case,   has a negative value according to the decreasing of the surface tension.Conversely,   has a positive value when  >  1 .In this case the free surface loses its temperature that is transferred to the core region of the flow jet and hence the free surface becomes stable.This phenomenon was illustrated in Figure 8 where the temperature coefficient   has a stabilizing effect when it is positive and the inverse occurs at   < 0.
In other words, heating the axial microsolid surface  1 has destabilizing effect at the interface and conversely for cooling.growth rate in the stability region.This result was obtained experimentally by D' Aubeterre et al. [31].The authors stated that the temperature has a destabilizing effect due to the fact that a substance is more volatile for increasing the temperate; this means that when a temperature gradient occurs, the substance tends to evaporate easily and provokes higher instabilities at the interface.Also, depending on the alcohol studied [31], turbulence begins at the interface at different temperature gradients.Instabilities begin near the interface, showing a sinuous movement.At higher temperatures and temperature gradients, instabilities increase and movements become continuous.The most evident effect occurs when temperature gradients are 15-17 ∘ C; at these values a sinuous wave was observed [31].Similar behavior can be noticed for the concentration as seen in Figure 10.In this figure we use   = ±(1/ 0 )(/) [30].The concentration coefficients of the surface tension   become negative if  <  1 .This implies that the concentration at the interface is less than the concentration at the core region.Because the surface tension is proportional with the concentration, the surface tension at the free surface accordingly reduces.So, the free surface loses its interfacial rigidity and its ability to movement increases.Therefore the instability of the free surface occurs.On the other hand, when the concentration of the interface is higher than that of the core region, the stability in the interface occurs.This means an increase of the interfacial rigidity, which in turn restricts any surface movements or Marangoni convection.Therefore, the surface becomes stable when   > 0 and the instability occurs for   < 0, as shown in Figure 10. Figure 11 confirms this result, where the concentration is plotted versus the growth rate.The region of the positive values of the growth rate , that is, the unstable region, occurs at  <  1 .Meanwhile, the stable region, that is, the negative growth rate values, occurs at  >  1 .Similar results were obtained by Agble and Mendes [32] and compared with experimental results.Finally, the instability criterion of the system can be discussed by studying the roots of the dispersion relation (49) to get the values of the growth rate , as shown in Table 1.The negative values of the roots of  indicate stability of the microcylindrical surface and the inverse occurs of the positive root values.It can be noticed from the roots of the growth rate  that the interface becomes stable when the microradius  0 increases which means an increase in the inertia with  0 .

Conclusion
In this paper, we have examined the influence of the existence of the porous structure on the EHD instability of the free microcylindrical fluid surface with thermosolutal Marangoni convection.The surface tension depends linearly on the heat and mass transfer.The main results of our study can be epitomized in the following points.
(1) The existence of the porous structure has a stabilizing effect for the free surface.Existence of the porous medium restricts the streaming velocity of the fluid.Similar behavior was obtained experimentally by [21] for the steam-liquid interface thorough a porous structure.
(2) As it is known, the electric field has a stabilizing effect, and this result coincides with the previously published results in the literature [7,20,24,25].
(3) The heating of the axial microsolid cylinder reduces the surface tension and makes the free surface volatile easily and hence the surface becomes unstable and the inverse occurs for cooling the axial microsolid cylinder.The same results were obtained experimentally by [31].
(4) The surface tension is proportional to the concentration variation.So, high concentration at the free surface inhibits the movement of the surface.This means stability of the free surface.The inverse behavior was obtained for a high concentration at the core region.Similar trend was obtained experimentally by [32].
Finally, our calculation of the free surface function was compared with the previous results of [3] in the absence of the porous medium and the electric field.Also, the present results were compared with the experimental results of [21,31,32].There was compatibility between the results.

Figure 1 :
Figure 1: Physical model and flow cylindrical geometry.
where the notation ‖ ‖ denotes the subtract at the fluid and surrounding vacuum.In other words ‖‖ =  1 −  2 .The subscript , in the electric potential   and the dielectric constant   , takes the values  = 1, 2 at the fluid and the surrounding vacuum, respectively.The continuity of the tangential electric filed components at the interface requires  × ‖‖ = 0 and in the linear form (where the linear terms only are retained and the nonlinear terms are omitted) becomes                      = 0, at  = ℎ.