The modelling of chlorine residual in water supply systems is of great importance in managing disinfectant concentrations throughout the network. First order decay kinetics are currently often used to describe both bulk and wall chlorine consumption. However, more complex approaches have been proposed, namely a parallel, two-reactant second order decay model (2R model) that has been reported as yielding better accuracy for simulating chlorine bulk decay in laboratory tests.