Numerical modeling of the flow and transport of radionuclides in heterogeneous porous media

This paper concerns numerical methods for the modeling of ow and transport of contaminant in porous media. The numerical methods feature the mixed nite element method over triangles as a solver to the Darcy ow equation, and a conser-vative nite volume scheme for the concentration equation. The convective term is approximated with the aid of a Godunov scheme considered over the nite volume mesh dual to a triangular grid, whereas the diusion-dispersion term is discretized by piecewise linear conforming triangular nite elements. It is shown that the scheme satis es a discrete maximum principle. Numerical examples demonstrate the eec-tiveness of the methodology for a coupled system which includes an elliptic equation and a diusion-convection-reaction equation arising in modeling of ow and trans-port in heterogeneous porous media. The proposed scheme is robust, conservative, ecient, and stable which is conrmed by numerical simulations