Mixed Convection Hybrid Nanofluid Flow Induced by an Inclined Cylinder with Lorentz Forces

Hybrid nanofluids may exhibit higher thermal conductivity, chemical stability, mechanical resistance and physical strength compared to regular nanofluids. Our aim in this study is to investigate the flow of a water-based alumina-copper hybrid nanofluid in an inclined cylinder with the impact of buoyancy force and a magnetic field. The governing partial differential equations (PDEs) are transformed into a set of similarity ordinary differential equations (ODEs) using a dimensionless set of variables, and then solved numerically using the bvp4c package from MATLAB software. Two solutions exist for both buoyancy opposing (λ < 0) and assisting (λ > 0) flows, whereas a unique solution is found when the buoyancy force is absent (λ = 0). In addition, the impacts of the dimensionless parameters, such as curvature parameter, volume fraction of nanoparticles, inclination angle, mixed convention parameter, and magnetic parameter are analyzed. The results of this study compare well with previously published results. Compared to pure base fluid and regular nanofluid, hybrid nanofluid reduces drag and transfers heat more efficiently.


Introduction
A nanofluid is a solid-liquid mixture containing nanoparticles and a base liquid. According to Daungthongsuk and Wongwises [1], the thermal conductivity of a fluid is boosted by mixing nanoparticles with a regular fluid. Because thermal conductivity is the most important issue in heat transfer problems, many mathematical models for nanofluids have been developed by researchers such as Buongiorno [2] and Tiwari and Das [3]. In addition, issues surrounding the improvement of heat transfer are also important for other cases such as employing microchannel heat sink (MCHS) with V-ribs and electrohydrodynamic (EHD) using water-based nanofluid. The heat transfer performance of a MCHS is enhanced due to the periodically arranged V-ribs that interrupt the thermal boundary, which produce chaotic convention and augment the heat transfer area [4]. Wang et al. [5] mentioned that the major processes involved in increasing heat transfer in microchannels are chaotic convection caused by the electrophoretic effect and thermophoretic effect. Although nanofluids have extremely efficient qualities and properties for industrial and engineering operations, researchers are still searching for fluid variants by studying several of their physical aspects such as the type, shape, and volume fraction of nanoparticles. Hybrid nanofluid exhibits significant developments in its thermal and rheological properties when compared to regular nanofluid, particularly in improving the heat conductivity of base fluids. The upgraded features of nanofluids result from the addition of two or more types of nanoparticles to a common fluid that is called a hybrid nanofluid and serves the purpose of increasing the heat transfer rate. Suresh et al. [6] mentioned that researchers have used a variety of nanoparticles including carbon nanotubes, metallic particles (Cu, Al, Fe, Au, and Ag), and non-metallic particles (Al 2 O 3 , CuO, Fe 3 O 4 , TiO 2 , and SiC). Furthermore, if the nanoparticle materials are properly chosen, they can positively enhance each other and significantly improve the thermal conductivity of the nanofluid. Based on the characteristics of metallic and non-metallic nanoparticles, it can be expected that the addition of metallic nanoparticles such as copper (Cu) into a nanofluid composed on a basis of Al 2 O 3 nanoparticles can enhance the thermophysical properties of this mixture. It has been determined that most hybrid nanofluids studied have a higher thermal conductivity than nanofluid, but Guan et al. [7] further explored the reason why hybrid nanofluids have high thermal conductivity. The nanolayer densities and diffusion coefficients were calculated for various hybrid nanofluids to explain the underlying mechanism of the enhancement in the thermal properties. The results showed that a hybrid nanofluid with Cu-Ag 50%/Ar can obtain the greatest improvement in thermal conductivity compared to liquid Ar, Cu/Ar and Ag/Ar nanofluids.
Waini et al. [8] reported that a hybrid nanofluid, also known as a base fluid with multiple types of nanoparticles, had improved thermal properties. Highly effective thermal conductivity, enhanced heat transfer, the synergistic effect of nanoparticles, stability, and good aspect ratio are the advantages of a hybrid nanofluid. Lower operating costs and higher energy efficiency are the best definitions of enhanced thermal conductivity of a hybrid nanofluid [9]. A hybrid nanofluid mixed convection flow has applications in various modern manufacturing processes, such as metal extrusion, wire drawing, welding, and heat pipe, etc. [10]. Suresh et al. [11] investigated how an Al 2 O 3 -Cu/H 2 O hybrid nanofluid affected heat transfer. They identified that the hybrid nanofluid had a greater friction coefficient than that of an Al 2 O 3 /H 2 O nanofluid. Alshuhail et al. [12] reviewed the thermal efficiency of hybrid nanofluids in solar energy applications. Ranga et al. [13] reported that for Cu-CuO/H 2 O hybrid nanofluid at 5% volume fraction, the collector efficiency increased by 2.175%. An experimental analysis conducted by Farajzadeh et al. [14] found that the collector efficiency increased by 26% on using Al 2 O 3 -TiO 2 /H 2 O hybrid nanofluids.
The mathematical analysis of the dual solutions of the boundary layer flow along moving surfaces has a practical application in the engineering scenario. It enables the determination of the most practical, stable, and physically acceptable solutions. Early research on Blasius flow for non-unique solutions by considering various flow situations and boundary conditions was conducted by Steinheuer [15]. He performed an investigation on boundary layer flow over a semi-infinite long flat plate with moving but impermeable soil to give a non-unique solution. Later, Klemp and Acrivos [16] reported that non-unique solutions exist in the case where the plate and the free stream move with the reverse flow. The flow field differs from the ordinary boundary layer in that an inviscid collision area exists around the site of detachment, where the reverse boundary-layer flow is turned in the direction of the mainstream. In addition, other scholars such as Bognár [17], Hussaini and Lakin [18], and Ishak and Bachok [19] reported similar problems with the moving flat plate for Newtonian fluid flow. It was discovered that the existence of solutions for a semi-infinite plate is dependent on the ratio of the plate surface velocity to the free stream velocity. It has been demonstrated that a solution exists only if this parameter does not exceed a specific critical value, and numerical computations have been performed to demonstrate that this solution is non-unique. However, the existence of non-unique solutions also occurs for the Blasius boundary conditions problem when the surface is at rest, as reported by Khashi'ie et al. [20], Waini et al. [21], Sohut et al. [22], etc.
A stagnation point (SP) is a point on a solid body which directly faces the stream and where the streamlines are separated. On the other hand, the flow of a fluid in the neighborhood of an SP is called stagnation point flow (SPF). The exploration of SPF is a hot topic in fluid mechanics because of its functions and applications in industries. The stagnation point flow, also commonly known as the Hiemenz flow, was established by Hiemenz [23] and can be accurately analyzed using the Navier-Stokes equation. He was the first researcher to find an outcome for the two-dimensional SPF model over a horizontal plate. Hiemenz [23] solved the Navier-Stokes equation numerically by transforming it into a regular differential equation using a suitable similarity transformation. Many cases that involved SPF have also been discovered by many scholars. For instance, Chiam [24] explored the SPF past a stretching sheet. Wang [25] examined both axisymmetric and twodimensional stagnation flows toward a shrinking sheet. He found that the flow structure is complicated due to the non-alignment of the stagnation flow on the shrinking sheet. A few years later, Lok et al. [26] continued the investigation by including the effect of MHD. However, due to the shrinkage effect, the streamlines are not always parallel, but a reverse flow is designed or originated near the surface. Furthermore, Awaludin et al. [27] investigated the SPF over a shrinking as well as a stretching sheet. They found that the solution for the stretching case is unique, but two solutions were found for the shrinking case.
Many industrial applications and natural phenomena involve mixed convection flows (MCFs), which involve a combination of free and forced convection flows, such as electronic equipment cooled by fans, flows in the ocean, underground cable systems, etc. In MCF, the existence of dual solutions has recently become a subject of discussion among scholars. For example, Ishak et al. [28] explored the magnetohydrodynamic MCF on a vertical continuous porous wall using a finite difference method. They reported that both opposing and assisting flows produce dual solutions. The existence of the solutions is dependent on the buoyancy parameter, where the range increases with suction effect at the boundary. Later, Aman et al. [29] continued the study by adding the slip effect at the boundary. Abbas et al. [30] studied the MHD mixed convection flow between concentric cylinders and reported that the mixed convection increases the velocity on the heated cylinder, but it decreases on the cold cylinder. The combined effects of MHD and slip velocity on the mixed convection flow around a circular cylinder was investigated by Ullah et al. [31]. They reported that an increase in slip factor enhances velocity as well as temperature. The MHD flow of a hybrid nanofluid past a vertical plate was considered by Zainal et al. [32], who found two branches of solutions in a certain range of the mixed convection strength. The temporal stability study verified that only one of the solutions was stable. Wang et al. [33] explored the influence of magnetic and electric fields on heat transfer to analyze the thermal conductivity and heat transfer enhancement of nanofluids, respectively, and the chaotic convention. They found that the applied magnetic and electric fields significantly affect the heat transfer performance of the nanofluid.
The study of fluid flow over an inclined cylindrical surface has significance for applications in industrial and engineering processes, e.g., magnetohydrodynamic (MHD) power generators, the polymer industry, and gas turbines [34]. Studies on the flow over a horizontal or a vertical cylinder have been explored widely; however, studies involving an inclined cylinder remain uncommon. An analysis of stratified stagnation flow induced by an inclined cylinder with mixed convection and a magnetic field was presented by Hayat et al. [35]. Bilal et al. [36] studied the Darcy Forchheimer flow of a hybrid nanofluid (consisting of carbon nanotubes) through the impermeable inclined cylinder using a homotopy analysis method. They found that the higher the inclination angle, the slower the fluid flow. Rehman et al. [37] extended the study to include logarithmic and parabolic curve fitting analysis and stated that analysis is the first step and can be of great help to previous studies. In addition, the study of MHD mixed convection on an inclined cylinder has also attracted the attention of researchers for the case of nanofluid flow, such as Dhanai et al. [38], who examined the effects of slip boundary conditions, Brownian motion, thermophoresis, and viscous dissipation over an inclined cylinder. They discovered that heat transfer increases with the augmentation of the mass transfer parameter but the opposite occurs for the thermal slip parameter. Moreover, Gupta and Sharma [39] studied a similar problem in the presence of thermal radiation. They used the differential transform method (DTM) with Padé approximation in their study and found that the DTM can overcome traditional perturbation limits, assumptions and restrictions.
Based on the previous investigations, hybrid nanofluid flow over an inclined cylinder has not been extensively studied. Binary hybrid nanofluid flow is a popular topic in fluid dynamics inspired by the above literature. The current investigation intends to discover and analyze the impact of magnetohydrodynamic on hybrid nanofluid flow over an inclined cylinder. This study is the extension of what has been investigated by researchers, such as Rehman et al. [37] and Devi and Devi [9], which considers binary hybrid Al 2 O 3 -Cu/H 2 O nanofluid. The stimulus of the governing parameters is generally represented via several graphs and tables, and the quantitative outcomes obtained are validated for the limited cases using the results of previous studies. Furthermore, the current study reveals that the stability investigation proves the physical outcome. These duality and stability results were also the main objectives of the current work.

