Heat and Mass Transfer in the Boundary Layer of Unsteady Viscous Nanofluid along a Vertical Stretching Sheet

Heat andmass transfer in the boundary-layer flow of unsteady viscous nanofluid along a vertical stretching sheet in the presence of magnetic field, thermal radiation, heat generation, and chemical reaction are presented in this paper.The sheet is situated in the xzplane and y is normal to the surface directing towards the positive y-axis. The sheet is continuously stretching in the positive x-axis and the externalmagnetic field is applied to the systemparallel to the positive y-axis.With the help of similarity transformations, the partial differential equations are transformed into a couple of nonlinear ordinary differential equations. The new problem is then solved numerically by a finite-difference scheme known as the Keller-box method. Effects of the necessary parameters in the flow field are explicitly studied and briefly explained graphically and in tabular form. For the selected values of the pertinent parameters appearing in the governing equations, numerical results of velocity, temperature, concentration, skin friction coefficient, Nusselt number, and Sherwood number are obtained.The results are compared to the works of others (from previously published journals) and they are found in excellent agreement.


Introduction
The flow over a stretching surface is an important problem in many engineering processes with applications in industries such as extrusion, melt-spinning, hot rolling, wire drawing, glass-fiber production, manufacture of plastic and rubber sheets, and cooling of a large metallic plate in a bath, which may be an electrolyte.In industry, polymer sheets and filaments are manufactured by continuous extrusion of the polymer from a die to a windup roller, which is located at a finite distance away.The thin polymer sheet constitutes a continuously moving surface with a nonuniform velocity through an ambient fluid [1].Bachok et al. [2] studied boundary layer flow of nanofluids over a moving surface in a flowing fluid and, moreover, a study on boundary layer flow of a nanofluid past a stretching sheet with a convective boundary condition was conducted by Makinde and Aziz [3].Olanrewaju et al. [4] examined boundary layer flow of nanofluids over a moving surface in a flowing fluid in the presence of radiation.An analysis of mixed convection heat transfer from a vertical continuously stretching sheet has been presented by Chen [5].
In many practical situations the material moving in a quiescent fluid is due to the fluid flow induced by the motion of the solid material and/or by the thermal buoyancy.Therefore, the resulting flow and the thermal field are determined by these two mechanisms, that is, surface motion and thermal buoyancy.It is well known that the buoyancy force stemming from the heating or cooling of the continuous stretching sheet alters the flow and the thermal fields and thereby the heat transfer characteristics of the manufacturing processes [6].Effects of thermal buoyancy on the flow and heat transfer over a stretching sheet were reported by many researchers.Chen and Strobel [7] investigated buoyancy effects in boundary layer adjacent to a continuous moving horizontal flat plate.Karwe and Jaluria [8] showed that the thermal buoyancy effects are more prominent when the plate moves vertically, that is, aligned with the gravity, than when it is horizontal.Ali [9] examined the buoyancy effect on the boundary layer induced by continuous surface stretched with rapidly decreasing velocities.Buoyancy driven heat and mass transfer over a stretching sheet in a porous medium with radiation and ohmic heating was studied by Dulal and Hiranmoy [10].Abo-Eldahab and El Aziz [11] presented the problem of steady, laminar, hydromagnetic heat transfer by mixed convection over an inclined stretching surface in the presence of spaceand temperature dependent heat generation or absorption effects.Ali and Al-Yousef [12,13] investigated the problem of laminar mixed convection adjacent to a moving vertical surface with suction or injection.On the other hand, Khan et al. [14] studied the unsteady free convection boundary layer flow of a nanofluid along a stretching sheet with thermal radiation and viscous dissipation effects in the presence of a magnetic field.
The study of MHD boundary layer flow on a continuous stretching sheet has attracted considerable attention during the last few decades due to its numerous applications in industrial manufacturing processes.In particular, the metallurgical processes such as drawing, annealing, and tinning of copper wires involve cooling of continuous strips or filaments by drawing them through a quiescent fluid.Controlling the rate of cooling in these processes can affect the properties of the final product.Thus, rate of cooling can be greatly controlled by the use of electrically conducting fluid and the application of the magnetic field [15].Magnetic field effects on free convection flow of a nanofluid past a vertical semi-infinite flat plate was studied by Hamad et al. [16].Effects of a thin gray fluid on MHD free convective flow near a vertical plate with ramped wall temperature under small magnetic Reynolds number [17] and free convective oscillatory flow and mass transfer past a porous plate in the presence of radiation of an optically thin fluid [18] have been studied.Moreover, MHD Flow and heat transfer over stretching/shrinking sheets with external magnetic field, viscous dissipation, and joule effects were studied by Jafar et al. [19].
Radiative heat transfer in which heat is transmitted from one point to another without heating the intervening medium has been found very important in the design of reliable equipment, nuclear plants, gas turbines, and various propulsion devices for aircraft, missiles, satellites, and space vehicles.Also, the effects of thermal radiation on the forced and free convection flows are important in the content of space technology and processes involving high temperature [10].Influence of thermal radiation, viscous dissipation, and hall current on MHD convection flow over a stretched vertical flat plate was studied by Gnaneswara Reddy [20].Vajravelu and Hadjinicolaou [21] examined heat transfer in a viscous fluid over a stretching sheet with viscous dissipation and internal heat generation.Hady et al. [22] analyzed the flow and heat transfer characteristics of a viscous nanofluid over a nonlinearly stretching sheet in the presence of thermal radiation.Postelnicu [23] studied the influence of chemical reaction on heat and mass transfer by natural convection from vertical surfaces in porous media considering Soret and Dufour effects.Shakhaoath Khan et al. [24] examined possessions of chemical reaction on MHD heat and mass transfer nanofluid flow on a continuously moving surface.Besides this, effects of chemical reactions, heat and mass transfer on nonlinear magnetohydrodynamic boundary layer flow over a wedge with a porous medium in the presence of ohmic heating, and viscous dissipation were studied [25].These days, because of the numerous applications of nanofluids in science and technology, a comprehensive study on heat and mass transfer in the boundary layer of unsteady viscous nanofluid in the presence of different fluid properties is indispensable.The paper entitled "Heat and mass transfer in the boundary layer of unsteady viscous nanofluid along a vertical stretching sheet" in the presence of thermal radiation, viscous dissipation, and chemical reaction is considered.Yet, the entitled paper has not been reported.Accordingly, we extended the works of Vendabai and Sarojamma [28] by incorporating viscous dissipation and chemical reaction terms in the energy and concentration equations, respectively, for more physical implications.The governing equations are reduced to a couple of nonlinear ODEs using similarity transformations; the resulting equations are solved numerically by using the Keller box.Effects of the pertinent parameters involved in the governing equations on velocity, temperature, concentration, skin friction, Nusselt number, and Sherwood number are briefly explained.

