Unit Operations of Chemical Engineering: January 2018
Unit Operations of Chemical Engineering: January 2018
Unit Operations of Chemical Engineering: January 2018
net/publication/322364940
CITATIONS READS
10 55,600
1 author:
Dennis Prieve
Carnegie Mellon University
129 PUBLICATIONS 5,602 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Dennis Prieve on 10 January 2018.
Unit Operations
of Chemical Engineering
by
Dennis C. Prieve
Department of Chemical Engineering
Carnegie Mellon University
Pittsburgh, PA 15213
An electronic version of this book in Adobe PDF® format was made available
to students of 06-361, Department of Chemical Engineering,
Carnegie Mellon University, Spring Semester, 2017.
Table of Contents
CHAPTER 1. COURSE INTRODUCTION .................................................................. 1
Syllabus Summary .........................................................................................................................1
Instructor ................................................................................................................................................ 1
Textbook.................................................................................................................................................. 1
Web Site .................................................................................................................................................. 1
Grading ................................................................................................................................................... 1
What Are “Unit Operations”? ......................................................................................................2
Objective of Course ................................................................................................................................ 2
Overview of ChE Curriculum ................................................................................................................. 3
CHAPTER 2. HEAT EXCHANGE EQUIPMENT ........................................................ 3
double-Pipe Hxer ...........................................................................................................................3
Shell-and-Tube Hxer .....................................................................................................................4
CHAPTER 3. FUNDAMENTALS OF HEAT TRANSFER: REVIEW .......................... 5
Log-Mean Driving Force ...............................................................................................................5
Co-Current versus Counter-Current Flow ........................................................................................... 12
Non-Linear Relation Between q(x) and ∆T(x)....................................................................................... 12
Resistances in Series ....................................................................................................................16
Overall Coefficients (U) vs. Film Coefficients (h) ................................................................................ 16
Fouling Factors .................................................................................................................................... 20
Correlations to Estimate h’s .......................................................................................................21
hi for Turbulent Flow In Long Pipes .................................................................................................... 22
Effect of Duct Shape on hi .................................................................................................................... 23
Nomenclature .............................................................................................................................237
Rectangular Diagrams ...............................................................................................................237
Flash Mixing of Two Liquids ....................................................................................................238
Single-Stage Extractor ...............................................................................................................240
Multi-Stage Extractor ................................................................................................................242
Operating Lines for Extractors ........................................................................................................... 243
Difference-Point Method .................................................................................................................... 244
Minimum Solvent Flowrate ................................................................................................................. 246
TABLE OF CONTENTS.................................................................................................. II
INDEX VI
SYLLABUS SUMMARY
Instructor
My name is Dennis Prieve and my office is DH 3120. I encourage each of you to come and
see me when you have any problems. Right after class is the best time; right before class is the
worst time (I’m usually preparing for class). If you don’t find me in my office, contact me by
email.
Textbook
Web Site
• lecture notes
• homeworks and solutions
• old exams from previous years
• Mathcad modules to help with homeworks
Grading
Your final letter grade will be assigned based primarily on your numerical average of the
following:
50% Two in-class exams (Feb. 22 and Apr. 5 – both are Wednesdays)
25% Final Exam (TBA)
10% Homework
10% Project Report
5% Attendance
Most of the grade will come from three exams: two in-class exams and a final exam given during
final-exam week. They will count 25% each. All of these exams are open-book and open-notes.
You will find examples of old exams on BlackBoard which can also be used as practice exams.
We will be having 11 weekly homework assignments usually due on Wednesday. The first
homework will be due a week from today and I’ll have more to say about homeworks in a later
lecture. Collectively, the homeworks will count 10% of the final average.
Attendance will count 5%. I consider regular attendance to be required for learning by the vast
majority of students. Of course, I also expect your full attention during class.
The remaining 10% comes from a term paper or project report which will be due near the end of
the semester. This is a 10-page written summary of some unit operation of your choosing, not
one of the unit ops covered in lectures. We will not have a homework assignment collected one
week so you can work on this project.
Objective of Course
The objective of this course is to learn the principles needed to size equipment used for
physical transformation of process streams. In particular, we will begin to appreciate the role of
mass transfer and phase equilibrium. Equipment which we will learn how to size includes:
In this course, we will apply the principles learned in the engineering science courses you
have had during the last three semesters; these will be applied to the design of equipment for
physical transformations. We might divide up the courses in the curriculum as follows:
Course Topics
intro conservation of mass
1) conservation of energy
thermo
2) criteria for phase/chemical equilibria
DOUBLE-PIPE HXER
One particularly simple heat-exchange device is the double-pipe heat exchanger, shown in
the figure at right. It is called this because the basic component is constructed from two pipes,
one inside the other. One fluid (Fluid B)
flows inside the inner pipe while the second
fluid (Fluid A) flows in the annular space
between the two pipes. Of course, one fluid is
hotter than the other and heat flows through
the wall of the inner pipe from the hot fluid to
the cold fluid. To obtain larger heat exchange
area, several pipes are arranged side-by-side
and fittings are attached to allow the fluids to
contact the pipes in series so the overall flow
is the same as if only one very long pipe were
used rather than several shorter ones.
SHELL-AND-TUBE HXER
If very large heat exchange areas are required, the double-pipe heat exchanger design
becomes impractical. Instead, the shell-and-tube heat exchanger design is used.
Shown above is a bundle of small-diameter tubes which are arranged parallel to each other and
reside inside a much larger-diameter tube called the “shell”, much like strands of uncooked
spaghetti come in a tube-shaped container. The tubes are connected to a common space (called a
“manifold”) at either end so that the “tube fluid” enters the left side and is distributed equally
among all the tubes. At the right side, the fluid exits from each tube, is mixed together in a
second manifold, then leaves as a single stream. The second fluid, called the “shell fluid”, flows
in the space in between the outside of tubes. Baffle plates inside the shell force the shell fluid to
flow across the tubes repeatedly as the fluid moves along the length of the shell.
In the design above (called a 1-1 hxer), the flow is from left to right inside the tubes and from
right to left for the shell-side fluid. Another common design is shown below (called a 1-2 hxer).
Half of the tubes have flow from left to right while the other half have flow in the opposite
direction.
=
q UA ∆T (1)
♠ An equal sign inside square brackets (i.e. “[=]”) should be read “has units of” rather than “equal to.”
* Mathcad documents and other files referred to in the Notes are usually available on the Blackboard website for this
course: under Handouts / Student Resources.
p j y p
R
Tho = 344.3 K <-----<-----<-----<-----<-----<-----<-----<-----<----- Thi = 377.6 K x
n This ratio should be nearly unity if integration is
400 = 0.99930
LT accurate
∆T1 = 6.402 K These two ∆T's must both be positive for the
simulation to be realistic; if one is negative,
350
Thi ∆T2 = 55.5 K change mc or mh; if both are negative check
the T's
Tci
300 3 kg
mh = 7.26 × 10 Thi = 377.6 K Tho = 344.3 K
hr
3 kg
mc = 2 × 10 Tci = 288.8 K Tco = 371.198 K
250 hr
0 0.2 0.4 0.6 0.8
Tci = 288.8 K ----->----->----->----->----->----->----->----->----->
xi Tco = 371.198 K watt D = 5 in L = 32.314 m
U = 653
2 ∆Tlm = 22.733 K
LT m ⋅K
Since ∆T is a continuous function of position inside the heat exchanger, we will have to
integrate a differential form of this equation. At any axial position x inside the heat exchange,
the local heat flux (rate per unit area) is given by
dq
heat
= = U ( x ) ∆T ( x )
flux
(2)
J
dA
m2 − s
where the heat transfer coefficient as well as driving force in general depend on x.
Mathematically, dq is
We are tempted to substitute this expression for dA into (2) and try to integrate both sides with
respect to x. This turns out to difficult because functions U(x) and especially ∆T(x) are unknown
functions. A better approach turns out to be solving the heat flux equation for dA:
dq
dA =
U ∆T
We then integrate over the entire heat exchanger (from x = 0 to x = L) to obtain the total area of
required for this heat duty, treating q as the independent (or integration) variable instead of x:
=A ( L ) A=
T q ( L ) qT
dq
=AT ∫ dA
= ∫ U ( q ) ∆T ( q )
=A ( 0 ) 0=q ( 0) 0
q ( L ) = qT
dq
or AT = ∫ U ( q ) ∆T ( q )
(3)
q ( 0) = 0
This serves as the general design equation for heat exchangers. To evaluate this integral, it’s
useful to think of q [= q(x)] as being the independent variable instead of x.
Example: To show how this done, let’s consider a simple example in which the flow of hot and
cold streams is countercurrent and no phase change is occurring in either fluid.
where m is the mass flowrate of the stream (units are x=0 x=x x=L
kg/sec), H is the enthalpy per unit mass (units are
J/kg), the subscripts and primes have the following meaning:
m H ( x ) − H1 = q ( x) = m′ H ′ ( x ) − H1′ (6)
rate of energy rate of energy
gained by lost by
cold stream hot stream
Now, let’s relate each enthalpy to the appropriate temperature. In general, enthalpy changes can
have both latent heat (caused by change in phase) and sensible heat (caused by change in
temperature) contributions.
=
dH c p dT + λ dy
sensible latent
heat heat
where
When the fluid consists of two phases (liquid and vapor), which have different heat capacities
cpL and cpV, we use weighed average of the two heat capacities*
cp = (1-y)cpL + y cpV
where y is the mass fraction of vapor present and 1-y is the mass fraction of liquid present. For
the time being, let’s assume:
Now we only have to worry about contributions from changes in the temperature (sensible
heat):
dH = cp dT
Integrating over the temperature change allows us to evaluate the change in enthalpy
* A vertical line drawn on the left side of a paragraph denotes material which was omitted in class
using using
assumption assumption
Hb #1 Tb #2
∫ dH
= ∫=
c p dT c p (Tb − Ta )
H
a
Ta δT
δH
In general, the specific heat (cp) depends on temperature. But if the temperature change is not
too large — say 10-20° — then the specific heat doesn't change too much. Let's assume we can
neglect these small changes:
Then we can factor cp out of the integral and the integration becomes trivial:
H – H1 = cp(T – T1)
and H ′ − H ′1= c ′p (T ′ − T ′1 )
where cp and c'p are the specific heats of the cold and hot streams. Thus (6) can be rewritten as:
at x=x: mc p T ( x=
) − T1 m′c ′p T ′ ( x )=
− T ′1 q ( x ) (7)
For example,
at x=0: ( )
mc p (T1 − T1 ) = m′c ′p T ′1 − T ′1 = q ( 0 ) = 0
In particular, we are interested in the total rate of heat transfer over the entire length L of the
exchanger:
at x=L: ( )
mc p (T2 − T1 )= m′c ′p T ′2 − T ′1 = q ( L ) ≡ qT
Using (7), I can show that the temperature of either stream [i.e. T(x) or T’(x)] is a linear
function of the amount of heat transferred, q(x):
T
1
T ( x=) T1 + q ( x) T' vs q
∆T2
mc p
T vs q 1
1 ∆T1
T ′ ( x=
) T1′ + q ( x)
mcp
m′c ′p ∆ T vs q
∆T1 ∆T2
Since both temperatures are linear functions of q, q qT
0
their difference ∆T is also a linear function of q. A
straight line has the very nice property that its slope
d ( ∆T ) ∆T2 − ∆T1
=
slope = const
= (8)
dq qT − 0
I can easily evaluate the driving force (∆T) at the two ends of the heat exchanger, and I know the
total heat duty (qT), so I can calculate this slope. The dq and d(∆T) which appear in (8) represent
the incremental rate of heat transfer and the change in driving force ∆T which occurred over a
differential length dx of the heat exchanger(from x to x+dx).
qT
1 dq
AT =
U ∫ ∆T ( q )
0
qT
1 dq d ( ∆T )
AT =
U ∫ d ( ∆T ) ∆T
0
Since the slope is a constant, so is the inverse of the slope dq/d(∆T) which can be factored out of
the integral. Its value is given by the inverse of (8). After substituting this constant value, we
have
∆T ∆T
1 dq 2
d ( ∆T ) 1 qT 2
d ( ∆T ) 1 qT ∆T
=AT = ∫ = ∫ ln ( ∆T ) ∆T2
U d ( ∆T ) ∆T U ∆T2 − ∆T1 ∆T U ∆T2 − ∆T1 1
∆T ∆T
1 1
∆T2
ln ∆T2 − ln ∆T1 =
ln
∆T1
In the above we have also changed the limits of integration to reflect the fact that we are now
integrating with respect to ∆T instead of q: on the left end of the heat exchanger where q(0) = 0,
the ∆T = ∆T1; whereas on the right end of the heat exchanger where q(L) = qT, the ∆T = ∆T2.
This leaves us with
∆T2
ln
q ∆T1
AT = T
U ∆T2 − ∆T1
∆T − ∆T1
Solving for qT: =
qT UAT ∆Tlm where ∆Tlm =2 (9)
∆T
ln 2
∆T1
is called the log-mean temperature difference (LMTD). It represents the average driving force
for heat transfer.
Summary: The general design equation for heat exchangers is (3). (9) is a simplified form which
is valid whenever 1) U(q) can be replaced by a suitable constant and 2) ∆T is a linear function of
q. For most of the design problems in this course, we will use (9) as the design equation.
We have just shown that these two conditions are met by making the following assumptions:
In the homework, you will show that (9) is also correct if the hot and cold streams flow co-
currently; you will also show that δH can also be all latent heat and the LMTD is still the
appropriate average driving force. For example, if you use condensing steam as your hot fluid,
as long as the steam remains saturated:
One last detail: in (1) we said A is the area of the inner pipe. Owing Di
to the finite thickness of the pipe wall, the inner and outer diameter of Do
the inner pipe are different (Di < Do), which implies that the outer
surface Ao of the pipe wall is larger than the inner surface Ai of the pipe
wall. Either area can be used to calculate the heat duty, but a different
heat transfer coefficient must also be used to obtain the same heat duty: Ui and Uo are defined
such that:
qT = U o Ao ∆Tlm = U i Ai ∆Tlm
From the figure, it is clear that Do>Di. This implies that Ao> Ai so that Uo<Ui.
Given a choice between cocurrent or countercurrent flow, which is preferred? It turns out
that countercurrent flow is usually more efficient. To see why, let’s look at an example:
Example: A cold fluid is being heated from 15.8ºC to 62.9ºC while the hot fluid is cooled from
104.6ºC to 71.3ºC. Compare ∆Tlm for co-current flow to that for countercurrent flow.
104.6ºC → 71.3ºC
Solution: For co-current flow, we have the adjacent temperatures
shown at right. This corresponds to the following ∆T’s at the two ends
of the heat exchanger: 15.8ºC → 62.9ºC
8.4°C − 88.8°C
( ∆Tlm )cocurrent =
8.4°C
=
34.1°C
ln
88.8°C
71.3ºC ← 104.6ºC
For counter-current flow, we the adjacent temperatures shown at right.
This corresponds to the following ∆T’s at the two ends of the heat
exchanger: 15.8ºC → 62.9ºC
41.7°C − 55.5°C
( ∆Tlm )countercurrent =
41.7°C
=
48.3°C
ln
55.5°C
Notice that, for exactly the same heat duty, the average driving for countercurrent flow is higher
than for co-current flow. This means that less heat exchanger area is needed if the flow is
countercurrent.
The general design equation (3) can be replaced by the much simpler form of (9), provided 1)
U(x) can be taken to be a constant over the entire exchanger and 2) ∆T is a linear function of q.
In the analysis at the beginning of this chapter, we showed that the latter is a reasonable
approximation if the flow is counter-current and no phase change occurs. In the homework, you
will show that this can be extended to 1) co-current flow with no phase change and 2) isothermal
phase change in one fluid. So you might well ask “when can’t (9) be used?
To give one example, suppose the hot fluid is superheated steam. Over part of the heat
exchanger surface, we have cooling of the hot vapor to its saturation temperature. Over the rest
of the surface, we have condensation of the saturated steam. The result is a nonlinear
relationship between ∆T(x) and q(x).
q(x)/qT
Tbp ( ps)
400 0.5
Hot fluid T
350
Cold fluid T
q(x)/qT
300 0
Tci = 3000 K 0.2 0.4 0.6 0.8 Tco = 400 K
----->----->----->----->----->----->----->----->----->
x/L
Figure 1
Figure 1 above shows the hot and cold temperature profiles and the cumulative heat duty
q(x). Figure 2 below shows the corresponding ∆T(x) and ∆T(q) functions. Notice that neither
function is linear over the entire length of the exchanger. But note that ∆T(q) consists of two
linear segments: one for the region in which superheated vapor is cooling to its saturation
temperature and a second line for the region in which saturated vapor is condensing.
150
deltaT vs q/qT
deltaT vs x/L∆T2
delta T (degC)
100 ∆T1
50
0
0 0.2 0.4 0.6 0.8
q/qT or x/L
Figure 2
Since ∆T(q) is piecewise linear, we might evaluate (3) by breaking the integral into two parts
and taking advantage of the linearity over each region:
qT q1 qT
dq dq dq
=AT ∫=
U ( q ) ∆T ( q ) ∫ U ( q ) ∆T ( q ) + ∫ U ( q ) ∆T ( q )
0 0 q1
q1 q − q1
or =AT + T (10)
U1∆Tlm1 U 2 ∆Tlm 2
where q1, U1 and ∆Tlm1 are the heat duty, overall heat-transfer coefficient and log-mean driving
force for the region of the hxer surface in which saturated vapor is condensing; qT – q1, U2 and
∆Tlm2 are the same quantities evaluated for the other region in which superheated vapor is being
cooled to its saturation temperature.
Example: evaluate the area required for the situation shown in Figure 1 in which 1,000 lb/hr of
water is being heated from 300 K to 400 K by superheated steam at 500 K and 4 atm pressure
using U = 700 btu/ft2-hr-ºF.♣
Hhi = 1223 btu/lb (4 atm, 500K) and Hvap = 1178 btu/lb (4 atm, 418K)
♣ Condensing steam has a different U compared to cooling of a vapor. We should really use two different U’s for
the two different regions but in this example we will use the same U for both regions.
♥ Appendices A.2-9 (saturated steam) and A.2-10 (superheated steam).
By taking the outlet enthalpy to equal the enthalpy of saturated liquid, we are assuming that the
steam leaves the shell side as condensate at its boiling temperature and at the same pressure at
the inlet (4 atm). The flowrate of steam can be calculated from an overall energy balance:
Substituting mc = 1,000 lb/hr, cpc = 1 btu/lb-ºF, Tco = 400 K and Tci = 300 K along with the
values of the enthalpies above, we obtain qT = 52.75 kW and mh = 187 lb/hr. The part of the heat
duty associated with condensing steam is
To calculate the ∆Τ’s appearing in (10), we need to evaluate the hot and cold temperatures at
the point in the heat exchanger where the superheated steam becomes saturated steam. We will
denote the temperatures at this point Th1 and Tc1. Th1 = 417.6 K is the temperature of saturated
steam at the inlet pressure (4 atm). The other temperature can be calculated from
Comment #1: If instead, we blindly applied (9) using the two end-point ∆T’s, we would have
calculated
which considerably overestimates the average driving force, leading to an large underestimate of
the area required: 1.32 ft2 instead of 2.51 ft2 (almost a factor of 2 error).
Comment #2: A much better approximation would be to ignore any superheating or subcooling
of the steam: just assume the steam enters and leaves at the saturation temperature.
This leads to ∆Tlm = (117.6 K, 17.6 K)lm = 52.7 K which lies in between the two exact ∆T’s. In
the next chapter (Evaporation), this is exactly the approximation we shall make in preliminary
designs.
RESISTANCES IN SERIES
=
qT UAT ∆Tlm (9)
We have just discussed how to calculate ∆T in this equation. Now we will turn our attention to
calculating U. Heat transfer coefficients are usually estimated from empirical correlations of
experimental data. However the correlations usually do not give U directly; instead, the
correlations are usually for the “film coefficients” hi and ho, which are related to U by the
following formula:
1 1 D x D 1
= o + w o +
U o hi Di kw Dlm ho
We will now review the theory of heat transfer which leads to this result.
Consider heat transfer through the inner pipe of a double-pipe heat exchanger. Let’s take a
thin axial slice (of length dL) out of the heat exchanger and examine the temperature profile
along a radius to the pipe. The labelled temperatures in Figure 3 below assume (somewhat
arbitrarily) that the hot fluid flows in the inner pipe while the cold fluid flows in the annular
space between the two pipes.
Figure 3
One can identify three regions in which the temperature varies sharply. These correspond to the
three resistances to heat transfer: 1) the pipe wall itself, 2) the fluid film inside the pipe and 3)
the fluid film outside the pipe.
Ui and Uo are called overall heat-transfer coefficients because when each is multiplied by
the overall driving force and the appropriate area, they give the heat flux: ♣
definition of U=
i: dq U i dAi (Th − Tc ) (11)
πDi dL
definition of=
Uo: dq U o dAo (Th − Tc ) (12)
πDo dL
♣ The temperature measured at the center of the pipe (T in Figure 3) differs slightly from the “mixing cup”
i
temperature Th which instead is used to characterize the average temperature of the fluid inside the pipe. Similarly
the temperature measured at the outer edge of the annulus (To in Figure 3) differs slightly from the mixing cup
temperature Tc of the cold fluid. The mixing-cup temperature is what would be measured if a sample of the fluid
were collected in a cup and allowed to come to a uniform temperature without any addition or loss of heat. The
mixing cup temperature is also what is used in evaluating the enthalpy of a stream for purpose of an energy balance.
The temperature drop across the fluid film on the inside of the pipe wall is Th – Twh. The
film heat-transfer coefficient hi is defined as the proportionality constant between heat flux
dq/dAi and this driving force:
Since we are calculating the flux through the inner film, we use the area dAi of the inner surface
of the pipe wall. Similarly, the film coefficient ho is defined as the proportionality constant
between heat flux dq/dAo and the driving force across the outer film:
Film coefficients and their definitions (14) and (15) are used to describe convective heat-
transfer rates across fluid films undergoing turbulent fluid flow. By contrast, heat transfer
through the pipe wall occurs by pure conduction. Nonetheless, the relationship between rate of
heat transfer, driving force and heat transfer area can be written in a similar form. For steady-
state conduction through an annulus, the relationship is:
kw
=
conduction through pipe wall: dq dAlm (Twh − Twc ) (16)
xw
where kw is the thermal conductivity of the pipe wall and xw is the wall thickness
Do − Di
xw = ro − ri =
2
dAo − dAi π ( Do − Di ) dL
and where dAlm = = = πDlm dL
ln
dAo Do
dAi
ln
Di
is the log-mean of dAo and dAi. A relationship similar to (16) (including the log-mean area) was
probably derived in your heat-and-mass transfer class.*
Notice that the same dq is used as the rate of heat transfer in (11) through (16). This
sameness is a consequence of assuming steady-state heat transfer. At steady state, there can be
no accumulation of heat at any location.
The two areas dAo and dAi do not differ by very much so neither do the values Uo and Ui.
The choice of which area to use is arbitrary. We will proceed using the outer area. Then (12)
can be written as
dq
= Th − Tc
U o dAo
dq
=Th − Thw + Thw − Tcw + Tcw − Tc
U o dAo
dq dq dq
hi dAi kw ho dAo
dAlm
xw
1 1 D x D 1
= o + w o + (17)
U o hi Di kw Dlm ho
after substituting the expressions relating the dA’s to D’s and dL. Starting with (11) instead of
(12), the analogous expression for Ui is
1 1 xw Di 1 Di
= + +
U i hi kw Dlm ho Do
If we divide (17) by the outer area (Ao = πDoL) we obtain what Geankoplis calls the total
resistance to heat transfer:
1 1 xw 1 1
= + + =
U o Ao hi Ai kw Alm ho Ao U i Ai
total resistance of resistance of resistance of total
resistance inner film pipe wall outer film resistance
To see why 1/UA might be called the resistance, recall again the design equation (1):
1
=
q U A ∆T or ∆
T=q
I
V UA
R
The heat duty q is like the electric current I in Ohm’s law (V = IR) and the driving force ∆T is
like the voltage drop V across a resistor having resistance 1/UA. In this expression the total
resistance is either 1/UiAi or 1/UoAo. Since the same current is passing through all three
resistors, the resistors must be arranged in series.
Fouling Factors
The usefulness of analogies is that we can extend our knowledge in one area using
knowledge in other. For example, if there are additional layers through which heat must be
conducted, then we can just add their resistances to get the total. A situation involving additional
layers is:
fouling — deposition of dirt, scale or any solid deposit on one or both sides of the tube
Fouling generally reduces heat transfer by imposing an additional resistance. It's effect is
accounted for by use of
k
fouling factors: hd = d
xd
which are calculated from simple conduction theory, assuming the fouling layer is planar film of
thickness xd; compare with coefficient in (16).
There may be a fouling factor for both the inside and the outside of the pipe. Adding these
resistances to the others:
1 1 1 x D 1 Do 1 Do
=+ + w o + +
U o ho hdo kw Dlm hdi Di hi Di
fouling fouling
outside pipe inside pipe
By means of the electrical analogy, we are able to guess the form of the correction without
having to go through the analysis.
1. What is the geometry? (e.g flow inside pipe, flow across bank of tubes, flow around sphere)
2. Is the flow turbulent or laminar?
3. Does a phase change occur?
Our book presents a number of common correlations which cover the most common situations:
For flow inside long tubes, the criteria for turbulent flow is simple: just calculate the
Reynolds number:
ρ v Di GDi
=
N Re =
µ µ
Q
The latter quantity is defined as v =
Si
Another term which is sometimes used to calculate Reynolds number is the mass velocity:
ρQ m
G≡ρ =
v = (units of mass per area per time)
Si Si
The mass velocity is often used in problems involving gas flow where both density and velocity
depend on pressure, but their product is independent of pressure (since it represents mass
♣ You should make a special effort to read these pages before working the HWK problems, since I don’t have time
in class to present all the details.
flowrate and at steady state, a mass balance requires the mass flowrate to be the same at all axial
positions). However Reynolds number is calculated, the criteria for turbulent flow is as follows:
For turbulent heat transfer (NRe > 6000) of high Prandtl number (0.7 < NPr < 16,000) fluids
in long (L/D > 60) pipes, the individual film coefficient inside a pipe can be estimated using the
Sieder-Tate equation [eq. (4.5-8) in G]:♦
0.8 1/ 3 0.14
hi Di GDi c pi µi µi
= 0.023 (18)
ki µi
ki µiw
N Nu N Re N Pr
ρQ m
and where G≡ρ =
v =
Si Si
is called the mass velocity (units of mass per area per time). Unless specifically indicated, all
fluid properties are evaluated at the mixing-cup temperature of the fluid:
The subscript “iw” indicates the property is to evaluated at the inner wall temperature:
Experimentally, one finds that heat transfer coefficients can be different for heating and cooling
over the same temperature range. The reason is that the viscosity at the wall is very important in
determining transport and µw can be very different from µ. This is accounted for in the µ/µw
factor.
♦ The coefficient of 0.023 in this equation is the value which appears in McCabe, Smith & Harriott. Geankoplis
gives a value of 0.027 instead. In homeworks and exam, you can use either value. The Notes use 0.023 because
most of the homework problems using this equation were taken from MS&H.
When the inlet and outlet T’s of a heat-transfer fluid differ, at what T do you use to evaluate
the properties? Strictly speaking, you should use the local value of T(x), which then leads to
different U at each x. In turn, (3) must be used as the design equation. This is significantly more
work than using (9).
In most cases, U does not vary significantly with local temperature. Then a good
approximation is to treat U as a constant evaluated using the arithmetic mean values for T and
Tw:
T + To T + Two
T = i and Tw = wi
2 2
where the subscript “i” refers to the inlet and “o” refers to the outlet.* Since the wall temperature
Tw depends on the heat transfer coefficient, the determination of h is often a trial-and-error (see
Example 4.5-2).
Water is often used in heat-transfer equipment. The following simplified equation contains
the variation in physical properties with temperature and can be used in the temperature range 4
< T < 105°C:
v 0.8
for liquid water: hL =1429 (1 + 0.0146T °C ) (19)
D 0.2
The two coefficients in this equation have some weird units (involving fractional exponents).
Rather than giving those weird units, we just point out that the constants were chosen to work
with SI units: the equation yields hL in W/(m2-°K) if v is expressed in m/s and D is expressed in
meters. Similar simplified correlations exist for air at 1 atm or for organic liquids:
v 0.8
for air at 1 atm: hL = 3.52
D 0.2
v 0.8
for organic liquids: hL = 423
D 0.2
* This represents a change in notation: previously “i” stood for inner surface while “o” stood for outer surface.
pipes (or ducts) having noncircular shape: for example, the annular space between the inner and
outer pipes of a double-pipe heat exchanger (see blue region in figure at right).
The main question is: what is used in (18) for Di? Note that an annulus has two diameters,
which are labelled DA and DB in the figure above. The answer is to use the hydraulic diameter
which is defined as
S
DH ≡ 4
P
where S is the cross-sectional area of the duct and P is its wetted perimeter. For a circular pipe
of inside diameter Di, the area and wetted perimeter are
π 2
Scircle = Di and Pcircle = πDi
4
Note that the hydraulic diameter is just the usual diameter for a circular pipe. For the annulus in
the figure above, the area and wetter perimeter are
=
Sannulus
π 2
4
( 2
DB − DA ) and π ( DB + DA )
Pannulus =
DB2 − DA2
so the hydraulic diameter is D=
H = DB − DA
DB + DA
TUBE PITCH
The bundle of tubes in a shell-and-tube heat exchanger can be stacked in one of two different
ways
The use of triangular pitch allows the tubes to be more tightly packed -- more tubes and therefore
more area per unit volume of shell. This makes the shell cheaper. On the other hand, the square
pitch has the advantage that it is easier to clean. You can get a long brush in between two rows
or between two columns. Of course, the shell has to be pulled off of the tube bundle before you
can get at it in either case.
BAFFLES
To some extent, it is advantageous to reduce the cross-sectional area available for flow. For
the same volumetric flowrate, this will increase the mass flux or the mass velocity of the fluid
flow, which will increase ho. Of course, this also increases the pressure drop and the pumping
cost.
Ds = diameter of shell
A number of these plates is placed along the length of the exchanger with the location of the
baffle window alternating between the top and bottom of the plate.
A second advantage of these plates is that they prevent the formation of large stagnant regions
(having low heat transfer rates) which usually accompany “channelling” of the fluid in a
direction path connecting inlet to outlet. The spacing between the baffle plates is called the
baffle pitch:
P = baffle pitch
Besides reducing the area available for flow on the shell side, the baffles also force the fluid
across all of the tubes, which normally will enhance the heat transfer coefficient.
An approximate but generally useful correlation for estimating the shell-side coefficient in a
shell-and-tube heat exchanger is the Donohue equation:
♠ A PDF file containing a scan of these pages can be found under “Handouts” in Blackboard.
and coefficients are different. Another difference is that the mass velocity Ge is a weighted
average
m Kg
=
Gb =[ ] (22)
Sb m2s
m
Gc = (23)
Sc
where m is the mass flow rate (Kg/s), Sb is the open area available for flow of the shell-side
fluid through the baffle window and Sc is the open area available for crossflow perpendicular to
the bank of tubes at the widest point in the shell. The formulas for calculating these two areas
are given in the Appendix starting on page 37.
MULTIPASS CONSTRUCTION
Just as we can improve the heat transfer coefficient on the shell side by reducing the baffle
pitch, we can improve the heat transfer coefficient on the tube side by restricting the flow to half
of the tubes. Both have the effect of reducing the area available for flow, thereby increasing the
fluid velocity and the heat transfer coefficient.
Suppose we partition the manifold in such a way that the inlet fluid can only enter half of the
tubes. The other half of the tubes are used for return flow. We say that such an exchanger has
two tube passes.
♣ An equal sign enclosed in square brackets (i.e. [=]) should be read “has units of …”
Although two tube passes improves hi, one of the passes will be co-current. A 2-2 exchanger
improves both hi and ho, and in addition, the flow is more nearly countercurrent. We could also
partition the shell side with a horizontal baffle in addition to the vertical baffles, producing two
shell passes.
• 1-2 exchanger
• 2-4 exchanger
where the first number is the number of shell passes and the second number is the number of
tube passes.
An even number of tube passes is desirable because then input and output of the tube-side
fluid occur on the same end of the exchanger. This allows the other end of the exchanger to
"float" in response to thermal expansion of the metal, which can otherwise be a serious problem
when the shell changes temperature during startup or shutdown.
qT = Uo Ao ∆T (24)
For double-pipe heat exchangers, the appropriate average driving force ∆T is usually the log-
mean ∆Tlm of the two ∆T evaluated at either end of the heat exchanger. For multipass shell-and-
tube heat exchangers, the flow is neither counter-current nor co-current. Instead, the flow is
some combination of counter-current, co-current and cross flow. As a consequence, the average
driving force is no longer the log mean. Instead the appropriate average driving force is called
the true mean temperature difference (MTD) and is somewhat less than the log-mean:
For pure crossflow to a bank of tubes, or for a mixture of crossflow and counter-current flow, the
factor is less than unity:
1-2 shell-and-tube
2-4 shell-and-tube
cross flow to a bank of tubes
T −T T − Tho
Y = co ci and Z = hi Tho Thi
Thi − Tci Tco − Tci
where the subscript “c” refers to the cold fluid Tci Tco
while the subscript “h” refers to the hot fluid and
the subscript “i” refers to the inlet while “o”
refers to the outlet.* Using this same notation,
the log-mean ∆T is calculated from
This is the log-mean ∆T for counter-current flow: it is used for all shell-and-tube heat
exchangers, regardless of type (1-1, 1-2 or 2-4).
* This represents a change in notation: previously “i” stood for inner surface while “o” stood for outer surface.
so that Z = 0 and FT = 1
for condensing steam as the hot fluid. Similarly if we have boiling water as the “cold” fluid, its
temperature doesn’t change (i.e. Tci = Tco = Tbp) and then we have
Y = 0 and FT = 1
Although the flow in the exchanger is complex, we calculate ∆Tlm in (25) as if the flow were
counter-current. The deviation caused by not being counter-current is included in the values of
the empirical correction factor FT.
An Example:*
A steel mill needs to cool 3,000 gal/hr of a dilute acid from 200°F to 130°F using 5,000 gal/hr of
80°F river water. A shell-and-tube heat exchanger is available (see description below).
Determine if this heat exchanger is sufficient for this task.
type: 1-2
tubes: 158 1-inch, 14BWG stainless tubes, 16 ft long, stacked in 1¼ inch square pitch
baffles: 25% with 6 inch pitch
shell: 20 inch inside diameter, mild steel
Solution: The strategy is to compare the required heat duty to the heat duty which can be
expected under these conditions. If the former is larger (smaller), the heat exchanger is not (is)
capable of performing this task.
First we calculate the required heat duty: the rate at which heat must be removed from the
dilute acid is
In making this calculation, we have assumed that the properties (density and heat capacity) of
dilute acid are essentially those of pure water. This is to be compared with the available heat
duty, which is calculated from the design equation:
True Mean ∆T. One of the four temperatures is not given: we have to determine the outlet
temperature of the cold stream. The rate of energy lost by the hot stream is the rate of energy
gained by the cold stream
qT m h c ph (Thi − T
= =ho ) m
c c pc (Tco − Tci ) (21)
Assuming the physical properties of dilute acid and river water are the same as pure water, we
can take cph = cpc so that heat capacity can be cancelled out of the second equation above,
leaving
Tco = 122°F♦
Next, we calculate the two quantities needed to evaluate the true-mean driving force:
♦ Such a high temperature in the river could be harmful to fish and other wildlife. Engineers need to be mindful of
the impact of their designs on the environment as well as on the health and safety of people in the vicinity. Faced
with this result, the engineer should suggest altering the design (e.g. increase the flowrate of cooling water to obtain
a lower increase in temperature — most states allow no more than a 10ºF increase). For the purpose of this
academic example, we will ignore this objection and proceed through the analysis.
0.86
0.35
FT = 0.86
Although the flow in the exchanger is complex, we calculate ∆Tlm as if the flow were counter-
current:
∆T − ∆T2
∆Tlm = 1 =63°F
∆T1
ln
∆T2
Overall heat transfer coefficient Uo. We will obtain the overall heat transfer coefficient by first
estimating the two film coefficients and then using the resistances in series formula (17):
1 1 D x D 1
= o + w o +
U o hi Di kw Dlm ho
The shell-side coefficient (ho) can be estimated using Donohue’s correlation (20):
Since acid is corrosive, we would be well-advised to place it inside the stainless tubes, leaving
river water on the shell side (shell is mild steel). With this choice, the mass flowrate on the shell
side is
gal lb 4 lb
m 5000
= 8.34 = 4.17 × 10
hr gal hr
The mass flux Ge appearing in Donohue’s correlation (20) is the geometric mean [defined by
(21)] of Gb and Gc defined by (22) and (23). To calculate these last two mass fluxes, we need to
divide the mass flowrate m by either of two different areas Sb and Sc whose formulas are
derived in the Appendix (starting on page 37).
πD 2 πD 2
= Sb fb s − Nb o (28)
4 4
area of
blue region
fb = 0.1955
Let’s assume that the fraction of tubes which pass through the baffle window equals the fraction
of the shell cross-section represented by the baffle window. In other words:
where 158 was given as the total number of tubes in the shell. Now everything in (28) is known:
πD 2 πD 2
Sb = fb s − Nb o =0.258 ft 2
4 4
♦ A 25% baffle plate has a window whose height is 25% of the diameter of the shell (as shown in the figure above).
m lb
G=
b = 1.62 × 105
Sb 2
ft -hr
For crossflow, the maximum available area for flow can be calculated from (29):
D
=Sc PDs 1 − o
p
p = tube pitch
P = baffle pitch
using P = 6 inch, Ds = 20 inch, Do = 1 inch (see Appendix A.5-2), and p = 1.25 inch. The cross-
flow area is
D
Sc = PDs 1 − o = 0.167 ft 2
p
m lb
G=
c = 2.5 × 105
Sc ft 2 -hr
lb
=
Ge =
Gb Gc 2.01 × 105
ft 2 -hr
Before we can estimate ho from (20), we need some fluid properties. Since we have no idea
what the wall temperature is, let’s first get a rough estimate of ho by setting µ/µw = 1. The other
properties (besides µw) are evaluated at the average temperature of the shell-side fluid [T =
0.5(80°F + 122°F) = 101°F (38.3°C)]. Using Appendix A.2-11:
Since the thermal properties of dilute acid are very nearly those of pure water, we will estimate
the inside film coefficient using (19):
v 0.8
hL =1429 (1 + 0.0146T °C )
D 0.2
The tubes are one-inch 14 BWG. Properties can be found in table A.5-2:
Since there are two tube passes, the entering acid is distributed among half of the 158 tubes. The
average velocity in any one tube is then
gal ft 3
3000
hr 7.48 gal ft m
=v = 0.372 = 0.113
158
2
(
0.00378 ft 2 )
s s
(19) can then be evaluated at the average acid temperature: (130+200)/2 = 165°F = 74ºC.
hi = 1125 W/m2-°C
The overall heat transfer coefficient can be evaluated from the sum-of-resistances formula (17).
We just a thermal conductivity. The tubes are constructed of stainless steel. Table A.3-16 gives
thermal conductivities of a number of metals. For 304 stainless at 100°C:
kw = 16.3 W/(m-°K)
Total Heat Transfer Area. The heat duty which can be supplied by the available heat exchanger
is given by (24). The outside area of all 158 tubes (L = 16 ft) is
The required heat duty was evaluated at the beginning as 510 kW. Since the available heat duty
exceeds what’s needed, the available heat exchanger is capable of performing the job. Had the
two heat duties been more nearly equal, we would have had to go back and include the correction
µ/µw in both hi and ho. This is a trial-and-error.
In case you are curious as to how these formulas we obtained, I have put their derivation in
this appendix.
widest
row of
tubes
P
p
Ds Do
Sc is the area open to fluid flow across the widest row of tubes in the shell and between two
adjacent baffle plates. First the total area (ignoring tubes) is
where P is the baffle pitch. To calculate the fraction of PDs which is open. let's subdivide Ds
into repeating unit cells:
Multiplying this fraction by the total area gives the area open:
D
=Sc PDs 1 − o (29)
p
Of course, we have ignored some end effects which could be significant if the number of tubes is
not large.
Likewise the open area of a baffle window is the total area of the window minus the area
blocked off by tubes:
πD 2 πD 2
=Sb fb s − Nb o (30)
4 4
To get this value of fb, we have to do a little geometry. First, recall that the area of a circle of
radius r is πr2.
Acircle = πr 2
= 1=
Atriangle bh 1 xy
2 2
For a 25% baffle, the height of the segment, r-x, will be 25% of the diameter -- or 50% of the
radius:
x π 3
= 0.5 ⇒ θ= cos −1 ( 0.5 )= 60°= ⇒ y= r 2 − x 2= r
r 3 2
1 r 3 3 2
=
Atriangle = r = r 0.217 r 2
2 2 2 8
π
Asec tor =θr 2 = r 2
3
π 2 3 2 π 3 2
Asegment = r − 2 r =−
0.61418r 2
r =
3 8 3 4
Now fb is the fraction of the total area of the circle which is represented by this segment:
π 3 2
− r
Asegment 3 4 0.61418
=fb = = = 0.1955 (31)
Acircle πr 2 3.1416
Chapter 5. Evaporation
evaporation -- concentrating a solution of a nonvolatile solute by boiling away the solvent
(usually water)
Examples: vapor
Evaporation is not the only unit operation that involves boiling. Some others that involve boiling
include:
• distillation - involves two or more volatile components (evaporation has just one)
• drying - product is a solid (in evaporation, product is usually a more concentrated liquid
solution)
• crystallization - product is a slurry of crystals precipitated from solution (not a simple
solution)
Usually these differences are sufficient to require different equipment to perform the operation.
Part #1: HX (for adding the latent heat of vaporization to the liquid feed). Usually HX is a
shell-and-tube with condensing steam as the hot fluid on the shell side (because steam
is free of minerals which form scale).
Part #2: vapor space or vapor head (larger chamber in which liquid, entrained in the vapor as
droplets or foam, can be separated from the vapor, usually by impinging the two-phase
mixture onto a plate.)
You also need something to move the liquid through the heat exchanger. Sometimes a pump is
used and sometimes it's just left to gravity, either in the form of a falling film on the outside of
the HX tubes or the entrainment of liquid by rising bubbles.
• pump
• gravity - falling film
• gravity - entrainment by rising bubbles
steam flow is allowed to bleed out of the shell. This is done to allow air (which is often present
with the steam) to escape. Otherwise non-condensable gases like air will accumulate in the shell,
causing the pressure to rise and decreasing the steam flowrate.
Long-tube vertical evaporators are especially effective in concentrating liquids that tend to
foam. Foam is broken when the two-phase mixture hits the baffle plate.
Forced-Circulation Evaporator
EVAPORATOR PERFORMANCE
Regardless of which type of evaporator we use, there are two standard measures of performance:
economy – lb water evaporated per lb steam used (usually somewhat less than unity)
p = 4 inHg
If the feed is 55,000 lb/hr at 70°F, what area of heat exchanger is required and what flowrate of
15 psig steam?
Find: m s , A
Solution: First we do a mass balance to determine the flowrates. A mass balance on the solids
yields:
is the flowrate of the liquor. The corresponding flowrate of the vapor must be (from an total
mass balance):
The heat duty of the HXer is just the heat lost by the condensing stream, or the heat gained by
the feed. In terms of enthalpy per mass of the streams, the heat gained by the feed is:
For high-molecular-weight substances like orange juice and other “organic colloids”♥ we can
usually neglect boiling-point elevation and heat of dilution. Then the enthalpies of the feed and
the liquor are essentially the same as pure water under those conditions. The enthalpy of liquid
water and water vapor can be read from the Steam Tables found in Appendices A.2-9 (saturated
steam) and A.2-10 (superheated steam).♦
But what is the temperature of the liquor? Well, since the liquor is in equilibrium with the vapor,
their temperature and pressure should be the same. Since the vapor is just saturated steam at 4
inHg, we look up its properties in the steam table:
♥ “Organic colloids” are solutions of high molecular-weight polymers or dispersions of fine particles. Then the
molar concentration of these solutions is negligibly small and their properties are close to that of pure water: see
example on page 48. An example of an organic colloid is orange juice.
♦ For consistency with later examples in which heat of dilution effects are important, we read all the enthalpies
from the steam tables, which have a reference state (where H is defined to be zero btu/lb) corresponding to saturated
liquid water at its freezing point. In Example 8.4-1, Geankoplis works a very similar problem in which all
enthalpies are calculated using the conditions of the liquor as the reference state. Then (by definition) H in (32) is
zero; the enthalpy of the feed Hf differs from that of the liquor just in sensible heat, which can be estimated as (8.4-
9); and the enthalpy of the vapor Hv differs from that of the liquor just in latent heat (i.e. Hv is just the heat of
vaporization at the T,p of the vapor space). See the discussion on “steam economy” on page 45. Try re-working
Example 8.4-1 evaluating the enthalpies appearing in (8.4-7) using the steam tables (as above): you should get very
nearly the same answer.
♠ Strictly speaking, the enthalpies tabulated in A.2-9 are for saturated liquid, meaning that the pressure on the liquid
equals the vapor pressure. Here the pressure on the feed is probably 1 atm (or 14.7 psia) rather than 0.3622 psia (the
vapor pressure). Pressure has a very weak effect on the enthalpy of liquid water and the effect of pressure is usually
neglected. From thermodynamics (see 4.2-24 on p169 of Sandler):
∂H ∂V
∂P = V − T ∂T
T P
which identically vanishes for an ideal gas ( PV = RT ). For most liquids like water, the second term is negligible
compared to the first. Then the effect of a change on pressure equal to 1atm on the enthapy of liquid water can be
estimated as
cm3 btu
∆p V ≈ (1 atm ) 1
∆H T = =0.044
g lb m
which is equivalent to a change in temperature of only 0.044 °F. As long as we are dealing with much larger
sensible temperature changes, we can neglect the effect of pressure on the enthalpy of liquid water.
Now, using the design equation, we can calculate the required area:
q = UA∆Tm (33)
In a typical evaporator problem, steam is condensing on the shell side while water is evaporating
on the tube side. When a phase change occurs with a single component (water in this case), the
temperature stays at the boiling point. Thus the driving force ∆T is practically constant along the
length of the heat exchanger.
it’s temperature is above the boiling point. Until the Steam leaves
temperature drops to the boiling point, condensation subcooled
does not occur. Similarly, so much heat might be T
removed from the steam, that it completely condenses Feed solution
and is then subcooled below the boiling point. Finally, enters subcooled
the feed solution has to be heated up to its boiling point 0 L
before evaporation occurs.♥ The temperature profile Position in Hxer
at right illustrates the most general situation.
♣ In working the homework assignments, you will have frequent need to perform this interpolation. A Mathcad
document has been prepared to assist you: Steamtab.mcdx. It can be found at our Blackboard website: see the
“Student Resources” folder under “Handouts”.
♥ The heat required to preheat the feed from 70°F (the feed conditions) to 125°F (the boiling point at the pressure in
the vapor space) is
lb btu btu btu
55000 93.1 − 38.09 =0.30 × 107
hr lb lb hr
This is only 6.3% of the heat duty q. We are willing to tolerate this size of error in our preliminary design.
♠ For the alternative to making this assumption, see the example beginning on page 14.
∆Tm ≈ Ts - T
“Proof”: Assuming that the sensible heat associated 1) with superheating the steam, 2) with sub-cooling the
condensate and 3) with pre-heating the feed up to its boiling point are all negligible compared with the latent heat of
vaporization of either stream, then the fraction of the heat exchanger area experiencing a driving force different from
Ts – T will be negligible; thus the entire heat exchanger can be approximated as having Ts – T as the driving force: in
particular,
∆T − ∆T2
∆T1 = ∆T2 = Ts – T and ∆Tlm =1 = Ts – T
∆T
ln 1
∆T2
When ∆T1 and ∆T2 have identical values, it seems reasonable (almost obvious) that their mean should have the same
value; however the expression above for log-mean is mathematically indeterminate: ∆Tlm = 0/0. This indeterminacy
can be resolved if we consider the asymptotic limit of this expression as the two quantities being averaged approach
each other. Define
∆T1 = ∆T2 (1 + ε )
∆T1 1 1 1 ε→0 1
and ln =ln (1 + ε ) =ε − ε 2 + ε3 − ε 4 + → ε − ε 2
∆T2 2 3 4 2
where second equation represents a Taylor series expansion. Finally, we take the limit as ε→0:
When the two numbers being averaged are nearly the same, the log-mean and the simple mean yield the same result.
Back to our simple evaporator problem: Knowing the pressure of the steam, we can find its
temperature in the steam tables:
The heat of vaporization λs is the difference between the enthalpy of saturated vapor Hs and the
enthalpy of the condensate Hc leaving (assumed here to be saturated liquid):
The driving force for heat transfer is the difference in temperature of saturated steam and
saturated vapor:
Btu Btu
4.77 × 107
= 500 2 ( A )(125°F )
h h − ft − °F
or A = 768 ft2
The amount of steam required is just the heat duty divided by the latent heat of the steam:
Btu
4.77 × 107
q h= 50,800 lb
s =
m=
λs Btu hr
945.6
lb
Recall that our two main performance indicators were “capacity” and “economy”. Capacity is
just the rate of evaporation, which will be the flowrate of the vapor. For this example,
lb
capacity = m f − m
=44,000
hr
The economy is the mass of vapor produced per mass of steam consumed:
lb
m f − m 44, 000 hr
= =
economy = 0.87
m s lb
50,800
hr
Steam Economy
Why is the economy less than one? The answer becomes clear, if (instead of looking up the
enthalpies of water and its vapor in a table) we estimate them as sensible and latent heat. First,
let’s select the liquor as the reference state for the calculations of the enthalpies. Then
The vapor is at the same temperature as the liquor so there are no sensible heat changes, but the
vapor has a much higher enthalpy owing to its heat of vaporization:
Hv = λ
Since the feed is also a liquid (like the liquor), its enthalpy differs from the liquor only in
sensible heat:
Hf = cp(Tf - T)
( m f ) ( )
− m λ − m f c p T f − T = m s λ s
or
m f − m λ s m f c p T − T f
= −
( ) (34)
m s λ m s
λ
<1
Generally the heat of vaporization of water decreases with increases in temperature. This is a
consequence of the specific heat of vapor being less than that of liquid water. To show this,
recall that the heat of vaporization λ is the difference between the enthalpies of saturated vapor
and liquid:
=
λ HV − H L
The enthalpy derivatives are heat capacities, which have been substituted in the second equation
above. Now the heat capacity of steam is about half of that for liquid water, which leads to the
numerical values in the third and fourth equations above. In the above calculation, we have held
pressure constant whereas the heat of vaporization is usually evaluated at the vapor pressure of
the liquid which increases significantly with temperature. However the effect of pressure on
enthalpy is generally negligible (see footnote ♠ on page 43). Thus we have shown that
Alternatively, We can use the steam table A.2-9 on p964 of Geankoplis to directly evaluate
the changes in λ with temperature.
The drop in latent heat in this example is about 62.5 BTU/lbm which is close to 58 BTU/lbm
obtained by multiplying cpV-cpL by 100ºF (or 58 BTU/lbm).
Since the steam must be hotter than the water being boiled, this means that λs < λ. So the
first term in (34) is already less than one. The second term decreases it further; this second term
represents the energy required to preheat the feed up to the temperature of the liquor where
boiling occurs. So a pound of steam boils less than one pound of water because some of the
energy of the steam is required to preheat the feed. However, even if the feed were at the same
temperature as the liquor, a pound of steam would boil less than one pound of water because the
latent heat of the hot steam is less than the latent heat of the water being boiled.
∆Tm ≈ const. = Ts - T
where T is the boiling point of the liquor at the pressure in the vapor space. This is approximate
when there is significant BPE because T will vary along the length of the heat exchanger tube as
the concentration of solute also changes.
• Evaluation of enthalpies with negligible BPE and negligible heat of dilution (example on
page 42 and also Example 8.4-1). Assuming the solvent is water, we can read the enthalpies
of all three streams from the steam tables.
♣ Since enthalpies are being taken from more than one source, we should take care that the reference state is the
same for all tables and charts used in a given problem. The steam tables in Geankoplis use saturated liquid water at
32°F and 1 atm as the reference state. The enthalpy chart for NaOH solutions in Geankoplis also use 0 wt% NaOH
at 32°F and 1 atm as the reference state.
COMPLICATIONS
Boiling-Point Elevation
If the solution is more concentrated (in a mole-fraction sense), complications arise because
the enthalpy and partial pressure of a liquid more generally depend on the solute concentration as
well as the temperature. Whenever a solute is added to water, the partial pressure of the water is
reduced. Recall Raoult's law from thermodynamics.
Note that adding any solute with water causes xw<1 so that
pw < pwo
In general, this implies that I must heat the solution up to a higher temperature (i.e. higher than
the b.p. of pure solvent) before it will boil. This increase in boiling point due to addition of
solute is called
In effect, the boiling point of a solution depends, not only on the total pressure, but also on the
concentration of solute:
Tbp = Tbp(p, c)
This effect can be quite significant for concentrated solutions of inorganic salts. For
example, for at a pressure of 1 atm, the boiling point of pure water is
If the pressure is other than 1 atm, the temperature can be read from the steam tables. Given this
boiling point for pure water, we can look up the boiling point of solutions of NaOH from a
Dühring chart (see Fig. 8.4-2 of Geankoplis):
Notice that the lines (which correspond to a particular mass fraction) are nearly parallel. If they
were perfectly straight and parallel, then BPE would be independent of pressure. From the
Dühring chart, we can read the boiling point of a concentrated NaOH solution:
Note that the boiling point of the solution is significantly greater than the boiling point of pure
solvent. BPE is a colligative property, which means
Example: Use Raoult’s law to estimate the partial pressure of water at 100°C of a) a 50wt%
NaOH solution and b) a 50wt% suspension of “organic colloids” (use MWaverage = 1000 g/mol).
Estimate the BPE for each.
Solution: To calculate the mole fraction of water in the solution, take 1 gram of solution as a
basis. This 1 gram mixture consists of 0.5 g of water (MW = 18) and 0.5 g of NaOH (MW =
40). We calculate the corresponding number of moles of each and the mole fraction:
NaOH (MW=40) Organic Colloids
(MW=1000)
0.5g 0.5g
g g
18 18
mol = 0.690 mol = 0.982
xw
0.5g 0.5g 0.5g 0.5g
+ +
g g g g
18 40 18 1000
mol mol mol mol
At 100°C (the b.p. of pure water), the vapor pressure of pure water is 1 atm (= pwo). Using
Raoult’s law, the partial pressure is estimated as
To cause the solution to boil, the partial pressure must be raised to 1 atm. For the same mole
fraction of NaOH, this means raising the vapor pressure to
1 atm 1 atm
pwo = 1.450 atm = 1.018 atm
0.690 0.982
Using the steam tables (Appendix A.2-9), we find that the temperature has to be
T 110.0°C 100.8°C
Of course, this is only an estimate. Raoult’s law is only asymptotically correct as xw→1. The
actual BPE for NaOH is 39°C.
Although BPE is significant for a low-molecular-weight solute like NaOH, the same 50wt%
concentration of an "organic colloid" (e.g. orange juice) may have a negligible BPE because the
molecular weight is much higher. Thus the same wt% gives a much smaller molar concentration
of a high-molecular-weight solute and therefore a much smaller BPE.
Heat of Dilution
If you dilute concentrated sulfuric acid by adding water, large quantities of heat can be
released. If the mixing is adiabatic and both fluids are initially at room temperature, the final
temperature after mixing can by high enough to causing boiling. You were probably warned to
use great caution when diluting sulfuric acid in chemistry lab. The problem is that the enthalpy
of a solution depends not only on its temperature, but also on its concentration:
H = H(c,T)
which means you then need more than just a heat capacity to calculate sensible heat changes.
Figure 6 gives the enthalpy-concentration diagram for NaOH solutions.
Example: Suppose we dilute 1-lb of 50wt% NaOH solution at 70°F with 1-lb of pure water, also
at 70°F. What is the final temperature if the mixing is adiabatic?
Multiplying each by the corresponding weight and adding, we obtain the final enthalpy of the
mixture (“adiabatic” just means that no energy is lost during mixing):
Dividing by the mass of the final mixture, we obtain the enthalpy per pound:
This represents the specific enthalpy of the final mixture, which we also know is 25wt%. Again
using Figure 6 (Fig. 8.4-3 in Geankoplis), to get this H, we need to choose the temperature as
T = 123°F
Thus the final solution is warmer than the two initial solutions, due to the “heat of dilution.”
Question: when diluting sulfuric acid (which also exhibits a significant heat of dilution), why is
is wiser to add (concentrated) acid to water, rather than water to (concentrated) acid?
M1
2m = (36)
m1 M2
where 2m and m1 are the lengths of the two lines drawn in the schematic of Figure 7: 2m is the
length of the line connecting point “2” with point “m” while m1 is the length of the line
connecting point “m” with point “1”. Note that decreasing the amount of M1 relative to M2
shifts point ‘m’ closer to ‘2’ (i.e. it decreases 2m ). In the limit that M1 = 0, point ‘m’ coincides
with ‘2’, which makes sense. Armed with (36), we can solve such problems graphically without
doing any balances.
By collecting together terms multiplied by either M1 or M2, these two relationships can be
rearranged to show
H 2 − H m x2 − xm M1 2m
= = = (37)
H m − H1 xm − x1 M 2 m1
H 2 − H m H m − H1 rise
= = = slope
x2 − xm xm − x1 run
where the ratio represents the slope of either of the two line segments in Figure 7. That the two
line segments have equal slopes and share a common point means that they are part of the same
straight line. This serves a proof that the mixing point “m” lies on a straight line connecting “1”
and “2”.
To find the exact location of the mixing point on this straight line, consider the two shaded
right triangles shown in Figure 7. Since the slope of the two hypotenuses is the same, and the
other sides are either horizontal or vertical, the corresponding angles of the two triangles must be
the same for all three angles defining either triangle. This means the two triangles are similar in
the geometric sense. For two similar triangles, the ratio of lengths of corresponding sides are the
same for all three sides. (37) gives the ratio of the heights (i.e. ∆H) as M1/M2; therefore the ratio
of the lengths of the two hypotenuses must also equal M1/M2, which completes the proof of (36).
M1 l2
M1l1 = M 2 l2 or =
M 2 l1
Notice the inverse relationship between the mass and the distance from the pivot.
MULTIPLE EFFECTS
We will now look at two methods to improve steam economy of an evaporator:
1. multiple “effects”
2. vapor recompression
?
the product streams is water vapor — also known as steam.
Steam has a utility cost associated with it. Instead of
discarding this vapor product, why not use this vapor in place
of steam (which would save this utility cost)? The answer is
V
simple: if the vapor replaced the steam, there would be no
driving force for heat transfer because the “steam” and the L
liquid would have the same temperature: Ts = T and ∆Tm = 0.
You can have any number of evaporators connected together in this way. Each evaporator in this
sequence is called an effect, and the entire sequence is called a multiple-effect evaporator. In
the figure above, we show a double-effect evaporator. In order to maintain a driving force for
heat transfer in each effect, we must have
T2 < T1 < Ts
♣ Assuming we have saturated liquid and vapor in each effect, the pressure in any effect is just the vapor pressure at
the temperature for that effect. Recall that vapor pressure increases with temperature.
This requires some device to lower the pressure in the last effect. While a mechanical vacuum
pump (a compressor applied in reverse) could be used, it is usually much cheaper to use a
barometric condenser (see description in section 8.6C of Geankoplis).
q2 = U2A2∆T 2 (38)
∆T 1 = Ts - T1 > 0
q1 = U1A1∆T 1 (39)
Approximation Solution
For simple designs where BPE and ∆Hdil are negligible and where the sensible heat needed
to preheat the feed up to its boiling point is negligible, the same latent heat added to the vapor in
the first effect is recovered by condensing in the second effect:
Usually the pressure in the last effect and the pressure of the steam are known. Using the steam
tables, we can convert these pressures into temperatures. The unknowns are the temperatures
and pressure in the other effects.
Example: given Ts and T2, find T1. Solution: (38) and (39) into (40) yields
This represents one equation in one unknown (T1), assuming that the U’s and A’s are known.
If instead, we have N effects, the analog of (38), (39) and (40) yield N-1 equations in N-1
unknowns: equating the heat duties
q1 ≈ q2 ≈ ... ≈ qN
♦ The expression for ∆T which follows assumes that the vapor produced in the first effect condenses in the steam
chest of the second effect at the same temperature at which it boiled. If BPE is significant, the vapor will condense
at a lower temperature than it boiled (even if the pressure is the same). See discusion of BPE starting on page 57.
Now Ts and TN are known, but T1, T2 … TN-1 represent N-1 unknowns. These equations are
strongly coupled. A simpler approach is to try and solve for the ∆Ti’s instead of the Ti’s. We
will illustrate this by re-solving the last example. From (38) and (39), we have
q
∆T1 = 1 (41)
U1 A1
q
and ∆T2 = 2 (42)
U 2 A2
2 q1 q
Adding (41) and (42): ∑ ∆T=i + 2
U1 A1 U 2 A2
(43)
i =1
q1 1
q1 = q2
∆T1 U1 A1 U1 A1
= ⇒ (44)
∑ ∆T q1 U1 A1
q
+ 2
U 2 A2
1
U1 A1
+ 1
U 2 A2
Since Σ∆T can be easily calculated from the info given in the problem statement:
1
U1 A1
=∆T1 (Ts − T2 )
1 + 1
U1 A1 U 2 A2
q
(41): ∆Ti = i for i = 1,2, ... N (46)
U i Ai
qi
∆Ti U i Ai
=
∑ ∆T q1 U1 A1
q
+ 2
U 2 A2
q
+ + N
U N AN
Applying the analog of (40) (q1 = q2 = ... = qN ≡ q), we can cancel out the qi, leaving:
1
∆Ti U i Ai
= (47)
∑ ∆T 1U1 A1 + 1U 2 A2 + + 1U N AN
When BPE is negligible, the sum of the driving forces can be calculated explicitly without
knowing the temperatures in the intermediate effects (i.e. without knowing Ti for i = 1, 2 … N–
1):
Notice that the all of the unknown temperatures cancel out. Substituting (48) into (47):
1
U i Ai
∆Ti = (Ts − TN ) for i = 1,2, ... N
1 + 1 + + 1
U1 A1 U 2 A2 U N AN
In the simple case in which all the A’s are equal and all the U ’s are equal, this last result reduces
to
1
∆=
Ti (Ts − TN ) for i = 1,2, ... N (49)
N
In short, the total driving force is divided equally among each of the N effects.
Performance
Each effect results in the evaporation of water at a rate equal to q/λi. Summing the rate of
evaporation in each effect, we obtain the total rate of evaporation:
N q
mV = ∑ i
λ
i =1 i
Steam is used only in the first effect. The rate of steam consumption in the first effect is
q
m s = 1
λs
N q
∑ λi
m V i =1 i q1= q2= = q N N
1
= ⇒ λs ∑ ≈N
m s q1
i =1
λ i
λs
If we neglect changes in the latent heat with temperature and pressure, then all the λ’s are equal
and the sum reduces to N. If changes in the latent heat with temperature are considered, the sum
will be less than N.
Summary: For the same capacity, steam requirements are significantly less for multiple-effect
evaporators, than for single-effect evaporators: generally, the economy increases with the
number of effects. The price that is paid for this increased steam economy is greater equipment
costs: generally the cost of the equipment is proportional to the number of effects.
The reason that the economy increased is that we need to supply steam only to the first effect.
The source of heat for 2nd, 3rd, ... Nth effects is supplied by condensing the vapor produced in the
previous effect.
For simplicity, let’s assume the heat-exchanger area and heat-transfer coefficient are equal for
each effect:
BTU
U=
1 U=
2 U=
3 500 2
ft -°F
and A=
1 A=
2 A=
3 ?
T − T3 250 − 125
∆T1 =∆T2 =∆T3 =s = =41.7 °F
3 3
Ts − T1 T1 − T2 T2 − T3
Although the heat duty in each effect is the same, the amount evaporated is not quite equal
because the heat of vaporization depends on temperature. The total of evaporation rates in each
effect is
q q q lbm
+ + =
44, 000
λ1 λ 2 λ3 hr
btu
q 3.66 × 106
=
hr
Using the same U as in the previous example, we can estimate the hxer area:
q
=A = 704 ft 2
U ∆Ti
Of course, we have three heat exchangers, each of this area. The steam requirements can be
estimated as
q lbm
=
S = 1.547 × 104
λs hr
F − LN
and the steam economy is = 2.84
S
Comment: The economy is slightly less than the number of effects, but significantly greater than
unity. While the improved steam economy is the main advantage of using multiple effects, the
tradeoff is greater heat exchanger costs: we need 3×704 = 2112 ft2 compared to only 768 ft2 for
a single effect evaporator. Of course, the heat-exchanger cost is a capitol expense, which you
pay only once, whereas steam is an operating cost which is ongoing. As the cost of energy rises,
the use of multiple effects becomes more attractive.
Consider what happens to the temperature of the vapor leaving the first effect as it condenses
in the hxer of the second effect. When it leaves the first effect, the vapor is at the boiling point
of the solution in the first effect
But this same vapor condenses at the boiling of pure water (not the boiling point of the solution
from whence it came). The difference between the two temperatures is the BPE for the first
effect BPE1.
This means that the driving force for heat transfer in the second effect is not T1 – T2. Instead:
∆T2 = T1 – BPE1 – T2
The same correction must be made for each of the subsequent effects. Thus the effect of boiling
point elevation is to reduce the driving force for heat transfer in the downstream effects
To the extent that sensible heat is still negligible compared to latent heat, we can still claim the
heat duties of each effect is the same: q1 = q2 = ... = qN. (46) and (47) still hold but (48) must be
replaced by (50). The final result is
N −1 1
U i Ai
∆Ti = Ts − TN − ∑ BPEi (51)
1
U1 A1 + U 2 A2 + + U N AN
1 1
i =1
Boiling point elevation thus lowers the driving force for each effect and lowers the amount of
evaporation which occurs in each effect. Unfortunately, BPEi cannot be evaluated until Ti is
known. So the problem is now trial-and-error.
Equality of the heat duties (i.e. q1 = q2 = ... = qN) represents a first approximation. Sensible
heat is generally not entirely negligible and this causes the heat duties to vary from one effect to
another. If the heat duties are not equal, they will not cancel out as in (47). This leads to the
following generalization of (51):
qi
N −1 U i Ai
∆Ti = Ts − TN − ∑ BPEi (52)
q1 q q
i =1 U A + 2 U A + + N U A
1 1 2 2 N N
Unfortunately, the individual qi’s cannot be calculated until the temperatures are known, so the
problem requires a trial-and-error solution, which is described in the following recipe.
Example: A triple-effect evaporator is being used to concentrate 50,000 lb/hr of a 10wt% sugar
solution at 80ºF to 50wt% sugar. BPE (in ºF) and heat capacity cp (in btu/lb) can estimated
from♥
ps = 29.8 psia
p3 = 1.94 psia
♥ Unlike orange juice, sugar is generally a low-molecular weight compound. For example, sucrose has a molecular
weight of 342 g/mol.
A1 = A2 = A3 = ?
Calculate the area of each effect and the steam requirements. The flow diagram containing some
of the input data follows:
Find: S and A1 = A2 = … = AN
Step 1) Knowing the pressures pN and ps and the composition xN in the last effect, use a Dühring
chart or other thermo data provided to determine the boiling point of the solution TN and then use
the steam tables to determine the temperature of the steam Ts.
From the steam table, the boiling point of pure water at 1.94 psia (or 29.8 psia) is found as
124.921ºF (or 250ºF). Add the boiling-point elevation BPE(0.5) = 4.400ºF gives the boiling
temperature of the liquor in the last effect as
Step 2) calculate final liquor flowrate and total capacity from overall mass balances
x F N
LN = F
xN
and ∑ V=i F − LN
i =1
3
0.1 lb lb lb
=L3 = 50, 000 10, 000
0.5 hr hr
and ∑ Vi = 50, 000 − 10, 000 = 40, 000 hr
i =1
Step 3) As a first guess, assume the vapor flowrates from each effect are equal:
F − LN
Vi = for i = 1,2 … N (53)
N
lb
40, 000
V= hr= 13,333 lb
1 V=
2 V=
3
3 hr
Step 4) calculate the liquid flowrate and composition leaving each effect using total mass and
solute mass balances:
x F
Step 4a) L1= F − V1 and x1 = F for i = 1
L1
x L
Step 4b) =
Li Li −1 − Vi and xi = i −1 i −1 for i = 2, 3 … N–1
Li
lb 50, 000
L1 = 50, 000 − 13,333 = 36.667 = = 0.1364
and x1 0.5
hr 36, 667
lb 36, 667
L2 = 36, 667 − 13,333 = 23.333 = = 0.2143
and x2 0.1364
hr 23,333
We don’t need to do the calculations for i = N (except possibly as a check) because the results for
L3 and x3 are already known.
Step 5) Knowing all the concentrations xi we can now estimate the boiling-point elevation for
each effect. Although we don’t know the pressures (except in the last effect), BPE is practically
independent of pressure. This is why the lines of a Dühring chart are nearly parallel.
BPEi ≈ BPE(xi)
Using the known xis and the BPE formula provided in the problem statement, the BPEs can be
calculated:
Step 6) Estimate the driving force in each effect using (52) in which the equal areas cancel out:
qi
N −1 Ui
∆Ti = Ts − TN − ∑ BPEi for i = 1,2 … N (54)
q1 q q
i =1 U + 2 U + + N U
1 2 N
Although the q’s will turn out to be unequal in Step 10), the first time through this loop, we will
assume them to be equal (in order to get started):
N −1 1
Ui
∆Ti = Ts − TN − ∑ BPEi (55)
1 + 1 + + 1
i =1 U1 U2 UN
Comment: Geankoplis corrects the value of ∆T1 “for the cold feed”. This correction (if done)
takes the form:
Fc pF (T1 − TF )
( ∆T1 )corrected =
1 + ∆T1
V1λ s1
where the heat capacity cpF, heat of vaporization λs1 and Ti are evaluated below. This increase
in ∆T1 is then subtracted equally from the remaining values so that the sum of the ∆T’s is
unchanged.
At the end of the first full iteration, we will calculate each of the q’s including sensible heat
effects; the ∆T’s will be re-calculated from (54) instead of (55). If such rigorous calculations are
planned, the main effect of making this correction is perhaps to reduce the number of iterations
by one (is it worth the effort?). On the other hand, if we are willing to settle for an approximate
solution without iteration, then this correction is probably worthwhile. In the example which
follows, I will not make this correction.
With an eye to eventually evaluating (55): First, the sum of the driving forces is calculated as
2
Ts − T3 − ∑ BPE j = 250 − 129.3 − ( 0.6 + 1.2 ) = 118.8ºF
j =1
Notice that the BPE in the last effect is not included in the above summation. Next, we will
calculate the sum of 1/Ui:
1 1 1 1 1 1 ft 2 ⋅ hr ⋅ °F
+ + = + + = 0.00968
U1 U 2 U 3 500 350 200 BTU
Using the given overall heat transfer coefficients we can now estimate the driving force in each
effect from (55):
1
∆T1= 118.8 = 22.3°F
500 × 0.00968
In this example we will pass on making the correction for the cold feed.
Step 7) Calculate temperatures of the hot and cold streams in each effect. The hot fluid
temperature is the temperature at which the vapor from the previous effect condenses. If there
were no BPE (as with saturated steam used in the first effect), then the hot fluid temperature in
the ith effect is just be the boiling temperature from the previous effect. But if there is BPE, then
the hot fluid temperature in the ith effect is the condensation temperature of the vapor from the (i-
1)th effect, which in turn is the boiling temperature of the solution in the (i-1)th effect minus the
BPE:
The cold fluid temperature is the hot-fluid temperature minus the driving force:
T=
i Thi − ∆Ti for i = 1, 2 … N–1 (57)
There is no need to evaluate again TN in the above calculation, because its calculation in Step 1)
is exact; however repeating the calculation might be worthwhile as a check on other calculations.
The hot- and cold-stream temperatures in each effect can be estimated from (56) and (57) as
Step 8) Knowing the temperature and composition of all the streams, we now determine their
enthalpies from given thermo data. All vapor streams are pure water, so we use the steam tables
to evaluate their enthalpies. We look for the enthalpy of saturated vapor at the boiling
temperature:
Remember the temperature of the vapor leaving the ith effect is Ti.
Since steam tables will be employed in most calculations, we adapt the reference state used
there (H ≡ 0 for pure liquid water at 0ºC). If enthalpy data is available for the concentrated
solution (for example, H for NaOH solutions is given in Figure 6 or Fig. 8.4-3 in Geankoplis),
then use that:
Important: If you use more than one source of enthalpy data, make sure all sources employ the
same reference state for enthalpy. We have already noted that the steam tables from Geankoplis
define H ≡ 0 for saturated liquid water at 0ºC. By inspection of Figure 6 we see that H = 0 for 0
wt% NaOH at 32ºF. So this second source is OK.
When heat of dilution is negligible, enthalpy data for solutions can be estimated by
specifying the heat capacity as a function of concentration. Then enthalpy can be calculated
from
where the reference temperature for enthalpy (that employed in the steam tables) is 0ºC or Tref =
273ºK.
We will also need the heat of vaporization of the “steam” used in the ith effect. Since some
of the vapors are superheated, we need to correct as follows:
Knowing their temperatures and compositions, we can now estimate the enthalpy of each stream.
The liquid heat capacities and enthalpies♣ are found to be
♦ Owing to BPE, the vapor produced by boiling the solution is superheated (the pressure of the vapor is lower than
the vapor pressure of pure water at this temperature); we are neglecting the effect of pressure on vapor enthalpies
(see justification in footnote ♠ on page 43).
♣ Liquid enthalpies were calculated as c (T - 32°F). This assumes that all sugar solutions (regardless of sugar
pi i
concentration) have zero enthalpy at the freezing point of pure water. Actually the reference condition for enthalpy
is pure water at 32°F. There can be a slight effect of concentration: note that the isotherms in Figure 6 for NaOH
solutions at 40°F and 50°F are nearly horizonal, but not exactly.
cpi = (0.924, 0.880, 0.720) btu/lb-ºF and Hi = (181.3, 141.3, 70.5) btu/lb
The vapor enthalpies are read directly from the steam table using the boiling temperature Ti:
To evaluate the enthalpy of the condensate leaving the ith effect, we need the condensation
temperatures calculated from (56):
The heat given up by condensing the superheated vapor from the previous effect can be
estimated from (58) and (59):
Step 9) Recalculate all the flowrates using mass and energy balances:
F = L1 + V1
Li-1 = Li + Vi for i = 2, 3 … N
These are 2N equations in 2N unknowns: L1, L2 … LN-1, S, V1, V2 … VN. The steam flowrate S
replaces LN in this list of unknowns since LN is known. Collecting terms involving the unknown
flows on the left-hand side and know flows on the right-hand side, gives the following matrix
equation to solve (shown for N=3):
H1 0 −λ s1 H v1 0 0 L1 FH F
− H1 H2 0 −λ s1 Hv2 0 L2 0
0 −H2 0 0 −λ s3 H v3 S − L3 H 3
= (61)
1 0 0 1 0 0 V1 F
−1 1 0 0 1 0 V2 0
0 −1 0 0 0 1 V3 − L3
The first three rows are the three energy balance equations and the last three rows are the mass
balances. While this could be solved by hand using elimination, it can be very easily solved in
Mathcad (or other software) using the LSOLVE function.
In preparation for doing mass and energy balances we construct the following two matrices
M and B in (61):
This matrix equation MX = B was then solved for X using the LSOLVE function in Mathcad.
The result for X is
L1 37, 611
L
2 24, 222
S 19,934
=
X = lb/hr
V1 12,389
V2 13,389
V3 14, 222
Step 10) Calculate the heat duties for each effect from
The vapor flowrates V1, V2 … VN calculated in Step 9) will generally differ from those
approximated by (53). If this difference is only a few percent, you can stop here and skip to Step
12). If the difference is significant (and it almost always is significant the first time you reach
this step), you need to iterate in Step 11) until the difference from the previous step is negligible.
The heat duties of each effect are then calculated from (62) using the modified latent heats λsi
from (60):
Notice that the heat duties are not all the same [as we had assumed in Step 6 – see (55)]. This
means we have to iterate at least once more. We will iterate until the vapor flowrate converge
within some tolerance.
Step 11) Repeat Step 4b) through Step 10) using the updated ∆Ts until the flowrates converge.
repeat Step 5&6) We then recalculate the driving forces, using the qi’s:
repeat Step 7-10) Calculate a new set of flowrates (all three sets are shown):
Notice that the changes are relatively small after the 2nd iteration, so we stop.
Step 12) Once converged, the areas can be calculated from the Hxer design equation:
qi
Ai = (63)
U i ∆Ti
The area should be the same for each effect, because that was what was assumed in (54).
The areas should be equal, but are not quite equal because the mass and energy balances are not
perfectly converged. Additional iterations would make them more nearly equal. The area to use
in design is the average of the three:
A = 1129 ft2
40, 000
The economy is = 2.03
19, 722
Although it was not asked for, you might be interested to know how the results depend on the
number of effects. Here is a summary:♦
N
S VN ∑ Ai (ft2) economyAs
(lb/hr) (lb/hr) i Economy Total Area
1 44,843 40,000 703 0.89 expected,
steam 4.5
2 25,575 20,450 1,454 1.56
3 18,984 14,000 2,208 2.11 consumpti 4
4 15,717 10,800 2,978 2.55 on S 3.5
decreases
3
with N (or economy increases). Somewhat surprising is the
large increase with N in total area of all N heat exchangers: 2.5
the total area is essentially proportional to N. The reason 2
for this large increase can be rationalized as follows. Recall 1.5
the design equation for any of the heat exchangers: 1
0.5
qi
Ai = 0
U i ∆Ti
1 2 3 4
Let’s again neglect sensible heat effects and BPE as well as
the T-dependence of the heat of vaporization. Then
F − LN
∑ qi ≈ λ
i
F − LN T − TN
and qi ≈ and ∆Ti ≈ s
Nλ N
Substituting these results back into the expression for Ai and summing:
F − LN
qi Nλ = F − LN F − LN
Ai = ≈
(Ts − TN ) λ
and ∑ Ai ≈ N (T ∝N
U i ∆Ti Ts − TN
i s − TN ) λ
N
This leads to the prediction that the total area is proportional to N, as in our more rigorous
calculation in the example. Thus the substantial increase in steam economy comes at the price of
a substantial increase in heat exchanger area. Is it worth it? The answer depends on the relative
cost of heat exchanger area and steam.
♦The economy for N=3 is slightly different from what was calculated in the detailed example above. This is
because in varying N I have used U = 500 for all hxers, whereas in the above example, each hxer had a different
value.
without compression: ∆T = Ts - T = 0
Adiabatic compression of water vapor (treated as an ideal gas) increases its temperature
according to♣
1− 1
Tout pout γ
= (64)
Tin pin
where “out” and “in” refer to the outlet and inlet of the compressor and where γ is the heat
capacity ratio, which is
ideal
c p gas cp water
γ≡ ≈ = 1.324 or 1− 1 =0.245
cv R γ
cp −
M
psia
4.7
Temperature (degR)
=1
r p in
n fo
1000 ssio
pr e
C om
ab atic
t Adi
Boiling Poin
1 10 100
Figure 8: Red curve show how temperature increases with adiabatic compression, starting with
saturated vapor at 14.7 psia. For comparison, the boiling point is shown by the black curve.
This increase in compressor outlet temperature is plotted in Figure 8. Note that the outlet
temperature rises faster with increases in outlet pressure than the boiling point increases with
pressure. Thus adiabatic compression of a saturated vapor nearly always produces a superheated
vapor. When heat is removed from this compressed vapor, it will eventually condense at its
boiling point for that pressure.
Increasing the pressure of the vapor also increases its boiling point, so that as the vapor
condenses, it stays at a high temperature Ts so we always have some driving force.
Example: Suppose we compress the vapor above a solution of organic colloids (negligible BPE)
from pin = 1 atm to pout = 2 atm. Evaluate ∆T for use in the Hxer design equation.
Solution: The inlet temperature for the compressor is the boiling point for 1 atm: Tin = 212ºF.
The outlet temperature for the compressor can be calculated from (64): Tout = 339ºF. But this is
not the temperature of the hot fluid which should be used. The compressed vapor is superheated.
The boiling point at pout = 2 atm is Ts = 250ºF. It is this second temperature which should be
used in preliminary designs:
Unfortunately, raising the temperature and pressure of the vapor actually lowers its latent
heat per pound
so the vapor no longer has enough latent heat to boil the same mass of liquid. Thus make-up
steam is required, although generally much less steam is needed than if the vapor’s latent heat is
not recycled.
The reduction in steam is not without cost. We still have to pay the utilities (electricity or
steam turbine) to operate the compressor. Generally, compressors are expensive to buy,
maintain, and operate; so you have to do some economic calculations to see if using a
compressor is economical. At the usual conditions of an evaporator, the vapor can be treated as
an ideal gas. Then the work required for adiabatic compression can be estimated from♦
1− 1 γ
γ RTin p
= Btu
=W pr out
− 1 [ ] (65)
γ −1 M pin lb
where the subscripts “in” and “out” refer to conditions in the inlet or outlet to the compressor and
R is the universal gas constant. Now (65) gives the work needed to compress a unit mass of gas,
under ideal conditions (no irreversible losses of energy). To get the power requirements, we
must multiply by the mass flow rate and divide by the mechanical efficiency (η) of the
compressor [typically η ≈ 0.7]:
=P
(=
m f − m ) W pr
[=]
Btu
[ ] kW
η hr
By compressing the vapor to raise its temperature, we are lowering its latent heat of vaporization.
This means that, even if we completely condense all the vapor, we still won’t have enough latent
heat to boil the liquor at its lower temperature. The difference must come from make-up steam.
To estimate the make-up steam requirements, let’s return to the equation (32) for the total
heat duty added to the cold fluid:
(
)
q = m f − m H v + mH
−
liquor
m f H f
vapor feed
In the simple case where BPE and ∆Hdil are negligible (i.e. for inorganic colloids), Hf and H
differ from the enthalpy of saturated liquid only in sensible heat. To the extent that sensible heat
is negligible,♣ we can substitute the enthalpy of saturated liquid for Hf and H:
( ) (
= m f − m H v − m f − m H L)
( m f
= − m ) ( H v − H L )
This heat duty is partially offset by the latent heat recovered by condensing the compressed
vapor. The remainder comes from make up steam:
=
q (
m f − m ) λ s
+ m s λ s
(67)
compressed vapor make-up steam
where we are assuming we are compressing the vapor to the same pressure as the make-up
steam. Eliminating q between (66) and (67), we can solve for the steam requirements:
λ
(
m s ≈ m f − m − 1)
λs
Notice that any compression allows us to recover all the latent heat of the vapor, thus
significantly reducing the steam consumption which would be otherwise be required without
compression. However greater compression actually increases steam consumption (because λs
decreases with pressure). Rearranging this expression, we can calculate the economy:
m f − m λs
≈ >> 1
m s λ − λs
which is now very much greater than unity. However, this way of computing economy is unfair
because the reduction in steam consumption is at the cost of operating the compressor. In
computing the “economy” for an evaporator with vapor recompression, it is customary to add an
♣ Sensible heat is often negligible compared to latent heat. To see why, note that the latent heat of vaporization of
water at 1atm, 212°F is 970.3 btu/lb, where 1 btu is the amount of sensible heat required to raise the temperature of
one lb of water 1°F. This means that, for sensible heat to equal latent heat, I need to raise the temperature of water
970°F. A temperature rise of only 30°F, say, represents a sensible heat change of 30/970 or 3% of the latent heat of
s.
water. This is an estimate of the relative error made in this estimate of m
equivalent steam rate for the compressor to the make-up steam rate before dividing into the
capacity.
• If the plant uses a steam-powered turbine or if the plant uses steam generators to make its
own electricity, then the apparent steam usage is
P btu hr lb
( m s )app =
m s + [=
] [=
]
λ ′s ηs btu lb hr
where P is the electrical power supplied to the compressor, λ ′s is the latent heat of the steam
used (the prime is added because the pressure of the steam used here might be different from that
used in the evaporator) and ηs is the efficiency of electricity generation from steam (typically
0.35).
• On the other hand, if electricity is purchased from a utility to drive the compressor, then it
makes sense to convert the electricity usage into an equivalent amount of steam (on a
financial basis):
$
( m s )app = m s + P ×
$
kW-hr
lb of stm
If P is in kW then the last term yields lb/hr of steam. Regardless of which method of converting
the power usage is appropriate, the economy is computed as:
m f − m
economy ≡
( m s )app
which decreases the latent heat of vaporization, requiring more make-up steam. Only when the
cost of the heat exchangers is considered does an optimum exist. The cost of the heat exchanger
decreases with pressure because a larger driving force for heat transfer leading to smaller area.
Comment: It might appear surprising that the amount of steam increases with the pressure of the
compressed vapor. After all, we are arguing that compression reduces steam requirements; you
might expect that more compression reduces steam requirements more – but it doesn’t. The
reason is simple: the higher the pressure of the compressed vapor, the more its latent heat of
vaporization λs is reduced (see example on page 47) and the more make-up steam is required
from (67) for the same vapor flowrate mf –
m and the same heat duty q. Cos ts of Simple Ev aporato r
2 .10
6
14.7 23.5
At the minimum, the total costs with
vapor recompression are about $159,000
1.5 .10
6
per year. If we use all steam instead of
Annual Costs ($)
vapor recompression plus make-up steam,
the cost of the steam alone would exceed 1 .10
6
0
15 20 25 30
The next unit operation we will consider is distillation in which we attempt to separate a
mixture of two miscible components by creating a second phase: either partially vaporizing a
liquid or partially condensing a vapor. This is quite different from evaporation because both
components are present in both liquid and vapor phases (whereas evaporation involves an
nonvolatile solute).
In general, the liquid and vapor phases (created by partially vaporizing a liquid feed or
partially condensing a vapor feed) will have different concentrations of the two components.
This is the basis of separation by distillation. We will start by reviewing the principals of vapor-
liquid equilibrium, which limits the amount of separation that can be accomplished.
z
One Component
T
The phase diagram for water is shown below:
105
piston
104
Liquid
103
Pressure (kPa)
102
Solid
101
Vapor
100
10-1
10-2
-50 0 50 100 150 200
Temperature (oC)
To comprehend the meaning of such phase diagrams, it’s helpful to recall Gibbs phase rule:♦
P+F=C+2
where P is the number of phases present, F is the number of degrees of freedom and C is the
number of components. For example, in the phase diagram above, we have only one component
present (water♥), so C = 1 and Gibbs phase rule says
for C = 1: P+F=3
Suppose we are interested in portions of the diagram having only one phase present at
equilibrium: P = 1. Then Gibbs predicts F = 2; i.e. two degrees of freedom. This means that we
can vary both pressure and temperature over some range and still have only one phase present.
In other words, single phase regions are represented on this P-T chart as 2-D areas, which of
course are marked with the words “liquid”, “vapor” or “solid” representing the single phase
found under those conditions.
Two Components
♦ The somewhat silly mnemonic I use to remember this is “Police Force = Cops + 2”.
♥ It’s important to exclude air which might otherwise be present. Air (itself a mixture of nitrogen and oxygen)
behaves like a second or third component and will affect phase equilibria.
(C8). The vapor-liquid phase behavior of two-component systems is a little more complicated
because I have an additional variable — composition — as well as T and P. Thus we need a 3-D
graph to represent the phase diagram for two component mixtures. The vapor-liquid portion of
the phase diagram for two components is shown in Figure 10.
200
Vapor 200
L&V
Temperature (degC)
Temperature (degC)
150 150
100 100
50 50
0 L&V 0
Liquid
-50 -50
1
(a tm) 1
(a tm)
0.2 sure 0.8 sure
0.4
0.6
0.8 0.1
Pres 0.6
0.4
0.2 10
Pres
1.0 0.0
Mole Fraction Mole Fraction
Pentane Pentane
Figure 10: Phase diagram for mixtures of pentane and octane. The two views are from the low-
pressure, pentane-rich side (left) or from the high-pressure, pentane-lean side (right).
Suppose we fix the pressure inside our piston. Instead of a single temperature at which both
liquid and vapor can co-exist (as in Figure 9), there is a range of temperatures. This can more
easily be visualized by taking a vertical slice though the 3-D phase diagram of Figure 10 (a
P=const plane).
The composition (y) of this first bubble can also be read from the Txy diagram. In general, the
vapor will be richer than the liquid in one component and leaner in the other. This difference in
composition of the vapor and liquid can serve as the basis for separating two or more
components of a mixture. By completely condensing the vapor, I produce a liquid which is
richer in one component than the original liquid. Thus I have achieved a degree of separation of
the two components.
For a binary system having two phases, we denote the composition of the liquid and vapor
using x and y:
Keep in mind that liquid and vapor must have the same temperature at equilibrium. The line
connecting the two points representing vapor and liquid is called a tie line. It’s horizontal on a
Txy diagram.
y
Ki ≡ i
xi
The subscripts i and j refer to components in the two-phase mixture. Generally the most volatile
component will have the largest distribution coefficient and will have its concentration in the
vapor increased (compared to the liquid) by the largest percentage. The ratio of the distribution
coefficients of two components is the relative volatility.
For a binary mixture we don’t need subscripts because x and y denote the mole fractions of
the more volatile component and 1-x and 1-y are the mole fractions of the other component.
Combining the two equations above for this binary mixture:
K more vol y x
α≡ =
K less vol (1 − y ) (1 − x )
αx
Solving for y(x): y= (68)
( α − 1) x + 1
Although relative volatility is generally a function of temperature [this implies that α = α(x) in
(68)], for ideal solutions of components having nearly equal molar heats of vaporization, the
relative volatility is nearly constant (see “Constant Relative Volatility Approx.” on page 84).
Then the VLE for an N-component mixtures can be entirely represented using only N-1 relative
volatilities. For a binary mixture, only one constant is needed and the xy diagram can be
represented by (68).
To see why relative volatility might be independent of temperature, let’s look at ideal
solutions whose partial pressure pi is given by Raoult’s law:
P = total pressure
Substituting xi from Raoult’s law into the definition of partition coefficient, we have
yi pi P Pi′
Ki ≡ = = (69)
xi pi Pi′ P
Ki P′
and αij= = i (70)
K j Pj′
d ln Pio λ d ln Pio λ
= i or = − i
dT RT 2 d (1 T ) R
where λi is the molar heat of vaporization (in the past we have used λ to denote the heat of
vaporization per unit mass, rather than per mole) of the pure component. Integrating this with
respect to T (ignoring the weak temperature dependence of λi) yields♣
λi
i′ Ai −
ln P= (71)
RT
♣ A somewhat more accurate but empirical approximation is often used to summarize vapor pressure data in
handbooks (e.g. Lange’s Handbook of Chemistry, 13th edition, pages 10-29 to 10-54):
Bi
i ′ Ai −
ln P=
T + Ci
It is called Antoine’s equation and it accounts for the variation of λi with T. Other handbooks (e.g. CRC Handbook
of Chemistry & Physics, 57th edition, pages D-183 to D-215) give vapor pressures at several temperatures and leave
you to develop your own interpolation scheme.
λ
or Pi′(T ) Ai′ exp − i
= (72)
RT
where Ai is just some integration constant. Since vapor pressure depends strongly on
temperature, so does the partition coefficient Ki given by (69). However, the relative volatility α
is much less dependent on temperature. The reason becomes clear when we substitute (72) into
(69):
P ′ Ai′ λi − λ j
αij =
≡ i exp −
Pj′ A′j RT
Note that when the two molar heats of vaporization are equal, the temperature-dependence
vanishes. Even if they are only approximately equal, the temperature-dependence of the relative
volatility will be much less than of the vapor pressures. For this reason, the assumption of
constant relative volatility is frequently used.
xB = 1-xA
yB = 1-yA
Adding the partial pressures from (73) and (74) we obtain the total pressure:
P = pA + pB = xAPA'(T) + (1-xA)PB'(T)
Keep in might that the vapor pressures in these equations depend on temperature. So if you
specify a T and P, I can determine the vapor pressures and calculate the composition of the vapor
and liquid which will be at equilibrium.
Given: T,P
Thus I can construct a complete phase diagram for the mixture using only vapor-pressure data
for the pure components.
Example.♥ n-heptane and n-octane form a nearly ideal mixture. Use Raoult's law to construct a
phase diagram at constant pressure (1 atm).
Solution. We first look up the normal boiling points of the two pure components.
♥ Illustration 9.1 from Treybal. A similar example is 10.1 from S&VN (p305).
126
124
122
120 Vapor
Temperature (deg C)
118
116
114
112
L&V
110
108
106 Liquid
104
102
100
98
0 0.2 0.4 0.6 0.8 1
mole fraction of heptane
NONIDEAL BEHAVIOR
Fig. 1. Txy and xy diagrams for nonideal systems. On left is carbon disulfide in acetone at one
atmosphere pressure. On right is acetone in chloroform. Taken from Treybal, 2nd ed., pages 290 and
295.
Examples Txy and xy diagrams for binary mixtures which display significant nonideal
behavior are illustrated in Fig. 1. In particular notice that the curve on the xy-diagram can cross
the 45° line. This point of crossing is called an azeotrope:
azeotrope - vapor and liquid have the exactly the same composition at equilibrium (i.e. y = x)
(other than the trivial case of pure components)
Notice also that the boiling temperature of the azeotrope is not bounded by the boiling points
of the two pure components. The example on the left is called a minimum-boiling azeotrope
while the example on the right is called a maximum-boiling azeotrope. The occurrence of
azeotropes can pose difficulty for distillation since producing a new phase by condensation or by
evaporation produces no separation of components.
The solution is to evaporate only partially, because the first vapors which come off will be
the richest in the more volatile component. Vapors produced later will only dilute the
condensate. This process is called batch distillation because we process one batch at a time
rather than continuously processing as with flash distillation.
Let's now try to predict how the composition of the liquid and vapor change as we evaporate
more and more of the charge.
n (t )
x (t ) = A = mole fraction of A in liq. remaining (77)
n (t )
Suppose we evaporate a small amount of liquid at time t over a time interval dt:
ndx = (y-x)dn
dn dx
=
n y−x
Integrating from t=0, at which time x,n have their initial values:
n (t ) x (t )
dn dx
we obtain ∫ n
= ∫ y−x
n0 x0
x0
n dx
Integrating: ln 0 =
n (t ) ∫ y−x
x (t )
x0
dx
Solving for n(t): =(
n t ) n0 exp −
y − x ∫ (79)
x (t )
at equilibrium: y = f (x)
For each x between x(0) and x(t), we can look up the corresponding y(t)
and compute the integrand in Rayleigh's equation 1/(y-x). The integral is the area under the
curve between these two limits.
Step 1: Assign a value to x(t) in the interval 0.2 - 0.5, say x=0.2.
Step 2: Calculate the corresponding y. We expect the vapor to be richer in the more volatile
component. For constant relative volatility, the equilibrium curve is given by (68):
αx
y= 1
( α − 1) x + 1 α = 1.7
0.5
1 0
= 10.2 0 0.5 1
y−x x
0.2 0.5
n (t ) 0.5 dx
= exp −
n0 ∫
y − x 1
10
0.2 y
i
x
i
5
The integrand is plotted in the figure at right. The area
under the curve between the two dotted lines is
0
0 0.5 1
x
i
0.5
dx
∫ y−x
= 2.45
0.2
n (t ) 0.5 dx
Rayleigh's equation yields
n0
exp −
=
∫ =
y − x
0.086
0.2
This represents the fraction of the initial charge which remains in the still as liquid. The fraction
of the initial charge which has been distilled is just unity minus this
n (t )
1− =
0.914
n0
0.8
Step 6: Calculate the condensate concentration
Find: D, B, yD and xB
F=D+B (80)
where F,D,B are total molar flowrates of each stream and x and y are mole fractions of the more
volatile component. The remaining equation comes from requiring the vapor and liquid products
to be at equilibrium:
(80)-(83) represent four equations in the four unknowns. Of course, HF, HD and HB are also
needed. In general, for a given pressure (which we are assuming to be held constant), the
enthalpy of any stream depends on its temperature and composition. For a 2-component system,
having two phases at equilibrium, temperature and composition are not independent. For a given
value of yD (for example), there is a single temperature at which liquid coexists and its
♠ We use x for mole fraction of the liquid and y for mole fraction of the vapor. More generally, when we have a
two-phase mixture (or if we are talking more generally and don’t wish to specify which phase), we will use z for
mole fraction of the mixture (i.e. total moles of the more volative component in both phases divided by the total
number of moles in the two phases).
composition is xB. Thus for a particular pressure and yD, the values of HD, HB, xB and
temperature are unique and can be found from the given “thermo data.” The details will become
clearer in upcoming examples.
D HF − HB Q F
f=
≡ + (84)
F HD − HB HD − HB
The ratio D/F represents the (molar) fraction of the feed which is vaporized. This fraction is
defined as f. Dividing (80) by F and substituting (84):
B D
=1 − =1 − f (85)
F F
z F = yD f + xB (1 − f ) (86)
f −1 z
Solving for yD: =yD xB + F (87)
f f
<0
APPROXIMATE SOLUTION
The solution to the four equations (80)-(83) can be obtained in a number of ways. In some
mixtures (especially of homologous series* like the alkanes), the molar heat of vaporization is
nearly the same for both components. For this reason, we tend to solve all distillation problems
using mole fraction rather than mass fractions, and express enthalpies per unit mole rather than
enthalpy per unit mass. If you can assume that the molar heat of vaporization of the mixture is
completely independent of its composition, then solution is quite easy:
* A “homologous series” is a set of very similar items, which differ in only one property. In the case of alkanes, the
series is a set of components of similar molecular structure which differ only in the number of methyl groups in the
chain.
Let's choose saturated liquid having the composition of the feed as our reference state for
enthalpies and neglect sensible heat effects. Then we can re-state our example problem:
Find: D, B, yD and xB
Approx. Solution: if we neglect sensible heat, the two liquid streams have the same enthalpy. If
we choose the saturated liquid as the reference state for enthalpies, then that enthalpy is zero:
HF = HB = 0
The vapor stream has a higher enthalpy equal to the heat of vaporization
HD = λ
D 0 Q F Q
When HF = HB, (84) reduces to: f ≡ ≈ + =
F λ λ λF
(87) is a straight line (relating the two unknowns xB and yD) which can now be
constructed on the xy diagram. Knowing the value of f we first calculate the slope:
f −1
slope = <0
f
We just need one point to finish the construction. Of course, we could use the y-intercept: yD =
zF/f = xF/f for xB = 0. But this value is frequently off-scale (e.g. yD = ∞ for f = 0 and xB = 0). A
better choice for an arbitrary point is to substitute xB = xF. Then (87) yields
f −1 xF f −1 + 1
=
yD xF + = xF = xF
f f f
operating line: relationship between yD and xB which is imposed by the mass and
energy balances
To find which point on this line is the solution for a given Q, we need also to satisfy (83), which
represents phase equilibrium. This relationship is a curve on the xy diagram, known as the
The intersection of operating line and the equilibrium curve satisfies all four of the original
equations and yields the answer. Note that the vapor mole fraction yD actually exceeds the feed
xF which in turn exceeds the liquid product mole fraction xB:
yD > xF > xB
So we have achieved some degree of separation of the two components in the feed.
In the previous example we assumed the feed was only liquid and we added heat to create a
vapor. Suppose instead, the feed is only vapor and we remove heat (i.e. Q < 0).
Find: D, B, yD and xB
Approx. Solution: In terms of the enthalpies, the only differences compared to Example #1 are
1) that Q<0 (we are removing heat) and 2) the feed is now a vapor so
HF = λ
D HF − HB Q F λ Q F Q
(84) gives: = + =+ =
1+ <1
F HD − HB HD − HB λ λ Fλ
which will be <1 for Q<0. But if we stick with the original definition of f as
D Q
f ≡ =+
1 <1
F Fλ
then (85) – (87) still hold. In particular, the construction of the operating line on the xy-diagram
is the same as in the previous example. Its intersection with the equilibrium curve gives the
answer to this flash problem.
There is also a simple graphical procedure for estimating the fraction vaporized on a Txy
diagram. Suppose you are given a feed composition and (instead of specifying the fraction
evaporated) the final temperature of the mixture at equilibrium.
D z F − xB B yD − z F
=
f = and = (88)
F y D − xB F y D − xB
where D, B and F are proportional to the distances indicated on the Txy diagram above. These
relations are another example of the lever-arm rule (see page 54).
Proof:
D z F − xB
Solving (86) for f: f ≡ = (89)
F y D − xB
which confirms the first equation in (88). Substituting (89) into (85):
B D z − xB
=1 − =1 − F =
y D − xB − z F − xB( ) = yD − z F (90)
F F y D − xB y D − xB y D − xB
determine the corresponding point (xB, yD) on the xy-diagram; we didn’t need to evaluate f first.
The downside of this method of solution is that I don’t know what Q is required; to estimate Q, I
need to evaluate enthalpies and then I am forced to make approximations like those in Examples
1&2, unless I have full enthalpy data (see next section).
RIGOROUS SOLUTION
Our solution to Example #1 (the interaction of the operating line and the equilibrium curve in
Error! Reference source not found. on page Error! Bookmark not defined.) was approximate
because of the two assumptions we made about the enthalpies. Although this approximate
solution works reasonable well for some mixtures, for others it might not be good enough. Then
you will need a lot more thermodynamic data besides latent heat and heat capacity: you will need
an Hxy diagram.
As we shall soon see, we will need to determine the tie line through a particular point. No matter
how many tie lines your graph has, it would be a real stroke of luck if one passes through the
particular point specified. So we need some way to interpolate tie lines. This is accomplished
by means of:
If we form a right triangle out of each tie line, using the tie line as the hypotenuse, then the locus
of the points where the right angle is located is the auxiliary line. Of course, the ends of the
auxiliary line coincides with the ends of the saturated vapor line. By reversing this procedure,
we can locate more tie lines.
Now let's see how this diagram can be use do solve our flash distillation problems. First, I will
give you the recipe and then later I will show why the recipe works.
Example #4:
Solution (recipe):
The main idea which we need to prove is that the point (zF, HF+Q/F) [call this point “F”] lies on
the tie line. We know that (xB, HB) [call this point “B”] and (xD, HD) [call this point “D”] lie on
the tie line since the two streams leaving an equilibrium flash unit are assumed to be in phase
equilibrium with each other. Our approach will be show that the slope of the line drawn from B
to F is the same as the slope of the tie line drawn from F to D.
B yD − z F
Dividing by D and solve for B/D: = (91)
D z F − xB
Q D B
Divide (82) by F: HF + = HD + HB
F F F
Q B B B
HF + = 1 − H D + H B = H D − ( H D − H B )
F F F F
Q
HD − HF +
Solving for B/F:
B
= F
(92)
F HD − HB
Q
HF + − HB
D F
= (93)
F HD − HB
Q
HD − HF +
(92) ÷ (93):
B
= F
(94)
D q
HF + − HB
F
Q Q
HD − HF + HF + − HB
F
= F
yD − z F xF − xB
Note that both sides of this equation have the form of rise/run on an Hxy diagram; thus each side
represents the slope of a line drawn between two points on an Hxy diagram. These two slopes
involve three distinct points:
The equality of these two slopes means that these three points must lie on the same line. That
line connecting B and D must be a tie line. The conclusion is that
This helps a little but it's still a long way from 0.95 and you
don't generate much product with such a small fraction
vaporized.
One solution is to take the distillate product of the first flash and partially condense it. This
will produce a liquid phase which will be leaner in the more volatile component, leaving the
remaining vapor richer:
♦ H -H is not the “heat of vaporization” (since the composition of these two streams is different). Instead heat of
D B
vaporization is the heat to totally vaporize one mole of feed. On the Hxy diagram, heat of vaporization is the length
of a vertical line between the two curves.
In the schematic above, we have repeated this step N-1 times which is required to obtain
some arbitrary enrichment of the more volatile component. We have also separated the heat
exchanger from the vapor space of the flash (for reasons which will become apparent in the next
section).
Another problem is that each flash unit requires a heat exchanger and then we will need one
additional heat exchanger to condense the product from the final flash unit.
COUNTER-CURRENT CASCADE
Both of these drawbacks (very low recovery and very high heat-exchanger costs) can be
removed by a slight modification: we recycle the bottoms of one flash unit into the feed of the
previous flash unit.
To condense the vapor arriving from the i-1th flash unit, we need to transfer some of its latent
heat to the liquid recycled from the i+1th flash unit. For this to work, the temperature of the
vapor must be greater than that of the liquid. Let’s see if this inequality is satisfied: generally,
the vapor will get richer in the more volatile component as we move left to right:
Assuming that the vapor and liquid streams leaving each flash are
in equilibrium, then the corresponding inequalities for tempera-
ture can be deduced from the typical Txy diagram at right:
So the recycled liquid will be cooler than the vapor it's mixed
with. Thus the liquid will absorb some of the heat of the vapor. This condenses some of the
vapor and vaporizes some of the liquid.
Although they are not shown in the schematic of the counter-current cascade shown above,
we still need two heat-exchangers: 1) to vaporize a liquid stream to produce the feed stream
labelled V0 and 2) to partial condense VN to generate the stream labelled LN+1. Now that we
have shown that most of the heat exchangers might be eliminated, let’s return to a discussion of
the countercurrent cascade.
The product recovery for a countercurrent cascade is much better than for the series of flash
units without recycle. We have only one waste stream (the bottoms product from the first unit,
labelled L1) -- instead of one from each flash.
stage -- one unit of a cascade (e.g. one flash unit)
ideal stage -- a stage whose exit streams are at equilibrium
The contacting of vapor and liquid on any given tray is facilitated by bubble caps, one of
which is shown in more detail in Fig. 3. Each tray contains many bubble caps arranged in a
array covering the central region of the tray (see Fig. 4).
Fig. 4. Top view of arrangement of bubble caps on a tray. Taken from Fig. 6.6 of Treybal.
In the particular design shown in Fig. 3, vapor penetrates the tray through a circular hole 2.47
inches in diameter, upward through a 3-inch length of pipe before being diverted downward by
the cap and forced through the liquid outside the pipe.
Fig. 5. Flow patterns under normal and abnormal operating of a bubble cap. Taken from Fig. 6.7 of
Treybal.
The flow pattern under both normal and abnormal operation is shown in Fig. 5.
An alternative tray design is the sieve-plate tray (see Fig. 6). Sieve-plate trays are much
cheaper to fabricate compared to bubble-cap trays since the tray consists of a circular plate with
holes drilled through for the passage of vapor. These holes are generally much smaller in
diameter than the bubble caps.
Fig. 6. Side view of a “sieve plate” tray. Taken from Fig. 6.8 of Treybal.
A disadvantage of sieve plates is that they are much more prone to weeping and flooding
(conditions corresponding to poor contact of liquid and vapor). “Weeping” is the leaking of
liquid down through the small holes of the sieve plate. “Flooding” is the failure of the gas phase
to flow through the small holes; instead gas flow may stop completely.
To convey some idea of the size of distillation tower, see the photo on the next page. This
photo is of a DuPont plant (in Texas) which purifies ethylene. Ethylene is used (for example) in
the production of polyethylene sheet.
The photo shows six separate towers: five are distillation columns and one is an gas-
absorption column (we will cover gas absorption in Chapter 13). The five distillation towers
range in height from 67-197 feet and are 3 or 4 feet in diameter.
• Li, Vi are molar flowrates of total (A+B) liquid or total vapor leaving the ith stage.
• L0 and VN+1 are the feed flows at the top and bottom of the cascade
• xi, yi are the mole fraction of the more volatile component in the liquid and vapor leaving the
ith stage.
Suppose we know everything about the two input streams: their flowrate, composition, and
enthalpy. Suppose we also have all the thermodynamic data: Hxy diagram and xy diagram.
We want to determine the flowrate, composition and enthalpy of the two output streams and of
all the intermediate streams:
I haven't bothered to list the enthalpies or temperatures as unknowns, because if I know the mole
fraction and pressure, I can look up the enthalpy and temperature of saturated liquid or vapor.
As our system, let's select the top "i" stages. An overall mole balance for this system
requires:
Do I have enough equations to calculate the 4N unknowns? How many independent equations
do I have?
The additional equation is obtained by assuming that the two streams leaving any given stage are
in thermodynamic equilibrium with each other (same temperature and same chemical potential
for each component). When this assumption applies, we say the stage is ideal. Mathematically,
this means
So, at least when I am dealing with ideal stages, I have now as many equations as unknowns.
Many distillation towers have N=20 stages or more, so I could easily be talking about 80
equations -- a formidable problem. Fortunately, there are assumptions which can make this
problem tractable so that it can be solved easily with pencil and paper.
Assumptions:
Substituting this result into the total mole balance about the top i stages (95) yields:
L0 = Li for all i
In other words, all of the vapor streams have the same molar flowrate VN+1 within the cascade
(although their composition and mass flowrate will change) and all the liquid streams have the
same molar flowrate L0:
L0 = L1 = L2 = ... = LN ( = L, say)
This state is called constant molal overflow. Since the input flowrates L0 and VN+1 are known,
we know all the flowrates. Thus we have just halved the number of unknowns from 4N to 2N.
We can then drop the subscripts on the flowrates and forget the enthalpy balance (97) and the
total mole balance (95). This leaves just the component mole balance (96):
L L
which can be rearranged to: yi +1 = xi + y1 − x0 (99)
V V
m b
operating line — yi+1 vs. xi; a graphical representation of the mass and energy balances.
In the next section, we will see how these two curves can be used to determine the number of
equilibrium stages required to achieve some desired degree of separation.
L
R≡ (100)
D
xD = x0 = y1 (101)
L L
yi +1= xi + 1 − xD (102)
V V
The ratio of the flowrates can be expressed in terms of the reflux ratio R using (100):
L L RD R
= = = (103)
V L + D RD + D R + 1
L R R +1− R 1
1− =1 − = =
V R +1 R +1 R +1
R 1
=
yi +1 xi + xD (104)
R +1
+
R
1
slope intercept
Thus the slope of the operating line (a plot of yi+1 versus xi) is R/(R+1). One point on the
operating line is xi = yi+1 = xD. In other words, if we substitute xD for both xi and yi+1, (104) is
automatically satisfied:
♠ Alternatively, we could “start over” and perform a component mole balance on the top i stages of the cascade
(indicated by the blue line in the sketch).
R 1
=
yi +1 xD + xD
R +1 R +1
Proof:
R +1
= xD = xD
R +1
Armed with the slope and one point we could plot the operating
line on an xy diagram (see figure at right).
xD
second point: 0,
R +1
R 1 xD
Proof: y ( 0=
) ×0+ x=
D
R +1 R +1 R +1
Example #1. Suppose I have a feed stream containing 50% alcohol in water which I want to
enrich to 95% alcohol. How many stages do I need if I use a reflux ratio of unity and what will
be the flowrate of this product?
Find: N and D
Solution: The first step is to draw the operating line on the xy-diagram. One point is (xD,xD) or
x
(0.95, 0.95). The second point is 0, D or (0, 0.475). Knowing two points, we can draw the
R +1
operating line (as in figure above left).* Once the operating line is plotted on the xy diagram
along with the equilibrium curve, we can determine the composition of other streams as follows.
*The operating line is not drawn to scale in these sketches. For that matter, neither is the equilibrium curve.
However, given the equilibrium curve and operating line as drawn, the method of stepping off the number of stages
is correct.
answer: N=4
Comment #2: After taking 4 steps (or any integer number), we don’t end up exactly at the feed
composition (except by some extreme stroke of luck). Instead we end up stepping below the
feed composition. If the trays really are perfectly ideal, and the feed composition is 0.5, the
distillate composition would turn out to be slightly richer than 0.95. Although any real column
has an integer number of trays, it is customary to report the answer as a fraction; for example, to
go from a feed of 0.50 to a distillate of 0.95, we require 3.7 ideal stages (or 3 + 0.7, where 0.7 is
the fraction of the last step which corresponds to the feed composition) in the current problem.
Comment #3: Suppose that our feed composition had been a little lower, say yb'. No matter how
many steps we took, we could never get a vapor composition lower than that at the intersection
of the operating and equilibrium curves. In that case, any solution for the specified reflux ratio is
impossible. Instead the reflux ratio R needs to be increased, which will rotate the operating line
counter-clockwise about the point (xD, xD). This will lower the composition at the intersection
point and make the separation possible. It turns out there exists a minimum value of the reflux
ratio which allows a particular distillate composition xD to be produced from a given feed
composition xF. This called the minimum reflux ratio.
Let’s finish this problem by calculating the flowrates. Since the vapor flowrate leaving each
tray is the same, and because we are creating the vapor to the first tray by completely vaporizing
the feed, we can say that
V=F
The ratio of liquid to vapor flowrate inside the cascade is given in terms of the known reflux
ratio by (103):
L R R =1 1 1 1
= = or =
L = V F
V R +1 2 2 2
1
B= L= F
2
Finally, the distillate flowrate can be calculated from the definition of the reflux ratio: recall
(100):
L R =1 1
=
D = = L F
R 2
A STRIPPING CASCADE
Use of a single cascade with 3.7 stages allows us to
convert a feed of 50mol% alcohol into one containing
95mol%. But a significant amount of the alcohol is being lost
in the second product stream. It’s possible to strip the alcohol
from this stream by feeding it to a second cascade.
is recycled into the bottom of the cascade. The vapor stream at the top represents the second
product stream. The analog of the reflux ratio for a stripping cascade is the boil-up ratio
V
RB ≡
B
Because the bottoms product is only partially vaporized to form the vapor stream which is
recycled to the cascade, its composition is not the same as the bottoms product. In other words
xN ≠ yN+1 ≠ xB
So no analogy can be made to (101). However we still have a operating-line relationship like
(99) or (102) which is imposed by mole balances. From a balance over the entire stripping
cascade on the more volatile component, we know that
=
Fx = BxB + Vy1
Lx
F 0
in out
Substituting this result into (99), the operating line for a stripping cascade becomes♠
L B
y=
i +1 xi − xB (105)
V V
The ratio of the flowrates can be expressed in terms of the boil-up ratio:
B B 1
= =
V RB B RB
L V + B RB B + B RB + 1
= = =
V V RB B RB
RB + 1 x
The operating line (105) becomes:=
yi +1 xi − B (106)
RB RB
slope
Like all the operating lines we’ve encountered so far, this one
has a slope of
L RB + 1
slope: = >1
V RB
which is always greater than unity for a stripping cascade. One point on this operating line is
Proof: substituting yi+1 = xi = xB automatically satisfies (106). Having one point and the slope,
we can draw the operating line. Alternatively, a second point is
RB + xB
second point: R + 1 ,1
B
R +1 xB
Proof: Substituting y = 1 into =
(106): 1 B x−
RB RB
RB x RB RB + xB RB + xB
= 1 + B= = R +1
RB + 1 RB RB + 1 RB
then solving for x: x
B
Example #2. Suppose I have a feed stream containing 50% alcohol in water which I want to
strip the alcohol from to produce a stream containing only 5% alcohol. How many stages do I
need if I boil half of the liquid leaving the bottom of the cascade and return it as vapor?
Find: N
Solution: Once again, vapor and liquid streams leaving the same stage are assumed to at
equilibrium; thus (xi,yi) lies on the equilibrium curve. On the other hand, the composition of
adjacent streams (xi,yi+1) lies on an operating line for this cascade. One point on this operating
line is xi = yi+1 = xB. A second point is:
RB + xB 1 + 0.05
= ,1 = ,1 ( 0.525,1)
RB + 1 1 + 1
♣ Clearly the figure at right is not drawn to scale (since the y=1 intercept is located well above 0.5 on the x-axis).
On the other hand, (xN, yN+1) are adjacent streams in the cascade and so their compositions must
lie on the operating line:
Similarly, (xN, yN) lies on the equilibrium curve, while (xN-1, yN) lies on the operating line.
We continue to step off stages in this manner until we step past a liquid composition equal to
the feed xF; in other words, we continue stepping until we cross the blue vertical feed line in
Error! Reference source not found.. In this case, three steps are required (actually only 2.3).
But one of these steps is the reboiler; only two additional equilibrium trays are required:
Comment #1: Notice that the “feed line” for a stripping cascade is vertical, whereas the feed line
for a rectifying cascade was horizontal (see page 114). This is because the feed into the bottom
of a rectifying cascade is a vapor, whereas the feed to a stripping cascade is a liquid.
Comment #2: Suppose that our feed composition had been a little higher, say x'F, which is above
the intersection of the operating and equilibrium curves. No matter how many steps we took, we
could never get a liquid composition higher than that at the intersection of the operating and
equilibrium curves. In that case, any solution for the specified boil-up ratio RB is impossible.
Instead RB would have to be increased, which will rotate the operating line clockwise about
point (xB, xB). This will raise the liquid composition at the intersection point and make the
separation possible. It turns out that, just like the minimum reflux ratio for a rectifying cascade,
there exists a minimum boil-up ratio RB which allows a particular bottoms composition xB to be
produced from a given feed composition xF.
Let's finish the problem by determining the flowrates. Once again, all the flowrates can be
deduced from the equimolal overflow♠ assumption and the boil-up ratio. All the liquid flowrates
correspond to the flowrate of the liquid fed to the top:
L=F
Since I’m boiling half of this liquid to return to the cascade as vapor, all the vapor flowrates are
V = L/2 = F/2
The bottoms flowrate is what remains after half of the liquid to the reboiler is boiled:
B = L – V = F/2
♠ “Equimolal overflow” is another name for the “constant molal overflow” approximation: see page 110.
STEAM DISTILLATION
Instead of using a heat exchanger to partially boil the liquid
stream leaving the bottom stage to produce the vapor entering the
bottom stage, we could directly inject steam as the vapor. This has
the advantage that it avoids the cost of a heat exchanger and you are
getting direct contact between the hot and cold fluids.
Find: N and N
where N and N are the number of stages in enriching and stripping sections.
Solution: Having two cascades is a little trickier than one, so let me first tell you how to solve
the problem by giving you the recipe, and later we will rationalize some of the steps in this
recipe.
This recipe is called the McCabe-Thiele Method and it works when you have only two
components and you can assume equal molal overflow.
We use “zF” to denote the mole fraction of the more volatile component (A) in the feed, ignoring
which phase the feed is in:
HV − H F
q≡
HV − H L
Step 4) Plot q-line, which is a line given by the following equation (QL):
q zF
=y x+ (107)
q −1 1− q
q
slope:
q −1
which is zero (horizontal) for a saturated vapor (q = 0) or infinite (vertical) for a saturated liquid
(q = 1). Otherwise, it’s easier to plot two points. A second point is where the q-line crosses the
y-axis (i.e. at x=0). Substituting x=0 into (107), we calculate the y-coordinate of the second
point:
Should this y-coordinate be larger than one or less than zero, an alternative second point is where
the q-line crosses the y=1 line. Substituting y=1 into (107), we solve for the x-coordinate of the
alternative second point:
Of course, the name "q-line" is given because this line depends on q — the thermal condition
of the feed. For example, suppose the feed is ...
q
superheated vapor: HF>HV → q<0 → 0< <1
q −1
sat'd vapor:
q
HF=HV → q=0 → = 0 (horizontal line)
q −1
2-phase mixture:
q
HL<HF<HV → 0<q<1 → −∞ < <0
q −1
Figure 11: an xy diagram
In this case, q represents the fraction of the feed which is liquid. showing q-lines for feeds having
different thermal conditions (i.e.
q different HF).
sat'd liquid: HF=HL → q=1 → = ∞ (vertical line)
q −1
q
subcooled liquid: HF<HL → q>1 → 1< < +∞
q −1
As it turns out, the two operating lines for the two cascades
intersect somewhere along this line: for a proof, see Role of q-
Line below on page 137.
R x
=y x+ D (108)
R +1
R +1
L
V
xD
another point: 0,
R +1
Step 6) Plot the operating line for the stripping section (SOL):
RB + 1 x
=y x− B (109)
RB RB
L
V
Comment #1: equations (108) and (109) for the two operating lines share a common feature also
found in the generic form of the operating line given by (99): the slope of any operating line is
the ratio of liquid to vapor flowrates in that cascade: L V for the rectifying cascade and L V for
the stripping cascade.
Comment #2: The ROL was plotted in Step 5) using the same two points used in Example #1
(for an isolated rectifying cascade on page 113). But for the SOL, we used only one of the two
points used in Example #2 (for an isolated stripping cascade on page 117): namely (xB, xB).
Instead of using the second point in Example #2, Step 6) of the recipe above used the
intersection of the QL and the ROL.♦ While the second point used in Example #2 is also
correct here, it would have required us to determine the boil-up ratio RB. The procedure
suggested in this recipe avoids this determination.
Note that the boil-up ratio RB and the reflux ratio R cannot both be specified in the problem
statement (we specified only R). Instead these two values are related by the constraint that the
ROL and SOL must intersection along the QL.
You could start at either end (x=y=xD or x=y=xB), although it is customary to start at the top.
♦ A proof that ROL and SOL intersect along the QL is offered below in Role of q-Line on page 136.
Either of the above two stepping schemes is OK, but one scheme leads to fewer total trays than
the other.
optimum feed plate: to get the fewest number of total stages, you switch operating lines when
that allows you to take larger steps
In other words, if you have a choice of which operating line to use, use that one which takes you
furthest from the equilibrium curve. This choice gives the biggest step and fewest total stages.
Rule: an infinite number of steps is required to get to an intersection of the operating line and
equilibrium curve.
Proof: when the size of your steps are no longer large compared to the width of the pen or pencil
you are using to draw them, you should "blow up" the drawing to reduce drawing errors.
Figure 12
To demonstrate the infinity of steps required, let me present the results of stepping off stages in a
particular example:
Example: given the xy-diagram for acetone and methanol at 1 atm pressure (has an azeotrope at
x = y = 0.8000), step off stages for the following conditions:
Solution: This is a particularly difficult separation because the desired distillation composition
xD = 0.7999 is very close to the azeotrope. We anticipate a large number of steps. Available
online is a Mathcad document♣ which steps off stages automatically using the optimum feed
location.
The results are shown in the figure above. The program reports the total number of steps and the
feed location. In this figure, there are a total of 64 steps with the feed introduced on the 60th
from the top. In the upper right corner of the xy-diagram, the steps are becoming so small as to
be invisible.
Below on the left is a blow-up of the region from 0.7 to 0.8. Even on this blow-up the steps
are becoming invisible in the upper right corner. So the second graph below on the right is
another blow-up: this from 0.79 to 0.80. Still the steps are becoming vanishingly small. Notice
that in these blow-ups, both the equilibrium curve and the operating line are very linear.
♣ McCabe-Thiele.mcdx can be found in the “Student Resources” folder under the “Handouts” heading.
0.78 0.798
0.76 0.796
0.74 0.794
0.72 0.792
0.7 0.79
0.7 0.72 0.74 0.76 0.78 0.79 0.792 0.794 0.796 0.798
log
( ya* − ya )
N=
( yb* − yb )
(110)
log
( ya* − yb* )
( ya − yb )
which is called Kremser's equation.♥ The subscripts a and b denote streams at the top and
bottom of the cascade, respectively (i.e. “b” is for “bottom”). The y*s are vapor compositions in
equilibrium with the liquid at the top and bottom of the cascade.
yi = mxi + b (111)
On the other hand, (xi, yi+1) is a point on the OL, which is always linear (assuming equimolal
overflow) [recall (99) on page 111]:
L L
yi +1 = xi + ya − xa
V V
where we have substituted y1=ya and x0=xa. Using (111) to eliminate xi the OL equation
becomes:
L yi − b L
=
yi +1 + ya − xa (112)
m
V
V
mA mA
Now we define a new parameter, which represents the ratio of the slope of the operating line
(L/V) to the slope of the equilibrium curve (m):
L
A≡
mV
From (111), we see that mxa+b is just the vapor concentration in equilibrium with liquid at
concentration xa, which we denote at y*a. Making this substitution, the above equation becomes
yi +1 =Ayi − Ay *a + ya (113)
This serves as a recursion formula, which we can apply repeatedly to obtain the concentration on
successive stages, starting at the top. For i=1, recalling y1 = ya, (113) yields
i = 1: y2 = A y1 − Ay *a + ya = ya (1 + A ) − Ay *a
ya
y3 =Ay2 − Ay *a + ya =A ya (1 + A ) − Ay *a − Ay *a + ya
( ) ( )
i = 2:
= ya 1 + A + A2 − y *a A + A2
where the second equality above results from substituting y2 from the previous step. Now we
begin to see a pattern developing:
i = n: ( ) (
yn +1= ya 1 + A + A2 + + An − y *a A + A2 + + An
) (114)
A (1+ A + A2 + + An −1 )
Each sum is a partial sum of a geometric series, whose total is known analytically:♣
1 − An + 1
1 + A + A2 + + An =
1− A
1 − An + 1 1 − An
=yn + 1 ya − y *a A
1− A 1− A
For the bottom stage of the cascade, n=N and yN+1 = yb:
1 − A N +1 1 − AN
=yb ya − y *a A
1− A 1− A
yb (1 − A )= ya (1 − A N +1 ) − y *a A (1 − A N )
yb − Ayb =ya − A N +1 ya − Ay *a + A N +1 y *a
N +1
yb − ya + A ( y *=
a − yb ) A ( y *a − ya ) (115)
This can be further simplified. Recall that A is defined as the ratio of the slopes of the operating
line to the equilibrium curve. This ratio can be calculated from the concentration defined on the
xy-diagram above:
( ya − yb )
slope of OL ( xa − xb ) ya − yb
A≡ = = (116)
slope of EC ( y *a − y *b ) y *a − y *b
( xa − xb )
We can use this result to obtain ya – yb = A(y*a – y*b), which will now be substituted into (115):
♣ see (2.6) in M.D. Greenberg, Foundations of Applied Mathematics, Prentice-Hall, 1978. Available electronically
since 2013 as a Dover Edition.
A ( y *b − y *a ) + A ( y=
*a − yb ) A N +1 ( y *a − ya )
( )
A y *b − y *a + y *a − yb= A ( y *b − yb )
y * −y
log b b
AN =
y *b − yb
or N= y *a − ya
y *a − ya log A
Taking the log of both sides, recalling that log(AN) = N logA, then solving for N we obtain the
second expression above. To obtain the form of Kremser's equation given as (110), we substitute
(116) for A.
Comment #1: In the xy-diagram shown next to (110), the equilibrium curve lies above the
operating line. This is always the arrangement when we are dealing with gas absorption
(transport from gas to liquid). In order to obtain N>0, the following inequalities must be obeyed:
If you inadvertently (but consistently) interchange the subscripts a and b, you will obtain a
negative value for N.
Comment #2: If the operating line and equilibrium curve should cross (which is physically
impossible), when evaluating (110) you will find yourself trying evaluate the log of a negative
number.
Comment #3: The situation in which the operating line and equilibrium curve are parallel (i.e.
A = 1) is a special case: (110) becomes indeterminant. However, since the driving force y*-y is
constant, the number of steps is just the total change in concentration ya-yb divided by the height
of one step y*a-ya:
ya − yb
for A = 1: N=
y *a − ya
Now we will return to discuss the motivation for Kremser's equation: approaching the
intersection of OL and EC:
Notice that, as xb approaches the point of intersection (see Figure 12 on page 124), yb* → yb
and
( )
log
= 0 log [ ∞ ]
=N →∞
log [ ] log [ ]
Thus you can never reach the point of intersection for any finite number of steps. In particular, a
very large number of stages is required as we increase the purity of either of our two product
streams. For example, as xD→1, the end-point (xD, xD) of our operating line approaches (1,1),
which is always one point on the equilibrium curve: as xD→1, N→∞.
Example: in the previous example in which xD is close to the azeotrope, both the OL and EC
appear to be straight in the interval 0.7900 ≤ x ≤ 0.7999 (see right-hand figure on page 126). Use
this fact and Kremser’s equation to calculate the number of stages in this interval.
Solution: we need to calculate the four y’s which appear in (110). Two of these y’s are given by
the operating line. One of the points on the operating line is
xa = ya = xD = 0.7999
R 1
yb = xb + xD =0.791029
R +1 R +1
Substituting xb = 0.79 and R = 8.618, we obtain the value given in the second equation above.
To evaluate points on the equilibrium curve, we linearly interpolate two points at each end of this
interval: (0.790, 0.792146)♠ and (0.800, 0.800).
0.800 − 0.792146
y * −0.792146
= ( x − 0.790 )
0.800 − 0.790
For xa = 0.7999, we obtain ya* = 0.799921 and for xb = 0.79, we obtain yb* = 0.792146.
Substituting these numbers into Kremser’s equation, we obtain N = 30.1.
♠ The value of y* = 0.792146 was obtained using the cubic-spline interpolation [i.e. yeq(0.79)] of the tabular data
for this methanol-acetone system found in McCabe-Thiele.mcdx, except one point was omitted from the table: (1,
1). Values reported above were calculated by McCabe-Thiele near pinch (in-class example).mcd. Omission of this
point – which is outside the interval of interest – alters slightly the result of the interpolation: including (1, 1) gives
yeq(0.79) = 0.791899; excluding (1, 1) gives yeq(0.79) = 0.792146.
y* − y
log
a a ( )
log 0.799921 − 0.799900
=N =
yb
*
− yb
( )
0.792146 − 0.791029
log
y* − y*
a b ( )
0.799921 − 0.792146
log
0.799900 − 0.791029
( y − y )
a b
0.000021
log
= = 0.001117 log ( 0.0019 ) −1.72
= = 30.1
0.007775 log ( 0.876 ) −0.057
log
0.008871
Stepping off the stages with Mathcad, we require 33 steps to get the liquid mole fraction below
0.79.
1) Existence of Rmin
R
slope of ROL=
R +1
This point of tangency is called a pinch point. In stepping off the stages, it will take an
infinite number just to reach the pinch, but if the reflux ratio is increased ever so slightly, then
the separation can be achieved, although the number of stages might be large.
minimum reflux ratio -- smallest value of R for which the desired separation can be achieved
with a finite number of stages.
The pinch point might also occur on the SOL as shown at right.
Regardless of where the pinch occurs, once you have drawn the q-
line and operating lines, the minimum reflux ratio (thus determined
graphically) can be calculated from:
x − y′
R= D (117)
y ′ − x′
where (x', y' ) is any point on the ROL except (xD, xD). For
maximum accuracy in a graphical solution, choose a point as far as
possible from (xD, xD); usually this would be the y-intercept of the ROL.
Proof: (117) was derived by setting xi = x' and yi+1 = y' in (104) on page 112:
R x
=
yi +1 xi + D (104)
R +1 R +1
y′ x′
and then solving for R. The point (x', y' ) can be chosen for convenience: for example, you might
want to choose the y-intercept.
An important thing to notice about the minimum reflux ratio is that, as we approach it, the
number of stages approaches infinity:
as R → Rmin, N → ∞
Clearly you don't want to operate near the minimum reflux ratio because it would cost too much
to build a column with so many stages. Generally as you increase the reflux ratio, both operating
lines move away from the equilibrium curve so that the steps become larger. Then fewer stages
will be needed.
as R↑, N↓
In the limit that R approaches infinity, the y-intercept of the ROL becomes zero. The ROL
becomes the 45° line. Of course, the SOL and ROL always intersect along the q-line, which
means that the SOL also becomes the 45° line. This is the furthest the operating lines can ever
be from the equilibrium curve, so this condition gives the minimum number of stages:
and N → Nmin
Although it might appear that this is the best value of the reflux ratio, it turns out that there is a
"dark side" to this solution. The problem is what is happening to flowrates inside the column.
Recall R ≡ L/D
Now to get infinite reflux ratio, one of two things has to happen:
R → ∞: L→∞ or D→0
In most designs, you are interested in processing a fixed amount of material, this implies that
F,D, and B are fixed, which eliminates the second alternative above. If V = (R+1)D, L = RD and
D remains fixed as R→∞, we see that both L→∞ and V→∞:
L → ∞: V → ∞, qc → ∞, qr → ∞, tower diameter → ∞
♣ The term “fixed costs” is used in accounting because the depreciation of construction costs is the same every
month regardless of how the equipment is operated. In contrast, the cost of utilities (electricity, steam, cooling
water) are termed “variable costs” because they vary month by month, depending on how the equipment is used.
reflux (where the diameter of the column becomes large). The costs of utilities increase about
linearly with reflux ratio.
The total costs will experience a minimum at some reflux ratio, which we call the
However most chemical plants actually operate at slightly larger R since the additional cost is not
very much and it buys you flexibility in case you decide later you want to change xD or in case
the feed composition changes, you have some room to maneuver:
Find: qr and qc
F=D+B
zFF = xDD + xBB
These are the same pair of equations we solved earlier for flash distillation [see (80) and (81) on
page 93]. The solution is given by (89) and (90):
x − xB x − xF
D= F F and B= D F (118)
xD − xB xD − xB
Now we’ll determine the internal flows by starting at the top of the column and work our
way down. First we use the definition of reflux ratio:
In a total condenser, we are condensing all of the overhead vapor. Assuming the condenser
output is a saturated liquid, the condenser heat duty is just
qc = V λ
x − xB
qc= ( R + 1) Dλ= ( R + 1) F Fλ
xD − xB
where in the second and third equations, we have substituted (119) and (118), respectively.
V +F +L=V +L (120)
HL = 0, HV = λ
Recall the definition of q from page 121. This allows us to express the enthalpy of the feed in
terms of q:
HV − H F λ − H F
q≡ = or HF = (1 – q)λ
HV − H L λ
Substituting these values for enthalpy into (121) and solving for V gives:
F(1-q)λ + V λ = Vλ or V = V − (1 − q ) F (122)
L= L + qF (123)
Thus introducing the feed onto a plate causes a discontinuity in the flowrates of liquid and vapor
in the column. Although the flowrates are the same for each stage in the cascade above the feed,
and the same for each stage in the cascade below the feed, the flowrates are different in the two
cascades.
For example, for a saturated liquid feed, q=1, the liquid leaving the feed tray has a higher
flowrate than the liquid entering by an amount equal to the feed flow rate:
On the hand, if the feed is saturated vapor, then the vapor flowrate leaving
the stage will be higher than entering by an amount equal to the feed rate.
The flows for other thermal conditions of the feed are summarized graphically
below:
When the feed is a superheated vapor (q<0), not only does all of its flow join the vapor stream,
but it also boils up some of the liquid. When the feed is a two-phase mixture (0<q<1), then the
liquid fraction joins the liquid stream and the vapor fraction joins the vapor stream. When the
feed is a subcooled liquid (q>1), all of its flow enters the liquid and some of the vapor is also
condensed and joins the liquid stream too.
In any case, given F, L and V, we can calculate L and V from (122) and (123). Knowing
all the flowrates, we can now calculate the reboiler heat duty. The reboiler has to boil all the
vapor sent into the bottom cascade. For a partial reboiler (where the bottoms product is
withdrawn from the liquid remaining in the kettle), there can be no superheating of the vapor (it
must be a saturated vapor). Then the reboiler heat duty is
qr = V λ
x − xB x − xB
q=
r V − (1 − q ) F =
λ ( R + 1) F
xD − xB
F − (1 − q ) F = λ ( R + 1) F
xD − xB
+ q − 1 F λ
V
ROLE OF Q-LINE
Recall (see “Step 4” on p121) that we said one point of the
SOL is the intersection of the q-line (QL) and the ROL. In this
section, we present the arguments supporting this claim.
We claim that the q-line represents the intersection of these two lines. Think of these as two
equations in two unknowns (x, y). The point of intersection of these two lines is that point (x, y)
which simultaneously satisfies both equations.
Now we are not going to try to solve these two equations for (x, y). But we can more easily
find another linear relationships between x and y by combining these two equations. One linear
combination is particularly convenient in that it yields an equation which we can plot without
knowing any of the flowrates that appear in (124) and (125):
Notice that the last term in this equation represents the total rate of
component A leaving the tower. This must equal the total rate at
which A enters the tower. So the above becomes:
Dividing by F gives:
V −V L−L
= y x − zF
F F
q −1 q
q zF
or =y x+
q −1 1− q
which is the equation of our q-line. This equation represents a linear combination of the two
equations which are satisfied at the intersection of the two operating lines. Thus this equation
must also be satisfied at that point. It turns out that the QL represents the locus of intersections
of the ROL and SOL for different reflux ratios.
PLATE EFFICIENCY
So far we have always assumed that the vapor and liquid
streams leaving a given stage are at equilibrium. In stepping off the
stages, we have always gone all the way from the operating line to
the equilibrium curve. If the contacting of the two phases on the
tray is good (that is, if the area of contact is large enough and the
time of contact is long enough) we can approach this ideal state, but
we will never quite reach it. For a non-equilibrium stage, our
horizontal step will be a little shorter than for an ideal stage. Of
course, we still end up at a point on the operating line, because this
represents a mole balance which must still be satisfied.
Because the real steps are shorter, it will take more of them to achieve the same separation.
The ratio of the number of ideal to real stages required will be less than one and is called the:
While this efficiency is easy to use, the value of the efficiency is hard to
determine. A more fundamental quantity is the
See equilibrium curve and operating line sketched at right. Basically the definition is the size of
a real step from the operating line toward the equilibrium curve divided by the size of an ideal
step (which would take us all the way to the equilibrium curve).
The corresponding definition for the efficiency based on the vapor mole fractions is
This definition is based on the height of a step (involves changes in y) rather than based on the
width of a step (changes in x). Unless the equilibrium curve and operating lines are parallel,
these two definitions are not equivalent, but the differences are probably within experimental
error.
Some typical values of plate efficiencies are given in Table 14-12 (taken from Perry’s
Chemical Engineering Handbook, 8th edition, 2008).♦
We locate the two operating lines and the q-line, just as before. But before we step off the stages
we use the efficiency to determine a new curve between the operating lines and the equilibrium
curve.
curve depends both on the vapor-liquid equilibrium and on the operating lines. Thus if we
change the reflux ratio, for example, we will have to replot the pseudo-equilibrium curve.
Nc components
Nc-1 columns
This assumes that each column can achieve a sharp separation between two components of
adjacent volatility. “Alternative #1” represents one scheme for separating a 3-component feed
into 3 nearly pure products.
SHARP SPLIT
The first step in the design of a single column is to specify the product distribution. This can
be estimated easily by assuming the column achieves a “sharp split” between two components.
To some extent, the choice of which components to split is up to the designer who can tune the
reflux ratio and feed location to obtain different sharp splits.
Comment #2: Sharp separations are not always possible, but we will focus on them in this
introduction, for the sake of simplicity.
In the description above, components M and N play a special role and are called
The more volatile of the two is called the “light key,” while the other is the “heavy key”:♣
light key = M
heavy key = N
given: zF6 = 0.33, zF7 = 0.37, zF8 = 0.30, xD7 = 0.01, and xB6 = 0.01
find: D/F, B/F, and xDi, xBi for the other components
Solution:
♣ Here “light” and “heavy” suggest weight, as in molecular weight: in a homologous series, lower molecular weight
compounds are more volatile so the light key is more volatile than the heavy key.
If we had a more complex mixture (i.e. not all components members of the same homologous
series), then we would look up the boiling points of the pure components at the pressure in the
distillation column. The ranking then is from lowest boiling point (i.e. most volatile) to highest
boiling point (least volatile).
Let's take a basis of 100 moles of feed (F=100) and try to fill out the following table:
Component FzFi zFi DxDi xDi BxBi xBi
hexane 33 0.33 -- -- -- 0.01 light key
heptane 37 0.37 -- 0.01 -- -- heavy key
octane 30 0.30 -- -- -- --
totals = 100 1 D 1 B 1
1. Components heavier than the HK are assumed to end up in entirely in the bottoms
2. Components lighter than the LK are assumed to end up in entirely in the distillate
In essence, we are then making a sharp split between hexane and heptane with the “1%” values
specifying the sharpness of the split desired. In other words, hexane has been chosen as the light
key and heptane as the heavy key. If we assume that there is virtually no octane in the distillate
(octane is heavier than the heavy key), then we can finish specifying the composition of the
distillate stream. In this case, this means 1) that the mole fraction of octane in the distillate is
zero and 2) the molar flowrate of octane in the bottoms is equal to that of the feed: 30. Finally,
to have the sum of the mole fractions of the distillate equal unity, the mole fraction of hexane can
be deduced to be 1 – 0.01 = 0.99:
Knowing the concentration of hexane in both of the product streams allows us to determine the
flowrates from mole balances:
total: D + B = 100
Now using the fact that the balance of the distillate and bottoms streams must be heptane, we can
fill out the table:
Comment #2: The only additional information that is needed if the number of components is
larger than 3 is the mole fraction of those components in the feed. Although there will be more
rows in the table (if there are more components), the additional unknown mole fractions in the
distillate and bottoms streams equals additional constraints obtained by applying the “rules.”
Example #2: A mixture with 25% butane, 25% pentane, 25% hexane and 25% heptane is to be
distilled to recover 95% of the butane in the distillate and 95% of the pentane in the bottoms.
Estimate the remaining compositions of the distillate and bottoms.
Given: zF4 = zF5 = zF6 = zF7 = 0.25, DxD4 = 0.95×FzF4 and BxB5 = 0.95×FzF5
Find: D/F, B/F, and xDi, xBi for the other components
Solution:
Recall that 95% recovery refers to the component flowrate in the distillate or bottoms as a
fraction of the component flowrate in the feed:
Let's take a basis of 100 moles of feed (F=100) and try to fill out the following table:
Components heavier than the heavy key end up entirely in the bottoms. This gives us the
flowrate of the hexane and pentane components.
Knowing the component flowrate of any component in one of the product streams allows us to
determine that components flowrate in the other product stream:
Knowing the flowrate of each component in each stream, we can divide by the total flowrate to
component the mole fraction of each component in each stream:
Comment: Notice that the extra unknowns introduced by an additional component were matched
by an equal number of extra constraints arising from applying the “rules.”
x x
ln Di Bi
xDj xBj
N min =
ln αij ( )
To obtain such a simple analytical result requires some gross assumptions:
In case you have forgotten, the definition of relative volatility can be found in the next
section and the “proof” that it is nearly constant for ideal solutions of components having the
same molar heat of vaporization can be found on page 83.
Proof:
Ki
Recall αij ≡ (128)
Kj
y
Ki ≡ i (129)
xi
♥ We have one rule for each component lighter than the light key and heavier than the heavy key, which is one real
for each component less the two keys.
is the distribution coefficient for component i (or j). With these assumptions, we can start at the
top of the column and work out the compositions of streams leaving every stage.
y1i
Ki x
αij ≡ = 1i
K j y1 j
x1 j
x1i 1 y1i y
or = = α ji 1i (131)
x1 j αij y1 j y1 j
x1i x
Substituting (130): = α ji Di (132)
x1 j xDj
x1i x
for i=j: = αii Di or 1=1
x1i xDi
1
1 1
which doesn’t tell us anything new; thus we have only Nc-1 independent relations of this type
instead of Nc. Choosing one element as the reference (say j=Nc), we can write Nc-1 relations for
the Nc-1 independent mole fractions. The missing equation is Σx1i = 1. Thus the composition of
the liquid leaving tray 1 has been completely determined.
Next we determine y2i by mole balances (system is enclosed by red rectangle in figure
above).
for i=1,...,Nc.
Knowing the composition of the vapor leaving stage 2, we can calculate the composition of
the liquid leaving stage 2 from the known volatilities. Rewriting (131), replacing the stage
number 1 by 2:
x2i y
= α ji 2i
x2 j y2 j
x2i x
From (133): = α ji 1i
x2 j x1 j
x2i x
Using (132): = α 2ji Di (134)
x2 j xDj
From (132) and (134), we begin to see a trend. The generalization of this trend yields:
xni x
= α nji Di
xnj xDj
In particular, we are interested in the number of steps required to get from the distillate to the
bottoms:
xBi N x
= α jimin Di
xBj xDj
The exponent represents the minimum number of steps on the xy diagram. One of these steps is
the partial reboiler. Solving for Nmin:
x x
ln Di Bi
xDj xBj
N min =
( )
ln αij
Now in any real system, the relative volatility will not be constant. One way to patch things
up is to use some kind of average volatility. The usual choice is the geometric mean of the
volatilities evaluated at the two extreme temperatures in the column:
Ki
αij ≡ (135)
Kj
y
where Ki ≡ i (136)
xi
is the distribution coefficient for component i (or j). The latter quantity describes how a given
component distributes itself between vapor (y) and liquid phases (x). For mixtures of
components which form ideal solutions, Raoult's law applies:
pi = xi Pio (T ) (137)
where Pio (T ) is the vapor pressure of pure component i at the same temperature as the mixture.
Recall that partial pressure of a component is defined as the mole fraction times the total
pressure:
where P is the total pressure, which also equals the sum (over all Nc components) of the partial
pressures:
Nc
P = ∑ pi
i =1
Substituting xi from (137) and yi from (138) into the definition of partition coefficient (136), we
have
pi
yi P Pio (T )
Ki ≡ = = (139)
xi pi P
Pio
Ki P o (T )
αij= = i (140)
K j Pjo (T )
Using these constants, we then can calculate the vapor pressure at any temperature from:
Bi
log10 Pio (T )=
A −
i (141)
T + Ci
Recall that, as temperature drops, the “dew point” is the temperature at which the first drop
of dew is formed in the vapor. Usually, the composition of the vapor is specified as yi for i = 1, 2
… Nc. At the dew point, the vapor having composition yi is in equilibrium with that first drop of
dew having composition xi. Knowing yi, we can calculate xi at any T (for which two phases
coexist) from (136):
yi
xi =
Ki (T )
Nc
y
∑ K (Ti ) = 1
i =1 i dew
This is a trial-and-error for T. For a particular guess of T, the K-values can either be read from a
graph (like Fig. 11.7-2 in Geankoplis) or they can be calculated from vapor pressure data using
(139). Substituting Ki from (139) into the equation above and dividing through by P:
Nc
yi 1
∑ =
o
(
i =1 Pi Tdew ) P
To evaluate the dew point of the distillate, we choose yi = xDi and P = 1.2 atm. The dew
point is found to be 75.2°C [see Bubble point and dew point of multicomponent mix.xmcd].
This Mathcad document uses vapor pressures calculated from (141). At this temperature, the
vapor pressures are evaluated and the relative volatility is calculated from (140). The results are
summarized by
As temperature rises, the bubble point is the temperature at which the first bubble of vapor is
formed. Usually, the composition of the liquid is specified as xi for i = 1, 2 … Nc. Knowing xi,
we can calculate yi at any T from (136):
yi = Ki (T ) xi
Nc
∑ Ki (Tbub ) xi = 1
i =1
Multiplying this equation through by P and substituting (139), we obtain an alternate equation:
Nc
∑ Pio (Tbub ) xi = P
i =1
Once again, solving either equation for Tbub is a trial-and-error. In this example, the bubble
point of the bottoms is calculated taking xi = xBi and P = 1.2 atm [see Bubble point and dew
point of multicomponent mix.xmcd]. The result is 114.2°C and the results are summarized by
the following table:
To calculate Nmin, Fenske’s equation is applied to the heavy key and the light key. So we are
interested in αLK,HK
α=
67 2.547 × 2.213
= 2.374
x x
ln Di Bi ln 0.99 0.01
xDj xBj
=N min = =0.01 0.544
9.82
ln αij ( )
ln ( 2.374 )
To the extent that the assumptions made during its derivation are valid, Fenske’s equation
should yield the same value for Nmin using any choice of component i and j. Once the value of
Nmin is known, we can choose an i or j corresponding to a component whose mole fraction was
assumed to be zero for the purpose of estimating the product distribution. Then Fenske’s
equation can be used to calculate the actual mole fraction of that component as a check to see
whether the mole fraction is indeed small enough to be neglected.
For example, in Example #1, we assumed that the octane fraction in the distillate is
negligible. Let’s now check to see what the octane mole fraction actually is. We use i=8
(octane) in Fenske’s equation rather then i=6 (hexane, the light key):
α=
87 0.401 × 0.463
= 0.431
x 0.445
ln D8
= 9.82
N min = 0.01 0.545
ln ( 0.431)
which indeed is negligible compared to the other mole fractions in the distillate.
Assuming a saturated liquid feed (q=1) and constant relative volatility, the minimum reflux
ratio turns out to be
xD 1 − xD
− α AB
z 1 − zF
Rmin = F (142)
α AB − 1
where xD and zF are the mole fraction of the light key in the pseudo-binary mixture and 1-xD and
1-zF are the mole fraction of the heavy key in the pseudo-binary mixture:
z F , LK z F , HK
zF = and 1 − z F =
z F , LK + z F , HK z F , LK + z F , HK
where zF,LK and zF,HK are the mole fractions in the multi-
component mixture.*
Proof: We can obtain the minimum reflux ratio from the point of
intersection of the q-line and the equilibrium curve. Recall (117)
from page 132:
x − y′
Rmin = D (143)
y ′ − x′
which was derived using two points on the ROL to determine the
slope and knowing that the slope is R/(R+1). For a saturated
liquid feed, the q-line is vertical so that
x′ = zF (144)
The y-coordinate of this point is determined knowing that ( x ′, y ′ ) lies on the equilibrium curve.
For constant relative volatility, we can express the equilibrium curve as a simple equation using
the definition of relative volatility:
K A y A x A y (1 − y )
α AB ≡ = =
K B y B xB x (1 − x )
αx
Solving for y: y=
1 + ( α − 1) x
α AB z F
y′ = (145)
1 + ( α AB − 1) z F
This method makes assumptions similar to those used in Fenske’s equation. In particular, the
assumptions made by Underwood’s method are:
1) equimolal overflow
2) constant relative volatility
3) a sharp split with only one “pinch”
First, we will just state the method and then we will take a quick look at its proof.
Recipe:
Step 1) Find the root (denoted φ) in the interval α HK < φ < α LK which satisfies the following
polynomial (order of polynomial = number of components)
Nc
αz
∑ α i −Fiφ = 1− q (146)
i =1 i
Step 2) Once we have φ, we can compute the minimum reflux ratio from
αx
Rmin + 1 =∑ αi −Diφ (147)
i i
where the sum is just over those components present in the distillate (i.e. we ignore components
heavier than the heavy key).
yn , i
Ki = yi/xi or K n, i =
xn, i
yn , i
Vn +1=
yn +1, i Ln + DxD, i
K n, i
For constant molal overflow, the flowrates will be the same for each
tray in the rectifying cascade, just as with binary mixtures. Thus we can drop the subscripts on
flowrate.
yn , i
Vy=
n +1, i L + DxD, i (149)
K n, i
When we approach the minimum reflux ratio, we will encounter a "pinch" somewhere in the
column. At the pinch, we take an infinite number of steps without changing the mole fractions:
D xDi
yi = (151)
V∞ L∞
1−
V∞ K ∞, i
We have added the subscript “∞” to each flow rate to emphasize that this corresponds to the
minimum reflux, which requires an infinite number of trays. If we require (151) to be satisfied
for each component and sum over all components:
D xDi
∑ yi = 1= ∑
V∞ i L∞
(152)
i 1−
V∞ K ∞, i
K
αi = i
K ref
While the flowrate ratios in the rectifying section are related to the reflux ratio by (103):
L∞ Rmin D L∞ R 1
= and= =
V∞ Rmin + 1 V∞ V∞ Rmin + 1
αx
(152) becomes: ∑ αi=
Di
−φ
Rmin + 1 (147)
i i
Rmin
where φ≡
K ref ( Rmin + 1)
If Kref were known, then Rmin could be determined as the root of this equation. Unfortunately, to
calculate Kref we need to know the temperature at the pinch, which is not known.
Generally there is also a pinch in the stripping section of the tower. For sharp splits, the
temperature at this second pinch will be the same. We can then derive a second equation like
(147) by performing our balances about the bottom n stages of the column to obtain a second
equation.
B xBi
∑ yi = 1= ∑
V∞ i L∞
i −1
V ∞ K ∞, i
where the overbars denotes flowrates in the stripping section. Then we will have two equations
and two unknowns: Kref and Rmin. Their solution is the value of φ which satisfies the following
equation:
αi z F , i
∑α = 1− q
i i −φ
where zF,i is the composition of the feed and q is its thermal condition. There are several φ's
which satisfy this equation. It turns out that we want:
Once φ is determined, we can calculate Rmin from (147). This is called Underwood's method.
Example #4. Calculate the Rmin for the separation in Example #1 if the feed is 60% vapor (q =
0.4).
α z
∑ α i −Fiφ = 1− q
i i
which is a cubic equation. The left-hand side of this equation [defined as f (φ)] is plotted below.
This equation has three singular points at which the function diverges to ±∞:
φ = α1, α2 and α3
These are shown as the 3 vertical lines in the plot. Notice that there are also 3 different φ’s for
which
f (φ) = 1 – q = 0.6
which turns to be φ = 1.739. Once we have the value of φ, we calculate the reflux ratio from
(147):
αx
Rmin + 1 =∑ αi −Diφ
i i
where the sum only includes those components present in significant quantities in the distillate.
In this case, octane (which is heavier than the heavy key) was neglected in the distillate for
Example #1. It should also be neglected in this sum, although in this case, it’s inclusion does not
change the result.
Underwood's method can be extended to calculate the number of plates analytically when the
equal molal overflow assumption and constant relativity are appropriate. But even then, the
calculations are sufficiently involved that a computer would be helpful. If we are going to use a
computer, we might as well do the rigorous tray-by-tray calculations (see Method 3); thus we
will not extend Underwood’s method further in these Notes.
An empirical but much easier to use method is the Erbar-Maddox correlation. This is an
empirical relation between the number of ideal stages and the operating reflux ratio. You also
need to know the minimum number of ideal stages (at total reflux) and the minimum reflux ratio.
Because of its simplicity, this correlation is widely used for preliminary estimates.
R N R
= f min , min
R +1 N Rmin + 1
The function is presented in Geankoplis at Fig. 11.7-3 (except subscript “min” has been replaced
by “m”).
N − N min R − Rmin
= f
N +1 R +1
lim N →∞
R → Rmin
FEED-PLATE LOCATION
The Kirkbride equation provide an approximate method for estimating the feed plate
location once the total number of trays Ne+Ns has been estimated for the operating reflux ratio:
2 0.206
N e z F , HK B xB, LK
= × × (153)
N s z F , LK D xD, HK
where Ne is the number of ideal stages in the enriching section (above the feed) and Ns is the
number of ideal stages in the stripping section (below the feed).
In previous section, we presented a method for estimating the total number N of stages,
which also equals
N = Ne + Ns (154)
Thus (153) and (154) represent two equations in two unknowns: Ne and Ns. Of course, Ne
represents the number of the feed tray (counting from the top of the column).
• M Equations are Mole balances for each component i (Nc equations for each stage):
Note that the first subscript is the component i and the second subscript for the stage n.
• E Equations describe the liquid-vapor Equilibrium relations among the two outlet
streams, which are assumed to be at equilibrium (Nc equations for each stage):
• S Equations are the mole fraction Summations (2 equations for each stage):
• H Equations is the energy balance (H stands for enthalpy) for each stage (1 equation for
each stage):
Ln −1 H L, n −1 + Vn +1 HV , n +1 + Fn H F , n = ( Ln + U n ) H L, n + (Vn + Wn ) HV , n + Qn
If there are a total of Ns stages in the column and Nc components, then the number of equations
is
where Tn is the temperature of the liquid and vapor streams leaving stage n.
Of course, there are more variables involved in the equations. We assume that the feeds and any
side draw-off flows are specified:
Also, we assume that the thermodynamic relations exists so that we can compute the enthalpy of
saturated liquid and saturated vapor, given the composition, pressure and temperature:
known functions: (
H L, n = H L, n Tn , Pn , x1, n xNc , n)
HV , n = HV , n (Tn , Pn , y1, n y N , n )
c
In this section we will use the McCabe-Thiele* solution for the number of stages required to
meet a specified distillate and bottoms compositions and solve the MESH equations to see if the
specified distillate and bottoms compositions can actually be obtained using the number of trays
predicted.
♣ Actually, the problem whose solution is presented here is slightly different from Hwk6 Prob 3. The
thermodynamic model used by Chemsep produces a slightly different equilibrium curve. In particular, the azeotrope
occurs at a mole fraction of 0.76 instead of 0.80. Thus the distillate composition of 0.78 specified in the original
problem cannot be achieved. So we adjusted the distillate composition to 0.715 and have imported Chemsep’s VLE
data into Mathcad.
♦ Our solution is obtained using CHEMSEP, a process simulator available through the CACHE of AIChE. For any
student who is interested, this is available at www.chemsep.org.
The discontinuity of flowrates on the 17th stage occurs because this is the feed tray. Because the
feed is 30% liquid, the liquid flowrate in the bottom cascade is higher than in the top; because
the feed is 70% vapor, the vapor flowrate in the top cascade is higher than in the bottom.
Beyond the expected changes in flowrates at the feed tray, there is a continuous drop in flowrate
of both streams as we move down either the rectifying cascade or the stripping cascade. These
changes in flowrates within a cascade are NOT consistent with the equimolal overflow
assumption.
These variations in flowrate within a cascade stem from the fact that the enthalpy of saturated
liquid is not independent of temperature or acetone content as is required for "equimolal
overflow."
How important are these flowrate differences? On the McCabe-Thiele diagram, these
differences show up as curvature in the operating lines. In the graph below we compare the
slopes of the operating lines.
While the variations in slope look significant, actually plotting the operating lines themselves
shows the difference is not too great:
The plot above shows points (the red diamonds) taken from the operating lines in the
rigorous ChemSep separation and compares them with the ROL and SOL calculated using the
McCabe-Thiele method. These calculated operating lines are the straight, solid lines on the plot.
Connecting the red diamonds with a smooth curve would yield the actually operating lines.
There is very little difference; indeed all the points appear to lie on one of the calculated
operating lines.
The main prediction for this problem is what the distillate composition and flowrate will be
xD D/F
Model
McCabe 0.715 0.3404
Chemsep 0.7131 0.3414
The difference between the two is insignificant. This insignificant difference is typically the
case when the molar heats of vaporization differ by less than 20%.
The main error made by the equimolar overflow assumption is estimating the heat duties
which are proportional to the vapor flowrates. Recalling the flowrates graph on page 166, we
find (for a feed flowrate of 100):
Model V V
McCabe 327.5 257.5
Chemsep 328.3 210.5
Chemsep predicts about the same flowrate through the condenser, but predicts a reboil rate about
18% below that obtained using the equimolar overflow assumption.
Also shown in the table above are the average relative volatilities. The values shown were
determined from the thermodynamic model of Chemsep, rather than the vapor pressure
correlations presented in the problem statement (Hwk #8, Prob. #1), in order to provide a fair
comparison between Fenske’s equation, Underwood’s method and the more rigorous predictions
of Chemsep. However as a consequence of using a different thermodynamic model here, the
remaining results will differ slightly from your solution to the homework problem.
Using Fenske’s equation, we predict Nmin = 6.42. Using this value of Nmin we can then use
Fenske again to check on the two ‘zeros’ of mole fractions in the table. The results are
Rmin = 0.403
R = 0.524
Chemsep was supplied with most of this info and asked to calculate the product compositions
for a reflux ratio of 0.524 and a recovery of 95 of the pentane in the bottoms. The main results
are compared below with our predictions:
Distillate Bottoms
Component Hwk7 Chemsep Hwk7 Chemsep
butane 0.6178 0.6183 1.2×10-4 4×10-5
pentane 0.3668 0.3663 0.0355 0.0370
hexane 0.0154 0.0155 0.5390 0.5382
heptane 5×10-5 1×10-5 0.4255 0.4248
total flowrate 64.75 64.69 35.25 35.31
The agreement is very good. Let’s take a closer look at some of the assumptions. A major
assumption in both Fenske’s equation and Underwood’s method is equimolar overflow. This is
an assumption about the enthalpies. Below is a plot of the enthalpies of the liquid and vapor
streams
Both streams show variations in enthalpies, especially below the feed. Most of the changes in
temperature and composition occur below the feed:
situations in which the tray efficiency is low (perhaps only 5% instead of 50%). Because of this
very low efficiency, very large numbers of trays would be required -- perhaps 100's or 1000's.
Fabrication costs just become prohibitively expensive. Fortunately, a viable alternative exists:
the packed tower.
A packed tower is simply a tube or pipe, which is filled with some sort of "packing." The
packing typically consists of particles around an inch in diameter. In commercial packed towers,
the usual choice are particles with one of three different shapes:
• Raschig ring (which is just a piece of pipe which has been cut into segments, whose length
and diameter are about the same)
L ≈ D ≈ 1 to 1 1 inches
2 2
• Berl saddle
• Pall ring
Although, in a pinch, almost anything you have laying around would do -- ping-pong balls, golf
balls, etc. The purpose of the packing is to promote good contact between the liquid and vapor
streams which are being brought together to permit interfacial mass transfer. In particular, what
is desired is a large interfacial area per unit volume.
The liquid stream is usually fed into the top of the tower while the vapor is fed into the
bottom. Thus we have countercurrent flow of the two streams, which has the same advantages
for mass transfer as it did for heat transfer.
The packing promotes good contact between the phases by dividing the two feed streams into
many parallel interconnected paths. Ideally, you would like the liquid to flow downward as a
thin film over the surface of the packing. This would give the maximum surface area of contact
between the gas and liquid.
If you just pour the liquid from the end of the pipe onto the top of the packing in tower
having much larger diameter than the pipe, most of the packing will not even be wet. Only some
of the channels will be carrying flow. This is called:
So some sort of device to distribute the flow over the entire cross section of the tower is needed.
This device is called a distributor.
Even if the flow is evenly distributed at the top of tower, channeling might still develop as
the fluid trickles down. When two thin films converge they tend to form a thick film and a dry
patch, which results in a reduction in contact area. So redistributors are placed every 10-15 feet
along the length of the tower.
• flowrate of water
• diameter of tower
• height of packing
Solution Overview:
Just as with distillation, there is a minimum L/V required to achieve the desired separation.
In typical operation of an absorber, the liquid rate is chosen to be just above the minimum.
The diameter of the tower is usually chosen on the basis of gas mass velocity GG. Generally,
for a particular L and V, the smaller the column diameter, the larger the mass velocities will be,
and the larger the pressure drop. Generally, large mass velocities are desirable because they give
high mass transfer coefficients, but too large mass velocities cause “flooding” which severely
decreases mass transfer rates. The gas mass velocity above is about as close to flooding as you
dare get.
b y
V S dy
Tower height:* ZT =
Kya ∫ y − y*
ya
H Oy
NOy
The height of the tower is determined by mass transfer rates. Basically, the gas and liquid
phases need to be in contact for a certain time for the acetone to have time to diffuse from the gas
phase into the liquid. The equation above is called the "design equation". The integral
yb
dy yb − ya
NOy ≡ ∫ y − y*
=
∆y
ya
is called the number of transfer units, where y is the mole fraction of the transferable
component, y – y* is the local driving force for mass transfer and ∆y is some average driving
force (perhaps the log-mean of the driving force at the top and bottom). The number of transfer
units is loosely analogous to the number of ideal trays required: both can be thought of as
measures of the difficulty of the separation.
Since NOy is dimensionless, to get units of height, the remaining factor in the design equation
(defined to be HOy) must have units of height and is called the height of one transfer unit:
V S
H Oy ≡
Kya
* This particular formula is for dilute solutions, which is probably not really applicable to the current example, but it
is applicable in the homework problems.
where Kya is the overall mass transfer coefficient times the interfacial area per unit volume of
packing.
Detailed Solution:
First of all, what is “SCFM”? This stands for “standard cubic feet per minute”
where “standard” means that the volume is evaluated at standard conditions. We can convert
SCFM into moles. Starting with the ideal gas law,
PV = nRT
V RT R ( 273° K ) liters ft 3
or = = = 22.4 = 359
n P (1 atm ) gmol lbmol
depending of which set of units are used to express R. Recall that a gmol is the quantity of
material whose mass equals the molecular weight in grams. Similarly, a lbmol is the quantity of
material whose mass equals the molecular weight in pounds. For example, air has a molecular
weight of 29:
g lb
MW= = 29
of air 29
gmol lbmol
Now let’s move on to determine the liquid flowrate required. Recall that in distillation the
slope of the ROL is given by
L R
=
slope of ROL =
V R +1
You should also recall that there exists a minimum value of the reflux ratio for any given
combination of product and feed specifications. This minimum value of R translates into a
minimum allowable L/V. There is a similar minimum L/V for gas absorption. To find the
minimum reflux ratio in distillation, we needed to plot an operating line and an equilibrium
curve.
EQUILIBRIUM CURVE
In Prob. 22-1, McCabe gives us equilibrium data in the form of partial pressure of acetone in
the gas phase, pA, as a function of the mole fraction of acetone in the liquid, x:
o
p=
A PA γ A x (155)
Activity coefficient is a measure of nonideality of the liquid phase. For ideal solutions the
activity coefficient is unity:
γA = 1 (ideal solution)
For ideal solutions and γA = 1, (155) reduces to Raoult’s law. Our solution of acetone in water is
not ideal, but McCabe kindly gives us a model
A 1.95 (1 − x )
2
ln γ=
where we use P = 1 atm. Repeating this for different x’s to obtain y’s up to 0.14 (the feed
concentration), we obtain the curve at right, where
OPERATING LINE
Suppose we are trying to absorb acetone from a mixture of air and acetone by contacting the
air mixture with water. Let
L V y − La xa
=
y x+ a a
V V
This relationship y(x) between the mole fractions in the liquid and gas streams is again known as
the operating line, since the relationship is imposed by a component mole balance.
Unlike distillation (where the equimolal overflow assumption makes the ratio L/V constant
within any cascade), L and V are not constant over the height of the absorber: acetone is
transferred from the gas to liquid, making L and V change. In particular, their ratio changes.
Thus the operating line is not straight. This makes it difficult to determine the minimum water
flowrate.
There is one limiting case where we can easily predict the change in L and V:
It is often possible to assume that only one component is undergoing transfer between the
liquid and gas streams. For example, in our problem acetone is being transferred from the air to
the water:
transferable: acetone
Although some of the water will evaporate when it contacts the air and some of the air will
dissolve in the water, the molar rates of transfer of these components can often (but not always)
be neglected compared to rate of acetone transfer. When you can neglect the molar transfer rate
of all but one component, then significant simplification can be made. If we can neglect
evaporation of water and dissolution of air, then the moles of water (the nontransferable
component) in the liquid stream must be the same at all elevations:
L′
(1-x)L ≡ L' = const so L= (157)
1− x
In the above equation and what follows, a prime (' ) is used to denote a flowrate or mole
fraction of the nontransferable component. Similarly, the moles of air in the gas stream must be
the same at all elevations:
V′
(1-y)V ≡ V' = const so V= (158)
1− y
L ′xa V ′y L ′x V ′ya
(157) and (158) into (156): + = + (159)
1 − xa 1 − y 1 − x 1 − ya
(1 − ya )( x − xa ) R + ya (1 − xa ) (1 − x )
y ( x) =
(1 − ya )( x − xa ) R + (1 − xa ) (1 − x )
L'
where R≡
V'
x moles of A in liquid
X ≡ =
1 − x moles of non-A in liquid
♦ To confirm that the lines are not straight, align a straight-edge (e.g. a ruler or the just the edge of a paper
envelope) with any of the red lines.
y moles of A in gas
Y≡ =
1 − y moles of non-A in gas
L′ L′
or Y= X + Ya − X a (160)
V′ V′
0≤x≤1
but 0≤X≤∞
xa
xa = 0.0002 →=
Xa = 0.0002
1 − xa
V′ V′
ya = 0.05 yb
1 − ya 1 − yb
= =
Ya 0.05 Yb 0.00814
0.14
= 0.1628
1− 0.14
Ya
or =
ya = 0.00807
1 + Ya
Y = Yb = 0.1628
Extending this straight line from the pinch up to Y = Yb, we locate a second point on the line:
(Xb,max,Yb) = (0.081, 0.1628). Using the two points (Xa,Ya) and (Xb,max,Yb), we can calculate
the slope corresponding to the minimum water flowrate:
L′ 0.1628 − 0.00814
=
′ = 1.91
V min 0.081 − 0.0002
′ =
Lmin 1.91V ′ = {
1.91 × 1.39
lbmol
min
(1 − 0.14 ) =}
2.29
lbmol
min
=4.94
gal
min
Multiplying this by the molecular weight of water (18 lb/lbmol) and dividing by the density of
water (8.33 lb/gal) yields a minimum water flowrate of 4.94 gal/min.
yb
V S dy
ZT =
Kya ∫ y − y*
ya
where V is the molar flowrate of gas, S is the cross-sectional area of the tower and Kya is a mass
transfer coefficient. The denominator of the integrand, y-y*, is the local driving force for mass
transfer; on an xy diagram, it is also the vertical distance between the operating line y(x) and the
equilibrium curve y*(x).
In distillation, the tower height was determined by the number of plates required times the
plate separation (which is usually 1-2 feet). The number of plates required for a given separation
is determined by the operating and equilibrium lines. For packed towers, the height also depends
on the operating and equilibrium lines, but in addition it is inversely proportional to the mass
transfer coefficient, which in distillation plays only a minor role in determining plate efficiency.
Now that we have established the importance of interfacial mass transfer in packed towers,
let's talk about modelling mass transfer across a phase boundary. There are two main differences
between interfacial heat transfer (as employed in Hxer design) and interfacial mass transfer:
yi ≠ xi
T yi = T xi = T i
The reason for this discontinuity in concentration across the interface has to do with
thermodynamic criteria for phase equilibrium. Recall:
thermal equilibrium: TV = TL
To illustrate this, recall the simplest case of VLE: an ideal gas mixture in equilibrium with an
ideal solution. This leads to Raoult's law:
o
=
p j y=
j P x j Pj
♣ Instead of component i, the superscript “i ” will be used to denote quantities evaluated at the interface between
two phases. In place of i, we will use j to denote components.
y j Pjo
Note that: = ≠1
xj P
even in the simplest case of vapor-liquid equilibrium, the mole fractions are not equal, except in
the trivial case when you have only one component and the total pressure is the vapor pressure.
Recall from the first part of this course, the local heat flux through the interface can be
related to the local temperatures using any one of three types of local heat transfer coefficients:
Basically these equations say that the rate of heat transfer is proportional to the driving force
(which is the departure from equilibrium) and the proportionality constant is the heat transfer
coefficient.
Comparing the mass transport expressions in (162) with their heat-transfer analogues in (161),
there is a good deal of similarity — especially in the one-phase relations. For a single phase, the
driving force is the difference between the concentration or temperature in the bulk and in the
concentration or temperature at the interface (but in the same phase).
But the driving force for the overall coefficients look a little different. The overall driving
force for heat transfer is just the difference between the bulk temperatures of the hot and cold
fluid
Th - Tc
By simple analogy, you might expect the overall driving force for mass transfer to be the
difference in concentration of the two phases:
Instead y - y* and x* - x
Example #1
Find: xi and yi
If the fluxes were not equal, we would have acetone building up at the
interface. We will denote this interfacial flux by NA.
Now the bulk compositions — x, y — are known, together with the two single-phase mass
transfer coefficients. Think of the two interfacial concentrations as two unknowns. We need
two equations. One equation is provided by the requirement that (xi, yi) must lie on the
equilibrium curve. The second relation is (163), which can be re-written as a linear relationship
between yi and xi:
k k
yi =− x xi + y + x x (164)
ky ky
slope intercept
Example #2: Next, let’s determine the value of the overall coefficient which leads to the same
flux for these interfacial concentrations.
Given: kx and ky
Find: Ky
The relationship among the concentrations is shown in the figure above. Adding and subtracting
yi from the overall driving force y-y*:
y − y * = y − y + yi − y * (166)
i
NA NA
Ky ky
Using (165) we can assign meaning to two of the three differences appear above, leaving:
NA NA
= + yi − y * (167)
Ky ky
( xi − x ) m
The remaining difference yi-y* can be related to xi-x if we multiply and divide the last term by
xi-x:
y − y * NA
yi − y* = ( xi − x ) i = m (168)
xi − x kx
NA
m
kx
1 1 m
= + (169)
K y k y kx
y − y*
where m≡ i = avg. slope of EC
xi − x
is the average slope of the equilibrium curve in the region of the xy diagram between x and xi.
Of course, if the equilibrium curve is straight (as it will be in dilute solutions), then m is its true
slope over the entire range.
Comment: This is similar to the expression for overall heat transfer coefficients for a double-pipe
heat exchanger: ♣
1 1 1
Di = Do: = + (170)
U hx hy
We said that 1/U represented the total resistance to heat transfer through the two phases, which is
just the sum of the resistances of each phase. The main difference between (169) and (170) is
the appearance of m in (169). For heat transfer, the slope of the equilibrium line is unity (m = 1)
because at thermal equilibrium Ty = m Tx = Tx.
In short, (169) states that the total resistance to mass transfer equals the sum of the
resistances of the gas and liquid phases. We could also have showed, in a similar fashion, that
1 1 1
= + (171)
K x m′k y k x
♣ In our earlier treatment of heat transfer, we also distinguished between inner and outer coefficients (i.e. U and
i
Uo). This was because the inside area and outside area of a pipe were somewhat different. We will drop this
distinction here because the inside area and outside area are equal for an interface. In other words, the ‘pipe wall’
thickness is zero.
y − yi
where now: m′ ≡
x * − xi
1 m
straight EC: =
K y Kx
Diffusion-Induced Convection
We have just seen that the overall driving force for interfacial mass transfer is y-y* or x*-x (it
is not y-x). This is one of two major differences between interfacial mass transfer and interfacial
heat transfer. The second major difference arises because diffusion of mass usually induces
convection of the fluid (i.e. motion of the center of mass).
dT energy
qz =
−k = ]
heat flux [ =
dz area × time
If there is also bulk fluid motion (due to pressure-driven flow, for example), an additional term is
added to the heat flux to account for convection of heat:
energy vol
qz = −k
dT
(
+ ρc p T ( z ) − Tref vz
)
dz
conduction convection
where vz is the local z-component of velocity of the fluid and ρcp(T–Tref) is the enthalpy per unit
volume. Usually heat conduction does not induce flow.*
By contrast, diffusion of mass nearly always induces flow. Fick’s law of diffusion must be
modified by adding a contribution from convection of mass:
*An important exception is “natural convection” which arises because the nonuniform temperature of the fluid
causes nonuniform density of fluid which in a gravitational field can cause flow.
mol vol
dc A N +N
N Az =
− DAB + cA Az Bz (172)
dz
c
diffusion v = mol average
velocity
convection
where c = cA+cB is the total molar concentration. The second term represents convection of
matter which arises because the center of mass moves with time.
n P
=
c = = const
V RT
Actually there is second and an even stronger argument for equimolar counter-diffusion
when the two reservoirs are rigid and finite in volume. Any departure from equimolar counter-
diffusion of A and B would cause the pressure in one reservoir to increase with time and the
pressure in the other reservoir to decrease. Any difference in pressure between the two
reservoirs would cause pressure-driven flow of the gas mixture toward the lower pressure (i.e.
convection). Unless the diffusion length is very short, convection generally dominates any
diffusion. This is why diffusion boundary layers tend to be very thin. Convection then would
quickly erode any pressure-difference.
In this way, the steady mole-average velocity v = 0 and the total pressure in both reservoirs
remains constant. For NAz = –NBz, (172) gives
dc A
N Az = − DAB (173)
dz
See Example 6.2-1 in Geankoplis for details leading to a linear profile of partial pressure or of
molar concentration of either gas. For 1-D steady equimolar counter-diffusion (from a reservoir
having cA1 to a second reservoir having cA2), taking the flux NAz in (173) to be independent of z,
(173) integrates to yield cA varying linearly with z. Substituting this result into (173) yields
DAB cDAB
=
N Az ( c A1 −=
c A2 ) ( y A1 − y A2 )
z2 − z1 z2 − z1 (174)
≡ k y′ ( y A1 − y A2 )
where z2-z1 is the length of the capillary through which diffusion is occurring. The last equation
is the definition of k'y, the single-phase (gas) mass transfer coefficient to be used with mole
fractions and equimolar-counter diffusion.
If instead, we have diffusion of benzene vapor above its liquid through air (which is virtually
insoluble in the liquid), the air must remain stagnant (i.e. NBz = 0) since it cannot enter the liquid
(see Fig. 6.2-2a). This is diffusion of A through stagnant B. For NBz = 0, (172) gives
dc cA dc dy
N Az =
− DAB A + N Az or (1 − y A ) N Az =
− DAB A =−cDAB A (175)
dz c dz dz
yA
cDAB
=N Az (y − y )
( z2 − z1 )(1 − y A )M A1 A2 (176) pB
≡ k y ( y A1 − y A2 )
z1 z2
(1 − y A1 ) − (1 − y A2 )
where (1 − y A )M ≡
(1 − y A1 )
ln
(1 − y A2 )
represents the log-mean of 1-yA evaluated at either end of the diffusion path. Once again, the
second equation in (176) is the definition of ky, the single-phase (gas) mass transfer coefficient to
be used with mole fractions and diffusion through stagnant film of B.
∂c A ∂N Az
+ =
0
∂t ∂z
then requires the flux NAz to be constant with respect to z. Taking NAz and cDAB as constants,
(175) is a first-order, ordinary, differential equation in yA(z), with constant coefficients, whose
general solution is
N
y A ( z ) = 1 + A exp Az z (177)
cDAB
where A is an integration constant. Applying the boundary condition yA(z1) = yA1, we could
evaluate the integration constant A:
N
A =− (1 − y A1 ) exp − Az z1
cDAB
N
(1 − y A1 ) exp Az ( z − z1 )
1 − yA ( z ) = (178)
cDAB
cDAB 1 − y A2
N Az = ln (179)
L 1 − y A1
1 − y A ( z ) 1 − y A2
z L
=
1 − y A1 1 − y A1
Finally, if we multiply and divide (179) by the driving force yA1–yA2 across the capillary tube,
then we obtain (176).
k y′
=ky > k y′ (180)
(1 − y A )M
For the same driving force (yA1 – yA2), the fluxes NAz are different. Since 1-yA is always less
than one, we see that equimolar counter-diffusion (174) is slower than diffusion through a
stagnant fluid (176). This can be qualitatively understood as follows. Suppose that to get to
class, you need to walk down a corridor that's crowded with other students. If everyone else was
standing almost still (stagnant fluid), it would be easier to walk around them than if everyone is
walking toward you (counter diffusion).
Although (180) was derived for one particularly simple case (no pressure-driven or gravity-
driven convection), it applies much more generally. Similar relations exist for diffusion in the
liquid phase:
k x′ (1 − x A1 ) − (1 − x A2 )
kx = where (1 − x A )M ≡
(1 − x A )M (1 − x A1 )
ln
(1 − x A2 )
and between the overall mass-transfer coefficients for equimolar counter-diffusion and diffusion
through a stagnant film:
K y′ K x′
Ky = and Kx =
(1 − y A )*M (1 − x A )*M
(
(1 − y A )*M =
1 − y*A ) − (1 − y AG ) (1 − x*A ) − (1 − x AL )
(1 − x A )*M = *
(1 − y*A ) (1 − xA )
where and
ln ln
(1 − y AG ) (1 − x AL )
and where yAG is the mole fraction in the bulk of the gas, y*A is the mole fraction which would
be in equilibrium with the bulk liquid having a mole fraction of xAL, and x*A is mole fraction in
the bulk of the liquid which would be in equilibrium with the bulk gas having a mole fraction of
yAG.
We could have used any one of the other expressions in (162). Just like the overall heat transfer
coefficient Uo depends on the flowrates of both fluids being contacted in the heat exchanger, the
overall mass transfer coefficient Ky depends on both flowrates
K y = K y ( L, V )
To obtain the rate of transfer across the interface in our slice of column, we need to multiply
(181) by the interfacial area in that slice. In a heat exchanger, the heat transfer surface is fixed
by the geometry of the equipment selected: it is just the area of the pipe wall or the tubes. In
particular, the heat transfer area does not depend on the flowrates of the hot and cold streams.
On the other hand, the boundary between liquid and gas in a packed bed is very complex. Most
importantly, the area also depends on the flowrates of the gas and liquid streams (L,V). The
interfacial area is usually expressed as a, the interfacial area per unit volume of packing:
interfacial area
≡a=a ( L, V ) (182)
volume of packing
One empirical correlation relating area to flowrates is the Schulman Correlation (see Table 6.3
of Treybal above):
a = mGxn G yp
where Gx and Gy are the mass flowrates of the liquid or gas stream divided by the cross-
sectional area of the tower
mass time
G≡ =
mass velocity
πD 2
4
and where m, n and p are constants which depend on the type and size of packing used.
Since Ky and a usually cannot be measured independently, what is usually found in correlations
of mass transfer rates is their product Kya.
rate of transfer = r S ∆z
πDT2
S = empty tower cross-section = (184)
4
S ∆z = volume of slice
At steady state, the rate of transport of acetone into the vapor must equal the rate out:
in = out
(Vy ) z + ∆z − (Vy ) z
or = −rS
∆z
d (Vy )
= −rS (185)
dz
Now the total molar flowrate (V) will change as the gas phase looses acetone, but the
flowrate of acetone-free air (V' ) does not change with z:
V′ V′y
V= → Vy =
1− y 1− y
d d V′y d y dy dz V dy
(Vy=) = V′ = = V′
=
dz dz 1 − y dz 1 − y (1 − y )2 1 − y dz
V dy
(185) becomes: =−rS =− K y a ( y − y *) S (186)
1 − y dz
The second equation is obtained by substituting our expression for the rate of absorption r from
(183). We can now solve this equation for the height dz of the slice needed for a certain change
dy in concentration:
−Vdy
dz =
(1 − y )( y − y *) ( K y a ) S
So if we know the local mole fractions y and y* and the change in the mole fraction dy of the
gas phase which occurred between the top and bottom of this slice, we can calculate the height of
the slice. The total height of packing is the sum of the height of each slice:
∫ ∫
ZT ya
−Vdy
=
ZT =dz
0 yb (1 − y )( y − y *) ( K y a ) S
Factoring out the S, which is constant along the entire height
of the packing, and changing the order of integration (since
yb > ya):
∫
yb
1 Vdy
ZT =
S ya (1 − y )( y − y *) ( K y a )
The equation above represents the most general form of
the design equation. Unfortunately, it hard to give any
meaning to the integral in this general form. So let’s examine some limiting cases. For the case
of a dilute gas stream, we know that
V S
∫
yb
dy
ZT = (187)
Kya
y y − y*
avg
a
H Oy NOy
which is a little simpler. We have added the subscript “avg” to suggest that, in case V and Kya
are not quite constant, we might still evaluate this quotient at each end of the packed height (i.e.
at y = ya and at y = yb) and average the two values.
EVALUATION OF NTU'S
To better understand how to
evaluate this integral, we will continue
the example we started a couple of
lectures back (see page 174). The
problem is summarized in the figure at
right. The water flowrate given in the
figure is 1.08 times the minimum,
where the minimum water flowrate L'
was determined on page 181:
L′ L′
= 1.08
= ′ 2.06
V′
V min
1.91
x x y y
L′ − a = V ′ − a (188)
1 − x 1 − xa 1 − y 1 − ya
On the other hand, y*(x) is the mole fraction of acetone in the gas, which would be in
equilibrium with a liquid having mole fraction x. Recall (155) and the text which follows:
′
pace Pace
y * ( x=
) x 0.33 x exp 1.95 (1 − x )
2
= γ= (189)
P P
yb
∫y
dy
= 21.0
a y − y*
dy
=
dy = dx y ′ ( x ) dx
yb dx xb y′ ( x )
∫ ∫
dy
= dx
ya y − y * xa y ( x ) − y * ( x )
The maximum value of the integrand occurred at the point in the column where the operating
line comes closest to the equilibrium curve. The denominator y-y* is a measure of this distance.
Because the equilibrium curve turns downward in this case, the "pinch" point occurs in the
middle of the range. In other cases (e.g. NH3 in water), the equilibrium curve turns upward
instead.
The figure at left below shows the OL and EC obtained in the solution of MS&H5 Prob. 22.2
involving the absorption of NH3 in water. The problem calls for removing 99% of the ammonia
from the air by contacting the ammonia-air mixture with pure water. This corresponds to (xa, ya)
= (0.0000, 0.0025) and (xb, yb) = (0.090, 0.200). Notice also that the equilibrium curve (red) is
concave up (unlike the previous example). The diluteness of both streams at the top of the tower
plus the concave-up shape of the EC cause a very tight “pinch” to occur at the left end of the OL
and EC
From the corresponding plot of the integrand vs y (plot on right side above), the denominator of
the integrand becomes very small near the lower limit; thus the integral has much of its area
under a “spike” located at the lower limit.
To integrate such a function accurately, you will need a very small grid spacing — which
means a lot of points — or at least very small spacing under the spike — which means a
nonuniform grid. One solution is to transform the integral by multiplying and dividing the
integrand by y:
yb yb yb
∫ ∫ ∫
dy y dy y
= = × d ln y
ya y − y * ya y − y * y ya y − y *
TRANSFER UNIT
Even if the integrand varies significantly over
the domain of integration, we can still pull it
outside the integral, provided we replace it with
some appropriate constant, which can be thought of as the mean value for this range. Recall the
Mean Value Theorem of Calculus:
yb dy yb yb − ya
∫ ∫
1
= = dy (190)
ya y − y * ( y − y *)avg ya ( y − y *)avg
Of course, it is not usually known what type of average to use, but for the present purpose, the
precise value is not important. The integral can be thought of as the total change in gas mole
fraction divided by the average driving force.
(Gas Phase) Transfer Unit -- a slice (not necessarily thin) of an absorber in which the gas
undergoes a change in y equal to ( y − y *)avg which represents the average driving force
over the entire absorber.
No. of TU's -- the number of these it takes to make the entire absorber (not necessarily an
integer)
From (190), we see that the integral clearly represents the number of transfer units:
yb
∫
dy
NOy = (191)
ya y − y *
More specifically, this integral is called the number of overall gas-phase transfer units since it is
based on the overall driving force expressed as the gas-phase mole fraction. Since the product in
(187) represents the total height of packing, the coefficient of this integral must represent the
height of a transfer unit:
V S
H Oy =
K y a avg
ZT = HOyNOy
which is called the HTU method for sizing an absorber. Throughout the derivation above we
made use of one type of mass transfer coefficient:
Of course, other types of mass transfer coefficient are commonly encountered and can also be
used to compute the height of packing required:
to mention a few. Of course, there are still more definitions based on driving forces expressed as
differences in molar concentration or partial pressure (see Table 7.2-1 in Geankoplis). Any one
of these mass transfer coefficients can be used in the HTU method:
where
xa dx L S
=NOx
∫ =
xb x * − x
H Ox
K x a avg
xa dx L S
=Nx
∫
=
xb xi − x
Hx
k x a avg
yb dy V S
=Ny
∫
=
ya y − yi
Hy
kya
avg
Note that the Ns are not equal to each other; thus the number of transfer units depends on which
driving force you use.
Ny ≠ NOy
It’s still not apparent what “average” to use to compute ( y − y *)avg . There is one case in
which we can provide a simple answer: For very dilute solutions,* this integral in (190) can be
evaluated analytically (i.e. numerical integration is not required). If both x<<1 and y<<1 hold,
* Up to now, “dilute” has meant y and x are small enough (compared to 1) so that V and L are virtually constant
throughout the tower. This makes the OL straight, but not necessarily the EC. When x and y are very small
compared to unity, the EC will also be straight: it should have the form of Henry’s law: y* = mx.
then the molar flowrates of the two streams won't change much as the solute is transferred across
the phase boundary. In particular, the ratio of the flowrates will be constant:
This means that the operating line will be straight, even on mole fraction coordinates:
Also, for sufficiently dilute solutions, Henry's law always applies so that
y* = mx
will also be a straight line. When both operating and equilibrium curves are straight, then the
integration can be done analytically instead of numerically. The result is:
yb − ya
straight OL & EC: NOy = (192)
( y − y *) L
where ( y − y *) L is the log-mean of yb-yb* and ya-ya* which represent the driving force at the
bottom and top of the tower, respectively.
( y − y *)avg =
( y − y *) L
(
=
yb − yb* ) − ( ya − ya* )
yb − yb*
ln
ya − ya*
This result is quite analogous to our design equation for double-pipe heat exchangers. To see
this analogy, first recall the design equation for double-pipe heat exchangers: see (9) on page 11:
qT = UAT∆TL
where qT is the total heat duty — in other words, the total rate of heat transfer. The analogous
quantity for gas absorption is the total rate of ammonia transfer over the entire packed tower,
which can be calculated from a mole balance as the difference between the rate of ammonia into
the tower (via the gas phase) minus the rate of ammonia leaving the tower:
qT = V(yb-ya) (193)
Recall the design equation (187). After substituting (192) into (187) we have:
V S yb − ya
= =
ZT H Oy NOy
K y a ( y − y *)
L
where SZT is the total volume of packing in the absorber and a is the interfacial area per unit
volume of packing. Thus
Thus, for dilute solutions, the design equation for an absorber is identical in form to that for heat
exchangers: with the overall mass transfer coefficient Ky analogous to the overall heat transfer
coefficient U and the log-mean of y-y* at the two ends of the absorber analogous to the log-mean
∆T.
Dividing both sides of (169) by a (the interfacial area per volume of packing), we obtain:
1 1 m
= + (194)
K y a k y a kx a
Of course, if the equilibrium curve is straight (as it will be in dilute solutions), then m is its slope.
If we now multiply (194) by V/S we have the height of an overall gas-phase transfer-unit on the
left-hand side:
V S V S mV S V S V LS
= + = +m
K ya kya kx a kya L kx a
H Oy Hy Hx
Recalling the definitions for various heights of a transfer unit we end up with
V
= Hy + m
H Oy Hx
L
Had instead we started with (171) and multiplied by L/S we would obtained
L
=
H Ox H y + Hx
m ′V
In our original motivating example (see page 174), we said the designer of an absorber needs
to size three quantities: the liquid flowrate, the tower height and the tower diameter. We have
already discussed the basis for choosing the first two. In this section, we address the third
quantity: tower diameter. The selection of this quantity involves controlling the pressure drop
across the tower. A typical plot of pressure-drop vs. the mass velocities of the gas and liquid
phases is shown in Figure 7:
Fig. 7. Pressure drop across a packed bed of 1-inch ceramic Intalox saddles as a function of air and
water flowrates. Taken from Fig. 22-4 of MS&H5.
lb/hr of liq M L
Gx = mass velocity of liq = = x
cross-section of tower S
lb/hr of gas M yV
Gy = mass velocity of gas = =
cross-section of tower S
where Mx and My are the average molecular weights of the liquid and gas, respectively, and S is
the cross-sectional area of the tower given by (184). The tower diameter DT exerts its influence
on packed-tower design through S in these equations for the mass velocities.
Notice in Figure 7 that the curves with nonzero liquid flow (i.e. “wet” packing) curl up at
high gas flow rates. This can be explained as follows:
The upward flow of the gas exerts a upward shear force on the liquid, retarding its downward
motion. To get the same flowrate with a slower velocity, you need a thicker film. Thus liquid
holdup increases with increased gas flowrate. As the gas flowrate increases so does this shear
force. When the shear force becomes comparable to gravity, the liquid flow might slow to zero
(rightmost figure above) and flooding occurs in the tower.
intersticial volume -- "empty" space between and inside packing particles; that volume
occupied by liquid or gas.
For a given liquid flowrate, there is a maximum flowrate of gas which can be forced through the
column. If you exceed this maximum, then the shear forces on the liquid are so high that they
exceed the weight of the liquid. Then the net force on the liquid is upward and you blow the
liquid back out the top of the column. This is called:
Since you can't get the liquid through the column, flooding is clearly a condition to be
avoided. Generally, to get a high area of contact between liquid and gas, you want to operate
near but comfortably below the flooding velocity:
Normal operating pressure drops are just below those for which loading begins.
where Fp is the packing factor for the particular packing used in the column (see Geankoplis
Table 10.6-1 reproduced as Fig. 10 below). This formula gives ∆Pflood in inches of water per
foot of packed height. Once Gyat flooding is determined then
Notice (from Fig. 7) that when we increase the liquid flowrate, flooding occurs at lower gas
flows. If we increase the size of the packing particle, we can tolerate higher gas flows. While
higher gas flows can be obtained with larger packing particles, the amount of interfacial area is
generally less (i.e. a is decreased). These effects are summarized by Fig. 8 (below).
Fig. 8. Air flowrate at which flooding occurs for ceramic Intalox saddles at different water flowrates.
Taken from Fig. 22.5 in MS&H5.
Fig. 9. General correlation of pressure drop across packed beds for random packings. Taken from
Fig. 10.6-5 of Geankoplis. Units of ordinate assume units given in table below.
Fig. 7 and Fig. 8 are for a particular packing (1” Intalox saddles) and for a particular
temperature and pressure (20 °C and 1 atm). A more general correlation of pressure drop is
provided by Fig. 9 where
mG
GG = = mass velocity of gas [=] kg-s-1-m-2
S
G
vG = G = superficial velocity of gas [=] m/s
ρG
µ
ν = L = kinematic viscosity of liquid [=] m2/s*
ρL
♦ 1 poise = g-cm-1-s-1 and 1 centipoise = 1 cp = 10-2 g-cm-1-s-1 (approximate viscosity of water at room
temperature).
* 1 stokes = 1 cm2/s and 1 centistokes = 10-2 cm2/s. Yes, the “s” belongs at the end here: it does not denote the
plural but is part of a man’s name (George Gabriel Stokes, who gave us “Stokes law”).
Fig. 10. Characteristics of common packings. Taken from Table 10.6-1 of Geankoplis.
TOWER DIAMETER
Once the total molar liquid and gas flowrates (L and V) are known, we can choose the
diameter of tower we need. The diameter of the tower is usually chosen such that the pressure
drop is some prescribed value below flooding:
where the flooding velocity GG,flood is found from Fig. 9 such that the pressure drop is given by
(195).
M GV
Similarly, GG = (196)
S
So to get a certain pressure drop for a certain set of flowrates, we must choose a particular
value for the tower diameter. Although it looks like you have to guess the diameter DT in order
to calculate GL and GG to get ∆p, a trial-and-error can be avoided by noting that the abscissa
(x-coordinate) of Fig. 9 does not depend on the diameter:
MLL
GL S M L
= = L
GG M GV M GV
S
Note that the unknown cross-sectional area S conveniently cancels out. What remains is known.
The procedure is illustrated by the following example:
Given: L, V
Find: DT
L′ lbmol V′ lbmol
=
Lb = 3.09 and =
Vb = 1.619
1 − xb min 1 − yb min
We will also need the average molecular weights of the two streams at the bottom (calculated by
taking a basis of one mole of feed):
gm
M G yb M ace + (1 − yb ) M
= = 33.1
air
mol
58.1 29
gm
and M L xb M ace + (1 − xb ) M=
= 20.8
H
2 O
mol
58.1 18
Finally we need the mass densities of the two streams. Water’s density we can look up (we
ignore the small amount of acetone dissolved in it).
gm
ρL =
1.00
cm3
The density of gases can be calculated using the ideal gas law:
P gm lb
=
ρG M G= 0.001343= 0.084
RT cm3 ft 3
GL ρG M L Lb ρG
= = 0.044
GG ρ L M GVb ρ L
We have 1-inch Raschig rings♠ as packing material (see original problem statement on p174).
According to Table 10.6-1, this packing has a packing factor of
Fp = 179 ft −1
inch water
∆Pflood =
2
ft of packing
Step 3) locate the point on the curve of Fig. 9 corresponding to ∆Pflood which also has
the calculated value of the flow parameter (abscissa)
♠ Named for its inventor, the German chemist Friedrich Raschig (1863 – 1928). In 1891, he started his own
chemical company which exists today as Raschig GmbH (www.raschig.de).
Step 4) read the corresponding capacity parameter (ordinate of this point) from Fig. 9
ρG
vG Fp ν 0.05 =
1.9 (197)
ρ L − ρG
where ν is the kinematic viscosity of the liquid (in centistokes). The kinematic viscosity is
gm
0.0085
µL cm-sec cm 2
=
ν = = 0.0085 = 0.0085 stokes
= 0.85 centistokes
ρL gm sec
1.00
cm3
Step 5) calculate GG which gives this value for the capacity parameter. This gives
GG,flood
ρ L − ρG −0.5 −0.05
=vG 1.9
ρG
( Fp ) ν
0.5
1.00 ( −0.5
( 0.85 )−0.05 3.9 ft
=
1.9 179 )
0.001343 sec
The final units were assigned using the table located below Fig. 9. Multiplying this velocity
times the density gives the corresponding mass velocity of the gas phase:
lb ft lb
GG ,flood =
ρG vG =
0.084 3 3.9 sec =
0.327
2
ft ft -sec
Step 6) calculate the operating value of GG using GG/GG,flood ≈ 0.5 – 0.7. We will take
the lower end of this range (to be on the safe side).
lb lb
= =
GG 0.5GG ,flood 0.5 0.327= 0.164
ft 2 -sec ft 2 -sec
πDT2 M GV
= S=
4 GG
First calculate the desired cross-sectional area of the tower from (196):
lb lbmol
33.1 1.619
M GVb lbmol min
=S = = 5.5 ft 2
GG lb sec
0.164 2 60 min
ft -sec
4S
=
DT = 2.6 ft
π
Another important unit operation for mass-transfer is the membrane separator. Here the two
phases are miscible — usually they are both gases or both are aqueous solutions. Because of
this, they cannot be intimately contacted. Instead, a membrane is placed between them. This
membrane serves two functions:
2. controls selectivity
Because of differences in molecular size or solubility of the gas in the membrane material,
different species diffuse through the membrane at different rates. This difference is the basis of
separation.
EXAMPLES
With the proper choice of membrane material, a membrane separator can achieve significant
separations in a single stage. Indeed, as we shall see, it usually turns out that multistage
operations with membranes are not economical (compared to alternatives like cryogenic
distillation). Most industrial applications of membranes do not require very high purity products.
Gases
• separation of 238UF6 from 235UF6. This was a very important part of the atom bomb
project in World War II.♠
• separation of He from CH4 (natural gas). Because of its very small molecular size, half
of the helium in a feed stream containing <1% He can be recovered in a single stage,
resulting in a permeate stream contain 30 times greater concentration of He.
Liquids
• dialysis (or liquid permeation). Separation of urea and other small molecules from protein
and cells in whole blood is accomplished in artificial kidneys. Recovery of caustic
(NaOH) from hemicellulose produced in the manufacture of rayon.
• electro-dialysis. Brackish water (containing less salt than sea water) can be made potable
by applying electric field across membrane.
• reverse osmosis. Brackish water can also be made potable by applying pressure to force
water through a membrane which is impermeable to ions.
In the 1950’s, President Eisenhower was convinced that the easy availability of drinking water
was the key to world peace in many parts of the world. Eisenhower’s successor as President was
John Kennedy. In a famous speech, President Kennedy expressed his intention to “land a man
on the moon and make the deserts bloom.” The second part of this quote has been obscured by
the first part. Nonetheless, Kennedy did follow through by creating the Office of Saline Water,
which provided financial support for much of the early research on membranes.
EQUIPMENT
In our schematics of membranes, the membrane will be drawn as if it is a simple planar
sheet. The problem with this geometry is that you are very limited about how much area you can
get in a given volume, sort of like the double-pipe heat exchanger. One way to get a greater area
♠ See Henry De Wolf Smyth, Atomic Energy for Military Purposes, “Chapter X: The Separation Of The Uranium
Isotopes By Gaseous Diffusion” (available at http://www.atomicarchive.com/Docs/SmythReport/smyth_x.shtml ).
For a less rigorous summary, see https://en.wikipedia.org/wiki/Gaseous_diffusion under “History”
in a given volume is to stack the sheets using spacers which allow for inflow and outflow to each
sheet in the stack. If you roll-up the stack of sheets, you have the spiral-wound membrane (see
below).
This is the first of two geometries used. The second is a variation on the shell-and-tube heat
exchanger, namely a bundle of tubes. Suppose instead of sheets, we can construct the membrane
as a hollow tube. Since we want to operate the membrane under very high pressure drops, the
membrane “tubes” are actually fibers with a diameter of 50-200 microns (1000 microns = 1
mm):
Above is a photomicrograph of one such hollow fiber. The actual membrane of this fiber is
only 2 microns thick. The thickness of the annulus representing the actual membrane is marked
by the two red arrows. The rest of the fiber (the spongy annulus) is a porous support for the
membrane which allows very large pressure drop (perhaps 1000 psi) to occur without exploding
the membrane.
tube sheets
Hundreds of these fibers are then arranged in a bundle with each of their ends sealed in a disk
which will separate the tube-side flow from the shell-side flow much like a tube-sheet does in a
shell-and-tube heat exchanger (see schematic above). In heat exchangers, the metal tubes would
be welded to the tube sheet, one at a time, usually by hand. These hollow fibers are made of a
polymer (plastic) as is the tube sheet. In order to make these hollow-fiber bundles economically,
a technology was developed to seal all the fibers in the tube sheet at once. This is perhaps the
most proprietary part of making these devices.
The figure above shows the radial-flow pattern inside bundle (end view). Shown below are
two possible axial-flow patterns inside bundle (side views).
♣ DuPont patented the first commercial reverse osmosis membranes for treating brackish water in 1969 and
improved Permasep to the point that it was capable of desalinizing seawater in 1974. See
http://www2.dupont.com/Phoenix_Heritage/en_US/1969_a_detail.html.
♥ http://www.airproducts.com/products/Gases/supply-options/prism-membranes/prism-membrane-products.aspx
(a) shows hollow fibers with right ends sealed and left ends open and exiting to the manifold. V
denotes the flowrate of all gas permeating through the membrane whereas L2 denotes the
remainder of the feed which does not permeate (i.e. the “retentate”). (b) shows hollow fibers
with both ends open and exiting to the manifolds.
SEPARATION OF GASES
Porous Membranes
Perhaps the simplest basis for separation of two gases is molecular size. Gases of different
molecular weights have different diffusion coefficients and will diffuse at different rates through
the membrane.
For equimolar counter-diffusion of gases A and B, the Chapman-Enskog theory predicts the
diffusion coefficient is given by*
1 1
T3 +
M A MB
DAB = 0.0018583
pσ2AB Ω D, AB
where Mi is the molecular weight of species i, p is pressure, σAB is the collision diameter (sum of
radii of molecules A and B), ΩD,AB is a function of kT/εAB, and εAB is the depth of the Lennard-
Jones well for the interaction of molecules of A and B. Note that larger molecules diffuse more
slowly.
For membranes having pores whose diameter is comparable or smaller than the mean-free
path of gas molecules, the diffusion coefficient is inversely proportional to the square-root of
molecular weight:
T
for Knudsen diffusion:♦ DA = 9700 r (198)
MA
Notice once again that larger molecules (i.e. larger MA) have smaller diffusion coefficients.
dC A C (0) − CA ( L )
J Az =
− DA DA A
=
dz L
For membranes having fine pores through which the gas can
diffuse, the fluxes of two gases, A and B, through a porous
membrane are given by
C A (0) − C A ()
J Az = DAe (199)
τ
CB (0) − CB ()
J Bz = DBe
τ
where l is the thickness of the membrane and τ (which is greater than unity) is a dimensionless
factor called the tortuosity which accounts for the fact that the pore length is generally larger
than the thickness of the membrane. Also the flux is the rate of diffusion per unit of total
membrane area; if only a fraction ε of this membrane area is open to diffusion, then the effective
diffusion coefficients are a fraction of the true molecular diffusion coefficient:
DAe = ε DA
where ε is called the porosity of the membrane. Even if there are equal driving forces for the two
components to diffuse across the membrane, they will have unequal concentrations one the
downstream side because of their different diffusion coefficient. This is the basis for separation
in a porous membrane.
The best known example of gas separation by a porous membrane is separation of 238UF6
from 235UF6. Most gases have much smaller molecules than these and polymer membranes are
used instead of porous membranes
C A = S A p A and CB = S B pB (200)
Again, at steady state the concentration profile inside the membrane will be linear and Fick’s law
yields:
For the purpose of comparison, the flux through a porous membrane [given by (199)] can
also be written in terms of the partial pressures. For an ideal gas, the molar concentration CA can
be related to the partial pressure of the gas PA by
p AV = n A RT
nA p A
C=A =
V RT
DA p A (0) − p A ()
(199) becomes: J Az = porous membrane
RT
p (0) − p A ()
(201): J Az = DA S A A polymer membrane
Although these two expressions have the same form, the coefficients are different. In particular,
the solubility constant SA varies much more (between one gaseous solute and another) than does
the diffusion coefficient DA. This implies that polymer membranes tend to be much more
selective in separation various gas species than porous membranes. Unfortunately, diffusion
coefficients in solids and liquids are much smaller than for gases so this increase in selectivity is
often at the expense of lower permeation rates.
which has units of flux per unit pressure gradient. Since these complex units occur frequently in
such problems, a new unit is defined, which is known as the Barrer:
Often the thickness of the commercial membranes is not accurately known. Then it is easier to
measure the flux per unit pressure drop, rather than per unit pressure gradient. The flux per unit
pressure drop is known as the
DA S A
permeance: QA ≡
J A= QA ( p A1 − p A2 )= QA ( xp1 − yp2 )
Because we are using a vacuum pump to pull the gas through the membrane, the pressure on the
permeate side of the membrane will be much below atmospheric, thus:
J A ≈ xQA p1 (202)
J B ≈ (1 − x ) QB p1
Since the entire permeate stream must come through the membrane, the mole fraction of A in
the permeate stream is determined by the relative amounts of the two species which come
through. In time t the number of moles of species A which comes through a membrane of area A
is JAAt. Similarly, the number of moles of species B which comes through is JBAt. The mole
fraction of A is the number of moles of A divided by the total number of moles. The area and
time cancel leaving
J A At JA
=y = (203)
J A At + J B At J A + J B
xQA p1 xQA
=y =
xQA p1 + (1 − x ) QB p1 xQA + (1 − x ) QB
αx αx
=y = (204)
αx + (1 − x ) ( α − 1) x + 1
QA
where α≡
QB
is called the separation factor. It's role is similar to relative volatility in distillation [see (68) on
page 82] so we give it the same symbol.
Example: a fairly typical value for the selectivity for the separation of oxygen and nitrogen in a
polymer membrane is α = 7.
Find: y
If we did not make the assumption of a vacuum on the downstream side, the analysis
becomes a little more complicated, but it is still tractable. More generally, (202) becomes
=J A QA p1 ( x − Ry )
(205)
J B QB p1 (1 − x ) − R (1 − y )
=
p
where R ≡ 2 <1
p1
and (203) yields a value for y which is the root of the following quadratic equation:
R ( α − 1) y 2 − {( α − 1) ( x + R ) + 1} y + αx = 0 (206)
1
R=0
0.1
0.9 0.2
0.3
0.4
0.8
0.5
1
Permeate Mole Fraction, y
0.7
0.6
0.5
0.4
0.3
0.2
1 .10 1 .10
3 4
1 10 100
Note that if R = 0, (206) reduces to (204). On the other hand, when the pressure ratio
approaches unity (i.e. R→1), nonextraneous root of the quadratic corresponds to the permeate
having the same concentration as the feed:
lim y = x
R →1
αQB p1 ( x − y )
JA =
=J B QB p1 ( y − x )
JA α QB p1 ( x − y ) α( x − y)
=y = = (207)
J A + J B α QB p1 ( x − y ) + QB p1 ( y − x ) ( α − 1) ( x − y )
While it is tempting to cancel out the common factor of x-y in the numerator and
denominator, this leads to y = α/(α-1) which is a mole fraction larger than unity (recall α>1).
Actually, by cancelling out the common factor, we obtain the extraneous root of the quadratic [in
y(x,α)] which (207) represents before cancelling the common factor. Instead of cancelling the
common factor x-y, suppose we multiply both sides by the denominator: (α–1)(x–y). This results
in a quadratic for y(x,α). The nonextraneous root of this quadratic is y=x and then (207) is
indeterminant ( 0/0 ).
Above is a plot of y versus α for several different R’s. A couple of trends are apparent. First,
lim y = x
α→1
α=1 means there is no difference in permeability of the membrane to oxygen or nitrogen: then
the permeate has the same concentration as the feed. Second, note that the permeate tends to
pure oxygen if the permeability ratio is high enough:
1 for R ≤ x
lim y = x
α→∞ R for R > x
A closer look reveals that pure oxygen is only obtained if the pressure ratio is small enough. In
particular, the ratio must be smaller than the feed mole fraction. When R<x, then y→1 as α→∞.
This behavior can be rationalized as follows: if oxygen is to permeate the membrane, the
partial pressure of oxygen in the permeate can never exceed that of the feed:
yp2 ≤ xp (208)
1
pO2 in permeate pO2 in feed
x
(209) into (208): y≤
R
If this inequality is not met, we would have diffusion uphill. The upper bound in oxygen partial
pressure is reached when α → ∞:
x
lim y = min 1,
α→∞ R
Of course, the mole fraction must always be less than unity as well as meeting this constraint.
Stage Cut
Although we did not state it as an explicit assumption, the analysis in the previous section
only applies if the fraction of the feed permeating the membrane is negligibly small. If the
fraction permeating the membrane is significant, a number of complications arise. For example,
the concentration of the gas remaining (in the “residue” or “concentrate” or “reject” or
“retentate” stream)♦ will be different from the feed. This requires the gas on the high pressure
side of the membrane to vary continuously in concentration from the feed end of the membrane
to the residue end. It is no longer clear what average should be used as x in the formulas above.
The fraction of the feed which permeates through the membrane is known as the stage cut:
This also means the composition of the feed doesn’t change much as a result of passing though
the separator, so that
If the stage cut is not small, we will partially deplete the feed in the more permeable
component: oxygen and x then varies locally along the surface of the membrane.
♦ These are all names applied to the same stream: the outlet stream
for gases which do not permeate the membrane.
♣ This inequality assumes x is the mole fraction of the more
permeable component.
residue (a.k.a. “reject”). The simplest case (see figure at right) is to assume the gas composition
is uniform on both sides of the membrane despite the reject composition being different from the
feed. The stirrer in the schematic suggests that the gas is well mixed like a well-stirred tank.
The concentrations on either side of the membrane were also assumed to be uniform in the
derivation of (206), which still applies. The only remaining question is what is x in this situation
with significant stage cut? xin or xout? In a well-stirred tank, the steady-state concentration in
the tank is the same as the concentration xout of the output stream, which might differ from that
of the feed xin. This suggests that x should be replaced by xout:
R ( α − 1) yout
2
− {( α − 1) ( xout + R ) + 1} yout + αxout = 0 (210)
Now let’s look at the effect of stage cut on the permeate concentration.
Solution: xout, xin and yout are also related by a mass balance:
=
Lin xin Lout xout + Vout yout
Dividing through by Lin and relating the flowrate ratios to the stage cut θ:
Lout V
xin= xout + out yout= (1 − θ ) xout + θyout
Lin Lin
1−θ θ
x − θyout
Solving for xout or yout: xout = in (211)
1− θ
x − (1 − θ ) xout
yout = in (212)
θ
Mole Fraction
also approaches xin as the pressure
0.5
ratio R approaches 1. 0.4
Of course, this makes sense because if the entire feed stream permeates the membrane; then the
permeate must have the same composition as the feed (by a mass balance). To obtain the
minimum value of xout, we substitute yout = xin into (210) and solve for xout:
Example: If air is used as feed in a membrane having α=7, what is the richest concentration of
nitrogen that can be obtained?
Solution: we want the smallest value of xout predicted by (213). Since α>1, the smallest xout
corresponds to R=0 (see Figure 15). Substituting R=0, α=7 and xin=0.21, (213) predicts xout =
0.0365 or 1-xout = 0.964 or 96.4% nitrogen.
liquid-liquid extraction — separation of two components of a liquid (the “feed”) by contact with
a second immiscible liquid (the “solvent”)
Extraction is used primarily when direct distillation is not economical (since the solvent usually
has to be recovered by distillation). This might be the case with mixtures having very similar
boiling points, or substances that cannot withstand the temperatures of distillation, even under
vacuum. Examples include:
• recovery of acetic acid from dilute aqueous solutions by contact with ethyl acetate
[CH3CH2COOCH3] (fingernail polish solvent) or ethyl ether. Extraction greatly
reduces the amount of water which has to be distilled.
• separation of high-MW fatty acids from vegetable oil by contact with liquid propane
EQUIPMENT
Contacting of the two liquids might be
performed in a countercurrent cascade of stages.
Sometimes this can be accomplished in a sieve-
tray tower like that used in distillation (see figure
at right), where the difference in density drives
the heavy fluid down while the lighter one floats
upward.
To understand why these names are used, consider the most common application: to extract a
solute from an aqueous solution by contacting the solution with some oil phase. Since the role of
the oil phase is to extract the solute, the oil-rich product stream is called the “extract.” What
remains of the aqueous phase might be considered to be “refined” so the water-rich product
stream is called the “raffinate.” The word raffinate comes from the French word raffiner which
means to refine. Since oils are generally less dense than aqueous solutions, the extract will likely
rise to the top of the settler.
Mixers-settlers can also be arranged as a countercurrent cascade of stages, as shown in the figure
above.
TRIANGLE DIAGRAMS
Many extraction systems involve three components which form two immiscible liquid
phases. For example, acetone and water are miscible in all proportions. When an organic
solvent, methyl isobutyl ketone (MIK) is added to a water-rich acetone solution, two phases are
formed: a water rich phase and an MIK-rich phase, both phases containing all three components.
The composition of the two phases can be graphically represented by means of tie lines on a
“triangle diagram.”
Although gas absorption also involved three components, only two of them appeared in
either phase (e.g. acetone in air or acetone in water).♣ Then we could represent the equilibrium
curve on an xy-diagram as in distillation of binary mixtures. By contrast, L/L extraction has
three components in both phases. This means that two independent mole fractions must be
specified in order to know the composition of either phase. A simple xy diagram is no longer
sufficient to represent the phase diagram.
♣ Most of our analysis of gas absorption was based on the assumption of “one-transferrable component” which was
present in both product phases. The gas and liquid phases each contained a second component which was different
for the two phases, for a total of three components. For example, liquid mixtures of NH3 and H2O were separated
by contact with air. The assumption of one-transferrable component means that the mole fraction of H2O in the air
phase and the mole fraction of air in H2O phase are negligible; NH3 is the transferrable component present in both
phases.
x A + xB + xC =
1
Thus each point inside the triangle represents the composition of a possible three-component
mixture. A point at one of the three corners represents a pure component, while the normal
distance between an interior point and the side opposite the corner for pure A represents the mole
fraction or mass fraction of component A in the mixture.
The phase diagram for a mixture of acetone,♣ methyl isobutyl ketone (MIK) and water at
25°C is shown above. The region below the curve represents compositions for which two
immiscible liquid phases form, while the region above the curve are compositions for which only
one phase forms. The nearly horizontal lines inside the lower region represent tie lies connecting
the compositions of two liquid phases that co-exist at equilibrium.
In this system, acetone is completely miscible in both water and MIK, but MIK is only
slightly soluble in water and water is only slightly soluble in MIK. This is an example of a Type
I triangle diagram: the one-phase region completely includes two of the three sides of the
triangle.
Notice in the Type I triangle diagram that as the fraction of acetone gets higher, the tie lines
get shorter. Eventually, the tie line becomes infinitesimally short, such that the composition of
the two immiscible phases become equal at the
plait point — point on the triangle diagram at which the composition of the two immiscible
phases are equal.
A very different looking triangle diagram results when one of the liquids is not completely
miscible with either of other two. The diagram below is for the system aniline/heptane/methyl-
cyclohexane.♦ This is an example of a Type II triangle diagram: the one phase region
completely includes only one of the three sides of the triangle.
NOMENCLATURE
Under normal operation we contact the “feed” stream
with a “solvent” stream, producing two immiscible liquid Solvent Extract
phases which are separated to form two streams. The two MIK
product streams are called the “extract” and the “raffinate.”
raffinate — product stream which is leaner in the solvent (what remains of the feed after the
solute has been extracted)
RECTANGULAR DIAGRAMS
While triangle diagrams are
commonly used to present the
phase diagram for three-compo-
nent mixtures, the geometry of an
equilateral triangle complicates
construction. Three-component
phase behavior can also be
presented using a rectangular
diagram, an example of which is
shown at right for a ternary
system of Type I.
xA + xC = 1 or xB = 0
where the subscripts denote components. Lines parallel to the hypotenuse correspond to
if the line lies above the hypotenuse. Clearly, the latter situation is impossible so that points
outside the right triangle are extraneous.
M = total mass flowrate of the output stream (ignoring the phase distribution)
xAM, xCM = mass fraction of components A,C in the stream having flowrate M
Step 2: Locate point M on a straight line connecting L and V using the lever-arm rule:
L
VM = (214)
LM V
where LM or VM denote the length of the lines connecting points L and M or V and M,
respectively.
L+V=M (215)
L x AM − y A
= (218)
V x A − x AM
L xCM − yC
= (219)
V xC − xCM
x AM − y A xCM − yC
Equating (218) and (219): =
x A − x AM xC − xCM
♥ This proof is similar to that for Lever-Arm Rule for Heat of Dilution on page 53.
xCM − yC xC − xCM
= ≡ tan θ
x AM − y A x A − x AM
slope of line VM slope of line LM
where θ is the angle measured between line segment LV and the horizontal. Since line segments
LM or VM have the same slope and share a common point, the point M must lie on line segment
LV. Moreover, by dividing numerator and denominator of (218) by cosθ, we convert the lengths
of the adjacent sides of the two right triangles into the lengths of the hypotenuses:
( x AM − yA )
=
L
= cos θ VM
V ( x A − x AM ) LM
cos θ
where the overbar quantities denote the length of the line. This
last relationship proves (214). Notice that the two lengths could
be measured on an actual graph with a ruler. The units used for
length cancel out.
Comment: While we usually define the “lever-arm rule” as (214), in applications, it is often
more convenient to use (218) or (219). The three equations give the same value for L/V.
SINGLE-STAGE EXTRACTOR
If point M in the previous example lies in the 1-phase region of the phase diagram, then we
are done because we know the composition (xAM, xCM) of the single phase leaving the mixer. If
the point M lies in the two-phase region and the two phases are allowed to reach equilibrium and
separate in the settler, we will have (not one, but) two product streams. The additional equations
required to solve for the extra unknowns associated with the second product stream can be
obtained by assuming the two product streams are in phase-equilibrium with each other; in other
words, the mixer-settler unit performs as an ideal stage.
Let’s re-label the two input and two output streams as in the
figure at right. In a typical extraction problem, stream L0 is an stage
aqueous solution containing a solute we are trying to extract by #1
contact with an organic solvent. We will use the following
convention when referring to components:
A = solute
B = water
C = organic solvent
So let’s assume that the second feed stream (labelled V2) is pure solvent. Let’s try to determine
the composition of the two product streams.
EXAMPLE #2: Given: L0, xA0, xC0 = 0, V2, yA2 = 0 and yC2 = 1
Step 5: The flowrates L1 and V1 can be determined using the reverse of the procedure for
determining the mixing point M in Example #1. In particular, (214) is used
L1 V1M
= and L1 + V1 = L0 + V2
V1 L1M
Comment #1: Essentially the mixing point gives the composition of a hypothetical single-phase
mixture formed from the two feed streams. According to a mass balance, this is the same
mixture which would be formed by mixing the two product streams. In other words, M is the
mixing point for either L0 and V2 or L1 and V1.
MULTI-STAGE EXTRACTOR
If we are satisfied with the amount of solute we have been able to extract from the solution in
the single-stage extractor, then we are done. If we desire to extract more of the solute (i.e. we
wish to achieve a smaller value of xA1) then multiple stages might be used. Let’s specify an
arbitrarily small value for xAN, the concentration of solute remaining after N stages.
EXAMPLE #3: Given: L0, xA0, xC0 = 0, VN+1, yAN+1 = 0, yCN+1 = 1 and xAN
Solution: The two output streams LN and V1 involve a total of six unknowns (two mole fractions
plus the flowrate for each stream). One of those mole fractions is given (xAN) and we can
perform a total of three independent mole balances (when there are three components in the
mixture). The remaining two constraints arise from the fact that LN and V1 are products of a
equilibrium mixer-settler stage (but not the same stage, so they are not at opposite ends of the
same tie-line). Therefore the two points LN and V1 must lie somewhere on the equilibrium curve
(the curve separating the 1-phase region from the 2-phase region).
This will also be the mixing point for the two product streams LN and
V1. Even though these streams are not coming from the same stage,
they contain the same material as in the two feed streams and
therefore, they must have the same mixing point.
Assuming all the stages in the cascade are equilibrium stages, we also know that the points
LN and V1 lie somewhere on the equilibrium curve, but not necessarily on opposite ends of the
same tie line. Since these streams are not coming from the same stage, they don’t have to be in
equilibrium with each other.
Step 3: Using the known value of xAN, locate the point LN on the lower portion of the
equilibrium curve (i.e. on the raffinate portion of the curve). Read xCN.
Step 4: Draw a straight line connecting LN and M; then extend this line until it crosses the
equilibrium curve again.
Step 5: The intersection of this line with the equilibrium curve is point V1. Read the coordinates:
(yA1, yC1).
Step 6: Calculate the flowrates LN and V1 using the inverse lever-arm rule:
LN M V1
=
V1M LN
This gives us the ratio of the two unknown flowrates. We also know their sum:
Knowing the ratio and the sum, we can solve simultaneously for the two unknowns: LN and V1.
To determine the concentration and flowrates of the remaining streams internal to the
cascade (Li and Vi+1 for i = 1, 2 ... N-1), we need a convenient method for representing mass
balances graphically. The “mixing-point method” we have used above to handle mass balances
cannot be easily extended for this purpose: since we do not know both input flowrates to Stage
#1 or both output flowrates, the mixing point for Stage #1 cannot be determined using the lever-
arm rule. A trial-and-error would be required to determine the mixing point for each separate
stage (very tedious).
An easier graphical method for handling the mass balances is to use “operating lines”.
Operating lines on a triangle diagram connect two points representing adjacent streams. (e.g. Li
and Vi+1). It turns out that all such operating lines (when extended in either direction) share a
common point, called the difference point, which can be calculated in advance.
Difference-Point Method
Consider a total mass balance around some arbitrary subset of stages (denoted by the blue
line in the figure above) in the cascade. At steady state, the rate of total mass entering the system
equals the rate leaving, or:
L0 + Vi+1 = Li + V1
which can be rearranged to obtain difference between the flowrates of two “adjacent” streams in
the cascade
This can be repeated for any subset of stages in the cascade (i.e. i = 1,2 ... N) and the difference
in the flowrates of adjacent streams has exactly the same value:
L0 – V1 = L1 – V2 = … = LN – VN+1 ≡ ∆ (220)
Once again the difference the mass flowrates of component A in any two adjacent streams is
exactly the same, regardless of how many stages are included in the subset. For i = 1,2 ... N:
We denote this difference as ∆xA∆: the product of ∆ and a fictitious mass fraction which we call
xA∆.♦ We say “fictitious” because this mass fraction might be a negative number: it does not
have to lie between zero and unity.
The value of this fictitious mass fraction can be calculated by dividing (221) by (220) for n=1
(chosen because we know the flowrates and compositions of the two adjacent streams at this end
of the cascade from Example #3):
L x −V y L x −V y
x A∆ = 0 A0 1 A1 and xC ∆ = 0 C 0 1 C1 (222)
L0 − V1 L0 − V1
An equation similar to this can be obtained for component C, which leads to the second equation
above. The point (xA∆, xC∆) is called the difference point and is denoted by ∆.♠ All the
operating lines pass through this common point. In the example below, we see that the location
of this point can also be determined graphically without any calculations.
Solution: recipe
Step 6: Determine the tie line passing through V1. Locate L1 as the other end of this tie line.
Step 8: Locate V2 as the point where this operating line crosses the (upper part of) the
equilibrium curve
♦ We could just as easily have called this y . The main idea is that (x , x ) [or (y , y )] denotes a unique
A∆ A∆ C∆ A∆ C∆
point on the rectangle diagram.
♠ Geankoplis instead calls this the “operating point.”
Step 9: Determine the tie line passing through V2. Locate L2 as the other end of this tie line.
Step 10: Repeat steps 7-9 until the desired composition is obtained for LN.
Proof (of difference point method): We need to show that the slope of the line drawn from
L0 to V1 has the same slope as a line drawn from V1 to ∆. Only then can the points L0, V1 and ∆
lie on the same straight line. The equality of these two slopes requires
x A0 − y A1 x A0 − x A∆
=
xC 0 − yC1 xC 0 − xC ∆
L x −V y
x A0 − 0 A0 1 A1
x A0 − y A1 L0 − V1
=
xC 0 − yC1 L x −V y
xC 0 − 0 C 0 1 C1
L0 − V1
x A0 − y A1 L0 x A0 − V1 x A0 − L0 x A0 + V1 y A1
=
xC 0 − yC1 L0 xC 0 − V1 xC 0 − L0 xC 0 + V1 yC1
− V1 x A0 + V1 y A1
=
− V 1 xC 0 + V1 yC1
− x + y A1
= A0
− xC 0 + yC1
x − y A1
= A0
xC 0 − yC1
After cancellation, we end up with an identity. For two line segments having a common point to
have the same slope, these two segments must be part of the same straight line. Thus L0, V1 and
∆ lie on the same straight line. Recalling that (222) was
obtained by dividing (221) by (220), we could easily
generalize this proof to show that Ln, Vn+1 and ∆ lie on
the same straight line.
curve in the McCabe-Thiele method). As long as the slope of the tie line (dashed in the figure at
right) is more negative than the slope of the operating line for that step (OLs are solid lines), we
make progress in reducing the concentration from xAn-1 to xAn.
Should the slope of the operating line ever become equal to that of the tie line, we would step
back and forth endlessly without making progress. In effect, we would have reached a “pinch”.
Then even an infinite number of stages would not allow us to achieve our goal of reducing xA to
xAN.
This new ∆ point will make all the operating lines steeper
than in Example #4. If any of them becomes parallel to a
tie line, we have a pinch. This happens whenever the
solvent flowrate is at or below the minimum.
Recipe:
line connecting LN and VN+1, this intersection might occur either above VN+1 or below LN.
Step 5: The lowest point of intersection of these tie lines with the line in Step 2 is the ∆ point
corresponding to the minimum solvent flowrate.
Step 6: Connect ∆min with L0. Intersection with the extract portion
of the equilibrium curve is V1.
L M
VN + 1 = 0
VN+1M 0
L
Solution: drawing elements mentioned in the description below refer to the full page graph on
the next page.
Steps 1&2: The tie lines given as part of the equilibrium phase behavior are shown in cyan.
Other tie lines were interpolated using Mathcad. The heavy black line is the operating line
drawn through LN and VN+1.
Step 3: the heavy green line is the tie line which passes through L0.
Step 4: about ten other tie lines, equally spaced between L0 and LN are drawn in red.
240
230 0
∆ min
220
210
200
190
180
170
160
150
140
130
Ether wt%
120
110
VN + 1 100
100
90
V1
80
70
60
M
50
40
30
20
10
0 LN L0
10 0 10 20 30
Acetic acid wt%
Step 5: of the ten red tie lines and the green one, it is the green tie line which intersects the black
operating line closest to VN+1. This closest intersection point represents the critical ∆.
Typically (as in this problem) it is the tie line through L0 which determines the minimum solvent
flow rate.
Step 6: The intersection of the green line with equilibrium curve is V1 corresponding to
minimum solvent flowrate.
Step 7: Next we draw straight lines connecting the two feed points L0 and VN+1 and the two
product points LN and V1. Their intersection is the mixing point for either pair of streams.
Step 8: using the lever-arm rule, we can now calculate the minimum solvent flowrate by
measuring the lengths of the two line segments:
VN + 1 L0 M lb of ether
= = 1.283
L0 VN + 1 M lb of aqueous solution