Predictive Soave-Redlich-Kwong (PSRK)

Generic cubic equation of state (explicitly arranged for pressure):

<math> P = { { R \; T } \over { V - b } } - { { a(T) } \over { ( V + \epsilon * b ) * ( V + \sigma * b ) } } </math>

Parameters for Soave-Redlich-Kwong equation are:

<math> \epsilon = 0</math> and <math> \sigma = 1</math>

Thus:

<math> P = { { R \; T } \over { V - b } } - { { a(T) } \over { V * ( V + b ) } } </math>

PSRK mixing rule for calculating a(T) and b:

Cohesion pressure (attractive parameter):

<math> a(T) = b*R*T ( \sum x_i { {a_{ii}(T)} \over {b_i*R*T} } + { { { {g_0^E} \over {R*T} } + \sum x_i ln( b / b_i ) } \over { ln( u / (u + 1) )} } )</math>

Mathias-Copeman equation:

Gibbs-Excess energy:

modified UNIFAC:

Procedure for calculating vapour-liquid equilibria (VLE):