Application of ANFIS and LSSVM strategies for estimating thermal conductivity enhancement of metal and metal oxide based nanofluids

An extensive variety of chemical engineering processes include the transfer of heat energy. Since increasing the effective contact surface is known as one of the popular manners to improve the efficiency of heat transfer, the attention to the nanofluids has been attracted. Due to the difficulty and high cost of an experimental study, researchers have been attracted to fast computational methods. In this work, Adaptive neuro-fuzzy inference system and least square support vector machine algorithms have been applied as a comprehensive predictive tool to forecast the nanofluids thermal conductivity in terms of diameter, temperature, the thermal conductivity of the base fluid, the thermal conductivity of nanoparticle and volume fraction. To this end, a large and comprehensive experimental databank contains 1109 data points have been collected from reliable sources. The particle swarm optimization is utilized to reach the best structures of the proposed algorithms. A comprehensive statistical and graphical investigations are carried out to prove the accuracy and ability of proposed models. In addition, the comparisons outputs indicate that the least square support vector machine algorithm has the best performance among the existing correlations and Adaptive neuro-fuzzy inference system algorithms for forecasting thermal conductivity of different nanofluids.


Introduction
Owning to extensive applications of heat transfer phenomenon in different industrial instruments such as the cooling systems, compressor, evaporators, and boilers, it is well-known as a critical point. There is a direct relationship between the efficiency of such instruments and their requirement volume, which is significant for us to have economically process (Wen, Lin, Vafaei, & Zhang, 2009;Zalba, Marı´n, Cabeza, & Mehling, 2003). Furthermore, it should be noted that a kind of working fluid usually passes within such heat transfer devices through a pump or compressor and as the efficiency of heat transfer increases the required power consumption and thus processing cost decrease. Numerous studies have been implemented in order to improve heat transfer efficiency (Ebrahimnia-Bajestan, Moghadam, Niazmand, Daungthongsuk, & Wongwises, 2016;Mazumdar et al., 2016;Tomlinson, Manning, Schaefer, & Record, 2016). The improvement of effective surface, using microscale channels and vibration technique can be named as efforts that help to enhance the heat transfer efficiency. According to annually published articles, there are numerous attentions on nanofluid systems. Furthermore, the thermal conductivity property is known as a major factor in the improvement of heat transfer efficiency. the thermal conductivity property of some base fluids such as ethylene glycol (EG), water (H 2 O) and different types of oil, some addition of nanoparticles to these fluids can increase their thermal conductivity (Eiamsa-ard & Wongcharee, 2018;Zeeshan, Shehzad, Ellahi, & Alamri, 2017). These nanoparticles give us the advantage of having a simple fluidized process and they solve some major problems, for example, the precipitation of particles, the blockage of channels, and erosion problems because of their nanosize structures.
The suspension of particles in a base fluid was introduced by Chol as an innovative aspect of nanofluilds (Chol, 1995). Furthermore, an interesting application of them was suggested in microchannels (Chein & Chuang, 2007). In recent years, numerous attempts have been done to estimate the nanofluids thermal conductivity utilizing mathematical and computational methods. Ghalambaz et al. investigated heat transfer convection of Cu-water nanofluids through a porous media parallelogrammic enclosure numerically by utilizion Tiwari model. They studied the effect of volume fraction of nanoparticles, inclination angle and porous media matrix on heat transfer (Ghalambaz, Sheremet, & Pop, 2015). Sheremet et al. used Buongiorno' model to study free convection of a porous media contained a nanofluid (Sheremet, Groşan, & Pop, 2015). Their analysis was carried out for various parameters such as Brownian motion parameter, thermophoresis parameter, Lewis numbers, and Rayleigh number in the range using a finite difference method. Maxwell in 1904 proposed a simple correlation to estimate the effective thermal conductivity (Maxwell, 1881). This approach performed in an acceptable manner for forecasting thermal conductivity of the fluids containing microparticles, while the mentioned approach hasn't respectable forecasting thermal conductivity of nanofluids and the considerable amount of deviations are determined between its outcomes and actual values (Keblinski, Phillpot, Choi, & Eastman, 2002;Kleinstreuer & Feng, 2011;Lee, Lee, Choi, Jang, & Choi, 2011;Murshed, Leong, & Yang, 2008). The several experimentally studies regarding the determination of the nanofluids thermal conductivity and investigation of major parameters which have an impact on nanofluids thermal conductivity. The effective parameters are the nature of the particles and of base fluid, temperature, nanoparticles volume fraction. Furthermore, considerable efforts have been carried out to suggest a novel method for thermal conductivity improvement based on theoretical backgrounds. Vatani et al. carried out a study to forecast effective thermal conductivity of nanofluids based on a comprehensive databank collected from recently published papers and also they carried out a comparison between the results of their proposed model and seven other published correlations (Vatani, Woodfield, & Dao, 2015). It was concluded that all correlations can't be useful for different nanofluids with a high degree of accuracy, but these are known as popular methods regarding the prediction of nanofluid thermal conductivity. Due to this fact, the prediction of nanofluid thermal conductivity with high accuracy is a very stimulating subject and critical issue. In addition, experimental studies are time-consuming and costly to determine nanofluids thermal conductivity and also they can't investigate various working conditions (Ariana, Vaferi, & Karimi, 2015).
Moreover, the suggested existence methods can cover limited ranges of operational conditions. Hence applying an exact method for the identification of the relationships between related parameters to predict the thermal conductivity is a regularly problematic issue. To that end, utilization of computational methods and artificial intelligence methods such as the support vector machines, fuzzy inference systems and artificial neural networks (ANNs) which usually give us high accurate results, has been recommended by numerous researchers in different issues (Ahmadi, Ghahremannezhad, et al., 2019;Ahmadi, Mohseni-Gharyehsafa, et al., 2019;Ahmadi, Sadeghzadeh, Raffiee, & Chau, 2019;Ali Ghorbani, Kazempour, Chau, Shamshirband, & Taherei Ghazvinei, 2018;Baghban, Ahmadi, & Shahraki, 2015;Baghban, Ahmadi, Pouladi, & Amanna, 2015;Baghban, Bahadori, Mohammadi, & Behbahaninia, 2017;Baghban, Mohammadi, & Taleghani, 2017;Baghban et al., 2016;Bahadori et al., 2016;Chau, 2017;Cheng & Chau, 2002;Moazenzadeh, Mohammadi, Shamshirband, & Chau, 2018;Wu & Chau, 2011;Yaseen, Sulaiman, Deo, & Chau, 2019). ANN is a great satisfactory approach to reach optimal solutions for difficult problems, especially in chemical fields and can help us to reduce the time and cost considerably (Abdi-Khanghah, Bemani, Naserzadeh, & Zhang, 2018). Thanks to rapid process and broad applications of the ANN, numerous works have been done to predict the different properties of nanofluids such as viscosity, density and thermal conductivity (Atashrouz, Pazuki, & Alimoradi, 2014;Karimi & Yousefi, 2012;Karimi, Yousefi, & Rahimi, 2011;Sharifpur, Adio, & Meyer, 2015). Hojjat et al. used an artificial neural network model to estimate of non-Newtonian nanofluids thermal conductivity (Hojjat, Etemad, Bagheri, & Thibault, 2011). The input variables of the mentioned model are operational temperature, the concentration of nanoparticle, and the thermal conductivity of nanoparticle. moreover, another model as diffusion neural network algorithm was suggested by Papari et al. for predicting the thermal conductivity of some nanofluids containing multi-walled carbon nanotubes in four different base fluid such as distilled water, decene, oil, and ethylene glycol (Papari, Yousefi, Moghadasi, Karimi, & Campo, 2011). Longo et al. implemented two structures of ANN algorithms to forecast thermal conductivity of two nanoparticles of Al 2 O 3 and TiO 2 in water base fluid in terms of the average size, thermal conductivity and volume fraction of nanoparticle and temperature (Longo, Zilio, Ceseracciu, & Reggiani, 2012). An ANN-based method was proposed by Hemmati Esfe to forecast the thermal conductivity of Al 2 O 3 -Water nanofluid at various solid volume fractions and temperatures (Esfe, Afrand, Yan, & Akbari, 2015). In the recent decade, another types of artificial intelligence methods namely least square support vector machine algorithm and Adaptive neuro-fuzzy inference system have attracted the interests of researchers in different industries. These methods can be accurate predictors for different processes in chemical engineering (Baghban, Bahadori, et al., 2017;Baghban et al., 2016). The mentioned approach has been less used for prediction of nanofluids properties (Meybodi, Naseri, Shokrollahi, & Daryasafar, 2015).
Lack of enough experimental data, difficulty, time consuming and high cost of experimental studies and on the other hand limitations of all published approaches to number of nanofluids and low precision of them in prediction of thermal conductivity of nanofluids encourage the authors to do this study for overcoming these problems by using a general method for prediction of the thermal conductivity of nanofluids with satisfactory outputs.

