An Accurate and Computationally Efficient Explicit Friction Factor Model

The implicit Colebrook equation has been the standard for estimating pipe friction factor in a fully developed turbulent regime. Several alternative explicit models to the Colebrook equation have been proposed. To date, most of the accurate explicit models have been those with three logarithmic functions, but they require more computational time than the Colebrook equation. In this study, a new explicit non-linear regression model which has only two logarithmic functions is developed. The new model, when compared with the existing extremely accurate models, gives rise to the least average and maximum relative errors of 0.0025% and 0.0664%, respectively. Moreover, it requires far less computational time than the Colebrook equation. It is therefore concluded that the new explicit model provides a good trade-off between accuracy and relative computational efficiency for pipe friction factor estimation in the fully developed turbulent flow regime.


Introduction
Friction factor estimation is important for modeling flows in pipes and is relevant in most engineering disciplines, for example: chemical, civil and mechanical.Over the years, the Colebrook equation [1] [2] has been widely used for pipe friction factor estimation in the fully developed turbulent regime.The equation is expressed as: 1 2.51 2 log 3.71 Re The major drawback of Colebrook equation is that it is implicit in friction factor (f). Therefore, it requires it-eration to obtain its solution.For simulations of long pipes and network of pipes, the Colebrook equation must be solved a huge number of times [3].Therefore, an iterative solution to the Colebrook equation will be time consuming.The use of the Moody chart [4], as an alternative to the Colebrook equation, eliminates the requirement for iteration.However, it is a graphical tool and therefore not convenient for computer-based simulations.
The quest for a fast, non-iterative and accurate model, as an alternative to the Colebrook equation, has given rise to various explicit friction factor models.These explicit models differ in their accuracies and relative computational efficiencies, depending on their degree of complexity.
In this work, a new explicit model was developed for estimating friction factor in the range for which the Colebrook equation is valid.The trade-off between model accuracy and relative computational efficiency has been considered.
The remaining sections of this paper are organized as follows: Section 2 reviews the available explicit friction factor models based on accuracy, complexity and relative computational efficiency.In Section 3, the development of the proposed model is presented while Section 4 reports the performance of the proposed model in comparison with those of the selected existing explicit models.In the final section, relevant conclusions are drawn based on the results obtained in this study.

Accuracy
The accuracies of the existing explicit models have been reported using common criteria such as the mean square error (MSE), percentage relative error and absolute error [5]- [8].Model selection criteria (MSC) and Akaike information criterion (AIC) were used by Romeo, Royo and Monzon [9] for explicit model selection.These criteria were subsequently used by Genić et al. [7] and Yildrim [8] for comparison of several explicit models.Unfortunately, there is an apparent discrepancy in the MSC values reported [7] [9] for the same models.For example, the MSC values reported by Romeo, Royo and Monzon [9] and Genić et al. [7] for Moody [10] and Chen [11] models showed a wide contrast.
It has been shown that models with greater number of logarithmic functions are generally more accurate than those with lesser number of logarithmic functions, although the former require more computational time than the latter [6].For instance, it is observed from works of Brkić [12], Winning and Coole [6] and Fang, Xu and Zhou [13], that the most accurate approximations are those by Zigrang and Sylvester [14], Serghides [15], Romeo, Royo and Monzon [9] and Buzzelli [16].These models, with the exception of the model by Buzzelli [16], have three logarithmic functions (either natural logarithm or logarithm to base ten).
Brkić [12], based on maximum relative error criterion, classified the existing explicit models as extremely accurate (error ≤ 0.14%), very accurate (error up to 0.5%), moderately accurate (error up to 1.5%), less accurate (error up to 5%), non advisable (error up to 25%) and extremely inaccurate (error ≥ 80%).Based on this classification, the performances of several explicit models were evaluated and their accuracies are summarized in Table 1.Yildrim [8] conducted a comparative review of 16 explicit models.In his work, friction factor data were generated by digitizing the Moody chart.The turbulent portion of the Moody chart is a graphical solution of the Colebrook equation.Hence, digitizing the Moody chart [4] may have introduced secondary errors in the overall analysis [12].This view is supported by the error margin observed by Fang, Xu and Zhou [13].Ghanbari, Farshad and Rieke [17] also digitized the Moody chart [4] when developing their model.They claim that the model is valid for Reynolds number (Re) between 2100 ≤ Re ≤ 10 8 .It is not obvious how data was obtained for Reynolds number between 2100 and 3000 (critical zone), since the Moody chart does not contain Re values in this range.

