8.2 Fluid Dynamics and Flow

Laminar vs Turbulent Flow

Flow Characteristics and Reynolds Number

When fluid moves through a pipe or channel, it can behave in two fundamentally different ways. Understanding which type of flow you're dealing with is the first step in any fluid dynamics problem.

Laminar flow involves smooth, parallel layers of fluid moving in a predictable pattern with no mixing between layers. Think of honey slowly pouring off a spoon.

Turbulent flow features irregular fluctuations and mixing between fluid layers, resulting in chaotic motion. Think of rapids in a river.

The Reynolds number (Re) is the tool that predicts which flow regime you'll get. It's a dimensionless quantity that compares inertial forces (which promote turbulence) to viscous forces (which resist it):

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

where $\rho$ is fluid density, $V$ is velocity, $D$ is the characteristic length (like pipe diameter), and $\mu$ is dynamic viscosity.

For pipe flow, the critical thresholds are:

• $Re < 2300$: laminar flow
• $2300 < Re < 4000$: transitional (could go either way)
• $Re > 4000$: turbulent flow

Factors Influencing Flow Regime

Several factors push flow toward laminar or turbulent behavior:

• Fluid velocity: Higher speeds increase inertial forces, promoting turbulence
• Viscosity: Higher viscosity fluids (like oil) resist turbulent motion, so they stay laminar at higher velocities
• Density: Denser fluids have greater inertial forces, making turbulence more likely
• Pipe diameter: Larger diameters give turbulent eddies more room to develop
• Surface roughness: Rough pipe walls create disturbances that can trigger the transition to turbulence
• Obstructions and geometry changes: Valves, sudden expansions, and contractions disrupt smooth flow and generate turbulence

Open Channel Flow Considerations

Open channels (rivers, drainage ditches, canals) add another layer of complexity because there's a free surface exposed to the atmosphere.

The Froude number (Fr) is used alongside the Reynolds number to characterize open channel flow. While Reynolds number compares inertial to viscous forces, the Froude number compares inertial forces to gravitational forces:

$Fr = \frac{V}{\sqrt{gD_h}}$

where $D_h$ is the hydraulic depth and $g$ is gravitational acceleration.

Open channel flow regimes based on Froude number:

• $Fr < 1$: subcritical flow (slow, deep, tranquil)
• $Fr = 1$: critical flow
• $Fr > 1$: supercritical flow (fast, shallow, rapid)

Channel slope and hydraulic depth both influence which regime develops.

Continuity Equation for Fluid Flow

Fundamental Principles

The continuity equation comes from a simple but powerful idea: mass is conserved. Fluid doesn't just appear or disappear as it flows through a system.

For incompressible fluids (most liquids in civil engineering applications), this means the volumetric flow rate stays constant along a streamline:

$Q = A_1V_1 = A_2V_2$

where $Q$ is the volumetric flow rate, $A$ is the cross-sectional area, and $V$ is the fluid velocity.

In plain terms: if a pipe gets narrower, the fluid must speed up. If it gets wider, the fluid slows down. The volume passing any cross-section per unit time stays the same.

For steady-state flow, the mass flow rate entering a control volume equals the mass flow rate leaving it.

Applications and Variations

The continuity equation applies to both closed conduits (pipes) and open channels.

For compressible fluids (gases at high speeds or large pressure changes), density can vary, so the equation becomes:

$\rho_1 A_1 V_1 = \rho_2 A_2 V_2$

where $\rho$ is the fluid density at each section.

Common applications include:

• Fire hose nozzles: A nozzle narrows the flow area dramatically. If a hose has an internal diameter of 64 mm and the nozzle exit is 19 mm, the continuity equation tells you the exit velocity is roughly 11 times the hose velocity.
• Expanding pipe sections: When a pipe gradually widens, you can calculate how much the velocity decreases at the larger section.

Practical Considerations

• Temperature and pressure changes along the flow path can affect fluid density, which matters for long pipelines or systems with significant elevation changes
• For unsteady (time-varying) flow, the continuity equation must account for changes in stored volume over time
• Multi-phase flows (like water with air bubbles) require considering the volume fraction of each phase

Bernoulli's Equation for Fluid Dynamics

Fundamental Principles

Bernoulli's equation is derived from conservation of energy applied to a moving fluid. It relates three forms of energy along a streamline: pressure energy, kinetic energy, and potential energy.

$P_1 + \frac{1}{2}\rho V_1^2 + \rho g h_1 = P_2 + \frac{1}{2}\rho V_2^2 + \rho g h_2$

where:

• $P$ = pressure
• $\rho$ = fluid density
• $V$ = velocity
• $g$ = gravitational acceleration
• $h$ = elevation

The core idea: if one form of energy increases, at least one of the others must decrease. When fluid speeds up (higher $V$), pressure drops (lower $P$). This is why airplane wings generate lift and why water sprays faster from a narrowed garden hose.

