Please use this identifier to cite or link to this item:
Scopus Web of Science® Altmetric
Type: Conference paper
Title: Formulating the water distribution system equations in terms of head and velocity
Author: Simpson, A.
Elhay, S.
Citation: WDSA 2008: Proceedings of the 10th Annual Water Distribution Systems Analysis Conference, August 17-20, 2008, Kruger National Park, South Africa / K. Van Zyl (ed.): pp. 790-802
Publisher: American Society of Civil Engineers
Issue Date: 2008
ISBN: 9780784410240
Conference Name: Annual Symposium on Water Distribution Systems Analysis (10th : 2008 : Kruger National Park, South Africa)
Statement of
Angus R. Simpson and Sylvan Elhay
Abstract: The set of equations for solving for pressures and flows in water distribution systems are non‐linear due to the head loss‐velocity relationship for each of the pipes. The solution of these non‐linear equations for the heads and flows is usually based on the Todini and Pilati method. The method is an elegant way of formulating the equations. A Newton solution method is used to solve the equations whereby the special structure of the Jacobian is exploited to minimize the computations and this leads to an extremely fast algorithm. Each iteration firstly solves for the heads and then solves for the flows. In the EPANET implementation of the Todini and Pilati algorithm an initial guess of the flows is based on an assumed velocity of 1.0 fps (0.305 m/s) in each pipe in the network. Each flow is then determined from the continuity equation by multiplying the assumed velocity by the area. Usually velocities in pipes are in the range of 0.5 to 1.5 m/s (and perhaps sometimes higher up to 3 or 4 m/s). Thus the velocities to be solved for are all of the same order of magnitude. In contrast, the range of discharges may be quite large in a system — ranging from below 10 L/s up to above 700 L/s — thus possibly three orders of magnitude of difference. As an alternative to the usual formulation of the Todini and Pilati method in terms of flows and heads, this paper recasts the Todini and Pilati formulation in terms of heads and velocities to attempt to improve the convergence properties. Results are compared for the two formulations for a range of networks from 553 to 10,354 pipes. Convergence criteria for stopping the iterative solution process are discussed. The impact of the initial guess of the velocities in each of the pipes in the network on the convergence behavior is also investigated. Statistics on mean flows and velocities in the network and the minimum and maximum velocities for each of the example networks are given and finally operation counts are also provided for these networks.
Rights: © 2008 American Society of Civil Engineers
DOI: 10.1061/41024(340)68
Published version:
Appears in Collections:Aurora harvest 7
Civil and Environmental Engineering publications
Environment Institute publications

Files in This Item:
File Description SizeFormat 
hdl_66296.pdfAccepted version310.98 kBAdobe PDFView/Open

Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.