# Analysis of Peristaltic Motion of a Nanofluid with Wall Shear Stress, Microrotation, and Thermal Radiation Effects

Authors: C. Dhanapal; J. Kamalakkannan; J. Prakash; M. Kothandapani
Journal: Applied Bionics and Biomechanics
 (2016)

---

## Abstract

This paper analyzes the peristaltic flow of an incompressible micropolar nanofluid in a tapered asymmetric channel in the presence of thermal radiation and heat sources parameters. The rotation of the nanoparticles is incorporated in the flow model. The equations governing the nanofluid flow are modeled and exact solutions are managed under long wavelength and flow Reynolds number and long wavelength approximations. Explicit expressions of axial velocity, stream function, microrotation, nanoparticle temperature, and concentration have been derived. The phenomena of shear stress and trapping have also been discussed. Finally, the influences of various parameters of interest on flow variables have been discussed numerically and explained graphically. Besides, the results obtained in this paper will be helpful to those who are working on the development of various realms like fluid mechanics, the rotation, Brownian motion, thermophoresis, coupling number, micropolar parameter, and the nondimensional geometry parameters.


---

## Body

## 1. Introduction

Peristaltic pumping is one of the keystones for the development of science and engineering research in modern years. Peristalsis also plays an indispensable role in transporting physiological fluids inside living bodies, and many biomechanical and engineering devices have been designed on the basis of the principle of peristaltic pumping to transport fluids without internal moving parts. The problem of the mechanism of peristaltic transport has attracted the attention of many investigators since the first exploration of Latham [1]. A number of analytical, numerical, and experimental studies on peristaltic motion of different fluids have been described under various conditions with reference to physiological and mechanical environment [2–8].Micropolar fluids have been a subject of great interest to research workers and a number of research papers have been published on this flow model. Physically, micropolar fluids may represent fluids consisting of rigid, randomly oriented (or spherical) particles suspended in a viscous medium, where the deformation of fluid particles is ignored. This constitutes a substantial generalization of the Navier-Stokes model and opens a new field of potential applications including a large number of complex fluids. Animal bloods and liquid crystals (with dumbbell type molecules) are few examples of micropolar fluids. Local conservation laws of mass, linear, and angular momentum and the energy for polar fluids were received by Grad [9] by using the method of statistical thermodynamics. Eringen [10] proposed the theory of micropolar fluids in which the microscopic effect arises from local structure and fluid elements of micromotion are taken into account. Later, Eringen [11] generalized the micropolar fluids theory to include thermal effects. Using quasi-linearization finite difference technique, an impact of temperature dependent heat sources and frictional heating on the fully developed free convection micropolar fluid flow between two porous parallel plates was analyzed by Agarwal and Dhanapal [12]. Devi and Devanathan [13] premeditated the peristaltic motion of a micropolar fluid in a cylindrical tube with sinusoidal waves of small amplitude travelling down in its flexible wall for the case of low Reynolds number. Srinivasacharya et al. [14] recently examined the peristaltic transport of a micropolar fluid in a circular tube using low Reynolds number and long wavelength assumptions.Nowadays, there is a continuous focus of the researchers in the flow analysis of nanofluids because of its large number of applications in biomedical and industrial engineering. Choi [15] was the first who initiated this nanofluid technology. A detailed analysis of nanofluids was discussed by Buongiorno [16]. Sheikholeslami et al. [17] studied the natural convection in a concentric annulus between a cold outer square and heated inner circular cylinders in the presence of static radial magnetic field. After initiating a study of nanofluids flow under the effect of peristalsis by Akbar and Nadeem [18], Akbar et al. [19] discussed the slip effects on the peristaltic transport of nanofluid in an asymmetric channel. Recently, Mustafa et al. [20] examined the influence of wall properties on the peristaltic flow of a nanofluid. Mixed convection peristaltic flows of magnetohydrodynamic (MHD) nanofluids were analyzed by Hayat et al. [21]. The effects of wall properties on the peristaltic flow of an incompressible pseudoplastic fluid in a curved channel were investigated by Hina et al. [22]. Hayat et al. [23] studied the peristaltic transport of viscous nanofluid in an asymmetric channel. The channel walls satisfy the convective conditions and also effects of Brownian motion and thermophoresis have also been taken into account. The influence of nanofluid characteristic on peristaltic heat transfer in a two-dimensional axisymmetric channel was discussed analytically by Tripathi and Bég [24]. Moreover, the tremendous applications of nanofluids and the interaction of nanoparticles in peristaltic flows have obtained attentions of many researchers [25, 26].In the recent years, it is well known by physiologists [27, 28] that the intrauterine fluid flow due to myometrial contractions displays peristalsis and myometrial contractions may occur in both symmetric and asymmetric directions and also noted that blood behaves like as a non-Newtonian fluid in microcirculation [10–12]. Motivated from the above analysis and the importance of peristaltic flows, the purpose of the present paper is to investigate the effects of thermal radiation and heat source/sink on the peristaltic flow of micropolar nanofluids in the tapered asymmetric channel. Therefore, such a consideration of peristaltic transport may be used to evaluate intrauterine fluid flow in a nonpregnant uterus [29]. To the best of the author’s knowledge, no attempt is available in the literature which deals with the peristalsis flow of micropolar nanofluid in the tapered asymmetric channel. The present analysis of peristaltic flow is confined to large wavelength and low Reynolds number assumptions. Explicit solutions are developed for axial velocity, axial pressure gradient, stream function, microrotation of the nanofluids, nanofluid temperature, and nanoparticle concentration. The numerical discussion of the pressure rise, shear stresses, and trapping are also obtained and the results are discussed through graphs.