Description of the Mathematical Model
Consider the buoyancy effects on a magnetohydrodynamic SPF of a steady incompressible hybrid nanofluid past an inclined cylindrical surface, as shown in Figure 1. The coordinates (x, r) are measured in the corresponding axial and radial directions of the cylinder. It is assumed that the temperature of the cylindrical surface, T w , is higher than the surrounding temperature, T ∞ . In the current study, the thermophysical properties of the hybrid nanofluid are taken to be uniform. Meanwhile, copper nanoparticles are combined with the base fluid to produce Cu/H 2 O nanofluid, which is then combined with the alumina nanoparticles, Al 2 O 3 , to create the required hybrid nanofluid, Al 2 O 3 -Cu/H 2 O. transform method (DTM) with Padé approximation in their study and found that the DTM can overcome traditional perturbation limits, assumptions and restrictions. Based on the previous investigations, hybrid nanofluid flow over an inclined cylinder has not been extensively studied. Binary hybrid nanofluid flow is a popular topic in fluid dynamics inspired by the above literature. The current investigation intends to discover and analyze the impact of magnetohydrodynamic on hybrid nanofluid flow over an inclined cylinder. This study is the extension of what has been investigated by researchers, such as Rehman et al. [37] and Devi and Devi [9], which considers binary hybrid Al2O3-Cu/H2O nanofluid. The stimulus of the governing parameters is generally represented via several graphs and tables, and the quantitative outcomes obtained are validated for the limited cases using the results of previous studies. Furthermore, the current study reveals that the stability investigation proves the physical outcome. These duality and stability results were also the main objectives of the current work.

