VTPR (short for Volume-Translated Peng–Robinson) is an estimation method for the calculation of phase equilibria of mixtures of chemical components. The original goal for the development of this method was to enable the estimation of properties of mixtures which contain supercritical components. These class of substances couldn't be predicted with established models like UNIFAC.
Principle VTPR is a group contribution equation of state. This is class of prediction methods combine equations of state (mostly cubic) with activity coefficient models based on group contributions like UNIFAC. The activity coefficient model is used to adapt the equation of state parameters for mixtures by a so-called mixing rule. The usage of an equation of state introduces all thermodynamic relations defined for equations of state into the VTPR model. This allows the calculation of densities, enthalpies, heat capacities, and more.
Equations VTPR is based on a combination of the Peng–Robinson equation of state with a mixing rule whose parameters are determined by UNIFAC.
Equation of state The Peng–Robinson equation of state is defined as follows:
P = R T v − b − a α ( T ) v 2 + 2 b v − b 2 {\displaystyle P={\frac {R\;T}{v-b}}-{\frac {a\;\alpha (T)}{v^{2}+2bv-b^{2}}}}
The originally used α-function has been replaced by the function of Twu, Bluck, Cunningham and Coon .
α ( T r ) = T r N ( M − 1 ) e x p ( L ( 1 − T r M N ) ) {\displaystyle \alpha (T_{r})=T_{r}^{N\left(M-1\right)}exp\left(L\left(1-T_{r}^{MN}\right)\right)}
The parameters of the Twu equation are fitted to experimental vapor pressure data of pure components and guarantee therefore a better description of the vapor pressure than the original relation.
Mixing rule The VTPR mixing rule calculate the parameter a and b of the equation of state by
a ( T ) = b ⋅ ( ∑ i x i a i i ( T ) b i i + g r e s E − 0.53087 ) {\displaystyle a(T)=b\cdot \left(\sum _{i}{x_{i}}{\frac {a_{ii}(T)}{b_{ii}}}+{\frac {g_{res}^{E}}{-0.53087}}\right)}
with
P r e f = 1 a t m {\displaystyle P_{ref}=1\,atm}
and
b i j 3 / 4 = b i i 3 / 4 + b j j 3 / 4 2 {\displaystyle b_{ij}^{3/4}={\frac {b_{ii}^{3/4}+b_{jj}^{3/4}}{2}}}
b i i = 0.0778 ⋅ R ⋅ T c , i P c , i {\displaystyle b_{ii}=0.0778\cdot {\frac {R\cdot T_{c,i}}{P_{c,i}}}}
… excerpt ends here. Continue reading the full article.

