UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

I. OBJETIVOS

 Cálculo de la caída de presión de la columna empacada.

 Cálculo del NTOG y el HTOG de la columna.
 Cálculo del coeficiente “Cf” para el empaque de anillo Raschig de 1/2¨ de vidrio.

II. FUNDAMENTO TEÓRICO

Torres empacadas

Las columnas empacadas son los equipos de transferencia de masa más empleados en la
industria de transformación físico – química. Al inicio su uso era restringido mayormente a
laboratorios, plantas piloto y a columnas industriales de diámetro pequeño o servicios de
requerimientos de muy baja caída de presión y alta resistencia a la corrosión.
Actualmente, su aplicación se considera en todas las operaciones que implican
transferencia de masa, de calor y las que requieran el contacto directo entre fases.

En la siguiente tabla se presentan las operaciones en las que se pueden emplear torres
empacadas:

Operación Aplicaciones
Absorción física Remoción de contaminantes del aire,
eliminación de olores
Absorción química Absorción de CO2 en control de
emisiones y en el proceso Solvay.
Absorción en procesos de producción de
HNO3, H2SO4, Cl2 y HF
Desorción Decarbonatación de agua, despojamiento
de amoníaco.
Destilación Rectificación al vacío de crudo.
Separación de sustancias termolábiles.
Separación de isótopos.

Laboratorio de Operaciones Unitarias II (PI-136/A) 1

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

Extracción líquido - líquido Recuperación de compuestos orgánicos

de corrientes acuosas.
Deshumidificación
Desalinización de agua.

Enfriamiento Plantas térmicas y nucleares.

Enfriamiento de cloro gaseoso.

Propiedades de las torres de relleno

Las propiedades que debe poseer una torre de relleno son las siguientes:

Bajo peso por unidad de volumen: esto afecta no solamente al peso que ha de soportar la
torre, sino también a la envolvente de la torre. Un relleno que se vacía en la torre al azar,
puede ejercer un empuje local sobre las paredes, y si el relleno tiene un alto peso unitario,
esto puede afectar al coste de la construcción de la torre.

Gran sección transversal libre: esto tiene importancia porque afecta a la perdida de carga
por fricción a través de la torre y por consiguiente a la potencia necesaria para hacer
circular el gas, también una pequeña sección transversal libre, significa una velocidad
elevada para una cantidad dada de gas y por encima de ciertas velocidades libres existe
tendencia a soplar el líquido fuera de la torre.

Gran volumen libre: en algunos casos, tal como la absorción de dióxido de nitrógeno en
que hay que prevenir cierto tiempo para que se verifiquen las reacciones en la fase
gaseosa, este apartado puede tener gran importancia.
Pequeño peso del líquido retenido: esto generalmente, es una ventaja, puesto que
disminuye la carga sobre la torre tan rápidamente como es posible. En algunos casos
puede ser una gran desventaja, especialmente cuando la reacción entre el líquido y el gas
es lenta, o cuando la solubilidad del gas en el líquido no es grande. En torres acidas es de
desear tener pequeña cantidad de líquido retenido para disminuir el peligro de accidentes
cuando la torre se vacía.

Laboratorio de Operaciones Unitarias II (PI-136/A) 2

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

III. DATOS TEORICOS

Tabla N°1. Datos de la columna

Altura de empaque (m) 1.04

Diámetro de columna (pulg) 4
Tipo de empaque Anillos raschig de vidrio de ½”

Tabla N°2. Propiedad del gas y líquido

Densidad del aire (Kg/m3 ) 1.22

Densidad del agua (Kg/m3 ) 990
Viscosidad del agua (Kg/m. s) 0.00105

Laboratorio de Operaciones Unitarias II (PI-136/A) 3

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

IV. CUESTIONARIO

1. Tabla de resultados del experimento 1: Caída de presión

L = 30 lb/h Q = 10 pies3/min Promedio

∆P (mm H2O) 68 69 70 69 69
Volumen (ml) 250 240 250 245 245
Tiempo (s) 60 60 60 60 60
LSALIDA (ml/s) 4.08

L = 30 lb/h Q = 12 pies3/min Promedio

∆P (mm H2O) 97 96 96 95 96
Volumen (ml) 125 125 120 130 125
Tiempo (s) 30 30 30 30 30
LSALIDA (ml/s) 4.16