Description of the Mathematical Model
Consider the buoyancy effects on a magnetohydrodynamic SPF of a steady incompressible hybrid nanofluid past an inclined cylindrical surface, as shown in Figure 1. The x r are measured in the corresponding axial and radial directions of the cylinder. It is assumed that the temperature of the cylindrical surface, w T , is higher than the surrounding temperature, T ∞ . In the current study, the thermophysical properties of the hybrid nanofluid are taken to be uniform. Meanwhile, copper nanoparticles are combined with the base fluid to produce Cu/H2O nanofluid, which is then combined with the alumina nanoparticles, Al2O3, to create the required hybrid nanofluid, Al2O3-Cu/H2O. The boundary layer equations comprising hybrid nanofluid can take place as follows [9,37]: The boundary layer equations comprising hybrid nanofluid can take place as follows [9,37]: subject to the boundary conditions Here, v and u are the elements of velocity along the r and x axes, respectively, B 0 is the uniform magnetic field, g is the gravity acceleration, ω is the inclination angle, T is the temperature of the hybrid nanofluid, and a is a positive constant. Where In Table 1, φ hn f = φ 1 + φ 2 is the hybrid nanoparticle volume fraction, where φ f corresponds to the regular liquid. In addition, φ 1 and φ 2 indicate aluminum oxide (Al 2 O 3 ) and copper (Cu) nanoparticle volume fractions, and their solid components are symbolized by the subscripts Al and C, respectively. In addition, µ, ρ, k, C p , σ, and β refer to the absolute viscosity, density, TCN, specific heat capacity, electrical conductivity and thermal expansion coefficient, respectively. The subscripts hnf and f, respectively, stand for the requisite posited hybrid nanofluid and the regular fluid. Table 2 provides the values of the thermophysical characteristics of the nanoparticles and water (regular fluid).  [40], Ganguly et al. [41] and Waini et al. [42]).

