# MHD Natural Convection with Convective Surface Boundary Condition over a Flat Plate

Authors: Mohammad M. Rashidi; Mohammad Ferdows; Amir Basiri Parsa; Shirley Abelman
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

We apply the one parameter continuous group method to investigate similarity solutions of magnetohydrodynamic (MHD) heat and mass transfer flow of a steady viscous incompressible fluid over a flat plate. By using the one parameter group method, similarity transformations and corresponding similarity representations are presented. A convective boundary condition is applied instead of the usual boundary conditions of constant surface temperature or constant heat flux. In addition it is assumed that viscosity, thermal conductivity, and concentration diffusivity vary linearly. Our study indicates that a similarity solution is possible if the convective heat transfer related to the hot fluid on the lower surface of the plate is directly proportional to(
x
-
)
-
1
/
2 where x
- is the distance from the leading edge of the solid surface. Numerical solutions of the ordinary differential equations are obtained by the Keller Box method for different values of the controlling parameters associated with the problem.


---

## Body

## 1. Introduction

A review of the literature shows that to the best of our knowledge not much research has been reported on MHD flow over a flat plate with convective surface boundary conditions by applying the one parameter continuous group method. For this problem we apply similarity transformations on the partial differential equations. The transformed nonlinear coupled ordinary differential equations are solved numerically by the Keller Box method for different values of controlling parameters.Analysis of natural phenomena usually leads to partial differential equations and nonlinear ordinary differential equations. Nonlinear differential equations appear in physics, applied mathematics, and engineering sciences. In most cases for these problems exact solutions cannot be obtained. One of the most widely used applications of nonlinear differential equations is boundary-layer problems. Fluid flow and heat transfer are a relevant problem in many industrial processes such as metal and polymer extrusion processes, glass-fiber and paper production, manufacture and drawing of plastics and rubber sheets, and crystal growing. Magnetohydrodynamics (MHD) is the flow of an electrically conducting fluid in the presence of a magnetic field. This effect is of importance in various areas of technology and engineering such as MHD flow meters, MHD power generation, and MHD pumps [1–4]. The study of the interaction of conducting fluids with electromagnetic phenomena is important and such problems have received much attention from many researchers. Mukhopadhyay et al. [5], Andersson [6], Rashidi et al. [7], and Parsa et al. [8] investigated the effect of magnetic field over a stretching surface in various states. Numerical results for MHD free convection flow over a wedge in the presence of a magnetic field were presented by Watanabe and Pop [9]. Kumari and Nath [10] studied unsteady MHD viscous flow and heat transfer of Newtonian fluids induced by an impulsively stretched plane surface in two lateral directions by employing the homotopy analysis method. Rashidi et al. [1] solved the governing equations of suction and injection effects on the free convection boundary-layer flow over a vertical cylinder. In addition, a complete investigation of MHD studies and their technological applications was undertaken by Moreau [11]. Several interesting computational studies of reactive MHD boundary-layer flows with heat and mass transfer have appeared in recent years [12–15]. Effects of anisotropic scattering on steady nonsimilar free convective radiative hydromagnetic boundary-layer flow over a diffuse reflecting surface and solution of a separate equation for the magnetic field distribution were presented by Chen [16]. Ishak [17] studied steady laminar boundary-layer flow and heat transfer over a stationary permeable flat plate immersed in a uniform free stream with convective boundary condition. The problem of a vertical plate with convective boundary conditions was considered by Makinde [18]. Rashidi et al. [19] presented the first and second law analyses of an electrically conducting fluid past a rotating disk in the presence of a uniform vertical magnetic field by using the homotopy analysis method (HAM) and then applied artificial neural networks (ANN) and the particle swarm optimization (PSO) algorithm to minimize the entropy generation.The main objective in this paper is to investigate similarity solutions and scaling transformations of MHD heat and mass transfer flow of a steady viscous incompressible fluid over a flat plate with convective surface boundary conditions by using the one parameter continuous group method. A convective boundary condition instead of the commonly used constant surface temperature or constant heat flux boundary conditions is applied. The governing boundary-layer equations are transformed to a two-point boundary value problem in similarity variables, and the problem is solved numerically by the Keller Box method. The effects of governing parameters on fluid velocity, temperature, and particle concentration are investigated and shown graphically.


