Please use this identifier to cite or link to this item:
Scopus Web of Science® Altmetric
Type: Journal article
Title: Dealing with zero flows in solving the nonlinear equations for water distribution systems
Author: Elhay, S.
Simpson, A.
Citation: Journal of Hydraulic Engineering, 2011; 137(10):1216-1224
Publisher: ASCE-Amer Soc Civil Engineers
Issue Date: 2011
ISSN: 0733-9429
Statement of
Sylvan Elhay and Angus R. Simpson
Abstract: Three issues concerning the iterative solution of the nonlinear equations governing the flows and heads in a water distribution system network are considered. Zero flows cause a computation failure (division by zero) when the Global Gradient Algorithm of Todini and Pilati is used to solve for the steady state of a system in which the head loss is modeled by the Hazen-Williams formula. The proposed regularization technique overcomes this failure as a solution to this first issue. The second issue relates to zero flows in the Darcy-Weisbach formulation. This work explains for the first time why zero flows do not lead to a division by zero where the head loss is modeled by the Darcy-Weisbach formula. In this paper, the authors show how to handle the computation appropriately in the case of laminar flow (the only instance in which zero flows may occur). However, as is shown, a significant loss of accuracy can result if the Jacobian matrix, necessary for the solution process, becomes poorly conditioned, and so it is recommended that the regularization technique be used for the Darcy-Weisbach case also. Only a modest extra computational cost is incurred when the technique is applied. The third issue relates to a new convergence stopping criterion for the iterative process based on the infinity-norm of the vector of nodal head differences between one iteration and the next. This test is recommended because it has a more natural physical interpretation than the relative discharge stopping criterion that is currently used in standard software packages such as EPANET. In addition, it is recommended to check the infinity norms of the residuals once iteration has been stopped. The residuals test ensures that inaccurate solutions are not accepted. © 2011 American Society of Civil Engineers.
Keywords: Water distribution systems
Solving equations
Newton method
Todini and Pilati
Zero flows
Rights: © 2011 American Society of Civil Engineers
DOI: 10.1061/(ASCE)HY.1943-7900.0000411
Appears in Collections:Aurora harvest 7
Computer Science publications
Environment Institute publications

Files in This Item:
File Description SizeFormat 
hdl_69036.pdfAccepted version637.51 kBAdobe PDFView/Open

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