5.1 Laminar and turbulent flows

Laminar vs Turbulent Flow

Laminar and turbulent flow describe the two fundamental regimes of fluid motion. Laminar flow moves in smooth, orderly layers, while turbulent flow is chaotic with intense mixing. Distinguishing between these regimes, predicting when transition occurs, and quantifying their effects on pressure drop and drag are central problems in viscous flow analysis.

Laminar vs Turbulent Flow

Laminar flow consists of smooth, parallel layers (laminae) of fluid sliding past one another with no macroscopic mixing between layers. Momentum transfer occurs only through molecular viscosity.

Turbulent flow is irregular and chaotic. Fluid parcels mix across streamlines, velocity fluctuates in both magnitude and direction, and momentum transfer is dominated by turbulent (eddy) stresses rather than molecular viscosity alone.

The practical distinction matters because these two regimes produce very different pressure drops, heat transfer rates, and drag forces. A few rules of thumb:

• Laminar flow tends to occur at low velocities, in small-diameter conduits, or in highly viscous fluids.
• Turbulent flow tends to occur at high velocities, in large-diameter conduits, or in low-viscosity fluids.
• Turbulent flow has higher wall shear stress and friction losses but also much better mixing and heat transfer.

Reynolds Number

The Reynolds number is the dimensionless ratio of inertial forces to viscous forces in a flow:

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

where $\rho$ is fluid density, $v$ is a characteristic velocity, $D$ is a characteristic length (e.g., pipe diameter), and $\mu$ is dynamic viscosity. You can also write it as $Re = vD/\nu$, where $\nu = \mu/\rho$ is the kinematic viscosity.

A high $Re$ means inertial forces dominate and the flow is prone to turbulence. A low $Re$ means viscous forces dominate and the flow stays orderly.

Critical Reynolds Number

The critical Reynolds number is the value of $Re$ at which transition from laminar to turbulent flow begins.

• For internal flow in a circular pipe: $Re_{crit} \approx 2300$
• For flow over a flat plate: $Re_{x,crit} \approx 5 \times 10^5$ (based on distance $x$ from the leading edge)

These values are approximate. In very clean, disturbance-free laboratory conditions, laminar pipe flow has been maintained up to $Re \approx 100{,}000$. In practice, surface roughness, vibrations, and inlet disturbances cause transition near the standard critical values.

Predicting Flow Regime

1. Calculate $Re$ using the appropriate characteristic length and velocity for your geometry.
2. Compare to the critical value for that geometry.
3. If $Re < Re_{crit}$, expect laminar flow. If $Re > Re_{crit}$, expect turbulent flow.

Between roughly $2300$ and $4000$ in pipe flow, you're in a transitional zone where the flow can flicker between laminar and turbulent states. Design calculations in this range carry extra uncertainty.

Characteristics of Laminar Flow

Velocity Profile

In fully developed laminar pipe flow, the velocity profile is parabolic. The fluid at the wall satisfies the no-slip condition ($v = 0$), and the maximum velocity occurs at the centerline.

The profile is given by the Hagen-Poiseuille solution:

$v(r) = \frac{\Delta P}{4\mu L}(R^2 - r^2)$

where $\Delta P$ is the pressure drop over length $L$, $R$ is the pipe radius, and $r$ is the radial distance from the center. The average velocity across the cross-section is exactly half the centerline maximum: $v_{avg} = v_{max}/2$.

Pressure Drop

The Hagen-Poiseuille equation for volumetric flow rate $Q$ gives the pressure drop as:

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

This can also be written as $\Delta P = \frac{8\mu LQ}{\pi R^4}$. The key takeaway: pressure drop scales linearly with velocity (or flow rate) and is extremely sensitive to pipe diameter (inversely proportional to $D^4$). Halving the pipe diameter increases the pressure drop by a factor of 16.

Entrance Length

Near a pipe inlet, the velocity profile is still developing as the boundary layers grow inward from the walls. The entrance length $L_e$ is the distance required for the profile to become fully developed (parabolic).

For laminar flow:

$L_e \approx 0.06 \, Re \, D$

For example, at $Re = 1000$ in a 2 cm diameter pipe, $L_e \approx 0.06 \times 1000 \times 0.02 = 1.2$ m. Flow in the entrance region has a higher pressure drop than fully developed flow because the boundary layers are still accelerating the core fluid.

Characteristics of Turbulent Flow

Velocity Fluctuations

In turbulent flow, the instantaneous velocity at any point fluctuates around a time-averaged mean. You decompose the velocity using Reynolds decomposition:

$v = \bar{v} + v'$

