# Influence of Cross-Diffusion on Slip Flow and Heat Transfer of Chemically Reacting UCM Fluid between Porous Parallel Plates with Hall and Ion Slip Currents

Authors: Odelu Ojjela; Pravin Kashyap Kambhatla; N. Naresh Kumar; Samir Kumar Das
Journal: Journal of Computational Engineering
 (2016)

---

## Abstract

The present paper deals with the Hall and ion slip currents on an incompressible unsteady free convection flow and heat transfer of an upper convected Maxwell fluid between porous parallel plates with Soret and Dufour effects by considering the velocity slip and convective boundary conditions. Assume that there are periodic injection and suction at the lower and upper plates, respectively. The temperature and concentration at the lower and upper plates change periodically with time. The flow field equations are reduced to nonlinear ordinary differential equations by using similarity transformations and a semi-analytical-numerical solution has been obtained by the differential transform method. The velocity components, temperature distribution, and concentration with respect to different fluid and geometric parameters are discussed in detail and presented in the form of graphs. It is observed that the Biot number increases the temperature and concentration of the fluid. Further, the concentration of the fluid is enhanced whereas the temperature decreases with increasing slip. The present results are compared with the existing literature and are found to be in good agreement.


---

## Body

## 1. Introduction

The flow through porous channels is of great importance in both engineering and biological flows. Examples of this are found in aerodynamic heating, electrostatic precipitation, soil mechanics, food preservation, polymer technology, petroleum industry, the mechanics of the cochlea in the human ear, the flow of blood in the arteries, and artificial dialysis. The flows due to the periodic injection/suction at the boundary are of immense importance because of their wide range of applications. Several researchers have studied the incompressible laminar fluid flow between two porous parallel plates due to its mathematical simplicity and the flows can be analyzed easily in theoretical and experimental studies. Berman [1] examined the incompressible viscous fluid flow in a two-dimensional uniformly porous channel and a series solution was obtained for small values of the Reynolds number. Later, White Jr. et al. [2] considered the steady incompressible laminar viscous fluid flow between porous parallel plates with uniform suction or injection and the problem was analyzed for a wide range of suction Reynolds number. Walker and Davies [3] analyzed the mass transfer in an incompressible laminar viscous fluid flow between permeable parallel plates and a confluent hypergeometric function was used to obtain the solution. Hamza [4] investigated an incompressible laminar viscous fluid flow between two parallel rectangular as well as circular plates and obtained an approximate analytical solution by regular perturbation technique. The unsteady micropolar fluid flow between two parallel porous plates with periodic suction and injection was inspected by Srinivasacharya et al. [5]. The steady flow of chemically reacting micropolar fluid through a permeable channel was studied by Sheikholeslami et al. [6] and an analytical solution was obtained by using the homotopy perturbation method. Chen and Zhu [7] considered the slip effects on a Couette Poiseuille flow of a Bingham fluid between porous parallel plates and the reduced flow field equations are solved analytically. Singh and Pathak [8] analyzed the Hall currents and thermal radiation effects on unsteady MHD slip flow in a vertical rotating porous channel. Bhat and Mittal [9] inspected the Hall and ion slip currents on an electrically conducting viscous fluid flow through a parallel plate channel in the presence of uniform magnetic field and observed that the temperature of the fluid is reduced by Hall currents. Srinivasacharya and Kaladhar [10] examined analytically the Hall and ion slip currents on the mixed convective flow of couple stress fluid through a vertical parallel plate channel. Singh and Gorla [11] investigated the effects of thermal diffusion and heat source on free convection flow of a viscous fluid past an infinite vertical porous plate with Hall currents. Raptis et al. [12] studied the effects of Grashof number and permeability parameter on the free convection flow of a viscous fluid through a porous medium between two parallel plates. Abdulaziz and Hashim [13] investigated the heat and mass transfer of free convective flow of micropolar fluid in a vertical porous parallel plate channel and observed that the velocity decreases for increasing of Reynolds number. Makinde and Aziz [14] considered the convective boundary conditions for a MHD mixed convection fluid flow experiencing a first-order chemical reaction over a vertical plate embedded in a porous medium and a numerical solution is obtained by using shooting technique along with a sixth-order Runge-Kutta integration. Yao et al. [15] analyzed an incompressible laminar viscous fluid flow and heat transfer over a stretching/shrinking sheet with convective boundary conditions and have found that the convective boundary condition results in a temperature slip at the wall. RamReddy et al. [16] studied an incompressible laminar free convective micropolar fluid flow along a permeable vertical plate with the convective boundary condition and obtained a numerical solution using spectral quasilinearization method. Postelnicu [17] examined the effects of Soret and Dufour on a two-dimensional steady stagnation point flow of a free convective Darcian fluid with suction/injection and the governing partial differential equations are reduced to nonlinear ordinary differential equations using similarity transformations and then solved by the Keller box method. Hayat et al. [18] investigated the three-dimensional flow of a chemically reacting viscous fluid over exponentially stretching sheet with Soret and Dufour effects and observed that the thermal boundary layer increases for higher Dufour number while concentration boundary layer was increased for higher Soret values. Chamkha and Rashad [19] have numerically investigated the unsteady heat and mass transfer of an electrically conducting mixed convective chemically reacting viscous fluid flow over a rotating vertical cone with thermal-diffusion and diffusion-thermo effects. The MHD mixed convective heat and mass transfer of couple stress fluid through a porous channel with periodic suction and injection in the presence of cross-diffusion effects was considered by Ojjela and Naresh Kumar [20] and the reduced governing equations are solved numerically by the method of quasilinearization. The systematic framework of rate type viscoelastic fluids was developed by Oldroyd [21]. The unsteady three-dimensional flow of time dependent UCM fluid over a stretching sheet was studied by Awais et al. [22] and the series solution was obtained by HAM. Choi et al. [23] considered the two-dimensional steady incompressible laminar suction flow of a UCM fluid in a porous channel with the combined effects of inertia and viscoelasticity and obtained both analytical and numerical solutions. Hayat et al. [24] found an analytical solution for an electrically conducting Maxwell fluid peristaltic flow in a porous space. Abbas et al. [25] studied an incompressible magnetohydrodynamic two-dimensional boundary layer flow of UCM fluid through a rectangular porous channel using the homotopy analysis method. Mukhopadhyay and Gorla [26] analyzed the two-dimensional MHD flow and mass transfer of UCM fluid over an unsteady stretching sheet with first-order constructive/destructive chemical reaction and obtained a numerical solution by shooting method. Hayat and Abbas [27] inspected the effects of Deborah’s number, Schmidt number, and chemical reaction parameter on velocity and concentration fields in the flow of chemically reacting UCM fluid between porous parallel plates. The unsteady MHD flow of Maxwell fluid over an impulsively stretching sheet was analyzed by Alizadeh-Pahlavan and Sadeghy [28] and the flow field equations are solved by the homotopy analysis method. Differential transform method is first used by Zhou [29] in solving circuit problems. The theory of two-dimensional differential transform method was proposed by Chen and Ho [30]. Ayaz [31] studied the solutions of differential equations using differential transform method in a three-dimensional platform. Recently, several researchers have obtained solutions for various problems in fluid mechanics by differential transform method due to its simplicity and high accuracy. Mosayebidorcheh et al. [32] have investigated the effect of mass transfer on an incompressible laminar upper convected Maxwell fluid in a porous channel with high permeability medium and the reduced flow field equations are solved by differential transform method. Sheikholeslami et al. [33] have considered an incompressible steady nanofluid flow and heat transfer between two horizontal parallel plates with thermophoresis and Brownian motion under the influence of uniform magnetic field and the governing equations are reduced into nondimensional ordinary differential equations using similarity transformation and then solved by differential transform method. Sheikholeslami and Ganji [34] have studied the effect of Brownian motion on a nanofluid in the presence of variable magnetic field and the solution is obtained by differential transform method. Hatami and Jing [35] have analyzed two different problems, the unsteady two-dimensional squeezing nanofluid flow and heat transfer between two parallel plates and a non-Newtonian free convective nanofluid flow, and heat transfer between two vertical flat plates and obtained a solution for both by differential transform method and the results are compared with Runge-Kutta fourth-order method. Esmail and Taha [36] have found a solution to the atmospheric dispersion equation by using DTM and the results are shown to be in good agreement with field data.In the present paper, the effects of chemical reaction and Hall and ion slip currents on an unsteady incompressible free convective slip flow of an upper convected Maxwell fluid through a porous channel with Soret and Dufour in the presence of the convective boundary condition are considered. The reduced governing nonlinear ordinary differential equations are solved by differential transform method (DTM). The effects of various fluid and geometric parameters on nondimensional velocity components, temperature distribution, and concentration are discussed in detail and shown in the form of graphs. The present results are compared with Bujurke et al. [37] for the Newtonian fluid and presented in the form of a table.


