
Preliminaries Results on the Use of Non-Uniform Algebraic Dynamic Multilevel and Multiscale Method as a Preconditioner for the Numerical Simulation of Two-Phase Flows in Porous Media
Please login to view abstract download link
Nowadays, multiphase flow models often involve millions of degrees of freedom, resulting in large-scale problems. The systems of equations derived from their discretization are typically ill-conditioned, requiring significant computational resources and simulation times that can sometimes be prohibitive. Frequently, the oil industry demands the analysis of multiple scenarios, including the use of optimization strategies, history matching, and uncertainty quantification and propagation. A possible approach to address these challenges involves using multiscale techniques, either as reduced-order models or as preconditioners for high-fidelity models. In this work, we present the Algebraic Dynamic Multilevel method with Non-Uniform Resolution (NU-ADM), which leverages an Algebraic Multiscale (AMS) finite volume formulation with a classical Two-Point Flux Approximation (TPFA) for fully implicit simulation of two-phase flows in highly heterogeneous porous media. The pressure system is obtained using the Constrained Pressure Residual (CPR) method, which extracts the elliptic block from the Jacobian matrix. Multiscale operators are employed as preconditioners for solving this system. With this approach, the complete system of equations is solved by applying an iterative solver, such as the Generalized Minimum Residual Method (GMRES). To evaluate the performance of our strategy, we compare our results with other well-established preconditioning strategies, such as the Incomplete LU Factorization (ILU) and the Algebraic MultiGrid (AMG) method, which are widely used in industrial applications. For consistency, the same iterative solver is applied to both alternatives and our NU-ADM pressure approach. Our preliminary results are very promising, and in the near future, we also intend to explore the use of Multi-Stage Preconditioning (MSP) methods, incorporating multiple preconditioning stages and performing successive corrections. J. C. A. d. Santos, P. R. M. Lyra, J. P. R. d. Andrade, A. C. R. d. Souza, R. J. M. d. Lira Filho, D. K. E. d. Carvalho, An Algebraic Dynamic Multilevel and Multiscale Method with Non-Uniform Mesh Resolution and Adaptive Algebraic Multiscale Solver Operator for the Simulation of Two-Phase Flows in Highly Heterogeneous Petroleum Reservoirs, JCP, 462, 2022 M. Cusini, A. A. Lukyanov, J. Natvig, H. Hajibeygi, Constrained pressure residual multiscale (CPR-MS) method for fully implicit simulation of multiphase flow in porous media, JCP, 299, 2015