TY - GEN

T1 - Meshless Multi-Point Flux Approximation

AU - Lukyanov, Alexander A.

AU - Vuik, Cornelis

PY - 2017

Y1 - 2017

N2 - The reservoir simulation of the complex reservoirs with anisotropic permeability,which includes faults and non-orthogonal grids, with a fully discontinuous permeability tensor in the discretization is a major challenge. Several methods have already been developed and implemented within industry standard reservoir simulators for non-orthogonal grids (e.g., Multi-Point Flux Approximation (MPFA) “O” method). However, it has been noticed that some of the numerical methods for elliptic/parabolic equations may violate the maximum principle (i.e., lead to spurious oscillations), especially when the anisotropy is particularly strong. It has been found that the oscillations are closely related to the poor approximation of the pressure gradient in the flux computation. Therefore, proposed methods must correctly approximate underlying operators, satisfy a discrete maximum principle and have coercivity properties. Furthermore, the method must be robust and efficient. This paper presents the meshless multi-point flux approximation of second order elliptic operators containing a tensor coefficient. The method is based on a pressure gradient approximation commonly used in meshless methods (or Smoothed Particle Hydrodynamics method—SPH method). The proposed discretization schemes can be written as a sum of sparse positive semidefinite matrix and perturbation matrix. We show that convergence rates are retained as for finite difference methods O(hα), 1 ≤ α < 2, where h denotes the maximum particle spacing. The results are presented, discussed and future studies are outlined.

AB - The reservoir simulation of the complex reservoirs with anisotropic permeability,which includes faults and non-orthogonal grids, with a fully discontinuous permeability tensor in the discretization is a major challenge. Several methods have already been developed and implemented within industry standard reservoir simulators for non-orthogonal grids (e.g., Multi-Point Flux Approximation (MPFA) “O” method). However, it has been noticed that some of the numerical methods for elliptic/parabolic equations may violate the maximum principle (i.e., lead to spurious oscillations), especially when the anisotropy is particularly strong. It has been found that the oscillations are closely related to the poor approximation of the pressure gradient in the flux computation. Therefore, proposed methods must correctly approximate underlying operators, satisfy a discrete maximum principle and have coercivity properties. Furthermore, the method must be robust and efficient. This paper presents the meshless multi-point flux approximation of second order elliptic operators containing a tensor coefficient. The method is based on a pressure gradient approximation commonly used in meshless methods (or Smoothed Particle Hydrodynamics method—SPH method). The proposed discretization schemes can be written as a sum of sparse positive semidefinite matrix and perturbation matrix. We show that convergence rates are retained as for finite difference methods O(hα), 1 ≤ α < 2, where h denotes the maximum particle spacing. The results are presented, discussed and future studies are outlined.

UR - http://www.scopus.com/inward/record.url?scp=85018709342&partnerID=8YFLogxK

U2 - 10.1007/978-3-319-51954-8_5

DO - 10.1007/978-3-319-51954-8_5

M3 - Conference contribution

AN - SCOPUS:85018709342

SN - 978-3-319-51953-1

T3 - Lecture Notes in Computational Science and Engineering

SP - 67

EP - 84

BT - Meshfree Methods for Partial Differential Equations VIII

A2 - Griebel, Michael

A2 - Schweitzer, Marc Alexander

PB - Springer

CY - Cham

T2 - 8th International Workshop on Meshfree Methods for Partial Differential Equations, 2015

Y2 - 7 September 2015 through 9 September 2015

ER -