Numerical Treatment for the Three-Dimensional Eyring-Powell Fluid Flow over a Stretching Sheet with Velocity Slip and Activation Energy

In this manuscript, a computational paradigm of technique shooting is exploited for investigation of the three-dimensional EyringPowell fluid with activation energy over a stretching sheet with slip arising in the field of fluid dynamics. The problem is modeled and resulting nonlinear systemof PDEs is transformed into nonlinear systemofODEs usingwell-known similarity transformations. The strength of shooting based computing approach is employed to analyze the dynamics of the system. The proposed technique is well-designed for different scenarios of the system based on three-dimensional non-Newtonian fluid with activation energy over a stretching sheet. Slip condition is also incorporated to enhance the physical and dynamical analysis of the system. The proposed results are compared with the bvp4C method for the correctness of the solver. Graphical and numerical illustrations are used to envisage the behavior of different proficient physical parameters of interest includingmagnetic parameter, stretching rate parameter, velocity slip parameter, Biot number on velocity, and Lewis number on temperature and concentration.


Introduction
witnessed that, at moderate stresses, the fluid velocity is exponentially dependent on the stress, but as the stress level is enhanced, velocity behaves linearly in the flow.Patel et al. [2] reported that the Eyring-Powell model is useful and has significant benefits in comparison with the power law model.Fluid flow close to a dynamic plate for the Eyring-Powell with the help of three different techniques has been studied by Sirohi et al. [3].Considering the complexity of the fluid time scale with the aid of Powell-Eyring fluid model, the impact of the very small and large shear rate viscosities is reported by Yoon et al. [4].Nadeem et al. [5] explained the peristaltic flow by using Eyring-Powell fluid in their study.They concluded that, by enhancing the non-Newtonian Eyring-Powell flow parameter, increase in the peristaltic pumping region is observed.Time independent flow of incompressible boundary layer flow of Eyring-Powell nanofluid past shrinking sheet has been investigated by Motsa et al. [6].Jayachandra et al. [7] and later Hayat et al. [8] analyzed the double stratification effects using Eyring-Powell.
The field of nanotechnology has gotten much attention among the engineers and scientists due to its vast applications in industrial and medical sciences like ceramics and drug delivery, etc. Nanofluids are modeled because of the Advances in Mathematical Physics suspension of nanoshaped elements in the base liquid.Most commonly used nanoparticles are , , , , and their oxides.The nanofluid flow passing through stretching sheet near stagnation point was analyzed by Mustafa et al. [9].Makinde et al. [10] analyzed nonisothermal boundary layer nanofluid fluid because of linear stretching surface.Sheikholeslami et al. [11] reported the heat transfer properties for nanofluid flow in the channel.Goodarzi et al. [12] reported nanofluids flow for laminar and turbulent case in a channel of narrow cavity, by using two types mixture model.Malvandi et al. [13] presented the flow for multiple convection nanofluid underneath the transverse narrow channel by using modified Buongiorno's model.Kuznetsov et al. [14] reported natural convection nanofluid flow past a vertical plate.The 3D flow of second grade nanofluid was analyzed by Hayat et al. [15]; they used an exponentially stretching surface and studied the flow due to the thermal radiation effect and heat source/sink.Hedayati et al. [16] analyzed the convection enforced nanofluid flow in a narrow channels using asymmetrically provided heat.Using the heat source/sink, three-dimensional Maxwell nanofluid boundary layer flow was studied by Hayat et al. [17].Ellahi et al. [18] discussed entropy induced nanofluid flow using the behaviors of nanoparticles geometry.Latiff et al. [19] reported the bioconvection boundary layer micropolar nanofluid flow over a variable size sheet using diverse slip impacts.Uddin et al. [20] analyzed bioconvection nanofluid flow past a dynamic plate using the variations of Stefan injection along with diverse slip.
Several models for the slip at the boundary flow have been investigated in the span of previous few decades.Some intensively employed models are Maxwell's boundary-slip flow model of order first [21], slip boundary flow model of second order [22], Fukui-Kaneko (FK) slip boundary flow model [23], and slip boundary flow model of fractional order [24].All of these models have been originated with the help of kinetic theory for gases.The Fukui-Kaneko slip model is adequately precise; however, it very much relies on the shape of the flow and dependent upon the flow parameters involved.These velocity slip models of first grade, second grade, and fractional grade work well when the Kundsen number is taken between zero and one.The aforementioned models are extensively being employed in science and engineering problems.Currently, Lin Wu [25] studied enhanced slip model of second grade that is considered accurate and reliable in comparison with the FK model.This model is considered superior over previous slip models, because of its applicability for all Kundsen number.For more insightful understanding of the flow behaviors, several researchers worked on Wu's model using various physical parameters.Furthermore, Ibrahim and Shankar [26] introduced the slips of three types: velocity, solutes, and thermal, due to flow of MHD nanofluid through a stretching sheet.
Zhang et al. [27] developed a new scheme to optimize the GRNN parameter and also introduced a new nonlinear function fitting scheme.Duan et al. [28] presented incompressible SPH algorithm by soothing the constant invariant density scenario for fluid of constant density.Qing et al. [29] by using finite volume technique explored the impacts of water entry angle through cavitation flow region.
The activation energy is the least required energy that reactants have to obtain, so that a chemical reaction can be anticipated.The phenomenon of mass transfer happens because of concentration variance of types in a mixture.The types that can vary the concentration in a mixture flow from large concentration section to less concentration section.The activation energy has many advantages in chemical engineering industry, food processing, mechanics of water, oil emulsions and geothermal reservoirs, etc. Bestman [30] studied the flow using naturally convected binary mixture in a permeable medium along with activation energy.Makinde et al. [31] reported the flow of natural convection with the impacts of reaction, for order n and activation energy through numerical computations.Maleque [32] examined, with the aid of activation energy, the reactions of exothermic or endothermic on variant convective flows.Hongchun et al. [33] presented meshless barycentric interpolation collocation technique for finding the solution of PDEs which are encountered in many physical problems.Zhang et al. [34] presented numerical technique for fraction control problems using Chebyshev polynomials.Yao et al. [35] developed a mathematical model for investigating the loss of root stone due water flow in dam structure.Hemeda et al. [36] proposed an iterative technique along with integral iterative scheme for linear and nonlinear Fokker-Planck equations.Awad et al. [37] used the time dependent rotating binary fluid flow over a suddenly stretched sheet using the enhanced Arrhenius function.Abbas et al. [38] also studied flow of Casson nanofluid in the locality of stagnation point by using dual chemical reaction and activation energy effects.They modeled mathematically the flow system using spectralcollocation quasi-linearization method.Shafique et al. [39] used numerical approach and give the idea of rotating viscoelastic flow having species of chemically reaction with activation energy.
The novel contributions of the study are presented as follows: (i) A novel investigation is presented for modeling and analyzing numerically by shooting method for the three-dimensional Eyring-Powell fluid system with activation energy over a stretching sheet arising in the field of fluid dynamics (ii) The strength of similarity transformation is exploited to transform the nonlinear PDEs of fluidic model into system of nonlinear ODEs.(iii) The proposed results are compared with the stateof-the-art counterpart to prove its effectiveness and correctness.(iv) Numerical and graphical illustrations are utilized to ascertain the value of the scheme by visualizing and analyzing the performance of proficient physical quantities on velocity, temperature, and concentration profiles of the system.
The rest of the paper is organized as follows: mathematical formulations of the problem are described in Section 2; in Section 3, results of numerical experiments are provided, while the conclusion along with future research studies is provided in the last section.
In the current investigation, the three-dimensional Eyring-Powell fluid flow near a stretching surface with velocity slip and activation energy is discussed.To discover the full insight of the behavior of different parameters on velocities, energy, and concentration utilities, a well-known shooting technique along with fourth-order Runge-Kutta technique has been employed.