Objective of the study
The estimation of thermal conductivity of nanofluids is one of the important topics in heat transfer researches, so the proposing of an accurate and lowcost method for estimation of this parameter based on existing parameters will be a valuable research. In the present study, two artificial intelligence approaches based on least square support vector machine (LSSVM) and Adaptive neuro-fuzzy inference system (ANFIS) algorithms are proposed to predict the thermal conductivity of nanofluids in terms of nanoparticle volume fraction, diameter, temperature, the thermal conductivity of nanoparticle and base fluid. These input parameters are chosen based on the study of effective parameters in different existing correlations Also, our models are compared with some existing correlation based methods to identify an accurate tool for the estimation of thermal conductivity.
There are some operative medium models in order to estimate the effective thermal conductivity. The computational approaches such as Maxwell approach (Maxwell, 1873), Hamilton approach (Hamilton & Crosser, 1962) and Bruggeman approach (Bruggeman, 1935) estimated the effective thermal conductivity based on the thermal conductivity of particle, the volume fraction, and thermal conductivity of the base fluid. Furthermore some recent proposed and enhanced approaches are discussed here. The effect of the nanolayer were considered by these newly methods same as Yu-Choi's approach (Yu & Choi, 2003), Leong's approach (Leong, Yang, & Murshed, 2006), Xie's approach (Xie, Fujii, & Zhang, 2005) and Sohrabi's approach (Sohrabi, Masoumi, Behzadmehr, & Sarvari, 2010) and also some of them such as the Sohrabi's approach (Sohrabi et al., 2010), Koo's approach (Koo & Kleinstreuer, 2004) and Xu's approach (Xu, Yu, Zou, & Xu, 2006) considered the impact of the convective heat transfer comes from Brownian motion. Also, a study was done by Evans (Evans et al., 2008) in order to show that clustering and interfacial thermal resistance have an effect on the thermal conductivity of nanofluids. The literature papers were reported broad conducted studies by different scientists and engineers to estimate nanofluids thermal conductivity.

