Multipoint flux approximation (MPFA) methods were introduced to solve control-volume formulations on general grids. Although these methods are general in the sense that they may be applied to any grid, their convergence and monotonicity properties vary. We introduce a new MPFA method for quadrilateral grids termed the L-method. This method seeks to minimize the number of entries in the flux stencils, while honoring uniform flow fields. The methodology is valid for general media. For homogeneous media and uniform grids in two dimensions, this method has four-point flux stencils and seven-point cell stencils, whereas the MPFA O-methods have six-point flux stencils and nine-point cell stencils. The reduced stencil of the L-method appears as a consequence of adapting the method to the closest neighboring cells, or equivalently, to the dominating principal direction of anisotropy. We have tested the convergence and monotonicity properties for this method and compared it with the O-methods. For moderate grids, the convergence rates are the same, but for rough grids with large aspect ratios, the convergence of the O-methods is lost, while the Lmethod converges with a reduced convergence rate. Also, the L-method has a larger monotonicity range than the O-methods. For homogeneous media and uniform parallelogram grids, the matrix of coefficients is an M-matrix whenever the method is monotone. For strongly nonmonotone cases, the oscillations observed for the O-methods are almost removed for the L-method. Instead, extrema on no-flow boundaries are observed. These undesired solutions, which only occur for parameters not common in applications, should be avoided by requiring that the previously derived monotonicity conditions are satisfied. For local grid refinements, test runs indicate that the L-method yields almost optimal solutions, and that the solution is considerably better than the solutions obtained by the O-methods. The efficiency of the linear solver is in many cases better for the L-method than for the O-methods. This is due to lower condition number and a reduced number of entries in the matrix of coefficients.
Robustness of numerical methods for multiphase flow problems in porous media is important for development of methods to be used in a wide range of applications. Here, we discuss monotonicity for a simplified problem of single-phase flow, but where the simulation grids and media are allowed to be general, posing challenges to control-volume methods. We discuss discrete formulations of the maximum principle and derive sufficient criteria for discrete monotonicity for arbitrary nine-point control-volume discretizations for conforming quadrilateral grids in 2D. These criteria are less restrictive than the M-matrix property. It is shown that it is impossible to construct nine-point methods which unconditionally satisfy the monotonicity criteria when the discretization satisfies local conservation and exact reproduction of linear potential fields. Numerical examples are presented which show the validity of the criteria for monotonicity. Further, the impact of nonmonotonicity is studied. Different behavior for different discretization methods is illuminated, and simple ideas are presented for improvement in terms of monotonicity.
Mathematics Subject Classification (1991)
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.