Magnetohydrodynamic Boundary Layer Flow of Nanofluid over an Exponentially Stretching Permeable Sheet

Abstract

A mathematical model of the steady boundary layer flow of nanofluid due to an exponentially permeable stretching sheet with external magnetic field is presented. In the model, the effects of Brownian motion and thermophoresis on heat transfer and nanoparticle volume friction are considered. Using shooting technique with fourth-order Runge-Kutta method the transformed equations are solved. The study reveals that the governing parameters, namely, the magnetic parameter, the wall mass transfer parameter, the Prandtl number, the Lewis number, Brownian motion parameter, and thermophoresis parameter, have major effects on the flow field, the heat transfer, and the nanoparticle volume fraction. The magnetic field makes enhancement in temperature and nanoparticle volume fraction, whereas the wall mass transfer through the porous sheet causes reduction of both. For the Brownian motion, the temperature increases and the nanoparticle volume fraction decreases. Heat transfer rate becomes low with increase of Lewis number. For thermophoresis effect, the thermal boundary layer thickness becomes larger.

1. Introduction

The term “nanofluid” was proposed by Choi [1], referring to dispersions of nanoparticles in the base fluids such as water, ethylene glycol, and propylene glycol. The thermal conductivity enhancement characteristic of nanofluid was observed by Masuda et al. [2]. Buongiorno [3] discussed the reasons behind the enhancement in heat transfer for nanofluid and he found that Brownian diffusion and thermophoresis are the main causes. Later, Nield and Kuznetsov [4] and Kuznetsov and Nield [5] investigated the natural convective boundary layer flow of a nanofluid employing Buongiorno model.

The study of boundary layer flow and heat transfer due to stretching surface has numerous applications in industry and technology, such as in polymer extrusion, drawing of copper wires, artificial fibers, paper production, hot rolling, wire drawing, glass fiber, metal extrusion and metal spinning, and continuous stretching of plastic films. Crane [6] first studied the boundary layer flow due to linearly stretching sheet. Many researchers [7–17] extended the work of Crane, whereas Magyari and Keller [18] considered the boundary layer flow and heat transfer due to an exponentially stretching sheet. The flow and heat transfer over an exponentially stretching surface were investigated by Elbashbeshy [19] taking wall mass suction. Khan and Sanjayanand [20] presented the boundary layer flow of viscoelastic fluid and heat transfer over an exponentially stretching sheet with viscous dissipation effect and Partha et al. [21] reported a similarity solution for mixed convection flow past an exponentially stretching surface. Ishak [22] studied the magnetohydrodynamic (MHD) boundary layer flow over an exponentially shrinking sheet in presence of thermal radiation. Bhattacharyya [23] discussed the boundary layer flow and heat transfer caused due to an exponentially shrinking sheet and Bhattacharyya and Pop [24] showed the effect of external magnetic field on the flow over an exponentially shrinking sheet. Recently, Bhattacharyya and Vajravelu [25] described the stagnation-point boundary layer flow due to exponentially shrinking sheet for Newtonian fluid and Bachok et al. [26] investigated the same problem for nanofluid.

The boundary layer flow of nanofluid past a linearly stretching sheet was first studied by Khan and Pop [27] introducing the model of Nield and Kuznetsov [4]. The boundary layer flow induced in a nanofluid due to a linearly stretching sheet with convective boundary condition was described by Makinde and Aziz [28]. Kandasamy et al. [29] investigated the MHD boundary layer flow of a nanofluid past a vertical stretching permeable surface with suction/injection. Mustafa et al. [30] reported the flow of a nanofluid near a stagnation-point towards a stretching surface. Rana and Bhargava [31] illustrated the steady, laminar boundary layer flow due to the nonlinear stretching of a flat surface in a nanofluid. Later, Hady et al. [32] analysed the boundary layer flow and heat transfer characteristics of a viscous nanofluid over a nonlinearly stretching sheet in the presence of thermal radiation and variable wall temperature. Makinde et al. [33] discussed the combined effects of buoyancy force and magnetic field on stagnation-point flow and heat transfer in a nanofluid flow towards a stretching sheet. A theoretical study of unsteady boundary layer flow of a nanofluid over a permeable stretching/shrinking sheet was reported by Bachok et al. [34]. Recently, Nadeem and Lee [35] obtained analytic solutions of boundary layer flow of nanofluid over an exponentially stretching surface using homotopy analysis method (HAM).

