It is common that non-convergence problems occur for simulations under normal or abnormal operations. Among the main causes, there are the use of sublinear functions for calculation of the hydraulic steady-state, some derivative issues and an ill-conditioned iteration matrix in the solution method. The objective of this paper is to specify non-convergence cases for existing hydraulic software solutions. The damped Newton solution together with regularization techniques are proposed for ensuring global convergence, whatever the initial solution. Failures of convergence are shown to exist even on small case studies.