Magnetohydrodynamic hybrid nanofluid flow with the effect of Darcy–Forchheimer theory and slip conditions over an exponential stretchable sheet

The characteristics of water-based hybrid nanofluid flow when passing over an exponentially stretchable sheet with velocity and thermal slip factors are presented. This article provides a concept for a hybridized fluid comprising copper and cobalt iron oxide nanoparticles (NPs) dispersed into a base fluid (water). In addition, physical observations of the heat absorption behavior, the Darcy effect, the thermal radiation, and viscous dissipation are also taken into account. Because of their strong thermophysical properties, copper and cobalt iron oxide NPs are used in a wide range of applications in the engineering and medical fields. To study the dynamics of these NPs, a system of partial differential equations (PDEs) has been generated that forms a highly nonlinear coupled model. The PDE system is converted into nondimensional ordinary differential equations (ODEs) with the aid of similarity replacements. The semi-analytical homotopy analysis method (HAM) is applied to the set of dimensionless ODEs obtained to find the solution. Two engineering parameters, the Nusselt number and the skin friction, are plotted versus various parameters of the hybridized fluid using bar charts. It was observed that the no-slip condition, the suction parameter, and the Darcy–Forchheimer medium enhanced the thermal profile of the hybridized fluid.


Introduction
To reduce both costs and energy use significantly, intensification of thermal transfer is essential. Currently, advancements in science and technology are encouraging orders for compact devices with extraordinary features, the highest possible performance, precise operating characteristics, and long lifetimes. As a result, engineers and scientists have worked extensively on heat transmission analysis to upgrade the existing industrial work. A new class of nanotube suspensions created by Choi 1 demonstrated the enormous improvement in the thermal transport characteristics of these solid suspensions, which were called nanofluids, when compared with ordinary fluids. Nanofluid technology has proven to be a helpful tool in the development of oils and lubricants with higher performance. Nanofluids play vital roles in various industrial applications, including fabrication of paints and coatings, paper printing, power generation, cancer therapy, ceramics manufacture, drug delivery, and food production. Punith et al. 2 investigated nanofluid (ferromagnetic) flow with regard to the impact of a magnetic dipole over a stretching sheet and calculated that the thermal gradient was improved by augmentation of the Brownian motion and thermophoresis, with the Stefan blowing condition causing the nanofluid to display high heat transfer. Bhatti and Michaelides 3 analyzed nanofluid flow over a Riga plate in terms of the Arrhenius activation energy and found that the fluid's concentration profile was improved by augmentation of the activation energy; in this regard, the Brownian motion increased the temperature profile. Ullah 4 found that a reduction in the Nusselt number occurred as a function of the Brownian motion and thermophoresis parameter with a decreasing temperature gradient; however, a reduction in the Sherwood number also occurred as a function of the Brownian motion and the thermophoresis parameter under a similarly decreasing temperature gradient. Alsallami et al. 5 clearly explained the thermal transfer gradient and nanofluid flow behavior above a porous plate in terms of the impact of slip conditions; they claimed that addition of nanoscale particles with high thermal conductivity and a large volume fraction represents the best way to increase heat flux in nanofluids, but their simulation results contradicted the results reported previously by other researchers. [6][7][8][9][10][11][12][13][14] With regard to the development of these technologies, a new class of fluids for heat transfer has been introduced in the literature and has led to technological updates worldwide; these fluids are known as hybrid nanofluids. A nanofluid is the combination of a base fluid with a single type of added nanoparticle, while the hybrid nanofluid is an extension of the nanofluid that was invented by dispersion of dual or multiple types of nanoscale particles with sizes of less than 100 nm and excellent thermal conductivity when immersed into the base fluid. Algehyne et al. 15 investigated the attributes of a mixture of hybrid nanofluids over a heated plate exchanger and determined that the mixture provides the maximum coefficient of heat transfer. Al 2 O 3 + H 2 O has a higher pressure drop than the other hybrid nanofluid, with a value that is slightly higher than that of water. Ullah et al. 16 described an examination of the flow of a hybridized fluid across a curvy surface, and found that when injection and suction were compared, the hybrid nanofluid achieved a higher heat transfer rate than the simple nanofluid. Hussain et al. 17 investigated a computational model for hybridized fluid dynamics under convective conditions on a rotating sheet. Their results showed that as the stretching ratio increased, the velocity profiles of the CuO and TiO 2 hybrid base fluids decreased, while their velocity curves increased. Rashid et al. 18 provided insight into the flow behavior of hybrid nanofluids with a physical perception of mixed convection over a vertical surface fixed by a porous medium, revealing that there are dual solutions for opposing flows. Zhao et al. 19 studied the improvement in the thermal transmission of a hybridized fluid containing Ag + CuO + H 2 O NPs; their research showed that the chemical reaction, thermal gradient, and heat generation characteristics of the hybrid nanofluid were higher than those of nanofluids based on a single nanoparticle type, even in the presence of radiation. Tlili et al. 20 proposed a model of hybrid nanofluid (methanol) 3D flow based on magnetohydrodynamics (MHD) over a surface with uneven thickness and slip factors. Their results determined that the thermal transfer rate of the hybridized fluid was significantly greater than that of the corresponding nanofluid. The Lorentz force effect was also smaller on the hybridized fluid when compared with the effect on the nanofluid. Usman et al. 21 proposed a model of hybrid nanofluid MHD thin film flow over an unsteady spinning disk that considered the impact of couple stress, and the percentage improvement caused by heat transfer was calculated and shown graphically. Their investigation established that the hybridized fluid has a higher heat transmission capability when compared with that of conventional nanofluids. Ullah 22 also investigated nanofluid and hybridized fluid flow with entropy generation in thermal structures. Many more articles are available in the literature that specify that the hybrid nanofluid represents an important alternative solution to traditional thermal systems. [23][24][25][26][27][28][29][30][31][32] A stretchable sheet with the fluid motion characteristics described above in fluid mechanics terms is highly significant for both biomedical sciences and industrial processes. The ability to obtain a similar solution for fluid flow passing over a stretching sheet was first studied by Crane. 33 However, when compared with the stretching scenario itself, the flow passing over such a stretchy sheet received little attention. Alharbi et al. 34 stated that the flow over a shrinking sheet is a backward flow. With regard to this statement, Rashid et al. 35 explored the research on viscous fluid flow over a nonlinearly stretching sheet and found that the flow features revealed many interesting behavioral aspects that confirmed further study of the effects of nonlinear stretching on the flow characteristics is warranted. Alhowaity et al. 36 investigated the homotopy perturbation method and axisymmetric flow over a stretching sheet, and comparison of their results with the exact solutions showed that the homotopy perturbation method gives outstanding results. The analysis of Bilal et al. 37 of temperature field inflow over a stretching sheet with uniform heat flux showed that the temperature at a point declines with an upsurge in the Prandtl number. Zhao et al. 38 explained the effects of slip on MHD viscous flow over a stretching sheet and obtained an exact solution to show that there is only one physical solution for any combination of the slip, the magnetic parameters, and the mass transfer parameters. Chu et al. 39 explored the concept of a thermal gradient rate in viscoelastic fluid MHD flow over a stretchy sheet under the combined impact of variable thermal conductivity and a nonuniform heat source. Their results showed that in the boundary layer region, the combined effects of the radiation, the variable thermal conductivity, and the nonuniform heat source had a noteworthy effect on the control of the thermal gradient rate. The flow of fluid over a stretching sheet with respect to the effect of a magnetic field was studied using a power law model by Zhao et al. 40 and they concluded that in this case, a thinner boundary layer occurred because of the effect of the magnetic field, thus causing increased wall friction.
As a result of the wide range of its potential applications, hybridized fluid motion through permeable porous surfaces has now become the center of attention for many researchers. Possible application areas of flow through permeable porous surfaces include the environmental field, the solar sciences field, nuclear engineering, and materials science. To understand this type of porous surface flow, Darcy's law is erroneously but widely used. To use the permeability to calculate the effect of inertia, Forchheimer 41 devised a second-order polynomial for use in the momentum equation. Muskat and Wyckoff 42 undertook a study to identify the Forchheimer component. Darcy-Forchheimer concepts are used in many studies of the regularity of hybrid nanofluids, and several researchers have studied porous medium flow. A number of these studies are cited here. Ullah et al. 43 investigated the Darcy-Forchheimer hybrid nanofluid flow through a porous medium with variable features. Ullah et al. 44 analyzed the flow of a hybridized fluid under the influence of a Marangoni convective pattern through a Darcy-Forchheimer medium. Sun et al. 45 performed a study of a Darcy-Forchheimer hybrid nanofluid with MoS 2 and SiO 2 nanoparticles with extension of the entropy generation. Bilal et al. 46 proposed a model of heat and mass transfer through a Darcy-Forchheimer hybrid nanofluid flow over a stretching curved surface. Li et al. 47 investigated a dusty hybrid nanofluid with two phases via extension of the viscous dissipation over a cylinder. Nazeer et al. 48 proposed a model of the hybridized fluid dynamics due to solar radiation in a Darcy-Forchheimer permeable medium over a flat plate. Inspired by the approaches described above, this work presents the results of development of a model of a sophisticated form of hybridized fluid that contains nanoparticles of Cu (copper) and CoFe 2 O 4 (cobalt ferrite) flowing over an exponentially stretching sheet. The core of the proposed model is the analysis of the effects of different factors on the flow and heat transport behavior of the hybrid nanofluid. In addition, the effects of various factors on the dynamics of the hybridized fluid regimes are sketched both graphically and statistically.
Novelty is added to the recent work by extending the published work of Ahmed and Akbar 49 by addressing the following points: The recent work is based on a hybridized fluid composed of nanoparticles of Cu (copper) and