Mathematical Formulation
In this section, the detailed mathematical analysis for time independent incompressible three-dimensional Eyring-Powell fluid boundary layer flow, in the presence of activation energy, has been investigated.The region  > 0 that occupies the flow domain is considered.The sheet is extended in two directions with fixed origin.The velocities of the sheet are    =  and V   =  along the longitudinal and transverse directions.Coordinate scale and geometrical interpretation of the problem are exhibited in Figure 1.The velocity component  is along -axis, V is along -axis, and  is along -axis.  represents nanoparticles concentration at sheet and  ∞ is showing the ambient concentration.  and  ∞ are representing the temperature of the convective region and ambient surface.The equations of motions along energy and concentration of the Eyring-Powell model with the effects of activation energy can be expressed as ( The boundary conditions of the above-mentioned problem are defined as We define the below similarity transforms: By employing similarity transformations, the nondimensional forms of (1)-( 5) are transformed as () + Pr ( () +  ())   () +   ()   () With the help of similarity transformation (6) can be described as Various unitless parameters arising in the aforementioned equations are given as The local skin friction coefficients can be defined in  and  directions as Now we can define the skin frictions   and   in  and  directions; the thermal and mass heat fluxes are   and   defined as follows.
The dimensional local Nusselt number and Sherwood number, respectively, can be interpret in the form as Using ( 14) to (16), one can easily get Here   =   /] is the Reynolds number which is relied upon on the stretching velocity.

