§ 3.1 Module 3 — Head Loss: Friction
One dimensionless number decides which physics you are allowed to use. For everything a municipal station pumps, it is not a close call.
By the end of this lesson
Module 2 left you with total dynamic head as an accounting problem: static lift plus every loss between the two water surfaces. This module prices the largest single entry in that ledger for a long pipeline. But no friction equation is usable until you know which of two entirely different regimes the pipe is in — they have different resistance laws, different velocity profiles, and different lists of what matters, and in one of them the pipe's roughness does not appear at all. One dimensionless number decides. For a force main the honest answer is not "probably turbulent" but "turbulent by a factor of eleven, at the slowest velocity any standard permits, in the smallest pipe any standard permits".
Water driven along a pipe is resisted by two mechanisms that scale differently. Momentum carried bodily by the moving fluid — inertia — grows with ρV²; momentum handed sideways between layers by molecular friction — viscosity — grows with μV/D. Take the ratio and the units cancel completely:
Re = ρVD/μ = VD/ν V mean velocity, m/s
D inside diameter, m
ν kinematic viscosity, m²/s (= μ/ρ)
and because you are always given a FLOW, not a velocity,
substitute V = 4Q/πD² :
Re = 4Q/(πDν) Q flow, m³/s
The two forms you will actually need
That the result is dimensionless is the whole point. One number describes every pipe: a 6 mm laboratory tube and a 900 mm transmission main at the same Reynolds number have geometrically similar flow. That is why a friction chart drawn in the 1940s still sizes pipe today.
Q = 0.040 m³/s D = 0.250 m T = 15 °C
A = πD²/4 = π × 0.250² / 4 = 0.049087 m²
V = Q/A = 0.040 / 0.049087 = 0.8149 m/s
ν at 15 °C = 1.1395 × 10⁻⁶ m²/s
Re = VD/ν = 0.8149 × 0.250 / 1.1395e-6 = 178,780
in one step, from the flow:
Re = 4Q/(πDν) = 4(0.040) / (π × 0.250 × 1.1395e-6)
= 178,780 turbulent
Worked: 40 L/s in a 250 mm force main at 15 °C
The flow form is the one worth memorising, for the sign it exposes: at fixed flow the Reynolds number falls as the pipe gets larger. Velocity drops as D², the explicit D in the numerator recovers only one power of it, and the net dependence is 1/D. The instinct that a bigger pipe is "more turbulent" is exactly backwards — and it is the diameter ceiling of §1.3 wearing different clothes: oversize a force main and one decision lowers its velocity, lowers its Reynolds number and lengthens its detention.
Below Re of about 2,000 the fluid moves in sliding layers that do not trade parcels. All the resistance is viscous shear, and the problem has an exact solution — found forty years before anyone knew why it stopped working: hf = 32μLV/(ρgD²), which in Darcy-Weisbach's bookkeeping is f = 64/Re, exactly. Notice three features before they disappear: head loss is linear in velocity, there is nothing to look up, and the wall's roughness does not appear anywhere. A laminar pipe does not care what it is made of.
Above Re of about 4,000 the layers break down. Eddies carry parcels bodily across the flow, so momentum moves by bulk motion rather than by molecules and the resistance becomes mostly the cost of stirring. Every item above inverts: loss grows with roughly the square of velocity, the wall's texture matters because the eddies reach it, and there is no closed form — f comes from a correlation, which is what §3.2 is for.
| Laminar, Re < 2000 | Turbulent, Re > 4000 | |
|---|---|---|
| Friction factor | f = 64/Re, exact | a correlation in Re and e/D |
| Loss vs flow | hf ∝ Q1, exactly | hf ∝ Q1.8 to Q2.0 |
| Roughness | absent from the solution | first-order; often the design variable |
| Velocity profile | parabola, centreline = 2 × mean | blunt; centreline = 60/49 = 1.2245 × mean |
| Where the shear is | spread across the whole bore | crushed into a thin layer at the wall |
| Cross-stream mixing | molecular only; a dye thread survives | vigorous; a dye thread is gone in a pipe diameter |
| hf vs temperature | strong: hf ∝ μ | weak: 3.7% over 4–25 °C on unlined ductile iron, 11.7% on PVC |
The velocity profile is where the difference stops being notational. Laminar flow is a parabola — centreline exactly twice the mean, shear spread across the whole radius. Turbulent flow is far blunter: Prandtl's one-seventh power law puts the centreline at 60/49 = 1.2245 times the mean, with almost the entire velocity gradient compressed into a thin layer against the wall. That steep wall gradient is the friction.
Interactive 3D instrument
Reynolds\' dye filament, on a force main
A 3D instrument you drive yourself, one variable at a time. It needs JavaScript and WebGL, so it is not shown in this static copy of the page.
A calculated Reynolds number near transition does not justify interpolating confidently between laminar and turbulent formulae. Disturbance, fittings and surface condition can change the regime.
Expose the uncertainty and check both applicable bounds, or obtain service-specific evidence for the low-flow condition.
Between the two figures the pipe is not "somewhat laminar". It is intermittently both: turbulent puffs are born, travel with the bulk flow and die, separated by stretches of undisturbed laminar flow, so a pressure tapping at a fixed station sees resistance flickering between two laws rather than settling on one.
How far apart are those two laws? Take Re = 4,000, where both are still written down. The laminar law gives f = 64/4000 = 0.0160, no iteration and no roughness term. Colebrook-White for 250 mm unlined ductile iron — e/D = 0.26 mm / 250 mm = 1.040 × 10-3 — gives 0.04095: a factor of 2.559. The band is not a region of a few percent uncertainty. It contains two different answers and no basis for choosing.
Reynolds, 1883, and the number he never named
Osborne Reynolds (1842–1912) held the chair of engineering at Owens College, Manchester. His 1883 paper in the Philosophical Transactions of the Royal Society of London (vol. 174, pp. 935–982) carries a title that is itself the lesson: “An experimental investigation of the circumstances which determine whether the motion of water shall be direct or sinuous, and of the law of resistance in parallel channels.” He drew a filament of dye into a glass tube fed from a settling tank, and watched it.
What the retelling usually drops is that he did not find a number. Taking trouble to still the water in the supply tank let the filament survive to far higher speeds; disturbing the inlet pulled breakdown down sharply. He identified a lower value, near 2,000 in modern units, below which any disturbance decayed and order returned — that one is close to a property of pipe flow itself, the same figure whatever Newtonian liquid is in the pipe — and was explicit that the upper figure belonged to his apparatus. The band in every textbook since is Reynolds' own honesty, preserved.
Reynolds himself never wrote “Reynolds number”. Arnold Sommerfeld introduced the term in 1908, at the Fourth International Congress of Mathematicians in Rome; Nikolaus Rott traced the naming in “Note on the History of the Reynolds Number”, Annual Review of Fluid Mechanics 22 (1990), 1–11.
The laminar law was already old: Gotthilf Hagen in Annalen der Physik in 1839, and Jean Léonard Marie Poiseuille, forcing distilled water through glass capillaries a tenth of a millimetre across in pursuit of the circulation of blood, in the Comptes Rendus of 1840–41. Neither could say why so clean a linear law failed in bigger pipes. The turbulent side got its first usable fit from Heinrich Blasius (1883–1970), whose 1913 VDI paper (Forschungsarbeiten Heft 131) gave f = 0.3164·Re-1/4 for smooth pipe — used here as an independent check on our Colebrook solution, matching within 2.8% over Re 104–105.
Two working consequences. First, 2,000 and 4,000 are not the same kind of number: the lower is close to a physical constant, the upper a convention placed anywhere from 3,000 to 5,000. Second, the design rule — if a duty point lands inside the band, do not report a friction factor for it. Bound it with both laws and see whether the decision changes. src/core/hydraulics.js interpolates across the band to keep its curves continuous; every instrument here labels the band rather than quoting that interpolation as a result.
Two requirements already in force push the Reynolds number up from below, and they are requirements rather than tendencies. Ten State Standards asks for at least 0.6 m/s (2 ft/s) in a force main and forbids piping smaller than 100 mm (4 in) on raw wastewater. Put the two minima together — slowest permitted velocity, smallest permitted pipe — at a cold 10 °C, where ν is 1.3079 × 10-6 m²/s:
Re = VD/ν = 0.6 × 0.100 / 1.3079e-6 = 45,900 11.5 × the 4000 floor
Going the other way — the flow at which a pipe LEAVES turbulence,
from the closed-form inverse Q = πDνRe/4 at Re = 4000, 15 °C:
100 mm main 0.358 L/s
250 mm main 0.895 L/s = 3.22 m³/h, about a garden tap
600 mm main 2.148 L/s
Note the direction: at fixed Re, Q ∝ D. The bigger pipe needs MORE
flow to be laminar, because Re falls as 1/D at fixed flow.
The floor under Reynolds number, and the flow needed to escape it
Real duties sit far above that floor: the 40 L/s force main above at Re = 1.79 × 105, a 300 mm potable main at 1.2 m/s at 3.16 × 105, a 900 mm transmission main at 1.5 m/s at 1.18 × 106. Even the gravity sewer arriving at the station is turbulent: a 300 mm sewer half full at 0.4% grade with Manning n = 0.013 flows at 0.865 m/s, and on a hydraulic diameter of 0.300 m that is Re = 2.28 × 105. Reaching the band in the 250 mm main means throttling below 0.9 L/s — a flow at which self-cleansing and detention have been violated for hours.
Three things follow from being firmly turbulent, and they are the working content of this module.
tools/verify-hydraulics.mjs, the explicit Swamee-Jain fit's worst disagreement stays under 1% for Re ≥ 105 over e/D 10-5 to 10-2 — 0.80% at the grid points that file samples, and 0.90% at the true worst point of the range, Re = 105 with e/D ≈ 2.3 × 10-3, which sits between two of its roughness columns. It rises to 2.75% only at Re = 5 × 103 with very rough pipe — worst exactly where municipal flow never is.That last exponent is why a main sized for average flow is punished at peak, and why diameter is the decision that governs. The module epigraph says friction is levied on velocity squared; the exact local exponent is 1.9374 at this duty on unlined ductile iron and 1.8044 on smooth PVC, reaching 2 only in the fully-rough limit §3.2 arrives at. Squared is the right instinct and the wrong number.
None of this makes the Reynolds number a formality: a station contains more than its force main. Where the question stays live is small bore, low flow, or not water.
What this lesson has deliberately not done is name the friction factor for turbulent flow. We used the fact that one exists and quoted values from src/core/hydraulics.js; where 0.02133 comes from, and what the Moody chart's flat right-hand side means, is §3.2. Competing correlations and their error bars are §3.3 and §3.4. Fittings and valves are Module 4, and inside a station they are not minor. The pump is still a black box supplying whatever head the system asks for; what it can actually deliver is Module 6.
Check your understanding
Check your understanding
3 auto-graded questions with an explanation for every wrong answer. Requires JavaScript. (m3-l1-q1)
Lab 3.1
Classify the regime, and find the flow that would change it
Write regimeCheck(qLps, dMm, tC) returning an object with four properties. win.kinematicViscosity(tC) is provided and returns ν in m²/s, so you do not have to fit a viscosity curve; win.RE_LAMINAR_TOP (2000) and win.RE_TURBULENT_FLOOR (4000) are provided too. vMs — mean velocity, m/s. Flow over area, with A = πD²/4. Mind the units: the inputs are litres per second and millimetres. re — the Reynolds number. Either form will do, but 4Q/(πDν) avoids computing the area twice. regime — the string 'laminar' , 'transitional' or 'turbulent' . Below RE_LAMINAR_TOP is laminar; at or above RE_TURBULENT_FLOOR is turbulent; between them is transitional, and the whole point of the third name is that it is not a synonym for either of the other two. minTurbulentFlowLps — the flow, in L/s, at which this pipe at this temperature sits exactly at RE_TURBULENT_FLOOR . Invert Re = 4Q/(πDν); no iteration is needed. This is the number that tells you whether the regime question is live for a given main, and nobody has it to hand. Note what is not an argument: roughness. It has no place in either the Reynolds number or the classification. Graded in the browser against 7 assertions; the editor and harness require JavaScript.
Compute the Reynolds number once, at the start, and write it on the calculation. It is the sentence that says which of the following pages you are allowed to use.