A numerical model for the solute transport and dispersion in saturated porous media has been developed. The partial differential equations for water flow and solute transport are discretized using the finite difference technique and the resulting system of algebraic equations is solved using a dynamic programming method. The advantage of this method is that the problem is converted from solving an algebraic system of this method is that the problem is converted from solving an algebraic system of order NC(NL-1)XNC(NL-1) into that of solving a difference equation of order NCXNC over NL-1 steps and involving NL-1 matrix inversions of order NCxNC. The accuracy and precision of the solutions are shown by calculation of mass balance and comparing the results with MOC model developed by USGS.
展开▼