L = 30 lb/h Q = 14 pies3/min Promedio

∆P (mm H2O) 132 132 130 132 131.5
Volumen (ml) 135 120 120 125 125
Tiempo (s) 30 30 30 30 30
LSALIDA (ml/s) 4.16

L = 30 lb/h Q = 16 pies3/min Promedio

∆P (mm H2O) 187 187 190 188 188
Volumen (ml) 120 115 110 125 117.5
Tiempo (s) 30 30 30 30 30
LSALIDA (ml/s) 3.91

Para L = 30 lb/h constante:

Q (pies3/min) ∆P (mm H2O) Volumen (ml) Tiempo (s) LSALIDA (ml/s) LSALIDA (pies3/min)
10 69 245 60
12 96 125 30
14 131.5 125 30
16 188 117.5 30

Laboratorio de Operaciones Unitarias II (PI-136/A) 4

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

2. Calculo de Cf

Para realizar el cálculo se empleará el siguiente gráfico (Fuente: Operaciones de

transferencia de masa, Robert Treybal)

Figura N°2. Inundación y caída depresión en torres empacadas

L (lb/h) 30
L' 0.48
(kg/m2*s)
Q ∆P (mm G' ∆P/Z 𝑳′ 𝝆𝑮 𝟏/𝟐 𝑮′𝟐 𝑪𝒇 𝝁𝑳 𝟎.𝟏 𝑱 Cf
(pies3/min) H2O) (kg/m2*s) (Pa/m) ′
( )
𝑮 𝝆𝑳 − 𝝆𝑮 𝝆𝑮 (𝝆𝑳 − 𝝆𝑮 )𝒈𝒄
10 69 0.777 669.3 0.02177 0.13 516.18
12 96 0.932 931.2 0.01814 0.17 468.755
14 131.5 1.087 1275.6 0.01555 0.21 425.425
16 188 1.243 1823.6 0.01361 0.28* 431.109

(*) este dato se obtiene extrapolando con los datos de la gráfica para una abs=0.01361.

Cf (promedio) = 460.368

Laboratorio de Operaciones Unitarias II (PI-136/A) 5

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

Conclusiones:

 Cabe resaltar principalmente que estos datos obtenidos gráficamente, no son exactos.
Esto ocasiona error
 A medida que el flujo de gas comienza a incrementarse el Cf disminuye son inversamente
proporcionales, lo mismo ocurre cuando aumenta la caída de la presión respecto a Cf.
 Con los datos obtenidos en el laboratorio, solo pudimos calcular tres Cf el otro se hace
extrapolando.
Para la abs = 0.01361

Ord P
0.27795 1823.6
0.2 1200
0.15 800

3. Tabla de resultados del experimento 2: Corrida de desorción Gaseosa

Altura del Empaque (m) = 1.01

Flujo de Gas (pies3/min) = 14
L = 30 lb/h
Alimentación volumen de muestra (mL) = concentración del ácido (N) = 0.1
15
volumen usado de ácido (ml) 9.6 9.4
Líquidos de salida de los volumen de muestra (mL) = concentración del ácido (N) = 0.1
fondos 15
volumen usado de ácido (ml) 6.9 4 3.7 3.5 3.2

Volumen (ml) 45 35 50 35
Tiempo (s) 10 10 10 10

4. Porcentaje de desorción obtenido y porcentaje de desorción obtenido si la

columna fuera de altura infinita.

Calculo de las concentraciones del líquido

 Alimentación “se hace mediante una titulación con HCL 0.1N”

𝑁(HCl) ∗ V(HCl)
[𝐶𝑂𝑁. 𝑁𝐻3]𝑎𝑙𝑖𝑚𝑒𝑛𝑡𝑎𝑐𝑖𝑜𝑛 =
V(NH3)

Laboratorio de Operaciones Unitarias II (PI-136/A) 6

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

Muestra Volumen HCl (ml) Concentra (NH3) en la

alimentación
1 9.6 0.064
2 9.4 0.062
Concentración Promedio del NH3(N): alimentación 0.063

 Salida “ de igual manera se realiza mediante titulación con HCL, 0.1N”

