Analytical Investigation of Magnetic Field on Unsteady Boundary Layer Stagnation Point Flow of Water-Based Graphene Oxide-Water and Graphene Oxide-Ethylene Glycol Nanofluid over a Stretching Surface

This study explains the effect of magnetic field of the stagnation point flow of a water-based nanofluid graphene oxide-water (GO-W) and graphene oxide-ethylene glycol (GO-EG). Heat transfer analyses are discussed by converting the given partial differential equation into a nonlinear ordinary differential equation using the similarity transformation and solved using an approximate analytical method, namely, the optimal homotopy analysis method (OHAM), to obtain an approximate analytical solution of the nonlinear problem that analyzes the problem. The BVPh 2.0 package function of Mathematica is used to obtain the numerical results. The results of important parameters such as the magnetic field parameter, unsteady parameter, stretching parameter, Prandtl number, Eckert number, and kinematic parameter for both velocity and temperature profiles are plotted and discussed. The convergence control parameter of the approximate analytical method is obtained up to the 25th iteration using the BVPh 2.0 package. The skin friction coefficient and Nusselt number are explained in tabular form.


Introduction
In the history of fluid mechanics, the derivation of the boundary layer equation and its solution using similarity transformation is an important area for researchers. With the help of boundary layer theory, both Newtonian and non-Newtonian fluids can be modeled. e results obtained using the boundary layer theory are more similar to the experimental results. e industrial application of non-Newtonian fluids is higher than that of Newtonian fluids because of the use of non-Newtonian fluids in petroleum drilling polymer engineering and manufacturing of food and paper. Recently, the more active area of graphene oxide as a nanofluid is a rotating disk model problem. In a hexagonal lattice, graphene is a single layer of carbon atoms. Because graphene is a single layer of carbon and graphene oxide is a two-dimensional material, compared with graphene, researchers are interested in graphene oxide. In 1859, graphene oxide was introduced by Sir Second Baronent Benjamin Colline using the chemical method, Brodie, through the oxidation of bulk graphite with potassium chlorate and nitric acid. Lerf and Klinowsky used a modified Hummers' method to study the natural structure of graphene oxide. Graphene oxide has some key applications in engineering and industry, such as computer storage systems, turbomachinery hard disk jet motors, and turbine systems. It is known that the thermal conductivity of the base fluid is lower than that of the nanofluid, and to increase the heat transfer ratio, a different method is used; for example, the reduced heat transfer time and heat exchanger size can be minimized. It is known that the thermal conductivity of ethylene glycol is greater than that of the base fluid; therefore, ethylene glycol is used as a cooling fluid and antifreezing agent.
Many nonlinear relations were observed for the stress and the rate of strain for non-Newtonian fluids. It is difficult to express all the properties of a non-Newtonian fluid in a single equation. e flow owing to the stretching sheet of the boundary layer of non-Newtonian fluids has important applications in several manufacturing processes, such as molten polymers through a slit die for the production of plastic sheets, hot rolling wire and fiber coating, processing of food stuffs, metal spinning glass fiber production, and paper production. In mechanical applications, non-Newtonian fluid models play an important role. Numerous materials have incredible significance in ventures; for example, liquid plastic, nourishment items, divider paint, oils, lubricant oils, and penetrating mud have non-Newtonian fluid conduct. According to scientists, non-Newtonian liquids are the most pertinent to exploration. One of the precarious subclasses is Sisko fluid, which has more centrality as an event in many designing activities. e common sense application has delivered an enthusiasm for examining the feasibility of differential conditions overseeing streams of non-Newtonian fluids, which have various uses in the field of design, applied mathematics, and computer engineering. Munir et al. [1] used an extended sheet to study the bidirectional stream of the Sisko fluid. Olanrewaju et al. [2] used a level plate to study the Sisko precarious free convective liquid stream in the presence of heat exchange. Khan et al. [3] considered the Sisko fluid, which is consistent with the heat movement in an annular pipe. Khan et al. [4] used an extended surface to study the Sisko limit-layer liquid course. Patel et al. [5] studied a laminar Sisko fluid limit stream. Darji et al. [6] studied the Sisko fluid precarious common convective limit-layer stream. Siddiqui et al. [7] proposed a Sisko fluid film on a vertical belt for a waste framework.
Khan et al. [8] used an annular pipe to study the Sisko fluid. Sar et al. [9] reported the limit-layer conditions of a Sisko fluid. Marinca et al. [10] examined Maxwell fluid courses using permeable extending plates. Moallemi et al. [11] discovered a definite answer for the Sisko fluid stream in a pipe. Hayat et al. [12] investigated Maxwell nanofluids. Malik et al. [13] investigated an Erying-Powell nanofluid through an extension of the blended convective stream of magnetohydrodynamics (MHD). Nadeem et al. [14] investigated Maxwell fluid transfer through a vertical extending surface in the event of nanoparticles. Raju et al. [15] investigated an MHD nano non-Newtonian fluid through a cone with free convective heat and mass exchange. Rokni et al. [16] investigated a nanoliquid stream with heat exchange through plates. Nadeem et al. [13] investigated the numerical arrangement of a nanofluid stream on an extended sheet. Shehzad et al. [17] investigated the Jaffrey nanofluid MHD stream in the event of convective limit situations. Sheikholeslami et al. [18] used a nanofluid stream with an attractive field in the presence of heat exchange. Mahmoodi et al. [19] investigated a nanofluid stream with heat exchange for cooling applications. Recently, Shah et al. [20,21] investigated nanofluid streams with the effects of heat radiation and lobby current with a turning framework. e most recent theoretical study and investigation of Sheikholeslami on nanofluids utilizing disparate marvels with current application, assets, and properties with the use of different methodologies can be considered in [22][23][24][25].
Rehman et al. [26] studied analytically the Marangoni convection of a thin-film spray of CNT nanofluids using a stretching cylinder. In this study, the effect of a magnetic field on the unsteady boundary layer stagnation point flow of a water-based nanofluid graphene oxide-water (GO-W) and graphene oxide-ethylene glycol (GO-EG) with stretching surfaces is explained. Rehman et al. [27] studied analytically a thin-film unsteady nanofluid using a stretching sheet. Gul et al. [28] discussed the performance of graphene oxide nanofluids flowing in an upright channel through a permeable medium. Gul et al. [29] discussed the impact of Marangoni and thermal radiation convection on a graphene-oxide-water-based and ethylene-glycol-based nanofluid. Next, Gul et al. [30] discussed integer and noninteger order studies of GO-W and GO-EG nanofluids. Gul et al. [31] used rotating disks to study the stable dispersion of a graphene oxide nanofluid analytically.
In this paper, we will explain the approximate analytical solution of the unsteady boundary layer stagnation point flow of water-based GO-W and GO-EG nanofluids with stretching surfaces in the presence of a magnetic field, and variable thermal conductivity is obtained. e given partial differential equation is converted to a nonlinear ordinary differential equation using the similarity transformation and solved using the approximate analytical method (OHAM). Liao et al. [32] used this method for the first time and showed that the convergence control parameter of this method is faster than the approximate solution of the nonlinear problem as compared with other analytical methods. Using this method, we found a series solution of the velocity and temperature profiles.
e impacts of different parameters on both the velocity and temperature profiles were plotted and discussed.

Mathematical Formulation
Let us assume an unsteady two-dimensional incompressible flow of nanofluids over a stretching sheet. e stretching sheet velocity in x direction is denoted by U ϖ � bx/1 − ct, where b and c are constants, both of which are non-negative and t represents time. B 0 represents a uniform magnetic field, which was assumed from the Faradays. B(t) represents the magnetic field and is defined as represents the temperature distribution on the surface that changes in the x direction ( Figure 1).
Under the above conditions, the boundary layer equation for unsteady flow toward the stretching surface can be written as follows: e temperature distribution in the flow field is as follows: (3) e temperature is denoted by T, k represents the thermal conductivity, ρ is the fluid density, C p is the specific heat, and the dynamic viscosity of the nanofluid is denoted by μ s . u and v are the velocity components along the x and y directions, respectively, x denotes the distance along the sheet, y represents the distance perpendicular to the sheet, the velocity of the stagnation point is denoted by U s � ax with a > 0, where ] represents the kinematic viscosity of the fluid, the magnetic field is denoted by M and is defined as M � σB 2 0 /ρ nf , the unsteady parameter is denoted by S and is defined as S � c/b, Eckert number is defined as e boundary condition for velocity component is given as follows: e stretching velocity of the sheet is denoted by U w � cx, where c is the stretching constant. e boundary conditions for the temperature distribution were as follows: In equation (5), T w represents the temperature at the sheet and T ∞ represents the free stream temperature, which is assumed to be constant. e stream function is defined as follows: e similarity transformation is defined as follows: , and e similarity transformation from equation (7) is used in equations (2) and (3) to convert the PDE into a nonlinear ordinary differential equation.

