Duality Solutions in Hydromagnetic Flow of SWCNT-MWCNT/Water Hybrid Nanofluid over Vertical Moving Slender Needle

Recently, the topic of convection of heat transfer has created an interest among researchers because of its numerous applications in the daily life. The objective of this paper was to study theoretically the problem of mixed convection boundary layer flow and heat transfer of single-wall carbon nanotube (SWCNT) and multi-wall carbon nanotube (MWCNT) in presence of hydromagnetic effects. The problem was initiated by formulating a mathematical model in partial differential equation (PDE) for the hybrid nanofluid flow with appropriate boundary conditions. The similarity equation was used to transform the PDE into an ordinary differential equation (ODE) and solved using bvp4c in MATLAB. The graphical results on variation of skin friction coefficient, Cf, local Nusselt number, Nux, shear stress, f″c and local heat flux, −θ′c with the effects of magnetic, M, size of needle, c, mixed convection parameter, λ and volume fraction of nanoparticles, φ were presented and discussed in detail. The study revealed that duality of solutions appears when the buoyance force is in opposing flow of the fluid motion, λ<0. The presence of M in hybrid nanofluid reduced the skin friction coefficient and heat transfer. On the other hand, the Cf and Nux increased as different concentrations of φ1 and c were added. It gives an insight into the medical field, especially in treating cancer cells. By means, it reveals that CNTs hybrid nanofluid shows high potential in reaching the site of tumors faster compared with nanofluid. A stability analysis has to be carried out. It is noticed that the first solution was stable and physically realizable.


