Non-random two-liquid model

Non-random two-liquid model
VLE of the mixture of Chloroform and Methanol plus NRTL fit and extrapolation to different pressures

The non-random two-liquid model[1] (short NRTL equation) is an activity coefficient model that correlates the activity coefficients γi of a compound i with its mole fractions xi in the liquid phase concerned. It is frequently applied in the field of chemical engineering to calculate phase equilibria. The concept of NRTL is based on the hypothesis of Wilson that the local concentration around a molecule is different from the bulk concentration. This difference is due to a difference between the interaction energy of the central molecule with the molecules of its own kind Uii and that with the molecules of the other kind Uij. The energy difference also introduces a non-randomness at the local molecular level. The NRTL model belongs to the so-called local-composition models. Other models of this type are the Wilson model, the UNIQUAC model, and the group contribution model UNIFAC. These local-composition models are not thermodynamically consistent due to the assumption that the local composition around molecule i is independent of the local composition around molecule j. This assumption is not true, as was shown by Flemmer in 1976.[2]


Equations for a binary mixture

For a binary mixture the following equations[3] are used:

\left\{\begin{matrix} \ln\ \gamma_1=x^2_2\left[\tau_{21}\left(\frac{G_{21}}{x_1+x_2 G_{21}}\right)^2 +\frac{\tau_{12} G_{12}} {(x_2+x_1 G_{12})^2 }\right] \\
\\ \ln\ \gamma_2=x^2_1\left[\tau_{12}\left(\frac{G_{12}}{x_2+x_1 G_{12}}\right)^2 +\frac{\tau_{21} G_{21}} {(x_1+x_2 G_{21})^2 }\right]


\left\{\begin{matrix} \ln\ G_{12}=-\alpha_{12}\ \tau_{12}
\\ \ln\ G_{21}=-\alpha_{21}\ \tau_{21}

In here τ12 and τ21 are the dimensionless interaction parameters, which are related to the interaction energy parameters Δg12 and Δg21 by:

\left\{\begin{matrix} \tau_{12}=\frac{\Delta g_{12}}{RT}=\frac{U_{12}-U_{22}}{RT}
\\ \tau_{21}=\frac{\Delta g_{21}}{RT}=\frac{U_{21}-U_{11}}{RT}

Here R is the gas constant and T the absolute temperature, and Uij is the energy between molecular surface i and j. Uii is the energy of evaporation. Here Uij has to be equal to Uji, but Δgij is not necessary equal to Δgji.

The parameters α12 and α21 are the so-called non-randomness parameter, for which usually α12 is set equal to α21. For a liquid, in which the local distribution is random around the center molecule, the parameter α12 = 0. In that case the equations reduce to the one-parameter Margules activity model:

\left\{\begin{matrix} \ln\ \gamma_1=x^2_2\left[\tau_{21} +\tau_{12}  \right]=Ax^2_2
\\ \ln\ \gamma_2=x^2_1\left[\tau_{12}+\tau_{21} \right]=Ax^2_1

In practice, α12 is set to 0.2, 0.3 or 0.48. The latter value is frequently used for aqueous systems. The high value reflects the ordered structure caused by hydrogen bonds. However in the description of liquid-liquid equilibria the non-randomness parameter is set to 0.2 to avoid wrong liquid-liquid description. In some cases a better phase equilibria description is obtained by setting α12 = − 1[4]. However this mathematical solution is impossible from a physical point of view, since no system can be more random than random (α12 =0). In general NRTL offers more flexibility in the description of phase equilibria than other activity models due to the extra non-randomness parameters. However in practice this flexibility is reduced in order to avoid wrong equilibrium description outside the range of regressed data.

The limiting activity coefficients, aka the activity coefficients at infinite dilution, are calculated by:

\left\{\begin{matrix} \ln\ \gamma_1^\infty=\left[\tau_{21} +\tau_{12} exp{(-\alpha_{12}\ \tau_{12})} \right]
\\ \ln\ \gamma_2^\infty=\left[\tau_{12} +\tau_{21}exp{(-\alpha_{12}\ \tau_{21})}\right]

The expressions show that at α12 = 0 the limiting activity coefficients are equal. This situation that occurs for molecules of equal size, but of different polarities.
It also shows, since three parameters are available, that multiple sets of solutions are possible.

General equations

The general equation for ln(γi) for species i in a mixture of n components is[5]:

\ln(\gamma_i)=\frac{\displaystyle\sum_{j=1}^{n}{x_{j}\tau_{ji}G_{ji}}}{\displaystyle\sum_{k=1}^{n}{x_{k}G_{ki}}}+\sum_{j=1}^{n}{\frac{x_{j}G_{ij}}{\displaystyle\sum_{k=1}^{n}{x_{k}G_{kj}}}}{\left ({\tau_{ij}-\frac{\displaystyle\sum_{m=1}^{n}{x_{m}\tau_{mj}G_{mj}}}{\displaystyle\sum_{k=1}^{n}{x_{k}G_{kj}}}}\right )}


G_{ij}=\text{exp}\left ({-\alpha_{ij}\tau_{ij}}\right )


\tau_{ij}=A_{ij}+\frac{B_{ij}}{T}+\frac{C_{ij}}{T^{2}}+D_{ij}\ln{\left ({T}\right )}+E_{ij}T^{F_{ij}}

There are several different equations forms for αij and τij, the most general of which are shown above.

Temperature dependent parameters

To describe phase equilibria over a large temperature regime, i.e. larger than 50 K, the interaction parameter has to be made temperature dependent. Two formats are frequently used. The extended Antoine equation format:

\tau_{ij}=f(T)=a_{ij}+\frac{b_{ij}}{T}+c_{ij}\ \ln\ T+d_{ij}T

Here the logarithmic term is mainly used in the description of liquid-liquid equilibria (miscibility gap).

The other format is a third-order polynomial format:

\Delta g_{ij}=f(T)=a_{ij}+b_{ij}\cdot T +c_{ij}T^{2}

Parameter determination

The NRTL parameters are fitted to activity coefficients that have been derived from experimentally determined phase equilibrium data (vapor–liquid, liquid–liquid, solid–liquid) as well as from heats of mixing. The source of the experimental data are often factual data banks like the Dortmund Data Bank. Other options are direct experimental work and predicted activity coefficients with UNIFAC and similar models. Noteworthy is that for the same liquid mixture several NRTL parameter sets might exist. The NRTL parameter set to use depends on the kind of phase equilibrium (i.e. solid–liquid, liquid–liquid, vapor–liquid). In the case of the description of a vapor–liquid equilibria it is necessary to know which saturated vapor pressure of the pure components was used and whether the gas phase was treated as an ideal or a real gas. Accurate saturated vapor pressure values are important in the determination or the description of an azeotrope. The gas fugacity coefficients are mostly set to unity (ideal gas assumption), but for vapor-liquid equilibria at high pressures (i.e. > 10 bar) an equation of state is needed to calculate the gas fugacity coefficient for a real gas description.


  1. ^ Renon H., Prausnitz J. M., "Local Compositions in Thermodynamic Excess Functions for Liquid Mixtures", AIChE J., 14(1), S.135–144, 1968
  2. ^ McDermott (Fluid Phase Equilibrium 1(1977)33) and Flemr (Coll. Czech. Chem.Comm., 41 (1976) 3347)
  3. ^ Reid R. C., Prausnitz J. M., Poling B. E., The Properties of Gases & Liquids, 4th Edition, McGraw-Hill, 1988
  4. ^ Effective Local Compositions in Phase Equilibrium Correlations, J. M. Marina, D. P. Tassios Ind. Eng. Chem. Process Des. Dev., 1973, 12 (1), pp 67–71
  5. ^

Wikimedia Foundation. 2010.

Look at other dictionaries:

  • Non-Random Two Liquid model — The Non Random Two Liquid model [Renon H., Prausnitz J. M., Local Compositions in Thermodynamic Excess Functions for Liquid Mixtures , AIChE J., 14(1), S.135 144, 1968] (short NRTL equation) is an activity coefficient model that correlates the… …   Wikipedia

  • Random walk — A random walk, sometimes denoted RW, is a mathematical formalization of a trajectory that consists of taking successive random steps. The results of random walk analysis have been applied to computer science, physics, ecology, economics and a… …   Wikipedia

  • liquid — liquidly, adv. liquidness, n. /lik wid/, adj. 1. composed of molecules that move freely among themselves but do not tend to separate like those of gases; neither gaseous nor solid. 2. of, pertaining to, or consisting of liquids: a liquid diet. 3 …   Universalium

  • Ising model — The Ising model, named after the physicist Ernst Ising, is a mathematical model in statistical mechanics. It has since been used to model diverse phenomena in which bits of information, interacting in pairs, produce collectiveeffects.Definition… …   Wikipedia

  • NRTL — may refer to: Non Random Two Liquid model Nationally Recognized Testing Laboratories National Right To Life This disambiguation page lists articles associated with the same title. If an internal link led you he …   Wikipedia

  • gE-Modelle — Unter gE Modellen versteht man Methoden zur Vorhersage von Aktivitätskoeffizienten. Hierbei bedient man sich des Zusammenhangs zwischen freier Exzessenthalpie gE und dem Aktivitätskoeffizienten , der durch folgende Gleichung beschrieben werden… …   Deutsch Wikipedia

  • amorphous solid — ▪ physics Introduction       any noncrystalline solid in which the atoms and molecules are not organized in a definite lattice pattern. Such solids include glass, plastic, and gel.       Solids and liquids (liquid) are both forms of condensed… …   Universalium

  • Distillation — Distiller and Distillery redirect here. For other uses, see Distiller (disambiguation) and Distillery (disambiguation). For other uses, see Distillation (disambiguation). Laboratory display of distillation: 1: A heating device 2: Still pot 3:… …   Wikipedia

  • Chemical biology — is a scientific discipline spanning the fields of chemistry and biology that involves the application of chemical techniques and tools, often compounds produced through synthetic chemistry, to the study and manipulation of biological systems.… …   Wikipedia

  • Stock market — Financial markets Public market Exchange Securities Bond market Fixed income Corporate bond Government bond Municipal bond …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”

We are using cookies for the best presentation of our site. Continuing to use this site, you agree with this.