TORRES DE ABSORCIÓN POR ETAPAS - 101 -

44.
. TTORRES
ORRES DDE
E AABSORCIÓN
BSORCIÓN PPOR
OR EETAPAS
TAPAS

En esta operación dos fases: gas y líquido, se ponen en contacto y uno o

varios componentes de la fase gaseosa (soluto) son absorbidos por la fase
líquida. Se recurre a la absorción cuando: a) es necesario recuperar uno o
más componentes de una mezcla gaseosa y b) es necesario eliminar uno o más
componentes de una mezcla gaseosa.
El primer caso se presenta en la fabricación del amoniaco, donde del re-
actor sale una mezcla gaseosa compuesta por amoniaco y aire. El amoniaco se
recupera de esta mezcla poniéndola en contacto con agua pura en una torre de
absorción, como resultado de esta operación el amoniaco es absorbido por el
agua y el aire queda libre de amoniaco. Posteriormente la mezcla amoniaco –
agua es sometida a un proceso de destilación para obtener amoniaco relativa-
mente puro.
El segundo caso se presenta en la fabricación de ácido sulfúrico, donde
como producto intermedio se obtiene una mezcla gaseosa de SO2, CO2 y aire.
El SO2 y el aire deben pasar luego a los reactores catalíticos donde el SO2
se convierte en SO3, sin embargo se debe evitar el paso del CO2 pues este
compuesto ensucia los catalizadores. Con este fin se ponen en contacto la
mezcla gaseosa con una mezcla alcalina de agua e hidróxido de calcio (en una
torre de absorción). Como resultado de esta operación, el CO2 es absorbido
en la fase líquida y la fase gaseosa queda libre de este compuesto.
Otro ejemplo de absorción que está más a nuestro alcance se presenta en
la fabricación de bebidas gaseosas, donde el CO2 es absorbido en una mezcla
líquida (jarabe) a presiones moderadas y temperaturas bajas.
El proceso contrario a la absorción, es decir el paso de un componente de
la fase líquida a la fase gaseosa se conoce con el nombre de desorción.

444.1.
..11.. DDDatos
aattooss ddee ee
de qquuiilliibbrriioo
equilibrio

Para realizar el diseño y/o dimensionamiento de una torre de absorción

debemos contar con los datos de equilibrio gas – líquido para el sistema que
estemos tratando.
En esta materia el cálculo de los equipos de absorción es simplificado
pues sin importar cuan complejas sean las fases gaseosa, líquida y el solu-
to, siempre se asume que se trata de tres compuestos: uno para la fase
líquida, otro para la fase gaseosa y otro para el soluto.
Para los casos más comunes los datos de equilibrio están disponibles en
bibliografía. Sin embargo, para llevar a cabo un diseño confiable, estos
datos deben ser obtenidos experimentalmente, en condiciones similares a las
del equipo industrial. Para ello se colocan cantidades conocidas de las dos
fases en un recipiente cerrado a la temperatura a la que se operará en el
equipo industrial, entonces se agita repetidamente la mezcla hasta que se
alcanza el equilibrio (es decir hasta que no varía la presión).
Una vez alcanzado el equilibrio se efectúa un análisis químico de las dos
fases, con lo que se tiene un dato de equilibrio. El experimento se repite
empleando diferentes proporciones de gas y líquido llegando a conformar así
una tabla o gráfica de equilibrio para una temperatura dada.

..11..11.. LL
444.1.1. aa ll
La eey ddde
ley ee HHHenry
eennrryy

Cuando las concentraciones son bajas y no se cuenta con datos experimen-

tales de equilibrio, se puede emplear como primera aproximación la ley de
- 102 - Hernán Peñaranda V.

Henry, esta ley señala que la presión parcial de un compuesto “A” en la fase
gaseosa es directamente proporcional a su fracción molar en la fase líquida,
siendo la constante de proporcionalidad la constante de Henry, es decir:
PA = Hx A (1)
Donde “H” (la constante de Henry) está en unidades de presión, normalmen-
te atmósferas. Si ambos lados de la ecuación (1) se dividen entre la presión
total (P), se tiene:
PA H
= xA
P P
y A = H ' xA (2)
Donde H’ es la constante adimensional de Henry. Gráficamente las curvas
de equilibrio de los sistemas que siguen la ley de Henry son simplemente
líneas rectas que pasan por el origen, tal como se muestra en la siguiente
figura:

y o PA

Sistemas que obedecen

la ley de Henry

4.1.1.1. Ejemplo

Calcule los volúmenes de CO2 por volumen de agua a 293.15 K, en el siste-

ma CO2 – agua, que está en equilibrio a 2 atm de presión.
Solución:
La constante de Henry para el sistema CO2 – H2O a 293.15 K (obtenida del
apéndice A de Geankoplis) es 0.142x104 atm, entonces la fracción molar de
CO2 en el agua es:
xA = PA/H = 2 0.142E4 / = 1.40845070423E-3
Por lo tanto existen 0.00141 moles de CO2 por cada mol de solución. To-
mando como base 1 gmol de solución, el volumen ocupado por los moles de CO2
disueltos en condiciones normales de presión y temperatura es:
VCO2 = 1.40845070423E-3 R CONST 1_l*atm/(K*gmol) CONVERT UVAL * StdT CONST
UVAL * StdP CONST UVAL / = 3.15691510931E-2 litros
Debido a la pequeña cantidad de CO2 se puede considerar que el volumen
del líquido es básicamente el volumen de agua. El volumen del agua es
prácticamente 1 litro por kg, entonces el volumen de 1 gmol de agua (18.02
g) es:
VH2O = 18.02 1000 / = 0.01802 litros
TORRES DE ABSORCIÓN POR ETAPAS - 103 -

Por lo tanto los volúmenes de CO2 por volumen de agua son:

3.15691510931E-2 0.01802 / = 1.7518951772
Es decir existen 1.75 volúmenes de CO2 por volumen de agua.

4.1.1.2. Ejercicios

1. Una fábrica de gaseosas requiere que su producto tenga 5 volúmenes de

CO2 por volumen de agua a 4 °C de temperatura. Calcule la presión a la
que deberá trabajar el equipo de absorción para lograr este objetivo
(la constante de Henry a 4°C es 0.08528x104 atm).

2. Repita el ejercicio anterior si la temperatura es de 20°C.

444.2.
..22.. AAbbssoorrcciióónn ee
Absorción nn uu
en nnaa ee
una ttaappaa dd
etapa ee ee
de qquuiilliibbrriioo
equilibrio

En ocasiones el proceso de absorción se lleva a cabo en un equipo que

tiene cierta eficiencia con relación a una etapa de equilibrio. Para dimen-
sionar o diseñar tales equipos se realizan los cálculos en una etapa de
equilibrio y se emplea la eficiencia del equipo para convertir los resulta-
dos ideales en resultados reales.
En la operación de absorción se asume que nada del líquido inerte (el
solvente de la fase líquida) pasa a la fase gaseosa y que nada del gas iner-
te (el solvente de la fase gaseosa) pasa a la fase líquida. Esto por supues-
to no es del todo cierto, sin embargo, las proporciones de líquido inerte
que pasan a la fase gaseosa y de gas inerte que pasan a la fase líquida son
casi siempre muy pequeñas, razón por la cual esta simplificación es válida
en la práctica.
Esquemáticamente una etapa de equilibrio se representa como se muestra en
la siguiente figura:

V1 V0
1
L0 L1

Donde “L” es el flujo de líquido y “V” es el flujo de gas. En esta etapa

las corrientes de salida L1 y V1 están en equilibrio.
Denominaremos “x” a la fracción de soluto (molar o másica) en la fase
líquida y “y” a la fracción de soluto (molar o másica) en la fase gaseosa.
Igualmente denominaremos “A” al soluto, “B” al líquido inerte (solvente de
la fase líquida) y “C” al gas inerte (solvente de la fase gaseosa).
Realizando el balance general de materia en la etapa de equilibrio resul-
ta:
L0 + V0 = L1 + V1 = M (3)
Donde “M” es el flujo total de gas y líquido que ingresa y sale de la
etapa. Del balance de materia para el soluto se tiene:
L0 x0 + V0 y0 = L1 x1 + V1 y1 = A (4)
Donde “A” es el flujo total de soluto que ingresa y sale de la etapa.
En este sistema normalmente se conocen los flujos y composiciones de las
dos corrientes de entrada y el problema radica en calcular los flujos y
composiciones de las corrientes de salida. En consecuencia existen cuatro
- 104 - Hernán Peñaranda V.

incógnitas (L1, V1, x1, y1) y sólo dos ecuaciones, por lo tanto se requieren
dos ecuaciones más para resolver el sistema.
La primera de esas dos ecuaciones resulta de haber asumido que nada de
gas inerte pasa a la fase líquida y que igualmente nada de líquido inerte
pasa a la fase gaseosa, en consecuencia el flujo de líquido inerte (y por
supuesto el de gas inerte) es constante:
B = L0 (1 − x0 ) = L1 (1 − x1 ) = ctte. (5)
La cuarta ecuación es la relación de equilibrio. Sin importar que los da-
tos de equilibrio se encuentren en forma tabular, gráfica o como ecuación,
nos permite calcular el valor de “y” para un valor conocido de “x”:
y = f x ( x) (6)
Donde fx es la relación de equilibrio en la forma en que se encuentre. La
función de equilibrio contraria, es decir aquella que nos permite calcular
el valor de “x” para un valor conocido de “y” es:
x = f y ( y) (7)

En la mayoría de los casos prácticos, los datos de equilibrio se encuen-

tran en forma tabular o gráfica pues son obtenidos experimentalmente. Si los
datos de equilibrio se encuentran en forma gráfica pueden ser llevados a la
forma tabular leyendo los datos más significativos de la gráfica y formando
así una tabla.
Debido a que con frecuencia los problemas se resuelven de manera iterati-
va, es útil contar con un programa que empleando los datos tabulados, nos
devuelva el valor de la variable dependiente para un determinado valor de la
variable independiente. Ha este tipo de programas se conoce con el nombre de
generadores de funciones arbitrarias y básicamente son programas de interpo-
lación, siendo el método de interpolación más empleado el de interpolación
lineal.
Puesto que esta es una situación que se nos presentará con frecuencia en
las diferentes operaciones que estudiemos, elaboraremos primero un programa
para llevar a cabo la interpolación (o extrapolación) lineal en base a dos
puntos conocidos: (x1,y1) y (x2,y2).
En la interpolación lineal se emplea la ecuación de la línea recta:
y = a + bx (8)
Donde “x” es la variable independiente y “y” la variable dependiente. Las
constantes “b” y “a” se calculan con:
y1 − y 2
b= (9)
x1 − x 2