Mathematical Formulation
The electrically conducting and incompressible steadystate dynamics of a 2D Darcy-Forchheimer flow of H 2 O-based hybridized fluid past an exponentially stretching sheet are considered. The effects of partial slip and the thermal jump conditions on the dynamics are investigated. The nanoscale particles of copper (Cu) and cobalt ferrite (CoFe 2 O 4 ) are immersed in the base fluid of water (H 2 O). To develop a mathematical model, a Cartesian coordinate system is proposed in which the stretchable sheet surface is aligned with the x-axis and the y-axis represents the vertical direction, with the flow field being constrained to the area where y.0. The geometry of the problem is depicted in Figure 1. In establishing the governing equations for the flow of the hybridized fluid, the following assumptions are made: Because of the ambient velocity e 1 x ð Þ, the hybrid nanofluid flow permeates over the stretchy surface. We ignore a weak induced magnetic field generated by application of a magnetic field to the flow at an angle u. The nanoparticles (Cu and CoFe 2 O 4 ) immersed in the base fluid (H 2 O) are supposed to be in the state of thermal equilibrium, and therefore noslip phenomena occur between them. Combining the concepts of optically thick radiation, heat absorption, and viscous dissipation enhances the heat transfer rate. T l ðxÞ represents the temperature of the hybridized fluid at the stretchable sheet surface, while T ' represents the ambient temperature of the hybridized fluid. Because no externally applied electric field is present, the polarization effect can be overlooked.
The constitutive momentum flow equation and the energy equation for the hybrid nanofluid are then determined as follows using the assumptions above: 52 In these equations, the velocities of the hybrid nanofluid are represented by u and v along the x and y directions in the Cartesian coordinates, respectively, the temperature of the hybrid nanofluid is represented by T , and the exponential temperature of the hybrid nanofluid passing over the stretchy sheet is T l (x) = T ' + (T 0 À T ' )e ax=2l . The stretchable sheet velocity is denoted by e 1 (x) = e 0 e x=l . In the stretchable sheet, a, T 0 , and e 0 represent the reference temperature distribution parameter, temperature, and velocity, Using the Rosseland approximation, the radiation heat flux q r can be written as: The Stefan-Boltzmann constant and the heat absorption constant in equation (4) are denoted by s Ã and k Ã , respectively. In addition, a Taylor series expansion about T ' (the free stream temperature) is executed for linearization of the term T 4 while ignoring the higher degree factors, as follows: Implementation of equations (4) and (5) in equation (3) gives: The initial and boundary conditions are given by: at y = 0: The experimental values of the parameters for water, Cu, and CoFe 2 O 4 are reported in Table 1. In Table 2 The transformation variables are: By incorporating equation (8) into equations (2) and (3), we obtain: The transform conditions are: Here, M is the magnetic parameter, x is the dimensionless coordinate, Tb is the radiation parameter, Pr is the Prandtl number, Re is the Reynolds number, Ec is the Eckert number, and L represents the velocity slip. Ha represents the heat absorption, S is the injection/suction, d is the thermal slip parameter, Fr is the Forchheimer parameter, j denotes the local porosity parameter, and Gr is the Grashof number.
ð Þ n f m f Re 2 : ð12Þ Terms for the local Nusselt number Nu x and the skinfriction coefficient SF x are calculated in both dimensional and dimensionless forms. The dimensional forms of Nu x and SF x are given by:

HAM Solution
The HAM has been used to determine the solution for the proposed model. 57 For this purpose, the linear operators are given as follows: The extended forms of L and L Y are given as: By using Taylor's expansion, we obtain: Next, we obtain The equations above can also be written as: where the local Reynolds number is given by

Analysis of Results
Fluid flow phenomena over a stretching surface under the influence of various physical effects, including heat sources, a porous surface, and application of an inclined magnetic field, have been studied. The flow scenario has been modeled in the form of nonlinear partial differential equations, which have been handled analytically using the HAM. The results are presented graphically in the figures, while the relevant numerical values are tabulated in the tables.  Electrical Conductivity

Discussion of Results
This section evaluates the physical mechanism behind each figure in the velocity and temperature profiles. The following observations have been made.
Velocity and Energy Profiles. Figure 2(a) and (b) shows the relationships of the nondimensional magnetic parameter M with the velocity field F 0 h ð Þ and the temperature field Y h ð Þ, respectively, for both slip (L = 0:5) and no-slip (L = 0) conditions. Note that intensifying the value of M causes F 0 h ð Þ to decrease, while the value of Y h ð Þ is enhanced. Physically, a gradual improvement in M creates a resistance force to the flow of the hybrid nanofluid called the Lorentz force, through which the frictional drag force is incremented. Figure 3(a) and (b) are drawn to represent the effect of the suction parameter S.0 ð Þ being greater than zero on the velocity profile F 0 h ð Þ and the temperature profile Y h ð Þ, respectively, for the cases of the slip (L = 0:5) and no-slip (L = 0) conditions. When S.0, the values of F 0 h ð Þ and Y h ð Þ are reduced, but the corresponding profiles in Figure 3(c) and (d) show the opposite behavior when the injection parameter S\0 ð Þ is less than zero. Physically, the momentum boundary layer is fixed to the stretchy sheet surface in the scenario when suction occurs, and this breaks the flow momentum and reduces the profiles of both F 0 h ð Þ and Y h ð Þ for the hybrid nanofluid. Injection, in contrast, improves the fluid flow by causing a lateral mass flux to occur across the stretchy sheet, which adds momentum to the fluid flow. As a result, the hybrid nanofluid's F 0 h ð Þ and Y h ð Þ profiles improve in this case. Figure 4(a) and (b) illustrate the impact of the thermal radiation Tb and heat absorption Ha, respectively, on the temperature profile Y h ð Þ of the hybrid nanofluid in the cases of slip (L = 0:5) and no-slip (L = 0) conditions. Note that the value of Y h ð Þ is enhanced because of the increased value of Tb. The reason for this enhancement is that Tb is subject to the Y h ð Þ of the surroundings and is released in an electromagnetic wave form. Therefore, Tb is assumed to be a function of Y h ð Þ. The conduction phenomena of the hybrid nanofluid are improved by the thermal effect, and as a result, the boundary layer becomes thicker. Therefore, the Y h ð Þ of the hybridized fluid is noted to have an increasing trend. The Y h ð Þ is reduced by augmentation of Ha. Physically, enhancing the value of Ha increases the heat absorption capacity of the hybridized fluid significantly, and the Y h ð Þ profile thus vanishes. Figure 5(a) and (b) illustrate the relationships of the angle u to the velocity profile F 0 h ð Þ and the temperature profile Y h ð Þ of the hybrid nanofluid in the cases of slip (L = 0:5) and no-slip (L = 0) conditions. Increasing the value of u causes the F 0 h ð Þ to vanish while enhancing the Y h ð Þ. Physically, the magnetic field strength increases with increasing u, which generates a greater resistive force to F 0 h ð Þ, causing the fraction force to increase and Y h ð Þ to improve.     the cases of slip (L = 0:5) and no-slip (L = 0) conditions for the hybridized fluid. F 0 h ð Þ declines when l is increased because j is inversely related to the Darcian drag force; therefore, the hybrid nanofluid flow permeability improves, causing F 0 h ð Þ to vanish. Figure  8(b) show the relationship of j versus Y h ð Þ for the hybrid nanofluid in the cases of slip and no-slip conditions. In the presence of porous media, the resistance to hybrid nanofluid flow is enhanced, which induces higher temperatures and increased thermal layer thickness. Figure 9(a) illustrates the effect of the Grashof number Gr on the velocity profile F 0 h ð Þ in the cases of slip (L = 0:5) and no-slip (L = 0) conditions for the hybrid nanofluid, which shows that F 0 h ð Þ decreases as Gr increases. Figure 9(b) shows that an increase in Gr causes the Y h ð Þ to decline. This occurs because augmentation of Gr causes a convection current enhancement, thus improving Y h ð Þ while also reducing the hot surface area of the sheet. Figures 1 to 9 show that the temperature profile Y h ð Þ has a greater enhancement effect under the no-slip (L = 0:5) condition as compared with that under the slip condition (L = 0), while in the case of the velocity profile F 0 h ð Þ, the opposite phenomena can be observed. Physically, in the slip condition case, there is a fraction produced between the sheet surface and the fluid that enhances the temperature profile and reduces the velocity profile.
Skin Friction and Nusselt Number. Bar Chart 1(a) to (d) show the effect of increasing the value of the skin friction SF x versus the suction parameter S, the injection parameter ÀS, the magnetic parameter M, and the angle u in the cases of the slip and no-slip conditions. Similarly, Bar Chart 2(a) to (h) illustrate the relationships of the Nusselt number Nu x to the Eckert number Ec, the slip parameter d, the thermal radiation Tb, the angle u, the suction parameter S, the injection parameter ÀS, the magnetic parameter M, and the absorption parameter Ha. Note that increased values of Tb, S, and Ha enhance Nu x , while increased values of Ec, d, u, ÀS, and M reduce Nu x .

