MPFA methods were introduced to solve control-volume formulations on general simulation grids for porous media flow. While these methods are general in the sense that they may be applied to any grid, their convergence properties vary.
An important property for multiphase flow is the monotonicity of the numerical elliptic operator. In a recent paper , conditions for monotonicity on quadrilateral grids have been developed. These conditions indicate that MPFA formulations which lead to smaller flux stencils are desirable for grids with high aspect ratio or severe skewness and for media with strong anisotropy or strong heterogeneity. The ideas were recently pursued in  where the L-method was introduced for general media in 2D. For homogeneous media and uniform grid, this method has four-point flux stencils and seven-point cell stencils in two dimensions. The reduced stencils appear as a consequence of adapting the method to the closest neighboring cells
Here, we extend the ideas for discretization on 3D grids, and ideas and results are shown for both conforming and nonconforming grids. The ideas are particularly desirable for simulation grids which contain faults and local grid refinement.
We present numerical results herein which include convergence results for single phase flow on challenging grids in 2D and 3D, and some simple two-phase results. Also, we compare the L-method with the O-method.