DNAPL Infiltration in a Two-Dimensional Porous Medium—Influence of the Shape of the Solid Particles

Mustapha Hellou, Trong Dong Nguyen, Pascal Dupont
DOI: 10.4236/eng.2011.312148                

Abstract

The infiltration with atmospheric pressure of Dense Non Aqueous Phase Liquid (DNAPL) in a model of porous medium saturated by another liquid is studied when this DNAPL liquid has a contact angle characterizing wetting liquid. The model of the porous medium considered consists of an assembly of solid particles for various forms. The influence of the shape of the particles is studied. The results found show the retention capacity of such porous media in function of the shape of the solid particles.

Keywords

DNAPL, Multiphase Flow, Wetting Liquid, Solid Particle, VOF, Retention of Liquid Pollutant

1. Introduction

Multiphase flows in porous media involving immiscible fluids always interest scientists in terms of applied and fundamental point of view. For example, the transport of immiscible fluid with water known as NAPL (Non Aqueous Phase Liquid) is encountered in geomaterials, processes of remediation of the grounds, biofilters employed in waste water treatment, etc. These flows are complex and often unknown because of the large number of physical parameters involved to describe their properties. At the present time, influence of parameters as density, viscosity or nterfacial tension on the infiltration of NAPL into porous media is much studied by several authors [1-6]. Most of these studies are realized at a macroscopic scale. This macroscopic approach does not make it possible to understand the retention of small drops of DNAPL (Dense Non Aqueous Phase Liquid) in the pores of the medium. However, few studies carried out at the scale of the interstices showed the dispersion of the DNAPL as well in experiments and in numerical simulation [7-9].

In this work, we highlight the influence of a parameter which is poorly studied—geometry of the solid particles of porous media—on the infiltration and the dispersion of a drop of DNAPL. Numerical simulation is carried out for various shapes of the solid particles.

2. Problem Position

The fluids used in this study are water (carrier fluid) and the DNAPL (the drop). The density and the viscosity of water are of 998 kg·m^–3 and 1 × 10^–3 Pa·s. Those of the DNAPL are respectively equal to 1623 kg·m^–3 and 0.89 Pa·s. As an example, Perchlorethylen (PCE) possesses similar properties (INERIS, 2002). A contact angle of 65˚ is imposed in order to consider a wetting DNAPL because most DNAPL’s being in nature are wetting.

Figure 1 presents one of the three structures of porous

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	K. D. Pennell, Pope G. A. and L. M. Abriola, “Influence of Viscous and Buoyancy Forces on the Mobilization of Residual Tetrachloroethylene during Surfactant Flushing”, Environ. Sci. Technol., Vol. 30, NO. 4, 1996, pp. 1328-1335.
[2]	H. E. Dawson and P. V. Roberts, “Influence of Viscous, Gravitational, and Capillary Forces on DNAPL Saturation”, Ground Water, Vol. 35, N 2, 1997, pp. 261-269.
[3]	C. Hofstee, C. G. Ziegler, O. Tr?tschler and J. Braun, “Removal of DNAPL contamination from the saturated zone by the combined effect of vertical upward flushing and density reduction”, Journal of Contaminant Hydrology, Vol. 67, NO. 1, 2003, pp. 61-78.
[4]	S. W. Jeong and M-Y. Corapcioglu, “ Force analysis and visualization of NAPL removal during surfactant-related floods in a porous medium”, Journal of Hazardous Materials, A126, 2005, pp. 8-13.
[5]	A.M., Tartakovsky, A. L. Ward and P. Meakin,” Hetero- geneity Effects on Capillary Pressure-Saturation Relations Inferred from Pore-Scale Modeling”, Physics of Fluids, 19, 103301, DOI: 10.1063/1.2772529, 2007
[6]	P. Meakin, and A.M. Tartakovsky, “Modeling and simulation of pore scale multiphase fluid flow and reactive transport in fractured and porous media”, Reviews of Geophysics, 47, RG3002, doi:10.1029/2008RG00263, 2009.
[7]	R. Krishna and J. M.Van Baten, “Rise characteristics of gaz bubbles in a 2D rectangular column : VOF simulations vs experiments”, Int. Comm. Heat Mass Transfer, Vol. 66, NO.7, 1999, pp. 965-974.
[8]	G. J. Storr and M. Behnia, “ Comparisons between experiment and numerical simulation using a free surface technique of free-falling liquid jets”, Experimental Thermal and Fluid Science,Vol. 22, 2000, pp. 79-91.
[9]	D. J. E. Harvie, M. R. Davidson, J. J. Cooper-White and M. Rudman, “A parametric study of droplet deformation through a microfluidic contraction: Low viscosity Newtonian droplets”, Chemical Engineering Science, Vol. 61, 2006, pp. 5149-5158.
[10]	C. W. Hirt and B. D. Nichols, “Volume of Fluid (VOF) method for the dynamics of free boundaries”, J. Comp. phys., Vol. 39, 1981, pp. 201-225.
[11]	Nguyen T. D., “Infiltration de particules liquids ou solides dans un milieu poreux”, PhD Thesis, 2007, INSA Rennes France.
[12]	J.S. Hadamard, Mouvement permanent lent d'une sphère liquide et visqueuse dans un liquide visqueux. C. R. Acad. Sci. 152, Paris, 1735-1752, 1911.
[13]	W. Rybczynski, “On the translatory motion of a fluid sphere in a viscous medium”, Bull. Acad. Sci., Cracow, Series A, p. 40, 1911