Impact of NAPL architecture on interphase mass transfer: A pore network study

Abstract

Interphase mass transfer in porous media is commonly modeled using Sherwood number expressions that are developed in terms of fluid and porous medium properties averaged over some representative elementary volume (REV). In this work the influence of sub-grid scale properties on interphase mass transfer was investigated using a two-dimensional pore network model. The focus was on assessing the impact of (i) NAPL saturation, (ii) interfacial area (iii) NAPL spatial distribution at the pore scale, (iv) grain size heterogeneity, (v) REV or domain size and (vi) pore scale heterogeneity of the porous media on interphase mass transfer. Variability of both the mass transfer coefficient that explicitly accounts for the interfacial area and the mass transfer coefficient that lumps the interfacial area was examined. It was shown that pore scale NAPL distribution and its orientation relative to the flow direction have significant impact on flow bypassing and the interphase mass transfer coefficient. This results in a complex non-linear relationship between interfacial area and the REV-based interphase mass transfer rate. Hence, explicitly accounting for the interfacial area does not eliminate the uncertainty of the mass transfer coefficient. It was also shown that, even for explicitly defined flow patterns, changing the domain size over which the mass transfer process is defined influences the extent of NAPL bypassing and dilution and, consequently, the interphase mass transfer. It was also demonstrated that the spatial variability of pore scale parameters such as pore throat diameters may result in different rates of interphase mass transfer even for the same pore size distribution index.