# Three-Dimensional Flow and Heat Transfer Past a Permeable Exponentially Stretching/Shrinking Sheet in a Nanofluid

Authors: Syahira Mansur; Anuar Ishak; Ioan Pop
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

The three-dimensional flow and heat transfer of a nanofluid over a stretching/shrinking sheet is investigated. Numerical results are obtained using bvp4c in MATLAB. The results show nonunique solutions for the shrinking case. The effects of the stretching/shrinking parameter, suction parameter, Brownian motion parameter, thermophoresis parameter, and Lewis number on the local skin friction coefficient and the local Nusselt number are studied. Suction increases the solution domain. Furthermore, as the sheet is shrunk in thex-direction, suction increases the skin friction coefficient in the same direction while decreasing the skin friction coefficient in the y-direction. The local Nusselt number is consistently lower for higher values of thermophoresis parameter and Lewis number. On the other hand, the local Nusselt number increases as the Brownian motion parameter increases.


---

## Body

## 1. Introduction

Nanofluids are dispersions of nanometer-sized particles in a base fluid such as water, ethylene glycol, and propylene glycol, to increase their thermal conductivities. Choi and Eastman [1] showed that the addition of a small amount (less than 1% by volume) of nanoparticles to conventional heat transfer liquids increased the thermal conductivity of the fluid up to approximately two times. In his paper, Buongiorno [2] developed a model for convective transport in nanofluids which takes into account the Brownian diffusion and thermophoresis effects. Buongiorno’s nanofluid model was used in many recent papers, for example, Nield and Kuznetsov [3–5], Khan and Pop [6], Bachok et al. [7–9], Mansur and Ishak [10, 11], and Zaimi et al. [12] among others.The boundary layer flow over a stretching sheet is significant in applications such as extrusion, wire drawing, metal spinning, and hot rolling [13]. Wang [14, 15], Mandal and Mukhopadhyay [16], P. S. Gupta and A. S. Gupta [17], Andersson [18], Ishak et al. [19], and Makinde and Aziz [20] are among various names who published their papers on a stretching sheet. Miklavčič and Wang [21] studied flow over a shrinking sheet in which they observed that the vorticity is not confined within a boundary layer and the steady flow cannot exist without exerting adequate suction at the boundary. As the studies of shrinking sheet garner considerable attention, this finding proves to be crucial to these researches. In response to Miklavčič and Wang, numerous studies on these problems have been conducted by researches, namely, Wang [22], Fang et al. [23], Bachok et al. [24], Bhattacharyya et al. [25], Zaimi et al. [26], and Roşca and Pop [27] among others.All the above-mentioned studies dealt with problems involving linear stretching/shrinking sheet. The boundary layer flow induced by a stretching/shrinking sheet is very important in engineering processes [28] and has attracted many researchers to delve into this study such as Bachok et al. [29], Bhattacharyya and Vajravelu [30], and Rohni et al. [31]. However, most studies revolved around two-dimensional flows. Motivated by this, the objective of this paper is to solve the problem of three-dimensional flow and heat transfer past a permeable exponentially stretching/shrinking sheet in a nanofluid. The dependency of the local skin friction coefficient and the local Nusselt number on several parameters, namely, the stretching/shrinking, Brownian motion, and thermophoresis parameters is the main focus of the present investigation. Numerical solutions are presented graphically and in tabular forms to show the effects of these parameters on the local skin friction coefficient and the local Nusselt number.


## 2. Problem Formulation

