Liquid Pipeline Hydraulics Fundamentals

1. Fundamentals of Liquid Flow

Liquid pipeline hydraulics is governed by the principles of fluid mechanics applied to closed conduit flow. Unlike gas pipelines, where compressibility effects dominate, liquids are essentially incompressible, which simplifies the governing equations while introducing other challenges such as surge (water hammer) and cavitation.

Reynolds Number & Flow Regimes

The Reynolds number (Re) is the dimensionless ratio of inertial forces to viscous forces, and it determines the flow regime in the pipe:

Where:

• ρ = fluid density (lb/ft³ or kg/m³)
• V = average flow velocity (ft/s or m/s)
• D = pipe inside diameter (ft or m)
• μ = dynamic (absolute) viscosity (lb/(ft·s) or Pa·s)
• ν = kinematic viscosity (μ/ρ) (ft²/s or m²/s)

The Reynolds number defines three flow regimes:

Most liquid pipeline applications operate in the fully turbulent regime, with Reynolds numbers ranging from 10,000 to several million. Heavy crude oil pipelines and high-viscosity product lines may operate in laminar or transitional flow, requiring different friction factor correlations.

Continuity Equation & Velocity

For steady-state, incompressible flow the continuity equation reduces to:

Where Q is the volumetric flow rate, A is the cross-sectional flow area, and V is the mean velocity. This relationship is used to convert between flow rate (in BPD, GPM, or m³/h) and velocity (in ft/s or m/s) for a given pipe diameter. In US customary units with Q in gallons per minute and D in inches:

Friction Factor Correlations

The friction factor (f) is the dimensionless parameter that quantifies energy loss due to pipe wall shear. Several correlations are available:

• Colebrook-White (1939): The implicit equation that defines the Moody diagram. Requires iterative solution: 1/√f = −2 log(ε/3.7D + 2.51/(Re√f)). Accurate across the entire turbulent regime.
• Swamee-Jain (1976): An explicit approximation accurate to within ±1% of Colebrook-White for 5,000 ≤ Re ≤ 10&sup8; and 10²ε/D ≤ 0.05: f = 0.25 / [log(ε/3.7D + 5.74/Re&sup0;·&sup9;)]². Preferred for spreadsheet and calculator use.
• Churchill (1977): A single explicit equation valid for laminar, transitional, and turbulent flow. Useful for automated calculations that must handle all regimes without branching logic.
• Laminar flow: For Re < 2,100, the Hagen-Poiseuille solution gives f = 64/Re exactly, independent of pipe roughness.

Pipe Roughness Values

The absolute roughness (ε) represents the average height of surface irregularities on the pipe interior wall. It varies by pipe material and condition:

2. Darcy-Weisbach Equation

The Darcy-Weisbach equation is the most rigorous and universally applicable method for calculating frictional pressure drop in pipe flow. It is valid for any Newtonian fluid, any flow regime (laminar or turbulent), and any pipe material.

Pressure Drop Form

Where:

• ΔP = frictional pressure drop (psi)
• f = Darcy-Weisbach friction factor (dimensionless; note: some references use the Fanning friction factor, which is f/4)
• L = pipe length (ft)
• D = pipe inside diameter (ft)
• ρ = fluid density (lb/ft³)
• V = average flow velocity (ft/s)
• g_c = gravitational constant = 32.174 lb·ft/(lbf·s²)
• 144 = conversion factor (in²/ft²) to convert from lbf/ft² to psi

Head Loss Form

The same equation expressed as head loss (energy per unit weight of fluid) is often more convenient for pump system calculations:

Where h_f is the friction head loss in feet of liquid, and g = 32.174 ft/s² is the acceleration due to gravity. To convert head loss to pressure drop: ΔP = ρgh_f/144 (psi), or simply ΔP = h_f × SG / 2.31 where SG is the specific gravity.

When to Use Darcy-Weisbach