a = y 1 − bx 1 (10)
El algoritmo para llevar a cabo la interpolación lineal se presenta en el
diagrama de actividades de la siguiente página. El programa elaborado en
base a dicho algoritmo (añadido al archivo Metodos1.s) es el siguiente:
xNAME inlin ( Interpolación lineal )
:: ( Datos en la pila: x: valor a interpolar; x1,x2: valores )
CK5&Dispatch ( conocidos de la variable dependiente; y1,y2: valores )
# 11111 ( conocidos de la variable dependiente )
::
OVER SWAP %- 4PICK 4ROLL %- %/ ( b )
DUP 4ROLL %* ROT SWAP %- ( a )
TORRES DE ABSORCIÓN POR ETAPAS - 105 -

3UNROLL %* %+ ( y )
;
;
inlin: Interpolación lineal.

x: Valores a interpolar.
recibir x, x1, x2, y1, y2 x1, x2: Valores conocidos de la
variable independiente.
y1, y2: Valores conocidos de la
b = (y1-y2)/(x1-x2)
variable dependiente.

a = y1-bx1

devolver a+bx Valor interpolado

Por ejemplo dados los siguientes datos:

x Y
12 14.5
16 19.6
Podemos interpolar el valor de “y” para “x=15”:
y = 15 12 16 14.5 19.6 inlin = 18.325
Igualmente podemos extrapolar el valor de “y” para “x=18”:
y = 18 12 16 14.5 19.6 inlin = 22.15
Por supuesto también podemos interpolar el valor de “x” para un valor co-
nocido de “y”, así para “y=17” el valor de “x” es:
y = 17 14.5 19.6 12 16 inlin = 13.9607843137
Ahora podemos elaborar un programa que busque el segmento de interpola-
ción y una vez ubicado realice la interpolación lineal con el programa in-
lin. El algoritmo de este programa se presenta en el diagrama de actividades
de la siguiente página y el código del programa elaborado en base al mismo
(y escrito en el archivo metodos1.s) es el siguiente:
xNAME fun1 ( generador de funciones arbitrarias de una variable )
:: ( dato en la pila: x: valor de la variable independiente )
CK3&Dispatch ( vx: vector independiente; vy: vector dependiente)
# 144
:: ( 501 = dimensión inválida )
SWAP DUP MDIMSDROP ( n )
DUP #2 #<
IT :: 2DROP # 501 DO#EXIT ; ( si n<2 => error )
#2 ( j = 2; pila: x, vy, vx, n, j )
BEGIN
DUP 3PICK #< ( j < n )
6PICK 5ROLL 4PICK PULLREALEL ( x[j] )
SWAP 6UNROLL %> AND ( x > x[j] )
WHILE
#1+ ( j = j+1 )
REPEAT
SWAP DROP DUP #1- ( i = j-1; ; pila: x, vy, vx, j, i )
ROT OVER PULLREALEL ( x[i] )
SWAP 4PICK PULLREALEL SWAP DROP ( x[j] )
5ROLL 4ROLL PULLREALEL ( y[i] )
SWAP 5ROLL PULLREALEL SWAP DROP ( y[j] )
- 106 - Hernán Peñaranda V.

xinlin ( y )
;
;
fun1: Generador de funciones
arbitrarias de una variable.
x: Valor de la variable independiente.
recibir x, vx, vy
vx: Vector con los datos de la variable
independiente.
n = Nº de elementos en x vy:Vector con los datos de la variable
dependiente.
[n<2]

[else]
j=2
Dimensión inválida
[else]

[ (x>xj) y (j<n) ]
j = j+1 generar error

i = j-1

Valor de la variable
devolver inlin(x,xi,xj,yi,yj)
dependiente.

Por ejemplo si los datos son los siguientes:

X Y
1 1.1
5 7.2
7 9.8
9 14.5
12 16.7
16 19.6
20 25.5
La función arbitraria se crea simplemente guardando estos datos en varia-
bles, por ejemplo: vx y vy:
[ 1 5 7 9 12 16 20 ] ‘vx’ STO
[ 1.1 7.2 9.8 14.5 16.7 19.6 25.5 ] ‘vy’ STO
Y llamando al programa “fun1”, así por ejemplo los valores de “y” para
“x=8, 20, 21, 1.2 y 0.5” son:
y(8) = 8 vx vy fun1 = 12.15
y(20) = 20 vx vy fun1 = 25.5
y(21) = 21 vx vy fun1 = 26.975 (extrapolación)
y(1.2) = 1.2 vx vy fun1 = 1.405
y(0.5) = 0.5 vx vy fun1 = 0.3375 (extrapolación)
Igualmente los valores de “x” para “y=0.5, 10, 17, 23 y 30” son:
x(0.5) = 0.5 vy vx fun1 = .606557377049 (extrapolación)
x(0.5) = 10 vy vx fun1 = 7.08510638298
TORRES DE ABSORCIÓN POR ETAPAS - 107 -

x(0.5) = 17 vy vx fun1 = 12.4137931034

x(0.5) = 23 vy vx fun1 = 18.3050847458
x(0.5) = 30 vy vx fun1 = 23.0508474576 (extrapolación)
Otra alternativa consiste en ajustar los datos a una ecuación analítica y
emplear la ecuación resultante como la función de equilibrio: fx. El princi-
pal inconveniente de esta alternativa radica en encontrar la forma de la
ecuación que ajusta correctamente los datos tabulados. Además se debe veri-
ficar que el comportamiento de la ecuación resultante sea coherente, más
allá del intervalo de los datos, pues no es raro que la ecuación sea empleda
para obtener valores fuera del intervalo original (extrapolación).
Con el programa “fun1” estamos ahora en condiciones de resolver el siste-
ma de 4 ecuaciones con 4 incógnitas sin importar que los datos de equilibrio
se encuentren en forma tabular, gráfica o analítica.
Entonces para el caso más frecuentemente, donde se conocen los flujos y
composiciones de entrada, el sistema de 4 ecuaciones puede ser colocado en
función de una sola variable: x1, la fracción de soluto en la fase líquida:
f ( x1 ) = L1 x1 + V1 y1 − A = 0
y1 = f x ( x1 )
B
L1 =
1 − x1
V1 = M − L1 (11)
M = L0 + V0 = ctte.
B = L0 (1 − x0 ) = ctte.
A = L0 x0 + V0 y0 = ctte.
Que puede ser resuelta con uno de los métodos iterativos como el de Regu-
la – Falsi, empleando como valores límite 0 y la máxima concentración posi-
ble de soluto en la fase líquida, la cual corresponde al valor de “x” en
equilibrio con el valor inicial de soluto en la fase gaseosa (fy(y2). El
algoritmo para resolver esta ecuación ºes el siguiente:
ab1c1: Absorción en una etapa
de equilibrio.

recibir L0, X0, V0, Y0, fx, fy L0: Flujo inicial de líquido (kgmol/h).
x0: Fracción molar inicial de soluto en el líquido.
V0: Flujo inicial degas (kgmol/h).
M = L0+V0 y0: Fracción molar inicial de soluto en el gas.
fx, fy : Funciones de equilibrio: y = fx(x); x = fy(y)
B = L0(1-X0)

A = L0x0+V0y0

fx1 fx1:: Función de x1.

x1 = Refa(fx1,0,fy(y0))

devolver L1, x1, V1, y1

Donde el submódulo “fx1” tiene el siguiente diagrama de actividades:

- 108 - Hernán Peñaranda V.

fx1: Submódulo de ab1c1.

x1: Fracción molar final de soluto

recibir x1
en el líquido.

y1 = fx(x1)

L1 = B/(1-x1)

V1 = M-L1

devolver L1x1+V1y1-A

Para escribir los programas de este capítulo añada al proyecto “PRQ205”

el archivo “Absorción.s” (Project -> New Source File). El código elaborado
en base a estos algoritmos y escrito en el archivo “Absorción.s” es el si-
guiente:
xNAME ab1c1 ( absorción en una etapa de equilibrio: caso 1 )
::
%6 CKNNOLASTWD DROP
CK&DISPATCH1 ( Datos: L0,x0= Flujo y fracción inicial en el líquido )
# 11188 ( V2,y2= Flujo y fracción inicial en el gas )
:: ( fx = Función de equilibrio y=f[x] )
6ROLL ( fy = Función de equilibrio x=f[y] )
CK&DISPATCH1
ONE
::
6UNROLL
%0 SEVEN NDUPN DROP
{ LAM L0 LAM x0 LAM V0 LAM y0
LAM fx LAM fy LAM M LAM B LAM A
LAM x1 LAM y1 LAM L1 LAM V1 }
BIND
LAM L0 LAM V0 %+ ' LAM M STO
LAM L0 %1 LAM x0 %- %* ' LAM B STO
LAM L0 LAM x0 %* LAM V0 LAM y0 %* %+ ' LAM A STO
'
:: ( función de x1 )
' LAM x1 STO
LAM x1 LAM fx EVAL ' LAM y1 STO
LAM B %1 LAM x1 %- %/ ' LAM L1 STO
LAM M LAM L1 %- ' LAM V1 STO
LAM L1 LAM x1 %* LAM V1 LAM y1 %* %+ LAM A %-
;
%0 LAM y0 LAM fy EVAL xRefa DROP
LAM L1 LAM x1 LAM V1 LAM y1
ABND
;
ZERO
:: ( # 202 = Bad Argument Type )
6UNROLL # 202 DO#EXIT
;
;
;
TORRES DE ABSORCIÓN POR ETAPAS - 109 -

444.2.1.
..22..11.. EEEjemplos
jjeemmppllooss

4.2.1.1. Absorción de CO2 en una etapa de equilibrio

Se empleará un equipo de absorción equivalente a una etapa de equilibrio

para absorber CO2 de una mezcla de CO2-Aire empleando agua pura a 273.2 K.
El gas que ingresa al equipo tiene una presión de 3 atm, una fracción molar
de soluto igual a 0.45 y un flujo molar igual a 200 kgmol/h. El flujo molar
del agua es de 350 kgmol/h. Asumiendo que el sistema obedece la ley de Henry
calcule los flujos y composiciones de las corrientes de salida.
Solución
Los datos con los que se cuenta son los siguientes:
L0 = 350 kgmol/h
x0 = 0
V0 = 200 kgmol/h
y0 = 0.45
P = 3 atm
H = 0.0728x104 atm (del apéndice A3 de Geankoplis)
Entonces
H’ = H/P = 0.0728E4 3 / = 242.666666667
Como el sistema obedece a ley de Henry las funciones de equilibrio son:
fx = 242.666666667*x
fy = y/242.666666667
Se tiene entonces todos los datos para resolver el problema empleando el
programa ab1c1:
« 350. 0. 200. .45
« 242.666666667 * »
« 242.666666667 / »
ab1c1
4. →LIST
{ L1 x1 V1 y1} →TAG
LIST→ DROP
»
Haciendo correr el programa (con EVAL) se obtiene:
L1: 350.647657526
x1: 1.84703223414E-3
V1: 199.352342474
y1: .448213155485
Que son los flujos y composiciones de salida requeridos.

4.2.1.2. Absorción de metanol en una etapa de equilibrio

Se emplearán 1000 kgmol/h de agua pura en un equipo de absorción equiva-

lente a una etapa de equilibrio, para recuperar metanol de 400 kgmol/h de
una mezcla de metanol-aire que se encuentra a una presión de 720 mm de Hg y
una temperatura de 45°C. La presión parcial del metanol en esta mezcla es de
85 mm de Hg. Para este sistema se cuenta con los siguientes datos de equili-
brio:
XA PA (39.9°C) PA (59.4°C)
0.00 0 0
0.05 25.0 50
- 110 - Hernán Peñaranda V.

0.10 46.0 102

0.15 66.5 151
Calcule los flujos y composiciones de salida así como el porcentaje de
metanol recuperado.
Solución
Los datos con los que contamos son los siguientes:
L0 = 1000 kgmol/h
x0 = 0
V0 = 400 kgmol/h
PA0 = 85 mmHg
P = 720 mmHg
=> y0 = 85 720 / 5 RND = .11806
Se tienen además los datos de equilibrio aunque no a la temperatura del
problema, por lo que deben ser interpolados para 45°C:
« [ 0 25 46 66.5 ] [ 0 50 102 151 ]
→ PA1 PA2
« 1 4 FOR i
45 39.9 59.4 PA1 i GET PA2 i GET
inlin 1 RND
NEXT
4 →ARRY
»
»
Haciendo correr el programa se obtiene:
[ 0. 31.5 60.6 88.6 ]
O dividiendo entre la presión de operación:
YA = [ 0. 31.5 60.6 88.6 ] 720 / 5 RND = [ 0. .04375 .08417 .12306 ]
Por lo tanto los datos de equilibrio para la temperatura de operación
son:
XA PA (45°C) YA
0.00 0.00 0.00000
0.05 31.5 0.04375
0.10 60.6 0.08417
0.15 88.6 0.12306
Se cuenta entonces con todos los datos para resolver el problema emplean-
do el programa ab1c1. Para no escribir dos veces los datos de equilibrio es
conveniente guardarlos en dos variables:
[ 0 0.05 0.10 0.15 ] 'vx' STO
[ 0 .04375 .08417 .12306 ] 'vy' STO
Entonces las funciones de equilibrio son:
fx = « vx vy fun1 »
fy = « vy vx fun1 »
Las composiciones y flujos de salida son:
« 1000 0 400 0.11806
« vx vy fun1 »
« vy vx fun1 »
ab1c1
4 →LIST 5 RND
{ L1 x1 V1 y1 } →TAG
TORRES DE ABSORCIÓN POR ETAPAS - 111 -

LIST→ DROP
»
Haciendo correr el programa se obtiene:
L1: 1036.12354
x1: .03486
V1: 363.87646
y1: .03051
Entonces el porcentaje de metanol recuperado es:
1036.12354 .03486 * 1000 0.11806 * / 100 * 2 RND = 30.59%
_________________________________________________________
Si bien el caso donde se conocen todos los flujos y composiciones de en-
trada es el caso más frecuente, en una etapa de equilibrio existen 4 paráme-
tros de entrada y 4 de salida, por lo tanto se pueden presentar alrededor de
24 casos, no resulta práctico entonces analizar y programar cada uno de es-
tos casos, por lo que para los casos menos frecuentes se debe recurrir a las
ecuaciones (3) a (7) y resolverlas para el caso en particular.
Analicemos por ejemplo otro caso frecuente: aquel donde se conoce la com-
posición y flujo iniciales de la fase gaseosa (V0, y0), la composición ini-
cial de la fase líquida (x0) y el porcentaje de soluto a absorber (pa).
Con el porcentaje a absorber podemos calcular el flujo final de soluto en
la fase gaseosa (V1y1):

⎛ V y −V y ⎞ ⎛ Vy ⎞
pa = ⎜ 0 0 1 1 ⎟100 = ⎜1 − 1 1 ⎟100
⎝ V y
0 0 ⎠ ⎝ V0 y0 ⎠
V1 y1 pa
= 1− (12)
V0 y0 100
⎛ pa ⎞
V1 y1 = ⎜1 − ⎟ V0 y0
⎝ 100 ⎠
La composición final en la fase gaseosa puede ser calculada tomando en
cuenta que el flujo de gas inerte (C), al igual que de líquido inerte (B),
es constante:
C = V0 (1 − y0 ) = V1 (1 − y1 ) = ctte.
C = V1 − V1 y1 (13)
V1 = C + V1 y1
Entonces la fracción final de soluto en la fase gaseosa (y1) se calcula
con:
V1 y1
y1 = (14)
V1
Puesto que los flujos finales están en equilibrio x1 puede ser calculado
con la relación de equilibrio:
x1 = f y ( y1 ) (15)

Del balance de materia (ecuación 3) se tiene:

V0 − V1 = L1 − L0 = D (16)
Del flujo de líquido inerte constante (ecuación 5), se puede despejar el
flujo final de líquido (L1):
- 112 - Hernán Peñaranda V.

1 − x0
L1 = L0 (17)
1 − x1
Y reemplazando esta ecuación en la ecuación (16) obtenemos L0:

1 − x0 ⎛ 1 − x0 ⎞ ⎛x −x ⎞
L0 − L0 = L0 ⎜ − 1⎟ = L0 ⎜ 1 0 ⎟=D
1 − x1 ⎝ 1 − x1 ⎠ ⎝ 1 − x1 ⎠
(18)
⎛ 1 − x1 ⎞
L0 = ⎜ ⎟D
⎝ x1 − x0 ⎠
Con lo que el sistema queda resuelto y como se puede ver, en este caso,
la solución no es iterativa. El diagrama de actividades para automatizas el
cálculo de este caso es el siguiente:
ab1c2: Absorción en una etapa
de equilibrio.

recibir x0, V0, y0, pa, fy x0: Fracción molar inicial de soluto en el líquido.
V0: Flujo inicial degas (kgmol/h).
y0: Fracción molar inicial de soluto en el gas.
v1y1 = (1-pa/100) V0Y0 pa: Porcentaje de soluto a absorber.
fy : Función de equilibrio: x = fy(y).
C = V0(1-y0)

V1 = C+v1y1

y1 = v1y1/V1

x1 = fy(y1)

D = V0-V1

L0 = (1-x1)/(x1-x0)D

L1 = L0(1-x0)/(1-x1)

devolver L0, L1, x1, V1, y1

Y el programa elaborado en base al mismo es el siguiente:

xNAME ab1c2 ( absorción en una etapa de equilibrio: caso 2 )
::
CK5&Dispatch ( Datos: x0= Fracción inicial en el líquido )
# 11118 ( V2,y2= Flujo y fracción inicial en el gas )
:: ( pa= Porcentaje de soluto a absorber )
%0 EIGHT NDUPN DROP ( fy= Función de equilibrio x=fy[y] )
{ LAM x0 LAM V0 LAM y0 LAM pa LAM fy
LAM v1y1 LAM C LAM V1 LAM y1 LAM x1
LAM D LAM L0 LAM L1 }
BIND
%1 LAM pa %100 %/ %- LAM V0 LAM y0 %* %* ' LAM v1y1 STO
LAM V0 %1 LAM y0 %- %* ' LAM C STO
LAM C LAM v1y1 %+ ' LAM V1 STO
LAM v1y1 LAM V1 %/ ' LAM y1 STO
LAM y1 LAM fy EVAL ' LAM x1 STO
TORRES DE ABSORCIÓN POR ETAPAS - 113 -

LAM V0 LAM V1 %- ' LAM D STO

%1 LAM x1 %- LAM x1 LAM x0 %- %/ LAM D %* ' LAM L0 STO
LAM L0 %1 LAM x0 %- %* %1 LAM x1 %- %/ ' LAM L1 STO
LAM L0 LAM L1 LAM x1 LAM V1 LAM y1
ABND
;
;

444.2.2.
..22..22.. EEEjemplo
jjeemmpplloo

Se empleará agua pura en un equipo de absorción equivalente a una etapa

de equilibrio, para recuperar el 99% del metanol existente en 400 kgmol/h de
una mezcla de metanol-aire que se encuentra a una presión de 720 mm de Hg y
una temperatura de 45°C. La presión parcial del metanol en esta mezcla es de
85 mm de Hg. Los datos de equilibrio son los mismos que en el ejemplo
4.2.1.2. Calcule los flujos y composiciones de las corrientes faltantes.
Los datos con los que contamos son los siguientes:
x0 = 0
V0 = 400 kgmol/h
PA0 = 85 mmHg
P = 720 mmHg => y0 = 85 720 / 5 RND = .11806
pa = 99%
Y los datos de equilibrio:
[ 0 0.05 0.10 0.15 ] 'vx' STO
[ 0 .04375 .08417 .12306 ] 'vy' STO
Se tiene entonces todos los datos para resolver el problema empleando el
programa ab1c2:
« 0 400 0.11806 99 « vy vx fun1 »
ab1c2
5 →LIST
{ L0 L1 x1 V1 y1 } →TAG
LIST→ DROP
»
Haciendo correr el programa se obtiene:
L0: 30553.3770299
L1: 30600.1287899
x1: 1.52782886376E-3
V1: 353.24824
y1: 1.33685025579E-3
Como se puede observar el flujo de agua necesario para recuperar el 99%
del metanol en una etapa de equilibrio se requiere más de 76 veces del flujo
del gas, lo que es muy ineficiente.

444.2.3.
..22..33.. EEEjercicios
jjeerrcciicciiooss

3. Se empleará un equipo de absorción equivalente a una etapa de equili-

brio para absorber SO2 de una mezcla SO2-aire, empleando agua pura a
293.2 K. El gas se encuentra a una presión de 2.51x105 Pa y la presión
parcial inicial del SO2 en la mezcla es 1.76x104 Pa. El flujo inicial
del gas es de 6.4 kgmol/h y el del agua 330 kgmol/h. Calcule los flujos
y composiciones de las corrientes de salida así como el porcentaje de
soluto recuperado. Los datos de equilibrio son los siguientes:
- 114 - Hernán Peñaranda V.

xA PA (mm Hg)
0.0000000 0.0
0.0000562 0.5
0.0001403 1.2
0.0002800 3.2
0.0004220 5.8
0.0005640 8.5
0.0008420 14.1
0.0014030 26.0
0.0019650 39.0
0.0027900 59.0
0.0042000 92.0
0.0069800 161.0
0.0138500 336.0
0.0206000 517.0
0.0273000 698.0

4. Repita el ejercicio 3 asumiendo que el sistema obedece la ley de Henry.

La constante de Henry a 293.2 K es 18.09 atm.

5. Si en el sistema del ejercicio 3 se quiere recuperar el 70% del SO2 di-

suelto. Calcule los flujos inicial y final de agua así como las frac-
ciones molares en las fases gaseosa y líquida.

6. Si en el sistema del ejercicio 3 la fracción molar final de soluto en

la fase gaseosa es 7x10-3. Calcule los flujos inicial y final de agua,
las fracciones molares en las fases gaseosa y líquida y el porcentaje
de soluto recuperado.

7. Si en el sistema del ejercicio 3 se emplean 2 equipos de absorción,

cada uno equivalente a una etapa de equilibrio, donde en el primer
equipo se recupera el 30% de soluto, en el segundo el 25% y a ambos se
alimenta agua pura a 293.2 K. ¿Cuáles serían los flujos y composiciones
en cada uno de los equipos?.

8. Si en el sistema del ejercicio 3 se emplean 3 equipos de absorción,

donde en la fase líquida a cada equipo se alimenta 120 kgmol/h de agua
pura. En la fase gaseosa, el gas resultante de la primera etapa se ali-
menta a la segunda y el resultante de la segunda se alimenta a la ter-
cera. Calcule los flujos y composiciones faltantes así como el porcen-
taje de soluto recuperado.

9. Se empleará un equipo de absorción equivalente a una etapa de equili-

brio para absorber el 90% del amoniaco contenido en 65 kgmol/h de una
mezcla gaseosa de amoniaco – aire, que contiene 0.12 en fracción molar
de amoniaco. La presión del gas es de 1.723x105 Pa y para la absorción
se emplea agua pura a 293.2 K. Calcule los flujos y composiciones fal-
tantes. Los datos de equilibrio para este sistema son:

xA PA (mm Hg)
0.0000 0.0
0.0208 12.0
0.0258 15.0
0.0309 18.2
0.0405 24.9
0.0503 31.7
0.0737 50.0
TORRES DE ABSORCIÓN POR ETAPAS - 115 -

0.0960 69.6
0.1370 114.0
0.1750 166.0
0.2100 227.0
0.2410 298.0
0.2970 470.0

10. Se empleará un equipo de absorción equivalente a una etapa de equili-

brio para reducir el porcentaje de acetona contenido en 40 kgmol/h, de
una mezcla gaseosa de acetona – aire, desde 5% hasta 1%. Para ese fin
se empleará agua pura a 293.2 K. Calcule los flujos y composiciones
faltantes así como el porcentaje de acetona recuperado. Los datos de
equilibrio son los siguientes:

xA PA (mm Hg)
0.0000 0.0
0.0333 30.0
0.0720 62.8
0.1170 85.4
0.1710 103.0

11. Repita el ejercicio anterior asumiendo que el sistema obedece la ley de

Henry, siendo la constante de Henry igual a 2.53 atm.

444.3.
..33.. TToorrrreess dd
Torres ee aa
de bbssoorrcciióónn cc
absorción oonn mm
con úúllttiipplleess ee
múltiples ttaappaass ee
etapas nn cc
en oonnttrraaccoorrrriieennttee
contracorriente

En la mayoría de los casos prácticos no resulta económicamente factible

emplear una etapa de equilibrio para lograr la separación requerida. Una
disposición más eficiente se logra si se emplea un equipo con múltiples eta-
pas donde las fases gaseosa y líquida fluyen en contracorriente.
El esquema general de un equipo de múltiples etapas en contracorriente es
el siguiente:

V1 V2 Vi Vi+1 Vn V0
1 ... i ... n
L0 L1 Li-1 Li Ln-1 Ln
B
A
Donde se han identificado dos sistemas: “A” con las “n” etapas del equipo
y “B” con las primeras “i” etapas del equipo.
En un equipo con múltiples etapas de equilibrio se introduce una nueva
variable: el número de etapas “n”, por lo tanto se requiere un dato más para
resolver el sistema.
Analicemos primero el sistema “A”. En este sistema, dependiendo de los
datos conocidos, se presentan múltiples casos siendo el más frecuente aquel
donde se conocen los flujos y composiciones de entrada (L0, x0, V0, y0) y el
porcentaje de soluto absorbido en el equipo (pa).
El balance general de materia en este sistema es:
L0 + V0 = Ln + V1 = M (19)
Donde “M” es el flujo total de masa que ingresa (y sale) del equipo. El
balance de materia para el soluto es:
- 116 - Hernán Peñaranda V.

L0 x0 + V0 y0 = Ln xn + V1 y1 = A (20)
Donde “A” es el flujo total de soluto que ingresa y sale del equipo. Al
igual que en una etapa de equilibrio, se asumen constantes los flujos de
líquido (B) y gas inerte (C), por lo tanto:
B = L0 (1 − x0 ) = L1 (1 − x1 ) = ... = Ln (1 − xn ) = ctte. (21)

C = V0 (1 − y0 ) = V1 (1 − y1 ) = ... = Vn (1 − yn ) = ctte. (22)

Por supuesto se debe contar también con los datos de equilibrio:

y = f x ( x) (6)

x = f y ( y) (7)

Entonces, para el caso más frecuente, con el porcentaje de soluto absor-

bido (pa) se puede calcular el flujo final de soluto en la fase gaseo-
sa(V1y1):

⎛ pa ⎞
V1 y1 = ⎜1 − ⎟ V0 y0 (12)
⎝ 100 ⎠
Y con el flujo de gas inerte (C) el flujo final de la fase gaseosa (V1):
V1 = C + V1 y1 (13)
Entonces la fracción final de soluto en la fase gaseosa (y1) se calcula
con:
V1 y1
y1 = (14)
V1
Del balance general de materia obtenemos el flujo final de la fase líqui-
da (LN):
Ln = M − V1 (23)
Finalmente, del balance de materia para el soluto obtenemos la fracción
final de soluto en la fase líquida (xn):
A − V1 y1
xn = (24)
Ln
Entonces la solución del sistema “A” para el caso más frecuente no es
iterativa. El algoritmo que automatiza estos cálculos se presenta en el dia-
grama de actividades de la siguiente página y el código elaborado en base al
mismo es el siguiente:
xNAME abnc1 ( absorción en múltiples etapas de equilibrio: caso 1)
::
CK5&Dispatch
# 11111
::
%0 EIGHT NDUPN DROP
{ LAM L0 LAM x0 LAM V0 LAM y0 LAM pa
LAM M LAM A LAM v1y1 LAM C LAM V1
LAM y1 LAM Ln LAM xn }
BIND
LAM L0 LAM V0 %+ ' LAM M STO
LAM L0 LAM x0 %* LAM V0 LAM y0 %* %+ ' LAM A STO
%1 LAM pa %100 %/ %- LAM V0 LAM y0 %* %* ' LAM v1y1 STO
TORRES DE ABSORCIÓN POR ETAPAS - 117 -

LAM V0 %1 LAM y0 %- %* ' LAM C STO

LAM C LAM v1y1 %+ ' LAM V1 STO
LAM v1y1 LAM V1 %/ ' LAM y1 STO
LAM M LAM V1 %- ' LAM Ln STO
LAM A LAM v1y1 %- LAM Ln %/ ' LAM xn STO
LAM Ln LAM xn LAM V1 LAM y1
ABND
;
;
abnc1: Absorción en múltiples
etapas de equilibrio. Caso 1

recibir L0, x0, V0, y0, pa x0: Fracción molar inicial de soluto en el líquido.
V0: Flujo inicial degas (kgmol/h).
y0: Fracción molar inicial de soluto en el gas.
M = L0+V0 pa: Porcentaje de soluto a absorber.

C = V0(1-y0)

A = L0x0+V0y0

v1y1 = (1-pa/100)V0y0

V1 = C+v1y1

y1 = v1y1/V1

Ln = M-V1

xn = (A-v1y1)/Ln

devolver Ln, xn, V1, y1

444.3.1.
..33..11.. EEEjemplo
jjeemmpplloo

Se desea absorber el 90% de la acetona de un gas que contiene 1 mol% de

acetona en aire, empleando una torre con múltiples etapas en contracorrien-
te. El flujo inicial del gas es de 30 kgmol/h y el flujo de agua pura 90
kgmol/h. Calcule los flujos y composiciones de las corrientes que salen del
equipo.
Solución
Los datos con los que contamos son:
L0 = 90 kgmol/h
x0 = 0
V0 = 30 kgmol/h
y0 = 0.01
pa = 90
Se cuenta entonces con todos los datos para resolver el problema con el
programa abnc1:
- 118 - Hernán Peñaranda V.

« 90 0 0.01 30 90 pa x0 y0 V0 L0
« L0 x0 V0 y0 pa abnc1
4 LIST
{ Ln xn V1 y1 } TAG
LIST DROP
»
»
Haciendo correr el programa se obtiene:
Ln: 90.27
xn: 2.99102691924E-3
V1: 29.73
y1: 1.00908173562E-3
Que son los flujos y composiciones finales pedidos.

444.3.2.
..33..22.. EEEjercicios
jjeerrcciicciiooss

12. Se diseñará una torre de absorción, con múltiples etapas en contraco-

rriente, para absorber SO2 de una mezcla de SO2-aire, empleando 6000
kg/h.m2 de agua pura a 293 K. Al equipo ingresarán 150 kg/h.m2, de gas
con un 20% en peso de soluto y se desea que a la salida el porcentaje
de soluto sea reducido a 2% en peso. Calcule el porcentaje de soluto
absorbido en el equipo, los flujos de salida y el porcentaje final de
soluto en la fase líquida.

13. Se empleará una torre de absorción, con múltiples etapas en contraco-

rriente, para absorber el 95% del amoniaco contenido en 65 kgmol/h, de
una mezcla de amoniaco-aire, que contiene 20 mol% de amoniaco. Al equi-
po ingresan 80 kgmol/h de agua pura a 293.2 K. Calcule los flujos y
composiciones de las corrientes que salen del equipo.

14. Se empleará una torre de absorción, con múltiples etapas en contraco-

rriente, para absorber el 99% de la acetona contenida en 40 kgmol/h, de
una mezcla de acetona-aire que tiene un porcentaje inicial de acetona
del 5% (porcentaje molar). Calcule los flujos y composiciones de las
corrientes que salen del equipo.

Una vez resuelto el sistema “A”, podemos pasar a resolver el sistema “B”.
Del balance general de materia en este sistema se obtiene:
L0 + Vi +1 = Li + V1
(25)
L0 − V1 = Li − Vi +1 = DM
Donde “DM”, que es el diferencial entre la masa que ingresa y sale de una
etapa por un mismo lado, es constante y puede ser calculada pues los valores
de L0 y V1 se obtienen al resolver el sistema “A”. Del balance de materia
para el soluto resulta:
L0 x0 − V1 y1 = Li xi − Vi +1 yi +1 = DA (26)
Donde “DA”, que es la diferencia entre la masa de soluto que ingresa y
sale de una etapa por un mismo lado, es también una constante y puede ser
calculada pues los valores de L0, x0, V1 y y1 se obtienen al resolver el sis-
tema “A”. De esta ecuación podemos despejar la fracción de soluto que ingre-
sa a la etapa (yi+1):
Li xi − DA
yi +1 = (27)
Vi +1
TORRES DE ABSORCIÓN POR ETAPAS - 119 -

Esta ecuación es conocida como la línea de operación pues con ella se

puede trazar una línea en una gráfica x versus y, dándose valores de xi y
calculando los correspondientes valores de yi+1. Para ello se toma en cuenta
que el flujo de líquido inerte es constante (ecuación 21):
B
Li = (28)
(1 − xi )
Y que el flujo de gas que entra a la etapa (Vi+1) puede ser obtenido del
balance de materia (ecuación 25):
Vi +1 = Li − DM (29)
Con la curva de equilibrio y la línea de operación, trazadas en una
gráfica “x” versus “y”, como se muestra en la siguiente figura:
y

Línea de operación

Curva de equilibrio

Se puede determinar el número de etapas teóricas de una manera similar al

método de McCabe –Thiele. Comenzando con la fracción inicial de soluto en la
fase líquida “x0” se intersecta la línea de operación (ver la siguiente fi-
gura), lo que nos da el valor de y1, entonces con y1 se intersecta la curva
de equilibrio, lo que nos da x1, luego con x1 se intersecta la línea de ope-
ración obteniéndose y2 y con y2 se intersecta la curva de equilibrio obte-
niéndose x2. Se continúa de esta manera hasta que el valor de xi es mayor o
igual a la fracción final de soluto en la fase líquida “xn” (o hasta que el
valor de yi es mayor o igual a la fracción inicial de soluto en la fase ga-
seosa “y0”). Cuando se cumple una de estas condiciones el proceso concluye y
entonces se cuenta el número de veces que se ha repetido el proceso, lo que
corresponde al número de etapas teóricas.
- 120 - Hernán Peñaranda V.

y
y5
y0

Línea de operación
4
y4

3
y3

2
y2

1
y1 Curva de equilibrio

x
x0 x1 x2 x3 xn x4

Así por ejemplo en la figura se requieren casi 4 etapas teóricas para lo-
grar la separación.
Puesto que contamos con las ecuaciones analíticas para la línea de opera-
ción y con las funciones para los datos de equilibrio, podemos resolver el
problema empleando directamente las ecuaciones en lugar de seguir el proce-
dimiento gráfico.
Para ello comenzamos con la etapa 1, es decir i=1. En esta etapa, como
resultado del sistema “A”, se conoce la fracción final de soluto en la fase
gaseosa (y1) y como está en equilibrio con la fase líquida, podemos calcular
el valor de x1 con la función de equilibrio:
xi = f y ( yi )

Entonces con la ecuación (28) calculamos el valor de Li, con la ecuación

(29) el valor de Vi+1 y con la ecuación (27) el valor de yi+1 (y2). Ahora po-
demos repetir el proceso en la segunda etapa i=2, pues conocemos el valor de
y2, de esa manera obtendremos el valor de y3 y con este valor se puede repe-
tir el proceso para la tercera etapa, continuando de esa manera hasta que el
valor de xi sea mayor o igual a la fracción final de soluto en la fase
líquida (xn) o hasta que el valor de yi+1 sea mayor o igual a la fracción
inicial de soluto en la fase gaseosa (y0).
El algoritmo que automatiza el proceso de cálculo es el que se muestra en
el diagrama de actividades de la siguiente página, en el mismo el nuevo va-
lor de yi+1 se guardan en la variable y1 y no se conservan los valores de
yi+1, pues sólo se utilizan una vez. En Vi+1 se ha empleado la letra “j” en
lugar de “i+1”. La fracción útil de la última etapa se obtiene simplemente
como la razón entre la separación requerida (xn-xi-1) y la separación lograda
(xi-xi-1). El valor de xi-1 se almacena en la variable x0.
El código elaborado siguiendo el algoritmo es el siguiente:
xNAME abn ( Etapas de equilibrio en una torre )
:: ( de absorción en contracorriente )
%6 CKNNOLASTWD DROP ( Datos: )
CK&DISPATCH1 ( L0=13= flujo inicial de líquido; )
# 11118 ( x0=12= fracción inicial de soluto en el líquido )
:: ( xn=11= fracción final de soluto en el líquido )
6ROLL ( V1=10= flujo final de gas )
CK&DISPATCH1 ( y1= 9= fracción final de soluto en el gas )
ONE ( fy= 8= función de equilibrio x=fy[y] )
TORRES DE ABSORCIÓN POR ETAPAS - 121 -

:: ( variables locales: )
6UNROLL ( DM=7; DA=6; B=5; i=4; xi=3; Li=2; Vj=1 )
%0 SEVEN NDUPN DROP
' NULLLAM THIRTEEN NDUPN DOBIND
13GETLAM 10GETLAM %- 7PUTLAM ( DM = L0-V1 )
13GETLAM 12GETLAM %* 10GETLAM
9GETLAM %* %- 6PUTLAM ( DA = L0x0-V1y1 )
13GETLAM %1 12GETLAM %- %* 5PUTLAM ( B = L0[1-x0] )
::
BEGIN
9GETLAM 8GETLAM EVAL 3PUTLAM ( xi = fy[y1] )
3GETLAM 11GETLAM %>= IT :: 2RDROP ; ( si xi>=xn => salir )
3GETLAM 12PUTLAM ( x0=xi )
4GETLAM %1+ 4PUTLAM ( i= i+1 )
5GETLAM %1 3GETLAM %- %/ 2PUTLAM ( Li = B/[1-xi] )
2GETLAM 7GETLAM %- 1PUTLAM ( Vj = Li-DM )
2GETLAM 3GETLAM %* 6GETLAM %-
1GETLAM %/ 9PUTLAM ( y1= [Lixi-DA]/Vj )
AGAIN
;
4GETLAM 11GETLAM 12GETLAM %-
3GETLAM 12GETLAM %- %/ %+ ( n = i+[xn-x0]/[xi-x0] )
ABND
;
ZERO
::
6UNROLL # 202 DO#EXIT ( 202 = Bad Argument Type )
;
;
;
- 122 - Hernán Peñaranda V.

abn: Cálculo del número de etapas teóricas en una

torre de absorción en contracorriente.
recibir L0, x0, xn, V1, y1, fy L0 : Flujo inicial de la fase líquida.
x0 : Fracción inicial de soluto en la fase líquida.
xn : Fracción final de soluto en la fase líquida.
DM = L0-V1
V1 : Flujo final de la fase gaseosa.
y1 : Fracción final de soluto en la fase gaseosa.
DA = L0x0-V1y1 fy : Función de equilibrio x = fy(y)

B = L0(1-x0)

i=0

xi = fy(y1)

[xi>=xn]
[else]

x0 = xi

i = i+1

Li = B/(1-xi)

Vj = Li-DM

y1 = (Lixi-DA)/Vj

n = i+(xn-x0)/(xi-x0)

devolver n

Número de etapas de equilibrio.

444.3.3.
..33..33.. EEEjemplo
jjeemmpplloo

4.3.3.1. Absorción de acetona

Se desea absorber el 99% de la acetona de un gas que contiene 9 mol% de

acetona en aire, empleando una torre con múltiples etapas en contracorrien-
te. El flujo inicial del gas es de 30 kgmol/h y el flujo de agua pura 90
kgmol/h. La presión de operación de la torre es de 1.5x105 Pa y la tempera-
tura de 293.2 K. a) Calcule el número de etapas teóricas que deberá tener el
equipo para lograr la separación requerida; b) Calcule el número de etapas
reales si la eficiencia de cada etapa es del 60%. Los datos de equilibrio
son los del ejercicio 10.
Solución
Los datos con los que contamos son:
L0 = 90 kgmol/h
x0 = 0
V0 = 30 kgmol/h
y0 = 0.09
pa = 99
P = 7.2x104 Pa = 7.2E4_Pa 1_mmHg CONVERT 2 RND = 540.04_mmHg
TORRES DE ABSORCIÓN POR ETAPAS - 123 -

