# Slip Effect on MHD Chemically Reacting Convictive Boundary Layer Flow with Heat Absorption

Authors: Mekonnen Shifferaw Ayano; Negussie Tadege Demeke
Journal: Journal of Engineering
 (2016)

---

## Abstract

The aim of this paper is to investigate steady magneto micropolar fluid past a stretched semi-infinite vertical and permeable surface taking into account heat absorption, hall and ion-slip effect, first-order chemical reaction, and radiation effects. The system of coupled nonlinear equations is solved numerically and the effects of various parameters on the velocity, the microrotation, the temperature, and the concentration field are calculated. The following discovery was made: an increase in the hall parameter strongly enhances the fluid velocity, an increase in the heat absorption parameter increases the temperature, and an increase in the slip parameter decreases the velocity. Additionally, the local skin friction, Nusselt number, and Sherwood number are also analyzed for various parameters and presented in table form.


---

## Body

## 1. Introduction

The micropolar fluids are non-Newtonian fluids with microstructures such as polymeric additives, colloidal suspensions, and liquid crystals. Erigen [1, 2] first developed the theory of micropolar fluids and thermomicropolar fluids to describe such characteristics. A thorough review of the subject and applications of micropolar fluid mechanics was provided by Ariman et al. [3, 4]. These fluids are realistic and important from a technological point of view.Due to the vast applications of the heat and mass transfer over a vertical channel in geophysics and engineering, a lot of authors have studied and published about this subject. Considering the chemical reaction and radiation effects, the common applications that can be found are the combustion of fossil fuels, atmospheric reentry with suborbital velocities, plasma wind tunnels, electric spacecraft propulsion, hypersonic flight through planetary atmosphere photo-dissociation, and photo ionization. Articles have been published based on the study of heat generation or absorption, chemical reaction, and radiation effects. Studies have been done on the individual or collective effects. For instance, Muthucumaraswamy and Ganesan [5] studied numerically natural convection on a moving isothermal vertical plate with chemical reaction. In their study, they have shown that velocity increases during generative reaction and decreases in destructive reaction. Vajravelu and Hadjinicolaou [6] studied the heat transfer in viscous fluid over a stretching sheet with viscous dissipation and internal heat generation. Kelson and Desseaux [7] examined the effect of surface conditions on the flow of a micropolar fluid driven by a porous stretching sheet. Takhar et al. [8] have investigated the flow and mass diffusion of chemical species with first-order and higher order reactions over a continuously stretching sheet with the magnetic field effect. Influence of chemical reaction on heat and mass transfer by natural convection from vertical surfaces in porous media considering soret and dufour effects was analyzed by Postelnicu [9]. El-Arabawy [10] analyzed the problem of the effect of suction/injection on the flow of a micropolar fluid past a continuously moving plate in the presence of radiation.The effect of magnetic field with or without hall and ion-slip effects has been studied by several researchers. On uniformly stretched vertical permeable surface in the presence of heat generation/absorption and a chemical reaction, Chamkha [11] studied analytically the magnetic effects. Mahmoud [12] considered thermal radiation effects on MHD flow of a micropolar fluid over a stretching surface with variable thermal conductivity. Chemical reaction and heat and mass transfer on MHD flow over a vertical stretching surface with heat source and thermal stratification effects were examined by Kandasamy et al. [13]. Khedr et al. [14] examined numerically MHD flow of a micropolar fluid past a stretched permeable surface with heat generation or absorption neglecting the hall and ion-slip effects.The abovementioned studies addressed flows when one or more of the effects of chemical reaction, radiation, heat source, and MHD are absent. The objective of this paper is to analyze hall and ion-slip effects of micropolar fluid along a vertical semi-infinite permeable plate taking into account the slip boundary conditions with the combined effects of chemical reaction, radiation, and heat source with wall suction or injection effects. The governing systems of partial differential equations have been solved numerically by the Keller box method. The effects of the various parameters on the velocity, micrototation, temperature, and concentration field will be investigated.


