Tarea 2. Absorción de SO2 en Agua
Tarea 2. Absorción de SO2 en Agua
Tarea 2. Absorción de SO2 en Agua
Resumen
En la presente tarea se calculó el número de platos necesarios para obtener una absorción del 90% de
SO2 contenido en una mezcla con aire, esto mediante el uso de agua pura. Para hacer dicho cálculo
se empleó el uso del programa OCTAVE y PRO II, los cuales proporcionaron un número de platos
teóricos de 3, posteriormente mediante el uso de la eficiencia total de platos se calculó el número de
plato reales, los cuales fueron de 12.
1. Descripción del problema
2. Metodología
Lo primero que se tiene que hacer es pasar los flujos a kmol/h, de la siguiente manera:
𝑘𝑔 𝐴𝑖𝑟𝑒 1 𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙 𝐴𝑖𝑟𝑒
150 ( ) = 5.18
ℎ 28.96 𝑘𝑔 ℎ
𝑘𝑔 𝐻2 𝑂 1 𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙 𝐻2 𝑂
6000 ( ) = 333.33
ℎ 18 𝑘𝑔 ℎ
Después mediante balances de materia se calculan los flujos de entrada y de salida de la fase
gaseosa y líquida.
Estos balances se pueden programar o hacerlos a mano, para esta tarea yo haré ambos casos.
𝑦1 = (1 − 𝑥𝑅 )𝑦𝑛+1 = 0.02
𝑘𝑚𝑜𝑙
𝑉1 = 𝑉𝑎𝑖𝑟𝑒 + 𝑉𝑛+1 𝑦1 = 5.309
ℎ
- Flujos y composiciones de la fase líquida
𝑘𝑚𝑜𝑙
𝐿𝑛 = 𝐿0 + 𝑉𝑆𝑂2 (𝑥𝑅 ) = 334.49
ℎ
𝑉𝑆𝑂2 (𝑥𝑅 )
𝑋𝑛 = = 0.0035
𝐿𝑛
Una vez que se conocen todos los flujos y composiciones, mediante el programa OCTAVE se calcula
el número de platos teóricos, esto por el método gráfico y analítico (Ec. De Kremser).
En el método gráfico obtenemos 3 platos teóricos y en el analítico 2.2162 que sería prácticamente 3
ya que no puede haber mitades o pedazos de platos.
Ahora para calcular los platos reales, utilizamos la siguiente ecuación:
𝑝𝑙𝑎𝑡𝑜𝑠 𝑡𝑒ó𝑟𝑖𝑐𝑜𝑠
𝐸=
𝑝𝑙𝑎𝑡𝑜𝑠 𝑟𝑒𝑎𝑙𝑒𝑠
Como E=0.25, los platos reales son 12.
Ahora simulamos el proceso en PRO II, para ver si el número de platos calculados cumplen las
especificaciones necesarias.
El primero paso en PRO II es seleccionar las unidades y componentes, en este caso dejamos las
unidades por defecto y seleccionamos los componentes:
3
1 2
En el paso 4 es importante poner primero al soluto, que en este caso es el SO2 y después el aire.
Ahora ponemos la columna Destillation, deshabilitamos la opción del condensador y rehervido,
ponemos el número de platos teóricos calculados previamente y nos aseguramos de que la convención
sea “Vapor and liquid to be on the feed tray”. Posteriormente colocamos las flechas de alimentación
y de salida, quedando finalmente de la siguiente manera:
Una vez especificados los flujos de entrada, ahora especificamos la columna, de la siguiente manera:
Figura 8. Perfil de presión de la columna
Una vez realizado todo lo anterior corremos la simulación y generamos una tabla de propiedades.
Si hacemos lo mismo, pero ahora con N=12 o sea con el número de platos reales, obtenemos lo
siguiente:
Figura 11. Tabla de propiedades con flujos y composiciones para N=12
3. Discusión de resultados
Podemos observar en las figuras 10 y 11 que sí se obtienen los flujos y composiciones
especificados y calculados con el programa en OCTAVE, es decir que los 3 platos teóricos sí son
suficientes para absorber el SO2 requerido.
También es importante mencionar que, aunque con los platos reales se obtiene una mayor
absorción, realmente la diferencia no es tan grande y por lo tanto es conveniente plantearse si
realmente vale la pena invertir mucho más dinero en una columna con mayor número de platos.
De mi parte yo no lo haría ya que con 3 platos se cumple sin problema las especificaciones que
pide el problema.
4. Conclusiones
La metodología empleada en el programa OCTAVE es certera, ya que comprobamos con ayuda
de PRO II que los datos obtenidos en él son prácticamente los mismos que se obtienen con el modelo
Destillation. Es decir que con el programa elaborado en OCTAVE podemos obtener datos
suficientemente confiables que cumplan con las especificaciones requeridas.
Código utilizado
% Absorción de gases
clear all
close all
clc
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% ABSORCIÓN DE SO2 EN AGUA %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Alimentación de gas
Vaire=5.18; %kmol/hr
Vso2=1.295; %kmol/hr
Vnm= Vaire+Vso2; %kmol/hr
ynm= 0.2;
% Salida de gas
y1=(1-xR)*ynm;
V1=Vaire+Vnm*y1;
%Linea de equilibrio
M=1000;
dx=(ynm-0)/(M-1);
xeq=x0:dx:xn;
figure(2)
plot(xeq,yeq,'r'),hold on
xlabel('x_A'), ylabel('y_A')
% Linea de operación
yop=Ln*xeq/Vnm+(V1*y1-L0*x0)/Vnm;
plot(xeq,yop), hold on
legend('línea de equilibrio','Línea de operación')
%Dibujo de platos
N=1000;
xi=zeros(1,N);
yi=zeros(1,N);
xp=x0;
for i=1:N
yi(i)=Ln/Vnm*xp+(V1*y1-L0*x0)/Vnm;
xi(i)=yi(i)/H;
xp=xi(i);
if xi(i)>=xn
break
endif
endfor
plot(xi(1:i),yi(1:i),'o'), hold on
if i<N
disp(['Número de platos= ', num2str(i)])
elseif i>=N
disp('Platos infinitos')
endif
yz=yi(1:i);
xz=(yz-(V1*y1-L0*x0)/Vnm)/(Ln/Vnm)
plot(xz,yz,'ok')
xu=zeros(2*i,1);
yu=zeros(2*i,1);
cont=0;
for k=1:2:2*i
cont=cont+1;
xu(k)=xz(cont);
xu(k+1)=xi(cont);
yu(k)=yz(cont);
yu(k+1)=yi(cont);
endfor
plot(xu,yu,'-k')
%Ecuación de Kremser
m=H;
A=L0/(m*V1);