# Analysis of Heat Transfer in Berman Flow of Nanofluids with Navier Slip, Viscous Dissipation, and Convective Cooling

Authors: O. D. Makinde; S. Khamis; M. S. Tshehla; O. Franks
Journal: Advances in Mathematical Physics
 (2014)

---

## Abstract

Heat transfer characteristics of a Berman flow of water based nanofluids containing copper (Cu) and alumina (Al2O3) as nanoparticles in a porous channel with Navier slip, viscous dissipation, and convective cooling are investigated. It is assumed that the exchange of heat with the ambient surrounding takes place at the channel walls following Newton’s law of cooling. The governing partial differential equations and boundary conditions are converted into a set of nonlinear ordinary differential equations using appropriate similarity transformations. These equations are solved analytically by regular perturbation methods with series improvement technique and numerically using an efficient Runge-Kutta Fehlberg integration technique coupled with shooting scheme. The effects of the governing parameters on the dimensionless velocity, temperature, skin friction, pressure drop, and Nusselt numbers are presented graphically and discussed quantitatively.


---

## Body

## 1. Introduction

The study of fluid flow and heat transfer between two porous boundaries has gained tremendous attention of researchers due to its wide applications in engineering and industrial processes. Some of the practical interests include problems dealing with transpiration cooling where the walls of a channel containing heated fluid are protected from overheating by passing cooler fluid over the exterior surface of the channel; fluid flow occurring during the separation of isotopes of Uranium-235 and Uranium-238 by gaseous diffusion in order to produce fuel for nuclear reactors; controlling boundary layer flow over aircraft wings by injection or suction of fluid out of or into the wing; lubrication of porous bearings; petroleum technology; ground water hydrology; seepage of water in river beds; purification and filtration processes; methods of decreasing rates of heat transfer in combustion chambers exhaust nozzles and porous walled flow reactors, and so forth. In a pioneering work, Berman [1] presented an exact solution of the Navier-Stokes equations that describes the steady two-dimensional flow of an incompressible viscous fluid along a channel with parallel rigid porous walls, under the action of uniform suction or injection of fluid at the surface. Sellars [2] extended Berman’s work to high suction Reynolds number. Yuan [3] considered the flow in a channel with porous walls. He obtained solutions for small suction and injection values and asymptotic solution valid at large injection values. Terrill [4, 5] gave an exact series solution for the fully developed laminar flow in a pipe of circular cross section with porous wall driven by a spatially variable and time independent suction or injection at the surface. Studies on developing flow in porous-walled ducts with suction and injection effects were carried out by Sorour et al. [6] and Zaturska et al. [7]. In view of the above interests, several researchers have also investigated the heat transfer problems between two permeable parallel walls under different physical situations [8–10].Moreover, the applications of conventional heat transfer fluids such as water and glycol mixture in engineering flow processes are limited due to their low thermal properties. A potential solution to improve these thermal properties is to add nanoparticles into the conventional fluids, hence forming the so-called nanofluids as coined by Choi [11]. Nanofluids contain thermally conducting submicron solid particles and have great potential as a high-energy carrier. The nanoparticles such as copper, alumina, titania, and copper oxide, unlike larger-sized particles, can be suspended stably within the conventional fluids without settling out of suspension. Nanofluids are free from numerous problems such as abrasion, clogging, and high pressure loss and are considered to be the next-generation working fluids in modern heat transfer technologies [12]. Experimental results [13–15] have shown that, even with small solid volume fraction of nanoparticles (usually less than 5%), the thermal conductivity of heat transfer fluids can be enhanced by 10–50%. Several authors [16–18] have also theoretically investigated the heat transfer enhancement of nanofluids under different physical conditions.Meanwhile, advances in the manufacture of microdevices have enabled experimental investigation of fluid flow in nano- and microscale, and many experimental results have provided evidences to support the slip condition [19, 20]. In order to describe the slip characteristics of fluid on the solid surface, Navier [21] introduced a more general boundary condition, namely, the fluid velocity component tangential to the solid surface, relative to the solid surface, which is proportional to the shear stress on the fluid-solid interface. The proportionality is called the slip length which describes the slipperiness of the surface. However, from the literature survey, it is found that no study has been conducted on the heat transfer characteristics of Berman flow of nanofluids with Navier slip, viscous dissipation, or convective cooling at the walls. Hence, the present study is an attempt in this direction. The flow of water base nanofluids containing copper (Cu) and alumina (Al2O3) as nanoparticles in a uniformly porous wall channel with Navier slip, viscous dissipation, and convective heat exchange with the ambient surrounding is investigated. Sections 2–4 give in more details the model nonlinear governing equations together with the analytical and numerical solution techniques employed to tackle the problem. In Section 5, we present graphically and discuss the main features of the flow and heat transfer characteristics in a range of governing parameters. Final conclusions are drawn in Section 6. To the best of our knowledge, the results of this paper are new and they have not been published before.


## 2. Problem Formulation