Solution Methodology
The solutions for coupled nonlinear ODEs ( 8)- (11) with the aid of boundary conditions (12) are determined by employing the famous shooting method.For the determination of numerical solution, the domain of the problem has been taken as [0, 8] instead of [0, ∞].It is observed that the solution for  > 8 has disregarded variations.The initial value problem is obtained from the above boundary value problems.A system of initial value problem has been transformed, in which  by  1 ,  by  4 ,  by  7 , and  by  9 have been used.The resulting equations are written as where  > 0 is a small positive number.In the present study, the numerical outcomes are obtained for  = 10 −6 Table 1 shows the deep analysis of presently studied results of   (0) and   (0) against the Hartmann number  and stretching rate ratio parameter  with those of Freidoonimehr et al. [40] and Hayat et al. [41].The computed results in Table 1 show a very good comparison with the published numerical results.Table 2 shows the comparison of the numerical results of   (0) and   (0) against the stretching rate ratio parameter  with those of Wang [42] and Freidoonimehr et al. [40].A fabulous agreement is noticed between the computed results and the published numerical results in Table 2.For further validation and strengthening of proposed solver, the obtained numerical results have been compared with the MATLAB built-in solver bvp4c.
Ref. [35] Ref. [33] shooting bvp4c Ref. [35] Ref. [33] shooting bvp4c   magnetic parameter relies upon the Lorentz force which provides resistance to the flow.By enhancing the Hartmann number, the well-known Lorentz force increases, due to this reason decline in all the component of velocity is observed in the fluid.Moreover, along with this behavior, a small enhancement is observed in the temperature along with concentration distribution by rising the magnetic number .This demonstrates that it is very clear that the transport phenomena are opposed by the transverse magnetic field.It would be significant to describe here that the huge resistances on the fluid elements create heat generation in the flow field, due to this magnetic field increases vertically.Figures 6-9 narrates the velocity slip parameter  effects on both of the fluid velocity components, temperature, and concentration profiles.Here in this study, the variation of slip parameter is 0 <  < 1.By enhancing the values of , decrement is noticed in the fluid velocity components, while the increment is seen in the temperature and concentration profiles.In another sense, the effect of slip becomes stronger, when the smaller quantity of flow is taken outside and it moved ahead in all directions of the flow.By increasing,  = 0 creates a decline in the saturation of the fixed surface due to the boundary layer in both the axial and transverse directions, which results in decline in the boundary layer thickness of momentum, as reduction in the flow is observed with the enhancement in slip so that the skin friction is further decreased from the wall.As described earlier, like behavior is viewed for temperature profile along with slip effect for the nanoparticle concentration field.The effects stretching parameter  on the different components of fluid velocity, profiles of temperature, and concentration are demonstrated in Figures 10-13.By the definition of stretching parameter, when  = 0 shows the case of nonbidirectional stretching sheet, due to this fact,   tend to zero and improvement in the two-dimensional case is achieved.Moreover, in the case of  = 1, the stretching parameter becomes similar in the axial and transverse directions.It is worth mentioning here that, by increasing,  signifies the sheet velocity in  direction; that is,   increases; or decrements are noticed in  direction, that is,   .Generally, as  enhances from zero, then the lateral surface moves in the  direction.In addition to this, increasing in the values of  and decreases in the thermal boundary layer thickness are viewed, because boundary layer thickness and concentration profile are reduced.Figures 14 and 15  in Figure 16.In boundary layer, region  represents the comparative contribution of the rate of thermal diffusion to the species diffusion.Lewis number is increased due to the species boundary layer decline and the concentration distribution tends to zero.The behavior of the temperature profile against thermophoresis parameter is presented in Figure 17.It is noticed that if thermophoresis parameter  is increased, enhancement in the thermal and concentration boundary layer thicknesses is achieved.The reason is behind the fact that submersion of the nanoparticles improves the thermal conductivity of the fluid; due to this an increment is seen in the temperature profile.This is quite evident from Figure 18 that, by increasing the dimensionless activation energy , increment is seen in the nanoparticle concentration profile.The reason behind this fact is Arrhenius function.This function reduces when the activation energy increases, leading to start the generative chemical reaction.Figure 19 portrays the reaction rate constant  versus mass fraction field .It shows that, by increasing  results, the mass fraction field  reduces.The reason behind this fact is the presence of the huge concentration gradient at the boundary wall.

