Among the most intractable environmental remediation problems are those involving the release of dense non-aqueous phase liquids (DNAPLs), such as chlorinated solvents, to the subsurface. Research efforts have focused on the use of numerical models to investigate reductions in contaminant concentrations due to partial mass removal and improvements in the performance of complementary source zone remediation technologies. Previous numerical investigations, however, have been limited to two-dimensional systems. Furthermore, a lack of models capable of simulating the most promising complementary technology, metabolic reductive dechlorination, has limited its application. This work developed and applied compositional multiphase numerical simulators to examine the influence of dimensionality (two-dimensions versus three-dimensions) on DNAPL source zone simulations and to investigate the benefits of stimulating metabolic reductive dechlorination at a chlorinated ethene-DNAPL contaminated site. Results from the dimensionality investigation showed that the simulation of DNAPL migration, entrapment, and dissolution in two dimensions provided reasonable approximations to the behavior simulated in three dimensions. Commonly employed saturation distribution and mass recovery metrics were approximately equivalent. Flux- averaged concentrations simulated in two dimensions, however, tended to be three to four times higher than those simulated in three dimensions. This difference was attributed to dilution at the down gradient boundary. An alternative metric, mass flux reduction, however, yielded better agreement.