We consider the steady three-dimensional boundary layer flow of a viscous nanofluid past a permeable stretching/shrinking flat surface in a quiescent fluid. A locally orthogonal set of coordinates(
x
,
y
,
z
) is chosen with the origin O in the plane of the stretching/shrinking sheet. The x- and y-coordinates are in the plane of the sheet, while the coordinate z is measured in the perpendicular direction to the stretching/shrinking surface as shown in Figure 1. It is assumed that the flat surface is stretched/shrunk continuously in the both x- and y-directions with the velocities u
(
x
)
=
u
w
(
x
) and v
(
y
)
=
v
w
(
y
), respectively. It is also assumed that the mass flux velocity is w
w
(
x
,
y
), where w
w
(
x
,
y
)
<
0 is for suction and w
w
(
x
,
y
)
>
0 is for injection or withdrawal of the fluid. Further, we assume that the constant surface temperature and the constant surface volume fraction are T
w and C
w, while the constant temperature and the constant surface volume fraction of the ambient (inviscid) fluid are T
∞ and C
∞, respectively. Under these conditions, the boundary layer equations are
(1)
∂
u
∂
x
+
∂
v
∂
y
+
∂
w
∂
z
=
0
,
(2)
u
∂
u
∂
x
+
v
∂
u
∂
y
+
w
∂
u
∂
z
=
ν
∂
2
u
∂
z
2
,
(3)
u
∂
v
∂
x
+
v
∂
v
∂
y
+
w
∂
v
∂
z
=
ν
∂
2
v
∂
z
2
,
(4)
u
∂
w
∂
x
+
v
∂
w
∂
y
+
w
∂
w
∂
z
=
ν
∂
2
w
∂
z
2
,
(5)
u
∂
T
∂
x
+
v
∂
T
∂
y
+
w
∂
T
∂
z
=
α
∂
2
T
∂
z
2
+
τ
[
D
B
∂
C
∂
z
∂
T
∂
z
+
D
T
T
∞
(
∂
T
∂
z
)
2
]
,
(6)
u
∂
C
∂
x
+
v
∂
C
∂
y
+
w
∂
C
∂
z
=
D
B
∂
2
C
∂
z
2
+
D
T
T
∞
∂
2
T
∂
z
2
,
along with the boundary conditions
(7)
u
=
u
w
(
x
)
=
λ
1
U
w
(
x
)
,
v
=
v
w
(
y
)
=
λ
2
V
w
(
y
)
,
w
=
w
w
=
w
0
,
T
=
T
w
,
D
B
∂
C
∂
z
+
D
T
T
∞
∂
T
∂
z
=
0
at
z
=
0
u
⟶
0
,
w
⟶
0
,
T
⟶
T
∞
,
C
⟶
C
∞
as
z

⟶

