An energy stable evolution method for simulating two-phase equilibria of multi-component fluids at constant moles, volume and temperature

Jisheng Kou, Shuyu Sun, Xiuhua Wang

Research output: Contribution to journalArticlepeer-review

18 Scopus citations

Abstract

In this paper, we propose an energy-stable evolution method for the calculation of the phase equilibria under given volume, temperature, and moles (VT-flash). An evolution model for describing the dynamics of two-phase fluid system is based on Fick’s law of diffusion for multi-component fluids and the Peng-Robinson equation of state. The mobility is obtained from diffusion coefficients by relating the gradient of chemical potential to the gradient of molar density. The evolution equation for moles of each component is derived using the discretization of diffusion equations, while the volume evolution equation is constructed based on the mechanical mechanism and the Peng-Robinson equation of state. It is proven that the proposed evolution system can well model the VT-flash problem, and moreover, it possesses the property of total energy decay. By using the Euler time scheme to discretize this evolution system, we develop an energy stable algorithm with an adaptive choice strategy of time steps, which allows us to calculate the suitable time step size to guarantee the physical properties of moles and volumes, including positivity, maximum limits, and correct definition of the Helmhotz free energy function. The proposed evolution method is also proven to be energy-stable under the proposed time step choice. Numerical examples are tested to demonstrate efficiency and robustness of the proposed method.
Original languageEnglish (US)
Pages (from-to)283-295
Number of pages13
JournalComputational Geosciences
Volume20
Issue number1
DOIs
StatePublished - Feb 25 2016

Fingerprint

Dive into the research topics of 'An energy stable evolution method for simulating two-phase equilibria of multi-component fluids at constant moles, volume and temperature'. Together they form a unique fingerprint.

Cite this