Fluid Dynamics: Energy & Momentum

Learning Objectives

  • Relate fluid motion to pressure, gravity, viscous stress, and external work.
  • State the assumptions behind Euler's equation, Bernoulli's equation, and the Navier-Stokes equations.
  • Apply the extended mechanical-energy equation with pumps, turbines, friction, minor losses, and kinetic-energy correction factors.
  • Construct and interpret the hydraulic grade line and energy grade line.
  • Apply linear momentum to pipe bends, nozzles, jets, reducers, gates, and moving vanes using a consistent force sign convention.
  • Use momentum correction factors and angular momentum where velocity profiles or rotating machinery require them.
  • Recognize cavitation, vapor-pressure, absolute-pressure, and model-validity limits.

Fluid dynamics connects the kinematic description of motion to the forces and energy transfers that produce it. Civil-engineering applications combine conservation of mass, mechanical energy, and momentum to analyze pipes, meters, pumps, turbines, gates, bends, jets, spillways, and hydraulic structures.

Three Complementary Conservation Laws

  • Mass: Determines how discharge and velocity are related between sections.
  • Energy: Relates pressure, elevation, velocity, shaft work, and irreversible losses.
  • Momentum: Determines resultant forces and reactions from changes in velocity and pressure.

These equations answer different questions. An energy equation does not directly provide a support reaction, and a momentum equation does not by itself determine head loss.

Euler's Equation Along a Streamline

For steady, inviscid flow, Newton's second law applied tangentially to a differential fluid element gives

dpρ+V dV+g dz=0\frac{dp}{\rho}+V\,dV+g\,dz=0

Integrating for constant density produces Bernoulli's equation. The differential form makes clear that pressure gradients, acceleration, and gravity balance one another along the streamline.

Bernoulli Equation

Mechanical-energy relation for steady, incompressible, inviscid flow along a streamline with no shaft work.

pΞ³+V22g+z=H=constant\frac{p}{\gamma}+\frac{V^2}{2g}+z=H=\text{constant}

Variables

SymbolDescriptionUnit
ppStatic pressurePa
Ξ³\gammaSpecific weightN/m3N/m^3
VVLocal or appropriately averaged flow velocitym/s
zzElevation above a selected datumm
HHTotal mechanical headm

Bernoulli Assumptions

The basic Bernoulli equation requires steady flow, constant density, negligible viscous loss, no pump or turbine work between sections, and application along a streamline. It may be applied across streamlines only when the flow is irrotational. Real engineering systems normally require the extended energy equation.

Interactive Venturi Demonstration

The simulation applies continuity and ideal horizontal Bernoulli flow while using absolute pressure to identify vapor-pressure limits.

Bernoulli's Principle: Venturi Meter

V₁Vβ‚‚P₁Pβ‚‚
Inlet (1)
150.0 kPa abs
2.00 m/s
Throat (2)
120.0 kPa abs
8.00 m/s

This ideal horizontal-flow model neglects losses. Cavitation is assessed using absolute pressure and the approximate vapor pressure of water at 20Β°C. Once vapor pressure is reached, the single-phase Bernoulli model is no longer physically valid.

Kinetic-Energy Correction Coefficient

A cross section rarely has perfectly uniform velocity. Replacing the true velocity distribution with mean velocity requires the coefficient Ξ±\alpha in the kinetic-energy term.

Kinetic-Energy Correction Coefficient

Corrects mean-velocity kinetic energy to equal the actual energy flux.

Ξ±=1AV3∫Au3 dA\alpha=\frac{1}{A V^3}\int_A u^3\,dA

Variables

SymbolDescriptionUnit
Ξ±\alphaKinetic-energy correction coefficient-
uuLocal velocity normal to the sectionm/s
VVSection-average velocitym/s

Typical Values of $\alpha$

For fully developed laminar flow in a circular pipe, Ξ±=2\alpha=2. For ordinary turbulent pipe flow, the flatter profile commonly gives Ξ±\alpha near 1.03 to 1.10. Using Ξ±=1\alpha=1 is often adequate but should be stated, especially near entrances, outlets, contractions, or strongly nonuniform channels.

Extended Mechanical-Energy Equation

Relates two sections through shaft work and irreversible losses.

p1Ξ³+Ξ±1V122g+z1+hpβˆ’ht=p2Ξ³+Ξ±2V222g+z2+hL\frac{p_1}{\gamma} +\alpha_1\frac{V_1^2}{2g} +z_1 +h_p -h_t = \frac{p_2}{\gamma} +\alpha_2\frac{V_2^2}{2g} +z_2 +h_L

Variables

