Dynamical System of Three Magnetic Layers in the Presence of Porous Media

This paper concerns the linear stability of three viscous fluid layers in porous media. The system is composed of a middle fluid embedded between two semi-infinite fluids, in which the effect of the normal magnetic field is to introduce. The principle aim of this work is to investigate the influence of fluid viscosity and the porosity effect on the growth rate in the presence of normal magnetic field. The parameters governing the layers flow system, the magnetic properties and porosity effects strongly influence the wave forms and their amplitudes and hence the stability of the fluid. The stability criteria are discussed theoretically and numerically and stability diagrams are obtained, where regions of stability and instability are identified. It is found that the stabilizing role for the magnetic field is retarded when the flow is in porous media. Moreover, the increase in the values of permeability parameters plays a dual role, in stability behavior. It has been found that the phenomenon of the dual (to be either stabilizing or destabilizing) role is found for increasing the permeability parameter. It is established that both the viscosity coefficient and the magnetic permeability damps the growth rate, introducing stabilizing influence. The role of the magnetic field and Reynolds number is to increase the amplitude of the disturbance leading to the destabilization state of the flow system, promote the oscillatory behavior. Influence of the various parameters of the problem on the interface stability is thoroughly discussed.


Introduction
There has been a great deal of interest in magnetic fluids of the stability of hydrodynamic stability.A magnetic fluid, also known in the literature as a ferromagnetic fluid or simply a ferrofluid.In many previous researches has shown the importance of studying the hydrodynamic stability problems, for example capillary-gravity waves of permanent form at the interface between two unbounded magnetic fluids in porous media are investigated in paper [1].In which the system is influenced by the horizontal direction of the magnetic field to the separation face of two semi-infinite homogeneous and incompressible fluids.It is found that the stabilizing role for the magnetic field is retarded when the flow is in porous media.In paper [2] the influence of viscosity on the stability of the plane interface separating two incompressible superposed fluids of uniform densities, when the whole system is immersed in a uniform horizontal magnetic field, has been studied.The authors in the research [3] have carried out the instability of viscous potential flow in a horizontal rectangular channel.The analysis leads to an explicit dispersion relation in which the effects of surface tension and viscosity on the normal stress are not neglected but the effect of shear stresses is.The unsteady electrohydrodynamic stability has been investigated by Elhefnawy [4], where the stability analysis has been made of a basic flow of streaming fluids in the presence of an oblique periodic electric field.A series of studies for hydrodynamics stability have also been initiated by Drazin and Reid [5] and Joseph [6].
The authors in paper [7] have investigated the stability properties of time periodically streaming superposed magnetic fluids through porous media under the influence of an oblique alternating magnetic field.Also, Zakaria et al. [8] have analyzed the effect of an externally applied electric field on the stability of a thin fluid film over an inclined porous plane, using linear and non-linear stability analysis in the long wave limit.Wray et al. [9] have investigate the evolution and stability of a wetting viscous fluid layer flowing down the surface of a cylinder, and surrounded by a conductive gas.Liu et al. [10] have studied the instability properties of two-dimensional non-Newtonian liquid sheets moving in an inviscid gaseous environment.
Based on a modified Darcy's law for a viscoelastic fluid, Sirwah [11] has discussed the linear stability analysis of the electrified surface separating two coaxial Oldroyd-B fluid layers confined between two impermeable rigid cylinders in the presence of both interfacial insoluble surfactant and surface charge through porous media.Also, Tan and Masuoka [12] have extended Stokes?first problem to that for an Oldroyd-B fluid in a porous half space, where an exact solution was obtained by using Fourier sine transform.Zakaria [13] has investigated the time evolution of superposed layers of fluid flowing down inside an inclined permeable channel.Using the Kármán-Pohlhausen approximation, the problem is reduced to the study of the evolution equation for the liquid?liquid interface of the liquids film derived through a long wave approximation.Khan et al. [14] have demonstrated the analytical solutions for the magnetohydrodynamic flow of an Oldroyd-B fluid through a porous medium.They obtained the expressions for the velocity field and the tangential stress by means of the Fourier sine transform.Kumar and Singh [15] have investigated the stability of a plane interface separating two viscoelastic (Rivlin-Ericksen) superposed fluids in the presence of suspended particles.
In the present work we wish to consider an investigation of stability for flow in a porous medium under the effect of a magnetic field that is normal to the interface between the fluids.The considered system is composed of a viscous fluid layer of finite thickness embedded between two semi-infinite fluids.We have concentrated in this work to investigate the mechanisms of stability of three porous layers of fluids in the presence of normal magnetic field.This paper is organized as follows: This section has presented the motivation for the investigation in addition to relevant background information.In next section, we will give a formulation of the problem statement, including the basic equations of the fluid mechanics and Maxwell's equations governing the motion of our model.In the third section and its subsections, are concerned with the derivation of the characteristic equation and numerical estimations for stability configuration.The salient results of our analysis are discussed and some important conclusions are drawn in final section of this paper.

