Colebrook-White Equation

Computes the friction factor in pipes for given values of the Reynolds number (Re) and the relative roughness coefficient (epsilon).

MATLAB code to compute the friction factor in pipes for given values of the Reynolds number (Re) and the relative roughness coefficient (epsilon).

Syntax:
f = colebrook(Re,epsilon)

Example 1: Single Re, single epsilon
Re = 1e5;
epsilon = 1e-4;
f = colebrook(Re,epsilon)

Example 2: Multiple Re, single epsilon
Re = 5000:1000:100000;
epsilon = 1e-4;
f = colebrook(Re,epsilon);
plot(Re,f)

Example 3: Single Re, multiple epsilon
Re = 1e5;
epsilon = linspace(1e-4,1e-1,100);
f = colebrook(Re,epsilon);
plot(epsilon,f)

Example 4: Multiple Re, multiple epsilon
Re = logspace(4,8,100);
epsilon = linspace(1e-4,1e-1,100);
[RE,EPSILON] = meshgrid(Re,epsilon);
F = colebrook(RE,EPSILON);
surf(RE,EPSILON,F)

References:
[1] Colebrook, C. F., & White, C. M. (1937). Experiments with fluid
friction in roughened pipes. Proceedings of the Royal Society of
London. Series A - Mathematical and Physical Sciences, 161(906),
367-381.
[2] Colebrook, C. (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(4),
133-156.