Finite difference solutions of magneto hydrodynamic free
Transcript Of Finite difference solutions of magneto hydrodynamic free
IOSR Journal of Mathematics (IOSRJM) eISSN: 22785728, pISSN: 2319765X. Volume 11, Issue 1 Ver. II (Jan  Feb. 2015), PP 4258 www.iosrjournals.org
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction and variable thermal
conductivity in a DarcyForchheimer porous medium
1Ime Jimmy Uwanta, 2Halima Usman*
Department of Mathematics, Usmanu Danfodiyo University, Sokoto, Nigeria
Abstract: This paper presents the study of the effects of variable thermal conductivity and DarcyForchheimer on magnetohydrodynamic free convective flow in a vertical channel in the presence of constant suction. The resulting governing equations are nondimensionalised, simplified and solved using Crank Nicolson type of finite difference method. To check the accuracy of the numerical solution, steadystate solutions for velocity, temperature and concentration profiles are obtained by using perturbation method. Numerical results for the velocity, temperature and concentration profiles are illustrated graphically while the skin friction, Nusselt number and Sherwood number are tabulated and discussed for some selected controlling thermo physical parameters involved in the problem to show the behavior of the flow transport phenomena. It is found that the velocity and temperature increased due to increase in variable thermal conductivity parameter and there is decrease in concentration due to increase in DarcyForchheimer number. It is also observed that the numerical and analytical solutions are found to be in excellent agreement. Keywords: Variable thermal conductivity, Porous medium, Magnetohydrodynamic, Chemical reaction, DarcyForchheimer.
I. Introduction There has been a great deal of interest generated in the area of heat and mass transfer of free convective flow in porous medium due to it is numerous applications in a variety of industrial processes as well as in many natural circumstances. The subject of free convective flow in porous media has attracted considerable attention in the last few decades. Examples of such technological applications are ground water flows, thermal insulation engineering, food processing, soil pollution, fibrous insulation, geothermal extraction to name just a few. A comprehensive bibliography concerned with literature on convective flow in porous media can be found in the books of Nield and Bejan [12] and Ingham and Pop [6]. Processes involving heat and mass transfer effect have long been recognized as important principle in chemical processing equipment. The effect of chemical reaction, heat and mass transfer along a wedge with heat source and concentration in the presence of suction or injection have been examine by Kandasamy et al. [7]. Sattar [15] analyzed the effect of free and forced convection boundary layer flow through a porous medium with large suction. Mohammed et al. [11] investigated the effect of similarity solution for MHD flow through vertical porous plate with suction. In most of the studies on hydromagnetic heat transfer, thermal conductivity has been taken as constant. In metallurgical processing, the numerical value of thermal conductivity changes with temperature. Therefore the thermal conductivity is temperature dependent. In order to predict the accuracy of heat and mass transfer flow, mathematical models most consider the variation of thermal conductivity with temperature. Many papers were published with problem of thermal conductivity under different physical and geothermal conditions. Patowary and Dusmata [13] have presented the effects of thermal conductivity of micropolar fluid past a continuously moving plate with suction or injection in presence of magnetic field. Abdou [1] considered the effects of thermal radiation on unsteady boundary layer flow with temperature dependent viscosity and thermal conductivity due to a stretching sheet in porous medium. Similarly, Chiam [3] analyzed the effect of a variable thermal conductivity on the flow and heat transfer from a linearly stretching sheet. In another article, Seddek and Salema [18] have studied the effects of variable viscosity and thermal conductivity on an unsteady twodimensional laminar flow of a viscous incompressible conducting fluid past a vertical plate taking into account the effect of magnetic field in the presence of variable suction. Gitima [4] presented a numerical model on the effect of thermal radiation on unsteady boundary layer flow with variable viscosity and variable thermal conductivity due to stretching sheet in a porous channel in the presence of magnetic field. Additionally Hossain et al. [5] described the study of natural convection with variable viscosity and thermal conductivity from a vertical wavy cone. Sharma and Singh [19] have presented a theoretical study on the effects of variable thermal conductivity and heat source/sink on MHD flow near a stagnation point on a linearly stretching sheet. Very recently, Uwanta and Murtala [20] investigated the study of heat and mass transfer flow past an infinite vertical plate with variable thermal conductivity. Although considerable work has been reported on heat and mass transfer by free
DOI: 10.9790/572811124258
www.iosrjournals.org
42  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
convection in a porous medium, majority of porous studies have been on Darcy’s law which states that the volumeaveraged velocity is proportional to the pressure gradient which is valid only for slow flows through porous media with low permeability. However, Darcy’s law fails when the velocities and inertial effects are greater. Inertial porous phenomena can be accounted through the addition of a velocitysquared term in the momentum equation and the resulting model is known as DarcyForchheimer model Vafai and Tien [21]. The DarcyForchheimer model described the effect of inertia as well as viscous forces in porous media. Many investigators have reported work on Darcy and DarcyForchheimer. Numerical simulation of chemical reaction and viscous dissipation effects on DarcyForchheimer mixed convection in a fluid saturated media has been investigate by Mahdy and Chamkha [10]. Salem [17] has presented the radiation and mass transfer effects in DarcyForchheimer mixed convection from a vertical plate embedded in a fluid saturated porous medium. Lai [9] studied mixed convection heat and mass transfer from a vertical wall in a Darcian fluid saturated porous medium. Additionally, Kishan et al. [8] analyzed MHD free convective heat and mass transfer from a vertical surface embedded in a porous medium. A detailed numerical study has been carried out on the natural convection in a porous square cavity with an isoflux and isothermal discrete heater placed at the left wall by Saeid and Pop [16] using Darcy law. Rawat et al. [14] have investigated a theoretical study on free convection MHD micropolar flow, heat and species diffusion between vertical plates enclosing a nonDarcian porous medium with variable thermal conductivity and internal heat generation/absorption effects.
Inspite of all these studies, the aim of this work is to investigate the numerical solutions to the system of equations for magnetohydrodynamic free convective flow past a finite vertical channel with constant suction and thermal conductivity in a porous medium using finite difference method.
II. Mathematical analysis Consider an unsteady laminar onedimensional boundary layer heat and mass transfer flow in an electrically conducting fluid with thermal conductivity, chemical reaction and radiation in a DarcyForchheimer porous medium. Fluid suction or injection is imposed at the boundary of the cylinder. The x'  axis is taken
along the plate in the vertically upward direction and the y'  axis is taken normal to the plate. A uniform
transverse magnetic field B0 is applied normal to the flow direction. The magnetic Reynolds number is assumed
to be small so that the assumed magnetic field can be neglected. It is also assumed that the effects of joule
heating and viscous dissipation are neglected and the external electric field due to the effect of charge
polarization and the hall effect of MHD are neglected. Under these assumptions, the governing boundary layer
equations of continuity, momentum, energy and diffusion equations under Boussinesq’s approximations could
be written as follows:
v' 0 (1) y '
( ) ut ''  v0 uy'' = n 2yu'2'  æèççsrB02 + Kb'' öø÷÷u '2 Kn ' u' + gb(T '  T 0' ) + gb* C'  C0' (2)
T '  v T ' = k0
é ê1+a(T '  T '
)
T
'
ù ú

