DETERMINACIÓN DE NÚMERO DE PLATOS Y PLATO

DE ALIMENTACIÓN EN DESTILACIÓN CONTINUA

Morales Rodríguez, D.A. (20174159J); Timaná Lanfranco, J. L. (20152627K); Del Águila Rojas, P.
(20141262F); García Luyo, M. (20161153H); Castillo Trujillo, R. S. (20164133H); Gamarra Romero, J. E.
(20174153A)
Escuela de Ingeniería Petroquímica, FIP-UNI.
PI-136 / PQ-426 - Laboratorio de Operaciones Unitarias II, 8 de febrero del 2021
[email protected]; [email protected]; [email protected]; [email protected];
[email protected]; [email protected]

Resumen: En el presente laboratorio virtual se desarrolló un estudio para determinar el número de platos y el
plato de alimentación en un proceso de destilación continua de la mezcla cloroformo (1) / benceno (2). Se
obtuvieron los valores utilizando las curvas de equilibrio por medio de diferentes leyes termodinámicas mediante
el uso de códigos de programación, con el cual se modeló el comportamiento de la columna de destilación para
cada caso. Se determinó el plato de alimentación, el número de platos, así como sus propiedades de los flujos
del proceso. Finalmente se realizó una comparación con los resultados por los diferentes métodos aplicados.

1. INTRODUCCIÓN Vapor en el punto de rocío q=0

Vapor sobrecalentado q<0
Destilación
La destilación es un método para separar los Zona de enriquecimiento
componentes de una solución; depende de la Es la zona ubicada en la parte superior a la
distribución de las sustancias entre una fase gaseosa alimentación de la columna de destilación [1]. En
y una líquida, y se aplica a los casos en que todos esta sección tiene lugar el enriquecimiento del vapor
los componentes están presentes en las dos fases en componentes volátiles, por contacto con el
[1]. reflujo [4].
El balance de masa en esta zona queda definido
Método de McCabe – Thiele como:
Es un método gráfico basado en el método de 𝐺 =𝐿+𝐷 (5)
Lewis, mediante el cual se puede determinar el
número de platos o etapas teóricas necesarias para la 𝐺𝑦𝑛+1 = 𝐿𝑥𝑛 + 𝐷𝑥𝐷 (6)
separación de una mezcla binaria. Es adecuado en
aquellos casos donde los componentes de la mezcla De la cual la línea de operación de la zona de
binaria a separar son de tal naturaleza que sus enriquecimiento se define como:
mezclas tienen un comportamiento cercano a la 𝐿𝑛 𝐷
idealidad [2]. 𝑦𝑛+1 = 𝑥𝑛 + 𝑥𝐷 (7)
𝐺𝑛+1 𝐺𝑛+1
Condición de alimentación
Y tomando en cuenta que el reflujo externo se
Las condiciones de alimentación se pueden
define como:
representar por la fracción líquida, q, definida por el 𝐿𝑜
𝑅= (8)
cociente entre el calor necesario para convertir en
vapor saturado la alimentación y su calor latente de 𝐷
vaporización [3]. Se expresa por:
La línea de operación de la zona de enriquecimiento
ℎ𝑓𝑔 + 𝑐𝑝(𝑇𝑏 − 𝑇𝑓) en función del reflujo queda
𝐻𝑉 − 𝐻𝐹 (4) 𝑅 𝑥
𝑞= = 𝑦𝑛+1 = 𝑥𝑛 + 𝐷 (9)
𝐻𝑉 − 𝐻𝐿 ℎ𝑓𝑔 𝑅+1 𝑅+1

La fracción líquida puede adquirir los valores Zona de agotamiento

siguientes en función de las condiciones de la La sección del fraccionador abajo del punto de
alimentación: introducción de la mezcla de alimentación hace
referencia a la zona de agotamiento [1]. En esta
Liquido subenfriado q>1 zona
Líquidoelenagotamiento
ocurre el punto de burbuja q = 1 del
o empobrecimiento como:
Equilibrio Liquido-Vapor
líquido en componentes volátiles [4]. 0 < q<1 𝐿̅ = 𝐺̅ + 𝑊 (10)

El balance de masa en esta zona queda definido

 𝐿̅𝑥𝑚 = 𝐺̅ 𝑦𝑚+1 + 𝑊𝑥𝑊 (11) el conteo de platos para una eficiencia de Murphree
de 100% y se comparará con el indicado en la Tabla
De la cual la línea de operación de la sección de 1. Para la generación de lo indicado se usará el
agotamiento se define como: lenguaje de programación de Matlab para la
𝐿̅ 𝑊 creación del diagrama de McCabe Thiele y la
(12)
𝑦𝑚+1 = 𝑥𝑚 − 𝑥𝑊 obtención del número de platos correspondiente.
𝐺 ̅ 𝐺 ̅
La línea de operación en función del reflujo queda Tabla 1. Datos de obtenidos de experimentación para la mezcla
𝐿̅ 𝑊 (13) de cloroformo (1) y benceno (2)
𝑦𝑚+1 = 𝑥 − 𝑥
̅𝐿 − 𝑊̅ 𝑚 𝐿̅ − 𝑊̅ 𝑊 Valores
Condiciones
obtenidos
Eficacia de platos xF 0.44
El cálculo de los requerimientos de cada etapa
xD 0.92
teórica se ha realizado asumiendo que las fases se
encontraban en equilibrio [1]. Murphree define la xW 0.2
eficacia del plato en la fase líquida con la siguiente TF (°C) 22
ecuación. P (%) 28
𝑥 − 𝑥𝑛−1
𝐸𝑀𝐿 = ∗𝑛 (14) EML (%) 44
𝑥 −𝑦
𝑛 𝑛−1
3. RESULTADOS Y DISCUSION
OBJETIVOS 3.1 Resultados
 Determinar los flujos masicos de las En la Tabla 2 se muestran los valores obtenidos en
corrientes de salida en el tope y fondo de la la columna de destilación con las condiciones
columna por medio del balance de masa. mencionadas en la Tabla1:
 Determinar el estado de la alimentación por
Tabla 2 Valores obtenidos en la columna de destilación de para
la mezcla de cloroformo (1) y benceno (2)
medio de la línea q. Ley de
 Determinar las rectas de operación de Ley de Raoult
Valores Ley de
enriquecimiento y agotamiento de la Raoult modificada
obtenidos Raoult
columna. modificada volatilidad
 Determinar el número total de platos en la relativa
columna para una eficiencia de plato de Flujo de
100% y de 44%. 40.1925 40.1925 40.1925
tope (Kg/h)
2. MÉTODOS Y MATERIALES Flujo de
59.8075 59.8075 59.8075
En esta práctica se realizará el análisis de una fondo (Kg/h)
destilación continua para determinar el número de q 1.1941 1.2007 1.2007
platos para obtener una composición de destilado de
0.92 para la mezcla de cloroformo (1) / benceno (2), Reflujo
para lo cual se utilizará curvas de equilibrio interno 0.6734 0.6722 0.6893
obtenidas mediante la ley de Raoult, ley de Raoult mínimo
modificada y ley de Raoult modificada con
volatilidad relativa mediante el modelo de Margules Reflujo
simétrico. Para ello se trabajará con las condiciones externo 2.0622 2.051 2.2187
mostradas en la Tabla 1 y a una presión de 1 atm, mínimo
además se realizará Reflujo
0.862 0.8605 0.8823
interno
Reflujo
externo 6.2466 6.166 7.4978
Ley Numero
de RaoultdeModificada
platos 100% 9 9 9
EML
Plato de
alimentación 6° 5° 6°
100% EML
Numero de
platos 44% 19 19 19
EML
Plato de
alimentación 12° 11° 13°
44% EML Figura 3 Diagrama de McCabe Thiele para una eficiencia
de plato de 100% mediante la ley de Raoult Modificada
A continuación, se muestran las gráficas obtenidas
para la determinación de los platos en la columna
mediante con la construcción de la curva de
equilibrio mediante la ley de Raoult, ley de Raoult
modificada y la ley de Raoult modificada con
volatilidad relativa (los cálculos realizados se
encuentran en los apéndices):

Ley de Raoult

Figura 4 Diagrama de McCabe Thiele para una eficiencia de

plato de 44% mediante la ley de Raoult Modificada

Ley de Raoult Modificada alfa

Figura 1 Diagrama de McCabe Thiele para una eficiencia

de plato de 100% mediante la ley de Raoult

Figura 5 Diagrama de McCabe Thiele para una eficiencia

de plato de 100% mediante la ley de Raoult Modificada con
volatilidad relativa

Figura 2 Diagrama de McCabe Thiele para una eficiencia

de plato de 44% mediante la ley de Raoult
 con la alimentación en el plato 12°, 11° y
13° para las curvas realizadas mediante la
ley de Raoult, ley de Raoult modificada y
para la ley de Raoult modificada con alfa
respectivamente.

5. RECOMENDACIONES
 Se recomienda implementar un código que
utilice el método de Ponchon-Savarit en
lugar del método de McCabe-Thiele, para
poder comparar ambos métodos y medir las
Figura 6 Diagrama de McCabe Thiele para una eficiencia
de plato de 44% mediante la ley de Raoult Modificada con desviaciones.
volatilidad relativa  Se recomienda determinar el número de
platos para una mezcla que tenga un
Discusiones comportamiento más alejado de la idealidad
 El plato de alimentación varía según la Ley para evidenciar la importancia de la
Termodinámica utilizada, sin embargo, el determinación de la curva de equilibrio
número de platos totales es el mismo para mediante la ley de Raoult modificada.
cada ley.  Se recomienda calcular el calor especifico
 Se evidencia también una diferencia en los de la mezcla mediante el uso de una integral
valores obtenidos por cada Ley aplicada, que dependa de la temperatura.
esto se debe a que la Ley de Raoult
modificada y Ley de volatilidad con 6. BIBLIOGRAFÍA
volatilidad incluye más factores en su [1] D. Masa, R. E. Treybal, D. Francisco José
ecuación que permitirán obtener una mayor Lozano, and Bogota, “Operaciones De
cercanía a la realidad. Además, al tener un Transferencia.”
coeficiente de actividad cercano a 1, [2] L. Llanos, “Procesos de separación vapor,
muestra que la mezcla tiene un liquido-liquido ‘destilación.’”
comportamiento similar al de una solución [3] J. Ocon and G. Tojo, Problemas de
ideal. Ingeniería Química. 2006.
 Se obtuvo un valor del parámetro “q” [4] A. Marcilla Gomis, Introducción a las
mayor que 1 en todos los casos lo que operaciones de separación: Cálculo por
indica un líquido subenfriado en la etapas de equilibrio. 1998.
alimentación, a pesar de las pequeñas
diferencias por cada método.

4. CONCLUSIONES
 La alimentación en la columna se realiza en
condición de líquido subenfriado de la
mezcla cloroformo y benceno.
 Por el tope de la columna se obtuvo un flujo