In the present paper, the MHD boundary layer nanofluid flow dynamics due to an exponentially shrinking sheet is investigated. The obtained ordinary differential equations by similarity transformations are solved numerically using shooting technique with fourth-order Runge-Kutta method. Then computed results are plotted in graphs and discussed in detail.

2. Mathematical Formulation

Consider the steady boundary layer flow of nanofluid over an exponentially stretching sheet in presence of a transverse magnetic field. The governing equations of motion and the energy equation may be written in usual notation as [27, 35] 𝜕𝑢+𝜕𝑥𝜕𝑣𝑢𝜕𝑦=0,𝜕𝑢𝜕𝑥+𝑣𝜕𝑢𝜕𝜕𝑦=𝜐²𝑢𝜕𝑦²−𝜎𝐵²𝜌_𝑓𝑢𝑢,𝜕𝑇𝜕𝑥+𝑣𝜕𝑇𝜕𝜕𝑦=𝛼²𝑇𝜕𝑦²+𝜌𝑐_𝑝𝜌𝑐_𝑓𝐷_𝐵𝜕𝑁𝜕𝑦𝜕𝑇+𝐷𝜕𝑦𝑇^𝑇_∞𝜕𝑇𝜕𝑦²,𝑢𝜕𝑁𝜕𝑥+𝑣𝜕𝑁𝜕𝑦=𝐷_𝐵𝜕²𝑁𝜕𝑦²+𝐷_𝑇𝑇_∞𝜕²𝑇𝜕𝑦², (1) where 𝑢 and 𝑣 are the velocity components in 𝑥- and 𝑦-directions, respectively, 𝜐 is the kinematic viscosity, 𝜌_𝑓 is the density of the base fluid, 𝑇 is the temperature, 𝑇_∞ is constant temperature of the fluid in the inviscid free stream, 𝛼 is the thermal conductivity, (𝜌𝑐)_𝑝 is the effective heat capacity of nanoparticles, (𝜌𝑐)_𝑓 is heat capacity of the base fluid, 𝑁 is nanoparticle volume fraction, 𝐷_𝐵 is the Brownian diffusion coefficient, and 𝐷_𝑇 is the thermophoretic diffusion coefficient. Here, the variable magnetic field 𝐵(𝑥) is taken in the form [22, 24] 𝐵(𝑥)=𝐵₀𝑥exp2𝐿, (2) where 𝐵₀ is a constant.

The boundary conditions are given by 𝑢=𝑈_𝑤(𝑥),𝑣=𝑣_𝑤at𝑦=0,𝑢⟶0as𝑦⟶∞,𝑇=𝑇_𝑤=𝑇_∞+𝑇₀𝑥exp2𝐿at𝑦=0,𝑇⟶𝑇_∞as𝑦⟶∞,𝑁=𝑁_𝑤=𝑁_∞+𝑁₀𝑥exp2𝐿at𝑦=0,𝑁⟶𝑁_∞as𝑦⟶∞, (3) where 𝑇_𝑤 is the variable temperature at the sheet with 𝑇₀ being a constant which measures the rate of temperature increase along the sheet, 𝑁_𝑤 is the variable wall nanoparticle volume fraction with 𝑁₀ being a constant, and 𝑁_∞ is constant nanoparticle volume fraction in free stream. The stretching velocity 𝑈_𝑤 is given by 𝑈_𝑤(𝑥)𝑥=𝑐exp𝐿, (4) where 𝑐>0 is stretching constant. A physical model with the coordinate system of the problem is sketched in Figure 1.

Figure 1 

