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
