Abstract

A novel procedure is utilized to investigate the surface waves between two finite conducting fluids streaming through a porous medium in the presence of a horizontal electric field. Normal mode analysis is applied to study two- and three-dimension disturbances cases. The quadratic dispersion equation of complex coefficients representing the system is derived and discussed. It is noted that based on appropriate data selections, the stability criteria do not depend on the medium permeability. It is found that electrical conductivities, viscosities, medium porosity, and surface tension enhance the stability of the system while the dimension and the fluid velocities decrease the stability of the system. Finally, the fluid depths have a dual role (stabilizing as well as destabilizing effects) on the system.

1. Introduction

The Kelvin–Helmholtz instability (KHI) appears if two superposed fluids streaming with a relative horizontal velocity. This instability is very important in describing more kinds of astrophysical and space phenomena involving sheared stream. As a consequence of its connection to geophysical, astrophysical, and laboratory areas [1], it has been examined by numerous contributors. The linear KHI is treated in Chandrasekhar’s monogragh [2]. The impact of streaming is destabilized in the linear access in the KHI problem.

Electrohydrodynamics (EHD) can be treated as the branch of fluid mechanics interested in electrical force effects or as the part of EHD that is concerned with the influence of fluid motion on electric fields. Consequently, it is concerned with (electric fields polarized or free charges) interactions in fluids. The fluids may be very good insulators, slight conductors, or even high conductors. A survey on EHD subject with specific reference to several of the advancements in this area is provided by Melcher [3].

In the surface wave phenomena through dielectric fluids when both surface tension and gravity are essential, the capillary waves and the gravity couple to create one wave distinguished by the two influences. By contrast, the effect of an external electric field during conducting fluid flows has many several characteristics and various valuable engineering and physical applications. When surface charges share in the motion of fluids, the arising waves are mentioned as EHD surface waves. The effect of an electric field on the linear stability of a sharp interface dividing two nonconducting dielectric fluids of unlimited range has been investigated by Melcher and Schwarz [4]. A review of the role of interfacial stresses on EHD has been presented by Melcher and Taylor [5] and later the same characteristics have been reviewed by many authors (see Baygents and Baldessare [6] and Rudraiah et al. [7]).

In the linear EHD stability theory, it is well known that the horizontal electric field has a stabilizing effect while the vertical one always has a destabilizing influence [812]. EHD stability problems of the interface between conducting fluids in planar and cylindrical geometries have been as of late examined by many authors [1318].

Flows through a porous medium have been a topic of high interest for the latest many decades. This interest was driven by many engineering and geophysical applications in various majors [19]. Many recent works on this topic are investigated by many authors [2023]. In most earlier investigations on porous media, many researchers considered the treatments based on Darcy’s law and Forschheimer-extended Darcy’s law models. Darcy’s law is an experimental formula relating the bulk viscous resistance, the gravitational force, and the pressure gradient in a porous medium.

The KHI for flow in porous media has received great attention in the scientific literature. Sharma and Spanos [24] investigated the instability of the plane interface separating two uniform superposed fluids motion through porous media. El-Sayed [25] analyzed the instability of two superposed viscous streaming fluids through porous media in the presence of a horizontal electric field. For recent surveys regarding the developments of linear and nonlinear electromagnetic flows in porous media, see refs. [2636].

The goal of the present paper is to investigate the linear stability of two electrically conducting fluids of finite thickness moving through a porous medium under a uniform horizontal electric field. In addition, the uniform flow is taken through three dimensions to be more general and then using the normal modes analysis to obtain the solutions of system parameters. This problem, as far as we know, has not been studied yet. The stability results are presented in figures. All these flows are new and are presented for the first time in the literature.

2. Problem Description

We consider the three-dimensional finite amplitude surface waves propagating at the interface that represent the balance situation. The perturbed interface separates two conducting incompressible fluids that have a uniform thickness. Fluid (1) occupies the lower region while fluid (2) occupies the upper region . The two fluids are streaming with uniform velocities and along the -axis through a porous medium. The system is influenced by a constant horizontal electric field in the -direction. The gravitational acceleration directed normal to the interface. The surface tension forces are taken into account between the two fluids. We shall denote , , , , , , and , to the fluids densities, viscosities coefficients, dielectric constants, electrical conductivities, electric field elements, permeability, and porosity of the medium. A sketch of the physical problem is given in Figure 1. The disturbed interface may be written as

The unit outer perpendicular to the interface is

The equation of motion for an incompressible viscous fluid in a porous medium is

The continuity equation iswhere are the pressure and the velocity, respectively. Considering irrotational fluid flows, then the velocity potentials is such thatwhere is the unit vector in the -direction. The potentials satisfy the Laplace’s equations:

It is supposed that, in electrodynamics, the quasi-static approximation is valid. Hence, the electrical equations are

The charge conservation equation iswhere is the free current density (stationary current), is the electrical conductivity, and is the free charge density. Therefore, the electric field can be represented as the gradient of electrostatic potentials :

From equations (8) and (10), the scalar electric potentials to satisfy the Laplace’s equations are as follows:

3. The Boundary Conditions

(1)The vanishing of the normal fluid velocities on both the lower and upper boundaries is(2)The vanishing of the normal components of the electric fields also on these boundaries is(3)The kinematic boundary condition at the interface is given by(4)The tangential electric field component is supposed continuous at the interface :where describes the jump of the quantity over the interface.(5)When uniform conductivity of fluids is considered, the problem became more attractive yet difficult. Surface charges play an essential role in such cases; the conduction process in the interface region is a significant factor in many electrical systems. The continuity of stationary current normal to the interface must lead to charge accumulation on the interface [17]:(6)The interfacial normal stress component must be continuous at the interface :