Physical model and coordinate system.

Here 𝑣_𝑤 is the variable wall mass transfer velocity and is given by 𝑣_𝑤=𝑣₀𝑥exp2𝐿, (5) where 𝑣₀ is a constant with 𝑣₀<0 for mass suction and 𝑣₀>0 for mass injection.

Now, we introduce the similarity transformations: 𝜓=√2𝜐𝐿𝑐𝑓𝜂exp𝑥𝜂=2𝐿,𝜃𝑇−𝑇_∞𝑇_𝑤−𝑇_∞,𝜑𝜂=𝑁−𝑁_∞𝑁_𝑤−𝑁_∞,𝜂=𝑦𝑐𝑥2𝜐𝐿exp,2𝐿 (6) where Ψ is the stream function with 𝑢=𝜕𝜓/𝜕𝑦 and 𝑣=−𝜕𝜓/𝜕𝑥 and 𝜂 is the similarity variable.

In view of relations in (6) we finally obtain the following self-similar equations: 𝑓^′′′+𝑓𝑓^′′−2𝑓^′2−𝑀𝑓^′𝜃=0,_′′+Pr𝑓𝜃^′−𝑓^′𝜃+Nb𝜃^′𝜑^′+Nt𝜃^′2𝜑=0,_′′+Le𝑓𝜑^′−𝑓^′𝜑+Nt𝜃Nb^′′=0, (7) where 𝑀=2𝜎𝐵²₀𝐿/𝑐𝜌 is the magnetic parameter, Pr=𝜐/𝛼 is the Prandtl number, and Le=𝜐/𝐷_𝐵 is the Lewis number. The dimensionless parameters Nb (Brownian motion parameter) and Nt (thermophoresis parameter) are defined as Nb=𝐷_𝐵𝜌𝑐_𝑝𝜌𝑐_𝑓𝑁_𝑤−𝑁_∞𝜐,𝐷Nt=^𝑇𝑇_∞𝜌𝑐_𝑝𝜌𝑐_𝑓𝑇_𝑤−𝑇_∞𝜐. (8)

The boundary conditions (3) reduce to the following forms: 𝑓𝜂=𝑆,𝑓^′𝜂𝑓=1at𝜂=0,_′𝜂𝜃𝜂𝜂𝜑𝜂𝜂⟶0as𝜂⟶∞,=1at𝜂=0,𝜃⟶0as𝜂⟶∞,=1at𝜂=0,𝜑⟶0as𝜂⟶∞, (9) where 𝑆=−𝑣₀/𝜐𝑐/2𝐿 is the wall mass transfer parameter. 𝑆>0(𝑣₀<0) corresponds to mass suction and 𝑆<0(𝑣₀>0) corresponds to mass injection.

The quantities of physical interest for this problem are the local skin friction coefficient 𝐶_𝑓, the local Nusselt number Nu_𝑥, and the local Sherwood number Sh_𝑥, which are, respectively, defined as [35] 𝐶_𝑓=𝜐𝑈²_𝑤𝑒^2𝑥/𝐿𝜕𝑢𝜕𝑦_𝑦=0,Nu_𝑥𝑥=−𝑇_𝑤−𝑇_∞𝜕𝑇𝜕𝑦_𝑦=0,Sh_𝑥𝑥=−𝑁_𝑤−𝑁_∞𝜕𝑁𝜕𝑦_𝑦=0, (10) that is, 2Re_𝑥𝐶_𝑓=𝑓^′′(0),Nu_𝑥2Re_𝑥=−𝑥𝜃2𝐿^′(0),Sh_𝑥2Re_𝑥=−𝑥𝜑2𝐿^′(0), (11) where Re_𝑥=𝑈_𝑤𝑥/𝜐 is the local Reynolds number.

3. Numerical Method for Solution

