Block constrained pressure residual preconditioning for two-phase flow in porous media by mixed hybrid finite elements

Stefano Nardean,Massimiliano Ferronato,Ahmad Abushaikha

2023-03-24

Abstract:This work proposes an original preconditioner that couples the Constrained Pressure Residual (CPR) method with block preconditioning for the efficient solution of the linearized systems of equations arising from fully implicit multiphase flow models. This preconditioner, denoted as Block CPR (BCPR), is specifically designed for Lagrange multipliers-based flow models, such as those generated by Mixed Hybrid Finite Element (MHFE) approximations. An original MHFE-based formulation of the two-phase flow model is taken as a reference for the development of the BCPR preconditioner, in which the set of system unknowns comprises both element and face pressures, in addition to the cell saturations, resulting in a $3 \times 3$ block-structured Jacobian matrix with a $2 \times 2$ inner pressure problem. The CPR method is one of the most established techniques for reservoir simulations, but most research focused on solutions for Two-Point Flux Approximation (TPFA)-based discretizations that do not readily extend to our problem formulation. Therefore, we designed a dedicated two-stage strategy, inspired by the CPR algorithm, where a block preconditioner is used for the pressure part with the aim at exploiting the inner $2 \times 2$ structure. The proposed preconditioning framework is tested by an extensive experimentation, comprising both synthetic and realistic applications in Cartesian and non-Cartesian domains.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the inefficiency in solving linear equations during the two - phase flow simulation in porous media. Specifically, the paper proposes a new pre - processing method - Block Constrained Pressure Residual (BCPR) - to improve the efficiency of solving linearized system equations in the two - phase flow model discretized by the Mixed Hybrid Finite Element (MHFE) method.

The traditional Constrained Pressure Residual (CPR) pre - processing method is mainly designed for the discretization method based on Two - Point Flux Approximation (TPFA), and it has a poor effect on the Jacobian matrix with a complex block structure generated by MHFE. Therefore, the paper proposes a two - stage strategy specifically for the MHFE method, which improves the CPR algorithm by using block pre - processing of the pressure part of the internal 2×2 block structure, aiming to solve this problem.

The main contributions of the paper are as follows:
1. **Design of the BCPR pre - processing method**: This method combines block pre - processing technology and the CPR method, and is especially suitable for two - phase flow models containing element and face pressure variables.
2. **Proposal of an effective Schur complement approximation method**: Through the Explicit Decoupling Factor Approximation (EDFA) technology, the Schur complement is effectively approximated, which is crucial for dealing with complex 2×2 asymmetric block matrices.
3. **Provision of detailed numerical experiment results**: Through synthetic and practical application tests on Cartesian and non - Cartesian domains, the effectiveness and robustness of the BCPR pre - processing method are verified.

In conclusion, the goal of the paper is to improve the efficiency of solving linearized system equations in the two - phase flow model by developing and implementing the BCPR pre - processing method, thereby reducing the computational time and resource consumption while maintaining the accuracy of the solution.