Conclusion
A new numerical approach based on shooting method is designed effectively for the problem arising in the field of fluid dynamics represented with the three-dimensional  (iii) The activation energy strength enhances the mass fraction field, where as the wall mass flux declines.
(iv) Rise in the temperature and concentration profiles is witnessed with the increase of Biot number.
In future one may investigate stochastic numerical paradigm based on evolutionary and swarming optimized spline method for the superior numerical treatment of given fluidic system Kinematic viscosity : Brownian motion factor : A c t i v a t i o ne n e r g y : Material parameter : Thermophoresis factor : R e a c t i o nr a t ec o n s t a n t : F i t t e dr a t ec o n s t a n t (), (): Velocity similarity functions : T e m p e r a t u r e   : Convective surface temperature  ∞ : A m b i e n tt e m pe r a t u r e , V, : Componentsofvelocity , , : Cartesian coordinates : Thermophoresis factor : Thermal diffusivity : Slip parameter  0 : S l i pl e n g t h : Similarity parameter (): Temperature distribution : Stretching parameter : F l u i dd e n s i t y : Thermal conductivity , : Nondimensional fluid parameter (): Concentration.

From
the last few decades, considering the non-Newtonian nature of different physiological fluid and their use in the industry, researchers are paying their attention in the fluid mechanics field on non-Newtonian fluid models.Examples of this category of fluids are Reiner-Philippoff fluid, Casson fluid, micropolar fluid, Prandtl-Eyring fluid, Carreau fluid, power law fluid, Prandtl fluid, and Eyring-Powell fluid.Powell et al.

Figure 1 :
Figure 1: Geometrical representation of the problem.

Figures 2 - 1 M=0Figure 3 :
Figures 2-5 demonstrate the behavior of the velocity components   and   , temperature profile , and concentration distribution  against the Hartmann number .The drag force famous as Lorentz force is obtained by executing the magnetic field vertically through electrically conducting fluid.It is the capability of Lorentz force to make the flow slow which is passing through stretching sheet.With the enhancement in Hartmann number , velocity profile reduces in the directions of  and .As stated above, the

Figure 8 :Figure 9 :
Figure 8: Demonstrates the impact of  on temperature profile.

Table 2 :
Comparative analysis of shooting and bvp4c results with those of Wang