Impact of Smoluchowski Temperature and Maxwell Velocity Slip Conditions on Axisymmetric Rotated Flow of Hybrid Nanofluid past a Porous Moving Rotating Disk

Colloidal suspensions of regular fluids and nanoparticles are known as nanofluids. They have a variety of applications in the medical field, including cell separation, drug targeting, destruction of tumor tissue, and so on. On the other hand, the dispersion of multiple nanoparticles into a regular fluid is referred to as a hybrid nanofluid. It has a variety of innovative applications such as microfluidics, heat dissipation, dynamic sealing, damping, and so on. Because of these numerous applications of nanofluids in minds, therefore, the objective of the current exploration divulged the axisymmetric radiative flow and heat transfer induced by hybrid nanofluid impinging on a porous stretchable/shrinkable rotating disc. In addition, the impact of Smoluchowski temperature and Maxwell velocity slip boundary conditions are also invoked. The hybrid nanofluid was formed by mixing the copper (Cu) and alumina (Al2O3) nanoparticles scattered in the regular (viscous) base fluid (H2O). Similarity variables are used to procure the similarity equations, and the numerical outcomes are achieved using bvp4c in MATLAB software. According to the findings, double solutions are feasible for stretching (λ>0) and shrinking cases (λ<0). The heat transfer rate is accelerated as the hybrid nanoparticles increases. The suction parameter enhances the friction factors as well as heat transfer rate. Moreover, the friction factor in the radial direction and heat transfer enrich for the first solution and moderate for the second outcome due to the augmentation δ1, while the trend of the friction factor in the radial direction is changed only in the case of stretching for both branches.


