# Effect of Velocity Slip Boundary Condition on the Flow and Heat Transfer of Cu-Water and TiO2-Water Nanofluids in the Presence of a Magnetic Field

Authors: Abdelhalim Ebaid; Fahd Al Mutairi; S. M. Khaled
Journal: Advances in Mathematical Physics
 (2014)

---

## Abstract

In nanofluid mechanics, it has been proven recently that the no slip condition at the boundary is no longer valid which is the reason that we consider the effect of such slip condition on the flow and heat transfer of two types of nanofluids. The present paper considers the effect of the velocity slip condition on the flow and heat transfer of the Cu-water and the TiO2-water nanofluids over stretching/shrinking sheets in the presence of a magnetic field. The exact expression for the fluid velocity is obtained in terms of the exponential function, while an effective analytical procedure is suggested and successfully applied to obtain the exact temperature in terms of the generalized incomplete gamma function. It is found in this paper that the Cu-water nanofluid is slower than the TiO2-water nanofluid for both cases of the stretching/shrinking sheets. However, the temperature of the Cu-water nanofluid is always higher than the temperature of the TiO2-water nanofluid. In the case of shrinking sheet the dual solutions have been obtained at particular values of  the physical parameters. In addition, the effect of various physical parameters on such dual solutions is discussed through the graphs.


---

## Body

## 1. Introduction

Nanofluids are base fluids containing suspended nanoparticles. These nanoparticles are typically made of metals, oxides, carbides, or carbon nanotubes. The common base fluids include water, ethylene glycol, toluene, and oil. In the recent times, the study of nanofluid flow over a stretching sheet has gained much interest because of its numerous industrial applications. A few of these applications in engineering are polymer extrusion, wire drawing, continuous casting, manufacturing of foods and paper, glass fiber production, stretching of plastic films, glass blowing, manufacturing of plastic and rubber sheets, crystal growing, and continuous cooling and fiber spinning [1] and stretching of plastic films. In particular, during the manufacture of these sheets, the melt issues from a slit and is subsequently stretched to achieve the desired thickness. Hence, the final product with the desired characteristics strictly depends upon the stretching rate, the rate of cooling in the process, and the process of stretching [2]. Therefore, the choice of a proper cooling/heating liquid is essential as it has a direct impact on the rate of heat transfer. In addition, nanomaterials are of great interest in biology and medicine, owing to their numerous applications including DNA/protein detection, biomolecular regulators, cell imaging, and cancer cell diagnostics. Currently, the use of nanomaterials as drug delivery systems has become an emerging area in the field of nanomedicine. A wide variety of nanomaterials, such as nanotubes, nanorods, and nanoparticles, have been explored as carriers for delivering small molecule drugs, proteins, and genetic materials, exploiting their unique dimensions and specific physical and chemical properties [3].Choi [4] may be the first author to introduce the word nanofluid that represents the fluid in which nanoscale particles (diameter
<
50 nm) are suspended in the base fluid. The common nanoparticles that have been used are Aluminum, Copper, Silver, and Titanium or their oxides. Convective heat transfer of such nanoparticles in a base fluid is a topic of major interest both in applied sciences and engineering, where a very good review was presented by Saidur et al. [5] and Wang and Mujumdar [6, 7]. In this field, recent researches showed that nanoparticles changed the fluid characteristics because thermal conductivity of these particles was higher than convectional fluids. It was showed experimentally by Eastman et al. [8] and Xuan and Li [9] that the thermal conductivity of the base liquid can be enhanced by 10–20% with adding a small volumetric fraction of nanoparticles (usually
<
5
%). The enhanced thermal conductivity of nanofluids contributes to a remarkable improvement in the convective heat transfer coefficient. Recently, Hamad [10] studied the boundary-layer and heat transfer of nanofluids over impermeable isothermal stretching sheet for the metallic and metallic oxide nanoparticles. Noghrehabadi et al. [11] analyzed the heat transfer enhancement of water in the presence of SiO
2 nanoparticles over an isothermal stretching sheet. In addition, Noghrehabadi et al. [12] examined theoretically the flow and heat transfer of two types of nanofluids, namely Silver-water and Silicon Dioxide-water. They solved the governing equations by applying a combination of a symbolic power series and Padé approximation method. Further, Vajravelu et al. [13] studied the effect of variable viscosity on the flow and heat transfer of a viscous Ag-water and Cu-water nanofluids. They indicated that nanoparticle volume fraction is to increase the heat transfer and hence enhance the thermal boundary layer thickness. Moreover, an analysis has been performed by Alegn et al. [14] to study the heat transfer characteristics of steady two-dimensional boundary layer flow past a moving permeable flat plate in an alumina-water nanofluid. They have indicated that suction delays the boundary layer separation, while injection accelerates it.In [15], Das has investigated numerically the mixed convection stagnation point flow and heat transfer of Cu-water nanofluids impinging normally towards a shrinking sheet. Very recently, Aly and Ebaid [16] introduced a direct and effective approach to obtain the exact analytical solution for the nanoparticles-water flow over an isothermal stretching sheet with effect of the slip model. In particular, they examined and compared the effect of existence of five metallic and nonmetallic nanoparticles, namely, Silver, Copper, Alumina, Titania, and Silicon Dioxide, in a base of water. They introduced a very direct and effective approach to analytically obtain the exact solution which was found in a very good agreement with the previous published results.It has been shown in [17] that nanofluidic flow usually exhibits partial slip against the solid surface, which can be characterized by the so-called slip length (around 3.4–68 mm for different liquids). Accordingly, the authors in [18] discussed the effect of partial slip boundary condition on the flow and heat transfer of nanofluids past stretching sheet at constant wall temperature. Furthermore, the no-slip condition is no longer valid for fluid flows at the micro and nano scale and, instead, a certain degree of tangential slip must be allowed (see [19, 20]).Due to importance of the slip effect on modeling the boundary layer flows of nanofluids, a theoretical investigation has been introduced in this paper for studying the effect of the velocity slip boundary condition on the flow and heat transfer of Cu-water and TiO2-water nanofluids over a linearly stretching/shrinking sheets in the presence of magnetic field. Moreover, we introduce a direct and effective approach to analytically obtain the exact solutions for the governing nonlinear differential equation of the flow in terms of the exponential function while the temperature equation is solved exactly in terms of the generalized incomplete gamma function. The dual solutions have been also obtained in the case of the shrinking sheet. Further, we discuss the interested physical quantities, for example, the velocity and temperature, with comparison to results in literature. The structure of the paper is as follows. Description of the problem and basic equations are presented in Section 2. In Section 3, the exact analytical solutions of the flow and the temperature distribution are presented. In Section 4, we present the obtained numerical results. We finally conclude our paper in Section 5.


