9.1 Laminar and Turbulent Flow

Fluid flow in pipes falls into two main categories: laminar or turbulent. The distinction between them governs how you calculate pressure drops, friction losses, and velocity distributions. The Reynolds number is the single most important parameter for determining which regime you're dealing with, and it shows up constantly in pipe flow analysis.

This guide covers how to identify flow regimes, apply the Hagen-Poiseuille equation for laminar flow, work with friction factors for turbulent flow, and interpret velocity and shear stress profiles.

Laminar and Turbulent Flow in Pipes

Laminar vs turbulent flow

The Reynolds number ($Re$) is the dimensionless parameter that tells you which flow regime you're in:

$Re = \frac{\rho V D}{\mu}$

where $\rho$ is fluid density, $V$ is the average (mean) velocity, $D$ is the pipe's internal diameter, and $\mu$ is dynamic viscosity. Physically, $Re$ represents the ratio of inertial forces to viscous forces. When viscous forces dominate, flow stays orderly; when inertial forces dominate, flow becomes chaotic.

• Laminar flow ($Re < 2300$)
  • Fluid moves in smooth, parallel layers with no cross-stream mixing
  • The velocity profile is parabolic: maximum velocity at the pipe center, zero at the wall (no-slip condition)
  • Common in slow or highly viscous flows, such as oil in small-diameter hydraulic lines or flow through capillary tubes
• Turbulent flow ($Re > 4000$)
  • Fluid motion is irregular and chaotic, with significant mixing between layers
  • The velocity profile is much flatter than in laminar flow, meaning the velocity is more uniform across the cross-section, but it still drops to zero at the wall
  • Common in most practical engineering systems: water mains, HVAC ducts, fast-flowing rivers
• Transitional flow ($2300 < Re < 4000$)
  • The flow is unstable and can switch between laminar and turbulent behavior unpredictably
  • Design calculations in this range are unreliable, so engineers typically try to operate clearly in one regime or the other

Hagen-Poiseuille equation applications

The Hagen-Poiseuille equation gives the pressure drop for fully developed laminar flow in a circular pipe:

$\Delta P = \frac{128 \mu L Q}{\pi D^4}$

where $L$ is the pipe length and $Q$ is the volumetric flow rate. Rearranging to solve for flow rate:

$Q = \frac{\pi D^4 \Delta P}{128 \mu L}$

Notice the $D^4$ dependence. This is critical: if you double the pipe diameter, the flow rate increases by a factor of 16 for the same pressure drop. Conversely, even a small reduction in diameter (say from plaque buildup in an artery) causes a dramatic increase in resistance to flow.

The equation assumes:

• Steady, fully developed flow (not in the entrance region)
• Incompressible, Newtonian fluid
• Laminar regime ($Re < 2300$)

Typical applications include:

1. Microfluidic devices (lab-on-a-chip systems) where flow rates are tiny and $Re$ is well below 2300
2. Pressure drop calculations in small-diameter tubing such as hydraulic lines and capillary tubes
3. Biomedical analysis of blood flow in small vessels like arterioles and capillaries, where the laminar assumption is reasonable

Turbulent Flow and Pipe Flow Analysis

Friction factors in turbulent flow

For turbulent flow, the Hagen-Poiseuille equation no longer applies. Instead, you use the Darcy-Weisbach equation:

$\Delta P = f \frac{L}{D} \frac{\rho V^2}{2}$

where $f$ is the Darcy friction factor. This equation actually works for both laminar and turbulent flow (for laminar flow, $f = 64/Re$), but it's most commonly associated with turbulent analysis because that's where finding $f$ gets complicated.

The friction factor in turbulent flow depends on two things: the Reynolds number and the relative roughness $\varepsilon / D$, which is the ratio of the pipe wall's average roughness height ($\varepsilon$) to the pipe diameter. Smooth pipes like PVC or glass have very small $\varepsilon$; rough pipes like cast iron or concrete have much larger values.

Two main tools for finding $f$:

• Moody diagram: A graphical chart that plots $f$ against $Re$ for various values of $\varepsilon / D$. You look up your Reynolds number on the x-axis, follow the curve for your pipe's relative roughness, and read $f$ off the y-axis. This is the fastest approach for quick estimates.
• Colebrook equation: An implicit formula for $f$ in the turbulent regime:

$\frac{1}{\sqrt{f}} = -2.0 \log_{10}\left(\frac{\varepsilon/D}{3.7} + \frac{2.51}{Re\sqrt{f}}\right)$

Because $f$ appears on both sides, you can't solve it directly. You need an iterative approach:

1. Guess an initial value of $f$ (a common starting point is 0.02)
2. Plug it into the right side and compute a new $f$
3. Repeat until the value converges (typically 3-5 iterations)

Alternatively, the Haaland equation gives an explicit approximation that's accurate to within about 2% and avoids iteration entirely.

Velocity and shear stress profiles

Laminar flow velocity profile:

The velocity at any radial position $r$ from the center is:

$u(r) = 2V\left[1 - \left(\frac{r}{R}\right)^2\right]$

where $R$ is the pipe radius. This is a parabola: the maximum velocity at the centerline ($r = 0$) is exactly twice the mean velocity $V$, and the velocity drops to zero at the wall ($r = R$).

Laminar flow shear stress:

Shear stress varies linearly from zero at the center to a maximum at the wall:

$\tau_w = \frac{4\mu V}{R}$

This linear distribution makes sense because the velocity gradient is steepest right at the wall and zero at the centerline where the profile is flat.

Turbulent flow velocity profile:

The profile is much flatter than the parabolic laminar case. A common approximation is the power-law profile:

$\frac{u}{u_{\max}} = \left(1 - \frac{r}{R}\right)^{1/n}$

where $n$ depends on the Reynolds number ($n \approx 7$ is typical for many practical flows). Because the profile is flatter, the mean velocity is much closer to the maximum velocity than in laminar flow.

Turbulent flow shear stress:

The wall shear stress is related to the friction factor:

$\tau_w = \frac{1}{8} f \rho V^2$

Unlike the linear laminar distribution, turbulent shear stress varies non-linearly across the cross-section, with much higher values concentrated near the wall due to the steep velocity gradient in the thin viscous sublayer.

Why these profiles matter in practice:

• Piping system design: Knowing the velocity profile helps you minimize pressure drop and size pumps correctly
• Heat transfer analysis: Convective heat transfer rates depend directly on the near-wall velocity gradient and shear stress
• Erosion and corrosion: Regions of high wall shear stress experience more material wear, which affects pipe lifespan and maintenance scheduling