# 16.2: van der Waals and Redlich-Kwong Equations of State

- Page ID
- 13657

## The Van der Waals Equation of State

The van der Waals Equation of State is an equation relating the density of gases and liquids to the pressure, volume, and temperature conditions (i.e., it is a thermodynamic equation of state). It can be viewed as an adjustment to the ideal gas law that takes into account the non-zero volume of gas molecules, which are subject to an inter-particle attraction. It was derived in 1873 by Johannes Diderik van der Waals, who received the Nobel Prize in 1910 for this work.

Previously, we saw that for the pair potential

\[u(r) = \begin{cases} \infty & r \leq \sigma \\ -\dfrac{C_6}{r^6} & r > \sigma \end{cases} \label{1}\]

we could write the second virial coefficient as

\[B_2(T) = \dfrac{2}{3} \pi N_0 \sigma^3 \left[ 1 - \dfrac{C_6}{3 k_B T \sigma^6} \right] \label{2}\]

Let us introduce to simplifying variables

\[\begin{align} b &= \dfrac{2}{3} \pi N_0 \sigma^3 \\ a &= \dfrac{2 \pi N_0^2 C_6}{9 \sigma^3} \end{align} \label{3}\]

in terms of which

\[B_2(T) = b - \dfrac{a}{RT} \label{4}\]

With these definitions, the virial equation of state becomes

\[\begin{align} P &= \dfrac{nRT}{V} + \dfrac{n^2}{V^2} RT \left( b - \dfrac{a}{RT} \right) \\ &= \dfrac{nRT}{V} \left(1 + \dfrac{nb}{V} \right) - \dfrac{an^2}{V^2} \end{align} \label{5)}\]

If we assume \(nb/V\) is small, then we can also write

\[1 + \dfrac{nb}{V} \approx \dfrac{1}{1 - \dfrac{nb}{V}} \label{6)}\]

so that

\[P = \dfrac{nRT}{V - nb} - \dfrac{an^2}{V^2} \label{Eq7}\]

which is known as the *van der Waals* equation of state. Equation \(\ref{Eq7}\) can also be rewritten as

\[\left(P+{\dfrac {n^{2}a}{V^{2}}}\right)\left(V-nb\right)=nRT\]

The first term in this equation is easy to motivate. In fact, it looks very much like the equation of state for an ideal gas

\[PV=nRT\]

however, using \(V - nb\) fort the volume rather than \(V\). This part of the van der Waals equation is due entirely to the hard-wall potential \(u_o(r)\). Essentially, this potential energy term describes a system of “billiard balls” of diameter \(\sigma\). Figure \(\PageIndex{1}\) shows two of these billiard ball type particles at the point of contact (i.e, the distance of closest approach). At this point, they undergo a collision and separate, so they cannot be closer than that distance.

At this point, the distance between their centers is also \(\sigma\), as Figure \(\PageIndex{1}\) indicates. Because of this distance of closest approach, the total volume available to the particles is not \(V\) but some volume less than \(V\). This reduction in volume can be calculated as follows: Figure \(\PageIndex{1}\) shows a shaded sphere that just contains the pair of billiard ball particles. The volume of this sphere is the volume excluded from any two particles. The radius of the sphere is \(\sigma\) as the figure shows. Hence, the excluded volume for the two particles is \(4 \pi \sigma^3/3\), which is the volume of the shaded sphere. From this, we see that the excluded volume for any *one* particle is just half of this or \(2 \pi \sigma^3/3\). The excluded volume for a mole of such particles is just \(2 \pi \sigma^3 N_0/3\), which is the parameter \(b\):

\[b = \dfrac{2}{3} \pi \sigma^3 N_0 \label{\(\PageIndex{8}\)}\]

Given \(n\) moles of gas, the total excluded volume is then \(nb\), so that the total *available* volume is simply \(V - nb\).

## Redlich-Kwong Equation of State

The van der Waals cubic Equation of State had to wait almost 100 years before a real, successful improvement was introduced to it. This progress occurred once researchers committed themselves to finding the empirical temperature dependency of the attraction parameter \(a\) proposed by van der Waals. In contrast, very little attention has been paid to modifying the parameter \(b\) for co-volume. It makes a lot of sense that \(b\) would not be modified by temperature, because it represents the *volume* of the molecules, which should not be affected by their kinetic energy (measured in terms of temperature). The very first noteworthy successful modification to the attraction parameter came with the publication of the equation of state of Redlich-Kwong in 1949.

