On the Performance of the Node Control Volume Finite Element Method for Modeling Multi-phase Fluid Flow in Heterogeneous Porous Media

Abstract

In this paper, we critique the performance of the node control volume finite element (NCVFE) method for modeling multi-phase fluid flow in heterogeneous media. The NCVFE method solves for the pressure at the vertices of elements and a control volume mesh is constructed around them. Material properties are defined on elements, while transport is simulated on the control volumes. These two meshes are not aligned producing inaccurate results and artificial fluid smearing when modeling multi-phase fluid flow in heterogeneous media. We perform numerical tests to quantify and visualize the extent of this artificial fluid smearing in domains with different material properties. The domains are composed of tetrahedron finite elements. Large artificial fluid smearing is observed in coarse meshes; however, it decreases with the increase in mesh resolution. These findings prompt the use of high-resolution meshes for the method and the need for development of novel numerical methods to address this unphysical flow.

A novel block non-symmetric preconditioner for mixed-hybrid finite-element-based flow simulations

“…This analysis confirms the connection between the static and dynamic strategies, possibly inspiring the selection of better static connections from the the dynamic patterns. In particular, we focus on the pairs (n add , n ent ) equal to (1,4) and (2,6). Notice in Figure 7 that the newly added faces are primarily located following the main flow direction.…”

“…Finally, the effect of the time step size on the EDFA preconditioner in a full transient simulation is assessed in Table 3. The preconditioner is built by means of the dynamic technique, with the setting (1,4) for the pair (n add , n ent ). The range investigated spans the interval [0.01,1000] days, with the associated χ ∞ parameter up to almost 600.…”

“…The degree of approximation of the overall numerical model is defined, first of all, by the underlying mathematical model, but also by the selected discretization scheme. Discretization schemes should handle effectively non K-orthogonal unstructured grids, as well as highly heterogeneous and anisotropic rock/fluid properties, frequently introduced as full-tensors in the model (see, for instance, [1] about the numerical issues, related to abrupt changes in permeability, in node control volume finite element discretizations).…”

Recent Advance and Prospect of Enhanced Oil Recovery Mechanisms in Reservoir Sandstone

“…The hydrophobicity of the organo-montmorillonite increases with the lengthening of the Gemini surfactant's alkyl chain, but the spacer's length has only a modest influence on wettability [91]. Nonionic surfactants with shorter hydrophilic units have also been shown to improve wettability alteration, which is highly dependent on the modifier structure [92]. The structure-function relationships of both surfactants and polymers need to be studied for EOR application to utilize them in EOR application entirely [90].…”

“…According to Schlumberger, 70% of conventional oil reserves in the Middle East are in carbonate reservoirs. Fractures of various sizes and lengths, ranging from extremely tiny to kilometre-wide fissures, make up carbonate reservoirs [92]. Interaction between cracks and rock layers via capillary and gravity forces is the critical component that governs production in carbonate reservoirs.…”