A pore-throat model based on grain-size distribution to quantify gravity-dominated DNAPL instabilities in a water-saturated homogeneous porous medium

A pore-scale network model based on spherical pore bodies and cylindrical pore throats was developed to describe the displacement of water by DNAPL. The pore body size and the pore throat size were given by statistical distributions with user-specified values for the minimum, mean and maximum sizes. The numerical model was applied to a laboratory experiment conducted on a sand-filled glass column. The parameters relative to pore body size and pore throat size that were used in the construction of the equivalent network were derived from the discrete grain-size distribution of the real porous medium. The calculated arrival times of the DNAPL front were compared with those measured using optic fibre sensors placed at different points on the control section of the experimental device. Furthermore, the model simulated DNAPL pressure measured at the entrance section of the system. In general, the numerical results obtained with the model were in good agreement with the actual measurements.