An implicit time integration method for the simulation of steady and u
nsteady flow in pipes and channels is presented. It is based on the th
eory of Total Variation Diminishing (TVD) methods. A conservative line
arization leads to a block tridiagonal system of equations which can b
e cheaply solved by means of a non-iterative matrix decomposition meth
od. It keeps the advantages of classical implicit schemes, and properl
y deals with all kinds of flow. It is specially suited for steady flow
s including phenomena such as hydraulic jumps and also gives satisfact
ory results for time dependent problems. Several test cases are shown
to illustrate the performance of this implicit technique in single cha
nnels and networks.