Unsteady Stagnation Point Flow of Hybrid Nanofluid Past a Convectively Heated Stretching/Shrinking Sheet with Velocity Slip

Unsteady stagnation point flow in hybrid nanofluid (Al2O3-Cu/H2O) past a convectively heated stretching/shrinking sheet is examined. Apart from the conventional surface of the no-slip condition, the velocity slip condition is considered in this study. By incorporating verified similarity transformations, the differential equations together with their partial derivatives are changed into ordinary differential equations. Throughout the MATLAB operating system, the simplified mathematical model is clarified by employing the bvp4c procedure. The above-proposed approach is capable of producing non-uniqueness solutions when adequate initial assumptions are provided. The findings revealed that the skin friction coefficient intensifies in conjunction with the local Nusselt number by adding up the nanoparticles volume fraction. The occurrence of velocity slip at the boundary reduces the coefficient of skin friction; however, an upward trend is exemplified in the rate of heat transfer. The results also signified that, unlike the parameter of velocity slip, the increment in the unsteady parameter conclusively increases the coefficient of skin friction, and an upsurge attribution in the heat transfer rate is observed resulting from the increment of Biot number. The findings are evidenced to have dual solutions, which inevitably contribute to stability analysis, hence validating the feasibility of the first solution.


Introduction
The most common problem in boundary layer flow that had been treated so far is much focused on those for steady flows. Even though it is the steady flows that seem to have the utmost significance in real-world demands, some cases of time-varying in the boundary layer which is unsteady indicate an important role in several engineering problems. Some of the examples are start-up processes where the motions in rest are transits from a steady flow to another, and periodic motions of the working fluid [1]. The behavior of unsteady boundary layer flow describes an unusual pattern compared to the steady flow owing to the additional time-dependent terms in the governing equations, which exaggerated the separation of boundary layer and the fluid motion arrangement [2,3]. The thermal and mechanical properties of such an unsteady mechanism in the boundary layer approximation have been studied both analytically and numerically. Elbashbeshy and Bazid [4] presented the numerical investigation towards the unsteady stretching surface with heat transfer analysis. At the same time, Bhattacharya [5] are more relatable to real-world applications. For instance, Seth et al. [39] studied the unsteady hydromagnetic nanofluid flow and heat transfer past a non-linearly stretching surface with Navier's velocity slip and presented the analysis of entropy generation. Other valuable references regarding unsteady nanofluid flow with the velocity slip effect can be found in [40][41][42].
Ever since the outstanding inventions achieved by Choi and Eastman [43], who originated the brilliant idea of demonstrating the nanoparticle suspension in a base fluid and came out with the nanofluid term, a better type of working fluid is still being pursued. Acknowledging the sufficient improvement in the thermal conductivity of the conventional fluid is crucial, an advanced nanofluid form known as hybrid nanofluid is introduced, which intends to have highly developed heat conductivity. This modern type of fluid agent has fascinated numerous researchers owing to its reputation in the emergence and improvement of thermal characteristics in realistic applications, including micro-channel, heat pipes, heat exchangers, air conditioning systems, and mini-channel heat sink [44,45]. Gupta et al. [46] and Xian et al. [47] have reviewed the preparation method of hybrid nanoparticles along with the stabilization and its significance in industrial sectors. One of the critical elements in establishing a sustainable hybrid nanofluid suspension is selecting an appropriate combination of nanoparticles. The most widely used nanoparticles for the formation of hybrid nanofluid suspension are carbon materials (graphite, MWCNTs, CNTs), metals (Cu, Ag), metal oxides (Al 2 O 3 , CuO, Fe 2 O 3 ), metal carbide, and a metal nitride. Madhesh and Kalaiselvam [48] conducted an experimental study to examine the features of hybrid nanofluid as a coolant agent, Tahat and Benim [49] examined the efficiency of hybrid nanofluid on flat plate solar collector, and they verified that the viscosity, thermal conductivity, and density of the working fluid had increased together with the concentration of Al 2 O 3 /CuO concentration, thus enhancing the solar collector proficiency. Some early research on hybrid nanofluid that employed the numerical method was done by Labib et al. [50], who investigated the impact of base fluids and hybrid nanofluid using a two-phase mixture model in forced convective heat transfer. Moghadassi et al. [51] revealed that the heat transfer performance is enhanced by adding the nanoparticles of Al 2 O 3 -Cu hybrid nanofluid while creating a small pressure drop in the system regime, Devi and Devi [52] focused on the mathematical inspection towards a stretching sheet. In contrast, the evaluation of heat transfers in the natural convection of Al 2 O 3 /water nanofluid and Al 2 O 3 -Cu/water hybrid nanofluid with a discrete heat source was explored by Takabi and Salehi [53]. Additional details on this topic are well described in the literatures [54][55][56].
To the best of the authors' knowledge, the existing literature does not consider the unsteady stagnation point flow of hybrid nanofluids with the presence of velocity slip parameter and stability analysis in their models. Thus, the addressed issues above have inspired the authors to perform a numerical study in unsteady stagnation point flow towards a convectively heated stretching/shrinking sheet in alumina-copper/water (Al 2 O 3 -Cu/H 2 O) with the impact of velocity slip on heat transfer. The hybrid nanofluid is recognized by dispersing Al 2 O 3 nanoparticles into H 2 O, followed by Cu with different volume fractions and the thermophysical properties of the hybrid nanofluid are adopted from Ghalambaz et al. [57] and Takabi and Salehi [53], which were based on the feasible physical assumptions and are in agreement with the conservation of mass and energy. The present work also utilized the bvp4c approach, which can be accessed in the MATLAB programming system towards solving the formulated problem. The existence of more than one solution is predictable; thus, an analysis of solution stability is completed to confirm the steadiness of the solutions which has an actual physical interpretation. The explanation of the results and the convergence of the obtained solutions are deliberated on in detail. Particular cases of current findings are evaluated in accordance with those of Mahapatra and Nandy [34] and Wang [58]. Furthermore, the consensus between previous and current findings is outstanding, and the agreement is excellent.