## 2. Formulation of the Problem

Consider an unsteady incompressible slip flow of an electrically conducting upper convected Maxwell fluid between two parallel plates which are separated by “h.” Assume V
1
e
i
ω
t and V
2
e
i
ω
t are periodic injection and suction at the lower and upper plates, respectively. Also the temperature at the upper plate is maintained with T
2
e
i
ω
t, while the lower plate is governed by convective boundary condition and the concentrations at the lower and upper plates are C
1
e
i
ω
t and C
2
e
i
ω
t, respectively, as shown in Figure 1.Figure 1
Schematic diagram of fluid flow between porous parallel plates.The equations governing the heat and mass transfer of a UCM fluid flow in the presence of strong magnetic field and in the absence of body forces are [22, 28](1)
∂
u
∂
x
+
∂
v
∂
y
=
0
,
(2)
ρ
∂
u
∂
t
+
u
∂
u
∂
x
+
v
∂
u
∂
y
+
β
u
2
∂
2
u
∂
x
2
+
v
2
∂
2
u
∂
y
2
+
2
u
v
∂
2
u
∂
x
∂
y
=
-
∂
P
∂
x
+
μ
∇
2
u
-
σ
B
0
2
1
+
β
i
β
e
u
-
β
e
v
1
+
β
i
β
e
2
+
β
e
2
+
ρ
g
β
T
T
-
T
1
e
i
ω
t
+
ρ
g
β
C
C
-
C
1
e
i
ω
t
,
(3)
ρ
∂
v
∂
t
+
u
∂
v
∂
x
+
v
∂
v
∂
y
+
β
u
2
∂
2
v
∂
x
2
+
v
2
∂
2
v
∂
y
2
+
2
u
v
∂
2
v
∂
x
∂
y
=
-
∂
P
∂
y
+
μ
∇
2
v
-
σ
B
0
2
1
+
β
i
β
e
v
+
β
e
u
1
+
β
i
β
e
2
+
β
e
2
,
(4)
ρ
c
∂
T
∂
t
+
u
∂
T
∂
x
+
v
∂
T
∂
y
=
k
∂
2
T
∂
x
2
+
∂
2
T
∂
y
2
+
μ
2
∂
u
∂
x
2
+
2
∂
v
∂
y
2
+
∂
u
∂
y
+
∂
v
∂
x
2
+
σ
B
0
2
u
2
+
v
2
1
+
β
i
β
e
2
+
β
e
2
+
ρ
D
1
k
T
c
s
∂
2
C
∂
x
2
+
∂
2
C
∂
y
2
,
(5)
ρ
c
∂
C
∂
t
+
u
∂
C
∂
x
+
v
∂
C
∂
y
=
D
1
∂
2
C
∂
x
2
+
∂
2
C
∂
y
2
+
D
1
k
T
T
m
∂
2
T
∂
x
2
+
∂
2
T
∂
y
2
-
k
2
C
-
C
1
e
i
ω
t
.The associated boundary conditions are as in (6)(6)
u
x
,
λ
,
t
=
0
,
v
x
,
λ
,
t
=
V
1
e
i
ω
t
,
-
k
∂
T
∂
y
=
h
1
T
-
T
1
e
i
ω
t
,
C
x
,
λ
,
t
=
C
1
e
i
ω
t
a
t
y
=
0
,
u
x
,
λ
,
t
=
-
k
1
σ
1
∂
u
∂
y
,
v
x
,
λ
,
t
=
V
2
e
i
ω
t
,
T
x
,
λ
,
t
=
T
2
e
i
ω
t
,
C
x
,
λ
,
t
=
C
2
e
i
ω
t
a
t
y
=
h
.The governing equations are reduced to nondimensional equations using the following similarity transformations [20]:(7)
u
x
,
λ
,
t
=
U
0
a
-
V
2
x
h
f
′
λ
e
i
ω
t
,
v
x
,
λ
,
t
=
V
2
f
λ
e
i
ω
t
,
T
x
,
λ
,
t
=
T
1
+
μ
v
2
ρ
h
c
ϕ
1
λ
+
U
0
a
V
2
-
x
h
2
ϕ
2
λ
e
i
ω
t
,
C
x
,
λ
,
t
=
C
1
+
n
˙
A
h
υ
g
1
λ
+
U
0
a
V
2
-
x
h
2
g
2
λ
e
i
ω
t
,where f
(
λ
), ϕ
1
(
λ
), ϕ
2
(
λ
), g
1
λ, and g
2
λ are unknown functions to be determined.Substituting (7) in (2), (3), (4), and (5), we have(8)
f
i
v
1
-
Wi
⁡
Re
⁡
f
2
cos
⁡
ϕ
=
-
R
e
cos
⁡
ϕ
f
′
f
′
′
-
f
f
′
′
′
-
2
W
i
Re
⁡
f
′
2
f
′
′
+
f
f
′
′
2
cos
⁡
2
ϕ
+
H
a
2
1
+
β
i
β
e
1
+
β
i
β
e
2
+
β
e
2
f
′
′
-
E
c
G
r
ξ
ϕ
1
′
+
ξ
2
ϕ
2
′
-
S
h
G
c
ξ
g
1
′
+
ξ
2
g
2
′
,
ϕ
1
′
′
=
-
Re
⁡
P
r
cos
⁡
ϕ
1
-
D
u
S
c
S
r
4
f
′
2
+
H
a
2
1
+
β
i
β
e
2
+
β
e
2
f
2
-
f
ϕ
1
′
-
D
u
K
r
1
-
D
u
S
c
S
r
g
1
-
Re
⁡
Sc
Du
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
f
g
1
′
-
2
ϕ
2
,
ϕ
2
′
′
=
-
Re
⁡
P
r
cos
⁡
ϕ
1
-
D
u
S
c
S
r
f
′
′
2
+
H
a
2
1
+
β
i
β
e
2
+
β
e
2
f
′
2
+
2
f
′
ϕ
2
-
f
ϕ
2
′
-
D
u
K
r
1
-
D
u
S
c
S
r
g
2
-
Re
⁡
Sc
Du
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
f
g
2
′
-
2
f
′
g
2
,
g
1
′
′
=
-
2
g
2
+
K
r
1
-
D
u
S
c
S
r
g
1
+
Re
⁡
Sc
Du
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
f
g
1
′
+
Re
⁡
Sc
Pr
Sr
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
4
f
′
2
+
H
a
2
1
+
β
i
β
e
2
+
β
e
2
f
2
-
f
ϕ
1
′
,
g
2
′
′
=
-
K
r
1
-
D
u
S
c
S
r
g
1
+
Re
⁡
Sc
Pr
Sr
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
f
′
′
2
+
H
a
2
1
+
β
i
β
e
2
+
β
e
2
f
′
2
+
2
f
′
ϕ
2
-
f
ϕ
2
′
+
Re
⁡
Sc
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
f
g
2
′
-
2
f
′
g
2
,where prime denotes the differentiation with respect to λ and the nondimensional boundary conditions in terms of f, ϕ
1, ϕ
2, g
1, and g
2 are(9)
f
0
=
1
-
a
,
f
′
0
=
0
,
g
1
0
=
0
,
g
2
0
=
0
,
ϕ
1
′
0
=
-
γ
ϕ
1
0
,
ϕ
2
′
0
=
-
γ
ϕ
2
0
,
f
1
=
1
,
f
′
′
1
=
-
f
′
1
S
l
,
ϕ
1
1
=
1
E
c
,
ϕ
2
1
=
0
,
g
1
1
=
1
S
h
,
g
2
1
=
0
.