Formulation of the Problem
The system under consideration is composed of an infinite horizontal viscous fluid sheet of vertical height 2L confined between two semi-infinite superposed incompressible viscous fluids.The fluids are considered to be influenced by the gravity force in the negative y-direction.The x-axis of the co-ordinate system is parallel to the direction of the fluid sheet flow, and the y-axis is normal to the fluid sheet with its origin located at the middle plane of the fluid sheet as shown in

Scaling Variables for the Evolution Equations
First, we will use the dimensionless variables to provide improved insight into the physics and in order to understand hydrodynamic stability better.So we define the corresponding dimensionless variables using the half thickness of the middle fluid sheet L as a length scale.Thus the stream velocity and the time are made dimensionless using Lg and L g , while the applied magnetic field and the magnetic potential are made dimen- sionless by Lg ρ µ and ( ) ( ) L Lg ρ µ , respectively.In addition the viscosity ( ) , where j T is the surface tension coefficient.The dynamics of the problem are described by the simultaneous solution of three field equations: Maxwell's equations, Navier-Stokes equation, and the continuity equation.Assuming a quiescent initial state, therefore the base state velocity in the fluid layers is zero in which the flow is steady and fully developed.Fluid flow through a porous medium is often given by the phenomenological Darcy's equation.Thus, the equations governing twodimensional motion of a viscous incompressible fluid through porous medium are [3] and the equation of continuity will be where is the horizontal gradient operator.
In writing Maxwell's equations for the problem, we supposed that the electro-quasi-static approximation is valid for the problem, and hence the magnetic field equations read refers to the magnetic field intensity vector and ( ) r µ is the magnetic permeability.The construc- tion of a potential function ( ) r χ , can be representable as the gradient of the scalar potential such that and thus we have the Laplace equation in the form where 1 e and 2 e are unit vectors in x − and y − directions.

Boundary Conditions
Solution of the equations of motion cited before is accomplished by utilizing the convenient boundary conditions.The flow field solutions of the above governing equations have to satisfy the kinematic and dynamic boundary conditions at the two interfaces, which can be taken as 1 y ≈ ± .The normal component of the velocity vector in each of the phases of the system is continuous at dividing surface [16] [17]: where , e e is the outward normal unit vector to the interfaces.The condition that the inter- faces are moving with the fluids lead to In addition the jump in the shearing stresses is zero across the interfaces, this gives where, ( ) is the velocity vector due to disturbances and the notation ⋅ denotes the jump of a quan- tity across the interfaces.
Furthermore, the Maxwell's conditions on the magnetic field where no free surface charges are present on the interfaces.The continuity of the normal component of the magnetic displacement at the interfaces reads: ˆ0 at 1 and 0 at 1 y y y y y y The tangential component of the magnetic field is zero across the interfaces, this requires that from this equation, we have where, we use the zero order from continuity of the normal component of magnetic field to express both H .The completion of the mathematical description of the problem requires an additional interfacial condition determine the shape of the interface between the fluids, which is the dynamical equilibrium boundary condition in which the surface traction suffers a discontinuity due to the surface tension: These boundary conditions represented here are prescribed at the interface j y ξ = .As the interface is de- formed, all variables are slightly perturbed from their equilibrium values.

