# Numerical Study of the Influence of Heat Source on Stagnation Point Flow towards a Stretching Surface of a Jeffrey Nanoliquid

Authors: G. K. Ramesh
Journal: Journal of Engineering
 (2015)

---

## Abstract

An analysis is carried out to study the flow of Jeffrey fluid near a stagnation point towards a permeable stretching sheet. In particular, we investigate the effect of temperature dependent internal heat generation or absorption in the presence of nanoparticles. The governing system of partial differential equations is transformed into ordinary differential equations, which are then solved numerically using the fourth-fifth-order Runge-Kutta-Fehlberg method. Comparisons with previously published work on special cases of the problem are performed and found to be in excellent agreement. The results of the governing parametric study are shown graphically and the physical aspects of the problem are highlighted and discussed.


---

## Body

## 1. Introduction

Investigations on boundary layer flow and heat transfer of non-Newtonian fluids are increasing substantially due to the large number of practical applications in industrial and manufacturing processes. Examples of such applications are drilling muds, plastic polymers, optical fibers, hot rolling paper production, metal spinning, cooling of metallic plates in cooling baths, and many others. In the past, investigators proposed different non-Newtonian models because a single model cannot predict all the features of non-Newtonian materials. There is one subclass of non-Newtonian fluids known as Jeffrey fluid which has attracted much attention from the researchers in view of its simplicity. This fluid model is capable of describing the characteristics of relaxation and retardation times [1, 2]. The investigation of flow due to a stretching sheet has been intentional because of this flow’s various industrial applications, such as in the manufacturing of polymer sheets, filaments, and wires. During the manufacturing process, the moving sheet is assumed to stretch on its own plane, and the stretched surface interacts with the ambient fluid both mechanically and thermally. Stretching and shrinking can occur in a variety of materials each having a different strength, stretching transparency, and luster. Initially, Sakiadis [3] introduced the concept of a boundary layer flow over a stretching surface. Crane [4] modified the idea introduced by Sakiadis and extended this idea for both linear and exponentially stretching sheets.Flow in the neighborhood of a stagnation point in a plane was first studied by Hiemenz [5] and Mahapatra et al. [6–8] investigated the magnetohydrodynamic stagnation point flow towards a stretching sheet; they have shown that the velocity at a point decreases/increases with increase in the magnetic field when the free stream velocity is less/greater than the stretching velocity. Also they have studied the temperature distribution when the surface at constant temperature and constant heat fluxes. Further they have extended their work on power-law fluid and discussed the unique solutions of stagnation point flow of a power-law fluid towards a stretching surface. The study of heat source/sink effects on heat transfer is very important because its effects are crucial in controlling the heat transfer. Postelnicu et al. [9] examined the effect of variable viscosity on forced convection flow past a horizontal flat plate in a porous medium with internal heat generation, but in heat generation part they considered only space dependent heat source. Abel et al. [10] analyzed the non-Newtonian viscoelastic boundary layer flow of Walter’s liquid B past a stretching sheet, taking account of nonuniform heat source.Aforementioned studies were primarily concerned with the laminar flow of a clear fluid. In the recent past a new class of fluids, namely, nanofluids, has attracted the attention of the science and engineering community because of the many possible industrial applications of these fluids. Nanotechnology is an emerging science that is finding extensive use in industry due to the unique chemical and physical properties that the nanosized materials possess. These fluids are colloidal suspensions, typically metals, oxides, carbides, or carbon nanotubes in a base fluid. The term nanofluid was coined by Choi [11] in his seminal paper presented in 1995 at the ASME Winter Annual Meeting. It refers to fluids containing a dispersion of submicronic solid particles (nanoparticles) with typical length on the order of 1–50 nm. Kuznetsov and Nield [12] analytically studied the natural convective boundary layer flow of a nanofluid past a vertical plate. In a recent paper Khan and Pop [13] studied for the first time the problem of laminar fluid flow resulting from the stretching of a flat surface in a nanofluid. Mustafa et al. [14] investigated the stagnation point flow of a nanofluid towards a stretching surface using homotopy analysis method. Alsaedi et al. [15] examined the influence of heat generation/absorption on the stagnation point flow of nanofluid towards a linear stretching surface. Rahman et al. [16] investigated the dynamics of the natural convection boundary layer flow of water based nanofluids over a wedge in the presence of a transverse magnetic field with internal heat generation or absorption with the help of Matlab software. Nandy and Mahapatra [17] analyzed the effects of velocity slip and heat generation/absorption on magnetohydrodynamic stagnation point flow and heat transfer over a stretching/shrinking surface and obtained the solution numerically using fourth-order Runge-Kutta method with the help of shooting technique. Different from a stretching sheet, it was found that the solutions for a shrinking sheet are nonunique. Makinde et al. [18] studied the combined effects of buoyancy force, convective heating, Brownian motion, and thermophoresis on the stagnation point flow and heat transfer of an electrically conducting nanofluid towards a stretching sheet under the influence of magnetic field. Effect of magnetic field on stagnation point flow and heat transfer due to nanofluid towards a stretching sheet have been investigated by Ibrahim et al. [19]. Recently Ramesh and Gireesha [20] investigated the boundary layer flow of Maxwell fluid over a stretching surface with convective boundary condition in the presence of nanoparticles and heat source/sink effect. Nadeem et al. [21, 22] reported the numerical solutions of non-Newtonian nanofluid flow over a stretching sheet using the Jeffrey fluid model. Further they obtained the analytic solution for nonorthogonal stagnation point flow of a non-Newtonian nanofluid towards a stretching surface with heat transfer; here they use the second-grade model. Nadeem et al. [22] studied the natural convection boundary layer flow over a downward-pointing vertical cone in a porous medium saturated with a non-Newtonian nanofluid in the presence of heat generation or absorption and they used power-law model. Some interesting recent investigations related to the topic are presented in [23–25].In this paper, we study the behaviour of the stagnation point flow towards a stretching sheet with the effects of heat source/sink and suction in the presence of nanoparticles. Similarity transforms are presented for this problem, and nondimensionalized equations are addressed numerically. Graphical results for various values of the parameters are presented to gain thorough insight towards the physics of the problem. To the best of my knowledge, this problem has not been studied before.


