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 -