Analytic Simulation for Magnetohydrodynamic Unsteady Buongiorno Model Hybrid Nanofluid Flow over Stretching

In this work, we inspect and analyze a two dimensional, unsteady mixed convectional hybrid nanofluid hydromagnetic flow (Al2O3-Cu/H2O) over a convectional heated an extending/contracting surface with the influence of thermal radiation. Hybrid nanofluid (Al2O3-Cu/H2O) flows with magnetohydrodynamic and heat source or sink. Brownian motion and thermophoresis were incorporated using the Buongiorno model. Hybrid nanofluid with vol. fraction range limited to 1.5% and within the higher temperature range of 50°C to 70°C is considered for thermal conductivity and viscosity analysis. The proposed model is then converted into ODEs through similarity transformation with the help of homotopy analysis. The effect of embedded input factors on the temperature, velocity, and concentration profiles is visually demonstrated and explained. The magnetic field has inverse impact on velocity and temperature profiles. Velocity profile increases for both mixed convection and buoyancy ratio parameters. It has been noticed that the temperature profile increases with thermal radiation. For increasing values of Lewis number, the concentration of hybrid nanoparticles is considerably lowered. Moreover, we observed an increase in the concentration of hybrid nanoparticles through a destructive chemical reaction, whereas a generative chemical reaction has the reverse effect. It has been proved that skin friction is increasing function of ε,Mand decreasing function of λ1,Nr . On the other hand, Nusselt number increased with the increase of R,Q,Nb,Nt while Sherwood number is decreased, with the increase of Nb,Nt , Le.


Introduction
It is eminent that many studies have been carried out on two cases of nanoparticles floating in a base fluid known as "Hybrid Nanofluid" the forefront nanofluid. The key importance of the hybrid nanofluid is that by selecting the right mix of nanoparticles, favorable characteristics may be enhanced, and drawbacks can be mitigated owing to their interactive influence. It is reported that these hybrid nanofluids are new and have a number of application acoustics, defense, manufacturing, transportation, medical, microflui-dics, and naval constructions. Nanofluid flow, in particular, is long familiar for its high heat conveyance when equated to regular fluid. The hybrid nanofluid is applied to raise it even further. Jamshed and Aziz [1] studied the Cattaneo-Christov heat flux impact and discovered that spherical shaped nanoparticles transmit heat at the fastest rate when compared to hexagon, platelet type nanoparticles. The authors in [2] examined squeezing flow in a hybrid base fluid with nanoparticles suspended in it. The hybrid nanofluid stagnation point flow past an extending sheet was investigated by the researchers [3]. The influence on flow between two riga-plates was investigated by Ahmed et al. [4]. The authors in [5] studied the impact of changing viscosity on a stretched sheet containing hybrid nanoparticles using the Runge-Kutta fourth order technique. The researchers in [6] interrogated the synthesis of hybrid nanofluid while in [7] examined the effects of a hybrid nanofluid of (Al 2 O 3 -Cu/H 2 O) on heat transmission.
There are numerous numerical methods in the literature to inspect the boundary layer flow and heat transport of a hybrid nanofluid. These concepts were further explored by the researchers in their research [9] to conceptualize the main idea. The problem was then expanded to a three-dimensional flow capable to the Newtonian heating consideration by the authors [10]. In both the above studies, they obtained key results about the heat transfer. The difficulty of a hybrid nanofluid consisting of a rotating flow was described by Hayat and Nadeem [11]. Zainal et al. [12] inspected the unsteady stagnation point flow of a hybrid nanofluid past a convectional heated an extending/contracting surface while accounting for velocity slip's influence on heat transfer. Using the Buongiorno model, Daniel et al. [13] investigated an unstable blended convectional electrical magnetohydrodynamic (MHD) flow and heat conveyance generated by nanofluid across a permeable stretched surface. The friction coefficient and wear volume were studied by Xie et al. [14] to determine the tribological properties of hybrid nanoparticles. Devi and Anjali [10] investigated the three-dimensional flow of (Cu-Al 2 O 3 /water) copper alumina/water hybrid nanofluid using the RK-Fehlberg integration technique.
Thermal radiation's effect on MHD blood flow and heat conveyance in a permeable capillary in extending motion was inspected. It has been examined that at high temperatures, heat transfer and thermal radiation are recognized to have a significant influence on numerous physiological processes, technological, and engineering industry equipment. Solar collector performance, plume dynamics, rocket propulsion, high-dose cancer therapy, fire propagation, material processing, and combustion systems are just a few examples. With the advancement of analytical and computational techniques, thermal convection flows with high radiative flux have received more attention recently. The effect of heat radiation on peristaltic transport of ionic nanoliquids in biomicrofluid channels had analyzed by Prakash et al. [15]. In context from melting heat transfer and thermal radiation, the researcher [16] proposed the dynamics of stagnation point flow of carbon nanotubes. To get series solution, the optimal homotopy technique was used. Muhammad et al. [17] make another notable attempt to address the key aspects of heat radiation and viscous dissipation impacts in a viscous nanofluid. Mixed convectional and slide impacts in viscous nanofluid across an extended sheet were described by Hsiao [18]. The solution was numerically defined using the convective boundary conditions [19]. The reader can further study about heat transfer through nanofluid flow in refs. [20][21][22][23][24][25][26][27]. The relevant study has been seen in [28][29][30][31].
The main aim of this research is to see how magnetic fields, thermal radiation, and heat generation/absorption affect unstable hybrid nanofluid hydromagnetic flow (Al 2 O 3 -Cu/H 2 O) across a convectional heated an extend-ing/contracting surface. The thermophoresis and Brownian motion characteristics were included in the Buongiorno model. To minimize the autonomous variables in governance equations resulting from mathematical modelling, a suitable collection of dimensionless variables is employed. By implementing the homotopy approach, an analytical solution has been calculated. Comparison between HAM and ND solve has been shown in Table 1.