The Redlich–Kwong equation of state is an empirical, algebraic equation that relates temperature, pressure, and volume of gases. It is generally more accurate than the van der Waals equation and the ideal gas equation at temperatures above the critical temperature. It was formulated by Otto Redlich and Joseph Neng Shun Kwong in 1949, who showed that a simple two-parameter, cubic equation of state could well reflect reality in many situations. Redlich and Kwong revised the van der Waals Equation of State (Equation \(\ref{Eq7}\)) and proposed the following expressions:

\[ \left( P + \dfrac{a}{\sqrt{T} \bar{V} (\bar{V} + b)} \right) ( \bar{V}-b) = RT \label{10.1}\]

Notice that the fundamental change they introduced was to the functional form of \(δP_{attraction}\) (Equation 7.8). Additionally, they introduced the co-volume “b” into the denominator of this functional form. The important concept here is that the attraction parameter \(a\) of van der Waals needed to be made a function of temperature to do a better job of quantitatively matching experimental data. This was a realization that van der Waal had suggested, but no actual functional dependency had been introduced until the Redlich-Kwong equation.

We know what follows at this point. To come up with an expression for \(a\) and \(b\) of Equation \(\ref{10.1}\), we apply the criticality conditions to this EOS. As we recall, imposing the criticality conditions allows us to relate the coefficients “a” and “b” to the critical properties (P_{c}, T_{c}) of the substance. Once we have done that, we obtain the definition of “a” and “b” for the Redlich-Kwong EOS,

\[ a =0.42780 \dfrac{R^2T_c^{2.5}}{P_c} \label{10.2a}\]

\[ b =0.086640 \dfrac{RT_c}{P_c} \label{10.2b}\]

This EOS radically improved, in a quantitative sense, the predictions of vdW EOS. We now recall that vdW-type equations are cubic because they are cubic polynomials in molar volume and compressibility factor. It comes as no surprise then, that we can transform Equation \(\ref{10.1}\) into:

\[ \bar{v} ^3 - \left( \dfrac{RT}{P} \right)\bar{v}^2 + \dfrac{1}{P} \left( \dfrac{1}{T^{0.5}} - bRT - Pb^2 \right) \bar{v} - \dfrac{ab}{PT^{0.5}} = 0 \label{10.3}\]

and, by defining the following parameters,

\[ A = \dfrac{aP}{R^2T^{2.5}} \label{10.3a}\]

\[ B =\dfrac{bP}{RT} \label{10.3b}\]

and introducing the compressibility factor definition

\[Z = \dfrac{P\bar{v}}{RT}\]

we get:

\[ Z^3 -Z^2 + (A -B -B^2) Z -AB =0 \label{10.4}\]

We may also verify the two-parameter corresponding state theory by introducing Equations \(\ref{10.2a}\), \(\ref{10.2b}\), and \(\ref{10.3}\) into Equation \(\ref{10.4}\)),

\[Z^3−Z^2+\dfrac{P_r}{T_r} \left( \dfrac{0.42748}{T_r^{1.5}} −0.08664−0.007506 \dfrac{P_r}{T_r} \right)Z−0.03704 \dfrac{P^2_r}{T^{3.5}_r}=0 \label{10.5}\]

In Equation \(\ref{10.5}\), we can observe the same thing that we saw with vdW EOS: gases at corresponding states have the same properties. Equation \(\ref{10.5}\) is particularly clear about it: any two different gases at the same P_{r}, T_{r} condition __have the same compressibility factor.__

Just as any other cubic equation of state, Equations \(\ref{10.1}\) through \(\ref{10.5}\), as they stand, are to be applied to pure substances. For mixtures, however, we apply *the same equation,* __but__ we impose certain mixing rules to obtain “a” and “b”, which are functions of the properties of the pure components. Strictly speaking, we create a new “pseudo” pure substance that has the average properties of the mixture. Redlich-Kwong preserved the same mixing rules that vdW proposed for his EOS:

\[ a_m = \sum_i \sum_j y_iy_j a_{ij}\]

with

\[a_{ij} = \sqrt{a_ia_j}\]

and

\[ b_{m}=\sum_i y_ib_i\]

Naturally, Redlich and Kwong did not have the last word on possible improvements to the vdW Equation fo State. The Redlich-Kwong Equation fo State, as shown here, is no longer used in practical applications. Research continued and brought with it new attempts to improve the RK Equation fo State. After more than two decades, a modified RK EOS with very good potential was developed.

## Contributors and Attributions

Michael Adewumi (The Pennsylvania State University) Vice Provost for Global Program, Professor of Petroleum and Natural Gas Engineering