A New Approach for the Stability Analysis in Hydromagnetic Couette Flow

This paper analyses the effects of small injection/suction Reynolds number, Hartmann parameter, permeability parameter and wave number on a viscous incompressible electrically conducting fluid flow in a parallel porous plates forming a channel. The plates of the channel are parallel with the same constant temperature and subjected to a small injection/suction. The upper plate is allowed to move in flow direction and the lower plate is kept at rest. A uniform magnetic field is applied perpendicularly to the plates. The main objective of the paper is to study the effect of the above parameters on temporal linear stability analysis of the flow with a new approach based on modified Orr-Sommerfeld equation. It is obtained that the permeability parameter, the Hartmann parameter and the wave number contribute to the linear temporal stability while the small injection/suction Reynolds number has a negligible effect on the stability.


Introduction
The study of Couette flow in a rectangular channel of an electrically conducting viscous fluid under the action of a transversely applied magnetic field has immediate applications in many devices such as magnetohydrodynamic (MHD) power generators, MHD pumps, accelerators, aerodynamics heating, electrostatic precipitation, polymer technology, petroleum industry, purification of crude oil and fluid droplets spray.Channel flows of a Newtonian fluid with or without heat transfer were studied with or without Hall currents by many authors [1]- [8].
The effects of the injection/suction through the injection/suction parameter number, Hartmann parameter, permeabilty or Darcy parameter on the stability of the fluids flows were also studied with different approach [2] [5] [6] [7] [8].
Indeed, the heat source and the Soret effect on hydromagnetic oscillatory flow through a porous medium bounded by two vertical parallel porous plates have been analyzed by Chand [5], where one plate of the channel is kept stationary and the other is moving with uniform velocity.The plates of the channel are subjected to constant injection and suction velocities respectively.They remarked that the Lorentz force parameter i.e. the Hartmann number contributes to reduce the velocity and the skin friction profiles.The effect of the permeability parameter i.e.Darcy number is just opposite to that of Lorentz force parameter.Nayak et al. [8] studied an oscillatory effect on magnetohydrodynamic flow and heat transfer in rotating horizontal porous channel.They found that magnetic field fixed relative to the moving plate contributes more to the resultant velocity than the magnetic field fixed relative to the fluid in case of all the parameters.Another striking result is that frequency of oscillation has a distinct effect when the magnetic field is fixed relative to the fluid.The effect of all the pertinent parameters on phase angle is just opposite to that of resultant velocity owing to the relative positions of the magnetic fields.In [6], the authors have considered the unsteady hydromagnetic incompressible viscous fluid flow through a porous medium in a horizontal channel under prescribed discharge, under the influence of inclined magnetic field.Das [2] analyzed the effects of constant suction and sinusoidal injection on three dimensional Couette flow of a viscous incompressible electrically conducting fluid through a porous medium between two infinite horizontal parallel porous flat plates in presence of a transverse magnetic field.
The stationary plate and the plate in uniform motion are, respectively, subjected to a transverse sinusoidal injection and uniform suction of the fluid.It is observed that a growing magnetic parameter retards the main velocity and accelerates the cross flow velocity of the flow field and a growing permeability parameter or suction/injection parameter reverses the effect.Both Prandtl number and the suction/injection parameter have retarding effect on the temperature field.Further, a growing suction/injection parameter diminishes both the components of skin friction at the wall while the permeability parameter enhances the x-component and reduces the z-component of the skin friction at the wall.In [7], the authors studied the effects of variable viscosity and heat source on unsteady laminar flow of dusty conducting fluid between parallel porous plates through porous medium with temperature dependent viscosity.It is assumed that the parallel plates are porous and subjected to a uniform suction from above and injection from below.They found that the effect of the suction velocity on both the velocity and temperature of the fluid and particles is more pronounced for higher values of the porosity parameter.
In the present paper, we studied the effects of above parameters on the linear temporal stability of the fluid, in a Couette horizontal porous channel flow with the presence of a uniform transverse magnetic field fixed relative to the fluid.We

Modified Orr-Sommerfeld Equation
We considere a Couette viscous incompressilbe, electrically conducting fluid flow between two porous parallel plates of infinite lengh, distant h apart in the presence of uniform transverse constant magnetic field 0 B applied parallel to * y axis which is normal to the planes of the plates.We considered the simple case where, 0 B is fixed relative to the fluid.We work at constant temperature, the heat transfer aspect of the flow is not studied.We applied a small constant injection V ω , at the lower plate and a same small constant suction V ω , at the upper plate.The upper plate is allowed to move with non-zero uniform velocity 0 U U = in flow direction and the lower plate is kept at rest.We choose the origine on the plane ( ) and * x parallel to the direction of the motion of the upper plate.We assumed the magnetic Reynolds number very small for metallic liquids and neglected the induced magnetic field in comparison with the applied one [3] [7] [8].Initially, * 0 t = , both the fluid and plates are assumed to be at rest.When * 0 t > , the upper plate starts moving with a constant velocity U in coordinate system with the fluid.The equations of continuity, motion for the viscous incompressible electrically conducting fluid in vector form are: ( ) Figure 1.Physical model and coordinate system.
where ( 1)-( 7) are continuity, Newton's second law, Ampere's law, Faraday's law, Maxwell's law and Gauss law equations respectively, with ( ), and ( ) µ , e ε are the velocity, the magnetic field, the electric field, the current density vector, the fluid electrical conductivity, the magnetic permeability and absolute permittivity of the fluid respectively and * t denotes the time.The physical model of the problem is illustrated in Figure 1, where ( ) , , V u v w is the velocity vector in the * * * , , x y z directions respectively.
For the stability analysis, the flow is decomposed into the mean flow and the disturbance according to We take the dimensional basic flow for small suction and injection [2] [4] [5], ( ) By scaling these velocities as above, we obtain with To obtain the stability equations for the spatial evolution of three-dimensional, we take the dependent on time disturbances .
The pressure terms can be eliminated from Navier-Stokes equations.For such a mean profile (base flow), the divergence of Navier-Stokes equations and continuity, give The disturbances are taken to be periodic in the streamwise, spanwise directions and time, which allow us to assume solutions of the form i x z t f x y z t f y where f represents either one of the disturbances u , v , w or p and f the amplitude function; cos  ; with boundary conditions for all ( ) The first equation of system Equation ( 36) is a flow equation modified by the small injection/suction Reynolds number e R ω , the Hartmann parameter , and permeability parameter which is the so-called modified Orr-Sommerfeld equation, rewritten as an eigenvalue problem, where c is the eigenvalue and ˆp v the eigenfunction.

