Multiscale finite-volume method for density-driven flow in porous media

The multiscale finite-volume (MSFV) method has been developed to solve multiphase flow problems on large and highly heterogeneous domains efficiently. It employs an auxiliary coarse grid, together with its dual, to define and solve a coarse-scale pressure problem. A set of basis functions, which are local solutions on dual cells, is used to interpolate the coarse-grid pressure and obtain an approximate fine-scale pressure distribution. However, if flow takes place in presence of gravity (or capillarity), the basis functions are not good interpolators. To treat this case correctly, a correction function is added to the basis function interpolated pressure. This function, which is similar to a supplementary basis function independent of the coarse-scale pressure, allows for a very accurate fine-scale approximation. In the coarse-scale pressure equation, it appears as an additional source term and can be regarded as a local correction to the coarse-scale operator: It modifies the fluxes across the coarse-cell interfaces defined by the basis functions. Given the closure assumption that localizes the pressure problem in a dual cell, the derivation of the local problem that defines the correction function is exact, and no additional hypothesis is needed. Therefore, as in the original MSFV method, the only closure approximation is the localization assumption. The numerical experiments performed for density-driven flow problems (counter-current flow and lock exchange) demonstrate excellent agreement between the MSFV solutions and the corresponding fine-scale reference solutions.     

Node-centered finite volume discretizations for the numerical simulation of multiphase flow in heterogeneous porous media

Three node-centered finite volume discretizations for multiphase porous media flow are presented and compared. By combination of these methods two additional discretization methods are generated. The ability of these schemes to describe flows at textural interfaces of different geologic formations is investigated. It was found that models with nonzero-entry pressures for the capillary pressure-saturation relationship in conjunction with the Box discretization may give rise to spurious oscillations for flows around low permeable lenses. Furthermore, the applicability and sensitivity of the discretization methods with regard to the used computational grids is discussed. The schemes are used for the numerical study of two-phase flow in porous media with zones of different material properties. This revised version was published online in July 2006 with corrections to the Cover Date.     

A finite volume approach with local adaptation scheme for the simulation of free surface flow in porous media

We present a method for solving steady‐state flow with a free surface in porous media. This method is based on a finite volume approach and is halfway between a fixed and an adaptive mesh method, taking advantage of both approaches: computational efficiency and localization accuracy. Most of the mesh remains fixed during the iterative process, while the cells in contact with the free surface (free surface cells) are being reshaped. Based on this idea, we developed two methods. In the first one, only the volumes of the free surface cells are adapted. In the second one, the computational nodes of the free surface cells are relocated exactly at the free surface. Both adaptations are designed for a better application of the free surface boundary conditions. Implementation details are given on a regular finite volume mesh for the case of homogeneous and heterogeneous rectangular dams in 2D and 3D. Accuracy and convergence properties of the proposed approach are demonstrated by comparison with an analytical solution and with existing references. Copyright © 2011 John Wiley & Sons, Ltd.     

Construction of three-phase data to model multiphase flow in porous media: Comparing an optimization approach to the finite element approach

Multiphase flow modelling is a major issue in the assessment of groundwater pollution. Three-phase flows are commonly governed by mathematical models that associate a pressure equation with two saturation equations. These equations involve a number of secondary variables that reflect the fluid behaviour in a porous medium. To improve the computational efficiency of multiphase flow simulators, several simplified reformulations of three-phase flow equations have been proposed. However, they require the construction of new secondary variables adapted to the reformulated flow equations. In this article, two different approaches are compared to quantify these variables. A numerical example is given for a typical fine sand.     

Numerical modeling of multiphase fluid flow in deforming porous media: A comparison between two- and three-phase models for seismic analysis of earth and rockfill dams

In this paper, a fully coupled numerical model is presented for the finite element analysis of the deforming porous medium interacting with the flow of two immiscible compressible wetting and non-wetting pore fluids. The governing equations involving coupled fluid flow and deformation processes in unsaturated soils are derived within the framework of the generalized Biot theory. The displacements of the solid phase, the pressure of the wetting phase and the capillary pressure are taken as the primary unknowns of the present formulation. The other variables are incorporated into the model using the experimentally determined functions that define the relationship between the hydraulic properties of the porous medium, i.e. saturation, relative permeability and capillary pressure. It is worth mentioning that the imposition of various boundary conditions is feasible notwithstanding the choice of the primary variables. The modified Pastor–Zienkiewicz generalized constitutive model is introduced into the mathematical formulation to simulate the mechanical behavior of the unsaturated soil. The accuracy of the proposed mathematical model for analyzing coupled fluid flows in porous media is verified by the resolution of several numerical examples for which previous solutions are known. Finally, the performance of the computational algorithm in modeling of large-scale porous media problems including the large elasto-plastic deformations is demonstrated through the fully coupled analysis of the failure of two earth and rockfill dams. Furthermore, the three-phase model is compared to its simplified one which simulates the unsaturated porous medium as a two-phase one with static air phase. The paper illustrates the shortcomings of the commonly used simplified approach in the context of seismic analysis of two earth and rockfill dams. It is shown that accounting the pore air as an independent phase significantly influences the unsaturated soil behavior.     

Method of negative saturations for flow with variable number of phases in porous media: extension to three-phase multi-component case

Various EOR methods lead to the appearance of specific macroscopic surfaces called interfaces of phase transition (IPT) such that the number of phases on two sides of an IPT is different, and fluids separated by an IPT are in non-equilibrium. Therefore, the flow equations are also different on two sides of an IPT and cannot be deduced from each other by a continuous degeneration, which imposes difficulties in numerical modelling. To describe such systems, we developed a new conceptual mathematical method based on the replacement of the real single-phase fluid by an imaginary multi-phase multi-component continuum having fictitious properties. As the result, the fluid over all zones becomes multi-phase and can be described by uniform multi-phase hydro- and thermodynamic equations, which allows applying the direct numerical simulation. The equivalence principle determines the physical properties of the fictitious multi-phase fluid, as well as the structure of the uniform multi-phase equations. It also proves that the saturation of each phase becomes an extended function negative or higher than unity in non-equilibrium zones, which becomes the efficient method of tracking the interfaces, the number of phases at any point, and their degree of disequilibrium. The method was developed in [1, 2] for the two-phase case. In the present paper, the new version of the method is developed for the three-phase case with gravity, diffusion, and capillarity. We have obtained the new equivalent uniform multi-phase equations which contain additional non-classical terms responsible for the diffusion and gravity across an IPT. The comparison with classical method is presented. The presentation is illustrated by several examples of simulation by means of the code developed by the research group; their concern: EOR by miscible methods and CO ₂ bubble raising in aquifer.     

The influence of capillarity on multiphase flow within porous media: a new model for interpreting fluid levels in groundwater monitoring wells in dynamic aquifers

A new approach is presented to calculate the volume of oil in the underground at an oil spill site from fluid levels in monitoring wells. The approach includes the effects of hysteresis due to irregular pore geometry and to phase entrapment. It is possible to explain the drastic changes in the oil thickness in a monitoring well due to the decrease and increase in the groundwater table. A correct evaluation of the oil volume infiltrated underground from an oil spill and the effective control of remediation works can only be done by using the newly developed approach with a consideration of the dynamic changes in the groundwater table.