| Title |
Stable Multipoint Flux Approximation (MPFA) scheme for anisotropic porous media / |
| Authors |
Makauskas, Pijus ; Pal, Mayur |
| DOI |
10.1016/j.aej.2024.11.108 |
| Full Text |
|
| Is Part of |
Alexandria engineering journal.. Amsterdam : Elsevier. 2025, vol. 117, p. 418-429.. ISSN 1110-0168. eISSN 2090-2670 |
| Keywords [eng] |
numerical methods ; darcy flow ; anisotropy ; non K-orthogonal grids ; multipoint flux approximation ; monotonicity ; maximum principle ; spurious oscillations |
| Abstract [eng] |
This paper focuses on enhancing the multipoint flux approximation (MPFA) method for subsurface reservoir simulation with heterogeneous anisotropic permeability with discontinuities at grid block level. MPFA is a numerical method used in reservoir simulation to accurately model fluid flow through porous media by considering flux interactions at multiple points within each computational cell. MPFA aims to eliminate O ( h 1 ) inconsistency errors in the widely used two-point method (TPFA), however, it can suffer from spurious oscillations, especially in case of high anisotropy scenarios on non K-orthogonal grids. Various stability-improving schemes have been proposed, including optimal quadrature support schemes and non-linear flux-splitting techniques, which are computationally resource-intensive and do not always result in a fully oscillation free solution for a range of permeability anisotropies. This paper addresses the stability challenges in the standard MPFA-O method by enhancing the MPFA approach for highly heterogeneous, anisotropic domains represented by non K-orthogonal grids. Our new method implements a local rotation transformation, effectively reducing spurious oscillations. Results from the enhanced method demonstrate improved performance of the numerical solution, eliminating oscillatory behavior and maintaining the maximum principle and monotonicity. |
| Published |
Amsterdam : Elsevier |
| Type |
Journal article |
| Language |
English |
| Publication date |
2025 |
| CC license |
|