Mathematical Model
The unsteady two-dimensional stagnation-point flow of a hybrid Al 2 O 3 -Cu/H 2 O nanofluid over a convectively heated stretching/shrinking sheet with the influence of velocity slip is considered in this research work, as illustrated in Figure 1 (see Dzulkifli et al. [59]). The stretching/shrinking velocity is denoted by u w (x, t) = bx/(1 − ct), where b denotes a constant corresponds to stretching (b > 0) and shrinking (b < 0) cases while c signifies the unsteadiness problem and u e (x, t) = ax/(1 − ct) is the velocity of the free stream where a > 0 represents the strength of the stagnation flow. The ambient temperature and the reference temperature are T ∞ and T 0 , respectively. Now, we let the bottom of the sheet be heated by convection from a hot fluid at a specific temperature which supplies a coefficient of heat transfer, expressed by h f . From all of the assumptions above; the governing boundary layer equations can be acknowledged as [34].
where H 1 = H(1 − ct) 1/2 is the velocity slip factor, in which H refers to the initial value of the velocity slip factor. The copper (Cu) thermophysical properties, along with aluminum oxide (Al 2 O 3 ) and water (H 2 O) nanoparticles, are provided in Table 1, as demonstrated by [60]. In the meantime, Table 2 issued the thermophysical properties hybrid nanofluid as established by [53,57]. The nanoparticles solid volume fraction is represented by φ, ρ f indicates the H 2 O density, and ρ s is the density of the hybrid nanoparticle, C p is the constant pressure of heat capacity, while k f denotes the thermal conductivity of H 2 O and k s is the hybrid nanoparticles thermal conductivity.

Mathematical Model
The unsteady two-dimensional stagnation-point flow of a hybrid Al2O3-Cu/H2O nanofluid over a convectively heated stretching/shrinking sheet with the influence of velocity slip is considered in this research work, as illustrated in Figure 1 (see Dzulkifli et al. [59]). The stretching/shrinking velocity is denoted by     [57]).