## 2. Basic Equations

Here, we consider a laminar steady two-dimensional incompressible viscous nanofluid over linearly semi-infinite stretching/shrinking sheets. A constant magnetic field is applied normally to the sheet. The fluid is a water-based nanofluid containing two different volume fractions of the Copper (Cu) and the Titania (TiO2) nanoparticles. The thermophysical properties of the base fluid and nanoparticles are given (see Table  1 in [21]). The governing boundary-layer equations of the considered nanofluid can be written in the following form [10]:
(1)
f
′′′
η
+
1
-
ϕ
2.5
1
-
ϕ
+
ϕ
ρ
s
ρ
f
·
f
η
f
′′
η
-
f
′
η
2
-
M
1
-
ϕ
2.5
f
′
η
=
0
,
1
Pr
k
n
f
k
f
1
1
-
ϕ
+
ϕ
ρ
C
p
s
/
ρ
C
p
f
·
θ
′′
η
+
f
η
θ
′
η
=
0
,
where primes denote the differentiation with respect to a similarity variable η, f
(
η
) and θ
(
η
) are the dimensionless stream function and temperature, respectively, ϕ is the solid volume fraction, ρ
f and ρ
s are the densities, (
ρ
C
p
)
f and (
ρ
C
p
)
s are the heat capacitances, M is the magnetic parameter, Pr is the Prandtl number, and k
n
f is the thermal conductivity defined as follows [22, 23]:
(2)
k
n
f
k
f
=
k
s
+
2
k
f
-
2
ϕ
k
f
-
k
s
k
s
+
2
k
f
+
ϕ
k
f
-
k
s
,
where k
f and k
s are the thermal conductivities, where (
)
f and (
)
s denote the basic fluid and solid fractions, respectively. The flow is subjected to the boundary conditions:
(3)
f
0
=
0
,
f
′
0
=
1
+
λ
f
′′
0
,
f
′
∞
=
0
,
Stretching sheet
,
f
0
=
s
>
0
,
f
′
0
=
-
1
+
λ
f
′′
0
,
f
′
∞
=
0
,
Shrinking sheet
(4)
θ
0
=
1
,
θ
∞
=
0
,
where λ is a slip parameter. Equations (1) can be rewritten as
(5)
f
′′′
η
+
α
f
η
f
′′
η
-
f
′
η
2
-
γ
f
′
η
=
0
,
(6)
τ
θ
′′
η
+
f
η
θ
′
η
=
0
,
where
(7)
α
=
1
-
ϕ
2.5
1
-
ϕ
+
ϕ
ρ
s
ρ
f
,
γ
=
M
1
-
ϕ
2.5
,
τ
=
1
Pr
k
n
f
k
f
1
1
-
ϕ
+
ϕ
ρ
C
p
s
/
ρ
C
p
f
.
In the case of the stretching sheet, it should be noted that the first part of the current paper reduces to Hamad [10] in the absence of the slip effect; that is, λ
=
0.