∞
.
Here u
,
v, and w are the velocity components along x-,
y-, and z-axes, respectively; ν is the kinematic viscosity of the fluid, λ
1 is the constant stretching (
λ
1
>
0
) or shrinking (
λ
1
<
0
) parameter in the x-direction, and λ
2 is the constant stretching (
λ
2
>
0
) or shrinking (
λ
2
<
0
) parameter in the y- direction, respectively. Further, we assume that U
w
(
x
,
y
) and V
w
(
x
,
y
) are of the following form:
(8)
U
w
(
x
,
y
)
=
V
w
(
x
,
y
)
=
U
0
e
(
x
+
y
)
/
L
,
where L is the characteristic length and U
0 is the characteristic velocity of the stretching/shrinking sheet.Figure 1
Geometry of the problem.We introduce now the following similarity variables:(9)
u
=
U
0
e
(
x
+
y
)
/
L
f
′
(
η
)
,
v
=
U
0
e
(
x
+
y
)
/
L
g
′
(
η
)
,
w
=
-
(
ν
U
0
2
L
)
1
/
2
e
(
x
+
y
)
/
2
L
[
f
(
η
)
+
η
f
′
(
η
)
+
g
(
η
)
+
η
g
′
(
η
)
]
,
θ
(
η
)
=
(
T
-
T
∞
)
(
T
w
-
T
∞
)
,
ϕ
(
η
)
=
(
C
-
C
∞
)
(
C
w
-
C
∞
)
,
η
=
(
U
0
2
ν
L
)
1
/
2
e
(
x
+
y
)
/
2
L
z
,
where primes denote differentiation with respect to η. Next we take
(10)
w
w
(
x
,
y
)
=
-
(
ν
U
0
2
L
)
1
/
2
e
(
x
+
y
)
/
2
L
S
,
where S is the surface mass transfer parameter with S
>
0 for suction and S
<
0 for injection. Substituting the similarity variables (9) into (1) to (6), it is found that the continuity equation (1) is automatically satisfied, and (2) to (6) are reduced to the following ordinary (similarity) differential equations:
(11)
f
′′′
+
(
f
+
g
)
f
′′
-
2
(
f
′
+
g
′
)
f
′
=
0
,
g
′′′
+
(
f
+
g
)
g
′′
-
2
(
f
′
+
g
′
)
g
′
=
0
,
1
Pr
θ
′′
+
(
f
+
g
)
θ
′
+
Nb
ϕ
′
θ
′
+
Nt
θ
′
2
=
0
,
ϕ
′′
+
Le
(
f
+
g
)
ϕ
′
+
Nt
Nb
θ
′′
=
0
subject to the boundary conditions
(12)
f
(
0
)
=
S
,
g
(
0
)
=
0
,
f
′
(
0
)
=
λ
1
,
g
′
(
0
)
=
λ
2
,
θ
(
0
)
=
1
,
Nt
ϕ
′
(
0
)
+
Nb
θ
′
(
0
)
=
0
f
′
(
η
)
⟶
0
,
g
′
(
η
)
⟶
0
,
θ
(
η
)
⟶
0
,
ϕ
(
η
)
⟶
0
as
η
⟶
∞
,
where Pr is the Prandtl number, Le is the Lewis number, Nb is the Brownian motion parameter, and Nt is the thermophoresis parameter, which are defined as follows:
(13)
Pr
=
ν
α
,
Le
=
ν
D
B
,
Nb
=
β
D
B
(
C
w
-
C
∞
)
ν
,
Nt
=
β
D
T
(
T
w
-
T
∞
)
T
∞
ν
.
The physical quantities of practical interest are the local skin friction coefficients, C
f
x and C
f
y, and the local Nusselt number Nu
x, which are defined as follows:
(14)
C
f
x
=
2
τ
w
x
ρ
U
w
2
,
C
f
y
=
2
τ
w
y
ρ
U
w
2
,
Nu
x
=
q
w
k
(
T
w
-
T
∞
)
,
where τ
w
x and τ
w
y are the shear stresses in the x- and y-directions of the stretching/shrinking sheet and q
w is the heat flux from the surface of the sheet, which are given by
(15)
τ
w
x
=
μ
(
∂
u
∂
z
)
z
=
0
,
τ
w
y
=
μ
(
∂
v
∂
z
)
z
=
0
,
q
w
=
-
k
(
∂
T
∂
z
)
z
=
0
.
Substituting (9) into (14) and using (15), we obtain
(16)
R
e
x
1
/
2
C
f
x
=
2
f
′′
(
0
)
,
R
e
y
1
/
2
C
f
y
=
2
g
′′
(
0
)
,
R
e
x
-
1
/
2
Nu
x
=
-
θ
′
(
0
)
,
where R
e
x
=
U
w
L
/
ν and R
e
y
=
V
w
L
/
ν are the local Reynolds numbers.


## 3. Results and Discussions

The system of ordinary differential equations (11) subject to the boundary conditions (12) was solved numerically using the package bvp4c in MATLAB for different values of parameters: the stretching/shrinking parameter in x-direction λ
1, suction S, Brownian motion parameter Nb, thermophoresis parameter Nt, and Lewis number Le. We fixed the Prandtl number to be equal to 6.8 and the stretching/shrinking parameter in the y-direction λ
2 to be 1 (λ
2
=
1) throughout the paper. The relative tolerance is set to 10−10 and the boundary conditions (12) at η
=
∞ are replaced by η
=
10. This choice is sufficient for the velocity and the temperature profiles to reach the far field boundary conditions asymptotically.In this paper, we intend to study the three-dimensional flow and heat transfer of a nanofluid over a stretching/shrinking sheet. The analysis shows that the existence of solution depends on the suction parameterS and the stretching/shrinking parameter λ
1. Figures 2 and 3 show that the skin friction coefficient in the x-direction and the y-direction, respectively, decreases as λ
1 increases. From these figures, we can see that dual solutions exist for the problem. However, based on the previous studies [27, 32, 33], only the first solution is physically realizable and thus relevant to that studies. It is portrayed in Figures 2 and 3 that unique solution exists for λ
1
≥
-
1 and λ
1
=
λ
c, where λ
c is the critical values of λ
1. Furthermore, it is seen that the range of λ
1, where solutions exist, increases as S increases, as shown in Table 1. In addition, in Figure 2, it is shown that when the sheet is shrunk in the x-direction, the skin friction coefficient parallel to the direction increases as S increases. However, the skin friction coefficient in the y-direction decreases with increasing S as illustrated in Figure 3. Moreover, it is interesting to note that the shear stress in the x-direction is prominently higher than the shear stress in the y-direction.Table 1
Values ofλ
1
c.
S
λ
c
2.0
−1.7785
2.5
−2.2164
3.0
−2.7157Figure 2
Variation of the skin friction coefficient inx-direction with λ
1 for different values of S when λ
2
=
1.Figure 3
Variation of the skin friction coefficient iny-direction with λ
1 for different values of S when λ
2
=
1.Figure4 shows that the local Nusselt number increases with λ
1. However, the local Nusselt number decreases as thermophoresis parameter increases. This phenomenon may be caused by the thermal boundary layer that thickens as the thermophoresis parameter is increased. As opposed to this occurrence, the thermal boundary layer becomes thinner as the Brownian motion parameter increases. This leads to the increase of the local Nusselt number as Brownian motion parameter increases as shown in Table 2. The table also shows that the Lewis number lowers the local Nusselt number.Table 2
The local Nusselt number for different Nb and Le whenλ
1 = 1, Pr = 6.8, and Nt = 0.1.
Nb
Le
Nu
x
Re
x
-
1
/
2
0.10
5
11.689893
0.12
14.032429
0.14
19.445658
0.10
5
11.689893
8
11.260672
10
11.105908Figure 4
Variation of the local Nusselt number withλ
1 for different values of Nt when λ
2
=
1, Pr = 6.8, Nb = 0.1, and Le = 5.Figures5–7 show the velocity profiles for the flow in the x- and y-directions for different values S and λ
1. These profiles show that the far field boundary conditions are satisfied which validates the numerical result. Furthermore, these profiles also support the existence of dual solutions. The effect of S on both f
′
(
η
) and g
′
(
η
) is shown in Figures 5 and 6, respectively. From the two figures, it is noted that while S increases the velocity f
′
(
η
), it decreases the velocity g
′
(
η
). Figure 7 then shows the effect of λ on the velocity profiles f
′
(
η
) and g
′
(
η
). Increasing the stretching parameter in the x-direction causes f
′
(
η
) to increase. On the other hand, the velocity g
′
(
η
) is consistently lower for higher λ
1 although it is seen that the changes are minuscule.Figure 5
Velocity profiles inx-direction for different values of S when λ
1
=
-
1.7 and λ
2
=
1.Figure 6
Velocity profiles iny-direction for different values of S when λ
1
=
-
1.7 and λ
2
=
1.Figure 7
Velocity profiles for different values ofλ
1 when λ
2
=
1.


## 4. Conclusions

The three-dimensional flow and heat transfer of a nanofluid over a stretching/shrinking sheet was investigated numerically. The effects of various parameters on the skin friction coefficient and the local Nusselt number were discussed. The results showed that suction parameter increases the solution domain. Furthermore, as the sheet is shrunk in thex-direction, suction increases the skin friction coefficient in the same direction while decreasing the skin friction coefficient in the y-direction. As thermophoresis parameter and Lewis number increase, the local Nusselt number decreases. On the other hand, the local Nusselt number increases as Brownian motion parameter increases.



---

*Source: 517273-2014-08-24.xml*