SymbolDescriptionUnit
hph_pHead added by a pumpm
hth_tHead extracted by a turbinem
hLh_LTotal irreversible major and minor head loss from section 1 to 2m
Ξ±1,Ξ±2\alpha_1,\alpha_2Kinetic-energy correction coefficients-

Loss Head Is Directional and Nonnegative

With sections ordered in the actual direction of flow, hLβ‰₯0h_L\ge0. Reversing section labels without changing the assumed flow direction is a common sign error. Pump head is positive when energy is transferred to the fluid; turbine head is positive in the equation above when energy is extracted from the fluid.

Hydraulic Grade Line, HGL

The locus of piezometric head,

HGL=z+pΞ³HGL=z+\frac{p}{\gamma}

In a piezometer open to atmosphere, the liquid level indicates the local HGL for the connected fluid.

Energy Grade Line, EGL

The locus of total mechanical head,

EGL=z+pΞ³+Ξ±V22gEGL=z+\frac{p}{\gamma}+\alpha\frac{V^2}{2g}

The vertical separation between EGL and HGL is the corrected velocity head.

EGL and HGL Behavior

  • Both lines fall in the direction of flow through passive components because of loss.
  • The EGL rises across a pump and falls across a turbine.
  • The HGL can rise through a diffuser even while the EGL falls, because velocity head is converted to pressure head with some loss.
  • Where the HGL drops below the pipe centerline, gage pressure is negative.
  • Absolute pressure must still remain above vapor pressure and above any structural or air-entry limit.

Power Associated with Head

Converts pump, turbine, or loss head to a rate of energy transfer.

P=ρgQH=γQHP=\rho gQH=\gamma QH

Variables

SymbolDescriptionUnit
PPFluid power associated with head HW
QQVolumetric flow ratem3/sm^3/s

Navier-Stokes Equation

For a Newtonian fluid with constant viscosity, the differential linear-momentum equation is

ρ(βˆ‚Vβˆ‚t+(Vβ‹…βˆ‡)V)=βˆ’βˆ‡p+ΞΌβˆ‡2V+ρg\rho\left( \frac{\partial\mathbf{V}}{\partial t} +(\mathbf{V}\cdot\nabla)\mathbf{V} \right) = -\nabla p +\mu\nabla^2\mathbf{V} +\rho\mathbf{g}

The left side is mass density times material acceleration. The right side contains pressure force, viscous diffusion of momentum, and body force. Most practical hydraulic formulas are exact solutions, approximations, or empirical closures derived from this balance.

Turbulent Flow Requires Closure

Time-averaging turbulent Navier-Stokes equations introduces Reynolds stresses that are not known from the mean velocity alone. Friction factors, roughness formulas, turbulence models, and laboratory coefficients supply the required closure information.

Control-Volume Linear Momentum

For a fixed control volume, the sum of external forces equals the rate of momentum accumulation inside plus the net momentum flux through the control surface.

General Linear-Momentum Equation

Integral momentum balance for a fixed control volume.

βˆ‘Fext=ddt∫CVρV dV+∫CSρV(Vβ‹…n) dA\sum\mathbf{F}_{ext} = \frac{d}{dt}\int_{CV}\rho\mathbf{V}\,dV + \int_{CS}\rho\mathbf{V}(\mathbf{V}\cdot\mathbf{n})\,dA

Steady One-Dimensional Momentum

Section-average form for discrete inlets and outlets.

βˆ‘Fext=βˆ‘outβρQVβˆ’βˆ‘inβρQV\sum\mathbf{F}_{ext} = \sum_{out}\beta\rho Q\mathbf{V} - \sum_{in}\beta\rho Q\mathbf{V}

Variables

SymbolDescriptionUnit
Ξ²\betaMomentum correction coefficient-
Fext\mathbf{F}_{ext}All external forces acting on the fluid control volumeN

Momentum Correction Coefficient

Corrects mean-velocity momentum flux to equal the actual flux.

Ξ²=1AV2∫Au2 dA\beta=\frac{1}{A V^2}\int_A u^2\,dA

Typical Values of $\beta$

For fully developed laminar flow in a circular pipe, Ξ²=4/3\beta=4/3. For ordinary turbulent pipe flow, Ξ²\beta is commonly close to 1.01 to 1.05. Do not interchange Ξ±\alpha and Ξ²\beta: one corrects energy flux and the other momentum flux.

Control-Volume Force Analysis

  1. Draw a control volume around the fluid, not around the pipe material.
  2. Choose coordinate axes and assign signed velocity components.
  3. List all forces on the fluid: pressure, weight, wall/support reaction, shear when important, and atmospheric pressure if not canceled by consistent gage pressure.
  4. Apply continuity to obtain missing flow rates or mean velocities.
  5. Apply momentum separately in each coordinate direction.
  6. Solve for the force exerted by the pipe or support on the fluid.
  7. Reverse the sign to obtain the force exerted by the fluid on the pipe, bend, nozzle, gate, or support.