Properties Hybrid Nanofluid
Density In order to get similarity solutions, the dimensionless variables (6) are introduced [37,44]: where prime denotes differentiation with respect to η, and ψ is the stream function defined as u = r −1 ∂ψ/∂r and v = −r −1 ∂ψ/∂x, which identically satisfies the continuity Equation (1). The temperature at the surface of the cylinder is assumed to be T w = T ∞ + bx, where b is a constant. By utilizing the similarity variables (6), Equation (1) is identically satisfied, whereas Equations (2)-(4) are transformed into the following ODEs: where M stands for magnetic parameter, Pr is Prandtl number, K is the curvature parameter, and λ is the mixed convection parameter (constant). Further, the buoyancy term is equal to the ratio of the local Grashof number, Gr x , and the square of the local Reynolds number, Re x 2 . The dimensionless quantities are defined as: The skin friction coefficient and local Nusselt number (heat transfer rate) are given by: which, after applying (5) and (11), then become the following: where Re x = u e x/v f is the local Reynolds number.

Stability Analysis
Several researchers [45][46][47] have demonstrated the temporal stability of the numerical solutions. The purpose being to verify the stability of the solutions as time evolves. In the present study, two solutions for Equations (7)-(9) are found within a certain range of the physical parameters. It is worth investigating the time stability of the solutions and determining which one of the two solutions is reliable, and which one is not reliable, over time. For the working procedure, the following dimensionless time variable τ = at is introduced, as proposed by Merkin [48] and Weidman et al. [49]. The unsteady or time-dependent form of Equations (2) and (3) are considered, and then they are reduced to: ρC p hn f ∂T ∂t where t stands for time. Based on Equation (6), the dimensionless transformations are written as follows: Substituting Equation (15) into Equations (13) and (14) results in 1 Pr subjected to: To verify the temporal stability of the dual solutions, we follow Weidman et al. [49] by introducing the following perturbations: Here, the functions F and G are considered small compared to the steady solutions f 0 and θ 0 , and γ stands for an unknown eigenvalue. Substituting (19) into Equations (16)- (18) results in 1 Pr By rendering τ = 0, the stability of the steady flow solutions f 0 and θ 0 can be investigated, and hence the initial decay or growth of the solutions can be identified by F = F 0 (η) and G = G 0 (η), respectively. The linearized eigenvalue equations are: 1 Pr subject to the boundary conditions Following Harris et al. [50], the condition F 0 (η) → 0 as η → ∞ in Equation (25) is relaxed, and is substituted by the new boundary condition F 0 (0) = 1, without loss of generality. This replacement is performed in order to obtain the smallest eigenvalue γ 1 from the infinite set of eigenvalues γ 1 < γ 2 < γ 3 . . .. The flow solution is considered stable only when the generated smallest eigenvalue is positive (γ 1 > 0), and this positive smallest eigenvalue eventually approaches zero as λ → λ c , indicating that the solution is stable as the applied perturbation decreases with time.

