OPTIMIZED MUSKINGUM-CUNGE SOLUTION METHOD FOR SOLUTE TRANSPORT WITH EQUILIBRIUM FREUNDLICH REACTIONS

Because of its simplicity and efficiency, the Muskingum-Cunge (M-C) me thod is a widely used modelling tool applied to flood routing problems . However, in the M-C method, there is an infinite number of combinati ons of spatial and temporal step sizes for which the solution is feasi ble and not all of these lead to accurate results. This has been a maj or criticism of the method. The same approach has been adopted here to solve the problem of nonlinear reactive solute transport. An optimise d M-C scheme is proposed and tested which is very simple and efficient as compared with existing standard numerical solutions such as the Cr ank-Nicolson scheme. The optimisation principles are: (1) that the col umn holdup (area above the breakthrough curve normalised by the column length) is maintained; and (2) that the truncation error produced by the numerical scheme is minimised. The conditions derived are shown to be satisfactory for the mildly nonlinear reactions of the Freundlich- type typical of reactive chemical transport in soil.