## 2. Mathematical Formulation of Problem

Let us consider a steady, laminar, incompressible micropolar fluid flowing past a porous surface along a semi-infinite vertical plate in the presence of heat generation or absorption and thermal radiation effect. Further, it is assumed that there exists a first-order homogeneous chemical reaction between the species and the fluid, a magnetic field of a uniform strengthB
0 is applied transversely to the plate, and the magnetic Reynolds number is assumed to be small so that the induced magnetic field is neglected. When the fluid is assumed to be viscous and has constant properties retaining the hall and ion-slip effect, the equation governing the flow is given by(1)
∂
u
∂
x
+
∂
v
∂
y
=
0
,
(2)
ρ
u
∂
u
∂
x
+
v
∂
u
∂
y
=
μ
+
κ
∂
2
u
∂
y
2
+
κ
∂
ϖ
∂
y
-
σ
B
0
2
α
e
2
+
B
e
2
B
e
u
+
α
e
w
,
(3)
ρ
u
∂
w
∂
x
+
v
∂
w
∂
y
=
μ
+
κ
∂
2
w
∂
y
2
+
σ
B
0
2
α
e
2
+
B
e
2
α
e
u
-
B
e
w
,
(4)
-
2
κ
ϖ
-
κ
∂
u
∂
y
+
ω
∂
2
ϖ
∂
y
2
=
0
,
(5)
u
∂
T
∂
x
+
v
∂
T
∂
y
=
α
∂
2
T
∂
y
2
+
ν
C
p
∂
u
∂
y
2
-
1
C
p
∂
q
r
∂
y
+
Q
ρ
C
p
T
-
T
∞
,
(6)
u
∂
C
∂
x
+
v
∂
C
∂
y
=
D
∂
2
C
∂
y
2
-
K
1
C
-
C
∞
,where K
1 is chemical reaction parameter, B
0 is magnetic field intensity, D is the effective diffusive coefficient, ρ is density, u and w are fluid velocities, ϖ is microrotation component, T is temperature, C is concentration, μ is the dynamic viscosity, α
e
=
1
+
B
e
B
i is hall parameter, B
e is ion-slip parameter, α is the fluid thermal diffusivity, σ is fluid electrical conductivity, Q is dimensional heat generation (>0) or absorption (<0) coefficient, and q
r is the radiative heat flux.The boundary conditions are as follows:(7)
u
=
u
w
+
L
∂
u
∂
y
,
v
=
0
,
w
=
0
,
ϖ
=
-
n
∂
u
∂
y
,
T
=
T
w
,
C
=
C
-
w
at
y
=
0
,
u
⟶
0
,
w
⟶
0
,
ϖ
⟶
0
,
T
=
T
∞
,
C
=
C
∞
,
as
y
⟶
∞
,where L is the proportionality constant of the velocity. For microrotation boundary condition, n
=
0 corresponds to no slip, the concentrated particle where the microelements close to the plate surface are unable to rotate. And the case n
=
1
/
2 represents the fact that antisymmetric part of the stress is zero on the surface and weak concentration. The case n
=
1 represents turbulent boundary layer flow.In (5), the radiative heat flux under Rosseland approximation is given by (8)
q
r
=
-
σ
∗
3
k
∗
∂
T
4
∂
y
,where σ
∗ is Stephan-Boltzmann constant and k
∗ is mean absorption coefficient. If the temperature differences within the flow are sufficiently small, T
4 may be expanded by Taylor series. Hence, expanding T
4 about T
∞ and neglecting higher order terms, we get (9)
T
4
≅
T
∞
4
+
4
T
-
T
∞
T
∞
3
=
4
T
∞
3
T
-
3
T
∞
3
.Introduce the nondimensional variables through (10)
ψ
=
2
ν
x
U
0
f
η
,
u
=
∂
ψ
∂
y
=
U
0
f
′
η
,
v
=
-
∂
ψ
∂
x
=
-
U
0
ν
2
x
f
η
-
η
f
′
η
,
η
=
y
U
0
2
ν
x
,
w
=
U
0
h
η
,
ϖ
=
U
0
3
2
ν
x
g
η
,
θ
η
=
T
-
T
∞
T
w
-
T
∞
,
ϕ
η
=
C
-
C
∞
C
w
-
C
∞
.Substituting (8)–(10) into (1)–(6), we get the following nondimensional equations: (11)
f
′
′
′
+
f
f
′
+
N
1
g
′
-
M
2
α
e
2
+
B
e
2
α
e
f
′
+
B
e
h
=
0
,
(12)
h
′
′
+
1
2
h
′
f
+
M
2
α
e
2
+
B
e
2
B
e
f
′
-
α
e
h
=
0
,
(13)
λ
g
′
′
-
2
2
g
+
f
′
′
=
0
,
(14)
3
Ra
+
4
3
PrRa
θ
′
′
+
f
θ
′
+
A
θ
+
E
c
f
′
′
2
=
0
,
(15)
1
S
c
ϕ
′
′
+
f
ϕ
′
-
K
ϕ
=
0
,where primes denote differentiation with respect to the variable η, N
1
=
κ
/
ρ
ν is coupling number, A
=
2
Q
/
C
p
ρ
U
0 is internal heat generation/absorption parameter, E
c
=
U
0
2
/
c
p
(
T
w
-
T
∞
) is the Eckert number, Pr
⁡
=
ν
C
p
/
α is the Prandtl number, M
=
2
σ
B
0
2
/
ρ
U
0 is Hartmann number, Γ
=
γ
U
0
/
κ
ν is microrotation parameter, and Ra
=
k
k
∗
/
4
σ
∗
T
∞
3 is radiation parameter.The corresponding boundary conditions in dimensionless form are as follows:(16)
f
0
=
S
0
,
f
′
0
=
1
+
ε
f
′
′
0
,
h
0
=
0
,
g
0
=
-
n
f
′
′
0
,
θ
0
=
1
,
ϕ
0
=
1
,
(17)
η
⟶
∞
f
′
⟶
0
,
h
⟶
0
,
g
⟶
0
,
θ
⟶
0
,
ϕ
⟶
0
,where ε is slip parameter and S
0
=
-
V
0
ν
U
0
/
2
x is suction/injection velocity.The shearing stress component at the plate can be calculated in nondimensional form:(18)
C
f
=
μ
+
κ
∂
u
/
∂
y
y
=
0
1
/
2
ρ
U
0
2
=
-
2
R
e
x
1
/
2
f
′
′
0
.The local mass flux and the local Sherwood number are, respectively, given by(19)
Nu
x
=
-
x
∂
T
/
∂
y
y
=
0
T
w
-
T
∞
=
-
2
R
e
x
1
/
2
θ
′
0
,
Sh
x
=
-
x
∂
C
/
∂
y
y
=
0
C
w
-
T
∞
=
-
2
R
e
x
1
/
2
ϕ
′
0
,where R
e
x
=
U
0
x
/
ν is the Reynolds number.


