Nanofluids Thin Film Flow of Reiner-Philippoff Fluid over an Unstable Stretching Surface with Brownian Motion and Thermophoresis Effects

The current investigation is carried out on the thin film flow of Reiner-Philippoff fluid of boundary-layer type. We have analyzed the flow of thin films of Reiner-Philippoff fluid in the changeable heat transmission and radiation over a time-dependent stretching sheet in 2D. The time-dependent governing equations of Reiner-Philippoff fluid model are simplified with the help of transformation of similarity variables. To investigate the behavior of the Reiner-Philippoff fluid with variable stretching surface for different physical effects, we considered thermophoresis and Brownian motion parameters in the flow. The Homotopy Analysis Method is implemented in the reduced model to achieve a solution of the original problem. A numerical convergence of the implemented method is also analyzed. The behavior of temperature, velocity, and concentration profiles have been investigated with the variation of skin friction, Nusselt number, and Sherwood number. A comparative graphical survey is presented for the velocity gradient, under different parameters. An analytical analysis is presented for the time-dependent parameter over thin film flow. The results we obtained are better than the previously available results. For the survey, the physical representation of the embedded parameters, like, β depends on the stretching parameter ζ, and the Reiner-Philippoff fluid parameter e are discussed in detail and plotted graphically. Prandtl number Pr, Brownian motion parameter Nb, thermophoretic number Nt, and Schmidt number Sc are presented by graphs and discussed in detail.


Introduction:
In the last few decades, scientists have given great attention to thin film flow.The basic idea behind such an important concept is the applications and mechanism of thin film flow.Thin films are ubiquitous in nature and their mechanism is important to understand, because it has a wide range of practical uses.The traditional theory of Newtonian fluids is mainly focused on the linear relationship between stress and strain tensor or on their rates.Newtonian fluids are hardly pointing towards the doctrine of linear relation between the stress and strain tensor.In the same way, fluids that do not agree with the linear distribution between the stress and strain tensor are known as non-Newtonian fluids.There are two main classes of non-Newtonian fluids, visco-inelastic and visco-elastic fluids.The subsequent has great importance due its dual nature.Visco-elastic fluids, due to their viscosity, show an elastic behavior up to some degree.In such fluids, energy is stored in the form of strain energy.One cannot ignore the strain in these types of fluids, because it is responsible to recover the original state.These fluids came under the definition of Newtonian fluids, when there is a first order derivative of all the tensor formed from the velocity field, while fluids having higher order derivative tensors present come under the definition of non-Newtonian fluids.Some of these fluids are second order, third order, and fourth order fluids.Third order liquid, which is due to its sequential nonlinear constraints, has a variety of applications.
The study of non-Newtonian fluids from a theoretical point of view is too complicated and we need a mathematical relation to briefly illuminate the relation between shear stress and the shear rate.Therefore, a variety of non-linear relations have been suggested for the study of these fluids.As a consequence, to deliberate such a categorizing demeanor, various models have been initiated and discussed in [1].Non-Newtonian fluids are not so simple that a single mathematical relation can explain the whole scenario.Therefore, for this purpose, several models have been initiated and developed to briefly explain the nature of such fluids.Among all, Ellis fluid model, Sisko fluid model, the Carreau viscosity model, Ostwald-de Waele model, the cross-viscosity model, Carreau-Yasuda model, Powell-Eyring model, and Reiner-Philippoff fluid model are the most important models in explaining the nature of such fluids.These models have interesting properties in their own, in which the Ostwald-de Waele model is considered to be the basic model, normally known as the Power law model.The subsequent of all models discussed here is a time-independent three-parameter model, behaving like non-Newtonian fluids under intermediate shear rates and Newtonian fluids at extreme shear rates.Consequently, due to their dual nature, Reiner-Philippoff fluid has many applications in engineering sciences and other technologies.
A variety of fluids are available in nature, in which nanofluids are the most interesting fluids, due to their variety of applications.The most commonly uses of nanofluids are metal oxides, oxide ceramics and chemically stable metals, like Alumina, Silica, Zirconia, Titania, aluminum oxide, copper oxide, gold, copper and various forms of allotropes of carbon and metal carbides.Water, oils, polymeric solutions, lubricants, bio-fluids, and glycols are normally used as base fluids.Nanofluids are two-phased mixtures designed by spreading nanometer-sized particles, in which base fluid size ranging up to 100 nm.Nanoparticles play a key part in heat transfer analysis.On the other hand, applications of liquid film flow grow day by day.The most common uses of these flows are in heat exchange processes, techniques of coating, industrial and distillation processes, and many more.The applied usage of the liquid film flow is a fascinating interaction amongst fluid mechanics, structural mechanics, and technology.Some of the practical usages are polymer and metal extrusion, foodstuff processing, plastic sheets depiction, casting and fluidization of the reactor.
In the assessment of these applications, researchers have taken a keen interest in the study of liquid film flow on unstable surfaces.Stretching sheets at the beginning were treated as linear surfaces.Such a phenomenon encountered in many industrial processes, like, in cooling, extraction of polymer sheets and plastic sheets, etc.In these industrial processes the stretching sheet contacts with the fluid both mechanically and thermally.Sakiadis [2] work is considered to be the pioneering work in the study of boundary-layer flow over non-stationary and rigid surfaces.But in polymer industries stretching sheet play a key role, which is explained by Crane [3] in his famous work.After the work of Crane on stretching sheet and its numerous applications in the polymer industry, researchers have shown great interest in it.Stretching sheet problems have been investigated by different researchers under different physical parameters with their variations, like viscosity and thermal conductivity, magnetic and electric fields (MHD), thermal radiation, viscous dissipation, and chemical reactions etc.
Similar to the variations in the parameters for the stretching sheet, the same phenomenon of stretching problem is studied for different geometries and is developed from time to time.Siddiqui et al. [4,5] investigated non-Newtonian fluids on a moving built with a sloping plane in one direction for thin film flow.Tawade et al. [6] studied the effect of magnetic field upon a thin fluid stream passes over a temperamental stretching sheet with heat.They used two different numerical approaches, Newton-Raphson and RK-Felberg.They briefly explained and provided a detailed survey of different physical parameters.
Beside all these, the implemented techniques are not to be ignored.In literature, the boundarylayer equations obtained for stretching flow are solved by different approaches.Among all, numerical and perturbation approaches have been adopted and applied by many researchers successfully.In practice, numerical techniques are too difficult to apply due to the high non-linearity of the model equations, whereas perturbation techniques are not always applicable.Perturbation techniques need a small or large parameter to be presented in the equation, which is not always available to us.To overcome this situation, some new techniques were developed and implemented by the researchers successfully.Sajid and Hayat [7] used HAM and HPM to thin film flow of Sisko fluid and Oldroyed-6 constant.The effect of thermal radiations of blending convection stream over a steeping surface in a permeable channel is studied by Bakier [8].Nargis and Tahir [9] have given a more detailed survey of grade third fluid on a moving belt in the direction of a slanted plane.Stretching problem in permeable medium with thermal effects of a slanted plate is investigated by Moradi et al. [10].Chaudhary et al. [11] re-examined thermal radiation impacts of liquid on exponentially extending surfaces.Eldabe et al. [12] examined convection, radiation, and synthetic effects of MHD visco-elastic fluid flow in a permeable channel on a horizontal stretching sheet.Das [13] has investigated some important properties of thermal radiation and thermophoresis of MHD blended convective flow.Recently, Hsiao [14] has examined the heat and mass trasfer effects of Maxwell fluid.MHD flow of different models, like Powell-Eyring nanofluid and other non-Newtonian fluids on stretching surfaces are briefly explained in [15][16][17].Crane [3] for the first time studied the flow of gummy liquid in a stretched surface.The effect of heat exchange on an extending sheet for viscoelastic liquids is discussed by Dandapat [18].Wang [19] for the first time studied finite liquid film on an unsteady stretching domain.The problem discussed by Wang [19] has discussed by Usha and Sridharan [20] with a survey on different parameters.For a heat transfer analysis of liquid film fluid, numerical results for different parameters were obtained by Liu and Andersson [21].The repercussion in the thin liquid film on an unsteady stretching sheet due to the inner heat production was examined by Aziz et al. [22].
The thin liquid film flow of non-Newtonian fluids has a lot of practical features.Consequently, it becomes a common solute in engineering and other technologies.Andersson et al. [23,24] investigated the non-Newtonian thin liquid films at a time depending stretching sheet by taking the Power law model in consideration.After this pioneer investigation of Andersson, scientists have given more attention to stretching problems by using the Power law model, for more detail see [25,26].Other models also came in discussion during this era.Megahed et al. [27] examined the thin liquid film flow of Casson liquid for viscous promulgation with slip velocity and the transmission of variable heat transition.The same scenario was discussed by Abolbashari et al. [28] for nano particles with the generation of entropy.Buongiorno's model for nano fluid thin film on a temperamental extending stretched sheet was recently investigated by Qasim et al. [29].A steady flow of liquids through a porous medium is studied by Ariel [30].Ariel got a high non-linear coupled boundary value problem for the geometry under consideration and applied numerical methods to obtain an appropriate solution.Sahoo et al. [31] investigated heat exchange analysis with a uniform oblique magnetic field for non-Newtonian fluids.They successfully applied finite difference and Broyden's methods for the concatenation of the field of the velocity.Aiyesimi et al. [32,33] examined the thin liquid film flow of an MHD grade third fluid and obtained some interesting results by using perturbation techniques with a brief survey impact of slip parameters and magnetic parameters.Third grade fluid and its approximate analytical solution by using OHAM for three different kinds of flow has discussed by Islam and Shah et al. [34,35].Makinde [36] studied the same geometry with isothermal effects for hydro-dynamically third order liquid film flow.The approximate solution for velocity and temperature was obtained by a Hermitepade method.A brief discussion and explanation was given by Yao and Liu [37] of the second order fluid over flat plates for unsteady flows.Erdogan et al. [38] examined the properties of unsteady flow of the non-Newtonian fluids, with a brief description on Poiseuille, Couette, and generalized form of Couette flow.Abdulhameed et al. [39] successfully applied Laplace transformation, perturbation techniques, and separation of variable methods for the clarity of unsteady non-Newtonian fluids over an oscillating plate.Huan was the first to settle variational rules for nano thin film-lube [40] with the help of the method of semi-inverse [41][42][43][44][45]. Kapitza [46], Yih [47], Krishna and Lin [48], Anderson and Dahl [49], and Cheng et al. [50] considered thin film flow problems with distinct geometric expressions.As time passes, the thin film flow applications in engineering sciences increase day by day and as a result, the researchers extended the work to a new world.Recently, coating and fiber applications of thin film flow are described and discussed in [51][52][53][54].The geometry and other physical constraints have fixed, but some impurities have been introduced to the study as discussed in [55,56], to improve and enhance the heat transfer analysis.
An interesting and remarkable behavior of time-dependent non-Newtonian fluid is its pseoudo-plasticity, which vanishes with the expansion of shear rate.Many models have initiated as discussed earlier in the investigation of the behavior of such fluids.Among all, in Reiner-Philippoff fluid model [57], researchers have shown great interest.In 1965, Kapr and Gupta [58] have studied Reiner-Philippoff fluid two-dimensional flow in a linear channel.Different approaches have been adopted with different geometrical aspects to discuss this famous model.
In 1994, Tsung-Yen [59] investigated the boundary-layer flow problems by using this model.The basic idea of boundary-layer theory was introduced by L. Prandlt in Heidelberg, Germany, in August 1904, at the third International Congress of Mathematicians.It is the region in the fluid flow, developing at large Reynolds numbers.This region is strongly affected by inertial forces and viscous forces.Boundary layer theory is very important and has a variety of dimensions and visual perception of interest, and has been studied for a long interval of time.In 2009, this model was investigated by Yam [60] for the boundary flow past a stretching wedge.During this era Patel and Timol [61] used the technique of similarity solution for three dimensional boundary layer type equations for non-Newtonian fluids.Ahmad [62] examined the Reiner-Phillippoff fluid flow based nano-liquids past a stretching sheet.Recently, Ahmad et al. [63] discussed the same model with the same geometry with shifting and thickness in the stretching sheet.
Discussion has shown that different problems arise due to varying geometry as well as the fluid behavior, and different approaches are adopted by the researchers to meet their needs.Most of the problems that arise are highly non-linear and it is a difficult job to handle such problems with the usual available techniques in literature.Nowadays, perturbation techniques [64,65] are in the main stream for dealing with such problems.These methods work in the presence of small or large scale parameters.These parameters are not always available to us in applied sciences and we cannot apply these techniques to these types of problems.To deal with such problems, we use non-perturbative techniques like the "Lyapunov's artificial small parameter method" [66], the method of d-Expansion, and the adomians decomposition method (ADM) [67].Various approaches have been adopted by researchers to find the solution to their problems.An exact solution in literature is very rare.This is because of the complexity of the geometry of the problem.That is why we often see numerical approaches to find the approximate solution.Among all, homotopy analysis method (HAM) [68][69][70] is the one by virtue of which we can find the approximate solution.In 1992, Liao [71] for the first time developed and implemented this method and found solution in the form of series in a single variable.Liao also discussed the convergence of this proposed method and found a rapid convergence.HAM has some interesting points of interest.Most importantly, this method is independent of whether a given non-linear problem contains any small or large parameters or not.In HAM, we can modify and control the region of convergence, where necessary, and is helpful in selecting distinct sets of base operations, which approximate a non-linear problem with less effort.
The goal of our current investigations is to obtain the thin liquid film flow of Reiner-Philippoff fluid over a stretching sheet with heat transfer and thermal radiations.Boundary-layer equations are obtained from the physical demonstrated geometry.Thermophoresis effects and Brownian motion are also encircled with different physical parameters.A similarity solution is obtained with the help of new variables, due to which a complicated model is transformed into simple coupled ordinary differential equations.An analytical approach is adopted for the solution of the reduced system.HAM is implemented with initial guess as required for the implementation of the technique, due to its fast convergence.With the variation of different physical parameters, the results are plotted, tabulated, and discussed in detail.The physical significance of Sherwood number and skin friction is presented by tables.

Problem Formulation
Consider an electrically conducting and time-dependent thin film flow of Reiner-Philippoff fluid over spreading sheet.The elastic sheet start moving from fixed slit.The coordinates oxyz are adjusted in such way that ox and the plate are equal, and oy is along the sheet.The origin is at rest in the direction of the x-axis, due to the two equal and opposite forces of the stretching sheet flow.We take the x-axis in the direction spreading sheet and has the stress velocity U w (x, t) = γx(1 − ζt) −1 , where γ and ζ represent any fix numbers, which are vertical to y-axis as shown in Figure 1.The wall temperature [72,73] of the liquid is and the capacity of the nanoparticles is given by where ν f denotes the fluid kinematic viscosity, T 0 and C 0 denote the temperature of the slit and volume friction of the nanoparticles, while T r and C r represent the reference temperature and reference volume of the nanoparticles respectively.Assume that the effects of body forces are negligible in the field of flow.In light of the previous assumptions, the equation of continuity, the basic boundary governing equations, concentration, and heat transferring equations take the following forms.
The stress deformation behavior is well explained by Reiner-Philippoff in [74], and is considered one of the classical descriptions and is given by an implicit functional way: where the parameters µ 0 , µ ∞ , and τ 0 are greater than zero.This model gives interesting results due to its behavior, for large or small values of the model nearly agree with Newtonian fluids.Besides this, when the values of τ 0 are in between the extremes the model agrees with non-Newtonian fluids.
The Momentum equations for two dimensional flow together with the continuity equation takes the form: The components of stress, presented above, are difficult to present in a closed single explicit format.For this purpose, we assume a small τ 0 such that its higher powers greater than three vanishes.The stress components take the form: with the constitutive relations defined by (8) where τ 0 xx denotes the Newtonian stress component of the fluid with µ 0 , the coefficient of viscosity and τ xx is the residual contribution terms of the Reiner-Philippoff fluid, and prime should not be considered a derivative.
Using Equation (7) in Equation (3), we get Neglecting fourth and higher order terms, we get Hence, τ ij is explicitly related to the gradient of the velocity by the above relation.
Let us consider τ 0 as large, such that its exponent larger than three is negligible in the relation expressed below.
with the constitutive relations: (12) where τ ∞ xx denotes the Newtonian fluid stress component with the viscosity coefficient µ ∞ and τ xx , the residual contribution terms of the Reiner-Philippoff fluid, while the prime should not be considered a derivative sign.
In the absence of pressure gradient, the boundary-layer equations are simplified to ũ ∂y (16) where ∂T ∂t With the constraints defined at the boundaries ũ = U w , ṽ = 0, Here, ũ and ṽ represent the state variables, denotes the velocity components along x-axis and y-axis respectively, and ρ represents the density of the fluid, local temperature by T, and the fluid capacitance by C p .The ratio is the characteristic ratio of the base fluid to the nanoparticles heat capacitance; D B represents the direct Brownian diffusion constant; D T represents thermophoretic diffusion constant, K is the thermal conductivity, and T ∞ denotes the fluid temperature far away from the slit.Introducing the succeeding similarity transformations [75,76] where prime represents the change with respect to η, β = , for detail see [77,78].With the help of the newly introduced similarity transformations, Equations ( 14)-( 21) are reduced to the following equations, while the continuity equation is satisfied identically.
1 Pr The boundary constraints of the problem are: The generalized physical constraints obtained are defined as: S = γ is the non-dimensional measure of unsteadiness, = x υ , λ = µ 0 µ ∞ and γ = τ 2 or where Re x is known as the local Reynolds number and is defined as Re x = Ũw x ν and τ w is the value of τ on η = 0. Nu is the Nusselt number and is defined as , while Q w denotes the heat flux and represents the Sherwood number in which J w is the mass flux, where .
Sherwood number Sh and Nusselt number Nu take the dimensionless forms:

Solution by HAM
The approximate solution of the Equations ( 23)-( 26) corresponding to the Equations ( 27) and ( 28) are treated with Homotopy Analysis Method (HAM).The auxiliary parameters encircled the solution which normalize and switch to conjunction of the solutions.Let us take the initial guesses: Let us denote the linear operators by L f , L θ , and L φ defined as: with the property where j for j ∈ {1, 2, ..., 7} are the general solution coefficients.The fundamental procedure of the solution by using HAM is explained in [72,73,76].

Results and Discussion
The current analysis is carried out on the thin film flow of Reiner-Philippoff fluid of boundary layer type over a time-dependent stretching plate.The aim of this subsection is to study the velocity distribution, temperature distribution, and concentration profile with physical effects of different embedding parameters, which are discussed in Figures 2-16.
Figure 2 demonstrates the thin film thickness β during the fluid motion.In the performance of coating, thin film thickness play a key role.Physically, the thickness of the films is directly related to the velocity.Figure 2b  Figure 3 illustrates the effect of the stretching parameter .The velocity profile shows an increase with the increasing values of , because the lower plate always behaves directly to the flow fluid velocity.Physically, for > 0, the surface accelerating rises, < 0 decelerating the surface, while = 0 depict the random motion of the surface.Figure 3b  Figure 4 reveals the behavior of unsteady constraint S over f (η) for dissimilar values of the embedded parameters.It is observed that the velocity profile f (η) directly varies with unsteadiness parameter S. The velocity profile climes up with the increasing behavior of S. Furthermore, it is observed that the solution is possible only in the closed interval [0, 1] for S. Also, the increasing values of S increases the motion of the nanofluid.The effect of the thin film thickness β on temperature profile for different values of the embedded parameters shown in Figure 6.A similar effect in velocity profile is observed for β.The larger the thickness of the liquid film, the lesser the heat transfer.In other words, the flow of heat in the larger thickness film faces more more difficulty, as compared to a lesser thickness film.
Figure 7 illustrates the temperature distribution under Brownian motion parameter Nb.In general, due to the irregular motion of the particles, this causes a collision between these particles.An increase in heat of the fluid can be seen with the ascending order of the Brownian motion parameter Nb, consequently, free surface nanoparticle volume friction decreases.
The impact of the unsteadiness parameter S on the heat profile θ(η) is presented in Figure 8.It is observed that θ(η) varies directly with S.An increase in S increases the temperature of the fluid, which further increases the kinetic energy of the fluid, and results in increment of the liquid film motion.Figure 9 illustrates the effect of thermophoresis parameter Nt on temperature profile.The limitations thermophoresis helps in the increase of a surface temperature.The irregularity in motion (Brownian motion), causes a temperature increase due to the kinetic energy produced by nano suspended particles, which results in thermophoretic force generation.The intensity produced by this force compels the fluid to move away from the stretching sheet.As a result, larger values of Nt cause an increase in temperature, due to which the surface temperature also increases.
Figure 10 describes the Schmidt number effect over temperature profile.Schmidt number physically relates the boundary-layer of mass transfer to the hydrodynamics layer.Increasing rate of the viscous diffusion keeping the mass flux constant increases the Schmidt number, which as a result decreases the heat profile, as shown in the figure.Figure 12 describes the effect of thermophoresis parameter Nt on concentration field.It is clear from the figure that an increase in Nt increases the concentration field.This is because higher values of Nt increase the nanofluid molecules kinetic energy, and as a result the concentration increases.Figure 14 shows the concentration profile φ(η) behavior, under the effect of the unsteadiness parameter S. A direct relation has been observed between the unsteadiness parameter S and φ(η) the concentration profile .Increasing the unsteadiness parameter S, causes an increase in the temperature to be observed, that blows the kinetic energy off the fluid, which leads to an increase in the concentration of the liquid film.
Figure 15 reveals the opposite information as discussed in the temperature distribution under different parameters.The above diagram shows that the concentration profile decreases due to an increase in Schmidt number Sc, which as a result reduces the boundary-layer thickness.This is because of the physical significance of the Schmidt number, which relates both the mass and hydrodynamic layer.
Figure 16 illustrates the effect of thin film thickness β on φ(η) for the different values of the embedded parameters.It is clear that the concentration profile falls with higher values of β.The same effect has been observed for β in the velocity distribution as well as in temperature distribution.