The highly nonlinear coupled ODEs (7) along with the boundary conditions (9) form a two-point boundary value problem (BVP) and those are solved using shooting method [36–41]. The following first-order system is set: 𝑓^′=𝑝,𝑝^′=𝑞,𝑞^′=2𝑝²𝜃−𝑓𝑞+𝑀𝑝,_′=𝑟,𝑟^′=−Pr𝑓𝑟−𝑝𝜃+Nb𝑟𝑧+Nt𝑟²,𝜑^′=𝑧,𝑧^′=−Le𝑓𝑧−𝑝𝜑−Nt𝑟Nb^′, (12) with the boundary conditions 𝑓(0)=𝑆,𝑝(0)=1,𝜃(0)=1,𝜑(0)=1. (13) The set of nonlinear first-order ordinary differential equations (12) with boundary conditions (13) have been solved by shooting method using the fourth-order Runge-Kutta algorithm with a systematic guessing of 𝑞(0), that is, 𝑓^′′(0), 𝑟(0), that is, 𝜃^′(0), and 𝑧(0), that is, 𝜑^′(0). The step size is taken as Δ𝜂=0.01 and the suitable finite value of 𝜂→∞, 𝜂_∞, is taken as 20 in all cases. The guess values𝑓^′′(0), 𝜃^′(0), and 𝜑^′(0) are adjusted using “secant method” to give better approximation for the solution. An asymptotic convergence criterion of 10⁻⁵ level for the boundary conditions𝑓^′(𝜂_∞)=0, 𝜃(𝜂_∞)=0, and 𝜑(𝜂_∞)=0 is taken in the computation.

4. Results and Discussion

The numerical solutions are obtained using the above numerical scheme for some values of the governing parameters, namely, the magnetic parameter (𝑀), the wall mass transfer parameter (𝑆), the Prandtl number (Pr), the Lewis number (Le), Brownian motion parameter (Nb), and thermophoresis parameter (Nt). Effects of 𝑀, 𝑆, Pr, Le, Nb, and Nt on the steady boundary layer flow, heat transfer, and nanoparticle volume fraction over exponentially stretching sheet in nanofluid are discussed in detail.

To ensure the numerical accuracy, the values 𝑓^′′(0) and 𝑓(∞) are compared with the results of Magyari and Keller [18] in Table 1 without magnetic field (𝑀=0) and with nonporous stretching sheet (𝑆=0) and those are found in excellent agreement. Thus, we are very much confident that the present results are accurate.

Table 1 

The comparison of values of 𝑓^′′(0) and 𝑓(∞) for 𝑀=0 and 𝑆=0.

The nondimensional velocity 𝑓^′(𝜂), temperature 𝜃(𝜂), and nanoparticle volume fraction 𝜑(𝜂) for various values of magnetic parameter are shown in Figure 2. The velocity reduces with the increase of magnetic parameter. The magnetic field opposes the transport process. Actually, the increase of 𝑀 leads to the increase of the Lorentz force arising because of interaction of magnetic and electric fields for the motion of an electrically conducting fluid, and the stronger Lorentz force produces much more resistance to the transport phenomena. On the other hand, the temperature and the nanoparticle volume fraction increase with 𝑀. The Lorentz force has the tendency to increase the temperature and nanoparticle volume fraction in nanofluid motion. Consequently, the thermal boundary layer thickness and nanoparticle volume fraction boundary layer thickness become thicker for stronger magnetic field.

(a)

(b)

(c)

Figure 2 

Effect of magnetic parameter 𝑀 on (a) velocity 𝑓′(𝜂), (b) temperature 𝜃(𝜂), and (c) nanoparticle volume fraction 𝜑(𝜂).

In Figure 3, the velocity, temperature, and nanoparticle volume fraction are presented for variation in wall mass transfer parameter 𝑆. With increasing values of the mass suction parameter (𝑆>0), the velocity𝑓^′(𝜂), temperature 𝜃(𝜂), and nanoparticle volume fraction 𝜑(𝜂) in the boundary layer region decrease, whereas, due to the increase of mass injection (𝑆<0), all those increase. Due to mass suction, the fluid is brought closer to the sheet and it thins velocity boundary layer thickness as well as the thermal and nanoparticle volume boundary layer thicknesses. Opposite effect is found for mass injection case; that is, the fluid is taken away from the sheet. Consequently, the velocity, thermal, and nanoparticle volume boundary layer thicknesses become broader.