## 3. Exact Solutions

In this section, an effective analytical procedure shall be presented to exactly solve the system(5)-(6) with the boundary conditions (3)-(4) for both cases of the stretching and the shrinking sheets. In this regard, we observe from (5) that only unknown f
(
η
) is involved which is the reason that we begin with solving this f-differential equation in the next section.


### 3.1. Stretching Sheet

#### 3.1.1. Exact Solution for the Stream Functionf
(
η
)

In view of the boundary conditions given in(3), it is reasonable to assume f
(
η
) in the form:
(8)
f
η
=
a
+
b
e
-
β
η
,
in which the infinity boundary condition is automatically satisfied provided that β
>
0. On using (8) into (5), we have
(9)
β
b
β
2
-
α
a
β
-
γ
e
-
β
η
=
0
,
which leads to
(10)
β
2
-
α
a
β
-
γ
=
0
.
Applying the first two boundary conditions given in (2), we obtain
(11)
a
+
b
=
0
,
λ
b
β
2
+
b
β
+
1
=
0
.
Solving the system of the algebraic equations (10)-(11) for a, b, and β, we obtain
(12)
a
=
1
β
λ
β
+
1
,
b
=
-
1
β
λ
β
+
1
,
where β satisfies the following cubic algebraic equation:
(13)
λ
β
3
+
β
2
-
γ
λ
β
-
α
+
γ
=
0
.
It is also important here to refer to the fact that β is a positive real root of (13). Moreover, the exact solution given by (8) with a, b, and β given in (12) and (13) can be checked by direct substituting. Further, in the case of no slip, that is, λ
=
0, f
(
η
) reduces to that obtained by [10] (given by (14)-(15)). Regarding the heat transfer equation, we shall discuss in the next section an effective analytical procedure for obtaining the exact solution in terms of the generalized incomplete gamma function.


#### 3.1.2. Exact Solution of the Heat Transferθ
(
η
)

In order to solve the heat transfer equation, we first use(8) to rewrite (6) in terms of only θ
(
η
) as
(14)
τ
θ
′′
η
+
a
+
b
e
-
β
η
θ
′
η
=
0
.
On using separation technique and performing the integration once with respect to η from 0 to η, yields
(15)
θ
′
η
=
θ
′
0
e
-
a
/
τ
η
+
b
/
τ
β
e
-
β
η
-
b
/
τ
β
.
Integrating (15) once again with respect to η from 0 to η, we obtain
(16)
θ
η
=
θ
0
+
θ
′
0
∫
0
η
e
-
a
/
τ
σ
+
b
/
τ
β
e
-
β
σ
-
b
/
τ
β
d
σ
.
In view of the first condition in (4), we have
(17)
θ
η
=
1
+
θ
′
0
e
-
b
/
τ
β
∫
0
η
e
-
a
/
τ
σ
+
b
/
τ
β
e
-
β
σ
d
σ
.
The integration in the right-hand side can be analytically solved in terms of a well-known special function as declared by the following procedure. We first suppose that
(18)
μ
=
-
b
τ
β
e
-
β
σ
.
Accordingly
(19)
d
σ
=
-
d
μ
β
b
μ
,
e
-
σ
=
-
τ
β
b
1
/
β
μ
1
/
β
.
Inserting (18) and (19) into (17) gives
(20)
θ
η
=
1
+
1
β
b
-
τ
β
b
a
/
τ
β
θ
′
0
e
-
b
/
τ
β
·
∫
-
b
/
τ
β
e
-
β
η
-
b
/
τ
β
μ
a
/
τ
β
-
1
e
-
μ
d
μ
.
On using the definition of the generalized incomplete gamma function,
(21)
Γ
c
,
z
0
,
z
1
=
∫
z
0
z
1
μ
c
-
1
e
-
μ
d
μ
,
we can rewrite (20) as
(22)
θ
η
=
1
+
1
β
b
-
τ
β
b
a
/
τ
β
θ
′
0
e
-
b
/
τ
β
·
Γ
a
τ
β
,
-
b
τ
β
e
-
β
η
,
-
b
τ
β
.
Imposing the boundary condition θ
(
∞
)
=
0, we obtain the following explicit value for θ
′
(
0
):
(23)
θ
′
0
=
-
β
b
-
τ
β
/
b
-
a
/
τ
β
e
b
/
τ
β
Γ
a
/
τ
β
,
0
,
-
b
/
τ
β
.
Inserting this value of θ
′
(
0
) into (22), we get
(24)
θ
η
=
1
-
Γ
a
/
τ
β
,
-
b
/
τ
β
e
-
β
η
,
-
b
/
τ
β
Γ
a
/
τ
β
,
0
,
-
b
/
τ
β
.
This expression for θ
(
η
) can be further simplified by using the following identity for the generalized incomplete gamma function:
(25)
Γ
c
,
0
,
z
1
-
Γ
c
,
z
0
,
z
1
=
Γ
c
,
0
,
z
0
.
Hence, the exact solution for the temperature distribution θ
(
η
) is finally given by
(26)
θ
η
=
Γ
a
/
τ
β
,
0
,
-
b
/
τ
β
e
-
β
η
Γ
a
/
τ
β
,
0
,
-
b
/
τ
β
.


