Bruno Ramon Batista Fernandes

Investigation of Several Interpolation Functions for Unstructured Meshes in Conjunction with Compositional Reservoir Simulation

One of the key parameters in numerical simulation of fluid flow problems is the interpolation function used for adjusting the numerical value from the position where it is evaluated to the interface of the control volumes; for instance, if the finite-volume method is employed. In this article, we present an investigation of several interpolation functions in conjunction with the element based finite-volume method (EbFVM) using unstructured triangular meshes. We investigate the mass weighted upwind (MWU) and a modified version of this method, the upwind scheme, a streamline based upwind scheme, and a limited second order upwind scheme. These interpolation functions are implemented in a compositional simulator using vertex cell unstructured grids. The accuracy of these interpolation functions are evaluated for several case studies. For some of them, we also compare the results with the available analytical solutions.

A 3D Total Variation Diminishing Scheme for Compositional Reservoir Simulation Using the Element-Based Finite-Volume Method

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.