A functional relation for field-scale nonaqueous phase liquid dissolution developed using a pore network model

A pore network model with cubic chambers and rectangular tubes was used to estimate the nonaqueous phase liquid (NAPL) dissolution rate coefficient, K (diss)a(i), and NAPL/water total specific interfacial area, a(i), K(diss)a( i) was computed as a function of modified Peclet number (Pe') for various N APL saturations (S-N) and a(i) during drainage and imbibition and during di ssolution without displacement. The largest contributor to a(i) was the int erfacial area in the water-filled corners of chambers and tubes containing NAPL. When K(diss)a(i) was divided by a(i), the resulting curves of dissolu tion coefficient, K-diss versus Pe' suggested that an approximate value of K-diss could he obtained as a weak function of hysteresis or S-N. Spatially and temporally variable maps of K(diss)a(i) calculated using the network m odel were used in field-scale simulations of NAPL dissolution. These simula tions were compared to simulations using a constant value of K(diss)a(i) an d the empirical correlation of Powers et al. [Water Resour. Res. 30(2) (199 4b) 321]. Overall, a methodology was developed fut. incorporating pore-scal e processes into field-scale prediction of NAPL dissolution. (C) 2001 Publi shed by Elsevier Science B.V.