Example 05.17 Henry Constant for Methane in Benzene at 60 °C with the Help

of the Soave-Redlich-Kwong Equation of State


Calculate the Henry constant for methane(1) in benzene(2) at 60 °C with the help of the Soave-Redlich-Kwong equation of state (k12 = 0.08*)

*H. Knapp, R. Döring, L. Oellrich, U. Plöcker, J.M. Prausnitz, Vapor-Liquid Equilibria for Mixtures of Low Boiling Substances, DECHEMA Chemistry Data Series, Vol. VI, Frankfurt 1982.





input data

Critical data and acentric factor:

Binary parameter kSRK:



The Henry constant is defined as the ratio of gas fugacity to its mole fraction in the liquid in the limiting case of a gas partial pressure of zero. As in this case the liquid mole fraction of the gas is also zero, determination of Henry coefficients from experimental data has to use data for non-zero gas partial pressures and involves extrapolation to zero pressure.

In case the phase equilibrium is calculated from an equation of state, values for the fugacity coefficients in the two phases can also be found at zero concentration of a component. For the phase equilibrium it holds that:

The fugacity of the gas in the vapor phase is calculated as:

Using the definition of the Henry constant (for the limiting case of zero gas fugacity):

At a gas fugacity f1 = 0 it holds that:

In the binary case, the total pressure P is equal to the saturated vapor pressure of the solvent. In order to calculate the Henry constant at a temperature Temp, the phase equilibrium has to be calculated for the following conditions:

In the first step, the equation of state parameters for and between the pure components are calculated from the critical data and the acentric factor at the target temperature of 60°C. The parameters are identical for both phases and do not need to be updated as both phases have the same fixed composition. When calculating Henry coefficients in solvent mixtures, this simplification is in principle not valid, but at low partial pressures of non-associating solvents in the vapor phase, the vapor phase fugacity coefficients can asumed to be unity and for the calculation of the equilibrium pressure, only the fugacity coefficients in the (fixed concentration) liquid phase are required. The number of gases is irrelevant as their concentration in both phases is assumed to be zero.

As both phases are pure solvent at a gas partial pressure of zero, the mixture parameters for the phases are equal to the solvent pure component parameters:

In the next step, the pure component vapor pressure of the solvent (i.e. the total pressure) has to be determined. For this, the pressure explicit equation of state and functions to solve for the vapor and liquid volume are required:

With the help of the difference in fugacity coefficients the vapor pressure calculation function will later iterate the vapor pressure. The equation for mixtures is used (App. D, Equ. 0.215) here but a more simple pure component fugacity equation would be sufficient.

For the root function in the vapor pressure function a starting value for the pressure is required. The function employs a simple linear representation of the vapor pressure in the log(Pr) vs. 1/Tr diagram through the critical point and a decadic logarithm of the reduced pressure of -(w+1) at 1/Tr:

Verifying the equilibrium condition for this pressure leads to a sufficiently small deviation:

At this pressure the fugacity coefficient of the gas in the liquid is:

and the Henry coefficient of methane in benzene is:

The experimental value is approx. 513 bar [J. Horiuti, Sci. Papers Inst. Phys. Chem. Res. (Japan) 17 (341), 125-256 (1931)]. Although rather old, the experimental value is in good agreement with other data at different temperatures.

It should be noted that the result is not purely predictive, as a binary parameter kij has been used that was probably obtained from VLE data over the whole concentration range. Thus the calculation represents an extrapolation to the behavior at infinite dilution.

Using the functions defined above, the Henry coefficient can now be plotted as function of temperature: