Rotating Flow in a Nanofluid with CNT Nanoparticles over a Stretching/Shrinking Surface

The steady three-dimensional rotating flow past a stretching/shrinking surface in water and kerosene-based nanofluids containing single and multi-walled carbon nanotubes (CNTs) is investigated. The governing equations are converted to similarity equations, and then numerically solved using MATLAB software. The impacts of rotational, suction, and nanoparticle volume fraction on the flow and the thermal fields, as well as velocity and temperature gradients at the surface, are represented graphically and are analyzed. Further, the friction factor and the heat transfer rate for different parameters are presented in tables. It is found that the heat transfer rate increases with increasing nanoparticle volume fraction as well as suction parameter in water and kerosene-based nanofluids of single and multi-walled CNTs. However, the increment in the rotating flow parameter decreases the rate of heat transfer. Multi-walled carbon nanotubes and kerosene-based nanofluid contribute to heat transfer rates better than single-walled carbon nanotubes and water-based nanofluid, respectively. A unique solution exists for the stretching surface, while two solutions are obtained for the shrinking surface. Further analysis of their stabilities shows that only one of them is stable over time.


Introduction
Researchers have begun looking for fluids with high thermal conductivity in the last few decades since traditional fluids have low thermal conductivity. In 1873, Maxwell studied the dispersion of solid particles with a size of millimeters and micrometers in traditional fluids. However, the existence of these particles caused a drastic drop in pressure, clogged flow channels, caused sedimentation of the particles to occur faster and erosion on the channel, and certain particles are too large for the microsystem [1]. Choi and Eastman [1] proposed a new fluid called nanofluid, in which solid particles with sizes ranging from 1-100 nm are dispersed in traditional fluids such as water, oil, and ethylene glycol. This fluid is significant because of its unique chemical and physical properties, as well as its high thermal conductivity, which can improve the heat transfer rate [2]. Waqas et al. [3] investigated the impact of nanoparticle shape's factor in porous media and discovered that platelet and cylindrical shapes had the highest thermal conductivity compared to brick and spherical shapes. A comprehensive collection of nanofluids studies can be found in [4][5][6][7][8][9].
Recently, many studies have considered carbon nanotubes as nanoparticles since their thermal conductivity is better and higher than those of traditional fluids. Carbon nanotubes (CNTs) are a sheet of graphene rolled into a cylinder introduced by Iigima in 1991. Carbon nanotubes are classified as single-walled (SWCNT) and multi-walled (MWCNT). They have a nanometer diameter size which is much smaller than other advanced semiconductor devices, and the length of nanotubes can exceed 1 micrometer. In addition, carbon nanotubes also offer unique electrical characteristics that other semiconductors lack. Further studies on unique features and review of carbon nanotubes can be found through investigation by Popov [10], Patel et al. [11], and Taherian et al. [12].
Due to their performance in optimizing thermal conductivity, researchers carried out many works to study the influence of carbon nanotubes availability and other effects in the fluids on the flow behavior. Khan et al. [13] investigated the fluid flow and heat transmission of CNTs nanofluid over a flat plate using three types of base fluids: water, engine oil, and kerosene. In comparison to water and kerosene-CNTs, they discovered that engine oil-CNTs have the highest heat transfer rate. Using a homogeneous model, Akbar et al. [14] investigated the stagnation-point flow of carbon nanotubes towards a stretching sheet with a magnetic field, velocity slip effect, and convective boundary condition. It was discovered that a higher magnetic field reduces the heat transfer rate at the surface and that single-walled carbon nanotubes have a higher heat transfer rate than multi-walled carbon nanotubes. Later, using similar CNTs nanofluid, Hayat et al. [15] explored the stagnationpoint flow over a nonlinear stretching surface with homogeneous-heterogeneous reactions. Hussain et al. [16] studied a three-dimensional flow of carbon nanotubes nanofluids over a nonlinearly stretching sheet with heat generation/absorption effect in a porous medium. Sreedevi et al. [17] investigated the heat and mass transfer characteristics of nanofluids containing CNTs towards a vertical cone in a porous medium by considering the magnetic field, chemical reaction, and suction/injection parameters. The findings revealed that increasing the nanoparticle volume fraction improves the heat transfer rate, while employing multi-walled CNT. as nanoparticles enhance the heat transfer over single-walled CNT. The latter conclusion was similar to Anuar et al. [18], who analyzed the mixed convection flow towards a moving plate. However, Anuar et al. [18] obtained a different result when the nanoparticle volume fraction increased. Reddy and Sreedevi [19] studied the effects of thermal radiation and single-walled carbon nanotubes inside a square chamber. From their observation, increasing these two factors enhances the heat transfer rate. Ramzan et al. [20] applied the Cattaneo-Christov heat flux model to analyze the impact of nonlinear thermal radiation and thermal stratification on the unsteady magnetohydrodynamics (MHD) flow between two stretching rotatory discs in nanofluid with carbon nanotubes. They discovered that nonlinear thermal radiation gives rise to the temperature of both SWCNT and MWCNT. Relevant studies on CNT nanofluids can be found in [21][22][23].
Investigation of the rotating flow towards a stretching surface is vital due to its application in a variety of scientific, technical, and product applications, designing and modeling jet engines, pumps, and vacuum cleaners, as well as geophysical flows [24]. Numerical studies on rotating flow have been explored by many researchers such as Singh and Sathi [25], who obtained an exact solution using Laplace transform for an unsteady rotating flow over an infinite plate. Wang [26] was the first to study the rotating flow over a stretching sheet using a perturbation method to identify the velocity profiles of small parameter values. Later, Nazar et al. [27] extended the work of Wang [26] by considering an unsteady flow. This study was expanded by Ali et al. [28] to the shrinking surface, and they obtained dual solutions. Further, Sreelakshmi et al. [29] included the magnetic field and thermal radiation effects in their study. They reported that the magnetic field, thermal radiation parameter, and rotation parameter increase the heat transfer rate. Rosali et al. [30] investigated the rotating flow over an exponentially permeable shrinking surface and found dual solutions. Rana et al. [31] implemented the Boungiorno model to investigate the MHD flow across a stretching surface in a nanofluid and numerically solved using a finite element method. It was revealed that the local Nusselt number and the local Sherwood number drop with the growth of the rotational parameter. By combining the Boungiorno and Tiwari-Das models, Bakar et al. [32] solved a rotating flow towards a shrinking surface with three nanoparticles types: copper, alumina, and titania with water as the base fluid using a shooting method. Different from Rana et al. [31], they found dual solutions, and the 3 of 20 local Nusselt number and the local Sherwood number increase with rotational parameter. Further, Krishna [33] explored the MHD effect in a Jeffrey fluid on a porous surface with hall and ion slip.
Nadeem et al. [34] explored a rotating flow towards a stretching surface utilizing copper and titania as the nanoparticles with water as the base fluid by applying the Tiwari-Das model. They discovered that increasing the volume fraction of nanoparticles boosted the heat transfer rate while decreasing it with the rotation parameter. Hayat et al. [15] reported the analytical results of rotating flow with carbon nanotubes over a stretching porosity surface. Dzulkifli et al. [35] continued the study of Ali et al. [28] by considering the unsteady rotation flow on a stretched surface in a nanofluid with copper as nanoparticles.
According to the findings, as the rotating parameter and nanoparticle volume percentage grew, the heat transfer rate dropped. Later, Salleh et al. [36] expanded the study by Nadeem et al. [34] to a permeable shrinking surface and obtained dual solutions. A stability analysis was carried out to inspect the stability of the solutions in the long run. As stated by Merkin [37] and Harris et al. [38], the upper branch is always stable, and this discovery is similar to the works made by Salleh et al. [36], Mustafa et al. [39], and Tshivhi and Makinde [40]. Hafeez et al. [41] explored the MHD rotating flow in a hybrid nanofluid using a finite element method. Mehdipour et al. [42] designed a rotating liquid sheet (RLS) contactor to explore carbon dioxide absorption from a gas stream using water-based SiO 2 and ZnO nanofluids. CO 2 absorption is more significant in the presence of ZnO than in the presence of SiO 2 , and tube rotation can improve absorption performance.
Acharya et al. [43] conducted a comparative investigation of SWCNTs and MWCNTs in a rotating regime over a stretching sheet in the presence of a magnetic field using the RK-4 method. They noticed that increasing the volume fraction of SWCNTs and MWCNTs improved the nanofluid temperature. Later, Shah et al. [44] analyzed a similar investigation over a stretching sheet using an analytical method. Noranuar et al. [45] studied a noncoaxial rotation flow towards a permeable moving disk in an MHD Casson nanofluid. Hussain et al. [46] extended the work of Nadeem et al. [34] by considering engine oil as the base fluid with copper and aluminum dioxide as nanoparticles over a stretching surface with slip condition. They found that copper nanoparticles react as heat carriers better than aluminum oxide. Manjunatha et al. [47] recently performed a numerical analysis to study the effect of a magnetic field on the fluid flow over a rotating disc in a nanofluid using single and multi-walled carbon nanotubes as nanoparticles and water as the base fluid. They reported that a unique solution exists, and the increase in the magnetic field caused an increase in the radial velocity and temperature.
Based on the previous studies mentioned above, many researchers focused on a rotational flow/disk over a stretching surface rather than a permeable stretching/shrinking surface and only obtained a unique solution. A permeable stretching/shrinking surface is vital in industrial applications since the surface may not be static such as polymer extrusion process, liquid crystal, and ceramic production [48]. Therefore, the present study aims to extend the work of Nadeem et al. [34] and Salleh et al. [36] by exploring the impact of rotation flow parameter, suction and nanoparticle volume fraction on the heat and mass transfer characteristics of carbon nanotubes (SWCNTs and MWCNTs) over a permeable stretching/shrinking surface in a nanofluid. Since no studies have considered this mathematical model with carbon nanotubes nanoparticles, the results obtained are genuine and can contribute to the development of this field. The numerical results are obtained using the bvp4c solver in MATLAB software, and a stability analysis is performed to investigate the stability of the solutions as time evolves.

