An Accurate and Computationally Efficient Explicit Friction Factor Model

Show more

Received 23 March 2016; accepted 30 March 2016; published 24 May 2016

1. 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)

The major drawback of Colebrook equation is that it is implicit in friction factor (f). Therefore, it requires iteration 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.

2. Review of the Explicit Forms of the Colebrook Equation

2.1. 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.

2.2. 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

Table 1. Existing explicit friction factor models.

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.

3. The Proposed Nonlinear Model

3.1. 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.

3.2. Model Development

The plot of 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, 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:

(19)

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] .

3.3. Performance Criteria

1) (20)

2) (21)

3) (22)

4) (23)

Figure 1. Graph of Re against.

Table 2. Parameters of the new model.

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 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.

3.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.

Figure 2. Relative error distribution for this study, Zigrang and Sylvester [14] , Serghides [15] and Buzzelli [16] models when compared with the implicit Colebrook equation.

Table 3. Explicit models ordered by maximum relative error.

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

4. 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.

Acknowledgements

The authors are grateful to Dr. James F. Whidborne, a Reader at the School of Aerospace, Transport and Manufacturing, Cranfield University, United Kingdom, for his suggestions regarding the technical contents of this paper. Thanks to Emma Hughes, a Doctoral Candidate at La Trobe University, Australia, for proofreading this paper.

References

[1] Colebrook, C.F. (1939) Turbulent Flow in Pipes, with Particular Reference to the Transition Region between the Smooth and Rough Pipe Laws. Journal of the Institution of Civil Engineers, 11, 133-156.

http://dx.doi.org/10.1680/ijoti.1939.13150

[2] Colebrook, C.F. and White, C.M. (1937) Experiments with Fluid Friction Factor in Roughened Pipes. Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences, 161, 367-381.

http://dx.doi.org/10.1098/rspa.1937.0150

[3] Clamond, D. (2009) Efficient Resolution of the Colebrook Equation. Industrial & Engineering Chemistry Research, 48, 3665-3671.

http://dx.doi.org/10.1021/ie801626g

[4] Moody, L.F. (1944) Friction Factors for Pipe Flow. Transactions of the American Society of Mechanical Engineers, 66, 671-681.

[5] Zigrang, D.J. and Sylvester, N.D. (1985) A Review of Explicit Friction Factor Equations. Journal of Energy Resources Technology, 107, 280-283.

http://dx.doi.org/10.1115/1.3231190

[6] Winning, H.K. and Coole, T. (2013) Explicit Friction Factor Accuracy and Computational Efficiency for Turbulent Flows in Pipes. Flow, Turbulence and Combustion, 90, 1-27.

[7] Genic, S., et al. (2011) A Review of Explicit Approximations of Colebrook’s Equation. FME Transactions, 39, 67-71.

[8] Yildirim, G. (2009) Computer-Based Analysis of explicit approximations to the Implicit Colebrook—White Equations in Turbulent Flow Friction Calculation. Advances in Engineering Software, 40, 1183-1190.

http://dx.doi.org/10.1016/j.advengsoft.2009.04.004

[9] Romeo, E., Royo, C. and Monzon, A. (2002) Improved Explicit Equation for Estimation of the Friction Factor in Rough and Smooth Pipes. Chemical Engineering Journal, 86, 369-374.

http://dx.doi.org/10.1016/S1385-8947(01)00254-6

[10] Moody, L.F. (1947) An Approximate Formula for Pipe Friction Factors. Trans ASME, 69, 1005-1011.

[11] Chen, N.H. (1979) An Explicit Equation for Friction Factors in Pipes. Industrial & Engineering Chemistry Fundamentals, 18, 296-297.

http://dx.doi.org/10.1021/i160071a019

[12] Brkic, D. (2011) Review of Explicit Approximations to the Colebrook Relation for Flow Friction. Journal of Petroleum Science and Engineering, 77, 34-48.

http://dx.doi.org/10.1016/j.petrol.2011.02.006

[13] Fang, X., Xu, Y. and Zhou. Z. (2011) New Correlations of Single-Phase Friction Factor for Turbulent Pipe Flow and Evaluation of Existing Single-Phase Friction Factor Correlations. Nuclear Engineering and Design, 241, 897-902.

