Volume 53, Number 4, July-August 2019
|Page(s)||1125 - 1156|
|Published online||05 July 2019|
Combined face based and nodal based discretizations on hybrid meshes for non-isothermal two-phase Darcy flow problems⋆
Université Côte d’Azur, CNRS, Inria team Coffee, LJAD, Nice, France
2 BRGM, scientific and Technical Center, 3 avenue Claude Guillemin, BP 36009, 45060 Orléans Cedex 2, France
3 CSTJF, TOTAL S.A. – Avenue Larribau, 64018 Pau, France
* Corresponding author: firstname.lastname@example.org
Accepted: 18 February 2019
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.
Mathematics Subject Classification: 65M08 / 65M12 / 76S05
Key words: Finite Volume / gradient discretization / Darcy flow / two-phase Darcy flow / hybrid meshes
© The authors. Published by EDP Sciences, SMAI 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.