## 3. Numerical Scheme

The numerical scheme to solve (11)–(16) is outlined based on Keller [15] as follows: (11)–(15) are written as a first-order system of equation. For this, we define new variables U
,
V
,
P
,
Q
,
R
, and S as(20)
f
′
=
U
,
u
′
=
V
,
h
′
=
S
,
g
′
=
P
,
θ
′
=
Q
,
ϕ
′
=
R
,
V
′
+
f
V
-
n
1
P
-
M
2
α
e
2
+
B
e
2
α
e
U
+
B
e
h
=
0
,
h
′
′
+
f
S
+
M
2
α
e
2
+
B
e
2
B
e
U
-
α
e
h
=
0
,
λ
P
′
-
2
2
g
+
V
=
0
,
3
Ra
+
4
3
RaPr
Q
′
′
+
f
Q
+
A
θ
+
E
c
V
2
=
0
,
1
S
c
R
′
f
-
K
ϕ
=
0
.The boundary condition becomes(21)
η
=
0
,
f
=
0
,
U
=
0
,
h
=
0
,
g
=
0
,
θ
=
1
,
ϕ
=
1
,
U
⟶
0
,
h
⟶
0
,
g
⟶
0
,
θ
⟶
0
,
ϕ
⟶
0
.We write the difference equation: (22)
η
0
=
0
,
η
j
=
η
j
-
1
+
ξ
j
,
j
=
1,2
,
3
,
…
,
J
,
η
J
=
η
∞
,
η
j
-
1
/
2
=
1
2
η
j
+
η
j
-
1
,
j
-
1
/
2
=
1
2
j
+
j
-
1
,where ξ
j is Δ
η spacing. Therefore,(23)
ξ
j
-
1
f
j
-
f
j
-
1
=
U
j
-
1
/
2
,
ξ
j
-
1
U
j
-
U
j
-
1
=
V
j
-
1
/
2
,
ξ
j
-
1
g
j
-
g
j
-
1
=
P
j
-
1
/
2
,
ξ
j
-
1
θ
j
-
θ
j
-
1
=
Q
j
-
1
/
2
,
ξ
j
-
1
ϕ
j
-
ϕ
j
-
1
=
R
j
-
1
/
2
,
V
j
-
V
j
-
1
+
ξ
j
4
f
j
+
f
j
-
1
V
j
+
V
j
-
1
+
N
1
ξ
j
2
P
j
+
P
j
-
1
-
ξ
j
M
2
2
α
e
2
+
B
e
2
α
e
U
j
+
U
j
-
1
+
B
e
h
j
+
h
j
-
1
,
S
j
-
S
j
-
1
+
ξ
j
4
f
j
+
f
j
-
1
S
j
+
S
j
-
1
+
ξ
j
M
2
2
α
e
2
+
B
e
2
B
e
h
j
+
h
j
-
1
-
α
e
E
j
+
E
j
-
1
,
λ
P
j
-
P
j
-
1
-
2
ξ
V
j
+
V
j
-
1
+
2
g
j
+
g
j
-
1
=
0
,
3
Ra
+
4
3
RaPr
Q
j
-
Q
j
-
1
+
ξ
j
4
f
j
+
f
j
-
1
Q
j
+
Q
j
-
1
+
A
ξ
j
2
θ
j
+
θ
j
-
1
+
ξ
j
4
f
j
+
f
j
-
1
2
=
0
,
1
S
c
S
j
-
S
j
-
1
+
ξ
j
4
f
j
+
f
j
-
1
R
j
+
R
j
-
1
-
ξ
j
2
ϕ
j
+
ϕ
j
-
1
=
0
.Linearizing the resulting system of equations by Newton method, we will get a block tridiagonal structured matrix and then solve the system using Keller box method. This method has been proven to be adequate and give accurate results for boundary layer equations. The step size Δ
η
=
0.01 is chosen after testing with different step sizes. The calculations are repeated until some convergence criterion is satisfied and the calculations were stopped when δ
f
0
′
′
≤
1
0
-
8, δ
h
0
′
≤
1
0
-
8, δ
g
0
′
≤
1
0
-
8, δ
θ
0
′
≤
1
0
-
8,  and δ
ϕ
0
′
≤
1
0
-
8. In the present study, the boundary conditions for η at ∞ are replaced by a sufficiently large value of η where the velocity, microrotation, temperature, and concentration profiles approach zero. After some trials, we imposed a maximal value of η at ∞ of 10.