Overview of LSSVM strategy
However methods raised artificial neural network-based have good features like high performance, these models have some faults like they're not reproducible results due to alteration of optimization processes and unsystematic initialization of this form of algorithms. The support vector machine (SVM) approach utilizes a range of input parameters from non-linear functions to multidimensional mappings. Linear decision surface is used to connect the Input and output spaces. The SVM-based approaches require less adjustable parameters respect to ANN-based approaches; furthermore, SVM-based approaches don't require adjusting the number of hidden layers or their neurons so the mentioned methods would give more accurate generalization (Suykens & Vandewalle, 2000).
LSSVM was made by Suykens in 1999 which was more basic than SVM (Suykens & Vandewalle, 1999). In this method, the support vectors and linear equations have been utilized to make the quadratic programing systems simple and decrease the complexity of optimization. These types of models use the difference of determined parameters by the approaches and the measured data to obtain the regression error, however, in the SVM approach, this error reaches optimal value by calculation. the optimization process of LSSVM approach is done as below (Guo & Bai, 2009;Suykens, Van Gestel, & De Brabanter, 2002): Used in below equation where T is transposed matrix, w is weight vector, μ is a relation for total regression weight errors and single weights, g(x) denotes a mapping function, e i represents regression error and b is the bias terms. The regression weight coefficient (W) can be calculated as a function of an input vector (x i ) and Lagrangian multiplier (α i ) as The aforementioned equation can be altered by assuming the linear relation for dependent and independent parameters as below: Lagrange multiplier (α i ) will be obtained by using manipulation as The kernel function is constructed based on a modification of Equation (5) and making it applicable for nonlinear restrains so it concluded to: Where Kernel function can be defined as The usual form of Kernels which are ever applied in the evolvement of the LSSVM approach are Gaussian Radial Basis function Kernels and a sample of them is as follow: where σ 2 represents squared bandwidth and it is necessary to be determined by using an optimization method in the training phase.

