Boundary Layer Stagnation Point Flow and Heat Transfer over a Nonlinear Stretching/Shrinking Sheet in Hybrid Carbon Nanotubes: Numerical Analysis and Response Surface Methodology under the Influence of Magnetohydrodynamics

: The present study aims to offer new numerical solutions and optimisation strategies for the fluid flow and heat transfer behaviour at a stagnation point through a nonlinear sheet that is expanding or contracting in water-based hybrid nanofluids. Most hybrid nanofluids typically use metallic nanoparticles. However, we deliver a new approach by combining single-and multi-walled carbon nanotubes (SWCNTs-MWCNTs). The flow is presumptively steady, laminar, and surrounded by a constant temperature of the ambient and body walls. By using similarity variables, a model of partial differential equations (PDEs) with the magnetohydrodynamics (MHD) effect on the momentum equation is converted into a model of non-dimensional ordinary differential equations (ODEs). Then, the dimensionless first-order ODEs are solved numerically using the MATLAB R2022b bvp4C program. In order to explore the range of computational solutions and physical quantities, several dimensionless variables are manipulated, including the magnetic parameter, the stretching/shrinking parameter, and the volume fraction parameters of hybrid and mono carbon nanotubes. To enhance the originality and effectiveness of this study for practical applications, we optimise the heat transfer coefficient via the response surface methodology (RSM). We apply a face-centred central composite design (CCF) and perform the CCF using Minitab. All of our findings are presented and illustrated in tabular and graphic form. We have made notable contributions in the disciplines of mathematical analysis and fluid dynamics. From our observations, we find that multiple solutions appear when the magnetic parameter is less than 1. We also detect double solutions in the shrinking region. Furthermore, the increase in the magnetic parameter and SWCNTs-MWCNTs volume fraction parameter increases both the skin friction coefficient and the local Nusselt number. To compare the performance of hybrid nanofluids and mono nanofluids, we note that hybrid nanofluids work better than single nanofluids both in skin friction and heat transfer coefficients.