Conclusion
The analysis presented here reports the dynamic characteristics of a water-based hybridized fluid composed of copper (Cu) and cobalt ferrite (CoFe 2 O 4 ) NPs when passing over an exponentially stretchable sheet. The previously published work is extended by using the concepts of thermal and velocity slip factors, the heat absorption coefficient, the magnetic field effect, Darcy-Forchheimer parameters, and a hybridized nanomaterial. The core findings from these observations are presented as follows: The fluid motion F 0 h ð Þ of a hybrid nanofluid is slowed substantially by elevating the magnetic field angle and strength, but this elevation has a reverse effect on the energy profile Y h ð Þ: Raising the injection factor causes the fluid dynamics profile F 0 h ð Þ and energy profile Y h ð Þ to be enhanced, while the reverse effect is produced by raising the suction parameter and the Darcy-Forchheimer parameters. The thermal profile Y h ð Þ of the hybrid nanofluid can be increased by increasing the thermal slip factor d. The contributions of the heat absorption Ha and the thermal radiation Tb boost the heat transmission rate. The dispersion of the copper and cobalt ferrite NPs in the base fluid magnifies the velocity and the energy transmission rate remarkably because of the fluid's remarkable thermophysical properties.
The transformation of mechanical energy into the thermal energy of the nanoparticles occurs when the Eckert number Ec increases. As a result, the temperature field Y h ð Þ is enhanced by increasing Ec. In the case of the slip condition L = 0:5 ð Þ , the thermal profile curve is higher than in the case of the no slip L = 0 ð Þ condition, because a fraction is created between the sheet surface and the fluid as a result of the slip condition that helps to enhance the thermal profile Y h ð Þ.

Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.