Overview of ANFIS strategy
Jang suggested Adaptive neuro-fuzzy inference system or ANFIS which generates membership functions and fuzzy rule base (Jang, Sun, & Mizutani, 1997). The ANFIS algorithm consists of some connected nodes which affected by learning rules to reduce error. To reach the best structure of ANFIS network, different evolutionary approaches can be utilized (Afshar, Gholami, & Asoodeh, 2014;Baghban et al., 2016;Karkevandi-Talkhooncheh et al., 2017). In the present work, particle swarm optimization is utilized. In the purpose of better introducing ANFIS, the brief summary of ANFIS structure is shown in Figure 1. As demonstrated in this figure, there are two inputs and one output. The first layer of the network belongs to the linguistic terms come from input parameters. These linguistic terms can be organized by Gaussian membership function. The below formulation expresses the Gaussian function (Ahangari, Moeinossadat, & Behnia, 2015): here, Z and σ indicate the Gaussian parameters. The second layer belongs to the determination of weight terms based on rules: The third layer represents the process of determination of weight average values as following: In the fourth layer, a function in terms of average weight values is constructed as follows: Here, m, n, and r represent the resulting indexes.
In the last layer, the summation of previous layer outputs is determined such as the following:

Overview of PSO algorithm
Particle swarm optimization which was created based on existing populations in nature such as birds is known as one of the popular stochastic optimization methods. In this optimization algorithm, the promotion of initial population is optimization process; particles are known as the solutions and a swarm is a group of particles (Castillo, 2012;Eberhart & Kennedy, 1995;Onwunalu & Durlofsky, 2010;Panigrahi, Shi, & Lim, 2011). In this algorithm, the particle velocity V i (t) is updated by the following equation: p best,id and g best,id are known as the best position of particle and global position c, r and w denote the learning rate,random number and inertia weight. Also the position of particle can be determined as following (Chiou, Tsai, & Liu, 2012;Eberhart & Kennedy, 1995;Shi & Eberhart, 1998):