Introduction
The study of nanofluids piqued the interest of researchers because of different applications in science and technology, like hybrid-powered engines, fuel cells, measures of different heat transport with pharmaceutical, and micro-electronic dealings. Nanofluids are composed of nanoparticles, and the presence of ultrafine nanoparticles magnifies thermal conductivity. Choi and Eastman [1] first proposed the idea of nanofluids by dispersion of fluids comprising ultra-fine particles. Later, several researchers investigated the mechanism of nanofluid in different geometrical domains: such as Khan and Pop [2] addressed the flow problem past a stretchable sheet containing nanofluids. The effect of free buoyancy flow via a vertical plate induced by nanofluid saturated in a porous medium was inspected by Gorla and Chamkha [3]. Xu et al. [4] presented an exact solution of the time-dependent flow induced by thin fluid film subject to the stretchable sheet. Gireesha et al. [5] examined the impact of dust particles scattered in a nanofluid flow from a stretched surface. Krishnamurthy et al. [6] scattered fluid particles in a magnetic fluid flow induced by nanofluid past an exponentially stretchable sheet with viscous dissipation. Makarem et al. [7] numerically inspected the fluid flow and features of heat transfer with distinct nanofluids on a stretchable surface. Ghalambaz et al. [8] discussed the features of fluid flow and analysis of heat transfer through a vertical plate saturated in a porous media packed with nano-encapsulated suspensions. They observed that the heat transfer uplifts in the presence of nano-encapsulated materials. Wakif and Sehaqui [9] revised the two-phase model by utilizing the water-based nanofluid to investigate the features of metal nanoparticles with a magnetic effect. Hajjar et al. [10] inspected the free convective flow of time-periodic with heat transfer scattered nano-encapsulated materials in a cavity and utilized the finite element technique to get a solution. The impact of nanofluid on gravity-driven flow induced by thixotropic fluid comprising microorganisms through a vertical sheet was examined by Koriko et al. [11]. Gul et al. [12] scrutinized the influence of thermophoresis and Brownian motion on thin fluid film induced by Carreau fluid from a stretched sheet with couple stress. The impact of magnetic field on the unsteady thin film flow of nanofluid with irregular convective stream past an inclined stretched surface is inspected by Saeed et al. [13]. Mathew et al. [14] investigated the effects of nanoparticle shape and slips on the magneto stagnation-point flow induced by nanofluid with chemical reaction and thermal radiation. Gul et al. [15] investigated a 3D MHD steady flow of Casson nanofluid induced by gyrotactic microorganisms through the gap of a cone and disk and developed the homotopy method to find the result. Recently, Gul et al. [16] examined the dissipative flow of nanofluid past a time-dependent turning disc with the magnetic field. They observed that the magnetic factor declines the velocities in radial and transverse directions.
Generally, each nanofluid possesses only one nanoparticle, whereas the hybrid nanofluid consists of two different nanoparticles dispersed in a different element regular fluid or a combination of component regular fluids. Hybrid nanofluids are created by combining two distinct nanomaterials to improve thermal and rheological properties. The main objective of hybrid nanofluid advancement is to effectively manage heat transport phenomenon in the examination of flow characteristics field. It has a wide range of technological applications including microfluidics, damping, acoustics, naval, heat dissipation, dynamic sealing, and so on. Suresh et al. [17] investigated the effect of dissipation on the time-dependent free buoyancy consequences of a hybrid nanofluid flow through a circular tube. The analysis of mixed convection flow induced by hybrid nanofluid was experimentally inspected by Momin [18]. Takabi and Salehi [19] scrutinized the hybrid nanofluid within a sinusoidal enclosure. Devi and Devi [20] presented the numerical solution of hybrid liquids past stretchable surfaces. They discovered an enrich the heat transport, which was followed by an increase in nanoparticle concentration. Hayat and Nadeem [21] described an escalation in heat transfer provided by an Ag-CuO/water nano solution. Rostami et al. [22] investigated the mixed convective flow of a water-based silica-alumina hybrid nanofluid in conjunction with a vertical plate by imposing linearly unstable temperatures with analysis of heat transfer. Double solutions were presented for opposing flow as well as for assisting flow. A parallel problem was surveyed by Zainal et al. [23] by considering convective boundary conditions along with different nanomaterials. Acharya et al. [24] discussed the characteristics of Hall current through a moving disk by including the Cu/TiO 2 nanomaterials suspended into base fluid water along with subject to thermal radiation and magnetic effects. Recently, several researchers [25][26][27][28][29][30][31][32][33][34] explored the significance of hybrid nanofluid flow with different aspects.
Over the last decade or so, there was a lot of interest in the analysis of stagnation point flow past a rotating disc. Hiemenz [35] discovered an analytic solution to N-S equations describing 2D steady flow handled orthogonally by an extended flat plate. Later the aforesaid problem was extended by Homann [36] to include an axisymmetric flow case. Hannah [37] used the irrotational far-field flow to extend Homann's [36] classic stagnation-point flow on a flat surface to survey the flow dynamics versus a spinning disc. Alternatively, Agrawal [38] formed a novel axially symmetric stagnation point flow on rotating external flow and normal impact to an infinite or inestimable or vary far plane wall. Then, Weidman [39] prolonged Agrawal's problem numerically. The impact of magnetic field on a 3D rotated flow towards a stagnation-point through a stretched radial rotating disk was inspected by Turkyilmazoglu [40]. Weidman [41,42] scrutinized the Agrawal axisymmetric flow near a stagnation-point over a flat sheet and stretched sheet, respectively. Recently, Waini et al. [43] reconnoitered the Agrawal's flow and the features of heat transfer along with hybrid nanoparticles from a stretchable/shrinkable disc.
The component of tangential velocity at the surface is proportional to the shear stress wall is known as the slip boundary condition. Navier [44] coined the term slip boundary condition after discovering a relationship between shear stress and slip rate. Maxwell [45] proposed the simplest explanation of the phenomena of velocity slip, which depends on a gradient of velocity exerting normal to a creep term and the surface. Smoluchowski [46] later introduces the idea of temperature slip. Ramya et al. [47] investigated the effect of thermal and velocity slips on fluid dispersed in a nanofluid. Recently, Khashi et al. [48] scrutinized the impacts of convective boundary stipulations and slip on a 3D flow of a hybrid nanofluid past a stretchable/shrinkable sheet.
The goal of the current study is to look into axisymmetric rotating radiative flow and heat transfer induced by hybrid nanofluid through a stretchable/shrinkable sheet with Smoluchowski temperature and Maxwell velocity slip. The following are the key points of research of the current study:

•
The impacts of Smoluchowski temperature and Maxwell velocity slip on the axisymmetric rotating flow were not explored yet.

•
The radiation effect on the axisymmetric rotating flow induced by hybrid nanofluid was not yet considered. • Double solutions of axisymmetric rotating flow past a moving rotating disk in the presence of hybrid nanofluid were not presented before.

