4.7 Article

A hybrid finite element model for non-isothermal two-phase flow in deformable porous media

Journal

COMPUTERS AND GEOTECHNICS
Volume 135, Issue -, Pages -

Publisher

ELSEVIER SCI LTD
DOI: 10.1016/j.compgeo.2021.104199

Keywords

Multiphase flow; Geomechanic; Thermal flow; THM coupling; Hybrid solution

Funding

  1. Research Council of Norway through its Centers of Excellence funding scheme [262644]
  2. German Research Foundation (DFG) through the Collaborative Research Center [SFB 837]

Ask authors/readers for more resources

This paper presents a numerical method for modeling coupled thermo-hydro-mechanical processes in porous media saturated with two immiscible fluids. The method combines GFEM, CVFE, and SUCVFE approaches to ensure mass conservation and eliminate oscillations, improving accuracy and efficiency in simulation. The proposed method has been justified through numerical examples and comparison with analytical or numerical solutions.
This paper presents a numerical method to model the coupled thermo-hydro-mechanical (THM) processes in porous media saturated with two immiscible fluids. The basic equations of the system have been derived based on the averaging theory, considering skeleton deformation, two-phase fluid flow, and heat transport. As applying the standard Galerkin finite element method (GFEM) to solve this system of partial differential equations may lead to oscillatory results for saturation and temperature profiles, a hybrid numerical solution is proposed. In this frame, the GFEM is combined with a control volume based finite element (CVFE) approach, and a streamline upwind control volume finite element (SUCVFE) scheme, respectively for the mechanical, hydraulic and thermal part of the system. The CVFEM has been adopted to provide a smooth saturation profile by ensuring local mass conservation, while the streamline upwind scheme has been applied to remove the spurious temperature oscillation by adding stabilizing terms to the thermal part of the system. The CVFE and SUCVFE formulations have been derived using a similar approach as the standard FE practice in the context of weighted residual technique, but using different weighting functions. This will significantly facilitate the implementation of the proposed model in existing FE codes. Accuracy and efficiency of the proposed method have been justified using several numerical examples and comparing the results with available analytical or numerical solutions.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available