Water Distribution Network: Multi-Reservoir Model
| Pipe | Res Elev (m) | Diameter (mm) | Length (m) | Direction | Flow (L/s) | Velocity (m/s) | Head Loss (m) | Friction Slope (m/km) |
|---|
At any pipe junction, according to the conservation of mass, the total flow entering the node must equal the total flow leaving it. In this N-reservoir system, the solver iteratively finds the unique junction head ($H_j$) that satisfies:
Flow in each pipe depends on the head difference between the reservoir ($Z_i$) and the junction ($H_j$). Using the empirical Hazen-Williams formula (for metric units):
Where $H_f = |Z_i - H_j|$. Solving for Q gives the flow rate for a given junction head.
Because the flow equation is non-linear ($Q \propto H_f^{0.54}$), a direct algebraic solution is impossible. This tool employs the Bisection Method. It sets the initial bounds for $H_j$ between the highest and lowest reservoir elevations. It guesses the midpoint, calculates the net flow $\Sigma Q$, and iteratively halves the search interval until $\Sigma Q \approx 0$ (within a $10^{-6}$ error margin). This method guarantees absolute convergence for any number of pipes.
- Hazen, A., & Williams, G. S. (1920). Hydraulic Tables: The Elements of Gagings and the Friction of Water Flowing in Pipes (3rd ed.). John Wiley & Sons.
- Mott, R. L., & Untener, J. A. (2014). Applied Fluid Mechanics (7th ed.). Pearson.
- Bhave, P. R. (1991). Analysis of Water Distribution Networks. Technomic Publishing Co.
No comments:
Post a Comment