Model Complexity and Computational Efficiency
In the bid to develop accurate explicit models, sometimes, simplicity is sacrificed for accuracy without consideration for the actual computational power of such model for massive numerical requirements [12] [18].Generally, accuracy is obtained at the expense of additional computational complexity.To strike a balance between these two properties (complexity and accuracy), Zigrang and Sylvester [5] introduced the concept of complexity using friction factor models.Based on this concept, Brkić [12] computed the complexity and complexity index for 25 explicit models.His work revealed that models with three logarithmic functions or internal iterations were ( ) Not specified Extremely accurate Chen [11] 1.1098 ; where ( ) ( ) Extremely accurate Barr [19] 0.7 0 10 .52 1 e 1 Re 29 Less accurate 13 Churchill [24] ( ) ; where 16 0.9 Less accurate Continued 14 Round [25] ( ) Non-advisable 15 Brkić [26] 0.4343 Not specified Less accurate 16 Rao and Kumar [27] ( ) where 2 Re 0.33 ln 6.5 Swamee and Jain [28] 1.11 10 0.9 where ( ) ( ) Extremely accurate more complex than the models which have two and less internal iterations Winning and Coole [6] carried out a comparative review of 28 explicit friction factor models.They defined relative computational efficiency as the time taken by an explicit model to perform a task relative to the time taken by the Colebrook equation.The use of computational efficiency in their work clearly showed the impact of model complexity on the simulation time.They found that the models developed by Buzzelli [16] and Serghides [15] were the most accurate when ordered by absolute and relative errors, but when ordered by relative computational efficiencies, they ranked very low.The overall ranking reported was biased since it is not based on actual values of accuracy and relative computational efficiency.It was based on the number of available explicit models.If this number is altered, the values of the combined ranking may change.
Computational efficiency is observed to be dependent on the type of logarithmic function(s) contained in the reported models.The computation of the logarithm function in many computer languages is based on series expansion that requires several powers of arguments to be computed and added to each other [18].Glustolisi [18] and co-worker state that the natural logarithm function executes faster than the logarithmic function to base ten.This is based on the fact that the convergence function used for its computation is quite fast.Therefore, the computation of the logarithm function to base ten in many computer languages is based on the computation of the natural logarithm [18].It should be noted that an explicit equation which requires computational time longer than that of the Colebrook's equation defeats the aim of its development.An ideal explicit model should give a good trade-off between its accuracy and relative computational efficiency.

Data Generation
Using Microsoft Excel spread sheet, friction factor (f) data within an error limit of 10 −9 were obtained from Equation (1) for Re values in the range 4 × 10 3 ≤ Re ≤ 10 8 , using 1000 intervals in geometric order and (ε/D) value ranging from 10 −6 to 0.05 using 28 intervals in arithmetic order.Thus, producing a matrix of 28,000 datasets for f, Re and (ε/D) was obtained for model.

Model Development
The plot of Re f against Re for the different relative roughness gives straight lines, as shown in Figure 1.Therefore, the model presented in this work is derived from the implicit Colebrook equation by substitution of the parameter, Re f with the equation of a straight line having Re as the abscissa.Then natural logarithm function was introduced to enhance the computational efficiency of the model as noted by Glustolisi [18].After some rearrangements, the proposed new model was thus obtained as: Using surface-fitting function in the MATLAB curve-fitting toolbox, coefficients a, b, c, d, e, g, and h with their parameter bounds were obtained at 95% confidence level (Table 2).The uncertainties associated with the estimated parameters, which are a measure of the reliability of the parameters, and consequently, a measure of the adequacy of the model, are reported in Table 2.A model which has parameter estimates with low levels of uncertainties (narrow intervals) is deemed to be good and adequate [30].

