DA-2
NAME – RISHAV KUMAR
REG NO – 17BCM0026
SLOT – B1
COURSE- COMPUTATIONAL METHOD
IN PROCESS ENGINEERING
COURSE CODE- CHE 3001
DATE – 19-09-19
Q) A mixture of benzene and toluene containing 60% moles of benzene is to be
separated to give a product of 95% moles of benzene and a bottom product containing
of 10% moles benzene. The feed enters a column at its bubble point. It is proposed to
operate the column with reflux ratio 2.5. Find the number of theoretical plates needed.
Given equilibrium data:
x 0 0.1 0.3 0.5 0.6 0.8 1
y 0 0.21 0.5 0.7 0.77 0.9 1
MATLAB CODE:
global y
for i=1:11
y=0.1*(i-1);
ye(i)=0.1*(i-1);
xe(i)=fzero('equilib',0.5);
end
% Distillate, Bottom and Feed mole fractions are equal to 95 % mol, 10 % mol
% and 60 % mol respectively.
xd=0.95;
xb=0.1;
zf=0.6;
% Reflux ratio is equal to 2.5
R=2.5;
q=0.85;
yi=(zf+xd*q/R)/(1+q/R);
xi=(-(q-1)*(1-R/(R+1))*xd-zf)/((q-1)*R/(R+1)-q);
figure(1);
hold on;
axis([0 1 0 1]);
% plotting operating and feed lines and equilibrium curve
plot(xe,ye,'r');
set(line([0 1],[0 1]),'Color',[0 1 0]);
set(line([xd xi],[xdyi]),'Color',[1 0 1]);
set(line([zf xi],[zfyi]),'Color',[1 0 1]);
set(line([xb xi],[xbyi]),'Color',[1 0 1]);
% Stepping off stages
% Rectifying section
i=1;
xp(1)=xd;
yp(1)=xd;
y=xd;
while (xp(i)>xi),
xp(i+1)=fzero('equilib',0.5);
yp(i+1)=R/(R+1)*xp(i+1)+xd/(R+1);
y=yp(i+1);
set(line([xp(i) xp(i+1)],[yp(i) yp(i)]),'Color',[0 0 1]);
if (xp(i+1)>xi) set(line([xp(i+1) xp(i+1)],[yp(i) yp(i+1)]),'Color',[0 0 1]);
end
i=i+1;
end
% Stripping section
SS=(yi-xb)/(xi-xb);
yp(i)=SS*(xp(i)-xb)+xb;
y=yp(i);
set(line([xp(i) xp(i)],[yp(i-1) yp(i)]),'Color',[0 0 1]);
while (xp(i)>xb),
xp(i+1)=fzero('equilib',0.5);
yp(i+1)=SS*(xp(i+1)-xb)+xb;
y=yp(i+1);
set(line([xp(i) xp(i+1)],[yp(i) yp(i)]),'Color',[0 0 1]);
if (xp(i+1)>xb) set(line([xp(i+1) xp(i+1)],[yp(i) yp(i+1)]),'Color',[0 0 1]);
end
i=i+1;
end
hold off;
GRAPH:
RESULT:
Hence from the graph obtained we determined that the number of theoretical trays
required are 9.