Linear Perturbation and Solutions
The analysis of linear theory, as presented in Chandrasekhar book [18], depends on neglecting the nonlinear terms from equation of motion as well as from the boundary conditions.The solution of the above system of governing equations and boundary conditions can be facilitated by defining a stream function, ψ of the time and space coordinates, which automatically satisfies Equation ( 2), where To solve the equations for the fluid phases under consideration, the two-dimensional finite disturbances are introduced into the equation of motion and continuity equation as well as the boundary conditions.As a customary in hydrodynamic stability analysis [18], where all quantities have exponential time and a spatial dependence and in view of a standard Fourier decomposition, we may assume that the solutions are in the form ( ) , e ikx t j j x t c c where ˆl ξ is the initial amplitude of the disturbance, which is taken to be much smaller than the half-thickness L of the middle sheet, k is the wave number of the disturbance, which is assumed to be real and positive, ω is a complex frequency ( , where r ω represents the rate of growth of the disturbance, i ω is 2π times the disturbance frequency), the symbol i denotes 1 − , the imaginary number and c c ⋅ represents the complex conjugate of the preceding terms.
Eliminating the pressure term from Equations ( 1) and (2) and using ( 15) and ( 17), we obtain the following equation ˆ1 Using the normal mode solution (17) we can obtain the pressure from Equation (1): Substituting, the solution of the analytical solution of Equation ( 18) into Equation ( 15) we get ( ) ( ) Also, the solution of the magnetic potential, in view of Equation ( 5) may be taken the form

Derivation of the Characteristic Equation
In this section, we will derive the dispersion relation controlling the stability behavior of the system.When the obtained solutions of the stream function, magnetic potential and surface tension are inserted into Equations ( 6)-( 12), we have a linear homogeneous system of algebraic equations of the fourteen unknown coefficients pj A , lj B , ˆl ξ , ( ) 1, 2,3, 4 p = .These homogeneous system of equations can be expressed in matrix form as = MZ 0 (22) where 0 is a null vector, Z is a vector of unknown coefficients defined as ( ) where the superscript T indicates the matrix transpose.A non-trivial solutions of the unknown coefficients pr A , lr B , ˆj ξ , exists if and only if the determinant of the 14 × 14 matrix M must be equal to zero, which yields a dis- persion relation between the wave number k and the perturbation frequency ω for specified values of other parameters, given by which represents the linear dispersion equation for surface waves propagating through a viscous layer embedded between two other fluids with the influence of constant horizontal magnetic field.This dispersion relation controls the stability in the present problem.That is, each negative of the real part of ω corresponds to a stable mode of the interfacial disturbance.On the other hand, if the real part of ω is positive, the disturbance will grow in time and the flow becomes unstable.It is clear that the eigenvalue relation ( 24) is somewhat more general and quite complex, since r  involves square roots and so one can obtain other characteristic relation as limiting cases.For an inviscid fluid we get the characteristic equation as special case from Equation (24) when ( )  ˆ0 r η = .Thus by collecting the real and the imaginary terms in power order of ω with the help of symbolic computation software Mathematica, Equation (24) can be transformed into a polynomial algebraic equation of fourth order in the frequency ω.Zakaria et al.
[7] obtained a similar equation in their study of temporal stability of an inviscid fluids in porous media.Also, in the special case when the effect of the magnetic forces is absent and for the fluids flow through no-porous media, we get and in this case the dispersion relation ( 28) is reduced to a non-polynomial algebraic equation for the frequency ω which coincides with that obtained by Kwak and Pozrikidis [19].Another case is the limiting case of one interface between two continuum layers (non-porous medium), in which highly viscous fluids are considered.Thus we obtain a polynomial equation of fifth order in ω, which is obtained before by Kumar and Singh [15] and Sunil et al. [20].
In the following, numerical applications are carried out to demonstrate the effects of various physical parameters on the stability criteria of the system.In the present work, we will numerically solve the implicit dispersion relation by means of the Chebyshev spectral tau method [21].

