SlideShare a Scribd company logo
International Journal of Engineering Science Invention
ISSN (Online): 2319 – 6734, ISSN (Print): 2319 – 6726
www.ijesi.org ||Volume 4 Issue 10|| October 2015 || PP.33-42
www.ijesi.org 33 | Page
Analytical solution for Transient MHD flow through a Darcian
porous regime in a Rotating System
1
Hamida Khatun , ¥
Sahin Ahmed
1
Department of Mathematics, South Salmara College, Dhubri-783127, Assam, India,
2
Department of Mathematics, Rajiv Gandhi Central University, Rono Hills, Itanagar, Arunachal
Pradesh-791112, India,
Abstract: A rotating model is developed for a two-dimensional, unsteady, incompressible electrically
conducting, laminar free convection boundary layer flow of heat and mass transport in a saturated porous
medium, bounded by an infinite vertical porous surface in presence of an applied transverse magnetic field. The
porous plane surface and the porous medium are assumed to rotate in a solid body rotation. The vertical
surface is subjected to uniform constant suction perpendicular to it and the temperature at this surface
fluctuates in time about a non-zero constant mean. The basic equations governing the flow are in the form of
partial differential equations and have been reduced to a set of non-linear ordinary differential equations by
applying suitable similarity transformations. The problem is tackled analytically using classical perturbation
technique. Pertinent results with respect to embedded parameters are displayed through graphically and tables
for the velocity, concentration and skin friction profiles and were discussed quantitatively. Applications to the
flows of fluids through porous medium bounded by rotating porous systems find many industrial applications
particularly in the fields of centrifugation, filtration and purification processes.
Keywords: Coriolis and magnetic forces; Mass transport; free convection; porous medium; rotating system.
Mathematics Subject classification: 76(D, W05, 80A32, M45).
NOMENCLATURE
u, v, w velocity components in x, y and z directions respectively [ms1
]
w0 (>0) constant suction velocity of the liquid through the porous plane surface
[ms1
]
Z normal direction of vertical porous plane surface [m]
Z* dimensional normal distance [m]
C* dimensional species concentration [-]
Cp specific heat at constant pressure [J Kg1
K1
]
D chemical molecular diffusivity [m2
s1
]
g acceleration due to gravity [ms2
]
Gm modified Grashof number [-]
Gr Grashof number [-]
M Hartmann number (magnetic parameter) [-]
K permeability of the porous medium [m2
]
K*
permeability parameter [-]
Pr Prandlt number [-]
p pressure [mmHg]
 rotation parameter [-]
Sc Schmidt number [-]
 temperature [K]
T* dimensional temperature
t time [S]
t* dimensional time
U* dimensional velocity
Greek symbols
 volumetric coefficient of thermal expansion [K1
]
*
volumetric coefficient of expansion with concentration [Kg1
]
 (0 <  < 1) a constant [-]
 thermal conductivity, [J.m1
s1
K1
]
¥
Corresponding Author. Email: sahin.ahmed@rgu.ac.in
Analytical solution for Transient MHD…
www.ijesi.org 34 | Page
 kinematic viscosity [m2
s1
]
 density [Kgm3
]
 species concentration [Kg]
 frequency fo oscillation of the plate temperature, [s1
]
* dimensionless frequency
* angular velocity of the rotating frame of reference, [s1
]
Superscript
U derivative of U with respect to z
Subscripts
w conditions on the porous plane surface
 conditions away from the porous plane surface
I. Introduction
The study of flow problems, which involve the interaction of several phenomena, has a wide range of
applications in the field of science and technology. One such study is related to the effects of MHD free
convection flow, which plays an important role in agriculture, engineering and petroleum industries. The
problem of free convection under the influence of magnetic field has attracted the interest of many researchers
in view of its application in geophysics and astrophysics. In view of these applications, Eckert and Drake (1972)
have done pioneering work on heat and mass transfer. Elbashbeshy (1997) studied heat and mass transfer along
a vertical plate under the combined buoyancy effects of thermal and species diffusion, in the presence of the
magnetic field. Helmy (1998) presented the effects of magnetic field for an unsteady free convective flow past a
vertical porous plate. Soundalgekar (1982) analyzed the problem of free convection effects on flow past a
vertical uniformly accelerated vertical plate under the action of transversely applied magnetic field with mass
transfer. Kim (2000) investigated unsteady MHD convective heat transfer past a semi-infinite vertical porous
moving plate with variable suction by assuming that the free stream velocity follows the exponentially
increasing small perturb action law. The analytical solution of heat and mass transfer on the free convective
flow of a viscous incompressible fluid past an infinite vertical porous plate in presence of transverse sinusoidal
suction velocity and a constant free stream velocity was presented by Ahmed (2009). Also, Ahmed and Liu
(2010) analyzed the effects of mixed convection and mass transfer of three-dimensional oscillatory flow of a
viscous incompressible fluid past an infinite vertical porous plate in presence of transverse sinusoidal suction
velocity oscillating with time and a constant free stream velocity by use of classical perturbation technique.
Convective heat/mass transfer flow in a saturated porous medium has gained growing interest. This fact
has been motivated by its importance in many engineering applications such as building thermal insulation,
geothermal systems, food processing and grain storage, solar power collectors, contaminant transport in
groundwater, casting in manufacturing processes, drying processes, nuclear waste, just to name a few. A
theoretical and experimental work on this subject can be found in the recent monographs by Ingham and Pop
(1998) and Nield and Bejan (1998). Suction/blowing on convective heat transfer over a vertical permeable
surface embedded in a porous medium was analyzed by Cheng (1977). In that work an application to warm
water discharge along the well or fissure to an aquifer of infinite extent is discussed. Kim and Vafai (1989) have
analyzed the buoyancy driven flow about a vertical plate for constant wall temperature and heat flux. Raptis and
Singh (1985) studied flow past an impulsively started vertical plate in a porous medium by a finite difference
method. Ahmed (1983) investigated the effect of transverse periodic permeability oscillating with time on the
heat transfer flow of a viscous incompressible fluid through a highly porous medium bounded by an infinite
vertical porous plate, by means of series solution method. Raptis (1985) analyzed the unsteady flow through a
porous medium bounded by an infinite porous plate subjected to a constant suction and variable temperature.
Raptis and Perdikis (2008) further studied the problem of free convective flow through a porous medium
bounded by a vertical porous plate with constant suction when the free stream velocity oscillates in time about a
constant mean value. Ahmed (2010) studied the effect of transverse periodic permeability oscillating with time
on the free convective heat transfer flow of a viscous incompressible fluid through a highly porous medium
bounded by an infinite vertical porous plate subjected to a periodic suction velocity.
Rotating flow of electrically conducting viscous incompressible fluids has gained considerable
attention because of its numerous applications in physics and engineering which are directly governed by the
action of Coriolis and magnetic forces. In geophysics it is applied to measure and study the positions and
velocities with respect to a fixed frame of reference on the surface of earth which rotate with respect to an
inertial frame in the presence of its magnetic field. In the recent years a number of studies have appeared in the
literature involving rotation to a greater or lesser extent viz. Gupta (1972a), Jana and Datta (1977), Singh
[2000,1999]. Injection/suction effects have also been studied extensively for horizontal porous plate in rotating
frame of references by Ganapathy (1994), Mazumder (1991), Mazumder et al. (1976), Soundalgekar and Pop
(2005) for different physical situation. Singh and Sharma (2001) studied the effect of the permeability of the
Analytical solution for Transient MHD…
www.ijesi.org 35 | Page
porous medium on the three dimensional Couette flow and heat transfer. Recently, Ahmed and Joaquin (2011)
investigated the effects of Hall current, magnetic field, rotation of the channel and suction-injection on the
oscillatory free convective MHD flow in a rotating vertical porous channel when the entire system rotates about
an axis normal to the channel plates and a strong magnetic field of uniform strength is applied along the axis of
rotation.
In the present paper an attempt has been made to study the effects of the permeability of the porous
medium and rotation of the system on the free convective heat and mass transfer flow through a highly porous
medium when the temperature of the surface varies with time about a non-zero constant mean and the
temperature at the free stream is constant. The entire system rotates about an axis perpendicular to the planes of
the plates. Such flows are very important in geophysical and astrophysical problems.
II. Mathematical Formulation
The unsteady flow model of a viscous incompressible fluid through a porous medium occupying a semi-infinite
region of the space bounded by a vertical infinite porous surface in a rotating system under the action of a
uniform magnetic field applied normal to the direction of flow has been analyzed. The temperature of the
surface varies with time about a non-zero constant mean and the temperature at the free stream is constant. The
porous medium is, in fact, a non-homogenous medium which may be replaced by a homogenous fluid having
dynamical properties equal to those of a non-homogenous continuum. Also, we assume that the fluid properties
are not affected by the temperature and concentration differences except by the density  in the body force term;
the influence of the density variations in the momentum and energy equations is negligible.
Fig. 1: Sketch of the physical problem
The vertical infinite porous plate rotates in unison with a viscous fluid occupying the porous region with the
constant angular velocity * about an axis which is perpendicular to the vertical plane surface. The Cartesian
coordinate system is chosen such that X*, Y* axes respectively are in the vertical upward and perpendicular
directions on the plane of the vertical porous surface z = 0 while Z*-axis is normal to it as shown in Fig. 1 with
the above frame of reference and assumptions, the physical variables, except the pressure p, are functions of z
and time t only. Consequently, the equations expressing the conservation of mass, momentum and energy and
the equation of mass transfer, neglecting the heat due to viscous dissipation which is valid for small velocities,
are given by
0
z
W
*
*