Mathematical Model Description
A steady 3D rotating boundary layer flow of a carbon nanotube past a stretching or shrinking surface is examined, and its physical model is illustrated in Figure 1. In Figure 1, x, y, and z are the Cartesian coordinates with xand y-axes measured in the plane z = 0 where nanofluid is in the region z ≥ 0. The boundary z = 0 is elastic and linearly stretched or shrunk by two equal and opposite forces in the direction of the x-axis with the velocity U w = ax. The fluid is rotating at an angular velocity Ω in the z-direction. It is assumed that the nanofluid has a uniform wall and ambient temperatures, T w and T ∞ , respectively. Based on the above assumptions, the governing equations of mass, momentum, and energy are (see Nadeem et al. [34] and Salleh et al. [36]):

Mathematical Model Description
A steady 3D rotating boundary layer flow of a carbon nanotube past a stretching or shrinking surface is examined, and its physical model is illustrated in Figure 1. In Figure  1, , , and are the Cartesian coordinates with -and -axes measured in the plane = 0 where nanofluid is in the region 0. The boundary = 0 is elastic and linearly stretched or shrunk by two equal and opposite forces in the direction of the -axis with the velocity = . The fluid is rotating at an angular velocity Ω in the -direction. It is assumed that the nanofluid has a uniform wall and ambient temperatures, and , respectively. Based on the above assumptions, the governing equations of mass, momentum, and energy are (see Nadeem et al. [34] and Salleh et al. [36]): where , , and represent the velocity components in the -, -, and -directions, respectively, represents the fluid temperature, is the constant mass flux with > 0 indicates injection, while < 0 is for suction. Further, is the stretching or shrinking parameter with > 0 for stretching, < 0 for shrinking and = 0 for the static surface. In this study, we adopt the theoretical model proposed by Xue [49] as follows: The boundary conditions are where u, v, and w represent the velocity components in the x-, y-, and z-directions, respectively, T represents the fluid temperature, w 0 is the constant mass flux with w 0 > 0 indicates injection, while w 0 < 0 is for suction. Further, ε is the stretching or shrinking parameter with ε > 0 for stretching, ε < 0 for shrinking and ε = 0 for the static surface. In this study, we adopt the theoretical model proposed by Xue [49] as follows: In Equation (6), µ is dynamic viscosity, ρ is density, α is thermal diffusivity, k is thermal conductivity, C p is the specific heat capacity at constant pressure, ϕ is nanoparticle volume fraction parameter in which the subscripts nf, f, and cnt, respectively, refer to nanofluid, fluid, and carbon nanotube. The thermophysical properties of SWCNT and MWCNT for different base fluids namely, water (Pr = 6.2) and kerosene (Pr = 21) at the temperature 20-25 • C are given in Table 1 (Khan et al. [13]). Water and kerosene are chosen as the base fluids for comparison with the previous studies. Furthermore, water has the highest specific heat that can stabilize the temperature and is used as a cooling agent, while oil-based nanofluids are often utilized in industrial operations for lubrication or high-temperature applications [50]. The following similarity transformations (7) are introduced for solving Equations (2)-(4) together with the boundary conditions (5) that satisfy the continuity Equation (1). From Equation (7), the velocity mass flux can be defined as where S = f (0), S is the suction parameter in which S < 0 for injection and S > 0 for suction. Substituting Equations (6) and (7) into Equations (2)-(5) yields The BCs (5) become In Equations (9)-(12), prime denotes the differentiation w.r.t. similarity variable η, ω = Ω/a is the rotation parameter, and Pr is the Prandtl number.
The quantities of physical interest are the skin friction coefficients and the local Nusselt number (heat transfer rate at the surface). The skin friction coefficients in the xand ydirections, respectively, are given by: The local Nusselt number Nu x is defined as: Mathematics 2022, 10, 7 6 of 20 The shear stresses τ xz , τ yz and the heat flux q w at the surface are given by: Using Equation (15), Equations (13) and (14) become where Re x = Ux/ν f is the local Reynold number.

Results and Discussion
The system of nonlinear ordinary differential Equations (8)-(10) associated with the boundary conditions (11) was numerically computed for the diverse values of physical emerging parameters, namely, rotation parameter ω, stretching or shrinking parameter ε, nanoparticle volume fraction parameter ϕ, and suction parameter S. Two types of base fluids, which are water (Pr = 6.2) and kerosene (Pr = 21), were considered using both SWCNT and MWCNT. Numerical outcomes are obtained through the function bvp4c in MATLAB software are displayed in graphs and tables [52]. The numerical values of the velocity gradients at the surface of xcomponent f (0) and y-component h (0) were compared with the previous work of Wang [26] and Mustafa et al. [53] in the absence of stretching/shrinking parameter (ε = 0), nanoparticle volume fraction (ϕ = 0), and suction parameter (S = 0) for SWCNTs using water as the base fluid. Those values are listed in Table 2, and the acquired results were deemed adequate compared to previously obtained data by Wang [26] and Mustafa et al. [53]. Dual solutions are obtained as depicted in Figures 2-14, and it was worth performing a stability analysis as described in Section 2 to identify the stability of the solutions. The finding of stability analysis is discussed at the end of Section 3. Table 2. Comparison values of f (0) and h (0) for multiple values of rotation parameter ω when ε = ϕ = S = 0 and Pr = 6.2 (water) for SWCNT. ω Wang [26] Mustafa et al. [53] Present Study In this paper, the results in Figures 2-14 and Tables 3-6 are presented to discuss the impact of various physical parameters on the physical quantities of interest, such as the local skin friction coefficients Re 1/2 x C f x and Re 1/2 x C f y , and the local Nusselt number Re −1/2 x Nu x , which are proportional to f (0), h (0) and −θ (0), respectively, as well as the velocity and temperature profiles. Figures 2 and 3 present the variation of f (0) and h (0) for several values of rotation parameter ω versus stretching or shrinking parameter ε for SWCNT using water as the base fluid. These figures reveal that rising ω causes the gradients f (0) and h (0) to grow. However, it is noted that h (0) becomes zero when there is no rotation in the flow or when ω = 0. Further, a higher rotation rate causes fluid to move faster, which accelerates the velocity. The drag force between the fluid and the surface expedites and consequently increases the velocity gradient for both components. In addition, an increase in the values of f (0) and h (0) is due to the decrement in the momentum boundary layer thickness as ω increases (see Figure 7; Figure 8). Besides, as the stretching or shrinking parameter changes from the negative value to the positive one, both values of f (0) and h (0) seems to decrease. It happens because when the surface shrinks, less friction occurs on the surface than on the stretching surface. The rotating fluid forces the friction to appear in the flow significantly when the surface is being compressed. As a result, the velocity gradients rise for the shrinking surface compared to the stretching surface. Furthermore, the emergence of dual solutions is noticed when ε is in the region of ε c < ε ≤ −0.8. Meanwhile, unique solutions are obtained when ε > −0.8 and no solutions when ε < ε c .
The variation of f (0) and h (0) and −θ (0) for some values of nanoparticle volume fraction parameter ϕ versus suction parameter S are plotted in Figures 4-6 for MWCNT using kerosene as the base fluid. It is found from these plots that increasing the values of ϕ tends to diminish the f (0) and −θ (0). However, the reverse effect is noted for h (0) since adding more nanoparticles in the flow will slow down the collision between the molecules of the base fluid and nanoparticles, which in turn reduces f (0).
Since the surface is compressed or stretched in the x-direction, thus, the effect of ϕ on f (0) is more pronounced than h (0). As noticed in Figure 5, there is a slight increase in the value of h (0) when ϕ increases. Besides, the reduction in −θ (0), as displayed in Figure 6, is due to the thickening of the thermal boundary layer thickness as ϕ increases. This phenomenon prevents heat from the surface from being transferred to the surrounding fluid, and consequently, reduces −θ (0). It is also revealed that as S intensifies, f (0) and −θ (0) seem to increase, while the contrary phenomenon is noticed for h (0). The increment of −θ (0) is influenced by the surface texture where suction causes heat to be transferred from the wall to the fluid. It is observed that non-unique solutions occur for the particular  Figure 4 reveals that in the case of regular fluid (ϕ = 0), dual solutions exist in the range of S c = 1.9965 < S ≤ 2.31. Meanwhile, as the value of ϕ increases from 0 to 0.02 and 0.04, dual solutions are obtained in the range of S c = 2.0265 < S ≤ 2.34 and S c = 2.0584 < S ≤ 2.38, respectively. As can be seen in Figures 2-5, the range of solutions that exist are different based on the pertinent parameters involved. Figures 2 and 3 show that when different values of ω are plotted versus ε, the range of solutions increases with an increasing rotation parameter ω. However, Figures 4 and 5 depict that it decreases when the nanoparticle volume fraction ϕ increases. athematics 2021, 9, x FOR PEER REVIEW 9 for ℎ′(0) since adding more nanoparticles in the flow will slow down the collisio tween the molecules of the base fluid and nanoparticles, which in turn reduces ′′(0   athematics 2021, 9, x FOR PEER REVIEW 9 for ℎ′(0) since adding more nanoparticles in the flow will slow down the collision tween the molecules of the base fluid and nanoparticles, which in turn reduces ′′(0  Since the surface is compressed or stretched in the -direction, thus, the effect on ″(0) is more pronounced than ℎ′(0). As noticed in Figure 5, there is a slight incr        As can be seen in Figures 7 and 8, both velocities of xand y-components enhance with a higher rotation parameter for the first solution. An increase in the rotation rate will force the fluid to move rapidly and, in the meantime, accelerate both velocity components. This behavior will complicate the formation of the boundary layer thickness; thus, the boundary layer thickness becomes thinner. In contrast, the opposite phenomenon is reported for the second solution. Furthermore, it is found that the velocity of x-component decreases; meanwhile, the velocity of y-component and the temperature fields for the first solution show an increment with the increasing ϕ as displayed in Figures 9-11, respectively. enhance with a higher rotation parameter for the first solution. An increase in the rota rate will force the fluid to move rapidly and, in the meantime, accelerate both velo components. This behavior will complicate the formation of the boundary layer thickn thus, the boundary layer thickness becomes thinner. In contrast, the opposite pheno non is reported for the second solution. Furthermore, it is found that the velocity o component decreases; meanwhile, the velocity of -component and the tempera fields for the first solution show an increment with the increasing as displayed in ures 9-11, respectively.

First solution Second solution
Kerosene-MWCNT components. This behavior will complicate the formation of the boundary layer thick thus, the boundary layer thickness becomes thinner. In contrast, the opposite pheno non is reported for the second solution. Furthermore, it is found that the velocity o component decreases; meanwhile, the velocity of -component and the tempera fields for the first solution show an increment with the increasing as displayed in ures 9-11, respectively.    The presence of more suspended particles in the flow decelerates the velocity of xcomponent because the surface is only compressed in the x-direction. Since the surface is not compressed in the y-direction, the velocity of y-component increases due to less collision between the suspended molecules. The higher nanoparticle rate that thickens the thermal boundary layer thickness, as in Figure 11, will prevent heat diffusion from the surface to the fluid. As a result, the temperature of the flow rises. As the suction parameter enhances, the velocity of x-component for the first solution rises, while the opposite phenomenon is observed for both the velocity of y-component and temperature profiles as presented in Figures 12-14, respectively. Moreover, all the profiles obtained in this work satisfy the endpoint boundary conditions (11) asymptotically, thus supporting the numerical results gained for the current study. Further, the boundary layer thickness near the surface for the first solution is thinner than the second solution.
The presence of more suspended particles in the flow decelerates the velocity of component because the surface is only compressed in the -direction. Since the surface not compressed in the -direction, the velocity of -component increases due to less co lision between the suspended molecules. The higher nanoparticle rate that thickens th thermal boundary layer thickness, as in Figure 11, will prevent heat diffusion from th surface to the fluid. As a result, the temperature of the flow rises. As the suction paramet enhances, the velocity of -component for the first solution rises, while the opposite ph nomenon is observed for both the velocity of -component and temperature profiles a presented in Figures 12-14, respectively. Moreover, all the profiles obtained in this wor satisfy the endpoint boundary conditions (11) asymptotically, thus supporting the nume ical results gained for the current study. Further, the boundary layer thickness near th surface for the first solution is thinner than the second solution.   The presence of more suspended particles in the flow decelerates the velocity of component because the surface is only compressed in the -direction. Since the surface not compressed in the -direction, the velocity of -component increases due to less co lision between the suspended molecules. The higher nanoparticle rate that thickens th thermal boundary layer thickness, as in Figure 11, will prevent heat diffusion from th surface to the fluid. As a result, the temperature of the flow rises. As the suction paramete enhances, the velocity of -component for the first solution rises, while the opposite phe nomenon is observed for both the velocity of -component and temperature profiles a presented in Figures 12-14, respectively. Moreover, all the profiles obtained in this wor satisfy the endpoint boundary conditions (11) asymptotically, thus supporting the nume ical results gained for the current study. Further, the boundary layer thickness near th surface for the first solution is thinner than the second solution.   creases for both water-SWCNT and kerosene-MWCNT. This discovery demonstrate the suction effect stimulates particle collisions which causes more friction in the -d tion. Moreover, the suction impact on the skin friction in -component is more signi than -component since the surface is stretched in the -direction.   the suction effect stimulates particle collisions which causes more friction in the -d tion. Moreover, the suction impact on the skin friction in -component is more signifi than -component since the surface is stretched in the -direction.  x C f x and y-component Re 1/2 x C f y , for several values of ϕ, suction S, and ω for water and kerosene as the base fluids when ε = 1.0 (stretching case), respectively. It is observed from Tables 3 and 4, that an increase in ϕ tends to decrease Re 1/2 x C f x and Re 1/2 x C f y for both water-SWCNT and kerosene-MWCNT as reported by Anuar et al. [18]. This result differs from prior findings for copper, titanium, and alumina nanoparticles (see Nadeem et al. [34]; Mabood et al. [54]; Dinarvand et al. [55]). Furthermore, as the parameter S becomes more intense, the values of Re 1/2 x C f x increase, while Re 1/2 x C f y decreases for both water-SWCNT and kerosene-MWCNT. This discovery demonstrates that the suction effect stimulates particle collisions which causes more friction in the x-direction. Moreover, the suction impact on the skin friction in x-component is more significant than y-component since the surface is stretched in the x-direction.
In addition, both values of / and / increased with a higher tion rate. Physically, a high rotation rate causes the fluid to rotate continuously and erates more friction on the surface. The skin friction coefficients of the -and -com nents for SWCNTs are clearly higher than for MWCNTs. It is also worth noting that osene-based fluid has a greater -and -components skin friction coefficient than t of water-based fluid. Moreover, a positive value of the skin friction coefficient indi that the fluid exerts a drag force on the plate, whereas a negative sign suggests the o site.   tion rate. Physically, a high rotation rate causes the fluid to rotate continuously and erates more friction on the surface. The skin friction coefficients of the -and -co nents for SWCNTs are clearly higher than for MWCNTs. It is also worth noting tha osene-based fluid has a greater -and -components skin friction coefficient than of water-based fluid. Moreover, a positive value of the skin friction coefficient ind that the fluid exerts a drag force on the plate, whereas a negative sign suggests the o site.  In addition, both values of Re 1/2 x C f x and Re 1/2 x C f y increased with a higher rotation rate. Physically, a high rotation rate causes the fluid to rotate continuously and generates more friction on the surface. The skin friction coefficients of the xand y-components for SWCNTs are clearly higher than for MWCNTs. It is also worth noting that kerosene-based fluid has a greater xand y-components skin friction coefficient than those of water-based fluid. Moreover, a positive value of the skin friction coefficient indicates that the fluid exerts a drag force on the plate, whereas a negative sign suggests the opposite. when = 1.0 (stretching case). It is seen in Table 5 that the values of increase with increasing for both water/kerosene-SWCNTs and water/kerosene-MWCNTs However, the opposite observation was found for − ′(0), as displayed in Figure 6. I demonstrates that the thermal conductivity formula affects the heat transfer rate at the surface since it involves nanoparticle volume fraction (see Equations (6) and (15)). As a result of this circumstance, numerous studies are attempting to acquire the most effective thermal conductivity, as noted by Li et al. [56] and Fabre and Murshed [57]. However, a similar pattern was obtained for the values of      The numerical values of the local Nusselt number Re −1/2 x Nu x for some values of nanoparticle volume fraction ϕ, suction S and rotation parameter ω are given in Table 5 when ε = 1.0 (stretching case). It is seen in Table 5 that the values of Re −1/2 x Nu x increase with increasing ϕ for both water/kerosene-SWCNTs and water/kerosene-MWCNTs. However, the opposite observation was found for −θ (0), as displayed in Figure 6. It demonstrates that the thermal conductivity formula affects the heat transfer rate at the surface since it involves nanoparticle volume fraction (see Equations (6) and (15)). As a result of this circumstance, numerous studies are attempting to acquire the most effective thermal conductivity, as noted by Li et al. [56] and Fabre and Murshed [57]. However, a similar pattern was obtained for the values of Re −1/2 x Nu x and −θ (0), where a higher suction rate leading to an increase in Re −1/2 x Nu x and −θ (0). In contrast, the values of Re −1/2 x Nu x and −θ (0) decreased with the presence of the rotating fluid in the flow. This is due to the fluid rotation delaying heat transmission from the surface to the surrounding fluid, and the opposite trend occurs as suction increases. In the absence of suction and rotation parameters when ϕ = 0.01, SWCNT produces higher heat transfer rate than MWCNT for both water and kerosene base fluids. Without any disturbances like rotation in the fluid and suction at the surface, heat is quickly diffused through the thin wall (SWCNT) compared to the thick wall (MWCNT). In contrast, when there are suction and rotation parameters, a higher heat transfer rate is obtained for MWCNT. It is also found that kerosene has a higher heat transfer rate than water since it has a greater Prandtl number (Pr = 21) than water (Pr = 6.2).
Basically, the solution is stable if the trajectories do not deviate significantly in response to a slight perturbation and have the same values of f (0), h (0), and −θ (0) for η ∞ . In this study, the second solution is not stable since the solution is only valid for certain values of boundary layer thicknesses compared to the first solution. However, the existence of second solutions cannot be ignored since they may have physical sense in some circumstances. In this study, the stability analysis proposed by Merkin [37] was performed to determine which solution was stable as time evolved, by solving Equations (30)-(33) using bvp4c solver in MATLAB software. The values of the smallest eigenvalue γ for some values of rotation and stretching or shrinking parameters are generated and tabulated in Table 6 for single-wall carbon nanotubes. It was found that the values of the smallest eigenvalues γ for the first solutions were positive, whereas the second solutions had negative values. Referring to Equation (25), the flow was unstable since the initial growth of disturbance occurred for the negative value of γ. In contrast, the flow was stable when γ took the positive value, which implies an initial decay of disturbance in the system. It was noticed that as ε approached its critical value ε c , the smallest eigenvalue tended to zero for both first and second solutions, as shown in Figure 15. This criterion indicates that the transition occurs at the turning point. Figure 15 also shows the region for the stable and unstable solutions.

Conclusions
In the present paper, the numerical results of a steady rotational flow past a stretching/shrinking surface in a nanofluid with carbon nanotube as the nanoparticles (SWCNTs and MWCNTs) were obtained using bvp4c solver in MATLAB and reported in Tables and graphs. One of the most important discoveries in this paper is that dual solutions exist for the shrinking surface, and we found that only the first solution is stable in the long run. Moreover, the heat transfer enhancement is higher in water-MWCNTs than in water-SWCNTs. Kerosene-based carbon nanotubes have a higher heat transfer rate than water-based carbon nanotubes. Other findings of this study can be summarized as follows:

•
Increasing nanoparticle volume fraction tends to increase the heat transfer at the surface, while decreasing the temperature gradient in both nanofluids; • The heat transfer rate at the surface is increased by suction parameter, whereas it decreases by rotating flow parameter in both nanofluids; • In both nanofluids, increasing the nanoparticle volume fraction reduces the friction on the surface in the x and y directions, whereas rotational flow increases it; • The friction in x direction is higher than in y direction in both nanofluids.  Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: Not applicable.