Data preparation
Efficiency and consistency of any developed algorithm is extremely function of utilized real data points which need to be precise and easily accessible in order to carry out the algorithm ( Table S1. On the other hand, Table S2 summarized the thermal conductivity of particle and base fluid. The second step after preparation of data is choosing the input and output parameters for the estimating algorithms. The input parameters of proposed LSSVM and ANFIS models were the base fluid thermal conductivity (W/mK), the nanoparticle thermal conductivity (W/mK), temperature (K), void fraction (v) and particle diameter (nm). Also for the outputs of algorithms, the nanofluids thermal conductivity was assumed. The collected dataset has been divided into the testing and training data subsets to prepare and evaluate the algorithms. A quarter of total dataset was selected randomly as a testing dataset and other 75% of data points was applied for the training of the LSSVM and ANFIS approaches.

Model development and evaluation
Due to the fewer number of tuning terms and highspeed computations of radial basis function (RBF) kernel respect to other functions, the aforementioned function was selected as a useful kernel function in compliance with published papers. In the LSSVM approach combined with radial basis function, a remarkable problem is to find two terms (i.e. γ and σ 2 ) as mentioned before. These two tuning parameters have a dominant application in order to reach a high-performance LSSVM approach with an acceptable degree of globalizations and approximations, where σ 2 refers to the kernel sample variance γ expresses the regularization term. The current paper applies the particle swarm optimization (PSO) approach for the calculation of optimal values of γ and σ 2 . The fundamental of the optimization approach is to obtain less value of a cost function (i.e. mean absolute relative error for the testing phase in this study). The optimization process has been continued some iterations as efforts to reach an accesable global best condition based on the presented fitness function. The suggested ANFIS model includes 10 clusters and  120 membership function parameters, which were determined by the optimization algorithm. Figures 2 and 3 illustrates a diagrams of PSO-LSSVM and PSO-ANFIS models used in the present paper. Moreover, the characterizations of proposed models have been presented in Table 1. This table contains the number of utilized data in different steps of proposing algorithms, tuning parameters, number of iteration to reach the criteria and kernel function type.
The approximations by the suggested LSSVM and ANFIS models were assessed by using different methods such as the graphical and statistical confirmations. The cross diagrams of suggested models which are made through depicting predicted against measured values were shown in Figures 4 and 5 for LSSVM and ANFIS, respectively. Based on observations obtained from this figure, a satisfactory agreement exists between estimated and measured nanofluid thermal conductivity values for ANFIS and LSSVM algorithms due to the close distribution of data points around the identity line. As it is obvious in this figure the predicted and real thermal conductivity of nanofluid cover each other in an appropriate manner for both proposed models. Figure 6 illustrates the relative errors of predicted thermal conductivity of nanofluid and measured points for both models. As shown, the horizontal axis representes the measured thermal conductivity of nanofluid and the vertical axis denotes the absolute error percent. The outcomes of this figure express that the obtained error values by the LSSVM and ANFIS algorithms range mostly of 5 to −5 percentage which proves the remarkable performance capacity of suggested models. In order to indicate  deviations in a better way, relative deviations of two developed models are shown in Figure 7 versus dimensionless thermal conductivity which is a ratio of thermal conductivity of mixture and base fluid. As shown, both models have relatively similar deviations. Beside these graphical assessments, some popular statistical approaches were applied base on the below equations in order to confirm the satisfactory utilization of our suggested LSSVM and ANFIS models.
• The coefficient of determination (R 2 )  • average absolute relative deviation • Root Mean squared error (RMSE) • Standard deviations (STD) The abovementioned statistical analyses were summarized in Table 2 for the LSSVM and ANFIS models at three phases of train, test and overall.
The present table indicates that the LSSVM and ANFIS models display reasonably low values of RMSE, AARD, and STD and high values of R-squared. The  determined R 2 for LSSVM and ANFIS are 0.994 and 0.984 respectively. Also, the RMSE values for LSSVM and ANFIS were determined as 0.02 and 0.034 respectively. These analyses confirm the consistency and applicability of suggested models in order to predict the thermal conductivity of nanofluids. Also, these analyses show that the LSSVM algorithm has better performance respect to the ANFIS algorithm.
Lots of comparisons were implemented between the results of suggested ANFIS and LSSVM models and other published correlations. For Al 2 O 3 -water nanofluid has d p = 38.4 nm at temperature of 324 K, the predicted dimensionless thermal conductivities by Mintsa, Nan, Murshed, Yu-Choi, Xie, and Chon et al. correlations have been compared with suggested LSSVM and ANFIS models in Figure 8 (Chon et al., 2005;Mintsa et al., 2009;Murshed et al., 2008;Nan, Birringer, Clarke, & Gleiter, 1997;Xie et al., 2005;Yu & Choi, 2003). As can be seen, as the nanoparticle volume fraction increases the dimensionless thermal conductivity also increases.  Four other correlations such as the Godson, Wasp, Timofeeva, and Hamilton Correlations have been applied in order to compare with suggested LSSVM and ANFIS models for predicting Ag-water nanofluid dimensionless thermal conductivity with d p = 63 nm and temperature of 343 K (Godson, Raja et al., 2010;Hamilton & Crosser, 1962;Timofeeva, Moravek, & Singh, 2011;Wasp, Kenny, & Gandhi, 1977). This comparison was shown in Figure 9 and results obtained from the present model have great accuracy with corresponding experimental dimensionless thermal conductivity. Thang correlation has been compared with our suggested models to predict carbon nanotube-water nanofluid dimensionless thermal conductivity as shown in Figure 10 (Thang et al., 2015). Furthermore, another comparison of the mentioned correlation and predicting models at a temperature of 296.15 K for carbon nanotube-water nanofluid has been indicated in Figure 11. These above two figure can indicate a greatly satisfactory prediction of LSSVM model than Thang et al. correlation. Moreover, another comparison was carried out in Figure 12 between the outcomes of suggested LSSVM and ANFIS models and four other correlations namely the Corcione, Vajjha, Patel and Azmi correlations for determining the CuO-water nanofluid dimensionless thermal conductivity that has a particle diameter of 24 nm at T = 298.15 K (Azmi, Sharma, Mamat, Alias, & Misnon, 2012;Corcione, 2011;Vajjha & Das, 2009). Moreover, for TiO 2 -EG nanofluid at 298.15 K, the predicted dimensionless thermal conductivity by Jang's correlation and ANFIS and LSSVM models at different nanoparticle size and volume void fraction has been illustrated in Figure 13 (Jang & Choi, 2004).
As can be concluded, artificial intelligence based approaches can satisfactory applications in complicated cases such as nanofluid systems. Both developed models have great accuracy in order to estimate the thermal conductivity of nanofluids compared to other previously developed models. It is worth mentioning that the great advantage of these techniques is their simple application (as you can see the supplementary content) beside their accuracy.
The experimental data which are used in the comparisons in Figures 8-13 can be identified with blue color in the cross plot. As can be observed from all the above comparisons, the suggested LSSVM have the satisfactory predictive capacity to forecast the dimensionless thermal conductivity of nanofluids. The residual plot which shows a difference in residual variance versus experimental thermal conductivity was shown in Figure 14 for both models. Heteroscedasticity produces a distinctive fan or cone shape in residual plots. As can be seen, this pattern has a cone-shape and it is shown the errors are normally distributed.

Outlier detection and sensitivity analysis
The measured data applied in the computational investigation can highly influence the predictive ability of our proposed algorithms (Rousseeuw & Leroy, 2005). As there is a various source of measurement errors of actual datasets in the literature, identification of outlier and suspected data for utilized experimental data points can indicate greatly bad experimental data points.
William's plot as a graphical illustration can help us to identify accurately and suspected data points. This figure has been depicted in Figure 15 for LSSVM and ANFIS algorithm models to detect outliers.
Besides, the leverage boundary is defined base on the green line and the experimental data which their hat values (H) are more than the criterion value H* are known as outliers. H and H* can be determined as following: where X represents as the m×n matrix that n and m denote the model indexes number and number of samples (Mohammadi, Gharagheizi, Eslamimanesh, & Richon, 2012). The standardized residual values have the values of −3 and + 3 are shown by red color is shown in Figure 15 so the accurate experimental data are also located inside of the defined boundaries. One of the popular statistical approaches is Sensitivity analysis to illustrate the relationship between output and inputs. The impact of input parameters on the target parameter is identified by the relevancy factor (r) and in the current work, we are going find out the effect inputs on thermal conductivity of nanofluids. The relevancy factor can be determined as below: where Y i ,Ȳ,X k,i and X k are the 'i' th output, output average, k th of input and average of input (Hosseinzadeh & Hemmati-Sarapardeh, 2014;Razavi, Bemani, Baghban, Mohammadi, & Habibzadeh, 2019). The aforementioned index is varied between −1 and + 1 and as this term has the higher value is for a specific input variable, it can be understood that this input variable has much influence on the output. A negative value indicates a decreasing effect and Positive value has an increasing effect (Hosseinzadeh & Hemmati-Sarapardeh, 2014).
As it is obvious in Figure 16, the volume void fraction, thermal the conductivity of the base fluid and temperature has a straight relationship with nanofluid thermal conductivity values. Due to determined r values for different parameters, the particle diameter is the weakest parameter on target parameter determination and the most effective parameter is thermal conductivity of the base fluid.

Conclusions
In the current paper, different statistical and graphical analyses are carried out to find a high-performance method for estimation of thermal conductivity of nanofluids. This investigation was concluded to the following results: • In order to show the acceptable performance of LSSVM and ANFIS algorithms, a comprehensive dataset from existing experimental studies in literature has been gathered. • The LSSVM and ANFIS algorithms reach the best performance structure by the employment of the PSO algorithm. Furthermore, a high degree of accuracy was observed in the prediction of the thermal conductivity of different nanofluids by LSSVM and ANFIS algorithms. • In the purpose of illustrating the predictive ability of proposed models, the LSSVM and ANFIS algorithms have been compared with 15 different correlation based methods for estimation of nanofluids thermal conductivity. The comparisons expressed the amazing performance of LSSVM algorithm. • The current study has proposed two accurate methods which can be used as simple and user friend software in different projects. • In the current work, a great mathematical analysis has been used to improve the accuracy of predictive models by separating suspected experimental data. • A sensitivity analysis has been used to show the impact of different parameters on the thermal conductivity of nanofluids. This analysis expresses that thermal conductivity of base fluid is the most effective parameter in the determination of thermal conductivity of nanofluids. Also, the particle diameter is the least effective parameter in this process.
According to the aforementioned conclusions, the current work can be a comprehensive and helpful study about the heat transfer process in nanofluids which makes this work valuable for engineers and researchers. As a future suggestion, the current results can be compared with other modeling approaches such as numerical modeling and artificial intelligence based approaches (Genetic programing, committee machine intelligence system, and artificial neural network).

Disclosure statement
No potential conflict of interest was reported by the authors.