Non-isothermal compositional liquid gas Darcy flow: formulation, soil-atmosphere boundary condition and application to high-energy geothermal simulations

This article deals with the modelling and formulation of compositional gas liquid Darcy flow. Our model includes an advanced boundary condition at the interface between the porous medium and the atmosphere accounting for convective mass and energy transfer, liquid evaporation and liquid outflow. The formulation is based on a fixed set of unknowns whatever the set of present phases. The thermodynamic equilibrium is expressed as complementarity constraints. The model and its formulation are applied to the simulation of the Bouillante high-energy geothermal field in Guadeloupe characterised by a high temperature close to the surface.