parallel-pipe-systems

name	parallel-pipe-systems
description	Parallel pipe network analysis — Hardy Cross method, Newton-Raphson network solver, flow distribution in parallel branches, head loss balancing (Darcy-Weisbach, Hazen-Williams), looped networks, pipe network matrix formulation, pump-system curve intersection, flow control (balancing valves, orifice plates), manifold design, and AWWA/ASME pipe network standards.
metadata	{"priority":7,"promptSignals":{"phrases":["parallel pipe","pipe network","Hardy Cross","pipe flow distribution","loop pipe network","pipe manifold"],"minScore":3}}

Parallel Pipe Network Analysis — Complete Skill

Governing Equations

Head Loss in a Single Pipe

Darcy-Weisbach: h_f = f × (L/D) × V²/(2g) = f × (L/D) × Q²/(2g × A²) [m; f = Darcy friction factor; L = length; D = diameter; Q = flow; A = πD²/4]

Resistance coefficient form: h_f = K × Q² [K = resistance coefficient; K = f×L/(D×2g×A²) for a pipe] h_f = R × Qⁿ [general; n = 2 for Darcy-Weisbach; n = 1.852 for Hazen-Williams]

Hazen-Williams (for water distribution, D > 50 mm): V = 0.8492 × C_HW × R_h^0.63 × S^0.54 [m/s; C_HW = Hazen-Williams C; R_h = hydraulic radius = D/4; S = h_f/L] h_f = 10.67 × L × Q^1.852 / (C_HW^1.852 × D^4.87) [m; Q in m³/s]

Hazen-Williams C values: New smooth steel pipe: 140; Old steel with minor rust: 100–120; Cast iron (old): 80–100; PVC: 140–150; Concrete: 100–130

Colebrook-White (friction factor for turbulent flow): 1/√f = -2.0 × log₁₀(ε/(3.7D) + 2.51/(Re×√f)) [implicit; solve iteratively or use Swamee-Jain] Swamee-Jain (explicit): f = 0.25 / [log₁₀(ε/(3.7D) + 5.74/Re^0.9)]² [error < 3% for 10⁻⁶ ≤ ε/D ≤ 10⁻² and 5×10³ ≤ Re ≤ 10⁸]

Parallel Pipe Network Fundamentals

Conservation Laws

Node (junction) continuity — Kirchhoff's Current Law: Σ Q_in = Σ Q_out [net flow into any node = 0]

Loop head balance — Kirchhoff's Voltage Law: Σ h_f = 0 around any closed loop [sum of head losses around any loop = 0; sign convention: CW positive]

These two constraints fully determine flow distribution in a network.

Simple Parallel System (Two Branches)

Conditions for parallel branches: Same head loss across each branch: h₁ = h₂ [both connect same two nodes A and B] Total flow: Q_total = Q₁ + Q₂

Solution: h_f = K₁ × Q₁² = K₂ × Q₂² Q₁/Q₂ = √(K₂/K₁) Q₁ = Q_total × √(K₂) / (√(K₁) + √(K₂)) Q₂ = Q_total × √(K₁) / (√(K₁) + √(K₂))

Equivalent pipe resistance for parallel branches: 1/√K_eq = 1/√K₁ + 1/√K₂ + ... + 1/√Kₙ [for n branches in parallel; Darcy-Weisbach h = KQ²]

Example (two pipes in parallel): Pipe 1: D₁ = 100 mm, L₁ = 200 m, f₁ = 0.025 → K₁ = 0.025×200/(0.1×2×9.81×(π×0.01/4)²) = 4,315 Pipe 2: D₂ = 150 mm, L₂ = 200 m, f₂ = 0.022 → K₂ = 0.022×200/(0.15×2×9.81×(π×0.0225/4)²) = 498 Q₁/Q₂ = √(498/4315) = 0.340; if Q_total = 0.02 m³/s: Q₁ = 0.005 m³/s, Q₂ = 0.015 m³/s h_f = K₁×Q₁² = 4315×0.005² = 0.108 m (check: K₂×Q₂² = 498×0.015² = 0.112 m ≈ ✓)

Hardy Cross Method

Iterative Solution for Looped Networks

Algorithm (Hardy Cross):

1. Assume initial flows Q_i satisfying continuity at all nodes
2. For each loop L, compute correction ΔQ: ΔQ_L = -Σ(K_i × Q_i × |Q_i|^(n-1)) / (n × Σ(K_i × |Q_i|^(n-1))) For Darcy-Weisbach (n=2): ΔQ_L = -Σ(K_i × Q_i²) / (2 × Σ(K_i × |Q_i|)) For Hazen-Williams (n=1.852): ΔQ_L = -Σ(h_f,i) / (1.852 × Σ(h_f,i/|Q_i|))
3. Apply ΔQ to all pipes in loop L (+ for CW convention; − for CCW pipes)
4. Update Q for all pipes; a pipe in two loops gets two corrections
5. Repeat until |ΔQ| < tolerance (0.001 m³/s typical)

Convergence: Hardy Cross converges linearly; 5–20 iterations typical for simple networks

Manual Hardy Cross table: Loop → each pipe → Q_assumed, h_f = K×Q|Q|, h_f/|Q|, → ΔQ → Q_corrected