1
qr
(3)
t' 0 y' rCp y' ë
0 y' û rCp y'
C'  v C' = D 2C'  R*(C'  C' )
(4)
t' 0 y'
y'2
0
The corresponding initial and boundary conditions are prescribed as follows:
t 0, u' 0, T ' T 'w, C' C'w for all y'
t 0, u' 0, T ' T 'w , C' C 'w at y' 0 (5)
u' 0, T ' T0, C' C0 at y' = H
From continuity equation, it is clear that the suction velocity is either a constant or a function of time. Hence, on integrating equation (1), the suction velocity normal to the plate is assumed in the form,
' 0
where 0 is a scale of suction velocity which is nonzero positive constant. The negative sign indicates
that the suction is towards the plate and 0 0 corresponds to steady suction velocity normal at the surface.
DOI: 10.9790/572811124258
www.iosrjournals.org
43  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
The fourth and fifth terms on the right hand side of equation (2) denote the thermal and concentration buoyancy effects respectively. The last term of equation (3) represents the radiative heat flux term, while the last
term of equation (4) is the chemical reaction term. u' and v' are the Darcian velocity components in the x and y  directions respectively, t is the time, is the kinematic viscosity, g is the acceleration due to gravity,
is the coefficient of volume expansion, the density of the fluid, is the scalar electrical conductivity, *
is the volumetric coefficient of expansion with concentration, Cp is the specific heat capacity at constant
pressure, K ' is the Darcy permeability of the porous medium, b' is the empirical constant, k0 is the dimensionless thermal conductivity of the ambient fluid, is a constant depending on the nature of the fluid,
R* is the dimensionless chemical reaction, D is the coefficient of molecular diffusivity, v0 is the constant
suction parameter, qr is the radiative heat flux in the y direction, B0 is the magnetic induction of constant strength. T ' and T ' are the temperature of the fluid inside the thermal boundary layer and the fluid temperature
0
in the free stream respectively, while C ' and C0' are the corresponding concentrations.
To obtain the solutions of equations (2), (3) and (4) subject to the conditions (5) in nondimensional forms, we introduce the following nondimensional quantities:
u = u' ,t = t'u0 , y = y' ,q = T '  T 0' ,
u0 H 2
H
T 'w  T 0'
( ) ( ) C = C'  C'0 ,Pr = u0rCp , Sc = u0 , Da = K 'u0 ,
C 'w  C 0'
k0
D
H 2v
(6)
s B2H 2
H 2gb T 'w  T 0'
H 2gb * C 'w  C '0
M = 0 ,Gr =
,Gc =
,
ru0
u02
u02
( ) l = a T ' T '
, g = n0H ,K
=
R*H 2
,
R
=
16as
' 0
HT ' 0
3
, Fs =
b'u0
0
u ru
k 'u2
H 2v
0
0
0
Applying these nondimensionless quantities (6), the set of equations (2), (3), (4), and (5) reduces to the
following:
ut uy y2u2 M DFas u2 D1a u Gr GcC (7)
2 1
2 1
1 2 R
(8)
t y Pr y Pr
y Pr
C C 1 2C
2 KrC
(9)
t y Sc y
with the following initial and boundary conditions
t 0 u 0, 0, C 0, for all y,
t 0 u 0, 1, C 1, at y 0,
(10)
u 0, 0, C 0 at y 1.
Where Pr is the Prandtl number, Sc is the Schmidt number, M is the Magnetic field parameter, Gr is the thermal Grashof number, Gc is the Solutal or mass Grashof number, is the variable thermal
conductivity, is the variable suction parameter, R is the Radiation parameter, Kr is the chemical reaction
parameter, Da is the Darcy number, Fs is the DarcyForchheimer number, t is the dimensionless time while u and v are dimensionless velocity components in x  and y  directions respectively.
DOI: 10.9790/572811124258
www.iosrjournals.org
44  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
III. Analytical solutions In order to check the accuracy of the present numerical scheme of this model, there is need to compare numerical solutions with the analytical solutions. Since the governing equations are highly coupled and nonlinear, it is, therefore, of interest to reduce the governing equations of the present problem due to its nonlinearity into a form that can be solved analytically. At steady state the governing equations of the problem becomes:
2u u 1 u Gr GcC 0 (11) y2 y Da
2 Pr R 0 (12)
y2
y
2C
C
Sc KrScC 0
(13)
y2
y
The boundary conditions are
u 0, 1, C 1, at y 0,
u 0, 0, C 0 at y 1.
(14)
To find the approximate solution to equations (11)(13) subject to equation (14), we employed a regular perturbation method by taking a power series expansion and assume solution of the form:
u u0 y Ru1 y 0(R2 )
0 y R1 y 0(R2 ) (15) C C0 y RC1 y 0(R2 )
Using equation (15) into equations (11)(14) and equating the coefficient of like powers of R, the solutions to the governing equations are obtained as:
u F es1y F es2y F F e Pr y F eh1y F eh2y
9
10
11
12
13
14
(16)
R
F15e s1 y
F16es2 y
F17
F e 18
Pr y
F19 y F20
e Pr y
F1
F e Pr y 2
R F3
F e Pr y 4
F5
F6
ye Pr y
(17)
C F7eh1y F8eh2y
(18)
IV. Numerical solution procedure To solve the system of transformed coupled nonlinear partial differential equations (7)(9), under the initial and boundary conditions (10), an implicit finite difference scheme of CrankNicolson type which is unconditionally stable Carnahan et al, [2] has been employed. The equivalent finite difference approximations corresponding to equations (7)(9) are given as follows:
( ) æç ui, j+1  ui, j ö÷  g u  u =
ç è
Dt
÷ ø
2Dy
i+1, j
i1, j
( ) 1 u  2u + u + u  2u + u _
(19)
2 i+1, j+1
i, j+1 i1, j+1 i+1, j
i, j i1, j
( ) 2 Dy
( ) ( ) ( ) ( ) æçM + Fs ö÷ u
2