## 2. Mathematical Formulation

Let us consider the motion of peristaltic transport of an incompressible micropolar nanofluid through a tapered channel induced by sinusoidal wave trains propagating with constant speed but with different amplitudes and phases; see Figure1. The governing equations of motion for the present investigation are [13, 16, 30] (1)
∂
U
∂
X
+
∂
V
∂
Y
=
0
,
ρ
f
∂
U
∂
t
′
+
U
∂
U
∂
X
+
V
∂
U
∂
Y
=
-
∂
P
∂
X
+
μ
v
+
k
v
∇
2
U
-
k
v
∂
W
∂
Y
+
1
-
C
0
ρ
f
g
α
T
-
T
0
+
ρ
p
-
ρ
f
g
β
′
C
-
C
0
,
ρ
f
∂
V
∂
t
′
+
U
∂
V
∂
X
+
V
∂
V
∂
Y
=
-
∂
P
∂
Y
+
μ
v
+
k
v
∇
2
V
-
k
v
∂
W
∂
X
,
ρ
f
J
∂
W
∂
t
′
+
U
∂
W
∂
X
+
V
∂
W
∂
Y
=
γ
v
∇
2
W
+
k
v
∂
V
∂
X
-
∂
U
∂
Y
-
2
k
v
W
,
ρ
c
′
f
∂
T
∂
t
′
+
U
∂
T
∂
X
+
V
∂
T
∂
Y
=
κ
∂
2
T
∂
X
2
+
∂
2
T
∂
Y
2
-
∂
q
r
∂
X
+
∂
q
r
∂
Y
+
Q
0
T
-
T
0
+
ρ
c
′
p
D
B
∂
C
∂
X
∂
T
∂
X
+
∂
C
∂
Y
∂
T
∂
Y
+
D
T
ρ
c
p
T
m
∂
T
∂
X
2
+
∂
T
∂
Y
2
,
∂
C
∂
t
′
+
U
∂
C
∂
X
+
V
∂
C
∂
Y
=
D
B
∂
2
C
∂
X
2
+
∂
2
C
∂
Y
2
+
D
T
T
m
∂
2
T
∂
X
2
+
∂
2
T
∂
Y
2
,where U, V are the components of velocity along X and Y directions, respectively, t
′ is the dimensional time, the volumetric volume expansion coefficient is c
′, ρ
f is the density of the fluid, ρ
p is the density of the particle, g is the acceleration due to gravity, P is the pressure, μ
v, γ
v, and k
v are the martial parameters [9–12], T is the temperature, C is the nanoparticle concentration, α is the thermal expansion coefficient, ∂
/
∂
t
′ represents the material time derivative, β
′ is the coefficient of expansion with concentration, W is the microrotation of the nanofluid, j is the microgyration parameter, T
m is the fluid mean temperature, τ
=
(
ρ
c
′
)
p
/
(
ρ
c
′
)
f is the ratio of the effective heat capacity of nanoparticle material and heat capacity of the fluid with ρ being the density, κ is the thermal conductivity of the nanofluids, D
B is the Brownian diffusion coefficient, D
T is the thermophoretic diffusion coefficient, Q
0 is the constant heat addition/absorption, and the radioactive heat flux is q
r.Figure 1
Geometry of the generalized channel (tapered asymmetric channel) with peristaltic wave motion of wall.Hence, for the Rosseland approximation for thermal radiation, we have [31, 32](2)
q
r
=
-
4
σ
∗
3
k
∗
∂
T
4
∂
Y
,where σ
∗ and k
∗ are the Stefan-Boltzmann constant and the mean absorption coefficient.LetY
=
H
1 and Y
=
H
2 be, respectively, the left and right wall boundaries of the tapered asymmetric channel. Heat and mass transfer along with nanoparticle phenomena have been taken into account. The right wall of the channel is sustained at temperature T
1 and nanoparticle volume fraction C
1 while the left wall has temperature T
0 and nanoparticle volume fraction C
0. The geometry of the wall surface is defined as (3a)
H
1
X
,
t
′
=
-
d
-
m
′
X
-
a
1
s
i
n
2
π
λ
X
-
c
t
′
+
ϕ
,
(3b)
H
2
X
,
t
′
=
d
+
m
′
X
+
a
2
s
i
n
2
π
λ
X
-
c
t
′
,where a
1 and a
2are the amplitudes of left and right walls, respectively, λ is the wavelength, m
′ (m
′
≪
1) is the nonuniform parameter, the phase difference ϕ varies in the range 0
≤
ϕ
≤
π, ϕ
=
0corresponds to symmetric channel with waves out of the phase, and further a
1, a
2, d, and ϕ satisfy the condition for the divergent channel at the inlet of flow (4)
a
1
2
+
a
2
2
+
2
a
1
a
2
c
o
s
ϕ
≤
2
d
2
.The nondimensional parameters are as follows:(5)
x
=
X
λ
,
y
=
Y
d
,
t
=
c
t
′
λ
,
u
=
U
c
,
v
=
V
c
δ
,
δ
=
d
λ
,
h
1
=
H
1
d
,
h
2
=
H
2
d
,
Ω
=
d
W
c
,
j
=
J
d
2
m
=
m
′
λ
d
,
a
=
a
1
d
,
b
=
a
2
d
,
θ
=
T
-
T
0
T
1
-
T
0
,
σ
=
C
-
C
0
C
1
-
C
0
,
G
r
=
1
-
C
0
ρ
f
g
α
d
2
T
1
-
T
0
c
μ
v
,
β
=
Q
0
d
2
T
1
-
T
0
ν
c
p
,
P
r
=
μ
c
f
κ
,
N
b
=
τ
D
B
C
1
-
C
0
ν
,
N
t
=
τ
D
T
C
1
-
C
0
T
0
ν
,
R
n
=
16
σ
∗
T
0
3
3
k
∗
μ
v
c
f
,
B
r
=
ρ
p
-
ρ
f
g
β
′
d
2
C
1
-
C
0
c
μ
,
S
c
=
v
D
B
,
p
=
d
2
P
c
λ
μ
,
R
=
c
d
ρ
f
μ
v
.Using the above nondimensional quantities in (1)–(4), the resulting equations are(6)
R
δ
∂
∂
t
+
u
∂
∂
x
+
v
∂
∂
y
u
=
-
∂
p
∂
x
+
1
1
-
N
N
∂
Ω
∂
y
+
δ
2
∂
2
u
∂
x
2
+
∂
2
u
∂
y
2
+
G
r
θ
+
B
r
σ
,
R
δ
3
∂
∂
t
+
u
∂
∂
x
+
v
∂
∂
y
v
=
-
∂
p
∂
y
+
δ
2
1
-
N
-
N
∂
w
∂
x
+
δ
2
∂
2
v
∂
x
2
+
∂
2
v
∂
y
2
,
R
δ
j
1
-
N
N
∂
∂
t
+
u
∂
∂
x
+
v
∂
∂
y
Ω
=
-
2
Ω
+
δ
2
∂
v
∂
x
-
∂
u
∂
y
+
2
-
N
n
2
δ
2
∂
2
Ω
∂
x
2
+
∂
2
Ω
∂
y
2
,
R
δ
∂
θ
∂
t
+
u
∂
θ
∂
x
+
δ
v
∂
θ
∂
y
=
1
P
r
δ
2
∂
2
θ
∂
x
2
+
∂
2
θ
∂
y
2
+
R
n
∂
2
θ
∂
y
2
+
N
b
δ
2
∂
σ
∂
x
∂
θ
∂
x
+
∂
σ
∂
y
∂
θ
∂
y
+
β
θ
+
N
t
δ
2
∂
θ
∂
x
2
+
∂
θ
∂
y
2
,
R
δ
S
c
∂
σ
∂
t
+
u
∂
σ
∂
x
+
δ
v
∂
σ
∂
y
=
δ
2
∂
2
σ
∂
x
2
+
∂
2
σ
∂
y
2
+
N
t
N
b
δ
2
∂
2
θ
∂
x
2
+
∂
2
θ
∂
y
2
.in which N
=
k
v
/
(
μ
v
+
k
v
) is the coupling number (
0
≤
N
≤
1
) and n
2
=
d
2
k
v
(
2
μ
v
+
k
v
)
/
(
γ
v
(
μ
v
+
k
v
)
) is the micropolar parameter. We introduce the nondimensional variables such as p is dimensionless pressure, a and b are amplitudes of left and right walls, respectively, δ is wave number, m is the nonuniform parameter, R is the Reynolds number, ν is the nanofluid kinematic viscosity, Ω is the dimensionless microrotation, θ is the dimensionless temperature, σ is the dimensionless rescaled nanoparticle volume fraction, P
r is the Prandtl number, G
r is the local temperature Grashof number, B
r is the local nanoparticle Grashof number, Sc is the Schmidt number, N
b is the Brownian motion parameter, N
t is the thermophoresis parameter, and R
n is the radiation parameter as follows. The above equations can reduce to the classical Navier-Stokes equation when k
v
→
0.In several previous attempts [19–24], we employ the long wavelength and low Reynolds number approximations and thus (6) that (7)
∂
p
∂
x
=
1
1
-
N
∂
2
u
∂
y
2
+
N
1
-
N
∂
Ω
∂
y
+
G
r
θ
+
B
r
σ
,
(8)
∂
p
∂
y
=
0
,
(9)
-
2
Ω
-
∂
u
∂
y
+
2
-
N
n
2
∂
2
Ω
∂
y
2
=
0
,
(10)
1
+
R
n
P
r
P
r
∂
2
θ
∂
y
2
+
N
b
∂
σ
∂
y
∂
θ
∂
y
+
N
t
∂
θ
∂
y
2
+
β
θ
=
0
,
(11)
∂
2
σ
∂
y
2
+
N
t
N
b
∂
2
θ
∂
y
2
=
0
.The appropriate boundary conditions are (12a)
u
=
0
,
Ω
=
0
,
θ
=
0
,
σ
=
0
a
t
y
=
h
1
=
-
1
-
m
x
-
a
s
i
n
2
π
x
-
t
+
ϕ
,
(12b)
u
=
0
,
Ω
=
0
,
θ
=
1
,
σ
=
1
a
t
y
=
h
2
=
1
+
m
x
+
b
s
i
n
2
π
x
-
t
,which satisfy, at the inlet of channel, (13)
a
2
+
b
2
+
2
a
b
c
o
s
ϕ
≤
4
.


