Previous work has shown that Lagrangian methods are more efficient for modeling the transport of chemicals in a water distribution system. Two such methods, the Lagrangian Time-Driven Method (TDM) and Event-Driven Method (EDM) are compared for varying concentration tolerance and computational water quality time step. A new hybrid method (EDMNET) is developed which improves the accuracy of the Lagrangian methods. All the above methods are incorporated in an existing hydraulic simulation model. The integrated model is run for different network problems under varying conditions. The TDM-generated solutions are affected by both concentration tolerance and water quality time step, whereas EDM solutions are dependent on concentration tolerance. The EDMNET solutions are less sensitive to variations in these parameters. The threshold solutions are determined for all the methods and compared. The hybrid method simulates the nodal concentrations accurately with least maximum segmentation of network and reasonable computational effort as compared to the other Lagrangian methods.