Entonces contamos con todos los datos para resolver el sistema “A” emple-
ando el programa abnc1:
« 90 0 30 0.09 99
abnc1
4 →LIST 5 RND
{ Ln xn V1 y1 } →TAG
LIST→ DROP
»
Haciendo correr el programa se obtiene:
Ln: 92.673
xn: .02884
V1: 27.327
y1: .00099
Que son los flujos y composiciones a la salida del equipo redondeados en
el quinto decimal.
Para las funciones de equilibrio guardamos los datos del vector “x” en la
variable “vx”:
[ 0.0000 0.0333 0.0720 0.1170 0.1710 ] 'vx' STO
Para el vector “y” dividimos el vector de presiones entre la presión to-
tal y guardamos el vector resultante en la variable “vy”:
[ 0 30 62.8 85.4 103 ] 540.04 / 4 RND 'vy' STO
a) Contamos ahora con todos los datos para calcular el número de etapas
teóricas empleando el programa abn:
« 90 0 0.02884 27.327 0.00099
« vy vx fun1 »
abn 2 RND
"n" →TAG
»
Haciendo correr el programa se obtiene:
n: 5.75
Que es el número de etapas teóricas necesarias para lograr la separación
requerida.
b) El número de etapas reales se calcula simplemente dividiendo el número
de etapas teóricas entre la eficiencia:
Número de etapas reales = 5.75 0.6 / 2 RND = 9.58

