Numerical Study of Boundary Layer Flow and Heat Transfer of Oldroyd-B Nanofluid towards a Stretching Sheet

In the present article, we considered two-dimensional steady incompressible Oldroyd-B nanofluid flow past a stretching sheet. Using appropriate similarity variables, the partial differential equations are transformed to ordinary (similarity) equations, which are then solved numerically. The effects of various parameters, namely, Deborah numbers and , Prandtl parameter , Brownian motion , thermophoresis parameter and Lewis number , on flow and heat transfer are investigated. To see the validity of the present results, we have made the comparison of present results with the existing literature.


Introduction
The flow over a stretching sheet has been premeditated because of its numerous industrial applications such as industrialized of polymer sheet, filaments and wires. Through the mechanized process, the stirring sheet is assumed to extend on its own plane and the protracted surface interacts with ambient fluid both impulsively and thermally. Only Navier Stokes equations are deficient to explain the rheological properties of fluids. Therefore, rheological non-Newtonian fluid models have been proposed to overcome this deficiency. Sakiadis [1] was the first who discussed the boundary layer flow over a stretching surface. He discussed numerical solutions of laminar boundary-layer behavior on a moving continuous flat surface. Experimental and analytical behavior of this problem was presented by Tsou et al. [2] to show that such a flow is physically possible by validating Sakiadis [1] work. Crane [3] extended the work of Sakiadis [1] for both linear and exponentially stretching sheet considering steady twodimensional viscous flow. Free convection on a vertical stretching surface was discussed by Wang [4]. Heat transfer analysis over an exponentially stretching continuous surface with suction was presented by Elbashbeshy [5]. He obtained similarity solutions for the laminar boundary layer equations describing heat and flow in a quiescent fluid driven by an exponentially stretching surface subject to suction. Viscoelastic MHD flow heat and mass transfer over a stretching sheet with dissipation of energy and stress work was discussed by Khan et al. [6]. Ishak et al. [7] studied heat transfer over a stretching surface with variable heat flux in micropolar fluids. Nadeem et al. [8] coated boundary layer flow of a Jeffrey fluid over an exponentially stretching surface with radiation effects. Recently in another article Nadeem et al. [9] investigated the magnetohydrodynamic (MHD) boundary layer flow of a Casson fluid over an exponentially permeable shrinking sheet.
The term ''Nanofluids'' is used for the fluids having suspension of nano-sized metallic or non-metallic particles. The main idea of using nanoparticles is to enhance the thermal properties of a base fluid. Invokement of nanofluids with improved heat distinctiveness can be noteworthy in stipulations of more competent cooling systems, consequential in higher productivity and energy savings. Several prospective applications for nanofluids are heat exchangers, radiators for engines, process cooling systems, microelectronics, etc. Choi [10] was the first who have made the analysis on nanoparticles in 1995. Xuan and Roetzel [11] presented cautiously the flow of a nanofluid in a tube using a dispersal replica. Heat transfer enhancement in a two-dimensional flow utilizing nanofluids is presented by Khanafer et al. [12]. They discussed the problem physically for various flow parameters. The Cheng-Minkowycz problem of natural convection past a vertical plate, in a porous medium saturated by a nanofluid is studied analytically by Nield and Kuznetsov [13]. The use of nanofluid model incorporates the effects of Brownian motion and thermophoresis parameter. The natural convective boundary layer flow of a nanofluid over a vertical plate is studied analytically by Kuznetsov and Nield [14]. They found that the reduced Nusselt number is a decreasing function of thermophoresis number and Brownian motion number. The boundary-layer flow and heat transfer in a viscous fluid containing metallic nanoparticles over a nonlinear stretching sheet are analyzed by Hamad and Ferdows [15]. They studied different types of nanoparticles and found that the behavior of the fluid flow changes with the change of the nanoparticles type. Numerous recent studies on nanofluids can be found in Refs. [16][17][18][19][20][21][22][23][24][25].
Main objective of the present article is to discuss the Oldroyd B nanofluid flow model over a stretching sheet. Mathematical model of the proposed study has been constructed after applying the boundary layer approach. Then, invoking the similarity transformation, we reduce the system of nonlinear partial differential equations into the system of nonlinear ordinary differential equations. The reduced couple nonlinear ODEs are solved numerically. Excellent comparison of the present approach has presented with the previous literature. The effects of various flow controlling parameters on the velocity, temperature and mass fraction function profiles are discussed. Moreover, variation of the local Nusselt and Sherwood number for various nanoparticles parameters has been constructed. The formulation of the paper is designed as follow. The problem formulation is presented in section two. The numerical solutions graphically with physical interpretation are incorporated in section three. Section four contains the summary of the whole analysis.