Consider a two-dimensional steady flow of a viscous incompressible water base nanofluids containing copper (Cu) and alumina (Al2O3) as nanoparticles in a uniformly porous wall channel. The channel wall is subjected to Navier slip and convectively exchange heat with the ambient surrounding. We choose the Cartesian coordinates system in such a way that the x-axis is taken along the channel and the y-axis is normal to it as shown in Figure 1.Figure 1
Schematic diagram of the physical system.The governing equations which are those of conservation of mass, momentum, and energy are [1–5, 8–10]
(1)
∂
u
∂
x
+
∂
v
∂
y
=
0
,
(2)
ρ
nf
(
u
∂
u
∂
x
+
v
∂
u
∂
y
)
=
-
∂
p
∂
x
+
μ
nf
(
∂
2
u
∂
y
2
+
∂
2
u
∂
x
2
)
,
(3)
ρ
nf
(
u
∂
v
∂
x
+
v
∂
v
∂
y
)
=
-
∂
p
∂
y
+
μ
nf
(
∂
2
v
∂
y
2
+
∂
2
v
∂
x
2
)
,
(4)
u
∂
T
∂
x
+
v
∂
T
∂
y
=
k
nf
(
ρ
c
p
)
nf
(
∂
2
T
∂
x
2
+
∂
2
T
∂
y
2
)
+
2
μ
nf
(
ρ
c
p
)
nf
[
(
∂
u
∂
x
)
2
+
(
∂
v
∂
y
)
2
]
+
μ
nf
(
ρ
c
p
)
nf
(
∂
v
∂
x
+
∂
u
∂
y
)
2
,
where (
u
,
v
) are the velocity components of the nanofluid in the (
x
,
y
) directions, respectively, V
(
>
0
) uniform wall suction velocity, P is the pressure, T is the nanofluid temperature, μ
nf is the effective dynamic viscosity of the nanofluid [22], k
nf is the effective thermal conductivity of the nanofluid [23], and ρ
nf and (
ρ
c
p
)
nf are the nanofluid density and the heat capacitance of the nanofluid, respectively, which are given by [13–18]. 						(5)
μ
nf
=
μ
f
(
1
-
φ
)
2.5
,
ρ
nf
=
(
1
-
φ
)
ρ
f
+
φ
ρ
s
,
α
nf
=
k
nf
(
ρ
c
p
)
nf
,
k
nf
k
f
=
(
k
s
+
2
k
f
)
-
2
φ
(
k
f
-
k
s
)
(
k
s
+
2
k
f
)
+
φ
(
k
f
-
k
s
)
,
(
ρ
c
p
)
nf
=
(
1
-
φ
)
(
ρ
c
p
)
f
+
φ
(
ρ
c
p
)
s
,
where φ is the nanoparticles solid volume fraction, ρ
f is the reference density of the fluid fraction, ρ
s is the reference density of the solid fraction, μ
f is the viscosity of the fluid fraction, k
f is the thermal conductivity of the fluid fraction, c
p is the specific heat at constant pressure, and k
s is the thermal conductivity of the solid volume fraction.Due to the symmetric nature of the flow, the boundary conditions at the channel centerline and at the porous wall may be written as [1–4]
(6)
∂
u
∂
y
(
x
,
0
)
=
0
,
∂
T
∂
y
(
x
,
0
)
=
0
,
v
(
x
,
0
)
=
0
,
β
u
(
x
,
a
)
=
μ
nf
∂
u
∂
y
(
x
,
a
)
,
v
(
x
,
a
)
=
V
,
-
k
nf
∂
T
∂
y
T
(
x
,
a
)
=
h
[
T
(
x
,
a
)
-
T
w
]
,
where T
w is the ambient surrounding temperature, h is the coefficient of heat transfer, and β is the Navier slip coefficient. Introduce the stream function ψ and vorticity Ω into the governing equations (1)–(4) as follows:
(7)
u
=
∂
ψ
∂
y
,
v
=
-
∂
ψ
∂
x
,
Ω
=
∂
u
∂
y
-
∂
v
∂
x
=
∂
2
ψ
∂
y
2
+
∂
2
ψ
∂
x
2
.After eliminating the pressureP from (2) and (3), we obtain
(8)
ρ
nf
(
∂
ψ
∂
y
∂
Ω
∂
x
-
∂
ψ
∂
x
∂
Ω
∂
y
)
=
μ
nf
(
∂
2
Ω
∂
y
2
+
∂
2
Ω
∂
x
2
)
,
(
∂
ψ
∂
y
∂
T
∂
x
-
∂
ψ
∂
x
∂
T
∂
y
)
=
k
nf
(
ρ
c
p
)
nf
(
∂
2
T
∂
y
2
+
∂
2
T
∂
x
2
)
+
4
μ
nf
(
ρ
c
p
)
nf
(
∂
2
ψ
∂
x
∂
y
)
2
+
μ
nf
(
ρ
c
p
)
nf
(
-
∂
2
ψ
∂
x
2
+
∂
2
ψ
∂
y
2
)
2
.The following dimensionless variables and parameters are introduced into (8), together with their corresponding boundary conditions:
(9)
η
=
y
a
,
ψ
-
=
ψ
V
a
,
Ω
-
=
Ω
a
V
,
Φ
=
T
-
T
w
T
w
,
R
=
V
a
υ
f
,
X
=
x
a
,
υ
f
=
μ
f
ρ
f
,
Pr
=
μ
f
c
p
k
f
,
m
1
=
(
1
-
φ
+
φ
ρ
s
ρ
f
)
(
1
-
φ
)
2.5
,
λ
=
μ
f
β
a
,
P
-
=
a
P
V
μ
f
,
m
2
=
m
4
(
1
-
φ
+
φ
(
ρ
c
p
)
s
(
ρ
c
p
)
f
)
,
Bi
=
h
a
k
f
,
Ec
=
V
2
c
p
T
w
,
m
3
=
m
4
(
1
-
φ
)
2.5
,
m
4
=
(
k
s
+
2
k
f
)
+
φ
(
k
f
-
k
s
)
(
k
s
+
2
k
f
)
-
2
φ
(
k
f
-
k
s
)
,
and we obtain
(10)
m
1
R
(
∂
ψ
-
∂
η
∂
Ω
-
∂
X
-
∂
ψ
-
∂
X
∂
Ω
-
∂
η
)
=
(
∂
2
Ω
-
∂
η
2
+
∂
2
Ω
-
∂
X
2
)
,
(11)
m
2
R
Pr
(
∂
ψ
-
∂
η
∂
Φ
∂
X
-
∂
ψ
-
∂
X
∂
Φ
∂
η
)
=
(
∂
2
Φ
∂
η
2
+
∂
2
Φ
∂
X
2
)
+
4
m
3
Ec
Pr
(
∂
2
ψ
-
∂
X
∂
η
)
2
+
m
3
Ec
Pr
(
-
∂
2
ψ
-
∂
X
2
+
∂
2
ψ
-
∂
η
2
)
2
with
(12)
∂
2
ψ
-
∂
η
2
(
X
,
0
)
=
0
,
∂
Φ
∂
η
(
X
,
0
)
=
0
,
∂
ψ
-
∂
X
(
X
,
0
)
=
0
,
∂
ψ
-
∂
η
(
X
,
1
)
=
λ
(
1
-
φ
)
2.5
∂
2
ψ
-
∂
η
2
(
X
,
1
)
,
∂
ψ
-
∂
X
(
x
,
a
)
=
-
1
,
∂
Φ
∂
η
(
X
,
1
)
=
-
m
4
Bi
Φ
(
X
,
1
)
,
where R is the flow Reynolds number such that R
>
0 represents wall suction and R
<
0 represents wall injection; Ec is the Eckert number; Bi is the Biot number; Pr is the base fluid Prandtl number; m
1, m
2, m
3, and m
4 can be easily determined from the thermophysical properties of the base fluid and the nanoparticles; and λ is the Navier slip parameter such that λ
=
0 corresponds to no slip, while full lubrication is described in the limit λ
→
∞. We seek a similarity form of solution due to Berman [1]; that is,
(13)
ψ
-
(
X
,
η
)
=
X
F
(
η
)
,
W
(
X
,
η
)
=
X
d
F
d
η
,
Φ
(
X
,
η
)
=
H
(
η
)
+
X
2
θ
(
η
)
.Equations (10)-(11) together with the boundary conditions in (12) then become
(14)
d
4
F
d
η
4
=
m
1
R
(
d
F
d
η
d
2
F
d
η
2
-
F
d
3
F
d
η
3
)
,
(15)
d
2
θ
d
η
2
=
-
m
3
Ec
Pr
(
d
2
F
d
η
2
)
2
+
m
2
R
Pr
(
2
θ
d
F
d
η
-
F
d
θ
d
η
)
,
(16)
d
2
H
d
η
2
=
-
2
θ
-
4
m
3
Ec
Pr
(
d
F
d
η
)
2
-
m
2
R
Pr
F
d
H
d
η
,
(17)
d
2
F
d
η
2
(
0
)
=
0
,
d
θ
d
η
(
0
)
=
0
,
d
H
d
η
(
0
)
=
0
,
F
(
0
)
=
0
,
(18)
d
F
d
η
(
1
)
=
b
d
2
F
d
η
2
(
1
)
,
F
(
1
)
=
-
1
,
d
θ
d
η
(
1
)
=
-
m
4
Bi
θ
(
1
)
,
d
H
d
η
(
1
)
=
-
m
4
Bi
H
(
1
)
,
where b
=
λ
/
(
1
-
φ
)
2.5. The dimensionless fluid axial pressure gradient is given as follows:
(19)
-
∂
P
-
∂
X
=
X
A
,
where
(20)
(
1
-
φ
)
2.5
A
=
-
d
3
F
d
η
3
+
m
1
R
[
(
d
F
d
η
)
2
-
F
d
2
F
d
η
2
]
.Other quantities of practical interest in this study are the local skin friction coefficientC
f and the local Nusselt number Nu, which are defined as
(21)
C
f
=
a
τ
w
μ
f
V
,
Nu
=
a
q
w
k
f
T
w
,
where τ
w is the skin friction and q
w is the heat flux at the channel walls which are given by
(22)
τ
w
=
μ
nf
∂
u
∂
y
|
y
=
a
,
q
w
=
-
k
nf
∂
T
∂
y
|
y
=
a
.Using (9) and (13), we substitute (22) into (21) and obtain
(23)
C
f
=
X
(
1
-
φ
)
2.5
d
2
F
d
η
2
(
1
)
,
Nu
=
-
k
nf
k
f
[
d
H
d
η
(
1
)
+
X
2
d
θ
d
η
(
1
)
]
.In the following section, the boundary value problem (14)–(18) was solved analytically using regular perturbation method and numerically by the Runge-Kutta-Fehlberg method with shooting technique [24]. The results are utilised to compute the fluid pressure gradient, local skin friction, and local Nusselt number as highlighted in (20) and (23).


