Numerical modeling of three-dimensional contaminant migration from dissolution of multicomponent NAPL pools in saturated porous media

Abstract

A three-dimensional model for contaminant transport resulting from the dissolution of multicomponent nonaqueous phase liquid (NAPL) pools in three-dimensional saturated subsurface formations is developed. The solution is obtained numerically by a finite-difference scheme, and it is suitable for homogeneous porous media with unidirectional interstitial velocity. Each dissolved component may undergo first-order decay and may sorb under local equilibrium conditions. It is also assumed that the dissolution process is mass transfer limited. The nonaqueous phase activity coefficients of the NAPL pool components are evaluated at each time step. The model behavior is illustrated through a synthetic example with a NAPL pool consisting of a mixture of TCA (1,1,2-trichloroethane) and TCE (trichloroethylene). The numerical solution presented in this work is in good agreement with a recently developed analytical solution for the special case of a single component NAPL pool. The results indicate the importance of accounting for the necessary changes in the organic phase activity which significantly affects the equilibrium aqueous solubility.