Problem Formulation
Consider two-dimensional steady incompressible Oldroyd-B fluid past a stretching sheet. In addition, nanoparticles effects are saturated, while sheet is stretching along the plane y~0. The flow is assumed to be confined to yw0. Here we assumed that the sheet is stretched with the linear velocity u(x)~ax, where aw0 is constant and x{axis is measured along the stretching surface. The boundary layer equations of Oldroyd-B fluid model along with the thermal energy and nanoparticles equations for nanofluids are where u and v denote the respective velocities in the x{ and y{ directions respectively, r f is the density of the base fluid, n is the kinematic viscosity of the fluid, s is the electrical conductivity, L 1 and L 2 are the relaxation and retardation times, a is the thermal diffusivity, T the fluid temperature, C the nanoparticle fraction, Introducing the following similarity transformations where the stream function y is define as u~L y Ly and v~{ Ly Lx .
Making use of Eq. (6), Equation of continuity is identically satisfied and Eqs.
(2) to (4) along with (5) take the following form where q w and q m are the heat flux and mass flux, respectively.
Dimensionless form of Eq. (13) take the form where Re x~uw (x)x=n is local Reynolds number based on the stretching velocity u w (x).

Results and Discussion
The nonlinear coupled ordinary differential equations (7)-(9) subject to the boundary conditions (10)- (12) have been solved numerically using the fourth-fifth order Runge-Kutta-Fehlberg method. Figs. 1, 2, 3, 4, 5, and 6 illustrate the behavior of emerging parameters such relaxation time constant b 1 , retardation time constant b 2 , Prandtl parameter Pr, Brownian parameter N b , thermophoresis parameter N t and Lewis number L e on velocity profile f '(g), temperature profile h(g) and mass fraction function w(g). Fig. 1, depicts the variation of b 1 on f '(g), h(g) and w(g).   Since b 1 is a function of relaxation time L 1 and due to viscoelastic properties of fluid it always resist the motion of the fluid. As a result, the velocity profile f '(g) and boundary layer thickness are decreasing function of b 1 . On the other hand, both temperature profile h(g) and mass fraction function w(g) increases with an increase in Deborah number b 1 (see Fig. 1). Physical behavior of Fig. 2 is due to an increase in retardation time of any material enhances the flow. Consequently, with an increase of b 2 velocity profile increases and both temperature and mass fraction function decreases (see Fig. 2). Thus, it concluded that b 1 and b 2 have opposite results on f '(g), h(g) and w(g) due to relaxation and retardation times, respectively (see Fig. 1 and 2).
Physically it is observed that an increase in the elastic parameter, the resistance to fluid flow will increase. Table 1 illustrates an excellent agreement of the present results with Khan and Pop [17] in the absence of non-Newtonian parameters b 1 and b 2 . As expected, it is found from Fig. 3, that both temperature and nanoparticle concentration profiles exert the decreasing behavior with the influence of Pr. Fig. 4 shows that both temperature and nanoparticle concentration have the same behavior when it is compared with Fig. 3 for higher values of L e . Consequently, boundary layer thickness decreases indefinitely with an increase in Pr. Effects of Brownian motion and thermophoresis parameters on temperature profile h(g) and mass fraction function w(g) are shown in Figs. 5 and 6. It is observed that for higher values of both N b and N t , the temperature profile rises. On the other hand Fig. 5, shows opposite behavior for mass fraction function when it is compare with Fig. 6, for increasing values of both N b and Nt. In the absence of both nanoparticles and non-Newtonian effects there is an excellent agreement of the present results with Wang [4] (see Table 2). The effects of elastic parameter, Prandtl parameter, Brownian parameter, thermophoresis parameter and Lewis number on the Nusselt number and Sherwood number are presented in Figs. 7, 8, 9, and 10. It is seen from Fig. 7, 8 and Table 3 that the Nusselt number decreases with increasing N t for both cases when Pr is less or greater than L e for N b~0 :3,0:5,0:7. Figs. 9 and 10 and Table 3 show the variation in dimensionless mass transfer rates vs N t parameter for the selected values of other parameters. The dimensionless mass transfer rates decrease with the increase in N t . Finally, high Prandtl fluid has a low thermal conductivity reducing conduction which results in an increase in the heat transfer rate at the surface of sheet.

Conclusions
In this study we have presented the Oldroyd-B fluid model for nanofluid over a stretching sheet. The effects of elastic parameter, Brownian motion and thermophoresis parameters on flow and heat transfer are discussed numerically. The main results of present analysis are listed below.
N Effects of b 1 and b 2 have opposite behavior for velocity, temperature and mass fraction function. These phenomena    x Nu x and the local Sherwood number Re {1=2 x Sh x in the presence of nanoparticle with b 1~b2~0 :3,L e~1 and Pr = 6. N Both temperature and mass fraction function give same behavior for Pr and L e . Since Pr is the ratio of kinematic to dynamic viscosity. Indeed for higher values of Pr, temperature profile remains under control.
N Effects of N b and N t for temperature profile are similar. Since both N b and N t causes to enhance the temperature.
N Effects of N b and N t for mass fraction function are opposite.
Mathematically, it is seen that both N b and N t appeared in the function in Eqn. (9). Consequently, behavior of mass fraction function profile will be opposite for various values of both N b and N t . Author Contributions