## 3. Perturbation Method

Due to the nonlinear nature of the model equations (14)–(18), it is convenient to form a power series expansion in the parameter R; that is,
(24)
F
(
η
)
=
∑
i
=
0
∞
F
i
R
i
,
θ
(
η
)
=
∑
i
=
0
∞
θ
i
R
i
.Substituting the solution series in (24) into (14)–(18) and collecting the coefficients of like powers of R, we obtain the following.Zeroth Order. Consider
(25)
d
4
F
0
d
η
4
=
0
,
d
2
θ
0
d
η
2
=
-
m
3
Ec
Pr
(
d
2
F
0
d
η
2
)
2
,
d
2
H
0
d
η
2
=
-
2
θ
0
-
4
m
3
Ec
Pr
(
d
F
0
d
η
)
2
,
with
(26)
d
2
F
0
d
η
2
(
0
)
=
0
,
d
θ
0
d
η
(
0
)
=
0
,
d
H
0
d
η
(
0
)
=
0
,
F
0
(
0
)
=
0
,
d
F
0
d
η
(
1
)
=
b
d
2
F
0
d
η
2
(
1
)
,
F
0
(
1
)
=
-
1
,
d
θ
0
d
η
(
1
)
=
-
m
4
Bi
θ
0
(
1
)
,
d
H
0
d
η
(
1
)
=
-
m
4
Bi
H
0
(
1
)
.Higher Order
(
n
≥
1
). Consider
(27)
d
4
F
n
d
η
4
=
m
1
R
∑
i
=
0
n
-
1
(
d
F
i
d
η
d
2
F
n
-
i
-
1
d
η
2
-
F
i
d
3
F
n
-
i
-
1
d
η
3
)
,
d
2
θ
n
d
η
2
=
-
m
3
Ec
Pr
∑
i
=
0
n
(
d
2
F
i
d
η
2
d
2
F
n
-
i
d
η
2
)
+
m
2
R
Pr
∑
i
=
0
n
-
1
(
2
θ
i
d
F
n
-
i
-
1
d
η
-
F
i
d
θ
n
-
i
-
1
d
η
)
,
d
2
H
n
d
η
2
=
-
2
θ
n
-
4
m
3
Ec
Pr
∑
i
=
0
n
(
d
F
i
d
η
d
F
n
-
i
d
η
)
-
m
2
R
Pr
∑
i
=
0
n
-
1
(
F
i
d
H
n
-
i
-
1
d
η
)
,
with
(28)
d
2
F
n
d
η
2
(
0
)
=
0
,
d
θ
n
d
η
(
0
)
=
0
,
d
H
n
d
η
(
0
)
=
0
,
F
n
(
0
)
=
0
,
d
F
n
d
η
(
1
)
=
b
d
2
F
n
d
η
2
(
1
)
,
F
n
(
1
)
=
0
,
d
θ
n
d
η
(
1
)
=
-
m
4
Bi
θ
n
(
1
)
,
d
H
n
d
η
(
1
)
=
-
m
4
Bi
H
n
(
1
)
.The equations are solved iteratively and the series solutions for the velocity and temperature fields are given as(29)
F
(
η
)
=
η
(
η
2
+
6
b
-
3
)
2
-
6
b
+
η
R
m
1
(
η
2
-
1
)
280
(
3
b
-
1
)
3
×
(
3
η
4
b
-
η
4
+
3
η
2
b
-
η
2
-
18
b
+
2
)
+
O
(
R
2
)
,
(30)
θ
(
η
)
=
3
(
4
+
m
4
Bi
+
η
4
m
4
Bi
)
m
3
Ec
Pr
4
m
4
Bi
(
3
b
-
1
)
2
+
O
(
R
)
,
(31)
H
(
η
)
=
-
3
m
3
Ec
Pr
2
m
4
Bi
(
3
b
-
1
)
2
×
(
2
η
2
+
7
2
m
4
Bi
η
2
+
1
6
η
6
m
4
Bi
+
2
m
4
Bi
η
4
b
-
η
4
m
4
Bi
+
12
m
4
Bi
b
2
η
2
-
12
m
4
Bi
b
η
2
7
2
)
-
3
m
3
Ec
Pr
2
m
4
3
Bi
3
(
3
b
-
1
)
2
(
9
b
2
-
6
b
+
1
)
×
(
6
+
9
m
4
Bi
-
24
m
4
Bi
b
+
36
m
4
Bi
b
2
+
4
m
4
2
Bi
2
-
15
m
4
2
Bi
2
b
+
18
m
4
2
Bi
2
b
2
)
+
O
(
R
)
,
where b
=
λ
/
(
1
-
φ
)
2.5. We remark here that if we set the parameters φ
=
0 and b
=
0 in (29), we will recover the solution for the classical case of conventional fluid given in [1–4]. Using a computer symbolic algebra package (MAPLE), several terms of the above solution series in (29)–(31) are obtained. From (29)–(31) together with (20) and (23), we obtained the series solutions for the skin friction, Nusselt number, and axial pressure gradient as follows:
(32)
C
f
=
X
(
1
-
φ
)
2.5
[
R
2
m
1
2
(
3465
b
2
-
1743
b
+
394
)
40425
(
3
b
-
1
)
5
3
1
-
3
b
-
3
R
m
1
35
(
3
b
-
1
)
3
-
R
2
m
1
2
(
3465
b
2
-
1743
b
+
394
)
40425
(
3
b
-
1
)
5
+
O
(
R
3
)
R
2
m
1
2
(
3465
b
2
-
1743
b
+
394
)
40425
(
3
b
-
1
)
5
]
,
A
=
1
(
1
-
φ
)
2.5
×
[
3
1
-
3
b
-
9
R
m
1
(
105
b
3
-
140
b
2
+
63
b
-
9
)
35
(
3
b
-
1
)
3
+
R
2
m
1
2
(
1155
b
3
-
1309
b
2
+
364
b
-
26
)
13475
(
3
b
-
1
)
5
+
O
(
R
3
)
R
2
m
1
2
(
1155
b
3
-
1309
b
2
+
364
b
-
26
)
13475
(
3
b
-
1
)
5
]
,
k
f
k
nf
Nu
=
6
m
3
Ec
Pr
G
1
m
4
Bi
(
3
b
-
1
)
2
+
2
R
m
3
Ec
Pr
G
2
105
m
4
2
Bi
2
(
3
b
-
1
)
4
+
X
2
(
3
m
3
Ec
Pr
(
3
b
-
1
)
2
+
3
R
m
3
Ec
Pr
G
3
35
m
4
Bi
(
3
b
-
1
)
4
)
+
O
(
R
2
)
,
where the expression for G
i, i
=
1,2
,
3 is given in the appendix. We are aware that these power series solutions are valid for very small parameter values of R. However, using Hermite-Padé approximation technique (see Makinde [25]) that is based on the series summation and improvement method, the usability of the extended solution series is improved beyond small parameter values of R.


