Combined face based and nodal based discretizations on hybrid meshes for non-isothermal two-phase Darcy flow problems

Beaude, L.; Masson, R.; Lopez, S.; Samier, P.

ESAIM: Mathematical Modelling and Numerical Analysis , Tome 53 (2019) no. 4, pp. 1125-1156.

Suite au passage du modèle économique de la revue en S20, le texte intégral des articles des années 2019 et 2020 est accessible uniquement sur le site de l'éditeur et est réservé aux abonnés.

In the last 20 years many discretization schemes have been developed to approximate the Darcy fluxes on polyhedral cells in heterogeneous anisotropic porous media. Among them, we can distinguished cell based approaches like the Two Point Flux Approximation (TPFA) or the Multi Point Flux Approximation (MPFA) schemes, face based approaches like the Hybrid Finite Volume (HFV) scheme belonging to the family of Hybrid Mimetic Mixed methods and nodal based discretizations like the Vertex Approximate Gradient (VAG) scheme. They all have their own drawbacks and advantages which typically depend on the type of cells and on the anisotropy of the medium. In this work, we propose a new methodology to combine the VAG and HFV discretizations on arbitrary subsets of cells or faces in order to choose the best suited scheme in different parts of the mesh. In our approach the TPFA discretization is considered as an HFV discretization for which the face unknowns can be eliminated. The coupling strategy is based on a node to face interpolation operator at the interfaces which must be chosen to ensure the consistency, the coercivity and the limit conformity properties of the combined discretization. The convergence analysis is performed in the gradient discretization framework and convergence is proved for arbitrary cell or face partitions of the mesh. For face partitions, an additional stabilisation local to the cell is required to ensure the coercivity while for cell partitions no additional stabilisation is needed. The framework preserves at the interface the discrete conservation properties of the VAG and HFV schemes with fluxes based on local to each cell transmissibility matrices. This discrete conservative form allows to naturally extend the VAG and HFV discretizations of two-phase Darcy flow models to the combined VAG–HFV schemes. The efficiency of our approach is tested for single phase and immiscible two-phase Darcy flows on 3D meshes using a combination of the HFV and VAG discretizations as well as for non-isothermal compositional liquid gas Darcy flows on a vertical 2D cross-section of the Bouillante geothermal reservoir (Guadeloupe) using a combination of the TPFA and VAG discretizations.

[1] I. Aavatsmark, T. Barkve, O. Boe and T. Mannseth, Discretization on unstructured grids for inhomogeneous, anisotropic media. Part I: derivation of the methods. SIAM J. Sci. Comput. 19 (1998). | MR | Zbl

[2] A.S. Abushaikha, D.V. Voskov and H.A. Tchelepi, Fully implicit mixed-hybrid finite element discretization for general purpose subsurface reservoir simulation. J. Comput. Phys. 346 (2017) 514–538. | DOI | MR | Zbl

[3] O. Angelini, C. Chavant, E. Chénier, R. Eymard and S. Granet, Finite volume approximation of a diffusion-dissolution model and application to nuclear waste storage. Math. Comput. Simul. 81 (2011) 2001–2017. | DOI | MR | Zbl

[4] K. Aziz and A. Settari, Petroleum Reservoir Simulation. Applied Science Publishers (1979).

[5] L. Beaude, K. Brenner, S. Lopez, R. Masson and F. Smai, Non-isothermal compositional liquid gas Darcy flow: formulation, soil-atmosphere boundary condition and application to high energy geothermal simulations, Preprint https://hal.archives-ouvertes.fr/hal-01702391 (2018). | MR | Zbl

[6] R. Bertani, Geothermal power generation in the world 2010–2014 update report. Geothermics 60 (2016) 31–43. | DOI

[7] K. Brenner and R. Masson, Convergence of a vertex centered discretization of two-phase Darcy flows on general meshes. Int. J. Finite Volume Methods (2013). | MR | Zbl

[8] H. Class, R. Helmig and P. Bastian, Numerical simulation of non-isothermal multiphase multicomponent processes in porous media: 1. An efficient solution technique. Adv. Water Res. 25 (2002) 533–550. | DOI

