1. Reynolds Number & Flow Regimes
The Reynolds number (Re) is a dimensionless ratio that characterizes the flow regime in a pipe. It compares inertial forces to viscous forces and is the single most important parameter in determining the friction factor:
Where:
- Re = Reynolds number (dimensionless)
- ρ = fluid density (lb/ft³ or kg/m³)
- V = average fluid velocity (ft/s or m/s)
- D = pipe inside diameter (ft or m)
- μ = dynamic viscosity (lb/ft·s or Pa·s)
For liquid flow in practical units (GPM, inches, centipoise), the Reynolds number can be computed as:
Where Q is flow rate in GPM, SG is specific gravity (water = 1.0), d is pipe ID in inches, and μ is viscosity in centipoise. The constant 3,161 incorporates the unit conversions from GPM and inches to consistent units (equivalent to 50.66 × Q × ρ / (d × μ) when ρ is in lb/ft³ rather than SG).
Three Flow Regimes
The Reynolds number divides pipe flow into three distinct regimes, each with fundamentally different friction behavior:
- Laminar flow (Re < 2,100): Fluid moves in smooth, parallel layers. The friction factor is exact: f = 64/Re. Pipe roughness has NO effect. This regime is common in viscous crude oil pipelines and small-diameter instrument tubing.
- Transitional flow (2,100 < Re < 4,000): The flow alternates unpredictably between laminar and turbulent behavior. Friction factors in this zone are unreliable and should be avoided in design. If your operating point falls here, consider changing pipe size or flow rate.
- Turbulent flow (Re > 4,000): Chaotic, mixing flow with eddies and velocity fluctuations. The friction factor depends on BOTH Re and the relative pipe roughness ε/D. This is the most common regime for gas and liquid transmission pipelines.
2. The Colebrook-White Equation
The Colebrook-White equation (Colebrook, 1939) is the industry-standard relationship for determining the Darcy friction factor in turbulent pipe flow. It was developed by combining Prandtl's smooth-pipe law with von Karman's rough-pipe law into a single implicit equation:
Where:
- f = Darcy friction factor (dimensionless)
- ε = absolute pipe roughness (in or mm)
- D = pipe inside diameter (same units as ε)
- ε/D = relative roughness (dimensionless)
- Re = Reynolds number
Why It Is Implicit
The friction factor f appears on BOTH sides of the equation — inside the square root on the left and inside the logarithm on the right. This means you cannot algebraically isolate f. The equation must be solved iteratively.
Newton-Raphson Iterative Solution
The standard approach to solving the Colebrook-White equation is the Newton-Raphson method. Rearranging into the form F(f) = 0:
F(f) = 1/√f + 2 · log10( ε/(3.7D) + 2.51/(Re·√f) ) = 0
Starting from an initial guess (the Swamee-Jain approximation works well), the method converges rapidly — typically 3–5 iterations to reach machine precision (convergence criterion of 1 × 10−10 on the change in f).
Relative Roughness ε/D
The relative roughness is the ratio of the absolute surface roughness of the pipe wall to the inside diameter. It determines how far into the turbulent regime the friction factor remains dependent on Reynolds number versus becoming constant (complete turbulence). At high Re with large ε/D, the friction factor depends only on roughness — this is the "fully rough" or "complete turbulence" zone on the Moody diagram.
3. Explicit Approximations
Because the Colebrook-White equation requires iteration, several explicit (direct-solve) approximations have been developed. The two most widely used are the Swamee-Jain and Churchill equations.
Swamee-Jain Approximation (1976)
This elegant equation gives friction factors accurate to ±1% of the Colebrook-White solution over the ranges:
- Reynolds number: 5,000 ≤ Re ≤ 108
- Relative roughness: 10−6 ≤ ε/D ≤ 10−2
The Swamee-Jain equation is ideal for hand calculations, spreadsheets, and quick field estimates where iteration is inconvenient.
Churchill Equation (1977)
The Churchill correlation is a single equation that covers ALL three flow regimes — laminar, transitional, and turbulent — without requiring separate regime checks:
It uses intermediate terms A and B based on Re and ε/D, then combines them with the laminar solution (f = 64/Re) using a blending function. The result smoothly transitions through all flow regimes in one continuous function.
The Churchill equation is particularly useful in computer programs and process simulators where a single continuous function is preferred over branching logic for different flow regimes.
When to Use Each Method
- Swamee-Jain: Best for hand calculations, spreadsheets, and quick field estimates. Simple, explicit, accurate to ±1%.
- Churchill: Best for computer programs needing one continuous function across all regimes. Avoids if/else branching for laminar vs turbulent.
- Colebrook-White (iterative): The reference standard. Use when maximum accuracy is required or when validating other methods.
4. The Moody Diagram
The Moody diagram (Moody, 1944) is a graphical representation of the Colebrook-White equation. It plots the Darcy friction factor on the vertical axis (log scale) against the Reynolds number on the horizontal axis (log scale), with a family of parametric curves for different values of relative roughness ε/D.
Four Regions on the Diagram
- Laminar line: A single straight line (on the log-log plot) representing f = 64/Re. This line has a slope of −1 and applies for Re < 2,100 regardless of pipe roughness.
- Critical zone (Re = 2,100 – 4,000): A shaded or dashed region where the flow is transitional. Friction factors here are unpredictable, and the diagram intentionally shows no definitive curves.
- Transition region: The curved portion of each ε/D line where the friction factor depends on BOTH Reynolds number and relative roughness. As Re increases, the curves gradually flatten.
- Complete turbulence zone: The flat (horizontal) portion of each ε/D curve at high Reynolds numbers. Here the friction factor depends ONLY on ε/D, and the boundary layer is fully rough. A dashed line on the diagram marks the boundary where each curve becomes essentially flat.
How to Read the Moody Diagram
To find the friction factor: (1) Calculate Re and enter the diagram from the bottom axis. (2) Go up vertically to intersect your ε/D curve (read from the right axis). (3) Read the Darcy friction factor from the left axis. For intermediate ε/D values, interpolate between the printed curves.
Common Pipe Roughness Values
| Pipe Material | ε (inches) | ε (mm) |
|---|---|---|
| Commercial Steel / Wrought Iron | 0.0018 | 0.046 |
| Stainless Steel (new) | 0.000059 | 0.0015 |
| Cast Iron | 0.010 | 0.26 |
| Concrete | 0.012 – 0.12 | 0.3 – 3.0 |
| PVC / Plastic | 0.000059 | 0.0015 |
| Riveted Steel | 0.036 – 0.36 | 0.9 – 9.0 |
| Corroded Steel (10+ years) | 0.010 – 0.040 | 0.25 – 1.0 |
5. Practical Applications
Darcy-Weisbach Pressure Drop
The friction factor feeds directly into the Darcy-Weisbach equation for pressure drop in pipe flow:
Where ΔP is in psi, L is pipe length (ft), D is diameter (ft), ρ is density (lb/ft³), V is velocity (ft/s), and gc = 32.174 lbm·ft/(lbf·s²).
In head loss form: hf = f × (L/D) × (V²/2g), where hf is in feet of fluid.
Pipeline Sizing
The friction factor is central to pipeline design. For a given flow rate and allowable pressure drop, the engineer iterates on pipe diameter: assume a diameter, calculate Re, find f, compute ΔP, and check against the available pressure. The Darcy-Weisbach approach with Colebrook-White friction factors is the preferred method for petroleum liquids, refined products, and gas pipelines.
Hazen-Williams Alternative
For water distribution systems, the Hazen-Williams equation (using a C-factor instead of friction factor) remains widely used due to its simplicity. However, it is empirical, applies only to water near 60°F, and only in the turbulent regime. The Colebrook-White / Darcy-Weisbach approach is more fundamental and should be preferred for petroleum liquids, gas pipelines, and any fluid other than ambient-temperature water.
Common Design Pitfalls
- Darcy vs. Fanning confusion: Using a Fanning f in the Darcy-Weisbach equation (or vice versa) produces a factor-of-4 error in pressure drop. Always verify the friction factor convention.
- Ignoring roughness aging: New pipe roughness values underestimate friction in existing pipelines. Use aged roughness for in-service analysis.
- Not checking Reynolds number range: The Swamee-Jain approximation is not valid below Re = 5,000. Always verify that your operating Reynolds number falls within the valid range of your chosen friction factor method.
- Assuming turbulent flow: Viscous liquids (heavy crude, lube oils) often flow in the laminar regime where f = 64/Re and roughness is irrelevant. Check Re before applying turbulent friction correlations.
Standards & References
- Colebrook (1939): Turbulent flow in pipes with particular reference to the transition region between smooth and rough pipe laws
- Moody (1944): Friction factors for pipe flow (Transactions of ASME)
- Swamee & Jain (1976): Explicit equations for pipe-flow problems (Journal of the Hydraulics Division, ASCE)
- Churchill (1977): Friction-factor equation spans all fluid-flow regimes (Chemical Engineering)
- Crane TP-410: Flow of Fluids Through Valves, Fittings, and Pipe
Ready to calculate friction factors?
→ Launch the Friction Factor & Moody Diagram Calculator