Tables Discussion
Table 1 depicts the influence of Nusselt number.The effects of Nt, S, β, and Pr on Θ (0) are shown.It is clear that larger values of Nt and S decrease Θ (0), while the unsteadiness parameter S and thickness parameter β increase Θ (0).The influence of stretching parameters and λ of Reiner-Philippoff fluid and unsteadiness parameter S on skin friction C f is presented in

Conclusions
The goal of our current research is to investigate the effects of the embedded parameters on different state variables, like velocity, temperature, and the behaviour of the modeled equations under these parameters.The modeled equations are tackled with an analytical approach "homotopy analysis method" HAM.The embedded parameters and their effects are investigated and studied graphically.The motion of the suspended particles always has a great physical significance in heat transfer processes.Physically, the larger the Brownian motion parameter, the greater is the heat induction.Similar results have been observed for temperature profile with the larger values of Brownian motion.The effect of increasing values of Sc on boundary-layer thickness show a decrease in the thickness of the layer.The surface temperature increases with the increasing values of Pr, while an opposite effect is observed for unsteadiness parameter S, i.e., large values of S reduce the surface temperature.It is also observed that the temperature profile falls with large numbers of thermophoresis parameter Nt and vice versa.It is further observed that larger values of Nb reduce the mass flux, where Nt increases the mass flux.Higher values of Sc reduce the flux of mass, while it increases with increasing values of Sc.On the other hand, skin friction C f shows a decline in its behavior with larger values of the stretching parameters and λ, while the unsteadiness parameter S increases C f .In the solution procedure, the implemented technique convergence under the variation of physical parameters is observed numerically, which shows the reliability of our technique.
The central concluded points are as follows: • A comparative analysis for the stretching and unsteadiness parameters for the gradient of the velocity is discussed to observe the sensitivity of these parameters.

