Computation of Phase Equilibrium and Phase Envelopes

Research output: Book/ReportReport – Annual report year: 2017Researchpeer-review


View graph of relations

In this technical report, we describe the computation of phase equilibrium and phase envelopes based on expressions for the fugacity coefficients. We derive those expressions from the residual Gibbs energy. We consider 1) ideal gases and liquids modeled with correlations from the DIPPR database and 2) nonideal gases and liquids modeled with cubic equations of state. Next, we derive the equilibrium conditions for an isothermal-isobaric (constant temperature, constant pressure) vapor-liquid equilibrium process (PT flash), and we present a method for the computation of phase envelopes. We formulate the involved equations in terms of the fugacity coefficients. We present expressions for the first-order derivatives. Such derivatives are necessary in computationally efficient gradient-based methods for solving the vapor-liquid equilibrium equations and for computing phase envelopes. Finally, we describe a Matlab program that computes the phase envelope of a mixture. We present the source code and discuss practical details of the implementation.
Original languageEnglish
PublisherDTU Compute
Number of pages52
Publication statusPublished - 2017
SeriesDTU Compute-Technical Report-2017
Download as:
Download as PDF
Select render style:
Download as HTML
Select render style:
Download as Word
Select render style:

Download statistics

No data available

ID: 140838697