A New Approach for the Stability Analysis in Hydromagnetic Couette Flow ()
1. 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 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
, Hartmann parameter
, wave number
and permeability parameter
on the flow. The conclusion is presented in the last section.
2. Modified Orr-Sommerfeld Equation
We considere a Couette viscous incompressilbe, electrically conducting fluid flow between two porous parallel plates of infinite lengh, distant
apart in the presence of uniform transverse constant magnetic field
applied parallel to
axis which is normal to the planes of the plates. We considered the simple case where,
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
, at the lower plate and a same small constant suction
, at the upper plate. The upper plate is allowed to move with non-zero uniform velocity
in flow direction and the lower plate is kept at rest. We choose the origine on the plane
such as
and
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,
, both the fluid and plates are assumed to be at rest. When
, 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:
(1)
(2)
(3)
(4)
(5)
(6)
(7)
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
(8)
and
,
,
,
,
,
,
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
denotes the time. The physical model of the problem is illustrated in Figure 1, where
is the velocity vector in the
directions respectively.
(9)
(10)
(11)
where
is a constant. We assumed that no applied polarization voltage exists (i.e.,
). Then Equation (8) and Equation (11) give
(12)
and Equation (7) yields
(13)
We introduce the following non-dimensional quantities
,
,
,
,
,
,
,
,
(hydrodynamic Reynolds number),
(injection/suction Reynolds number),
(Hartmann parameter),
(permeability parameter).
So, Equation (1) and Equation (2) become
(14)
(15)
(16)
(17)
For the stability analysis, the flow is decomposed into the mean flow and the disturbance according to
(18)
(19)
We take the dimensional basic flow for small suction and injection [2] [4] [5] ,
(20)
(21)
(22)
By scaling these velocities as above, we obtain with
(
) the no-dimensional base flow
(23)
(24)
(25)
To obtain the stability equations for the spatial evolution of three-dimensional, we take the dependent on time disturbances
(26)
which are scaled in the same way as above.
Inserting Equations (18) (26) into Equations (15)-(17), we get
(27)
(28)
(29)
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
(30)
Taking the laplacian of Equation (28), we get after linearization with Equation (30)
(31)
The disturbances are taken to be periodic in the streamwise, spanwise directions and time, which allow us to assume solutions of the form
(32)
where
represents either one of the disturbances
,
,
or
and
the amplitude function;
and
are the wave numbers respectively on
and
axis directions;
is the frequency of the wave;
,
,
is the wave velocity which is taken to be complex,
and
are real because of temporal stability analysis consideration. Then Equation (31) becomes
(33)
where
; with boundary conditions for all
(34)
Taking
(35)
Equation (33) and the boundary conditions Equation (34), take the forms
(36)
The first equation of system Equation (36) is a flow equation modified by the small injection/suction Reynolds number
, the Hartmann parameter
#Math_100#, and permeability parameter
whichis the so-called
modified Orr-Sommerfeld equation, rewritten as an eigenvalue problem, where
is the eigenvalue and
the eigenfunction.
3. Linear Stability Analysis
We consider a three-dimensional disturbances. We use a temporal stability analysis as mentioned above. With
complex as defined above, when
, a stability mode takes place,
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
(37)
for
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
with
. 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
. Frame
correspond respectively to
,
and
.
Figure 2 presents the effect of injection/suction Reynolds number
on linear temporal stability of viscous incompressible electrically non conducting fluid (
) flow for different values of wave number. It is observed that for
and
(frames a and b), the stability is not affectd by
and the flow is unstable but for
and
(frames c and d),
affects the stability, and the flow stays stable, but increasing of
doesn’t contribute to the stability.
Figure 3 exhibits the effect of permeabilty parameter
on linear temporal stability of viscous incompressible electrically non conducting fluid flow for different values of wave number. It is observed that
affects the stability. For
and
, the frames
and
show that for
the flow is unstable (curves I) and stable for
(curves IV) but for
and
(curves II, III), we have the transition of the flow (See Table 1 for the criticals hydrodynamic Reynolds number values). For
and
(frames
and
) the flow is completely stable. On careful observation, we remark that for
, increasing of
contributes to the stability in frame
case, and the opposite is noticed in the frame
case, but when
, increasing of
contributes in the both cases. Thus, it
Figure 2. Ci vs. Re for
and
variable.
Figure 3. Ci vs. Re for
and
variable.
Figure 4. Ci vs. Re for
and
variable.
may be concluded that except the frame
case where, the increasing of
doesn’t contribute to stability for
, the increasing of permeabilty parameter contributes to the flow stability.
Figure 4 shows the effect of Hartmann parameter
on linear temporal stability of viscous incompressible electrically conducting fluid flow for different values of wave number. It is observed that
affects the stability. For
and
, frames
and
show that for
the flow is unstable (see curves I) and stable for
(see curves IV) but for
and
(curves II, III), we have the transition of the flow (See Table 1 for the critical Reynolds number values). For
and
(frames
and
) the flow is completely stable. Thus, we may concluded that the Hartmann parameter increasing contributes more to the flow stability.
Figure 5 depicts the effect of phase angle
for different values of the wave number on the flow stability. For
and
, frames
and
show that the flow is unstable and the instability increases when
increases. But, for
and
the flow is completely stable except the curve IV frame
, which presents a transition initialy and stays stable after.
Finally, Figures 6-8 (
, electrically conducting fluid) show that for
and
, the small injection/suction has no effect on the linear temporal stability of the flow. But for
and
, we remark a little influence of the
Figure 5. Ci vs. Re for
and
variable.
Figure 6. Ci vs. Re for
and
variable.
Figure 7. Ci vs. Re for
fixed and
variable.
Figure 8. Ci vs. Re for
and
variable.
small injection/suction on the stability only in a small range of
.
The critical Reynolds numbers
for which transition occurs are presented in Table 1.
4. 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
and
, the stabilizing effect appears.
Acknowledgements
The authors would like to thank very much the anonymous referees whose useful criticisms, comments and suggestions have helped strengthen the content and the quality of the paper.