Diagramas T-xy mediante el calculo de punto de burbuja usando ley de Raoult con MatLab

Se muestra como generar un diagrama de equilibrio liquido - vapor para mezclas binarias usando la ley de Raoult en Matlab.

% Este programa genera el diagrama T - xy usando el calculo de la temperatura
% del punto de burbuja empleando la Ley de Raoult para predecir el equilibrio

clear
clc

global P x
global A B C
global y

% Presion y composicion liquida
P = 2; % atm
x(1) = 0.0;
x(2) = 1 - x(1);

% Numero de datos calculados
n = 20;

% Parametros de la ecuacion de Antoine
A = [6.87987, 6.95087];
B = [1196.76, 1342.31];
C = [219.161, 219.187];

% Temperatura de ebullicion de los componentes
Tb = [80.1, 110.6];

for i = 1:n+1
% Temperatua supuesta
T = Tb(2);

% Soluci?n del Punto de Burbuja
T = fsolve('LeyRaoult',T)

% Asignaci?n de vectores para graficar
x_(i) = x(1);
y_(i) = y(1);
T_(i) = T;

% Reasignaci?n de nuevas xi
x(1) = x(1) + 1.0/n;
x(2) = 1.0 - x(1);
clc

end

plot(x_,T_,y_,T_)

xlabel('x - y')
ylabel('Temperatura, C')
axis square
axis normal
axis ([0 1 0.98*min(T_) 1.02*max(T_)])
legend('Linea de punto de burbuja, T-x','Linea de punto de rocio, T-y')

function Sy = LeyRaoult(T)

global P x
global A B C
global y

for i = 1:2
% Calculo de Ps (mmHg) con Ecc. de Antoine
Ps(i) = 10^(A(i)-B(i)/(T + C(i)));
% Calculo de K
k(i) = Ps(i)/(P*760);
% Calculo de la fraccion molar de vapor
y(i) = k(i)*x(i);
end

Sy = 1.0 - sum(y);

end