Mathematical Formulation
The configuration of the Agrawal flow problem is schematically shown in Figure 1, where the steady radiative axisymmetric flow of heat transfer along with rotational stagnationpoint induced by hybrid nanoparticles impinging radially a porous shrinkable/stretchable rotating disk is contemplated. The problem is initially expressed in terms of cylindrical coordinates (z, r, θ) considered in the following axial, radial, and azimuthal directions, respectively, with the associated component of velocities (w, u, v). The Agrawal flow is symmetric to the rθ−plane and also axisymmetric about the axial direction (z−axis), i.e., the variation along the coordinate θ (azimuthal direction) is ignored. The stagnation line is located at z = 0 and the region of the flow dynamics is in the upper half-plane. Therefore, the stretching/shrinking disk is rotating about the axial direction (z−axis) together with a fixed angular velocity ω. In addition, the hybrid nanofluid is composed of two dissimilar nanoparticles (i.e., alumina (Al 2 O 3 ) and copper (Cu)) along with viscous pure fluid (water). The physical properties of the binary hybrid nanomaterials are taken to be in equilibrium and no-slip occurs between them. It is assumed that the component of the free stream velocities is characterized by u e (r, z) = 2arz, v e (r, z) = 0 and w e (r, z) = −2az 2 , where a is a constant parameter measuring the strength of the Agrawal flow having units (LT) −1 , see Weidman [39]. Moreover, the surface velocity at the wall z = 0 is denoted by u w = a 2/3 υ 1/3 f r along with Maxwell slip velocity [45] having slip length 2−σ v σ v λ 0 ∂u ∂z were implemented to study the velocity slip effect and further v w = rω, here, ω is the unchanging rotational speed of the disc, however, the uniform mass flux velocity via the wall is equal to w 0 , where w 0 < 0 for suction and w 0 > 0 for blowing. The constant wall surface temperature T w along with the Smoluchowski slip temperature [46] is indicated by 2 ∂T ∂z , whereas the constant ambient temperature is signified by T ∞ with the case T w > T ∞ . Under the impact of the aforementioned assumptions, the modeled governing equations in form of PDEs are premeditated by (see Weidman [39]): ∂u ∂r u ∂v ∂r along with subject to the boundary conditions In the above governing equations, u, v and w are the component of velocities along r−, θ− and z−axes, P is the pressure, σ v is the tangential momentum accommodation coefficient, λ 0 is the coefficient of the main free path, T is the temperature, σ T is the thermal accommodation coefficient,γ is the specific heat ratio and λ is the constant stretching/shrinking parameters with λ > 0 for the stretching sheet, λ < 0 for the shrinking sheet, and λ = 0 for the static disk. Further, ρ hn f is the density, k hn f is the thermal conductivity, ρc p hn f is the heat capacity, µ hn f is the dynamic viscosity, and (σ) hn f is the electrical conductivity of the hybrid nanofluid, which is given as (see [19,20]).
where φ corresponding to the solid nanoparticle volume fraction and is equal to the sum of the two distinct nanoparticles such as φ = φ 1 + φ 2 , in which φ 1 signifies the copper (Cu) nanoparticles, φ 2 signifies to alumina (Al 2 O 3 ) nanoparticles, and φ = 0 signifies to a regular (viscous) fluid. Besides, ρ f ,k f , ρc p f , ρ 1 , ρ 2 , k 1 , k 2 , ρc p 1 and ρc p 2 are the densities, thermal conductivities, and heat capacitance of the (viscous) regular fluid and the hybrid nanoparticles, respectively. Thermophysical properties of the regular (viscous) fluid and both the distinct nanomaterials (copper (Cu) and alumina (Al 2 O 3 )) are written in Table 1. λ is the coefficient of the main free path, T is the temperature, T σ is the thermal accommodation coefficient, γ is the specific heat ratio and λ is the constant stretching/shrinking parameters with 0 λ > for the stretching sheet, 0 λ < for the shrinking sheet, and 0 λ = for the static disk. Further, hnf ρ is the density, hnf k is the thermal conductivity, ( ) ρ p hnf c is the heat capacity, μ hnf is the dynamic viscosity, and ( ) σ hnf is the electrical conductivity of the hybrid nanofluid, which is given as (see [19,20]).  In energy Equation (5) the last term is used for the radiative heat flux q r . Therefore, the simplified form of the term radiative heat flux using the Rosseland approximation for an optically (hybrid nanoliquid) thick layer one can write (see Hayat et al. [49]): where σ * is the Stefan Boltzman constant and k * is the mean absorption constant. Further, executing the well-known Taylor series about the point T ∞ , the corresponding fourth power of T 4 can be simplified as T 4 ∼ = 4T 3 ∞ T − 3T 4 ∞ by ignoring the higher-order term. To further, simplify the analysis of the given model, here, we introduce the following self-similarity dimensionless variables are (see Weidman [39]): where primes correspond the derivative with respect to the pseudo-similarity variable ξ and further the expression (9) leads us to yield: In the above Equation (10), the symbol S demonstrates the transparent factor with S > 0 and S < 0 correspond to suction and blowing, respectively, whereas S = 0 is for an impervious surface of the disk. Now utilizing the self-similarity variables in the governing equations, where Equation (1) is evidently satisfied, while the rest of the Eqs. are transmuted into following ODEs as follows: 1 with appropriate BCs are: The above similarity equations comprise some of the dimensionless influential parameters which are namely and symbolically premeditated as like Pr = µ f c p f /k f the Prandtl is the temperature slip parameter, and α A = ω 2 /a 4/3 υ 2/3 f corresponds the rotating disk factor, necessitating that α A is greater or equal to zero.
The important gradients of the problem are the following shear stress in the relative radial and azimuthal directions indicated by C f r and C f θ , respectively, and the local Nusselt number Nu r is defined as Using similarity transformation (9) in the above equations, we get where Re r = u w r υ f is called the local Reynolds number.

