Virial EoS

All virial models use the same argument list:

T : float
absolute temperature [K]
Tij: array
square matrix of critical temperatures [K]
Pij: array
square matrix of critical pressures [bar]
wij: array
square matrix of acentric factors
Tsonopoulos(T, Tij, Pij, wij)[source]

Returns array of virial coefficient for a mixture at given temperature with Tsonopoulos correlation for the first virial coefficient, B:

\[\frac{BP_c}{RT_c} = B^{(0)} + \omega B^{(1)}\]

Where \(B^{(0)}\) and \(B^{(1)}\) are obtained from:

\[\begin{split}B^{(0)} &= 0.1445 - \frac{0.33}{T_r} - \frac{0.1385}{T_r^2} - \frac{0.0121}{T_r^3} - \frac{0.000607}{T_r^8} \\ B^{(1)} &= 0.0637 + \frac{0.331}{T_r^2} - \frac{0.423}{T_r^3} - \frac{0.008}{T_r^8}\end{split}\]
Abbott(T, Tij, Pij, wij)[source]

Returns array of virial coefficients for a mixture at given temperature with Abbott-Van Ness correlation for the first virial coefficient, B:

\[\frac{BP_c}{RT_c} = B^{(0)} + \omega B^{(1)}\]

Where \(B^{(0)}\) and \(B^{(1)}\) are obtained from:

\[\begin{split}B^{(0)} &= 0.083 - \frac{0.422}{T_r^{1.6}}\\ B^{(1)} &= 0.139 + \frac{0.179}{T_r^{4.2}}\end{split}\]
ideal_gas(T, Tij, Pij, wij)[source]

Returns array of ideal virial coefficients (zeros). The model equation is

\[Z = \frac{Pv}{RT} = 1\]

Note: Ideal gas model uses only the shape of Tij to produce zeros.