## 4. Numerical Procedure

An efficient finite difference approach based on Runge-Kutta-Fehlberg method with shooting technique [24] has been employed to numerically solve the coupled nonlinear ordinary differential equations (14)–(16) subject to the boundary conditions (17)-(18) for different values of governing parameters. The boundary value problem is first transformed into an initial value problem (IVP). Let
(33)
z
1
=
F
,
z
2
=
d
F
d
η
,
z
3
=
d
2
F
d
η
2
,
z
4
=
d
3
F
d
η
3
,
z
5
=
θ
,
z
6
=
d
θ
d
η
,
z
7
=
H
,
z
8
=
d
H
d
η
.Substituting (33) into (14)–(18), we obtain a system of first-order differential equations, respectively, as follows:
(34)
d
z
1
d
η
=
z
2
,
d
z
2
d
η
=
z
3
,
d
z
3
d
η
=
z
4
,
d
z
4
d
η
=
m
1
R
(
z
2
z
3
-
z
1
z
4
)
,
d
z
5
d
η
=
z
6
,
d
z
6
d
η
=
-
m
3
Ec
Pr
z
3
2
+
m
2
R
Pr
(
2
z
5
z
2
-
z
1
z
6
)
,
d
z
7
d
η
=
z
8
,
d
z
8
d
η
=
-
2
z
5
-
4
m
3
Ec
Pr
z
2
2
-
m
2
R
Pr
(
z
1
z
8
)
,
subject to the initial conditions:
(35)
z
1
(
0
)
=
0
,
z
2
(
0
)
=
s
1
,
z
3
(
0
)
=
0
,
z
4
(
0
)
=
s
2
,
z
5
(
0
)
=
s
3
,
z
6
(
0
)
=
0
,
z
7
(
0
)
=
s
4
,
z
8
(
0
)
=
0
.By applying the shooting method with the Newton-Raphson algorithm to guess the unspecified conditionss
1, s
2, s
3, and s
4 in (35), the resulting initial value problem is then integrated numerically until the boundary conditions at η
=
1 are achieved. The step size is taken as Δ
η
=
0.001 and the convergence criteria were set to 10−7.


