Heat transfer in viscous free convective fluctuating mhd flow through porous media past a vertical porous plate with variable temperature

Mathematical Theory and Modeling www.iiste.org
ISSN 2224-5804 (Paper) ISSN 2225-0522 (Online)
Vol.2, No.6, 2012

Heat Transfer in Viscous Free Convective Fluctuating MHD Flow
Through Porous Media Past a Vertical Porous Plate with Variable
Temperature
S.R.Mishra, G.C.Dash and M.Acharya
Department of Mathematics,
Institute of Technical Education and Research,
Siksha ‘O’ Anusandhan University, Khandagiri,
Bhubaneswar-751030, Orissa, INDIA
E-mail: satyaranjan_mshr@yahoo.co.in
Abstract
In this paper, free convective Magnetohydrodynamics (MHD) flow of a viscous incompressible and
electrically conducting fluid past a hot vertical porous plate embedded in a porous medium has been studied. The
temperature of the plate varies both in space and time. The main objective of this paper is to study the effect of
porosity of the medium coupled with the variation of plate temperature with regards to space and time. The effect of
pertinent parameters characterizing the flow has been presented through the graph. The most interesting finding is
that presence of porous media has no significant contribution to the flow characteristics. Further, heating and cooling
of the plate due to convective current is compensated by the viscous dissipation.

Keywords: MHD flow/Span wise Co-sinusoidal/Free Convection /Heat transfer, porous medium.
1. Introduction
Many Industrial applications use Magnetohyrodynamics (MHD) effects to resolve the complex problems
very often occurred in industries. The available hydrodynamics solutions include the effects of magnetic field which
is possible as because the most of the industrial fluids are electrically conducting. For example, liquid metal MHD
takes its root in hydrodynamics of incompressible media which gains importance in the metallurgical industry,
Nuclear reactor, sodium cooling system ,every storage and electrical power generation[1974-76].Free convective
flows are of great interest in a number of industrial applications such as fiber and granular insulation, geothermal
system etc. Buoyancy is also of importance in an environment where difference between land and air temperature
can give rise to complicated flow patterns. The unsteady free convection flow past an infinite porous plate and
semi-infinite plate were studied by Nanda and Sharma [1962] .In their first paper they assumed the suction velocity
−1/2
at the plate varying in time as t ,where as in the second paper the plate temperature was assumed to oscillate in
time about a constant non-zero mean. Free convective flow past a vertical plate has been studied extensively by
Ostrach [1953] and many others. The free convective heat transfer on vertical semi-infinite plate was investigated by
Berezovsky [1977]. Martynenko et al. [1984] investigated the laminar free convection from a vertical plate.
The basic equations of incompressible MHD flow are non-linear. But there are many interesting cases where the
equations became linear in terms of the unknown quantities and may be solved easily. Linear MHD problems are
accessible to exact solutions and adopt the approximations that the density and transport properties be constant. No
fluid is incompressible but all may be treated as such whenever the pressure changes are small in comparison with
the bulk modulus. Mention may be made to the works in [1975, 72]. Ferdows et al. [2004] analysed free convection
flow with variable suction in presence of thermal radiation. Alam et al. [2006] studied Dufour and Soret effect with
variable suction on unsteady MHD free convection flow along a porous plate. Majumdar et al. [2007] gave an exact
solution for MHD flow past an impulsively started infinite vertical plate in the presence of thermal radiation.
Muthucumaraswamy et al. [2009] studied unsteady flow past an accelerated infinite vertical plate with variable
temperature and uniform mass diffusion. Recently, Dash et al[2009] have studied free convective MHD flow of a
visco-elastic fluid past an infinite vertical porous plate in a rotating frame of reference in the presence of
chemical reaction. In the present study we have set the flow through porous media with uniform porous matrix with
suction and blowing at the plate surface besides the free convective MHD effects and fluctuating surface temperature.

66

Vol.2, No.6, 2012

From the established result it is clear that the suction prevents the imposed non-torsional oscillations spreading away
from the oscillating surface (disk) by viscous diffusion for all values of frequency of oscillations. On the contrary the
blowing promotes the spreading of the oscillations far away from the disk and hence the boundary layer tends to be
infinitely thick when the disk is forced to oscillate with resonant frequency. In other words, in case of blowing and
resonance the oscillatory boundary layer flows are no longer possible.
Therefore, in the present study it aims at finding a meaningful solution for a non-linear coupled equation to bring
out the effects of suction/blowing with varying span- wise co-sinusoidal time dependent temperature in the presence of
uniform porous matrix in a free convective Magnetohydrodynamic flow past a vertical porous plate.