(a)

(b)

(c)

Figure 3 

Effect of mass transfer parameter 𝑆 on (a) velocity 𝑓′(𝜂), (b) temperature 𝜃(𝜂), and (c) nanoparticle volume fraction 𝜑(𝜂).

The influences of the Prandtl number Pr and the Lewis number Le on the temperature and nanoparticle volume fraction are depicted in Figures 4 and 5, respectively. The increment of Prandtl number results in major effects on temperature as well as on nanoparticle volume fraction. The thermal boundary layer thickness reduces with Prandtl number and it happens due to decrease of thermal diffusivity for the increment of Prandtl number. The nanoparticle volume fraction exhibits overshoot near the sheet for higher values of Pr, though the nanoparticle volume boundary layer thickness reduces. Hence, with uniform thermophoretic particle deposition, for larger values of Prandtl number the nanoparticle volume fraction 𝜑(𝜂) is higher in the fluid adjacent to the sheet than the value at the wall. Very minor variation (initially increasing near the sheet and then decreasing away from the sheet) is observed in the temperature with the increase in the Lewis number. For large values of Le, the nanoparticle volume fraction significantly decreases and also the nanoparticle volume boundary layer thickness reduces. But, for smaller value of Le the overshoot is found in 𝜑(𝜂) near the sheet (Figure 5(b)). The Brownian diffusion effect becomes nominal for larger values of Lewis number and for which the nanoparticle volume boundary layer thickness decreases.

(a)

(b)

Figure 4 

Effect of Prandtl number Pr on (a) temperature 𝜃(𝜂) and (b) nanoparticle volume fraction 𝜑(𝜂).

(a)

(b)

Figure 5 

Effect of Lewis number Le on (a) temperature 𝜃(𝜂) and (b) nanoparticle volume fraction 𝜑(𝜂).

The effect of Brownian motion parameter Nb on the dimensionless temperature and the dimensionless nanoparticle volume fraction is plotted in Figure 6 (Figures 6(a) and 6(b), respectively). The figure reveals that the temperature of the fluid increases and the nanoparticle volume fraction decreases with increasing values of Nb. In nanofluid system, due to the presence of nanoparticles, the Brownian motion takes place and for the increase in Nb the Brownian motion is affected and consequently the heat transfer characteristic of the fluid changes. Also, when the value of Nb increases, the nanoparticle volume boundary layer thickness decreases.

(a)

(b)

Figure 6 

Effect of Brownian motion parameter Nb on (a) temperature 𝜃(𝜂) and (b) nanoparticle volume fraction 𝜑(𝜂).

Thermophoresis parameter Nt is a key parameter for analysing the temperature distributions and nanoparticles volume fraction in nanofluid flow. The effect of thermophoresis parameter Nt on the temperature 𝜃(𝜂) and the nanoparticle volume fraction 𝜑(𝜂) is presented in Figure 7. With increase of Nt, the temperature of the fluid increases. The nanoparticle volume fraction increases with Nt and the overshoot near the wall is found. Increase in Nt causes increment in the thermophoresis force which tends to move nanoparticles from hot to cold areas and consequently it increases the magnitude of temperature profiles and nanoparticle volume fraction profiles. Ultimately, the thickness of nanoparticle volume boundary layer becomes significantly large for slightly increased value of thermophoresis parameter.

(a)

(b)

Figure 7 

Effect of thermophoresis parameter Nt on (a) temperature 𝜃(𝜂) and (b) nanoparticle volume fraction 𝜑(𝜂).