### 3.2. Shrinking Sheet

Following the same analysis as in the stretching sheet we can get the velocity exact solution as(27)
f
η
=
c
+
d
e
-
δ
η
,
in which the infinity boundary condition is automatically satisfied provided that δ
>
0, where c and d are given by
(28)
c
=
s
-
1
δ
λ
δ
+
1
,
d
=
1
δ
λ
δ
+
1
,
and δ should be at least one positive root of the cubic algebraic equation:
(29)
λ
δ
3
+
1
-
α
λ
s
δ
2
-
α
s
+
γ
λ
δ
+
α
-
γ
=
0
.
The temperature distribution is given by the same formula in the case of the stretching sheet with replacing a by c, replacing b by d, and replacing β by δ; that is,
(30)
θ
η
=
Γ
c
/
τ
δ
,
0
,
-
d
/
τ
δ
e
-
δ
η
Γ
c
/
τ
δ
,
0
,
-
d
/
τ
δ
.
Here, it is important to refer to that at particular values for the physical parameters we can get two positive roots for (29); this leads to the dual solutions for the velocity and the temperature as clarified latter in the next section.


## 3.1. Stretching Sheet

### 3.1.1. Exact Solution for the Stream Functionf
(
η
)

In view of the boundary conditions given in(3), it is reasonable to assume f
(
η
) in the form:
(8)
f
η
=
a
+
b
e
-
β
η
,
in which the infinity boundary condition is automatically satisfied provided that β
>
0. On using (8) into (5), we have
(9)
β
b
β
2
-
α
a
β
-
γ
e
-
β
η
=
0
,
which leads to
(10)
β
2
-
α
a
β
-
γ
=
0
.
Applying the first two boundary conditions given in (2), we obtain
(11)
a
+
b
=
0
,
λ
b
β
2
+
b
β
+
1
=
0
.
Solving the system of the algebraic equations (10)-(11) for a, b, and β, we obtain
(12)
a
=
1
β
λ
β
+
1
,
b
=
-
1
β
λ
β
+
1
,
where β satisfies the following cubic algebraic equation:
(13)
λ
β
3
+
β
2
-
γ
λ
β
-
α
+
γ
=
0
.
It is also important here to refer to the fact that β is a positive real root of (13). Moreover, the exact solution given by (8) with a, b, and β given in (12) and (13) can be checked by direct substituting. Further, in the case of no slip, that is, λ
=
0, f
(
η
) reduces to that obtained by [10] (given by (14)-(15)). Regarding the heat transfer equation, we shall discuss in the next section an effective analytical procedure for obtaining the exact solution in terms of the generalized incomplete gamma function.


### 3.1.2. Exact Solution of the Heat Transferθ
(
η
)