where $\bar{v}$ is the time-averaged velocity and $v'$ is the fluctuating component. The time average of $v'$ is zero by definition, but $\overline{v'^2}$ is not. The turbulence intensity is defined as the ratio of the root-mean-square of the fluctuations to the mean velocity:

$I = \frac{\sqrt{\overline{v'^2}}}{\bar{v}}$

Typical turbulence intensities in pipe flow range from about 1% to 10%.

Turbulent Eddies and the Energy Cascade

Turbulent flow contains rotating fluid structures called eddies spanning a wide range of sizes. The largest eddies are set by the geometry (e.g., pipe diameter), while the smallest are set by viscosity.

Energy flows from large eddies to progressively smaller ones in a process called the energy cascade. At the smallest scales (the Kolmogorov microscale), kinetic energy is finally dissipated into heat by viscosity. This cascade is why turbulence requires a continuous energy input to sustain itself.

Turbulent Boundary Layer Structure

A turbulent boundary layer over a surface has three distinct regions:

• Viscous sublayer: A very thin layer right at the wall where viscous effects dominate and the velocity profile is nearly linear.
• Buffer layer: A transitional region where both viscous and turbulent stresses are significant.
• Fully turbulent (log-law) region: The bulk of the boundary layer, where turbulent mixing dominates and the velocity follows a logarithmic profile.

The turbulent boundary layer grows faster with downstream distance than a laminar one, and it produces higher wall shear stress.

Transition from Laminar to Turbulent

Transition is not a sudden switch. It occurs over a transition region where the flow intermittently alternates between laminar and turbulent patches.

Transition Mechanisms

• Natural (Tollmien-Schlichting) transition: Small disturbances in the laminar flow are amplified through linear instability, eventually breaking down into turbulence. This is the classical route.
• Bypass transition: Large disturbances (high freestream turbulence, surface roughness) skip the linear instability stage and trigger turbulence directly.
• Separation-induced transition: The laminar boundary layer separates from the surface (e.g., at a sharp edge), and the separated shear layer transitions to turbulence. The turbulent flow may then reattach downstream.

Factors Affecting Transition

• Surface roughness introduces disturbances that promote earlier transition.
• Pressure gradient: An adverse pressure gradient (decelerating flow) destabilizes the boundary layer and promotes transition. A favorable pressure gradient (accelerating flow) stabilizes it and delays transition.
• Freestream turbulence: Higher turbulence levels in the external flow trigger bypass transition at lower $Re$.
• Wall suction or blowing: Suction stabilizes the boundary layer; blowing destabilizes it.

Friction Factors

The Darcy friction factor $f$ relates pressure loss to flow conditions in a pipe. It appears in the Darcy-Weisbach equation (covered below) and depends on whether the flow is laminar or turbulent.

Laminar Friction Factor

For fully developed laminar flow in a circular pipe:

$f = \frac{64}{Re}$

This is an exact analytical result derived from the Hagen-Poiseuille solution. The friction factor decreases as $Re$ increases because viscous effects become relatively less dominant (though the flow eventually transitions to turbulence before $f$ drops too far).

Turbulent Friction Factor

In turbulent flow, $f$ depends on both $Re$ and the relative roughness $\epsilon/D$, where $\epsilon$ is the average roughness height of the pipe wall.

The Colebrook equation is the standard implicit relation:

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

Because $f$ appears on both sides, you need to solve iteratively (or use an explicit approximation like the Swamee-Jain equation). For smooth pipes ($\epsilon/D \to 0$), the roughness term drops out and $f$ depends only on $Re$.

Moody Diagram

The Moody diagram plots $f$ versus $Re$ on a log-log scale, with curves for different values of $\epsilon/D$. It's the graphical equivalent of the Colebrook equation and is extremely useful for quick lookups.

How to read it:

1. Find your $Re$ on the horizontal axis.
2. Locate the curve corresponding to your pipe's relative roughness.
3. Read $f$ from the vertical axis.

The diagram clearly shows the laminar region (a single straight line with slope $-1$), the transition zone, and the family of turbulent curves. At very high $Re$, the curves flatten out, meaning $f$ becomes independent of $Re$ and depends only on roughness. This is the fully rough regime.

Flow in Pipes

Laminar Pipe Flow

Fully developed laminar pipe flow is one of the few cases with an exact analytical solution. The parabolic velocity profile and the Hagen-Poiseuille pressure drop relation (covered above) apply. Pressure drop scales linearly with flow rate.

Turbulent Pipe Flow

The turbulent velocity profile is much flatter than the parabolic laminar profile. Most of the velocity variation is concentrated near the wall. The ratio $v_{avg}/v_{max}$ is typically around 0.8 to 0.85 in turbulent flow, compared to exactly 0.5 in laminar flow.

Pressure drop in turbulent flow scales approximately with the square of the velocity, making it significantly higher than laminar pressure drop at the same flow rate.

Pressure Losses

Total pressure loss in a pipe system has two components:

Major (friction) losses from wall shear along straight pipe sections, calculated with the Darcy-Weisbach equation:

$h_f = f \frac{L}{D} \frac{v^2}{2g}$

where $h_f$ is the head loss, $f$ is the Darcy friction factor, $L$ is pipe length, $D$ is diameter, $v$ is mean velocity, and $g$ is gravitational acceleration.

Minor losses from fittings, bends, valves, expansions, and contractions:

$h_m = K \frac{v^2}{2g}$

where $K$ is a loss coefficient specific to each fitting. Despite the name "minor," these losses can dominate in systems with many fittings and short pipe runs.

Flow over Surfaces

Laminar Boundary Layer

When a uniform flow encounters a flat plate, a laminar boundary layer grows from the leading edge. The Blasius solution gives the boundary layer thickness as:

$\delta \approx \frac{5x}{\sqrt{Re_x}}$

where $Re_x = \rho v_\infty x / \mu$ is the local Reynolds number based on distance $x$ from the leading edge. The boundary layer grows as $\sqrt{x}$, so it thickens slowly.

Turbulent Boundary Layer

Once $Re_x$ exceeds approximately $5 \times 10^5$, the boundary layer transitions to turbulent. The turbulent boundary layer grows faster (roughly as $x^{4/5}$) and has higher skin friction than the laminar layer. However, it is also more resistant to separation because the turbulent mixing brings high-momentum fluid from the outer flow down toward the wall.

Separation and Reattachment

Flow separation occurs when the boundary layer encounters a strong enough adverse pressure gradient that the near-wall fluid decelerates to zero velocity and reverses direction. At that point, the boundary layer detaches from the surface.

Consequences of separation include:

• A large increase in pressure drag (form drag)
• Reduced lift on airfoils (stall occurs when separation covers most of the upper surface)
• Unsteady phenomena such as vortex shedding, which can cause structural vibrations

After separation, the flow may reattach downstream, enclosing a recirculation zone called a separation bubble. Turbulent boundary layers separate later than laminar ones because their higher near-wall momentum resists the adverse pressure gradient more effectively. This is why golf balls have dimples: the roughness triggers turbulence, delays separation, and reduces the wake, cutting overall drag.

Turbulence Modeling

Directly resolving every scale of turbulence in a simulation is usually impractical, so computational fluid dynamics (CFD) relies on turbulence models of varying fidelity and cost.

Reynolds-Averaged Navier-Stokes (RANS)

RANS models decompose every flow variable into a time-averaged part and a fluctuating part, then solve the time-averaged equations. This averaging introduces unknown terms called Reynolds stresses ($\overline{u_i' u_j'}$), which must be modeled.

Common RANS closure models include:

• $k$-$\epsilon$ model: Solves transport equations for turbulent kinetic energy $k$ and its dissipation rate $\epsilon$. Robust and widely used for general-purpose industrial flows, but struggles with strong separation and swirling flows.
• $k$-$\omega$ model: Solves for $k$ and the specific dissipation rate $\omega$. Better near walls than $k$-$\epsilon$, making it popular for boundary layer flows.
• SST (Shear Stress Transport): A hybrid that blends $k$-$\omega$ near walls with $k$-$\epsilon$ in the freestream, combining the strengths of both.

RANS is computationally cheap and suitable for most engineering design work, but it only provides time-averaged results and can miss important unsteady features.

Large Eddy Simulation (LES)

LES directly resolves the large, energy-carrying eddies and models only the small-scale (subgrid-scale) eddies using a subgrid-scale model. The Navier-Stokes equations are spatially filtered to remove scales smaller than the grid size.

LES captures much more flow detail than RANS, including unsteady vortex dynamics and large-scale mixing. The cost is significantly higher because it requires fine grids and time-resolved simulations, especially near walls where the turbulent scales become very small.

Direct Numerical Simulation (DNS)

DNS resolves every scale of turbulent motion, from the largest energy-containing eddies down to the Kolmogorov dissipation scale. No modeling is involved; the Navier-Stokes equations are solved exactly (within numerical discretization error).

The computational cost scales roughly as $Re^3$ (in 3D), which limits DNS to relatively low Reynolds numbers and simple geometries. DNS is primarily a research tool used to generate benchmark data and to develop and validate turbulence models.

Applications of Laminar and Turbulent Flows

Heat Transfer

Turbulent flow dramatically enhances convective heat transfer because eddies transport thermal energy across the flow much faster than molecular conduction alone. Heat exchangers are designed to operate in the turbulent regime to maximize the heat transfer coefficient.

Laminar flow is preferred when uniform, gentle heating or cooling is needed, such as in some electronic cooling applications or in processes where thermal stress must be minimized.

Mass Transfer

The same mixing that enhances heat transfer also enhances mass transfer. Chemical reactors often operate in turbulent flow to ensure rapid mixing of reactants and uniform concentration fields, which increases reaction rates.

Laminar flow is used where precise, controlled transport is needed. Microfluidic devices for biomedical applications, for example, rely on laminar flow to keep fluid streams separate and allow diffusion-controlled mixing.

Mixing and Dispersion

In environmental fluid mechanics, turbulent mixing governs the dispersion of pollutants in rivers, the atmosphere, and the ocean. Predicting how a contaminant plume spreads requires understanding turbulent diffusion, which is orders of magnitude faster than molecular diffusion.

Conversely, laminar flow is exploited in manufacturing processes (e.g., layered composite fabrication) where mixing must be minimized to maintain distinct material boundaries.