Abstract: An apparatus and method are provided for solving a non-linear S-shaped function F=f(S) which is representative of a property S in a physical system, such saturation in a reservoir simulation. A Newton iteration (T) is performed on the function f(S) at Sv to determine a next iterative value Sv+1. It is then determined whether Sv+1 is located on the opposite side of the inflection point Sc from Sv. If Sv+1 is located on the opposite side of the inflection point from Sv, then Sv+1 is set to Sl, a modified new estimate. The modified new estimate, Sl, is preferably set to either the inflection point, Sc, or to an average value between Sv and Sv+1, i.e., Sl=0.5(Sv+Sv+1). The above steps are repeated until Sv+1 is within the predetermined convergence criteria. Also, solution algorithms are described for two-phase and three-phase flow with gravity and capillary pressure.

Abstract: A multi-scale finite-volume (MSFV) method to solve elliptic problems with a plurality of spatial scales arising from single or multi-phase flows in porous media is provided. The method efficiently captures the effects of small scales on a coarse grid, is conservative, and treats tensor permeabilities correctly. The underlying idea is to construct transmissibilities that capture the local properties of a differential operator. This leads to a multi-point discretization scheme for a finite-volume solution algorithm. Transmissibilities for the MSFV method are preferably constructed only once as a preprocessing step and can be computed locally.

Abstract: A multi-scale finite-volume (MSFV) method to solve elliptic problems with a plurality of spatial scales arising from single or multi-phase flows in porous media is provided. Two sets of locally computed basis functions are employed. A first set of basis functions captures the small-scale heterogeneity of the underlying permeability field, and it is computed to construct the effective coarse-scale transmissibilities. A second set of basis functions is required to construct a conservative fine-scale velocity field. The method efficiently captures the effects of small scales on a coarse grid, is conservative, and treats tensor permeabilities correctly. The underlying idea is to construct transmissibilities that capture the local properties of a differential operator. This leads to a multi-point discretization scheme for a finite-volume solution algorithm. Transmissibilities for the MSFV method are preferably constructed only once as a preprocessing step and can be computed locally.

Abstract: A multi-scale finite-volume (MSFV) method to solve elliptic problems with a plurality of spatial scales arising from single or multi-phase flows in porous media is provided. Two sets of locally computed basis functions are employed. A first set of basis functions captures the small-scale heterogeneity of the underlying permeability field, and it is computed to construct the effective coarse-scale transmissibilities. A second set of basis functions is required to construct a conservative fine-scale velocity field. The method efficiently captures the effects of small scales on a coarse grid, is conservative, and treats tensor permeabilities correctly. The underlying idea is to construct transmissibilities that capture the local properties of a differential operator. This leads to a multi-point discretization scheme for a finite-volume solution algorithm. Transmissibilities for the MSFV method are preferably constructed only once as a preprocessing step and can be computed locally.

Abstract: A multi-scale finite-volume (MSFV) method to solve elliptic problems with a plurality of spatial scales arising from single or multi-phase flows in porous media is provided. Two sets of locally computed basis functions are employed. A first set of basis functions captures the small-scale heterogeneity of the underlying permeability field, and it is computed to construct the effective coarse-scale transmissibilities. A second set of basis functions is required to construct a conservative fine-scale velocity field. The method efficiently captures the effects of small scales on a coarse grid, is conservative, and treats tensor permeabilities correctly. The underlying idea is to construct transmissibilities that capture the local properties of a differential operator. This leads to a multi-point discretization scheme for a finite-volume solution algorithm. Transmissibilities for the MSFV method are preferably constructed only once as a preprocessing step and can be computed locally.