## 2. Mathematical Formulation of the Problem

The problem of two-dimensional steady MHD heat and mass transfer laminar flow of a viscous incompressible and electrically conducting fluid past a flat plate is considered. Thex
- axis is taken along the plate and the y
- axis is normal to the plate. The gravitation acceleration vector is parallel to plate. A magnetic field of uniform strength B
0 is applied perpendicular to the direction of the plate. The viscosity, thermal conductivity, and concentration diffusivity of fluid are assumed to vary linearly. The top surface of the plate is kept at uniform temperature T
w which is assumed to be greater than the full stream temperature T
∞. The species concentration C
w at the surface is uniform and the full stream concentration is C
∞. The bottom surface of the plate is heated by convection from a hot fluid of temperature T
f which provides a heat transfer coefficient h
f. The induced magnetic field due to the motion of the electrically conducting fluid is negligible. This assumption is valid for small magnetic Reynolds numbers. It is also assumed that the external electric field is zero and the electric field due to polarization of charges is negligible. It is also assumed that the pressure gradient and viscous and electrical dissipation are neglected. The physical configuration and schematic of the problem are shown in Figure 1. It is known that this is a type of Falkner-Skan flow.Figure 1
Physical configuration and the coordinate system.Furthermore the following assumptions are considered: (i) fluid has constant kinematic viscosity and thermal diffusivity and the Boussinesq approximation may be adopted for steady laminar flow, (ii) the particle diffusivity is constant, (iii) the concentration of particles is sufficiently dilute that particle coagulation in the boundary layer is negligible, and (iv) the magnetic Reynolds number is small so that the induced magnetic field is negligible in comparison with the applied magnetic field. Under these assumptions the governing Prandtl boundary-layer equations in dimensional form are as follows (see Kays et al. [20] and White [21]):
(1)
∂
u
-
∂
x
-
+
∂
v
-
∂
y
-
=
0
,
(2)
u
-
∂
u
-
∂
x
-
+
v
-
∂
u
-
∂
y
-
=
u
-
e
d
u
-
e
d
x
-
+
1
ρ
∞
∂
∂
y
-
[
μ
(
T
)
∂
u
-
∂
y
-
]
+
g
β
T
(
T
-
T
∞
)
+
g
β
C
(
C
-
C
∞
)
-
σ
B
0
2
ρ
(
u
-
-
u
-
e
)
,
(3)
u
-
∂
T
∂
x
-
+
v
-
∂
T
∂
y
-
=
1
ρ
∞
c
p
∂
∂
y
-
[
κ
(
T
)
∂
T
∂
y
-
]
,
(4)
u
-
∂
C
∂
x
-
+
v
-
∂
C
∂
y
-
=
∂
∂
y
-
(
D
(
C
)
∂
C
∂
y
-
)
,
where u
- and v
- are the velocities in the x
- and y
- directions, respectively, T is the temperature within the boundary layer, T
∞ is the temperature far away from the plate, C is the species concentration, g is the acceleration due to gravity, β
T is the volumetric coefficient of thermal expansion, β
C is the volumetric coefficient of concentration expansion, α is the thermal conductivity, and D is the molecular diffusivity. The respective boundary conditions are
(5)
v
-
=
0
,
u
-
=
N
1
ν
∂
u
-
∂
y
-
,
T
=
T
w
+
D
1
∂
T
∂
y
-
,
C
=
C
w
at
y
-
=
0
,
u
-
⟶
u
-
e
,
T
⟶
T
∞
,
g
g
g
g
g
g
g
g
g
g
i
C
⟶
C
∞
as
y
-
⟶
∞
,
where u
-
e is the velocity over the plate that should be in the form u
-
e
=
c
x
m. This condition will be imposed later. μ
(
T
), κ
(
T
), and D
(
C
) are variable viscosity, thermal conductivity, and molecular diffusivity, respectively; the dimensions of N
1 are (
velocity
)
-
1 and the dimension of D
1 is length. It is assumed that the temperature dependent viscosity and thermal conductivity vary linearly and are given by (see Seddeek and Salem [22])
(6)
μ
(
T
)
=
μ
∞
[
1
+
b
1
(
T
f
-
T
)
]
,
κ
(
T
)
=
κ
∞
[
1
+
c
(
T
-
T
∞
)
]
,
where μ
∞ and κ
∞ are the constant undisturbed viscosity and undisturbed thermal conductivity, b
1 is a constant with b
1
>
0, and c is a constant which depends on the fluid. It is assumed that the concentration diffusivity varies linearly and is given by (see Seddeek and Salem [22])
(7)
D
(
C
)
=
D
m
[
1
+
c
(
C
-
C
∞
)
]
=
D
m
[
1
+
D
c
ϕ
]
,
where D
m is the constant concentration diffusivity.The following dimensionless variables are introduced:(8)
x
=
x
-
L
,
y
=
y
-
L
Re
1
/
2
,
u
=
u
-
U
∞
,
v
=
v
-
U
∞
Re
1
/
2
,
θ
=
T
-
T
∞
T
f
-
T
∞
,
ϕ
=
C
-
C
∞
C
w
-
C
∞
,
where Re is the Reynolds number, L is the characteristic length, θ is the dimensionless temperature variable, and ϕ is the dimensionless concentration variable. Introducing the stream function ψ such thatu
=
∂
ψ
/
∂
y and v
=
-
∂
ψ
/
∂
x, continuity equation (1) is satisfied identically and (2)–(4) now yield
(9)
Δ
1
≡
∂
ψ
∂
y
∂
2
ψ
∂
x
∂
y
-
∂
ψ
∂
x
∂
2
ψ
∂
y
2
-
u
e
d
u
e
d
x
-
-
(
a
+
A
(
1
-
θ
)
)
∂
3
ψ
∂
y
3
+
A
∂
2
ψ
∂
y
2
∂
θ
∂
y
-
g
(
T
f
-
T
∞
)
L
U
∞
2
β
T
θ
-
g
(
C
w
-
C
∞
)
L
U
∞
2
β
C
ϕ
+
M
(
∂
ψ
∂
y
-
u
e
)
=
0
,
Δ
2
≡
∂
ψ
∂
y
∂
θ
∂
x
-
∂
ψ
∂
x
∂
θ
∂
y
-
1
Pr
[
1
+
S
θ
]
∂
2
θ
∂
y
2
-
1
Pr
S
(
∂
θ
∂
y
)
2
=
0
,
Δ
3
≡
∂
ψ
∂
y
∂
ϕ
∂
x
-
∂
ψ
∂
x
∂
ϕ
∂
y
-
1
Sc
[
1
+
D
c
ϕ
]
∂
2
ϕ
∂
y
2
-
1
Sc
D
c
(
∂
ϕ
∂
y
)
2
=
0
.The boundary conditions are(10)
∂
ψ
∂
x
=
0
,
∂
ψ
∂
y
=
a
∂
2
ψ
∂
y
2
,
θ
=
1
+
b
∂
θ
∂
y
,
h
h
h
h
h
h
h
h
h
h
h
h
h
i
h
h
h
h
h
h
h
ϕ
=
1
at
y
=
0
,
∂
ψ
∂
y
⟶
u
e
(
x
)
,
θ
⟶
0
,
ϕ
⟶
0
as
y
⟶
∞
.In the above equations the parameters are defined as(11)
Re
=
U
∞
L
ν
,
Sc
=
ν
D
,
M
=
σ
B
0
2
L
ρ
U
∞
,
Pr
=
μ
c
p
κ
,
A
=
b
1
(
T
f
-
T
∞
)
,
S
=
c
(
T
f
-
T
w
)
,
a
=
N
1
ν
L
Re
,
b
=
D
1
Re
L
,
D
c
=
c
(
C
w
-
C
∞
)
,
where Re is the Reynolds number, Sc is the Schmidt number, M is the magnetic parameter, Pr is the Prandtl number of the fluid, S is the thermal conductivity parameter, A is the viscosity parameter, D
c is the concentration diffusivity parameter, a is the velocity slip parameter, and b is the thermal slip parameter.