4. Linear Theory and Characteristic Equation

To study the stability of the problem, first we shall analyze the perturbation quantities into normal modes, and we suppose that all the perturbed quantities take the formwhere is any function of , is the carrier wave phase, , and represent the wave number elements in and directions, respectively, while , , , and denote the angular frequency, the complex amplitude of the surface wave, the imaginary unit, and the complex conjugate of the previous terms, respectively.

Here, using equation (19) together with equations (6), (7), (11), and (12) with the suitable boundary conditions, the potential function solutions may be written as

Now, , and to obtain a nontrivial starting solution, the frequency and the wave number should satisfy the next characteristic equation:

Equation (21) is simplified in the dispersion relationwhere

We have the following special cases:(1)In two semiinfinite nonconducting fluids, equation (22) reduces to the corresponding dispersion relation obtained first by El-Sayed [25](2)In two-dimensional nonporous case, equation (22) reduces to the same equation derived by Melcher [8](3)In the absence of the electric field and the porosity, equation (22) reduces to the corresponding relation presented by Chandrasekhar [2]

5. Stability Analysis and Discussion

Applying here the Routh–Hurwitz stability criterion [37] to equation (22), the system will be linearly stable according to

Since and are positive parameters, the first condition in (24) is trivially satisfied. The second inequality is satisfied ifwhere

From the previous stability criteria, it is obvious that the permeability of the meduim has no effect on the linear stability of the system. The horizontal electric field has a dual role according to the sign of the term in contrast with the pure dielectric fluids case, and this was well studied by Moatimid [38].The surface tension has stabilizing effects, while the uniform streaming velocity is strictly destabilizing. Next, we shall make a numerical discussion to clarify the stability of the problem by graphing the transition curves represented by equation (26) in plane when . These transition curves that represent the marginal stability state separate the upper stable S regions from the lower U unstable ones. The calculations are made in that case, for a system having  gr/,  gr/,  farad/cm,  farad/cm,  mho/cm,  mho/cm,  cm,  cm, /sec, /sec,  cm/sec,  cm/sec,  cm/,  dyn/cm, and  sec/cm. In Figures 28, the first dimension expresses disturbance wave number , the second expresses the critical electric field , and the third expresses one of the previous parameters.

Figure 2 shows the variation of the electric field versus for different electrical conductivity values in the two-dimension and three-dimension cases. It was clear that, for , the stable region increases by increasing , while it decreases by increasing the dimension for small wave number range after that the stability does not depend on since the curves coincide. Also, for any value and small values, the system is always unstable, while for high , the electric field has a stabilizing effect since stable region increases by increasing the electric field and wave number value . Then, both the electrical conductivites and the electric field have stabilizing effects, while the dimension has a destabilizing effect only for small wave number value .

Also the effects of the porosity of the medium (see Figure 3) as well as the fluid viscosities and (see Figure 4) on the stability are found to be similar to the effect of the electrical conductivites given in Figure 2. Thus, each of the porosity of the medium and the fluid viscosities and has stabilizing effects on the system.

In Figure 5, we draw versus for different surface tension values in two- and three-dimension disturbances cases. It was clear that, when and for small values of , the stability increases by increasing while it decreases by increasing the dimension for small wave numbers just like the previous figures. In contrast, for high values of , t is seen that the stability increases by increasing and the dimension for wave numbers range . We conclude that, the dimension plays a dual role (stablizing and destablizing) when there are large variation in surface tension force between the fluids.

Figure 6 shows the variation of against for different fluid velocities values and in two- and three-dimension cases. It is shown that, when , the instability increases by increasing the fluid velocities and , and for any fixed values of and , the instability also increases by increasing for small wave number values after which has no effect on the stability. Therefore, both the fluid velocities and and the dimension have destabilizing effects, while for the high wave number value , the dimension has no effect on the stability.

Finally, Figures 7 and 8 show the effects of the fluid depths and accompanied with the dimension on the stability, stabilizing the system for small values of and (see Figure 7) while destabilizing it for large values of and (see Figure 8). Thus, the fluid depths and the dimension have a dual role in the system.

6. Conclusions

In this paper, new results of bounded flows through porous medium have been presented and analyzed for weakly electrically conducting fluids. Exact analytical solutions have been given for all parameters.

Using the normal mode analysis, we obtain quadratic dispersion relation of complex coefficients characterizing the behaviour of the disturbed system. Based on appropriate data selections, we conclude the following:(1)The permeability of the medium has no effect on the linear stability of the system(2)The horizontal electric field , porosity of the medium , surface tension , fluid viscosities and , and electrical conductivites have stabilizing effects or enhance the stability(3)The fluid velocities and (including the case of absence of fluid velocities and or Rayleigh–Taylor instability) have destabilizing effects or tend to reduce the stability(4)Both the fluid depth and and the dimension have a dual role in the system, i.e., stabilizing the system for small values of and , while destabilizing it for large values of and (including the case of two semiinfinite fluids).

According to the importance of nonlinear effects in EHD phenomena, the describing equations of EHD flow are nonlinear. Therefore, we will discuss the nonlinear stability characteristics for this problem in a subsequent article.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This project was supported by the Deanship of Scientific Research at Prince Sattam Bin Abdulaziz University under the research project 2017/01/7102.