The friction factor x for turbulent flow in a pipe of diameter D and interior roughness characterized by a roughness coefficient σ is given by Colebrook’s formula:
where Re is the Reynolds number indicating laminar or turbulent flow. Implement OCTAVE program to find the root of the function f(x) in (2), that is, the value of the friction factor x which solves the nonlinear equation, f(x) = 0, using following methods:
- Bisection method
- Regular Falsi method
- Secant method
Read the input data from keyboard or from a file. The input data are D, σ, Re, accuracy level of your program, the left endpoint of the selected initial and the right endpoint of the selected interval. You can assume that D = 0.1 m, σ = 0.0025 m and Re = 3 × 104