(1)
Analytical solution for Transient MHD…
www.ijesi.org 36 | Page
ρ
uσB
u
K
ν
z
u
ν)C(Cgβ)T(TgβvΩ2
z
u
W
t
u
*2
0*
*2*
*2
*******
*
*
*
*
*










(2)
ρ
vσB
v
K
ν
z
v
νuΩ2
z
v
W
t
v
*2
0*
*2*
*2
**
*
*
*
*
*









(3)
*
**
*
W
K
ν
z
p
ρ
1
0 


 (4)
2*
*2
p
*
*
*
*
*
z
T
ρC
k
z
T
W
t
T








(5)
2*
*2
*
*
*
*
*
z
C
D
z
C
W
t
C








(6)
with the boundary conditions
*
w
*ωti**
w
*
w
***
CC,)eT(TεTT0,v0,u  
at z*
= 0
 
*******
zasCC,TT0,v,u (7)
where all the symbols are defined in the Nomenclature section.
In a physically realistic situation, we cannot ensure perfect insulation in any experimental setup. There
will always be some fluctuations in the temperature. The plate temperature is assumed to vary harmonically
with time. It varies from )T(TεT
**
w
*
w 
 as t varies from 0 to 2/. Since  is small, the plate temperature
varies only slightly from the mean value
*
w
T .
For constant suction, we have from Eq. (1) in view of (7)
0
*
WW  (8)
Considering
***
Uivu  and taking into account Eq. (8), then Eqs. (2) and (3) can be written as
*
*2*
*2
*******
*
*
0*
*
U
K
ν
z
U
ν)C(Cg β)Tβ(TgUiΩ2
z
U
W
t
U










(9)
Let us introduce the following non-dimensional quantities:
2
0
*
3
0
**
w
*
3
0
**
w
2
*2
0p
**
w
**
**
w
**
2
0
**2
0
0
**
0
W
νΩ
Ω,
W
)C(Cνgβ
Gm,
W
)Tνgβ(T
Gr,
ν
KW
K,
ρνC
Pr
,
D
ν
Sc,
CC
CC
,
TT
TT
,
ω
νω
ω,
ν
ωt
t,
W
U
U,
ν
ZW
Z



















In view of the above non-dimensional quantities, eqs. (9), (5) and (6) reduce, respectively, to
 UMK
z
U
GmGrUi2
z
U
t
U 21
2
2








 
 (10)
2
2
zPr
1
t 






 
z
(11)
2
2
zSc
1
zt 






 
(12)
and the boundary conditions (7) become


zas0C0,T0,U
0zat1C,eε1T0,U
ωti
(13)
III. Method of Solution
In order to reduce the system of partial differential equations (10)(12) under their boundary conditions (13), to
a system of ordinary differential equations in the non-dimensional form, we assume the following for velocity,
temperature and concentration of the flow field as the amplitude  (<< 1) of the permeability variations is very
small:
Analytical solution for Transient MHD…
www.ijesi.org 37 | Page
(z)Ueε(z)Ut)U(z, 1
ωti
0

(z)Teε(z)Tt)T(z, 1
ωti
0
 (14)
(z)Ceε(z)Ct)C(z, 1
ωti
0

Substituting (14) into the system (10)(12) and equating harmonic and non-harmonic terms we get:
  )Gm(GrUMK2iRUUU 000
21
000
 

(15)
  )Gm(GrU2R)i( ωMKUU 111
21
11
 

(16)
0Pr 00
  (17)
0PriωPr 111
  (18)
0S 00
  c (19)
0SciωSc 111
  (20)
The appropriate boundary conditions reduce to
1(0)1,(0)0,(0)U 000
 
0(0)1,(0)0,(0)U 111
 
0)(0,)(0,)(U 000
  (21)
0)(0,)(0,)(U 111
 
The solutions of the eqs. (15) to (20) subject to the boundary conditions (21) are:
))((
)(
GrεAAAt)U(z,
5141
zz
tiωz
3
Scz
2
Prz
1
15
3






 ee
eeee (22)
ztiωPrz 1
eeεet)(z,

 

(23)
Scz
et)(z,

 (24)
Now, it is convenient to write the primary and secondary velocity fields, in terms of the fluctuating parts,
separating the real and imaginary part from Eqs. (22) and (23) and taking only the real parts as they have
physical significance, the velocity and temperature distribution of the flow field can be expressed in fluctuating
parts as given below.
t)ωsinNtωcos(Nεu
w
u
ir0
0
 (25)
t)ωcosNtωsin(Nε
w
v
ir0
0
 v (26)
where 000
Uivu  and 1ir
UiNN  .
Hence, the expressions for the transient velocity profiles for
2
π
tω  are given by
(z)εN(z)u
2 ω
π
z,
w
u
i0
0






and
(z)εN(z)v
2 ω
π
z,
w
v
r0
0






.
Skin Friction
The skin friction at the plate z = 0 is given by
)(
Gr
εe)AScAA(
dZ
dU
eε
dZ
dU
dZ
dU
41
tωi
3321
0
1ti
0
0
0






Pr
zzz
(27)
Rate of Heat Transfer
The heat transfer coefficient in terms of the Nusselt number at the plate z = 0 is given by
Analytical solution for Transient MHD…
www.ijesi.org 38 | Page
1
ti
0
1ti
0
0
0
eε
dZ
d
eε
dZ
d
dZ
d

 


PrNu
zzz
(28)
The constants A1 to A3, 1 to 5, Ni, Nr have been presented in the Appendix section.
IV. Results and Discussion
The problem of unsteady MHD free convective flow with heat and mass transfer effects in a rotating porous
medium has been considered. The solutions for primary and secondary velocity field, temperature field and
concentration profiles are obtained using the perturbation technique. The effects of flow parameters such as the
magnetic parameter M, porosity parameter K, and the rotation parameter  on the transient velocity profiles
u/w0 and v/w0 have been studied analytically and presented with the help of Figs. 2a-2c and 3a-3b. The effects
of flow parameter on concentration profiles have been presented with the help of Fig. 4. Further, the effects of
flow parameters on the rate of heat transfer have been discussed with the help of Table 1.
Primary Velocity Profile (u/w0)
From equations (22), (23) and (24), it is observed that the steady part of the velocity field has a three layer
character. These layers may be identified as the thermal layer arising due to interaction of the thermal field and
the velocity field and is controlled by the Prandtl number; the concentration layer arising due to the interaction
of the concentration field and the velocity field, and the suction layer as modified by the rotation and the
porosity of the medium. On the other hand, the oscillatory part of the velocity field exhibits a two layer
character. These layers may be identified as the modified suction layers, arising as a result of the triangular
interaction of the Coriolis force and the unsteady convective forces with the porosity of medium.
With a rise in M, from 0, 5, 10 through 15 to 20 there is a strong reduction in magnitude of the primary
velocity across the boundary layer. A velocity peak arises close to the wall for all profiles. The magnetic
parameter is found to decelerate the primary velocity of the flow field to a significant amount due to the
magnetic pull of the Lorentizian hydromagnetic drag force acting on the flow field. With constant M value (=
5.0), as the rotation parameter  increases from 0.0 through 5.0, 10.0, 15.0 to 20.0, there is a distinct depression
in velocity field i.e. the flow is decelerated. Conversely, an increase in porosity parameter K from 0.1, 0.3
through 0.5 to 1.0 the flow velocity accelerated near the plate and then steadily reduces to zero in the free
stream. Irrespective of M or  or K value, it is important to highlight that there is a flow reversal i.e. backflow
has been seen in the boundary layer regime. The velocity, u, sustains negative values throughout motion in
boundary layer.
Figure 2a: Effects of M on primary velocity profiles against Z
-0.08
-0.07
-0.06
-0.05
-0.04
-0.03
-0.02
-0.01
0
0 2 4 6 8 10
u/w0
z
M=0
M=5
M=10
M=15
M=20
Pr=0.71, Gr=5=Gm, =5,
t=2, K=0.5, =10, =0.005
Analytical solution for Transient MHD…
www.ijesi.org 39 | Page
Figure 2b: Effects of  on primary velocity profiles against Z
Figure 2c: Effects of K on primary velocity profiles against Z
Secondary Velocity Profile (v/w0)
The Secondary velocity profiles are shown in Fig. 3a-3b for various values of the magnetic field
(Hartmann number) and rotation parameter. This is indeed the trend observed in Fig. 3a where a strong
deceleration in the flow is achieved with an increase in M from 0.0 (non-conducting case i.e. Lorentz force
vanishes) through 1.0, 2.0, 3.0 to 5.0. It is observed from Fig. 3b that the secondary velocity profiles increases
when there is an increase in rotation parameter of the system. Significantly, the velocity profiles show the
opposite trend to the primary velocity (Fig. 2b) when there is an increase in the rotation parameter. In no case
however there is flow reversal i.e. velocities remain positive through the boundary layer.
Figure 3a: Effects of M on secondary velocity profiles against Z
-0.09
-0.06
-0.03
0
0 2 4 6 8 10
u/w0
z
0
5
10
15
20
Pr=0.71, Gr=5=Gm, M=5,
t=2, K=0.5, =10, =0.005
-0.08
-0.07
-0.06
-0.05
-0.04
-0.03
-0.02
-0.01
0
0 2 4 6 8 10
u/w0
z
K=0.1
K=0.3
K=0.5
K=1.0
Pr=0.71, Gr=5=Gm, =5,
t=2, M=5, =10, =0.005
0
0.02
0.04
0.06
0.08
0.1
0.12
0.14
0.16
0.18
0.2
0.22
0 2 4 6 8 10
v/w0
z
M=0
M=1
M=2
M=3
M=5
Pr=0.71, Gr=5=Gm, =5,
t=2, K=0.5, =10, =0.005
Analytical solution for Transient MHD…
www.ijesi.org 40 | Page
Figure 3b: Effects of  on secondary velocity profiles against Z
Concentration Profiles ()
In Fig. 4, the profiles of concentration () for various values of the Schmidt number are presented. The values
of the Schmidt number Sc are chosen in such a way that they represent the diffusing chemical species of most
common interest in air. For example, the values of Sc for H2, H20, NH3, propyl benzene and helium in air are
0.22, 0.60, 0.78, 2.62 and 0.30, respectively as reported by Perry (1963). It is seen that the concentration profiles
decreases with an increase in the Schmidt number. It is noted that for heavier diffusing foreign speices, i.e.,
increasing the Schmidt number reduces the concentration in both magnitude and extent and thining of
concentration boundary layer occurs.
Figure 4: Effects of Sc on concentration distribution
Skin Friction coefficient
The distribution of skin friction coefficients are presented in Figs. 4a-4b for different values of the magnetic
field, rotation parameter and porosity parameter. Inspection shows that increasing rotation () reduces skin
friction. Increasing porosity (K) is found to increase the skin friction, in consistency with earlier discussion for
the velocity response (Fig. 2c). Flow reversal is observed i.e. skin friction becomes negative in both the figures
4a and 4b. Clearly all profiles decay as M increases since larger Hartmann number corresponds to greater
magnetohydrodynamic drag which decelerates the flow.
Figure 4a: Effects of  on skin friction against M
-6.25E-1
0.015
0.03
0.045
0 2 4 6 8 10
v/w0
z
5
10
20
15
Pr=0.71, Gr=5=Gm, =5,
t=2, K=0.5, =10, =0.005
0
0.2
0.4
0.6
0.8
1
0 1 2 3 4 5

Z
Sc=0.60
Sc=0.30
Sc=0.78
Sc=2.62
-6
-5
-4
-3
-2
-1
0
0 5 10 15 20

M
0
5
10Pr=0.71, Gr=5=Gm, K=0.5,
t=2, =10, =0.005
Analytical solution for Transient MHD…
www.ijesi.org 41 | Page
Figure 4b: Effects of K on skin friction against M
Rate of Heat Transfer (Nu)
Table 1: Heat transfer coefficient when  = 0.005, in the case air and water:
Pr = 0.71
(Air)
Pr = 7.0
(Water)
 t M Nu Nu
2 2 2 0.34410 3.71514
2 2 3 0.37135 3.80426
3 2 3 0.80371 4.37293
3 3 3 0.03704 2.10853
3 3 4 0.02472 2.05830
4 3 4 0.90514 4.43715
4 4 4 0.07632 2.35304
Table 3 shows the magnitude of the heat transfer coefficient for different values of , t and M. An
increase either in the Hartmann number (M) or time (t) the magnitude of the heat transfer coefficient reduces for
both the cases of air and water. However, the heat transfer coefficient rises due to an increase in the frequency
parameter (). Significantly, it is observed that the heat transfer coefficient in the case of water for any
particular values of , t and M is much higher than that of air.
V. Conclusion
The present theoretical analysis brings out the following results of physical interest on the primary and
secondary velocities, concentration profiles, skin friction and rate of heat transfer of the flow field for a rotating
system in a saturated porous regime:
 It is noticed that all the velocity profiles increase steadily near the plate and thereafter they show a
constant decrease and reach the value zero at the free stream.
 The magnetic / rotation parameter is found to decelerate the primary velocity of the flow field.
 The porosity drag force has an effect accelerating the primary velocity profiles and reduces the skin
friction.
 There is a clear back flow for all the primary velocity profiles.
 The secondary velocity profiles and skin friction coefficient are opposite to each other due to the effect
of rotation.
 The magnetic parameter has the effect of decreasing the heat transfer profiles whereas it has the effect
of increasing the skin friction.
 The presence of foreign species reduces the concentration boundary layer and further reduction occurs
with increasing values of the Schmidt number.
 The Prandtl number and the frequency parameter have the effect of increasing the heat transfer
coefficient.
Appendix:
))(Pr(Pr
Gr
A
32
1
 
 ,
))(Sc(Sc
Gm
A
32
2
 
 , )A(AA 213
 ,
-2.5
-2
-1.5
-1
-0.5
0 5 10 15 20

M
K=0.1
K=0.3
K=0.5
K=1.0
Pr=0.71, Gr=5=Gm, =5,
t=2, =10, =0.005
Analytical solution for Transient MHD…
www.ijesi.org 42 | Page
2
)Pr4i ωPrPr(
2
1

 ,
2
K
1
M2i411
2
2















 ,
2
K
1
M2i411
2
3















 ,
2
iω
K
1
M2i411
2
4















 ,
2
iω
K
1
M2i411
2
5















 ,  










zz
15
tiω
r
15
ee
Gr
epartRealN


,
  .ee
Gr
epartImaginaryN
zz
15
tiω
i
15













References
[1] Eckert, E.R., and Drake, R.M., 1972, “Analysis of Heat and Mass Transfer,” McGraw-Hill, pp. 171.
[2] Elbashbeshy, E.M.A., 1997, “Heat and Mass Transfer along a vertical plate with variable surface tension and concentration in the
presence of the magnetic field,” Int. J. Engineering Science, 34(5), pp. 515-522.
[3] Helmy, K.A., 1998, “MHD Unsteady free convection flow past a vertical porous plate,” ZAMM, 78, pp. 255-270.
[4] Soundalgekar, V.M., 1982, “Effects of mass transfer on flow past a uniformly accelerated vertical plate,” Letters in heat and
mass transfer, 9, pp.65-72.
[5] Kim, Y.J., 2000, “Unsteady MHD convective heat transfer past semi-infinite vertical porous moving plate with variable suction,”
International Journal of Engineering Science, 38, pp. 833-845.
[6] Ahmed, S., 2009, “The study of heat and mass transfer on free convective three-dimensional unsteady flows over a porous
vertical plate,” Journal of Energy Heat and Mass Transfer, 31, pp. 89-110.
[7] Ahmed, S., and Liu, I.C., 2010, Mixed convective three-dimensional heat and mass transfer flow with transversely periodic
suction velocity,” Int. J. Applied Mathematics and Mechanics, 6(1), pp. 58-73.
[8] Ingham, D.B., and Pop, I., 1998, Transport Phenomena in Porous Media, Pergamon, Oxford.
[9] Nield, D.A. and Bejan, A., 1998, Convection in Porous Media, 2. Ed., Springer-Verlag, Berlin.
[10] Cheng, P., 1977, “The influence of lateral mass flux on free convection boundary layers in a saturated porous medium,” Int.J.
Heat Mass Transfer, 20, pp. 201- 206.
[11] Kim, S., and Vafai K. 1989. Analysis of natural convection about a vertical plate embedded in a porous medium, Int. J. Heat
Mass Transfer, vol. 32, pp. 665-677.
[12] Raptis, A., and Singh, A.K., 1985, “Free convection flow past an impulsively started vertical plate in a porous medium by finite
difference method,” Astrophys. Space Sci., 112, pp. 259-265.
[13] Raptis, A.A., 1983, “Unsteady free convection through porous medium,” Int. J. Engg. sci., 21, pp. 345-348.
[14] Raptis, A.A., and Perdikis, C.P., 1985, “Oscillatory flow through a porous medium by the presence of free convective flow,” Int.
J. Engg. sci., 23, pp.51-55.
[15] Ahmed, S., 2008, Transient three dimensional flow through a porous medium with transverse permeability oscillating with
time,” Emirate Journal for Engineering research, 13(3), pp. 11-17.
[16] Ahmed, S., 2010, “Free convective transient three-dimensional flow through a porous medium oscillating with time in presence
of periodic suction velocity,” Int. J. Applied Mathematics and Mechanics, 6(11), pp. 1-16.
[17] Gupta, A.S., 1972a, “Ekman layer on a porous plate,” Phys. Fluids, 15, pp. 930-931.
[18] Jana, R.N., and Datta, N., 1977, “Couette flow and heat transfer in a rotating system,” Acta Mech., 26, pp. 301-306.
[19] Singh, K.D., 2000, “An oscillatory hydromagnetic Couette flow in a rotating system,” J. Appl. Math and Mech., 80, pp. 429-432.
[20] Singh, K.D., 1999, “Three dimensional Coutte flow with transpiration cooling,” Z. Angew. Math. Phys., 50, pp. 661-668.
[21] Ganapathy, R., 1994, “A note on oscillatory Couette flow in a rotating system,” ASME J. Appl. Mech., 61, pp. 208-209.
[22] Mazumder, B.S., 1991, “An exact solution of oscillatory Couette flow in a rotating system,” ASME J. Appl. Mech., 58, pp. 1104-
1107.
[23] Mazumder, B.S., Gupta, A.S., and Dutta, N., 1976, “Flow and heat transfer in the hydromagnetic Ekman layer on a porous plate
with Hall effects,” Int. J. Heat Mass Transfer, 19, pp. 523-527.
[24] Singh, K.D., Gorla, M.G., and Hans, Raj., 2005, “A periodic solution of oscillatory Couette flow through porous medium in
rotating system,” Indian J. pure appl. Math., 36, no. 3, pp. 151-159.
[25] Singh, K.D., and Sharma, R., 2001, “Three dimensional Couette flow through a porous medium with heat transfer,” Indian J.
pure appl. Math., 32(12), pp. 1819-1829.
[26] Ahmed, S., and Zueco, J., 2011, “Modelling of heat and mass transfer in a rotating vertical porous channel with hall current,”
Chemical Engineering Communications, 198(10), pp. 1294-1308.
[27] Perry, E.D., 1963, Chemical Engineers Handbook, (4th
Edn.), McGraw-Hill Book Company, New York.

More Related Content

PDF
D04702040055
PDF
Effects of Thermal Radiation and Chemical Reaction on MHD Free Convection Flo...
PDF
B0412011026
PDF
MHD Stagnation Point Flow of A MICROPOLAR Fluid over a Stretching Surface wit...
PDF
Radiation and Soret Effect on Visco-Elastic MHD Oscillatory Horizontal Channe...
PDF
Effect of Radiation on Mixed Convection Flow of a Non-Newtonian Nan fluid ove...
PDF
Do36689702
PDF
Soret Effect And Effect Of Radiation On Transient Mhd Free Convective Flow Ov...
D04702040055
Effects of Thermal Radiation and Chemical Reaction on MHD Free Convection Flo...
B0412011026
MHD Stagnation Point Flow of A MICROPOLAR Fluid over a Stretching Surface wit...
Radiation and Soret Effect on Visco-Elastic MHD Oscillatory Horizontal Channe...
Effect of Radiation on Mixed Convection Flow of a Non-Newtonian Nan fluid ove...
Do36689702
Soret Effect And Effect Of Radiation On Transient Mhd Free Convective Flow Ov...

What's hot (20)

PDF
Unsteady MHD Flow Past A Semi-Infinite Vertical Plate With Heat Source/ Sink:...
PDF
L044047379
PDF
Heat Transfer on Steady MHD rotating flow through porous medium in a parallel...
PDF
first order chemical reaction and hall effect on mhd flow past an infinite ve...
PDF
F0733038
PDF
Gv3311971210
PDF
Effect of Thermo-Diffusion and Chemical Reaction on Mixed Convective Heat And...
PDF
Effect of Chemical Reaction and Radiation Absorption on Unsteady Convective H...
PDF
NUMERICAL INVESTIGATION OF AN UNSTEADY MIXED CONVECTIVE MASS AND HEAT TRANSFE...
PDF
C027013025
PDF
Ac32195204
PDF
Double Diffusive Convection and the Improvement of Flow in Square Porous Annulus
PDF
A045001011
PDF
E0362232
PDF
Effects of Variable Fluid Properties and MHD on Mixed Convection Heat Transfe...
PDF
28Dufour effect on radiative effect flow and heat transfer over a vertically ...
PDF
Effect of Radiation and Thermo-Diffusion on Convective Heat and Mass Transfer...
PDF
MHD Flow past a Vertical Oscillating Plate with Radiation and Chemical Reacti...
PDF
Hydromagnetic Flow of Nanofluids over a Nonlinear Shrinking Surface with Heat...
PDF
The International Journal of Engineering and Science (IJES)
Unsteady MHD Flow Past A Semi-Infinite Vertical Plate With Heat Source/ Sink:...
L044047379
Heat Transfer on Steady MHD rotating flow through porous medium in a parallel...
first order chemical reaction and hall effect on mhd flow past an infinite ve...
F0733038
Gv3311971210
Effect of Thermo-Diffusion and Chemical Reaction on Mixed Convective Heat And...
Effect of Chemical Reaction and Radiation Absorption on Unsteady Convective H...
NUMERICAL INVESTIGATION OF AN UNSTEADY MIXED CONVECTIVE MASS AND HEAT TRANSFE...
C027013025
Ac32195204
Double Diffusive Convection and the Improvement of Flow in Square Porous Annulus
A045001011
E0362232
Effects of Variable Fluid Properties and MHD on Mixed Convection Heat Transfe...
28Dufour effect on radiative effect flow and heat transfer over a vertically ...
Effect of Radiation and Thermo-Diffusion on Convective Heat and Mass Transfer...
MHD Flow past a Vertical Oscillating Plate with Radiation and Chemical Reacti...
Hydromagnetic Flow of Nanofluids over a Nonlinear Shrinking Surface with Heat...
The International Journal of Engineering and Science (IJES)
Ad

Similar to Analytical solution for Transient MHD flow through a Darcian porous regime in a Rotating System (20)

PDF
The Effect of Radiation on the Convective Heat and Mass Transfer Flow of a Vi...
PDF
Minimal M-gs Open and Maximal M-gs Closed Sets In Interior Minimal Space
PDF
E0325022040
PDF
Finite difference solutions of magneto hydrodynamic free convective flow with...
PDF
I05844759
PDF
Ac25149159
PDF
MHD Free Convection and Mass Transfer Flow past a Vertical Flat Plate
PDF
Soret Effect And Effect Of Radiation On Transient Mhd Free Convective Flow Ov...
PDF
Dufour and Soret Effects on Convective Heat and Mass Transfer in Non-Darcy Do...
PDF
D04841832
PDF
MHD Convection Flow of Kuvshinski Fluid past an Infinite Vertical Porous Plat...
PDF
Laplace Analysis of Periodic Heat and Mass Transport on a Parabolic Started S...
PDF
asianjournal.pdf
PDF
Radiation Effects on MHD Free Convective Rotating Flow with Hall Effects
PDF
International Journal of Computational Engineering Research(IJCER)
PDF
A0350300108
PDF
Effects on Study MHD Free Convection Flow Past a Vertical Porous Plate with H...
PDF
P0748292
PDF
OSCILLATORY FLOW OF MHD POLAR FLUID WITH HEAT AND MASS TRANSFER PAST A VERTIC...
The Effect of Radiation on the Convective Heat and Mass Transfer Flow of a Vi...
Minimal M-gs Open and Maximal M-gs Closed Sets In Interior Minimal Space
E0325022040
Finite difference solutions of magneto hydrodynamic free convective flow with...
I05844759
Ac25149159
MHD Free Convection and Mass Transfer Flow past a Vertical Flat Plate
Soret Effect And Effect Of Radiation On Transient Mhd Free Convective Flow Ov...
Dufour and Soret Effects on Convective Heat and Mass Transfer in Non-Darcy Do...
D04841832
MHD Convection Flow of Kuvshinski Fluid past an Infinite Vertical Porous Plat...
Laplace Analysis of Periodic Heat and Mass Transport on a Parabolic Started S...
asianjournal.pdf
Radiation Effects on MHD Free Convective Rotating Flow with Hall Effects
International Journal of Computational Engineering Research(IJCER)
A0350300108
Effects on Study MHD Free Convection Flow Past a Vertical Porous Plate with H...
P0748292
OSCILLATORY FLOW OF MHD POLAR FLUID WITH HEAT AND MASS TRANSFER PAST A VERTIC...
Ad

Recently uploaded (20)

PPTX
Digital-Transformation-Roadmap-for-Companies.pptx
PDF
Encapsulation theory and applications.pdf
PDF
Architecting across the Boundaries of two Complex Domains - Healthcare & Tech...
PDF
Blue Purple Modern Animated Computer Science Presentation.pdf.pdf
PPTX
Programs and apps: productivity, graphics, security and other tools
PDF
Advanced methodologies resolving dimensionality complications for autism neur...
PDF
gpt5_lecture_notes_comprehensive_20250812015547.pdf
PDF
Machine learning based COVID-19 study performance prediction
PPTX
MYSQL Presentation for SQL database connectivity
PDF
Diabetes mellitus diagnosis method based random forest with bat algorithm
PDF
Dropbox Q2 2025 Financial Results & Investor Presentation
PPTX
ACSFv1EN-58255 AWS Academy Cloud Security Foundations.pptx
PDF
Reach Out and Touch Someone: Haptics and Empathic Computing
PDF
Approach and Philosophy of On baking technology
PDF
Mobile App Security Testing_ A Comprehensive Guide.pdf
PDF
Profit Center Accounting in SAP S/4HANA, S4F28 Col11
PDF
Agricultural_Statistics_at_a_Glance_2022_0.pdf
PPTX
20250228 LYD VKU AI Blended-Learning.pptx
PDF
Empathic Computing: Creating Shared Understanding
PPTX
Spectroscopy.pptx food analysis technology
Digital-Transformation-Roadmap-for-Companies.pptx
Encapsulation theory and applications.pdf
Architecting across the Boundaries of two Complex Domains - Healthcare & Tech...
Blue Purple Modern Animated Computer Science Presentation.pdf.pdf
Programs and apps: productivity, graphics, security and other tools
Advanced methodologies resolving dimensionality complications for autism neur...
gpt5_lecture_notes_comprehensive_20250812015547.pdf
Machine learning based COVID-19 study performance prediction
MYSQL Presentation for SQL database connectivity
Diabetes mellitus diagnosis method based random forest with bat algorithm
Dropbox Q2 2025 Financial Results & Investor Presentation
ACSFv1EN-58255 AWS Academy Cloud Security Foundations.pptx
Reach Out and Touch Someone: Haptics and Empathic Computing
Approach and Philosophy of On baking technology
Mobile App Security Testing_ A Comprehensive Guide.pdf
Profit Center Accounting in SAP S/4HANA, S4F28 Col11
Agricultural_Statistics_at_a_Glance_2022_0.pdf
20250228 LYD VKU AI Blended-Learning.pptx
Empathic Computing: Creating Shared Understanding
Spectroscopy.pptx food analysis technology

Analytical solution for Transient MHD flow through a Darcian porous regime in a Rotating System

  • 1. International Journal of Engineering Science Invention ISSN (Online): 2319 – 6734, ISSN (Print): 2319 – 6726 www.ijesi.org ||Volume 4 Issue 10|| October 2015 || PP.33-42 www.ijesi.org 33 | Page Analytical solution for Transient MHD flow through a Darcian porous regime in a Rotating System 1 Hamida Khatun , ¥ Sahin Ahmed 1 Department of Mathematics, South Salmara College, Dhubri-783127, Assam, India, 2 Department of Mathematics, Rajiv Gandhi Central University, Rono Hills, Itanagar, Arunachal Pradesh-791112, India, Abstract: A rotating model is developed for a two-dimensional, unsteady, incompressible electrically conducting, laminar free convection boundary layer flow of heat and mass transport in a saturated porous medium, bounded by an infinite vertical porous surface in presence of an applied transverse magnetic field. The porous plane surface and the porous medium are assumed to rotate in a solid body rotation. The vertical surface is subjected to uniform constant suction perpendicular to it and the temperature at this surface fluctuates in time about a non-zero constant mean. The basic equations governing the flow are in the form of partial differential equations and have been reduced to a set of non-linear ordinary differential equations by applying suitable similarity transformations. The problem is tackled analytically using classical perturbation technique. Pertinent results with respect to embedded parameters are displayed through graphically and tables for the velocity, concentration and skin friction profiles and were discussed quantitatively. Applications to the flows of fluids through porous medium bounded by rotating porous systems find many industrial applications particularly in the fields of centrifugation, filtration and purification processes. Keywords: Coriolis and magnetic forces; Mass transport; free convection; porous medium; rotating system. Mathematics Subject classification: 76(D, W05, 80A32, M45). NOMENCLATURE u, v, w velocity components in x, y and z directions respectively [ms1 ] w0 (>0) constant suction velocity of the liquid through the porous plane surface [ms1 ] Z normal direction of vertical porous plane surface [m] Z* dimensional normal distance [m] C* dimensional species concentration [-] Cp specific heat at constant pressure [J Kg1 K1 ] D chemical molecular diffusivity [m2 s1 ] g acceleration due to gravity [ms2 ] Gm modified Grashof number [-] Gr Grashof number [-] M Hartmann number (magnetic parameter) [-] K permeability of the porous medium [m2 ] K* permeability parameter [-] Pr Prandlt number [-] p pressure [mmHg]  rotation parameter [-] Sc Schmidt number [-]  temperature [K] T* dimensional temperature t time [S] t* dimensional time U* dimensional velocity Greek symbols  volumetric coefficient of thermal expansion [K1 ] * volumetric coefficient of expansion with concentration [Kg1 ]  (0 <  < 1) a constant [-]  thermal conductivity, [J.m1 s1 K1 ] ¥ Corresponding Author. Email: sahin.ahmed@rgu.ac.in
  • 2. Analytical solution for Transient MHD… www.ijesi.org 34 | Page  kinematic viscosity [m2 s1 ]  density [Kgm3 ]  species concentration [Kg]  frequency fo oscillation of the plate temperature, [s1 ] * dimensionless frequency * angular velocity of the rotating frame of reference, [s1 ] Superscript U derivative of U with respect to z Subscripts w conditions on the porous plane surface  conditions away from the porous plane surface I. Introduction The study of flow problems, which involve the interaction of several phenomena, has a wide range of applications in the field of science and technology. One such study is related to the effects of MHD free convection flow, which plays an important role in agriculture, engineering and petroleum industries. The problem of free convection under the influence of magnetic field has attracted the interest of many researchers in view of its application in geophysics and astrophysics. In view of these applications, Eckert and Drake (1972) have done pioneering work on heat and mass transfer. Elbashbeshy (1997) studied heat and mass transfer along a vertical plate under the combined buoyancy effects of thermal and species diffusion, in the presence of the magnetic field. Helmy (1998) presented the effects of magnetic field for an unsteady free convective flow past a vertical porous plate. Soundalgekar (1982) analyzed the problem of free convection effects on flow past a vertical uniformly accelerated vertical plate under the action of transversely applied magnetic field with mass transfer. Kim (2000) investigated unsteady MHD convective heat transfer past a semi-infinite vertical porous moving plate with variable suction by assuming that the free stream velocity follows the exponentially increasing small perturb action law. The analytical solution of heat and mass transfer on the free convective flow of a viscous incompressible fluid past an infinite vertical porous plate in presence of transverse sinusoidal suction velocity and a constant free stream velocity was presented by Ahmed (2009). Also, Ahmed and Liu (2010) analyzed the effects of mixed convection and mass transfer of three-dimensional oscillatory flow of a viscous incompressible fluid past an infinite vertical porous plate in presence of transverse sinusoidal suction velocity oscillating with time and a constant free stream velocity by use of classical perturbation technique. Convective heat/mass transfer flow in a saturated porous medium has gained growing interest. This fact has been motivated by its importance in many engineering applications such as building thermal insulation, geothermal systems, food processing and grain storage, solar power collectors, contaminant transport in groundwater, casting in manufacturing processes, drying processes, nuclear waste, just to name a few. A theoretical and experimental work on this subject can be found in the recent monographs by Ingham and Pop (1998) and Nield and Bejan (1998). Suction/blowing on convective heat transfer over a vertical permeable surface embedded in a porous medium was analyzed by Cheng (1977). In that work an application to warm water discharge along the well or fissure to an aquifer of infinite extent is discussed. Kim and Vafai (1989) have analyzed the buoyancy driven flow about a vertical plate for constant wall temperature and heat flux. Raptis and Singh (1985) studied flow past an impulsively started vertical plate in a porous medium by a finite difference method. Ahmed (1983) investigated the effect of transverse periodic permeability oscillating with time on the heat transfer flow of a viscous incompressible fluid through a highly porous medium bounded by an infinite vertical porous plate, by means of series solution method. Raptis (1985) analyzed the unsteady flow through a porous medium bounded by an infinite porous plate subjected to a constant suction and variable temperature. Raptis and Perdikis (2008) further studied the problem of free convective flow through a porous medium bounded by a vertical porous plate with constant suction when the free stream velocity oscillates in time about a constant mean value. Ahmed (2010) studied the effect of transverse periodic permeability oscillating with time on the free convective heat transfer flow of a viscous incompressible fluid through a highly porous medium bounded by an infinite vertical porous plate subjected to a periodic suction velocity. Rotating flow of electrically conducting viscous incompressible fluids has gained considerable attention because of its numerous applications in physics and engineering which are directly governed by the action of Coriolis and magnetic forces. In geophysics it is applied to measure and study the positions and velocities with respect to a fixed frame of reference on the surface of earth which rotate with respect to an inertial frame in the presence of its magnetic field. In the recent years a number of studies have appeared in the literature involving rotation to a greater or lesser extent viz. Gupta (1972a), Jana and Datta (1977), Singh [2000,1999]. Injection/suction effects have also been studied extensively for horizontal porous plate in rotating frame of references by Ganapathy (1994), Mazumder (1991), Mazumder et al. (1976), Soundalgekar and Pop (2005) for different physical situation. Singh and Sharma (2001) studied the effect of the permeability of the
  • 3. Analytical solution for Transient MHD… www.ijesi.org 35 | Page porous medium on the three dimensional Couette flow and heat transfer. Recently, Ahmed and Joaquin (2011) investigated the effects of Hall current, magnetic field, rotation of the channel and suction-injection on the oscillatory free convective MHD flow in a rotating vertical porous channel when the entire system rotates about an axis normal to the channel plates and a strong magnetic field of uniform strength is applied along the axis of rotation. In the present paper an attempt has been made to study the effects of the permeability of the porous medium and rotation of the system on the free convective heat and mass transfer flow through a highly porous medium when the temperature of the surface varies with time about a non-zero constant mean and the temperature at the free stream is constant. The entire system rotates about an axis perpendicular to the planes of the plates. Such flows are very important in geophysical and astrophysical problems. II. Mathematical Formulation The unsteady flow model of a viscous incompressible fluid through a porous medium occupying a semi-infinite region of the space bounded by a vertical infinite porous surface in a rotating system under the action of a uniform magnetic field applied normal to the direction of flow has been analyzed. The temperature of the surface varies with time about a non-zero constant mean and the temperature at the free stream is constant. The porous medium is, in fact, a non-homogenous medium which may be replaced by a homogenous fluid having dynamical properties equal to those of a non-homogenous continuum. Also, we assume that the fluid properties are not affected by the temperature and concentration differences except by the density  in the body force term; the influence of the density variations in the momentum and energy equations is negligible. Fig. 1: Sketch of the physical problem The vertical infinite porous plate rotates in unison with a viscous fluid occupying the porous region with the constant angular velocity * about an axis which is perpendicular to the vertical plane surface. The Cartesian coordinate system is chosen such that X*, Y* axes respectively are in the vertical upward and perpendicular directions on the plane of the vertical porous surface z = 0 while Z*-axis is normal to it as shown in Fig. 1 with the above frame of reference and assumptions, the physical variables, except the pressure p, are functions of z and time t only. Consequently, the equations expressing the conservation of mass, momentum and energy and the equation of mass transfer, neglecting the heat due to viscous dissipation which is valid for small velocities, are given by 0 z W * *    (1)
  • 4. Analytical solution for Transient MHD… www.ijesi.org 36 | Page ρ uσB u K ν z u ν)C(Cgβ)T(TgβvΩ2 z u W t u *2 0* *2* *2 ******* * * * * *           (2) ρ vσB v K ν z v νuΩ2 z v W t v *2 0* *2* *2 ** * * * * *          (3) * ** * W K ν z p ρ 1 0     (4) 2* *2 p * * * * * z T ρC k z T W t T         (5) 2* *2 * * * * * z C D z C W t C         (6) with the boundary conditions * w *ωti** w * w *** CC,)eT(TεTT0,v0,u   at z* = 0   ******* zasCC,TT0,v,u (7) where all the symbols are defined in the Nomenclature section. In a physically realistic situation, we cannot ensure perfect insulation in any experimental setup. There will always be some fluctuations in the temperature. The plate temperature is assumed to vary harmonically with time. It varies from )T(TεT ** w * w   as t varies from 0 to 2/. Since  is small, the plate temperature varies only slightly from the mean value * w T . For constant suction, we have from Eq. (1) in view of (7) 0 * WW  (8) Considering *** Uivu  and taking into account Eq. (8), then Eqs. (2) and (3) can be written as * *2* *2 ******* * * 0* * U K ν z U ν)C(Cg β)Tβ(TgUiΩ2 z U W t U           (9) Let us introduce the following non-dimensional quantities: 2 0 * 3 0 ** w * 3 0 ** w 2 *2 0p ** w ** ** w ** 2 0 **2 0 0 ** 0 W νΩ Ω, W )C(Cνgβ Gm, W )Tνgβ(T Gr, ν KW K, ρνC Pr , D ν Sc, CC CC , TT TT , ω νω ω, ν ωt t, W U U, ν ZW Z                    In view of the above non-dimensional quantities, eqs. (9), (5) and (6) reduce, respectively, to  UMK z U GmGrUi2 z U t U 21 2 2            (10) 2 2 zPr 1 t          z (11) 2 2 zSc 1 zt          (12) and the boundary conditions (7) become   zas0C0,T0,U 0zat1C,eε1T0,U ωti (13) III. Method of Solution In order to reduce the system of partial differential equations (10)(12) under their boundary conditions (13), to a system of ordinary differential equations in the non-dimensional form, we assume the following for velocity, temperature and concentration of the flow field as the amplitude  (<< 1) of the permeability variations is very small:
  • 5. Analytical solution for Transient MHD… www.ijesi.org 37 | Page (z)Ueε(z)Ut)U(z, 1 ωti 0  (z)Teε(z)Tt)T(z, 1 ωti 0  (14) (z)Ceε(z)Ct)C(z, 1 ωti 0  Substituting (14) into the system (10)(12) and equating harmonic and non-harmonic terms we get:   )Gm(GrUMK2iRUUU 000 21 000    (15)   )Gm(GrU2R)i( ωMKUU 111 21 11    (16) 0Pr 00   (17) 0PriωPr 111   (18) 0S 00   c (19) 0SciωSc 111   (20) The appropriate boundary conditions reduce to 1(0)1,(0)0,(0)U 000   0(0)1,(0)0,(0)U 111   0)(0,)(0,)(U 000   (21) 0)(0,)(0,)(U 111   The solutions of the eqs. (15) to (20) subject to the boundary conditions (21) are: ))(( )( GrεAAAt)U(z, 5141 zz tiωz 3 Scz 2 Prz 1 15 3        ee eeee (22) ztiωPrz 1 eeεet)(z,     (23) Scz et)(z,   (24) Now, it is convenient to write the primary and secondary velocity fields, in terms of the fluctuating parts, separating the real and imaginary part from Eqs. (22) and (23) and taking only the real parts as they have physical significance, the velocity and temperature distribution of the flow field can be expressed in fluctuating parts as given below. t)ωsinNtωcos(Nεu w u ir0 0  (25) t)ωcosNtωsin(Nε w v ir0 0  v (26) where 000 Uivu  and 1ir UiNN  . Hence, the expressions for the transient velocity profiles for 2 π tω  are given by (z)εN(z)u 2 ω π z, w u i0 0       and (z)εN(z)v 2 ω π z, w v r0 0       . Skin Friction The skin friction at the plate z = 0 is given by )( Gr εe)AScAA( dZ dU eε dZ dU dZ dU 41 tωi 3321 0 1ti 0 0 0       Pr zzz (27) Rate of Heat Transfer The heat transfer coefficient in terms of the Nusselt number at the plate z = 0 is given by
  • 6. Analytical solution for Transient MHD… www.ijesi.org 38 | Page 1 ti 0 1ti 0 0 0 eε dZ d eε dZ d dZ d      PrNu zzz (28) The constants A1 to A3, 1 to 5, Ni, Nr have been presented in the Appendix section. IV. Results and Discussion The problem of unsteady MHD free convective flow with heat and mass transfer effects in a rotating porous medium has been considered. The solutions for primary and secondary velocity field, temperature field and concentration profiles are obtained using the perturbation technique. The effects of flow parameters such as the magnetic parameter M, porosity parameter K, and the rotation parameter  on the transient velocity profiles u/w0 and v/w0 have been studied analytically and presented with the help of Figs. 2a-2c and 3a-3b. The effects of flow parameter on concentration profiles have been presented with the help of Fig. 4. Further, the effects of flow parameters on the rate of heat transfer have been discussed with the help of Table 1. Primary Velocity Profile (u/w0) From equations (22), (23) and (24), it is observed that the steady part of the velocity field has a three layer character. These layers may be identified as the thermal layer arising due to interaction of the thermal field and the velocity field and is controlled by the Prandtl number; the concentration layer arising due to the interaction of the concentration field and the velocity field, and the suction layer as modified by the rotation and the porosity of the medium. On the other hand, the oscillatory part of the velocity field exhibits a two layer character. These layers may be identified as the modified suction layers, arising as a result of the triangular interaction of the Coriolis force and the unsteady convective forces with the porosity of medium. With a rise in M, from 0, 5, 10 through 15 to 20 there is a strong reduction in magnitude of the primary velocity across the boundary layer. A velocity peak arises close to the wall for all profiles. The magnetic parameter is found to decelerate the primary velocity of the flow field to a significant amount due to the magnetic pull of the Lorentizian hydromagnetic drag force acting on the flow field. With constant M value (= 5.0), as the rotation parameter  increases from 0.0 through 5.0, 10.0, 15.0 to 20.0, there is a distinct depression in velocity field i.e. the flow is decelerated. Conversely, an increase in porosity parameter K from 0.1, 0.3 through 0.5 to 1.0 the flow velocity accelerated near the plate and then steadily reduces to zero in the free stream. Irrespective of M or  or K value, it is important to highlight that there is a flow reversal i.e. backflow has been seen in the boundary layer regime. The velocity, u, sustains negative values throughout motion in boundary layer. Figure 2a: Effects of M on primary velocity profiles against Z -0.08 -0.07 -0.06 -0.05 -0.04 -0.03 -0.02 -0.01 0 0 2 4 6 8 10 u/w0 z M=0 M=5 M=10 M=15 M=20 Pr=0.71, Gr=5=Gm, =5, t=2, K=0.5, =10, =0.005
  • 7. Analytical solution for Transient MHD… www.ijesi.org 39 | Page Figure 2b: Effects of  on primary velocity profiles against Z Figure 2c: Effects of K on primary velocity profiles against Z Secondary Velocity Profile (v/w0) The Secondary velocity profiles are shown in Fig. 3a-3b for various values of the magnetic field (Hartmann number) and rotation parameter. This is indeed the trend observed in Fig. 3a where a strong deceleration in the flow is achieved with an increase in M from 0.0 (non-conducting case i.e. Lorentz force vanishes) through 1.0, 2.0, 3.0 to 5.0. It is observed from Fig. 3b that the secondary velocity profiles increases when there is an increase in rotation parameter of the system. Significantly, the velocity profiles show the opposite trend to the primary velocity (Fig. 2b) when there is an increase in the rotation parameter. In no case however there is flow reversal i.e. velocities remain positive through the boundary layer. Figure 3a: Effects of M on secondary velocity profiles against Z -0.09 -0.06 -0.03 0 0 2 4 6 8 10 u/w0 z 0 5 10 15 20 Pr=0.71, Gr=5=Gm, M=5, t=2, K=0.5, =10, =0.005 -0.08 -0.07 -0.06 -0.05 -0.04 -0.03 -0.02 -0.01 0 0 2 4 6 8 10 u/w0 z K=0.1 K=0.3 K=0.5 K=1.0 Pr=0.71, Gr=5=Gm, =5, t=2, M=5, =10, =0.005 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 0.22 0 2 4 6 8 10 v/w0 z M=0 M=1 M=2 M=3 M=5 Pr=0.71, Gr=5=Gm, =5, t=2, K=0.5, =10, =0.005
  • 8. Analytical solution for Transient MHD… www.ijesi.org 40 | Page Figure 3b: Effects of  on secondary velocity profiles against Z Concentration Profiles () In Fig. 4, the profiles of concentration () for various values of the Schmidt number are presented. The values of the Schmidt number Sc are chosen in such a way that they represent the diffusing chemical species of most common interest in air. For example, the values of Sc for H2, H20, NH3, propyl benzene and helium in air are 0.22, 0.60, 0.78, 2.62 and 0.30, respectively as reported by Perry (1963). It is seen that the concentration profiles decreases with an increase in the Schmidt number. It is noted that for heavier diffusing foreign speices, i.e., increasing the Schmidt number reduces the concentration in both magnitude and extent and thining of concentration boundary layer occurs. Figure 4: Effects of Sc on concentration distribution Skin Friction coefficient The distribution of skin friction coefficients are presented in Figs. 4a-4b for different values of the magnetic field, rotation parameter and porosity parameter. Inspection shows that increasing rotation () reduces skin friction. Increasing porosity (K) is found to increase the skin friction, in consistency with earlier discussion for the velocity response (Fig. 2c). Flow reversal is observed i.e. skin friction becomes negative in both the figures 4a and 4b. Clearly all profiles decay as M increases since larger Hartmann number corresponds to greater magnetohydrodynamic drag which decelerates the flow. Figure 4a: Effects of  on skin friction against M -6.25E-1 0.015 0.03 0.045 0 2 4 6 8 10 v/w0 z 5 10 20 15 Pr=0.71, Gr=5=Gm, =5, t=2, K=0.5, =10, =0.005 0 0.2 0.4 0.6 0.8 1 0 1 2 3 4 5  Z Sc=0.60 Sc=0.30 Sc=0.78 Sc=2.62 -6 -5 -4 -3 -2 -1 0 0 5 10 15 20  M 0 5 10Pr=0.71, Gr=5=Gm, K=0.5, t=2, =10, =0.005
  • 9. Analytical solution for Transient MHD… www.ijesi.org 41 | Page Figure 4b: Effects of K on skin friction against M Rate of Heat Transfer (Nu) Table 1: Heat transfer coefficient when  = 0.005, in the case air and water: Pr = 0.71 (Air) Pr = 7.0 (Water)  t M Nu Nu 2 2 2 0.34410 3.71514 2 2 3 0.37135 3.80426 3 2 3 0.80371 4.37293 3 3 3 0.03704 2.10853 3 3 4 0.02472 2.05830 4 3 4 0.90514 4.43715 4 4 4 0.07632 2.35304 Table 3 shows the magnitude of the heat transfer coefficient for different values of , t and M. An increase either in the Hartmann number (M) or time (t) the magnitude of the heat transfer coefficient reduces for both the cases of air and water. However, the heat transfer coefficient rises due to an increase in the frequency parameter (). Significantly, it is observed that the heat transfer coefficient in the case of water for any particular values of , t and M is much higher than that of air. V. Conclusion The present theoretical analysis brings out the following results of physical interest on the primary and secondary velocities, concentration profiles, skin friction and rate of heat transfer of the flow field for a rotating system in a saturated porous regime:  It is noticed that all the velocity profiles increase steadily near the plate and thereafter they show a constant decrease and reach the value zero at the free stream.  The magnetic / rotation parameter is found to decelerate the primary velocity of the flow field.  The porosity drag force has an effect accelerating the primary velocity profiles and reduces the skin friction.  There is a clear back flow for all the primary velocity profiles.  The secondary velocity profiles and skin friction coefficient are opposite to each other due to the effect of rotation.  The magnetic parameter has the effect of decreasing the heat transfer profiles whereas it has the effect of increasing the skin friction.  The presence of foreign species reduces the concentration boundary layer and further reduction occurs with increasing values of the Schmidt number.  The Prandtl number and the frequency parameter have the effect of increasing the heat transfer coefficient. Appendix: ))(Pr(Pr Gr A 32 1    , ))(Sc(Sc Gm A 32 2    , )A(AA 213  , -2.5 -2 -1.5 -1 -0.5 0 5 10 15 20  M K=0.1 K=0.3 K=0.5 K=1.0 Pr=0.71, Gr=5=Gm, =5, t=2, =10, =0.005
  • 10. Analytical solution for Transient MHD… www.ijesi.org 42 | Page 2 )Pr4i ωPrPr( 2 1   , 2 K 1 M2i411 2 2                 , 2 K 1 M2i411 2 3                 , 2 iω K 1 M2i411 2 4                 , 2 iω K 1 M2i411 2 5                 ,             zz 15 tiω r 15 ee Gr epartRealN   ,   .ee Gr epartImaginaryN zz 15 tiω i 15              References [1] Eckert, E.R., and Drake, R.M., 1972, “Analysis of Heat and Mass Transfer,” McGraw-Hill, pp. 171. [2] Elbashbeshy, E.M.A., 1997, “Heat and Mass Transfer along a vertical plate with variable surface tension and concentration in the presence of the magnetic field,” Int. J. Engineering Science, 34(5), pp. 515-522. [3] Helmy, K.A., 1998, “MHD Unsteady free convection flow past a vertical porous plate,” ZAMM, 78, pp. 255-270. [4] Soundalgekar, V.M., 1982, “Effects of mass transfer on flow past a uniformly accelerated vertical plate,” Letters in heat and mass transfer, 9, pp.65-72. [5] Kim, Y.J., 2000, “Unsteady MHD convective heat transfer past semi-infinite vertical porous moving plate with variable suction,” International Journal of Engineering Science, 38, pp. 833-845. [6] Ahmed, S., 2009, “The study of heat and mass transfer on free convective three-dimensional unsteady flows over a porous vertical plate,” Journal of Energy Heat and Mass Transfer, 31, pp. 89-110. [7] Ahmed, S., and Liu, I.C., 2010, Mixed convective three-dimensional heat and mass transfer flow with transversely periodic suction velocity,” Int. J. Applied Mathematics and Mechanics, 6(1), pp. 58-73. [8] Ingham, D.B., and Pop, I., 1998, Transport Phenomena in Porous Media, Pergamon, Oxford. [9] Nield, D.A. and Bejan, A., 1998, Convection in Porous Media, 2. Ed., Springer-Verlag, Berlin. [10] Cheng, P., 1977, “The influence of lateral mass flux on free convection boundary layers in a saturated porous medium,” Int.J. Heat Mass Transfer, 20, pp. 201- 206. [11] Kim, S., and Vafai K. 1989. Analysis of natural convection about a vertical plate embedded in a porous medium, Int. J. Heat Mass Transfer, vol. 32, pp. 665-677. [12] Raptis, A., and Singh, A.K., 1985, “Free convection flow past an impulsively started vertical plate in a porous medium by finite difference method,” Astrophys. Space Sci., 112, pp. 259-265. [13] Raptis, A.A., 1983, “Unsteady free convection through porous medium,” Int. J. Engg. sci., 21, pp. 345-348. [14] Raptis, A.A., and Perdikis, C.P., 1985, “Oscillatory flow through a porous medium by the presence of free convective flow,” Int. J. Engg. sci., 23, pp.51-55. [15] Ahmed, S., 2008, Transient three dimensional flow through a porous medium with transverse permeability oscillating with time,” Emirate Journal for Engineering research, 13(3), pp. 11-17. [16] Ahmed, S., 2010, “Free convective transient three-dimensional flow through a porous medium oscillating with time in presence of periodic suction velocity,” Int. J. Applied Mathematics and Mechanics, 6(11), pp. 1-16. [17] Gupta, A.S., 1972a, “Ekman layer on a porous plate,” Phys. Fluids, 15, pp. 930-931. [18] Jana, R.N., and Datta, N., 1977, “Couette flow and heat transfer in a rotating system,” Acta Mech., 26, pp. 301-306. [19] Singh, K.D., 2000, “An oscillatory hydromagnetic Couette flow in a rotating system,” J. Appl. Math and Mech., 80, pp. 429-432. [20] Singh, K.D., 1999, “Three dimensional Coutte flow with transpiration cooling,” Z. Angew. Math. Phys., 50, pp. 661-668. [21] Ganapathy, R., 1994, “A note on oscillatory Couette flow in a rotating system,” ASME J. Appl. Mech., 61, pp. 208-209. [22] Mazumder, B.S., 1991, “An exact solution of oscillatory Couette flow in a rotating system,” ASME J. Appl. Mech., 58, pp. 1104- 1107. [23] Mazumder, B.S., Gupta, A.S., and Dutta, N., 1976, “Flow and heat transfer in the hydromagnetic Ekman layer on a porous plate with Hall effects,” Int. J. Heat Mass Transfer, 19, pp. 523-527. [24] Singh, K.D., Gorla, M.G., and Hans, Raj., 2005, “A periodic solution of oscillatory Couette flow through porous medium in rotating system,” Indian J. pure appl. Math., 36, no. 3, pp. 151-159. [25] Singh, K.D., and Sharma, R., 2001, “Three dimensional Couette flow through a porous medium with heat transfer,” Indian J. pure appl. Math., 32(12), pp. 1819-1829. [26] Ahmed, S., and Zueco, J., 2011, “Modelling of heat and mass transfer in a rotating vertical porous channel with hall current,” Chemical Engineering Communications, 198(10), pp. 1294-1308. [27] Perry, E.D., 1963, Chemical Engineers Handbook, (4th Edn.), McGraw-Hill Book Company, New York.