Properties Hybrid Nanofluid
Dynamic viscosity Density Thermal capacity Thermal conductivity In order to express the governing Equations (1)-(3) concerning the boundary conditions (4) in a much simpler form, the subsequent similarity transformations are presented [34] where ψ is the stream function that can be specified as u = ∂ψ/∂y, v = −∂ψ/∂y and η is the similarity variable. Thus, we attain In view of the above relations, by employing the similarity variables (5) and (6), Equations (2) and (3) reduce to the following set of nonlinear similarity differential equations 1 Pr Here, ε measures the unsteadiness parameter with ε = c/a, Pr represents the Prandtl number where Pr = ν f /α f . Next, the initial and boundary conditions (4) now transform into From Equation (9), λ symbolises as the ratio of velocity parameter, γ and Bi are the dimensionless velocity slip parameter and Biot number, respectively, which are described as Next, we define the skin friction coefficient C f and the local Nusselt number (Nu x ) as The shear stress along the x− direction is represented by τ w , while q w signifies the surface heat flux that accentuated by By applying (5) and (12) into (11), we acquire is the local Reynolds number in xaxis.

Analysis of Solution Stability
By obeying the efforts of Merkin [61] and Merill et al. [62] from their outstanding discoveries of stability analysis scheme, the unsteady equations need to be deliberated in order to ultimately identify the reliable and stable solution since we notice the appearance of non-uniqueness solutions in the boundary value problem (7)- (9). Now, in accordance with the unsteady-state problem, a new similarity conversion is proposed Employing the similarity variables of Equation (14) to Equations (7) and (8), we now obtain the following converted differential equations 1 Pr with respect to In accordance with Weidman et al. [63], to test the stability of the steady flow f (η) = f 0 (η) and θ(η) = θ 0 (η) which fulfil the boundary value problem and boundary conditions (refer to (7)-(9)), we write by which ω is the eigenvalue of unidentified variables, while functions F(η) and G(η) are relatively small to f 0 (η) and θ 0 (η). The eigenvalue problems (15) and (16) result in an infinite group of eigenvalues ω 1 < ω 2 < ω 3 ...... that detect an early decay when ω 1 is positive, while an early growth of disruptions is observed when ω 1 is negative, which exposes the unstable flow. Substituting (18) into (15)- (17), 1 Pr and the boundary conditions are The heat transfer stability and steady-state flow solutions f 0 (η) and θ 0 (η) was implemented via τ → 0 , therefore F = F 0 (η) and G = G 0 (η) in (19)- (21). As a consequence, an early growth of Equation (18) is detected, and the subsequent generalized eigenvalue problem is recognized 1 Pr subject to The range of possible eigenvalues can be calculated by resting a boundary condition [64] of the present problem. In this study, we choose to repose F 0 (η) → 0 , and the linear eigenvalue problems (22)-(24) are disclosed as F 0 (0) = 1 for a fixed value of ω 1 . It is worth mentioning that the values of ω 1 are proficient in measuring the stability of the corresponding solutions f 0 (η) and θ 0 (η).

