Proyecto Listado Txy Xy

2
% Calcular las temperaturas de ebullición de los componentes puros y determinar cuál es el compuesto más volátil clear all %Ecuación de Antoine para Benceno y Hexano Ab= 6.89745; Bb=1206.350; Cb=220.237; Ah= 6.87773; Bh=1171.530; Ch=224.366; P=input('Introduzca la presión a una atmósfera o a la presión del laboratorio de la en mm de Hg') Tbb=(Bb./(Ab-(log10(P))))-Cb; Tbh=(Bh./(Ah-(log10(P))))-Ch; %Generar un vector de temperaturas desde la temperatura de ebullición del más volátil hasta la del menos volátil. %Para cada temperatura del vector generado calcular la presión de vapor de cada componente, la fracción mol en la fase líquida (x) y la fracción mol en la fase vapor (y), éstas últimas del compuesto más volátil, según sea el caso condicion1=Tbb<Tbh condicion2=Tbb>Tbh if condicion1 T=Tbb:0.1:Tbh Pb= 10.^(Ab-(Bb./(T+Cb))); Ph= 10.^(Ah-(Bh./(T+Cb))); x=(P-Ph)./(Pb-Ph) y=(x.*Pb)./P elseif condicion2 T=Tbh:0.1:Tbb Pb=10.^(Ab-(Bb./(T+Cb))); Ph=10.^(Ah-(Bh./(T+Ch))); x=(P-Pb)./(Ph-Pb)

description

Matlab

Transcript of Proyecto Listado Txy Xy

% Calcular las temperaturas de ebullicin de los componentes puros y determinar cul es el compuesto ms voltilclear all%Ecuacin de Antoine para Benceno y Hexano Ab= 6.89745; Bb=1206.350; Cb=220.237;Ah= 6.87773; Bh=1171.530; Ch=224.366;P=input('Introduzca la presin a una atmsfera o a la presin del laboratorio de la en mm de Hg')Tbb=(Bb./(Ab-(log10(P))))-Cb;Tbh=(Bh./(Ah-(log10(P))))-Ch;%Generar un vector de temperaturas desde la temperatura de ebullicin del ms voltil hasta la del menos voltil.%Para cada temperatura del vector generado calcular la presin de vapor de cada componente, la fraccin mol en la fase lquida (x) y la fraccin mol en la fase vapor (y), stas ltimas del compuesto ms voltil, segn sea el casocondicion1=TbbTbhif condicion1T=Tbb:0.1:TbhPb= 10.^(Ab-(Bb./(T+Cb)));Ph= 10.^(Ah-(Bh./(T+Cb)));x=(P-Ph)./(Pb-Ph)y=(x.*Pb)./Pelseif condicion2T=Tbh:0.1:TbbPb=10.^(Ab-(Bb./(T+Cb)));Ph=10.^(Ah-(Bh./(T+Ch)));x=(P-Pb)./(Ph-Pb)y=(x.*Ph)./Pend%Graficar: Temperatura vs. x, y (diagrama de fases) plot(x,T,y,T)% Graficar y vs. x (curva de equilibrio lquido-vapor)plot(x,y)