Numerical Estimation for Stability Configuration
In this section, the goal is to determine the numerical assess for the stability pictures for surface waves propagating through porous media.In order to present this examination, Equation ( 24) is used to control the stability behavior, which requires specification of the parameters: the magnetic field, the magnetic permeability, the porosity effect, the density, the viscosity.In the calculations given below all the parameters are sought in the dimensionless form as defined above.The stability of fluid sheets corresponds to negative values of the disturbance growth rate (i.e.0 r ω < ), and the disturbance growth rates of different fluids can be gained through solving the above corresponding dispersion relation numerically.
To show the effect of changes of the magnetic permeability ratio ( ) The graph displayed in this plane is evaluated for a system having the parameters given in the caption of Figure 2, while the magnetic permeability ( ) 3 μ has some variation for the sake of comparison.In the graph 2 the values 0.5, 1.5 and 1.5 are selected for ( ) correspond to the continuous, dashed and dotted curves respectively.Before, we discuss the stability of this graph we firstly define the critical wave number (also called the cutoff wave number) as given in [10] the value of the wave number at the point where the growth rate curve crosses the wave number axis in the plots of wave growth rate versus wave number.In other words the critical wave number is the value of the wave number, which separates the stable motions from the unstable ones and conversely, and can be obtained from the corresponding dispersion relations by setting 0 r ω = .It is apparent from the inspection of Figure 2, under the influence of the magnetic permeability ( ) 3 μ , the growth rates with different magnetic permeability ratio keep almost identical for the wave numbers less than 0.4, but increase  correspondingly at higher values of the wave number, further the plane ( ) The first is 0 0.11 k < < , which represent a stabilizing effect for increasing the parameter ( ) 3 μ .The second region lie in the range 1.5 < < 0.11 k , since in this range, we notice that, when the magnetic permeability is increased, both the growth rates and the cutoff wave numbers of fluid sheets decrease.A general conclusion of the graph 2 reveals that, the stabilizing influence is found for increasing magnetic permeability ratio.
The examination of the influence of the magnetic field H on the stability picture is illustrated in Figure 3, for the same system considered in Figure 2. The results for calculations are displayed in the plane ( ) r k ω − .Since the stability arises according to the negative sign of the real part of the complex frequency ω.Thus when the wave number is over the cutoff wave number, the fluid sheet is unstable.In this figure the solid curve is plotted at the value 3 H = .The inspection of Figure 3 indicates that as the magnetic field is increased the growth rates increased, on other meaning the stable regions under the curves are decreased.Therefore, it is concluded that the magnetic field effects has destabilizing influence in the fluid sheets.
The examination of change of the lower to the middle fluid viscosity ratio ( ) 1 η in the stability criteria is illu- strated in Figure 4, where the values 0.2, 2.5 and 3 are choosing for the quantity ( ) 1 η .It is obvious from this graph, for every value of the ( ) 1 η , the corresponding curve crosses the wave number axis at three points (the cutoff wave numbers) and formed areas of stability and instability regions.Inspection of Figure 4 revels that the increasing of the viscosity ratio ( ) 1 η tends to a reduction in the width of the unstable regions, whereas the sta- bility region extended under the influence of the increasing ( ) η has a stabilizing influence on the stability of the movement of the waves.Ozen et al. [22] have been ob- tained a similar conclusion in their studies of electrohydrodynamic linear stability of two immiscible fluids in channel flow.In Figure 5, in which the real part of the frequency ω is plotted against the wave number, the Reynolds number e R has the values 0.1, 0.3 and 0.7 to show its effect on the stability picture.Having checked the stability diagrams of this figure, it is discovered that the increasing of the Reynolds number leads to a contraction in the stability areas under the wave number axis, and consequently the Reynolds number has a destabilizing role on the stability behavior.Similar results were reported by Liu et al. [10] in their studies of the instability of two-dimensional non-Newtonian liquid sheets., leads to an contraction in the width of the instability regions (the regions under the curves and above the wave number axis correspond to the positive sing of the disturbance growth rate).On the other hand, through the interval 0.7 2 e R < < the growth rates of insta- bilities with different ( )

Q
are increased.In general view of the graph 6, it is noticed that there are two roles of the variation of the the porous parameter ( ) 2 Q , the first one is a stabilizing when the Reynolds number e R less than the value 0.7, and the other role is a destabilizing when e R lies between the values 0.7 and 2. Hence the phenomenon of the dual (irregular) role is found for increasing the permeability parameter ( ) 2 Q .The influence of magnetic field