## 3. Exact Solution

By integration of (11) with respect to y and implementation in (10) and boundary conditions of (12a) and (12b), the nanoparticles temperature field is obtained as(14)
θ
=
s
i
n
h
h
1
θ
1
c
o
s
h
θ
1
y
-
c
o
s
h
h
1
θ
1
s
i
n
h
θ
2
y
s
i
n
h
h
1
θ
1
c
o
s
h
θ
1
h
2
-
c
o
s
h
h
1
θ
1
s
i
n
h
θ
2
h
2
.Substituting (14) into (11), moreover integrating (11) with respect to y and using proper boundary conditions of (12a) and (12b), the nanoparticle concentration field is received as (15)
σ
=
h
-
y
N
b
+
B
N
t
s
i
n
h
h
2
θ
2
+
A
N
t
c
o
s
h
h
2
θ
1
+
y
-
h
2
N
t
B
s
i
n
h
h
1
θ
2
+
A
c
o
s
h
h
1
θ
1
N
b
h
1
-
h
2
-
N
t
N
b
A
c
o
s
h
θ
1
y
+
B
s
i
n
h
θ
2
y
. Equation (7) can be written in the following form:(16)
∂
2
u
∂
y
2
=
∂
∂
y
1
-
N
∂
p
∂
x
y
-
N
ϖ
-
1
-
N
G
r
θ
-
1
-
N
B
r
σ
.Integration of the above equation yields(17)
∂
u
∂
y
=
1
-
N
∂
p
∂
x
y
-
N
Ω
-
1
-
N
G
r
θ
1
θ
2
θ
2
A
s
i
n
h
θ
1
y
+
θ
1
B
c
o
s
h
θ
2
y
-
1
-
N
B
r
2
c
y
+
D
y
2
2
-
N
t
N
b
θ
1
θ
2
A
θ
2
s
i
n
h
θ
1
y
+
B
θ
1
c
o
s
h
θ
2
y
+
G
x
.From (9) and (17), one can write(18)
∂
2
Ω
∂
y
2
-
n
2
Ω
=
n
2
1
-
N
2
-
N
∂
p
∂
x
y
+
A
4
s
i
n
h
θ
1
y
+
A
5
c
o
s
h
θ
2
y
+
A
6
y
+
A
7
y
2
+
m
2
G
x
2
-
N
.The general solution of the above equation can be written as (19)
Ω
=
E
c
o
s
h
n
y
+
F
s
i
n
h
n
y
-
1
-
N
2
-
N
∂
p
∂
x
y
+
A
4
s
i
n
h
θ
1
y
θ
1
2
-
n
2
+
A
5
c
o
s
h
θ
2
y
θ
2
2
-
n
2
-
A
6
y
n
2
-
A
7
n
2
y
2
+
2
n
4
-
G
x
2
-
N
.Making the above equation into (17), one obtains(20)
u
=
∂
p
∂
x
y
2
2
-
N
-
A
8
N
s
i
n
h
n
y
n
-
A
11
N
c
o
s
h
n
y
n
-
A
9
N
s
i
n
h
n
y
n
-
A
12
N
c
o
s
h
n
y
n
+
c
o
s
h
θ
1
y
1
-
N
B
r
N
t
A
N
b
θ
1
2
-
1
-
N
G
r
A
θ
1
2
-
N
A
4
θ
1
θ
1
2
-
n
2
+
N
A
7
3
n
2
-
1
-
N
B
r
D
6
y
3
+
s
i
n
h
θ
2
y
1
-
N
B
r
N
t
B
N
b
θ
2
2
-
1
-
N
G
r
B
θ
2
2
-
N
A
5
θ
2
θ
2
2
-
n
2
+
2
A
7
y
N
n
4
+
N
A
6
y
2
2
n
2
+
G
x
2
y
2
-
N
+
A
13
N
c
o
s
h
n
y
n
+
A
10
N
s
i
n
h
n
y
n
+
H
x
.The volume flux through each cross section in the wave frame is given by (21)
F
=
∫
h
1
h
2
u
d
y
.Using (21), we find that(22)
∂
p
∂
x
=
F
-
A
21
θ
1
sinh
⁡
θ
1
h
2
-
sinh
⁡
θ
1
h
1
-
A
22
θ
2
c
o
s
h
θ
2
h
2
-
c
o
s
h
θ
2
h
1
-
A
32
h
2
3
-
h
1
3
3
-
A
29
n
s
i
n
h
n
h
2
-
s
i
n
h
n
h
1
-
A
30
n
c
o
s
h
n
h
2
-
c
o
s
h
n
h
1
-
A
24
h
2
4
-
h
1
4
4
-
A
31
h
2
2
-
h
1
2
4
-
A
25
h
2
-
h
1
·
A
20
h
2
2
-
h
1
2
2
-
N
+
h
2
3
-
h
1
3
6
-
3
N
+
A
27
n
c
o
s
h
n
h
2
-
c
o
s
h
n
h
1
+
A
28
n
s
i
n
h
n
h
2
-
c
o
s
h
n
h
1
+
A
26
h
2
-
h
1
-
1
.The corresponding stream function from (20) is(23)
ψ
=
∂
p
∂
x
y
2
2
-
N
-
A
8
N
c
o
s
h
n
y
-
1
n
2
-
A
11
N
s
i
n
h
n
y
n
2
-
A
9
N
c
o
s
h
n
y
n
2
-
A
12
N
s
i
n
h
n
y
n
2
+
sinh
⁡
θ
1
y
θ
1
1
-
N
B
r
N
t
A
N
b
θ
1
2
-
1
-
N
G
r
A
θ
1
2
-
N
A
4
θ
1
θ
1
2
-
n
2
+
N
A
7
12
n
2
-
1
-
N
B
r
D
24
y
4
.The constant values appeared are listed out in Appendix.The nondimensional expression for the pressure rise per wavelengthΔ
p
λ is(24)
Δ
P
λ
=
∫
0
1
∫
0
1
∂
p
∂
x
y
=
0
d
x
d
t
.Interestingly, we note that the stress tensor in micropolar fluid is not symmetric in behavior. For that reason, the dimensionless form of the shear stress implicated in the present problem under consideration is given by [30] (25)
τ
x
y
=
∂
u
∂
y
-
N
1
-
N
Ω
,
(26)
τ
x
y
=
3
y
2
B
r
D
N
-
1
6
+
A
7
N
3
n
2
-
∂
p
∂
x
y
N
-
1
+
A
8
N
c
o
s
h
n
y
+
A
11
N
s
i
n
h
n
y
-
2
N
y
N
-
1
2
N
-
2
+
G
A
10
N
c
o
s
h
n
y
-
N
N
-
2
+
A
13
N
s
i
n
h
n
y
+
1
+
c
o
s
h
θ
2
y
B
G
r
N
-
1
θ
2
+
A
5
N
n
2
-
θ
2
2
-
B
B
r
N
t
N
-
11
N
b
θ
2
-
N
N
-
1
A
5
c
o
s
h
θ
2
y
n
2
-
θ
2
2
-
F
s
i
n
h
n
y
-
G
N
-
2
-
E
c
o
s
h
n
y
+
A
4
s
i
n
h
θ
1
y
n
2
-
θ
1
2
+
A
7
n
2
y
2
+
2
n
4
+
A
6
y
n
2
+
∂
p
/
∂
x
N
-
1
N
-
2
+
θ
1
s
i
n
h
θ
1
y
A
G
r
N
-
1
θ
1
2
+
A
4
N
θ
1
n
2
-
θ
1
2
-
A
B
r
N
t
N
-
1
N
b
θ
1
2
-
A
9
N
c
o
s
h
n
y
-
A
12
N
s
i
n
h
n
y
+
2
A
7
N
n
4
+
A
6
N
y
n
2
,
τ
y
x
=
1
1
-
N
∂
u
∂
y
+
N
1
-
N
Ω
,
τ
y
x
=
-
1
N
-
1
3
y
2
B
r
D
N
-
1
6
+
A
7
N
3
n
2
-
∂
p
∂
x
y
N
-
1
+
A
8
N
c
o
s
h
n
y
+
A
11
N
s
i
n
h
n
y
-
2
N
y
N
-
1
2
N
-
2
+
G
A
10
N
c
o
s
h
n
y
-
N
N
-
2
+
A
13
N
s
i
n
h
n
y
+
1
+
c
o
s
h
θ
2
y
B
G
r
N
-
1
θ
2
+
A
5
N
n
2
-
θ
2
2
-
B
B
r
N
t
N
-
11
N
b
θ
2
+
θ
1
s
i
n
h
θ
1
y
A
G
r
N
-
1
θ
1
2
+
A
4
N
θ
1
n
2
-
θ
1
2
-
A
B
r
N
t
N
-
1
N
b
θ
1
2
+
2
A
7
N
n
4
+
A
6
N
y
n
2
-
A
9
N
c
o
s
h
n
y
-
A
12
N
s
i
n
h
n
y
-
N
N
-
1
A
5
c
o
s
h
θ
2
y
n
2
-
θ
2
2
-
F
s
i
n
h
n
y
-
G
N
-
2
-
E
c
o
s
h
n
y
+
A
4
s
i
n
h
θ
1
y
n
2
-
θ
1
2
+
A
7
n
2
y
2
+
2
n
4
+
A
6
y
n
2
+
∂
p
/
∂
x
N
-
1
N
-
2
.The numerical computation for the shear stress τ
x
y is obtained at left wall of the channel and whose graphical representation is presented in the next section.