## 5. Results and Discussion

The flow of water base nanofluids containing Cu and Al2O3 as nanoparticles and their heat transfer characteristics in a symmetrically porous channel with Navier slip and convective cooling at the surface are investigated. The governing partial differential equations and the corresponding boundary conditions are converted into a set of nonlinear ordinary differential equations and tackled both analytically using the perturbation method coupled with series improvement technique and numerically using Runge-Kutta Fehlberg integration technique coupled with shooting scheme. Thermophysical properties of base fluid and nanoparticles are presented in Table 1. For pure water, the momentum diffusivity is dominant and convection is very effective in transferring within the fluid in comparison to pure conduction. Following [16–18], we take Pr
=
6.2 in the numerical computation. Note that when φ
=
0, no nanoparticle is present in the based fluid (water). The solid volume fraction in the base fluid is taken as φ
=
0 to 0.3 (i.e., ranging from 0 to 30 percent). In order to get a clear insight into the entire flow structure and thermal development, we have assigned numerical values to other parameters encountered in the problem. Numerical solutions are displayed in Tables 2 and 3 together with Figures 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, and 23. The numerical values of suction driven normal velocity profile (
R
=
1
) are displayed in Table 2. In the absence of nanoparticles (
φ
=
0
) and Navier slip (
λ
=
0
) at the channel walls, the results agreed well with the one already in the literature (see [1–3]) for the flow of conventional fluid in a symmetrical porous wall channel. However, Cu-water nanofluid is affected more by the combined effects of suction and Navier slip at wall in comparison to Al2O3-water nanofluid. Table 3 shows the perfect agreement between the series solution and the numerical solution for the axial velocity profiles along the channel centerline with increasing concentration of nanoparticles for both Cu-water and Al2O3-water nanofluids.Table 1
Thermophysical properties of the fluid phase (water) and nanoparticles [13–18].
Physical properties
Fluid phase (Water)
Cu
Al2O3
c
p (J/kgK)
4179
385
765
ρ (kg/m3)
997.1
8933
3970
k (W/mK)
0.613
400
40Table 2
Computation showing the normal velocity profiles forR
=
1.
η
-
F
(
η
)
λ
=
0, φ
=
0  Reference [1–3]
-
F
(
η
)
λ
=
0, φ
=
0  Present
-
F
(
η
)
λ
=
0.05, φ
=
0.1  Cu-water
-
F
(
η
)
λ
=
0.05, φ
=
0.1  Al2O3-water
0.0
0.00000
0.000000
0.000000
0.000000
0.1
0.14874
0.148739
0.160680
0.160920
0.2
0.29455
0.294548
0.317732
0.318184
0.3
0.43449
0.434493
0.467527
0.468132
0.4
0.56564
0.565639
0.606428
0.607103
0.5
0.68504
0.685042
0.730781
0.731419
0.6
0.78974
0.789739
0.836887
0.837373
0.7
0.87672
0.876724
0.920958
0.921192
0.8
0.94292
0.942920
0.979045
0.978994
0.9
0.98514
0.985137
1.006937
1.006713
1.0
1.00000
1.000000
1.000000
1.000000Table 3
Computation showing the axial velocity profiles forR
=
1, λ
=
0.05, X
=
1, and  η
=
0.
φ
W
(
1,0
) Cu-water(series)
W
(
1,0
) Cu-water(numerical)
W
(
1,0
) Al2O3-water(series)
W
(
1,0
) Al2O3-water(numerical)
0.00
1.581405
1.581405
1.581405
1.581405
0.01
1.583733
1.583733
1.584084
1.584084
0.05
1.594688
1.594688
1.596218
1.596218
0.10
1.612850
1.612850
1.615303
1.615303
0.15
1.637660
1.637660
1.640343
1.640343
0.20
1.672037
1.672037
1.674097
1.674097
0.25
1.721196
1.721196
1.721395
1.721395
0.30
1.795420
1.795420
1.791485
1.791485Figure 2
Nanofluids axial velocity profiles.Figure 3
Axial velocity profiles with increasingφ.Figure 4
Axial velocity profiles with increasingR.Figure 5
Axial velocity profiles with increasingλ.Figure 6
Nanofluids normal velocity profiles.Figure 7
Normal velocity profiles with increasingφ.Figure 8
Normal velocity profiles with increasingR.Figure 9
Normal velocity profiles with increasingλ.Figure 10
Nanofluids temperature profiles.Figure 11
Temperature profiles with increasingφ.Figure 12
Temperature profiles with increasingBi.Figure 13
Temperature profiles with increasingR.Figure 14
Temperature profiles with increasingλ.Figure 15
Temperature profiles with increasingEc.Figure 16
Temperature profiles with increasingX.Figure 17
Skin friction with increasingφ.Figure 18
Skin friction with increasingR and λ.Figure 19
Axial pressure gradient with increasingφ.Figure 20
Axial pressure gradient with increasingR and λ.Figure 21
Nusselt number with increasingφ.Figure 22
Nusselt number with increasingR and λ.Figure 23
Nusselt number with increasingX, Bi, and Ec.