de 40.1925 Kg/h, mientras que por el fondo
un flujo de 59.8075 Kg/h.
 Para las curvas realizadas mediante la ley
de Raoult, ley de Raoult modificada y para
la ley de Raoult modificada con alfa, se
obtuvieron un total de 9 platos teóricos para
una eficiencia de 100% por plato en la
columna, con la alimentación en el plato 6°,
5° y 6° respectivamente
 Por otro lado, para una eficiencia de 44%
por plato se obtuvo un total de 19 platos
reales
APÉNDICES
Ley de Raoult
Códigos usados: Algoritmo para el número de platos
%constantes de Antoine
Antoine=[3.96288 1106.904 218.552; 3.98523 yi=(xf+xd*q/R)/(1+q/R);
1184.240 217.572]; %cloroformo(1)/ xi=(-(q-1)*(1-R/(R+1))*xd-xf)/((q-
benceno(2)% 1)*R/(R+1)-q);
x=[0 0.125 0.25 0.375 0.5 0.625 figure(1);
0.75 0.875 1]; hold on;
%definir y equilibrio de acuerdo al modelo axis([0 1 0 1]);
utilizado%
yeq= [0 0.203567951 0.37507046 %Dibujando las curvas de operación y la
0.520790839 0.64557909 0.753211323 curva de equilibrio
0.846661367 0.92829306 1]; plot(x,yeq,'r');
%condiciones de operación% set(line([0 1],[0 1]),'Color',[0 1 0]);
P=1.01325; set(line([xd xi],[xd yi]),'Color',[1 0
Porcetanje=0.28; 1]);
EML=0.44; set(line([xf xi],[xf yi]),'Color',[1 0
F=100; %flujo en kg/h 1]);
xf=0.44; xd=0.92; xw=0.2; %composiciones set(line([xw xi],[xw yi]),'Color',[1 0
molares% 1]);
Tf=22; %°C%
%datos de los componentes% %Eficiencia EML=0.44
M1=119.38; M2=78.11; %pesos moleculares
del cloroformo(1)/ benceno(2) g/mol% %Zona de enriquecimiento
hfg1=29.24; hfg2=30.72; %entalpias de i=1;
vaporización a 25°C del cloroformo(1)/ xp(1)=xd;
benceno(2) KJ/mol% yp(1)=xd;
Cp1=0.8648; Cp2=1.72; %Capacidades while (xp(i)>xi),
caloríficas del cloroformo(1)/
benceno(2)KJ/kg.K% [Tr,xeq(i+1)]=TrRaoult(yp(i),P,Antoine);
%Probar eficiencia teórica xp(i+1)=xp(i)-(xp(i)-xeq(i+1))*EML;
%eficienciaL=[1 1 1 1 1 1 1 1 1]; yp(i+1)=R/(R+1)*xp(i+1)+xd/(R+1);
%Balance de masa% y=yp(i+1);
FF=F/(xf*M1+(1-xf)*M2)*1000; %mol/h% set(line([xp(i) xp(i+1)],[yp(i)
WW=FF*(xf-xd)/(xw-xd); %mol/h% yp(i)]),'Color',[0 0 1]);
DD=FF-WW; %mol/h% if (xp(i+1)>xi) set(line([xp(i+1)
W=WW*(xw*M1+(1-xw)*M2)/1000; xp(i+1)],[yp(i) yp(i+1)]),'Color',[0 0
D=DD*(xd*M1+(1-xd)*M2)/1000; 1]);
[Tb,y]=TbRaoult(Antoine,xf,P); end
Cp=(Cp1*xf*M1+Cp2*(1-xf)*M2)/1000; %Cp de i=i+1;
la mezcla en la alimentación Kj/mol.K) end
hfg=xf*hfg1+(1-xf)*hfg2;
q=(hfg+Cp*(Tb-Tf))/hfg; %Zona de agotamiento
xm=0; SS=(yi-xw)/(xi-xw);
error=1; yp(i)=SS*(xp(i)-xw)+xw;
while error>10^(-6) y=yp(i);
[T,y]=TbRaoult(Antoine,xm,P); set(line([xp(i) xp(i)],[yp(i-1)
ym=(q/(q-1))*xm-(1/(q-1))*xf; yp(i)]),'Color',[0 0 1]);
error=y-ym;
xm=xm+0.00001; while (xp(i)>xw),
end
[Tr,xeq(i+1)]=TrRaoult(yp(i),P,Antoine);
Rmin=(xd-xm)/(ym-xm)-1; %(L/D)min xp(i+1)=xp(i)-(xp(i)-xeq(i+1))*EML;
rmin=Rmin/(Rmin+1); %(L/V)min yp(i+1)=SS*(xp(i+1)-xw)+xw;
r=(Porcetanje+1)*rmin; % y=yp(i+1);
(L/V)operacion R=r/(1-r); % set(line([xp(i) xp(i+1)],[yp(i)
(L/D)operacion yp(i)]),'Color',[0 0 1]);
if (xp(i+1)>xw) set(line([xp(i+1)
%Intersección entre la linea de xp(i+1)],[yp(i) yp(i+1)]),'Color',[0 0
enriquecimiento y agotamiento 1]);
end
i=i+1;
end
hold off;

Algoritmo para presión de saturación

