Coupling of a two phase gas liquid compositional 3D Darcy flow with a 1D compositional free gas flow∗

A model coupling a three dimensional gas liquid compositional Darcy flow and a one dimensional compositional free gas flow is presented. The coupling conditions at the interface between the gallery and the porous medium account for the molar normal fluxes continuity for each component, the gas liquid thermodynamical equilibrium, the gas pressure continuity and the gas and liquid molar fractions continuity. This model is applied to the simulation of the mass exchanges at the interface between the repository and the ventilation excavated gallery in a nuclear waste geological repository. The spatial discretization is essentially nodal and based on the vertex approximate gradient (VAG) scheme. Compared with classical nodal approaches such as the Control Volume Finite Element method, the VAG scheme has the advantage to avoid the mixture of different material properties and models in the control volumes located at the interfaces. The discrete model is validated using a quasi analytical solution for the stationary state, and the convergence of the VAG discretization is analysed for a simplified model coupling the Richards approximation in the porous medium and the gas pressure equation in the gallery.