Key assumptions: steady flow, incompressible fluid, inviscid (frictionless) flow, and along a single streamline.

Applications and Problem Solving

Bernoulli's equation lets you solve for an unknown variable when the others are known. Here's a typical approach:

1. Identify two points along the streamline
2. Write Bernoulli's equation for those two points
3. Cancel terms that are equal or zero at both points
4. Solve for the unknown

Example: To find the exit velocity of water draining from a large open tank, set point 1 at the water surface and point 2 at the outlet. The surface velocity is approximately zero (large tank), and both points are at atmospheric pressure. The equation simplifies to Torricelli's theorem: $V = \sqrt{2gh}$, where $h$ is the height of water above the outlet.

Two useful visualization tools come from Bernoulli's equation:

• Hydraulic Grade Line (HGL): represents the sum of pressure head and elevation head at each point
• Energy Grade Line (EGL): adds velocity head on top of the HGL, representing total energy

The EGL is always above the HGL, and the vertical distance between them equals the velocity head $\frac{V^2}{2g}$.

Limitations and Considerations

Bernoulli's equation assumes no energy losses, which never perfectly holds in real systems. For real-world applications:

• Friction losses must be added separately (covered in the next section)
• Compressible flows or flows with significant density changes need modified forms
• Unsteady flow conditions and flows with significant energy exchanges (like pumps or turbines) require additional terms

Pressure Losses in Pipe Systems

Major Losses Due to Friction

In real pipe systems, fluid loses energy to friction along the pipe walls. These are called major losses because in long pipelines, they account for most of the total energy loss.

The Darcy-Weisbach equation calculates friction head loss:

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

where:

• $f$ = Darcy friction factor (dimensionless)
• $L$ = pipe length
• $D$ = pipe diameter
• $V$ = fluid velocity
• $g$ = gravitational acceleration

Finding the friction factor $f$ depends on the flow regime:

• Laminar flow: $f = \frac{64}{Re}$ (straightforward calculation)
• Turbulent flow: Use the Moody diagram or the Colebrook-White equation, which requires both the Reynolds number and the pipe's relative roughness ($\varepsilon / D$)

For non-circular conduits and open channels, replace $D$ with $4R_h$, where $R_h$ is the hydraulic radius (cross-sectional area divided by wetted perimeter).

Minor Losses Due to Fittings

Valves, elbows, tees, expansions, and contractions all cause additional energy losses called minor losses. (The name is misleading; in short pipe systems with many fittings, these can actually dominate.)

Minor losses are calculated as:

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

where $K$ is a loss coefficient specific to each fitting type. These $K$ values are found in reference tables.

Total head loss in a system is the sum of all major and minor losses:

$h_{total} = h_f + \sum h_m$

The equivalent length method is a handy shortcut: it converts each fitting's minor loss into an equivalent length of straight pipe, so you can lump everything into one Darcy-Weisbach calculation.

Practical Considerations

• Pipe material and age significantly affect surface roughness. Old cast iron pipes can have roughness values many times higher than new pipes due to corrosion and deposits.
• Temperature changes affect fluid viscosity, which in turn changes the Reynolds number and friction factor
• Accurate pressure drop calculations are essential for sizing pumps and ensuring the system delivers adequate flow

Flow Analysis in Pipe Networks

Hardy Cross Method

Real water distribution systems aren't single pipes; they're interconnected networks with loops and junctions. The Hardy Cross method is a classic iterative technique for finding flow rates and pressures throughout these networks.

It relies on two principles:

• Conservation of mass: Flow into each junction must equal flow out
• Conservation of energy: The net head loss around any closed loop must equal zero

Steps for the Hardy Cross method:

1. Assume initial flow rates in each pipe (these are guesses, but they must satisfy continuity at every junction)
2. Calculate the head loss in each pipe using the Darcy-Weisbach equation
3. For each loop, compute the net head loss (sum of head losses around the loop, with sign conventions for flow direction)
4. Calculate a correction factor for each loop and adjust the flow rates
5. Repeat steps 2-4 until the net head loss around every loop is close enough to zero (convergence)

A typical application is analyzing flow distribution in a municipal water supply network to ensure adequate pressure at all service points.

Alternative Techniques and Considerations

• The Newton-Raphson method converges faster than Hardy Cross for large networks but requires more complex setup
• The linear theory method offers another approach by linearizing the head loss equations
• For large-scale networks, computer software (like EPANET) handles the computational complexity

Realistic models also incorporate:

• Pump characteristic curves
• Tank water levels and their variation over time
• Pressure-dependent demands (e.g., simulating the effect of fire hydrant usage on system pressures)

Advanced Network Analysis

• Extended period simulation accounts for how demands change throughout the day (peak morning usage, low overnight demand)
• Water quality modeling can be integrated with hydraulic analysis to track contaminant movement or chlorine residual decay through the network
• Optimization techniques help engineers design networks that meet performance requirements at minimum cost