Scribd
Upload
57 views2 pages

Packed Column Ex

This document contains properties and flowrates for a packed column separation process. It provides temperature, pressure, component flowrates and compositions for the feed and outlets of the column. It also lists parameters for the packing material and calculations of various stage efficiencies including NOG (Number of Transfer Units), using a MATLAB code to model the column and determine the NOG.

Uploaded by

Danielle Barker
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
57 views2 pages

Packed Column Ex

This document contains properties and flowrates for a packed column separation process. It provides temperature, pressure, component flowrates and compositions for the feed and outlets of the column. It also lists parameters for the packing material and calculations of various stage efficiencies including NOG (Number of Transfer Units), using a MATLAB code to model the column and determine the NOG.

Uploaded by

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

Cabanas | Delfin | Imperial

ChE 133 Packed Column Exercise, Oct. 05, 2016

Properties Flowrates

Temperature (K) 293.15 In Out

Pressure (atm) 1 Feed (kg/h) 6000


Feed
ky 0.0005 188.6792453
(kmol/h)
SO2
kx 0.0139 15.09433962
(kmol/h)
packing factor (m-1) 170 V' (kmol/h) 173.5849057
gas phase density
1.21 Compositions
(kg/m3)
liquid phase density
1000 In Out
(kg/m3)
liquid phase
0.001 SO2 (v/v) 0.08
viscosity (N-s/m2)
aw (m2/m3) 138 SO2 transferred (kmol/hr) 14.4
SO2
Molecular Weights (kg/mol) 0.954
(%recovery)
SO2 64

Air 29

H2O 18

Feed Gas 31.8

Mole Fractions

y2 0.08 y1 0.003984064

Y2 0.086956522 Y1 0.004

x1 0.002930403 x0 0

X1 0.002939015 X0 0

Y* 0.08 X*max 0.002930403

EC slope 27.3 HOG

(L'/V')min 28.22595652 (1-y)lm 0.957505115

(L'/V') 39.51633913 g 1.450318765

L' 6859.44 S 0.689503593

L'min 4899.6 HL 1.507736315

HG 1.150521934

HOG 37.15451181

(from Matlab
NOG 6.664856292
code)

h (m) 247.6294818
Cabanas | Delfin | Imperial
ChE 133 Packed Column Exercise, Oct. 05, 2016

NOG Matlab Code

clear while abs(ac-ap)>0.00001


clc N=N+3;
MWgas=64*0.08+29*(1-0.08); yv=linspace(a,b,N+1);
V=6000/MWgas;
SO2in=V*0.08; for i=1:N+1
Vp=V-SO2in;
ysv(i)=EC_xtoys(Xtox(OL_YtoX(ytoY(yv(i)))));
end
y2=0.08;
y1=(.046*SO2in)/(.046*SO2in+Vp); h=-(b-a)/N;
x0=0; I=0;
for m=1:N/3
Y2=y2/(1-y2); I=I+(f(yv(3*m-2), ysv(3*m-
Y1=y1/(1-y1); 2))+3*f(yv(3*m-1), ysv(3*m-1))+3*f(yv(3*m),
X0=x0/(1-x0); ysv(3*m))+f(yv(3*m+1), ysv(3*m+1)))*(3*h/8);
end
Ystar2=Y2; ap=ac;
x1max=y2/27.3; ac=I;
X1max=x1max/(1-x1max); end
N_OG=ac
LVmin=(Ystar2-Y1)/(X1max-X0);
LV=1.4*LVmin;
X1=Y2/LV;

x1=X1/(1+X1);

ytoY=@(y) y/(1-y);
OL_YtoX=@(Y) (Y-Y1)/LV;
Xtox=@(X) X/(1+X);
EC_xtoys=@(x) x*27.3;

num=@(y, ys) ((1-y)-(1-ys))/(log((1-y)/(1-


ys)));
den=@(y, ys) (1-y)*(ys-y);
f=@(y, ys) num(y, ys)/den(y,ys);
a=y1;
b=y2;

N=0;
ap=0;
ac=0;
for n=1:2
N=N+3;
yv=linspace(a,b,N+1);
ysv=[];

for i=1:N+1

ysv(i)=EC_xtoys(Xtox(OL_YtoX(ytoY(yv(i)))));
end

h=-(b-a)/N;
I=0;
for m=1:N/3
I=I+(f(yv(3*m-2), ysv(3*m-
2))+3*f(yv(3*m-1), ysv(3*m-1))+3*f(yv(3*m),
ysv(3*m))+f(yv(3*m+1), ysv(3*m+1)))*(3*h/8);
end
ap=ac;
ac=sum(I);
end

You might also like