Formulation of the Problem
Unsteady two-dimensional boundary layer flow and heat transfer of a nanofluid along a stretching sheet coinciding with the plane  = 0 are considered.The Cartesian coordinate system has its origin located at the leading edge of the sheet with the positive -axis extending along the sheet in the upward direction, while the -axis is measured normal to the surface of the sheet and is positive in the direction from the sheet to the fluid.A schematic representation of the physical model and coordinates system is shown in Figure 1.For  < 0, it is assumed that the fluid and heat flows are steady but the unsteady scenario starts at  = 0. We assume that the sheet is being stretched with the velocity   (, ) along the -axis, keeping the origin fixed.Let   (, ) and   (, ) be the temperature and concentration of the sheet whereas let  ∞ (, ) and  ∞ (, ) be the ambient temperature and concentration, respectively.An external variable magnetic field  =  0 / √ 1 −  is applied along the positive direction.The induced magnetic field is sufficiently weak to ignore magnetic induction effects; that is, magnetic Reynolds number is small.The charge density, external electrical field effects, and polarization voltage are ignored.Using Oberbeck Boussinesq approximation for buoyancy-driven flows, following Buongiorno [6,15,[28][29][30], the governing equations for the flow become The boundary conditions associated to the differential equations are where   (, ) = − 0 / √ 1 −  is suction/injection velocity ( 0 > 0 corresponds to suction velocity),  = (  )  /(  )  ,  = /(  )  ,   is the radiative heat flux, and the heat generation coefficient  is defined by  =  0 /(1 − ), where  0 represents the heat source if  0 > 0 and the heat sink if  0 < 0.
The continuous sheet moves in its own plane with the nonuniform velocity   (, ) = /(1 − ) where  and  are positive constants with dimensions(time) −1 ,  is the initial stretching rate, and /(1 − ) is the effective stretching rate which is increasing with time.  (, ) =  ∞ + (/(1 − ) 2 ) is the temperature distribution of the sheet and the concentration distribution on the wall is given by   (, ) =  ∞ + (/(1 − ) 2 ) where  is a constant and has a dimension (temperature/length) with  > 0 and  < 0 corresponds to the assisting and opposing flows and  = 0 is forced convection limit (absence of buoyancy force).The expressions for   (, ), (), (), and   (, ) are valid only for time  < 1/ unless  = 0.
Using the Rosseland diffusion approximation, the radiative heat flux   is given by We assume that the temperature differences within the flow are sufficiently small such that  4 may be expressed as a linear function of temperature.Giving Taylor series expansion about  ∞ and neglecting higher order terms we get Using ( 6) and ( 7) we obtain the expression   / as follows: In order to transform the governing equations into a system of ordinary differential equations, we introduce the following nondimensional quantities into (1)-( 4): where the stream function  is defined by  = / and V = −/ which identically satisfies the continuity equation (1).Substitution of the similarity variables into (2)-(4) gives and the corresponding boundary conditions become () → 0,  () → 0, ℎ () → 0, as  → ∞.
We are interested to study the skin-friction coefficient   , the local Nusselt number Nu  , and the local Sherwood number Sh  .These quantities are defined as , where   ,   , and   are the skin friction, heat flux, and mass flux at the surface and these quantities are, respectively, defined by Using ( 14) and ( 15), the dimensionless skin friction coefficient (surface drag), wall heat, and mass transfer rates become where Re  =   / is the local Reynolds number.
It is important to mention that if  = 0, the problem under consideration is reduced to the steady state flow scenario; the absence of the viscous and chemical reaction terms in this paper reduces it to the works of Vendabai and Sarojamma [28].It is also worth mentioning that the absence of , Nr, Nt, Nb, , Ec, and  from ( 10) and ( 11) together with the impermeability condition of the sheet corresponds to the works of Ishak et al. [27] whereas the absence of these parameters except Nr and letting  = 0 in their paper corresponds to the works of Vajravelu et al. [6].Assigning  = 0, Λ = 0 and the absence of , , Nt, Nb, , and Ec in (10) and ( 11) reduces these equations to those of Grubka and Bobba [26] with  = 1 in their paper; in this case the problem has the closed form solution as given in the cited journal.

Numerical Solution
As ( 10)-( 12) are nonlinear, it is impossible to get the closed form solutions.Consequently, the equations with the boundary conditions ( 13) are solved numerically by means of a finite-difference scheme known as the Keller-box method, as mentioned by Cebeci and Bradshaw [31].According to Vajravelu et al. [6], the principal steps in the Keller box method to get the numerical solutions are the following: (i) reduce the given ODEs to a system of first order equations; (ii) write the reduced ODEs to finite differences; (iii) linearize the algebraic equations by using Newton's method and write them in vector form; (iv) solve the linear system by the block tridiagonal elimination technique.
One of the factors affecting the accuracy of the method is the appropriateness of the initial guesses.The following initial guesses are chosen: In this study a uniform grid of size Δ = 0.006 is taken and the solutions are obtained with an error of tolerance 10 −5 in all cases, which gives about four decimal places accurate for most of the prescribed quantities as shown in all tables.

Results and Discussion
Our main objective in this flow scenario is to investigate the effects of unsteadiness on the flow quantities of the nanofluid.In order to investigate the flow quantities like velocity, temperature, concentration, and so forth a parametric study has taken place to illustrate effects of the various parameters like magnetic parameter, unsteady parameter, suction parameter, viscous dissipation parameter, buoyancy ratio number, heat source parameter, radiation parameter, Prandtl number, free convection parameter, chemical reaction parameter, Lewis number, thermophoresis parameter, and Brownian motion parameter upon the nature of flow   and transport; the numerical results are displayed graphically (from Figure 2 to Figure 31) and in tables (from Table 1 to Table 5).The present results for temperature gradient are compared to the works of Grubka and Bobba [26], Ishak et al. [27], Vendabai and Sarojamma [28], and Vajravelu et al. [6] and they are found in good agreement as shown in Table 1.
The velocity, temperature, concentration, skin friction coefficient, wall heat, and mass transfer rates for some prescribed values of the various parameters Pr, Nb, Nt, Λ, , Ec, Nr, Le, , , , , , and Re  are briefly presented,  explained, and interpreted graphically and in table form.The velocity, temperature, and concentration profiles are clearly illustrated graphically whereas the skin friction coefficient, wall heat, and mass transfer rates for the pertinent parameters are tabulated.First, let us see effects of these parameters on velocity, temperature, and concentration graphically.
Figure 2 illustrates effects of suction and unsteady parameters on velocity profile.The suction parameter retards the velocity profile both in the steady and in the unsteady cases.It is certain that the role of unsteady parameter is to reduce the velocity profile.It is observed that when the surface is impermeable (  = 0), the velocity in the steady flow case reduces to zero faster (see Figures 3-10).As we move away from the boundary (as the distance  from the boundary increases), the velocity of the nanofluid continuously declines to zero.It is found that as the free convection parameter Λ increases, the velocity profile in both the steady and the unsteady cases increase significantly (Figure 3).It is worth mentioning that the free convection parameter Λ > 0 corresponds to the heating of the fluid (assisting flow), Λ < 0 corresponds to cooling of the fluid (opposing flow), and Λ = 0 means the absence of free convection currents [28].An increase in Λ brings about the enhancement of velocity due to the enhancement of convection currents and thus the momentum boundary layer thickness increases.
Figure 4 depicts effects of buoyancy ratio number Nr and unsteady parameter  on velocity profile.The buoyancy parameter opposes the velocity profile both in the steady and in the unsteady cases whereas the Brownian motion parameter Nb enhances the velocity profile significantly both in the steady and in the unsteady cases.Figure 5 shows the effect of Brownian motion parameter and unsteady parameter on velocity profile.For Nt = 0.5, Nb = 0.5, For Nt = 0.5, Nr = 0.1,  = 0.5, For Nb = 0.5, Nr = 0.1,  = 0.5, Le = 1, Ec = 0.1,  = 0.1, Λ = 1,  = 0.1, Re  = 0.3,   = 0 For Nt = 0.5, Nb = 0.5, Nr = 0.1,  = 0.5, Figures 6 and 7 show the effects of thermophoresis parameter Nt and thermal radiation parameter , respectively, on velocity profile.The increments of both parameters assist the velocity profile to grow both in the steady and in the unsteady cases.
Figures 8 and 13 show effects of magnetic field on velocity and temperature profiles, respectively.The presence of magnetic field reduces the velocity throughout the boundary layer which is in conformity with the fact that the Lorentz force (magnetic force) acts as a retarding force and, consequently, it reduces the momentum boundary layer thickness significantly both in the steady and in the unsteady cases whereas it increases thermal boundary layer thickness.This happens because as the strength of the applied magnetic field increases in an electrically conducting fluid, it produces the resistive Lorentz force ( 2 /).This force decelerates the motion of the fluid in the boundary layer.On the other hand, we can define thermal energy as the additional work done required to drag the nanofluid against the action of the magnetic field .The work done heats up the conducting nanofluid and upgrades the temperature profile.Thus, the presence of magnetic field in the flow regime decreases the momentum boundary layer thickness and enhances the thermal boundary layer thickness [32].
The presence of Prandit number also retards the velocity profile significantly both in the steady and in the unsteady flows as shown in Figure 9.
Figure 10 depicts the effect of heat source parameter  on velocity profile.This parameter enhances velocity of the flow significantly in the steady state situation whereas the enhancement of the parameter in the unsteady state situation is very minimal (see more on explanations for Figure 18).
Figure 11 shows the effect of suction parameter   on temperature profile.The presence of suction in the flow decreases   the temperature profile both in the steady and in the unsteady situations.It is obviously illustrated that the presence of unsteady parameter in a flow is to reduce the temperature profile.In addition to this, the free convection parameter Λ significantly reduces the temperature profile both in the steady and in the unsteady state situations (see Figure 12).On the other hand, in the steady state flow scenario, the thermal boundary layer thickness grows considerably with increasing value of the Lorentz force.But the increment of the thermal boundary layer in the case of unsteady flow situations is nominal (see Figure 13).Figures 14 and 15 show the effects of thermal radiation and Brownian motion parameters, respectively, on temperature.Both parameters significantly enhance the thickening of the thermal boundary layer thickness in the steady and unsteady state situations.As we move away from the boundary farther, the Brownian motion parameter makes the velocity decay to zero faster.This is due to the fact that thermal radiation inspires in thickening the thermal boundary layer at the expense of releasing heat energy from the flow region and it causes the system to cool.In reality this is true because   temperature increases as a result of increasing the Rosseland diffusion approximation for radiation   .
Figure 16 shows the effects of thermophoresis and unsteady parameters on temperature profile.In the steady state situation, thermophoresis parameter enhances the thickening of thermal boundary layer thickness considerably whereas the thickening of the boundary layer in the unsteady state situation is not significant.Figure 17 illustrates the effect of Eckert number Ec on temperature profile.As the viscous dissipation parameter increases, the temperature profiles of both the steady and the unsteady state flows also increase but not significantly.
Figure 18 shows effects of the heat source parameter on temperature.Increasing the heat source parameter from 0.1 to 0.4 increases the temperature profile and hence it thickens the thermal boundary layer moderately in both the steady and the unsteady state flows.An increase in the values of heat source parameter  increases both the velocity and the temperature profiles.Because the presence of source of heat (in the flow regime) enhances thermal energy, as a result of this, temperature profile rises.The rise in temperature allows the fluid to increase the velocity profile due to the effect of buoyancy; effects of heat sink on velocity and temperature profiles play oppositely.
As the buoyancy ratio number increases from 0 to 0.6, the thermal boundary layer thickness increases in the case of steady state flows.For unsteady state situations, the effect of buoyancy ratio number on temperature is not significant as shown in Figure 19.
The effect of Prandtl number Pr on temperature is shown in Figure 20.This parameter significantly reduces the thermal boundary layer thickness in the cases of both steady and  unsteady state situations.By definition, Prandtl number is a dimensionless number which is the ratio of momentum diffusivity to thermal diffusivity; that is, Pr =   /.An increase in the values of Pr is equivalent to momentum diffusivity which dominates thermal diffusivity.Hence, thermal boundary layer thickness reduces as Pr increases.This is due to the fact that the larger the Prandtl number Pr is, the higher the viscosity (sticker) of the fluid and the thicker the momentum boundary layer will be compared to the thermal boundary layer.Consequently, heat transfer will be less convective.On the other hand, the presence of the free convection parameter Λ minimizes nanoparticles volume fraction moderately in the cases of both steady and unsteady state flows (see Figure 21).
The effect of Brownian motion parameter on nanoparticles concentration is shown in Figure 22.As this parameter increases, nanoparticles volume fraction decreases both in the case of steady and in the case of unsteady state situations.Figure 23 shows the effect of thermophoresis parameter on nanoparticles volume fraction.It is observed from the figure that the concentration decays monotonically to zero as the   distance  increases from the boundary.For larger values of Nt in the steady state situation, the concentration profile attains its maximum in the boundary layer and then declines to zero faster when we go farther.Increment of this parameter enhances the nanoparticles volume fraction strongly in both the steady and the unsteady state cases.
Effects of Prandtl number on concentration profile is depicted in Figure 24.Increasing the value of Prandtl number leads to an increase in the concentration profile of both the steady and the unsteady state cases.This increment ceases as we go farther from the boundary.The effect of Lewis number  on concentration profile is shown in Figure 25.It is observed that increasing Lewis number significantly decreases the nanoparticles volume fraction in the case of both steady and unsteady state situations.Lewis number is a dimensionless number which is defined as the ratio of thermal diffusivity to mass diffusivity.Increasing the value of Le means increasing thermal boundary layer thickness at the expense of reducing the concentration boundary layer thickness; this leads to an increase in mass transfer rate.Figure 26 illustrates the influence of chemical reaction parameter on concentration profile.In the steady state situation, it is observed that the chemical reaction parameter enhances nanoparticles volume fraction but in the case of unsteady state situations, the effect of this parameter is minimal.It is understandable that the effect of unsteady parameter is to reduce the nanoparticle volume fraction.Moreover, the concentration decreases slowly as the suction parameter increases in the case of both the steady and unsteady state flow conditions (see Figure 27).profile.In the steady state flow, both parameters moderately enhance the nanoparticles volume fraction.On the other hand, in the case of unsteady flow, the effect of magnetic field on the concentration profile is minimal whereas the effect of buoyancy parameter is almost nil.
Figure 30 shows the effect of radiation parameter on concentration profile.In the case of both steady and unsteady state situations, increment of the radiation parameter minimizes the nanoparticles volume fraction and soon after a certain distance from the boundary the concentration  profile remains constant.Moreover, in the case of steady state situation, the heat source parameter significantly decreases the concentration profile near the boundary layer while the effect of this parameter far from the boundary in both steady and unsteady cases has negligible effect on concentration (see Figure 31).
Table 2 shows the effect of unsteady, free convection, buoyancy ratio, and magnetic parameters on skin friction coefficient-  (0), wall heat transfer rate-  (0), and wall mass transfer rate-ℎ  (0).Increasing the value of unsteady  parameter increases all the skin friction coefficient, wall heat transfer rate, and wall mass transfer rate.On the other hand, free convection parameter enhances both wall mass and wall heat transfer rates while it reduces the skin friction coefficient.Both buoyancy ratio and magnetic parameters retard heat transfer and mass transfer rates whereas these parameters enhance the skin friction coefficient.
Effects of suction, thermal radiation, and Brownian motion parameters on skin friction, heat transfer, and mass transfer rates are shown in Table 3.It is observed that the enhanced values of suction parameter lead to increase skin friction, heat transfer, and mass transfer rates.On the other hand, both radiation, and Brownian motion parameters enhance mass transfer rate while these parameters reduce both skin friction and heat transfer rate.
Table 4 shows influences of thermophoresis parameter, Prandtl number, and heat source parameters on skin friction, Nusselt number, and Sherwood number.Thermophoresis parameter reduces skin friction, heat transfer, and mass transfer rates.However, increasing the values of Prandtl number leads to increase both skin friction coefficient and heat transfer rate while it decreases mass transfer rate.This is due to the fact that the higher the Prandtl number the thinner the thermal boundary layer and the thicker the nanoparticle volume fraction boundary layer (refer to Figures 20 and  24).As a result of this, the rate of heat diffusion increases while the rate of mass diffusion decreases with increasing values of Pr.On the other hand, increasing the values of heat source parameter results in reducing both the skin friction coefficient and heat transfer rate whereas it makes mass transfer rate increase more.
Table 5 shows the effect of Lewis number, viscous dissipation parameter, and chemical reaction parameter on skin friction, wall heat, and wall mass transfer rates.As it is clearly shown, Eckert number and Lewis number reduce both skin friction and heat transfer rate while both parameters increases mass transfer rate.However, chemical reaction parameter retards both the skin friction coefficient and the heat transfer rate whereas it enhances mass transfer rate.

Conclusion
The problem of heat and mass transfer in the boundary layer of unsteady viscous nanofluid along a vertical stretching sheet has been studied.The nonlinear governing equations associated to the boundary conditions were transformed into coupled ODEs with the help of similarity transformation equations.The solutions of these problems are solved numerically with the help of the Keller box method.The following are some of important results among many: (i) The velocity, temperature, and concentration profiles of the unsteady flow are less than the corresponding parts of the steady state flow scenario.
(ii) The velocity profile decreases with an increase in the buoyancy ratio number, magnetic parameter, suction parameter, and Prandtl number whereas it increases with the increment of free convection, radiation, Brownian motion, thermophoresis, and heat source parameters.
(iii) Viscous dissipation, thermal radiation, Brownian motion, buoyancy ratio, heat source, and magnetic and thermophoresis parameters enhance the temperature profile whereas Prandtl number, suction parameter, and free convection parameter reduce the temperature profile.
(iv) Nanoparticles volume fraction is enhanced by Prandtl number, thermophoresis parameter, magnetic parameter, and buoyancy ratio number whereas Lewis number, heat source parameter, Brownian motion parameter, radiation parameter, chemical reaction parameter, suction parameter, and free convection parameter reduce the concentration profile.
(v) The skin friction coefficient, heat transfer, and mass transfer rates of unsteady flow are greater than the corresponding parts of the quantities in the steady state situations.
(vi) Increasing the values of the buoyancy ratio number, magnetic parameter, suction parameter, and Prandtl number enhances the skin friction coefficient whereas it decreases with increasing values of free convection, Brownian motion, radiation, thermophoresis, chemical reaction, heat source, viscous dissipation parameters, and Lewis number.
(vii) The presence of Lewis number, thermal radiation, Brownian motion, heat source, magnetic, buoyancy ratio, chemical reaction, thermophoresis, and viscous dissipation parameters in the flow field is to reduce the rate of thermal boundary layer thickness whereas Prandtl number, suction parameter, and free convection parameter maximize the thermal boundary layer thickness.
(viii) The wall mass transfer rate is an increasing function of Lewis number, free convection, chemical reaction, viscous dissipation, suction, Brownian motion, heat source, and radiation parameters while Prandtl number, thermophoresis, and magnetic and buoyancy ratio parameters reduce the mass transfer rate at the plate surface.

Figure 1 :
Figure 1: Physical model and coordinate system.

Figure 2 :
Figure 2: Effects of   and  on velocity profile.

Figure 3 :
Figure 3: Effects of Λ and  on velocity profile.

Figure 4 :
Figure 4: Effects of Nr and  on velocity profile.

Figure 5 :Pr = 1 ,
Figure 5: Effects of Nb and  on velocity profile.

Figure 6 :
Figure 6: Effects of Nt and  on velocity profile.

Figure 7 :Pr = 1 ,
Figure 7: Effects of  and  on velocity profile.

Figure 8 :
Figure 8: Effects of  and  on velocity profile.

Figure 9 : 5 , 1 R
Figure 9: Effects of Pr and  on velocity profile.

Figure 10 :
Figure 10: Effects of  and  on velocity profile.

Figure 16 :
Figure 16: Effects of Nt and  on temperature profile.

Figure 18 :
Figure 18: Effects of  and  on temperature profile.

Figure 19 :
Figure 19: Effects of Nr and  on temperature profile.

Figure 20 :
Figure 20: Effects of Pr and  on temperature profile.

Figure 22 :
Figure 22: Effects of Nb and  on concentration profile.

Figure 23 : 1 𝜙
Figure 23: Effects of Nt and  on concentration profile.

Figure 24 :
Figure 24: Effects of Pr and  on concentration profile.
Figure 26  illustrates the influence of chemical reaction parameter on concentration profile.In the steady state situation, it is observed that the chemical reaction parameter enhances nanoparticles volume fraction but in the case of unsteady state situations, the effect of this parameter is minimal.It is understandable that the effect of unsteady parameter is to reduce the nanoparticle volume fraction.Moreover, the concentration decreases slowly as the suction parameter increases in the case of both the steady and unsteady state flow conditions (see Figure27).Figures 28 and 29 illustrate effects of magnetic parameter and buoyancy ratio number, respectively, on concentration

Figure 25 : 1 0
Figure 25: Effects of Le and  on concentration profile.

Figure 26 :
Figure 26: Effects of  and  on concentration profile.

Figure 28 :
Figure 28: Effects of  and  on concentration profile.

Figure 30 :
Figure 30: Effects of  and  on concentration profile.

Nomenclature𝐴 :
Unsteady parameter : Magnetic field strength : Nanoparticle concentration   : Skin-friction coefficient   : Nanoparticles concentration at the stretching surface  ∞ : Nanoparticle concentration far from the sheet   : Specific heat capacity at constant pressure   : Brownian diffusion coefficient   : Thermophoresis diffusion coefficient : Dimensional heat generation coefficient Ec: Eckert number : Dimensionless stream function  * : Acceleration due to gravity : Dimensionless temperature ℎ: Dimensionless concentration  * : Rosseland mean absorption coefficient   : Thermal conductivity of the base fluid   : Thermal conductivity of the nanoparticles   : Chemical reaction parameter Le: Lewis number : Magnetic parameter Nb: Brownian motion parameter Nr: Buoyancy ratio number Nt: Thermophoresis parameter Nu  : Nusselt number Pr: Prandtl number : Radiation parameter Re  : Local Reynolds number Sh  : Sherwood number : Fluidtemperature   : Temperature at the surface  ∞ : Temperature of the fluid far away from the stretching sheet , V: Velocity components in the and -axes, respectively   : Velocity of the wall along the -axis ; : Cartesian coordinates measured along the stretching sheet.Greek Symbols : Thermal diffusivity of the base fluid : Thermal expansion coefficient : Electrical conductivity  * : Stefan-Boltzmann constant : Stream function : Dimensionless similarity variable Λ: Free convection parameter : Scaled chemical reaction parameter : Heat source parameter : Dynamic viscosity of the base fluid : Kinetic viscosity of the base fluid   : Density of the base fluid   : Density of the nanoparticle

Table 1 :
Comparison of some of the values of wall temperature gradient −  (0) obtained by the cited researchers and the present results for  =  = Nb = Nt = Nr = Ec =  =   = 0.

Table 2 :
Effects of free convection parameter Λ, unsteady parameter , buoyancy ratio number Nr, and magnetic parameter  on skin friction coefficient, heat transfer, and mass transfer rates.

Table 3 :
Effects of suction parameter   , unsteady parameter , radiation parameter , and Brownian motion parameter Nb on skin friction coefficient, heat transfer, and mass transfer rates.

Table 4 :
Effects of thermophoresis parameter Nt, unsteady parameter , Prandtl number Pr and heat source parameter  on skin friction coefficient, heat transfer and mass transfer rates.

Table 5 :
Effects of viscous dissipation parameter Ec, unsteady parameter , Lewis number Le, and chemical reaction parameter  on skin friction coefficient, heat transfer, and mass transfer rates.