PSRK (short for Predictive Soave–Redlich–Kwong) 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 containing supercritical components. This class of substances cannot be predicted with established models, for example UNIFAC.
Principle PSRK is a group-contribution equation of state. This is a class of prediction methods that combines equations of state (mostly cubic) with activity coefficient models based on group contributions, such as UNIFAC. The activity coefficient model is used to adapt the equation-of-state parameters for mixtures by a so-called mixing rule. The use of an equation of state introduces all thermodynamic relations defined for equations of state into the PRSK model. This allows the calculation of densities, enthalpies, heat capacities, and other properties.
Equations As stated previously, the PSRK model is based on a combination of the Soave–Redlich–Kwong equation of state with a mixing rule whose parameters are determined by the UNIFAC method.
Equation of state The equation of state of Soave is defined as follows:
P = R T v − b − a α ( T ) v ( v + b ) . {\displaystyle P={\frac {RT}{v-b}}-{\frac {a\alpha (T)}{v(v+b)}}.}
The original α-function has been replaced by the function of Mathias–Copeman:
α ( T r ) = [ 1 + c 1 ( 1 − T r ) + c 2 ( 1 − T r ) 2 + c 3 ( 1 − T r ) 3 ] 2 . {\displaystyle \alpha (T_{r})=\left[1+c_{1}\left(1-{\sqrt {T_{r}}}\right)+c_{2}\left(1-{\sqrt {T_{r}}}\right)^{2}+c_{3}\left(1-{\sqrt {T_{r}}}\right)^{3}\right]^{2}.}
The parameters of the Mathias–Copeman equation are fitted to experimental vapor-pressure data of pure components and provide a better description of the vapor pressure than the original relation. The form of the equation is chosen as it can be reduced to the original Soave form by setting the parameters c2 and c3 to zero. Additionally, the parameter c1 can be obtained from the acentric factor, using the relation
c 1 = 0.48 + 1.574 ω − 0.176 ω 2 . {\displaystyle c_{1}=0.48+1.574\,\omega -0.176\,\omega ^{2}.}
This may be performed if no fitted Mathias–Copeman parameter is available.
Mixing rule The PSRK mixing rule calculates the parameters a and b of the equation of state by
a b R T = ∑ i x i a i b i R T − g 0 E R T + ∑ x i ln b b i 0.64663 {\displaystyle {\frac {a}{bRT}}=\sum _{i}x_{i}{\frac {a_{i}}{b_{i}RT}}-{\frac {{\frac {g_{0}^{E}}{RT}}+\sum x_{i}\ln {\frac {b}{b_{i}}}}{0.64663}}}
and
… excerpt ends here. Continue reading the full article.