Temporal Stability Analysis
The present Agrawal hybrid nanofluid flow problem was solved numerically and admits double solutions (upper branch solution and lower branch solution). Therefore, the multiple branch outcomes for a structure of a different geometry of the problems were documented by Weidman et al. [50], and Chu et al. [51], where the upper branch solution is stable and physically trustworthy, while the lower branch solution is unstable and not physically acceptable. Keeping these published available work in mind, we should test the properties of these stated efforts by letting the two-point problem of Equations (11)- (13). Hence, we define here the new non-dimensional time variable Π = a 2/3 υ 1/3 f t. To test and ease the procedure of the stability analysis, rewrite Equation (9) with new nondimensional variables as follows: Now executing Equation (17), the governing Equations (11)-(13) along with BCs (14), take the following form: 1 and the altered BCs become To check the working method of the temporal analysis of the time-independent Agrawal hybrid nanofluid flow outcome F(ξ) = F 0 (ξ) , G(ξ) = G 0 (ξ) and θ(ξ) = θ 0 (ξ) satisfy the two-point problem (11) to (13), we write (see Weidman et al. [50], and Chu et al. [51]): Here, the notation Σ is an unknown eigenvalue parameter, and functions f (ξ), g(ξ) and q(ξ) are comparatively small to F 0 (ξ),G 0 (ξ) and θ 0 (ξ), respectively. Plugging Equation (22) into Equations (18)- (20) along with the boundary conditions (21), the following linear eigenvalue problem is obtained: and the BCs (21) becomes Solving the linear eigenvalue problems (23)(24)(25) one observes an infinite number of eigenvalues Σ 1 < Σ 2 < Σ 3 . . .. If we found the Σ (smallest eigenvalue) positive, then the flow is stable and if Σ negative then the flow is unstable.

Numerical Procedure of the Considered Scheme
In this segment of the paper, we showed the mathematical technique of the studied numerical scheme in comprehensive form, as well as the accuracy of the code for the supplied flow problem. The modified similarity set of ordinary differential Equations (11)-(13) along with the border condition (14) are formed in the extremely nonlinear form, which is quite difficult to solve analytically after applying the similarity variables (9). As a result, these aforementioned numbers of highlighted equations are solved approximately by employing a finite difference scheme, which is known as asbvp4c. This package is a legitimate built-in code for MATLAB software (MATLAB R2020b) that is based on the three-stage Lobatto IIIA formula and was well-described in detail by Shampine et al. [52] and Khan et al. [53]. According to these available articles, the set of higher-order similarity equations were converted into a first-order set of ODEs by incorporating new symbols or notations. To continue our working process of the scheme, let the variables are: Now applying the Equation (27) into Equations (11)-(13) with the subject BCs (14), we get the following set of ODEs in the first-order form, which can be written as follows: with appropriate ICs are In ongoing to this process of the problem, we recommended here the essential early estimates at the specific mesh point to grip the system of first-order dimensional forms of Equation (28) with the subjected proper ICs (29), respectively. The mesh size ∆ξ in ξ, and the thickness of the boundary layer ξ ∞ have to be adjusted for dissimilar values of the influential constraints to endure the certain goal of accuracy. This problem has two distinct branch solutions. Therefore, the guess for the solution branch of upper is straightforward, while the solution branch of lower required an appropriate guess, which is quite hard to find out but trying until to satisfy the far-field BCs.

Validation of the Dual Solutions Code
To check the validation, confirmation or accuracy of the given numerical scheme, a comparison was constructed and shown in this subsection of the work. Table 2 highlights the comparison of wall drag force in the radial direction for several values of α A (without the impact of hybrid nanoparticles, stretching/shrinking parameter, radiation parameter, velocity and temperature slip parameters, and mass flux parameter) with published work of Lok et al. [54]. From this table, we have established that the existing work is completely mapping with exceptional accuracy along with the subject available reported work for the outcome of the first branch solution. Therefore, it revealed an outstanding agreement or sound and give exceptional power to catch or discover the unavailable outcomes using the considered scheme.