• Any liquid: Water, crude oil, NGL, refined products, glycol, amine, brine — all Newtonian fluids.
• All flow regimes: Works for laminar (f = 64/Re), transitional, and fully turbulent flow.
• Any pipe material: The friction factor accounts for roughness through the Colebrook-White equation or its approximations.
• Elevated temperatures and pressures: Valid as long as fluid properties (density, viscosity) are evaluated at the correct operating conditions.

Advantages

• Theoretically rigorous — derived from dimensional analysis and validated against extensive experimental data.
• Applicable to all fluids and all flow regimes without empirical restrictions.
• Separates the friction factor from the geometry and velocity terms, making it easy to evaluate sensitivities.

Limitations

• Requires iterative calculation of f when using the Colebrook-White equation (though explicit approximations like Swamee-Jain eliminate this).
• Requires accurate fluid properties (density and viscosity) at operating temperature and pressure.
• Does not directly account for minor losses (fittings, valves); these are added separately using equivalent lengths or K-factors per Crane TP-410.

3. Hazen-Williams Method

The Hazen-Williams equation is a widely used empirical formula for calculating head loss in water distribution and fire protection piping. Its simplicity and the availability of tabulated C-factors make it popular for systems carrying water or water-like fluids at ambient temperatures.

The Equation

In US customary units with head loss per unit length:

Where:

• h_f = friction head loss (ft)
• L = pipe length (ft)
• Q = volumetric flow rate (GPM)
• C = Hazen-Williams C-factor (dimensionless roughness coefficient)
• d = pipe inside diameter (inches)

Unlike the Darcy-Weisbach friction factor, the C-factor is not dimensionless in the strict sense — the equation is empirical, and the coefficient 10.67 absorbs the unit conversions for the specific set of units listed above.

C-Factor Table by Pipe Material & Age

A higher C-factor indicates a smoother pipe with lower friction loss. The C-factor decreases over time as the pipe interior corrodes, scales, or accumulates deposits. Use the aged C-factor for design to account for long-term performance degradation.

When to Use Hazen-Williams

• Water and water-like fluids: The equation was empirically derived from water flow data and is accurate only for fluids with kinematic viscosity near that of water (approximately 1.0 cSt at 68°F).
• Turbulent flow only: The C-factor does not account for laminar flow. Hazen-Williams should only be applied when Re > 4,000.
• Temperature range: Best accuracy between 40°F and 75°F (4°C to 24°C). Outside this range, viscosity changes that are not captured by the constant C-factor can introduce significant errors.
• Common applications: Municipal water distribution, fire protection systems, cooling water piping, and irrigation systems.

Comparison with Darcy-Weisbach

4. Pipe Sizing & Selection

Pipe sizing for liquid pipelines balances three competing requirements: the allowable pressure drop budget, the maximum and minimum velocity constraints, and the capital cost of the pipe. The designer must find the smallest standard pipe size that satisfies all three constraints simultaneously.

Three Calculation Modes

• Pressure drop calculation: Given the pipe diameter, flow rate, and fluid properties, calculate the friction pressure drop and total head loss along the pipeline. This is the most common mode for evaluating an existing or proposed pipe size.
• Flow capacity calculation: Given the pipe diameter and available pressure differential (or pump head), calculate the maximum flow rate the pipeline can deliver. Used for debottlenecking studies and pump selection.
• Required diameter calculation: Given the flow rate and allowable pressure drop, solve for the minimum pipe inside diameter. Round up to the next standard NPS (Nominal Pipe Size) per ASME B36.10.

Erosional Velocity Limits per API RP 14E

API RP 14E provides the erosional velocity formula for mixed-phase flow, but it is also applied to single-phase liquids as an upper velocity bound:

Where V_e is the erosional velocity (ft/s), C is an empirical constant (typically 100 for continuous service, up to 150 for intermittent service), and ρ is the fluid density (lb/ft³). For clean, non-corrosive liquid service in carbon steel pipe, the erosional velocity provides a conservative upper limit that prevents accelerated wall thinning from fluid impingement.

Velocity Guidelines for Liquid Pipelines

