Che 312 Notebook Distillation
Che 312 Notebook Distillation
Che 312 Notebook Distillation
CHE 312
SEPARATION PROCESSES 1
DISTILLATION
Distillation can be defined as the process of separating, concentrating or purifying liquid by boiling it and then
condensing the resulting vapour. It is the most widely used method of separating liquid mixtures into their vapour
components. It is one of the major operations in the process industries and the key operation in any oil refinery.
Distillation, in the industry, takes place in a distillation column. A distillation column is equipment for separating
liquid mixtures. It is a hollow vertical column, fitted with trays inside or packing material, in which liquid
mixtures are separated into their components by boiling the mixture and condensing the resulting vapour. The
vertical cylindrical column provides, in a compact form and with the minimum of ground requirements, a large
number of separate stages of vaporization and condensation. In engineering terms, distillation columns have to be
designed with a larger range in capacity than any other types of processing equipment, with single columns 0.3–
10 m in diameter and 3–75 m in height. Designers are required to achieve the desired product quality at minimum
cost and also to provide constant purity of product even though there may be variations in feed composition.
The separation of benzene from a mixture with toluene, for example, requires only a simple single unit as shown
in Figure 1 and virtually pure products may be obtained. A more complex arrangement is shown in Figure 2
where the columns for the purification of crude styrene formed by the dehydrogenation of ethyl benzene. In this
case, several columns are required and it is necessary to recycle some of the streams to the reactor.
VAPOUR–LIQUID EQUILIBRIUM
The composition of the vapour in equilibrium with a liquid of given composition is determined experimentally
using an equilibrium still. The results are conveniently shown on a temperature–composition diagram as shown in
Figure 3.In the normal case shown in Figure 3a, the curve ABC shows the composition of the liquid which boils
at any given temperature and the curve ADC the corresponding composition of the vapour at that temperature.
Thus, a liquid of composition x1 will boil at temperature T1, and the vapour in equilibrium is indicated by point D
of composition y1. It is seen that for any liquid composition x the vapour formed will be richer in the more volatile
component, where x is the mole fraction of the more volatile component in the liquid, and y in the vapour.
Examples of mixtures giving this type of curve are benzene–toluene, n-heptane–toluene, and carbon disulphide–
carbon tetrachloride. In Figures 3b and c, there is a critical composition xgwhere the vapour has the same
composition as the liquid, so that no change occurs on boiling. Such critical mixtures are called azeotropes.
For distillation purposes, it is more convenient to plot y against x at a constant pressure, since the majority of
industrial distillations take place at substantially constant pressure. This is shown in Figure 4 where it should be
noted that the temperature varies along each of the curves.
is obtained as shown by ABCJ in Figure 5. The corresponding dew point curve ADEJ shows the temperature at
which a vapour of composition y starts to condense.
If a mixture of composition x2 is at a temperature T3 below its boiling point, T2, as shown by point G on the
diagram, then on heating at constant pressure the following changes will occur:
When the temperature reaches T2, the liquid will boil, as shown by point B, and some vapour of
composition y2, shown by point E, is formed.
On further heating, the composition of the liquid will change because of the loss of the more volatile
component to the vapour and the boiling point will therefore rise to some temperature T1. At this
temperature, the liquid will have a composition represented by point L, and the vapour a composition
represented by point N. Since no material is lost from the system, there will be a change in the proportion
of liquid to vapour, where the ratio is:
On further heating to a temperature T1, all of the liquid is vaporized to give vapour D of the same
composition y1 as the original liquid. It may be seen that partial vaporization of the liquid gives a vapour
richer in the more volatile component than the liquid. If the vapour initially formed, as for instance at
point E, is at once removed by condensation, then a liquid of composition x3 is obtained, represented by
point C. The step BEC may be regarded as representing an ideal stage, since the liquid passes from
composition x2 to a liquid of composition x3, which represents a greater enrichment in the more volatile
component than can be obtained by any other single stage of vaporization.
Starting with superheated vapour represented by point H, on cooling to D condensation commences, and the first
drop of liquid has a composition K. Further cooling to T1 gives liquid L and vapour N. Thus, partial condensation
brings about enrichment of the vapour in the more volatile component in the same manner as partial vaporization.
The industrial distillation column is, in essence, a series of units in which these two processes of partial
vaporization and partial condensation are effected simultaneously.
4
Since in an ideal gas or vapour the partial pressure is proportional to the mole fraction of the constituent, then:
= ………………………………………………………………………………………………….. (1)
For an ideal mixture, the partial pressure is related to the concentration in the liquid phase by Raoult’s law which
may be written as:
= ………………………………………………………………………………………………… (2)
This relation is usually found to be true only for high values of , or correspondingly low values of , although
mixtures of organic isomers and some hydrocarbons follow the law closely. For low values of , a linear relation
between and again exists, although the proportionality factor is Henry’s constant H, and not the vapour
pressure of the pure material.
For a liquid solute A in a solvent liquid B, Henry’s law takes the form:
= ………………………………………………………………………………………………….. (3)
If the mixture follows Raoult’s law, then the vapour pressure of a mixture may be obtained graphically from
knowledge of the vapour pressure of the two components. Thus, in Figure 6, OA represents the partial pressure
of A in a mixture, and CB the partial pressure of B, with the total pressure being shown by the line BA. In a
mixture of composition D, the partial pressure is given by DE, PB by DF, and the total pressure P by DG,
from the geometry of Figure 6.
Figure 7 shows the partial pressure of one component A plotted against the mole fraction for a mixture that is not
ideal. It is found that over the range OC the mixture follows Henry’s law, and over BA it follows Raoult’s law.
Although most mixtures show wide divergences from ideality, one of the laws is usually followed at very high
and very low concentrations. If the mixture follows Raoult’s law, then the values of for various values of
may be calculated from knowledge of the vapour pressures of the two components at various temperatures.
Thus,
= , = ,
∴ = , and = ………………………………………………………………………….(4)
But + = 1, + =1
∴ + =1
(1 − )
+ =1
+ (1 − )=
− + =
− = −
−
∴ = …………………………………………………………………………………………….(5)
−
Example 1: The vapour pressures of n-heptane and toluene at 373 K are 106 and 73.7 kN/m2 respectively. What
are the mole fractions of n-heptane in the vapour and in the liquid phase at 373 K if the total pressure is 101.3
kN/m2?
106 × 0.856
= = = 0.896
101.3
6
One of the most widely used correlations of saturated vapour pressure is that proposed by Antoine. This takes the
form:
ln = 1 − 2 /( + 3 )………………………………………………………………………………(6)
This is only valid over limited ranges of both temperature and pressures.
Similarly,
= 102.22742 = 168.82 22.51 / 2
2
= = 22.51 × 0.5 = 11.255 /
Example 3:What is the boiling point of an equimolar mixture of benzene and toluene at 101.3 kN/m2?
Solution: The saturation vapour pressures are calculated as a function of temperature using the Antoine equation
and the constants for benzene and toluene. Using Raoult’s law, the actual vapour pressures are given by:
= and = .
2
Then, the temperature at which ( + ) = 101.3 / is determined by a process of trial and error. The
data with pressures in / 2 , are:
T (K) ( + )
373 180.006 74.152 90.003 37.076 127.079
353 100.988 38.815 50.494 19.408 69.902
363 136.087 54.213 68.044 27.106 95.150
7
T (K) ( + )
373.2 180.006 0.2813 74.152 0.6831 0.9644
371.2 170.451 0.2872 69.760 0.7261 1.0233
371.7 172.803 0.2931 70.838 0.7150 1.0081
371.9 173.751 0.29 71.273 0.7107 1.0021
372.0 174.225 0.2907 71.491 0.7085 0.9992
As 0.9992 is near enough to 1.0000, the dew point may be taken as 372.0 K.
Relative Volatility
The relationship between the composition of the vapour yA and of the liquid xA in equilibrium may also be
expressed in a way, which is particularly useful in distillation calculations. If the ratio of the partial pressure to the
mole fraction in the liquid is defined as the volatility, then:
= and = …………………………………………………….(7)
The ratio of these two volatilities is known as the relative volatility αgiven by: = =
But = and =
∴ = …………………………………………………………………………………………….(8)
=∝ ……………………………………………………………………………………………...(9)
This gives a relation between the ratio of A and B in the vapour to that in the liquid.
But + = 1 and + = 1, such that =1− and = 1 − , then:
1−
∝=
1−
∝
∴ = 1+ ∝ −1 ……………………………………………………………………………………(10)
Similarly,
= ∝− ∝ −1 ………………………………………………………………………………………(11)
For separation to be achieved, must not equal 1 and, considering the more volatile component, as increases and
the separation becomes much easier.
From the definition of the volatility of a component, it is seen that for an ideal system, the volatility is
numerically equal to the vapour pressure of the pure component. Thus, the relative volatility may be expressed as:
∝= ………………………………………………………………………………………………..(12)
METHODS OF DISTILLATION
8
For a binary mixture with a normal y-x curve, the vapour is always richer in the more volatile component than the
liquid from which it is formed. There are three main methods used in distillation practice which all rely on this
basic fact. These are:
(a) Differential distillation: In this process, the vapour formed, when a still pot initially full, is heated at constant
rate, on boiling the liquid, is removed at once from the system. At the end of the process, the liquid which has
not been vaporized is removed as the bottom product.
(b) Flash or equilibrium distillation: This is frequently carried out as a continuous process. It consists of
vaporizing a definite fraction of the liquid feed in such a way that the vapour evolved is in equilibrium with
the residual liquid. The feed is usually pumped through a fired heater and enters the still through a valve
where the pressure is reduced. The still is essentially a separation in which the liquid and vapour produced by
the reduction in pressure have sufficient time to reach equilibrium. The vapour is removed from the top of the
separator and is then usually condensed, while the liquid leaves from the bottom.
(c) Rectification: In the two processes considered, the vapour leaving the still at any time is in equilibrium with
the liquid remaining, and normally, there will be only a small increase in concentration of the more volatile
component. The essentially merit of rectification is that it enables a vapour to be obtained that is substantially
richer in the more volatile component (MVC) than is the liquid left in the still. This is achieved by an
arrangement known as a ‘fractionating column’ which enables successive vaporization and condensation to be
established in one unit.
+1 −1 + −1 −1 = + + + ………………….(13)
This equation is difficult to handle for the majority of mixtures, and some simplifying assumptions are
usually made. Thus, with good lagging, the heat losses will be small and may be neglected, and for an
ideal system, the heat of mixing is zero. For such mixtures, the molar heat of vaporization may be taken
as constant and independent of the composition. Thus, one mole of vapour Vn-1 on condensing releases
sufficient heat to liberate one3 mole of vapour Vn. It follows that Vn = Vn-1, so that the molar vapour flow
is constant up the column unless material enters or is withdrawn from the section. The temperature
change from one plate to the next will be small and may be taken as equal to +1 . Thus, = +1 ,
so that the moles of liquid reflux are also constant in this section of the column. Thus, Vn and Ln are
constant over the rectifying section, and Vm and Lm are constant over the stripping section. For these
conditions, there are two basic methods for determining the number of plates required.
= +1 +1 +
11
+1
= +1 + ……………………………………………………………………………(15)
This equation relates the composition of the vapour rising to the plate to the composition of the liquid on
any plate above the feed plate. Since the molar liquid overflow is constant, may be replaced by +1
and:
= +1 + ………………………………………………………………………………(16)
Similarly, taking a material balance for the total streams and for the MVC from the bottom to above plate
m, as indicated by the loop II in Figure 10, and noting that = +1 gives:
+1 = + ……………………………………………………………………………………..(17)
= +1 −
∴ = +1 − …………………………………………………………………………(18)
This equation, which is similar to equation (16), gives the corresponding relation between the
composition of the vapour rising to a plate and the liquid on the plate, for the section below the feed plate.
These two equations are the equations of the operating lines.
Example 4: A mixture of benzene and toluene containing 40 mole percent benzene is to be separated to
give a product containing 90 mole percent benzene at the top, and a bottom product containing not more
than 10 mole percent benzene. The feed enters the column at its boiling point, and the vapour leaving the
column which is condensed but not cooled, provides reflux and product. It is proposed to operate the unit
with a reflux ratio of 3 kmol/kmol of product. It is required to find the number of theoretical plates
needed and the position of entry for the feed. Pressure is at 100 kN/m2.
Since the feed is all liquid at its boiling point, this will all run down as increased reflux to the plate below.
Thus: = + = 112.5 + 100 = 212.5
Also, = − = 212.5 − 62.5 = 150
∴ = +1 −
212.5 62.5
+1 − = × 0.1
150 150
= 1.415 +1 − 0.042……………………………………………………………………....(iv)
Mole fraction of benzene in liquid 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
Mole fraction of benzene in vapour 0 0.22 0.38 0.51 0.63 0.70 0.78 0.85 0.91 0.96
With the two equations (iii) and (iv) and the equilibrium curve, the composition on the various
plates may be calculated by working either from the still up to the condenser, or in the reverse
direction. Since all the vapour from the column is condensed, the composition of the vapour
from the top plate must equal that of the product , and that of the liquid returned as reflux .
The composition of the liquid on the top plate is found from the equilibrium curve and, since
it is in equilibrium with vapour of composition, = 0.90, = 0.79.
13
Figure 11: Calculation of the Number of Plates by the Lewis–Sorel Method for Example 4.
−3 = (0.75 × 0.492) + 0.225 = (0.369 + 0.225) = 0.594, −3 from the equilibrium curve = 0.382
This last value of composition is sufficiently near to that of the feed for the feed to be introduced on plate
(t − 3). For the lower part of the column, the operating line equation (iv) will be used.
Thus: −4 = (1.415 × 0.382) − 0.042 = (0.540 − 0.042) = 0.498, −4 from the equilibrium curve = 0.298
−5 = (1.415 × 0.298) − 0.042 = (0.421 − 0.042) = 0.379, −5 from the equilibrium curve = 0.208
−6 = (1.415 × 0.208) − 0.042 = (0.294 − 0.042) = 0.252, −6 from the equilibrium curve = 0.120
−7 = (1.415 × 0.120) − 0.042 = (0.169 − 0.042) = 0.127, −7 from the equilibrium curve = 0.048
This liquid −7 is slightly weaker than the minimum required and it may be withdrawn as the bottom product.
Thus, −7 will correspond to the reboiler, and there will be seven plates in the column.
= +1 − ………………………………………………………………………….(18)
These equations are used in the Lewis–Sorel method to calculate the relation between the composition of the
liquid on a plate and the composition of the vapour rising to that plate. McCabe and Thiele pointed out that, since
these equations represent straight lines connecting with +1 and with +1 , they can be drawn on the
same diagram as the equilibrium curve to give a simple graphical solution for the number of stages required. Thus,
14
the line of equation (16) will pass through the points 2, 4 and 6 shown in Figure 11, and similarly the line of
equation (18) will pass through points 8, 10, 12 and 14.
and this equation represents a line passing through the point = +1 = . If +1 is put equal to zero, then,
= / , giving a second easily determined point. The top operating line is therefore drawn through two
points of coordinates ( , ) and (0, ).
For the bottom operating line, equation (18), if +1 = , then:
= − ……………………………………………………………………………………(20)
Since = − , it follows that = . Thus the bottom operating line passes through the point C, that
is ( ), and has a slope . When the two operating lines have been drawn in, the number of stages
required may be found by drawing steps between the operating line and the equilibrium curve starting from point
A. This method is one of the most important concepts in chemical engineering and is an invaluable tool for the
solution of distillation problems. The assumption of constant molar overflow is not limiting since in very few
systems do the molar heats of vaporization differ by more than 10 per cent. The method does have limitations,
however, and should not be employed when the relative volatility is less than 1.3 or greater than 5, when the
reflux ratio is less than 1.1 times the minimum, or when more than twenty-five theoretical trays are required. In
these circumstances, the Ponchon–Savarit method should be used.
Solution
From a material balance for a feed of 100 kmol,
= = 150; = 112.5; = 212.5; D = 37.5 and W = 62.5 kmol
The equilibrium curve and the diagonal line are drawn in as shown in Figure 12.
The equation of the top operating line is: = 0.75 +1 + 0.225
Thus, the line AB is drawn through the two points A (0.9, 0.9) and B (0, 0.225).
The equation of the bottom operating line is: = 1.415 +1 − 0.042
This equation is represented by the line CD drawn through C (0.1, 0.1) at a slope of 1.415.
Starting at point A, the horizontal line is drawn to cut the equilibrium line at point 1. The vertical line is dropped
through 1 to the operating line at point 2 and this procedure is repeated to obtain points 3–6.
A horizontal line is drawn through point 6 to cut the equilibrium line at point 7 and a vertical line is drawn
through point 7 to the lower enrichment line at point 8. This procedure is repeated in order to obtain points 9–16.
The number of stages is then counted, that is points 2, 4, 6, 8, 10, 12, and 14 which gives the number of plates
required as 7.
15
Figure 12: Determination of Number of Plates by the McCabe–Thiele Method for Example 5
= + ………………………………………………………………………………………(21)
and: = − …………………………………………………………………………….…(22)
or − = − − ( )………………………………………………………………..(23)
A material balance over the feed plate gives,
+ + = +
or:
− = − − ………………………………………………………………………………(24)
16
To obtain a relation between and , it is necessary to make an enthalpy balance over the feed plate, and to
consider what happens when the feed enters the column. If the feed is all in the form of liquid at its boiling point,
the reflux Lm overflowing to the plate below will be + . If however the feed is a liquid at a temperature ,
that is less than the boiling point, some vapour rising from the plate below will condense to provide sufficient
heat to bring the feed liquor to the boiling point. If is the enthalpy per mole of feed, and is the enthalpy of
one mole of feed at its boiling point, then the heat to be supplied to bring feed to the boiling point is ( - ),
and the number of moles of vapour to be condensed to provide this heat is ( - )/ , where λ is the molar
latent heat of the vapour. The reflux liquor is then:
= + + ( - )/ = + ( + - )/ = + ……………………………….(25)
ℎ 1
where = ℎ ℎ
Thus, − = − ……………………………………………………………………………….(26)
A material balance of the more volatile component over the whole column gives:
= +
Thus, from equation 23: −1 = −
or: = − ………………………………………………………………………….......(27)
−1 −1
This equation is commonly known as the equation of the q-line. If xq = xf, then yq = xf . Thus, the point of
intersection of the two operating lines lies on the straight line of slope q/(q − 1) passing through the point (xf , xf ).
When yq = 0, xq = xf /q. The line may thus be drawn through two easily determined points. From the definition of
q, it follows that the slope of the q-line is governed by the nature of the feed as follows.
Altering the slope of the q-line will alter the liquid concentration at which the two operating lines cut each other
for a given reflux ratio. This will mean a slight alteration in the number of plates required for the given separation.
Whilst the change in the number of plates is usually rather small, if the feed is cold, there will be an increase in
reflux flow below the feed plate, and hence an increased heat consumption from the boiler per mole of distillate.
Figure 13: Effect of the Condition of the Feed on the Intersection of the Operating Lines
for a Fixed Reflux Ratio
= +1 + …………………………………………………………………………..(28)
+ +
= +1 + ………………………………………………………………………….……..(29)
+1 +1
Any change in the reflux ratio R will therefore modify the slope of the operating line and, as may be seen from
Figure 11.15, this will alter the number of plates required for a given separation. If R is known, the top line is
most easily drawn by joining point A (xd, xd ) to B (0, xd/(R + 1)) as shown in Figure 14. This method avoids the
calculation of the actual flow rates Ln and Vn, when the number of plates only is to be estimated.
If no product is withdrawn from the still, that is D = 0, then the column is said to operate under conditions of total
reflux and, as seen from equation 11.47, the top operating line has its maximum slope of unity, and coincides with
the line x = y. If the reflux ratio is reduced, the slope of the operating line is reduced and more stages are required
to pass from xf to xd, as shown by the line AK in Figure 14. Further reduction in R will eventually bring the
operating line to AE, where an infinite number of stages is needed to pass from xd to xf . This arises from the fact
that under these conditions the steps become very close together at liquid compositions near to xf , and no
enrichment occurs from the feed plate to the plate above. These conditions are known as minimum reflux, and the
reflux ratio is denoted by Rm. Any small increase in R beyond Rm will give a workable system, although a large
number of plates will be required. It is important to note that any line such as AG, which is equivalent to a smaller
value of R than Rm, represents an impossible condition, since it is impossible to pass beyond point G towards xf.
Two important deductions may be made. Firstly that the minimum number of plates is required for a given
separation at conditions of total reflux, and secondly that there is a minimum reflux ratio below which it is
impossible to obtain the desired enrichment, however many plates are used.
18
Figure 14: Influence of Reflux Ratio on the Number of Plates Required for a Given Separation
If the q -line is horizontal as shown in Figure 11.18, the enrichment line for minimum reflux is given by AC,
where C has coordinates (xc, yc). Thus:
−
=
+1 −