In order to solve the heat transfer equation, we first use(8) to rewrite (6) in terms of only θ
(
η
) as
(14)
τ
θ
′′
η
+
a
+
b
e
-
β
η
θ
′
η
=
0
.
On using separation technique and performing the integration once with respect to η from 0 to η, yields
(15)
θ
′
η
=
θ
′
0
e
-
a
/
τ
η
+
b
/
τ
β
e
-
β
η
-
b
/
τ
β
.
Integrating (15) once again with respect to η from 0 to η, we obtain
(16)
θ
η
=
θ
0
+
θ
′
0
∫
0
η
e
-
a
/
τ
σ
+
b
/
τ
β
e
-
β
σ
-
b
/
τ
β
d
σ
.
In view of the first condition in (4), we have
(17)
θ
η
=
1
+
θ
′
0
e
-
b
/
τ
β
∫
0
η
e
-
a
/
τ
σ
+
b
/
τ
β
e
-
β
σ
d
σ
.
The integration in the right-hand side can be analytically solved in terms of a well-known special function as declared by the following procedure. We first suppose that
(18)
μ
=
-
b
τ
β
e
-
β
σ
.
Accordingly
(19)
d
σ
=
-
d
μ
β
b
μ
,
e
-
σ
=
-
τ
β
b
1
/
β
μ
1
/
β
.
Inserting (18) and (19) into (17) gives
(20)
θ
η
=
1
+
1
β
b
-
τ
β
b
a
/
τ
β
θ
′
0
e
-
b
/
τ
β
·
∫
-
b
/
τ
β
e
-
β
η
-
b
/
τ
β
μ
a
/
τ
β
-
1
e
-
μ
d
μ
.
On using the definition of the generalized incomplete gamma function,
(21)
Γ
c
,
z
0
,
z
1
=
∫
z
0
z
1
μ
c
-
1
e
-
μ
d
μ
,
we can rewrite (20) as
(22)
θ
η
=
1
+
1
β
b
-
τ
β
b
a
/
τ
β
θ
′
0
e
-
b
/
τ
β
·
Γ
a
τ
β
,
-
b
τ
β
e
-
β
η
,
-
b
τ
β
.
Imposing the boundary condition θ
(
∞
)
=
0, we obtain the following explicit value for θ
′
(
0
):
(23)
θ
′
0
=
-
β
b
-
τ
β
/
b
-
a
/
τ
β
e
b
/
τ
β
Γ
a
/
τ
β
,
0
,
-
b
/
τ
β
.
Inserting this value of θ
′
(
0
) into (22), we get
(24)
θ
η
=
1
-
Γ
a
/
τ
β
,
-
b
/
τ
β
e
-
β
η
,
-
b
/
τ
β
Γ
a
/
τ
β
,
0
,
-
b
/
τ
β
.
This expression for θ
(
η
) can be further simplified by using the following identity for the generalized incomplete gamma function:
(25)
Γ
c
,
0
,
z
1
-
Γ
c
,
z
0
,
z
1
=
Γ
c
,
0
,
z
0
.
Hence, the exact solution for the temperature distribution θ
(
η
) is finally given by
(26)
θ
η
=
Γ
a
/
τ
β
,
0
,
-
b
/
τ
β
e
-
β
η
Γ
a
/
τ
β
,
0
,
-
b
/
τ
β
.


## 3.1.1. Exact Solution for the Stream Functionf
(
η
)

In view of the boundary conditions given in(3), it is reasonable to assume f
(
η
) in the form:
(8)
f
η
=
a
+
b
e
-
β
η
,
in which the infinity boundary condition is automatically satisfied provided that β
>
0. On using (8) into (5), we have
(9)
β
b
β
2
-
α
a
β
-
γ
e
-
β
η
=
0
,
which leads to
(10)
β
2
-
α
a
β
-
γ
=
0
.
Applying the first two boundary conditions given in (2), we obtain
(11)
a
+
b
=
0
,
λ
b
β
2
+
b
β
+
1
=
0
.
Solving the system of the algebraic equations (10)-(11) for a, b, and β, we obtain
(12)
a
=
1
β
λ
β
+
1
,
b
=
-
1
β
λ
β
+
1
,
where β satisfies the following cubic algebraic equation:
(13)
λ
β
3
+
β
2
-
γ
λ
β
-
α
+
γ
=
0
.
It is also important here to refer to the fact that β is a positive real root of (13). Moreover, the exact solution given by (8) with a, b, and β given in (12) and (13) can be checked by direct substituting. Further, in the case of no slip, that is, λ
=
0, f
(
η
) reduces to that obtained by [10] (given by (14)-(15)). Regarding the heat transfer equation, we shall discuss in the next section an effective analytical procedure for obtaining the exact solution in terms of the generalized incomplete gamma function.


## 3.1.2. Exact Solution of the Heat Transferθ
(
η
)

In order to solve the heat transfer equation, we first use(8) to rewrite (6) in terms of only θ
(
η
) as
(14)
τ
θ
′′
η
+
a
+
b
e
-
β
η
θ
′
η
=
0
.
On using separation technique and performing the integration once with respect to η from 0 to η, yields
(15)
θ
′
η
=
θ
′
0
e
-
a
/
τ
η
+
b
/
τ
β
e
-
β
η
-
b
/
τ
β
.
Integrating (15) once again with respect to η from 0 to η, we obtain
(16)
θ
η
=
θ
0
+
θ
′
0
∫
0
η
e
-
a
/
τ
σ
+
b
/
τ
β
e
-
β
σ
-
b
/
τ
β
d
σ
.
In view of the first condition in (4), we have
(17)
θ
η
=
1
+
θ
′
0
e
-
b
/
τ
β
∫
0
η
e
-
a
/
τ
σ
+
b
/
τ
β
e
-
β
σ
d
σ
.
The integration in the right-hand side can be analytically solved in terms of a well-known special function as declared by the following procedure. We first suppose that
(18)
μ
=
-
b
τ
β
e
-
β
σ
.
Accordingly
(19)
d
σ
=
-
d
μ
β
b
μ
,
e
-
σ
=
-
τ
β
b
1
/
β
μ
1
/
β
.
Inserting (18) and (19) into (17) gives
(20)
θ
η
=
1
+
1
β
b
-
τ
β
b
a
/
τ
β
θ
′
0
e
-
b
/
τ
β
·
∫
-
b
/
τ
β
e
-
β
η
-
b
/
τ
β
μ
a
/
τ
β
-
1
e
-
μ
d
μ
.
On using the definition of the generalized incomplete gamma function,
(21)
Γ
c
,
z
0
,
z
1
=
∫
z
0
z
1
μ
c
-
1
e
-
μ
d
μ
,
we can rewrite (20) as
(22)
θ
η
=
1
+
1
β
b
-
τ
β
b
a
/
τ
β
θ
′
0
e
-
b
/
τ
β
·
Γ
a
τ
β
,
-
b
τ
β
e
-
β
η
,
-
b
τ
β
.
Imposing the boundary condition θ
(
∞
)
=
0, we obtain the following explicit value for θ
′
(
0
):
(23)
θ
′
0
=
-
β
b
-
τ
β
/
b
-
a
/
τ
β
e
b
/
τ
β
Γ
a
/
τ
β
,
0
,
-
b
/
τ
β
.
Inserting this value of θ
′
(
0
) into (22), we get
(24)
θ
η
=
1
-
Γ
a
/
τ
β
,
-
b
/
τ
β
e
-
β
η
,
-
b
/
τ
β
Γ
a
/
τ
β
,
0
,
-
b
/
τ
β
.
This expression for θ
(
η
) can be further simplified by using the following identity for the generalized incomplete gamma function:
(25)
Γ
c
,
0
,
z
1
-
Γ
c
,
z
0
,
z
1
=
Γ
c
,
0
,
z
0
.
Hence, the exact solution for the temperature distribution θ
(
η
) is finally given by
(26)
θ
η
=
Γ
a
/
τ
β
,
0
,
-
b
/
τ
β
e
-
β
η
Γ
a
/
τ
β
,
0
,
-
b
/
τ
β
.


## 3.2. Shrinking Sheet

Following the same analysis as in the stretching sheet we can get the velocity exact solution as(27)
f
η
=
c
+
d
e
-
δ
η
,
in which the infinity boundary condition is automatically satisfied provided that δ
>
0, where c and d are given by
(28)
c
=
s
-
1
δ
λ
δ
+
1
,
d
=
1
δ
λ
δ
+
1
,
and δ should be at least one positive root of the cubic algebraic equation:
(29)
λ
δ
3
+
1
-
α
λ
s
δ
2
-
α
s
+
γ
λ
δ
+
α
-
γ
=
0
.
The temperature distribution is given by the same formula in the case of the stretching sheet with replacing a by c, replacing b by d, and replacing β by δ; that is,
(30)
θ
η
=
Γ
c
/
τ
δ
,
0
,
-
d
/
τ
δ
e
-
δ
η
Γ
c
/
τ
δ
,
0
,
-
d
/
τ
δ
.
Here, it is important to refer to that at particular values for the physical parameters we can get two positive roots for (29); this leads to the dual solutions for the velocity and the temperature as clarified latter in the next section.


## 4. Results and Discussion

In the previous section the exact analytical solutions for the slip effect on the flow and heat transfer of nanofluids past a stretching sheet/shrinking sheet (with suction) have been calculated. The present paper is an extension of the study made by the author [10] in which the slip boundary condition and the shrinking case were not considered. Therefore, we aim in this section not only to compare the current results (when λ
=
0) with those obtained in [10] at the same values of the physical parameters but also to study the effect of the slip parameter on the behaviour of the fluid velocity f
′
(
η
) and the temperature distribution θ
(
η
) in the case of shrinking sheet. For purpose of comparison, we present in Figure 1 the effect of the magnetic parameter M on the distribution of the velocity in the case of no slip, that is, at λ
=
0. The results showed good agreement with those obtained by [10] (Figure 1).Figure 1
Effect ofM on velocity distribution f
′
η for Pr
=
6.2 and ϕ
=
0.1.Regarding the temperature distribution, the exact solution forθ
(
η
) was obtained in [10] in the case of no slip in terms of the Whittaker M function and the hypergeometric function which are more complicated special functions than the generalized incomplete gamma function of the current analysis. It may be reasonable to refer here to that the analytical procedure followed in this paper is easier than the one considered in [10]. Figure 2 shows the effect of the magnetic parameter M on the temperature distribution in the case of no slip. The results also showed good agreement with those obtained by [10] (Figure 2). Although different special functions have been used in [10] to express θ
(
η
), the same results have been obtained with the impression that our procedure is easier. In order to confirm this point additional plots are depicted in Figure 3 for the effect of the volume fraction of the Cu-nanoparticles on temperature distribution, where also good agreement has been achieved.Figure 2
Effect ofM on temperature distribution θ
(
η
) for Pr
=
6.2 and ϕ
=
0.1.Figure 3
Effect ofϕ on temperature distribution θ
(
η
) for Pr
=
6.2 and M
=
1.The effect of the slip parameterλ on the velocity distribution f
′
(
η
) is depicted in Figures 4 and 5 for the Cu-water and the TiO2-water nanofluids, respectively, when Pr
=
6.2 (water), M
=
1, and ϕ
=
0.1. It is observed form these figures that the velocities of the Cu-water and the TiO2-water nanofluids decrease with increasing the slip parameter. It seems from Figures 4 and 5 that the Cu-water and the TiO
2-water nanofluids have the same magnitude of velocity at all the values of the slip parameter. In order to declare this point, these two figures are collected in Figure 6 which shows that f
′
(
η
)
|
Cu-water
<
f
′
(
η
)
|
TiO
2
-water for the small values of the slip parameter, that is, λ
≤
3, while f
′
(
η
)
|
Cu-water
≈
f
′
(
η
)
|
TiO
2
-water for the higher values of the slip parameter, that is, λ
>
3. The results reveal that the Cu-water nanofluid is slower than the TiO2-water nanofluid when the slip parameter is in the range 0
≤
λ
≤
3; otherwise, the two nanofluids have the same magnitude of velocity. Figures 7 and 8 display the effect of the slip parameter on the temperature distribution of the Cu-water and the TiO2-water nanofluids, respectively. It is clear from these figures that the temperature of the two nanofluids increases with the increase of the slip parameter. Figure 9 shows that the temperature of the Cu-water nanofluid is always higher than the temperature of the TiO2-water nanofluid whatever the values of the slip parameter.Figure 4
Effect ofλ on velocity distribution f
′
η for Pr
=
6.2, M
=
1, and ϕ
=
0.1.Figure 5
Effect ofλ on velocity distribution f
′
η for Pr
=
6.2, M
=
1, and ϕ
=
0.1.Figure 6
Effect ofλ on the velocity distribution f
′
η for the Cu-water and the TiO2-water nanofluids at Pr
=
6.2, M
=
1, and ϕ
=
0.1.Figure 7
Effect ofλ on temperature distribution θ
(
η
) for Pr
=
6.2 and M
=
1.Figure 8
Effect ofλ on temperature distribution θ
(
η
) for Pr
=
6.2 and M
=
1.Figure 9
Effect ofλ on the temperature distribution θ
(
η
) for the Cu-water and the TiO2-water nanofluids at Pr
=
6.2, M
=
1, and ϕ
=
0.1.In the case of the shrinking sheet it is has been shown in Figures10 to 14 that dual solutions have been obtained for the velocity and the temperature at the indicated parameters. In Figure 10(a) the effect of the volume fraction of the Cu-nanoparticles on the dual velocity is plotted. For the upper solution branch it is shown that f
′
(
η
) decreases with the increase of ϕ; however, a converse behaviour is observed for the lower solution branch of the velocity where slight increasing is observed. The situation for the effect of ϕ on the dual temperature is unique as depicted in Figure 10(b), where increasing ϕ increases the temperature for both the upper and the lower branches of the temperature. At different values of the selected parameters, it is shown in Figure 11(a) that both the upper and the lower branches of the dual velocity decreases with increasing the the volume fraction of the Cu-nanoparticles. However, the situation for the effect of ϕ on the dual temperature in Figure 11(b) is also unique and this agrees with the results obtained in Figure 10(b). The effect of the slip parameter λ on the dual velocity distributions f
′
(
η
) is depicted in Figures 12(a) and 12(b) for the Cu-water nanofluid. For the upper solution there is an increase in the velocity with increasing the slip parameter, while a converse behaviour is found for the lower solution branch. Figure 13 shows the effect of the small values of the slip parameter on the variation of the dual temperature solutions. It may be concluded from Figure 13 that the small variation in λ has no effect on the variation of the temperature. The effect of such slip parameter at higher values on the dual solutions of both the Cu-water and the TiO2-water nanofluids is depicted in Figures 14(a) and 14(b). The upper solution of both the two nanofluids increases with the increase of the slip parameter. For the lower solutions, the increase in the slip parameter decreases the velocities of the two nanofluids. One further observation here is that in the case of shrinking sheet the velocity of the Cu-water nanofluid is also slower than the TiO2-water nanofluid and these results are similar to those obtained in the case of the stretching sheet.(a) Dual velocity for several values ofϕ at s
=
1, λ
=
1, and M
=
1. (b) Dual temperature for several values of ϕ at s
=
1, λ
=
1, and M
=
1.
(a)
(b)(a) Dual velocity for several values ofϕ at s
=
2, λ
=
3, and M
=
1. (b) Dual temperature for several values of ϕ at s
=
2, λ
=
3, and M
=
1.
(a)
(b)(a) Upper velocity solution for several values ofλ at s
=
2, ϕ
=
0.1, and M
=
1.5. (b) Lower velocity solution for several values of λ at s
=
2, ϕ
=
0.1, and M
=
1.5.
(a)
(b)Figure 13
Dual temperature for several values ofλ at s
=
2, ϕ
=
0.1, and M
=
1.5.(a) Upper velocity solutions for several values ofλ for the Cu-water and the TiO2-water nanofluids at s
=
2, ϕ
=
0.1, and M
=
1. (b) Lower velocity solutions for several values of λ for the Cu-water and the TiO2-water nanofluids at s
=
2, ϕ
=
0.1, and M
=
1.
(a)
(b)