## 4. Numerical Results and Discussion

In view of the fact that the constant value of rate of volume flowF gives the pressure rise Δ
P
λ always negative and hence pumping action cannot be perceived. To discuss the found results quantitatively, we assume the form of the instantaneous volume rate of the flow F
(
x
,
t
) [30, 33–37],(27)
F
x
,
t
=
Θ
+
a
s
i
n
2
π
x
-
t
+
b
s
i
n
2
π
x
-
t
+
ϕin which F
=
Q
¯
/
c
d, Θ
=
q
¯
/
c
d, F
=
∫
h
1
h
2
u
d
y
=
ψ
(
h
2
)
-
ψ
(
h
1
), and Θ is the time-average of flow over one period of wave.In order to get nearby into the given substantial problem, we observe physical characteristics of average rise in pressure, axial velocity, microrotation velocity, shear stress, nanoparticles temperature, and concentration with respect of various values to the parameters appearing in the problem by sketching Figures2–20 with the constant values x
=
0.5 and t
=
0.2. It is found that, in the absence of nonuniform parameter, local temperature Grashof number, and local nanoparticle Grashof number, the present analysis reduces to approximate analytical solution of peristaltic flow of a Newtonian fluid in asymmetric channel [32]. Thus, the results introduced may be applied in the real life problems associated with nanoparticles movement in the gastrointestinal tract, intrauterine fluid motion induced by uterine contraction, and flow through small blood vessels and intrapleural membranes.Figure 2
Variation ofm on the velocity u with respect to y.Figure 3
Variation ofb on the velocity u with respect to y.Figure 4
Variation of Gr on the velocityu with respect to y.Figure 5
Variation ofN
b on the velocity u with respect to y.Figure 6
Variation ofN on the velocity u with respect to y.Figure 7
Variation ofn on the velocity u with respect to y.Figure 8
Temperature profileθ
(
y
) for β.Figure 9
Nanoparticle volume fractionσ
(
y
) for β.Figure 10
Temperature profileθ
(
y
) for N
b.Figure 11
Nanoparticle volume fractionσ
(
y
) for N
b.Figure 12
Temperature profileθ
(
y
) for R
n.Figure 13
Nanoparticle volume fractionσ
(
y
) for R
n.Figure 14
Microrotation velocity profileΩ
(
y
) for different values of N.Figure 15
Microrotation velocity profileΩ
(
y
) for different values of N
b.Figure 16
Microrotation velocity profileΩ
(
y
) for different values of m.Figure 17
Shear stress profileτ
x
y
(
x
) for different values of m.Figure 18
Shear stress profileτ
x
y
(
x
) for different values of n.Figure 19
Shear stress profileτ
x
y
(
x
) for different values of N.Figure 20
Streamlines form
=
0, N
=
0.01, n
=
0.1 (a), m
=
0.3, N
=
0.01, n
=
0.1 (b), m
=
0.3, N
→
1, n
=
0.1 (c), and m
=
0.3, N
→
1, n
=
0.9 (d). The other parameters chosen are a
=
0.2, b
=
0.3, m
=
0.3, Θ
=
1.5, ϕ
=
π
/
2, N
t
=
0.2, N
b
=
0.3, R
n
=
0.9, P
r
=
0.5, G
r
=
0.6,  B
r
=
0.2, β
=
0.2, and   t
=
0.2.
(a)
(b)
(c)
(d)


