Non-isothermal two-phase flow in deformable porous media: systematic open-source implementation and verification procedure

We propose a formulation for non-isothermal two-component two-phase flow through deformable porous media. The approach covers phase transitions among both phases, i.e. liquid phase components evaporate into the gas phase while gas phase components dissolve or condensate into the liquid phase. These phase transitions always take place in thermodynamic equilibrium. The set of model equations is thereby largely independent of the specific constitutive relations. Starting from general equilibrium equations, we show the evolution of the system of weak formulations of all governing equations, which are then discretised with Taylor-Hood elements in a standard finite element approach. The model equations and the construction of the constitutive equilibria are implemented in the open-source simulator OpenGeoSys, which can be freely used and modified. To verify the implementation, we have selected a number of complementary test cases covering a wide range of process couplings. The numerical model is compared with analytical and semi-analytical solutions of these problems as well as with experimental results. It is shown in the paper that by including thermodynamic effects, solid mechanics, and phase transition processes, the proposed numerical model covers many characteristic features of unsaturated geomaterials and can be employed for the description of a broad range of problems encountered in geotechnical engineering.

Automated summary

This paper proposes a formulation for non-isothermal two-component two-phase flow through deformable porous media, and demonstrates the implementation of the model equations and construction of constitutive equilibria. Through verification experiments and comparative analysis, it is shown that the proposed numerical model can be employed for various problems encountered in geotechnical engineering.