Volumen Concentración
Flujo Muestra
HCl mol/L

45ml 15 ml 6 0.04
(10s)

35ml 15 ml 4 0.027
(10s)

50ml 15 ml 3.7 0.025

(10 s)

35ml 15 ml 3.5 0.023

(10s)

38ml 15ml 3.2 0.021

(10s)

Para este caso consideraremos la concentración estacionaria 0.021 N

 Para el cálculo de las fracciones molares de entrada y salida

Se considera para sistemas agua y amoniaco como una solución diluida y las
características del fluido serán netamente del agua

 Haciendo el cálculo “ALIMENTACION”

n(NH3)
X2 = n(NH3)+𝑛(H2O) = 0.00114

X2=( X2 )/ (1 - X2) =0,00114

Laboratorio de Operaciones Unitarias II (PI-136/A) 7

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

 Haciendo el cálculo “SALIDA”

molNH3
X1 = 0.000385 moles totales

X1=( X1 )/ (1 – X1) =0,000385

Hallando el porcentaje de desorción:
% de desorción = (X2 – X1) / (X2)

% de desorción =66.3 %

Conclusión:

 El porcentaje es considerable, un 66. 3 % puede ser medianamente optimo, pero

se espera tener un mayor porcentaje como por ejemplo 90 %
 Para los cálculos de las fracciones parciales se uso netamente del agua esto porque
la solubilidad del NH3 en agua es muy pequeña a temperaturas en las cuales
hemos trabajado.

Hallando el porcentaje de desorción si la columna fura de altura infinita:

Graficaremos nuestra curva de equilibrio de 16°C y 20°C esto se explicara en su debido
momento en la pregunta 5, “usando los datos del Perry” el cuadro de solubilidad está en
la pregunta 6.
Interpolamos para las temperaturas de 16 y 20 °C
Para ello se hallaron los polinomios de cada curva:

Laboratorio de Operaciones Unitarias II (PI-136/A) 8

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

Curvas de equilibrio 16°C y 20 ° C

0.45
0.4 Y* = 0.7133X3 + 0.6743X2 + 0.6157X - 2E-05
0.35 Y* = 1.2952X3 + 0.7473X2 + 0.7462X - 2E-05
0.3
0.25
0.2
0.15
0.1
0.05
0
-0.05 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4

 Obtenemos los polinomios

o 20°C; Y* = 1.295X3 + 0.7473X2 + 0.7462X - 0.00002

Como de tenemos una concentración en el tope de X2= 0.00114, Obtenemos en el

equilibrio Y2* = 0.000831

o 16°C; Y* = 0.7133X3 + 0.6743X2 + 0.6157X - 2E-05

De igual manera con X1= 0.000385, obtenemos Y1* = 0.000217

Ahora si obtenemos 2 puntos y podemos graficar nuestra curva de equilibrio

Puntos a graficar serán:

X Y*

0 0

0.000385 0.000217

0.00114 0.000831

Laboratorio de Operaciones Unitarias II (PI-136/A) 9

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

Curva de equilibrio de nuestro proceso

Curva de equilibrio para nuestro caso

0.0009
0.0008
0.0007
0.0006
0.0005
0.0004
0.0003
0.0002
0.0001
0
0 0.0002 0.0004 0.0006 0.0008 0.001 0.0012
-0.0001

Se realiza un ajuste polinómico y a partir de ello extendemos la curva para

concentraciones diluida de NH3

Procedemos al cálculo de la fracción (y2) salida del gas:

Calculo de Flujo de Gas:
G1=14 pies3 /min = 6,61 L/S
P(aire) ∗ V(aire) (1atm) ∗ (6,61)
n(aire) = =
R .T L
(0,082057 atm. ) ∗ (300,15ºK)
mol. K
𝐺1 = 0,2683mol/seg
𝑮𝑺 = 𝑮𝟏 = 𝟎, 𝟐𝟔𝟖𝟑𝒎𝒐𝒍/𝒔𝒆𝒈
Calculo del Flujo de Liquido
L2=30 lb/h = 3.83gr/s ; (MH2O)=18
Pero la masa molecular de la solución:
3.833gr/s
Entonces L2= = 0.2129 mol/s
18 mol/g

Ls= L2*(1-X2); X2=0.00114

