Scribd
Upload
2 views2 pages

Flash Distillation Calculation Tool

The document contains code for calculating flash distillation parameters including vapor and liquid flow rates. Input variables like feed composition and operating conditions are requested. Equations are defined to calculate vapor to liquid ratio, component flow rates in vapor and liquid, and overall material balances. The code then displays output like realized vapor to liquid ratio, individual and total component flow rates in vapor and liquid phases.

Uploaded by

Shandy Adesya Rahman
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
2 views2 pages

Flash Distillation Calculation Tool

The document contains code for calculating flash distillation parameters including vapor and liquid flow rates. Input variables like feed composition and operating conditions are requested. Equations are defined to calculate vapor to liquid ratio, component flow rates in vapor and liquid, and overall material balances. The code then displays output like realized vapor to liquid ratio, individual and total component flow rates in vapor and liquid phases.

Uploaded by

Shandy Adesya Rahman
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

12/11/17 1:02 AM C:\Users\SHANDY\Do...\flashdestillation.

m 1 of 2

clear
clc
display ('Diketahui data:');
z1=input ('Fraksi Umpan Etana = ');
z2=input ('Fraksi Umpan Propana = ');
z3=input ('Fraksi Umpan Isobutana = ');
z4=input ('Fraksi Umpan n-Pentana = ');
F=input ('Umpan (kmol/jam) = ');
display (' ');
display ('Nilai K dari grafik de priester pada suhu 60C dan tekanan 14bar
(1400kPa):');
K1=input ('K1 (Etana) = ');
K2=input ('K2 (Propana) = ');
K3=input ('K3 (Isobutana) = ');
K4=input ('K4 (n-Pentana) = ');
atebak=input ('V/L tebak = ');

%Rumus

L = @(a) (F*z1/((a*K1)+1))+(F*z2/((a*K2)+1))+(F*z3/((a*K3)+1))+(F*z4/((a*K4)+1));
V = @(a) (F*z1/((1/(a*K1))+1))+(F*z2/((1/(a*K2))+1))+(F*z3/((1/(a*K3))+1))+(F*z4/((1/
(a*K4))+1));

x1 = @(a) (F*z1)/(L(a)*((a*K1)+1));
x2 = @(a) (F*z2)/(L(a)*((a*K2)+1));
x3 = @(a) (F*z3)/(L(a)*((a*K3)+1));
x4 = @(a) (F*z4)/(L(a)*((a*K4)+1));

O = @(a) (V(a)*K1*x1(a))+(L(a)*x1(a))+(V(a)*K2*x2(a))+(L(a)*x2(a))+(V(a)*K3*x3(a))+(L(a)
*x3(a))+(V(a)*K4*x4(a))+(L(a)*x4(a))-F;
a=fzero(O,atebak);

display(' ');
disp(['V/L real : ',num2str(a)]);

L = (F*z1/((a*K1)+1))+(F*z2/((a*K2)+1))+(F*z3/((a*K3)+1))+(F*z4/((a*K4)+1));
V = (F*z1/((1/(a*K1))+1))+(F*z2/((1/(a*K2))+1))+(F*z3/((1/(a*K3))+1))+(F*z4/((1/(a*K4))
+1));

display(' ');
disp('Laju alir uap dan cair');
disp(['Laju Alir Vapour/Uap (V) : ',num2str(V),'kmol/jam']);
disp(['Laju Alir Liquid/Cair (L) : ',num2str(L),'kmol/jam']);

%data tabel
etana=1;
propana=2;
isobutana=3;
npentana=4;
sigma=1234;
12/11/17 1:02 AM C:\Users\SHANDY\Do...\flashdestillation.m 2 of 2

sigmazi=z1+z2+z3+z4;

x1 = (F*z1)/(L*((a*K1)+1));
x2 = (F*z2)/(L*((a*K2)+1));
x3 = (F*z3)/(L*((a*K3)+1));
x4 = (F*z4)/(L*((a*K4)+1));
sigmaxi=x1+x2+x3+x4;

y1=K1*x1;
y2=K2*x2;
y3=K3*x3;
y4=K4*x4;
sigmayi=y1+y2+y3+y4;

L1=L*x1;
L2=L*x2;
L3=L*x3;
L4=L*x4;
sigmaL=L1+L2+L3+L4;

V1=V*y1;
V2=V*y2;
V3=V*y3;
V4=V*y4;
sigmaV=V1+V2+V3+V4;

disp([sprintf
('------------------------------------------------------------------------')]);
display ('Komponen zi xi yi L(Liquid) V(Vapour)'
)
disp([sprintf
('------------------------------------------------------------------------')]);
disp([sprintf('%4.0f %13.4f %12.4f %12.4f %12.4f %12.4f',[etana z1 x1 y1 L1 V1])]);
disp([sprintf('%4.0f %13.4f %12.4f %12.4f %12.4f %12.4f',[propana z2 x2 y2 L2 V2])]);
disp([sprintf('%4.0f %13.4f %12.4f %12.4f %12.4f %12.4f',[isobutana z3 x3 y3 L3 V3])]);
disp([sprintf('%4.0f %13.4f %12.4f %12.4f %12.4f %12.4f',[npentana z4 x4 y4 L4 V4])]);
disp([sprintf
('------------------------------------------------------------------------')]);
disp([sprintf('%4.0f %13.4f %12.4f %12.4f %12.4f %12.4f',[sigma sigmazi sigmaxi sigmayi
sigmaL sigmaV])]);
disp([sprintf
('------------------------------------------------------------------------')]);
display(' ');
display('Keterangan:');
display('1=etana');
display('2=propana');
display('3=isobutana');
display('4=n-pentana');
display('1234=TOTAL');

You might also like