### 5.1. Velocity Profiles with Parameter Variation

The effects of parameter variation on both the axial and normal velocity components are displayed in Figures2–9. Generally, it is interesting to note that the effect of Navier slip is to cause flow reversal at the channel walls. In Figure 2, it is observed that Cu-water moves faster along the channel centerline region and subjected to higher flow reversal at the wall as compared to Al2O3-water in the presence of suction. With Cu-water as the working nanofluid, it is observed that increasing nanoparticles volume fraction concentration from 0 to 30% increases both the axial velocity along the centerline region and the flow reversal at the channel walls as shown in Figure 3. Similar trend is observed in Figures 4-5 with a growing suction parameter and Navier slip parameter. Figures 6–9 show the nanofluids normal velocity profiles. With suction and Navier slip, the Cu-water moves faster to the wall as compared to Al2O3-water as illustrated in Figure 6. A further increase in normal velocity towards the walls is observed with a growing in suction parameter, Navier slip parameter, and nanoparticles volume fraction as shown in Figures 7–9.


### 5.2. Temperature Profiles with Parameter Variation

Figures10–16 illustrate the nanofluids temperature profiles across the channel with different parameter variation. Generally, a decrease in the fluid temperature near the channel walls is observed due to convective heat loss to ambient surrounding. It is noteworthy that the temperature of Cu-water nanofluid is generally higher than that of Al2O3-water nanofluid under the same flow condition as shown in Figure 10. In Figure 11, it is observed that the nanofluid temperature decreases with growing in nanoparticles volume fraction. Similar effect of a decrease in nanofluid temperature is observed in Figure 12 with Cu-water as working nanofluid as the Biot number increases. This is expected, since an increase in Biot number indicates a rise in convective cooling due to heat loss to the ambient surrounding from the channel walls. Meanwhile, a combine increase in the suction, Navier slip, and viscous dissipation as shown in Figures 13–15 causes an increase in the nanofluid temperature. This may be attributed to the fact that R, λ, and Ec increase the internal heat generation within the fluid due velocity gradient increases, leading to a rise in temperature. Figure 16 elucidates the temperature profiles with increasing axial distance along the channel. The nanofluid temperature decreases within the channel centerline region and increases near the wall region with increasing axial distance. Moreover, it is interesting to note that at point η
=
0.5 within the channel, the nanofluid temperature is independent of the axial distance.