Introduction
The need to enhance the transmission of heat over various surfaces of devices has continued to escalate in thermal management systems.Traditionally, base fluids such as water and kerosene work as heat transfer media (HTM) in order to transport heat to the environment.With rising industrial requirements to transfer heat at a faster rate, the innovative idea of mixing nanoparticles (NPs) with a base fluid has encouraged researchers to explore the most effective combinations for the preparation of nanofluids (NFs).
Because NPs have a large surface area with a nanometre scale of size, the investigations into the manipulation of NPs using experimental, theoretical, and numerical methods have still been ongoing since Choi's pioneering work [1] in this area.The fruitful development in nanotechnology has led to scientists not just suspending single types of NPs with a base fluid but successfully inventing two different types of NPs working with the base fluid.This process produces a new HTM called hybrid nanofluids (HNFs).The addition of HNFs to the flow has significantly improved the flow by promising a high thermal conductivity and, hence, increasing the heat-transmitting process in many applications.Turcu et al. [2] were reported as the first team to produce HNFs by combining NPs from carbon-based and non-metal materials.Reddy et al. [3] enumerated the possible uses of HNFs in numerous fields, including heat exchangers, engine cooling, and solar collectors.
The excellent achievement of HNFs in their narratives can be debated when it comes to determining the most effective amalgamation of NPs.Despite the inherent advantages of HNFs over traditional fluids in terms of heat transfer enhancement, the selection of an appropriate nanoparticle match for heat transfer enhancement in HNFs has been a subject of ongoing assessment and discussion.Concerning this challenge, Devi and Devi, Khashi'ie et al. and Khashi'ie et al. [4][5][6] recommended that researchers vary the mixture of NPs in order to produce the most productive HNFs that could potentially improve the flow and heat transfer behaviour.
Among the best ways to synthesise HNFs is by applying carbon-based materials.Sajid et al. [7] reported that Turcu and his team [2] became the first group to use multi-walled carbon nanotubes (MWCNTs), which are one of the nanoparticles from carbon-based particles that act as the main nanoparticles in HNFs.As CNT materials were identified as suitable materials for environmental applications by Navrotskaya et al. [8], demand for producing HNFs from these materials has increased significantly in nanotechnology.The increase in applying CNTs in HNFs is also supported, notably, by other relevant explorations of CNTs, which found them to contain excellent physical and thermal characteristics.However, there has been less discussion on the experimental and theoretical study of hybrid-based carbon nanotubes, especially in suspending single and multi-walled carbon nanotubes (SWCNTs and MWCNTs) as potential HNFs.The studies conducted by Hanaya et al., Sulochana et al.,  were identified in the literature as actively investigating the feasibility of manipulating hybrid SWCNTs and MWCNTs in order to study flow and heat transfer behaviour over different geometrical surfaces.From their research, all of them concluded that hybrid SWCNTs and MWCNTs performed better than single nanofluids, either SWCNTs or MWCNTs.
The boundary layer analysis near bodies of different geometrical shapes plays a vital role in investigating the flow behaviour of certain fluids.Svorcan et al.'s findings [13] suggest that studying the boundary layer has the potential to aid in the development of highly efficient boundary layer control devices.Due to the high application of shrinking and stretching sheets in cooling systems, researchers have focused attention on these surfaces.In a considerable proportion of the literature, most studies have explored these bodies, taking into account that the boundary velocity varies linearly with the fluid flow.Shateyi [14][15][16][17][18] were among the groups that studied the flow characteristics of a sheet linearly moved that was either stretching or shrinking using various mathematical models.In their models, it was assumed that the sheet moved at a linear velocity at the speed of the boundary layer.From our survey, the majority of researchers have been motivated by Crane et al.'s [19] study on the linear velocity along the stretched surface.However, realistically, it is important to explore the nonlinear variation of velocity over these surfaces, as emphasised by Vajravelu et al. [20].According to Nayakar et al. [21], a stretching or shrinking sheet that moves at a nonlinear velocity will employ at least a quadratic function velocity along the x-axis in a two-dimensional Cartesian coordinate system.The current understanding of the nonlinear velocity change occurring in an expanding or contracting sheet is considered insufficient for comprehensive discourse.This topic needs to be looked into more in the future, especially when it comes to investigating non-Newtonian fluid flow behaviour, like HNFs.
A search of the literature has revealed that the current works of flow over shrinking and stretching sheets have put NFs and HNFs in the spotlight.One of the numerous studies on flow in nanofluids was conducted by Rahman et al., .By taking the effect of magnetohydrodynamics into account in nanofluid flow over a nonlinear shrinking and stretching sheet, they analysed the possible region of non-unique solutions.They reported that the shrinking case contributed to producing multiple solutions while stretching only executed a single solution.The other attempt performed by Mahabaleshwar et al. [25] explained the thermal radiation and mass transpiration reactions in the boundary layer flow of nanofluids containing carbon nanotubes (CNTs) past linear shrinking and stretching sheets, and they found the existence of dual solutions in the skin friction coefficient.
Aly and his team [26] chose to look at the stagnation point flow of HNFs across a sheet that is either stretching or shrinking at a linear velocity because they had some interesting data from HNFs studies.They chose copper and alumina oxide (Cu-Al 2 O 3 ) as the main composition of HNFs, cooperating with water.They found that the duality solution appeared for defined parameters, and HNFs worked better than single nanofluids both in heat transfer and skin friction coefficients.The different positions of the stretching or shrinking sheet were scrutinised numerically by Khan et al. [27] by dealing with the vertical stretching or shrinking sheet.They observed that the process of separating the boundary layer became slower when the suction parameter increased.Although multidisciplinary research has been carried out on stretching/shrinking sheets, to date, the exploration of the flow of hybrid CNTs over a nonlinear sheet has not yet been organised by any researchers.Additionally, up to now, more work has been conducted on making models of how HNFs flow over a sheet that is linearly stretching or shrinking.Examples include the work by Zainal et al., .
The examination of stagnation point flow is currently the most extensively researched area pertaining to the flow characteristics across a surface that is undergoing either stretching or shrinking due to a wide range of industrial processes.In Merkin et al.'s [31] study, they found that the stagnation point was present on all body surfaces, regardless of whether they were subjected to stretching or shrinking.In the stagnation area, as seen in Figure 1, the fluid particles encountered a state of rest velocity in the vicinity of the stretching or shrinking sheet.According to Merkin et al. [31], this phenomenon offers a high transmission rate of heat in this region.
Magnetohydrodynamics (MHD) plays a crucial role in the improvement of many applications, including nuclear reactors, accelerators, and generators.Alekseev et al. and Ragupathi et al. [32,33] describe that MHD is used to explain the motion of electrically conducting and incompressible fluids in the presence of a magnetic field.To the best of our knowledge, not much research has been performed on computing the contribution of the numeric value and properties of the electrical conductivity of nanoparticles that are used in numerical studies that look into how MHD affects fluid dynamics.The comprehensive depiction of the MHD impact on flow motion may be hindered if computational analysis fails to include the electrical conductivity of nanoparticles in the model.The investigation of boundary layer CNTs with an MHD effect was carried out by Mahabaleshwar et al. [34].This team discovered that MHD produced a positive impact on the fluid flow when the flow passed through a linearly stretching or shrinking sheet.A recent study was undertaken by Mahesh and his colleagues [35], whereby they analysed the flow of carbon nanotubes (CNTs) in the boundary layer over a sheet that is either linearly extending or contracting.The research was performed under the influence of MHD.When the performance of SWCNTs and MWCNTs was compared, it was found that both types of nanotubes had better velocity profiles when the MHD level went up.However, to date, no investigation has yet been conducted on the boundary layer flow of hybrid CNTs (SWCNTs-MWSCNTs) across a nonlinear stretching or shrinking sheet in the presence of the MHD effect.
better velocity profiles when the MHD level went up.However, to date, no investigation has yet been conducted on the boundary layer flow of hybrid CNTs (SWCNTs-MWSCNTs) across a nonlinear stretching or shrinking sheet in the presence of the MHD effect.Motivated by the above studies, our team is planning to extend the original idea proposed by Anuar et al. [36].The previous work by Anuar et al. [36] presented the analysis of HNFs on different forms of nonlinear velocity over stretching or shrinking sheets.They reported that the greater the value of the polynomial function , the greater the contribution of detaching the boundary layer at a faster pace.However, this team did not use CNTs and performed optimisation procedures.In this study, we improve Anuar et al.'s model [36] by examining the impact of the MHD effect on the fluid flow characteristics and conducting an optimisation of the heat transfer.To build a different model from others, the MHD effect in that equation is counted together with the effect of the electrical Motivated by the above studies, our team is planning to extend the original idea proposed by Anuar et al. [36].The previous work by Anuar et al. [36] presented the analysis of HNFs on different forms of nonlinear velocity over stretching or shrinking sheets.They reported that the greater the value of the polynomial function n, the greater the contribution of detaching the boundary layer at a faster pace.However, this team did not use CNTs and performed optimisation procedures.In this study, we improve Anuar et al.'s model [36] by examining the impact of the MHD effect on the fluid flow characteristics and conducting an optimisation of the heat transfer.To build a different model from others, the MHD effect in that equation is counted together with the effect of the electrical conductivity of CNTs and water.Specifically, we use HNFs differently from the previous model by considering HNFs based on CNTs as the basis of our investigation.The decision made to enhance the prior model has a substantial impact on the new numerical outcome and the flow behaviour of HNFs.The use of MHD has found widespread application in several fields.Reddy et al. [3] have shown that this mechanism significantly influences the magnetic and electrical conductivity characteristics of fluid behaviour.In developing our model, we are taking a different approach from Jaafar et al.'s work [37] in order to investigate the MHD effect on the momentum equation and employ similarity variables that were suggested by Anuar et al. [36] before.In order to develop a highly effective model, our study focuses on the concept of MHD, as discussed in the article by Jaafar et al. and Khashi'ie et al. [37,38].Because Jaafar et al.'s model [37] did not take the flow at stagnation points into account and Khashi'ie et al. [38] studied a moving plate, our study, which is MHD-influenced, differs from that of Jaafar et al. and Khashi'ie et al. [37,38].
In creating the model, we adhere to the suggestion of Zafar et al. [39] since it is based on the assumption of flow in the laminar phase.Hence, the mathematical study of this model is conducted based on the Tiwari and Das model.Zafar et al. [39] concluded that utilising the Tiwari and Das model was the best option for illustrating the uniform fluid flow of nanofluids compared to other models.
For the purpose of ascertaining the most ideal magnetic parameter for this model, we provide a novel approach to quantifying the maximum heat transfer rate.To achieve optimal heat transfer, we use the response surface methodology (RSM) with the desirability function approach to design the numerical experiment.The use of the RSM approach has not been implemented in the models previously mentioned.The face-centred composite design (CCF) is employed to build a design of experiment (DOE) that incorporates other possible controllable factors.The objective of this DOE is to maximise the heat transmission capacity using a quadratic regression model.Matsui's study [40] indicates that the quadratic regression model is more adaptable than the linear regression model for assessing the interplay of several factors and responses to predict the best values of these interactions.To conduct the numerical analysis, we use the bvp4c function built in MATLAB, whereas the RSM analysis is carried out using Minitab.Since the integration of numerical approaches with RSM has garnered very little attention from scholars, this is an opportunity for further exploration of this combination.Both the numerical and RSM strategies that are built into this model may contribute to a comprehensive knowledge of the relationship between computational and experimental studies.

