We found a match
Your institution may have rights to this item. Sign in to continue.
- Title
Numerical modeling of two-phase binary fluid mixing using mixed finite elements.
- Authors
Sun, Shuyu; Firoozabadi, Abbas; Kou, Jisheng
- Abstract
Diffusion coefficients of dense gases in liquids can be measured by considering two-phase binary nonequilibrium fluid mixing in a closed cell with a fixed volume. This process is based on convection and diffusion in each phase. Numerical simulation of the mixing often requires accurate algorithms. In this paper, we design two efficient numerical methods for simulating the mixing of two-phase binary fluids in one-dimensional, highly permeable media. Mathematical model for isothermal compositional two-phase flow in porous media is established based on Darcy's law, material balance, local thermodynamic equilibrium for the phases, and diffusion across the phases. The time-lag and operator-splitting techniques are used to decompose each convection-diffusion equation into two steps: diffusion step and convection step. The Mixed finite element (MFE) method is used for diffusion equation because it can achieve a high-order and stable approximation of both the scalar variable and the diffusive fluxes across grid-cell interfaces. We employ the characteristic finite element method with moving mesh to track the liquid-gas interface. Based on the above schemes, we propose two methods: single-domain and two-domain methods. The main difference between two methods is that the two-domain method utilizes the assumption of sharp interface between two fluid phases, while the single-domain method allows fractional saturation level. Two-domain method treats the gas domain and the liquid domain separately. Because liquid-gas interface moves with time, the two-domain method needs work with a moving mesh. On the other hand, the single-domain method allows the use of a fixed mesh. We derive the formulas to compute the diffusive flux for MFE in both methods. The single-domain method is extended to multiple dimensions. Numerical results indicate that both methods can accurately describe the evolution of the pressure and liquid level.
- Subjects
DIFFUSION; GAS-liquid interfaces; FLUID dynamics; THERMODYNAMICS; FINITE element method; COMPUTER simulation; GRID computing
- Publication
Computational Geosciences, 2012, Vol 16, Issue 4, p1101
- ISSN
1420-0597
- Publication type
Article
- DOI
10.1007/s10596-012-9306-2