0% found this document useful (0 votes)
565 views47 pages

Ponchon Savarit Method

The document describes how to construct an enthalpy-composition diagram for distillation systems involving non-ideal mixtures, which requires data on heat capacities, heats of mixing and vaporization, and bubble points as functions of temperature, composition, and pressure. It presents equations to calculate liquid and vapor enthalpies and develops an example enthalpy-concentration diagram for a heptane-ethyl benzene system. The diagram can be used to evaluate streams in the distillation column graphically.

Uploaded by

Bench Gueco
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
565 views47 pages

Ponchon Savarit Method

The document describes how to construct an enthalpy-composition diagram for distillation systems involving non-ideal mixtures, which requires data on heat capacities, heats of mixing and vaporization, and bubble points as functions of temperature, composition, and pressure. It presents equations to calculate liquid and vapor enthalpies and develops an example enthalpy-concentration diagram for a heptane-ethyl benzene system. The diagram can be used to evaluate streams in the distillation column graphically.

Uploaded by

Bench Gueco
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 47

 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 RUHEMANN, 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
Hm = hi + y i i (3)
i i
Enthalpy-composition diagram
The equations used to calculate enthalpy of liquid and
vapor are:

hmix = x A hA + xB 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 + x B B (6)
EXAMPLE 2
Devise an enthalpy-concentration diagram for the
heptane-ethyl benzene system at 760 mm Hg, using the
pure liquid at 0C 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.2C
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,920 + 8,600 = 14,511 cal/mole


t = 129.5C
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 h H
(cal/mole) (cal/mole)
136.2 0.000 5,911 14,511
129.5 0.080 5,708 14,226
122.9 0.185 5,527 13,937
119.7 0.251 5,450 13,793
116.0 0.335 5,365 13,621
110.8 0.487 5,267 13,368
106.2 0.651 5,197 13,129
103.0 0.788 5,160 12,952
100.2 0.914 5,127 12,790
98.5 1.000 5,112 12,687
20,000
18,000
Vapor
16,000 Saturated vapor
14,000
12,000
H

10,000 2 Phase
8,000
6,000
4,000
2,000 Liquid Saturated liquid

0
0 0.2 0.4 0.6 0.8 1
x
The enthalpy-concentration diagram may be used to
evaluate graphically the enthalpy and composition of
streams added or separated.
Over-all material balance:
V
F=V+L (7)
Component material balance
F q F xF = V y + L x (8)

Enthalpy balance:
F hF = V H + L h (9)
L
Steady-state flow system
with phase separation
and heat added
For adiabatic process, q = 0:
Substituting eq. (7) to (9) gives:
(V + L ) hF = V H + L h
V (H hF) = L (hF h) (10)
L H hF
= (11)
V hF h

Substituting eq. (7) to (8) gives:

V (y xF) = L (xF x) (12)


L y xF
= (13)
V xF x
Substituting eq. (12) to (13) gives:

H hF y x F (14)
=
hF h x F x

Eq. (14) can be rearranged:

H hF hF h
= (15)
y xF xF x
___ H hF
The slope of line VF is :
y xF
V
H
F
hF
___ hF h
The slope of line FL is :
h
L xF x

According to eq. (15), the


slopes of both lines are the
x xF y
same.
Enthalpy-
concentration lines Since both lines go through
adiabatic, q = 0 the same point (F), the lines
lie on the same straight line.
LEVER-ARM RULE PRINCIPLE
V L H hF
H =
V hF h
F A hF
Consider triangle LBV
___ ___
L FV AV H hF L
___ = ___ = =
B h
LF BA hF h V
___
x xF y
L FV
= ___
V LF
___ ___
V LF L FV
Similarly: = ___ = ___
F LV F LV
Over-all material balance:
V1 qC
F=D+B (16)
D, xD, hD
L0
x0 Component material balance:
h0
F xF = D xD + B xB (17)
F
xF
hF F xF = D xD + (F D) xB (18)
F (x F x B )
D= (19)
xD xB
qB
B, xB, hB Enthalpy balance:
F hF + qB = D hD + qD + B hB (20)
A
qC Material balance around condenser:
V1 = L 0 + D (21)
V1 D
L0 xD
Component material balance:
V1 y1 = L0 x0 + D xD (22)
L1

