This paper describes the modifications applied to EPANET, a public-domain water distribution system modelling software package, that does not correctly compute the hydraulics of a water distribution network (WDN) with variable tank heads in (slow) unsteady flow conditions. Firstly the methodology adopted to extend the GGA algorithm implemented in the original EPANET source code to the EPS-GGA is described. Then the convergence and stability conditions of the theta method, used for the discretisation in time of the set of differential equations describing the hydraulic behaviour of a WDN, are discussed. The reasons for EPS-GGA numerical stability are demonstrated and a fully implicit discretisation of differential equations (i.e. theta = 1) is suggested as the optimal choice as implicitly proposed in Giustolisi et al. but without theoretical justification. Both the modified and original versions of EPANET are applied to a particularly severe test case of a WDN. Moreover, the procedures for the correct numerical representation of the tanks' maximum and minimum level boundary conditions are developed and compared with previously proposed procedures. The modified version of EPANET source code does not show the significant instabilities which are evident in the original version, nor the lack of consistency due to the improper maximum and minimum level boundary condition schematisations formerly proposed in the scientific literature.

The extension of EPANET source code to simulate unsteady flow in water distribution networks with variable head tanks

Avesani, Diego;Righetti, Maurizio;Bertola, Paolo
2012

Abstract

This paper describes the modifications applied to EPANET, a public-domain water distribution system modelling software package, that does not correctly compute the hydraulics of a water distribution network (WDN) with variable tank heads in (slow) unsteady flow conditions. Firstly the methodology adopted to extend the GGA algorithm implemented in the original EPANET source code to the EPS-GGA is described. Then the convergence and stability conditions of the theta method, used for the discretisation in time of the set of differential equations describing the hydraulic behaviour of a WDN, are discussed. The reasons for EPS-GGA numerical stability are demonstrated and a fully implicit discretisation of differential equations (i.e. theta = 1) is suggested as the optimal choice as implicitly proposed in Giustolisi et al. but without theoretical justification. Both the modified and original versions of EPANET are applied to a particularly severe test case of a WDN. Moreover, the procedures for the correct numerical representation of the tanks' maximum and minimum level boundary conditions are developed and compared with previously proposed procedures. The modified version of EPANET source code does not show the significant instabilities which are evident in the original version, nor the lack of consistency due to the improper maximum and minimum level boundary condition schematisations formerly proposed in the scientific literature.
4
Avesani, Diego; Righetti, Maurizio; D., Righetti; Bertola, Paolo
File in questo prodotto:
File Dimensione Formato  
Avesani Righetti et al_2012 Journal of Hydroinformatic_ HYDRO-D-11-00013-4.pdf

Solo gestori archivio

Tipologia: Versione editoriale (Publisher’s layout)
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 719.87 kB
Formato Adobe PDF
719.87 kB Adobe PDF   Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11572/93849
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 11
  • ???jsp.display-item.citation.isi??? 12
social impact