EPANET is one of the most commonly used open-source programs in hydraulic modelling water distribution networks (WDNs), based on steady-state and extended period simulation approaches. These approaches effectively estimate flow capacity and average pressures in networks; however, EPANET is not yet fully effective in modelling incompressible unsteady flows in WDNs. In this study, the hydraulic solver capacity of EPANET 3 is extended with the Rigid Water Column Global Gradient Algorithm (RWC-GGA) to model incompressible unsteady flow hydraulics in WDNs. Moreover, we incorporated dynamically more accurate valve expressions than the existing ones in the default EPANET code and introduced a new global convergence algorithm, Convergence Tracking Control Method (CTCM), in the solver code. The RWC-GGA, CTCM, and valve expressions are tested and validated in three different WDNs varying from simple to sophisticated set-ups. The results show that incompressible unsteady flows can be modelled with RWC-CGA and dynamic valve representations. Finally, the convergence problem due to the valve motion and the pressure-dependent algorithm (PDA) is solved by the implemented global convergence algorithm, i.e. CTCM.