Equilibrium Diagrams

In order to solve equilibrium stage problems, you must have a model or correlation for the vapor-liquid equilibrium (VLE) and physical properties of the system. When two phases are in equilibrium, VLE data enables us to relate the composition of a liquid phase to the composition of the vapor phase.

Equilibrium data may be obtained by experiment, by thermodynamic calculation, or in published sources. It is typically presented either in tabular form or as an equilbrium diagram. Diagrams may take several forms:

Equilibrium Diagram T vs. x,y

The figure (also see Fig 21.3 in MSH6) shows one common way of plotting equilibrium data -- the Txy diagram. It represents a binary mixture, and all compositions are expressed as mole fractions of the more volatile component; x in the liquid phase or y in the vapor phase. (The less volatile component is thus 1-x and 1- y.)

The lower curve plots the bubble point of the binary mixture as a function of composition. The upper curve is the dew point. For a given temperature and composition, this diagram tells us the nature and composition of each phase of the mixture that is present.

Similar drawings can be constructed at constant temperature while allowing pressure to vary. You should have seen some of these diagrams in your material balance and stoichiometry class (for instance, see Felder & Rousseau (3rd ed.), Section 6.4d).

yx Diagrams

In this class, we will often work with plots of vapor phase composition vs. liquid phase composition, y vs. x (also see Figures 21.2 and 21.20 in MSH6).

Equilibrium Diagram y vs. x
Usually we'll include a "45 degree line" (x=y) on the diagram for reference. These diagrams are typically made at constant pressure, so each point represents a different temperature.

An xy diagram like this may be constructed from a Txy diagram by picking a temperature, reading the corresponding y and x and plotting them against each other.

Tip: When the envelope enclosed by the equilibrium curve and the 45 degree line is "fat", distillation will probably be an easy way to make separations of the mixture.

Ideal VLE

Some equations in this document are being displayed using MINSE, a browser independent approach to displaying equations on the web. If the equations are not properly formatted by your brower, you need to RENDER EQUATIONS (select this link) by invoking Ping's MINSE polymediator. This will run a special program that should cause them to be formatted for viewing by your browser.

Problems can be greatly simplified in cases where the VLE behavior of the system is ideal or can be represented in a simple form. This usually means we can work with one of three main alternatives:

  1. Raoult's Law and Henry's Law
  2. Equilibrium Coefficients ("K-values")
  3. Relative Volatilities

Roult's and Henry's Laws

Raoult's Law (review your notes from material balances and from thermodynamics) relates the vapor and liquid equilibrium compositions of a mixture to the pure component vapor pressure and system pressure:

y*P=x*p;sat

There are few theoretical tools for predicting the equilibria of gases dissolved in liquid solvents. One way of correlating the data is to use Henry's Law, which takes the same form as Raoult's Law but replaces the vapor pressure with a Henry's Law coefficient (also a function of temperature). Henry's Law applies to low solubility gases dissolved in liquids (it would be a bad choice for ammonia, say ...).

Equilibrium "K-values"

A "K-value" or equilibrium coefficient is the ratio of the vapor composition to the liquid composition. These are often correlated as functions of temperature and can be found in published data, as nomographs, or as equations, such as:

K;i='quot(y;i,x;i)=A;i+B;i*T+C;i*'exp(T,2)+'ellipsis
You likely encountered equilibrium coefficients in your multicomponent thermodynamics course. Smith, Van Ness, and Abbott (5th ed., pp. 494-95) have a set of "DePriester Charts", which are plots of correlated K values for hydrocarbons. It might be smart to review using these.

A related concept, the "distribution coefficient", can be used to relate the equilibrium compositions of two liquid phases.

Relative Volatility

The "relative volatility" is the ratio of the K values for two components:

?alpha?;ij='quot(K;i,K;j)='quot('quot(y;i,x;i),'quot(y;j,x;j))
which in the case of a binary system can be rearranged to:
?alpha?;ij='quot(y*(1-x),(1-y)*x)
Better yet, if the relative volatility is constant, the expression can be rearranged into the form needed to plot an equilibrium curve for a set of x-values:
equilibrium form of alpha
(the math can be found at the end of this document).

Relative volatility is generally a much less strong function of temperature than the component vapor pressures; in many systems, it is acceptable to assume that the relative volatility is constant over a range of temperatures and compositions.

High relative volatilities produce xy diagrams with a great deal of separation between the equilibrium curve and the 45 degree line. A relative volatility less than one is probably "upside down" -- the more volatile component is in the denominator.

Making Equilibrium Diagrams

Computer calculation aids immensely in construction of equilibrium diagrams. I've prepared examples of two approaches on MATHCAD. You may down load them by following the links:

NonIdeal Systems

In cases where the VLE cannot be treated as ideal, more rigorous models based on multicomponent thermodynamics are required. We won't worry about these for awhile.

Azeotropes

