 For a non-ideal system, where the molar latent heat is

no longer constant and where there is a substantial

heat of mixing, the calculations become much more
tedious.
 For binary mixtures of this kind a graphical model has
been developed by , PONCHON, and SAVARIT, based
on the use of an enthalpy-composition chart.
 It is necessary to construct an enthalpy-composition
diagram for particular binary system over a
temperature range covering the two-phase vapor-
liquid region at the pressure of the distillation.
The following data are needed:

1. Heat capacity as a function of temperature,

composition and pressure.
2. Heat of mixing and dilution as a function of
temperature and composition.
3. Latent heats of vaporization as a function of
composition and pressure or temperature.
4. Bubble-point temperature as a function of
composition and pressure.
Enthalpy of liquid:
n
hm   hi (1)
i

In “regular” / ideal mixtures:

hi  x i hio
(2)

For gaseous / vapor mixtures at normal T and P:

n n
H m   hi   xi i (3)
i i
Enthalpy-composition diagram
The equations used to calculate enthalpy of liquid and
vapor are:

hmix  x A hA  x B hB  Hsol (3)

hmix  x A CP ,A T  Tref   x B CP ,B T  Tref   Hsol (4)

Hmix  hmix   mix (5)

mix  x A  A  xB B (6)
EXAMPLE 1
Devise an enthalpy-concentration diagram for the
heptane-ethyl benzene system at 760 mm Hg, using the
pure liquid at 0C as the reference state and assuming
zero heat of mixing.

SOLUTION
heptane ethyl benzene

TB (C) 136.2 98.5

CP (cal/mole K) 51.9 43.4
 (cal/mole) 7575 8600
 t, C xH yH
136.2 0.000 0.000
129.5 0.080 0.233
122.9 0.185 0.428
119.7 0.251 0.514
116.0 0.335 0.608
110.8 0.487 0.729
106.2 0.651 0.834
103.0 0.788 0.904
100.2 0.914 0.963
98.5 1.000 1.000
t = 136.2C
xH = 0.0 xEB = 1.0

h  x H CP ,H t  tref   x EB CP ,EB t  tref   Hsol

= 0 + 1.0 (43.4) (136.2 – 0) + 0 = 5,911 cal/mole

mix = xH H + xEB EB = 0 + 1.0 (8,600) = 8,600

H = h + mix = 5,911 + 8,600 = 14,511 cal/mole

t = 129.5C
xH = 0.08 xEB = 0.92

h  x H CP ,H t  tref   x EB CP ,EB t  tref   Hsol

= 0.08 (51.9) (129.5) + 0.92 (43.4) (129.5)

= 5,708 cal/mole
mix = xH H + xEB EB = 0.08 (7,575) + 0.92 (8,600) = 8,518

H = h + mix = 5,708 + 8,518 = 14,226 cal/mole

The computations are continued until the last point

where xH = 1.0 and xEB = 0.0
t, C xH yH h H
(cal/mole) (cal/mole)

136.2 0.000 0.000 5,911 14,511

129.5 0.080 0.233 5,708 14,226

122.9 0.185 0.428 5,527 13,937

119.7 0.251 0.514 5,450 13,793

116.0 0.335 0.608 5,365 13,621

110.8 0.487 0.729 5,267 13,368

106.2 0.651 0.834 5,197 13,129

103.0 0.788 0.904 5,160 12,952

100.2 0.914 0.963 5,127 12,790

98.5 1.000 1.000 5,112 12,687

 20,000
18,000
Vapor
16,000 H vs. yA (Saturated vapor)
14,000
12,000
H

10,000 2 Phase
8,000 tie line

6,000
4,000
2,000 Liquid h vs. xA (Saturated liquid)

0
0 0.2 0.4 0.6 0.8 1
x
 The material balance equation
QD
maybe rearranged in the from
of difference:

V1 D
L0 xD
L 0 – V1 = L 1 – V2 = L 2 – V3
V2
L1
V3  = . . . . = Ln – Vn+1
L2
=–D =

n L 0 – V1 = – D =  (35)
Vn+1
Ln

LF
For the component material balance:

L0 x0 – V1 y1 = L1 x1 – V2 y2
= L2 x2 – V3 y3 = . . . .
= Ln xn – Vn+1 yn+1
= – D x D =  x
L0 x0 – V1 y1 = – D xD =  x (36)

Combining eqs. (35) and (36):

xD = x (37)
For the enthalpy balance:

