# 10.4: Soave-Redlich-Kwong EOS (1972)

• • Vice Provost Emeritus for Global Program & Professor Emeritus (Petroleum and Natural Gas Engineering) at Pennsylvania State University
• Sourced from John A. Dutton: e-Education Institute

In 1972, Soave proposed an important modification to the RK EOS — or shall we say, a modification to vdW EOS. Between the time of vdW EOS and Redlich-Kwong’s, a new concept for fluid characterization was being discussed. Pitzer had introduced the concept of acentric factor in 1955.

All modifications to the vdW EOS had focused on the temperature dependency of the attractive parameter. Soave expanded this by proposing a two-variable dependency for “a”:

$a=a(T, \omega) \label{10.7}$

It was the first time that “a” was expressed not only as a function of temperature, but also as a function of the shape (sphericity) of the molecules (through w, Pitzer’s acentric factor). As we recall, Pitzer’s acentric factor is a measure of the configuration and sphericity of the molecule. It can also be seen as a measure of the deformity of the molecule.

The Soave-Redlich-Kwong EOS is given by the expression:

$\left(P+\frac{\alpha a}{\bar{v}(\bar{v}+b)}\right)(\bar{v}-b)=R T \label{10.8a}$

Like all cubic equations of state, the SRK EOS is also explicit in pressure. Notice, for example, how the SRK EOS readily becomes:

$P=\frac{R T}{\bar{v}-b}-\frac{\alpha a}{\bar{v}(\bar{v}+b)} \label{10.8b}$

where,

$\alpha=\left[1+\left(0.48508+1.55171 \omega-0.15613 \omega^{2}\left(1-\sqrt{T}_{r}\right)\right]\right) \label{10.8c}$

The influence of acentric factor and temperature on the attractive term is introduced now through “a”. What do we do next? We apply the criticality conditions to Equation \ref{10.8b}. Notice that expression \ref{10.8c} becomes unity at $$Tr=1$$, throughout the critical isotherm. We obtain:

$a=0.427480 \frac{R^{2} T_{c}^{2}}{P_{c}} \label{10.9a}$

$b=0.086640 \frac{R T_{c}}{P_{c}} \label{10.9b}$

Now we show the cubic form (in compressibility factor) of Soave-Redlich-Kwong EOS. Defining,

$A=\frac{\alpha a P}{R^{2} T^{2}} \label{10.10a}$

$B=\frac{b P}{R T} \label{10.10b}$

we are able to obtain:

$z^{3}-Z^{2}+\left(A-B-B^{2}\right) Z-A B=0 \label{10.11}$

For mixtures, Soave proposed a “little” modification to the mixing rules with which we have dealt with so far by introducing the use of “binary interaction parameters” ($$k_{ij}$$):

$(\alpha a)_{m}=\sum \sum y_{i j_{j}(\alpha a)_{i j}} \label{10.12a1}$

$(\alpha a)_{i j}=\sqrt{(\alpha a)_{i}(\alpha a)}_{j}\left(1-k_{i j}\right) \label{10.12a2}$

$b_{m}=\sum_{i} y_{i} b_{i} \label{10.12b}$

The use of binary interaction parameters ($$k_{ij}$$) generated a lot of resistance upon their first introduction. This is because there is no analytical, science-based derivation that justifies their existence. Nowadays, they are regarded just as they are, empirical factors used to tune equations of state and make them match experimental data for mixtures. This has become the heuristic justification for their existence: with them, EOS can do a better job of matching experimental data. Heuristically speaking, they are a measure of interaction between a pair of dislike molecules. Based on this “definition,” their value is zero for pairs of molecules that are alike. Actually, this is no more than a mathematical requirement in order for Equations \ref{10.12a1} and \ref{10.12a2} to give $$(\alpha a)_{i j}=(\alpha a)_{i}$$ when $$j=i$$. The determination of kij is based on experimental data from binary systems; “kij” results from the value that allows the given equation of state (through the expression in 10.12a) to yield the closest match. These values are assumed to be constant (and so are used) when the same two components are part of a more complex multi-component mixture.