Minimum velocity constraints ensure that solids (sand, scale, wax) remain in suspension and do not settle in the pipe. Maximum velocity constraints prevent erosion, excessive vibration, and water hammer severity.

NPS Selection Process

After calculating the minimum required inside diameter, select the next standard NPS from ASME B36.10 pipe schedules. The process is:

• Calculate the theoretical minimum inside diameter from the flow rate and maximum allowable velocity.
• Select the next standard NPS whose inside diameter (based on the specified wall thickness schedule) equals or exceeds the calculated minimum.
• Verify that the pressure drop at the selected NPS does not exceed the available differential pressure.
• Verify that the velocity at the selected NPS falls within the acceptable range (not below the minimum and not above the erosional limit).
• If constraints conflict, iterate up to the next NPS until all requirements are met.

Economic Diameter Considerations

The economic pipe diameter balances the capital cost of the pipe (which increases with diameter) against the operating cost of pumping (which decreases with diameter as friction losses drop). For a given flow rate, a larger pipe has lower pressure drop and requires less pump power, but costs more to purchase, transport, and install. The optimum is typically found by performing a net present value (NPV) analysis over the pipeline’s operating life, comparing total installed cost plus discounted pumping energy cost for each candidate pipe size.

5. Practical Design Considerations

Beyond the friction loss equations, several real-world factors must be incorporated into a liquid pipeline hydraulic design to ensure reliable and safe operation.

Elevation Effects & Static Head

Unlike gas pipelines, liquid pipelines are strongly affected by elevation changes. The total pressure at any point in the pipeline is the sum of the friction pressure loss and the static (elevation) head:

Where the elevation pressure change is: ΔP_elevation = 0.433 × SG × Δh (psi), with Δh in feet and SG as the fluid specific gravity. An increase in elevation adds to the required pump discharge pressure, while a decrease in elevation reduces it. On hilly terrain, the hydraulic gradient line must be plotted to verify that the pipeline pressure remains above the fluid vapor pressure at every high point to prevent vapor lock or column separation.

Temperature Effects on Viscosity & Density

Liquid viscosity is strongly temperature-dependent. For crude oil and heavy hydrocarbons, viscosity can decrease by an order of magnitude with a 50°F temperature increase. This dramatically affects the friction factor and pressure drop. Design calculations should evaluate pressure drop at both the coldest expected operating temperature (maximum viscosity, highest friction) and the warmest (minimum viscosity, lowest friction) to bound the operating envelope.

Density also varies with temperature, though less dramatically. For petroleum products, the density correction follows the API gravity tables in API MPMS Chapter 11 (formerly API 2540). For water, density is approximately 62.4 lb/ft³ at 60°F and decreases to about 60.0 lb/ft³ at 200°F.

Multi-Phase Considerations

When free gas is present in a liquid pipeline (as may occur at pressure drops below the bubble point, or in well flowlines with gas-liquid mixtures), single-phase liquid hydraulic methods no longer apply. Multi-phase flow requires specialized correlations such as Beggs-Brill, Dukler, or mechanistic models that account for flow pattern transitions (bubble, slug, annular, stratified). If any segment of a nominally liquid pipeline may experience two-phase flow, use a multi-phase flow simulator for that segment rather than the single-phase Darcy-Weisbach equation.

Standards & References

• Crane TP-410: Flow of Fluids Through Valves, Fittings, and Pipe — the primary reference for single-phase liquid and gas friction loss calculations, equivalent lengths, and K-factors
• API RP 14E: Recommended Practice for Design and Installation of Offshore Production Platform Piping Systems — erosional velocity formula and velocity guidelines
• ASME B31.4: Pipeline Transportation Systems for Liquids and Slurries — design pressure, wall thickness, and testing requirements for liquid pipelines
• ASME B36.10: Welded and Seamless Wrought Steel Pipe — standard NPS dimensions, wall thickness schedules
• 49 CFR 195: Transportation of Hazardous Liquids by Pipeline — federal safety regulations for liquid pipelines