Scribd
Upload
46 views

AMT Assignment 2 Anup

This document describes the calculations to determine the bubble point and distillation column design for a mixture containing methanol, ethanol, and propanol. It provides the process conditions, assumptions, and steps to calculate vapor-liquid equilibrium using the Antoine equation and generate tridiagonal matrices to solve for mole fractions in each stage of the column. The key outputs are the normalized mole fractions of each component in each stage.

Uploaded by

PRATEEK KUMAR DAS
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
46 views

AMT Assignment 2 Anup

This document describes the calculations to determine the bubble point and distillation column design for a mixture containing methanol, ethanol, and propanol. It provides the process conditions, assumptions, and steps to calculate vapor-liquid equilibrium using the Antoine equation and generate tridiagonal matrices to solve for mole fractions in each stage of the column. The key outputs are the normalized mole fractions of each component in each stage.

Uploaded by

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

Bubble and Dew Point Calculations

fprintf('GIVEN CONDITIONS: \n');


F=input('\n Feed input in the middle stage of the bubble point
mixture(kmol/h):');
m=input('Percentage Composition of Methanol:');
M=(m/100)*F;
fprintf('Methanol input(kmol/h): %d',M)
e=input('\n Percentage Composition of Ethanol:');
E=(e/100)*F;
fprintf('Ethanol input(kmol/h): %d',E)
p=input('\n Percentage Composition of Propanol:');
P=(p/100)*F;
fprintf('Propanol input(kmol/h): %d',P)
Rv=input('\nVapour flow rate leaving Reboiler and at each stage(kmol/h):');
fprintf('Three equilibrium stages in the column \n')
fprintf('A total Condenser \n')
fprintf('A partial reboiler')
fprintf('\n Feed tray is the 3rd tray \n')
D=input('Distillate rate(kmol/h):');
R=input('Reflux rate(kmol/h):');
fprintf('\n ASSUMPTIONS: \n')
fprintf('\n1.Constant Molar Overflow \n');
fprintf('2.Ideal Solution with Raoult law K-values from vapour pressure \n')
fprintf('3.Linear temperature profile for tear variables and liquid flow
rates ')
B=F-D;
fprintf('\n Bottom flow rate:%d',B)
fprintf('\n Methanol in the Distillate \n')
fprintf('Ethanol and Propanol in the Bottoms \n')
LRaF=input('Liquid Rate above the feed(kmol/h):');
LRbF=input('Liquid Rate below the feed(kmol/h):');
Mbp=64.7;
fprintf('\n Boiling point of Methanol(in deg C): %.2f',Mbp)
Ebp=78.2;
fprintf('\n Boiling point of Ethanol(in deg C): %.2f',Ebp)
Pbp=97.2;
fprintf('\n Boiling point of Propanol(in deg C): %.2f',Ebp)
fprintf('\n Temperature Approximations: \n')
T1=64.7;
fprintf('Distillate Temperature(in deg C),T1: %.2f',T1)
T5=(Ebp+Pbp)/2;
fprintf('\n Bottom Temperature(in deg C),T5: %.2f',T5)
d=(T5-T1)/4;
T2=T1+d;
fprintf('\n T2(in deg C): %.2f',T2)
T3=T1+2*d;
fprintf('\n T3(in deg C): %.2f',T3)
T4=T1+3*d;
fprintf('\n T4(in deg C): %.2f \n',T4)
j=[1;2;3;4;5];
Vj=[0;Rv;Rv;Rv;Rv];
Tj=[64.7;T2;T3;T4;T5];
Lj=[R;R;LRbF;LRbF;B];
fprintf('Stage Vj Tj Lj \n')
for i=1:5
fprintf('%d %.2f %.2f %.2f \n',j(i),Vj(i),Tj(i),Lj(i))
end
fprintf('\n CONSTANTS for ANTOINE EQUATION: \n')
fprintf('For Methanol \n')
mC1=16.5785;
mC2=3638.27;
mC3=239.500;

fprintf('mC1=%f \n mC2=%f \n mC3=%f',mC1,mC2,mC3)


fprintf('\n For Ethanol \n')
eC1=16.8958;
eC2=3795.17;
eC3=230.918;

fprintf('eC1=%f \n eC2=%f \n eC3=%f',eC1,eC2,eC3)


fprintf('\n For Propanol \n')
pC1=16.1154;
pC2=3483.67;
pC3=205.807;

fprintf('pC1=%f \n pC2=%f \n pC3=%f',pC1,pC2,pC3)


P=101.300;
fprintf('\nPressure taken in kPa: %.4f',P)
t=[T1;T2;T3;T4;T5];
fprintf('\nK-values and Vapour pressure for Methanol:')
for i=1:5
T=t(i);
mPs=exp(mC1-(mC2/(T+mC3)));
mK=mPs/P;
fprintf('\n mK%d=%f and mPs%d=%f \n',i,mK,i,mPs)
end
fprintf('\nK-values and Vapour pressure for Ethanol:')
for i=1:5
T=t(i);
ePs=exp(eC1-(eC2/(T+eC3)));
eK=ePs/P;
fprintf('\n eK%d=%f and ePs%d=%f \n',i,eK,i,ePs)
end
fprintf('\nK-values and Vapour pressure for Propanol:')
for i=1:5
T=t(i);
pPs=exp(pC1-(pC2/(T+pC3)));
pK=pPs/P;
fprintf('\n pK%d=%f and pPs%d=%f \n',i,pK,i,pPs)
end
%K-values and Vapour pressure values for METHANOL
mK1=1.000289;mPs1=101.329247;
mK2=1.248778;mPs2=126.501213;
mK3=1.546447;mPs3=156.655074;
mK4=1.900479;mPs4=192.518515;
mK5=2.318700;mPs5=234.884266;
mStg1=[-(D+R),Rv*mK2,0,0,0];
mStg2=[R,-(Rv*mK2+R),(Rv*mK3),0,0];
mStg3=[0,R,-((Rv*mK3)+LRbF),Rv*mK4,0];
mStg4=[0,0,LRbF,-((Rv*mK4)+3000),(Rv*mK5)];
mStg5=[0,0,0,3000,-((Rv*mK5)+B)];
fprintf('TRIDIAGONAL FORM FOR METHANOL')
mStage=[mStg1;mStg2;mStg3;mStg4;mStg5]
mD=[0;0;-(F*(m/100));0;0]
mA=[0,R,R,LRbF,3000];
mB=[-(D+R),-(Rv*mK2+R),-((Rv*mK3)+LRbF),-((Rv*mK4)+3000),-((Rv*mK5)+B)];
mC=[Rv*mK2,Rv*mK3,Rv*mK4,Rv*mK5,0];
mD=[0,0,-(F*(m/100)),0,0];
n=length(mD);
mz=zeros(1,n-1);
mz(1)=mC(1)/mB(1);

for k = 2:n-1
mz(k)=mC(k)/(mB(k)-mA(k)*mz(k-1));
end
mw(1)=mD(1)/mB(1);
for k=2:n
mw(k)=(mD(k)-mA(k)*mw(k-1))/(mB(k)-mA(k)*mz(k-1));

end
fprintf('\nmX%d=%.4f',5,mw(5))
mX(n)=mw(n);
for k=n-1:-1:1
mX(k)=mw(k)-mz(k)*mX(k+1);
fprintf('\nmX%d=%.4f',k,mX(k))
end

%K-values and Vapour pressure values for ETHANOL


eK1=0.570995;ePs1=57.841798;
eK2=0.729476;ePs2=73.895953;
eK3=0.923436;ePs3=93.544052;
eK4=1.158880;ePs4=117.394563;
eK5=1.442482;ePs5=146.123421;
eStg1=[-2600,2600*eK2,0,0,0];
eStg2=[2000,-(2600*eK2+2000),(2600*eK3),0,0];
eStg3=[0,2000,-((2600*eK3)+3000),2600*eK4,0];
eStg4=[0,0,3000,-((2600*eK4)+3000),(2600*eK5)];
eStg5=[0,0,0,3000,-((2600*eK5)+400)];
fprintf('\nTRIDIAGONAL FORM FOR ETHANOL')
eStage=[eStg1;eStg2;eStg3;eStg4;eStg5]
eD=[0;0;-(F*(e/100));0;0]
eA=[0,R,R,LRbF,3000];
eB=[-(D+R),-(Rv*eK2+R),-((Rv*eK3)+LRbF),-((Rv*eK4)+3000),-((Rv*eK5)+B)];
eC=[Rv*eK2,Rv*eK3,Rv*eK4,Rv*eK5,0];
eD=[0,0,-(F*(e/100)),0,0];
n=length(eD);
ez=zeros(1,n-1);
ez(1)=eC(1)/eB(1);

for k = 2:n-1
ez(k)=eC(k)/(eB(k)-eA(k)*ez(k-1));
end
ew(1)=eD(1)/eB(1);
for k=2:n
ew(k)=(eD(k)-eA(k)*ew(k-1))/(eB(k)-eA(k)*ez(k-1));

end
fprintf('\neX%d=%.4f',5,ew(5))
eX(n)=ew(n);
for k=n-1:-1:1
eX(k)=ew(k)-ez(k)*eX(k+1);
fprintf('\neX%d=%.4f',k,eX(k))
end

%K-values and Vapour pressure for PROPANOL


pK1=0.251331;pPs1=25.459809;
pK2=0.328592;pPs2=33.286416;
pK3=0.424935;pPs3=43.045897;
pK4=0.543907;pPs4=55.097753;
pK5=0.689487;pPs5=69.845058;
pStg1=[-2600,2600*pK2,0,0,0];
pStg2=[2000,-(2600*pK2+2000),(2600*pK3),0,0];
pStg3=[0,2000,-((2600*pK3)+3000),2600*pK4,0];
pStg4=[0,0,3000,-((2600*pK4)+3000),(2600*pK5)];
pStg5=[0,0,0,3000,-((2600*pK5)+400)];
fprintf('\nTRIDIAGONAL FORM FOR PROPANOL')
pStage=[pStg1;pStg2;pStg3;pStg4;pStg5]
pD=[0;0;-(F*(p/100));0;0]
pA=[0,R,R,LRbF,3000];
pB=[-(D+R),-(Rv*pK2+R),-((Rv*pK3)+LRbF),-((Rv*pK4)+3000),-((Rv*pK5)+B)];
pC=[Rv*pK2,Rv*pK3,Rv*pK4,Rv*pK5,0];
pD=[0,0,-(F*(p/100)),0,0];
n=length(pD);
pz=zeros(1,n-1);
pz(1)=pC(1)/pB(1);

for k = 2:n-1
pz(k)=pC(k)/(pB(k)-pA(k)*pz(k-1));
end
pw(1)=pD(1)/pB(1);
for k=2:n
pw(k)=(pD(k)-pA(k)*pw(k-1))/(pB(k)-pA(k)*pz(k-1));

end
fprintf('\npX%d=%.4f',5,pw(5))
pX(n)=pw(n);
for k=n-1:-1:1
pX(k)=pw(k)-pz(k)*pX(k+1);
fprintf('\npX%d=%.4f\n',k,pX(k))
end
fprintf('Normalization of mole fractions\n')
mx=zeros(1,5);
ex=zeros(1,5);
px=zeros(1,5);

mX5=0.1359; eX5=0.1893; pX5=0.4471;


mX4=0.2913; eX4=0.2619; pX4=0.3267;
mX3=0.4979; eX3=0.2883; pX3=0.2136;
mX2=0.7282; eX2=0.2839; pX2=0.1074;
mX1=0.9094; eX1=0.2071; pX1=0.0353;

mx=[mX1,mX2,mX3,mX4,mX5];
ex=[eX1,eX2,eX3,eX4,eX5];
px=[pX1,pX2,pX3,pX4,pX5];
for i=1:5
sum(i)=mx(i)+ex(i)+px(i);
fprintf('\nsum%d=%f\n',i,sum(i))
if sum(i)~=0
n_mx(i)=mx(i)/sum(i);
n_ex(i)=ex(i)/sum(i);
n_px(i)=px(i)/sum(i);
fprintf('Normalized values: ')
fprintf('n_mx%d=%.4f n_ex%d=%.4f
n_px%d=%.4f\n',i,n_mx(i),i,n_ex(i),i,n_px(i))
end
end

You might also like