Introduction
Fluid dynamics field have shown great growth over the years, especially in nanotechnology industries, starting with viscous fluid, and now the nanofluid, which was introduced by [1]. The nanofluid is classified as dispersion of small-dimension nanoparticles into the base fluid, which gives good thermophysical properties. In, [2] colloidal stability, different phases of nanofluid, and their rheological behavior towards the applications were reviewed. It can be widely employed in various aspects, especially in electronics and nanotechnologies. Over the last decades, a lot of researchers have contributed their knowledge in conducting experimentally or theoretically studying using nanofluid with different amount concentrations, sizes of particles, and surfaces under a variety of effects [3][4][5][6][7][8][9]. Findings reported that nanoparticles in the fluids enhance the performance of the heat transfer. The types of nanoparticles used were characteristically from oxides, carbides, metals, or carbon nanotubes. Besides, the exposure of the goodness in using carbon nanotubes (CNTs) has perhaps greatly contributed to the revolution in nanotechnology industries due to its thermophysical properties and has eventually evolved into one of the most widely studied nanomaterials. Carbon nanotubes are classified into two categories, i.e., multi-wall and single-wall carbon nanotubes (MWCNTs and SWC-NTs). The pioneered CNTs were first discovered by Russian scientists Radushkevich and Lukyanovich [10] in 1951 on MWCNTs. Later in 1991, a male Japanese scientist, Iijima [11], structured a helical microtubules of graphitic carbon model. Two years later, Iijima and Ichihashi [12] exposed their findings on CNT with single wall, also known as single-walled carbon nanotube (SWCNT). To date, numerous researchers seeme to be engrossed in understanding the characteristics of CNTs and their novelty in boundary layer flow and heat transfer. Two different types of CNTs, i.e., SWCNT and MWCNT, were incorporated in water for the flow model. It was observed that an excellent result was achieved when MWCNT was diffused in the base fluid. It gives lowest thermal capacity and highest heat transfer in contrast with other nanofluids. Moreover, Hayat et al. [13] explored carbon nanotubes surface thickness with heat transfer in peristaltic flow. Besides, Shafiq et al. [14] studied the boost implications of heat transfer radiation on Darcy-Forchheimer (DF) flow of nanotubes of carbon along a stretched rotating surface by the methodology of responding surfaces. Through their findings, it is noted that the sensitivity of SFC via SWCNT-water towards the permeability number is higher than the solid volume fraction for medium and higher permeability levels. Besides, [15] scrutinized experimentally on the thermophysical properties and stability of carbon nanostructures and metallic oxides nanofluid. The mixed convective flow over a heated curve surface by influence of magnetic and velocity slip were scrutinized by Acharya et al. [16]. A recent article by Kamis et al. [17] presented their work in thin-film hybrid carbon nanotubes over an unsteady stretching sheet.
Problems on convection have caught the attention of researchers and academics for the past few decades. Convection happens when the heat is transmitted, hence exchanging the heat from the hot to cool position. There are three modes of convection, namely forced, free (natural), and mixed convection. Forced convection can be described as an external source mechanism or heat mode in which fluid movement is produced. Few literacies on forced convection can be found in [18][19][20][21]. Natural convection, on the other hand, is the process by which heat is transmitted by variations in density in the fluid because of temperature gradients. The temperature difference affects the density of the fluid. A collection of research papers on natural convection can be seen in [22][23][24][25].
The most favorable mode in convection to discuss is mixed convection. Mixed convection is a combination of forced and natural convections. The application of mixed convection usually in the systems with small scale such as in most of the engineering devices. Numerical investigation on mixed convection was carried out by [26]. They investigated non-isothermal bodies subjected to a nonuniform free-stream velocity. In their paper, a flow was considered effectively pure (either force or free) if the heat transfer or the friction coefficient deviates by no more than 5% from the value associated with the completely pure flow. Later, [27] studied over a vertical plate. In their findings, they concluded that the increment of convective parameter enhances the skin friction coefficient and the heat transfer rate at the surface; [28] examined a vertical plate in a porous surface. Both papers also found that two branches of solution appeared in high opposing flow. [29] extended the problem from [26] by adding the effect of heat generation/absorption through his study.
Besides, the application on thin needle also drawn an interest among researchers. It has many utilizations in the boundary layer and heat-transfer processes that are associated with our daily lives. For instance, calculating tools test, hot wire of the anemometer, or protected thermocouple. The term "thin needle" refers to a body of revolution whose indistinguishable order of diameter develops as velocity or thermal boundary layer. The mixed convection boundary layer flow over a vertical needle has been studied by many researchers. [30] studied the natural convection from needles with variable wall heat flux. In [31], the mixed convection flow about slender bodies of revolution was studied; [32] studied the problem of free and mixed convection boundary layer flows on a vertical adiabatic thin needle with a concentrated heat source at the tip of the needle. Later, [33] investigated steady laminar viscous fluid for both assisting and opposing flow. The results highlighted the dual solution obtained at opposing and assisting flows. The study was continued by [34] in nanofluids. In their conclusions, they perceived a favorable pressure gradient in assisting flow due to the buoyancy force. Consequently, the flow was accelerated, leading to a larger skin fraction coefficient compared with non-buoyant or opposing flow case. In 2018, [35] contributed to the idea by adding the effect of magnetohydrodynamics. In their study, they chose Buongiorno's model to observe the effect of Brownian motion and thermophoresis towards the vertical thin needle. They revealed that increasing the value of the magnetic, size of the needle, and thermophoresis parameter reduced the heat transfer on the needle surface. Besides, the local Sherwood number decreased as the Brownian motion rate increased. The above studies showed that a significant impact on flow and heat transfer characteristics was observed on shape, size, and wall heat flux variation for the thin needles.
Technologists found that the mixture of two different types of nanoparticles (hybrid nanofluid), scattered in a base fluid, improves the heat capacity of the fluid. Many researchers have shown interest due to the numerous applications of hybrid nanofluids in various engineering field such as in nuclear safety, manufacturing, transportation, and cooling electronic heaters in modern electronic and computer devices. Experimental work was conducted by [36] on enhancing the thermophysical of the nanoparticles. In advance, [37] reviewed on hybrid nanofluid over the backward and forward step. They presented the result experimentally and numerically on the usage of hybrid nanofluids. [38] explored the significance of thermal radiation on mixed convective boundary layer flow of a hybrid (SiO 2 -MoS 2 /H 2 O) nanofluid. Other detailed publications for hybrid nanofluid can be seen in [39][40][41][42][43].
To the authors' concern, the study about the concurrent effects of hydromagnetic and mixed convection over a moving thin needle of hybrid carbon nanotubes has not been publicized previously in the literature. In view of this fact, the importance of the magnetic field appended to the CNTs can be seen in medical applications, for example, in treating cancerous cells. Moreover, in the tumor's treatment, one of the effective approaches is by injecting it with magnetic nanoparticles in the nearest blood vessel to the tumor while putting a magnet close to the tumor. Ironically, the particles react like heat sources in the existence of alternating nature of magnetic field [44]. Note also, the stability analysis is executed to demonstrate the steadiness of the first and second solutions obtained throughout our research. Hence, the findings throughout this study can give benefit to other fields, not only in medical applications.

