Testing linear solvers for global gradient algorithm

Steady-state Water Distribution Network models compute pipe flows and nodal heads for assumed nodal demands, pipe hydraulic resistances, etc. The nonlinear mathematical problem is based on energy and mass conservation laws which is solved by using global linearization techniques, such as global gradient algorithm (GGA). The matrix of coefficients of the linear system inside GGA belongs to the class of sparse, symmetric and positive definite. Therefore a fast solver for the linear system is important in order to achieve the computational efficiency, especially when multiple runs are required. This work aims at testing three main strategies for the solution of linear systems inside GGA. The tests are performed on eight real networks by sampling nodal demands, considering the pressure-driven and demand-driven modelling to evaluate the robustness of solvers. The results show that there exists a robust specialized direct method which is superior to all the other alternatives. Furthermore, it is found that the number of times the linear system is solved inside the GGA does not depend on the specific solver, if a small regularization to the linear problem is applied, and that pressure-driven modelling requires a greater number which depends on the size and topology of the network and not only on the level of pressure deficiency.

[1]  Daniele B. Laucelli,et al.  Deterministic versus Stochastic Design of Water Distribution Networks , 2009 .

[2]  Zoran Kapelan,et al.  Using high performance techniques to accelerate demand-driven hydraulic solvers , 2013 .

[3]  Pramod R. Bhave,et al.  Comparison of Methods for Predicting Deficient-Network Performance , 1996 .

[4]  Timothy A. Davis,et al.  Row Modifications of a Sparse Cholesky Factorization , 2005, SIAM J. Matrix Anal. Appl..

[5]  Aaron C. Zecchin,et al.  Steady-State Behavior of Large Water Distribution Systems: Algebraic Multigrid Method for the Fast Solution of the Linear Step , 2012 .

[6]  Olivier Piller,et al.  Least Action Principles Appropriate to Pressure Driven Models of Pipe Networks , 2003 .

[7]  J. E. van Zyl,et al.  The potential of graphical processing units to solve hydraulic network equations , 2012 .

[8]  Avi Ostfeld,et al.  The Battle of the Water Sensor Networks (BWSN): A Design Challenge for Engineers and Algorithms , 2008 .

[9]  Thomas M. Walski,et al.  Extended Global-Gradient Algorithm for Pressure-Dependent Water Distribution Analysis , 2009 .

[10]  Orazio Giustolisi,et al.  Pressure-Driven Demand and Leakage Simulation for Water Distribution Networks , 2008 .

[11]  Istvan Lippai.  Colorado Springs Utilities Case Study: Water System Calibration / Optimization , 2005 .

[12]  Luigi Berardi,et al.  Testing linear solvers for WDN models , 2011 .

[13]  E Todini,et al.  A more realistic approach to the “extended period simulation” of water distribution networks , 2003 .

[14]  Richard Burd,et al.  BATTLE OF THE WATER CALIBRATION NETWORKS (BWCN) , 2011 .

[15]  Y. Notay.  An aggregation-based algebraic multigrid method , 2010 .

[16]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[17]  J. E. Van Zyl,et al.  EXTENSION OF EPANET FOR PRESSURE DRIVEN DEMAND MODELING IN WATER DISTRIBUTION SYSTEM , 2005 .

[18]  Orazio Giustolisi,et al.  Demand Components in Water Distribution Network Analysis , 2012 .

[19]  E. Todini,et al.  Unified Framework for Deriving Simultaneous Equation Algorithms for Water Distribution Networks , 2013 .

[20]  J. Pasciak,et al.  Computer solution of large sparse positive definite systems , 1982 .

[21]  Luigi Berardi,et al.  A computationally efficient modeling method for large size water network analysis , 2012 .

[22]  Long Chen.  INTRODUCTION TO MULTIGRID METHODS , 2005 .

[23]  Orazio Giustolisi,et al.  Considering Actual Pipe Connections in WDN Analysis , 2010 .

[24]  Timothy A. Davis,et al.  Modifying a Sparse Cholesky Factorization , 1999, SIAM J. Matrix Anal. Appl..

[25]  Timothy A. Davis,et al.  Multiple-Rank Modifications of a Sparse Cholesky Factorization , 2000, SIAM J. Matrix Anal. Appl..

[26]  Timothy A. Davis,et al.  Algorithm 8 xx : a concise sparse Cholesky factorization package , 2004 .

[27]  Yvan Notay,et al.  Aggregation-Based Algebraic Multilevel Preconditioning , 2005, SIAM J. Matrix Anal. Appl..

[28]  Jakobus E. van Zyl,et al.  Operational Optimization of Water Distribution Systems using a Hybrid Genetic Algorithm , 2004 .