The presentation will begin with a derivation of a model of nuclear contaminant transport in naturally-fractured porous media. Then, a Locally-Conservative Eulerian Lagrangian Method (LCELM) will be implemented to numerically solve the system of equations that arises. Several numerical results will be shown. In particular, the effects of the fractures on the flow will be illustrated.