Enthalpy balance:
qC + V1 H1 = L0 h0 + D hD (23)
qC
Designating: Q CD =
D
V1 H1 = L0 h0 + D (hD QCD) (24)

Combining eqs. (21) and (24):

L0 hD QCD H1 (25)
=
D H1 hD

Enternal reflux is shown as:


L0 h D QCD H1
= (26)
V1 h D QCD h 0
A
qC
Internal reflux between
each plate, until a point
in the column is reached D
where a stream is L0 xD
added or removed, can Lm
m
be shown as: Vm+1

Lm hD QCD Hm+1
=
Vm+1 hD QCD hm
(27)
(hD QD), xD

hD QCD H1

H1

H or h

V1
H1 hD

h0, hD
L0, D

y1, x0, xD
D, yD Material balance:
A qC V1 = L 0 + D (28)
Component material balance:
v1 L0 y1 V1 = x0 L0 + yD D (29)
v2
L1 Enthalpy balance:
v3
L2 qCD + V1 H1 = L0 h0 + D HD (30)

Designating: Q CD =
q CD
vF
D
F LF-1
xF V1 H1 = L0 h0 + D (hD QCD) (31)
LF
Combining eqs. (28) and (30):
L0 HD QCD H1
= (32)
D H1 h0

Internal reflux is shown as:

L0 HD QCD H1 (33)
=
V1 HD QCD h0

Internal reflux between each plate, until a point in the


column is reached where a stream is added or removed,
can be shown as:
Lm HD QCD Hm+1
= (34)
Vm+1 HD QCD hm
(HD QCD), xD

HD QCD H1
V1 D


HD, yD
H or h

H1 h0


h0, x0

y1, x0, yD
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 xD = 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 QCD) = h

Combining eqs. (23) and (35):


h = hD QCD (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 QCD)]
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, HD or hD to establish the location of point with x = xD


and h = hD QCD or h = HD QCD
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
qS
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:

LM hM VM+1 HM+1 = B(hB QSB) = LM1 hM1 VM 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 QSB (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 QSB)].
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 .
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 VM1

hB LM LM-1
H or h

x or y

( , x B )
The construction may start from
V1 qC
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
qS difference point is used to
B complete the construction.

H or h

2 1
9 8 7 6 4 3
5

F

xF xD
xB
V1

H or h L1

y1

x1
EXAMPLE 3
Using the enthalpy-concentration diagram from Example 2,
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.
e. Reboiler duty, Btu/hr.

SOLUTION
(a) From the graph: hD = h0 = 5,117 cal/mole
H1 = 12,723 cal/mole
L 0 h D Q CD H1
=
D H1 h D

5,117 QD 12,723
2 .5 = QCD = 26,621 cal/mole
12,723 5,117
The coordinate of point is:
x = xD = 0.97

h = hD QCD = 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
35,000

30,000

25,000

20,000

H
15,000

10,000

5,000
F

1
0
0.9 0 0.2 0.4 0.6 0.8 1

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
= 21,700 cal/mole
20,000 (b)
18,000

16,000 L0 hD QCD 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
D min 12,723 5,117
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 5 4
14,000 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
(d) hD QCD = h = 31,738 cal/mole

hD = 5,117 cal/mole

QCD = 26,621 cal/mole

cal 1.8 Btu lb mole mole D 10,000 lb F hr


QCD = 26,621 0.426
mole cal mole mole F 103 lb mole F
= 1,981,843 Btu/hr

(e) hB QSB = 14,350 cal/mole

hB = 5,886 cal/mole

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

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


QCD = 20,236 0.574
mole cal mole mole F 103 lb mole F

= 2,631,751 BTU/hr

You might also like