We now discuss the variations of the physical quantities of engineering importance, that is, the local skin friction coefficient 𝐶_𝑓, the local Nusselt number Nu_𝑥, and the local Sherwood number Sh_𝑥 for different values of 𝑀, 𝑆, Pr, Le, Nb, and Nt. And for which the quantities −𝑓^′′(0), −𝜃^′(0), and −𝜑^′(0) related to local skin friction coefficient, the local Nusselt number and the local Sherwood number, respectively, are plotted in Figures 8, 9, 10, 11, 12, 13, 14, 15, and 16 for various values of parameters. The values of −𝑓^′′(0), −𝜃^′(0), and −𝜑^′(0) versus the mass transfer parameter 𝑆 are plotted in Figure 8 for different values magnetic parameter. For stronger magnetic field the value of −𝑓^′′(0), that is, the local skin friction coefficient, increases, whereas the values of −𝜃^′(0), that is, local Nusselt number, and −𝜑^′(0), that is, local Sherwood number, decreases. Also, the values of −𝜃^′(0) and −𝜑^′(0) are depicted against Brownian motion parameter Nb and thermophoresis parameter Nt for various 𝑀 in Figures 9 and 10. Similar results are also observed in those graphs. In physical viewpoint, it can be noticed that the powerful Lorentz force that arose in flow field for larger magnetic field reduces the values of local Nusselt number and local Sherwood number. On the other hand, the skin friction coefficient increases (decreases) with mass suction (injection) (Figure 8(a)). −𝜃^′(0) and −𝜑^′(0) are displayed against Nb and Nt for several 𝑆 in Figures 11 and 12. From all graphical results it should be noted that, similar to that of skin friction, the local Nusselt number and the local Sherwood number also decrease with injection and increase with suction. But for greater values of thermophoresis parameter Nt, the contrary effect of mass transfer parameter on the local Sherwood number is obtained, namely, the value of−𝜑^′(0) increases with mass injection and decreases with mass suction (Figure 12(b)). The effects of Prandtl number Pr and Lewis number Le on the local Nusselt number Nu_𝑥 and the local Sherwood number Sh_𝑥 are illustrated in Figures 13–16. Due to higher values of Prandtl number the local Nusselt number increases and local Sherwood number reduces. The heat transfer rate is strengthened when the Prandtl number increases, whereas it is worth noting that the Lewis number affects the Nusselt number (reduction) and the Sherwood number (increment) in completely reverse manner compared to that of Prandtl number. For small values of Lewis number, the Brownian diffusion effect is large and accordingly increased heat transfer rate is found (greater Nusselt number, Figure 15(a)). From Figures 9 to 16, the variations in local Nusselt number and the local Sherwood number for Brownian motion parameter Nb and thermophoresis parameter Nt can be understood. For high rate of Brownian motion, that is, for increasing Nb, the values of −𝜃^′(0) reduce and those of −𝜑^′(0) increase. But the increase in thermophoresis parameter Nt causes reduction in both, the local Nusselt number and the local Sherwood number.

(a)

(b)

(c)

Figure 8 

Variations of (a) −𝑓′′(𝜂), (b) −𝜃′(0), and (c) −𝜑′(0) for different values of magnetic parameter 𝑀 and mass transfer parameter 𝑆.

(a)

(b)

Figure 9 

Variation of −𝜃′(0) with (a) Brownian motion parameter Nb and (b) thermophoresis parameter Nt, for different values of magnetic parameter 𝑀.

(a)

(b)

Figure 10 

Variation of −𝜑′(0) with (a) Brownian motion parameter Nb and (b) thermophoresis parameter Nt, for different values of magnetic parameter 𝑀.

(a)

(b)

Figure 11 

Variation of −𝜃′(0) with (a) Brownian motion parameter Nb and (b) thermophoresis parameter Nt, for different values of mass transfer parameter 𝑆.

(a)

(b)

Figure 12 

Variation of −𝜑′(0) with (a) Brownian motion parameter Nb and (b) thermophoresis parameter Nt, for different values of mass transfer parameter 𝑆.

(a)

(b)

Figure 13 

Variation of −𝜃′(0) with (a) Brownian motion parameter Nb and (b) thermophoresis parameter Nt, for different values of Prandtl number Pr.

