Please use this identifier to cite or link to this item:
https://hdl.handle.net/2440/76746
Citations | ||
Scopus | Web of Science® | Altmetric |
---|---|---|
?
|
?
|
Type: | Journal article |
Title: | Steady-state behavior of large water distribution systems: Algebraic multigrid method for the fast solution of the linear step |
Author: | Zecchin, A. Thum, P. Simpson, A. Tischendorf, C. |
Citation: | Journal of Water Resources Planning and Management, 2012; 138(6):639-650 |
Publisher: | ASCE-Amer Soc Civil Engineers |
Issue Date: | 2012 |
ISSN: | 0733-9496 1943-5452 |
Statement of Responsibility: | A. C. Zecchin; P. Thum; A. R. Simpson; and C. Tischendorf |
Abstract: | The 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(105) 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. |
Keywords: | Water distribution systems Global gradient algorithm Steady-state Algebraic multigrid methods |
Rights: | © 2012 American Society of Civil Engineers |
DOI: | 10.1061/(ASCE)WR.1943-5452.0000226 |
Published version: | http://dx.doi.org/10.1061/(asce)wr.1943-5452.0000226 |
Appears in Collections: | Aurora harvest 4 Civil and Environmental Engineering publications Environment Institute publications |
Files in This Item:
File | Description | Size | Format | |
---|---|---|---|---|
hdl_76746.pdf | Accepted version | 457.71 kB | Adobe PDF | View/Open |
Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.