Steady-state behavior of large water distribution systems: Algebraic multigrid method for the fast solution of the linear step

dc.contributor.authorZecchin, A.
dc.contributor.authorThum, P.
dc.contributor.authorSimpson, A.
dc.contributor.authorTischendorf, C.
dc.date.issued2012
dc.description.abstractThe Newton-based global gradient algorithm (GGA) (also known as the Todini and Pilati method) is a widely used method for computing the steady-state solution of the hydraulic variables within a water distribution system (WDS). The Newton-based computation involves solving a linear system of equations arising from the Jacobian of the WDS equations. This step is the most computationally expensive process within the GGA, particularly for large networks involving up to O(10<sup>5</sup>) variables. An increasingly popular solver for large linear systems of the M-matrix class is the algebraic multigrid (AMG) method, a hierarchical-based method that uses a sequence of smaller dimensional systems to approximate the original system. This paper studies the application of AMG to the steady-state solution of WDSs through its incorporation as the linear solver within the GGA. The form of the Jacobian within the GGA is proved to be an M-matrix (under specific criteria on the pipe resistance functions), and thus able to be solved using AMG. A new interpretation of the Jacobian from the GGA is derived, enabling physically based interpretations of the AMG's automatically created hierarchy. Finally, extensive numerical studies are undertaken where it is seen that AMG outperforms the sparse Cholesky method with node reordering (the solver used in EPANET2), incomplete LU factorization (ILU), and PARDISO, which are standard iterative and direct sparse linear solvers. © 2012 American Society of Civil Engineers.
dc.description.statementofresponsibilityA. C. Zecchin; P. Thum; A. R. Simpson; and C. Tischendorf
dc.identifier.citationJournal of Water Resources Planning and Management, 2012; 138(6):639-650
dc.identifier.doi10.1061/(ASCE)WR.1943-5452.0000226
dc.identifier.issn0733-9496
dc.identifier.issn1943-5452
dc.identifier.orcidZecchin, A. [0000-0001-8908-7023]
dc.identifier.orcidSimpson, A. [0000-0003-1633-0111]
dc.identifier.urihttp://hdl.handle.net/2440/76746
dc.language.isoen
dc.publisherASCE-Amer Soc Civil Engineers
dc.rights© 2012 American Society of Civil Engineers
dc.source.urihttps://doi.org/10.1061/(asce)wr.1943-5452.0000226
dc.subjectWater distribution systems
dc.subjectGlobal gradient algorithm
dc.subjectSteady-state
dc.subjectAlgebraic multigrid methods
dc.titleSteady-state behavior of large water distribution systems: Algebraic multigrid method for the fast solution of the linear step
dc.typeJournal article
pubs.publication-statusPublished

Files

Original bundle
Now showing 1 - 1 of 1
No Thumbnail Available
Name:
hdl_76746.pdf
Size:
457.71 KB
Format:
Adobe Portable Document Format
Description:
Accepted version