Unsteady MHD Mixed Convection Flow in Hybrid Nanofluid at Three-Dimensional Stagnation Point

There has been significant interest in exploring a stagnation point flow due to its numerous potential uses in engineering applications such as cooling of nuclear reactors. Hence, this study proposed a numerical analysis on the unsteady magnetohydrodynamic (MHD) mixed convection at three-dimensional stagnation point flow in Al2O3–Cu/H2O hybrid nanofluid over a permeable sheet. The ordinary differential equations are accomplished by simplifying the governing partial differential equations through suitable similarity transformation. The numerical computation is established by the MATLAB system software using the bvp4c technique. The bvp4c procedure is excellent in providing more than one solution once sufficient predictions are visible. The influence of certain functioning parameters is inspected, and notable results exposed that the rate of heat transfer is exaggerated along with the skin friction coefficient while the suction/injection and magnetic parameters are intensified. The results also signified that the rise in the volume fraction of the nanoparticle and the decline of the unsteadiness parameter demonstrates a downward attribution towards the heat transfer performance and skin friction coefficient. Conclusively, the observations are confirmed to have multiple solutions, which eventually contribute to an investigation of the analysis of the solution stability, thereby justifying the viability of the first solution.


Introduction
The observation of magnetohydrodynamic (MHD) flow behavior is essential in diverse fields of engineering and has attracted considerable attention due to its significance in industrial applications, for instance in fossil-fueled power generation [1]. The existence of MHD in a fluid that conducts electricity gives rise to a resistive type force, which causes the fluid particle's motion resistance knows as Lorentz force. The Lorentz force intensifies the fluid temperature and concentration significantly, thereby slowing down the boundary layer's separation. The analysis throughout unsteady MHD flow at the forward stagnation point was first triggered by Katagiri [2]. Pavlov [3] developed the study of an electrically conducting MHD fluid in the boundary layer, including a transverse magnetic field related to a stretching sheet. The research was expanded by Takhar and Gupta [4] who testified the stabilizing effect detection on Taylor-Görtler three-dimensional disturbances and reviewed the solution stability in the magnetic field. Ever since a substantial number of studies with the consideration of MHD have been performed including [5][6][7][8]. Such research, however, has exempted the mixed convection flow.
Convective heat transfer, also known as convection, is a phenomenon where the flow of fluids transfers heat from one position to the next. The process under which an external source induces fluid motion is called forced convection. Free or normal convection, alternatively, is a process in which buoyancy forces alone produce fluid motion, arising from density differences. When both natural and forced convection systems function together, mixed convection mechanism occurs. The topic of mixed convection flow has drawn much interest from the researchers because of its prominence in manufacturing industries, for example, solar and nuclear collectors, heat exchangers, and atmospheric boundary layer flows [9]. The groundbreaking work on numerical analysis of mixed convection stagnation point flow past heated vertical flat surfaces was performed by Ramachandran et al. [10]. They extended the work done by Merkin [11] who found a non-uniqueness solution within a specific range of mixed convection parameter. The research was then broadened by Merkin [12] in his next exploration, where the stability of the results has been identified. Since then, a significant number of publications on mixed convection with the presence of MHD have been produced. Oztop et al. [13] performed MHD mixed convection laminar flow in a lid-driven cavity, Daniel and Daniel [14] investigated the MHD mixed convection flow with thermal radiation effect towards a stretching porous surface by utilizing the homotopy analysis method while Jamaludin et al. [15] examined the influence of the heat source/sink in MHD mixed convection stagnation point in a hybrid nanofluid. It is found that the heat transfer of the conventional alumina/water nanofluid is greater than then copper-alumina water nanofluid with the increment in the heat source/sink parameter.
To provide better-evolved heat conductivity, an innovative type of nanofluid introduced as hybrid nanofluid is invented. This alternative form of working fluid has now captivated many scientists due to its popularity in thermal properties advancement [16][17][18]. Specifying a good nanoparticles mixture is part of the main components in maintaining a robust nanofluid hybrid suspension. Xian et al. [19] and Gupta et al. [20] has studied the hybrid nanoparticle preparation method and the stabilization mechanism as well as its importance in the industrial sectors. Suresh et al. [21] conducted an exploratory practice to analyze the thermophysical characteristics of Al 2 O 3 -Cu/H 2 O hybrid nanofluid. Their discovery reveals that the stability of hybrid nanofluid depends mostly on the concentration of nanoparticles and the stability of higher-concentration nanofluid is incompetent. Consequently, the experimental results disclose Al 2 O 3 -Cu/H 2 O is capable of increasing thermal conductivity efficiency and reliability. In another report, Suresh et al. [22] used different concentrations of nanocomposite powder in Al 2 O 3 -Cu/H 2 O to conduct the synthesis, characterization and stability. The stability of prepared nanofluid is observed to decrease when the volume concentration is increased. Since hybrid nanofluid is claimed to have a range of beneficial characteristics in improving thermal conductivity, many researchers performed an investigation on the hybrid nanofluid by considering diverse aspect and conditions [23][24][25][26][27].
Researchers have widely studied the stagnation point flow focused on its uses in engineering areas and sectors such as wire drawing, paper making, hot rolling, and several others. The stagnation point flow field and heat transfer can be used to verify the consistency of certain goods. The flow of stagnation points was initially proposed by Hiemenz [28] in 1911. To tackle the two-dimensional stagnation point flow, he employed the similarity variables and achieved the exact solution to the problems. Afterwards, Libby [29] performed a boundary layer analysis on heat and mass transfer in three-dimensional stagnation point flow. Chiam [30] revealed the stagnation point flow analysis against a stretching sheet and expanded his study on heat transport analysis using the regular perturbation technique in stagnation point flow. Earlier stages of mixed convection flow in stagnation point flow have been studied by Takhar [31] in an incompressible fluid. Chamkha [32] studied continuous, two-dimensional, MHD mixed convection flows near a stagnation point of an electrically conducting and heat-absorbing fluid on a semi-infinite vertical permeable surface with arbitrary variations of surface heat flux. Meanwhile, Abdelkhalek [33] performed a numerical analysis to examine the impact of mass transfer in MHD mixed convection in view of stagnation point flow over a heated vertical permeable sheet. Jamaludin et al. [34] scrutinized the stability analysis of the mixed convection phenomenon toward a nonlinearly permeable surface. Researchers reviewed some recent inspections based on stagnation point flow through [35][36][37][38][39].
The overarching focus of this study is to conduct an analysis on unsteady MHD mixed convection stagnation point in the three-dimensional flow of alumina-copper/water hybrid nanofluid with stability analysis. According to the works described above, this propose problem remains briefly addressed in the literature. The idea of the present article was motivated by Noor et al. [40] and the evaluation was carried out by applying the bvp4c feature to obtain non-uniqueness solutions in the opposing flow past a permeable surface. The combination of alumina and copper nanoparticle were chosen in this study based on the outstanding works of Suresh et al. [21,22] as discussed earlier. The correlations properties of hybrid nanofluid are employed inspired by Takabi and Salehi [41] and Ghalambaz et al. [42]. The present work applied the bvp4c tool in the MATLAB platforms to address the constructed problem. More than one solution has been productively recognized by the stated approach method. In addition, analysis of solution stability is conducted to verify the solutions constancy for a valid physical explanation. This major participation will lead to stimulating industrial progress, particularly in the engineering and manufacturing sectors.