1
u
+ Gr q
+ Gc C
è
Da ø i, j Da i, j
i, j
i, j
DOI: 10.9790/572811124258
www.iosrjournals.org
45  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
( ) æçqi,
j+1
qi,
j
ö ÷

g
q
q
=
ç è
Dt
÷ ø
2Dy
i+1, j
i1, j
( ) H 2 qi+1, j+1  2qi, j+1 +qi1, j+1 +qi+1, j  2qi, j +qi1, j + (20) ( ) 2Pr Dy
( ) ( ) l
q
q
2

R
q
2 i+1, j i, j
( ) Pr Dy
Pr i, j
( ) æç Ci, j+1  Ci, j ö÷  g C  C =
ç è
Dt
÷ ø
2Dy
i+1, j
i1, j
( ) 1 C  2C +C + C  2C + C _
(21)
2 i+1, j+1
i, j+1
i1, j+1
i+1, j
i, j
i1, j
( ) 2Sc Dy
Kr (Ci, j )
The initial and boundary conditions take the following forms
ui, j = 0, qi, j = 0, Ci, j = 0
u0, j = 0, q0, j = 1, C0, j = 1
(22)
uH , j = 0, qH , j = 0, CH , j = 0 where H corresponds to 1.
The index i corresponds to space y and j corresponds to time t. y and t are the mesh sizes along y
direction and time tdirection respectively. The analytical solutions obtained in this work are used to check on the accuracy and effectiveness of the numerical scheme. Computations are carried out for different values of physical parameters involved in the problem.
The skin friction coefficient, Nusselt number and Sherwood number are important physical parameters for this type of boundary layer flow, which respectively are given by:
u
C
C f y y0 , Nu y y0 , Sh y y0 (23)
V. Results and discussion Numerical computations have been carried out for various values of thermo physical parameters
involved in this work such as thermal Grashof number ( Gr ), solutal Grashof number ( Gc ), Magnetic
parameter ( M ), Darcy number Da , suction parameter ( ), radiation parameter ( R ), chemical reaction
parameter ( Kr ), DarcyForchheimer number Fs , variable thermal conductivity ( ), Prandtl number ( Pr )
and Schmidt number ( Sc ) using an implicit method of CrankNicolson type. In order to illustrate the results graphically, the numerical values are plotted in Figs. (125). Values of local skin friction C f , Nusselt number
Nu and Sherwood number Sh are presented in Tables 1 and 2 for various values of physical parameters.
Solutions are obtained for fluids with Prandtl number (Pr = 0.71, 1.0 and 7.0) corresponding to air, salt water and water respectively. The values of Schmidt number are chosen in such a way as they represent the diffusing chemical species of most common interest in air and are taken for (Sc = 0.22 0.60, 0.78, 2.62) corresponding to
hydrogen, water, NH3 and Propyl benzene respectively. Default values of the physical parameters for velocity,
temperature and concentration profiles are specified as follows:
DOI: 10.9790/572811124258
www.iosrjournals.org
46  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
Gr = 5.0,Gc = 5.0, M = 2.0, Kr = 1.0,Pr = 0.71,Sc = 0.22, R = 3.0, Da = 0.01,
Fs = 3.0,g = 5.0,l = 0.05.
All graphs therefore correspond to these values unless specifically indicated on the appropriate graphs. The effects of various thermo physical parameters on the fluid velocity are illustrated in Figs.111 at
time t 0.1. The influence of thermal Grashof number and solutal Grashof number are presented in Figs. 1 and
2. As expected, it is observed that there is a rise in the velocity profiles due to the enhancement of thermal and mass buoyancy force. The velocity distribution increases rapidly near the porous plate and then decrease smoothly to the free stream velocity. Fig. 3 illustrates the effect of magnetic parameter on the velocity profile. It
is interesting to note that increase in M leads to decreasing the flow velocity. This is because the presence of
magnetic field in an electrically conducting fluid tends to produce a body force, which reduces the velocity of the fluid. Fig. 4 depicts typical velocity distribution across the boundary layer for various values of Darcy number. It is seen that the dimensionless velocity increases with increase in Darcy number. This is due to the
fact that larger values of Da parameter correspond to higher permeability porous media, which implies that,
less porous fiber resistance to the flow and therefore acceleration in transport. The effect of DarcyForchheimer
parameter is demonstrated in Fig. 5. It is noticed that an increase in Fs parameter increases the resistance to the
flow and so a decrease in the velocity profile. Figs. 68 graphically illustrate the behavior of Prandtl number, suction and Schmidt number respectively. It is observed that the momentum boundary layer decreases with the
increase in Pr this is due to the fact that convection current becomes weak and hence results in a decrease in the velocity. Also, it is seen that the velocity decreases with an increase in g parameter. This is because suction
parameter decelerates fluid particles through porous wall there by reducing the growth of the fluid boundary
layers. Figure 8 illustrates the effect of Sc on the velocity profiles. Clearly, it is observed that the velocity
distribution decreases with an increase in Sc. This is due to the fact that as Sc increases, the concentration buoyancy effect decreases yielding a reduction in the velocity profiles. The influence of variable thermal conductivity parameter, radiation parameter and dimensionless time are displayed in Figs. 911 respectively. It
is observed that as and t increases the momentum boundary layer increases, while an increase in radiation
parameter obviously decreases the velocity within the boundary layer. This is because larger values of
R correspond to an increased dominance of conduction over radiation, thereby decreasing the buoyancy force
and the thickness of the momentum boundary layer. The influences of various thermo physical parameters on the temperature distribution are respectively
presented in Figs. 1217 at time t 0.2 . Fig. 12 depicts the behavior of Prandtl number that is inversely
proportional to the thermal diffusivity of the working fluid on temperature distribution. It is noticed that the
thermal boundary layer decreases with increasing Pr . This is due to the fact that as values of Pr increases, the
thermal diffusivity of the fluid increases which results in a corresponding decrease in the temperature within the fluid. The temperature profile decreases with increasing suction parameter in Fig.13. This is because suction decelerates the fluid properties through porous wall thereby reducing the thermal boundary layer. The effect of variable thermal conductivity parameter on temperature profile is shown in Fig. 14. Here it is observed that the temperature profile decreases with increase in the variable thermal conductivity parameter. Fig. 15 described the
influence of radiation parameter on the temperature profile. An increase in R leads to increasing the
temperature boundary layer. This is obvious because the effect of radiation is to increase the temperature distribution in the thermal boundary layer. Figs.16 and 17 graphically illustrate the effects of DarcyForchheimer parameter and dimensionless time. It is observed that there is no significant change in the temperature profile as DarcyForchheimer parameter is increased, but there is a rise in the thermal boundary layer as the time increases.
Figs. 1821 illustrate the concentration profiles for various values of Schmidt number, chemical
reaction parameter, suction parameter and DarcyForchheimer parameter respectively at time t 0.4 . A
decrease in concentration profile with increasing these parameters is observed from these figures. This is due to the fact that an increase of Schmidt number means decrease in molecular diffusivity, while suction parameter stabilizes the boundary layer growth and therefore the concentration profile decreased. The effect of dimensionless time on the concentration profile is presented in Fig. 22. It is seen that the concentration profile increases with increase in the dimensionless time. The validity of the present model has been verified by comparing the numerical solutions and the analytical solutions through Fig. 2325 for velocity, temperature and concentration profiles respectively.These results are presented to illustrate the accuracy of the numerical solution. It is observed that the agreement between the results is excellent because the curves corresponding to analytical and numerical solutions coincide with one another. This has established confidence in the numerical results reported in this paper.
DOI: 10.9790/572811124258
www.iosrjournals.org
47  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
0.08 0.06
Gr = 5, 10, 15, 20
Vel(U)
0.04
0.02
0
0
0.2
0.4
0.6
0.8
1
y
Fig.1. Effect of Gr on the velocity profiles
Fig.2. Effect of Gc on the velocity profiles
Fig.3. Effect of M on the velocity profiles
DOI: 10.9790/572811124258
www.iosrjournals.org
48  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
Fig.4. Effect of Da on the velocity profiles
Fig.5. Effect of Fs on the velocity profiles
Fig.6. Effect of Pr on the velocity profiles
DOI: 10.9790/572811124258
www.iosrjournals.org
49  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
Fig.7. Effect of ( = S) on the velocity profiles Fig.8. Effect of Sc on the velocity profiles
Fig.9. Effect of ( = V) on the velocity profiles
DOI: 10.9790/572811124258
www.iosrjournals.org
50  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
Vel(U)
0.12
0.1
R = 0, 5, 10, 15
0.08
0.06
0.04
0.02
0
0
0.2
0.4
0.6
0.8
1
y
Fig.10. Effect of R on the velocity profiles
0.05
0.04
t = 0.1, 0.2, 0.3, 0.4
0.03
Vel(U)
0.02 0.01
0
0
0.2
0.4
0.6
0.8
1
y
Fig.11. Effect of t on the velocity profiles
1
0.8
Pr = 0.71, 1.0, 3.0, 7.0
0.6
0.4
0.2
Tem(T)
0
0
0.2
0.4
0.6
0.8
1
y
Fig.12. Effect of Pr on the temperature profiles
DOI: 10.9790/572811124258
www.iosrjournals.org
51  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction and variable thermal
conductivity in a DarcyForchheimer porous medium
1Ime Jimmy Uwanta, 2Halima Usman*
Department of Mathematics, Usmanu Danfodiyo University, Sokoto, Nigeria
Abstract: This paper presents the study of the effects of variable thermal conductivity and DarcyForchheimer on magnetohydrodynamic free convective flow in a vertical channel in the presence of constant suction. The resulting governing equations are nondimensionalised, simplified and solved using Crank Nicolson type of finite difference method. To check the accuracy of the numerical solution, steadystate solutions for velocity, temperature and concentration profiles are obtained by using perturbation method. Numerical results for the velocity, temperature and concentration profiles are illustrated graphically while the skin friction, Nusselt number and Sherwood number are tabulated and discussed for some selected controlling thermo physical parameters involved in the problem to show the behavior of the flow transport phenomena. It is found that the velocity and temperature increased due to increase in variable thermal conductivity parameter and there is decrease in concentration due to increase in DarcyForchheimer number. It is also observed that the numerical and analytical solutions are found to be in excellent agreement. Keywords: Variable thermal conductivity, Porous medium, Magnetohydrodynamic, Chemical reaction, DarcyForchheimer.
I. Introduction There has been a great deal of interest generated in the area of heat and mass transfer of free convective flow in porous medium due to it is numerous applications in a variety of industrial processes as well as in many natural circumstances. The subject of free convective flow in porous media has attracted considerable attention in the last few decades. Examples of such technological applications are ground water flows, thermal insulation engineering, food processing, soil pollution, fibrous insulation, geothermal extraction to name just a few. A comprehensive bibliography concerned with literature on convective flow in porous media can be found in the books of Nield and Bejan [12] and Ingham and Pop [6]. Processes involving heat and mass transfer effect have long been recognized as important principle in chemical processing equipment. The effect of chemical reaction, heat and mass transfer along a wedge with heat source and concentration in the presence of suction or injection have been examine by Kandasamy et al. [7]. Sattar [15] analyzed the effect of free and forced convection boundary layer flow through a porous medium with large suction. Mohammed et al. [11] investigated the effect of similarity solution for MHD flow through vertical porous plate with suction. In most of the studies on hydromagnetic heat transfer, thermal conductivity has been taken as constant. In metallurgical processing, the numerical value of thermal conductivity changes with temperature. Therefore the thermal conductivity is temperature dependent. In order to predict the accuracy of heat and mass transfer flow, mathematical models most consider the variation of thermal conductivity with temperature. Many papers were published with problem of thermal conductivity under different physical and geothermal conditions. Patowary and Dusmata [13] have presented the effects of thermal conductivity of micropolar fluid past a continuously moving plate with suction or injection in presence of magnetic field. Abdou [1] considered the effects of thermal radiation on unsteady boundary layer flow with temperature dependent viscosity and thermal conductivity due to a stretching sheet in porous medium. Similarly, Chiam [3] analyzed the effect of a variable thermal conductivity on the flow and heat transfer from a linearly stretching sheet. In another article, Seddek and Salema [18] have studied the effects of variable viscosity and thermal conductivity on an unsteady twodimensional laminar flow of a viscous incompressible conducting fluid past a vertical plate taking into account the effect of magnetic field in the presence of variable suction. Gitima [4] presented a numerical model on the effect of thermal radiation on unsteady boundary layer flow with variable viscosity and variable thermal conductivity due to stretching sheet in a porous channel in the presence of magnetic field. Additionally Hossain et al. [5] described the study of natural convection with variable viscosity and thermal conductivity from a vertical wavy cone. Sharma and Singh [19] have presented a theoretical study on the effects of variable thermal conductivity and heat source/sink on MHD flow near a stagnation point on a linearly stretching sheet. Very recently, Uwanta and Murtala [20] investigated the study of heat and mass transfer flow past an infinite vertical plate with variable thermal conductivity. Although considerable work has been reported on heat and mass transfer by free
DOI: 10.9790/572811124258
www.iosrjournals.org
42  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
convection in a porous medium, majority of porous studies have been on Darcy’s law which states that the volumeaveraged velocity is proportional to the pressure gradient which is valid only for slow flows through porous media with low permeability. However, Darcy’s law fails when the velocities and inertial effects are greater. Inertial porous phenomena can be accounted through the addition of a velocitysquared term in the momentum equation and the resulting model is known as DarcyForchheimer model Vafai and Tien [21]. The DarcyForchheimer model described the effect of inertia as well as viscous forces in porous media. Many investigators have reported work on Darcy and DarcyForchheimer. Numerical simulation of chemical reaction and viscous dissipation effects on DarcyForchheimer mixed convection in a fluid saturated media has been investigate by Mahdy and Chamkha [10]. Salem [17] has presented the radiation and mass transfer effects in DarcyForchheimer mixed convection from a vertical plate embedded in a fluid saturated porous medium. Lai [9] studied mixed convection heat and mass transfer from a vertical wall in a Darcian fluid saturated porous medium. Additionally, Kishan et al. [8] analyzed MHD free convective heat and mass transfer from a vertical surface embedded in a porous medium. A detailed numerical study has been carried out on the natural convection in a porous square cavity with an isoflux and isothermal discrete heater placed at the left wall by Saeid and Pop [16] using Darcy law. Rawat et al. [14] have investigated a theoretical study on free convection MHD micropolar flow, heat and species diffusion between vertical plates enclosing a nonDarcian porous medium with variable thermal conductivity and internal heat generation/absorption effects.
Inspite of all these studies, the aim of this work is to investigate the numerical solutions to the system of equations for magnetohydrodynamic free convective flow past a finite vertical channel with constant suction and thermal conductivity in a porous medium using finite difference method.
II. Mathematical analysis Consider an unsteady laminar onedimensional boundary layer heat and mass transfer flow in an electrically conducting fluid with thermal conductivity, chemical reaction and radiation in a DarcyForchheimer porous medium. Fluid suction or injection is imposed at the boundary of the cylinder. The x'  axis is taken
along the plate in the vertically upward direction and the y'  axis is taken normal to the plate. A uniform
transverse magnetic field B0 is applied normal to the flow direction. The magnetic Reynolds number is assumed
to be small so that the assumed magnetic field can be neglected. It is also assumed that the effects of joule
heating and viscous dissipation are neglected and the external electric field due to the effect of charge
polarization and the hall effect of MHD are neglected. Under these assumptions, the governing boundary layer
equations of continuity, momentum, energy and diffusion equations under Boussinesq’s approximations could
be written as follows:
v' 0 (1) y '
( ) ut ''  v0 uy'' = n 2yu'2'  æèççsrB02 + Kb'' öø÷÷u '2 Kn ' u' + gb(T '  T 0' ) + gb* C'  C0' (2)
T '  v T ' = k0
é ê1+a(T '  T '
)
T
'
ù ú

1
qr
(3)
t' 0 y' rCp y' ë
0 y' û rCp y'
C'  v C' = D 2C'  R*(C'  C' )
(4)
t' 0 y'
y'2
0
The corresponding initial and boundary conditions are prescribed as follows:
t 0, u' 0, T ' T 'w, C' C'w for all y'
t 0, u' 0, T ' T 'w , C' C 'w at y' 0 (5)
u' 0, T ' T0, C' C0 at y' = H
From continuity equation, it is clear that the suction velocity is either a constant or a function of time. Hence, on integrating equation (1), the suction velocity normal to the plate is assumed in the form,
' 0
where 0 is a scale of suction velocity which is nonzero positive constant. The negative sign indicates
that the suction is towards the plate and 0 0 corresponds to steady suction velocity normal at the surface.
DOI: 10.9790/572811124258
www.iosrjournals.org
43  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
The fourth and fifth terms on the right hand side of equation (2) denote the thermal and concentration buoyancy effects respectively. The last term of equation (3) represents the radiative heat flux term, while the last
term of equation (4) is the chemical reaction term. u' and v' are the Darcian velocity components in the x and y  directions respectively, t is the time, is the kinematic viscosity, g is the acceleration due to gravity,
is the coefficient of volume expansion, the density of the fluid, is the scalar electrical conductivity, *
is the volumetric coefficient of expansion with concentration, Cp is the specific heat capacity at constant
pressure, K ' is the Darcy permeability of the porous medium, b' is the empirical constant, k0 is the dimensionless thermal conductivity of the ambient fluid, is a constant depending on the nature of the fluid,
R* is the dimensionless chemical reaction, D is the coefficient of molecular diffusivity, v0 is the constant
suction parameter, qr is the radiative heat flux in the y direction, B0 is the magnetic induction of constant strength. T ' and T ' are the temperature of the fluid inside the thermal boundary layer and the fluid temperature
0
in the free stream respectively, while C ' and C0' are the corresponding concentrations.
To obtain the solutions of equations (2), (3) and (4) subject to the conditions (5) in nondimensional forms, we introduce the following nondimensional quantities:
u = u' ,t = t'u0 , y = y' ,q = T '  T 0' ,
u0 H 2
H
T 'w  T 0'
( ) ( ) C = C'  C'0 ,Pr = u0rCp , Sc = u0 , Da = K 'u0 ,
C 'w  C 0'
k0
D
H 2v
(6)
s B2H 2
H 2gb T 'w  T 0'
H 2gb * C 'w  C '0
M = 0 ,Gr =
,Gc =
,
ru0
u02
u02
( ) l = a T ' T '
, g = n0H ,K
=
R*H 2
,
R
=
16as
' 0
HT ' 0
3
, Fs =
b'u0
0
u ru
k 'u2
H 2v
0
0
0
Applying these nondimensionless quantities (6), the set of equations (2), (3), (4), and (5) reduces to the
following:
ut uy y2u2 M DFas u2 D1a u Gr GcC (7)
2 1
2 1
1 2 R
(8)
t y Pr y Pr
y Pr
C C 1 2C
2 KrC
(9)
t y Sc y
with the following initial and boundary conditions
t 0 u 0, 0, C 0, for all y,
t 0 u 0, 1, C 1, at y 0,
(10)
u 0, 0, C 0 at y 1.
Where Pr is the Prandtl number, Sc is the Schmidt number, M is the Magnetic field parameter, Gr is the thermal Grashof number, Gc is the Solutal or mass Grashof number, is the variable thermal
conductivity, is the variable suction parameter, R is the Radiation parameter, Kr is the chemical reaction
parameter, Da is the Darcy number, Fs is the DarcyForchheimer number, t is the dimensionless time while u and v are dimensionless velocity components in x  and y  directions respectively.
DOI: 10.9790/572811124258
www.iosrjournals.org
44  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
III. Analytical solutions In order to check the accuracy of the present numerical scheme of this model, there is need to compare numerical solutions with the analytical solutions. Since the governing equations are highly coupled and nonlinear, it is, therefore, of interest to reduce the governing equations of the present problem due to its nonlinearity into a form that can be solved analytically. At steady state the governing equations of the problem becomes:
2u u 1 u Gr GcC 0 (11) y2 y Da
2 Pr R 0 (12)
y2
y
2C
C
Sc KrScC 0
(13)
y2
y
The boundary conditions are
u 0, 1, C 1, at y 0,
u 0, 0, C 0 at y 1.
(14)
To find the approximate solution to equations (11)(13) subject to equation (14), we employed a regular perturbation method by taking a power series expansion and assume solution of the form:
u u0 y Ru1 y 0(R2 )
0 y R1 y 0(R2 ) (15) C C0 y RC1 y 0(R2 )
Using equation (15) into equations (11)(14) and equating the coefficient of like powers of R, the solutions to the governing equations are obtained as:
u F es1y F es2y F F e Pr y F eh1y F eh2y
9
10
11
12
13
14
(16)
R
F15e s1 y
F16es2 y
F17
F e 18
Pr y
F19 y F20
e Pr y
F1
F e Pr y 2
R F3
F e Pr y 4
F5
F6
ye Pr y
(17)
C F7eh1y F8eh2y
(18)
IV. Numerical solution procedure To solve the system of transformed coupled nonlinear partial differential equations (7)(9), under the initial and boundary conditions (10), an implicit finite difference scheme of CrankNicolson type which is unconditionally stable Carnahan et al, [2] has been employed. The equivalent finite difference approximations corresponding to equations (7)(9) are given as follows:
( ) æç ui, j+1  ui, j ö÷  g u  u =
ç è
Dt
÷ ø
2Dy
i+1, j
i1, j
( ) 1 u  2u + u + u  2u + u _
(19)
2 i+1, j+1
i, j+1 i1, j+1 i+1, j
i, j i1, j
( ) 2 Dy
( ) ( ) ( ) ( ) æçM + Fs ö÷ u
2

1
u
+ Gr q
+ Gc C
è
Da ø i, j Da i, j
i, j
i, j
DOI: 10.9790/572811124258
www.iosrjournals.org
45  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
( ) æçqi,
j+1
qi,
j
ö ÷

g
q
q
=
ç è
Dt
÷ ø
2Dy
i+1, j
i1, j
( ) H 2 qi+1, j+1  2qi, j+1 +qi1, j+1 +qi+1, j  2qi, j +qi1, j + (20) ( ) 2Pr Dy
( ) ( ) l
q
q
2

R
q
2 i+1, j i, j
( ) Pr Dy
Pr i, j
( ) æç Ci, j+1  Ci, j ö÷  g C  C =
ç è
Dt
÷ ø
2Dy
i+1, j
i1, j
( ) 1 C  2C +C + C  2C + C _
(21)
2 i+1, j+1
i, j+1
i1, j+1
i+1, j
i, j
i1, j
( ) 2Sc Dy
Kr (Ci, j )
The initial and boundary conditions take the following forms
ui, j = 0, qi, j = 0, Ci, j = 0
u0, j = 0, q0, j = 1, C0, j = 1
(22)
uH , j = 0, qH , j = 0, CH , j = 0 where H corresponds to 1.
The index i corresponds to space y and j corresponds to time t. y and t are the mesh sizes along y
direction and time tdirection respectively. The analytical solutions obtained in this work are used to check on the accuracy and effectiveness of the numerical scheme. Computations are carried out for different values of physical parameters involved in the problem.
The skin friction coefficient, Nusselt number and Sherwood number are important physical parameters for this type of boundary layer flow, which respectively are given by:
u
C
C f y y0 , Nu y y0 , Sh y y0 (23)
V. Results and discussion Numerical computations have been carried out for various values of thermo physical parameters
involved in this work such as thermal Grashof number ( Gr ), solutal Grashof number ( Gc ), Magnetic
parameter ( M ), Darcy number Da , suction parameter ( ), radiation parameter ( R ), chemical reaction
parameter ( Kr ), DarcyForchheimer number Fs , variable thermal conductivity ( ), Prandtl number ( Pr )
and Schmidt number ( Sc ) using an implicit method of CrankNicolson type. In order to illustrate the results graphically, the numerical values are plotted in Figs. (125). Values of local skin friction C f , Nusselt number
Nu and Sherwood number Sh are presented in Tables 1 and 2 for various values of physical parameters.
Solutions are obtained for fluids with Prandtl number (Pr = 0.71, 1.0 and 7.0) corresponding to air, salt water and water respectively. The values of Schmidt number are chosen in such a way as they represent the diffusing chemical species of most common interest in air and are taken for (Sc = 0.22 0.60, 0.78, 2.62) corresponding to
hydrogen, water, NH3 and Propyl benzene respectively. Default values of the physical parameters for velocity,
temperature and concentration profiles are specified as follows:
DOI: 10.9790/572811124258
www.iosrjournals.org
46  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
Gr = 5.0,Gc = 5.0, M = 2.0, Kr = 1.0,Pr = 0.71,Sc = 0.22, R = 3.0, Da = 0.01,
Fs = 3.0,g = 5.0,l = 0.05.
All graphs therefore correspond to these values unless specifically indicated on the appropriate graphs. The effects of various thermo physical parameters on the fluid velocity are illustrated in Figs.111 at
time t 0.1. The influence of thermal Grashof number and solutal Grashof number are presented in Figs. 1 and
2. As expected, it is observed that there is a rise in the velocity profiles due to the enhancement of thermal and mass buoyancy force. The velocity distribution increases rapidly near the porous plate and then decrease smoothly to the free stream velocity. Fig. 3 illustrates the effect of magnetic parameter on the velocity profile. It
is interesting to note that increase in M leads to decreasing the flow velocity. This is because the presence of
magnetic field in an electrically conducting fluid tends to produce a body force, which reduces the velocity of the fluid. Fig. 4 depicts typical velocity distribution across the boundary layer for various values of Darcy number. It is seen that the dimensionless velocity increases with increase in Darcy number. This is due to the
fact that larger values of Da parameter correspond to higher permeability porous media, which implies that,
less porous fiber resistance to the flow and therefore acceleration in transport. The effect of DarcyForchheimer
parameter is demonstrated in Fig. 5. It is noticed that an increase in Fs parameter increases the resistance to the
flow and so a decrease in the velocity profile. Figs. 68 graphically illustrate the behavior of Prandtl number, suction and Schmidt number respectively. It is observed that the momentum boundary layer decreases with the
increase in Pr this is due to the fact that convection current becomes weak and hence results in a decrease in the velocity. Also, it is seen that the velocity decreases with an increase in g parameter. This is because suction
parameter decelerates fluid particles through porous wall there by reducing the growth of the fluid boundary
layers. Figure 8 illustrates the effect of Sc on the velocity profiles. Clearly, it is observed that the velocity
distribution decreases with an increase in Sc. This is due to the fact that as Sc increases, the concentration buoyancy effect decreases yielding a reduction in the velocity profiles. The influence of variable thermal conductivity parameter, radiation parameter and dimensionless time are displayed in Figs. 911 respectively. It
is observed that as and t increases the momentum boundary layer increases, while an increase in radiation
parameter obviously decreases the velocity within the boundary layer. This is because larger values of
R correspond to an increased dominance of conduction over radiation, thereby decreasing the buoyancy force
and the thickness of the momentum boundary layer. The influences of various thermo physical parameters on the temperature distribution are respectively
presented in Figs. 1217 at time t 0.2 . Fig. 12 depicts the behavior of Prandtl number that is inversely
proportional to the thermal diffusivity of the working fluid on temperature distribution. It is noticed that the
thermal boundary layer decreases with increasing Pr . This is due to the fact that as values of Pr increases, the
thermal diffusivity of the fluid increases which results in a corresponding decrease in the temperature within the fluid. The temperature profile decreases with increasing suction parameter in Fig.13. This is because suction decelerates the fluid properties through porous wall thereby reducing the thermal boundary layer. The effect of variable thermal conductivity parameter on temperature profile is shown in Fig. 14. Here it is observed that the temperature profile decreases with increase in the variable thermal conductivity parameter. Fig. 15 described the
influence of radiation parameter on the temperature profile. An increase in R leads to increasing the
temperature boundary layer. This is obvious because the effect of radiation is to increase the temperature distribution in the thermal boundary layer. Figs.16 and 17 graphically illustrate the effects of DarcyForchheimer parameter and dimensionless time. It is observed that there is no significant change in the temperature profile as DarcyForchheimer parameter is increased, but there is a rise in the thermal boundary layer as the time increases.
Figs. 1821 illustrate the concentration profiles for various values of Schmidt number, chemical
reaction parameter, suction parameter and DarcyForchheimer parameter respectively at time t 0.4 . A
decrease in concentration profile with increasing these parameters is observed from these figures. This is due to the fact that an increase of Schmidt number means decrease in molecular diffusivity, while suction parameter stabilizes the boundary layer growth and therefore the concentration profile decreased. The effect of dimensionless time on the concentration profile is presented in Fig. 22. It is seen that the concentration profile increases with increase in the dimensionless time. The validity of the present model has been verified by comparing the numerical solutions and the analytical solutions through Fig. 2325 for velocity, temperature and concentration profiles respectively.These results are presented to illustrate the accuracy of the numerical solution. It is observed that the agreement between the results is excellent because the curves corresponding to analytical and numerical solutions coincide with one another. This has established confidence in the numerical results reported in this paper.
DOI: 10.9790/572811124258
www.iosrjournals.org
47  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
0.08 0.06
Gr = 5, 10, 15, 20
Vel(U)
0.04
0.02
0
0
0.2
0.4
0.6
0.8
1
y
Fig.1. Effect of Gr on the velocity profiles
Fig.2. Effect of Gc on the velocity profiles
Fig.3. Effect of M on the velocity profiles
DOI: 10.9790/572811124258
www.iosrjournals.org
48  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
Fig.4. Effect of Da on the velocity profiles
Fig.5. Effect of Fs on the velocity profiles
Fig.6. Effect of Pr on the velocity profiles
DOI: 10.9790/572811124258
www.iosrjournals.org
49  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
Fig.7. Effect of ( = S) on the velocity profiles Fig.8. Effect of Sc on the velocity profiles
Fig.9. Effect of ( = V) on the velocity profiles
DOI: 10.9790/572811124258
www.iosrjournals.org
50  Page
Finite difference solutions of magneto hydrodynamic free convective flow with constant suction ….
Vel(U)
0.12
0.1
R = 0, 5, 10, 15
0.08
0.06
0.04
0.02
0
0
0.2
0.4
0.6
0.8
1
y
Fig.10. Effect of R on the velocity profiles
0.05
0.04
t = 0.1, 0.2, 0.3, 0.4
0.03
Vel(U)
0.02 0.01
0
0
0.2
0.4
0.6
0.8
1
y
Fig.11. Effect of t on the velocity profiles
1
0.8
Pr = 0.71, 1.0, 3.0, 7.0
0.6
0.4
0.2
Tem(T)
0
0
0.2
0.4
0.6
0.8
1
y
Fig.12. Effect of Pr on the temperature profiles
DOI: 10.9790/572811124258
www.iosrjournals.org
51  Page