Results and Discussion
The mathematical computations of this research work were achieved by employing the bvp4c function in the MATLAB programming system subject to the governing ordinary differential Equations (7) and (8) together with the boundary conditions (9). A bvp4c method is a notable tool for solving the boundary value problem that has been extensively established by various researchers to clarify the boundary value concern. In order to ensure the obtainment of desired solutions, early estimation of the primary mesh point and variations step size is crucial. Also, a reasonable assumption of the thickness of the boundary layer along with effective preliminary approximation is essential for defining the non-uniqueness solutions. The consistency of the results generated in the present study is evaluated with those in [34,58], as accessible in Table 3, which is in excellent agreement.
In this study, we recognized the hybrid Al 2 O 3 -Cu/H 2 O nanofluid by dispersing the first nanoparticle φ 1 (alumina) into the base fluid (water) followed by the second nanoparticle φ 2 (copper) with various amounts of volume fractions. It is important to declare that the dispersion of discrete nanoparticle Al 2 O 3 /Cu is capable of developing Al 2 O 3 -Cu/H 2 O and Cu-H 2 O nanofluids. Apart from that, the Prandtl (Pr) number is set to be fixed at Pr = 6.2 corresponds to water as the reference-based fluid, while as for the hybrid Al 2 O 3 -Cu/H 2 O nanofluid, the size of the nanoparticles is assumed to be standardized, and the thermophysical properties effect of nanoparticles agglomeration is ignored. The main purpose of the present study is to examine the influence of the control parameter such as the nanoparticles volume fraction (φ 1 , φ 2 ), the unsteadiness parameter (ε), the velocity slip parameter (γ), and the Biot number (Bi) towards the coefficients of skin friction variations ( f (0)) and the heat transfer rate (θ (0)). The alumina and copper nanoparticles volume fraction in the present work is chosen within the range of 0.00 ≤ φ 2 ≤ 0.04 which motivated by the experimental work done by Suresh et al. [65] who conducted the synthesis, characterization of Al 2 O 3 − Cu/H 2 O nanocomposite powder for different volume concentrations 0.1%, 0.33%, 0.75%, 1%, and 2%. In their valuable study, the stability of the prepared nanofluids was determined by measuring the pH of nanofluids, and the nanofluid stability was found to diminish with increasing volume concentration. On another note, various values of the controlling parameter that has been used, are set within the following extent; 0.1 ≤ ε ≤ 0.2, 0.1 ≤ γ ≤ 0.4, and 0.2 ≤ Bi ≤ 0.7 to ensure the certainty of the obtained solutions. The present study also interested in witnessing the non-uniqueness solutions that appear in the ordinary governing differential of the ensuing problem, hence confirm the real and valid solutions through the stability analysis by utilizing the bvp4c approach in MATLAB operating system (MATLAB R2019b, MathWorks, Natick, MA, USA). Table 3. Evaluation values of f (0) when ε = γ = Bi = 0, by certain values of λ. The generated results of Equations (7) and (8) perceive the non-uniqueness (dual) solutions together with the boundary conditions (9) to a specific scope of λ c where λ c manifests the non-uniqueness solutions meeting point, and this too is regarded as a critical point. There is a single solution at a critical point, and that is the distinct line. The flow separation happens to occur after the critical point, and the flow is no longer laminar and does not obey the principle of boundary layer theory. According to the numerical results attained in this research work, it is proven that a non-uniqueness solutions appear, namely first and second solutions, as depicted in Figures 2-15. The existence of the non-uniqueness solutions contributes to the analysis of solution stability so that one may be able to verify the theoretically relevant solution. The first solution, however, is predicted to be reliable and fundamentally exist in practice. The smallest eigenvalues ω 1 for some values of λ when ε = 0.1, φ 1 = φ 2 = 0.02, Bi = 0.2, and Pr = 6.2 with γ = 0.1, 0.2, 0.4, are tabulated in Table 4. The values of ω 1 is noted approaching zero in the first and second solutions for λ → λ c in selected cases of γ. Hence, the authors may assure that the stability formulation and practices of the current problem are accurate and reliable. The flow is measured unstable if ω 1 is negative because it implies an early development of disruptions that leads to flow separation. The smallest eigenvalue ω 1 with positive value implies the flow is practically attainable and sustainable, which elucidates the solution stabilizing property to overwhelm the allowing disruptions. Also, it denotes an initial deterioration of disruptions that appear. The positive value of ω 1 describes the stable mode of the flow as in the first solution. The significance of dual solutions that contribute to stability analysis is in revealing the other possibilities of flow behavior, which may be useful for future references in the extrusion process. For example, in the present study, we have presented that the second solutions are unstable solutions with negative eigenvalues. These unstable solutions with the negative eigenvalues infer the growth of the disturbance in the solutions, especially at some range of the stretching/shrinking parameter (λ < −1.0). When the shrinking rate increases, it limits other external forces' effect at the sheet and eventually showed the opposite behavior of the transport phenomena than the first solutions in the same range of λ < −1.0, as the respective parameter varies. Moreover, any boundary value problem can generate more than one solution because of the nonlinearity in the boundary value problem; for example, see the mathematical model (7)- (9). Also, changes in the governing parameter values cause bifurcations in solutions that yield the existence of dual solutions [1]. Therefore, the mathematical model (7)-(9), which obeyed the boundary layer assumptions, managed to exhibit the variety in the fluid flow and heat transfer behavior through the uniqueness and existence of the solutions.
of disruptions that appear. The positive value of 1 ω describes the stable mode of the flow as in the first solution. The significance of dual solutions that contribute to stability analysis is in revealing the other possibilities of flow behavior, which may be useful for future references in the extrusion process. For example, in the present study, we have presented that the second solutions are unstable solutions with negative eigenvalues. These unstable solutions with the negative eigenvalues infer the growth of the disturbance in the solutions, especially at some range of the stretching/shrinking parameter ( ) When the shrinking rate increases, it limits other external forces' effect at the sheet and eventually showed the opposite behavior of the transport phenomena than the first solutions in the same range of 1.0, λ < − as the respective parameter varies. Moreover, any boundary value problem can generate more than one solution because of the nonlinearity in the boundary value problem; for example, see the mathematical model (7)- (9). Also, changes in the governing parameter values cause bifurcations in solutions that yield the existence of dual solutions [1]. Therefore, the mathematical model (7)- (9), which obeyed the boundary layer assumptions, managed to exhibit the variety in the fluid flow and heat transfer behavior through the uniqueness and existence of the solutions.              Figure 2 presents the coefficient of skin friction ( f (0)) towards the stretching/shrinking parameter λ, which revealed a decline of f (0) with the occurrence of velocity slip effect (γ = 0.1, 0.2, 0.4) at the boundary. The same results are obtained in the previous literature, as reported by Dzulkifli et al. [59]. Based on a physical perspective, an improvement in the slip parameter reflects the fact that the vorticity produced by the stretching/shrinking velocity is gradually decreased, thus the vorticity stays restricted within the boundary layer for greater stretching/shrinking velocity with the same straining velocity of the stagnation flow, and subsequently, the steady solution is achievable for some broad values of λ, as stated by Mahapatra and Nandy [34]. Figure 2 emphasizes that the solution for a certain value of γ persist prior to a critical value λ = (λ c < 0) at which the boundary layer splits from the convectively heated stretching/shrinking sheet and the solution on the basis of the boundary layer approximations is not feasible. In short, no solution exists for λ c < −1.6851. Moreover, the expansion of γ results in the increment of |λ c | suggesting that the velocity slip parameter is efficiently influenced by raising the range of dual solutions. Thus, it is proven that the existence of a slip velocity impact may prolong the separation of the boundary layer. Figure 2 also highlights as the sheet is stretching at the rate of λ = 1, the value of f (0) = 0 which describes no frictional drag exerted at the convectively heated stretching/shrinking sheet. Figure 3 exposes an upward trend in −θ (0) when the velocity slip parameter arises on the convectively heated stretching/shrinking sheet, which is proportionate to the heat transfer rate. Those findings are in line with the remarkable work reviewed by Mahapatra and Nandy [34] and Khashi'ie et al. [66]. Evidently, in the case of the first solution, the rate of heat transfer presents an ascendant trend with an increase of γ while the second solution permits a top-down direction when the velocity slip γ enlarges. From the existing and current findings, the authors can conclude that the velocity slip contributes to the improvement of the heat transfer rate significantly, prior to this case study. However, the authors would also like to declare that such results may vary if different control parameters are taken into consideration. The distribution of velocity ( f (η)) and temperature (θ(η)) profiles over several values of γ in the case of a convectively heated shrinking sheet is certified in Figures 4 and 5. As demonstrated in Figure 4, it is confirmed that as γ increases, f (η) is subsequently decreased in the first solution and revealed an upward trend in the velocity profile in the second solution with the presence of velocity slip condition. As the slip occurs, the velocity flow near the sheet in no longer equal to stretching velocity. This is due to the slip condition where the pulling force of the stretching sheet is partly shifted to the fluid, resulting in the decrement of the velocity profiles in the first solution.
In short, the slip condition reduces the momentum transfer from the sheet to the fluid. A contradict results are obtained in θ(η) as the velocity slip effect is exaggerated, as shown in Figure 5. The first solution in the temperature profile distribution increases as γ intensifies, while it diminishes in the second solution. This may occur due to the appearance of slip boundary conditions on the wall, which triggers the hybrid nanofluid to retain the velocity on the walls, and such slip may prevent from total heat exchange of the hybrid nanofluid. It is therefore noticed when the slip coefficient is applied, the difference in temperature rate is increased. Additionally, the effect of convective heating progress that has been reflected in the present study, which controls the temperature of the stretching/shrinking surface also might contribute to this phenomenon. As a matter of fact, greater convection leads to increased surface temperatures, which permit the thermal impact to penetrate deeper within the hybrid Al 2 O 3 -Cu/H 2 O nanofluid. Figures 6 and 7 expose the coefficient of skin friction ( f (0)) and the rate of heat transfer (−θ (0)) of a conventional nanofluid (φ 1 = 0.02, φ 2 = 0) and hybrid nanofluid (φ 1 = 0.02, φ 2 = 0.02, 0.04) past a convectively heated stretching/shrinking sheet when φ 2 varies from 0.00 to 0.02 where ε = 0.1, γ = 0.2, Bi = 0.2, and Pr = 6.2. Figure 6 Figure 8. The velocity profile in Figure 8 clarifies that the momentum boundary layer thickness was diminished in response to the rise of φ 2 , thereby raising the velocity of the fluid and boosting the gradient of velocity. In fact, the thinner momentum boundary layer continues to evolve the wall shear stress as well as the convectively heated shrinking sheet, leading to an enhancement of f (0). As f (0) increases, the result implies the increase of the frictional drag exerted on the convectively heated shrinking surface, which may delay the boundary layer flow separation. Besides, Figure 6 also highlights when the sheet is stretching at the rate of λ = 1, the value of f (0) = 0 which explains no frictional drag exerted at the sheet surface. Meanwhile, Figure 7 illustrates an increasing trend of the heat transfer characteristic or −θ (0) when the values of φ 2 past a convectively heated shrinking sheet, and this trend holds true to the first solution but is in dispute with the results in the second solution. In essence, as the conventional Al 2 O 3 -H 2 O nanofluid becomes the hybrid Al 2 O 3 -Cu/H 2 O nanofluid, the heat transfer efficiency improves. This finding upholds the assumption of the convective heat transfer system can be improved by optimizing the nanoparticle concentration when φ 2 increased. The results obtained in Figures 6 and 7 are consistent with Waini et al. [20] and Zainal et al. [21], whereby adding the concentrations of hybrid nanoparticles may contribute to the improvement of the heat transfer rate, accordingly. The temperature profile in Figure 9 describes the temperature variations when the conventional Al 2 O 3 -H 2 O nanofluid becomes the hybrid Al 2 O 3 -Cu/H 2 O nanofluid in both first and second solutions. The incline in the temperature of hybrid nanofluid proliferates the thermal conductivity, which may be triggered by the extra energy dispersed through the increment of the nanoparticles volume fraction over the state of convectively heated stretching/shrinking sheet. Figures 10-13 show the impact of the unsteadiness parameter (ε) towards a convectively heated stretching/shrinking sheet when ε shifts from 0.1 to 0.2. The hybrid Al 2 O 3 -Cu/H 2 O nanofluid characteristic with regard to the coefficient of skin friction ( f (0)) is depicted as in Figure 10. Figure 10 captures that when the sheet shrinks, the increment in ε conclusively increases the trend of f (0) in the first solution. An increment in the unsteadiness parameter results in the reduction of the boundary layer thickness, as depicted in Figure 12 and consequently upsurge the velocity gradient on the convectively heated stretching/shrinking sheet, thus f (0) improves. Furthermore, the existence of nanoparticle volume fraction in the working fluid (Al 2 O 3 -Cu/H 2 O) might also trigger the increment of f (0) owing to an uplift of the hybrid nanofluid viscosity. This result is aligned with the study done by Ismail et al. [67]. Figure 11 presents the unsteadiness parameter effect towards the rate of heat transfer −θ (0). According to the generated results, the heat transfer rate increases when ε increases when the convectively heated sheet is shrinking. The forced convective heat transfer is indeed proportional to the effectiveness of nanofluid heat conductivity, hence raising the reduced local Nusselt number notably. The dimensionless velocity profiles f (η) with a different value of ε are depicted in Figure 12, where the presence of dual velocity profiles is also observed. As illustrated in Figure 12, the first solution increases proportionally to the increment of ε values while the second solution displayed contradict results of the first solution, possibly because of the enhancement in the unsteadiness of the flow. Meanwhile, the same trend of the solution in Figure 12 also reflected the graph of temperature distribution θ(η) in the convectively heated stretching/shrinking sheet with the existence of unsteadiness parameter, as portrayed in Figure 13. Apart from that, it is proven that the second solution in both profiles, i.e., velocity and temperature distribution showed larger boundary layer thickness in each solution of the unsteady cases than those of the first solution. Figures 14 and 15 depict the variants of heat transfer rate −θ (0) and temperature distribution profile θ(η) with a different value of the Biot number (Bi = 0.2, 0.5, 0.7) towards the convectively heated stretching/shrinking sheet. In Figure 14, the heat transfer rate shows an upsurge trend in the dual solutions along with the augmentation of Bi values. The Biot number signifies the conduction resistance ratio within the sheet to convection resistance at the sheet. The critical values of the different usage of Biot numbers suggest no substantial impact on the magnitude of −θ (0), as highlighted by Jusoh et al. [68]. An increment in the Biot number related to the improvement of convective heating is observed to decrease the fluid temperature proficiently in the first and second solutions, as displayed in Figure 15. This result is in contrast with the idea of a large Biot number representing larger internal thermal resistance of the sheet compare to the boundary layer thermal resistance since the temperature distribution profile spike as the value of Bi increases. However, note that the Biot number is specifically correlated to the coefficient of heat transfer h f ; therefore, it is conversely related to the thermal resistance of the current problem. Consequently, the heat resistance decreases as the Biot number increases, thereby increasing the heat transfer rate at the stretching/shrinking sheet and decreasing the temperature distribution (see Figure 15).