## 2. Mathematical Analysis

Consider the flow of an incompressible Jeffrey fluid in the regiony
>
0 driven by a stretching surface located at y
=
0 with a fixed stagnation point at x
=
0. The stretching velocity U
w
(
x
) and the ambient fluid velocity U
∞
(
x
) are assumed to vary linearly from the stagnation point; that is, U
w
x
=
c
x and U
∞
x
=
b
x, where c and b are constant as shown in Figure 1. The problem under consideration is governed by the following boundary layer equations of Jeffrey fluid and nanoparticles and heat generation or absorption are given by Nadeem et al. [21]:(1)
∂
u
∂
x
+
∂
v
∂
y
=
0
,
(2)
1
+
λ
ν
u
∂
u
∂
x
+
v
∂
u
∂
y
=
∂
2
u
∂
y
2
+
U
∞
d
U
∞
d
x
+
λ
1
u
∂
3
u
∂
x
∂
y
2
+
v
∂
3
u
∂
y
3
-
∂
u
∂
x
∂
2
u
∂
y
2
+
∂
u
∂
y
∂
2
u
∂
x
∂
y
2
,
(3)
u
∂
T
∂
x
+
v
∂
T
∂
y
=
α
∂
2
T
∂
x
2
+
∂
2
T
∂
y
2
+
τ
D
B
∂
C
∂
x
∂
T
∂
x
+
∂
C
∂
y
∂
T
∂
y
∂
T
∂
x
2
+
∂
T
∂
y
2
h
h
h
h
h
h
h
+
D
T
T
∞
∂
T
∂
x
2
+
∂
T
∂
y
2
+
Q
0
ρ
f
c
p
T
-
T
∞
,
(4)
u
∂
C
∂
x
+
v
∂
C
∂
y
=
D
B
∂
2
C
∂
x
2
+
∂
2
C
∂
y
2
+
D
T
T
∞
∂
2
T
∂
x
2
+
∂
2
T
∂
y
2
,where u and v are the velocity components along the x- and y-axes, respectively. Further, α, ρ
f, ρ
p, ν, T, and T
∞ are, respectively, the thermal diffusivity, density of the base fluid, density of the particles, kinematic viscosity of the fluid, fluid temperature, and ambient fluid temperature. λ and λ
1 are ratios of relaxation to retardation times and retardation time, Q
0 is the dimensional heat generation/absorption coefficient, D
B is the Brownian diffusion coefficient, D
T is the thermophoresis diffusion coefficient, and c
p is the specific heat at constant pressure. Here τ is the ratio of the effective heat capacity of the nanoparticle material to the heat capacity of the ordinary fluid and C is the nanoparticle volume fraction.Figure 1
Schematic diagram of the flow.The associated boundary conditions for the present problem are(5)
u
=
U
w
x
,
v
=
-
V
w
x
,
T
=
T
w
,
C
=
C
w
at
y
=
0
,
u
⟶
U
∞
x
,
v
⟶
0
,
∂
u
∂
y
⟶
0
,
T
⟶
T
∞
,
C
⟶
C
∞
as
y
⟶
∞
,where U
w
x
=
c
x is the stretching sheet velocity, c
>
0; this is known as stretching rate. T
w and C
w are the temperature of fluid and nanoparticles fraction at wall, and C
∞ is the ambient nanoparticle volume fraction.The specific forms of the stretching velocity and the surface temperature and concentration are chosen to allow the coupled nonlinear partial differential equations (1)–(4) to be converted to a set of coupled, nonlinear ordinary differential equations by the similarity transformation(6)
η
=
U
w
ν
x
1
/
2
y
,
f
η
=
ψ
x
ν
U
w
1
/
2
,
θ
η
=
T
-
T
∞
T
w
-
T
∞
,
φ
η
=
C
-
C
∞
C
w
-
C
∞
,where η is the similarity variable; f, θ, and ϕ are the dimensionless stream function, temperature, and concentration, respectively. The velocity components u and v in (6) automatically satisfy the continuity equation (1). In terms of f
η, θ
η, and ϕ
η the momentum equation (2), energy equation (3), and concentration equation (4) can be written as(7)
f
η
η
η
+
β
f
η
η
2
-
f
f
η
η
η
η
+
(
1
+
λ
)
(
f
f
η
η
-
f
η
2
)
+
A
2
=
0
,
(8)
θ
η
η
+
Pr
⁡
⁡
f
θ
η
+
N
b
θ
η
ϕ
η
+
N
t
θ
η
2
+
S
θ
=
0
,
(9)
ϕ
η
η
+
Le
f
ϕ
η
+
N
t
N
b
θ
η
η
=
0
.Here f
η
(
η
) represents an ordinary derivative with respect to η and the corresponding boundary conditions in the nondimensional form are(10)
f
=
f
0
,
f
η
=
1
,
θ
=
1
,
ϕ
=
1
at
η
=
0
,
f
η
⟶
A
,
f
η
η
⟶
0
,
θ
⟶
0
,
ϕ
⟶
0
as
η
⟶
∞
.The dimensionless parameters in (7), (8), and (9) are A
=
c
/
b is the velocity ratio, β
=
c
λ
1 is the Deborah number, f
0
=
V
w
/
(
ν
c
)
1
/
2 is suction parameter, N
b
=
τ
D
B
φ
w
-
φ
∞
/
ν is the Brownian motion, N
t
=
τ
D
T
(
T
w
-
T
∞
)
/
ν
T
∞ is the thermophoresis parameter, and Le
=
ν
/
D
B is the Lewis number. Pr
⁡
=
ν
/
α is the Prandtl number; S
=
Q
0
/
c
ρ
f
c
p is the heat generation (
A
>
0
) or absorption (
A
<
0
) parameter.


