Methods and systems to reduce or eliminate numerical oscillations in solutions that occur when using conventional MPFA when modeling flow in a reservoir are provided. The technique may be referred to as enriched multi-point flux approximation (EMPFA) and may be used to improve the consistency and accuracy in constructing pressure interpolations in cells for the purpose of determining flux equations used in predicting flow in a reservoir.