## 3. Application of Group Transformations

Determining similarity solutions of (9)-(10) is equivalent to determining invariant solutions of these equations under a particular continuous one parameter group (Hamad et al. [23] and Kandasamy et al. [24]). Thus we search for a transformation group from the elementary set of one-parameter scaling transformations as one of the techniques that are defined by the following group which is called G1:
(12)
G
1
:
x
*
=
x
e
ε
α
1
,
y
*
=
y
e
ε
α
2
,
ψ
*
=
ψ
e
ε
α
3
,
θ
*
=
θ
e
ε
α
4
,
ϕ
*
=
ϕ
e
ε
α
5
,
β
T
*
=
β
T
e
ε
α
6
,
β
C
*
=
β
C
e
ε
α
7
,
u
e
*
=
u
e
e
ε
α
8
.Hereε(≠
0) is a parameter of the group and the α’s are arbitrary real numbers whose connection will be determined by our analysis. The transformations listed in (12) may be treated as point transformations which transform the coordinates
(13)
(
x
,
y
,
ψ
,
θ
,
ϕ
,
β
T
,
β
C
,
u
e
)
to
(
x
*
,
y
*
,
ψ
*
,
θ
*
,
ϕ
*
,
β
T
*
,
β
C
*
,
u
e
*
)
.The system (9)-(10) remains invariant under the group transformation G1. Hence we have the following relationships among the parameters, namely,
(14)
2
α
3
-
α
1
-
2
α
2
=
2
α
8
-
α
1
=
α
3
-
3
α
2
=
α
3
-
α
2
=
α
3
+
α
4
-
3
α
2
=
α
8
=
α
4
+
α
6
=
α
5
+
α
7
,
α
3
+
α
4
-
α
1
-
α
2
=
α
4
-
2
α
2
=
2
α
4
-
2
α
2
,
α
3
+
α
5
-
α
1
-
α
2
=
α
5
-
2
α
2
=
2
α
5
-
2
α
2
.From boundary conditions (10), these will be invariant if
(15)
α
2
-
α
3
=
2
α
2
-
α
3
,
-
α
4
=
α
2
-
α
4
,
α
5
=
0
,
α
2
-
α
3
=
-
α
8
.Solving (14) and (15), we obtain
(16)
α
2
=
α
4
=
α
5
=
0
,
α
1
=
α
3
=
α
6
=
α
7
=
α
8
.With these relations the boundary conditions remain invariant.The set of transformations G1 in (12) then reduces to
(17)
x
*
=
x
e
ε
α
1
,
y
*
=
y
,
ψ
*
=
ψ
e
ε
α
1
,
θ
*
=
θ
,
ϕ
*
=
ϕ
,
β
T
*
=
β
T
e
ε
α
1
,
β
C
*
=
β
C
e
ε
α
1
,
u
e
*
=
u
e
e
ε
α
1
.Using a Taylor series expansion in powers ofε, retaining terms up to first order, and neglecting higher powers of ε results in
(18)
x
*
-
x
=
ε
α
1
x
,
y
*
-
y
=
0
,
ψ
*
-
ψ
=
ε
α
1
ψ
,
θ
*
-
θ
=
0
,
ϕ
*
-
ϕ
=
0
,
β
T
*
-
β
T
=
ε
α
1
β
T
,
β
C
*
-
β
C
=
ε
α
1
β
C
,
u
e
*
-
u
e
=
ε
α
1
u
e
.The characteristic equations are(19)
d
x
α
1
x
=
d
y
0
=
d
ψ
α
1
ψ
=
d
θ
0
=
d
ϕ
0
=
d
β
T
α
1
β
T
=
d
β
C
α
1
β
C
=
d
u
e
α
1
u
e
.Solving the above characteristic equations gives(20)
η
=
y
,
ψ
=
x
f
(
η
)
,
θ
=
θ
(
η
)
,
ϕ
=
ϕ
(
η
)
,
β
T
=
β
T
0
x
,
β
C
=
β
C
0
x
,
u
e
=
U
∞
x
.Substituting (20) into (9)-(10) yields
(21)
[
1
+
A
(
1
-
θ
)
]
f
′′′
+
(
f
-
A
θ
′
)
f
′′
-
f
′
2
-
M
(
f
′
-
1
)
+
1
+
Gr
θ
+
Gc
ϕ
=
0
,
[
1
+
S
θ
]
θ
′′
+
S
θ
′
2
+
Pr
f
θ
′
=
0
,
[
1
+
D
c
ϕ
]
ϕ
′′
+
D
c
ϕ
′
2
+
Sc
f
ϕ
′
=
0
.HereGr
=
g
(
T
f
-
T
∞
)
L
β
T
0
/
U
∞
2and Gc
=
g
(
C
w
-
C
∞
)
L
β
c
0
/
U
∞
2 are Grashof numbers based on temperature and on concentration, respectively.The corresponding boundary conditions are(22)
f
(
0
)
=
0
,
f
′
(
0
)
=
a
f
′
′
(
0
)
,
θ
(
0
)
=
1
+
b
θ
′
(
0
)
,
ϕ
(
0
)
=
1
,
f
′
(
∞
)
⟶
1
,
θ
(
∞
)
⟶
0
,
ϕ
(
∞
)
⟶
0
.To obtain a similarity solution for the energy equation, the quantityb must be independent of x
- and for this to occur the heat transfer coefficient h
f must be directly proportional to (
x
-
)
-
1
/
2.