Mathematical Formulation
In order to govern the boundary layer equations, which are organised into nonlinear partial differential equations (PDEs), we establish the following assumptions: • The flow experiences an incompressible, steady, and laminar flow without a slip effect.• The flow is described in two-dimensional (2D) space, with x and y acting as the axes on the cartesian coordinates.• The sheet is extended and contracted with the free stream velocity U ∞ (x) and the stretching or shrinking velocity U w (x).• The terms U ∞ (x) and U w (x) can be written as U ∞ (x) = ax n and U w (x) = bx n , where a, b and n are positive constants, such that a, b > 0 and n > 1.
• The sheet is surrounded by the constant surface temperature T w and the constant ambient hybrid CNT nanofluids temperature T ∞ .• Temperature does not affect the thermophysical properties of hybrid CNTs.
• The size of SWCNTs and MWCNTs is not varied, and there is no issue of agglomeration in the flow.
To design our model, Anuar et al.'s study [36] is adopted, and the effect of the magnetic field on the strength B 0 in the dimensionless unit is taken into consideration in the flow, where B 2 0 x n−1 = B 2 .However, the influence of an induced magnetic field is neglected.By setting u and v as the velocity components that are dependent on x and y, our model with specific boundary conditions is defined as follows: The variables u and v in Equations ( 1)-( 4) represent the component velocities that are dependent on the Cartesian coordinates x and y.These components are written in non-dimensional units using the following equations: The terms µ hn f , ρ hn f and σ hn f in Equation ( 2) indicate the dynamic viscosity, density and electrical conductivity of hybrid CNT nanofluids, while α hn f in Equation ( 3) illustrates the thermal diffusivity of hybrid CNT nanofluids.

Non-Dimensional Similarity Ordinary Equations
To reduce the complexity of solving the above model, we introduce several similarity variables that were employed by Anuar et al. [36].These variables play a vital role in transforming the PDEs system into the non-dimensional ordinary differential equations (ODEs) model.Crank [41] asserts that by converting PDEs into ODEs, it is possible to examine the physical system of fluid flow and heat transfer using a less complex system.This claim was supported by Burden and Faires [42] who said that ODEs are more numerically stable than PDEs.The introduction of these similarity variables can be deployed in the following expressions: By inserting Equation (6) into Equations ( 2)-( 5), we discover the following system of ordinary differential equations (ODEs).
subject to boundary conditions: The term A i , for i = 1, 2, 3, 4, 5 can be recorded as below: where n, M and ε = a/b are the nonlinear, magnetic and stretching/shrinking velocity parameters, respectively, in dimensionless units.f and θ are functions of η.Hence, we can differentiate f with respect to η to obtain f ′ , f ′′ and f ′′′ , and differentiate θ with respect to η to obtain θ ′ and θ ′′ .Referring to Jaafar et al. and Khashi'ie et al. [37,38], the term M in Equation ( 10) can be denoted as To measure the physical quantities that can be applied in studying the application of the model, we define the local skin friction coefficient, C f , and the local Nusselt number, N ux , in the following equations: and where τ w and q w represent the surface shear stress and the surface heat flux, respectively.The terms τ w and q w can be written as and Applying the similarity variables from Equation ( 6), we obtain the reduced skin friction and heat transfer coefficients, respectively, in Equations ( 15) and ( 16) as follows: and where  [43].The mathematical correlation of the thermal conductivity of hybrid CNTs in this model is also drawn from the most up-to-date research on hybrid CNTs conducted by Aladdin et al. [11] and Haider et al. [44].The correlations and specific numeric values for the physical properties can be viewed in Tables 1 and 2.
Table 1.The thermophysical characteristics and correlations of hybrid CNTs dispersed in water.

Characteristics Correlations
Viscosity

Numerical Procedure
For the purpose of being able to numerically solve Equations ( 7)-( 9), the higher-order dimensionless ODEs are altered into first-order ODEs.This approach allows for the model to be entered into the bvp4c function in MATLAB.The transformation process is initiated by establishing the following equations: Therefore, the higher orders of dimensionless f ′′′ and θ ′′ are expressed as follows: The boundary conditions pointed out in Equation ( 9) turn into these mathematical expressions:

Optimisation Procedure
To conduct the optimisation procedures, we apply the Response Surface Methodology (RSM).According to Myres et al. [45], the adoption of RSM can boost the efficiency of experimental and numerical operations by offering optimum solutions based on statistical and mathematical concepts.Researchers such as Benim et al. [46] successfully manipulated the RSM to optimise the aerofoil profiles in wind turbines.
To begin the optimisation process, we first identify three factors or parameters that have the potential to impact the heat transfer performance.According to Equations ( 7)- (10), the volume fraction, magnetic and nonlinear parameters are selected as the elements in the input set.These parameters are recognised as the most significant inputs that may contribute to the model's performance.Additionally, the Tiwari and Das model emphasises that the volume parameter is the key factor in determining the fluid flow and heat transmission characteristics.Adhering to this model in developing the mathematical model has become the main reason for choosing this parameter as part of the input set.The matrix of these factors with different levels is illustrated in Table 3. Table 3 is organised according to the relative limitations of each parameter.[+1] = [1] stands for the highest level, [0] for the middle level and [−1] for the lowest level.Numerical simulations are used to produce the output responses of hybrid (y hCNT ) and mono (y CNT ) CNTs nanofluids by systematically altering the values of the factors.We also utilise a face-centred central composite design (CCF) to create an experimental design, with the number of trials N calculated using the formula below: where F and C are the number of factors and the centre point, respectively.
Note: The symbol [] denotes a coded component.
After screening the factors and their output variables, we develop the multivariate quadratic regression equations by modifying the quadratic model expressed previously by Myres et al. [45] in order to investigate the relationship between factors and their response variables.These equations are also formed to optimise the heat transfer for both hybrid and mono CNT nanofluids.The regression models of (y hCNT ) and (y CNT ) are written as follows: and The terms β 0 , β i , i = 1, 2, 3, β j , j = 4, 5, 6, β k , k = 7, 8, 9 and ϵ represent the mean value, the linear coefficients, the interaction between factors coefficients, the quadratic coefficients and the error term, respectively.The mathematical formulation of Equations ( 21) and ( 22), which include quadratic components, can be helpful in our analysis to determine the optimal amount of the heat transfer coefficient.According to Myres et al. [45], a quadratic model is suitable for fitting the process of an optimal issue.To determine the value of coefficients in Equations ( 21) and ( 22), we use the analysis of variance (ANOVA).The ANOVA data for the y hCNT and y CNT models are computed with the aid of the Minitab program.

Numerical Solutions
To validate the programming code developed in the bvp4c function (MATLAB), we compare the data with Anuar et al.'s model [36].By setting ϕ SWCNT = ϕ MWCNT = M = 0, n = 1, ε = −1.1, and Pr= 0.7, the model produces the reduced skin friction f ′′ (0) and the reduced heat transfer −θ ′ (0) as displayed in Table 4.The results presented in Table 4 indicate that the current model and the previous model are in good agreement.The result of this agreement gives us the opportunity to commence our investigation by examining the effects of magnetic fields on the heat transfer and flow characteristics of the hybrid CNTs.[36].

Anuar et al.'s Model [36]
The Present Solution

Anuar et al.'s Model [36]
The Present Solution In order to examine the effect of the magnetic parameter M on the flow separation, it is necessary to identify the region in which the duality solution occurs.It is possible to predict the point at which the boundary layer separates from the surface using the dual solution.We assume that the upper solution (the first solution) and the lower solution (the second solution) are the attached and separated flow solutions, respectively.To investigate this region, we set M to vary from 0 to 0.2 and −1.5 < ε < 1.5.The Waqas et al. model [47] serves as a guide for changing the variable M. The other parameters are kept at constant values, where Pr = 6.2, ϕ 1 = ϕ SWCNT = 0.01, ϕ 2 = ϕ MWCNT = 0.01 and n = 2.The first solution (solid line) and the second solution (dash line) appear when ε c ≤ ε < 0, where ε is in shrinking mode and ε c represents the critical point of producing the dual solutions.These phenomena are visually represented in Figure 2.An additional observation that can be made from this figure is that the boundary layer delays being separated as M increases.Increasing n also results in a noticeable lag in the detachment of boundary layers, as displayed in Figure 3.The aforementioned outcomes are produced through the configuration of various values for n and ε while holding all other variables constant.From Figures 2 and 3, the rise of M and n contributes to the widening range of f ′′ (0) and −θ ′ (0).The variations in velocity f ′ (η) and temperature θ(η) profiles for hybrid CNTs nanoflu- ids resulting from the modifications in M and n can be demonstrated in Figures 4 and 5, respectively.In Figure 4, the value of M is modified from 0 to 0.2, while in Figure 5, the value of n is adjusted from 2 to 4. We noticed that both figures (4 and 5) generate the dual solution when the sheet moves at stretching velocity (ε < 0), where the first solution is thinner than the second solution.Figures 4 and 5 also provide strong evidence that the profiles asymptotically meet the boundary constraints specified in Equation ( 9).The velocity profile f ′ (η) for the first solution has an increasing trend as the values of M and n grow.But the rise in the values of M and n leads to a distinct pattern in the profile of f ′ (η) for the second solution.To illustrate the thermal boundary layer of the first solution, we discover that an increase in M and n results in a drop in θ(η).However, the value of θ(η) for the second solution gets higher as a consequence of the increase in M.   The performance of hybrid CNTs both for the skin friction and heat transfer rate under the impact of M, ϕ 1 = ϕ SWCNT , and ϕ 2 = ϕ MWCNT are discussed in Figure 6.The values of M are set at 0, 0.1 and 0.2, while ϕ 1 and ϕ 2 vary from 0 to 0.1.The values of ϕ 1 and ϕ 2 undergo modifications at low concentrations, as shown by Kumar et al. [48].They believe that the performance of the transfer is effectively enhanced under low-concentration conditions.The other parameters are held constant, with Pr = 6.2, ε = 0.5 and n = 2. Based on this figure, we identified that the rise in M, ϕ 1 and ϕ 2 triggers a boost in both the skin friction and heat transfer coefficients.Prior research conducted by Sun et al. [49] demonstrated the enhancement of heat transfer resulting from the application of a magnetic field.According to Shah et al. [50], they stated that the magnetic field induces an increase in heat transfer via the action of Lorentz forces.With the goal of evaluating the performance of hybrid carbon nanotubes (CNTs) nanofluids with mono-CNTs nanofluids in terms of skin friction and heat transfer, we conduct investigations utilising various values of ϕ 1 , ϕ 2 and M. To choose the most suitable nanoparticles that effectively act as mono-CNTs, we adhere to Samat et al.'s model [51].They found that SWCNTs were better at transferring heat than MWCNTs.This means that SWCNTs represent mono-CNTs.By leaving Pr at 6.2 and n = 2, we can discover from Figures 7 and 8 that the skin friction and heat transfer coefficients for hybrid CNTs nanofluids and mono-CNTs nanofluids both rise.Furthermore, hybrid carbon nanotubes (CNTs) nanofluids exhibit superior performance compared to mono-CNTs nanofluids in terms of both skin friction and heat transfer coefficients.

RSM
In order to compute the numerical values of the responses  and  , we employ a mixture of several parameter levels, including  ,  ,  and .The distribution of various values of these independent variables and their corresponding responses is shown in Table 5. Employing Equation (20) with  = 3,  = 6, we perform 20 trials.