Linear Stability Analysis
We consider a three-dimensional disturbances.We use a temporal stability analysis as mentioned above.With c complex as defined above, when 0 for e R ω small (i.e.small suction) is considered.The eigenvalue problem Equation ( 36) is solved numerically with the suitable boundary conditions.The solutions are found in a layer bounded at 1 y = ± with ( ) ( )          small injection/suction on the stability only in a small range of e R .The critical Reynolds numbers ec R for which transition occurs are presented in Table 1.

Conclusion
In this paper, we have investigated the effects of small injection/suction Reynolds number, Hartmann parameter, permeability parameter and wave number on a viscous incompressilbe electrically conducting fluid flow, in a porous parallel plates forming a channel.We have derived the appropriate equation named modified Orr-Sommerfeld equation in order to make the stability analysis of the flow.Through this approach, we have found that the small injection/suction has a negligible effect on the linear temporal stability of hydromagnetic Couette flow.We noticed that the permeability parameter (Darcy number), the Lorentz force (the Hartmann parameter) and the wave number contribute to the temporal linear stability of hydromagnetic Couette flow.We remarked also that at low wave numbers, the phase angle θ doesn't contribute to the stability of the fluid flow, but for

A
. V. Monwanou et al. used a new approach based on a derived equation named modified Orr-Sommerfeld equation.The corresponding eigenvalue problem is resolved in order to study the linear stability of the flow.The plates of the channel are considered porous and flow within the channel is due to the uniform motion of the upper plate.Such linear temporal stability analysis through the so-called modified Orr-Sommerfeld equation has been made earlier [4] in a Poiseuille flow without injection/suction.The authors showed that the magnetic field has a stabilizing effect on the electrically conducting fluid flows.The paper is organized as follows: Sect. 2 addresses the so-called modified Orr-Sommerfeld equation governing the stability analysis in the hydromagnetic Couette horizontal porous plates flow.Section 3 deals with analysis of the effects of small injection/suction Reynolds number e R ω , Hartmann parameter M , wave number k and permeability parameter p K on the flow.The conclusion is presented in the last section.
Taking the laplacian of Equation (28), we get after linearization with Equation velocity which is taken to be complex, α and β are real because of temporal stability analysis consideration.Then Equation (31) becomes mode takes place, 0 i c = corresponds to neutral stability and elsewhere corresponds to instability.We employ Matlab 7.8.0.(R2009a) version in all our numerical computations to find the eigenvalues.The Couette horizontal porous plates flow with the basic velocity profile 1 ,1, 0 2

.
The results of calculations are presented in the figures below.The black, red, green and blue colors are respectively, curves I, II, III, IV and the yellow color figure corresponds to the neutral mode 0

Figure 2
Figure 2 presents the effect of injection/suction Reynolds number e R ω on linear temporal stability of viscous incompressible electrically non conducting fluid ( 0 M = ) flow for different values of wave number.It is observed that for 1 k = and 1.02 k = (frames a and b), the stability is not affectd by e R ω and the flow is unstable but for 2 k = and 3 k = (frames c and d), e R ω affects the stability, and the flow stays stable, but increasing of e R ω doesn't contribute to the stability.

Figure 3 KFigure 2 .
Figure 3 exhibits the effect of permeabilty parameter

Figure 4
Figure 4 shows the effect of Hartmann parameter M on linear temporal stability of viscous incompressible electrically conducting fluid flow for different values of wave number.It is observed that M affects the stability.For 1 k = and 1.02 k = , frames a and b show that for 10 M = the flow is unstable (see curves I) and stable for 100 M = (see curves IV) but for 50 M = and 80 M = (curves II, III), we have the transition of the flow (See Table 1 for the critical Reynolds number values).For 2 k = and 3 k = (frames c and d )

Figure 5 = and 3 kFigure 5 .
Figure 5 depicts the effect of phase angle θ for different values of the wave number on the flow stability.For 1 k = and 1.02 k = , frames a and b show that the flow is unstable and the instability increases when θ increases.But, for 2 k = and 3 k = the flow is completely stable except the curve IV frame c , which presents a transition initialy and stays stable after.Finally, Figures 6-8(

Table 1 .
Table of critical values.