Conclusions
An analysis of the unsteady stagnation point flow of hybrid nanofluid over a convectively heated stretching/shrinking sheet incorporating the velocity slip impact on heat transfer was verified in this study. The results were generated by employing the bvp4c features in the MATLAB programming platform. The effect of diverse controlling parameters-namely, the nanoparticle volume fraction, the velocity slip, the Biot number, and the unsteadiness parameter-were reviewed. Our discoveries happen to determine that the presence of non-uniqueness solutions (first and second solutions) is verifiable within the hybrid Al 2 O 3 -Cu/H 2 O nanofluid for a specific range of control parameters, and the stability analysis authorizes the reliability of our first solution. The nanoparticle volume fraction increment improved both the skin friction coefficient and local Nusselt number in the hybrid nanofluid (Al 2 O 3 -Cu/H 2 O). From this, it is proven that the heat transfer rate improves when the ordinary Al 2 O 3 -H 2 O nanofluid becomes the hybrid nanofluid (Al 2 O 3 -Cu/H 2 O) by expanding the nanoparticles concentration. The recent addition of the velocity slip parameter at the boundary had encouraged a reduction in the skin friction coefficient and velocity profile. However, it increased the rate of heat transfer significantly. The temperature profiles escalate as the magnitude of velocity slip upsurge because of such slip may prevent from exchange of total heat in the hybrid nanofluid. An increase in the unsteadiness parameter consequently raises the velocity gradient on the convectively heated stretching/shrinking sheet, thus improving the skin friction coefficient. Besides, an upsurge in the Biot number intensity boosts the heat transfer rate since the Biot number is directly associated with the heat transfer rate coefficient. Thus, it is conversely interrelated to the thermal resistance of the designated problem. Apart from that, the critical values of the different practices of the Biot number imply no significant outcome towards the magnitude of the heat transfer coefficient.