4.3.3.2. Reutilización de un equipo existente

Una industria cuenta con una torre de absorción con múltiples etapas en
contracorriente equivalentes a 12 etapas teóricas. Se quiere emplear este
equipo para absorber el SO2 contenido en una mezcla de SO2 – aire. La con-
centración inicial de SO2 en el gas es del 20 % (molar) y el flujo de aire
inerte es de 160 kg/h.m2. Para la absorción se empleará agua pura con un
flujo inicial igual a 3785 kg/h.m2. El equipo opera a una presión de 1x105
Pa. Calcule el porcentaje de SO2 de azufre que puede ser absorbido en el
equipo. Los datos de equilibrio son los del ejercicio 3.
Solución
Los datos con los que se cuenta son los siguientes:
n = 12
- 124 - Hernán Peñaranda V.

y0 = 0.2
C = 160 kg/h.m2 = 160 28.97 / 2 RND = 5.52 kgmol/h.m2
L0 = 3785 kg/h.m2 = 3785 18.02 / 2 RND = 210.04 kgmol/h.m2
x0 = 0
P = 1x105 Pa = 1E5_Pa 1_mmHg CONVERT 2 RND = 750.06_mmHg
Se tiene además los datos de equilibrio que por comodidad se guardan en
variables globales:
[ 0 0.0000562 0.0001403 0.00028 0.000422 0.000564 0.000842 0.001403 0.001965
0.00279 0.0042 0.00698 0.01385 0.0206 0.0273 ] 'vx' STO
[ 0 0.5 1.2 3.2 5.8 8.5 14.1 26 39 59 92 161 336 517 698 ] 750.06 / 7 RND
'vy' STO
Puesto que el flujo de gas inerte es constante, con “C” podemos calcular
V0:
V0 = 5.52 1 0.2 - / 2 RND = 6.9 kg/h.m2
En este caso no contamos con el porcentaje absorbido o la fracción final
de soluto en la fase gaseosa. Entonces la solución del problema es iterati-
vo: a) se asume un porcentaje de absorción (pa); b) con el valor asumido y
el programa “abnc1” se calculan los flujos y composiciones de salida (Ln,
xn, V1, y1); c) con los flujos y composiciones de salida calculados en el
inciso anterior y el programa “abn” se calcula el número de etapas teóricas;
c) se compara el número etapas calculadas con el número de etapas conocidas,
si son iguales el proceso concluye, siendo el porcentaje absorbido el valor
asumido (pa), caso contrario el proceso se repite desde el inciso (a).
Puesto que el problema iterativo es función de una sola variable: “pa”,
puede ser resuelto empleando uno de los métodos iterativos como el método de
la Secante o Regula Falsi. La función a resolver es:
f(pa) = Nº de etapas calculadas – Nº de etapas conocidas = 0
El programa que resuelve esta función, empleando el método de Regula Fal-
si, buscando la solución entre 1% y 99% es el siguiente:
« 12 0.2 160 28.97 / 3785 18.02 / 0
1E5_Pa 1_mmHg CONVERT UVAL
[ 0 0.0000562 0.0001403 0.00028 0.000422 0.000564
0.000842 0.001403 0.001965 0.00279 0.0042
0.00698 0.01385 0.0206 0.0273 ]
[ 0 0.5 1.2 3.2 5.8 8.5 14.1 26 39 59 92 161 336
517 698 ] 0 0 0 0 0 0
n y0 C L0 x0 P vx vy fy V0 Ln xn V1 Y1
« vy P / 'vy' STO « vy vx fun1 » 'fy' STO
C 1 y0 - / 'V0' STO
« 'pa' STO
L0 x0 V0 y0 pa abnc1 4 LIST
{ Ln xn V1 y1 } STO
L0 x0 xn V1 y1 fy abn n -
»
1 99 Refa "pa" TAG
»
»
Haciendo correr el programa se obtiene:
pa: 97.4632670011
Que es el porcentaje de soluto que puede ser absorbido en el equipo.
TORRES DE ABSORCIÓN POR ETAPAS - 125 -