### 4.1. Flow Characteristics

Figures2–7 express that the variation of axial velocity u with respect to y for different values of the nonuniform parameter (
m
), amplitude of right wall (
b
), local temperature Grashof number (
G
r
), Brownian motion parameter (
N
b
), micropolar parameter (
n
), and coupling number (
N
). The dimensionless axial velocity profiles (
u
) satisfy the boundary conditions and are varied by a smooth curve for different (
m
) values in Figure 2. This nonuniform parameter effect has a tendency to slow down the motion of the fluid which results in decreasing the axial velocity profiles at the core of the channel. Figure 3 represents the variation of u with y for various values of the amplitude of right wall (
b
). It is found that an increase of b results in increase of the velocity of the fluid near the right part of the channel. The effects of G
r and N
b on the axial velocity pattern are shown in Figures 4 and 5. When G
r and N
b values are larger and u curves are larger too at the central part of the channel, so the velocity profiles become larger with right wall. This is due to the fact that buoyancy force gives rise to fluid flow. This force has a tendency to accelerate the motion of the fluid which results in increasing the axial velocity profiles at right wall; otherwise it gets decreased. The effect of coupling number N on axial velocity is depicted in Figure 6. It is observed that the flow reversal near the left wall of the channel increases with the increase of coupling number, while the reversal trend occurs in the vicinity of the right channel wall. But similar significant change is found for the variation of micropolar parameter n in Figure 7.