## 3. Solution of the Problem

The nonlinear equations (8) along with the boundary conditions (9) are transformed into recurrence equations (11)–(15) and initial values (16) using the transformations in (10)(10)
f
⟶
F
1
k
,
ϕ
1
⟶
F
2
k
,
ϕ
2
⟶
F
3
k
,
g
1
⟶
F
4
k
,
g
2
⟶
F
5
k
,
f
ϕ
1
⟶
∑
l
=
0
k
F
1
l
F
2
k
-
l
,
f
ϕ
1
ϕ
2
⟶
∑
s
=
0
k
∑
m
=
0
k
-
s
F
1
s
F
2
m
F
3
k
-
s
-
m
,
F
1
k
=
1
k
!
d
k
f
d
x
k
,
f
x
=
∑
j
=
0
k
F
1
j
x
j
,
d
n
f
d
x
n
⟶
k
+
1
k
!
F
1
k
+
n
,
F
2
k
=
1
k
!
d
k
ϕ
1
d
x
k
,
ϕ
1
x
=
∑
j
=
0
k
F
2
j
x
j
,
F
3
k
=
1
k
!
d
k
ϕ
2
d
x
k
,
ϕ
2
x
=
∑
j
=
0
k
F
3
j
x
j
,
F
4
k
=
1
k
!
d
k
g
1
d
x
k
,
g
1
x
=
∑
j
=
0
k
F
4
j
x
j
,
F
5
k
=
1
k
!
d
k
g
2
d
x
k
,
g
2
x
=
∑
j
=
0
k
F
5
j
x
j
,
(11)
F
1
k
+
4
=
R
e
cos
⁡
ϕ
∑
l
=
0
k
l
+
3
l
+
2
l
+
1
F
1
l
+
3
F
1
k
-
l
-
l
+
2
l
+
1
F
1
l
+
2
k
-
l
+
1
F
1
k
-
l
+
1
+
Re
⁡
W
i
cos
⁡
2
ϕ
∑
s
=
0
k
-
1
∑
m
=
0
k
-
s
s
+
4
s
+
3
s
+
2
s
+
1
F
1
m
F
1
s
+
4
F
1
k
-
m
-
s
-
2
Re
⁡
W
i
cos
⁡
2
ϕ
∑
s
=
0
k
∑
m
=
0
k
-
s
m
+
1
s
+
1
F
1
m
+
1
F
1
s
+
1
F
1
k
-
m
-
s
+
m
+
2
m
+
1
s
+
2
s
+
1
F
1
m
+
2
F
1
s
+
2
F
1
k
-
s
-
m
+
H
a
2
1
+
β
i
β
e
1
+
β
i
β
e
2
+
β
e
2
k
+
2
k
+
1
F
1
k
+
2
-
E
c
G
r
ξ
k
+
1
F
2
k
+
1
+
ξ
2
F
3
k
+
1
-
S
h
G
c
ξ
k
+
1
F
4
k
+
1
+
ξ
2
F
5
k
+
1
1
-
F
1
0
2
Re
⁡
W
i
cos
⁡
2
ϕ
k
+
4
k
+
3
k
+
2
k
+
1
-
1
,
(12)
F
2
k
+
2
=
-
Re
⁡
P
r
cos
⁡
ϕ
1
-
D
u
S
c
S
r
∑
l
=
0
k
4
l
+
1
F
1
l
+
1
k
-
l
+
1
F
1
k
-
l
+
1
-
l
+
1
F
2
l
+
1
F
1
k
-
l
-
Re
⁡
P
r
cos
⁡
ϕ
1
-
D
u
S
c
S
r
∑
l
=
0
k
H
a
2
1
+
β
i
β
e
2
+
β
e
2
F
1
l
F
1
k
-
l
-
D
u
K
r
1
-
D
u
S
c
S
r
F
4
k
-
Re
⁡
Sc
Du
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
∑
l
=
0
k
l
+
1
F
4
l
+
1
F
1
k
-
l
-
2
F
3
k
k
+
2
k
+
1
-
1
,
(13)
F
3
k
+
2
=
-
Re
⁡
P
r
cos
⁡
ϕ
1
-
D
u
S
c
S
r
∑
l
=
0
k
l
+
2
l
+
1
F
1
l
+
2
k
-
l
+
2
k
-
l
+
1
F
1
k
-
l
+
2
+
∑
l
=
0
k
H
a
2
1
+
β
i
β
e
2
+
β
e
2
l
+
1
F
1
l
+
1
k
-
l
+
1
F
1
k
-
l
+
1
+
∑
l
=
0
k
2
l
+
1
F
1
l
+
1
F
3
k
-
l
-
l
+
1
F
3
l
+
1
F
1
k
-
l
-
D
u
K
r
1
-
D
u
S
c
S
r
F
5
k
-
Re
⁡
Sc
Du
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
∑
l
=
0
k
l
+
1
F
5
l
+
1
F
1
k
-
l
-
2
l
+
1
F
1
l
+
1
F
5
k
-
l
k
+
2
k
+
1
-
1
,
(14)
F
4
k
+
2
=
-
2
F
5
k
-
K
r
1
-
D
u
S
c
S
r
F
4
k
-
Re
⁡
Sc
Du
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
∑
l
=
0
k
l
+
1
F
4
l
+
1
F
1
k
-
l
+
Re
⁡
Sc
Pr
S
r
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
∑
l
=
0
k
4
l
+
1
F
1
l
+
1
k
-
l
+
1
F
1
k
-
l
+
1
+
H
a
2
1
+
β
i
β
e
2
+
β
e
2
F
1
l
F
1
k
-
l
-
l
+
1
F
2
l
+
1
F
1
k
-
l
k
+
2
k
+
1
-
1
,
(15)
F
5
k
+
2
=
Re
⁡
S
c
cos
⁡
ϕ
1
-
D
u
S
c
S
r
∑
l
=
0
k
l
+
1
F
5
l
+
1
F
1
k
-
l
-
2
l
+
1
F
1
l
+
1
F
5
k
-
l
+
Re
⁡
Sc
Du
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
∑
l
=
0
k
l
+
2
l
+
1
F
1
l
+
2
k
-
l
+
2
k
-
l
+
1
F
1
k
-
l
+
2
+
Re
⁡
S
c
Pr
Sr
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
∑
l
=
0
k
H
a
2
1
+
β
i
β
e
2
+
β
e
2
l
+
1
F
1
l
+
1
k
-
l
+
1
F
1
k
-
l
+
1
+
Re
⁡
Sc
Pr
Sr
⁡
cos
⁡
ϕ
1
-
D
u
S
c
S
r
∑
l
=
0
k
2
l
+
1
F
1
l
+
1
F
3
k
-
l
-
l
+
1
F
3
l
+
1
F
1
k
-
l
-
K
r
1
-
D
u
S
c
S
r
F
5
k
k
+
2
k
+
1
-
1
.The transformed equations of the boundary conditions are (16)
F
1
0
=
1
-
a
,
F
1
1
=
0
,
F
2
0
=
n
3
,
F
3
0
=
0
,
F
4
0
=
0
,
F
5
0
=
0
,
F
1
2
=
n
1
,
F
1
3
=
n
2
,
F
2
1
=
γ
E
c
,
F
3
1
=
n
4
,
F
4
1
=
n
5
,
F
5
0
=
n
6
.The values in (16) are used in the recurrence relations (11)–(15) and n
1, n
2, n
3, n
4, n
5, and n
6 are calculated such that they satisfy the boundary conditions at λ
=
1; then we obtain nondimensional velocity components, temperature distribution, and concentration for various parameters by substituting these constants into the respective Taylor series.