Mathematical Formulation
Assume the unsteady hybrid nanofluid hydromagnetic flow (Al 2 O 3 -Cu/H 2 O) past a convectional heated an extending/ contracting surface in the bearing of the magnetic field, chemical reaction, heat generation/absorption, and thermal radiation. The flow problem plot is depicted in Figure 1, where u w ðx, tÞ = bx/ð1 − ctÞ, the extending/contracting is surface velocity, b indicates a constant that agree to extending ðb > 0 Þ and contractingðb < 0Þ instances, and c denotes the problem of unsteadiness. The velocity of the free stream is notified by u e ðx, tÞ = ax/ð1 − ctÞ, where a > 0 is the strength of the stagnation flow. T 1 and T 0 are used to represent the ambient and reference temperature separately. Then, we consider that the bottom of the surface is heated by convectional from a hot fluid at a particular temperature T f ðx, tÞ = T 1 − T 0 ðax 2 /2υÞ ð1 − ctÞ −3/2 and coefficient of heat transfer has noted by h f while the mass transfer coefficient is indicated by h s . The governing boundary layer equations may be recognized as [12,13] based on all of the assumptions mentioned.
∂T ∂t Applying concept of Rosseland approximation q r is as follows: 2 Advances in Mathematical Physics where the parameters σ 1 and k * are explained in nomenclature Further simplification leads us to the following: Here, we rewrite Equation (3) as follows: where u denotes the factor of velocity in x-axis, v is the velocity factor in y-axis, μ hnf is the where H 1 = Hð1 − ctÞ 1/2 is the velocity slip term, in which H stands for the starting merit of the velocity slip term. Table 2 lists the thermophysical characteristics of copper (Cu), as well as aluminium oxide (Al 2 O 3 ) and water (H 2 O) nanoparticles. The thermophysical characteristics of hybrid nanofluid are shown in Table 3. The solid volume fraction of nanoparticles is described as ϕ, where ρ f denotes the density, ρ s is the density of the hybrid nanoparticle, c p is the constant pressure of heat capacity, k f represents the thermal conductivity of H 2 O, and k s is the thermal conductivity of the hybrid nanoparticle.

Advances in Mathematical Physics
where ψ is the stream function that can be specified as u = ∂ψ/∂y and v = −∂ψ/∂x and η is the similarity variable. Thus, we attain Equations (2), (7), and (8) are converted into the following set of nonlinear similarity ODEs by using the similarity variables (10) and (11) in light of the aforementioned relationships.
Now the starting and final constraints (9) are also changed into the following: where ε, M, R, Pr, N r , Bi, Re x , λ, Q, λ 1 , Sc, N t , and N b are unsteadiness parameter, magnetic parameter, thermal radiation parameter, Prandtl number, buoyancy ratio parameter, Biot number, local Reynolds number in x-axis, ratio of velocity parameter, heat source/sink parameter, mixed convection parameter, Lewis number, thermophoretic parameter, and Brownian motion parameter.   Advances in Mathematical Physics Table  3: Applied models for thermophysical properties of the hybrid nanofluid [12].

Property
Hybrid nanofluid Thermal capacity Thermal conductivity where In dimensional form, we have from above as 2.2. Idea of HAM for the Model. Here, we used HAM to solve Equations (12)-(14) using boundary conditions (15). The Mathematica software is utilized for this, and the proposed model can be solved using HAM in the following manner: where operators are specified as L f _ and L θ _; moreover, we have We also define the operators, such that For Equations (12)- (14), the 0th-order system is shown as follows: with the BCs are While the embedding constraint is ζ ∈ ½0, 1, to regulate for the solution convergence ℏ f _, ℏ θ _, and ℏ Φ _ are used. At specific values ζ = 0 and ζ = 1, the following is obtained: ð26Þ While BCs are While