http://dx.doi.org/10.1016/j.nucengdes.2010.12.019

[14] Zigrang, D.J. and Sylvester. N.D. (1982) Explicit Approximations to the Solution of Colebrook’s Friction Factor Equation. AIChE Journal, 28, 514-515.

http://dx.doi.org/10.1002/aic.690280323

[15] Buzzelli, D. (2008) Calculating Friction in One Step. Machine Design, 80, 54-55.

[16] Seghides, T.K. (1984) Estimate Friction Factor Accurately. Chemical Engineering Journal, 91, 63-64.

[17] Ghanbari, A., Farshad, F. and Rieke, H. (2011) Newly Developed Friction Factor Correlation for Pipe Flow and Flow Assurance. Journal of Chemical Engineering and Materials Science, 2, 83-86.

[18] Glustolisi, O., Berardi, L. and Walski, T.M. (2011) Some Explicit Formulations of Colebrook-White Friction Factor Considering Accuracy vs. Computational Speed. Journal of Hydroinformatics, 13, 401-418.

http://dx.doi.org/10.2166/hydro.2010.098

[19] Barr, D.I.H. (1981) Solutions of the Colebrook-White Function for Resistance to Uniform Turbulent Flow. Proceedings of the Institution of Civil Engineers, Part 2, 71, 529-535.

[20] Schorle, B.J., Churchill, S.W. and Shacham, M. (1980) Comments on: “An Explicit Equation for Friction Factor in Pipe”. Industrial & Engineering Chemistry Fundamentals, 19, 228-229.

http://dx.doi.org/10.1021/i160074a019

[21] Sonnad, J.R. and Goudar, C.T. (2006) Turbulent Flow Friction Factor Calculation Using a Mathematically Exact Alternative to the Colebrook-White Equation. Journal of Hydraulic Engineering, 132, 863-867.

http://dx.doi.org/10.1061/(ASCE)0733-9429(2006)132:8(863)

[22] Haaland, S.E. (1983) Simple and Explicit Formulas for the Friction Factor in Turbulent Pipe Flow. Journal of Fluids Engineering, 105, 89-90.

http://dx.doi.org/10.1115/1.3240948

[23] Manadilli, G. (1997) Replace Implicit Equations with Signomial Functions. Chemical Engineering Journal, 104, 129-130.

[24] Churchill, S.W. (1977) Friction Factor Equations Spans All Fluid-Flow Regimes. Chemical Engineering Journal, 84, 91-92.

[25] Round, G.F. (1980) An Explicit Approximation for the Friction Factor-Reynolds Number Relation for Rough and Smooth Pipes. The Canadian Journal of Chemical Engineering, 58, 122-123.

http://dx.doi.org/10.1002/cjce.5450580119

[26] Brkic, D. (2011) New Explicit Correlations for Turbulent Flow Friction Factor. Nuclear Engineering and Design, 241, 4055-4059.

http://dx.doi.org/10.1016/j.nucengdes.2011.07.042

[27] Rao, A.R. and Kumar, B. (2007) Friction Factor for Turbulent Pipe Flow. Division of Mechanical Sciences, Civil Engineering Indian Institute of Science, Bangalore, ID Code 9587.

[28] Swamee, D.K. and Jain, A.K. (1976) Explicit Equations for Pipe Flow Problems. Journal of the Hydraulics Division, 102, 657-664.

[29] Vantankhah, A.R. and Kouchakzadeh, S. (2008) Discussion of “Turbulent Flow Friction Factor Calculation Using a Mathematically Exact Alternative to Colebrook-White Equation” by Jagadeesh R. Sonnad and Chetan T. Goudar. Journal of Hydraulic Engineering, 134, 1187.

http://dx.doi.org/10.1061/(ASCE)0733-9429(2008)134:8(1187)

[30] Motulsky, H.J. and Christopoulos, A. (2003) Fitting Models to Biological Data Using Linear and Nonlinear Regression: A Practical Guide to Curve Fitting. GraphPad Software Inc., San Diego.