A family of multi-point flux approximation schemes for general element types in two and three dimensions with convergence performance

SUMMARY A family of flux-continuous, locally conservative, control-volume-distributed multi-point flux approximation (CVD-MPFA) schemes has been developed for solving the general geometry-permeability tensor pressure equation on structured and unstructured grids. These schemes are applicable to the full-tensor pressure equation with generally discontinuous coefficients and remove the O(1) errors introduced by standard reservoir simulation schemes when applied to full-tensor flow approximation. The family of flux-continuous schemes is characterized by a quadrature parameterization. Improved numerical convergence for the family of CVD-MPFA schemes using the quadrature parameterization has been observed for structured and unstructured grids in two dimensions. The CVD-MPFA family cell-vertex formulation is extended to classical general element types in 3-D including prisms, pyramids, hexahedra and tetrahedra. A numerical convergence study of the CVD-MPFA schemes on general unstructured grids comprising of triangular elements in 2-D and prismatic, pyramidal, hexahedral and tetrahedral shape elements in 3-D is presented. Copyright © 2011 John Wiley & Sons, Ltd.

[1]  Todd Arbogast,et al.  Logically rectangular mixed methods for Darcy flow on general geometry , 1995 .

[2]  Bradley T. Mallison,et al.  A compact multipoint flux approximation method with improved robustness , 2008 .

[3]  Michael G. Edwards,et al.  Higher‐resolution hyperbolic‐coupled‐elliptic flux‐continuous CVD schemes on structured and unstructured grids in 3‐D , 2006 .

[4]  Ivar Aavatsmark,et al.  Discretization on Unstructured Grids for Inhomogeneous, Anisotropic Media. Part I: Derivation of the Methods , 1998, SIAM J. Sci. Comput..

[5]  Michael G. Edwards,et al.  A family of MPFA finite-volume schemes with full pressure support for the general tensor pressure equation on cell-centered triangular grids , 2011, J. Comput. Phys..

[6]  Patrick Jenny,et al.  A Finite-Volume Method with Hexahedral Multiblock Grids for Modeling Flow in Porous Media , 2002 .

[7]  Michael G. Edwards,et al.  Quasi M-Matrix Multifamily Continuous Darcy-Flux Approximations with Full Pressure Support on Structured and Unstructured Grids in Three Dimensions , 2011, SIAM J. Sci. Comput..

[8]  I. Aavatsmark,et al.  Numerical convergence of the MPFA O‐method and U‐method for general quadrilateral grids , 2006 .

[9]  Michael G. Edwards,et al.  A quasi-positive family of continuous Darcy-flux finite-volume schemes with full pressure support , 2008, J. Comput. Phys..

[10]  Louis J. Durlofsky,et al.  Unstructured grid optimization for improved monotonicity of discrete solutions of elliptic equations with highly anisotropic coefficients , 2006, J. Comput. Phys..

[11]  Michael G. Edwards,et al.  Higher‐resolution hyperbolic‐coupled‐elliptic flux‐continuous CVD schemes on structured and unstructured grids in 2‐D , 2006 .

[12]  Michael G. Edwards.  Split Full Tensor Discretization Operators for Structured and Unstructured Grids in Three Dimensions , 2001 .

[13]  Mayur Pal,et al.  Positive-definite q-families of continuous subcell Darcy-flux CVD(MPFA) finite-volume schemes and the mixed finite element method , 2008 .

[14]  Daniil Svyatskiy,et al.  Monotone finite volume schemes for diffusion equations on unstructured triangular and shape-regular polygonal meshes , 2007, J. Comput. Phys..

[15]  G. T. Eigestad,et al.  On the convergence of the multi-point flux approximation O-method: Numerical experiments for discontinuous permeability , 2005 .

[16]  Mayur Pal,et al.  Quasimonotonic Continuous Darcy-Flux Approximation for General 3D Grids of Any Element Type , 2007 .

[17]  Johannes Mykkeltveit,et al.  Symmetric Positive Definite Flux-Continuous Full-Tensor Finite-Volume Schemes on Unstructured Cell-Centered Triangular Grids , 2008, SIAM J. Sci. Comput..

[18]  Michael G. Edwards,et al.  Double-families of quasi-positive Darcy-flux approximations with highly anisotropic tensors on structured and unstructured grids , 2010, J. Comput. Phys..

[19]  Mayur Pal,et al.  Convergence study of a family of flux‐continuous, finite‐volume schemes for the general tensor pressure equation , 2006 .

[20]  Christophe Le Potier,et al.  Schéma volumes finis pour des opérateurs de diffusion fortement anisotropes sur des maillages non structurés , 2005 .

[21]  Michael G. Edwards,et al.  Unstructured, Control-Volume Distributed, Full-Tensor Finite-Volume Schemes with Flow Based Grids , 2002 .

[22]  Ragnar Winther,et al.  Convergence of multipoint flux approximations on quadrilateral grids , 2006 .

[23]  Mayur Pal,et al.  Non‐linear flux‐splitting schemes with imposed discrete maximum principle for elliptic equations with highly anisotropic coefficients , 2011 .

[24]  Michael G. Edwards,et al.  M-Matrix Flux Splitting for General Full Tensor Discretization Operators on Structured and Unstructured Grids , 2000 .

[25]  Michael G. Edwards,et al.  Finite volume discretization with imposed flux continuity for the general tensor pressure equation , 1998 .

[26]  Khalid Aziz,et al.  A Control Volume Scheme for Flexible Grids in Reservoir Simulation , 1997 .

[27]  Ivar Aavatsmark,et al.  Numerical Convergence of the MPFA O-Method for General Quadrilateral Grids in Two and Three Dimensions , 2006 .