### 3.1. Parameters of Physical Interest

We are interested in the friction factorC
f, Nusselt number Nu, and Sherwood number Sh, respectively. Physically, C
f indicates wall shear stress and Nu indicates the rate of heat transfer whilst Sh indicates the rate of mass transfer. These quantities may be conveniently determined from
(23)
C
f
=
μ
ρ
U
∞
2
(
∂
u
-
∂
y
-
)
y
-
=
0
,
Nu
=
-
x
-
T
w
-
T
∞
(
∂
T
∂
y
-
)
y
-
=
0
,
Sh
=
-
x
-
C
w
-
C
∞
(
∂
C
∂
y
-
)
y
-
=
0
.By substituting (10) and (19) into (23), we obtain
(24)
Re
1
/
2
C
f
=
[
1
+
A
(
1
-
θ
(
0
)
)
]
f
′
′
(
0
)
,
Re
-
1
/
2
Nu
=
-
θ
′
(
0
)
,
Re
-
1
/
2
Sh
=
-
ϕ
′
(
0
)
.From (24) it can be shown that the skin friction factor C
f, the Nusselt number Nu, and the Sherwood number Sh are proportional to the numerical values f
′
′
(
0
), -
θ
′
(
0
), and -
ϕ
′
(
0
), respectively.


## 3.1. Parameters of Physical Interest

We are interested in the friction factorC
f, Nusselt number Nu, and Sherwood number Sh, respectively. Physically, C
f indicates wall shear stress and Nu indicates the rate of heat transfer whilst Sh indicates the rate of mass transfer. These quantities may be conveniently determined from
(23)
C
f
=
μ
ρ
U
∞
2
(
∂
u
-
∂
y
-
)
y
-
=
0
,
Nu
=
-
x
-
T
w
-
T
∞
(
∂
T
∂
y
-
)
y
-
=
0
,
Sh
=
-
x
-
C
w
-
C
∞
(
∂
C
∂
y
-
)
y
-
=
0
.By substituting (10) and (19) into (23), we obtain
(24)
Re
1
/
2
C
f
=
[
1
+
A
(
1
-
θ
(
0
)
)
]
f
′
′
(
0
)
,
Re
-
1
/
2
Nu
=
-
θ
′
(
0
)
,
Re
-
1
/
2
Sh
=
-
ϕ
′
(
0
)
.From (24) it can be shown that the skin friction factor C
f, the Nusselt number Nu, and the Sherwood number Sh are proportional to the numerical values f
′
′
(
0
), -
θ
′
(
0
), and -
ϕ
′
(
0
), respectively.