### 4.2. Heat Transfer and Nanoparticle Mass Transfer Distributions

The dimensionless nanoparticle temperature and concentration profilesθ
,
σ satisfy the boundary conditions and are varied by a smooth curve for different β, N
b, and R
n values. In the absence of amplitude of left wall a
=
0, heat source/sink β
=
0, Prandtl number P
r
=
1, thermal radiation parameter R
n
=
0, and nonuniform parameter m
=
0 our nanoparticle temperature and concentration distributions results are in close agreement with earlier works of Tripathi and Bég [24]. The effects of heat source/sink parameter β on the nanoparticle temperature and concentration distributions are displayed in Figures 8 and 9. That is, increase in the heat source strength amounts to increase in energy supply to the tapered asymmetric channel and opposite behavior is noticed form nanoparticle mass transfer. The effects of Brownian motion parameter on the nanoparticle temperature and concentration distributions are considered in Figures 10 and 11. When N
b is larger and σ curve is lower, so the mass transfer effect is higher for a larger N
b and opposite behavior found in heat transfer. In nanofluid system, the size of the nanoparticle generates Brownian motion which affects the heat and mass transfer properties. As the particle size scale approaches to the nanometer scale, the Brownian motion of particles and its diffusion effect on the liquid play a significant role in heat transfer. For prominent values of N
b, the Brownian diffusion effect is large compared to the thermal diffusion effect. Figure 12 shows the variation of thermal radiation parameter over the temperature field. Therefore higher values of thermal radiation parameter imply higher surface heat flux and so, it decreases the temperature within the tapered asymmetric channel, whereas opposite effects are observed in a nanoparticle mass transfer as shown in Figure 13.