Pressure Forces Act Inward on the Control Surface

Pressure is compressive. At an inlet or outlet, the pressure-force vector acting on the fluid is opposite the outward control-surface normal: Fp=βˆ’pAn\mathbf{F}_p=-pA\mathbf{n}. Memorizing β€œinlet positive, outlet negative” without drawing normals often fails for inclined or multiport systems.

Forces on Bends and Reducers

A bend changes the momentum-vector direction and may also change its magnitude. The support reaction must balance momentum-flux change, pressure forces at every section, fluid weight when significant, and wall shear when it cannot be neglected. Three-dimensional bends require vector components in all relevant directions.

Jets Striking Stationary Surfaces

For a free jet at atmospheric pressure, pressure forces at exposed jet sections cancel in gage form. The force on a stationary plate or vane follows from the change in momentum. A normal flat plate that stops the jet's normal velocity component experiences approximately

F=ρQVF=\rho QV

when splash leaves with negligible normal velocity.

Moving Vanes and Relative Velocity

For a vane moving at speed uu, the mass flow that actually enters the vane is based on relative speed. For a jet of area AA moving in the same direction as the vane,

mΛ™=ρA(Vβˆ’u)\dot m=\rho A(V-u)

Momentum analysis must use the absolute inlet and outlet velocities, while the vane geometry commonly controls relative velocity. Neglecting this distinction produces major force and power errors.

Angular-Momentum Equation

Relates torque to the change in moment of momentum, central to pumps and turbines.

βˆ‘MO=ddt∫CVρ(rΓ—V) dV+∫CSρ(rΓ—V)(Vβ‹…n) dA\sum M_O = \frac{d}{dt}\int_{CV}\rho(\mathbf{r}\times\mathbf{V})\,dV + \int_{CS}\rho(\mathbf{r}\times\mathbf{V})(\mathbf{V}\cdot\mathbf{n})\,dA

For steady one-inlet/one-outlet turbomachinery,

T=mΛ™(r2VΞΈ2βˆ’r1VΞΈ1)T=\dot m(r_2V_{\theta2}-r_1V_{\theta1})

Variables

SymbolDescriptionUnit
TTTorque exerted on the fluidNβ‹…mN\cdot m
VΞΈV_\thetaTangential or whirl component of absolute velocitym/s

Torricelli, Pitot, and Venturi Applications

  • Torricelli: For a large reservoir discharging ideally to atmosphere, Vβ‰ˆ2gHV\approx\sqrt{2gH}.
  • Pitot tube: Stagnating a streamline converts velocity head to pressure head, so V=C2gΞ”hV=C\sqrt{2g\Delta h} after calibration.
  • Venturi meter: Continuity and pressure-head difference determine discharge, modified by a discharge coefficient.

Each device inherits Bernoulli assumptions and requires appropriate coefficients, pressure references, tap elevations, and fluid-column conversions.

Cavitation

Formation and collapse of vapor cavities when local absolute pressure falls to or below the liquid vapor pressure at the operating temperature.

Use Absolute Pressure for Cavitation

Gage pressure cannot be compared directly with vapor pressure unless atmospheric pressure is first added. Once vapor forms, a single-phase incompressible Bernoulli calculation no longer represents the actual flow. Cavitation can cause noise, vibration, erosion, loss of capacity, and unstable operation.

Energy–Momentum Problem Strategy

  1. Establish the actual flow direction and pressure reference: absolute or gage.
  2. Apply continuity and calculate section-average velocities.
  3. Use the extended energy equation to find head, pressure, discharge, losses, or shaft work.
  4. Construct HGL/EGL checks and confirm that absolute pressure remains physically acceptable.
  5. Use momentum, not energy alone, to obtain reactions or forces.
  6. Include Ξ±\alpha and Ξ²\beta when profiles are strongly nonuniform or laminar.
  7. Check units, signs, and whether empirical coefficients match the device and flow regime.
Key Takeaways
  • Bernoulli's equation is an inviscid special case; real systems use the extended energy equation with pumps, turbines, losses, and often Ξ±\alpha.
  • HGL represents piezometric head, while EGL includes corrected velocity head.
  • Momentum analysis requires a clearly drawn fluid control volume, inward pressure forces, vector velocity components, and a final action–reaction sign reversal for structural loads.
  • Ξ±\alpha corrects kinetic-energy flux and Ξ²\beta corrects momentum flux; they are not interchangeable.
  • Navier-Stokes provides the local momentum balance, while engineering friction factors and coefficients supply practical closures.
  • Cavitation checks always use absolute pressure and vapor pressure at the actual temperature.