Interpolation function is a key parameter for numerical simulation using finite-difference, finite-element, and finite-volume methods, especially when the advective terms of the conservation equations are considered. Due to the flow orientation, a first-order interpolation scheme such as upwind introduces a considerable degree of numerical diffusion in the numerical solution. In this work, we present a second-order total variation diminishing scheme in conjunction with 3D compositional reservoir simulation using the element-based finite-volume method (EbFVM). The results of several case studies using the hexahedron element are shown in terms of oil, water, and gas production, as well as saturation field.

In reducing the grid orientation effect for the numerical solution of partial differential equations, interpolation functions play an important role when the advective transport of the governing equations is considered. This is due to the fact that, in general, the unknowns are evaluated in the vertices of the elements and such properties must be extrapolated to inner parts of the elements. First-order schemes, such as upwind, are the easiest methods to use for performing the extrapolation of the properties. However, such methods introduce a large amount of numerical diffusion in the solution. A few higher-order interpolation schemes, on the other hand, are capable of providing solutions free of numerical diffusion, increasing the accuracy of the method and reducing the computational efforts required. In this work, we investigate the TVD interpolation scheme for three-dimensional unstructured grids in conjunction with Element-based Finite Volume Method (EbFVM) using four types of elements: hexahedron, tetrahedron, prism and pyramid.

No abstract

No abstract