## 3. Solutions for Some Special Cases

In the limiting case ofA
=
0, β
=
0, λ
=
0, and S
=
0 (i.e., in the absence of velocity ratio, Deborah number, ratio of relaxation to retardation times, and heat source/sink parameter ) our system of (7) and (9) reduces to those of Khan and Pop [13] (Newtonian nanofluid).In the presence of velocity ratio parameter, when there is no Deborah number, ratio of relaxation to retardation times, and heat source/sink parameter, system of (7) and (9) reduces to those of Mustafa et al. [14] (A
≠
0, β
=
0, λ
=
0, and S
=
0). Further in the absence of melting effect, heat source/sink, and nanoparticles volume fraction, the equations are similar to the ones studied by Hayat [27].In the absence of nanoparticles, the stagnation point flow of Jeffrey fluid and heat transfer problem degenerates. In this case, the approximate numerical solutions for the velocity field, temperature, and concentration fields are obtained.


## 4. Result and Discussion

The nonlinear coupled ordinary differential equations (7)–(9) that are subject to the boundary conditions (10) have been solved numerically using a fourth-fifth-order Runge-Kutta-Fehlberg method to obtain the missing values of f
′
′
0, -
θ
′
(
0
), and -
ϕ
′
0 for some values of the governing parameters, namely, the velocity ratio parameter A, heat source/sink parameter S, Deborah number β, ratio of relaxation to retardation times parameter λ, Brownian parameter N
b, thermophoresis parameter N
t, Lewis number Le, and Prandtl number Pr
⁡. In order to validate the numerical results obtained, we compare our results with those reported by Mahapatra and Gupta [7], Hayat et al. [28], and Ibrahim et al. [19] for various values of A as shown in Tables 1 and 2, and they are found to be in a favorable agreement.Table 1
Comparison of the values of skin friction coefficientf
′
′
(
0
) for various values of velocity ratio A.
A
Mahapatra and Gupta [7]
Hayat et al. [26]
Ibrahim et al. [19]
Present study
0.01
—
−0.9982
−0.9980
−0.99910
0.1
−0.9694
−0.9695
−0.9694
−0.96941
0.2
−0.9181
−0.9181
−0.9181
−0.91812
0.5
−0.6673
−0.6673
−0.6673
−0.66729
2.0
2.0175
2.0176
2.0175
2.01750
3.0
4.7293
4.7296
4.7292
4.72923Table 2
Comparison of the values of Nusselt number-
θ
′
(
0
) for various values of velocity ratio A.
Pr
A
Mahapatra and Gupta [7]
Hayat et al. [26]
Ibrahim et al. [19]
Present study
1
0.1
0.603
0.6021
0.6022
0.60223
0.2
0.625
0.6244
0.6245
0.62452
0.5
0.692
0.6924
0.6924
0.69249
1.5
0.1
0.777
0.7768
0.7768
0.77687
0.2
0.797
0.7971
0.7971
0.79712
0.5
0.863
0.8647
0.8648
0.86475The velocity profiles for different values ofA are presented in Figure 2. It is found that when the stretching velocity is less than the free stream velocity (
A
>
1
), the flow has a boundary layer structure, physically saying that the straining motion near the stagnation region increases so the acceleration of the external stream increases which leads to decrease in the thickness of the boundary layer with increase in A. When the stretching velocity c
x of the surface exceeds the free stream velocity b
x
(
A
<
1
) inverted boundary layer structure is formed and A
=
1; there is no boundary layer formation because the stretching velocity is equal to the free stream velocity. The temperature and concentration profiles for different values of A with other parameters being fixed are presented in Figures 3 and 4, respectively. It can be seen from these figures that both θ
(
η
) and ϕ
(
η
) decrease with increase in A.Figure 2
Velocity profile for different values ofA.Figure 3
Temperature profile for different values ofA.Figure 4
Concentration profile for different values ofA.Figure5 is plotted for velocity distribution for different values of β. When A
<
1 velocity increases with the increase of β, whereas when A
>
1 the velocity decreases with the increase of β. Similar effects can be seen in Figure 6 for various values of λ when A
>
1, while the opposite effect can be seen when A
<
1. This is because β is dependent upon λ (retardation time); physically larger retardation time of any material makes it less viscous resulting in an increase in its motion. It is worth mentioning that the liquid-like behavior is associated with small Deborah number. However, large Deborah number signifies the solid-like behavior.Figure 5
Velocity profile for different values ofβ.Figure 6
Velocity profile for different values ofλ.The graph of velocity profilef
′ versus η for different values of f
0 is plotted in Figure 7. It is found that for a fixed value of A
<
1 the velocity decreases with the increase of f
0. The velocity profiles tend asymptotically to the horizontal axis; the nondimensional velocities absorb maximum at the wall. It is a fact that suction stabilizes the boundary layer growth. At A
>
1 velocity increases with the increase of f
0. Figure 8 presents temperature profiles for different values of S. For S
>
0 (heat source), it can be observed that the thermal boundary layer generates the energy, and this causes the temperature in the thermal boundary layer to increase with increase in S, whereas S
<
0 (heat sink) leads to a decrease in the thermal boundary layer.Figure 7
Velocity profile for different values off
0.Figure 8
Temperature profile for different values ofS.Figures9 and 10 illustrate the variation of θ
,
ϕ with η for various values of N
b. It is found that the increase in the value of N
b is to increase the θ
(
η
) in the boundary layer, whereas in concentration boundary layer reduces as N
b increases which thereby enhances the nanoparticles concentration at the sheet. The graphs of N
t on the θ and ϕ profile are shown in Figures 11 and 12. From these plots, it is observed that the effect of increasing values of N
t is to increase the temperature and concentration profiles. Figure 13 displays the effect of Le on concentration profiles. It is noted that the concentration of the fluid decreases with increase of Le, physically due to the fact that mass transfer rate increases as Le increases. It also reveals that the concentration gradient at surface of the plate increases.Figure 9
Temperature profile for different values ofN
b.Figure 10
Concentration profile for different values ofN
b.Figure 11
Temperature profile for different values ofN
t.Figure 12
Concentration profile for different values ofN
t.Figure 13
Concentration profile for different values ofLe.Temperature and concentration profiles for the selected values ofPr
⁡ are plotted in Figures 14 and 15. The graph depicts that the temperature decreases when the values of Pr
⁡ increase. This is due to the fact that a higher Pr
⁡ fluid has relatively low thermal conductivity, which reduces conduction and thereby the thermal boundary layer thickness, and as a result temperature decreases. We note that N
b, N
t, Pr
⁡, and S have no influence on the flow field, which is clear from (7). From Table 3 one can see that the skin friction coefficient f
′
′
(
0
) is negative at A
<
1. Physically, negative value of f
′
′
(
0
) means the surface exerts a drag force on the fluid, and positive value means the opposite. This is not surprising since, in the present problem, we consider the case of a stretching sheet, which induces the flow. Negative value of θ
′
(
0
) means the heat flows from the fluid to the solid surface. This is not surprising since the fluid is hotter than the solid surface.Table 3
Values of skin friction coefficient-
f
′
′
(
0
), Nusselt number -
θ
′
(
0
), and Sherwood number -
ϕ
′
(
0
) for different values of the physical parameters with A
=
0.2.
β
λ
f
0
Le
N
b
N
t
Pr
S
-
f
′
′
(
0
)
-
ϕ
′
(
0
)
-
θ
′
(
0
)
0.05
0.05
0.02
1
0.1
0.1
3
0.1
0.9299
0.8437
0.8342
1
0.6836
0.8918
0.8785
2
0.5639
0.9156
0.8993
0.05
0.05
0.02
1
0.1
0.1
3
0.1
0.9299
0.8437
0.8342
1
1.3202
0.7776
0.7640
2
1.6397
0.7348
0.7077
0.05
0.05
0
1
0.1
0.1
3
0.1
0.9219
0.8354
0.8019
1
1.3479
1.1947
2.7517
2
1.7526
1.6380
4.9918
0.05
0.05
0.02
1
0.1
0.1
3
0.1
0.9299
0.8437
0.8342
2
0.9299
1.5778
0.7897
3
0.9299
2.1199
0.7695
0.05
0.05
0.02
1
0.1
0.1
3
0.1
0.9299
0.8437
0.8342
1
0.9299
1.2511
0.0874
2
0.9299
1.2436
0.0298
0.05
0.05
0.02
1
0.1
0.1
3
0.1
0.9299
0.8437
0.8342
1
0.9299
1.1108
0.3028
2
0.9299
3.1595
0.0844
0.05
0.05
0.02
1
0.1
0.1
1
0.1
0.9299
0.3671
0.5026
2
0.9299
0.6057
0.7138
3
0.9299
0.8437
0.8342
0.05
0.05
0.02
1
0.1
0.1
3
−0.1
0.9299
0.6449
1.0726
0
0.9299
0.7404
0.9586
0.1
0.9299
0.8437
0.8342Figure 14
Temperature profile for different values ofPr
⁡.Figure 15
Concentration profile for different values ofPr
⁡.


