An Adaptive Sequential Fully Implicit Domain-Decomposition Solver
- 17 September 2021
- journal article
- research article
- Published by Society of Petroleum Engineers (SPE) in SPE Journal
- Vol. 27 (01), 566-578
- https://doi.org/10.2118/203991-pa
Abstract
Modern reservoir simulation must handle complex compositional fluid behavior, orders-of-magnitude variations in rock properties, and large velocity contrasts. We investigate how one can use nonlinear domain-decomposition preconditioning to combine sequential and fully implicit (FI) solution strategies to devise robust and highly efficient nonlinear solvers. A full simulation model can be split into smaller subdomains that each can be solved independently, treating variables in all other subdomains as fixed. In subdomains with weaker coupling between flow and transport, we use a sequential fully implicit (SFI) solution strategy, whereas regions with stronger coupling are solved with an FI method. Convergence to the FI solution is ensured by a global update that efficiently resolves long-range interactions across subdomains. The result is a solution strategy that combines the efficiency of SFI and its ability to use specialized solvers for flow and transport with the robustness and correctness of FI. We demonstrate the efficacy of the proposed method through a range of test cases, including both contrived setups to test nonlinear solver performance and realistic field models with complex geology and fluid physics. For each case, we compare the results with those obtained using standard FI and SFI solvers. This paper is published as part of the 2021 Reservoir Simulation Conference Special Issue.Keywords
This publication has 42 references indexed in Scilit:
- Trust-region based solver for nonlinear transport in heterogeneous porous mediaJournal of Computational Physics, 2013
- Domain decomposition strategies for nonlinear flow problems in porous mediaJournal of Computational Physics, 2013
- Results of the Brugge Benchmark Study for Flooding Optimization and History MatchingSPE Reservoir Evaluation & Engineering, 2010
- Unconditionally convergent nonlinear solver for hyperbolic conservation laws with S-shaped flux functionsJournal of Computational Physics, 2009
- Adaptive fully implicit multi-scale finite-volume method for multi-phase flow and transport in heterogeneous porous mediaJournal of Computational Physics, 2006
- Algorithm 832ACM Transactions on Mathematical Software, 2004
- Nonlinearly Preconditioned Inexact Newton AlgorithmsSIAM Journal on Scientific Computing, 2002
- Upstream Differencing for Multiphase Flow in Reservoir SimulationSIAM Journal on Numerical Analysis, 1991
- Mathematical Structure of Compositional Reservoir SimulationSIAM Journal on Scientific and Statistical Computing, 1989
- 4. Algebraic MultigridPublished by Society for Industrial & Applied Mathematics (SIAM) ,1987