Analysis of Results
The similarity equations contained distinct controlling parameters like radiation parameter R d , rotating disk parameter α A , velocity slip parameter δ 1 , stretching/shrinking parameter λ, mass suction parameter S, temperature slip parameter δ 2 , and the hybrid nanoparticles φ 1 and φ 2 . For the persistence of the code simulations, we guaranteed the following range values of the selected distinguished constraints like 0.0 < R d < 5.0, 0.0 < α A ≤ 0.5, 0.01 < δ 1 < 2.0, 0.0 < δ 2 < 2.0, 0.0 < S < 5.0, −5.0 < λ < 2.0,0.02 < φ 1 < 0.04 and 0.02 < φ 2 < 0.04, whereas, the Prandtl number Pr = 6.2. Furthermore, the consequence of these selected comprised parameters on the wall drag forces along the radial and azimuthal directions, eigenvalues, and heat transfer of the hybrid (Cu-Al 2 O 3 /water) nanoparticles for the two distinct (first and second) branch solutionsversus λ are portrayed in Figures 2-11, respectively, whereas the Tables 3 and 4 are also prepared for the computational values of the gradients for several distinct values of the pertinent included influential parameters. More accurately, the numerical values of the wall drag forces along the radial and azimuthal directions for the sundry values of α A , S, δ 1 , φ 1 and φ 2 are constructed in Table 3 when λ = −1.4 (shrinking sheet), R d = 1.5 and δ 2 = 0.5. Meanwhile, Table 4 exhibitions the heat transfer numerical values for the enormous distinct varying values of α A , S, δ 1 ,R d ,δ 2 , φ 1 and φ 2 when λ = −1.4 (shrinking sheet)and Pr = 6.2. Tables 5-7 are prepared to compute values of normal nanofluid, and pure fluid as well as the presence and absence of slip factors. From the outcome of the tables, it is observed that the friction factor in the radial direction and heat transfer enriches for the first and second branch solutions with higher values of α A ,S, φ 1 and φ 2 . Meanwhile, the friction factor reduces in both branches with higher values of S but it is decreases and increases with hybrid nanoparticles and then increases and reduces in the first and second branch solution with α A . Moreover, the wall drag force in the radial direction and heat transfer rises for the first branch and decline for the second branch with larger values of δ 2 whereas the shear stress along the azimuthal direction abruptly shrinkages for both branches. Besides, the heat transfer decays and develops for the first and second branch outcomes with superior values of δ 2 and R d , respectively. Physically, the fluid particles captivate additional heat transfer by booming the parameter R d as a consequence, the Nusselt number upsurges. Moreover, the shear stress and heat transfer rate uplift as the velocity slip as well as temperature slip parameters enhance. Furthermore, with the presence of hybrid nanofluids, the heat transfer enhances further. Therefore, the performance of the cooling will be efficient for the inclusion of hybrid nanoparticles. Furthermore, the heat transfers and the magnitude of the friction factors in the radial and azimuthal directions in the stable branch are higher for the nanofluid as compared to the working base fluid. and 3 respectively. Similarly, here the second branch outcome breaks at some specific point against the stretching/shrinking parameter for the larger values of 1 δ due to the non-rotating disk parameter.
Finally, Figure 11exposes the smallest eigenvalues of Σ for the numerous values of the stretchable/shrinkable parameter λ . From the constructed graph, the negative value of Σ specifies an initial growth of disturbance, and the flow is an unstable mode. On the other hand, the flow is said to be in a stable mode when the value of Σ are positive which creates the initial decay of disturbance. In addition, we can noticed that the value of Σ tends to zero either from the lower branch or the upper branch solution as the values of λ are approaching, λ C . This shows that the transitions from positive (stable) to negative (unstable) of Σ occur at the turning points.

Conclusions
In this idea of the manuscript, we studied the impact of Smoluchowski tempe and Maxwell velocity slip boundary conditions on thermal radiative axisymmetri tional stagnation-point flow induced by hybrid (Cu-Al2O3/water) nanoparticles im ing on a porous stretchable/shrinkable rotating disk with heat transfer analysis. Th ilarity Equations (11-13) along with BCs (14) were reduced from the governing PD exercising the pertinent self-similarity variables. The obtained set of similarity equ comprised distinct dimensionless controlling parameters. Therefore, the multipl  Table 3. Values of shear stress along radial and azimuthal directions for several values of selected parameters when λ = −1.4 (shrinking sheet), R d = 1.5 and δ 2 = 0.5.   Table 6. Values of shear stress along azimuthal direction for several values of selected parameters when λ = −1.1 (shrinking sheet), δ 1 = 0, φ 2 = 0, δ 2 = 0 and R d = 1.5. The deviations of the rotating disk parameter α A on the wall drag forces along the radial and the azimuthal directions and heat transfer − k hn f /k f + (4/3)R d θ (0) of the hybrid (Cu-Al 2 O 3 /water) nanoparticles for the two distinct (first and second) branch solutions against λ are bounded in Figures 2-4, respectively. From the graphs, we have observed that the wall drag force in the radial direction and rate of heat transfer upsurge for the first outcome and reduces for the second outcome with improving values of the parameter α A , while the trend of the shear stress in the azimuthal direction continuously progresses for both distinct branch outcomes. Physically, the deceleration of velocity in the azimuthal direction is due to the fact of higher values of the rotating disk parameter which gets moderately transmitted to the upcoming adjacent fluid layers. As a result, the wall drag force in the azimuthal direction elevates. The branch of first outcomes is continuously smooth and no breaking exists for higher α A (see Figures 2-4) while the line was breaking at some finite values of λ 0 of λ due to the rotating disk parameter. Further, the first and second branch outcome was constructed in these graphs owing to larger α A for the phenomena of stretching and shrinking parameter. In other words, the nonunique outcomes were possible to find for the case of λ < 0 and λ > 0, whereas the outcome was unique for λ = λ C and no solution was accessible to found for λ < λ C . For higher values of α A , the following three distinct critical points are obtained and were written in each window of Figures 2-4, respectively. These points are also emphasized in the zoom window of the pictures by a sequence number such as 1, 2, and 3, while it is also spotted by the solid black balls. For mounting values of α A , the absolute values of the bifurcation points elevate, which shows that the trend of the separation of the boundary layer declines. Figures 5-7 elucidate the influence of S on the wall drag forces along the radial and azimuthal directions and heat transport of the hybrid nanomaterials for the two individual (first and second) branch solutions against λ, respectively. The wall drag force in the radial direction enlarges for the first outcome and declines in the unstable outcome owing to the higher impacts of S. Physically, an improvement in S conveys the flow of hybrid nanoparticles near the surface of the rotating disk which can bring velocity profile down, and consequently, the drag coefficient along with the radial direction upsurges. On the other hand, the − k hn f /k f + (4/3)R d θ (0) develops in both solution branches for the higher consequences of the parameter S while the wall drag force along with the azimuthal directions abruptly declines for both branches and escalates near the bifurcation points for the unstable branch outcomes with S. Moreover, the second branch results of the key physical quantities of interest are breaking or disconnected through some finite values of λ 0 of λ when the values of S is increasing. The gap in the first solution curves is finer as related to the distance in the second solution curves for the distinct selected values of S. With the variant values of S, we have perceived the following bifurcation values −1.47182, −1.98727 and −2.57303, respectively. From these critical points, it is seen that higher values of S the solution domain expands faster towards the larger negative values of λ and therefore, the magnitude of the bifurcation values is also boosted up.

