0% found this document useful (0 votes)
234 views10 pages

Distillation Column Design With Matlab

This project is worked by students of Chemical Engineering Department at Addis Ababa Science and Technology University, Ethiopia

Uploaded by

Theødřøš Äb
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
234 views10 pages

Distillation Column Design With Matlab

This project is worked by students of Chemical Engineering Department at Addis Ababa Science and Technology University, Ethiopia

Uploaded by

Theødřøš Äb
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/ 10

Addis Ababa Science and Technology University

College of Biological and chemical Engineering


Department of Chemical Engineering

Course Name: Computer Aided Process Design and


Simulation
Course Code: Cheg5192

Prepared By
Name ID No.
Yishak Gizaw……………...TETS152/07
Eskedar Workineh…………TETS113/07
Rediet Mekonen……………TETS141/07
Kalkidan Mosse…………….TETS121/07

Submitted to: Baru D.(MSc)


Submission date: 18/12/2018
% of Hexanol and Pentanol Mixtures
EDm=[0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9
0.95 1;0 0.0844619 0.1632757 0.2369315 0.3058709 0.37049 0.4311458 0.4881602
0.5418238 0.5923998 0.6401265 0.6852204 0.727878 0.8065858 0.842787 0.8429481
0.8775032 0.9103758 0.9416804 0.9715228 1];[r,c]=size(EDm);
ED=ones(r,c);
MtP=88.15;
MtH=102.162;
for i =1:r;
for j= 1:c;
ED(i,j)=(EDm(i,j)/MtP)/((EDm(i,j)/MtP)+((1-EDm(i,j))/MtH));
end
end
ys=ED(2,1:length(ED));
xs=ED(1,1:length(ED));
c= polyfit(xs,ys,6);
figure('Name', 'McCabe-Thiele Method');
hold on
eq=plot(0:0.1:1,polyval(c,0:0.1:1));
op=line([0 1],[0 1]);
set(eq,'color','b');
set(op,'color','r');
xlabel('Pentanol in liquid form');ylabel('Pentanol vapour form');
title('McCabe-Thiele Method method');
xd=0.95;
xb=0.05;
xs=[];ys=[];
xs(1)=xd;
ys(1)=xd;
f=@(x,y)(c(1)*x^6+c(2)*x^5+c(3)*x^4+c(4)*x^3+c(5)*x^2+c(6)*x-y);
i=1;
while xs(i)>xb
xs(i+1)=fzero(f,0.4,[],ys(i));
line([xs(i) xs(i+1)],[ys(i) ys(i)],'color', 'g');
ys(i+1)=xs(i+1);
if xs(i+1)>xb
line([xs(i+1) xs(i+1)],[ys(i) ys(i+1)],'color', 'g');
end
i=i+1;
end
disp(i) % displays the value of variables i without printing the variable name.

Results from Command Window

 The number of Stages are Twelve (12) with reboiler


%Carbod disulfide and Carbon tetrachloride ED Mixtures

EDm=[0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9
0.95 1;0 0.13005 0.240264 0.334331 0.415193 0.485199 0.54624 0.599845 0.647261
0.689519 0.727474 0.761849 0.793267 0.822276 0.849376 0.875038 0.899731 0.92395
0.948251 0.97331 1];[r,c]=size(EDm);

ED=ones(r,c);

MtCD=76.145;

MtCTC=153.82;

for i =1:r;

for j= 1:c;

ED(i,j)=(EDm(i,j)/MtCD)/((EDm(i,j)/MtCD)+((1-EDm(i,j))/MtCTC));

end

end

ys=ED(2,1:length(ED));

xs=ED(1,1:length(ED));

c= polyfit(xs,ys,6);

figure('Name', 'McCabe-Thiele');

hold on

eq=plot(0:0.1:1,polyval(c,0:0.1:1));

op=line([0 1],[0 1]);

set(eq,'color','b');

set(op,'color','r');

xlabel('Carbon disulfide in liquid form');ylabel('Carbon disulfide in vapor form');

title('McCabe-Thiele Method');

xd=0.95;

xb=0.05;
xs=[];ys=[];

xs(1)=xd;

ys(1)=xd;

f=@(x,y)(c(1)*x^6+c(2)*x^5+c(3)*x^4+c(4)*x^3+c(5)*x^2+c(6)*x-y);

i=1;

while xs(i)>xb

xs(i+1)=fzero(f,0.4,[],ys(i));

line([xs(i) xs(i+1)],[ys(i) ys(i)],'color', 'g');

ys(i+1)=xs(i+1);

if xs(i+1)>xb

line([xs(i+1) xs(i+1)],[ys(i) ys(i+1)],'color', 'g');

end

i=i+1;

end Results from Command Window

disp(i)

 The number of Stages are Eight (8) with reboiler


%Diphenyl methane and P-tetra butyl phenol mixture

EDm=[0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9
0.95 1;0 0.081464 0.157874 0.229647 0.297159 0.390748 0.420723 0.477361 0.530917
0.581619 0.629679 0.675287 0.718617 0.759831 0.799074 0.836482 0.872178 0.906277
0.938883 0.970094 1];[r,c]=size(EDm);

ED=ones(r,c);

MtBP=150.211;

MtDP=186.21;

for i =1:r;

for j= 1:c;

ED(i,j)=(EDm(i,j)/MtBP)/((EDm(i,j)/MtBP)+((1-EDm(i,j))/MtDP));

end

end

ys=ED(2,1:length(ED));

xs=ED(1,1:length(ED));

c= polyfit(xs,ys,6);

figure('Name', 'McCabe-Thiele');

hold on

eq=plot(0:0.1:1,polyval(c,0:0.1:1));

op=line([0 1],[0 1]);

set(eq,'color','b');

set(op,'color','r');

xlabel('P-Tetrabutyl phenol in liquid form');ylabel('P-Tetrabutyl phenol in vapor form');

title('McCabe-Thiele Method');

xd=0.95;

xb=0.05;
xs=[];ys=[];

xs(1)=xd;

ys(1)=xd;

f=@(x,y)(c(1)*x^6+c(2)*x^5+c(3)*x^4+c(4)*x^3+c(5)*x^2+c(6)*x-y);

i=1;

while xs(i)>xb

xs(i+1)=fzero(f,0.4,[],ys(i));

line([xs(i) xs(i+1)],[ys(i) ys(i)],'color', 'g');

ys(i+1)=xs(i+1);

if xs(i+1)>xb

line([xs(i+1) xs(i+1)],[ys(i) ys(i+1)],'color', 'g');

end

i=i+1;

end

disp(i)
Results from Command Window

 The number of Stages are Thirten (13) with reboiler


 Equilibrium data for Hexanol and Pentanol (Excel)
xs 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
ys 0 0.08446 0.16328 0.23693 0.30587 0.37049 0.43115 0.48816 0.54182 0.5924 0.64013
xs 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1
ys 0.68522 0.72788 0.76828 0.80659 0.84295 0.8775 0.91038 0.94168 0.97152 1

 X and Y values for theoretical number of stage calculations


x 0.95 0.95 0.915 0.915 0.855 0.855 0.7675 0.7675 0.649 0.649 0.5075 0.5075 0.3675
value
y 0 0.95 0.95 0.915 0.915 0.855 0.855 0.7675 0.7675 0.649 0.649 0.5075 0.5075
value
x 0.3675 0.2475 0.2475 0.158 0.158 0.0975 0.0975 0.095 0.095 0.05 0.05 0.05
value
y 0.3675 0.3675 0.2475 0.2475 0.158 0.158 0.0975 0.0975 0.095 0.095 0.05 0
value

Result: The number of Stages are Eleven (11) with reboiler

 It has some difference with the matlab result

0.9

0.8

0.7
Pentanol In vapour form

0.6

0.5

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Pentanol in liquid form


Equilibrium data for Tetra butyl phenol and Diphenyl (Excel)

xs 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
ys 0 0.081464 0.157874 0.229647 0.297159 0.360748 0.420723 0.473736 0.530917 0.581619 0.629679
xs 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1
ys 0.675287 0.718617 0.759831 0.799074 0.836482 0.872178 0.906277 0.938883 0.970094 1

Theoretical number of stage calculation values


x value 0.95 0.95 0.92 0.92 0.87 0.87 0.8 0.8 0.7 0.7 0.58 0.58 0.45
y value 0 0.95 0.95 0.92 0.92 0.87 0.87 0.8 0.8 0.7 0.7 0.58 0.58
x value 0.45 0.33 0.33 0.23 0.23 0.15 0.15 0.098 0.098 0.056 0.056 0.039
y value 0.45 0.45 0.33 0.33 0.23 0.23 0.15 0.15 0.098 0.098 0.056 0.056

Result: The number of Stages are Thirteen (13) with reboiler

Theoretical stage calculation for Tetrabutyl phenol and Di


1
phenol
0.9

0.8

0.7
% Tetrabutyl phnol vapour form

0.6

0.5

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

%Tetrabutyl phenol liquid form


Equilibrium data for Carbon disulfide and Carbon tetrachloride(Excel)
xs 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45

ys 0 0.13005 0.240264 0.334331 0.415193 0.485199 0.54624 0.599845 0.647261 0.689519

xs 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1

ys 0.727474 0.761849 0.793267 0.822276 0.849376 0.87503 0.899731 0.92395 0.948251 0.97331 1

Number of stage drawing values


x value 0.95 0.95 0.9025 0.9025 0.805 0.805 0.62
y value 0 0.95 0.95 0.9025 0.9025 0.805 0.805
x value 0.62 0.37 0.37 0.17 0.17 0.0675 0.0675 0.025
y value 0.62 0.62 0.37 0.37 0.17 0.17 0.0675 0.0675

Result: The number of Stages are Eight (8) with reboiler

0.9

0.8

0.7
Carbon disulfide vapour form

0.6

0.5

0.4

0.3

0.2

0.1

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Carbon disulfide vapour form

You might also like