Method of Solution
e nondimensionless equations (8) and (9) are solved analytically by an approximate analytical method, namely, the OHAM, which is given as follows: where L represents a linear operator, x represents an independent variable, g(x) represents an unknown function, N represents a nonlinear operator, and B(u) represents a boundary operator. Using this method, we first obtain a family of equations: where p represents an embedding parameter and lies in [0, 1], H(p) represents a nonzero auxiliary function for p ≠ 0 and H(0) ≠ 0, and ϕ(x, p) is an unknown function. e initial estimates obtained from equations (8) and (9) are as follows: with constant properties

Mathematical Problems in Engineering
For equations (7) and (8), the average squared residual error is obtained from the following equation:

Result and Discussion
is section explains the consequences of different nondimensional parameters such as kinematic viscosity, magnetic field, stretching parameter, Prandtl, Eckert number, and unsteady parameter on velocity and temperature distribution (], M, A, Pr, Ec, and S). e thermophysical properties of the base fluid and nanofluid are presented in Tables 1 and  2 Tables 3  and 4. Table 3 shows that the skin friction coefficient is a decreasing function of the magnetic parameter and unsteady parameter; that is, by increasing the value of the magnetic parameter and the unsteady parameter, the skin friction coefficient decreases in both GO-W and GO-EG. Table 4 shows that the Nusselt number coefficient is a decreasing function of the Prandtl number and Eckert number; that is, higher values of the Prandtl number and Eckert number decrease the Nusselt number coefficient in both GO-W and GO-EG. Tables 5 and 6 show the convergence of the given approximate analytical method for both GO-W and GO-EG nanofluids up to the 25th iteration, and it is observed from Tables 5 and 6 that with an increase in the number of iterations, a strong convergence is attained. Figure 2 shows the influence of the magnetic field parameters on the velocity profile. e relationship between f′(η) and M is inversely related. e velocity profile is a decreasing function of the magnetic field parameter, that is, a higher value of the magnetic field parameter decreases the velocity profile. With the increase in magnetic field resistance, certain types of forces known as Lorentz force are produced, which oppose the motion of the fluid particles in both GO-EG and GO-W; therefore, the velocity profile decreases, as presented in Figure 2. Figure 3 shows the influence of the unsteady parameters on the velocity profile. e values of f ′ (η) and M were inversely related. e velocity profile is a decreasing function of the unsteady parameter; that is, a higher value of the unsteady parameter decreases the velocity profile. By increasing the unsteady parameter, the motion of the fluid particles in both GO-EG and GO-W decreases, resulting in a decrease in the velocity profile, as shown in Figure 3. e influence of the stretching parameter on the velocity profile is shown in Figure 4. e relation between the velocity field f ′ (η) and the stretching parameter A is direct. e velocity is the increasing function of the stretching parameter; that is, the increasing values of the stretching parameter A increase      the velocity field f ′ (η). Physically, by increasing the stretching parameter A, the position of the fluid particle changes; thus, the moment of the particle increases, and the velocity field is enhanced by increasing the stretching parameter A, as shown in Figure 4. e influence of the kinematic viscosity on the velocity profile is shown in Figure 5. e relationship between the velocity profile and kinematic viscosity is inverse, or the velocity profile is a decreasing function of the kinematic viscosity; that is, higher values of the kinematic viscosity α decrease the velocity profile f ′ (η), as shown in Figure 5. By increasing the kinematic viscosity, viscous forces that oppose the motion of fluid particles in both GO-EG and GO-W are produced; therefore, the velocity profile decreases with increasing kinematic viscosity. Figure 6 shows the effect of the magnetic field on the temperature profile. It can be observed that by increasing the magnetic field, the temperature profile increases or the temperature profile is a decreasing function of the magnetic parameter. Such a state occurs because of the increase in resistance forces, which opposes the motion of fluid particles in both GO-EG and GO-W, resulting in friction forces that increase the temperature profile. e influence of the Prandtl number on the temperature profile is shown in Figure 7. From Figure 7, it can be observed that the temperature profile is a decreasing function of the Prandtl number; that is, a higher value of Pr decreases the   Figure 7.