Upper Solution Lower Solution
Furthermore, the power of the velocity slip parameter δ 1 on the wall drag forces along the radial and azimuthal directions and heat transfer of the hybrid nanoparticles are depicted in Figures 8-10, respectively. The wall drag force in the radial direction and heat transfer enrich for the first solution and moderates for the second outcome due to the augmentation δ 1 , while the trend of the wall drag force in the radial direction is completely changed only in the case of stretching for both branches with higher values of δ 1 (see Figure 8). Generally, in the phenomenon of Maxwell slip BCs, the speed of the disk and the fluid particles are not the same at the disk surface which depreciates the fluid speed and generates a reduction in the velocity profile as a conclusion, the wall drag force in the relevant radial direction escalates. Alternatively, the wall drag force along the azimuthal direction decelerates for both branches of outcome with superior impacts of δ 1 while the tendency or behavior of the flow near the bifurcation values for the second branch rises owing to the increment in δ 1 . In addition, the impacts of the velocity slip parameter δ 1 on the gradients are very small, and therefore we can zoom a specific part of the graphs to show the difference between the first curves and second curves. The zooming portion of the graphs also comprised the associated bifurcation values for the corresponding values of δ 1 , which can also highlight in these figures by positions 1, 2, and 3 respectively. Similarly, here the second branch outcome breaks at some specific point against the stretching/shrinking parameter for the larger values of δ 1 due to the non-rotating disk parameter.
Finally, Figure 11 exposes the smallest eigenvalues of Σ for the numerous values of the stretchable/shrinkable parameter λ. From the constructed graph, the negative value of Σ specifies an initial growth of disturbance, and the flow is an unstable mode. On the other hand, the flow is said to be in a stable mode when the value of Σ are positive which creates the initial decay of disturbance. In addition, we can noticed that the value of Σ tends to zero either from the lower branch or the upper branch solution as the values of λ are approaching, λ C . This shows that the transitions from positive (stable) to negative (unstable) of Σ occur at the turning points.

Conclusions
In this idea of the manuscript, we studied the impact of Smoluchowski temperature and Maxwell velocity slip boundary conditions on thermal radiative axisymmetric rotational stagnation-point flow induced by hybrid (Cu-Al 2 O 3 /water) nanoparticles impinging on a porous stretchable/shrinkable rotating disk with heat transfer analysis. The similarity Equations (11)-(13) along with BCs (14) were reduced from the governing PDEs via exercising the pertinent self-similarity variables. The obtained set of similarity equations comprised distinct dimensionless controlling parameters. Therefore, the multiple outcomes for various involved influential parameters on engineering physical quantities of interest for the first and second branch outcomes were demonstrated through various dissimilar graphs and as well as in a tabular form. The main findings of the considered simulations are summarized as follows: Figure 11.Smallest eigenvalues  for different values of  .