Ls= 0.2126 moles de agua/s

Laboratorio de Operaciones Unitarias II (PI-136/A) 10

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

Para el cálculo de Y2 haremos un balance de masa

GS*(Y2 – Y1) = LS*( X2 - X1)

0.2683*(Y2 – 0) = 0.2126*(0.00114-0.000385)

Y2= 0.000598
Ahora tenemos los puntos: (X1,Y1) y (X2,Y2) graficamosla curva de operación y equilibrio

Curva de equilibrio para nuestro caso

0.0009
0.0008
0.0007
0.0006
0.0005
0.0004
0.0003
0.0002
0.0001
0
-0.0001 0 0.0002 0.0004 0.0006 0.0008 0.001 0.0012

Para lograr una altura infinita solo corremos la curva de operación paralelamente hacia la
izquierda hasta chocar con la curva de equilibrio.

Obtenemos resolviendo el polinomio de la curva de equilibrio X2 infinito = 0.00049

Curva de equilibrio para nuestro caso

0.0009
0.0008
y = 219.2x2 + 0.4796x - 9E-19
0.0007
0.0006
0.0005
0.0004
0.0003 y = 0.7921x - 0.0001
0.0002
y = 0.7921x - 0.0003
0.0001
0
-0.0001 0 0.0002 0.0004 0.0006 0.0008 0.001 0.0012

Laboratorio de Operaciones Unitarias II (PI-136/A) 11

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

Para el cálculo del porcentaje con la altura infinita:

% de desorción = (X2infinito – X1) / (X2infinito)

% de desorción =38.7%
Conclusión:

 Notamos que el % de desorción a disminuido si es que tiene la misma concentración de

entrada
 Para caso que tenga una concentración de 0.0001 el % aumentaría.

5. Estimación del rango de temperatura de operación en la columna.

T (1) y T (2), las medimos con termómetro, cercanas a las temperaturas del borde superior
y borde inferior
T (Borde inferior del relleno) = 21°C
T (1) = 16 ° C
Como se observa las temperaturas en la parte inferior se diferencias por 5ºC. Entonces
esto quiere decir que el calor de disolución es poca, ya que la solución es diluida, entonces
podemos asumir que en la parte superior tiene una misma variación de temperatura
entonces:
T (Borde superior del relleno) – T (2) = 5ºC, siendo T(2)= 25ºC,
Entonces concluimos que T (Borde superior del relleno)= 20ºC

Laboratorio de Operaciones Unitarias II (PI-136/A) 12

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

6. Datos y fuente consultadas para la curva de equilibrio

Manual de ingenieros perry tabla 12.3 solubilidad del NH3 en agua

x y*0°C y*10°C y*20°C y*40°C

0,1059 0,0330 0,0550 0,0916 0,2197
0,1588 0,0562 0,0922 0,1500 0,3592
0,2118 0,0842 0,1362 0,2184 0,5197
0,2647 0,1178 0,1895 0,2987 0,7026
0,3176 0,1566 0,2500 0,3921 0,9105
0,4235 0,2500 0,3961 0,6184 1,4013
0,5294 0,3618 0,5776 0,9026 2,0000
0,6353 0,5000 0,7895 1,2434 2,8026
0,7412 0,6579 1,0263 1,5395 3,6316
0,8471 0,8368 1,2987 1,9079 4,3421

Laboratorio de Operaciones Unitarias II (PI-136/A) 13

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

7. Ecuación y grafica de la curva de equilibrio y línea de operación

Curva de equilibrio para nuestro caso

0.0009
0.0008
y = 219.2x2 + 0.4796x - 9E-19
0.0007
0.0006
0.0005
0.0004
0.0003 y = 0.7921x - 0.0001
0.0002
y = 0.7921x - 0.0003
0.0001
0
-0.0001 0 0.0002 0.0004 0.0006 0.0008 0.001 0.0012

Para la ecuación de la curva de equilibrio hacemos un ajuste polinómico y observamos que

nuestra curva de operación es curva
Y* = 219.2X2 + 0.4796X - 9E-19
Para la curva de operación observamos que es una línea recta con pendiente
(Y2 – Y1) / (X2 - X1), entonces la ecuación será:
Y*= 0.7921X - 0.0003