### 4.3. Spin Velocity Distribution

The variation of microrotation velocity(
Ω
) with respect to y for various values of coupling number (
N
), Brownian motion (
N
b
), and nonuniform parameter (
m
) are shown in Figures 14–16. Parabolic microrotation velocity profile is observed for the present flow problem. The velocity is maximum or minimum near the center of channel. The cause of coupling number on the microrotation velocity is illustrated in Figure 14. It is seen that in the microrotation velocity increases by increasing the coupling number N (i.e., microrotation velocity for the micropolar nanofluid is wider than that of the region for Newtonian nanofluid). In the case of N
=
0, there is no appreciable difference between Newtonian nanofluids and micropolar nanofluids. The microrotation velocity for the Brownian motion parameter N
b is plotted in Figure 15. It is seen that, with the increase in the Brownian motion parameter, microrotation velocity profile decreases and maximum microrotation velocity occurs at N
b
→
0. Figure 16 depicts the microrotation velocity field for different values of nonuniform parameter. It is also viewed that the microrotation velocity for a divergent channel (m
>
0) is higher compared to its value for a uniform channel (m
=
0).


### 4.4. Shear Stress Distribution

It is well known that the stress tensor is not symmetric in micropolar nanofluid. In Figures17–19, we have plotted the shear stresses τ
x
y at the left wall for values of the nonuniform parameter, micropolar parameter, and coupling number. One can observe from Figure 17 that the wall shear stress decreases with an increase in the nonuniform parameter. From Figure 18, we notice that shear stress is in oscillation behavior, which may be due to the creation of contraction and expansion walls. It is indicated that the shear stress decreases with an increase in the micropolar parameter n, while it increases as the coupling number N increases in Figure 19. Hence, we observed that the shear stress for a Newtonian nanofluid is less than that for a micropolar nanofluid.


### 4.5. Trapping Phenomena

The effects of nonuniform parameter, coupling parameter, and micropolar parameter on the streamlines are shown in Figure20. Besides, a comparison between uniform channel and the nonuniform channel is made in Figures 20(a) and 20(b). We note that the size of the trapping bolus increases with increasing nonuniform parameter and symmetry nature is also noticed with respect to uniform symmetric channel. The effects of N on trapping are presented in Figures 20(b) and 20(c). This figure reveals that the size of lower bolus decreases with an increase in N. From Figures 20(c) and 20(d), it is clear that the trapped bolus increases in size as n increases.


## 4.1. Flow Characteristics

Figures2–7 express that the variation of axial velocity u with respect to y for different values of the nonuniform parameter (
m
), amplitude of right wall (
b
), local temperature Grashof number (
G
r
), Brownian motion parameter (
N
b
), micropolar parameter (
n
), and coupling number (
N
). The dimensionless axial velocity profiles (
u
) satisfy the boundary conditions and are varied by a smooth curve for different (
m
) values in Figure 2. This nonuniform parameter effect has a tendency to slow down the motion of the fluid which results in decreasing the axial velocity profiles at the core of the channel. Figure 3 represents the variation of u with y for various values of the amplitude of right wall (
b
). It is found that an increase of b results in increase of the velocity of the fluid near the right part of the channel. The effects of G
r and N
b on the axial velocity pattern are shown in Figures 4 and 5. When G
r and N
b values are larger and u curves are larger too at the central part of the channel, so the velocity profiles become larger with right wall. This is due to the fact that buoyancy force gives rise to fluid flow. This force has a tendency to accelerate the motion of the fluid which results in increasing the axial velocity profiles at right wall; otherwise it gets decreased. The effect of coupling number N on axial velocity is depicted in Figure 6. It is observed that the flow reversal near the left wall of the channel increases with the increase of coupling number, while the reversal trend occurs in the vicinity of the right channel wall. But similar significant change is found for the variation of micropolar parameter n in Figure 7.


## 4.2. Heat Transfer and Nanoparticle Mass Transfer Distributions