Outcomes with Discussion
In the following discussion section, the outcomes of various parametric quantity on velocity, concentration profiles, and temperature are shown graphically and discussed. These parameters are unsteadiness parameter, buoyancy ratio parameter, magnetic parameter, thermal radiation parameter, Lewis number, mixed convection parameter, Prandtl number and thermophorotic parameter, and Brownian motion parameter. The geometry of the problem is shown in Figure 1. It can be seen that when the magnetic parameter M assigns maximum values, the velocity distribution retains its declining tendency as shown in Figure 2. There is a periodic oscillation velocity with diminishing amplitude while the magnetic force produces a resistive nature force known as the Lorentz force which controls the flow of fluid particles. As a result, to demonstrate the value of reducing the velocity of moving particles, the effect of the unsteadiness parameter on ε the nanofluid velocity profile has seen in Figure 3. Because of the acceleration situation, ðε > 0Þ brings down rate of fluid flow and a narrower momentum boundary barrier thickness; the behavior occurs. The velocity profile reduces for the higher acceleration. Figure 4 illustrates the influence of the buoyancy ratio parameter on velocity drawings and also increasing levels of N r ; there is a noticeable increase in velocity. This increase in velocity drawings is more noticeable in stable condition than in the unsteadness situation. Figure 5 shows the effect of λ 1 on f ′ðηÞ. From Figure 5, it can be expected that f ′ðη Þ has greater values of λ 1 . Physically, this is because of the larger values of buoyancy force. The impact of heat generation/absorption on the temperature is plotted in Figure 6. The temperature of the fluid and the thickness of the thermal boundary level are both increased when a heat source ðQ > 0Þ is used. The heat sink ðQ < 0Þ, on the other hand, lowers the fluid's temperature and thins the thermal boundary level thickness.   To denote ðQ = 0Þ, the lack of heat generation or absorption, the increase in thermal radiation raises the temperature but lowers the concentration profiles, as seen in Figure 7. This is because an increase in thermal radiation gives to a greater extent heat to the hybrid nanofluid, resulting in a rise in temperature and the thickness of the thermal boundary layer. Physically, radiative aspect creates the Brownian motion of minute ingredients faster than normal; thus, random migrated particles strike with one another, and the caused frictional energy transforms it to thermal energy.
The effect from the Prandtl number along the temperature profile has seen in Figure 8. By increasing the Prandtl number, the fluid's thermal conductivity decreases. This fact is due to inverse relation of Pr with thermal diffusivity, and it is well-known fact that the fluid with higher values of Pr has weaker thermal diffusion so that the temperature declines. Figures 9 and 10 demonstrate the consequence of the thermophoresis parametric quantity ðN t Þ along the thermal profile. When the thermophoresis parametric quantity is increased which enhance the behavior of the system, thermophoresis effect is associated with movement of nanoparticles from a hot wall to a cold wall, and since it is generated by temperature gradients, it creates a fast flow away from the moving plate; thus, more fluid is heated away from the surface, and this leads to an increase in the temperature within the thermal boundary layer. Figure 10 demonstrates the consequence of the Brownian diffusion parameter ðN b Þ on the thermal profile. Physically, the imperfect nature of the Brownian motion parametric quantity enables heating the physical setup. This heating causes nanoparticles to be transferred from the cooler stretched sheet area to the quiescent fluid zone. Because the particles migrate from a high up to a low concentration part, an increase in the Brownian motion parameter N b and ΦðηÞ leads decrease in the fluid concentration gradient, as seen in Figure 11. The Brownian motion is the movement of fluid particles from higher to lower concentration so increases in the thermophoresis parameter ðN t Þ produce a huge rise in the concentration dispersion of the fluid flow ΦðηÞ, which meet to zero at the boundary level, as shown in Figure 12. A minor 13 Advances in Mathematical Physics change in the thermophoresis parameter causes the fluid particles to move quickly, releasing surplus heat energy and causing a huge rise in the concentration dispersion. As a result, with an increase in the value of ðN t Þ, Figure 12 indicates a considerable rise in the concentration dispersion. Figure 13 also shows the influence of the Schmidt number on the concentration dispersion. It is noticed that for larger values of Sc, the fluid concentration decreases. This fact is quite similar with that of Prandtl number impact on temperature. Actually, Sc has inverse relation with mass diffusivity; therefore, larger values of Sc correspond thinner concentration boundary layer. Table 4 shows that C f is increased when the values of ε, M are increased. The C f is decreased, when the values of λ 1 , N r are increased. Table 5 shows that Nu x is increased when the values of R, Q, N b , N t are increased. The Nu x is decreased, when the values of Pr are increased. Table 6 shows that Sh is decreased, when the values of N b , N t , Le are increased. Table 1 shows the very excellent agreement of the HAM and ND solve solution by computer-based package Mathematica 11.1.0.

Conclusions
In this research work, the uniform mixed convectional with the combined consequences of thermal radiation and heat generation/absorption along the magnetohydrodynamic (MHD) flow of hybrid nanofluid via a convectional heated the stretching/shrinking surface is studied. Brownian motion and thermophoresis were incorporated using the Buongiorno model. We came at the following key conclusions based on our findings with various factors in this study: (i) Increasing the combined convection and buoyancy ratio parameters improved nanofluid flow properties (ii) It is well established that an increase in the unsteadiness and magnetic parameters result in a decrease in the velocity profile (iii) Consequently, a rise in the thermophorotic parameter, radiation parameter, or Brownian motion parameter means an increase in the temperature profile (iv) It is noticed that the increase of the Prandtl number decreases the temperature profile