(a)

(b)

Figure 14 

Variation of −𝜑′(0) with (a) Brownian motion parameter Nb and (b) thermophoresis parameter Nt, for different values of Prandtl number Pr.

(a)

(b)

Figure 15 

Variation of −𝜃′(0) with (a) Brownian motion parameter Nb and (b) thermophoresis parameter Nt, for different values of Lewis number Le.

(a)

(b)

Figure 16 

Variation of −𝜑′(0) with (a) Brownian motion parameter Nb and (b) thermophoresis parameter Nt, for different values of Lewis number Le.

5. Concluding Remarks

The laminar MHD boundary layer flow of nanofluid past an exponentially permeable stretching sheet has been studied. The spatial focus is offered to the effects of Brownian motion and thermophoresis on the heat transfer in boundary layer flow of nanofluid due to an exponentially stretching sheet in presence of magnetic field. The findings of the analysis can be summarized as follows.(a)Due to thermophoresis particle deposition and Brownian motion the heat transfer from the surface to the ambient fluid significantly increases.(b)The temperature and the nanoparticle volume fraction are enhanced with thermophoresis.(c)For Brownian motion the temperature of the fluid increases and the nanoparticle volume fraction decreases.(d)Due to the presence of the nanoparticles in the fluid, the Brownian motion and thermophoresis take place and those cause significant change in the heat transfer efficiency of the fluid.(e)The magnetic field reduces the heat transfer rate, though it causes the increment in the temperature and the nanoparticle volume fraction inside the boundary layer due to exponential stretching of sheet.(f)Due to mass suction, the velocity boundary layer thickness as well as the thermal and nanoparticle volume boundary layer thicknesses becomes thinner, whereas mass injection makes those thicker.

Nomenclature

𝐵∶	Variable magnetic field
𝐵0:	A constant
𝑐:	Stretching constant
𝐶𝑓:	Local skin friction coefficient
𝐷𝐵:	Brownian diffusion coefficient
𝐷𝑇:	Thermophoretic diffusion coefficient
𝑓:	Dimensionless stream function
𝑓′:	Dimensionless velocity
Le:	Lewis number
𝑀:	Magnetic parameter
𝑁:	Nanoparticle volume fraction
𝑁∞:	Constant nanoparticle volume fraction in free stream
𝑁𝑤:	Variable wall nanoparticle volume fraction
𝑁0:	A constant
Nb:	Brownian motion parameter
Nt:	Thermophoresis parameter
Nu𝑥:	Local Nusselt number
𝑝:	A variable
Pr:	Prandtl number
𝑞:	A variable
Re𝑥:	Local Reynolds number
𝑆:	Wall mass transfer parameter
Sh𝑥:	Local Sherwood number
𝑇:	Temperature
𝑇∞:	Constant temperature of the fluid in the free stream
𝑇𝑤:	Variable temperature at the sheet
𝑇0:	A constant
𝑢:	Velocity component in 𝑥-direction
𝑈𝑤:	Stretching velocity
𝑣:	Velocity component in 𝑦-direction
𝑣𝑤:	Variable wall mass transfer velocity
𝑣0:	A constant
𝑥:	Distance along the sheet
𝑦:	Distance perpendicular to the sheet
𝑧:	A variable.

Greek Symbols

𝛼:	Thermal conductivity
𝜂:	Similarity variable
𝜐:	Kinematic viscosity
𝜑:	Dimensionless nanoparticle volume fraction
𝜌𝑓:	Density of the base fluid
(𝜌𝑐)𝑝:	Effective heat capacity of nanoparticles
(𝜌𝑐)𝑓:	Heat capacity of the base fluid
Ψ:	Stream function
𝜃:	Dimensionless temperature.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

One of the authors (Krishnendu Bhattacharyya) is grateful to the National Board for Higher Mathematics (NBHM), Department of Atomic Energy, Government of India, for the financial support in pursuing this work. The authors are thankful for the suggestions and comments of the referees, which have led to improvement of the paper.