Conclusions
In this idea of the manuscript, we studied the impact of Smoluchowski temperature and Maxwell velocity slip boundary conditions on thermal radiative axisymmetric rotational stagnation-point flow induced by hybrid (Cu-Al2O3/water) nanoparticles impinging on a porous stretchable/shrinkable rotating disk with heat transfer analysis. The similarity Equations (11)(12)(13) along with BCs (14) were reduced from the governing PDEs via exercising the pertinent self-similarity variables. The obtained set of similarity equations comprised distinct dimensionless controlling parameters. Therefore, the multiple outcomes for various involved influential parameters on engineering physical quantities of interest for the first and second branch outcomes were demonstrated through various dissimilar graphs and as well as in a tabular form. The main findings of the considered simulations are summarized as follows:  The wall drag force coefficient in the radial direction and heat transfer upsurges for the branch of the first outcome and reduces for the branch of the second outcome owing to higher values of A  . However, the shear stress along the azimuthal direction elevates for both branches with A  , while the trends or behaviors are inverted near the bifurcation values for the branch of second solution curves.  Improvement in the mass suction parameter S displays an enhancement in the heat transfer, but the reduction in the shear stress along the azimuthal direction for both solution branches while the wall drag force in the radial direction rises for the first branch and declines for the second branch curves.  The heat transfer develops with radiation parameter for both solution branches, while shrinkages due to the temperature slip parameter 2  .
 The shear stress along the radial direction upturns for the first solution and drops for the second solution with an enhancement in 1  , while the wall drag force in the azimuthal direction diminished for both solution branches, but heat transfer upsurges for the first branch and decays for the second branch with higher values of 1  .
 An improvement in the heat transfer is also remarked with the larger impacts of hybrid nanoparticles. On the other hand, the shear stress in the radial direction augments for the first branch and reduces for the second branch with a development The wall drag force coefficient in the radial direction and heat transfer upsurges for the branch of the first outcome and reduces for the branch of the second outcome owing to higher values of α A . However, the shear stress along the azimuthal direction elevates for both branches with α A , while the trends or behaviors are inverted near the bifurcation values for the branch of second solution curves.

Conclusions
In this idea of the manuscript, we studied the impact of Smoluchowski temperature and Maxwell velocity slip boundary conditions on thermal radiative axisymmetric rotational stagnation-point flow induced by hybrid (Cu-Al2O3/water) nanoparticles impinging on a porous stretchable/shrinkable rotating disk with heat transfer analysis. The similarity Equations (11)(12)(13) along with BCs (14) were reduced from the governing PDEs via exercising the pertinent self-similarity variables. The obtained set of similarity equations comprised distinct dimensionless controlling parameters. Therefore, the multiple outcomes for various involved influential parameters on engineering physical quantities of interest for the first and second branch outcomes were demonstrated through various dissimilar graphs and as well as in a tabular form. The main findings of the considered simulations are summarized as follows:  The wall drag force coefficient in the radial direction and heat transfer upsurges for the branch of the first outcome and reduces for the branch of the second outcome owing to higher values of A  . However, the shear stress along the azimuthal direction elevates for both branches with A  , while the trends or behaviors are inverted near the bifurcation values for the branch of second solution curves.  Improvement in the mass suction parameter S displays an enhancement in the heat transfer, but the reduction in the shear stress along the azimuthal direction for both solution branches while the wall drag force in the radial direction rises for the first branch and declines for the second branch curves.  An improvement in the heat transfer is also remarked with the larger impacts of hybrid nanoparticles. On the other hand, the shear stress in the radial direction augments for the first branch and reduces for the second branch with a development Improvement in the mass suction parameter S displays an enhancement in the heat transfer, but the reduction in the shear stress along the azimuthal direction for both solution branches while the wall drag force in the radial direction rises for the first branch and declines for the second branch curves.

Conclusions
In this idea of the manuscript, we studied the impact of Smoluchowski temperature and Maxwell velocity slip boundary conditions on thermal radiative axisymmetric rotational stagnation-point flow induced by hybrid (Cu-Al2O3/water) nanoparticles impinging on a porous stretchable/shrinkable rotating disk with heat transfer analysis. The similarity Equations (11)(12)(13) along with BCs (14) were reduced from the governing PDEs via exercising the pertinent self-similarity variables. The obtained set of similarity equations comprised distinct dimensionless controlling parameters. Therefore, the multiple outcomes for various involved influential parameters on engineering physical quantities of interest for the first and second branch outcomes were demonstrated through various dissimilar graphs and as well as in a tabular form. The main findings of the considered simulations are summarized as follows:  The wall drag force coefficient in the radial direction and heat transfer upsurges for the branch of the first outcome and reduces for the branch of the second outcome owing to higher values of A  . However, the shear stress along the azimuthal direction elevates for both branches with A  , while the trends or behaviors are inverted near the bifurcation values for the branch of second solution curves.  Improvement in the mass suction parameter S displays an enhancement in the heat transfer, but the reduction in the shear stress along the azimuthal direction for both solution branches while the wall drag force in the radial direction rises for the first branch and declines for the second branch curves.