Mathematical Modeling
The current work considers unsteady MHD mixed convection near the stagnation point in the three-dimensional flow of Al 2 O 3 -Cu/H 2 O hybrid nanofluid. In this problem, u, v and w are the velocity component through x−, y− and z− axes with the origin at the nodal stagnation point N, as displayed in Figure 1. T wx (x, t) and T wy (y,t) are the variable temperatures, and T ∞ is the surrounding fluid temperature, where T wx (x, t), T wy (y, t) < T ∞ denote the opposing flow and T wx (x, t), T wy (y, t) > T ∞ apply to the assisting flow. We assume that where T 0 is the characteristic temperature of the surface of the sheet, with T 0 > 0 for assisting flow and T 0 < 0 denotes the opposing flow, while L is the characteristic length of the surface of the sheet. Here, δ is a parameter showing the unsteadiness of the problem towards t, time. Generally, δ = 0 signifies the steady inviscid flow, δ > 0 accelerates the outer potential flow, while δ < 0 corresponds to the reverse flow. The outer flow is assumed to be u e (x) = ax/(1 − δt) and v e (x) = by/(1 − δt) in x− and y− axes, respectively. Noticeably, c = b/a, where c represents the three-dimensional stagnation points parameter or ratio of velocity gradients at the edge of the boundary layer, with a and b are the principal curvatures parameter at N or the velocity gradients at the edge of the boundary layer along the x− and y− axes, respectively [43,44]. Both a, b are positive constants. Further, b = a corresponds to the axisymmetric case, while b = 0 is the plane stagnation flow problem. If a, b are positive, the solution of the corresponding equations results in nodal points of attachment, that is, 0 ≤ c ≤ 1. On the contrary, the saddle points of attachment are declared if a and b are negative, that is, −1 ≤ c ≤ 0. Also, it should be noted that when c = 0, the problem will convert to a two-dimensional case, while the axisymmetric case can be regained as c = 1. The plane of the body remains stagnant and there is a mass flux velocity denoted by w 0 as the surface is permeable, with w 0 > 0 referring to injection and w 0 < 0 representing the suction condition. The transverse magnetic field B 2 (x) = B 0 2 /(1 − δt) is considered normal, where B 0 is the magnetic field strength. From the above assumptions, the governing boundary layer equations can be defined as [9,43]: ∂T ∂t together with: together with:  Table 1 provides the correlation properties of Al2O3-Cu established by [41,42] while Table 2 presents the thermophysical properties [4 working fluid.   Table 1 provides the correlation properties of Al 2 O 3 -Cu/H 2 O as established by [41,42] while Table 2 presents the thermophysical properties [45] of the working fluid.

Properties Al 2 O 3 -Cu/H 2 O
Density Thermal capacity Dynamic viscosity Electrical conductivity Substituting Equation (6) into Equations (2)-(5), a series of similarity differential equations may be interpreted as: 1 Pr considering that: Here ε = δ/a refers to the unsteadiness parameter with ε = 0 signifies steady-state flow, ε < 0 implies the decelerating flow and ε > 0 denotes an accelerating flow and Pr stands for Prandtl number. The solution is obtained numerically and not using perturbation. The effect of ε is very well and in detail presented in Section 4: Discussion and Results. The mixed convection parameter symbolizes by Ω, where Ω > 0 signifies the assisting flow and Ω < 0 suggests an opposing flow, S represents the steady mass flux parameter (S > 0 for suction and S < 0 for injection) and M is the magnetic coefficient, which are described as where is the local Reynolds number. Subsequently, the physical quantities of interest are: Note that C f x , C f y is the skin friction coefficient along x− and y− axes, respectively, and Nu x is the local Nusselt number. The surface heat flux is identified as q w , whereas τ wx , τ wy are the shear stresses illustrated by: By exerting Equations (6) and (13) into Equation (12), we earn: provided that Re x = ax 2 /(1 − δt)ν f and Re y = ay 2 /(1 − δt)ν f .

Stability Analysis
A stability analysis is essential to verify the reliability of the obtained solutions since there exist more than one solution in the problem Equations (7)- (10). Following the contributions of [11,47], we introduce a dimensionless time variable τ, associated with the initial value problem. Now, a new conversion of similarity is proposed in accordance with the unsteady-state query as follows: Employing Equation (15) into Equations (7)-(9), the subsequent equations are guaranteed: 1 Pr subject to: As expressed by [48], to investigate the steady flow consistency f (η) = f 0 (η), h(η) = h 0 (η) and θ(η) = θ 0 (η), we write: where ω is the undetermined parameter of eigenvalue, as F(η), H(η) and I(η) are comparatively small to f 0 (η), h 0 (η) and θ 0 (η). The eigenvalue problem in Equations (16)- (18) leads to an infinite set of eigenvalues ω 1 < ω 2 < ω 3 . . . that trace a steady flow movement and primary deterioration while ω 1 is positive. However, when ω 1 is negative, the initial development of delays is observed, which exposes the erratic flow. Replacing Equation (20) into Equations (16)- (19), we have: 1 Pr and the boundary conditions are as follows: The steady-state flow solutions f 0 (η) and θ 0 (η) were implemented via τ → 0 . Subsequently, the corresponding linearized eigenvalue problem is defined: 1 Pr together with: By relaxing a boundary condition, the potential eigenvalues could be estimated [49]. Now, we assume F (η) → 0 , consequently, the eigenvalue problems in Equations (25)- (27) are discovered when F (0) = 1, where ω 1 is fixed.

Discussion and Results
The nonlinear ordinary differential equations presented in Equations (7)-(10) were solved by using the bvp4c feature in the MATLAB program [50]. The problem of the expressed boundary value is simplified to an ordinary first-order differential equation system, initially. The bvp4c feature is a noteworthy approach commonly exercised by numerous researchers to explain the justification of the reference value. A preliminary forecast of the variations step size and primary mesh point is required beneficial to the necessary response confirmation. In order to find more than one solution, the accurate estimation of boundary layer thickness, together with an early intervention guess is important.
The comparisons of results in Tables 3 and 4 are presented to validate the numerical procedure of the current study with the steady (ε = 0) numerical results from Noor et al. [40] and Eswara and Nath [43] for a different type of fluid. The previous studies by Noor et al. [40] and Eswara and Nath [43] tackled the viscous fluid problem, while the present study implemented the hybrid nanofluid. Further, Eswara and Nath [43] used an implicit finitedifference scheme with a quasilinearization technique while this study applied the bvp4c procedure in the MATLAB programming. It is observed that the present results are in good agreement with the solutions obtained for the steady regular fluid case; thus, this gives us confidence that the computational structure to analyze the hybrid nanofluid flow behaviors and heat transfer in this study can be employed with significant assurance. The key component for determining nanofluid flow behaviors and the efficiency of heat transfer is the organization of compatible single/hybrid nanofluids. Suresh et al. [22] performed the synthesis of Al 2 O 3 -Cu/H 2 O nanocomposite powder and its characteristics for various volume concentrations. In their noteworthy study, nanofluid stability is observed to decrease as volume concentration increases. Since the Al 2 O 3 -Cu/H 2 O hybrid nanofluid is supported in this current work, a different set of φ values fraction is limited in between 0.005 ≤ φ ≤ 0.02, corresponds to the work of [51]. The Al 2 O 3 -Cu/H 2 O hybrid nanofluid is selected in this study because of the outstanding work of Suresh et al. [21] in developing an exploration practice to scrutinize the Al 2 O 3 -Cu/H 2 O thermophysical properties. By diffusing the alumina (Al 2 O 3 ) nanoparticle followed by copper (Cu) into H 2 O, the Al 2 O 3 -Cu/H 2 O hybrid nanofluid is established with various sums of volume fractions [51,52].  Table 3. Approximation values of f (0) and h (0) by certain values of c when φ 1 = φ 2 = ε = Ω = M = S = 0, and Pr = 0.7.

c Present Result Noor et al. [40] Eswara and Nath [43]
f " (0) h " (0) f " (0) h " (0) f " (0) h " (0) that the scope of mixed convection parameters for which the solution occurs decreases with the inclusion of nanoparticle volume fraction in the first solution while the sheet is shrinking. The same findings have been found in the preceding literature, as stated by Waini et al. [24]. The buoyant force behaves in the same direction as the fluid movement in the assisting flow. This suggests the velocity of fluid flow to enhance thus capable of supporting the buoyant force, thereby, improves shear stress of the permeable surface. In contrast, the consequence of the buoyant force's opposing flow may affect the fluid velocity to become weak, hence affecting the fluid flow to slow down and minimizing the surface shear stress. Further, Figure 2 stresses that when  varied is denoted in Figure 3 and this phenomenon is in contrast with the second solutions. Concisely, the rate of heat transfer decreases as the nanoparticles volume fraction improves in a hybrid nanofluid. This observation suggests that the addition of nanoparticles volume fractions in the boundary layer may minimize the thickness of the thermal boundary layer thus improving the heat flux. Hence, prior to this subsequent case, we may infer that the incorporation of the nanoparticle volume fraction leads to the acceleration of the boundary layer separation.  Figure 5. Realize that the suction impact permits the molecules of Al2O3-Cu/H2O hybrid nanofluid to occupy the surface and then physically improve the heat transfer rate at the permeable sheet.      Figure 8, where dual velocity profiles are noted. As exemplified in Figure 8, the first solution decreases in proportion to the deteriorating of  , whereas the second approach revealed contradictory outcomes. In the meantime, the diverse progress of the solution in Figure 8 reflected as well the temperature profile    with the presence of unsteadiness parameter, which can be seen in Figure 9. Overall, both profiles asymptotically fulfilled the far-field boundary conditions (10) when         toward  , as portrayed in Figure 10. The magnetic field emergence across the electrically conducting fluid contributes to the Lorentz force appearance which prompts endurance to the motion of the fluid particle and therefore increases the fluid velocity (see Figure 12). Furthermore, an escalation of M leads to an intensification of the heat transfer rate as promoted in Figure 11. In short, improves in conjunction with the heat transfer rate along the permeable surface. More nanoparticles are drawn to the surface by the Lorentz force ensuing in greater temperature near the permeable sheet. In addition, by expanding the magnetic effect of the working fluid system, the thickness of the boundary layer is increased, thereby reducing the convection mechanism dramatically over the permeable wall surface, as shown clearly in the first solution of Figure 13.       Table 5. The flow represents an erratic flow when 1  appears negative because an initial extension of interruptions is proposed. The smallest eigenvalue, 1  clarifies the solution stability property to fix the authorizing disturbances, the flow is therefore steady ( 1  remains positive). It also suggests an early deterioration in the appearance of disruptions.

Conclusions
A numerical assessment of the unsteady MHD mixed convection stagnation point in Al2O3-Cu/H2O hybrid nanofluid at three-dimensional flow was established in the present work. The engagement of the bvp4c features in the MATLAB programming platform is employed to perform the numerical computation. The result of different regulating parameters, for example, the suction/injection parameter, the nanoparticle volume fraction, the unsteadiness and magnetic parameter were examined. Our analyses suggest that the occurrence of dual solutions is demonstrable for a wide variety of operating parameters, besides the stability analysis permits the first solution reliability. The  Figure 2 expresses the trend of f (0) and h (0) toward Ω when φ 2 varied. The findings obtained indicate that the scope of mixed convection parameters for which the solution occurs decreases with the inclusion of nanoparticle volume fraction in the first solution while the sheet is shrinking. The same findings have been found in the preceding literature, as stated by Waini et al. [24]. The buoyant force behaves in the same direction as the fluid movement in the assisting flow. This suggests the velocity of fluid flow to enhance thus capable of supporting the buoyant force, thereby, improves shear stress of the permeable surface. In contrast, the consequence of the buoyant force's opposing flow may affect the fluid velocity to become weak, hence affecting the fluid flow to slow down and minimizing the surface shear stress. Further, Figure 2 stresses that when Ω = 0 (static surface), f (0), h (0) = 1 which explains the lack of frictional drag on the sheet. A diminishing behavior over the first solution of the heat transfer performance or −θ (0) when φ 2 varied is denoted in Figure 3 and this phenomenon is in contrast with the second solutions. Concisely, the rate of heat transfer decreases as the nanoparticles volume fraction improves in a hybrid nanofluid. This observation suggests that the addition of nanoparticles volume fractions in the boundary layer may minimize the thickness of the thermal boundary layer thus improving the heat flux. Hence, prior to this subsequent case, we may infer that the incorporation of the nanoparticle volume fraction leads to the acceleration of the boundary layer separation. Figures 4 and 5 expose the effects of various value in S toward Ω past a permeable sheet. The characteristics of f (0) and h (0) in Al 2 O 3 -Cu/H 2 O is described in Figure 4. Figure 4 proves an improvement in S will decisively upsurge f (0) and h (0) in the first solution. In reality, the suction event may facilitate the boundary layer steadiness. In addition, the suction diminishes the friction of the external flow on the bodies, thereby reducing the boundary layer thickness and magnifying the velocity differential of the permeable sheet by removing the fluid across the low momentum surface. Both solutions convey an increase in −θ (0) as S escalates across the permeable sheet, as clarified in Figure 5. Realize that the suction impact permits the molecules of Al 2 O 3 -Cu/H 2 O hybrid nanofluid to occupy the surface and then physically improve the heat transfer rate at the permeable sheet.
The influence of the unsteady parameter ε toward Ω when ε shifts from −0.2 to −0.6 are demonstrated in Figures 6-9. The Al 2 O 3 -Cu/H 2 O hybrid nanofluid characteristic is demonstrated in Figure 6 with respect to the skin friction coefficient f (0), h (0) when ε varied in the unsteady case. Figure 6 captures that as ε reduced, the first solution has decreased in f (0), h (0) and the second solution has demonstrated a reverse effect. The reduction in ε leads to the expansion of the boundary layer thickness and subsequently declines the velocity gradient of the permeable sheet, hence f (0), h (0) diminished. The presence of nanoparticle volume fraction could also initiate the reduction of f (0), h (0) because of the rise in Al 2 O 3 -Cu/H 2 O hybrid nanofluid viscosity in the permeable surface. Moreover, according to the generated results in Figure 7, −θ (0) is decreased in the first solution which is proportional to the rate of heat transfer, when Ω < 0 (opposing flow) in the permeable sheet as ε reduces. In contrast, the second solution demonstrated an upward trend of −θ (0) the values of ε declines. From the current and existing evidence, the authors can infer that the unsteadiness parameter promotes significantly to the degradation of heat transfer. Even so, if multiple control parameters are taken into account, the authors would also like to claim that those effects may vary. The dimensionless profiles of velocity f (η), h (η) with various ε are accessible in Figure 8, where dual velocity profiles are noted. As exemplified in Figure 8, the first solution decreases in proportion to the deteriorating of ε, whereas the second approach revealed contradictory outcomes. In the meantime, the diverse progress of the solution in Figure 8 reflected as well the temperature profile θ(η) with the presence of unsteadiness parameter, which can be seen in Figure 9. Overall, both profiles asymptotically fulfilled the far-field boundary conditions (10) when η ∞ = 4 is implemented. Figures 10 and 11 illustrate the magnetic properties impact toward f (0), h (0) besides −θ (0). Apparently, we recognize that the first solution of f (0), h (0) has improved when the values of M rises in the Al 2 O 3 -Cu/H 2 O hybrid nanofluid flow toward Ω, as portrayed in Figure 10. The magnetic field emergence across the electrically conducting fluid contributes to the Lorentz force appearance which prompts endurance to the motion of the fluid particle and therefore increases the fluid velocity (see Figure 12). Furthermore, an escalation of M leads to an intensification of the heat transfer rate as promoted in Figure 11. In short, −θ (0) improves in conjunction with the heat transfer rate along the permeable surface. More nanoparticles are drawn to the surface by the Lorentz force ensuing in greater temperature near the permeable sheet. In addition, by expanding the magnetic effect of the working fluid system, the thickness of the boundary layer is increased, thereby reducing the convection mechanism dramatically over the permeable wall surface, as shown clearly in the first solution of Figure 13.
A stability analysis was further carried out by employing the bvp4c application in the MATLAB systems software. The smallest eigenvalues, ω 1 for certain values of Ω when φ 1 = 0.01, φ 2 = 0.02, S = 2.2, M = 0.02, c = 0.5, ε = −0.2 are listed in Table 5. The flow represents an erratic flow when ω 1 appears negative because an initial extension of interruptions is proposed. The smallest eigenvalue, ω 1 clarifies the solution stability property to fix the authorizing disturbances, the flow is therefore steady (ω 1 remains positive). It also suggests an early deterioration in the appearance of disruptions.

Conclusions
A numerical assessment of the unsteady MHD mixed convection stagnation point in Al 2 O 3 -Cu/H 2 O hybrid nanofluid at three-dimensional flow was established in the present work. The engagement of the bvp4c features in the MATLAB programming platform is employed to perform the numerical computation. The result of different regulating parameters, for example, the suction/injection parameter, the nanoparticle volume fraction, the unsteadiness and magnetic parameter were examined. Our analyses suggest that the occurrence of dual solutions is demonstrable for a wide variety of operating parameters, besides the stability analysis permits the first solution reliability. The augmentation in nanoparticle volume concentration surprisingly reduced the coefficient of skin friction and local Nusselt number. Thus, this leads to the conclusion that as the concentration of nanoparticles expands, the heat transfer rate decreases in Al 2 O 3 -Cu/H 2 O hybrid nanofluid for this particular problem. Meanwhile, an upsurge in the suction parameter intensity is capable of boosting the skin friction coefficient and the heat transfer rate of the Al 2 O 3 -Cu/H 2 O hybrid nanofluid. Consequently, a decrement in the unsteadiness parameter decreases the coefficient of skin friction with opposing flow over the permeable surface. On the contrary, it is also reported that improved magnetic control in Al 2 O 3 -Cu/H 2 O hybrid nanofluid escalates the rate of heat transfer. The magnetic fields escalation continues to interrupt the fluid development of the current study. Finally, the stability analysis is executed since the dual solutions are perceived to exist. The first solution's consistency and steadiness were verified by stability analysis, while the second solution is unconvincing and unstable.