Computationally efficient modeling method for large water network analysis