Distillation 1
Distillation 1
UNIT 2
DISTILLATION
OUTLINE
Introduction to distillation- basics and theory
Vapor liquid equilibrium
Phase Diagram
X-Y Diagram
Azetrope
Batch distillation
Flash distillation
Continuous distillation
Tray distillation column
Analysis of McCabe and thiele method for binary system
Analysis of Ponchon-savariat method for binary system
Tray efficiency
Packed distillation column
Approximate column sizing
DEPARTMENT OF CHEMICAL AND
PROCESS SYSTEM ENGINEERING/
2
HIT/ HARARE
INTRODUCTION
Distillation is a method of separation of components
from a liquid mixture which depends on the
differences in boiling points of the individual
components.
CLAUSIUS-CLAPEYRON EQUATION:
Where pV and p1V are the vapour pressures
in Pascal at absolute temperature T and T1 in
K. λ is the molar latent heat of vaporization
which is independent of temperature.
ANTOINE EQUATION:
The rate of species i in the distillate flowing out of the batch distillation
system, xD dD, is equal to the decreasing rate of species i the liquid in the
still pot -(d Wxw)
Graphically determine the area under the curve from xo, to xw,final
or perform the integration numerically
After the numerical value of the integration is obtained, the value of Wfinal
(i.e. the amount of liquid remained in the still pot) can be obtained
Boil the mixture, condense the vapour and collect the distillate.
Repeat the procedure until the desired purity is obtained.
EACH STILL IS A STEP ON THE X-Y CURVE+
McCabe-Thiele (1925)
in the present:
complex mathematical matrix methods allow to find exact solutions of nonlinear
equation systems: Wang-Henke (1966),
The component heat capacity changes and the heat of mixing are
negligible compared to the heat of vaporization (considering ideal
behaviour of binary mixtures).
The distillation column, the condenser and the reboiler are well
insulated so that heat losses to environment are negligible.
That means, that all liquid and vapor molar flow rates in the
rectifying section are constant and that all liquid and vapor molar
flow rates in the stripping section are constant but not the same
as those in the rectifying section.
DEPARTMENT OF CHEMICAL AND
PROCESS SYSTEM ENGINEERING/
54
HIT/ HARARE
DESIGN PROCEDURE McCabe-Thiele Method
The steps to be followed to determine the number of theoretical stages by
McCabe-Thiele Method:
1. Start with VLE data and draw an x-y diagram
2. On the x-y diagram draw the line y=x or the 45oC line
3. Determine the Rectifying section operating line (ROL) and draw it on the
x-y diagram
4. Determine the feed condition and use it for the determination of the feed
section operating line (q-line) and draw it on the x-y diagram
5. Determinate the stripping section operating line (SOL) and draw it on
the x-y diagram
6. Using the ROL and the equilibrium curve draw steps from the distillate
composition to the feed point.
7. when you reach point of intersection of q-line, ROL and SOL start stepping
using equilibrium curve and SOL until you reach the required composition
of the bottoms.
8. Count the number of steps. Each one equals an equilibrium or theoretical
stage
9. Using the stripping operating line and the equilibrium curve draw steps
from the bottoms composition to the feed point
10. Convert this to a number of plates by dividing by the plate efficiency.
The stripping operating line (SOL) which is a straight line with slope ( L' / V' ) and
intercept ( B𝑥𝐵 / V'). When 𝑥 = 𝑥𝐵 , y =𝑥𝐵 , the SOL passes through (𝑥𝐵 , 𝑥𝐵 ) on the
45 o diagonal line.
DEPARTMENT OF CHEMICAL AND
PROCESS SYSTEM ENGINEERING/
59
HIT/ HARARE
STRIPPING LINE ON X-Y DIAGRAM
SOL is constructed starting from desired bottom product ( xB )
composition. A vertical line is drawn from this point to the
diagonal line, and a line of slope L/V is drawn as illustrated in
the diagram below.
368.2 − 363.2
𝑇𝑏.𝑝 𝑚𝑖𝑥𝑡𝑢𝑟𝑒 𝑎𝑡 𝑥𝐹 = 368.2𝐾 − 0.411 − 0.45 = 367𝐾
0.411 − 0.581
Thus feed enters at a temperature below boiling point of mixture that’s
our feed condition is sub-cooled liquid.
𝐶𝑝,𝐿 𝑇𝑏𝑝 − 𝑇𝐹 + 𝜆
For cold feed 𝑞=
𝜆
Where Tb= 367K and Tf= 327.6K inlet feed temperature
159(367 − 327.6)
𝑞 =1+ = 1.195
13800 DEPARTMENT OF CHEMICAL AND
PROCESS SYSTEM ENGINEERING/
77
HIT/ HARARE
SOLUTION CONT’:
Slope of q line is:
𝑞 1.195
= = −6.12
𝑞;1 1.195;1
DEPARTMENT OF CHEMICALAND
PROCESS SYSTEM ENGINEERING/
85
HIT/ HARARE
MINIMUM REFLUX
using definition of slope of line = R/(R+1)
Where x’, y’ is the intersection
of line with curve
xA yA
1.000 1.000
0.780 0.900
0.581 0.777
0.411 0.632
0.258 0.456
0.130 0.261
0 0
Material balance
overall mass balance: F = D +W
204 = D + W……………1
= 1 − 0.611 = 0.389
𝑞 0.389
The slope of the q-line is: = = −0.637
𝑞;1 0.389;1
𝑉𝑛 − 𝐿𝑛:1 = 𝐷 ……………..2.4.1
Component balance:
𝑉𝑛 𝑦𝑛 − 𝐿𝑛:1 𝑥𝑛:1 = 𝐷𝑥𝐷 ……….….2.4.2
The left hand side net flow upwards of the more volatile component. Since for a
given distillation the right-hand side is constant it follows that the difference or net
rate of flow upward is constant, independent of tray number in
this section of tower and equal to that
DEPARTMENT OF CHEMICALAND
permanently withdrawn at the top. PROCESS SYSTEM ENGINEERING/
113
HIT/ HARARE
Enthalpy balance:
𝑉𝑛 𝐻𝑉𝑛 + 𝐿𝑛:1 𝐻𝐿𝑛:1 = 𝑄𝐶 + 𝐷𝐻𝐷 ……….2.4.3
𝑄𝐶
𝑉𝑛 𝐻𝑉𝑛 − 𝐿𝑛:1 𝐻𝐿𝑛:1 = 𝐷 + 𝐻𝐷 ……..2.4.4
𝐷
𝑄𝐶
Let , 𝑄′ = + 𝐻𝐷 ……………….2.4.5
𝐷
𝑉𝑛 𝐻𝑉𝑛 − 𝐿𝑛:1 𝐻𝐿𝑛:1 = 𝐷𝑄′ ……….2.4.6
Left-hand side is the net flow of heat upward. Since for a given distillation the right-
hand side is constant it follows that the difference or net rate of flow of heat upward is
constant, independent of tray number in this section of tower and equal to that
permanently taken out at the top with distillate and at condenser.
Elimination of D between Eqs. 2.41 and 2.4.2 and between 2.4.1 and
2.4.6 yields
𝐿𝑛+1 𝑥𝐷 ;𝑦𝑛 𝑄′ ;𝐻𝑉𝑛
= = …………………………..2.4.7
𝐷 𝑦𝑛 ;𝑥𝑛+1 𝐻𝑉𝑛 ;𝐻𝐿𝑛+1
A stream ∆𝐷 may be defined with mass equal to the difference between
the vapour and liquid streams between two plates and it has composition
𝑥𝐷 and enthalpy 𝑄′ . Then the three quantities
𝑉𝑛 ,𝐿𝑛:1 , and 𝑄′ on a straight line passing through
DEPARTMENT OF CHEMICALAND
𝑄′ as shown in figure below. PROCESS SYSTEM ENGINEERING/
114
HIT/ HARARE
Loop II
Below the feed plate a similar series of equations for material
and heat balance can be written
Total material balance:
𝐿𝑚:1 − 𝑉𝑚 = 𝐵 ……………………..2.4.8
𝐿𝑚:1 𝑥𝑚:1 − 𝑉𝑚 𝑦𝑚 = 𝐵𝑥𝐵 ………….2.4.9
𝐿𝑚:1 𝐻𝐿𝑚:1 − 𝑉𝑚 𝐻𝑉𝑚 = 𝐵𝐻𝐵 − 𝑄𝐵 ……..2.4.10
𝑄𝐵
𝐿𝑚:1 𝐻𝐿𝑚:1 − 𝑉𝑚 𝐻𝑉𝑚 = 𝐵 𝐻𝐵 − …….2.4.11
𝐵
′ 𝑄𝐵
Let , 𝑄 ′ = 𝐻𝐵 − ……….….2.4.12
𝐵
𝐿𝑚:1 𝐻𝐿𝑚:1 − 𝑉𝑚 𝐻𝑉𝑚 = 𝐵𝑄 ′′ ………2.4.13
𝐿𝑚+1 ;𝑥𝐵 : 𝑦𝑚 ;𝑄 ′′ : 𝐻𝑉𝑚
= = ……..2.4.14
𝐵 𝑦𝑚 ; 𝑥𝑚+1 𝐻𝑉𝑚 ;𝐻𝐿𝑚+1
DEPARTMENT OF CHEMICAL AND
PROCESS SYSTEM ENGINEERING/
115
HIT/ HARARE
Eqn. 2.4.14 represents any operating line below the
feed plate, and shows that all such lines pass through
a common pole ΔB of coordinates xB and 𝑄′′ . As with
the rectifying section a stream
ΔB may be defined by mass 𝐿𝑚:1 − 𝑉𝑚 and its
composition xB and enthalpy𝑄′′ .
Then
𝐹 = ∆𝐷 + ∆𝐵 …………………………..2.4.15
𝐹𝑥𝐹 = ∆𝐷 𝑥𝐷 + ∆𝐵 𝑥𝐵 …………………2.4.16
DEPARTMENT OF CHEMICALAND
PROCESS SYSTEM ENGINEERING/
129
HIT/ HARARE
PROBLEM 2
A feed containing 40 mole percent n-hexane and 60 percent n-octane is fed
to a distillation column. A reflux ratio of 1.2 is maintained. The overhead
product is 95 percent hexane and the bottoms 10 percent hexane. Find the
number of theoretical stages and the optimum feed stage. Assume that a
total condenser is used. The column is to operate at 1 atm.
Equilibrium data is collected.
VLE Data, Mole Fraction Hexane, 1 atm
x (liquid) 0.0 0.1 0.3 0.5 0.55 0.7 1.0
y (vapor) 0.0 0.36 0.70 0.85 0.90 0.95 1.0
Enthalpy-Concentration Data
Mole Fraction Hexane Enthalpy cal/gmol
Sat. Liquid Sat. Vapor
0.0 7000 15,700
0.1 6300 15,400
0.3 5000 14,700
0.5 4100 13,900
0.7 3400 12,900
0.9 3100 11,600
1.0 3000 10,000
DEPARTMENT OF CHEMICAL AND
PROCESS SYSTEM ENGINEERING/
HIT/ HARARE 130
PROBLEM 2
1000 kg/hr of a mixture containing 42 mole percent heptane and 58 mole
percent ethyl benzene is to be fractionated to a distillate containing 97
mole percent heptane and a residue containing 99 mole percent ethyl
benzene using a total condenser and feed at its saturated liquid
condition. The enthalpy-concentration data for the heptane-ethyl
benzene at 1 atm pressure are as follows:
xheptane 0 0.08 0.18 0.25 0.49 0.65 0.79 0.91 1.0
Hl (kJ/kmol) x 10-3 24.3 24.1 23.2 22.8 22.05 21.75 21.7 21.6 21.4
Hv (kJ/kmol) x 10-3 61.2 59.6 58.5 58.1 56.5 55.2 54.4 53.8 53.3
Typical values for tray efficiency are 0.5 to 0.75, i.e. they are
0.5 to 0.75 times as effective as an ideal stage.
Divide the number of ideal stages by the efficiency, add 10%
extra trays.
The overall efficiency is extremely easy to measure and use;
thus, it is the most commonly used efficiency value in the
plant.
However, the overall efficiency is not representative of
column operation because the different compositions on the
various trays result in different tray efficiencies.
References
DEPARTMENT OF CHEMICALAND
PROCESS SYSTEM ENGINEERING/
142
HIT/ HARARE
Determination of the column height
1. The theoretical stage method
2. The transfer unit method
𝑘𝑥 𝑎 𝑦1 ;𝑦
Rearranging − =
𝑘𝑦 𝑎 𝑥1 ;𝑥
DEPARTMENT OF CHEMICALAND
PROCESS SYSTEM ENGINEERING/
148
HIT/ HARARE
Material balance over an incremental column
height, h with constant molar overflow:
𝑉𝑑𝑦 = 𝑘𝑦 𝑎 𝑦1 − 𝑦 𝑠𝑑 Where s is the cross- sectional area
over an incremental column height, h
L𝑑𝑥 = 𝑘𝑥 𝑎 𝑥 − 𝑥1 𝑠𝑑 with constant molar overflow
Integrating over the rectifying section
𝑦2 𝑥𝐷
𝑉𝑑𝑦 𝑉𝑑𝑦
𝑅 = 𝑑 = =
𝑦𝐹 𝑘𝑦 𝑎𝑆 𝑦1 − 𝑦 𝑥𝐹 𝑘𝑥 𝑎𝑆 𝑥 − 𝑥1
𝑦2 𝑥𝐷
𝐻𝐺 𝑑𝑦 𝐻𝐿 𝑑𝑦
𝑅 = 𝑑 = =
𝑦𝐹 𝑦1 − 𝑦 𝑥𝐹 𝑥 − 𝑥1
Integrating over the stripping section
𝑦𝐹 𝑥𝐹
𝑉𝑑𝑦 𝑉𝑑𝑦
𝑆 = 𝑑 = =
𝑦1 𝑘𝑦 𝑎𝑆 𝑦1 − 𝑦 𝑥1 𝑘𝑥 𝑎𝑆 𝑥 − 𝑥1
𝑦𝐹 𝑥𝐹
𝐻𝐺 𝑑𝑦 𝐻𝐿 𝑑𝑦
𝑅 = 𝑑 = =
𝑦1 𝑦1 − 𝑦 𝑥1 𝑥 − 𝑥1
DEPARTMENT OF CHEMICAL AND
PROCESS SYSTEM ENGINEERING/
HIT/ HARARE 149
Values of ky and kx vary over the packed height,
causing the slope (kxa/kya) to vary.
DEPARTMENT OF CHEMICALAND
PROCESS SYSTEM ENGINEERING/
152
HIT/ HARARE
Column Height
The bottom of the tower must be tall enough to serve as a liquid
reservoir. Depending on your boss's feelings about keeping
inventory in the column, you will probably design the base for about
5 minutes of holdup, so that the total material entering the base can
be contained for at least 5 minutes before reaching the bottom tray.
The total of height added to the top and bottom will usually amount
to about 15% or so added to that required by the trays.
You rarely will see a real tower that is more than about 175 ft. tall.
Tall, skinny towers are not a good idea, so watch the
height/diameter ratio.
You generally want to keep it less than 20 or 30. If your tower ends
up exceeding these values, you probably want to look at a redesign,
maybe by reducing the tray spacing, or splitting the tower into two
parts.
DEPARTMENT OF CHEMICALAND
PROCESS SYSTEM ENGINEERING/
153
HIT/ HARARE
Diameter of tray towers
The flooding condition fixes the upper limit of vapour
velocity. A high vapour velocity is needed for high plate
efficiencies, and the velocity will normally be between 70
to 90 per cent of that which would cause flooding.
For design, a value of 80 to 85 per cent of the flooding
velocity should be used. The flooding velocity can be
estimated from the correlation given by Fair (1961):
𝜍 0.2 𝜌𝐿 − 𝜌𝑉
𝑢𝑓𝑙𝑜𝑜𝑑𝑖𝑛𝑔 = 𝐾𝑉
20 𝜌𝑉
DEPARTMENT OF CHEMICALAND
PROCESS SYSTEM ENGINEERING/
154
HIT/ HARARE
DEPARTMENT OF CHEMICAL AND
PROCESS SYSTEM ENGINEERING/
155
HIT/ HARARE
4𝑉𝑚𝑎𝑠𝑠
𝐷𝑇 =
𝐴𝑑
𝑢𝑓𝑙𝑜𝑜𝑑𝑖𝑛𝑔 𝑓𝑑𝑒𝑠𝑖𝑔𝑛 𝜋𝜌𝑣 1 − 𝐴
DEPARTMENT OF CHEMICALAND
PROCESS SYSTEM ENGINEERING/
156
HIT/ HARARE
STUDY QUESTIONS
1. What equipment is included in a typical distillation operation?
2. What determines the operating pressure of a distillation column?
3. Under what conditions does a distillation column have to operate under vacuum?
4. Why are distillation columns arranged for counter-current flow of liquid and vapour?
5. Why is the McCabe–Thiele graphical method useful in this era of more rigorous,
computer-aided algebraic methods used in process simulators?
6. Under what conditions does the McCabe–Thiele assumption of constant molar overflow
hold?
7. In the McCabe–Thiele method, between which two lines is the staircase constructed?
8. What is meant by the reflux ratio? What is meant by the boil-up ratio?
9. What is the q-line and how is it related to the feed condition?
10. What are the five possible feed conditions?
11. In the McCabe–Thiele method, are the stages stepped off from the top down or the
bottom up? In either case, when is it best, during the stepping, to switch from one
operating line to the other? Why?
12. Can a column be operated at total reflux? How?
13. How many stages are necessary for operation at minimum reflux ratio?
14. What is meant by a pinch point? Is it always located at the feed stage?
15. What is meant by sub-cooled reflux? How does it affect the amount of reflux inside the
column?
16. Is it worthwhile to preheat the feed to a distillation column?
17. Why is the stage efficiency in distillation higher than that in absorption?
18. What kind of a small laboratory column is useful for obtaining plate efficiency data?
DEPARTMENT OF CHEMICALAND
PROCESS SYSTEM ENGINEERING/
157
HIT/ HARARE