444.3.4.
..33..44.. EEEjercicios
jjeerrcciicciiooss

15. Se empleará una torre de absorción, con múltiples etapas en contraco-

rriente, para absorber el 97% del metanol contenido en 400 kgmol/h de
una mezcla de metanol-aire que se encuentra a una presión de 1 atm y
una temperatura de 55°C. La presión parcial inicial del metanol en esta
mezcla es de 120 mm de Hg y para la absorción se emplearán 700 kgmol/h
de agua pura. a) Calcule el número de etapas teóricas para lograr la
separación requerida. b) Calcule el número de etapas reales si el ren-
dimiento de cada etapa es del 35%. Utilice los datos de equilibrio del
ejemplo 4.2.1.2.

16. Se empleará agua pura, en una torre de absorción, con múltiples etapas
en contracorriente, para reducir el porcentaje de CO2 de una mezcla de
CO2 – aire que se encuentra a 10 atm y una temperatura de 0°C. El por-
centaje inicial (molar) de CO2 en la mezcla es del 70% y se quiere re-
ducir el mismo a 20%. El flujo inicial del aire es de 250 kgmol/h y el
de agua 16100 kgmol/h. Calcule el número de etapas teóricas que deberá
tener el equipo para lograr la absorción requerida. Asuma que el siste-
ma obedece la ley de Henry, siendo la constante de Henry 0.0728x104
atm.

17. Una industria cuenta con una torre de absorción equivalente a 4.3 eta-
pas de equilibrio. Se quiere emplear esta torre para absorber amoniaco,
con agua pura, de una mezcla de amoniaco – aire que se encuentra a
1.4x105 Pa y 293.2 K. El flujo inicial de gas es de 80 kgmol/h y tiene
un porcentaje (molar) de amoniaco igual al 23%. El flujo inicial de
agua pura es de 90 kgmol/h. Calcule el porcentaje de soluto que puede
ser absorbido en esta torre. Utilice los datos de equilibrio del ejer-
cicio 9.

18. Una industria cuenta con una torre de absorción equivalente a 5.6 eta-
pas teóricas que piensa emplear para absorber el 99.4% del amoniaco
contenido en una mezcla de amoniaco – aire, que se encuentra a 0.89 atm
y 293.2 K. El flujo inicial de gas es de 78 kgmol/h y tiene un porcen-
taje (molar) de amoniaco igual a 25%. Para la absorción se empleará
agua pura a 293.2 K. Calcule el flujo de agua que deberá emplearse para
lograr la absorción requerida. Utilice los datos de equilibrio del
ejercicio 9.

444.4.
..44.. FFlluujjoo mm
Flujo íínniimmoo dd
mínimo ee ll
de ííqquuiiddoo
líquido

El flujo mínimo de líquido que puede ser empleado para lograr una absor-
ción determinada, corresponde a la máxima concentración final de soluto que
puede alcanzar la fase líquida (xmax) y ocurre cuando la línea de operación
se hace tangente a la curva de equilibrio, tal como se muestra en la si-
guiente figura:
- 126 - Hernán Peñaranda V.

Línea de operación

y0

y1

Curva de equilibrio

x0 xmax x

Para estas condiciones el número de etapas teóricas es infinito.

El flujo mínimo puede ser encontrado por tanteo, haciendo correr los pro-
gramas “abnc1” y “abn”, para diferentes valores de L0, hasta que se encuen-
tra el flujo mínimo para el cual todavía es posible calcular el número de
etapas.

444.4.1.
..44..11.. EEjjeemmpplloo:: DD
Ejemplo: eetteerrmmiinnaacciióónn dd
Determinación eell ff
del lluujjoo mm
flujo íínniimmoo pp
mínimo oorr tt
por aanntteeoo
tanteo

Calcule el flujo mínimo de agua pura que puede ser empleada para recupe-
rar el 99% de la acetona de un gas que contiene 9 mol% de acetona en aire.
El flujo inicial del gas es de 30 kgmol/h y la presión de operación de la
torre es de 7.2x104 Pa a una la temperatura de 293.2 K. Los datos de equili-
brio son los del ejercicio 10.
Solución
Los datos con los que contamos son:
x0 = 0
V0 = 30 kgmol/h
y0 = 0.09
pa = 99
P = 7.2x104 Pa = 7.2E4_Pa 1_mmHg CONVERT 2 RND = 540.04_mmHg
Se tienen también los datos de equilibrio:
[ 0.0000 0.0333 0.0720 0.1170 0.1710 ] 'vx' STO
[ 0 30 62.8 85.4 103 ] 540.04 / 4 RND 'vy' STO
La búsqueda del flujo mínimo la llevaremos a cabo con el siguiente pro-
grama, donde se recibe el flujo inicial de líquido “L0” y con el mismo se
calculan las composiciones y flujos finales empleando el programa “abnc1” y
con estos resultados se calcula el número de etapas teóricas con el programa
“abn”:
« 0 30 0.09 99 0 0 0 0
→ L0 x0 V0 y0 pa Ln xn V1 y1
«
L0 x0 V0 y0 pa
abnc1
4 →LIST
{ Ln xn V1 y1 } STO
L0 x0 xn V1 y1
TORRES DE ABSORCIÓN POR ETAPAS - 127 -

