CELE Hydraulics & Fluid Mechanics — Fundamentals of Fluid FlowSummary

For steady, incompressible flow through a stream tube, the volume flow rate Q (measured in m³/s) remains constant along the flow path. This is expressed as Q = A₁v₁ = A₂v₂, where A is cross-sectional area (m²) and v is velocity (m/s). The equation reflects the physical principle that mass cannot be created or destroyed: if a pipe narrows, water must accelerate to maintain the same discharge. For a circular pipe, the velocity relationship simplifies to v₂ = v₁(D₁/D₂)², showing that velocity is inversely proportional to the square of the diameter. This principle is fundamental to all flow problems and is the first equation applied in any hydraulic analysis.

The total mechanical energy per unit weight (total head H, in meters) at any point in a flowing fluid comprises three components: pressure head p/γ (where p is absolute pressure in Pa and γ is specific weight = 9.81 kN/m³ for water), velocity head v²/2g (where g = 9.81 m/s²), and elevation head z above a reference datum. Between two points, with a pump adding head h_A, a turbine extracting head h_E, and friction removing head h_L, the equation is: p₁/γ + v₁²/2g + z₁ + h_A = p₂/γ + v₂²/2g + z₂ + h_E + h_L. When no machines and no losses are present (ideal case), the total head is conserved: H₁ = H₂. The left side represents available energy input; the right side, energy output plus losses. This is the master equation for pipe and channel flow design.

The Energy Grade Line is a graphical representation of total head H = p/γ + v²/2g + z plotted along the flow path. It always slopes downward in the direction of flow (except where a pump adds energy). The Hydraulic Grade Line is the same plot but with velocity head subtracted: HGL = p/γ + z. In an open channel, the HGL coincides with the free water surface. In a closed pipe under pressure, the HGL lies below the EGL by exactly the velocity head. At a large-diameter section (low velocity), the HGL is relatively high; at a narrow section (high velocity), the HGL drops. The vertical distance between EGL and HGL at any point equals the velocity head v²/2g. These lines provide immediate visual insight into pressure distribution and identify low-pressure zones (cavitation risk) or high-pressure zones (burst risk) in a piping system.

The net force acting on a fluid control volume equals the rate of change of momentum: ΣF = ρQ(v₂ − v₁), where ρ is fluid density (kg/m³, typically 1000 for water), Q is discharge (m³/s), and v₂ − v₁ is the change in velocity (m/s). Applied component-wise (x, y, z directions), this equation determines reaction forces on pipe elbows, tees, nozzles, and guide vanes. For example, a 90° elbow with constant velocity magnitude but changing direction experiences a reaction force perpendicular to both inlet and outlet flows. The momentum equation is derived from Newton's second law (F = ma) applied to a control mass of fluid, making it a direct consequence of fundamental mechanics. It is independent of viscous losses—the force depends only on mass flow rate and velocity change, not on head loss.

The power available in a flowing stream of water is P = γQH (in watts), where γ = 9.81 kN/m³ is the specific weight of water, Q is discharge (m³/s), and H is the head (m) across which energy is available. A pump requires input power P_input = γQH/η to raise water a height H against resistance, where η (pump efficiency, 0 < η ≤ 1) accounts for mechanical and fluid losses. A turbine produces output power P_output = η·γQH from a head H, the product being mechanical power available at the shaft. For example, a pump lifting 0.1 m³/s through 30 m head at 85% efficiency requires P_input = (9.81 × 0.1 × 30)/0.85 ≈ 34.6 kW of electrical input. Conversely, a turbine with the same discharge and head at 88% efficiency produces P_output = 0.88 × 9.81 × 0.1 × 30 ≈ 26.0 kW. Power calculations are direct, but the sign and placement of efficiency (denominator for pump input, numerator for turbine output) must be correct.

Head loss h_L (in meters of water column) represents the irreversible conversion of mechanical energy into thermal energy due to friction between fluid layers and between fluid and pipe wall. Total head loss comprises major losses (friction in straight pipe) and minor losses (at entrance, exit, elbows, valves, etc.). The major loss is calculated using the Darcy-Weisbach equation: h_f = f(L/D)(v²/2g), where f is the friction factor (dimensionless, depends on Reynolds number and pipe roughness), L is pipe length (m), D is diameter (m), and v is velocity (m/s). The friction factor is determined from the Moody diagram or Colebrook equation. Minor losses are expressed as h_m = K(v²/2g), where K is a minor-loss coefficient (K ≈ 0.5 for pipe entrance, 1.0 for exit, 0.9 for 90° elbow, etc.). Total loss: h_L = h_f + Σh_m. Head loss always opposes the flow and is always subtracted in the Bernoulli equation.

Absolute pressure is measured from perfect vacuum (zero reference); gauge pressure is measured from atmospheric pressure (101.325 kPa). In the Bernoulli equation, pressure must be absolute. At the free surface of an open channel, gauge pressure is zero but absolute pressure is 101.325 kPa. In a closed pipe, if a gauge pressure of +50 kPa is measured, the absolute pressure is 101.325 + 50 = 151.325 kPa. Conversely, a negative gauge pressure (vacuum or partial vacuum) of −20 kPa means absolute pressure is 101.325 − 20 = 81.325 kPa. This distinction is critical for identifying cavitation risk: cavitation occurs when absolute pressure falls below the vapor pressure of water (≈2.3 kPa at 20°C). A gauge pressure of −95 kPa would produce an absolute pressure of 6.325 kPa, risking cavitation. Confusion between gauge and absolute pressure is a leading source of error in student calculations.

Continuity equation determines velocity at each point (Q = Av). Higher velocity means higher velocity head (v²/2g) and higher friction losses (∝ v²). Thus, a pipe reducer accelerates the flow, raising velocity head and friction losses. This justifies energy loss in the Bernoulli equation and sets boundary conditions for momentum calculations.

The Bernoulli equation with h_A (pump head) and losses allows calculation of required pump head: h_A = (p₂ − p₁)/γ + (v₂² − v₁²)/2g + (z₂ − z₁) + h_L. Once h_A is known, power is P = γQh_A/η, which guides pump selection. Conversely, a chosen pump's characteristic curve constrains the system operating point. Turbine head h_E is similarly derived and power is P = η·γQh_E.

The momentum equation gives the reaction force on a pipe bend: F = ρQ(v₂ − v₁). This force must be resisted by anchors, guide blocks, or tie-rods in the supporting structure. For a 90° bend at constant speed (|v₂| = |v₁| = v), the resultant force is F = ρQv√2, directed at 45° away from the pipe axis. Knowledge of this force is essential for structural design of pipe supports and hydraulic machines (nozzle reaction force, turbine blade loading).

The rate of head loss per unit length is dh_L/dx = f(v²/2gD). A steep slope on the EGL indicates high losses; a gentle slope indicates smooth, efficient flow. The HGL, lying below EGL by velocity head, also slopes downward but by a smaller amount if velocity is constant. Changes in diameter (and thus velocity) produce 'jumps' in the HGL (positive at enlargements, negative at contractions) while the EGL smoothly reflects energy loss accumulation.

At high-velocity sections (small pipes or nozzles), velocity head v²/2g is large. If this occurs in a region of low static pressure p/γ, the total head may still be adequate, but the absolute pressure p_abs = p + p_atm can fall dangerously low. When p_abs < p_vapor (≈2.3 kPa), cavitation bubbles form, causing erosion and performance loss. Bernoulli identifies low-pressure points; continuity identifies where velocity (and thus velocity head) is high; absolute pressure assessment determines cavitation onset.

A nozzle converts pressure head into velocity head. Continuity shows that area reduction causes v₂ = v₁(A₁/A₂). Bernoulli (frictionless) gives the pressure drop: p₁/γ + v₁²/2g = p₂/γ + v₂²/2g, so p₂ < p₁. Momentum shows the reaction force: F = ρQ(v₂ − v₁), directed opposite to the jet. All three equations work in concert: continuity sets the geometry, Bernoulli predicts the pressure and velocity, and momentum predicts the reaction force.

A municipal water supply system from a reservoir to a residential district must deliver Q = 0.25 m³/s over a distance of 5 km with elevation rise of 80 m. Using continuity, the pipe diameter is selected to keep velocity below 1.2 m/s (to minimize losses). Using Bernoulli with estimated head loss (Darcy-Weisbach), the required pump head is computed: h_pump = Δz + Δ(p/γ) + Δ(v²/2g) + h_L. For typical conditions, h_pump might be 85–95 m. With efficiency η = 0.82, the pump power is P = 9.81 × 0.25 × 90 / 0.82 ≈ 27 kW. This determines the pump motor rating and operating cost. Head loss along the line is also mapped using the HGL, ensuring that pressure does not drop below safe minimum (e.g., 10 m gauge) at any point to prevent cavitation and maintain service.

In Philippine agricultural regions, an irrigation network distributes water from a central pump station to multiple field laterals. Main-line velocity is set at 0.5–1.0 m/s to minimize loss and pumping cost. At branch points, continuity is applied to divide discharge among laterals. Friction head loss in the main line (often 2–5 km long) is calculated; minor losses at tees and laterals are added. The hydraulic gradient (HGL) is plotted: it must remain above minimum pressure (typically 5–10 m) to ensure adequate pressure at the most remote field inlet. If the HGL drops too low, a booster pump is installed. Power consumption for the main pump is P = γQH/η; annual energy cost drives the pipe-sizing decision: larger diameter → higher initial cost but lower pumping cost.