## 4. The Keller Box Method

Equation (21) subject to boundary conditions (22) is solved numerically using a very efficient finite difference scheme known as the Keller Box method. The details of this method are described in Cebeci and Bradshaw [25] and Na [26]. For more information refer to Keller [27, 28].


## 5. Results and Discussion

Applying scaling group transformations to analyze the governing equations and the boundary conditions, the two independent variables are reduced by one. Consequently the governing equations reduce to a system of nonlinear ordinary differential equations with the appropriate boundary conditions. The transformed momentum, energy, and concentration equation (21) subject to the boundary conditions (22) were solved numerically by using the Keller Box method. We obtained velocity, temperature, and concentration profile graphs for different values of governing parameters.Figures2, 3, and 4 show the effects of the viscosity parameter and the magnetic parameter on the velocity, temperature, and concentration distributions, respectively. The velocity distribution decreases with increasing A and M whereas they have no significant effect on the temperature and concentration distributions. This behavior can be predicted from (21) and also the physical definition of parameters A and M, since the viscosity and magnetic parameters only appear in the momentum equation. Figures 5, 6, and 7 show the effect of the thermal and mass Grashof numbers on the velocity, temperature, and concentration distributions, respectively. Physically, since the thermal Grashof number (
Gr
) is the ratio of buoyancy to viscous forces in the boundary layer, increasing its value suggests an increase in the buoyancy forces relative to the viscous forces and this is clearly reflected in the progressive increase in the velocity of the flow. Increase in the mass transfer Grashof number (
Gc
) yields a similar effect on the velocity of the flow. Moreover, the reverse trend is seen for the temperature and concentration distributions.Figure 2
Effects of the viscosity parameterA and the magnetic parameter M on the dimensionless velocity.Figure 3
Effects of the viscosity parameterA and the magnetic parameter M on the dimensionless temperature.Figure 4
Effects of the viscosity parameterA and the magnetic parameter M on the dimensionless concentration.Figure 5
Effects of the thermal Grashof numberGr and the mass Grashof number Gc on the dimensionless velocity.Figure 6
Effects of the thermal Grashof numberGr and the mass Grashof number Gc on the dimensionless temperature.Figure 7
Effects of the thermal Grashof numberGr and the mass Grashof number Gc on the dimensionless concentration.Figures8, 9, and 10 illustrate the influence of the thermal conductivity parameter S and Prandtl number Pr on the velocity, temperature, and concentration distributions, respectively. It is observed that the velocity and temperature distributions increase with increasing thermal conductivity parameter and decrease with increasing Prandtl number. This is in agreement physically since the thermal boundary-layer thickness decreases with increasing Pr. The thermal conductivity parameter S and the Prandtl number Pr have no significant effect on the concentration distribution and this can be predicted from (21). The effects of the concentration diffusivity parameter D
c and the Schmidt number Sc on the velocity, temperature, and concentration distributions are shown in Figures 11–13. The velocity and concentration distributions increase with increasing concentration diffusivity parameter whereas they decrease with increasing Schmidt number. Since Schmidt number is the ratio of viscosity to diffusivity, this behavior can be predicted. From Figure 12, the concentration diffusivity parameter and the Schmidt number have no significant effect on the temperature distribution. In Figures 14–16 effects of the velocity slip parameter (
a
) and the thermal slip parameter (
b
) are depicted. In Figure 14 it is observed that velocity distribution increases with increasing velocity slip parameter and decreases with increasing thermal slip parameter. From Figure 15 we observe that as the velocity slip parameter and the thermal slip parameter increase the temperature distribution decreases. Figure 16 shows that the concentration distribution decreases with increasing velocity slip parameter and increases with increasing thermal slip parameter. In some of the velocity profiles an overshoot of the velocity profile is observed. This depends on the boundary conditions. In other words, since there is not a no-slip condition on the plate, a larger velocity rather than free stream velocity can exist inside the boundary layer. With attention to boundary conditions (5) an overshoot of the velocity is likely.Figure 8
Effects of the thermal conductivity parameterS and the Prandtl number Pr on the dimensionless velocity.Figure 9
Effects of the thermal conductivity parameterS and the Prandtl number Pr on the dimensionless temperature.Figure 10
Effects of the thermal conductivity parameterS and the Prandtl number Pr on the dimensionless concentration.Figure 11
Effects of the concentration diffusivity parameterD
c and the Schmidt number Sc on the dimensionless velocity.Figure 12
Effects of the concentration diffusivity parameterD
c and the Schmidt number Sc on the dimensionless temperature.Figure 13
Effects of the concentration diffusivity parameterD
c and the Schmidt number Sc on the dimensionless concentration.Figure 14
Effects of the velocity slip parametera and the thermal slip parameter b on the dimensionless velocity.Figure 15
Effects of the velocity slip parametera and the thermal slip parameter b on the dimensionless temperature.Figure 16
Effects of the velocity slip parametera and the thermal slip parameter b on the dimensionless concentration.In Table1 choosing A
=
1.0, S
=
1.0, D
c
=
1.0, Sc
=
1.0; a
=
1.0, and b
=
1.0, numerical values of f
′′
(
0
), θ
′
(
0
), and φ
′
(
0
) are shown for different values of the parameters Gr, Gc, M, and Pr. Results of Figures 2–16 are verified.Table 1
Numerical results off
′
′(0), θ
′(0), and φ
′(0)  for different values of the parameters Gr, Gc, M, and Pr when A
=
1.0, S
=
1.0, D
c
=
1.0, Sc
=
1.0;  a
=
1.0, and b
=
1.0.
Gr
Gc
M
Pr
f
′
′(0)
-
θ
′(0)
-
φ
′(0)
1
1
1
0.72
0.835798
0.33
0.51803
0
1
1
0.72
0.751505
0.32094
0.49805
2
1
1
0.72
0.912888
0.33783
0.53563
1
0
1
0.72
0.713589
0.31708
0.48958
1
2
1
0.72
0.947207
0.34088
0.54264
1
1
0
0.72
0.842129
0.33159
0.52146
1
1
2
0.72
0.833216
0.32895
0.51581
1
1
1
0.3
0.87512
0.23957
0.52956
1
1
1
1
0.820823
0.36801
0.51389