RSM
In order to compute the numerical values of the responses y hCNT and y CNT , we employ a mixture of several parameter levels, including ϕ hn f , ϕ n f , M and n.The distribution of various values of these independent variables and their corresponding responses is shown in Table 5. Employing Equation (20) with F = 3, C = 6, we perform 20 trials.The numerical findings arranged in Table 5 are tested using the analysis of variance (ANOVA) in Minitab.The ANOVA results in analysing the heat transfer coefficients for hybrid CNTs and mono CNTs are displayed in Tables 6 and 7 and Figures 9-14.Both models, y hCNT and y CNT , are considered well-fitting models as their residuals are normally distributed.According to Myres et al. [45], the model is statistically significant if the P-value of the sources is less than 0.05.Since the P-values of the several inputs in Tables 6 and 7 exceed 0.05, we exclude these factors when establishing the reduced models for hybrid CNTs and mono CNTs.The reduced models for hybrid CNTs and mono CNTs are constructed in Tables 8 and 9, respectively.After removing these terms, the ANOVA results in Tables 8 and 9 show that all the sources have p-values less than 0.05.As a consequence, the reduced models fit statistically to produce the optimal solution for heat transfer  and  .Therefore, we have productively formulated the final quadratic regression equations of the heat transfer for hybrid CNTs in Equation ( 23) and mono CNTs in Equation ( 24).The equations are written   After removing these terms, the ANOVA results in Tables 8 and 9 show that all the sources have p-values less than 0.05.As a consequence, the reduced models fit statistically to produce the optimal solution for heat transfer y hCNT and y CNT .Therefore, we have productively formulated the final quadratic regression equations of the heat transfer for hybrid CNTs in Equation ( 23) and mono CNTs in Equation (24).The equations are written as follows: y hCNT = 1.81639 + 0.129914ϕ hn f + 0.001125M + 0.313135n + 0.004701ϕ 2 hn f − 0.027299n 2 + 0.024337ϕ hn f n, ( y CNT = 1.67915 + 0.064819ϕ n f + 0.000988M + 0.289030n + 0.000907ϕ 2 n f − 0.025219n 2 + 0.011281ϕ n f n. The satisfactory outcomes shown in Tables 8 and 9 also enable the determination of the ideal solutions for heat transmission in both kinds of nanofluids.Utilising the response optimizer, the rate of heat transfer for hybrid CNTs and mono CNTs is maximised when ϕ hn f , ϕ n f , M and n are at their highest values.With a composite desirability value of 100%, Tables 10 and 11 reveal that hybrid CNTs have excellent heat transmission performance compared to individual CNTs.In terms of heat transfer, this result shows that the RSM has successfully confirmed the previous computational findings.These findings show that hybrid carbon nanotubes (CNTs) are better than mono CNTs.

Conclusions
The investigation of hybrid CNTs flow and heat transfer past a nonlinear stretching or shrinking sheet is successfully performed using the numerical and RSM approaches under the influence of magnetohydrodynamic effects.The model is developed by taking the magnetic and volume fractions of hybrid and mono CNTs, nonlinear, and shrinking or stretching parameters as the leading variables in non-dimensionless units.The similarity equations in the form of non-dimensional ordinary differential equations are solved using the similarity variables method.The bvp4C MATLAB function is employed to numerically compute the first order of the ODEs system.By referring to the mathematical model and using the numerical experiment data, the design of the heat transfer coefficients both for hybrid CNTs and mono CNTs is statistically optimised in Minitab.The design of this optimisation procedure involves three parameters and twenty simulation runs.The findings show that: • The multiple solution is generated when the sheet moves nonlinearly in the shrinking region where ε c ≤ ε < 0. • The boundary layer of the first solution is thinner than that of the second solution.

•
The model may provide a couple of solutions when the magnetic parameter is below 1. • The magnetic and non-linear parameters delay the occurrence of boundary layer separation and expand the range of potential solutions for the reduced skin friction f ′′ (0) and the reduced heat transfer −θ ′ (0).• The magnetic and hybrid CNTs volume fraction parameters enhance the favourable effect on skin friction and heat transfer coefficients.• To maximise the heat transfer rate, the magnetic, nonlinear, and hybrid CNTs volume fractions are set at their greatest values.

•
Hybrid CNTs are superior to mono CNTs for the skin friction coefficient.• The numerical analysis and response surface methodology demonstrate that hybrid CNTs are better than mono CNTs with regard to the heat transfer rate.

Future Works
Due to some limitations in this model, there are several recommendations to extend the model in the future.They are as follows: • The investigation of hybrid nanofluids flow-based carbon nanotubes over a nonlinear stretching/shrinking sheet under an unsteady flow case.• The study of hybrid carbon nanotube flow across a nonlinear stretching/shrinking sheet using the 3-Dimension (3D) flow approach.• The examination of hybrid nanofluids flow-based carbon nanotubes over a nonlinear stretching/shrinking sheet with distinct effects, including velocity and thermal slip effects.• The exploration of hybrid carbon nanotube flow past different geometric configurations, such as a Riga plate.
The suggestions mentioned above for future research could potentially improve the comprehension of the fluid motion and heat transfer characteristics of hybrid nanofluids, particularly for hybrid carbon nanotubes.

Figure 1 .
Figure 1.The schematic model of stagnation point flow over nonlinear (a) stretching sheet and (b) shrinking sheet.

Figure 1 .
Figure 1.The schematic model of stagnation point flow over nonlinear (a) stretching sheet and (b) shrinking sheet.

Figure 4 .
Figure 4. Variation of solutions (a) velocity; (b) temperature profiles using different values of  in hybrid SWCNTs-MWCNTs water-based nanofluids.

Figure 6 .Figure 6 .
Figure 6.Variation of solutions (a) skin friction; (b) heat transfer coefficients using different values of  and  in hybrid SWCNTs-MWCNTs water-based nanofluids.

Figure 7 .Figure 7 .Figure 8 .
Figure 7.Comparison performance of hybrid CNTs and mono CNTs in (a) skin friction; (b) heat transfer coefficients using different values of  =  and  =  .

2 Figure 8 .
Figure 8.Comparison performance of hybrid CNTs and mono CNTs in (a) skin friction; (b) heat transfer coefficients using different values of M.

Figure 13 .
Figure 13.Contour plots for the interaction of (a) x 1 and x 2 ; (b) x 1 and x 3 ; (c) x 2 and x 3 for mono CNTs.

Figure 14 .
Figure 14.Surface plots for the interaction of (a)  and  ; (b)  and  ; (c)  and  for mono CNTs.

Figure 14 .
Figure 14.Surface plots for the interaction of (a) x 1 and x 2 ; (b) x 1 and x 3 ; (c) x 2 and x 3 for mono CNTs.

Table 2 .
The thermophysical characteristics of hybrid CNTs and water at room temperature.

Table 3 .
Different levels of values of independent variables (factors).
Computation 2024, 12, x FOR PEER REVIEW 11 of 32 conduct investigations utilising various values of  ,  and  .To choose the most suitable nanoparticles that effectively act as mono-CNTs, we adhere to Samat et al.'s model [51].They found that SWCNTs were better at transferring heat than MWCNTs.This means that SWCNTs represent mono-CNTs.By leaving  at 6.2 and  = 2 , we can discover from Figures 7 and 8 that the skin friction and heat transfer coefficients for hybrid CNTs nanofluids and mono-CNTs nanofluids both rise.Furthermore, hybrid carbon nanotubes (CNTs) nanofluids exhibit superior performance compared to mono-CNTs nanofluids in terms of both skin friction and heat transfer coefficients.

Table 5 .
Design of experiment of heat transfer in hybrid CNTs and mono CNTs with waterbased nanofluids.

Table 6 .
Analysis of variance for the full model of heat transfer in hybrid CNTs with waterbased nanofluids.
Note: The symbol * indicates that the data is too small.

Table 7 .
Analysis of variance for the full model of heat transfer in mono CNTs with waterbased nanofluids.The symbol * indicates that the data is too small.: The symbol * indicates that the data is too small. Note:Note

Table 7 .
Analysis of variance for the full model of heat transfer in mono CNTs with water-based nanofluids.
Note: The symbol * indicates that the data is too small.

Table 8 .
Analysis of variance for the reduced model of heat transfer in hybrid CNTs with waterbased nanofluids.
Note: The symbol * indicates that the data is too small.

Table 9 .
Analysis of variance for the reduced model of heat transfer in mono CNTs with water-based nanofluids.
Note: The symbol * indicates that the data is too small.

Table 10 .
The optimal solution of heat transfer in hybrid CNTs with water-based nanofluids.

Table 11 .
The optimal solution of heat transfer in mono CNTs with water-based nanofluids.