Conclusión

 La curva de equilibrio la tomamos de las ecuaciones polinómicas formadas de las

curvas de 20°C y 16°C, usando sus ecuaciones polinómicas
 Y la curva de equilibrio se obtiene con los cuatro datos hallados
experimentalmente.
 Note que para una desorción la línea de operación esta debajo de la curva de
equilibrio.

Laboratorio de Operaciones Unitarias II (PI-136/A) 14

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

8. Calculó del Ntog en forma rigurosa y en forma simplificada

𝑌1
𝑑𝑌 1
𝑁𝑇𝑂𝐺 = ∫ + ∗𝐿𝑁[(1+𝑌2)/(1+𝑌1)]
𝑌−𝑌∗ 2
𝑌2

 Forma rigurosa soluciones diluidas Ntl =Ntog

Escogemos puntos:

X Y Y*
0.000385 0 0.000217
0.00114 0.000598 0.000831

Usando método numéricos 𝑁𝑇𝑂𝐺 = 2.65

 Calculo simplificado como es una solución diluida:

Ecuación del treybal 8.50

NTOG = Ln ((y1-mx2)(1-1/A)/(y2-mx2)+1/A)
1-1/A
Donde:
m es la pendiente de la curva de equilibrio = 0.769
𝐿𝑠
𝐴=
𝑚𝐺𝑠
𝑁𝑇𝑂𝐺 = 2.08
Conclusiones

 Se podría usar tranquilamente el método simplificado ya que nuestras soluciones

son muy diluidas y trabajamos por eso con los Y en vez de y (fracciones molares)

9. Cálculo de Htog en m y del Kya en Kmol/sm3 , discusiones y conclusiones

𝑮′
𝑯𝒕𝑶𝑮 =
𝑲𝒚 𝒂(𝟏 − 𝒚)∗𝑴

Tenemos anillos rashing de ½’’

Para un flujo L’= 0.48

Según la tabla 6.4 del Treybal

m = 28.01

Laboratorio de Operaciones Unitarias II (PI-136/A) 15

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

n = 0.2323L’- 0.30
p = -1.04
808G′ n
Entonces aAW = m( ρ 0.5
) L′p
G
φLo
aA = aAW
φLoW
Donde ϕ se consigue de la tabla 6.5 del Treybal

ds = 0.01774

φLs =0.0045

β = 0.331

φLsW =0.032

φLtW =0.0080

H = 1.008
φLoW = 0.0244
φLo = 0.0246
φLt = 0.0291
aAW =17.05
aA =17.19
De la tabal 6.3 del Treybal
ϵ = 0.63
ϵLo = 0.63-0.0291 = 0.6
Sc = 0.63
FG = 0.044 Kmol/m2.s
FGa = 0.7634 kmol/m3.s
Kya = FGa/(1-y)*M
Kya = 0.7639 kmol/sm3
HtOG = 0.3232 m

Laboratorio de Operaciones Unitarias II (PI-136/A) 16

UNIVERSIDAD NACIONAL DE INGENIERÍA 2015-1

10. Utilizar dos correlaciones de la bibliografía para el cálculo de HtOG y Zcal.

Comparar Zcal con Zreal.

Zcal = HtOG*NtOG

Zcal= 0.85 m
Zreal = 1.01 m

% ERROR = 15%

Discusión y Conclusiones

 Es posible notar que los cálculos realizados se acercan bastante al Zreal y que el
error creció por la manera de operar los datos obtenidos, ya que no se tenía datos
exactos como en el caso de las curvas de caída de presión.

V. Bibliografía

 Robert E. Treybal- “Operaciones de transferencia de masa” -2° Edición-Mac Graw

Hill - Págs. 157-235
 Héctor Jr De La Hoz Siegler – Universidad Nacional de Colombia - “Diseño de un
módulo para el dimensionamiento hidráulico de torres empacadas”
http://www.revistavirtualpro.com/files-bv/20051001/20051001-015.pdf
 Bach. Misari Cuestas, Marlon y Bach. Orihuela Solier, Elena - Tesis de la
Universidad del Centro del Perú – Huancayo
http://cip.org.pe/imagenes/temp/tesis/42960616.pdf

Laboratorio de Operaciones Unitarias II (PI-136/A) 17