### 5.3. Skin Friction, Pressure Gradient, and Nusselt Number

Figures17-18 depict the skin friction profiles for both Cu-water and Al2O3-water nanofluids at the channel walls. The skin friction generally increases with an increase in nanoparticles volume fraction; however, it is noticed that the skin friction produced by Cu-water is more intense than the one produced by Al2O3-water as shown in Figure 17. This is expected since the velocity gradient of Cu-water near the channel walls is higher than that of Al2O3-water. In Figure 18, it is observed that the skin friction generally increases with an increase in Navier slip. Meanwhile, a growing in suction (
R
>
0
) increases the skin friction, while a growing in injection (
R
<
0
) decreases the skin friction. The pressure drop along the channel is illustrated in Figures 19 and 20. For both Cu-water and Al2O3-water nanofluids, the pressure drop increases with increasing nanoparticles volume fraction. Interestingly, the pressure drop produced by Al2O3-water is slightly higher than that of Cu-water as shown in Figure 19. Figure 20 shows a general increase in pressure drop with a rise in Navier slip. A growing in suction (
R
>
0
) decreases the pressure drop, while a growing in injection (
R
<
0
) increases the pressure drop along the channel. Figures 21–23 elucidate the heat transfer rate at the channel walls with different parameter variation. It is observed that the wall heat transfer rate (Nu) decreases with an increase in the nanoparticles volume fraction as shown in Figure 21. A slight increase in Al2O3-water Nusselt number is noticed as compared to Cu-water Nusselt number. In Figure 22, it is seen that the increase in the Navier slip parameter results in increase in Nusselt number. This may be due to a rise in the nanofluid temperature gradient at the channel walls. Meanwhile, as the suction increases, the heat flux at the wall increases, while a decrease in wall heat is observed with a rise in injection. The strength of the wall heat flux is enhanced with increasing axial distance, viscous dissipation, and convective cooling as illustrated in Figure 23. This can be attributed to a rise in the temperature gradient due to convective heat exchange with the ambient along the channel walls.


## 5.1. Velocity Profiles with Parameter Variation

The effects of parameter variation on both the axial and normal velocity components are displayed in Figures2–9. Generally, it is interesting to note that the effect of Navier slip is to cause flow reversal at the channel walls. In Figure 2, it is observed that Cu-water moves faster along the channel centerline region and subjected to higher flow reversal at the wall as compared to Al2O3-water in the presence of suction. With Cu-water as the working nanofluid, it is observed that increasing nanoparticles volume fraction concentration from 0 to 30% increases both the axial velocity along the centerline region and the flow reversal at the channel walls as shown in Figure 3. Similar trend is observed in Figures 4-5 with a growing suction parameter and Navier slip parameter. Figures 6–9 show the nanofluids normal velocity profiles. With suction and Navier slip, the Cu-water moves faster to the wall as compared to Al2O3-water as illustrated in Figure 6. A further increase in normal velocity towards the walls is observed with a growing in suction parameter, Navier slip parameter, and nanoparticles volume fraction as shown in Figures 7–9.