Conclusions
In this idea of the manuscript, we studied the impact of Smoluchowski temperature and Maxwell velocity slip boundary conditions on thermal radiative axisymmetric rotational stagnation-point flow induced by hybrid (Cu-Al2O3/water) nanoparticles impinging on a porous stretchable/shrinkable rotating disk with heat transfer analysis. The similarity Equations (11)(12)(13) along with BCs (14) were reduced from the governing PDEs via exercising the pertinent self-similarity variables. The obtained set of similarity equations comprised distinct dimensionless controlling parameters. Therefore, the multiple outcomes for various involved influential parameters on engineering physical quantities of interest for the first and second branch outcomes were demonstrated through various dissimilar graphs and as well as in a tabular form. The main findings of the considered simulations are summarized as follows:  The wall drag force coefficient in the radial direction and heat transfer upsurges for the branch of the first outcome and reduces for the branch of the second outcome owing to higher values of A  . However, the shear stress along the azimuthal direction elevates for both branches with A  , while the trends or behaviors are inverted near the bifurcation values for the branch of second solution curves.  Improvement in the mass suction parameter S displays an enhancement in the heat transfer, but the reduction in the shear stress along the azimuthal direction for both solution branches while the wall drag force in the radial direction rises for the first branch and declines for the second branch curves.  The heat transfer develops with radiation parameter for both solution branches, while shrinkages due to the temperature slip parameter 2  .
 The shear stress along the radial direction upturns for the first solution and drops for the second solution with an enhancement in 1  , while the wall drag force in the azimuthal direction diminished for both solution branches, but heat transfer upsurges for the first branch and decays for the second branch with higher values of 1  .
 An improvement in the heat transfer is also remarked with the larger impacts of hybrid nanoparticles. On the other hand, the shear stress in the radial direction augments for the first branch and reduces for the second branch with a development The shear stress along the radial direction upturns for the first solution and drops for the second solution with an enhancement in δ 1 , while the wall drag force in the azimuthal direction diminished for both solution branches, but heat transfer upsurges for the first branch and decays for the second branch with higher values of δ 1 .

Conclusions
In this idea of the manuscript, we studied the impact of Smoluchowski temperature and Maxwell velocity slip boundary conditions on thermal radiative axisymmetric rotational stagnation-point flow induced by hybrid (Cu-Al2O3/water) nanoparticles impinging on a porous stretchable/shrinkable rotating disk with heat transfer analysis. The similarity Equations (11)(12)(13) along with BCs (14) were reduced from the governing PDEs via exercising the pertinent self-similarity variables. The obtained set of similarity equations comprised distinct dimensionless controlling parameters. Therefore, the multiple outcomes for various involved influential parameters on engineering physical quantities of interest for the first and second branch outcomes were demonstrated through various dissimilar graphs and as well as in a tabular form. The main findings of the considered simulations are summarized as follows:  The wall drag force coefficient in the radial direction and heat transfer upsurges for the branch of the first outcome and reduces for the branch of the second outcome owing to higher values of A  . However, the shear stress along the azimuthal direction elevates for both branches with A  , while the trends or behaviors are inverted near the bifurcation values for the branch of second solution curves.  Improvement in the mass suction parameter S displays an enhancement in the heat transfer, but the reduction in the shear stress along the azimuthal direction for both solution branches while the wall drag force in the radial direction rises for the first branch and declines for the second branch curves.  The heat transfer develops with radiation parameter for both solution branches, while shrinkages due to the temperature slip parameter 2  .  An improvement in the heat transfer is also remarked with the larger impacts of hybrid nanoparticles. On the other hand, the shear stress in the radial direction augments for the first branch and reduces for the second branch with a development An improvement in the heat transfer is also remarked with the larger impacts of hybrid nanoparticles. On the other hand, the shear stress in the radial direction augments for the first branch and reduces for the second branch with a development in the hybrid nanoparticles, but the impact of the parameter is contrary for the shear stress in the azimuthal direction.

Conclusions
In this idea of the manuscript, we studied the impact of Smoluchowski temperature and Maxwell velocity slip boundary conditions on thermal radiative axisymmetric rotational stagnation-point flow induced by hybrid (Cu-Al2O3/water) nanoparticles impinging on a porous stretchable/shrinkable rotating disk with heat transfer analysis. The similarity Equations (11)(12)(13) along with BCs (14) were reduced from the governing PDEs via exercising the pertinent self-similarity variables. The obtained set of similarity equations comprised distinct dimensionless controlling parameters. Therefore, the multiple outcomes for various involved influential parameters on engineering physical quantities of interest for the first and second branch outcomes were demonstrated through various dissimilar graphs and as well as in a tabular form. The main findings of the considered simulations are summarized as follows:  The wall drag force coefficient in the radial direction and heat transfer upsurges for the branch of the first outcome and reduces for the branch of the second outcome owing to higher values of A  . However, the shear stress along the azimuthal direction elevates for both branches with A  , while the trends or behaviors are inverted near the bifurcation values for the branch of second solution curves.  Improvement in the mass suction parameter S displays an enhancement in the heat transfer, but the reduction in the shear stress along the azimuthal direction for both solution branches while the wall drag force in the radial direction rises for the first branch and declines for the second branch curves.  The heat transfer develops with radiation parameter for both solution branches, while shrinkages due to the temperature slip parameter 2  .