The dimensionless nanoparticle temperature and concentration profilesθ
,
σ satisfy the boundary conditions and are varied by a smooth curve for different β, N
b, and R
n values. In the absence of amplitude of left wall a
=
0, heat source/sink β
=
0, Prandtl number P
r
=
1, thermal radiation parameter R
n
=
0, and nonuniform parameter m
=
0 our nanoparticle temperature and concentration distributions results are in close agreement with earlier works of Tripathi and Bég [24]. The effects of heat source/sink parameter β on the nanoparticle temperature and concentration distributions are displayed in Figures 8 and 9. That is, increase in the heat source strength amounts to increase in energy supply to the tapered asymmetric channel and opposite behavior is noticed form nanoparticle mass transfer. The effects of Brownian motion parameter on the nanoparticle temperature and concentration distributions are considered in Figures 10 and 11. When N
b is larger and σ curve is lower, so the mass transfer effect is higher for a larger N
b and opposite behavior found in heat transfer. In nanofluid system, the size of the nanoparticle generates Brownian motion which affects the heat and mass transfer properties. As the particle size scale approaches to the nanometer scale, the Brownian motion of particles and its diffusion effect on the liquid play a significant role in heat transfer. For prominent values of N
b, the Brownian diffusion effect is large compared to the thermal diffusion effect. Figure 12 shows the variation of thermal radiation parameter over the temperature field. Therefore higher values of thermal radiation parameter imply higher surface heat flux and so, it decreases the temperature within the tapered asymmetric channel, whereas opposite effects are observed in a nanoparticle mass transfer as shown in Figure 13.


## 4.3. Spin Velocity Distribution

The variation of microrotation velocity(
Ω
) with respect to y for various values of coupling number (
N
), Brownian motion (
N
b
), and nonuniform parameter (
m
) are shown in Figures 14–16. Parabolic microrotation velocity profile is observed for the present flow problem. The velocity is maximum or minimum near the center of channel. The cause of coupling number on the microrotation velocity is illustrated in Figure 14. It is seen that in the microrotation velocity increases by increasing the coupling number N (i.e., microrotation velocity for the micropolar nanofluid is wider than that of the region for Newtonian nanofluid). In the case of N
=
0, there is no appreciable difference between Newtonian nanofluids and micropolar nanofluids. The microrotation velocity for the Brownian motion parameter N
b is plotted in Figure 15. It is seen that, with the increase in the Brownian motion parameter, microrotation velocity profile decreases and maximum microrotation velocity occurs at N
b
→
0. Figure 16 depicts the microrotation velocity field for different values of nonuniform parameter. It is also viewed that the microrotation velocity for a divergent channel (m
>
0) is higher compared to its value for a uniform channel (m
=
0).


## 4.4. Shear Stress Distribution

It is well known that the stress tensor is not symmetric in micropolar nanofluid. In Figures17–19, we have plotted the shear stresses τ
x
y at the left wall for values of the nonuniform parameter, micropolar parameter, and coupling number. One can observe from Figure 17 that the wall shear stress decreases with an increase in the nonuniform parameter. From Figure 18, we notice that shear stress is in oscillation behavior, which may be due to the creation of contraction and expansion walls. It is indicated that the shear stress decreases with an increase in the micropolar parameter n, while it increases as the coupling number N increases in Figure 19. Hence, we observed that the shear stress for a Newtonian nanofluid is less than that for a micropolar nanofluid.


## 4.5. Trapping Phenomena

The effects of nonuniform parameter, coupling parameter, and micropolar parameter on the streamlines are shown in Figure20. Besides, a comparison between uniform channel and the nonuniform channel is made in Figures 20(a) and 20(b). We note that the size of the trapping bolus increases with increasing nonuniform parameter and symmetry nature is also noticed with respect to uniform symmetric channel. The effects of N on trapping are presented in Figures 20(b) and 20(c). This figure reveals that the size of lower bolus decreases with an increase in N. From Figures 20(c) and 20(d), it is clear that the trapped bolus increases in size as n increases.


## 5. Concluding Remarks

A mathematical model is presented to study wall induced flow of a micropolar nanofluid in the most generalized (tapered asymmetric) channel with the presence of heat source and thermal radiation parameters. The flow model in the rotation of nanoparticles was included. Long wavelength and low Reynolds number assumptions are used in the mathematical modeling. In this investigation, special emphasis has been paid to study the flow features, the axial velocity, the nanoparticles temperature and concentration, the shear stress, and the trapping phenomena. The study leads to the following conclusions:(i)
The axial velocity of fluid decreases at the core part of channel whenm is increased as anticipated.
(ii)
The axial velocity increases near the left wall of channel and decreases near the right wall of the channel with increase of the coupling number and micropolar parameter.
(iii)
The nanoparticles mass transferσ has reverse behavior when compared to heat transfer.
(iv)
Coupling number and Brownian motion parameters have opposite effects on the microrotation velocity.
(v)
The walls shear stressτ
x
y decreases with the increase of nonuniform parameter m at the left wall but the opposite behavior is identified for coupling number.
(vi)
It may be interesting to note that the size of trapped bolus gets increased with increasing of nonuniform parameter.We ultimately conclude that our theoretical analysis bears the potential to be useful in the field of biomedical and industrial engineering.



---

*Source: 4123741-2016-09-05.xml*