2. Formulation of the problem
An unsteady flow of a viscous incompressible electrically conducting fluid through a porous medium past an
insulated, infinite, hot, porous plate lying vertically on the x − z plane is considered. The x -axis is oriented in the
* * *
*
direction of the buoyancy force and y -axis is taken perpendicular to the plane of the plate. A uniform magnetic field
* * *
of strength B0 is applied along the y*-axis. Let (u , v , w ) be the component of velocity in the direction
( x , y , z ) respectively. The plate being considered infinite in x* direction, hence all the physical quantities are
* * *

independent of x . Thus, following Acharya and Padhi [1983], ω is independent of z and the equation of continuity
* * *

gives v = −V (constant) throughout.
*

We assume the span wise co-sinusoidal temperature of the form
T = T * + ε (T0* − T∞ ) cos(π z * / l − ω *t * ),
*

(2.1)
*
The mean temperature T of the plate is supplemented by the secondary temperature
ε (T0 − T∞ ) cos(π z / l − ω *t * ) varying with space and time. Under the usual Boussineques approximation the free
* * *

convective flow through porous media is governed by the following equations:

ut* + v*u * = υ (u * + u * ) + g β (T * − T∞ ) − σ B02 u * / ρ − υ u * / K P* ,
y yy zz
*
(2.2)

K µ
Tt * + v*Ty * = (Tyy * + Tzz * ) + ((u * )2 + (u z * ) 2 ) (2.3)
ρ CP ρ CP y

The boundary conditions are given by

y* = 0 : u * = 0, v* = V (Const.), T * = T0* + ε (T0* − T∞* ) cos(π z * / l − ω *t * ) (2.4)

y → ∞ : u → 0, T → T∞ .
* * * *

Introducing the non dimensional quantities defined in the nomenclature, we get,

ω 1 1 1
ut − u y = (u yy + u zz ) + Grθ − (M 2 + )u , (2.5)
Re Re Re KP
ω 1
θt − θ y = (θ yy + θ zz ) + Re Ec(u y 2 + u z 2 ), (2.6)
Re Pr Re
with corresponding boundary conditions:

67

Vol.2, No.6, 2012

y = 0 : u = 0,θ = 1 + ε cos(π z − t ),
(2.7)
y → ∞ : u → 0,θ → 0.
Since the amplitude, ε (<< 1), of the plate temperature is very small, we represent the velocity and temperature in the
neighborhood of the plate as

u ( y, z , t ) = u0 ( y ) + ε u1 ( y, z , t ) + o(ε 2 ),
(2.8)
θ ( y, z , t ) = θ 0 ( y ) + εθ1 ( y, z, t ) + o(ε 2 ),
Comparing the coefficient of like powers of ε after substituting (2.8) in (2.5) and (2.6), we get the following zeroth
order equations:

1
u0 yy + Re u0 y − ( M 2 + )u0 = − Re Grθ 0 , (2.9)
KP

θ0 yy + Pr Re θ 0 y = − Pr Re2 Ecu0 y 2 (2.10)

For solving the above coupled equations we use the following perturbation equations with perturbation parameter Ec,
the Eckert number,

u0 = u01 + Ecu02 + o( Ec 2 ) 

2 
(2.11)
θ0 = θ01 + Ecθ02 + o( Ec )  
Substituting (2.11) into (2.9) and (2.10) we get the following zeroth and first order equations of Ec .

1
′′ ′
u01 + Re u0 − ( M 2 + )u01 = − Re Grθ 01 , (2.12)
KP

θ 01 + Pr Re θ 01 = 0,
′′ ′ (2.13)

1
′′ ′
u02 + Re u02 − ( M 2 + )u02 = − Re Grθ 02 , (2.14)
KP

′′ ′ ′
θ 02 + Pr Re θ 02 = − Pr Re 2 u012 , (2.15)

The corresponding boundary conditions are:

y = 0; u01 = 0, θ 01 = 1, u02 = 0,θ 02 = 0,
(2.16)
y → ∞; u01 = 0,θ 01 = 0, u02 = 0,θ 02 = 0.

The solution of the equations (2.12) to (2.15) under the boundary conditions (2.16) are

u01 ( y ) = C3 (e− Pr Re y − e − m1 y ), (2.17)

68

Vol.2, No.6, 2012

θ 01 ( y ) = e − Pr Re y , (2.18)

u02 ( y ) = C4e− m1 y + A7e− Pr Re y + A8e−2Pr Re y + A9e −2 m1 y + A10e− (Pr Re + m1 ) y , (2.19)

θ 02 ( y ) = A3e − Pr Re y + A4e−2Pr Re y + A5e−2 m y + A6e − (Pr Re +m ) y ,
1 1
(2.20)

The terms of the coefficient of ε give the following first order equations:

1
ωu1t − Re u1 y = (u1 yy + u1zz ) + Gr Re θ1 − ( M 2 + )u1 , (2.21)
KP
1
ωθ1t − Re θ1 y = (θ1 yy + θ1zz ) + Re2 u0 y u1 y . (2.22)
Pr
In order to solve (2.21) and (2.22),it is convenient to adopt complex notations for velocity and temperature profile as,

u1 ( y, z , t ) = φ ( y )ei (π z −t ) . (2.23)

θ1 ( y, z, t ) = φ ( y )ei (π z −t ) .

The solutions obtained in terms of complex notations, the real part of which have physical significance.

Now, substituting (2.23) into (2.21) and (2.22) we get the following coupled equations:

1
φ ′′( y ) + Re φ ′( y ) + [ωi − π 2 − ( M 2 + )]φ = − Re Grψ ( y ), (2.24)
KP

ψ ′′( y ) + Pr Reψ ′( y ) + (ω Pr i − π 2 )ψ = −2 Re Pr uoyφ ′, (2.25)

Again, to uncouple above equations we assume the following perturbed forms:

φ = φ0 + Ecφ1 + o( Ec 2 ),
(2.26)
ψ = ψ 0 + Ecψ 1 + o( Ec 2 ),
Substituting (2.26) into (2.24) and (2.25) and equating the coefficient of like powers of Ec we get the subsequent
equations:

1
φ0′′ + Re φ0′ + [ωi − π 2 − ( M 2 + )]φ0 = − Re Grψ 0 , (2.27)
KP

′′ ′
ψ 0 + Pr Reψ 0 + (ω Pr i − π 2 )ψ 0 = 0, (2.28)

1
φ1′′+ Re φ1′ + [ωi − π 2 − ( M 2 + )]φ1 = − Re Grψ 1 , (2.29)
K

69

Vol.2, No.6, 2012

ψ 1′′ + Pr Reψ 1′ + (ω Pr i − π 2 )ψ 1 = −2 Re2 Pr u01φ0′ ,
′ / (2.30)

Where the prime denotes the differentiation with respect to y.
The corresponding boundary conditions are:
y = 0; φ0 = 0, φ1 = 0,ψ 0 = 1,ψ 1 = 0,
(2.31)
y → ∞;φ0 = 0, φ1 = 0,ψ 0 = 0,ψ 1 = 0,

The solutions of equations (2.27) to (2.30) under the boundary conditions (2.31) are:

φ0 ( y ) = A11 (e− m y − e− m y ),
2 3
(2.32)

ψ 0 ( y) = e−m y ,
2
(2.33)

ψ 1 ( y ) = A12e− (Pr Re +m ) y + A13e− (Pr Re + m ) y + A14e − ( m + m ) y + A15e − ( m + m ) y + A16e − m y ,
2 3 1 2 1 2 4
(2.34)

φ1 ( y ) = A17 e− (Pr Re + m ) y + A18e − (Pr Re + m ) y + A19e − ( m + m ) y + A20e − ( m + m ) y + A21e − m y + A22e − m y . (2.35)
2 3 1 2 1 3 4 3

The important flow characteristics of the problem are the plate shear stress and the rate of heat transfer at the plate.
The expressions for shear stress ( τ ) and Nusselt number ( Nu ) are given by

du
τ = τ *l / µV = Re al at y = 0
dy

= C3 (− Pr Re+ m1 ) − Ec[C4 m1 + A7 Pr Re+ 2 A8 Pr Re+ 2m1 A9 + A10 (Pr Re+ m1 )]
(2.36)
+ε F Cos (π z − t + α )

Where

F = Fr + iFi