Results and Discussion
The system of Equations (7)-(9) was numerically solved using the bvp4c package available from MATLAB R2021a software. It is necessary to set an initial guess at the initial mesh point, and to set the step size, ∆η, in order to obtain the specified accuracy. In the present study, we set the boundary layer thickness as η ∞ = 10 and the step size as ∆η = 0.01 with the relative error tolerance set to 10 −10 . The initial guesses are dependent on the parameter values applied for computation. The validity of the numerical results was checked by examining the related profiles. Moreover, the numerical solutions were compared with the available results for the particular case of this study.
Several representing values of parameters were selected in such a way as to examine the behaviors of the heat transport and flow dynamics characteristics. For instance, the present problem comprised several influential parameters such as curvature parameter K, mixed convection parameter λ, magnetic parameter M, inclination angle ω, and the solid volume fraction of nanoparticles φ 1 , φ 2 . For clarity, subscripts 1 and 2, which represent Al 2 O 3 and Cu, respectively, were introduced. In addition, the effects of these parameters on fluid motion, skin friction coefficient, temperature and heat transfer profiles are shown in several graphs and tables. For the purpose of the computations, we used the following default values: M = 0.1, ω = 30 • , Pr = 6.2, K = 0.2 and φ 1 = φ 2 = 0.01. Furthermore, the values of f (0) and −θ (0), which signify the skin friction coefficient and the heat transfer rate, were compared with those obtained previously by Lok et al. [51] and Ishak et al. [28], which show a good agreement. This validation is highlighted for the case of assisting flow, λ = 1 in the absence of hybrid nanoparticles φ 1 = φ 2 = 0, curvature parameter K = 0, inclination angle ω = 0, and magnetic parameter M = 0, for different values of Prandtl number Pr, as presented in Tables 3 and 4.           the skin friction increases, raising the velocity gradient and the temperature ent (0) θ ′ − , which are proportional to the skin friction and the heat transfer rate, r tively (see Equation (12)). In addition, in the presence of a magnetic field (M > 0), th transfer rate is enhanced as the inclination angle increases. This is because the ma force counteracts the gravitational buoyancy force at <90° but no longer counterac ω = 90°.       increases when the curvature parameter increases. This tendency is seen be of the larger impacts of the curvature parameter; the more it decreases the radius cylinder and then accelerates the fluid motions, the less contact there is with the cyl surface, which causes a reduction in resistance to fluid flow. However, the opposite  increase for the first solution. This is due to the decrement in the radius of the cylinder and the increment of the fluid motions. Less contact between the fluid and the cylinder surface causes a reduction in resistance to the fluid flow. Moreover, increasing K will delay the boundary layer separation. Using cylindrical surfaces in boundary layer problems is one approach to obtain a high skin friction and a better heat transfer rate. However, the second solution demonstrates the contrary, with f (0) values that reduce as K rises. The second solution presented in Figure 3 produces the opposite results for both flow regions. The heat transfer rate −θ (0) increases for the phenomenon of buoyancy assisting flow (λ > 0); however, it decreases for the phenomenon of buoyancy opposing flow (λ < 0) due to the superior impacts of K. Figures 4 and 5 demonstrate the influences of the hybrid nanoparticle volume fractions φ 1 , φ 2 on the skin friction and heat transfer rate, respectively, when M = 0.1, Pr = 6.2, ω = 30 0 and K = 0.2. It was discovered that the bifurcation value λ c for water (φ 1 = φ 2 = 0.0) was −4.6636; for Cu/H 2 O nanofluid (φ 1 = 0.0, φ 2 = 0.01) it was −4.9360; and for Al 2 O 3 -Cu/H 2 O hybrid nanofluid (φ 1 = φ 2 = 0.01), it was −5.0731. It was realized that the bifurcation values λ c (in the opposing region) and the values of f (0) increase when the hybrid nanoparticle volume fractions increase. The results prove that the augmentation of the hybrid nanoparticle volume fraction improves fluid viscosity and thus enhances skin friction, leading to an increase in drag forces. Figure 5 illustrates the same dual solutions pattern as Figure 3. As the nanoparticle volume fractions φ 1 , φ 2 increase, the heat transfer rate for the first solution also increases, which implies that the Al 2 O 3 -Cu/H 2 O hybrid nanofluid transfers heat more efficiently than the Cu/H 2 O nanofluid and pure base fluid. However, for the second solution, the heat transfer rate decreases for the opposing flow but rises for the assisting flow. It is worth noting that a hybrid nanofluid can have a higher heat transfer rate than a regular nanofluid as well as a pure fluid, as has been proved experimentally and numerically by numerous researchers, see for example [52,53].
Moreover, Figures 6 and 7 display the changes in the values of the skin friction coefficient and the heat transfer rate of the dual solutions for the numerous values of an inclination angle ω when Pr = 6.2, M = 0.1, φ 1 = φ 2 = 0.01 and K = 0.2, respectively. The value of f (0) for opposing flow (λ < 0) increases as the inclination angle rises. Physically, as the value of ω increases relative to the posited x-axis and when the curvature parameter is present, the effect of buoyancy force due to gravity decreases, causing f (0) to increase but −θ (0) to decline. The critical values λ c = −5.0731, −6.2111, and −8.7722 for ω = 30 • , 45 • , and 60 • , respectively, are shown in this illustration as well. The larger values of inclination angle augment the drag force, thus lowering the velocity within the boundary layer. The reason for this is the decreasing effect of gravity (see Rehman et al. [54]). Figure 7 illustrates how the increasing value of the inclination angle reduces the rate of heat transmission. As the temperature rises, the rate of heat transmission decreases. The influences of the drag force and heat transfer rate against λ for some values of magnetic parameter M when ω = 30 • , Pr = 6.2, K = 0.2 and φ 1 = φ 2 = 0.01 are demonstrated in Figures 8 and 9. Because of the Lorentz force, the velocity and temperature gradients escalate with rising values of the magnetic parameter. This force is generated by the presence of a magnetic field in the fluid flow, which forms resistance to the movement of the fluid. As a result, the skin friction increases, raising the velocity gradient f (0) and the temperature gradient −θ (0), which are proportional to the skin friction and the heat transfer rate, respectively (see Equation (12)). In addition, in the presence of a magnetic field (M > 0), the heat transfer rate is enhanced as the inclination angle increases. This is because the magnetic force counteracts the gravitational buoyancy force at <90 • but no longer counteracts it at ω = 90 • . Figures 10-19 reveal the variations in velocity and temperature profiles for assisting flow (λ > 0) with several values of non-dimensional parameters. Figures 10 and 11 demonstrate the influences of K on the fluid velocity as well as the fluid temperature. At first, f (η) increases when the curvature parameter increases. This tendency is seen because of the larger impacts of the curvature parameter; the more it decreases the radius of the cylinder and then accelerates the fluid motions, the less contact there is with the cylinder surface, which causes a reduction in resistance to fluid flow. However, the opposite effect is seen when η ∞ is approximately greater than one. In contrast, the temperature field shows the opposite behavior to the curves of velocity. The temperature initially decreases (η ∞ < 1) and then rises when η ∞ is roughly greater than one.
in Figure 8. Additionally, Figure 17 shows that the temperature distribution reduc the magnetic parameter rises. Deceleration in the thermal boundary layer thickness l to the increment of the heat flux. This shows a similar trend to that shown in Figur This finding agrees with the illustration in Figure 9. Figures 18 and 19 show the influences of the buoyancy force on the velocity and perature for numerous values of the mixed convection parameter λ when ω = 30°, Pr = K = 0.2, ϕ1 = ϕ2 = 0.01 and M = 0.1. In Figure 18, the velocity increases for the first solu whereas the second solution shows the opposite trend with the increment of λ (oppo to assisting flow). Nevertheless, the temperature for the first solution decreases when mixed convection parameter increases, as shown in Figure 19. The second solution sents inconsistent behavior in the temperature profile.                   Table 5 presents the smallest eigenvalues, γ, of the Al2O3-Cu/H2O hybrid nano for several values of K when M = 0.1, Pr = 6.2, ω = 30° and ϕ1 = ϕ2 = 0.01. The obta minimum eigenvalues are positive for the first solution, but negative for the second s tion. Referring to Equation (19), positive values of γ indicate a decay of disturban time evolves τ → ∞ . On the other hand, negative values of γ show an initial grow disturbance, and thus the solutions obtained for this case are unstable in the long Moreover, for both solutions, the smallest eigenvalue γ approaches zero as λ approa the bifurcation value λc. This means that, at turning points (from the first to the se Further, the effects of the inclination angle ω on the profiles are depicted in Figures 12 and 13. Physically, an increase in the inclination angle reduces the buoyancy effects, thus decreasing the velocity but enhancing the temperature in the boundary layers. The results in these figures correspond to the results shown in Figures 6 and 7. As can be seen in Figure 12, when the inclination angle rises, the velocity decreases (for both solutions). However, the second solution augments when the boundary layer thickness η ∞ is roughly greater than two. Meanwhile, Figure 13 shows a similar trend where the second solution displays the opposite behavior. We note that in the present study, there is a reduction to a vertical cylinder when ω = 0 • and to a horizontal cylinder, when ω = 90 • . The buoyancy effects for the latter case are very small and thus can be neglected. Figures 14 and 15 portray the influences of the hybrid nanoparticle volume fraction φ 1 , φ 2 for the base fluid (φ 1 = φ 2 = 0.0), Cu/H 2 O nanofluid (φ 1 = 0.0, φ 2 = 0.01) and Al 2 O 3 -Cu/H 2 O hybrid nanofluid (φ 1 = φ 2 = 0.01). It was found that the velocity elevates when the volume of the nanoparticles augments, as presented in Figure 14. The thickness of the momentum boundary layer reduces as the concentration of nanoparticles uplifts. Hence, the velocity of the fluid flow increases. Figure 15 shows that as the hybrid nanoparticle rises, the temperature curves and the thermal boundary layer thickness also upsurges. The surface heat flux lowers as the thermal boundary layer thickens. These discoveries are according with the results illustrated in Figures 4 and 5. However, the second solutions for both figures show contradictory behavior, which is associated with the assisting flow (λ = 1).
Apart from that, Figures 16 and 17 present the effects of the magnetic parameter M on the velocity and temperature profiles when ω = 30 • , Pr = 6.2, K = 0.2, φ 1 = φ 2 = 0.01 and λ = 1. By increasing the values of M, the velocity increases, as portrayed in Figure 16. The momentum boundary layer becomes thinner as the impact of M is boosted, causing an increase in the velocity gradient at the surface f (0). Thus, the skin friction coefficient increases with increasing M for the present case. This result is consistent with that shown in Figure 8. Additionally, Figure 17 shows that the temperature distribution reduces as the magnetic parameter rises. Deceleration in the thermal boundary layer thickness leads to the increment of the heat flux. This shows a similar trend to that shown in Figure 15. This finding agrees with the illustration in Figure 9. Figures 18 and 19 show the influences of the buoyancy force on the velocity and temperature for numerous values of the mixed convection parameter λ when ω = 30 • , Pr = 6.2, K = 0.2, φ 1 = φ 2 = 0.01 and M = 0.1. In Figure 18, the velocity increases for the first solution, whereas the second solution shows the opposite trend with the increment of λ (opposing to assisting flow). Nevertheless, the temperature for the first solution decreases when the mixed convection parameter increases, as shown in Figure 19. The second solution presents inconsistent behavior in the temperature profile. Table 5 presents the smallest eigenvalues, γ, of the Al 2 O 3 -Cu/H 2 O hybrid nanofluid for several values of K when M = 0.1, Pr = 6.2, ω = 30 • and φ 1 = φ 2 = 0.01. The obtained minimum eigenvalues are positive for the first solution, but negative for the second solution. Referring to Equation (19), positive values of γ indicate a decay of disturbance as time evolves τ → ∞ . On the other hand, negative values of γ show an initial growth in disturbance, and thus the solutions obtained for this case are unstable in the long run. Moreover, for both solutions, the smallest eigenvalue γ approaches zero as λ approaches the bifurcation value λ c . This means that, at turning points (from the first to the second solution), there is a change in the values of λ from stable (positive) to unstable (negative). Thus, this finding supports the stability and physical reliability of the first branch solutions. As shown by the temporal stability analysis, the second solution is unstable in the long run. Although such solutions may lack physical significance (in the long run), they are still of mathematical interest as far as the differential equations are concerned. The second solution is also a solution to the differential equations. Similar equations may arise in other situations where the corresponding solution can have a more realistic meaning.

Conclusions
The present study investigated the behavior of the mixed convection flow of an Al 2 O 3 -Cu/H 2 O hybrid nanofluid over an inclined cylindrical surface. The effects of the dimensionless parameters on the skin friction coefficient and the heat transfer rate at the surface were numerically computed and graphically presented. The summarized discoveries are listed as follows: • Two solutions were obtained for both buoyancy assisting and opposing flows, whereas a unique solution was found in the absence of buoyancy force (λ = 0); • Larger inclination angles ω lead to a lower gravitational force, increasing the friction factor f (0) but decreasing the heat transfer −θ (0); • The skin friction coefficient and the heat transfer rate of the Al 2 O 3 -Cu/H 2 O hybrid nanofluid (φ 1 = φ 2 = 0.01) is higher than that of the nanofluid (φ 1 = 0%, φ 2 = 1%) and water (φ 1 = φ 2 = 0.0); • Larger values of the curvature parameter K and mixed convection parameter λ result in a slower detachment of the boundary layer; • Velocity increases with increasing values of the buoyancy parameter λ, magnetic parameter M, and hybrid nanofluid nanoparticle volume fraction φ 1 , φ 2 , whereas it declines for curvature parameter K and the inclined angle ω; • The velocity and temperature gradients upsurge with a higher impact of the magnetic parameter; • According to the temporal stability analysis, the first solution is physically stable as time evolves, whereas the second solution is not reliable in the long run. Data Availability Statement: Not applicable.

Conflicts of Interest:
The authors declare no conflict of interest.