Most of the time, we deal with systems where the equilibrium curves have no inflection points; that is, as the concentration of the less volatile component increases, so do the dew point and bubble point.

If, however, there are strong physical or chemical interactions between components, diagrams may look different.

In some such systems, there is a critical composition where the liquid and vapor compositions are identical. Once this composition is reached, separation cannot continue without changing pressure. These mixtures are called azeotropes -- they will have minimum or maximum boiling points. Minimum boilers are more common. Composition of the vapor produced from an azeotrope is the same as the liquid, so an azeotrope can be boiled at constant pressure without changing composition.

Enthalpy Concentration Diagrams

For a binary mixture, cp, ?Delta?H, ?lambda?, etc. are functions of composition. All of this information can be used to construct an enthalpy concentration diagram (see MSH6 Fig 21.24 for another example):

Equilibrium Diagram H vs. x,y
These are constructed at constant pressure and each point represents a different temperature (but beware! T(x=0.5) is not equal to T(y=0.5)). The upper curve is a plot of saturated vapor enthalpy vs. vapor composition, the lower curve saturated liquid enthalpy vs. liquid composition; consequently, for any given composition, the difference between the two curves is the latent heat of vaporization of the mixture. Points between the two curves ("inside the envelope") represent two phase mixtures. Points outside the envelope represent superheated vapor or subcooled liquid. Constant temperature tie-lines between the vapor and liquid curves are sometimes provided (although rarely as many as one might want).

Developing Hxy diagrams isn't particularly difficult, but it does require attention to detail. To begin with, reference temperatures must be selected and applied with consistency. It is also important to remember that liquid solution enthalpies must include both sensible heat and heats of mixing:

L = c;p*(T;sys - T;ref) +?Delta?H;soln
It may be most convenient to set the reference temperature at that corresponding to the heat of solution data (which will be negative if heat is evolved on mixing).

Remember that enthalpy is a state function, so you may calculate it along any "path" that you choose; and usually want to choose the one that makes calculation easiest. Vapor enthalpies are typically calculated using heat capacities and the latent heat of vaporization:

H;V = c;pL*(T;bp - T;ref) + ?lambda?(T;bp) + c;pV*(T;sys - T;bp)

Mixture values are usually hard to find, so as an approximation, assume the unmixed liquids are heated separately to the dew point, vaporized at that temperature, and then mixed:

H;V = y;1*(c;p1*(T;dp - T;ref) + ?lambda?;1) + y;2*(c;p2*(T;dp - T;ref) + ?lambda?;2)

There is a direct relationship between Hxy and xy diagrams. If you have tie lines connecting equilibrium liquid and vapor compositions on an Hx diagram, you can use it to construct an xy diagram. The endpoints of the tie lines are at the equilibrium compositions needed/used on the xy diagram.

Constructing an equilibrium curve from an Hxy diagram (or adding tie- lines to an Hxy diagram from data on an Hxy diagram) is a nice construction. Place the equilibrium diagram below the Hxy and line up the scales. Starting from a tie line liquid endpoint on the Hxy, drop a vertical line onto the equilibrium plot. This will be an x-coordinate. Next drop a line from a vertical endpoint on the Hxy. When you hit the 45 degree line, turn at a right angle. The intersection between this line and the previous line will give you one (x,y) point on the equilibrium curve.

You can reverse the construction and use an xy diagram to get the tie lines on your Hx plot. Take x coordinates straight up and mark the intersection with the bubble point curve. Take y-coordinates over to the 45-degree line, turn and carry them up and mark the intersection with the dew point curve. Connect the two intersections to get a tie line.


References:

  1. Kister, H.Z., Distillation Design, McGraw Hill, p. 14.
  2. McCabe, W.L. and J.C. Smith, Unit Operations of Chemical Engineering (3rd Edition), McGraw-Hill, 1976, pp. 486-488.
  3. McCabe, W.L., J.C. Smith, and P. Harriott, Unit Operations of Chemical Engineering (5th Edition), McGraw-Hill, 1993, pp. 508, 553-554.
  4. McCabe, W.L., J.C. Smith, and P. Harriott, Unit Operations of Chemical Engineering (6th Edition), McGraw-Hill, 2001, pp. 628, 672-674.
  5. Seader, J.D. and E.J. Henley, Separation Process Principles, John Wiley, 1998, pp. 36-42, 173-76, 207-11.
  6. Treybal, R.E., Mass-Transfer Operations, 3rd Edition (Reissue), McGraw-Hill, 1987, pp. 343-450,357-360.


Rearrangement of Relative Volatility Expression

Rearrangement of Alpha
Function Form

R.M. Price
Original: 3/6/97, 1/5/98; 8/13/2002
Modified: 1/8/98, 1/8/99; 2/10/2003

Copyright 1998, 1999, 2002, 2003 by R.M. Price -- All Rights Reserved