Problem Formation
An incompressible laminar of two-dimensional steady boundary layer with a combination of SWCNT and MWCNT hybrid nanofluid passing through a moving vertical slander needle was used. The needle was considered slender if the thickness of the needle did not exceed that of the boundary layer above it. Figure 1 shows the schematic representation of the physical model and coordinate system where x and r are the axial and radial coordinates, respectively, in which the radius is given by r = R(x). The x-axis measures the leading edge of the needle in the vertical direction, while r is always normal to the x-axis. In this system, we assumed the temperature of needle surface is T w and the ambient temperature T ∞ where when T w > T ∞ corresponds to a heated needle (assisting flow) and when T w < T ∞ corresponds to a cooled needle (opposing flow). Meanwhile, the needle moves with a constant velocity, U w , in the same or opposite way to the free stream velocity, U ∞ . The magnetic field, B 0 , is imposed parallel to the direction of the flow. Since the needle is assumed slender, the pressure gradient along the body is neglected but the importance is on the effect of its transverse curvature [45]. Using the above assumption, the governing equations for continuity, momentum, and energy for this problem are written in the cylindrical coordinates of x and r as follows [34,46]: with boundary conditions (BCs) in which u and υ are the component of velocity for x axes and r axes, respectively α hn f = k hn f / ρc p hn f is the thermal diffusivity of the hybrid nanofluid. Meanwhile, g is the gravity acceleration and B = √ 2B 0 /x refers to the magnetic field applied to the fluid flow. The induced magnetic field was disregarded since the magnetic number of Reynolds, R m , is assumed to be very small. The hybrid nanofluid's thermophysical attributes are described in Table 1. Following [48,49], a set of thermophysical properties that relate to these parameters are simplified in Table 2. The formation of hybrid nanofluid was formed by mixing the MWCNT nanoparticles into 0.01 volume of SWCNT/water. It is worth highlighting that 0.01 volume of solid fraction SWCNT (ϕ 1 = 0.01) was added constantly to the water throughout most of the problem. Meanwhile, to produce a SWCNT-MWCNT/water, various volumes of solid fraction of MWCNT (ϕ 2 ) were added. Table 2. Hybrid nanofluid's thermophysical properties [48,49].

Thermophysical Hybrid Nanofluids
Density Thermal Conductivity The subscripts of hnf, nf, and f representi the three categories of fluids, which are hybrid nanofluids, nanofluids, and fluids, respectively. While, ϕ 1 and ϕ 2 for two different nanoparticles of solid volume fractions which ϕ 1 is for SWCNT and ϕ 2 is for MWCNT. Besides, s1 was definitely for SWCNT s2 for MWCNT nanoparticles and S = s1 + s2.
To solve the governing equations, similarity resolution was acquired to diminish the independent variable of x and r. Thus, the above partial differential equations were reduced into ordinary differential equations with η as the new independent variable. Introducing a stream function, ψ, where u = 1 r ∂ψ ∂r and υ = − 1 r ∂ψ ∂x , the composite velocity between the needle and the free stream flow was define as U = U w + U ∞ . Hence, the obtaining similarities are: Thus, Equation (1) was satisfied identically. Equation (6) signifies the size of the needle where Invoke Equation (5) into Equations (2) and (3), and set c = η, the reduction of momentum and energy equations are in the following form: 2 Pr where λ = Gr x /Re 2 x is the constant mixed convection parameter or Richardson number with λ > 0, λ = 0 and λ < 0 represent the assisting flow, forced convection flow, and opposing flow, respectively. In this context, f is the local Grashof number and Re x = Ux/ν f is the local Reynolds number. While ε = U w U is velocity ratio parameter between the needle and the free stream. Furthermore, parameter M = σB 0 /ρ hn f U is the magnetic field and Pr = ν f α f is a Prandtl number. Our interests in physical quantities are the skin friction coefficient, C f , and local Nusselt number, Nu x , which can be easily expressed as Using Equations (5) and (6) to solve Equation (10), we have

Stability Flow Solution
Ref. [50] conveyed the solution for stability analysis in his work. To test the fact that the two solutions are either stable or unstable, a stability analysis has to be conducted. Consider the unsteady form by adding ∂u ∂t and ∂T ∂t at the beginning of Equations (2) and (3), respectively, where t denotes the dimensionless time. Thus, we have ∂T ∂t And the new similarity transformations are in the following form; Note that a new dimensionless similarity variable, τ = 2Ut x , is introduced in order to model the problem for temporal stability analysis. It is worth addressing that the use of τ in relation to the initial value problem consistent with the solution that will be obtained in practice (stable and physically realizable). Apply Equation (14) into Equation (12) and Equation (13), we have 2 Pr subjected to BCs, To determine the stability of the steady flow solution f (η) = f 0 (η) and θ(η) = θ 0 (η) that satisfied Equations (7) and (8), we introduced these linearized eigenvalues equations such as in [51]: where F(η, τ) and G(η, τ) are set as the small relative to f 0 (η) and θ 0 (η) respectively while γ is an unknown parameter of eigenvalue. In order to identify the growth or decomposition of Equation (18), we need to set τ = 0 , so that the function of F and G can be written as F 0 and G 0 . Differentiating Equation (18) and equating it to Equations (15) and (16), we obtain the following linearized eigenvalue equations: 2 Pr with BCs; After the meticulous derivation of Equations (19)- (21), the stability of the corresponding steady flow solution of f 0 and G 0 was determined by the least eigenvalue γ. According to [52], after relaxation of the BC on F 0 (η) or G 0 (η), the range of possible eigenvalues can be resolute. In this problem, F 0 (∞) → 0 was chosen to be relaxed and a new boundary condition was presented, which is F 0 (0) = 1, in order to replace the relaxing condition. Results on stability were obtained using bvp4c MATLAB.

Numerical Computation
There is a lot of computational programming to solve the boundary layer problem. Some researchers are familiar with MAPLE, Fortran, or MATLAB. However, generally, researchers would like to use a higher-order method that is robust and capable of solving systems of equations. So, bvp4c in MATLAB was chosen because it offers a convenient and proficient in solving fairly sophisticated problems, especially for solving nonlinear systems of equations [53,54]. The algorithm relies on an iteration structure. BVP4C is, in particular, a coding that implements the Lobatto IIIa three-stage formula. This is a collocation formula which is formed by a polynomial collocation. It provides a C 1 -continuous solution that is fourth-order accurate uniformly in x ∈ [a, b]. The crucial part in solving the bvp4c is the variation step and early guessing of the mesh point. Besides, the efficiency will eventually depend on one's ability to provide the algorithm with an initial guess for the solution. We created two folders, namely, code a and code b, for the trial-and-error initial guess and continuous iterations that nearly to the initial guess, respectively. To initiate the computation, the reduction of Equations (7) and (8) together with the BCs (9) to a set of first order is essential. Let, The BCs in Equation (9) become, Some commands in handling the function such as @odeBVP and @odeBC are from the syntax of the solver sol = bvp4c (@OdeBVP, @OdeBC, solinit, options). These commands are specified by Equations (22) and (23) and boundary conditions of Equation (24). The numerical results obtained from the solver are transformed to graphs. Duality solutions exist to the problem and also satisfy the boundary conditions. Further information regarding on bvp4c can be found in [55].

