The example explains how to use the MPFA-O solver to solve a problem with homogeneous, anisotropic permeability. The resulting pressure solutions are compared and contrasted to similar solutions computed by the TPFA and a mimetic method.
The MPFA-O scheme is consistent on grids that are not necessarily K-orthogonal. This example compares and contrasts grid-orientation effects for the MPFA-O, mimetic, and TPFA solvers for a single-phase flow problem with isotropic permeability represented on a curvilinear grid in which most cells are not K-orthogonal.
To get a consistent discretization for anisotropic permeablities and general polyhedral grids, multipoint flux-approximation method use, as the name says, more than two points to approximate the flux across each inter-cell face. The mpfa module implements one such method, the MPFA-O method . The implementation utilizes the fact that some variants of the method can be formulated as a mimetic method. In the resulting local-flux mimetic formulation [2,3] of the MPFA-O method, each face in the grid is subdivided into a set of subfaces, one subface per node that makes up the face. The inner product of the local-flux mimetic method gives exact result for linear flow and is block diagonal with respect to the faces corresponding to each node of the cell, but it is not symmetric. The block-diagonal property makes it possible to reduce the system into a cell-centered discretisation for the cell pressures. This naturally leads to a method for calculating the MPFA transmissibilities.
The module offers two basic routines
The incompMPFA solver is compatible with the incompressible fluid models and transport solvers from the incomp module.
See also: two-point and mimetic discretizations
Published December 8, 2016