•
The temperature profile climbs up with larger values of Brownian motion parameter Nb.

•
Wth larger values of S, the thermal boundary layer thickness reduces.

•
Higher values of Pr increase the surface temperature, where an opposite effect is observed for unsteady parameter S, i.e., large values of S reduce the temperature of the surface.

•
It is examined that the heat profile decreases with increasing values of thermophoresis parameter Nt, and increases with small numbers.

•
The increasing values of Nb reduce the mass flux, where Nt increases the mass flux, while it rises with rising values of Sc.

•
The effect of Prandtl number Pr on concentration and temperature profile is analyzed and a similar decline is observed in both the profiles.

•
The convergence of the HAM method with the variation of the physical parameters is observed, and found its convergence more rapid as compared with other techniques.

Abbreviations
The following abbreviations and parameters with their possible dimensions stated here are used in this article:

Figure 1 .
Figure 1.Geometry of the physical model.

K
is the Prandtl number, Nt = τD w (T w −T ∞ ) νT ∞ represents thermophoresis constraint, Nb = τD B (C w −C ∞ ) ν represents the limitation of the Brownian motion, and Sc = ν D B denotes Schmidt number.All these parameters and numbers are well defined and explained briefly in literature.C f x and Nu x represent the local skin-friction coefficient and local Nusselt number respectively, and are defined as:

Figure 3 .
Figure3illustrates the effect of the stretching parameter .The velocity profile shows an increase with the increasing values of , because the lower plate always behaves directly to the flow fluid velocity.Physically, for > 0, the surface accelerating rises, < 0 decelerating the surface, while = 0 depict the random motion of the surface.Figure3breflects the sensitivity of , under smaller values of β and S.

Figure 4 . 8 Pr = 0. 6 Figure 5 .
Figure 4. (a) Impact of S on f (η), when β = 0.1, = 0.2 and (b) Impact of β on f (η), when β = 0.01, = 0.01.The impact of Pr on θ(η) is shown in Figure 5.An inverse relation has been observed between the temperature and the Prandlt number.Physically, for small values of Pr these fluids have larger thermal conductivity and vice versa.As a result, for larger values of Pr the thermal boundary layer declines.

4 Figure 12 . 4 Figure 13 .
Figure12.Impact of Nt on φ(η), when Sc = 0.7, = 0.3, β = 0.9, S = 0.5, Nb = 0.5, Pr = 0.6.The impact of Pr on φ is shown in Figure13.Larger values of the Prandtl number Pr cause the concentration to falls down.The information from the figure reveals that large values of Pr cause the concentration profile to fall down.Physically, the thermal boundary-layer vanishes with greater values of Pr and as a result the concentration profile falls.The same phenomenon is observed for the heat profile.

Table 2 .
It is observed that the increasing values of stretching parameters of Reiner-Philippoff fluid and λ decrease C f , while unsteadiness parameter S increases C f .The effects of Nb, Nt, Pr, Sc, and S on the Sherwood number Φ (0) are demonstrated inTable 3. It is observed that local Sherwood number values increases due to an increase in thermoporetic parameter Nt.Increasing values of Schmidt number Sc decreases the Sherwood number, while increasing unsteady parameter and Prandtl number decreases the Sherwood number.

Table 1 .
Variation in Nusselt number with different values of the parameters Nt, Pr, β, and S.

Table 2 .
Variation in skin friction with different values of the parameters , λ, and S.

Table 3 .
Variation in Sherwood number with different values of the parameters Nb, Nt, Sc, S and Pr.