## 4. Results and Discussion

Computations have been performed using the following parameters:S
c
=
0.2, P
r
=
0.72, λ
=
2, E
c
=
0.02, Ra very large number, and N
1
=
0.5, which represents the presence of microelements; the fluid is micropolar (in (2), when κ
→
0 and when N
1
→
0 in (11), equations (11) and (13) decouple and the fluid behaves as Newtonian). To validate the code, results are compared with those of El-Arabawy [10] and it was found that they are in good agreement as shown in Table 1.Table 1
A comparison of values of Pr whenS
0
=
0, E
c
=
0.02, λ
=
2, A
=
0, R
a
=
1000, N
1
=
0.2, and K
=
0.
Pr
El-Arabawy [10]
Present
0.733
0.501327
0.50138
7
1.931150
1.93117
10
2.33700
2.33714
20.0
3.360750
3.36076The effects of the suction/injection parameter on velocities, microrotation, temperature, and concentration are presented in Figures1–5. It is observed that injecting a fluid at the plate (
S
0
<
0
) increases all the velocities, microrotation, temperature, and concentrations. On the other hand, the opposite is observed on all velocities, microrotation, temperature, and concentrations in the case of suction (
S
0
>
0
). Similar trend is observed on the boundary layer thicknesses.Figure 1
The effect of suction/injection onf
′.Figure 2
The effect of suction/injection onh.Figure 3
The effect of suction/injection ong.Figure 4
The effect of suction/injection onθ.Figure 5
The effect of suction/injection onϕ.Imposition of a magnetic field to an electrically conducting fluid creates a drag like force called the Lorentz force. As it is known, this force has the tendency to slow down the flow. The effect of the hall parameterα
e on the velocity profiles studied is shown in Figure 6. Figure 6 shows that velocity increases with increasing α
e and the effective conductivity decreases with increasing α
e which reduces the magnetic damping force on velocities. In general, increasing the hall and ion-slip parameters increases primary and secondary velocities as seen in Figures 6 and 7. When α
2
=
0, it is seen from Figure 7 that velocity component h is zero. The secondary velocity profiles h plotted in Figure 7 reveal that, increasing the values of the hall parameter α
e, a sharp rise in the velocity profile is observed near the plate. The effect of hall parameter on microrotation depicted in Figure 8 shows similar effect as the velocity component h is shown. The effect of hall and ion-slip however shows a reverse effect on the temperature and concentration as it is depicted in Figures 9 and 10, respectively.Figure 6
The effect ofB
e and α
e on f
′.Figure 7
The effect ofB
e and α
e on h.Figure 8
The effect ofB
e and α
e on g.Figure 9
The effect ofB
e and α
e on θ.Figure 10
The effect ofB
e and α
e on ϕ.The effect of the chemical reaction parameterK on the dimensionless concentration is shown in Figure 11. The chemical reaction parameter has a significant effect on the concentration distribution. Increasing the chemical reaction parameter decreases the concentration and the concentration boundary layer thickness.Figure 11
The effect ofK on ϕ.Our study focused on the caseE
c
=
0.02, which corresponds to plate cooling, that is, loss of heat from the plate to the fluid. The presence of heat absorption decreases in the thermal state of the fluid. Hence, from Figure 12, it is observed that the fluid temperature is noticeably enhanced with an increased A which also shows the increase in the thermal boundary layer thickness.Figure 12
The effect ofA on θ.The effects of the slip parameterε on velocity f
′ are shown in Figure 13. Increasing the sip parameter decreases the velocity, showing that the skin friction decreases at the surface. Figure 14 shows the nondimensional microrotation increases with increase in n. It is revealed that microrotation close to the boundary is unable to rotate for the case n
=
0.Figure 13
The effect ofε on f
′.Figure 14
The effect ofn on g.The boundary layer produces a drag on the plate due to the viscous stresses, which are developed at the wall, skin friction. The quantity of heat exchanged between the body and the fluid is given by the temperature gradientN
u (Nusselt number) and the rate of mass transfer S
h (Sherwood number) presented in Tables 2–7 for different values of the parameters. Increasing the suction/injection parameter increases all the skin friction, Nusselt number, and Sherwood numbers. Table 3 shows the effect of the ion-slip parameter. As the ion-slip parameter increases, the Nusselt number and Sherwood numbers increase. Table 4 shows similar effects for the hall parameter. The effect of the heat absorption parameter is shown in Table 5. It is observed that the Nusselt number decreases when the absorption parameter increases; however, there is no effect on the skin friction and Sherwood numbers as the equations are decoupled. The slip effect presented in Table 6 shows that as the slip parameter increases skin-friction, Nusselt number, and Sherwood numbers decrease. Increasing n, both N
u and S
h decrease whereas the skin friction C
f increases as shown in Table 7.Table 2
The effect of suction/injection on skin friction and heat and mass transfer coefficients.
S
0
f
′
′
(
0
)
θ
′
(
0
)
ϕ
′
(
0
)
−1
0.86861
0.46342
0.07382
−0.5
1.06557
0.59918
0.11458
0
1.30799
0.77337
0.16943
0.5
1.59560
0.98546
0.23824
1
1.92459
1.23034
0.31910Table 3
The effect ofB
e on skin friction and heat and mass transfer coefficients.
B
e
f
′
′
(
0
)
θ
′
(
0
)
ϕ
′
(
0
)
0
1.68842
0.74200
0.14730
1
1.52515
0.76014
0.16020
2
1.30799
0.77337
0.16943
3
1.16179
0.78229
0.17587
4
1.06222
0.78849
0.18054Table 4
The effect ofα
e on skin friction and heat and mass transfer coefficients.
α
e
f
′
′
(
0
)
θ
′
(
0
)
ϕ
′
(
0
)
0
2.97849
0.71014
0.13581
1
1.67501
0.75506
0.15781
2
1.30799
0.77337
0.16943
3
1.12814
0.78386
0.17680
4
1.01978
0.79066
0.18190Table 5
The effect ofA on skin friction and heat and mass transfer coefficients.
A
f
′
′
(
0
)
θ
′
(
0
)
ϕ
′
(
0
)
0
1.30799
0.35464
0.16943
−1
1.30799
0.84147
0.16943
−2
1.30799
1.11637
0.16943
−3
1.30799
1.33259
0.16943
−4
1.30799
1.51709
0.16943Table 6
The effect ofε on skin friction and heat and mass transfer coefficients.
ε
f
′
′
(
0
)
θ
′
(
0
)
ϕ
′
(
0
)
0
0.94209
0.51326
0.08851
0.2
0.85908
0.50550
0.08394
0.4
0.77743
0.49771
0.07942
0.7
0.65752
0.48599
0.07273
1
0.54063
0.47425
0.06618Table 7
The effect ofn on skin friction and heat and mass transfer coefficients.
n
f
′
′
(
0
)
θ
′
(
0
)
ϕ
′
(
0
)
0
0.94209
0.51326
0.08851
0.2
1.03079
0.50886
0.08490
0.4
1.11950
0.50444
0.08132
0.7
1.25257
0.49777
0.07602
1
1.38564
0.49107
0.07081


## 5. Conclusion

The steady incompressible chemically reacting micropolar fluid flow along a semi-infinite vertical plate in the presence of heat generation or absorption with heat source under uniform transverse magnetic field has been studied numerically. The effects of hall, ion, chemical reaction, slip, and suction/injection parameters on velocities, microrotation, temperature, and concentration have been investigated. The hall and ion-slip parameters increase the main velocity whereas the slip parameter decreases. The microrotation increases with increase in suction/injection, slip, hall, and ion-slip parameters. Temperature decreases with the increases in hall, ion-slip, and suction/injection parameters whereas it decreases with heat absorption parameters. The chemical reaction profile is highly influenced by the chemical reaction parameter. The skin friction at the wall decreases by increasing the hall, ion-slip, and slip parameters.



---

*Source: 2689493-2016-04-11.xml*