Colebrook
Relative error 100 ( )   5) Relative Computational efficiency: According to Winning and Coole [6], relative computational efficiency is the ratio of the time required by the explicit model to perform a task to the time required by the Colebrook equation to perform the same task.It means that a model with relative computational efficiency value greater than one (1.0) will require more time than the Colebrook equation to perform a particular task and vice-versa for a model with a value less than one (1.0).
Ten million friction factor calculations were performed using the available explicit models in the ranges of Re and D ε for which the Colebrook equation is valid.These calculations were performed four times and the average was recorded for each of the explicit model.For this analysis, f values for the Colebrook equation were determined using the method developed by Clamond [3] because of its speed of convergence.The relative computational efficiency was thereafter determined based on the approach proposed by Winning and Coole [6].The results are as shown in Table 4.

Model Accuracy, Adequacy and Computational Efficiency
It is observed from Table 3 that the new model (for this study), having the least mean relative and maximum relative errors of 0.0025% and 0.0664%, respectively, is more accurate than the selected extremely accurate models.In addition to the high accuracy of the new model from this study, its parameters are observed to have very low uncertainties ≤ 0.076% (see Table 2).This indicates that the parameters are known precisely.Consequently, the model is deemed very accurate and adequate for predicting friction factor.
It is observed from Table 4 that all the existing extremely accurate models, with the exception of Buzzelli [16] equation, have relative computational efficiencies greater than one (1.0).This is not unexpected, given their complexity with respect to the number of logarithmic functions contained in the models.On the contrary, relative computational efficiency values of less than one have been reported in the work of Winning and Coole [6] for all the extremely accurate models.These values are disputable considering the complexity of these models (in terms of the numbers of logarithmic functions).Our findings show that the Buzzelli [16] model is almost two times faster than the Serghides [15], Romeo, Royo and Monzon [9], Zigrang and Sylvester [14] models.The Buzzelli [16] model has only two logarithmic functions, a combination of logarithm to base ten and the natural logarithm functions.The Buzzelli's [16] model, based on the analysis in this study, is the best existing model in terms of accuracy and relative computational efficiency.However, it is found that that the new model is 39 and 1.9 times (in terms of mean and maximum relative errors, respectively) more accurate than the Buzzelli [16] model (see Table 3).Interestingly, the new model has two logarithmic functions and a higher accuracy (see Figure 2 for error distribution).It has approximately the same relative computational efficiency as the Buzzelli [16] model, which has only two logarithmic functions.Thus, the new model is regarded as a superior model to the existing extremely accurate explicit models.

Conclusion
A new explicit model is developed for predicting friction factor in the range for which the Colebrook equation is valid.Until now, the best predictions are obtained with models having three logarithmic functions.The new simple model having only two logarithmic functions and maximum relative error of 0.0664% in this study is found to be more accurate than the selected existing extremely accurate models.Moreover, the relative computational efficiency (0.3883) of the new model is in close agreement with that (0.3836) of the Buzzelli [16] which was adjudged as the best existing model in this work.Therefore, the new model provides a good trade-off between accuracy and relative computational efficiency.Thus it is superior model to the existing explicit models for estimating pipe friction factor in the fully developed turbulent flow regime.

Figure 1 .
Figure 1.Graph of Re against Re f .

Table 1 .
Existing explicit friction factor models.

Table 2 .
Parameters of the new model.

Table 3 .
Explicit models ordered by maximum relative error.

Table 4 .
Computational efficiencies of the proposed and existing explicit models.