[9] K.H. Coats, Implicit compositional simulation of single-porosity and dual-porosity reservoirs. In: SPE Symposium on Reservoir Simulation. Society of Petroleum Engineers (1989). | DOI

[10] J. Droniou, R. Eymard, T. Gallouët and C. Guichard, R. Herbin, The gradient discretization method: a framework for the discretization of linear and nonlinear elliptic and parabolic problems. Preprint https://hal.archives-ouvertes.fr/hal-01382358 (2016).

[11] J. Droniou, R. Eymard, T. Gallouët and R. Herbin, A unified approach to mimetic finite difference, hybrid finite volume and mixed finite volume methods. Math. Models Methods Appl. Sci. 20 (2010) 265–295. | DOI | MR | Zbl

[12] J. Droniou, R. Eymard, T. Gallouët and R. Herbin, Gradient schemes: a generic framework for the discretisation of linear, nonlinear and nonlocal elliptic and parabolic equations. Math. Models Methods Appl. Sci. 23 (2013) 2395–2432. | DOI | MR | Zbl

[13] J. Droniou, Finite volume schemes for diffusion equations: introduction and review of modern methods. Math. Models Methods Appl. Sci. 24 (2014) 1575–1619. | DOI | MR | Zbl

[14] M.G. Edwards and C.F. Rogers, Finite volume discretization with imposed flux continuity for the general tensor pressure equation. Comput. Geosci. 2 (1998) 250–290. | DOI | MR | Zbl

[15] R. Eymard, C. Guichard and R. Herbin, Small-stencil 3D schemes for diffusive flows in porous media, ESAIM: M2AN 46 (2010) 265–290. | DOI | Numdam | MR | Zbl

[16] R. Eymard, T. Gallouët and R. Herbin, Discretization of heterogeneous and anisotropic diffusion problems on general nonconforming meshes SUSHI: a scheme using stabilisation and hybrid interfaces. IMA J. Numer. Anal. 30 (2010) 1009–1043. | DOI | MR | Zbl

[17] R. Eymard, C. Guichard and R. Herbin, Benchmark 3D: the VAG scheme, in Finite Volumes for Complex Applications VI – Problems and Perspectives, edited by J. Fort, J. Furst, J. Halama, R. Herbin and F. Hubert. In Vol 2 of Springer Proceedings in Mathematics (2011) 213–222.

[18] R. Eymard, C. Guichard, R. Herbin and R. Masson, Vertex centred discretization of compositional multiphase darcy flows on general meshes. Comput. Geosci. 16 (2012) 987–1005. | DOI | MR

[19] R. Huber and R. Helmig, Node-centered finite volume discretizations for the numerical simulation of multi-phase flow in heterogeneous porous media. Comput. Geosci. 4 (2000) 141–164. | DOI | MR | Zbl

[20] S. Lacroix, Y.V. Vassilevski and M.F. Wheeler, Decoupling preconditioners in the implicit parallel accurate reservoir simulator (IPARS). Numer. Linear Algebra Appl. 8 (2001) 537–549. | DOI | MR | Zbl

[21] A. Lauser, C. Hager, R. Helmig and B. Wohlmuth, A new approach for phase transitions in miscible multi-phase flow in porous media. Adv. Water Res. 34 (2011) 957–966. | DOI

[22] D.W. Peaceman, Fundamentals of Numerical Reservoir Simulations. Elsevier (1977).

[23] P. Samier and R. Masson, Implementation of a vertex centered method inside an industrial reservoir simulator – Practical issues and comprehensive comparison with CPG and PEBI grid models on a field case. SPE Reservoir Simulation Symposium (2016).

[24] R. Scheichl, R. Masson and J. Wendebourg, Decoupling and block preconditioning for sedimentary basin simulations. Comput. Geosci. 7 (2003) 295–318. | DOI | MR | Zbl

[25] E. Schmidt, Properties of Water and Steam in S.I. units. Springer-Verlag (1969).

[26] F. Xing, R. Masson and S. Lopez, Parallel numerical modeling of hybrid-dimensional compositional non-isothermal Darcy flows in fractured porous media. J. Comput. Phys. 345 (2017) 637–664. | DOI | MR | Zbl