## 5. Conclusions

In the present investigation, the influence of the different parameters on the velocity, temperature, and concentration profiles is illustrated and discussed. The numerical results give a view towards understanding the response characteristics of the Jeffrey fluid in the presence of nanoparticles and heat generation/absorption. It is found that boundary layer is formed whenA
>
1; on the other hand inverted boundary layer is formed when A
<
1. Some results of thermal characteristics at the wall are usually analyzed from the numerical results and the same are documented in Tables 3 and 4. Analyzing this table, it reveals that the effect of increasing the values of f
0 and Pr
⁡ is to decrease the -
θ
′
(
0
) and -
ϕ
′
(
0
) and the effect of increasing the values of λ, Le, N
b, N
t, and S is to increase the -
θ
′
(
0
) and -
ϕ
′
(
0
). Also one can observe that there is no change in -
f
′′
(
0
), when Le, N
b, N
t, Pr
⁡, and S vary.Table 4
Values of skin friction coefficientf
′
′
(
0
), Nusselt number -
θ
′
(
0
), and Sherwood number -
ϕ
′
(
0
) for different values of the physical parameters with A
=
2.0.
β
λ
f
0
Le
N
b
N
t
Pr
S
f
′
′
(
0
)
-
ϕ
′
(
0
)
-
θ
′
(
0
)
0.05
0.05
0.02
1
0.1
0.1
3
0.1
1.9625
1.2545
1.1880
1
1.7009
1.2204
1.1565
2
1.4597
1.2022
1.1385
0.05
0.05
0.02
1
0.1
0.1
3
0.1
1.9625
1.2545
1.1880
1
1.0659
1.1603
1.1056
2
0.4741
1.0993
1.0491
0.05
0.05
0
1
0.1
0.1
3
0.1
1.9530
1.2453
1.1579
1
2.3962
1.5971
2.9677
2
2.7461
1.9483
5.1165
0.05
0.05
0.02
1
0.1
0.1
3
0.1
1.9625
1.2545
1.1880
2
1.9625
2.0475
1.1350
3
1.9625
2.6191
1.1091
0.05
0.05
0.02
1
0.1
0.1
3
0.1
1.9625
1.2545
1.1880
1
1.9625
1.6833
0.1720
2
1.9625
1.6664
0.0108
0.05
0.05
0.02
1
0.1
0.1
3
0.1
1.9625
1.2545
1.1880
1
1.9625
3.0176
0.4412
2
1.9625
6.5682
0.1676
0.05
0.05
0.02
1
0.1
0.1
1
0.1
1.9625
0.6307
0.8591
2
1.9625
0.9589
1.0804
3
1.9625
1.2545
1.1880
0.05
0.05
0.02
1
0.1
0.1
3
−0.1
1.9625
1.1225
1.3437
0
1.9625
1.1876
1.2671
0.1
1.9625
1.2545
1.1880



---

*Source: 382061-2015-03-23.xml*