A run-of-river hydropower project exploits a river drop (head H, e.g., 50 m) and discharge (Q, e.g., 8 m³/s). The available power is P_available = γQH = 9.81 × 8 × 50 ≈ 3.93 MW. The turbine efficiency (typically 88–92% for modern machines) gives P_output = 0.90 × 3.93 ≈ 3.54 MW. The penstock (pressure pipe) is sized using continuity: if v ≤ 2.5 m/s is desired, A = Q/v = 8/2.5 = 3.2 m², giving D ≈ 2.02 m. Friction loss in the penstock (estimated 1–2 m for a 500 m length) is subtracted from H in the Bernoulli equation when calculating actual turbine head. Force on turbine blades (momentum equation) determines blade shape and stresses. Cavitation risk at the turbine exit is assessed by ensuring absolute pressure remains above vapor pressure. These calculations govern both feasibility and final design of Philippine hydropower sites.

A high-velocity water main (v = 3 m/s, diameter D = 600 mm, discharge Q ≈ 0.85 m³/s) encounters a 90° elbow. By continuity, velocity remains 3 m/s (same pipe size). By momentum, the change in momentum vector (from entering along x-axis to exiting along y-axis) produces a resultant reaction force. For equal velocity magnitudes at 90° angle: F_resultant = ρQv√(2) = 1000 × 0.85 × 3 × √2 ≈ 3,606 N, directed at 45° to both arms. The elbow support must be designed to resist this force via anchor blocks, guide rods, or tie-rods. Neglecting this force leads to pipe movement, joint separation, or support failure—a common failure mode in aging Philippine water systems. Momentum calculations are thus essential for infrastructure durability.

At a pump intake on a high-altitude or deep-suction application, absolute pressure may approach vapor pressure, risking cavitation. For example, at 1,500 m elevation, atmospheric pressure ≈ 84.5 kPa. A suction lift of 5 m (intake 5 m below pump center) reduces absolute pressure further: p_abs ≈ 84.5 − 9.81 × 5 ≈ 35.1 kPa. If the pump inlet has a velocity v ≈ 2 m/s (from continuity: A = Q/v), velocity head ≈ 0.2 m, which consumes another 1.96 kPa. The absolute pressure available is thus ~33 kPa. If this falls below vapor pressure (~2.3 kPa), cavitation occurs. Prevention measures include: (1) submerging the intake to increase p_atm contribution, (2) using a bell-mouth inlet to reduce velocity head (lower K value), (3) reducing suction lift, or (4) choosing a low-NPSHR (net positive suction head required) pump. Bernoulli and absolute pressure analysis guide all these design decisions.

During heavy rains, a flood channel must safely pass a design discharge (e.g., Q = 500 m³/s from a 100-year storm). The channel is trapezoidal with bed slope S_0 = 0.002. Using continuity with uniform flow assumptions, the normal depth is found (typically iterative, using Manning equation). At the normal depth, velocity v ≈ 1.8 m/s and hydraulic radius R_h ≈ 8 m. The HGL is essentially the water surface elevation, which must remain below the channel crown (freeboard requirement). If an obstruction (e.g., bridge) locally narrows the channel, continuity demands higher velocity at that section; Bernoulli predicts a rise in water level (backwater) upstream and a drawdown downstream. EGL and HGL plots show whether flooding will occur and where mitigation (wider section, raise crown) is needed. In Philippine rivers prone to flooding, such analysis is essential for public safety and infrastructure protection.

A 300 mm supply pipe (v₁ = 2 m/s, p₁ = 200 kPa gauge) narrows to a 150 mm nozzle. By continuity: v₂ = v₁(D₁/D₂)² = 2 × 4 = 8 m/s. By Bernoulli (neglecting loss): p₂/γ + v₂²/2g = p₁/γ + v₁²/2g, so p₂/9.81 + 8²/19.62 = 200/9.81 + 2²/19.62, giving p₂ ≈ 121 kPa gauge. The pressure drops from 200 kPa to 121 kPa as kinetic energy increases. By momentum: F = ρQ(v₂ − v₁) acts as a reaction force pulling backward on the nozzle. Knowledge of these quantities is essential for: (1) pressure rating of the nozzle (121 kPa is safe; ~20 MPa is a fail point), (2) nozzle support design (reaction force), and (3) cavitation risk (p₂ is still safe, but in a smaller or longer nozzle, p₂ might drop below vapor pressure). This example ties all three principles together.