## 4. Results and Discussion

In order to understand the flow characteristics in a better way, the numerical results for the nondimensional velocity components, temperature, and concentration with respect to various fluid and geometric parameters such as Weissenberg number Wi, slip parameter Sl, Biot numberγ, Soret number Sr, Dufour number Du, Hall parameter β
e, ion slip parameter β
i, and Prandtl number Pr are calculated in the domain [
0,1
] and presented in the form of graphs. A comparative study for the case of Newtonian fluid has been studied and the skin friction values using differential transform method and long series with polynomial coefficients as reported by Bujurke et al. [37] are presented in Table 1.Table 1
Comparison of the present nondimensional skin friction values with Bujurke et al. [37] for Newtonian case.
Serial number
a
Re
f
′
′
(
0
)
Bujurke et al. [37]
Present
1
0.04106
13.8833
0.0953
0.0954
2
0.15420
31.3577
0.3190
0.3395
3
0.00790
41.0960
0.0166
0.0176
4
0.16480
52.9262
0.3300
0.3593
5
0.01310
67.3621
0.0261
0.0298The effect of Weissenberg number Wi on velocity components, temperature distribution, and concentration is presented in Figure2. It is observed that as Wi increases the radial velocity and concentration also increase towards the upper plate and the velocity in X-direction increases towards the center of the plates and then decreases. However, the temperature decreases near the lower plate and then increases towards the upper plate. It is a known fact that for low Reynolds number flows (
0
<
R
e
<
5
) the viscoelastic effect dominates the inertial effect.Figure 2
Effect of Wi on (a) axial velocity, (b) radial velocity, (c) temperature, and (d) concentration ata
=
0.2; R
e
=
3; ϕ
=
0.4; H
a
=
4; β
i
=
8; β
e
=
0.2; G
r
=
10; G
c
=
10; γ
=
0.2; D
u
=
0.2; S
r
=
0.2; S
c
=
0.22; K
r
=
0.2; E
c
=
1; S
h
=
1; P
r
=
0.71; S
l
=
0.2.
(a)
(b)
(c)
(d)Figure3 gives the effect of slip parameter Sl on velocity components, temperature, and concentration. It is evident that as Sl increases, the velocity components and concentration are following the similar trend of Wi whereas the temperature is decreasing. It is due to the fact that the specific permeability decreases the velocity components of the fluid.Figure 3
Effect of slip parameter on (a) axial velocity, (b) radial velocity, (c) temperature, and (d) concentration ata
=
0.2; R
e
=
3; W
i
=
7; ϕ
=
0.4; H
a
=
4; β
i
=
8; β
e
=
0.2; G
r
=
10; G
c
=
10; γ
=
0.2; D
u
=
0.2; S
r
=
0.2; S
c
=
0.22; K
r
=
0.2; E
c
=
1; S
h
=
1; P
r
=
0.71.
(a)
(b)
(c)
(d)The effect of Biot numberγ on velocity components, temperature, and concentration is displayed in Figure 4. It is noticed that the radial velocity, temperature, and concentration profiles are decreasing with increasing γ and the axial velocity is decreasing towards the center of the plates and then increases. Since γ is inversely proportional to the thermal conductivity of the fluid, the heat transfer between the plates is decreased.Figure 4
Effect of Biot number on (a) axial velocity, (b) radial velocity, (c) temperature, and (d) concentration ata
=
0.2; R
e
=
3; W
i
=
1; ϕ
=
0.4; H
a
=
10; β
i
=
2; β
e
=
0.4; G
r
=
10; G
c
=
10; D
u
=
0.2; S
r
=
0.2; S
c
=
0.22; K
r
=
0.2; E
c
=
1; S
h
=
1; P
r
=
3; S
l
=
0.2.
(a)
(b)
(c)
(d)The effect of Soret number Sr on velocity components, temperature distribution, and concentration is studied in Figure5. As Sr increases, the axial velocity decreases up to the center of the channel and then increases towards the upper plate, whereas the radial velocity and concentration of the fluid decrease towards the upper plate. However, the temperature distribution increases initially and then decreases. This is because of the fact that the mean heat transfer between the plates decreases with the increasing of thermal diffusion.Figure 5
Effect of Sr on (a) axial velocity, (b) radial velocity, (c) temperature, and (d) concentration ata
=
0.2; R
e
=
3; W
i
=
7; ϕ
=
0.4; H
a
=
4; β
i
=
8; β
e
=
0.2; G
r
=
10; G
c
=
10; γ
=
0.2; D
u
=
0.2; S
c
=
0.22; K
r
=
0.2; E
c
=
1; S
h
=
1; P
r
=
0.71; S
l
=
0.2.
(a)
(b)
(c)
(d)The effect of Dufour number Du on the velocity components, temperature, and concentration is shown in Figure6. It is inferred that as Du increases, the radial velocity and concentration also increase whereas the temperature decreases. However the axial velocity increases towards the center of the channel and then decreases towards the upper plate. This is due to the fact that the concentration of the fluid is enhanced with diffusion in the temperature.Figure 6
Effect of Du on (a) axial velocity, (b) radial velocity, (c) temperature, and (d) concentration ata
=
0.2; R
e
=
3; W
i
=
7; ϕ
=
0.4; H
a
=
4; β
i
=
8; β
e
=
0.2; G
r
=
10; G
c
=
10; γ
=
0.2; S
r
=
0.2; S
c
=
0.22; K
r
=
0.2; E
c
=
1; S
h
=
1; P
r
=
0.71; S
l
=
0.2.
(a)
(b)
(c)
(d)The effect of Hall parameterβ
e on velocity components, temperature distribution, and concentration is shown in Figure 7. β
e increases the radial velocity and concentration also increases towards the upper plate, whereas the axial velocity increases towards the middle of the channel and then decreases. However, the temperature decreases near the lower plate and then increases towards the upper plate. This is because the effective conductivity decreases the damping force on the flow field.Figure 7
Effect ofβ
e on (a) axial velocity, (b) radial velocity, (c) temperature, and (d) concentration at a
=
0.2; R
e
=
3; W
i
=
1; ϕ
=
0.4; H
a
=
4; β
i
=
0.2; G
r
=
10; G
c
=
10; γ
=
0.2; Du = 0.2; S
r
=
0.2; S
c
=
0.22; K
r
=
0.2; E
c
=
1; S
h
=
1; P
r
=
0.71; S
l
=
0.2.
(a)
(b)
(c)
(d)Figure8 represents the effect of ion slip parameter β
i on velocity components, temperature, and concentration. It is observed that as β
i increases, the velocity components, temperature, and concentration are following the opposite trend of β
e. This is because the increase in the ion slip scales down the effect of Lorentz force [38].Figure 8
Effect ofβ
i on (a) axial velocity, (b) radial velocity, (c) temperature, and (d) concentration at a
=
0.2; R
e
=
3; W
i
=
1; ϕ
=
0.4; H
a
=
10; β
e
=
0.4; G
r
=
10; Gc = 10; γ
=
0.2; Du = 0.2; S
r
=
0.2; S
c
=
0.22; K
r
=
0.2; E
c
=
1; S
h
=
1; P
r
=
0.71; S
l
=
0.2.
(a)
(b)
(c)
(d)The effect of Prandtl number Pr on velocity components, temperature, and concentration is presented in Figure9. It is analyzed that as Pr increases the axial velocity decreases towards the center of the plates and then increases, whereas the radial velocity and concentration are decreasing from the lower plate to the upper plate. However, the temperature increases near the lower plate and then decreases towards the upper plate. Physically, if Pr increases the thermal diffusivity decreases and this leads to the reduction in the heat transfer ability at the thermal boundary layer.Figure 9
Effect of Pr on (a) axial velocity, (b) radial velocity, (c) temperature, and (d) concentration ata
=
0.2; R
e
=
3; W
i
=
1; ϕ
=
0.4; H
a
=
10; β
i
=
2; β
e
=
0.4; G
r
=
10; G
c
=
10; γ
=
0.2; Du = 0.2; S
r
=
0.2; S
c
=
0.22; K
r
=
0.2; E
c
=
1; S
h
=
1; S
l
=
0.2.
(a)
(b)
(c)
(d)


## 5. Conclusions

The effects of thermal-diffusion and diffusion-thermo on an unsteady free convective chemically reacting upper convected Maxwell fluid flow between two porous parallel plates with Hall and ion slip currents under the influence of slip and convective boundary conditions are considered. The flow field equations are reduced to nonlinear ordinary differential equations by using similarity transformations and then solved by differential transform method. The results are obtained for the nondimensional velocity components, temperature, and concentration distributions with respect to various fluid and geometric parameters such as Wi, Sl,γ, Sr, Du, β
e, β
i, and Pr and we observe the following:(i)
The concentration of the fluid is enhanced with slip parameter, whereas the temperature of the fluid is decreased with increasing of Biot number.
(ii)
Soret and Dufour parameters exhibit opposite trend on concentration distribution.
(iii)
The Weissenberg number and Hall parameter show similar effects on velocity components.
(iv)
The Hall and ion slip parameters exhibit reverse trend on temperature and concentration distributions.
(v)
The present results are compared with the previously published work [37] and are found to be in good agreement.These results have possible applications in science and technology such as fermentation, biorheology, geophysics, polymer industry, composite processing, food processing, and petroleum industries.



---

*Source: 9160956-2016-06-29.xml*