Separation - Absorption
Separation - Absorption
STRIPPING/DESORPTION
• solute is removed from a liquid by contacting it
with a gas
GAS-LIQUID EQUILIBRIUM
Henry’s laws: p = Hx or yA = H’xA
A A
where
pA = partial pressure of component A (atm)
H= Henry’s law constant (atm/mol fraction) in Appendix A.3-18
H’ = Henry’s law constant (mol fraction gas/mol fraction liquid) = H/P
xA = mole fraction of component A in liquid
yA = mole fraction of component A in gas = pA/P
P= total pressure (atm)
SINGLE-STAGE EQUILIBRIUM CONTACT
yA1 yA2
xA0 xA1
Liquid & gas phases are brought into contact and separated
Long enough for equilibrium
Gas phase – solute A & inert gas B
Liquid phase – solute A & inert liquid/solvent C
Total material balance: L0 + V2 = L1 + V1
Balance on A: L0xA0 + V2yA2 = L1xA1 + V1yA1
Balance on A can also be written as
x y x y
A0 A2 A1 A1
L’ 1- x + V’ 1- y = L’ 1- x + V’ 1- y
A0 A2 A1 A1
where L’ = moles inert C
V’ = moles inert B
Example
Example 10.3-1
V1 V2 = 100 kmol/h
1 atm
yA1 yA2=0.2
L0 = 300 kmol/h 293K L1
xA0 = 0 xA1
At 293K, Henry’s law constant from App. A.3-18 = 0.142 x 104 atm/mol frac.
yA1 = H’xA1 = (H/P)xA1 = (0.142 x 104/1.0)xA1 (2)
Substitue yA1 = 0.142 x 104xA1 into eq. (1):
0 0.2 x 0.142 x 104 x
A1 A1
300 1- 0 + 80 1- 0.2 = 300 1- x + 80 1- 0.142 x 104 x
A1 A1
L0 + VN+1 = LN + V1
Balance on A:
operating line
COUNTERCURRENT MULTIPLE-CONTACT STAGES
Graphical determination of N:
Dilute Concentrated
Example 10.3-2
VN+1 = 30 kmol/h
yAN+1 = 0.01
L0 = 90 kmol/h
xA0 =0
Example 10.3-2
VN+1 = 30 kmol/h
yAN+1 = 0.01
L0 = 90 kmol/h
xA0 =0
300K, 101.3 kPa A = acetone, B = air, C = water 90% acetone absorb
Acetone in VN+1 = 0.01(30) = 0.3 kmol/h
Acetone in LN = 0.9(0.3) = 0.27 kmol/h
Balance of Acetone in V1 = 0.03 kmol/h
Entering air = 30 – 0.3 = 29.7 kmol/h
Entering water= 90 kmol/h
V1 = 29.7 + 0.03 = 29.73 kmol/h
yA1 = 0.03/29.73 = 0.00101
LN = 90 + 0.27 = 90.27 kmol/h
xAN = 0.27/90.27 = 0.003
Example 10.3-2
V1 = 29.73 kmol/h VN+1 = 30 kmol/h
yA1 = 0.00101 yAN+1 = 0.01
L0 = 90 kmol/h LN = 90.27 kmol/h
xA0 =0 xAN = 0.003
Kremser equations – valid only when operating & equilibrium lines are straight
Absorption: y - mx 1 1
ln N 1 0 1
y - mx
A A
N 1 0
ln A
where
m = slope of equilibrium line
A = absorption factor = Aav. = √(A1AN)
A1=L0/(mV1)
AN = LN/(mVN+1)
y -y
N N 1 1
When A = 1 y 1 - mx 0
ANALYTICAL/ KREMSER EQUATION
V1 VN+1
yA1 yAN+1
L0 LN
xA0 xAN
yN 1
Stripping: x0 -
m
ln y N 1 1 A A
xN -
N m
1
ln
where A
m = slope of equilibrium line
A = absorption factor = Aav. = √(A1AN)
A1=L0/(mV1)
AN = LN/(mVN+1)
x0 - xN
N y
When A = 1 xN - N 1
m
Example 10.3-3
xNA = 0.0218
Example 10.6-3 Number of theoretical trays = 4
xNA = 0.0218
yN+1A y - mx 1 1
ln N 1 0 1
y - mx
A A
N 1 0
ln A
A = absorption factor = Aav. = √(A1AN)
A1=L0/(mV1) & AN = LN/(mVN+1)
V1 = 98.02 kmol/h, VN+1 = 100 kmol/h,
L0 = 88.86 kmol/h, LN = 90.84 kmol/h
y1A
A1=L0/(mV1) = 88.86/[(0.68)(98.02)] = 1.333
x0A
xNA = 0.0218 AN = LN/(mVN+1) = 90.84/[(0.68)(100)] = 1.336
Equilibrium line: y = 0.68x A = √(A1AN) = √[(1.333)(1.336)] = 1.335
xA0 = 0, yAN+1 = 0.022, yA1 = 0.002244, xNA = 0.0218
0.022- (0.68)0 1 1
1
ln
0.00244- (0.68)0 1.335 1.335
N
4.04
ln 1.335
PACKED TOWER
PRESSURE DROP & FLOODING IN PACKED TOWERS
where:
G = superficial gas velocity (ft/s) = GG/G
GL = liquid mass velocity (Ibm/s.ft2)
G = gas density (Ibm/ft )
3
Accuracy = ± 11%
Fp = packing factor (ft-1) in Table 10.6-1
v = kinematic viscosity (centstokes) = L/(L/62.4)
L= liquid viscosity (cp)
PRESSURE DROP IN STRUCTURED PACKING
FLOODING PRESSURE DROP IN PACKED & STRUCTURED
PACKINGS
Pflood = 0.115FP0.7
where:
Pflood = pressure drop at flooding (in. H2O/ft height of packing)
Fp = packing factor (ft-1) in Table 10.6-1
Conversion: 1 in H2O/ft height packing= 83.33 mm H2O/m height of packing
FP from 9 - 60
Accuracy = ± 10-15 %
FP 60 or higher, Pflood = 2.0 in. H2O/ft (166.7 mm H2O/m)
PRESSURE DROP & TOWER DIAMETER IN PACKED &
STRUCTURED PACKINGS
Procedure:
1. From the type of packing used, get FP from Table 10.6-1
2. Determine Pflood from Pflood = 0.115FP0.7 or Pflood =2 in H2O/ft packing
height when FP ≥ 60
3. Calculate flow parameter using the gas and liquid flows in the bottom of
the tower. Using Fig. 10.6-5 or 10.6-6, read off the capacity parameter
4. Calculate G from the capacity parameter which is equal to GG/G where
GG = maximum value of gas mass velocity at flooding
5. Using a given % of the floodingGG , obtain new GL and GG based on the
given liquid-to-gas ratio GL/GG
6. Calculate the cross-sectional area of the tower (D2/4) from the given gas
flow rate and hence, the diameter of the tower
7. Calculate the total flow rates of the outlet and inlet water assuming all the
solute is absorbed
PACKED TOWERS FOR ABSORPTION
Balance on A :
x 2 y
1
x
1
y
2
L’ 1 - x + V’ 1 - y = L’ 1 - x + V’ 1 - y
2 1 1 2
Operating line y x y
x 1 1
L’ 1 - x + V’ 1 - y = L’ 1 - x + V’ 1 - y
1 1
Absorption Stripping
LIMITING & OPTIMUM L’/V’ RATIOS
Balance on A: x
2
y
1
x
1
y
2
L’ 1 - x + V’ 1 - y = L’ 1 - x + V’ 1 - y
2 1 1 2
Absorption Stripping
At point P, liquid flow L’ = L’min & x1 = x1max
x y x y
2 1 1max 2
L’min 1 - x + V’ 1 - y =L’min 1 - x + V’
1- y
2 1 1max 2
L x1 x 2
L x1 x 2
z z
k'x aS x i x
K'x aS x * x
M M
y - y i
y -y* 1
ln y y *1 / y 2 y * 2
ln y 1 y i1 / y 2 y i2
M
M 1
x i - x
xi1 x1 xi2 x2 x *1 x1 x *2 x 2
x * -x
ln x i1 x1 / x i2 x 2
ln x *1 x1 / x * 2 x 2
M M
av
∫ y2 y y *
L1 - x x1
dx L1 - x *M
∫
x1
∫ dx
z iM
x x z
k'x aS 1 - x x2 i
K'x aS1 - x x2 x * x
av av
1 y 1 y
1 y * 1 y
1- y *
ln1 y * /1 y
1 - y i
M
iM
ln1 y i /1 y
1 x 1 xi
1 x 1 x *
1- x 1- x *
iM
ln1 x /1 x i
ln1 x /1 x *
M
V’ = 13.65 kmol/h
y1A = 0.026 L1
V’ = 13.65 kmol/h
y1A = 0.026 L1
Material balance on A:
x A2 y
A1
x
A1
y
A2
L’ 1 - x + V’ 1 - y = L’ 1 - x + V’ 1 - y
A2 A1 A1 A2
0.026 x 0.005
0
1
45.36 1 - 0 + 13.65 1- 0.026 = 45.36 1 - x + 13.65 1- 0.005
1
x1 = 0.00648
Example 10.6-4 V2
L’ =45.36 kmol/h
y2A = 0.005
x2=0
T = 293K
P = 101.3 kPa
k’ya = 3.78 x 10-2 kmol A/s.m3.mol frac.
S = 0.186m2
k’x a= 6.16 x 10-2 kmol A/s.m3.mol frac.
V’ = 13.65 kmol/h L1
y1A = 0.026 x1 = 0.00648
1. Plot operating line & equilibrium line
Example 10.6-4 k’ya = 3.78 x 10-2 kmol A/s.m3.mol frac.
k’x a= 6.16 x 10-2 kmol A/s.m3.mol frac.
V’ = 13.65 kmol/h, L’ = 45.36 kmol/h , S = 0.186 m2
y1 = 0.026 , x1 = 0.00648, y2 = 0.005, x2 = 0
k'x a/1 - x
2. By trial-and-error, determine yi, xi or y*, x*: 1st trial, Slope = -
k'y a/1 - y
For point 1: (y1 = 0.026 , x1 = 0.00648) 6.16 x 10- 3 / 1 - 0.00648
Slope = - -2
1.60
3.78 x 10 / 1 - 0.026
Example 10.6-4
For point 1:
From the plot, yi1 = 0.0154 , xi1 = 0.013, y*1 = 0.0077
k'x a/1 - x
iM
For 2nd trial, Slope = -
k'y a/1 - y
iM
1 y 1 y
1 0.0154 1 0.026
1 - y i
0.979
iM
ln1 y i /1 y ln1 0.0154 /1 0.026
1 - y 0.979
1 - x 0.993
iM
iM
k'x a/1 - x
iM
For 2nd trial, Slope = -
k'y a/1 - y
iM
-3
6.16 x 10 /0.993
Slope = - -2
1.61
3.78 x 10 /0.979
Since the latest slope and the former slope is approximate close, the values
yi1 = 0.0154 , xi1 = 0.013, y*1 = 0.0077 are accurate enough.
Example 10.6-4
Since the slope at point 2 and point 1 changes little in the tower, the value
of the slope -1.62 from the 1st trial is acceptable. Plotting the slope at point
2 gives yi2 = 0.002, xi2 = 0.0018 and y*2 = 0
Example 10.6-4 k’ya = 3.78 x 10-2 kmol A/s.m3.mol frac.
k’x a= 6.16 x 10-2 kmol A/s.m3.mol frac.
V’ = 13.65 kmol/h, L’ = 45.36 kmol/h , S = 0.186 m2
y1 = 0.026 , x1 = 0.00648, y2 = 0.005, x2 = 0
yi1 = 0.0154 , xi1 = 0.013, y*1 = 0.0077
yi2 = 0.002, xi2 = 0.0018 and y*2 = 0
where
y - y
y1 y i1 y 2 y i2 0.026
0.0154
0.005 0.002
0.00602
ln y 1 y i1 / y 2 y i2 ln 0.026 0.0154 / 0.005 0.002
i M
Example 10.6-4 k’ya = 3.78 x 10-2 kmol A/s.m3.mol frac.
k’x a= 6.16 x 10-2 kmol A/s.m3.mol frac.
V’ = 13.65 kmol/h, L’ = 45.36 kmol/h , S = 0.186 m2
y1 = 0.026 , x1 = 0.00648, y2 = 0.005, x2 = 0
yi1 = 0.0154 , xi1 = 0.013, y*1 = 0.0077
yi2 = 0.002, xi2 = 0.0018 and y*2 = 0
y - y 0.00602
i M
V y1 y 2 (V1 V2 )/2 y 1 y 2
z
k'y aS y y i
k'y aS
y y i
M M
V y1 y 2 (V1 V2 ) y 1 y 2
z
k'y aS y y i
k'y aS
y y i
M M
(3.893x10- 3 3.811x10 3 )/2 0.026 0.005
3.78x10- 2 (0.186) 0.00602
z 1.911 m
HEIGHT & NUMBER OF TRANSFER (HTU & NTU)
Equilibrium & Operating lines = straight
V y1 y 2 V y1 y 2
z
[HG ]NG z
[HOG ]NOG
k'y aS y y i
K'y aS y y *
M M
L x1 x 2
L x1 x 2
z [H L ]N L z [HOL ]NOL
k'x aS x i x
K'x aS x * x
M M
y1
V1 - y y1 dy [H ]N V1 - y *M dy [H ]N
z
k'y aS 1 - y
iM
y2∫
av
y yi G G z
K'y aS1 - y ∫
y2
av
y y* OG OG
L1 - x x1 L1 - x *M x1
dx [H ]N
∫ dx [H ]N
∫
z iM z
k'x aS 1 - x
x2 xi x L L
K'x aS1 - x x2 x * x OL OL
av av
HEIGHT & NUMBER OF TRANSFER (HTU & NTU)
When major resistance to mass transfer is in gas phase, NOG
or NG should be used (absorption)
When major resistance to mass transfer is in liquid phase,
NOL or NL should be used (stripping)
1 y - mx
N ln 1 1 1 2 1
Operating line = straight OG 1 A y 2 - mx 2 A
1
A
N ln A
Equilibrium & Operating lines = straight & not parallel NOG
1 - 1/A
y 1
x2 -
1
m A
Operating line = straight NOL ln 1 A
y
1 A
x1 -
1
m
L
HOL = HL + mV HG
where
Height of tower, z
At 303K and 101.32kPa, from A.3-3, air= 1.666 kg/m3 and = 1.866 x 10-5 kg/m.s
From Table 6.2-1, for Air-CO2 system, at 276K DAB = 0.142 x 10-4 m2/s
1.75
T P
2
Correcting for 303K, DAB303K = DAB276K T P1
1 2
1.75
-4 303
101.32 0.167x10 4
DAB303K = 0.142 x 10
276
101.32
Example 10.8-1
μ 1.866x10 5
NSc = Schmidt number = ρD 0.958
4
AB 1.666 0.167x10
0.5 0.5 0.35
0.226 0.958
4.069 0.5424
HG
0.2426m
1.34 0.660
6.782
0.678
Example 10.8-1 Gx = 4.069 kg/s.m2 and Gy = 0.5424 kg/s.m2.
0.3
G X
0.5
0.357 N
μ
HL
Sc
f
372
6.782
P -3
0.8937x10
0.3
4.069
0.5
0.357 354.3
0.8007x10 -3
H L
0.2306 m
1.34 372
6.782
-3
0.8937x10
Example 10.8-1
Gx = 4.069 kg/s.m2 and Gy = 0.5424 kg/s.m2.
HG = 0.2426 m & HL = 0.2306 m
L
HOL = HL + mV HG
where
m = slope of equilibrium line
From A.3-18, for CO2 at 1 atm, pA = 0.186 x 104xA,
yA = (pA/P)xA = (0.186 x 104/1)xA
L
0.2261
0.2426 0.2322 m
HOL = HL + mV HG
0.2306
0.186x104 (0.01872)