, respectively. e influence of skin friction and Nusselt number is presented for both GO-W and GO-EG in
is effect is because the thickness of the momentum boundary layer is greater than that of the thermal boundary layer, and as a result, a higher value of the Prandtl number decreases the thermal boundary layer. e influence of the Eckert number on the temperature profile is shown in Figure 8. From Figure 8, it can be seen that the temperature profile is an increasing function of the Eckert number; that is, the increasing value of Eckert number increases the temperature field, and thermal diffusion is higher viscous diffusion; therefore, the improved value of Ec increases the temperature profile. e relationship between Nusselt number, Prandtl number, and Eckert number is shown in Figure 9. We see that Nusselt number is the decreasing function of Prandtl number and Eckert number in both nanofluids. Figure 10 shows the relationship between skin friction coefficient, unsteady parameter, and magnetic field parameter. We see that skin friction coefficients are the decreasing function of unsteady parameter and M (magnetic field parameter) in both nanofluids.

Conclusion
In this paper, we explain the approximate analytical solution of the unsteady boundary layer stagnation point flow of water-based nanofluids GO-W and GO-EG with stretching surfaces in the presence of a magnetic field, and variable thermal conductivity is obtained. e given partial differential equation is converted to a nonlinear ordinary differential equation by using the similarity transformation and solved by an approximate analytical method, namely, the OHAM. e practical application of this study is to produce cooling and heating effects, which also have significant applications in engineering and industrial machinery; the key objective of which is the circular pressure gradient and stability of centrifugal forces. e numerical results were obtained using the BVPh 2.0 package function of Mathematica.
e results of important parameters such as the magnetic field parameter, unsteady parameter, stretching parameter, Prandtl number, Eckert number, and kinematic parameter were obtained. e effects of the skin friction coefficient and Nusselt number coefficient are explained in table form. e outputs obtained were as follows: (i) e velocity profile is the decreasing function of the unsteady parameter S; that is, it decreases the velocity field by increasing the unsteady parameter (ii) e velocity profile is an increasing function of the stretching parameter A; that is, the velocity profile   Mathematical Problems in Engineering increases with an increase in the stretching parameter (iii) e temperature profile is the decreasing function of the Prandtl number; that is, the temperature profile decreases with the increasing value of the Prandtl number (iv) e temperature profile is an increasing function of the Eckert number; that is, the temperature profile increases with increasing Eckert number (v) e velocity profile is a decreasing function of the magnetic field M; that is, the velocity profile decreases with an increase in the magnetic field M [33][34][35] Nomenclature S: Unsteady parameter C p : Specific heat μ nf : Dynamic viscosity of the nanofluid axial directions W w : Stretching velocity ρ nf : Density of the nanofluids η: Similarity variable Ec: Eckert number φ: Solid particle volume fraction Q: Heat generation/absorption parameter σ * : Stefan-Boltzmann constant ] nf : Kinematic viscosity of the nanofluids C f : Skin friction T: Local temperature δ: ickness of the liquid film Pr: Prandtl number U w : Suction, injection velocity T δ : Temperature at the free surface k * : Mean absorption coefficient k * : Mean absorption coefficient q r : ermal radiative heat flux β: Nondimensional thickness of the film liquid k nf : ermal conductivity of the nanoparticles Nu: Nusselt number β nf : ermal expansion coefficient of solid particles (u, w): Velocity components.

Data Availability
No data were used to support this study.

Conflicts of Interest
e authors declare that they have no conflicts of interest for the publication of this paper.