## 5. Conclusion

The effect of the velocity slip boundary condition on the flow and heat transfer of Cu-water and TiO2-water nanofluids in the presence of a magnetic field is investigated in this paper. The system of nonlinear differential equations governing the flow is exactly solved in terms of the exponential function for the fluid velocity and in terms of the generalized incomplete gamma function for the temperature distribution. In comparison with the previously published results, an easier analytical procedure has been suggested in this work to obtain the exact solutions. The main results of the present paper can be summarized as follows.(I)
In the case of the stretching sheet,
(i)
the velocities of the Cu-water and the TiO2-water nanofluids decrease with increasing the slip parameter λ;
(ii)
velocity
|
Cu-water
<
velocity
|
TiO
2
-water, for λ
≤
3, while velocity
|
Cu-water
≈
velocity
|
TiO
2
-water, when λ
>
3;
(iii)
the Cu-water nanofluid is slower than the TiO2-water nanofluid when the slip parameter is in the range 0
≤
λ
≤
3; otherwise, the two nanofluids have the same magnitude of velocity;
(iv)
the temperature of the two nanofluids increases with the increase of the slip parameter;
(v)
the temperature of the Cu-water nanofluid is always higher than the temperature of the TiO2-water nanofluid whatever the values of the slip parameter.
(II)
In the case of the shrinking sheet,
(i)
there are dual solutions for the velocity and the temperature at particular values of the physical parameters;
(ii)
the Cu-water nanofluid is slower than the TiO2-water nanofluid as in the case of the stretching sheet;
(iii)
the effect ofϕ on the dual temperature is unique, where increasing ϕ increases the temperature for both the upper and the lower branches of the temperature.



---

*Source: 538950-2014-08-05.xml*