= A11 (m3 − m2 ) − Ec[(m3 A22 + m4 A21 + A17 (Pr Re+ m2 ) + A18 (Pr Re+ m3 ) + A19 (m1 + m2 ) + A20 (m3 + m1 )]

The amplitude and the phase angle are given by

Fi
F = Fr2 + Fi 2 , α = tan −1
Fr

q*l dθ
Nu = − =Real at y = 0
k (T0 − T∞ )
* *
dy

= − Pr Re− Ec[Pr Re( A3 + 2 A4 ) + 2m1 A5 + A6 (Pr Re+ m1 + β )]
(2.37)
−ε G Cos (π z − t + β )

70

Vol.2, No.6, 2012

Where G = Gr + iGi

= m2 + Ec[(m14 A16 + A12 (Pr Re+ m2 ) + A13 (Pr Re+ m3 ) + A14 (m1 + m2 ) + A14 (m1 + m3 )]

The amplitude and the phase angle are given by

Gi
G = Gr2 + Gi2 , β = tan −1 .
Gr

3. Result and Discussions
This section analyses the velocity, temperature, and amplitude, phase angle of shear stress and rate of heat
transfer. The present discussion brings the following cases as particular case.
• M = 0, represents the case of non-conducting fluid without magnetic field.
• K P → ∞, represents without porous medium.
The most important part of the discussion is due to the presence of the sinusoidal variation of surface
temperature with space and time and the forcing forces such as Lorentz force, resistive force, due to porosity of the
medium, thermal buoyancy and cross flow due to permeable surface.
From equation (2.5) the following results follows.

 Vl 
In the absence of cross flow, V = 0 ⇒ Re  =  = 0, the u component of velocity remains unaffected by
 υ 

convective acceleration and thermal buoyancy force and the equation (2.5) reduces to

1  2 1  
ut = ( u yy + u zz ) −  M + u
ω  KP  
More over, the viscosity contributes significantly with a combined retarding effect caused by magnetic force and

resistance due to porous medium with an inverse multiplicity of the frequency of the temperature function. Further,

M = 0, K P → ∞, reduces the problem to a simple unsteady motion given by ut =
ω
1
(u yy + u zz ) .

Thus, in the absence of cross flow, that is, in case of a non-permeable surface, frequency of the fluctuating

 ω *l 2 
temperature parameter ω =  acts as a scaling factor.In case of impermeable surface, that is, in absence of
 υ 

cross flow, equation (6) reduces to θt =
1
(θ + θ ) .
ω Pr yy zz
This shows that the unsteady temperature gradient is reduced by high prandtle number fluid as well as
increasing frequency of the fluctuating temperature and viscous dissipation fails to effect the temperature
distribution.
From fig.1 it is observed that maximum velocity occurs in case of cooling of the plate ( Gr > 0 ) without
magnetic field (curve IX) and reverse effect is observed exclusively due to heating of the plate ( Gr < 0 ) with a

71

Vol.2, No.6, 2012

back flow (curve X). Thermal buoyancy effect has a significant contribution over the flow field.
From curves (VIII & IX), it is seen that the Lorentz force has a retarding effect which is in conformity with the
earlier reported result [14].
In the present study, velocity decrease is 7 times and there by steady state is reached within a few layers of the
flow domain. Further, it is seen that flow reversal occurs only in case of heating of the plate (curve X& XX) both in
porous and non porous medium.
Moreover, it is interesting to note that increase in Pr leads to significant decrease of the velocity (curve I (air)
and IV (water)) throughout the flow field but an increase in cross flow Reynolds number, increases the velocity in
the vicinity of the plate. Afterwards, it decreases rapidly. From curves, I & VIII it is observed that an increase in
viscous dissipation increases the velocity at all points. Thus, the energy loss sets a cooling current vis-à-vis
accelerating the velocity to reach the steady state.
These observations are clearly indicated by the additive and subtractive terms in the equation discussed earlier.
Fig.2 exhibits asymptotically decreasing behavior of the temperature distribution. For higher Pr fluid (Pr =7.0),
the temperature reduces drastically (curves I&IV).
It is interesting to remark that the cooling and heating of the plate (curves I & X) makes no difference in temperature
distribution which is compensated due to increase in dissipative energy loss.
Fig.3 shows the variation of amplitude of shear stress in both the cases i.e. presence and absence of porous
medium. Shear stress is almost linear for all values of frequency of fluctuation. For high Reynolds number as well as
greater buoyancy force shearing stress increases (curves II, XII, VII & XVII).
From fig.4 a steady increasing behavior is marked with an increasing frequency parameter but when ω
exceed 15.0 (approx.), slight decrease is marked. There is always a phase lead for all the parameters and for all
values of ω . A decrease in phase angle is marked due to increase in Reynolds number and magnetic parameter but
the reverse effect is observed in case of Prandtl number .Phase angle remains unaffected due to buoyancy effect,
dissipative loss and presence of porous medium.
Fig.5 and fig.6 exhibits the variation of amplitude G and phase angle tan β in case of heat transfer.
Further, G exhibits three layer characters due to high, medium and low value of Pr.

4. Conclusion
• Heating of the plate leads to back flow.
• Lorentz force has retarding effects.
• Presence of porous medium has no significant contribution.
• Viscous dissipation generates a cooling current which accelerates the velocity.
• There are three layer variation of amplitude of heat transfer with a phase lead.
• A phase lead is marked in case of shear stress and a phase lag for heat transfer which remains
unaffected by magnetic parameter, Grashoff number and Eckert number
Reference:
Stangeby,P.C(1974), A review of the status of MHD power generation technology including generation for a
Canadian MHD research programme, UTIAS rev, vol.39.
Lielausis,O.A,(1975), Liquid metal magnetohydrodynamics, Atomic Energy rev.,vol.13, pp.527
Ham, J.C.R and Moreu, R, (1976), Liquid metal magnetohydrodynamics with strong magnetic field”, a report on
Euromech, J.Fluid mech. vol.78, pp.261.

Nanda,R.S. and Sharma,V.P,(1962), Possibility similarity solutions of unsteady free convection flow past a vertical
plate with suction, J. phys. Soc. Japan ; vol.17(10),pp.1651-1656.

Ostrach,S.(1953), New aspects of natural convection heat transfer, Trans.Am.Soc.Mec.Engrs.; vol.75,pp.1287-1290.

Berezovsky, A. A., Martynenko, O. G. and Sokovishin,Yu. A, (1977) Free convective heat transfer on a vertical semi
infinite plate, J. Engng. Phys., vol.33, pp.32-39.

72

Vol.2, No.6, 2012

Martynenko,O. G., Berezovsky, A. A. and Sokovishin ,Yu. A,(1984), Laminar free convection from a vertical plate”,
Int. J. Heat Mass Transfer, vol.27, pp.869-881.

Murty ,S.N,(1975),MHD Ekman layer on a porous plate, Nuow cimento, vol. 29 ,pp.296.

Debnath ,L,(1972) “On unsteady MHD boundary layers in a rotating flow” ZAAM ,vol.52,pp.623.

Acharya,B.P. and Padhy,S,(1983), Free convective viscous flow along a hot vertical porous plate with periodic
temperature. Indian J. Pure appl.Math., vol.14(7),pp.838-849.

Ferdows, M., Sattar,M.A and Siddiki,M.N.A,(2004), Numerical approach on parameters of the thermal radiation
interaction with convection in boundary layer flow at a vertical plate with variable suction.
Thammasat International Journal of Science and Technology, vol.9, pp. 19-28.

Alam,M.S., Rahman ,M.M. and Samad ,M.A.,(2006), Dufour and Soret effects on unsteady MHD free convection
and mass transfer flow past a vertical plate in a porous medium. Nonlinear Analysis: Modeling and Control ,vol.11,
pp.271-226.

Majumder,M.K. and Deka,R.K,(2007), MHD flow past an impulsively started infinite vertical plate in the presence
of thermal radiation,Rom. Jour. Phys., vol.52 (5-7), pp.565-573.

Dash,G.C., Rath,P.K., Mohapatra,N. and Dash,P.K,(2009), Free convective MHD flow through porous media of a
rotating visco-elastic fluid past an infinite vertical porous plate with heat and mass transfer in the presence of a
chemical reaction, AMSE Modeling B,vol.78,No.4,pp.21-37.

Muthucumaraswamy R, Sundarraj M and Subhramanian V.S.A,(2009) Unsteady flow past an accelerated infinite

vertical plate with variable temperature and uniform mass diffusion. Int. J. of Appl. Math. and Mech. 5(6): 51-56.

Nomenclature:

( x* , y * , z * ) Cartesian Co-ordinate system V Suction velocity

l Wave length ε Amplitude of the span
wise co-sinusoidal temperature
υ Kinematics coefficient of viscosity µ Coefficient of viscosity
y = y* / l 

 Dimensionless space variable θ Dimensionless temperature
z = z /l 
*


(T * − T∞ ) (T0* − T∞ )
* *

β Coefficient of volumetric expansion B0 Uniform magnetic field strength

g Acceleration due to gravity β Coefficient of volume approximation

T0* Temperature of the plate *
T∞ ambient temperature

73

Vol.2, No.6, 2012

ρ Density CP Specific heat at constant pressure

K Thermal conductivity ω Frequency of temperature fluctuation-

ω *l 2 / υ

Re Reynolds number- Vl / υ Pr Prandtl number- µCP / K

Gr Grashoff number- υ g β (T0* − T∞ ) / V 3
*
Ec Eckert number- V 2 CP (T0* − T∞ )
*

M Hartmann number- B0l σ
µ τ Dimensionless skin friction

Nu Nusselt number K P Permeability of the medium

Appendix
1
Re+ Re 2 + 4( M 2 + )
KP Pr Re+ Pr 2 Re2 − 4(Pr iω − π 2 )
m1 = , m2 =
2 2

1
Re+ Re2 − 4(iω − π 2 − ( M 2 + ))
KP Pr Re+ Pr 2 Re2 − 4(Pr iω − π 2 )
m3 = , m4 =
2 2

− Re Gr 6
Pr Re2 C32
C3 = , A3 = −∑ Ai , A4 = − ,
1 2
Pr Re − Pr Re − ( M +
2 2 2 2
) i =4
KP

Pr Re 2 C32 m12 2 Pr 2 Re3 C32 m1 − Re GrA3
A5 = − , A6 = , A7 =
4m12 − 2m12 Pe Re (Pr Re+ m1 ) 2 − 2m12 Pr Re Pr 2 Re 2 − Pr Re2 − ( M 2 +
1
)
KP

− Re GrA4 − Re GrA5
A8 = , A9 =
1 1
4 Pr 2 Re2 − 2 Pr Re2 − ( M 2 + ) 4m12 Re2 − 2m1 Re− ( M 2 + )
KP KP

− Re GrA6 10
A10 = , C4 = −∑ Ai ,
1
(Pr Re+ m1 ) 2 − Re(Pr Re+ m1 ) − ( M 2 + ) i =7
KP

74

Vol.2, No.6, 2012

− Re Gr 2 Re3 Ec Pr 2 m2C3 A11
A11 = , A12 = −
m22 − Re m2 + [iω − π 2 − M 2 +
1
)] (Pr Re+ m2 )2 − Pr Re(Pr Re+ m2 ) + (iω Pr − π 2 )
KP

2
2 Re3 Ec Pr 2 m3C3 A11 2 Re Ec Pr m1m2 A11
A13 = , A14 =
(Pr Re+ m3 ) − Pr Re(Pr Re+ m3 ) + (iω Pr − π )
2 2
(m1 + m2 ) + Pr Re(m1 + m2 ) + (iω Pr − π 2 )
2

2 Re2 Ec Pr m1m3 A11 15
A15 = − , A16 = − ∑ Ai
(m1 + m3 ) + Pr Re(m1 + m3 ) + (iω Pr − π )
2 2
i =12

Re GrA13
A17 = −
1
(Pr Re+ m3 ) 2 − Re(Pr Re+ m3 ) + [iω − π 2 − ( M 2 + )]
KP

Re GrA12
A18 = −
1
(Pr Re+ m2 )2 − Re(Pr Re+ m2 ) + [iω − π 2 − ( M 2 + )]
KP

Re GrA14
A19 = −
1
(m1 + m2 ) 2 − Re(m1 + m2 ) + [iω − π 2 − ( M 2 + )]
KP

Re GrA15
A20 = −
1
(m1 + m3 ) 2 − Re(m1 + m3 ) + [iω − π 2 − ( M 2 + )]
KP

21
Re GrA16
A21 = − , A22 = −∑ Ai
1
m4 2 − Re m4 + [iω − π 2 − ( M 2 + )] i =17
KP

75

Vol.2, No.6, 2012

Re Pr M Gr Ec
3.5 I 2.0 0.71 2.0 5.0 0.01
IX
II 4.0 0.71 2.0 5.0 0.01
XIX III 2.0 2.00 2.0 5.0 0.01
3.0 VII IV 2.0 7.00 2.0 5.0 0.01
VIII,XVIII V 2.0 0.71 4.0 5.0 0.01
VI 2.0 7.00 4.0 5.0 0.01
IV,XIV VII 2.0 0.71 2.0 10.0 0.01
2.5 I,XI VIII 2.0 0.71 2.0 5.0 0.02
XVII IX 2.0 0.71 0.0 5.0 0.02
X 2.0 0.71 2.0 -5.0 0.02
V,XV
2.0 II,XII XI 2.0 0.71 2.0 5.0 0.01
XII 4.0 0.71 2.0 5.0 0.01
XIII 2.0 2.00 2.0 5.0 0.01
XIV 2.0 7.00 2.0 5.0 0.01
1.5 III,XIII XV 2.0 0.71 4.0 5.0 0.01
XVI 2.0 7.00 4.0 5.0 0.01
X,XX
u

XVII 2.0 0.71 2.0 10.0 0.01
XVIII 2.0 0.71 2.0 5.0 0.02
1.0 XIX 2.0 0.71 0.0 5.0 0.02
XX 2.0 0.71 2.0 -5.0 0.02

0.5

0.0
0.0 0.5 y 1.0 1.5 2.0
-0.5

POROUS MEDIUM
-1.0
WITOUT POROUS
MEDIUM

Fig.1 Velocity disrtibution.Kp=1(I-X),Kp=1000(XI-XX)

1.0 Re Pr M Gr Ec
I. 2.0 0.71 2.0 5.0 0.01
II. 4.0 0.71 2.0 5.0 0.01
IV,VI,XIV,XVI III. 2.0 2.00 2.0 5.0 0.01
0.8 IV. 2.0 7.00 2.0 5.0 0.01
III,XIII V. 2.0 0.71 4.0 5.0 0.01
VI. 2.0 7.00 4.0 5.0 0.01
II,XII VII. 2.0 0.71 2.0 10.0 0.01
0.6 VIII. 2.0 0.71 2.0 5.0 0.02
I,V,VIII,X,XI, IX. 2.0 0.71 0.0 5.0 0.02
XV,XVIII,XX X. 2.0 0.71 2.0 -5.0 0.02
θ

IX,XIX Re Pr M Gr Ec
0.4 XI. 2.0 0.71 2.0 5.0 0.01
VII,XVII XII. 4.0 0.71 2.0 5.0 0.01
XIII. 2.0 2.00 2.0 5.0 0.01
XIV. 2.0 7.00 2.0 5.0 0.01
XV. 2.0 0.71 4.0 5.0 0.01
0.2
XVI. 2.0 7.00 4.0 5.0 0.01
XVII. 2.0 0.71 2.0 10.0 0.01
XVIII. 2.0 0.71 2.0 5.0 0.02
XIX. 2.0 0.71 0.0 5.0 0.02
0.0
XX. 2.0 0.71 2.0 -5.0 0.02
0.0 0.5 1.0 1.5 2.0
y porous medium
Non-porous
Fig.2.Temperature distribution.Kp=1(I-X),Kp=1000(XI-XX)
medium

76

Vol.2, No.6, 2012

II,XII
3.0 VII,XVII Re Pr M Gr Ec
I. 2.0 0.71 2.0 5.0 0.01
II. 4.0 0.71 2.0 5.0 0.01
III. 2.0 2.00 2.0 5.0 0.01
IV. 2.0 7.00 2.0 5.0 0.01
2.5 V. 2.0 0.71 4.0 5.0 0.01
I,VIII,X,XI, VI. 2.0 7.00 4.0 5.0 0.01
XVIII,XX VII. 2.0 0.71 2.0 10.0 0.01
VIII. 2.0 0.71 2.0 5.0 0.02
2.0 V,XV IX. 2.0 0.71 0.0 5.0 0.02
VI,XVI X. 2.0 0.71 2.0 -5.0 0.02
IX,XIX
IV,XIV III,XIII
|F|

1.5
Re Pr M Gr Ec
XI. 2.0 0.71 2.0 5.0 0.01
XII. 4.0 0.71 2.0 5.0 0.01
XIII. 2.0 2.00 2.0 5.0 0.01
1.0 XIV. 2.0 7.00 2.0 5.0 0.01
XV. 2.0 0.71 4.0 5.0 0.01
XVI. 2.0 7.00 4.0 5.0 0.01
XVII. 2.0 0.71 2.0 10.0 0.01
XVIII. 2.0 0.71 2.0 5.0 0.02
0.5
XIX. 2.0 0.71 0.0 5.0 0.02
XX. 2.0 0.71 2.0 -5.0 0.02

0.0
5 10 15 20 25 30

ω
Fig. 3. The amplitude |F| of the shear stress.Kp=1(I-X),Kp=1000(XI-XX)

IX XIX
0.7 Re Pr M Gr Ec
I,VII,VIII,X,XVII, I. 2.0 0.71 2.0 5.0 0.01
XVIII,XI,XX III,XIII II. 4.0 0.71 2.0 5.0 0.01
III. 2.0 2.00 2.0 5.0 0.01
IV. 2.0 7.00 2.0 5.0 0.01
0.6
IV,XIV V. 2.0 0.71 4.0 5.0 0.01
VI. 2.0 7.00 4.0 5.0 0.01
VII. 2.0 0.71 2.0 10.0 0.01
VI,XVI II,XII VIII. 2.0 0.71 2.0 5.0 0.02
0.5 IX. 2.0 0.71 0.0 5.0 0.02
X. 2.0 0.71 2.0 -5.0 0.02

V,XV
Tan α

0.4 Re Pr M Gr Ec
XI. 2.0 0.71 2.0 5.0 0.01
XII. 4.0 0.71 2.0 5.0 0.01
XIII. 2.0 2.00 2.0 5.0 0.01
XIV. 2.0 7.00 2.0 5.0 0.01
0.3
XV. 2.0 0.71 4.0 5.0 0.01
XVI. 2.0 7.00 4.0 5.0 0.01
XVII. 2.0 0.71 2.0 10.0 0.01
XVIII. 2.0 0.71 2.0 5.0 0.02
0.2 XIX. 2.0 0.71 0.0 5.0 0.02
XX. 2.0 0.71 2.0 -5.0 0.02

0.1
5 10 15 20 25 30

ω
Fig. 4. The phase angle Tanα of shear stress,Kp=1(I-X),Kp=1000(XI-XX)

77

Vol.2, No.6, 2012

21.0 Re Pr M Gr Ec
I. 2.0 0.71 2.0 5.0 0.01
II. 4.0 0.71 2.0 5.0 0.01
19.0 III. 2.0 2.00 2.0 5.0 0.01
IV. 2.0 7.00 2.0 5.0 0.01
17.0 V. 2.0 0.71 4.0 5.0 0.01
VI. 2.0 7.00 4.0 5.0 0.01
VII. 2.0 0.71 2.0 10.0 0.01
15.0 I,V,VII,VIII,IX,X,XI,XV,XVII, VIII. 2.0 0.71 2.0 5.0 0.02
IX. 2.0 0.71 0.0 5.0 0.02
XVIII,XIX,XX X. 2.0 0.71 2.0 -5.0 0.02
13.0
|G|

III,XIII II,XII
Re Pr M Gr Ec
11.0 IV,VI,XIV,X XI. 2.0 0.71 2.0 5.0 0.01
VI XII. 4.0 0.71 2.0 5.0 0.01
9.0 XIII. 2.0 2.00 2.0 5.0 0.01
XIV. 2.0 7.00 2.0 5.0 0.01
XV. 2.0 0.71 4.0 5.0 0.01
7.0 XVI. 2.0 7.00 4.0 5.0 0.01
XVII. 2.0 0.71 2.0 10.0 0.01
XVIII. 2.0 0.71 2.0 5.0 0.02
5.0
XIX. 2.0 0.71 0.0 5.0 0.02
XX. 2.0 0.71 2.0 -5.0 0.02
3.0
5 10 15 20 25 30

ω
Fig.5. The amplitude |G| of the ate of heat transfer.
Kp=1(I-X),Kp=1000(XI-XX)

Re Pr M Gr Ec
-0.1
I. 2.0 0.71 2.0 5.0 0.01
5 10 15 20 25 30 II. 4.0 0.71 2.0 5.0 0.01
III. 2.0 2.00 2.0 5.0 0.01
IV. 2.0 7.00 2.0 5.0 0.01
-0.2 V. 2.0 0.71 4.0 5.0 0.01
VI. 2.0 7.00 4.0 5.0 0.01
VII. 2.0 0.71 2.0 10.0 0.01
VIII. 2.0 0.71 2.0 5.0 0.02
-0.3 IX. 2.0 0.71 0.0 5.0 0.02
X. 2.0 0.71 2.0 -5.0 0.02
Tan β

Re Pr M Gr Ec
-0.4 XI. 2.0 0.71 2.0 5.0 0.01
XII. 4.0 0.71 2.0 5.0 0.01
XIII. 2.0 2.00 2.0 5.0 0.01
XIV. 2.0 7.00 2.0 5.0 0.01
-0.5 III,XIII XV. 2.0 0.71 4.0 5.0 0.01
XVI. 2.0 7.00 4.0 5.0 0.01
XVII. 2.0 0.71 2.0 10.0 0.01
XVIII. 2.0 0.71 2.0 5.0 0.02
I,V,VII,VIII,IX,X,XI, XIX. 2.0 0.71 0.0 5.0 0.02
-0.6 XV,XVII,XVIII,XIX,XX IV,VI,XIV,XVI XX. 2.0 0.71 2.0 -5.0 0.02

II,XII

-0.7

ω
Fig.6.The phase angle Tanβ of the rate of heat transfer.
Kp=1(I-X),Kp=1000(XI-XX)

78

This academic article was published by The International Institute for Science,
Technology and Education (IISTE). The IISTE is a pioneer in the Open Access
Publishing service based in the U.S. and Europe. The aim of the institute is
Accelerating Global Knowledge Sharing.

More information about the publisher can be found in the IISTE’s homepage:
http://www.iiste.org

The IISTE is currently hosting more than 30 peer-reviewed academic journals and
collaborating with academic institutions around the world. Prospective authors of
IISTE journals can find the submission instruction on the following page:
http://www.iiste.org/Journals/

The IISTE editorial team promises to the review and publish all the qualified
submissions in a fast manner. All the journals articles are available online to the
readers all over the world without financial, legal, or technical barriers other than
those inseparable from gaining access to the internet itself. Printed version of the
journals is also available upon request of readers and authors.

IISTE Knowledge Sharing Partners

EBSCO, Index Copernicus, Ulrich's Periodicals Directory, JournalTOCS, PKP Open
Archives Harvester, Bielefeld Academic Search Engine, Elektronische
Zeitschriftenbibliothek EZB, Open J-Gate, OCLC WorldCat, Universe Digtial
Library , NewJour, Google Scholar

Heat transfer in viscous free convective fluctuating mhd flow through porous media past a vertical porous plate with variable temperature

More Related Content

What's hot (18)

Viewers also liked (8)

Similar to Heat transfer in viscous free convective fluctuating mhd flow through porous media past a vertical porous plate with variable temperature (20)

More from Alexander Decker (20)

Recently uploaded (20)

Heat transfer in viscous free convective fluctuating mhd flow through porous media past a vertical porous plate with variable temperature