An attempt to improve numerical stability of EPANET 2 hydraulic solver