Pore-scale simulation of entrapped non-aqueous phase liquid dissolution

We investigated the dissolution of non-aqueous phase liquids (NAPLs) in a three-dimensional random sphere-pack medium using a pore-scale modeling approach to advance fundamental understanding and connect rigorously to microscale processes. Residual NAPL distributions were generated using a morphological approach and the entrapped non-wetting phase was quantitatively characterized by calculating volume, orientation, interfacial area, and shape of isolated NAPL regions. With a detailed aqueous-phase flow field obtained by a multiple-relaxation time lattice Boltzmann approach, we solved the advective–diffusive equation in the pore space using a high-resolution, adaptive-stencil finite-volume scheme and an operator-splitting algorithm. We show good agreement between the mass transfer rates predicted in the computational approach and previously published experimental observations. The pore-scale simulations presented in this work provide the first three-dimensional comparison to the considerable experimental work that has been performed to derive constitutive relations to quantify mass transfer from a residual NAPL to a flowing aqueous phase.