function [Pv]=Psat(T,Antoine) Algoritmo para la temperatura de burbuja
Pv=10.^(Antoine(:,1)- function [T,y]=TbRaoult(Antoine,x,P)
Antoine(:,2)./(T+Antoine(:,3))); Tv=Tsat(P,Antoine);
end n=length(x);
T=zeros(1,n);
y=zeros(1,n);
for i=1:n
function [T]=Tsat(P,Antoine)
z=x(i); T=Antoine(:,2)./(Antoine(:,1)-
error=1; log10(P))-Antoine(:,3);
emin=1E-4; end
T1=Tv(1,1)*z+Tv(2,1)*(1-z);
while error>emin
Pv=Psat(T1,Antoine);
Ley de Raoult Modificada
Pb=Pv(1,1)*z+Pv(2,1)*(1-z); Algoritmo para el número de platos
Pv1=Pv(1,1)*P/Pb; %mas liviano %constantes de Antoine
cloroformo(1)% Antoine=[3.96288 1106.904 218.552; 3.98523
Tb=Antoine(1,2)/(Antoine(1,1)- 1184.240 217.572]; %cloroformo(1)/
log10(Pv1))-Antoine(1,3); benceno(2)%
error=abs(T1-Tb); x=[0;0.125;0.25;0.375;0.5;0.625;0.75;0.875
T1=Tb; ;1];
end yeq= [0 0.178037028 0.351882387
T(i)=T1; 0.510702466 0.648698092 0.764486121
y(i)=x(i)*Pv(1,1)/P; 0.859592379 0.936977115 1];