H
on the flow behavior in terms of streamlines field is discussed through the parts of Figure 7 (a curve formed by the velocity vectors of each fluid particle at a certain time is called a streamline, in which the tangent at each point of this curve indicates the direction of fluid at that point).The streamlines in the physical domain are thus mapped into horizontal grid lines in the computational plane, thus resulting in a rectangular computational region.The streamlines show to be very effective tools to visualize a qualitative impression of the flow behavior during the motion.In the parts (a-c) of Figure 7, the streamlines picture is achieved by fixing the value of all the physical parameters as given in Figure 2

H
are increased to 1.2 and 1.8 respectively.
A conclusion that may be made from the comparison among the parts (a-c) of Figure 7 is that the magnetic field leads to crowd in the concentration of the streamlines in the movement of the fluids.In other words, in the light of stability configuration, we notice that corresponding the parts (a-c) of Figure 7 there are three different values of the disturbance growth rate ( ) r ω , which are 0.049, 0.66 and 0.502.Since the stability of fluid sheets arises to negative values of growth rate, thus it can be observed that the streamlines contours represent an unstable system.Hence the magnetic field has destabilizing influence, which coincide with the result given in Figure 3.

Conclusions
This work is concerned with the influence of the normal magnetic field on the gravitational stability of a viscous fluid sheet of finite thickness.The sheet is embedded between two semi-infinite fluids layers moving in porous media, under the influence of magnetic field.The solutions of the linearized equations of motion under the boundary conditions lead to an implicit dispersion relation between the growth rate and wave number.The parameters governing the layers flow system, the magnetic properties and porosity effects strongly influence the wave forms and their amplitudes and hence the stability of the fluid.The stability criteria have been performed theoretically and numerically in which the physical parameters are put in the dimensionless form.Some stability diagrams have been plotted and discussed, in which the influence of the various parameters of the problem on the interface stability is thoroughly analyzed.It has been found that the phenomenon of the dual (to be either stabilizing or destabilizing) role is found for increasing the permeability parameter.It is established that both the viscosity coefficient and the magnetic permeability damps the growth rate, introducing stabilizing influence, where a part of its kinetic energy may be absorbed.However, it is expected to be a more careful search would clarify that the motion of the interfacial waves will be more stable with the increase of the values of the viscosity as well as the magnetic permeability.In addition an increase of the lower to the middle fluid viscosity ratio decrease both the growth rate and the stability range of fluid sheet, which give a stabilizing influence on the stability behavior of the waves.This result confirmed the fact that when the lower fluid is more viscous than the upper, thus the system is stable.The role of the magnetic field and Reynolds number is to increase the amplitude of the disturbance leading to the destabilization state of the flow system, promote the oscillatory behavior.

Figure 1 .
The superscripts (1),(2) and(3) refer to the lower fluid, middle sheet and upper fluid respectively, where ( ) r ρ , ( ) 1, 2, 3 r = is the fluid densities, ( ) r µ is the magnetic per- meability and ( ) r η represents the viscosities of the fluids.The fluids are subjected to external magnetic field

Figure 1 .
Figure 1.Sketch for the system under consideration.The y-axis is taken vertically upwards, the x-axis is taken horizontally at the half of the middle sheet.

2 Lg ρ , the stream function 3 L
g and the permeability of the porous medium 2 L Q .Furthermore, in the equations of motion, we use the symbols: the magnetic permeability ratio( ) number of the middle layer and( ) the calculations are displayed in Figure 2, partitions (a), (b) and (c) indicate the plane ( )

2 H
= corresponds to the dashed line, while the dotted curve represents the value

1 η.
It is clear that from Figure 4 the viscosity ratio ( ) 1

Figure 3 .
Figure 3. Effects of the magnetic field ( ) 1 0 H in the plane ( ) r k ω − , on the wave growth rate with

Figure 5 .
Figure 5. Represents the stability diagrams in the plane ( ) r k ω −

Figure 6 2 Q 2 Q
Figure 6 exhibits the effects of the the permeability parameter ( ) 2 Q on the stability behavior of the fluid layers.In this graph the solid, dashed and dotted curves represent the values 0.2, 0.4 and 0.9 in the plane ( ) r e R ω −

H
has three value for comparison.Snapshots of instantaneous streamlines of the stream function, are shown in Figure 7(a) at of this graph reveals that the flow consists of cells (contours) consisting of clockwise (positive values of streamlines) and anti clockwise (negative values of streamlines) circulations.In parts (b) and (c) of this graph, the values of ( ) 1 0