« vy vx fun1 »
abn 2 RND
"n" →TAG
»
»
El cual por comodidad se guarda en la variable “lmin” ('lmin' STO)
Para 90 kgmol/h se obtiene:
90 lmin => n:5.75
Como es un número finito de etapas vamos bajando de 10 en 10:
80 lmin => n:6.74
70 lmin => n:8.41
60 lmin => n:12.
50 lmin => n:26.97
40 lmin => No se obtiene resultado en más de 1 minuto
Para detener el programa, que ha ingresado a un ciclo infinito, debe pul-
sar la tecla ON.
Entonces como para 40 kgmol/h se requiere un flujo infinito el flujo
mínimo debe estar comprendido entre 40 y 50 kgmol/h Continuamos la búsqueda
con 45 kgmol/h:
45 lmin => No se obtiene resultado en más de 1 minuto
Por lo tanto el flujo mínimo está entre 45 y 50 kgmol/h. Continuamos la
búsqueda con 47 kgmol/h:
47 lmin => n:60.1
Ahora sabemos que el flujo mínimo está entre 45 y 47 kgmol/h. Continuamos
la búsqueda con 46 kgmol/h:
46 lmin => n:183.82
Entonces el flujo mínimo está entre 45 y 46 kgmol/h. Continuamos la
búsqueda con 45.5 kgmol/h:
45.5 lmin => No se obtiene resultado en más de 1 minuto
Por lo tanto el flujo mínimo está entre 45.5 y 46 kgmol/h. Continuamos la
búsqueda con 45.7 kgmol/h:
45.7 lmin => No se obtiene resultado en más de 1 minuto
Ahora probamos con 45.9 kgmol/h:
45.9 lmin => n:305.03
Por lo tanto el flujo mínimo está comprendido entre 45.7 y 45.9 kgmol/h.
Probamos entonces con 45.8 kgmol/h
45.8 lmin => No se obtiene resultado en más de 1 minuto
En consecuencia el flujo mínimo que se puede emplear para lograr la sepa-
ración requerida (con un dígito de precisión después del punto) es 45.9
kgmol/h.
__________________________________________________________
El flujo mínimo de líquido puede ser calculado también recurriendo a la
ecuación del balance de materia para el soluto (ecuación 20), tomando en
cuenta que para el flujo mínimo de líquido la fracción molar final de soluto
xn es máxima xmax. La fracción máxima de soluto en la fase líquida se deter-
- 128 - Hernán Peñaranda V.

mina con la función de equilibrio, pues está en equilibrio con la fracción

inicial de soluto en la fase gaseosa (y0):
xmax = f y ( y0 )

El flujo final de soluto V1y1 se calcula con la (ecuación 12):

⎛ pa ⎞
V1 y1 = ⎜1 − ⎟ V0 y0
⎝ 100 ⎠
Entonces de la ecuación del balance de materia (ecuación 7), se puede
despejar el flujo mínimo de líquido inerte:

xmax x ⎛ x x ⎞
Ln xmax − L0 x0 = Bmin − Bmin 0 = Bmin ⎜ max − 0 ⎟ = V0 y0 − V1 y1
1 − xmax 1 − x0 ⎝ 1 − xmax 1 − x0 ⎠
V0 y0 − V1 y1
Bmin = (30)
xmax x
− 0
1 − xmax 1 − x0
Finalmente el flujo mínimo de líquido se calcula con:
Lmin = Bmin (1 − x0 ) (31)

El algoritmo que automatiza el cálculo del flujo mínimo de líquido se

presenta en el siguiente diagrama de actividades.
abLmin: Cálculo del flujo mínimo de líquido en una
torre de absorción.

recibir x0, V0, y0, pa, fy x0 : Fracción inicial de soluto en la fase líquida.
V0 : Flujo inicial de la fase gaseosa.
y0 : Fracción inicial de soluto en la fase gaseosa.
xmax = fy(y0)
pa: Porcentaje de soluto absorbido.
fy : Función de equilibrio x = fy(y)
V1y1 = (1-pa/100)V0y0

Bmin = (V0y0-V1y1)/(xmax/(1-xmax)-x0/(1-x0))

devolver Bmin(1-x0)

El código elaborado en base al algoritmo es el siguiente:

xNAME abLmin ( Flujo mínimo de líquido en una torre de absorción )
:: ( Datos: )
CK5&Dispatch ( x0=8= Fracción inicial de soluto en el líquido)
# 11118 ( V0=7= Flujo inicial de gas )
:: ( y0=6= Fracción inicial de soluto en el gas )
%0 THREE NDUPN DROP ( pa=5= Porcentaje de soluto absorbido )
' NULLLAM EIGHT ( fy=4= Función de equilibrio x=fy[y] )
NDUPN DOBIND ( xmax=3; V1y1=2; Bmin=1 )
6GETLAM 4GETLAM EVAL 3PUTLAM ( xmax= fy[y0] )
%1 5GETLAM %100 %/ %-
7GETLAM 6GETLAM %* %* 2PUTLAM ( V1y1= [1-pa/100]V0y0 )
7GETLAM 6GETLAM %* 2GETLAM %- ( [V0y0-V1y1] )
3GETLAM %1 3GETLAM %- %/ ( [xmax/[1-xmax] )
8GETLAM %1 8GETLAM %- %/ %- ( -x0/[1-x0]] )
%/ 1PUTLAM ( Bmin)
1GETLAM %1 8GETLAM %- %* ( Lmin )
TORRES DE ABSORCIÓN POR ETAPAS - 129 -

ABND
;
;

444.4.2.
..44..22.. EEjjeemmpplloo:: CC
Ejemplo: áállccuulloo dd
Cálculo eell ff
del lluujjoo mm
flujo íínniimmoo ee
mínimo mmpplleeaannddoo ee
empleando ccuuaacciioonneess aa
ecuaciones nnaallííttiiccaass
analíticas

Calcule el flujo mínimo de líquido que puede emplearse para las condicio-
nes del ejemplo 4.4.1, empleando las ecuaciones analíticas.
Solución
Tenemos los mismos datos que en el ejemplo 4.4.1, sólo que ahora el pro-
blema se resuelve directamente con el programa “abLmin”:
« 0 99 0.09 30 1.5E5_Pa 1_mmHg CONVERT UVAL
[ 0.0000 0.0333 0.0720 0.1170 0.1710 ]
[ 0 30 62.8 85.4 103 ] 0
x0 pa y0 V0 P vx vy fy
« vy P / 'vy' STO « vy vx fun1 » 'fy' STO
x0 V0 y0 pa fy abLmin "Lmin" TAG »
»
Haciendo correr el programa se obtiene:
Lmin: 45.7227477652
Que es el flujo mínimo de líquido que puede ser empleado en la torre. Es-
te resultado concuerda bien con el resultado obtenido por simple tanteo en
el ejemplo 4.4.1. Como se dijo, no se puede emplear este flujo en un equipo
real, porque para el mismo se requiere una torre de altura infinita, por lo
tanto, el flujo mínimo real deberá ser siempre mayor al valor calculado.
___________________________________________________________
En ocasiones, cuando el sistema se aleja mucho del comportamiento ideal,
puede suceder que, debido a su curvatura de la línea de equilibrio, la línea
de operación se haga tangente a la curva de equilibrio antes del punto
xmax,y0. En estos casos no queda otra alternativa que encontrar el flujo
mínimo por tanteo, como se hizo en el ejemplo 4.4.1.

444.4.3.
..44..33.. EEEjercicios
jjeerrcciicciiooss

19. Se empleará una torre de absorción en contracorriente para absorber el

90% del metano contenido en 400 kgmol/h de una mezcla de metanol-aire,
que se encuentra a una presión de 1 atm y una temperatura de 52°C. La
presión parcial inicial de metanol en esta mezcla es de 110 mm de Hg.
Calcule el número de etapas teóricas requeridas si se emplea un flujo
de agua pura igual a: a) 1.1 veces el flujo mínimo; b) 1.5 veces el
flujo mínimo; c) 2 veces el flujo mínimo. Los datos de equilibrio son
los del ejemplo 4.2.1.2.

20. Se empleará una torre de absorción en contracorriente para reducir el

porcentaje de CO2 de una mezcla de CO2 – aire que se encuentra a 1.4
atm y una temperatura de 0°C. El porcentaje inicial (molar) de CO2 en
la mezcla es del 12% y se quiere reducir el mismo a 0.9%. El flujo
inicial del aire es de 250 kgmol/h. Calcule el número de etapas teóri-
cas que deberá tener el equipo para lograr la absorción requerida si se
emplea influjo de agua pura igual a: a) 1.08 veces el flujo mínimo; b)
1.2 veces el flujo mínimo; c) 1.6 veces el flujo mínimo. Asuma que el
sistema obedece la ley de Henry, siendo la constante de Henry
0.0728x104 atm.
- 130 - Hernán Peñaranda V.

21. Se empleará una torre de absorción en contracorriente para absorber el

99.5% del amoniaco contenido en una mezcla de amoniaco – aire, que se
encuentra a 0.9 atm y 293.2 K. El flujo inicial del gas es de 78
kgmol/h y tiene un porcentaje (molar) de amoniaco igual a 25%. Para la
absorción se empleará agua pura a 293.2 K. Calcule el número de etapas
teóricas que deberá tener el equipo para lograr la absorción requerida
si se emplea influjo de agua pura igual a: a) 1.1 veces el flujo míni-
mo; b) 1.5 veces el flujo mínimo; c) 4 veces el flujo mínimo. Los datos
de equilibrio para este sistema son los del ejercicio 9.

444.5.
..55.. EEccuuaacciioonneess ss
Ecuaciones iimmpplliiffiiccaaddaass pp
simplificadas aarraa ee
para ll cc
el áállccuulloo dd
cálculo eell nn
del úúmmeerroo dd
número ee ee
de ttaappaass
etapas

Aun cuando las ecuaciones que hemos empleado para el cálculo del número
de etapas son simplificadas, pues se consideran constantes los flujos de gas
y líquido inerte. Pueden ser simplificadas más aún si se asume que la curva
de equilibrio es una línea recta, es decir se asume que el sistema obedece
la ley de Henry.
Las ecuaciones resultantes para estas condiciones han sido deducidas por
Kremser y para el caso de la absorción son las siguientes:

⎡ y − H ' x0 ⎛ 1 ⎞ 1 ⎤
log ⎢ 0 ⎜1 − ⎟ + ⎥
⎣ y1 − H ' x0 ⎝ A ⎠ A⎦
n= (32)
log ( A )

Donde “H'” es la constante adimensional de la ley de Henry y “A”, es el

factor de absorción, que está dada por la siguiente expresión:

⎛ L ⎞⎛ L ⎞
A = ⎜ '0 ⎟ ⎜ ' n ⎟ (33)
⎝ H V1 ⎠ ⎝ H V0 ⎠
Como se puede observar, para calcular el número de etapas con la ecuación
de Kremser, se requiere resolver el sistema “A”, pues se emplean los flujos
y composiciones a la entrada y salida del equipo. Por lo tanto, salvo el
hecho de que no es iterativa, esta ecuación no tiene mayor ventaja con rela-
ción al cálculo etapa por etapa (programa “abn”).
Automatizaremos la solución de esta ecuación para los casos donde sea
conveniente su uso. El algoritmo que resuelve la ecuación de Kremser es el
siguiente:
abnkr: Cálculo del número de etapas teóricas en una
torre de absorción con la ecuación de Kremser.
recibir L0, x0, Ln, V0, y0, V1, y1, H L0 : Flujo inicial de la fase líquida.
x0 : Fracción inicial de soluto en la fase líquida.
Ln : Flujo final de la fase líquida.
A = ((L0/(HV1)(Ln/(HV0)) )0.5
V0 : Flujo inicial de la fase gaseosa.
y0 : Fracción inicial de soluto en la fase gaseosa.
n = log((y0-Hx0)/(y1-Hx0)(1-1/A)+1/A)/log(A) V1 : Flujo final de la fase gaseosa.
y1 : Fracción final de soluto en la fase gaseosa.
H : Constante adimensional de Henry.
devolver n

Y el código elaborado en base al mismo es:

xNAME abnKr ( Número de etapas de absorción con la ecuación de Kremser )
:: ( datos: L0=9= Flujo inicial de líquido )
TORRES DE ABSORCIÓN POR ETAPAS - 131 -

%8 CKNNOLASTWD DROP ( x0=8= Fracción inicial de soluto en el líquido)

CK&DISPATCH1 ( Ln=7= Flujo final de líquido )
# 11111 ( V0=6= Flujo inicial de gas )
:: ( y0=5= Fracción inicial de soluto en el gas )
8ROLL 8ROLL 8ROLL ( V1=4= Flujo final de gas )
CK&DISPATCH1 ( y1=3= Fracción final de soluto en el gas )
# 111 ( H=2= Constante adimensional de Henry )
::
8UNROLL 8UNROLL 8UNROLL ( A=1 )
%0
' NULLLAM NINE NDUPN DOBIND
9GETLAM 2GETLAM 4GETLAM %* %/ ( L0/[H*V1] )
7GETLAM 2GETLAM 6GETLAM %* %/ ( Ln/[H*V0] )
%* %SQRT 1PUTLAM ( A=[...*...]^0.5 )
5GETLAM 2GETLAM 8GETLAM %* %- ( y0-H*x0 )
3GETLAM 2GETLAM 8GETLAM %* %- ( y1-H*x0 )
%/ %1 1GETLAM %1/ %- %* ( [.../...][1-1/A] )
1GETLAM %1/ %+ %LOG ( log[...+1/A] )
1GETLAM %LOG %/ ( n= .../log[A] )
ABND
;
# 000
::
8UNROLL 8UNROLL 8UNROLL
# 202 DO#EXIT
;
;
;

444.5.1.
..55..11.. EEEjemplo
jjeemmpplloo

Repita el ejemplo 4.3.3.1 empleando la ecuación de Kremser.

Solución
Los datos del mencionado ejemplo son los siguientes:
L0 = 90 kgmol/h
x0 = 0
V0 = 30 kgmol/h
y0 = 0.09
pa = 99
P = 7.2x104 Pa = 7.2E4_Pa 1_mmHg CONVERT 2 RND = 540.04_mmHg
Se cuenta además con los resultados el sistema “A”, obtenidos con el pro-
grama “abnc1”:
Ln = 92.673 kgmol/h
xn = 0.0288
V1 = 27.327 kgmol/h
y1 = 0.00099
Puesto que para aplicar la ecuación de Kremser se supone que el sistema
obedece la ley de Henry, ajustamos los 3 primeros datos de equilibrio a una
línea recta. Empleamos sólo los tres primeros datos porque sabemos que los
sistemas obedecen la ley de Henry sólo cuando las concentraciones son dilui-
das.
La pendiente de la recta ajustada es la constante de Henry. El programa
que realiza el ajuste es el siguiente:
« [ 0 0.0333 0.072 ]
- 132 - Hernán Peñaranda V.

[ 0 30 62.8 ] 540.04 / 4 RND

2 COL→ '∑DAT' STO
1 XCOL
2 YCOL
LINFIT
LR
DTAG 4 RND "H" →TAG
SWAP DROP
»
Haciendo correr el programa se obtiene:
H: 1.614
Que es la constante de Henry, resultante del ajuste, para el sistema ace-
tona-agua.
Ahora se tienen todos los datos para calcular el número de etapas con la
ecuación de Kremser (programa abnKr):
« 90 0 92.673 30 0.09 27.327
0.00099 1.614
abnKr
2 RND "n" →TAG
»
Haciendo correr el programa se obtiene:
n: 5.6
Que como vemos es ligeramente inferior al número de etapas obtenido en el
ejemplo 4.3.3.1 con el programa “abn”.

..55..22.. EEEjercicios
444.5.2. jjeerrcciicciiooss

22. Se empleará una torre de absorción, con múltiples etapas en contraco-

rriente, para absorber el 99% del metanol contenido en una mezcla de
metanol-aire que se encuentra a una presión de 1 atm y una temperatura
de 49°C. La presión parcial inicial del metanol en esta mezcla es de
130 mm de Hg y para la absorción se emplearán 900 kgmol/h de agua pura.
Empleando la ecuación de Kremser calcule el número de etapas teóricas
para lograr la separación requerida. Utilice los datos de equilibrio
del ejemplo 4.2.1.2.

23. Se empleará agua pura, en una torre de absorción, con múltiples etapas
en contracorriente, para reducir el porcentaje de CO2 de una mezcla de
CO2 – aire que se encuentra a 2.1 atm y una temperatura de 0°C. El por-
centaje inicial (molar) de CO2 en la mezcla es del 20% y se quiere re-
ducir el mismo a 0.9%. El flujo inicial del aire es de 270 kgmol/h y
el del agua 1.5 veces el flujo mínimo. Empleando la ecuación de Kremser
calcule el número de etapas teóricas que deberá tener el equipo para
lograr la absorción requerida. Asuma que el sistema obedece la ley de
Henry, siendo la constante de Henry 0.0728x104 atm.

24. Una industria cuenta con una torre de absorción equivalente a 5.1 eta-
pas de equilibrio. Se quiere emplear esta torre para absorber amoniaco,
con agua pura, de una mezcla de amoniaco – aire que se encuentra a
1.6x105 Pa y 293.2 K. El flujo inicial de gas es de 100 kgmol/h y tiene
un porcentaje (molar) de amoniaco igual al 25%. El flujo inicial de
agua pura es de 120 kgmol/h. Empleando la ecuación de Kremser calcule
el porcentaje de soluto que puede ser absorbido en esta torre. Utilice
los datos de equilibrio del ejercicio 9.
TORRES DE ABSORCIÓN POR ETAPAS - 133 -

Kremser ha deducido también una ecuación simplificada para el proceso

contrario al de la absorción: la desorción. La forma de la ecuación de Krem-
ser para la desorción es la siguiente:

⎡ y0 ⎤
⎢ x0 − H ' ⎥
log ⎢
y
(1 − A) + A⎥
⎢ xn − 0 ⎥
⎢⎣ m ⎥⎦ (34)
n=
⎛1⎞
log ⎜ ⎟
⎝ A⎠

Que como vemos es similar a la ecuación (32), por lo que su solución pue-
de ser automatizada siguiendo una lógica similar a la empleada en el progra-
ma “abnKr”.
- 134 - Hernán Peñaranda V.

4.5.2.1. Absorción una etapa, caso 3.

Datos conocidos: Composición inicial y final de la fase líquida, flujo y

composición inicial de la fase gaseosa.
El diagrama de actividades es el siguiente:
ab1c3: Absorción en una etapa de
equilibrio, caso 3.

recibir x0, x1, V0, y0, fx x0: Fracción molar inicial de soluto en el líquido.
x1: Fracción molar final de soluto en el líquido.
V0: Flujo inicial de gas (kgmol/h).
y1 = fx(x1) y0: Fracción molar inicial de soluto en el gas.
fx : Función de equilibrio: y = fx(x).
C = V0(1-y0)

V1 = C/(1-y1)
V1 V2
1
DV = V1-V0
L0 L1

x1 = fy(y1)

L0 = DV(1-x1)/(x0-x1)

L1 = L0-DV

devolver L0, L1, V1, y1

El código en System RPN es:

xNAME ab1c3 ( absorción en una etapa de equilibrio, caso 3)
::
CK5&Dispatch ( Datos: x0,x1= Fracciones inicial y final en el líquido )
# 11118 ( V0,y0= Flujo y fracción inicial en el gas )
:: ( fx= Función de equilibrio y=f[x] )
%0 SIX NDUPN DROP
{ LAM x0 LAM x1 LAM V0 LAM y0 LAM fx
LAM y1 LAM C LAM V1 LAM DV LAM L0
LAM L1 }
BIND
LAM x1 LAM fx EVAL ' LAM y1 STO
LAM V0 %1 LAM y0 %- %* ' LAM C STO
LAM C %1 LAM y1 %- %/ ' LAM V1 STO
LAM V1 LAM V0 %- ' LAM DV STO
LAM DV %1 LAM x1 %- %* LAM x0 LAM x1 %- %/ ' LAM L0 STO
LAM L0 LAM DV %- ' LAM L1 STO
LAM L0 LAM L1 LAM V1 LAM y1
ABND
;
;
Y un programa de prueba, con los datos del ejemplo 4.2.1.1 es:
« 0. 1.84703223414E-3 200. .45
« 242.666666667 * »
ab1c3 4. →LIST { L0 L1 V1 y1 }
TORRES DE ABSORCIÓN POR ETAPAS - 135 -

→TAG LIST→ DROP

»
Con el cual se obtienen los resultados:
L0: 350.000000391
L1: 350.647657918
V1: 199.352342473
y1: .448213155485
Que prácticamente son los mismos valores del mencionado ejemplo.
Con este programa se puede resolver también el caso 1, así para los datos
del ejemplo 4.2.1.1., resolvemos el problema con el siguiente programa:
« 350 0 200 0.45 « 242.666666667 * »
« 242.666666667 / » 0 0 0 0
L0 x0 V0 y0 fx fy LOC L1 V1 y1
« « 'x1' STO
x0 x1 V0 y0 fx ab1c3 4 LIST
{ L0C L1 V1 y1 } STO
L0C L0 -
»
1E-8 y0 fy EVAL Refa DROP
L1 x1 V1 y1 4 LIST { L1 x1 V1 y1 } TAG
LIST DROP
»
»
Con lo que se obtienen los resultados:
L1: 350.647657375
x1: 1.84703223415E-3
V1: 199.352342474
y1: .448213155488
Que son prácticamente los mismos resultados calculados en el mencionado
ejemplo.
Este es un caso que al no ser iterativo permite resolver de manera itera-
tiva prácticamente todos los otros casos siguiendo un procedimiento similar
al presentado en el anterior ejemplo.