Newton-Raphson Network Method (Matrix Formulation)

Nodal formulation (preferred for computer solution): Unknowns: node pressures (heads) H_j at each junction Pipe flow: Q_ij = sgn(H_i-H_j) × |H_i-H_j|^(1/n) / K_ij^(1/n) [from h = KQⁿ; n=2 DW; n=1.852 HW]

Residual equation at each node: R_j = Σ Q_ij (inflows positive) - q_j = 0 [q_j = external demand at node j]

Jacobian matrix assembly: ∂R_j/∂H_j = Σ (1/n) × |Q_ij| / |H_i-H_j| [diagonal term; sum over all pipes at node j] ∂R_j/∂H_i = -(1/n) × |Q_ij| / |H_i-H_j| [off-diagonal; from pipe i-j]

Newton-Raphson update: J × ΔH = -R [solve linear system; J = Jacobian matrix; ΔH = head correction vector] H_new = H + ΔH Converges quadratically; 4–6 iterations for typical networks

Manifold Design

Supply Manifold (One-to-Many)

Goal: equal flow through each branch (uniform distribution) Challenge: pressure decreases along header due to friction → downstream branches flow less

Bernoulli along header: P₁/ρg + V₁²/2g = P₂/ρg + V₂²/2g + h_f,1-2 [pressure drops along header] As flow is extracted: header velocity decreases → some pressure recovery

Design approaches:

1. Tapered header: diameter decreasing along header to maintain ~constant pressure
2. Flow control orifice at each branch: orifice controls flow independent of branch pressure
3. Manifold chamber (plenum): large diameter header → velocity pressure negligible → uniform static pressure

Manifold uniformity criterion: For uniform distribution within ±5%: ΔP_branch/ΔP_header ≥ 10 (orifice dominates) Rule of thumb: A_branch ≤ 0.05 × A_header per branch

Distribution coefficient: If all branches identical: Q_branch = Q_total / N_branches (ideal) Non-uniformity factor U = (Q_max - Q_min) / Q_mean × 100% [target: U < 5–10%]

Pump-System Interaction

System Curve

System curve H_system vs. Q: H_system = H_static + K_system × Q² [H_static = elevation difference; K_system = Σ pipe resistances] System curve: parabola passing through (0, H_static)

Parallel pump operation: Two identical pumps in parallel: same head H, flow doubles (Q_total = 2Q at same H) Combined pump curve: Q_combined(H) = 2 × Q_single(H) [add flows at same head] Operating point: intersection of combined pump curve with system curve

Parallel pumps — different pumps: For pumps A and B in parallel: cannot add flows at same H if H_shut-off differs If H_shut,A < H_shut,B: pump A doesn't contribute at high system head Selection: pumps for parallel operation should have similar shut-off heads

Series pumps: Combined head: H_combined(Q) = H_A(Q) + H_B(Q) [add heads at same flow] Series operation: higher head, same flow as single pump at system curve

Balancing Valves and Flow Control

Proportional Balancing (HVAC Networks)

Manual balancing valve: Resistance can be adjusted; set to achieve design flow at each terminal Procedure: start from index circuit (highest resistance); balance all other branches to match Balancing: K_valve adjusted until Q_measured = Q_design at each branch

Automatic balancing valves: Differential pressure independent flow controller: maintains flow regardless of network pressure variations Application: HVAC distribution; each coil gets exactly design flow

Orifice plate for flow restriction: β = D_orifice/D_pipe; Cd = 0.61–0.65 Q = Cd × A_orifice × √(2×ΔP/ρ) For a required ΔP_orifice: A_orifice = Q / (Cd × √(2×ΔP/ρ))

Standards and References

Standard	Scope
AWWA M32	Computer Modeling of Water Distribution Systems
ASME B36.10M	Welded and seamless wrought steel pipe dimensions
ASME B31.3	Process piping pressure and flow design
NFPA 13	Sprinkler system (parallel pipe network design)
ISO 4064	Measurement of water flow in distribution systems
Swamee-Jain (1976)	Explicit Colebrook-White friction factor approximation
Epanet 2.2 (US EPA)	Public-domain pipe network simulation software

Output

Provide: network description (nodes, pipes, loops, external demands [m³/s] and supply heads [m]), pipe data table (L [m], D [mm], ε [mm], K-value [m/(m³/s)²]), Hardy Cross solution (iteration table: Q_assumed, h_f, h_f/Q, ΔQ for each loop; converged Q [m³/s] per pipe), node pressures H_j [m] (from Newton-Raphson or derived from flows), parallel branch flow split (Q₁/Q_total [%]; head loss check h₁ = h₂ [m]), manifold uniformity (U [%]; branch-to-header area ratio; orifice ΔP if used), pump operating point (Q_op [m³/s]; H_op [m]; efficiency η_pump [%]; power P [kW]), balancing valve recommendations (branches requiring throttling; target ΔP_valve [kPa]), pressure drop summary (pipe friction [m], minor losses [m], control valves [m]; total system ΔP [m] or [kPa]), and applicable standard (AWWA M32, ASME B31.3, Epanet 2.2).