L0 h0 – V1 H1 = L1 h1 –V2 H2 = L2 h2 –V3 H3 = . . . .
(38)
= Ln hn – Vn+1 Hn+1 = – D (hD +QD) =  h

Combining eqs. (23) and (35):

h = h D + Q D QD=qD/D (39)

• These 3 independent equations [eqs. (35), (36), and (37)] can be

written for rectifying section of the column between each plate.
• On the enthalpy scale and on the composition scale, the
differences in enthalpy and in composition always pass through
the same point,  ([xD, (hD + QD ]
• This is designated as point , the difference point, and all lines
corresponding to the combined material and enthalpy balance
equations (operating line equations) for the rectifying section of
the column pass through this intersection.
PLATE-TO-PLATE GRAPHICAL PROCEDURE FOR DETERMINING
THE NUMBER OF EQUILIBRIUM STAGES:

1. Use R, xD, h  to establish the location of point  with x = xD

and h = (hD +QD)
2. Use Equilibrium data alone to establish the point L1 at (x1, h1).
Since L1 is assumed to be a saturated liquid, x1 must lie on the
saturated-liquid line.
3. Draw the operating line between L1 and . This line intersects
the saturated-vapor line at V2 (y2, H2).
4. Repeat steps 2 and 3 until the feed plate is reached.
   (x, h)

V4 V3 V2 V1
H or h

D
L2 L1
L3

x or y
m
V m 1 Lm


VM L M 1
N
V M 1 LM
qB
B
xB
The material balance equation maybe rearranged in
the from of difference:

L M  V M 1  B  L M 1  V M

 L M  2  L M 1  . . .

 L m  V m1   (40)
For the component material balance:
L M x M  V M  1 y M  1  B x B  L M 1 x M 1  V M y M

 L M  2 x M  2  V M 1 y M 1  . . .

 L m x m  V m 1 y m 1   x  (41)

Combining eqs. (40) and (42):

x   xB (42)
For the enthalpy balance: QB=qB/B

LM hM  V M 1 HM 1  B hB  QB   LM 1 hM 1  V M HM

 LM 2 hM 2  V M 1 HM 1  . . .

 L m hm  V m1 Hm1   h (43)

Combining eqs. (40) and (43):

h  hB  QB (44)
• These 3 independent equations [eqs. (40), (41), and
(43)] can be written for stripping section of the column
between each plate.
• On the enthalpy scale and on the composition scale,
the differences in enthalpy and in composition always
pass through the same point,  [xB, (hB – QB)].
• This is designated as point  , the difference point, and
all lines corresponding to the combined material and
enthalpy balance equations (operating line equations)
for the stripping section of the column pass through
this intersection.
QB is usually not known. It can be derived from over-all
material balance:

F=D+B (45)

Combining eq. (45) with eqs. (35) and (40) gives:

  F  (46)

Equation (46) implies that  lies on the extension of the

straight line passing through F and .
 Over-all material balance:
V1 qD
F=D+B (16)
D, xD, HLD
L0
x0
Component material balance:
HL0
F xF = D x D + B x B (17)
F
xF
HF F xF = D xD + (F – D) xB (18)
F x F  x B 
D (19)
xD  xB
qB
B, xB, HLB Enthalpy balance:
F hF  qB  D hD  qD  B hB (20)
PLATE-TO-PLATE GRAPHICAL PROCEDURE FOR
DETERMINING THE NUMBER OF EQUILIBRIUM STAGES:
1. Draw a straight line passing through F and .
2. Draw a vertical straight line at xB all the way down until it
intersects the extension of line F in 
3. Assuming the reboiler to be an equilibrium stage, the vapor
VM+1 is in equilibrium with the bottom stream.
4. Use equilibrium data alone to establish the value of ym+1 on
the saturated-vapor line.
5. Draw the operating line between Lm(xm, hm) and VM+1. This
line intersects the saturated-liquid line at 
6. Repeat steps 4 and 5 until the feed plate is reached.
 VM1 VM VM1

hB LM LM-1
H or h

x or y

  , xB 
 • The construction may start from
V1 qD
either side of the diagram,
D
indicating either the condition at
L0 the top or the bottom of the
column.
• Proceed as explained in previous

F slides.
• In either case, when an equilibrium
tie line crosses the line connecting
 the difference points through the
feed condition, the other difference
qB point is used to complete the
B construction.
 

H or h

2 1
9 8 7 6 4 3
5

F

 xF xD
xB
 V1

L1
H or h

y1

x1
EXAMPLE 2
Using the enthalpy-concentration diagram from Example 1,
determine the following for the conditions in Example 1, assuming a
saturated liquid feed.
a. The number of theoretical stages for an operating reflux ratio of
R = L0/D = 2.5
b. Minimum reflux ratio L0/D .
c. Minimum equilibrium stages at total reflux.
d. Condenser duty feeding 10,000 lb of feed/hr, Btu/hr. if xf =0.42,
xd =0.97 , xB =0.01 Reboiler duty, Btu/hr.
SOLUTION
(a) From the graph: hD = h0 = 5,117 cal/mole
H1 = 12,723 cal/mole

L0 hD  QD  H1

D H1  hD

5,117  QD  12,723
2.5 
12,723  5,117
QD = 26,621 cal/mole
The coordinate of point  is:
x = xD = 0.97

h = hD +QD = 5,117 + ( 26,621) = 31,738 cal/mole

• Draw a straight line passing through  and F.

• Extend the line until it intersects a vertical line passing
through xB, at 
• Draw operating lines and equilibrium lines in the whole
column using the method explained in the previous slides.

Number of stages = 11
t, C xH yH h H
(cal/mole) (cal/mole)

136.2 0.000 0.000 5,911 14,511

129.5 0.080 0.233 5,708 14,226

122.9 0.185 0.428 5,527 13,937

119.7 0.251 0.514 5,450 13,793

116.0 0.335 0.608 5,365 13,621

110.8 0.487 0.729 5,267 13,368

106.2 0.651 0.834 5,197 13,129

103.0 0.788 0.904 5,160 12,952

100.2 0.914 0.963 5,127 12,790

98.5 1.000 1.000 5,112 12,687

 35,000

30,000

25,000

20,000
H

15,000

10,000

5,000
F

0
0 0.2 0.4 0.6 0.8 1


  = 21,700 cal/mole
20,000 (b)
18,000
16,000
L0 hD  qD / D  H1
14,000 
D H1  hD
12,000
H

10,000
h  H1
8,000 
6,000 H1  hD
4,000
F
2,000  L0  21,700  12,723
  
0
1
0 0.2 0.4 0.6 0.8 1
 min 12,723  5,117
D
0.9

0.8
 = 1.18
0.7

0.6
y

0.5

0.4

0.3

0.2

0.1

0
0 0.2 0.4 0.6 0.8 1
x
 20,000
(c) 18,000
16,000
7 6
14,000 5 4 3 2 1
12,000
H 10,000
8,000
6,000
4,000 F
2,000
0
1 N=7
0 0.2 0.4 0.6 0.8 1
0.9

0.8

0.7

0.6
y

0.5

0.4

0.3

0.2

0.1

0
0 0.2 0.4 0.6 0.8 1
x
 hD + QD
(d)
= h = 31,738 cal/mole

hD = 5,117 cal/mole

QD = 26,621 cal/mole
cal 1.8 Btu lb mole  mole D  10,000 lb F hr
QD  26,621    0.426 
mole cal mole  mole F  103 lb mole F

= 1,981,843 Btu/hr

(e) hB – QB = 14,350 cal/mole

hB = 5,886 cal/mole

-QB = 14,350 + 5,886 = -20,236 cal/mole

cal 1.8 Btu lb mole  mole B  10,000 lb F hr

QD  20,236  
 0.574 
  103 lb mole F
mole cal mole  mole F 

= 2,631,751 hr
 A
qD Material balance around condenser:
V1 = L 0 + D (21)
V1 D
L0 xD
Component material balance:
L1
V1 y1 = L0 x0 + D x0 (22)

Enthalpy balance:
V1 H1 = L0 h0 + D hD + qD (23)
 qD
Designating: QD 
D
V1 H1 = L0 h0 + D (hD + QD) (24)

Combining eqs. (21) and (24):

L0 hD  QD  H1
 (25)
D H1  hD

Internal reflux is shown as:

L0 hD  QD  H1

V1 hD  QD  h0 (26)
 A
qD
Internal reflux between
each plate, until a point
in the column is reached D
where a stream is added L0 xD
or removed, can be Ln
n
shown as: Vn+1

Ln hD  QD  H n1

Vn1 hD  QD  hn (27)
 (hD + QD), xD


hD + QD – H1

H1

H or h

V1
H1 – hD

h0, hD
L0, D 

y1, x0, xD
Min Number of Stages
Min. Reflux Ratio
Thanks