%condiciones de operación%
end
P=1.01325; Porcetanje=0.28; EML=0.44;
T=T';
F=100; %flujo en kg/h
y=y';
xf=0.44; xd=0.92; xw=0.2; %composiciones
end
molares%
Tf=22; %°C%
Algoritmo para la temperatura de rocio %datos de los componentes%
function[Tr,x1]=TrRaoult(z1,P,Antoine) M1=119.38; M2=78.11; %pesos moleculares
Tsatu=Tsat(P,Antoine); del cloroformo(1)/ benceno(2) g/mol%
n=length(z1); hfg1=29.24; hfg2=30.72; %entalpias de
Tr=zeros(n,1); vaporización a 25°C del cloroformo(1)/
x1=zeros(n,1); benceno(2) KJ/mol%
Pva=zeros(n,2); Cp1=0.8648; Cp2=1.72; %Capacidades
for i=1:n caloríficas del cloroformo(1)/
error=1; benceno(2)KJ/kg.K%
emin=1E-4; %Probar eficiencia teórica
T00=Tsatu(1,1)*z1+Tsatu(2,1)*(1-z1); %eficienciaL=[1 1 1 1 1 1 1 1 1];
while error>emin %Balance de masa%
Pva(i,1)=10^(Antoine(1,1)- FF=F/(xf*M1+(1-xf)*M2)*1000; %mol/h%
Antoine(1,2)/(T00(i)'+Antoine(1,3))); WW=FF*(xf-xd)/(xw-xd); %mol/h%
Pva(i,2)=10^(Antoine(2,1)- DD=FF-WW; %mol/h%
Antoine(2,2)/(T00(i)'+Antoine(2,3))); W=WW*(xw*M1+(1-xw)*M2)/1000;
Pr=1/(z1(i)/(Pva(i,1))+(1- D=DD*(xd*M1+(1-xd)*M2)/1000;
z1(i))/(Pva(i,2))); [Tb,y]=TbRaoultM(Antoine,xf,P);
P2sat=Pva(i,2)*P./Pr; %mas Cp=(Cp1*xf*M1+Cp2*(1-xf)*M2)/1000; %Cp de
pesado:benceno(2) % la mezcla en la alimentación Kj/mol.K)
Tr(i)=Tsat(P2sat,Antoine(2,:)); hfg=xf*hfg1+(1-xf)*hfg2;
error=abs(Tr(i)'-T00(i)); q=(hfg+Cp*(Tb-Tf))/hfg;
T00(i)=Tr(i); xm=0;
end error=1;
x1(i)=z1(i)*P/(Pva(i,1)); while error>10^(-6)
end [T,y]=TbRaoultM(Antoine,xm,P);
end ym=(q/(q-1))*xm-(1/(q-1))*xf;
error=y-ym;
xm=xm+0.00001;
Algoritmo para la temperatura de saturación end

Rmin=(xd-xm)/(ym-xm)-1; %(L/D)min
rmin=Rmin/(Rmin+1); %(L/V)min
r=(Porcetanje+1)*rmin; %(L/V)operacion
R=r/(1-r); %(L/D)operacion

%Intersección entre la linea de

enriquecimiento y agotamiento

yi=(xf+xd*q/R)/(1+q/R);
xi=(-(q-1)*(1-R/(R+1))*xd-xf)/((q-
1)*R/(R+1)-q); %Dibujando las curvas de operación y la
figure(1); curva de equilibrio
hold on; plot(x,yeq,'r');
axis([0 1 0 1]); set(line([0 1],[0 1]),'Color',[0 1 0]);
set(line([xd xi],[xd yi]),'Color',[1 0 Pv=10.^(Antoine(:,1)-
1]); Antoine(:,2)./(T+Antoine(:,3)));
set(line([xf xi],[xf yi]),'Color',[1 0 end
1]);
set(line([xw xi],[xw yi]),'Color',[1 0
1]); Algoritmo para la temperatura de burbuja
function [T,y]=TbRaoultM(Antoine,x,P)
%Eficiencia EML=0.44 Tv=Tsat(P,Antoine);
n=length(x);
%Zona de enriquecimiento T=zeros(1,n);
i=1; y=zeros(1,n);
xp(1)=xd; A12=-0.236;
yp(1)=xd; A21=-0.167;
while (xp(i)>xi), for i=1:n
z=x(i);
[Tr,xeq(i+1)]=TrRaoultM(yp(i),P,Antoine); gamma1=exp(((1-x(i))^2)*(A12+2*(A21-
xp(i+1)=xp(i)-(xp(i)-xeq(i+1))*EML; A12)*x(i)));
yp(i+1)=R/(R+1)*xp(i+1)+xd/(R+1); gamma2=exp(((x(i))^2)*(A21+2*(A12-A21)*(1-
y=yp(i+1); x(i))));
set(line([xp(i) xp(i+1)],[yp(i) error=1;
yp(i)]),'Color',[0 0 1]); emin=1E-4;
if (xp(i+1)>xi) set(line([xp(i+1) T1=Tv(1,1)*z+Tv(2,1)*(1-z);
xp(i+1)],[yp(i) yp(i+1)]),'Color',[0 0 while error>emin
1]); Pv=Psat(T1,Antoine);
end Pb=gamma1*Pv(1,1)*z+gamma2*Pv(2,1)*(1-z);
i=i+1; Pv1=Pv(1,1)*P/Pb;
end Tb=Antoine(1,2)/(Antoine(1,1)-log10(Pv1))-
Antoine(1,3);
%Zona de error=abs(T1-Tb);
agotamiento SS=(yi- T1=Tb;
xw)/(xi-xw); end
yp(i)=SS*(xp(i)-xw)+xw; T(i)=T1;
y=yp(i); y(i)=x(i)*gamma1*Pv(1,1)/P;
set(line([xp(i) xp(i)],[yp(i-1) end
yp(i)]),'Color',[0 0 1]); T=T';
y=y';
while (xp(i)>xw), end

[Tr,xeq(i+1)]=TrRaoultM(yp(i),P,Antoine);
xp(i+1)=xp(i)-(xp(i)-xeq(i+1))*EML; Algoritmo para la temperatura de rocio
function[Tr,x1]=TrRaoultM(z1,P,Antoine)
yp(i+1)=SS*(xp(i+1)-xw)+xw;
Tsatu=Tsat(P,Antoine);
y=yp(i+1);
n=length(z1);
set(line([xp(i) xp(i+1)],[yp(i)
Tr=zeros(n,1);
yp(i)]),'Color',[0 0 1]);
x1=zeros(n,1);
if (xp(i+1)>xw) set(line([xp(i+1)
x2=zeros(n,1);
xp(i+1)],[yp(i) yp(i+1)]),'Color',[0 0
Pisatu=zeros(n,2);
1]);
for i=1:n
end
gamma1=1;
i=i+1;
gamma2=1;
end
error=1;
hold off;
emin=10^-3;
T00=Tsatu(1,1)*z1+Tsatu(2,1)*(1-z1);
Algoritmo para presión de saturación while error>emin
function [Pv]=Psat(T,Antoine) Pisatu(i,1)=10^(Antoine(1,1)-
Antoine(1,2)/(T00(i)'+Antoine(1,3)));
Pisatu(i,2)=10^(Antoine(2,1)-
Antoine(2,2)/(T00(i)'+Antoine(2,3)));
x1(i)=z1(i)*P./(gamma1*Pisatu(i,1));
x2(i)=(1-
z1(i))*P./(gamma2*Pisatu(i,2));
A12=-0.236;
A21=-0.167;
gamma1=exp(((1-
x1(i)).^2).*(A12+2*(A21-A12)*x1(i)));
gamma2=exp(((1-x2(i)).^2)*(A21+2*(A12- error=abs(Tr(i)'-T00(i));
A21).*x2(i))); T00(i)=Tr(i);
Pr=1/(z1(i)./(gamma1*Pisatu(i,1))+(1- end x1(i)=z1(i)*P/(gamma1*Pisatu(i,1));
z1(i))./(gamma2*Pisatu(i,2))); end
P2sat=Pisatu(i,2)*P./Pr; end
Tr(i)=Tsat(P2sat,Antoine(2,:));
 R=r/(1-r); %(L/D)operacion
Algoritmo para la temperatura de saturación
function [T]=Tsat(P,Antoine) %Intersección entre la linea de
T=Antoine(:,2)./(Antoine(:,1)-log10(P))- enriquecimiento y agotamiento
Antoine(:,3);
end yi=(xf+xd*q/R)/(1+q/R);
xi=(-(q-1)*(1-R/(R+1))*xd-xf)/((q-
1)*R/(R+1)-q);
Ley de Raoult Modificada con alfa figure(1);
Algoritmo para el número de platos hold on;
%constantes de Antoine axis([0 1 0 1]);
Antoine=[3.96288 1106.904 218.552; 3.98523
1184.240 217.572]; %cloroformo(1)/ %Dibujando las curvas de operación y la
benceno(2)% curva de equilibrio
x=[0:0.125:1]; plot(x,yeq,'r');
alf=1.76596; set(line([0 1],[0 1]),'Color',[0 1 0]);
yeq=alf.*x./(1+(alf-1).*x); set(line([xd xi],[xd yi]),'Color',[1 0
%condiciones de operación% 1]);
P=1.01325; Porcetanje=0.28; EML= 0.44; set(line([xf xi],[xf yi]),'Color',[1 0
F=100; %flujo en kg/h 1]);
xf=0.44; xd=0.92; xw=0.2; %composiciones set(line([xw xi],[xw yi]),'Color',[1 0
molares% 1]);
Tf=22; %°C%
%datos de los componentes% %Eficiencia EML=0.44
M1=119.38; M2=78.11; %pesos moleculares
del cloroformo(1)/ benceno(2) g/mol% %Zona de enriquecimiento
hfg1=29.24; hfg2=30.72; %entalpias de i=1;
vaporización a 25°C del cloroformo(1)/ xp(1)=xd;
benceno(2) KJ/mol% yp(1)=xd;
Cp1=0.8648; Cp2=1.72; %Capacidades while (xp(i)>xi),
caloríficas del cloroformo(1)/ xp(i+1)=xp(i)-(xp(i)-yp(i)/(alf-(alf-
benceno(2)KJ/kg.K% 1)*yp(i)))*EML;
%Probar eficiencia teórica yp(i+1)=R/(R+1)*xp(i+1)+xd/(R+1);
%eficienciaL=[1 1 1 1 1 1 1 1 1]; y=yp(i+1);
%Balance de masa% set(line([xp(i) xp(i+1)],[yp(i)
FF=F/(xf*M1+(1-xf)*M2)*1000; %mol/h yp(i)]),'Color',[0 0 1]);
% WW=FF*(xf-xd)/(xw-xd); %mol/h% if (xp(i+1)>xi) set(line([xp(i+1)
DD=FF-WW; %mol/h% xp(i+1)],[yp(i) yp(i+1)]),'Color',[0 0
W=WW*(xw*M1+(1-xw)*M2)/1000; 1]);
D=DD*(xd*M1+(1-xd)*M2)/1000; end
[Tb,y]=TbRaoultM(Antoine,xf,P); i=i+1;
Cp=(Cp1*xf*M1+Cp2*(1-xf)*M2)/1000; %Cp end
de la mezcla en la alimentación
Kj/mol.K) hfg=xf*hfg1+(1-xf)*hfg2; %Zona de agotamiento
q=(hfg+Cp*(Tb-Tf))/hfg; SS=(yi-xw)/(xi-xw);
xm=(-((xf+q)*(1-alf)+alf)+sqrt(((xf+q)*(1- yp(i)=SS*(xp(i)-xw)+xw;
alf)+alf)^2-4*(q*(alf-1))*-xf))/(2*q*(alf- y=yp(i);
1)); set(line([xp(i) xp(i)],[yp(i-1)
ym=(xm*q-xf)/(q-1); yp(i)]),'Color',[0 0 1]);
Rmin=(xd-xm)/(ym-xm)-1; %(L/D)min
rmin=Rmin/(Rmin+1); %(L/V)min while (xp(i)>xw),
r=(Porcetanje+1)*rmin; % xp(i+1)=xp(i)-(xp(i)-yp(i)/(alf-(alf-
(L/V)operacion 1)*yp(i)))*EML;
yp(i+1)=SS*(xp(i+1)-xw)+xw;
y=yp(i+1);
set(line([xp(i) xp(i+1)],[yp(i)
yp(i)]),'Color',[0 0 1]);
if (xp(i+1)>xw) set(line([xp(i+1)
xp(i+1)],[yp(i) yp(i+1)]),'Color',[0 0
1]);
end
i=i+1;
end
hold off;
Algoritmo para presión de saturación
function [Pv]=Psat(T,Antoine)
Pv=10.^(Antoine(:,1)-
Antoine(:,2)./(T+Antoine(:,3)));
end

Algoritmo para la temperatura de burbuja

function [T,y]=TbRaoultM(Antoine,x,P)
Tv=Tsat(P,Antoine);
n=length(x);
T=zeros(1,n);
y=zeros(1,n);
A12=-0.236;
A21=-0.167;
for i=1:n
z=x(i);
gamma1=exp(((1-x(i))^2)*(A12+2*(A21-
A12)*x(i)));
gamma2=exp(((x(i))^2)*(A21+2*(A12-A21)*(1-
x(i))));
error=1;
emin=1E-4;
T1=Tv(1,1)*z+Tv(2,1)*(1-z);
while error>emin
Pv=Psat(T1,Antoine);
Pb=gamma1*Pv(1,1)*z+gamma2*Pv(2,1)*(1-z);
Pv1=Pv(1,1)*P/Pb;
Tb=Antoine(1,2)/(Antoine(1,1)-log10(Pv1))-
Antoine(1,3);
error=abs(T1-Tb);
T1=Tb;
end
T(i)=T1;
y(i)=x(i)*gamma1*Pv(1,1)/P;
end
T=T';
y=y';
end

Algoritmo para la temperatura de saturación

function [T]=Tsat(P,Antoine)
T=Antoine(:,2)./(Antoine(:,1)-log10(P))-
Antoine(:,3);
end
Cálculos:  Luego hallamos h̅ fg
a) Balance de masa
h̅ fg = h̅ fg1 ∗ zF + h̅ fg2 ∗ (1 − zF )
Datos: h̅ fg = 29.24 ∗ 0.44 + 30.72 ∗
0.56
F = 100 kg/h zF = 0.44 h̅ fg = 30.0688 kJ/mol
xD = 0.92 xW = 0.20
 Por último, hallamos la Tb con
Cloroformo (1) → M̅ 1 = 119.38 g/mol
Benceno (2) → M̅ 2 = 78.11 g/mol LRM A la presión de 1.01325 bar

T1sat = 61.17ºC T2sat = 80.01ºC


̅ Primero hallamos la masa molar promedio
Tb = Tsat ∗ zF + Tsat ∗ (1 − zF)
Mprom = 0.44 ∗ 119.38 + 0.56 ∗ 78.11 1 2

M̅ prom = 96.2688 g/mol Tb = 71.72ºC

𝐴12 = −0.236
 Luego hallamos el valor de F en unidades de 𝐴21 = −0.167
 mol/h
𝛿1 = 𝑒𝑥𝑝((1 − 𝑥𝑖)2 ∗ (𝐴12 + 2 ∗ (𝐴21 − 𝐴12) ∗ 𝑥𝑖));
100 = 1.0386 kmol/h 𝛿1 = 0.9465
F= 𝛿2 = 𝑒𝑥𝑝((𝑥𝑖)2 ∗ (𝐴21 + 2 ∗ (𝐴12 − 𝐴21) ∗ (1 − 𝑥𝑖)));
96.2688
𝛿2 = 0.9538
 Aplicando el balance de materia, tenemos:
Psat = 1.411 bar Psat = 0.779 bar
1 2
sat sat
F= W+D Pb = P1 ∗ zF ∗ 𝛿1 + P 2 ∗ (1 − zF) ∗ 𝛿2
1.0386 = W + D Pb = 1.00371 bar

zF ∗ F = xW ∗ W + xD ∗ D P
0.45705 = 0.20 ∗ W + 0.92 ∗ D P1sat = P1sat ∗ = 1.4244 bar
Pb

𝐖 = 𝟓𝟗. 𝟖𝟎𝟕𝟓 Tb = 72.0316ºC

𝐤𝐠/𝐡
𝐃 = 𝟒𝟎. 𝟏𝟗𝟐𝟓 𝐤𝐠/𝐡 Luego se repite el proceso (iteraciones) hasta que los
valores no cambien, hallando:
b) Cálculo de q
Tb = 72.0186ºC
h̅𝐹g +Cp ∗(Tb −T𝐹 ) (1)
q= h̅𝐹
g
 Luego de hallar estos valores se reemplaza
en (1)
Datos:

h̅ fg1 = 29.24 kJ/mol h̅ fg2 = 30.72 30.0688 + 0.1207 ∗ (72.0186 − 22)

q=
kJ/mol Cp1 = 0.8648 kJ/kgK Cp2 = 1.72 30.0688
kJ/kgK
𝐪 = 𝟏. 𝟐𝟎𝟎𝟕
 Primero hallamos el valor de Cp
 Pendiente de la recta q
Cp = Cp1 ∗ M̅ 1 ∗ zF + Cp2 ∗ M̅ 2 ∗ (1 − zF )
q
Cp = 0.8648 ∗ 119.38 ∗ 0.44 + 1.72 ∗ 78.11 m = q − 1 = 5.9826
∗ 0.56
Cp = 0,1207 kJ/molK c) Cálculo del Reflujo interno y externo
mínimo
 Hallamos las coordenadas del punto que 𝐋
pasa por la recta q y la curva de equilibrio, 𝐑𝐞𝐱𝐭𝐞𝐫𝐧𝐨𝐦𝐢𝐧 = ( )𝐦𝐢𝐧 = 𝟐. 𝟐𝟏𝟖𝟕
𝐃
en donde la ecuación de la curva de
equilibrio tiene la siguiente forma: Para el reflujo interno mínimo, se calcula así:
Rmin
α∗x = 0.6893
y= Rmin + 1
1 + (α − 1) ∗ x
𝐋
donde α = 1.7659 para LRM, entonces: 𝐑𝐢𝐧𝐭𝐞𝐫𝐧𝐨𝐦𝐢𝐧 = ( )𝐦𝐢𝐧 = 𝟎. 𝟔𝟖𝟗𝟑
𝐕
1.7659x
y= d) Ecuaciones de las rectas de operación
1 + 0.7659 ∗ x (rectificación y agotamiento)
Hallando la pendiente de 2 puntos, e igualando con Si Rinterno = 1.28 ∗ Rinternomin
la pendiente de la recta q, obtenemos:
Rinterno = 1.28 ∗ 0.6893
1.7659X ) y (0,44; 0.44) 𝐑𝐢𝐧𝐭𝐞𝐫𝐧𝐨 = 𝟎. 𝟖𝟖𝟐𝟑
(x; y = 1+0.7659∗X
 1.7659x
Se sabe que:
0.44 − R
m = 5.9826 = 1 + 0.7659 ∗ x
= 0.8823
0.44 − x R+
1

x = 0.4682 𝐑𝐞𝐱𝐭𝐞𝐫𝐧𝐨 = 𝐑 = 𝟕. 𝟒𝟗𝟕𝟖

y = 0.6085
 Recta de enriquecimiento
 Ecuación de la recta q
R xD
y=( )∗𝑥+
y − 0.6085 R+ R+1
m = 5.9826 =
x − 0.4682 1
𝒚 = 𝟎. 𝟖𝟖𝟐𝟑𝒙 + 𝟎. 𝟏𝟎𝟖𝟑
y = 5.9826x − 2.1925
 Recta de agotamiento
 Hallamos la pendiente de la recta de
operación (Rmin) Teniendo la recta de enriquecimiento y la recta q,
vamos a igualarlas y hallar el valor de x e y
xD − 0.6085 0.92 − 0.6085
mRmin = = y = 0.8823𝑥 + 0.1083
xD − 0.4682 0.92 − 0.4682
y = 5.9826x − 2.1925
mRmin = 0.6895
XD
Obtenemos:
Se sabe que la recta pasa por el punto (0; )
Rmin+1
𝑥 = 0.4511 𝑦 = 0.5063
Hallando la pendiente de 2 puntos, e igualando con
Hallando la pendiente que pasa por el punto (xW, xW)
la pendiente de la recta de operación, obtenemos:
xD 0.5063 − 0.20
xD − 𝑚 = 0 − 5411 − 0,20 = 1.2199
R min + 1
m = 0.6895 =
xD − 0 Reemplazando en la ecuación y = mx + b, obtenemos
b:
Rmin = 2.2187
 𝟎. 𝟓𝟎𝟔𝟑 = 𝟏. 𝟐𝟏𝟗𝟗 ∗ 𝟎. 𝟒𝟓𝟏𝟏 + 𝐛
𝐛 = −𝟎. 𝟎𝟒𝟒𝟎

𝐲 = 𝟏. 𝟐𝟏𝟗𝟗𝒙 − 𝟎. 𝟎𝟒𝟒𝟎

e) Numero de platos y la posición del plato Tabla 3

de alimentación a una eficiencia de 100% Cálculo del número de etapas a 100%
Para hallar la cantidad de numero de platos se Etapas y x
tendrá que hallar el valor de “x” a través de la 1 0.9200 0.8669
siguiente ecuación: 2 0.8732 0.7958
y 3 0.8105 0.7077
x=
1.7659 − 0.7659 ∗ y 4 0.7327 0.6082
5 0.6449 0.5070
Luego de hallar x, hallaremos el valor de “y”, con la 6 0.5557 0.4146 < zF
ecuación de la zona de enriquecimiento hasta
7 0.4617 0.3270
encontrar el punto del platón de alimentación, el
cual luego seguiremos hallando el valor de “y”, pero
8 0.3549 0.2375
ahora con la ecuación de la zona de agotamiento. 9 0.2457 0.1558 < xW
Nota. Fuente: Elaboración propia
Ejemplo:
Se concluye que habrá 9 etapas (platos) y el plato de
Como nuestro primer valor es xD = 0.92, este será alimentación se encuentra en el plato 6.
nuestro valor de y
f) Numero de platos y la posición del plato
Es, decir y = 0.92 de alimentación a una eficiencia de 44%

0.92 Para este punto, el cálculo de xeq e y será igual al

x= punto e, pero para el cálculo de x será de la
1.7659 − 0.7659 ∗ 0.92
siguiente forma:
x = 0.8869
Ejemplo:
Luego hallamos y con la ecuación de la zona de
enriquecimiento. Tenemos como primer valor xD = 0.92, que será
nuestro valor inicial de x1 e “y”, es decir, usaremos
y = 0.8823x + 0.1083 la siguiente ecuación:
y = 0.8823 ∗ 0.8869 + 0.1083 𝑦𝑖
y = 0.8732 𝑥 =
𝑒𝑞(𝑖+1)
1.7659 − 0.7659 ∗ 𝑦𝑖
Y así sucesivamente hasta encontrar el plato de
alimentación. Para 𝑦1 = 0.92 e x1 = 0.92

0.92
Una vez encontrado el plato de alimentación la 𝑥𝑒𝑞(2 =
ecuación para hallar los valores de y hasta que )
1.7659 − 0.7659 ∗ 0.92
nuestro último valor de x sea menor a xW, será:
𝑥𝑒𝑞(2) = 0.8869
y = 1.2199x − 0.0440
y = 1.2199 ∗ 0.4146 − 0.0440 Ahora, con este valor, hallamos x con un EML = 44%
y = 0.4617
𝑥𝑖+1 = 𝑥𝑖 − (𝑥𝑖 − 𝑥𝑒𝑞(𝑖+1)) ∗ 𝐸𝑀𝐿
Para los siguientes valores se adjuntará la tabla 3
Para i = 1

𝑥2 = 𝑥1 − (𝑥1 − 𝑥𝑒𝑞(𝑖+1)) ∗ 𝐸𝑀𝐿 Luego hallado el valor de 𝑥2, con este valor hallo el
valor de 𝑦2, con la ecuación de la zona de
𝑥2 = 0.92 − (0.92 − 0.8869) ∗ 0.44
enriquecimiento como en el punto e.
𝑥2 = 0.8966
 Se concluye que habrá 19 etapas (platos) y el plato
𝑦2 = 0.8823𝑥2 + 0.1083 de alimentación se encuentra en el plato 13.
𝑦2 = 0.8823 ∗ 0.8966 + 0.1083
𝑦2 = 0.8994

Y así sucesivamente hasta encontrar el plato de

alimentación.

Una vez encontrado el plato de alimentación la

ecuación para hallar los valores de y hasta que
nuestro último valor de x sea menor a xW, será:

y = 1.2199x − 0.0440
y = 1.2199 ∗ 0.4271 − 0.0440
y = 0.4770

Tabla 4
Cálculo del número de etapas a 44%
Etapas y xeq x
1 0.9200 0.8669 0.8966
2 0.8994 0.8351 0.8695
3 0.8755 0.7993 0.8386
4 0.8482 0.7599 0.8040
5 0.8176 0.7174 0.7659
6 0.7841 0.6728 0.7249
7 0.7479 0.6269 0.6818
8 0.7098 0.5808 0.6373
9 0.6706 0.5355 0.5925
10 0.6311 0.4921 0.5483
11 0.5921 0.4511 0.5056
12 0.5544 0.4133 0.4650
13 0.5185 0.3788 0.4271 < zF
14 0.4770 0.3406 0.3890
15 0.4306 0.2998 0.3498
16 0.3827 0.2598 0.3102
17 0.3344 0.2215 0.2712
18 0.2868 0.1855 0.2335
19 0.2408 0.1523 0.1977 < xW
Nota. Fuente: Elaboración propia