UNIQUAC Model
% --------------------- %
% VLE Calculation %
% --------------------- %
comp = 2; %number of component
%Input data obtained from experiment
T = 323.15;
P =
[29.119;30.620;31.135;32.790;35.321;36.276;38.085;39.341;40.160;40.612;42.0
49;44.916;45.943;47.334];
xliq =
[0.2470;0.2842;0.294;0.3338;0.4028;0.4316;0.4872;0.5314;0.5513;0.5688;0.614
5;0.6989;0.7290;0.7730];
yvap =
[0.6710;0.7029;0.7109;0.7373;0.7772;0.7873;0.8100;0.8254;0.8326;0.8410;0.85
69;0.8890;0.8991;0.9160];
data = length(P);
for n=1:data
xliq(n,2) = 1-xliq(n,1);
yvap(n,2) = 1-yvap(n,1);
end
%Antoine constant
A = [7.25325;7.11572];
B = [1608.512;1684.123];
C = [-31.143;-43.568];
%experimental activity coefficient
for i=1:data
for j=1:comp
Psat(1,j) =10^(A(j,1)-(B(j,1)/(T+C(j,1))));
gamma_exp(i,j)=yvap(i,j)*P(i,1)/(xliq(i,j)*Psat(1,j));
end
end
%calculated activity coeeficient model using UNIQUAC model
% Volume parameters of pure components
%methanol
r_CH3 = 0.9011;
q_CH3 = 0.848;
r_OH = 1.0000;
q_OH = 1.200;
R1 = 1*r_CH3+1*r_OH;
Q1 = 1*q_CH3+1*q_OH;
%water
r_H2O = 0.9200;
q_H2O = 1.400;
R2 = 1*r_H2O;
Q2 = 1*q_H2O;
for i=1:data
sigma_r(i,1) = xliq(i,1)*R1+xliq(i,2)*R2;
phi(i,1) = xliq(i,1)*R1/sigma_r(i,1);
phi(i,2) = xliq(i,2)*R2/sigma_r(i,1);
sigma_q(i,1) = xliq(i,1)*Q1+xliq(i,2)*Q2;
theta(i,1) = xliq(i,1)*Q1/sigma_q(i,1);
theta(i,2) = xliq(i,2)*Q2/sigma_q(i,1);
sigma_q2(i,1) = xliq(i,1)*q1+xliq(i,2)*q2;
theta_2(i,1) = xliq(i,1)*q1/sigma_q2(i,1);
end
R = 1.987;