Absorcion de Gases MC Cabe
Absorcion de Gases MC Cabe
Absorcion de Gases MC Cabe
Absorción de gases
Este capítulo trata de las operaciones de transferencia de materia conocidas como ab-
sorción y eliminación o desorción de gases. En la absorción de gases un vapor soluble
se absorbe desde su mezcla con un gas inerte por medio de un líquido en el que el gas
(soluto) es más o menos soluble. Una aplicación primordial de la tecnología de absor-
ción es la eliminación de CO2 y H2S a partir del gas natural o del gas de síntesis por
absorción en disoluciones de sales de aminas o alcalinas. El lavado de amoniaco a partir
de una mezcla de amoniaco y aire por medio de agua líquida es otro ejemplo típico. El
soluto se recupera después del líquido por destilación, y el líquido absorbente se desecha
o se reutiliza. En ocasiones un soluto se remueve de un líquido poniendo este último en
contacto con un gas inerte; tal operación, que es inversa a la absorción, recibe el nombre
de desorción de gases o eliminación.
Entrada de líquido
Distribuidor de líquido
Sección
empacada
Sección
empacada
Entrada de gas
Salida de líquido
El gas que contiene el soluto, o gas rico, entra en el espacio de distribución situado
debajo del empaque y asciende a través de los intersticios del empaque en contracorrien-
te con el flujo del líquido. El empaque proporciona una gran área de contacto entre el
líquido y el gas, favoreciendo así un íntimo contacto entre las fases. El soluto contenido
en el gas rico es absorbido por el líquido fresco que entra en la torre, y el gas diluido o
agotado sale de la torre por la parte superior. El líquido se enriquece en soluto a medida
que desciende por la torre y el líquido concentrado, llamado licor concentrado, sale por
el fondo de la torre a través de la salida de líquido.
Los empaques de la torre se dividen en tres principales tipos: aquellos que son
cargados de forma aleatoria en la torre, los que son colocados a mano, y aquellos que se
conocen como empaques ordenados o estructurados. Los empaques aleatorios consisten
en unidades de 6 a 75 mm ( –14 o 3 in.) en su dimensión mayor; los empaques inferiores
a 25 mm se utilizan principalmente en columnas de laboratorio o de plantas piloto. Las
unidades de empaque ordenado son de tamaños comprendidos entre unidades de 50 a
200 mm (2 a 8 in.). Éstos se ocupan mucho menos que los emapaques aleatorios, pero
no se estudiarán aquí.
La mayoría de los empaques aleatorios de las torres se construyen con materiales
baratos e inertes, tales como arcilla, porcelana o diferentes plásticos. A veces se utilizan
anillos metálicos de pared delgada, de acero o aluminio. Se alcanzan altos espacios vacíos
(porosidad del lecho) y pasajes o pasos grandes para los fluidos haciendo las unidades
de empaque irregulares o huecas, de forma que se entrelazan para dar lugar a estructuras
abiertas con una porosidad de 60 a 90%.
En la figura 18.2 se ilustran empaques comunes de torres y sus características físicas
se presentan en la tabla 18.1. Las monturas de cerámicas Berl y los anillos Raschig son
los tipos de empaque más antiguos y no son muy usados en la actualidad, aunque repre-
sentaron una mejora importante respecto de las esferas de cerámica o la piedra triturada
que se introdujeron primero. Las monturas Intalox son similares a las monturas Berl, pero
la forma impide que las piezas queden demasiado juntas, y esto aumenta la porosidad del
lecho. Las monturas de Super Intalox tienen una pequeña variación con respecto al borde
escalopado; se encuentran disponibles en plástico o en cerámica. Los anillos Pall están
hechos de metal delgado con porciones de la pared inclinada hacia dentro, o de plástico
con ranuras en las paredes y costillas rígidas dentro. Los empaques Hy-pak metálicos
y Flexirings (no ilustrados) son similares en forma y funcionamiento a los anillos me-
tálicos Pall. Los lechos de anillos Pall tienen alrededor de 90% de fracción de huecos y
una ligera caída de presión que otros empaques de tamaño nominal parecido. El nuevo
IMTP de Norton (Instalox Metal Packing-empaque de torres metálico Intalox) tiene una
estructura muy abierta y una caída de presión más baja que la de los anillos Pall. Los
factores de empaque adicionales de caída de presión para algunos empaques comerciales
están dados por Robbins18 y, en unidades del SI, por Perry.16
Los empaques estructurados con orden geométrico han evolucionado desde los
empaques Stedman a finales de la década de 1930,26 pero se encontraron muy pocos
usos industriales hasta que se desarrollaron los empaques de Sulzer alrededor de 1965.24
Los primeros empaques estructurados se fabricaron de gasa alambre; los modelos más
recientes están hechos de láminas perforadas de metal corrugado, con láminas adyacentes
a) b) c)
d) e) f) g)
FIGURA 18.2
Empaques comunes en torres: a) anillos Raschig; b) anillo matálico Pall; c) anillo plástico Pall;
d) montura Berl; e) montura de cerámica Intalox; f) montura plástica Super Intalox;
g) montura metálica Intalox.
TABLA 18.1
Características de empaques para torres8,11b,21
Factores
Factoresdede
Tamaño, Densidad Área total, †
Porosidad empaque
empaque‡
‡
†
La densidad global y el área total están referidas a la unidad de volumen de la columna.
‡
El factor Fp es un factor de caída de presión y el factor ƒp es un coeficiente relativo de transferencia de materia. El factor ƒp se estudia en la página 632 en el
apartado “Funcionamiento de otros empaques”. Su uso se ilustra en el ejemplo 18.7.
§
Con base en datos de NH3–H2O; los demás factores se basan en datos de CO2–NaOH.
acomodadas de tal forma que el líquido se distribuye sobre sus superficies mientras que
el vapor fluye a través de los canales formados por los corrugados. Los canales están
colocados en ángulo de 45° respecto a la horizontal; el ángulo se alterna en dirección
en las capas sucesivas, como se muestra en la figura 18.3. Cada capa tiene unas cuantas
pulgadas de espesor. Varios empaques patentados difieren en el tamaño y distribución
de los corrugados y el tratamiento de las superficies de empaque.3 En general, las
corrugaciones triangulares tienen de 25 a 40 mm a lo largo de la base, de 17 a 25 mm
de lado, y de 10 a 15 mm de altura. Los intervalos de porosidad van de 0.93 a 0.97, y el
área específica de superficie de 60 a 76 ft2/ft3 (200 a 250 m2/m3). El empaque Sulzer BX,
fabricado de alambre de metal, tiene un área específica de superficie de 152 ft2/ft3 (500
m2/m3) con una porosidad de 0.90.
Dirección principal
de flujo del vapor
Dirección principal
de flujo del líquido
FIGURA 18.3
Trazado esquemático de un empaque estructurado.
La figura 18.4 muestra datos típicos de la caída de presión en una torre empacada. La caída
de presión por unidad de longitud (o profundidad) del empaque se debe a la fricción del
fluido; se grafica en coordenadas logarítmicas frente a la velocidad de flujo del gas Gy,
expresada en masa de gas por hora y por unidad de área de la sección transversal, conside-
rando que la torre está vacía. Por lo tanto, Gy está relacionada con la velocidad superficial
del gas por medio de la ecuación Gy = u0 ρy, donde ρy es la densidad del gas. Cuando el
empaque está seco, la línea que se obtiene es recta y tiene una pendiente del orden de 1.8.
Por consiguiente, la caída de presión aumenta con la velocidad elevada a una potencia
de 1.8. Si el empaque está irrigado con un flujo constante de líquido, la relación entre
la caída de presión y la velocidad de flujo del gas sigue al principio una línea paralela a la
del empaque seco. La caída de presión es mayor que en el empaque seco, debido a que
el líquido en la torre reduce el espacio disponible para el flujo de gas. Sin embargo, la
fracción de huecos no varía con el flujo de gas. Para velocidades moderadas del gas,
la línea para el empaque irrigado tiene una pendiente cada vez más pronunciada, debido
a que el gas impide el flujo descendente del líquido de forma que aumenta la retención de
éste con la velocidad de flujo del gas. El punto en el que la retención de líquido comienza
a aumentar, hecho que se aprecia por un cambio de la pendiente de la línea de la caída de
presión, recibe el nombre de punto de carga. Sin embargo, como se aprecia en la figura
18.4, no es fácil obtener un valor exacto para el punto de carga.
Al aumentar todavía más la velocidad del gas, la caída de presión se incrementa aún
más rápido, y las líneas se hacen casi verticales cuando la caída de presión es del orden de
2 a 3 in. de agua por pie de empaque (150 a 250 mm de agua por metro). En determinadas
regiones de la columna, el líquido se transforma en una fase continua y se dice que la co-
lumna está inundada. Temporalmente se utilizan flujos de gas más elevados, pero el líquido
se acumula con mayor rapidez, y la columna completa puede llenarse con líquido.
2.0
00
00
00
000
o
∆P, pulgadas de agua /ft de empaque
00
00
0
0
Sec
0
12 0
9 00
6 00
15 0
20 0
3 00
60
00
25
=3
=3
G =
G =
G =
G =
G =
G =
G =
x
x
x
x
x
x
x
G
1.0
0.8
0.6
Fp = 92
0.2
Velocidad del líquido, lb/ft2 − h,
como párametro FIGURA 18.4
Caída de presión en una torre empacada
0.1 para el sistema aire-agua con monturas
100 200 300 400 600 1 000 2 000 5 000
Intalox de 1 in. (1 000 lb/ft2 · h = 1.356
Velocidad másica del aire, Gy , lb/ft 2 − h kg/m2 · s; 1 in. H2O/ft = 817 Pa/m).
Es evidente que la velocidad del gas en la columna empacada en operación debe ser
inferior a la velocidad de inundación. Sin embargo, a medida que se aproxima la inunda-
ción, la mayor parte o toda la superficie de empaque se humedece, aumentando el área de
contacto entre el gas y el líquido. El diseñador debe escoger una velocidad suficientemente
distante de la velocidad de inundación para garantizar una operación segura, pero no tan
baja que se requiera de una columna mucho más grande. Bajar la velocidad del diseño hace
que se incremente el diámetro de la torre sin mucho cambio en la altura requerida, a partir
de que bajas velocidades de gas y líquido llevan a una reducción proporcional cercana a
la velocidad de la transferencia de masa. Uno de los beneficios de la baja velocidad del
gas es el decremento de la caída de presión, pero el costo de la energía consumida no es
por lo común un factor importante en la optimización del diseño. En algunas ocasiones, la
velocidad del gas es escogida como la mitad de la velocidad de inundación que se predice
a partir de una correlación generalizada. Esto quizá parezca muy conservador, pero hay
una dispersión considerable en los datos publicados para la velocidad de inundación, y las
correlaciones generalizadas no son muy exactas. Se emplea una aproximación más cercana
a la velocidad de inundación si se dispone de los datos detallados de funcionamiento para
la selección del empaque. También es posible diseñar torres empacadas con base en una
caída de presión definida por unidad de altura del empaque.
La velocidad de inundación depende en forma importante del tipo y tamaño del em-
paque y la velocidad másica de líquido. La figura 18.5 muestra información para monturas
Intalox tomados de la figura 18.4 y curvas similares para otros tamaños. Se supone que la
5 000
2 000
Velocidad de inundación Gy, lb/ft 2 − h
1 000
500
1
12 in.
2 in.
1
Tamaño del empaque 2
in.
1 in.
200
FIGURA 18.5
Velocidades de inundación
en monturas cerámicas
100 Intalox en un sistema aire-
1 000 2 000 5 000 10 000 20 000 50 000 agua. (1 000 lb/ft2 · h =
2
Velocidad másica del líquido Gx, lb/ft − h 1.356 kg/m2 · s).
velocidad de inundación ocurre a una caída de presión de 2.0 in. de H2O/ft de empaque,
puesto que las curvas de la caída de presión son verticales o cercanas a este punto. Para
bajas velocidades de líquido, la velocidad de inundación varía con la velocidad del líquido
elevado a la potencia de –0.2 a –0.3 y con el tamaño del empaque elevado a la potencia
de 0.6 a 0.7. Los efectos de la velocidad del líquido y el tamaño del empaque se vuelven
más pronunciados en altas velocidades másicas del líquido.
Se han propuesto varias correlaciones generalizadas para la caída de presión y la
velocidad de inundación en columnas empacadas. Muchas de ellas utilizan una gráfica log-
log con (Gx /Gy)(ρy / ρx)0.5 en la abscisa y una función que contiene Gy2 en la ordenada. Por
lo regular, la relación de flujo Gx /Gy se toma a partir del equilibrio y las consideraciones
económicas que se explicarán después en este capítulo, y Gy se determina directamente,
mientras que se requiere de la solución por prueba y error si Gy y Gx están en ejes sepa-
rados, como en la figura 18.5. Las características del empaque están dadas por un factor
de empaque Fp, el cual disminuye si se aumenta el tamaño del empaque o la fracción de
vacío (o fracción hueca). No es posible predecir los factores de empaque a partir de la
teoría utilizando la ecuación de Ergun [ecuación (7.22)] a causa de las formas complejas,
y por eso se determinan en forma empírica. Desafortunadamente, no hay correlaciones
únicas para caída de presión que den un buen ajuste para todos los empaques, y los
valores de Fp basados en el ajuste de datos para bajas caídas de presión tal vez difieran
significativamente de los valores obtenidos a partir del ajuste de datos para altas caídas
de presión o por ajuste de datos de la velocidad de inundación.
Una correlación muy utilizada para estimar caídas de presión en empaques colocados
de manera aleatoria se presenta en la figura 18.6 donde Gx y Gy están en lb/ft2 · s, µx está
en cP, ρx y ρy están en lb/ft3, y gc es 32.174 lb · ft/lbf · s2. Las primeras versiones de esta
0.60
0.40
0.040 0.25
0.020 0.10
0.010
0.006 0.05
0.004
0.002
FIGURA 18.6
0.001
Correlación generalizada para
0.01 0.02 0.04 0.06 0.1 0.2 0.4 0.6 1.0 2.0 4.0 6.0 10.0
la caída de presión en columnas
Gx ry
empacadas. (1 in. H2O/ft = 817
Gy rx − r y
Pa/m) (Según Eckert.3)
correlación incluían una línea de inundación por encima de la línea para ∆P = 1.5 in.
H2O/ft de empaque, pero estudios recientes estiman la inundación a caídas de presión de
sólo 0.7 a 1.5 in. H2O/ft para empaques de 2 o 3 in. Una ecuación empírica para el límite
de la caída de presión es9
La ecuación (18.1) se utiliza para factores de empaque desde 10 hasta 60. Para valores
mayores de Fp, la caída de presión en la inundación se toma como 2.0 in. H2O/ft.
Una correlación alternativa para la caída de presión en columnas empacadas fue pro-
puesta por Strigle27 y se presenta en la figura 18.7. La abscisa es esencialmente la misma que
para la figura 18.6, pero la ordenada incluye el factor de capacidad Cs = u0√ρy/(ρx – ρy), donde
u0 es la velocidad superficial en pies por segundo. La viscosidad cinemática del líquido ν
está en centistokes. La gráfica semilogarítmica permite que la interpolación sea más fácil
que la gráfica log-log, aunque ambas estén basadas en los mismos datos.
EJEMPLO 18.1 Se va a construir una torre empacada con monturas Intalox de cerámica de 1 in.
(25.4 mm) para tratar una alimentación de 25 000 ft3 (708 m3) de un gas por hora. El contenido
de amoniaco del gas que ingresa es de 2% en volumen. Agua libre de amoniaco se utiliza como
absorbente. La temperatura del gas y la del agua que ingresan es de 68 ºF (20 ºC); la presión es
de 1 atm. La relación entre flujo del líquido y flujo de gas es de 1.25 lb de líquido por libra de
gas. a) Si la caída de presión de diseño es de 5 in. de H2O por pie de empaque, ¿cuál debería ser
la masa velocidad del gas y el diámetro de la torre?; b) estime la relación de la velocidad del
2.4
∆P = (in. H2O/ft)
∆P
=1
2.0 1.0 .5
1.6 0.50
CsFp0.5n0.05
1.2 0.25
0.8 0.10
0.05
0.4
gas con la velocidad de inundación, usando una correlación generalizada y además usando
los datos específicos para las monturas Intalox.
Solución El peso molecular promedio del gas que ingresa es (29 0.98) (17 0.02) =
28.76. Entonces
En este caso, el valor de Gy se puede verificar usando los datos para las monturas Intalox de
1 in. de la figura 18.4. La interpolación entre las curvas para empaques secos y Gx = 3 000
lb/ft2 h da Gy 1 000 lb/ft h. Este valor se utiliza para el diseño.
Flujo total del gas: 25 000 0.0746 = 1 865 lb/h
Área transversal de la torre:
b) La correlación generalizada de la figura 18.7 no incluye una línea para inundación, pero
la caída de presión en la inundación es de alrededor de 2.0 in. de H2O/ft y la extrapolación de
este valor indica una ordenada de 1.95 comparada con el valor de diseño de 1.38. Por lo tanto,
Esta relación se aplica cuando Gx / Gy se mantiene constante. La relación es más baja con
un Gx más bajo. Si Gx se mantiene a 1 400 lb/ft2 h, la figura 18.5 indica inundación en Gx =
1.650 lb/ft2 h. Entonces,
Empaques estructurados
EJEMPLO 18.2 Sustituya el empaque Gempack 230 A2T, fabricado por Glitsch, por las
monturas Intalox de 1 in. del ejemplo 18.1. ¿Qué incremento se espera en la masa velocidad
del gas?
0.20
Montz B1-200
0.15 Gempak 230 A2T
Montz B1-100
u0,f r − r , m/s
0.08
ry
x
0.04
0.02
.01 .02 .03 .04 .06 .08 .10 .20 .40 .60
FIGURA 18.8
Gx ry
Velocidades de inundación
Gy rx
en empaques estructurados.4
FUNDAMENTOS DE LA ABSORCIÓN
Balances de materia
En una planta de contacto diferencial, tal como la torre empacada de absorción que se
ilustra en la figura 18.9, las variaciones de composición son continuas de un extremo a
otro del equipo. Los balances de materia para la porción de la columna por encima de
una sección arbitraria, tal como se representa con la línea de trazos discontinuos en la
figura 18.9, son los siguientes:
Balance total de material: La + V = L + Va (18.2)
L V y − La x a
y= x+ a a (18.6)
V V
Es posible graficar la línea de operación sobre una gráfica aritmética a lo largo de la curva
de equilibrio, tal como se muestra en la figura 18.10. La línea de operación debe estar
por arriba de la línea de equilibrio para que la absorción tenga lugar, puesto que esto
proporciona una fuerza impulsora positiva y – y* para la absorción.
En la ecuación (18.6), x y y representan las composiciones globales del líquido y el
gas, respectivamente, en contacto entre sí en cualquier sección de la columna. Se supone
que las composiciones para una altura determinada son independientes de la posición en
el empaque. La absorción de un componente soluble desde una mezcla gaseosa da lugar
Va FIGURA 18.9
Diagrama del balance de materia para una columna
ya Superficie empacada.
de control
La
xa
L V
x y
Vb
yb Lb
xb
a una disminución de la velocidad total del gas V a medida que el gas pasa a través de
la columna, mientras que el flujo del líquido L aumenta. Estos cambios provocan una
ligera curvatura de la línea de equilibrio, tal como se observa en la figura 18.10. En el
caso de las mezclas diluidas, que contienen menos de 10% de gas soluble, el efecto de las
variaciones del flujo total generalmente se ignora y el diseño se basa en las velocidades
de flujo promedio.
FIGURA 18.10
Relación gas-líquido limitante
(o mínima).
b b1
yb
in
)m
V
y (L/
n,
n
ió
ió
ac
ac
er
er
op
op
de yb*
de
ea
Lín
ea
rio
Lín
u ilib
q
ee
rvad
ya
a Cu
xa xb xb*
x
del líquido tiene que ser mayor que este mínimo para conseguir el cambio especificado en
la composición del gas.
En una columna de absorción en contracorriente, la relación L/V es importante desde el
punto de vista económico. La fuerza impulsora para la transferencia de materia es y – y*, que
es proporcional a la distancia vertical entre la línea de equilibrio y la línea de operación en un
diagrama como el de la figura 18.10. Al aumentar L/V, se incrementa la fuerza impulsora
en todas las partes de la columna excepto en la parte superior, y la columna de absorción
no necesita ser tan alta. Sin embargo, al utilizar más cantidad de líquido se obtiene una
disolución líquida más diluida, de la que es más díficil recuperar el soluto por desorción
o eliminación. El costo de la energía que se requiere para la eliminación suele representar
la parte más importante del costo total de la operación de absorción-eliminación. El flujo
óptimo de líquido para la absorción se obtiene equilibrando los costos de operación para
ambas unidades en función de los costos fijos del equipo. En general, la velocidad del
líquido para un absorbedor deberá estar comprendido en el intervalo de 1.1 a 1.5 veces la
velocidad mínima, a menos que el líquido se deseche y no se regenere.
Las condiciones en la parte superior del absorbedor son con frecuencia variables
de diseño que también han de establecerse teniendo en cuenta los costos del equipo y
los de operación. Por ejemplo, si unas especificaciones tentativas corresponden a 98%
de recuperación de un producto a partir de una corriente gaseosa, el diseñador debería
calcular qué altura adicional de columna se requeriría para lograr una recuperación de
99%. Si el valor del producto recuperado excede los costos adicionales, la recuperación
óptima debe ser al menos de 99% y habrá que repetir los cálculos para recuperaciones
aun mayores. Si el soluto no recuperado es contaminante, su concentración en el gas de
purga se determinará por omisiones estándar, y el porcentaje de recuperación requerido
excederá al valor óptimo basado en el valor del producto y los costos de operación.
Velocidad de absorción
La velocidad de absorción se expresa de cuatro formas diferentes utilizando coeficientes
individuales o globales basados en las fases gaseosa o líquida. Para la mayor parte de los
cálculos se utilizan coeficientes volumétricos debido a que es más difícil determinar
los coeficientes por unidad de área, ya que el propósito del cálculo del diseño consiste
por lo general en determinar el volumen total del absorbedor. En el tratamiento que
sigue se omiten, para simplificar, los factores de corrección para la difusión en una sola
dirección, y se desprecian las variaciones de las velocidades de flujo del gas y el líquido.
Las ecuaciones son estrictamente válidas sólo para gases diluidos, pero es posible em-
plearlas con poco error para mezclas con un contenido de hasta 10% de soluto. El caso
de la absorción a partir de gases ricos (con altas concentraciones) se trata más adelante
como un caso especial.
La velocidad de absorción por unidad de volumen de la columna empacada se de-
termina por cualquiera de las siguientes ecuaciones, donde y y x se refieren a la fracción
mol del componente que se absorbe:
r = k y a( y − yi ) (18.7)
r = k x a( xi − x ) (18.8)
(
r = K y a y − y∗ ) (18.9)
(
r = Kx a x∗ − x ) (18.10)
y − yi k x a
= (18.11)
xi − x k y a
Por lo tanto, una recta trazada desde la línea de operación con una pendiente –kx a/kya
intersectará a la línea de equilibrio en el punto (yi, xi), tal como se indica en la figura 18.11.
Por lo general no se necesitan conocer las composiciones en la superficie de contacto,
pero estos valores se utilizan para los cálculos cuando intervienen gases ricos o cuando
la línea de equilibrio presenta una curvatura pronunciada.
Las fuerzas impulsoras globales se determinan fácilmente como líneas verticales
u horizontales en el diagrama yx. Los coeficientes globales se obtienen a partir de ky a
y kxa utilizando la pendiente local de la curva de equilibrio m, tal como se indicó en el
capítulo 17 [ecuación (17.61)]:
1 1 m
= + (18.12)
K y a ky a kx a
1 1 1
= + (18.13)
K x a k x a mk y a
En la ecuación (18.12), los términos 1/(ky a) y m/(kxa) son las resistencias a la trans-
ferencia de masa en la película gaseosa y la película líquida, respectivamente. Cuando
los coeficientes kya y kxa son del mismo orden de magnitud,† y m es mucho mayor que
1.0, se dice que la resistencia de la película líquida está controlada. Esto quiere decir
que cualquier cambio en kxa tiene un cercano efecto proporcional tanto en Kya y Kxa
como en la velocidad de absorción, mientras que un cambio en kya sólo tiene un pequeño
efecto. Por ejemplo, el coeficiente de la ley de Henry para CO2 en agua a 20 °C es de
1 430 atm/fracción mol, lo cual corresponde a m = 1 430 para la absorción a 1 atm y
m = 143 para la absorción a 10 atm. Bajo estas condiciones, es evidente que la absorción
de CO2 en agua está controlada por la película líquida. Incrementar la velocidad del gas
aumentará kya pero tendrá un efecto despreciable en Kya. Al aumentar la velocidad del
líquido aumentará el área de la superficie de contacto a y probablemente también se
incrementará kx, lo cual lleva a un incremento en kxa y Kya.
Línea de
kx a operación
Pendiente = −
ky a
Línea de
y equilibrio
y
yi x*
xi
y*
FIGURA 18.11
x Identificación de las composiciones en la super-
x ficie de contacto.
†
kx a es por lo general de 3 a 20 veces kya en columnas empacadas.
Cuando la solubilidad del gas es muy alta, tal como ocurre con HCl en agua, m es muy
pequeña y la resistencia de la película gaseosa controla la velocidad de absorción. Con
gases de solubilidad intermedia, ambas resistencias son importantes, aunque el término
resistencia controlada se utiliza a veces para resistencias muy grandes. La absorción de
NH3 en agua con frecuencia se cita como un ejemplo del control de la película gaseosa,
puesto que la película gaseosa tiene de 80 a 90% de la resistencia total.
(
−Vdy = K y a y − y∗ SdZ ) (18.14)
Salida de gas
ya
dZ ZT
y x
Entrada de gas
Salida de líquido
yb
xb
FIGURA 18.13
yb Relación entre el número de unidades de
yb Línea de operación
Línea de operación transferencia (NTU) y el número de los
platos teóricos (NTP): a) NTU = NTP;
y y
b) NTU > NTP.
ya Línea de equilibrio
ya Línea de equilibrio
xb xb
x x
a) b)
L/S dx
Película líquida: Hx =
kx a
Nx = ∫x i −x
(18.22)
V /S dy
Gas global: HOy =
Kya
NOy = ∫ y−y ∗ (18.23)
L/S dx
Líquido global: HOx =
Kx a
NOx = ∫x ∗
−x
(18.24)
kga y Kga son por lo común mol/ft3 · h · atm. De manera análoga, los coeficientes de
película líquida se expresan como kLa y KLa, donde la fuerza impulsora es una diferencia
de concentración volumétrica; por lo tanto, kL es igual a kc definido por la ecuación (17.39).
Por consiguiente, kLa y KLa son iguales a kxa/ρM y Kxa/ρM, respectivamente, donde ρM es
la densidad molar del líquido. Las unidades de kLa y KLa son por lo general mol/ft3 · h ·
(mol/ft3) o h–1.
Si Gy /M o GM se sustituye por V/S en las ecuaciones (18.21) y (18.23), y Gx /M por
L/S en las ecuaciones (18.22) y (18.24), las ecuaciones para la altura de una unidad de
transferencia se escriben (puesto que MρM = ρx, la densidad del líquido)
GM GM
Hy = y HOy = (18.25)
kg aP K g aP
Gx / ρ x Gx / ρ x
Hx = y HOx = (18.26)
kL a KLa
Los términos HG, HL, NG y NL con frecuencia aparecen en la bibliografía en vez de Hy,
Hx, Ny y Nx, así como los correspondientes términos para los valores globales, pero en este
caso los distintos subíndices no significan ninguna diferencia en unidades o magnitud.
Si un diseño está basado en NOy, el valor de HOy se calcula a partir de Ky a, o bien
a partir de los valores de Hy y Hx, tal como se verá a continuación. Comenzando con
la ecuación para la resistencia global, ecuación (18.12), cada término se multiplica por
GM, y el último término se multiplica por LM /LM, donde LM = L/S = Gx /M, la velocidad
másica molar del líquido:
GM G mGM LM
= M +
K y a ky a k x a LM (18.27)
GM
HOy = H y + m Hx (18.28)
LM
LM
HOx = H x + Hy (18.29)
mGM
EJEMPLO 18.3 Una corriente gaseosa que contiene 3.0% de A se pasa a través de una columna
empacada para remover 99% de A por absorción en agua. El absorbedor operará a 25 °C y
1 atm, y las velocidades de gas y líquido están a 20 mol/h · ft2 y 100 mol/h · ft2, respectivamente.
Los coeficientes de transferencia de masa y los datos de equilibrio se dan a continuación:
y* = 3.1x a 25 °C
kxa = 60 mol/h · ft3 · unidad de fracción molar
ky a = 15 mol/h · ft3 · unidad de fracción molar
a) Encuentre NOy, HOy y ZT, asumiendo que se trata de una operación isotérmica y despreciando
los cambios en las velocidades de flujo de gas y líquido. ¿Qué porcentaje de la resistencia
total está en la fase gaseosa?
b) Calcule ZT, utilizando NOx y HOx.
Solución
a) Suponga que xa = 0. Puesto que GM ∆y = LM ∆x,
20 × 0.03 × 0.99
xb = = 0.00594
100
yb∗ = 3.1 × 0.0054 = 0.01841
En el fondo de la columna,
0.01159 − 0.0003
Entonces ∆y L = = 0.00309
ln(0.01159 / 0.0003)
∆y 0.03 × 0.99
NOy = = = 9.61
∆y L 0.00309
1 1 3.1
= + = 0.11833 K y a = 8.45
K y a 15 60
20
HOy = = 2.37 ft
8.45
ZT = 2.37 × 9.61 = 22.7 ft
0.00374 − 0.000097
∆x L = = 9.96 × 10 −4
[ (
ln 0.00374 / 9.677 × 10 −5 )]
0.00594
NOx = = 5.96
9.96 × 10 −4
1 1 1
= + = 0.03817 K x a = 26.2
K x a 60 3.1 × 15
100
HOx = = 3.817 ft
26.2
ZT = 5.96 × 3.817 = 22.7 ft
Obsérvese que aunque NOy y NOx son bastante diferentes, los valores de HOy y HOx son
también distintos, de modo que los valores de ZT son los mismos para los dos métodos
de cálculo.
Efecto de la presión
Las columnas de absorción con frecuencia son operadas bajo presión para incrementar
la capacidad y elevar las velocidades de transferencia de masa. La presión parcial de
equilibrio del soluto depende sólo de la composición del líquido y de la temperatura, así
que la fracción molar de equilibrio en el gas varía inversamente con la presión total
pA
yA = (18.30)
P
Si las velocidades de gas y líquido se mantienen constantes de modo que la línea de
operación no cambie hacia una presión mayor, se incrementa la fracción molar de la fuerza
impulsora, tal como se muestra en la figura 18.14, y se reduce el número de unidades
de transferencia. A mayor presión, la velocidad mínima del líquido es más pequeña, así
que la línea de operación podría cambiar para dar un producto más rico, como se indica
con la línea de trazos discontinuos en la figura 18.14, y alrededor del mismo número de
unidades de transferencia que antes.
El coeficiente global de transferencia de masa κya se incrementa con la presión porque
la resistencia de la película del líquido m/kxa decrece [véase ecuación (18.12)]. El coefi-
ciente de la película gaseosa κya no cambia mucho con la presión; la difusividad varía
inversamente con la presión, pero el gradiente de concentración para una determinada
∆y es proporcional a la presión. De este modo, HOy decaerá si la resistencia del líquido
es significativa, pero no cambiará si la resistencia de la película gaseosa controla.
1
=
P
b rio
u ili
eq
1
2
P=
de P =
=2
a ión oP
n
e ri
ció
Lín c ilib
e ra qu
era
y
op d ee
de
op
Lí nea
ea
de
Lín
ea
L ín
xb x′b
x
FIGURA 18.14
Efecto de la presión en la absorción.
ido
l líqu Tb
ra de
eratu uid
o
Temp Ti íq
el l
Ta rad Tb
yb atu
per
rio Tem
lib Ta
n
i
ió
u yb
eq Línea, eq.
ac
Línea, eq.
T, °C
er
y de Línea de operación
op
Tb o Tmáx
T, °C
ea y r i
de
Lín Línea, eq ib
Línea, eq. uil Línea, eq.
eq
ea
T +T
Ti Ti = a b e
Lín
d Tb
2
í nea
L
Línea, eq. Línea, eq.
ya Ta ya Ta
xa xi xb xa xb
x x
a) b)
FIGURA 18.15
Perfiles de temperatura y líneas de equilibrio para la absorción adiabática: a) sin evaporación del
solvente: b) evaporación del solvente apreciable o gas de alimentación frío.
EJEMPLO 18.4 Una corriente de gas con 6.0% de NH3 (base seca) y una velociad de flujo
de 4 500 SCFM (Standard Cubic foot per minute) (ft3/min a 0 °C, 1 atm) se depurará con agua
para bajar la concentración a 0.02%. El absorbedor operará a la presión atmosférica con tempe-
raturas de entrada de 20 y 25 °C para el gas y el líquido, respectivamente. El gas está saturado
con vapor de agua a la temperatura de entrada y se supone que saldrá como gas saturado a
25 °C. Calcule el valor de NOy si la velocidad del líquido es 1.25 veces el mínimo.
Bases 100 g mol de gas seco de entrada, contienen 94 moles de aire y 6 moles de NH3. La
salida de gas contiene 94 moles de aire.
Las moles de amoniaco en el gas que sale, a partir de que ya = 0.0002, son
0.0002
94 = 0.0188 mol NH 3
0.9998
Efectos del calor El calor de absorción es 5.98 × 8 310 = 49 690 cal. Llamemos a esto
Qa. Entonces
17.5
100 × = 2.36 mol
742.5
En la salida de gas es
23.7
94.02 × = 3.03 mol
736.3
La cantidad de agua vaporizada es por lo tanto 3.03 – 2.36 = 0.67 mol. A partir de que el calor
de evaporización es ∆Hυ = 583 cal/g,
Al resolver la ecuación (18.31) para Qsx, el cambio de calor sensible en el líquido, da como
resultado
La temperatura de salida del líquido Tb se encuentra por tanteo. Suponga que para el sol-
vente Cp = 18 cal/g mol · °C; Tb = 40 °C y xmáx = 0.031, como se estimó de la líneas de solubilidad
en equilibrio en la figura 18.16. Entonces el total de moles de líquido que sale Lb son
5.98
Lb = = 192.9 mol
0.031
Partiendo de que Ta = 25 °C
Para 1.25 veces el mínimo de la velocidad del agua, La = 1.25 × 175 = 219 moles, y Lb = 219
+ 6 = 225 moles. Entonces, el aumento de la temperatura del líquido es
39 260
Tb − 25 = = 9.7 oC
225 × 18
Absorción multicomponente
Cuando se absorbe más de un soluto de una mezcla de gas se necesita separar para cada
soluto las líneas de equilibrio y de operación, pero la pendiente de la línea de operación,
que es L/V, es la misma para todos los solutos. En la figura 18.17 se muestra un diagrama
yx típico para la absorción de dos solutos. En este ejemplo, B es un componente menor
0.08
40°C
0.07
0.06
ón
aci
35°C
30°C
per
0.05
eo
io
ad
libr
Líne
qui
yNH3
0.04
ee
ead
Lín
0.03
25°C
0.02
0.01
FIGURA 18.16
0 0.01 0.02 0.03 0.04 0.05 El diagrama yx para el
x NH3 ejemplo 18.4.
A
ra
pa
ón
ci
y*A
a
y A o yB
er
op
de
ea
yB,b
Lín
yB,a
yA,a
xB,b xA,a
xA o xB
del gas, y la velocidad del líquido se escogió de manera que permitiera la eliminación
de 95% de A con una altura empacada razonable. La pendiente de la línea de operación
es alrededor de 1.5 veces la pendiente de la línea de equilibrio para A, y NOy ≈ 5.5. La
línea de operación para B es paralela a la de A, y puesto que la línea de equilibrio para
B tiene una pendiente mayor que L/V, hay un acercamiento en el fondo de la columna, y
sólo es posible absorber una pequeña fracción de B. La línea de operación para B debe
trazarse a fin de dar un número de unidades de transferencia correcto para B, el cual es
por lo general igual que NOy para A. Sin embargo, en este ejemplo, xB,b es prácticamente el
mismo que xB* , el valor de equilibrio para yB,b y la eliminación fraccional de B se calcula
directamente a partir de un balance de materia:
( ) (
V y B, b − y B, a = L x ∗B, b − x B, a ) (18.32)
Desorción o eliminación
En muchos casos, un soluto que se absorbe a partir de una mezcla gaseosa, se desorbe a
partir del líquido para recuperar el soluto en una forma más concentrada y regenerar la
solución absorbente. Con el fin de favorecer las condiciones de desorción se incrementa
la temperatura o se reduce la presión total, o bien, se realizan ambas modificaciones de
manera simultánea. Si la absorción se realiza a alta presión, a veces es posible recuperar
una fracción grande de soluto simplemente efectuando una evaporación espontánea (flaw-
shing) a la presión atmosférica. Sin embargo, para una separación casi completa del soluto,
generalmente se requieren varias etapas y la desorción o eliminación se lleva a cabo en
una columna con flujos de gas y líquido en contracorriente. Como medio de eliminación
se puede utilizar un gas inerte o vapor de agua, pero la recuperación del soluto resulta más
fácil cuando se utiliza vapor de agua, ya que es posible condensar este componente.
En la figura 18.18 se muestran las líneas de operación y equilibrio para una elimi-
nación o desorción con vapor de agua. Cuando xa y xb se especifican, existe una relación
mínima de vapor a líquido correspondiente al punto en el que la línea de operación toca la
línea de equilibrio. Este punto se presenta en la parte intermedia de la línea de operación
si la línea de equilibrio presenta una curvatura hacia arriba, como en la figura 18.18, o
bien se presenta en la parte superior de la columna, para (ya, xa). Con el fin de simplificar,
la línea de operación se considera recta, aunque de una forma general tendrá una ligera
curvatura debido a las variaciones de las velocidades de vapor y líquido.
En un proceso global de absorción y desorción, el costo del vapor de agua es con fre-
cuencia el más importante; el proceso se diseña para utilizar la menor cantidad posible de
vapor. En la columna de desorción se opera con un flujo de vapor cercano al mínimo,
dejando que salga algo de soluto con la disolución agotada en vez de tratar de obtener
una recuperación completa. Cuando la línea de equilibrio presenta una curvatura hacia
arriba, como en la figura 18.18, el flujo mínimo de vapor de agua es mucho más alto a
medida que xb tiende hacia cero.
La altura de la columna de desorción se calcula a partir del número de unidades de
transferencia y de la altura de una unidad de transferencia, utilizando las mismas ecua-
ciones que para la absorción. Con frecuencia, la atención se fija en la concentración de
la fase líquida, y se utilizan NOx y HOx:
dx
ZT = HOx NOx = HOx ∫x
−x ∗
(18.33)
Línea de equilibrio,
P = 1 atm
ya, máx,
y Línea de operación,
ya para Vmín
Línea de operación,
Pendiente = L/V
FIGURA 18.18
xb xa Líneas de operación para una
x columna de desorción.
La desorción con aire se utiliza en algunos casos para eliminar pequeñas cantidades
de gases como amoniaco o solventes orgánicos contenidos en agua. Si no es necesario
recuperar el soluto en forma concentrada, la cantidad óptima de aire usado puede ser
mucho mayor que el mínimo, debido a que no resulta muy costoso el uso de más aire,
y la altura de la columna se reduce considerablemente. El siguiente ejemplo muestra el
efecto de la velocidad del aire en una operación de desorción.
EJEMPLO 18.5 Agua que contiene 6 ppm de tricloroetileno (TCE) va a ser purificada por
desorción con aire a 20 °C. El producto debe contener menos de 4.5 ppb de TCE para estar
dentro de los estándares de emisiones. Calcule la velocidad mínima de aire en metros cúbi-
cos estándar de aire por metro cúbico de agua y el número de unidades de transferencia si la
velocidad del aire es de 1.5 a 5 veces el valor mínimo.
atm ⋅ m 3 1 10 6 mol H 2 O
m = 0.0075 × × = 417
mol 1 atm 18 m3
Debido a que este valor de m es grande, la desorción está controlada en la fase líquida. A la
velocidad mínima del aire, la salida del gas estará en equilibrio con la solución a la entrada.
El peso molecular de TCE es 131.4, y
6 × 10 −6 mol TCE 18 g
xa = × = 8.22 × 10 −7
131.4 g H 2 O mol H 2 O
( )
ya = 417 8.22 × 10 −7 = 3.43 × 10 −4
VTCE =
(
10 6 6 × 10 −6 − 4.5 × 10 −9 )
131.4
= 4.56 × 10 −2 mol
4.56 × 10 −2
V= = 132.9 mol
3.43 × 10 −4
Debido a que 1 g mol = 0.0224 m3 std y puesto que la variación en la velocidad del flujo del
gas es muy pequeña,
La densidad de aire en condiciones estándar es 1.295 kg/m3, así que la velocidad mínima en
una base másica es
Gy 2.98 × 1.295
G = = 3.86 × 10 −3 kg aire/kg agua
x mín 1 000
3.43 × 10 −4
ya = = 2.29 × 10 −4
1.5
2.29 × 10 −4
x a∗ = = 5.49 × 10 −7
417
131.4
Ca∗ = 5.49 × 10 −7 × = 4.01 × 10 −6 g/g = 4.01 ppm
18
Ca − Ca∗ = ∆Ca = 6.0 − 4.01 = 1.99
En el fondo,
(C − C ∗ ) L =
1.99 − 0.0045
ln(1.99/0.0045)
= 0.3259 ppm
dC Ca − Cb
NOx =
∫ C − C = (C − C )
∗ ∗
L
6 − 0.0045
= = 18.4
0.3259
Cálculos similares para otros múltiplos de la velocidad mínima de flujo dan los siguientes
valores. La altura empacada se basa en el valor estimado de HOx = 3 ft; esto es un poco mayor
que los valores reportados para anillos plásticos Pall de 1 in.
flujo de líquido, a la vez que deberá incluirse el factor de corrección para la difusión en un
solo sentido. Por otra parte, los coeficientes de transferencia de masa no serán constantes
debido a las variaciones de las velocidades de flujo y de un gradiente de temperatura apre-
ciable en la columna, lo que dará lugar a un cambio en la línea de equilibrio.
La cantidad de soluto absorbido en una altura diferencial dZ es d(Vy), ya que tanto
V como y disminuyen a medida que el gas pasa a través de la torre.
Si solamente se transfiere A, dNA es igual a dV, de forma que la ecuación 18.34 se trans-
forma en
o Vdy
dN A = (18.36)
1− y
El efecto de la difusión en un solo sentido en la película gaseosa se traduce en un
aumento en la velocidad de transferencia de masa por el factor 1/(1 – y)L, tal como indica
la ecuación (17.49), de forma que el coeficiente global efectivo K′y a es algo mayor que
el valor normal de Kya:
1 (1 − y) L m (18.37)
= +
K y′ a ky a kx a
En este tratamiento se ha despreciado el efecto de la difusión en un solo sentido en la
película líquida. La ecuación básica de transferencia de masa es, por lo tanto,
dN A =
Vdy
1− y
= K y′ aSdZ y − y ∗ ( ) (18.38)
La altura de la columna se obtiene por integración gráfica, teniendo en cuenta las varia-
ciones de V, 1 – y, y – y* y K′y a:
1 b Vdy
ZT =
S ∫ (1 − y)( y − y )( K ′ a)
a ∗
y
(18.39)
1 b V (1 − y) L dy
ZT =
S ∫
a (
K y a(1 − y) y − y ∗ )
(18.40)
Puesto que ky a varía con aproximadamente V0.7, y Kya presentará casi la misma
variación cuando controla la película gaseosa, la razón V/Kya no variará mucho. Este
V /S a dy
ZT =
Kya
∫
b y − y∗
(18.41a)
Esta ecuación es la misma que la ecuación (18.16) para gases diluidos, excepto que el
primer término, HOy, es el valor medio para la columna en vez de una constante. Obser-
ve que en este caso ha de utilizarse Ky a de la ecuación (18.12), y no K′y a, puesto que el
término (1 – y)L se introduce en la derivación.
Si la película líquida es la resistencia que controla la transferencia de masa, se uti-
lizan los coeficientes de la película gaseosa para cálculos de diseño de acuerdo con la
ecuación (18.39). Si se utilizan los coeficientes de película líquida y si el factor (1 – x)L
se introduce para permitir una difusión en una sola dirección en el líquido, se deriva una
ecuación similar a la ecuación (18.41a):
L/S a dx (18.42)
ZT =
Kx a ∫ b x∗ − x
EJEMPLO 18.6 Se va a diseñar una torre empacada con anillos de 1 in. (25.4 mm) para
absorber dióxido de azufre desde aire depurando el gas con agua. El gas que entra contiene
20% de SO2 en volumen, mientras que el gas a la salida no debe contener más de 0.5% de
SO2 en volumen. El agua que entra está libre de SO2. La temperatura es de 80 °F y la presión
total es de 2 atm. El flujo de agua será el doble del valor mínimo, y la velocidad del flujo
de aire (sin considerar el SO2) será de 200 lb/ft2 · h (976 kg/m2 · h). ¿Qué altura de empaque se
necesitará?
Para la absorción de SO2 a 80 °F en torres empacadas con anillos de 1 in. se dispone de las
siguientes ecuaciones29 para los coeficientes de transferencia de masa:
k L a = 0.038Gx0.82
k g a = 0.028Gy0.7 Gx0.25
0.22
0.20 yb
ión
rac
0.18 xb
pe
eo
yi
ad
0.16 xi
e
Lín
280
0.14 260
io
240
br
li
ui
0.12 220
eq
de
y 200
a
rv
Cu
0.10 180
f
0.08 140
yi
xi 120
0.06 100
80
0.04 60
40
0.02 20
ya 0.2 0.1 y 0
0 1 2 3 4 5 6 7 8 9
FIGURA 18.19
x × 10 3 Diagrama para el ejemplo 18.6.
Puesto que los coeficientes se midieron a la presión atmosférica, y ky a = kgaP, el valor de kga es
igual a kya. Sin embargo, los coeficientes medidos no están corregidos para la difusión en un solo
sentido y, por lo tanto, corresponden a kya / (1 – y)L. La concentración de SO2 para los ensayos
publicados varía en el intervalo de 3 a 17%, y el valor promedio de (1 – y)L era de aproximada-
mente 0.9. Por consiguiente, la correlación publicada se multiplica por 0.9 para dar
k y a = 0.025Gy0.7 Gx0.25
La alimentación de SO2 es
8.62 × 0.2 = 1.724 mol/ft 2 ⋅ h
La alimentación de aire es
8.62 × 0.8 = 6.896 mol/ft 2 ⋅ h
El gas de salida contiene 0.5% de SO2; por lo tanto, la cantidad de SO2 que sale es
0.005
6.896 × = 0.035 mol/ft 2 ⋅ h
0.995
El SO2 absorbido es, por lo tanto, 1.724 − 0.035 = 1.689 mol/ft2 · h. Se absorbe 98% de SO2.
La alimentación del agua es
191.1 × 1.689 = 322.8 mol/ft 2 ⋅ h
x y 0.005
322.8 = 6.896 −
1− x 1 − y 0.995
1 0.818 20.1
= + = 0.192
K y′ a 12.15 161
1
K y′ a = = 5.21 mol/ft 3 ⋅ h
0.192
calcula a partir de la velocidad másica. Observe que V/S = GM. Cuando y > 0.05, la pendiente
m se toma como 20.1; cuando y < 0.05, m = 15.6.
La altura de la columna se obtiene en la forma que se indica a continuación. A partir de la
ecuación (18.39),
dZ =
(V /S)dy = f ′ dy
(1 − y)( y − y∗ )K y′ a
Los valores de f ′ que se presentan en la tabla 18.2 se han calculado a partir de diferentes
valores de y. La altura total se obtiene por integración gráfica de una gráfica de y versus f ′,
tal como se muestra en la figura 18.19, para dar una altura total ZT de 8.9 ft.
Puesto que la principal resistencia reside en la película de líquido, es posible utilizar la
ecuación (18.43) para una comprobación aproximada de la altura de la columna. La ecuación
análoga a la ecuación (18.37) es
1 1 (1 − y) L
= +
K x′ a k x a mk y a
Por lo tanto,
1 1 0.91
= + = 1.04 × 10 −2
K x′ a 161 20.1 × 11
K x′ a = 96.2 mol/ft 3 ⋅ h ⋅ unidad de diferencia de fracción mollar
De la ecuación (18.24)
322.8
HOx = = 3.36 ft
96.2
Despreciando la curvatura de las líneas de equilibrio y operación,
xb − xa
NOx =
( x∗ − x) L
0.00522
= = 3.0
0.00174
TABLA 18. 2
Integración para el ejemplo 18.6
V V/S
y y yy–−y*y yyii (1 − y)L Ky a f � Z∆Z
∗
S (1 − y)(y − y∗ )Ky a
xb xa
x
12
HOx
1.0 10
2
1
0.5 12 1
2
Tamaño del empaque, in.
0.2
0.1
100 1 000 10 000 100 000
2
Velocidad másica del líquido, Gx, Ib/ft h
aumento de kLa se debe en buena medida al aumento del área de la superficie de contacto
a, el resto procede de un aumento de kL. Para velocidades másicas elevadas, el empaque
está casi completamente húmedo y sólo se produce un ligero aumento de kLa con Gx, lo
que da lugar a que Hx sea aproximadamente proporcional a Gx. Observe que, por lo que se
refiere a la transferencia de masa, los empaques pequeños son sólo ligeramente mejores
que los grandes en el intervalo intermedio de flujos, aunque el área total varía en forma
inversamente proporcional con el tamaño del empaque. En las operaciones industriales
son preferibles los empaques más grandes debido a que presentan una capacidad muy
superior (mayor velocidad de inundación).
Los datos de la figura 18.21 se obtuvieron con velocidades de flujo de gas de 100 a
230 lb/ft2 · h, y en este intervalo no hubo efecto de Gy. Para velocidades de flujo de gas
comprendidas entre el punto de carga y la velocidad de inundación, Hx es ligeramente menor
debido al aumento de la retención del líquido. Sin embargo, para una columna diseñada para
operar a la mitad de la velocidad de inundación, el efecto de Gy sobre Hx es despreciable.
La resistencia de película líquida para otros sistemas se obtiene a partir de los datos
O2-H2O corrigiendo de las diferencias entre difusividad y viscosidad (como referencia,
a 25 °C, Dυ para oxígeno en agua es 2.41 × 10–5 cm2/s y Sc es 381):
n 0.5
1 G µ
Hx = x (18.44)
α µ ρDυ
donde α y n son constantes empíricas que están tabuladas en la bibliografía para alguno
de los antiguos tipos de empaque. El exponente de 0.5 en el número de Schmidt es consis-
tente con la teoría de la penetración, y es el que se esperaría aplicar al líquido fluyendo en
distancias cortas sobre piezas de empaque. El exponente n varía con el tamaño y tipo del
empaque, pero 0.3 se considera un valor típico. La ecuación (18.44) debe utilizarse con
precaución para otros líquidos diferentes al agua, puesto que los efectos de la densidad,
tensión superficial y viscosidad son inciertos.
Cuando un vapor es absorbido en un solvente de peso molecular alto, la velocidad
de flujo molar del líquido será mucho menor que si el agua se usara a la misma veloci-
dad másica. Sin embargo, el coeficiente kx a, que está basado en una fuerza impulsora de
fracción mol, varía en forma inversamente proporcional con respecto al peso molecular
promedio del líquido, y no hay efecto neto de M en Hx.
ρx
kx a = kL a (18.45)
M
Gx / M Gx / ρ x
Hx = = (18.46)
kx a kL a
El coeficiente kLa depende principalmente de la velocidad de flujo volumétrico, difusi-
vidad y viscosidad pero no del peso molecular, así que las correlaciones generales para
kLa o Hx son más simples que en el caso de kx a.
6 FIGURA 18.22
Carga Altura de una unidad de transferencia para
4
la absorción de amoniaco en agua con
anillos Raschig de cerámica de 1 –12 in.
Gx = 500
HOy y Hy, ft
2
1 500
1
0.8
4 500
0.6
HOy
0.4 Hy
de Hy que son alrededor de la mitad de los que se reportan en el caso del amoniaco a las
mismas velocidades másicas. La diferencia se atribuye a bolsas de líquido casi estancado
que contribuye de manera importante a la vaporización pero que se satura rápidamente
en una prueba de absorción de gas.22 Las bolsas estancadas corresponden a la retención
estática, es decir, el líquido que permanece en la columna mucho después de que el flujo
es expulsado. El resto del líquido constituye la retención dinámica, el cual se incrementa
con la velocidad de flujo del líquido. Se han desarrollado21 correlaciones para la retención
estática y dinámica y las correspondientes áreas de contacto y que conviene emplear para
relacionar estrechamente los resultados de absorción y vaporización de gas.
1.0
8.0
6.0
Tamaño de monturas Intalox, in.
4.0 1
2
3.0
2.0 1
1
12
2
1.0
1 000 2 000 3 000 6 000 10 000 20 000 30 000 60 000
Gx, lb/ft2 h
EJEMPLO 18.7 El gas procedente de un reactor tiene 3.0% de óxido de etileno (EO) y
10% de CO2; casi todo el resto está constituido por nitrógeno. Debe recuperarse 98% del EO
absorbiéndolo con agua. El absorbedor operarará a 20 atm, utilizando agua con 0.04% mol
Base de cálculo Entran 100 moles de gas seco y 140 moles de solución.
Entrada Salida
87 N2 87 N2
10 CO2 10 CO2
3 EO 0.06 EO (= 3 × 0.02)
100 97.06
Se supone que la absorción de CO2 es despreciable, así como el efecto del H2O sobre la
composición del gas.
0.06
En la parte superior: x = 0.0004 y= = 0.00062
97.06
0.06
40°C
0.05
0.04
30°C
Línea de
equilibrio
y EO
0.03
Línea de Línea de
operación operación
0.02
0.001
0.01 Línea de
equilibrio
0 0.001 0.002
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 FIGURA 18.24
x EO
Diagrama yx para el ejemplo 18.7.
2.94 + 0.056
En el fondo: x= = 0.0210
140 + 2.94
y = 0.030
∆y
∫ y − y = ∑ (y − y )
dy
NOy = ∗ ∗
L
y y – y* ∆NOy
0.03 0.008 ––
0.015 0.006 2.14
0.005 0.0024 2.55
0.0006 0.0003 4.36
NOy = 9.05 = 9.0 unidades de transferencia
A 40°C, ρy =
30.1
359
× 20 ×
273
313
(
= 1.46 lb/ft 3 0.0234 g/cm 3 )
Gx ρy 1.4 × 18 1.46
= = 0.130
Gy ρx − ρy 1 × 30.1 62.2 − 1.46
Gy2 Fp µ x0.1
= 0.045
(
ρ y ρ x − ρ y gc )
A partir de la tabla 18.1, Fp = 40. A 40°C, µ = 0.656 cP. Por lo tanto,
0.045(1.46)(62.2 − 1.46)(32.2)
Gy2 = = 3.35
40(0.656) 0.1
Gy = 1.83 lb/ft 2 ⋅ s = 6 590 lb/ft 2 ⋅ h
1.4 × 18
Gx = × 6 590 = 5 520 lb/ft 2 ⋅ h
1 × 30.1
Para una velocidad de alimentación de 10 000 moles/h × 30.1 = 3.01 × 105 lb/h,
3.01 × 10 5
S= = 45.7 ft 2 D = 7.6 ft
6 590
Se utiliza una columna de 8 ft de diámetro.
Dυ = 7.0 × 10 −3 cm 2 /s
µ 1.81 × 10 −4
Sc = = = 1.10
ρDυ ( )(
7.0 × 10 −3 2.34 × 10 −2 )
A partir de la tabla 18.1, fp = 1.36 para anillos Pall de 1.5 in. de forma que Hy es menor debido
a este factor. Se supone que Hy varía con Gy0.3 y Gx–0.4. Por lo tanto,
1/ 2 0.3 0.4
1.10 1 6 590 1 500
H y, EO = 1.4 = 1.71 ft
0.66 1.36 500 5 520
De la figura 18.21, Hx,O2 = 0.9 ft para Gx = 1 500. De la ecuación (17.33), Dυ = 2.15 × 10–5
cm2/s para EO en H2O a 40°C:
0.00656
Sc = = 305
1.0 × 2.15 × 10 −5
De acuerdo con la figura 18.24, el valor promedio de m es del orden de 1.0, y a partir de la
ecuación (18.28)
mGM 1.0 × 0.96
HOy = H y + H x = 1.71 + = 2.40 ft
LM 1.4
La absorción de gases puede realizarse en una columna equipada con platos perforados u
otro tipo de platos que normalmente se utilizan en destilación. En ocasiones se elige una
columna de platos perforados en vez de una columna empacada para evitar el problema de
la distribución del líquido en una torre de un gran diámetro y disminuir la incertidumbre
en el escalamiento. Los métodos para determinar el número de platos teóricos y estimar
la eficiencia media de los platos se estudiarán en el capítulo 21.
(18.48)
4 FIGURA 18.25
Absorción de CO2 en
una solución de mono-
etilamina (MEA). (Según
3
Freguia y Rochelle5.)
Línea de
operación
pCO2, kPa
1
Línea de
equilibrio
FIGURA 18.26
y entrada
Líneas de operación para flujos
paralelos y en contracorriente para
Líneas de operación para
la absorción con reacción química
flujos en contracorriente irreversible.
por la parte superior, el gas que sale por el fondo está expuesto al líquido enriquecido,
el que ha absorbido una gran cantidad de soluto; pero si y* = 0, la fuerza impulsora es
simplemente y y NOy se calcula de acuerdo con la ecuación (18.48), de igual manera que
en el caso de flujo en contracorriente.
La ventaja de la operación en paralelo es que no existe la inundación como una limi-
tación, y es factible utilizar un flujo de gas muy superior al normal. Esto permite reducir el
diámetro de columna requerido, y el correspondiente aumento de las velocidades másicas
del líquido y del gas da lugar a elevados coeficientes de transferencia de masa. Es posible
utilizar flujos de líquido tan elevados como 50 000 a 100 000 lb/ft2 · h (70 a 140 kg/m2).
Los valores de Kga o HOy se pueden estimar, aproximadamente, extrapolando los datos de
flujo en contracorriente para el mismo sistema a las velocidades másicas más elevadas.
Absorción en tanques
Los procesos en los cuales la absorción gaseosa va acompañada por una reacción química
normalmente se llevan a cabo en tanques donde las burbujas del gas se dispersan en la
fase líquida. El líquido se mezcla con un agitador o por medio de las mismas burbujas del
gas y también puede existir una apreciable mezcla del propio gas. Esto tiende a disminuir
la conversión. Sin embargo, cuando el gas que no ha reaccionado se puede reciclar o casi
no tiene costo (como el aire), no se requiere una alta conversión del gas. Un importante
ejemplo se da en el uso del aire para oxidar materia orgánica en aguas negras o lodo
residual de alcantarillas.
Uno de los métodos de aereación utiliza inyectores o difusores colocados cerca del
fondo del tanque. Los cabezales de tubería con múltiples orificios se denominan difusores
burdos y producen burbujas de 1 a 5 mm. Los tubos de cerámica porosa o de polímeros
generan burbujas más finas, pero la coalescencia puede aumentar el tamaño de las burbujas.
Los tanques grandes pueden contener más de 100 difusores, cada uno dispersando de 5 a
20 SCFM de aire y desahogando la mayoría del gas hacia la atmósfera. Otros dispositivos
para dispersión de gases incluyen turbinas cerca del fondo del tanque con un inyector en
la parte inferior e impulsores de flujo axial montados cerca de la superficie, para atrapar
el gas de el espacio del vapor, o para aspersar el líquido hacia el aire.23
En los equipos de aereación, la velocidad de transferencia de masa se controla
por medio de la resistencia de la película del líquido, debido a la baja solubilidad del
SÍMBOLOS
Letras griegas
α Constante en la ecuación (18.24)
γA Coeficiente de actividad del componente A
∆C Fuerza impulsora de concentración, g/g o ppm; ∆Ca, en el líquido que entra;
∆Cb, en el líquido de salida
∆H Calor de disolución, kcal/g mol; ∆Hυ, calor de vaporización
∆P Caída de presión, in. de agua/ft de empaque; ∆Pinund, en la inundación
∆xL Media logarítmica de x*b – xb y x*a – xa
∆yL Media logarítmica de yb – y*b y ya – y*a
ε Porosidad o fracción de huecos de la sección empacada
µ Viscosidad, cP o lb/f · h; µx, del líquido
ν Viscosidad cinemática, µ/ρ, m2/s o ft2/s
ρ Densidad, kg/m3 o lb/ft3; ρx, del líquido; ry, del gas
ρM Densidad molar, kg mol/m3 o lb mol/ft3; ρMx, del líquido; ρMy, del gas
φ Factor de avance en la absorción con reacción química, adimensional
PROBLEMAS
18.1. El diseño de una planta requiere un absorbedor para recuperar 95% de acetona desde una
corriente de aire utilizando agua como líquido absorbente. El aire que entra contiene 14%
mol de acetona. El absorbedor está provisto de enfriamiento y opera a 80 °F y 1 atm para
dar lugar a un producto que contiene 5.0% mol de acetona. La alimentación de agua a la
torre contiene 0.02% mol de acetona. La torre ha de diseñarse para operar al 50% de
la velocidad de inundación. a) ¿Cuántas libras de agua por hora deben introducirse como
alimentación de la torre, si el flujo de gas, medido a 1 atm y 32 °F, es de 500 ft3/min?
b) ¿Cuántas unidades de transferencia se necesitan con base en la fuerza impulsora global
de la fase gaseosa? c) Si la torre se empaca con anillos Raschig de 1 in., ¿cuál será la altura
del empaque?
Para el equilibrio. Suponga que pA = P′AγAx, donde ln γA = 1.95(1 – x)2. La presión de vapor
de acetona a 80 °F es 0.33 atm.
18.2. Un absorbedor debe recuperar 99% de amoniaco de una corriente de alimentación de aire-
amoniaco, utilizando agua como líquido absorbente. El contenido de amoniaco en el aire
es de 20% mol. La temperatura del absorbedor se mantiene a 30 °C mediante serpentines de
enfriamiento y la presión es de 1 atm. a) ¿Cuál es la velocidad de flujo mínimo de agua?
b) Para un flujo de agua 40% superior al valor mínimo, ¿cuántas unidades globales de
transferencia basadas en la fase gaseosa se necesitan?
18.3. Un gas soluble se absorbe en agua utilizando una torre empacada. La relación de equilibrio
puede tomarse como ye = 0.06xe. Las condiciones terminales son las siguientes:
x 0 0.08
y 0.001 0.009
18.4. Una columna de absorción se alimenta por el fondo con un gas que contiene 5% de benceno
y 95% de aire. Por la parte superior de la columna se introduce un aceite de absorción no
volátil que contiene 0.2% en peso de benceno. Otros datos son:
18.5. Una corriente de vapor que contiene 3.0% mol de benceno se va a purificar mediante lavado
con aceite en un absorbedor empacado con el fin de reducir la concentración de benceno
en el gas a 0.02%. El aceite tiene un peso molecular promedio de 250, una densidad de
54.6 lb/ft3 y contiene 0.015% de benceno. El flujo de gas es de 1 500 ft3/min a 25 °C y
1 atm. a) Si el sistema opera isotérmicamente a 25 °C con un flujo de líquido de 14 000 lb/h,
¿cuántas unidades de transferencia se requieren? b) ¿Cuántas unidades de transferencia se
necesitarían si operase adiabáticamente? c) ¿Cuál sería el principal efecto de operar con
un aceite de menor peso molecular, tal como M = 200?
18.6. Una corriente residual acuosa que contiene 1.0% en peso de NH3 se va a desabsorber con aire
en una columna empacada con el fin de separar 99% de NH3. ¿Cuál es el flujo mínimo de
aire, en kilogramos de aire por kilogramo de agua, si la columna opera a 20 °C? ¿Cuántas uni-
dades de transferencia se requieren para un flujo de aire igual al doble del valor mínimo?
18.7. Una columna de 8 ft de diámetro empacada con 20 ft de monturas Berl de 1 –12 in. opera
con flujo de aire a 1.5 atm y 40 °C. La torre está aparentemente próxima a la inundación,
puesto que ∆ρ = 24 in. de agua. La velocidad másica del líquido es 8.5 veces la del gas.
a) Si la torre se reempacara de nuevo con monturas Intalox de 1.5 in. ¿cuál sería la caída
de presión? b) Si la caída de presión fuese la misma que con las monturas Berl, ¿cuánto se
podrían aumentar las velocidades de flujo de operación?
18.8. Un absorbedor debe separar 99% del soluto A de una corriente gaseosa que contiene 4%
mol de A. Las soluciones de A en el solvente siguen la ley de Henry, y el aumento de tem-
peratura del líquido es despreciable. a) Calcule NOy para la operación a 1 atm, utilizando
un flujo de líquido libre de soluto 1.5 veces el valor mínimo. b) Para la misma velocidad
de líquido, calcule NOy para la operación a 2 atm y a 4 atm. c) ¿Sería el efecto de la presión
sobre NOy parcialmente compensado mediante una variación en HOy?
18.9. Un gas que contiene 2% de A y 1% de B se va a absorber con un solvente en el que A es 5
veces más soluble que B. Demuestre que utilizando dos columnas en serie con regeneración
separada del líquido de cada columna se podrían recuperar A y B en una forma relativa-
mente pura. Utilice un diagrama yx para representar las líneas de equilibrio y operación en
la absorción simultánea de A y B, y estime la relación entre A y B en el líquido procedente
del primer absorbedor.
18.10. Un absorbedor empacado con monturas Intalox de 1 in. opera a 50 °C y 10 atm con una
velocidad másica de líquido de 5 veces la velocidad másica del gas. Suponiendo que el gas
y el líquido son similares al aire y el agua, ¿qué velocidad másica del gas dará una caída
de presión de 0.5 in. de H2O/ft de empaque? Utilice las correlaciones generalizadas para
obtener el efecto del cambio de las propiedades físicas, y aplique una corrección a los datos
de la figura 18.4.
18.11. Una columna empacada con Mellapak 250Y es operada a 1 atm con aire y una solución
acuosa a Gx = 8Gy y 60% de la velocidad de inundación (véase figura 18.8). a) ¿Cuál es la
velocidad superficial del gas y la velocidad másica? b) Si la presión se cambiara a 2 atm,
¿qué velocidad másica y superficial debería utilizarse?
18.12. Una columna empacada se utiliza para remover tolueno a partir de agua por desorción
con aire. El agua contiene 600 ppm de tolueno (en peso), y la concentración debe ser re-
ducida a 2 ppm. La columna operará a 20 °C y una presión media de 1.1 atm. La relación
de equilibrio es
Ptol = 256 x
donde Ptol = presión parcial del tolueno, atm, en la fase gaseosa, y x = fracción molar del
tolueno en el líquido. a) Determine la velocidad mínima del aire en g mol/min para una
velocidad de agua de 100 L/min. b) Calcule NOx para una velocidad del aire del doble de la
mínima. c) Determine la altura empacada requerida si Hy = 0.7 m y Hx = 0.6 m.
18.13. Si una columna empacada con monturas Berl de 1 in. es reempacada con anillos metálicos
Pall, ¿cuáles serán los principales cambios en las características de operación?
18.14. Calcule la fracción de absorción de oxígeno de una burbuja de 2 mm que se eleva 15 ft
en agua pura a 20 ºC. Suponga que la burbuja se eleva a su velocidad terminal y tiene el
mismo coeficiente externo de transferencia de masa que una esfera rígida. La solubilidad
del oxígeno es de 9 mg/L y su difusividad es de 2 10–5 cm2/s.
18.15. Se empacará una columna de destilación de benceno-tolueno con un empaque IMTP de
1 –12 in. Para las condiciones en la parte más baja de la columna, donde L/V = 1.3, ¿cuál es
la caída de presión a esa relación?
18.16. Se diseñó un absorbedor de gas en columna empacada para extraer 98% de A del aire a un
flujo del líquido de 1.2 veces el mínimo. El solvente que ingresa está libre de A. a) ¿Cuántas
REFERENCIAS
1. Cheng, A.T.Y. y G.E. Storms. Trabajo presentado en la Pardue Industrial Waste Conference,
1996.
2. Danckwerts. P.V. Gas-Liquid Reactions. Nueva York: McGraw-Hill, 1970.
3. Eckert, J.S. Chem. Eng. Prog. 66(3): 39 (1970).
4. Fair, J.F. y J.L. Bravo. Chem. Eng. Prog. 86(1): 19 (1990).
5. Freguia, S. y G.T. Rochelle. AIChE J. 49: 1676 (2003)
6. Gmehling, J., U. Onken y W. Arlt. Vapor-Liquid Equilibria Data Collection, vol. 1. Dechema,
Frankfurt/Main, 1979.
7. Harriott, P. Environ. Sci. Tech. 23: 309 (1989).
8. Harriott, P. Chemical Reactor Design. Nueva York: Marcel Dekker, 2003.
9. Kister, H.Z. y D.R. Gill. Chem. Eng. Prog. 87(2): 32 (1991)
10. Kohl, A. y F. Reisenfeld. Gas purification, 2a. ed. Houston: Gulf Publishing, 1974.
11. Lincoff, A.H. y J.M. Gossett. International Symposium on Gas Transfer at Water Surfaces.
Ithaca, N.Y.: Cornell University, junio de 1983.
12. Metcalf y Eddy, Inc. Wastewater Engineering, 4a. ed. Wakefield, Mass.: Metcalf y Eddy, 2003,
pp. 437-54.
13. Norton Chemical Process Products Corp., Akron, Ohio, 1987.
14. Perry, D., D.E: Nutter y A. Hale. Chem. Eng. Prog. 86(1): 30 (1990).
15. Perry, J. H. Chemical Engineers’ Handbook, 6a. ed. Nueva York: McGraw-Hill, 1894; a) pp.
3-101; b) pp. 18-23; c) pp. 18-39.
16. Perry, R.H. y D.W. Green (eds.). Perry’s Chemical Engineers’ Handbook, 7a. ed. Nueva York:
McGraw-Hill, 1975, pp. 14-43.
17. Reynolds, T.D. y P. A. Richards. Unit Operations and Processes in Environmental Engineer-
ing. 2a. ed. Boston: PWS, 1996, p. 500.
18. Robbins, L.A., Chem. Eng. Prog. 87(5): 87 (1991).
19. Sherwood, T.K. y F.A. L. Holloway. Trans. AIChE 36: 21, 39 (1940).
20. Sherwood, T.K., R.L. Pigford y C. R. Wilke. Mass Transfer, Nueva York: McGraw-Hill, 1975,
p. 442.
21. Shulman, H.L., C.F. Ullrich, A.Z. Proulx y J.O. Zimmerman. AIChE J. 1: 253 (1955).
22. Shulman, H.L., C.F. Ullrich y N. Wells. AIChE J. 1: 247 (1955).
23. Sincero, A.P. y G.A. Sincero. Environmental Engineering; A Design Approach. Englewood
Cliffs, N.J.: Prentice-Hall, 1996, p. 331.
24. Sperandio, A., A. M. Richard y M. Huber. Chem. Ing. Tech. 37: 322 (1965).
25. Spiegel, L. y W. Meier. I. Chem. E. Symp. Ser. 104: A203 (1987).
26. Stedman, D.F. Trans. AIChE 33:153 (1937).
27. Strigle, L., Jr. Random Packings and Packed Towers. Houston, Tex.: Gulf Publishing, 1987.
28. VonStockar, U. y C.R. Wilke. Ind. Eng. Chem. Fund. 16: 88, 94 (1977).
29. Whitney, R.P. y J.E. Vivian. Chem. Eng. Prog. 45: 323 (1949).