Mojtaba Talebian, Rafid Al-Khoury and Lambertus J. Sluys
This paper aims to present a computationally efficient finite element model for the simulation of isothermal immiscible two-phase flow in a rigid porous media with a particular…
Abstract
Purpose
This paper aims to present a computationally efficient finite element model for the simulation of isothermal immiscible two-phase flow in a rigid porous media with a particular application to CO2 sequestration in underground formations. Focus is placed on developing a numerical procedure, which is effectively mesh-independent and suitable to problems at regional scales.
Design/methodology/approach
The averaging theory is utilized to describe the governing equations of the involved unsaturated multiphase flow. The level-set (LS) method and the extended finite element method (XFEM) are utilized to simulate flow of the CO2 plume. The LS is employed to trace the plume front. A streamline upwind Petrov-Galerkin method is adopted to stabilize possible occurrence of spurious oscillations due to advection. The XFEM is utilized to model the high gradient in the saturation field front, where the LS function is used for enhancing the weighting and the shape functions.
Findings
The capability of the proposed model and its features are evaluated by numerical examples, demonstrating its accuracy, stability and convergence, as well as its advantages over standard and upwind techniques. The study showed that a good combination between a mathematical model and a numerical model enables the simulation of complicated processes occurring in complicated and large geometry using minimal computational efforts.
Originality/value
A new computational model for two-phase flow in porous media is introduced with basic requirements for accuracy, stability, and convergence, which are met using relatively coarse meshes.