Analysis of Result
The numerical computations were carried out for several values of considered parameters which are M, λ, c, ϕ, and ε in order to furnish a clear understanding to the problem. The range of values of the parameters are M(0 ≤ M ≤ 0.05), c(0.1 ≤ c ≤ 0.2), and ϕ 1 (0 ≤ ϕ 1 ≤ 0.03) and the fixed values of Pr = 6.2. Equations (7) and (8) along the BCs (9) were computed numerically using bvp4c in MATLAB. The bvp4c solver uses the finite difference method. It requires a person to set an initial guess supplied at the initial mesh point and changing the step size in order to obtain the specified accuracy. Table 3 outlines the data produced by the present code and those reported by [56,57] when M = λ = ε = ϕ 1 = ϕ 2 = 0. It shows a great agreement, and hence we are confident to present our results. Table 3. Comparison values of f (c) when M = λ = ε = ϕ 1 = ϕ 2 = 0 for some thicknesses of needles, c when Pr = 1.  Figure 2 portrays the result for three types of fluid, which are viscous fluid (ϕ 1 = ϕ 2 = 0), SWCNT/water (ϕ 1 = 0.01, ϕ 2 = 0), MWCNT/water (ϕ 1 = 0, ϕ 2 = 0.010), and SWCNT-MWCNT/water (ϕ 1 = 0.01, ϕ 2 = 0.01) on variation of shear stress, f (c), and local heat flux, −θ (c). As we observed, we can see that hybrid nanofluid was highest compared with the viscous fluid and nanofluids in f (c) and −θ (c). Furthermore, we noticed that hybrid nanofluid delayed the separation of the boundary layer. Besides, the dual solutions seemed to exist at ε < 0, where the flow moves in opposite directions. In summary, the range of the dual solution was at ε c < ε ≤ −1.5 while no solution was found at ε < ε c . The result on the influence of size needle, c on the variation of shear stress, f (c) and local heat flux, −θ (c) is plotted in Figure 3. It can be seen that the shear stress for c = 0.1 was higher compared to c = 0.2. This is because the interaction between the needle surface and the fluid particles decreases when the size of needle decreases, which leads to the reduction of the drag force. Furthermore, the decreasing size of the needle gives a high thermal conductivity, allowing the heat being conducted from the surface through the needle to increase. It is also noted that the range of the solution widens as c increases, hence shortening the separation of the boundary layer. In particular, the dual solutions exist when ε c < ε ≤ −2 for both c = 0.1 and c = 0.2, respectively. Variation shear stress, f (c), and local heat flux, −θ (c), with mixed convection parameter, λ, for different values of M are presented in Figure 4. It is noticed that the values of f (c) and −θ (c) decrease when the imposed M is larger, following the fact that the rising strength of M will offer larger resistance on the fluid particles and cause the heat to be generated in the fluid. Therefore, a drag force known as Lorentz force is formed. This force may oppose and decelerate the motion of the hybrid nanofluid consequently, and reduce the wall friction and temperature, respectively. Besides, the stronger M makes the boundary layer become skinny and also reduces the range of boundary layer separation. Meaning, the boundary layer separation becomes faster as the value of M increases. It can be seen that the critical values of λ, λ c increase, i.e., M = 0( λ c = −3.1051), M = 0.01( λ c = −2.9221), and M = 0.05( λ c = −2.2850). Duality of solutions exist when λ c < λ ≤ 0, a unique solution appears at λ > 0, and no solution was found at λ < λ c . The buoyance force effect can be seen in assisting flow λ > 0 if the Prandtl number is set lower. According to [58], the physical explanation behind this is that the effect of buoyancy force is higher in low Prandtl number fluid. Due to the low fluid viscosity, the velocity in the boundary layer increases as the assisting buoyancy force functions as a favorable pressure gradient. For opposing flow (λ < 0), the buoyancy force is contrary to fluid movement, and hence the fluid is delayed, thus reducing the magnitude of velocity profiles within the boundary layer significantly. Thus, a reverse flow happens for λ = −1 near the surface. On the other note, the effect of λ is relatively less in −θ (c) because λ does not appear in Equation (8). Quite interesting result is on variation of shear stress, f (c), and local heat flux, −θ (c), with λ for different values of ϕ 1 . It is observed from Figure 5 that the increasing of ϕ 1 increased the f (c) near the wall and −θ (c) on the surface. The results showed good agreement with the experimental studies conducted by [59,60]. Increasing of volume fractions nanoparticles led to the increasing of the surface tension of the hybrid nanofluid. The number of nanoparticles suspended in the base fluid increased, which moved more nanoparticles to the liquid surface and adjacent to each other. As a result, an attractive force called Van der Waals forces occurred between the molecules, resulting in higher surface tension and consequently increased the shear stress. The strong force exerted between molecules led to the higher −θ (c). Duality of the solution exists at the opposing flow λ < 0, between the range λ c < λ ≤ −1, a unique solution appears at λ > 0, and no solution was found at λ < λ c . Besides, it can be seen clearly from the graph that the critical value, λ c , decreased as ϕ 1 increased, such as ϕ 1 = 0 (λ c = −2.8431), ϕ 1 = 0.01 (λ c = −2.9221), and ϕ 1 = 0.03 (λ c = −3.0880). Next, Figure 6 shows the influence of ϕ 2 on variation of C f and Nu x . Increasing ϕ 2 will consequently increase C f as ϕ 1 increases. Viscosity rate changes were barely observed at low fractions, likewise at high volume fractions. This is because the deposited of nanoparticles molecules increased the flow resistance, which resulted in increasing the shear stress and heat transfer. The behavior of C f and Nu x with ϕ for different values of c is depicted in Figure 7. It can be observed that C f and Nu x accelerated upwards with the increase of c as the value of ϕ increased. Note that both graphs are affected by the variation of c, as we can see in Equation (11). The friction force on the needle surface and fluid flow enhances due to the interaction between the needle's surface area and fluid particles when the size of the needle becomes bigger. The Nu x for c = 0.2 is higher compared with when c = 0.1. This follows the fact that heat transfer increases as the needle diameter increases because a large surface area requires higher temperature.  Figure 8 shows the influence of M on variation of C f and Nu x with ϕ. Both graphs portray that the value of C f and Nu x decreases with the increase in M and also increase with higher values of ϕ. It can be envisioned from the graph that the flow is flowing far from the surface instead of it became nearer to the surface when M is prominent. The reason for this behavior is because of the stronger interaction between the nanoparticles as ϕ increase but the motion was retarded by the Lorentz force that will gives disturbance to the flow. Thus, enhancing M tends to decrease both of C f and Nu x .   Figure  9 depicts the effect of M on both f (η) and θ(η). It can be observed that the Lorenz force retard the motion of the fluid, which decreased the velocity as M increased, but opposite trait happened on θ(η), where the force generated and consequently increased the heat as the value of M increased. However, Figure 10 clearly shows that the increasing of c slowed down the velocity but augmented the temperature of hybrid nanofluid for both branches. Furthermore, the effects of ϕ 2 on f (η) and θ(η) can be seen in Figure 11. In addition, all these profiles addressed and fulfilled the endpoint of BCs Equation (8) asymptotically.   Figure 9 depicts the effect of M on both f (η) and θ(η). It can be observed that the Lorenz force retard the motion of the fluid, which decreased the velocity as M increased, but opposite trait happened on θ(η), where the force generated and consequently increased the heat as the value of M increased. However, Figure 10 clearly shows that the increasing of c slowed down the velocity but augmented the temperature of hybrid nanofluid for both branches.
Furthermore, the effects of ϕ 2 on f (η) and θ(η) can be seen in Figure 11. In addition, all these profiles addressed and fulfilled the endpoint of BCs Equation (8) asymptotically.   The duality solutions obviously occurred at most in our numerical graph. Due to that, a stability analysis has to be conducted to test which explanations lead to a stable and reliable solution. To initiate them, we applied Equations (19)-(21) along with the new boundary condition into bvp4c in MATLAB software. In defining stability of the two solutions, the γ value is important. Numerical values obtained are written in Tables 4 and 5. Clearly, we can see that the first solution yielded a positive value of γ which is similar with the aforesaid [61][62][63]. However, an intriguing pattern was noticed at the second solution, where we also obtained positive values of γ for several λ → λ c . Figure 12 illustrates that for negative values of γ, the second solution specifies as an unstable solution; likewise for the positive value of γ, the second solution signifies as a stable solution. The explanation behind this is due to the rapid transition from turbulence flow to laminar flow. Generally, the boundary layer thickness for second solution is always thicker compared with the first solution. In this case, it is assumed that the laminar flow takes place in the thin boundary layer while the turbulence flow occurs in the thick boundary. For the present problem, the formation of the boundary layer was complicated with the presence of M, which gave an advantage for the turbulent flow to exchange to laminar flow.

Conclusions
This study was on mixed convection over a vertical moving slender needle which took into account the effect of M, c, λ and thermophysical properties (ϕ 1 and ϕ 2 ) on velocity, temperature, skin friction, and local Nusselt number. Therefore, the findings can be summarized in three important parts: 1.
Duality Solution

•
The presence of the duality solutions is prominent when in opposing flow λ < 0 • The range of duality of solutions is widely expanded under the influence of smaller M and c • Hybrid nanofluid possess faster boundary layer separation compared with viscous fluid and nanofluid.

•
The increment of c shortens the connection of first and second solution. • Stability analysis shows that first solution is favorable compared with the second solution.

Skin Friction
• Hybrid nanofluid reduced the value of skin friction as the value of M increased. • Hybrid nanofluid enhanced the skin friction as the concentration of the nanoparticles and c increased.

Heat Transfer
• Hybrid nanofluid reduced the value of heat transfer as the value of M increased. • Hybrid nanofluid enhanced the heat transfer as the concentration of the nanoparticles and c increased.

Acknowledgments:
The authors would like to express their deeply gratitude to the reviewers for all their meticulous, supportive also valuable comments and suggestions for this manuscript.

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