Cubic Equations of State¶
Contents
- Cubic Equations of State
- Introduction
- General Cubic Equation of State
- Property Package Options
- Required Parameters
- Calculation of Properties
- Mass Density by Phase
- Molar Density by Phase
- Molar Enthalpy by Phase
- Component Molar Enthalpy by Phase
- Molar Entropy by Phase
- Component Molar Entropy by Phase
- Component Fugacity by Phase
- Component Fugacity Coefficient by Phase
- Molar Gibbs Energy by Phase
- Component Gibbs Energy by Phase
Introduction¶
This module implements a general form of a cubic equation of state which can be used for most cubic-type equations of state. The following forms are currently supported:
- Peng-Robinson
- Soave-Redlich-Kwong
General Cubic Equation of State¶
All equations come from “The Properties of Gases and Liquids, 4th Edition” by Reid, Prausnitz and Poling. The general cubic equation of state is represented by the following equations:
where \(Z\) is the compressibility factor of the mixture, \(a_m\) and \(b_m\) are properties of the mixture and \(u\) and \(w\) are parameters which depend on the specific equation of state being used as show in the table below.
Equation | \(u\) | \(w\) | \(Omega_A\) | \(Omega_B\) | \(kappa_j\) |
---|---|---|---|---|---|
Peng-Robinson | 2 | -1 | 0.45724 | 0.07780 | \((1+(1-T_r^2)(0.37464+1.54226\omega_j-0.26992\omega_j^2))^2\) |
Soave-Redlich-Kwong | 1 | 0 | 0.42748 | 0.08664 | \((1+(1-T_r^2)(0.48+1.574\omega_j-0.176\omega_j^2))^2\) |
The properties \(a_m\) and \(b_m\) are calculated from component specific properties \(a_j\) and \(b_j\) as shown below:
where \(P_{c,j}\) and \(T_{c,j}\) are the component critical pressures and temperatures, \(y_j\) is the mole fraction of component :math`j`, \(k_{ij}\) are a set of binary interaction parameters which are specific to the equation of state and \(\Omega_A\), \(\Omega_B\) and \(\kappa_j\) are taken from the table above. \(\omega_j\) is the Pitzer acentric factor of each component.
The cubic equation of state is solved for each phase via a call to an external function which automatically identifies the correct root of the cubic and returns the value of \(Z\) as a function of \(A\) and \(B\) along with the first and second partial derivatives.
Property Package Options¶
When using the general cubic equation of state module, users must specify the type of cubic to use. This is done by providing a type option in the equation_of_state_options argument in the Phase definition, as shown in the example below.
from idaes.generic_models.properties.core.eos.ceos import Cubic, CubicType
configuration = {
"phases": {
"Liquid": {
"type": LiquidPhase,
"equation_of_state": Cubic,
"equation_of_state_options": {
"type": CubicType.PR}}}
Required Parameters¶
Cubic equations of state require the following parameters to be defined:
- omega (Pitzer acentricity factor) needs to be defined for each component (in the parameter_data for each component).
- kappa (binary interaction parameters) needs to be defined for each component pair in the system. This parameter needs to be defined in the general parameter_data argument for the overall property package (as it can be used in multiple phases).
Calculation of Properties¶
Many thermophysical properties are calculated using an ideal and residual term, such that:
The residual term is derived from the partial derivatives of the cubic equation of state, whilst the ideal term is determined using pure component properties for the ideal gas phase defined for each component.
Mass Density by Phase¶
The following equation is used for both liquid and vapor phases, where \(p\) indicates a given phase:
where \(MW_p\) is the mixture molecular weight of phase \(p\).
Molar Density by Phase¶
Molar density is calculated using the following equation
Molar Enthalpy by Phase¶
The residual enthalpy term is given by:
The ideal component is calculated from the weighted sum of the (ideal) component molar enthalpies.
Component Molar Enthalpy by Phase¶
Component molar enthalpies by phase are calculated using the pure component method provided by the users in the property package configuration arguments.
Molar Entropy by Phase¶
The residual entropy term is given by:
The ideal component is calculated from the weighted sum of the (ideal) components molar enthalpies.
Component Molar Entropy by Phase¶
Component molar entropies by phase are calculated using the pure component methods provided by the users in the property package configuration arguments.
Component Fugacity by Phase¶
Fugacity is calculated from the system pressure and fugacity coefficients as follows:
Component Fugacity Coefficient by Phase¶
The fugacity coefficient is calculated from the departure function of the cubic equation of state as shown below:
Molar Gibbs Energy by Phase¶
For both liquid and vapor phases, the molar Gibbs energy is calculated as the weighted sum of the component molar Gibbs energies for the given phase:
where \(x_{p, j}\) is the mole fraction of component \(j\) in the phase \(p\).
Component Gibbs Energy by Phase¶
Component molar Gibbs energies are calculated using the definition of Gibbs energy: