Impact of Magnetohydrodynamics on Stagnation Point Slip Flow due to Nonlinearly Propagating Sheet with Nonuniform Thermal Reservoir

Sarhad University of Science and Information Technology, Peshawar 25000, KPK, Pakistan Department of Mathematics, Abdul Wali Khan University Mardan, Mardan 23000, KPK, Pakistan Department of Mathematics, Kohat University of Science and Technology, Kohat 26000, KPK, Pakistan Department of Mathematics, College of Science Al-Zulfi, Majmaah University, Al-Majmaah 11952, Saudi Arabia Department of Mathematics and Computer Science, Faculty of Science, Beirut Arab University, Beirut, Lebanon Department of Computer Science and Engineering, Soonchunhyang University, Asan 31538, Republic of Korea Department of Mathematics, College of Arts and Sciences, Prince Sattam Bin Abdulaziz University, Wadi Al-Dawaser, Al-Kharj 11991, Saudi Arabia


Introduction
It is, to some extent, understood that the present generation depends on the achievements of physical sciences which are based on production industries. Magnetohydrodynamics (MHD) boundary-layer flow on an elongating surface is important because of its frequent uses in industrial engineering and various production processes such as the aerodynamic squeezing of polymers, rolling at high temperature, cooling control technology, and glass fiber production. e magnetohydrodynamics can display certain characteristics in heat conductivity as it has both fluid as well as magnetic features. Raju et al. [1] investigated an extensive study of the least squares finite element method over a varying boundary layer explored the stagnation flow. Gorla [2] studied a viscoelastic (a nonclassical) liquid of stagnation movement in pulsating magnetic field. Gorla concluded that the shear stress coefficients are directly proportional to the magnetic field. Takhar et al. [3] investigated magnetohydrodynamic nonlinear stagnation point interface movement. Besser et al. [4] presented a technique of annihilating the magnetic field in the limit of MHD equations for a noncompressible ionized fluid having constant viscosity and resistivity as exchange parameters. Massoudi and Ramezan [5] examined the heat convective features of fluid (i.e., viscoelastic) at stagnation point flow. Ariel [6] studied multidirectional stagnation point movement of viscoelastic material. Mahapatra and Gupta [7] derived an exactly matched solution of Navier-Stokes equations which characterizes a uniform axisymmetric stagnation point flux towards an elongating surface. It is observed that the flow shows a boundary-layer arrangement when the velocity of the elongating sheet is smaller than the free stream velocity and a counter boundary layer is made when the starching sheet velocity is greater than the free stream velocity. Abel et al. [8] studied the influence of varying temperature source on magnetohydrodynamics heat convection in the fluid thin layer on a nonuniform elongating surface.
Yazdi et al. [9] evaluated the magnetohydrodynamics slip streaming on the nonuniform porous elongating sheet in the existence of a chemical process. e slip stream takes place if the distinctive size of flow regime is much smaller or the flow pressure is weak. Hsiao [10] has worked on an incompressible uniform MHD stagnation point movement of a second-order viscoelastic liquid and thermal conduction caused by a flow oscillating surface, and it is observed that viscoelastic liquid flow thermal effect is better than nonviscoelastic liquid flow thermal effect. Rasheed et al. [11] analyzed numerical and analytical investigation of thin-film nanofluid flow over an angular surface. Roşca et al. [12] presented an analysis for the uniform boundary-layer movement and heat transformation of a non-Newtonian fluid in the stagnation point along with unsteady propagating plane sheet of the unrestricted flow slip velocity and presented a dual (upper and lower branch) solution for certain variables. Dessie and Kishan [13] studied the magnetohydrodynamics over boundary-layer stream and heat transformation of fluid material with varying viscosity, permeable channel, heat reservoir, and viscous dissipation, respectively. Rasheed et al. [14] investigated the two-dimensional viscoelastic fluid with nonuniform heat generation over permeable stretching sheet with slip condition. Hassan [15] described the charge carrying viscous interfacelayer flow and thermal transfer. Shen et al. [16] studied the MHD-varied heat transfer flow near a stagnation point flow on a nonuniform stretching surface with velocity slip. Tufail and Ali [17] investigated the effects of fluid flow and heat convection due to a nonlinear stretching sheet. Gireesha et al. [18] examined the hydromagnetic heat conduction in a mixed viscous fluid over uniform stretching thermovariant surface and observed the effects of different relevant parameters on flow and heat transfer. Zaidi and Mohyud-Din [19] studied the convective transformation of heat and significance of MHD effects in different technologies.
Khan et al. [20][21][22][23][24][25] discussed the impact of various non-Newtonian fluid materials for wire coating analysis in the presence of magnetic field. Kuman et al. [26] investigated the characteristics of entropy generation in radiative flow of CNTs Casson nanofluid in rotating channels with heat source. Khan et al. [27] studied the analytical and numerical solutions of an Oldroyd 8-constant fluid in double-layer optical fiber coating as a coating material. Nandeppanavar et al. [28] studied heat conduction over a nonuniform elongating surface with unsteady heat reservoir and flexible temperature at boundary [28], and fluid layer flow is caused by a nonuniform elongating surface [29], magnetohydrodynamic stagnation point stream movement and heat convection [30], and seconddegree slip effects over flow movement [31]. Ibrahim and Ul Haq [32] studied the heat transformation in fluid and MHD stagnation point movement and observed that the thickness of heat conducting layer rises when the thermophoresis factor reduces with Prandtl number Pr.
By studying the foresaid investigations, it is found that there is no research on the stagnation point movement of a non-Newtonian fluid and thermal convection with thermal jump and momentum in the existence of an unsteady heat reservoir as a result of the moving sheet. In this analysis, we studied the said constraints on flow and heat transfer. Figure 1 demonstrates physical configuration of the considered problem. Here, we measured two-dimensional steady, flow, and transfer of heat analysis in case of an incompressible fluid in the occurrence of transverse magnetic field strength B(x) affected normally on flow moment which give us a unique form given as follows:

Mathematical Formulation
where n is a constant and x is a coordinate along the plate measured from the leading edge. e plate is moving inside or outside the origin with the velocity u W (x) � ax n in an exterior (in viscid) flow of the velocity u e (x) � ax n , where u and v are the corresponding velocity components in the x and y directions, respectively. Here, T W (x) is assumed as a temperature of the plate and ambient fluid is T ∞ which is the constant temperature. e governing equations of continuity, momentum, and energy equations are as follows: where ρ is the conductivity (electrical) of the fluid, α is the thermal diffusivity constant, and u and v are velocity components in the xy − plane, respectively. e initial and boundary conditions are as follows: We guess the velocity slip factor N 1 and the temperature slip factor S 1 change with x in the form where T 0 is the characteristics temperature parameter and P is the wall temperature parameter. It is concluded from against physical point of view, n should vary in the range 0 ≤ n ≤ 1. If n > 1, then N 1 and S 1 become singular at x chose to the leading edge of the plate. It is remembered that the boundary layer does not start at x � 0 but starts in the vicinity of the leading edge of the plate [33]. erefore, the solution for n > 1 is realizable from the mathematical point of view. C P is the specific heat at constant pressure, σ is the electrical conductivity, B 0 is the applied magnetic field, μ is the viscosity, T is the temperature, k is the thermal conductivity of the fluid, and q‴ is the rate of nonuniform heat generation/absorption coefficient and defined as where A * and B * are the conditions of space-and temperature-dependent heat generation/absorption, respectively. We notice that A * > 0 and B * > 0 for the internal thermal generation and A * < 0 and B * < 0 for the internal thermal absorption [28].
Physical quantities for local skin-friction coefficient C f and local Nusselt number Nu x in this problem are defined as , where skin friction (shear stress) along the plate is τ w , and wall heat q w is given by From equations (8), (12), and (13), we acquired ��� where Re x � (u e (x)x/v) denotes the Reynolds number.

Solution by HAM (Homotopy Analysis Method)
e governing nonlinear partial differential equations (3) and (4) are converted into nonlinear ordinary differential equations (9) and (10). To find the solutions of equations (9) and (10), five boundary conditions are required: three on equation of motion and two on equations of temperature, respectively. But here f ″ (η) and θ ′ (η) are missing boundary conditions; hence, solving the boundary value problem of equations (9) and (10) is difficult. erefore, in the boundary conditions given in equation (11) we replace infinity to the finite value. Equations (9) and (10) with boundary condition (11) are solved analytically by HAM and numerical by BVPh2-midpoint methods. e optimal HAM [35][36][37] gives better results compared with perturbation techniques and other conventional investigative techniques. Firstly, the optimal HAM gives us a remarkable flexibility to pick the equation type of linear subproblems. Secondly, the optimal HAM works regardless of the possibility that there does not exists any small/large physical parameter in determining equations and boundary/initial conditions. Particularly, unlike perturbation and other analytic techniques, the optimal HAM gives us an advantageous approach to guarantee the convergence of a series solution by presenting the supposed convergence control parameter into the series solution.

Convergence of the Optimal HAM Method
e auxiliary parameters f and θ have a leading purpose of controlling the convergence of homotopic solutions. To get convergent solutions, we take the suggested values of these parameters. For this reason, residual errors are noticed for the momentum, and thermal energy equations by initiating the expressions are given as e convergence of the parametric values computed through optimal HAM is listed in Table 1 using the values of the parameters β � 1.2, λ � 0.5, Pr � 2.0, ϕ � 0.1, and c � 0.1, while the error decay for 10th-order approximation is shown in Figure 2.
Here, Δ t m � Δ f m + Δ θ m denotes the total discrete square residual error which is used to obtained optimal convergence control parameters. Table 2 shows the comparison of HAM, BVPh2-midpoint method, and published work reported by Bejan [34], and good agreement has been found.

Discussion on Results
In this paper, an analysis is introduced to investigate magnetohydrodynamic slip flow and thermal convection of the stagnate point of an elongating surface along with nonuniform oscillating plane along a free flow. e momentum and heat convection equations for the nearest neighbor layer of interface have been explained analytically and derived for various analytical equations for temperature profile. e outcomes of the semianalytical computations are displayed in different figures. e characteristic parameters which are investigated in this paper are power-law index n, boundary-layer temperature P, Prandtl number Pr, magnetic field parameter M, and spatial dependent heat reservoir parameters A * and B * , respectively. We now continue to discuss the results. Figures 3 and 4 demonstrate the stream velocity profiles for the various values of λ (oscillating parameter) and β (velocity slip parameter), respectively. It has been observed that the rising of λ and β causes reduction in the width of the boundary-layer movement. Figure 5 describes velocity and flow profiles of various input of power index variable n while the other parameters are kept fixed. It is observed that rising the power index variable, the stream movement and velocity rise, which cause an increment in the width of the momentum boundary wall layer. Figure 6 exhibits the velocity sketch for various input  of M. From here, it is detected that velocity along with the boundary-layer thickness reduces when M is increased. is may be due to the application of electromagnetic force (Lorentz force). Figures 7 and 8 depict the velocity profiles for the various inputs of λ and β. In these figures, we can see that as the value of λ advances, the thickness of velocity increases but we observed a reverse result when the velocity slip parameter β increases. Figures 9 and 10 show the skin-friction profile for different values of λ and β. In these profiles, it is clear that the        Figures 11 and 12 illustrate the temperature profile for various input values of power index n and some other physical variables which reveals that the thickness of the thermal-conducting boundary layer decreases with increasing the nonlinear stretching parameter n. Because of electromagnetic (Lorentz) force, the thickness of the heat convective boundary layer increases as the magnetic field strength increases; henceforth, Figure 13 is sketched for the temperature dependence upon the magnetic field. Figures 14 and 15 show the temperature profiles for different inputs of Pr. ese figures illustrate that, on increasing the Prandtl number Pr, the thickness of the thermal-conducting boundary layer decreases; hence, the temperature decreases.
Temperature profile for numerous values of positiondependent heat reservoir parameter A * and temperaturedependent heat reservoir parameter B * , respectively, is shown in Figures 16 and 17, respectively, while the wall temperature profile is given in Figure 18 which illustrates that the temperature on the boundary layer decreases as a result of rising the power index p.  parameters λ and β, respectively. It is observed that, by increasing both λ and β parameters, the thickness of the thermal boundary layer increases.    increasing value of the oscillating and velocity slip parameters, respectively, and the temperature gradient drops when the temperature slip parameter rises.

Conclusions
In this review, we have attempted to provide a glimpse of what we have studied from the stagnate point of magnetohydrodynamic flow movement and thermal convection in the presence of velocity slip, thermal slip, and nonlinear heat reservoir. We deduced some significant observations from these results given as follows: (i) By rising the power index, the velocity parameter rises, while in existence of magnetic field, a converse result is observed. (ii) Whether magnetic field exists or not, the velocity rises by increasing the power-law index. (iii) If the magnetic field value is zero, the temperature rises as the power index increases while falls when magnetic field is nonzero. (iv) As the Prandtl number increases, the temperature falls irrespective of the magnetic field. (v) In presence of uniform magnetic field, the temperature falls down as the parameters A * , B * , and p increase. Also, the temperature rises when there is a variable magnetic field.   (vi) e oscillating parameter λ decreases flow movement, velocity, variation of temperature in normal direction, skin friction, and heat transfer boundary layer. (vii) e velocity slip parameter increases the thermal boundary layer and gradient of temperature but decreases the flow movement, velocity (momentum), and skin friction of the boundary interface layer. (viii) e temperature slip parameter impacts merely on thermal interface layers which decreases the heatconvective boundary layer and so as reduces the gradient of temperature.