## 5.2. Temperature Profiles with Parameter Variation

Figures10–16 illustrate the nanofluids temperature profiles across the channel with different parameter variation. Generally, a decrease in the fluid temperature near the channel walls is observed due to convective heat loss to ambient surrounding. It is noteworthy that the temperature of Cu-water nanofluid is generally higher than that of Al2O3-water nanofluid under the same flow condition as shown in Figure 10. In Figure 11, it is observed that the nanofluid temperature decreases with growing in nanoparticles volume fraction. Similar effect of a decrease in nanofluid temperature is observed in Figure 12 with Cu-water as working nanofluid as the Biot number increases. This is expected, since an increase in Biot number indicates a rise in convective cooling due to heat loss to the ambient surrounding from the channel walls. Meanwhile, a combine increase in the suction, Navier slip, and viscous dissipation as shown in Figures 13–15 causes an increase in the nanofluid temperature. This may be attributed to the fact that R, λ, and Ec increase the internal heat generation within the fluid due velocity gradient increases, leading to a rise in temperature. Figure 16 elucidates the temperature profiles with increasing axial distance along the channel. The nanofluid temperature decreases within the channel centerline region and increases near the wall region with increasing axial distance. Moreover, it is interesting to note that at point η
=
0.5 within the channel, the nanofluid temperature is independent of the axial distance.


## 5.3. Skin Friction, Pressure Gradient, and Nusselt Number

Figures17-18 depict the skin friction profiles for both Cu-water and Al2O3-water nanofluids at the channel walls. The skin friction generally increases with an increase in nanoparticles volume fraction; however, it is noticed that the skin friction produced by Cu-water is more intense than the one produced by Al2O3-water as shown in Figure 17. This is expected since the velocity gradient of Cu-water near the channel walls is higher than that of Al2O3-water. In Figure 18, it is observed that the skin friction generally increases with an increase in Navier slip. Meanwhile, a growing in suction (
R
>
0
) increases the skin friction, while a growing in injection (
R
<
0
) decreases the skin friction. The pressure drop along the channel is illustrated in Figures 19 and 20. For both Cu-water and Al2O3-water nanofluids, the pressure drop increases with increasing nanoparticles volume fraction. Interestingly, the pressure drop produced by Al2O3-water is slightly higher than that of Cu-water as shown in Figure 19. Figure 20 shows a general increase in pressure drop with a rise in Navier slip. A growing in suction (
R
>
0
) decreases the pressure drop, while a growing in injection (
R
<
0
) increases the pressure drop along the channel. Figures 21–23 elucidate the heat transfer rate at the channel walls with different parameter variation. It is observed that the wall heat transfer rate (Nu) decreases with an increase in the nanoparticles volume fraction as shown in Figure 21. A slight increase in Al2O3-water Nusselt number is noticed as compared to Cu-water Nusselt number. In Figure 22, it is seen that the increase in the Navier slip parameter results in increase in Nusselt number. This may be due to a rise in the nanofluid temperature gradient at the channel walls. Meanwhile, as the suction increases, the heat flux at the wall increases, while a decrease in wall heat is observed with a rise in injection. The strength of the wall heat flux is enhanced with increasing axial distance, viscous dissipation, and convective cooling as illustrated in Figure 23. This can be attributed to a rise in the temperature gradient due to convective heat exchange with the ambient along the channel walls.


## 6. Conclusion

The combined effects of viscous dissipation, Navier slip, and convective cooling on Berman flow and heat transfer of water base nanofluids containing Cu and Al2O3 as nanoparticles are investigated. The nonlinear model problem is tackled both analytically using perturbation series method and numerically using Runge-Kutta Fehlberg integration technique coupled with shooting scheme. We summarize below some of the essential features of physical interest from the above analysis.(i)
Cu-water nanofluid moves faster with enhanced flow reversal at the walls as compared to Al2O3-water nanofluid.
(ii)
Nanofluids velocity and flow reversal at the walls increase with suction,λ, and φ.
(iii)
Cu-water produced higher temperature as compared to Al2O3-water. The nanofluids temperature increases with suction, λ, and Ec, but it decreases with Bi and φ.
(iv)
The skin friction produced by Cu-water is more intense than that of Al2O3-water. The skin friction increases with suction (
R
>
0
), λ, and φ but decreases with injection (
R
<
0
).
(v)
The pressure drop produced by Al2O3-water in more than that of Cu-water. The pressure drop is enhanced by injection, λ, and φ but decreases by suction.
(vi)
The Nusselt number increases with suction,λ, Bi, Ec, and X, but it decreases with injection and φ. A slight increase in Nu for Al2O3-water is noticed as compared to Cu-water.



---

*Source: 809367-2014-03-31.xml*