Abstract

The empirical logarithmic Colebrook equation for hydraulic resistance in pipes implicitly considers the unknown flow friction factor. Its explicit approximations, used to avoid iterative computations, should be accurate but also computationally efficient. We present a rational approximate procedure that completely avoids the use of transcendental functions, such as logarithm or non-integer power, which require execution of the additional number of floating-point operations in computer processor units. Instead of these, we use only rational expressions that are executed directly in the processor unit. The rational approximation was found using a combination of a Padé approximant and artificial intelligence (symbolic regression). Numerical experiments in Matlab using 2 million quasi-Monte Carlo samples indicate that the relative error of this new rational approximation does not exceed 0.866%. Moreover, these numerical experiments show that the novel rational approximation is approximately two times faster than the exact solution given by the Wright omega function.

Abstract

The logarithmic Colebrook flow friction equation is implicitly given in respect to an unknown flow friction factor. Traditionally, an explicit approximation of the Colebrook equation requires evaluation of computationally demanding transcendental functions, such as logarithmic, exponential, non-integer power, Lambert W and WrightΩ functions. Conversely, we herein present several computationally cheap explicit approximations of the Colebrook equation that require only one logarithmic function in the initial stage, whilst for the remaining iterations the cheap Padé approximant of the first order is used instead. Moreover, symbolic regression was used for the development of a novel starting point, which significantly reduces the error of internal iterations compared with the fixed value staring point. Despite the starting point using a simple rational function, it reduces the relative error of the approximation with one internal cycle from 1.81% to 0.156% (i.e., by a factor of 11.6), whereas the relative error of the approximation with two internal cycles is reduced from 0.317% to 0.0259% (i.e., by a factor of 12.24). This error analysis uses a sample with 2 million quasi-Monte Carlo points and the Sobol sequence.

Abstract

The unknown friction factor from the implicit Colebrook equation cannot be expressed explicitly in an analytical way, and therefore to simplify the calculation, many explicit approximations can be used instead. The accuracy of such approximations should be evaluated only throughout the domain of interest in engineering practice where the number of test points can be chosen in many different ways, using uniform, quasi-uniform, random, and quasi-random patterns. To avoid picking points with undetected errors, a sufficient minimal number of such points should be chosen, and they should be distributed using proper patterns. A properly chosen pattern can minimize the required number of testing points that are sufficient to detect maximums of the error. The ability of the Sobol quasi-random vs. random distribution of testing points to capture the maximal relative error using a sufficiently small number of samples is evaluated. Sobol testing points that are quasi-randomly distributed can cover the domain of interest more evenly, avoiding large gaps. Sobol sequences are quasi-random and are always the same, which allows the exact repetition of scientific results.