## 6. Conclusions

A numerical study based on the Keller Box method for MHD heat and mass transfer flow of a steady viscous incompressible fluid over a flat plate has been performed. We have investigated the effects of various governing parameters, namely, the viscosity parameterA, the magnetic field M, thermal Grashof number Gr, mass transfer Grashof number Gc, thermal conductivity parameter S, Prandtl number Pr, concentration diffusivity parameter D
c, Schmidt number Sc, velocity slip parameter a, and thermal slip parameter b on flow field and heat transfer characteristics. The following conclusions can be made.(1)
The thickness of the velocity boundary layer decreases with an increase in viscosity parameterA, magnetic field M, Schmidt number Sc, and thermal slip parameter b.
(2)
The thickness of the velocity boundary layer increases with an increase in thermal Grashof numberGr, mass transfer Grashof number Gc, thermal conductivity parameter S, concentration diffusivity parameter D
c, and velocity slip parameter a.
(3)
The thickness of the thermal boundary layer decreases with an increase in thermal Grashof numberGr, mass transfer Grashof number Gc, Prandtl number Pr, velocity slip parameter a, and thermal slip parameter b.
(4)
The thickness of the thermal boundary layer increases with an increase in thermal conductivity parameterS.
(5)
The thickness of the concentration boundary layer decreases with an increase in thermal Grashof numberGr, mass transfer Grashof number Gc, Schmidt number Sc, and velocity slip parameter a.
(6)
The thickness of the concentration boundary layer increases with an increase in concentration diffusivity parameterD
c and thermal slip parameter b.



---

*Source: 923487-2014-06-16.xml*