In the subsurface, small but significant quantities of chemicals will be dissolved by groundwater in contact with Non-aqueous Phase Liquids (NAPL) and can result in groundwater contamination. The rate of NAPL dissolution will influence the magnitude of dissolved chemical concentrations in groundwater and the length of time that NAPL sources may cause groundwater contamination. A suitable understanding of NAPL fate in the subsoil is required in order to evaluate the sources of existing dissolved chemical plumes, to allow predictions of future plume migration and to evaluate the effectiveness of dissolved plume purging in situations where NAPL remains in the subsurface. In this study the mathematical formulation necessary to represent the behavior of NAPL in the saturated zone is presented. The results of the simulation of a trichloroethene spill demonstrates that immiscible liquids can act as a continuous contamination source in an aquifer for many decades, even in small quantities.