Dense, non-aqueous phase liquids (DNAPLs) are common organic contaminants in subsurface environment. Once spilled or leaked underground, they slowly dissolved into groundwater
and generated a plume of contaminants. In order to manage the contaminated site and
predict the behavior of dissolved DNAPL in heterogeneous subsurface requires a comprehensive
numerical model. In this work, the Crank-Nicolson finite-element Galerkin (CN-FEG) numerical
scheme for solving a set coupled system of partial differential equations that describes fate and
transport of dissolved organic compounds in two-dimensional domain was developed and implemented.
Assumptions are made so that the code can be compared and verified with available