UNIVERSIDAD NACIONAL DE SAN CRISTÓBAL DE HUAMANGA

FACULTAD DE INGENIERÍA QUÍMICA Y METALURGIA

ESCUELA PROFESIONAL DE INGENIERÍA QUÍMICA

PRÁCTICA N°2

ASIGNATURA: TRANSFERENCIA DE MASA I

EVALUACIÓN DE COLUMNA DE RELLENO

PROFESOR DE PRÁCTICA: Ing. ARIAS JARA, Alfredo

PROFESOR DE PRÁCTICA: Ing. ZEGARRA VILA, Anna Elinor

ALUMNAS

 GUTIERREZ ORIUNDO, Ruth Karina

 YUPARI ASTO, Naysha Andrea

Día de practica: jueves Hora: 5:00-7:00 pm

Fecha de entrega: 13/05/2021

AYACUCHO- PERÚ

2021
 EVALUACIÓN DE COLUMNA DE RELLENO

I. OBJETIVOS
Estudiar la respuesta de la columna a los caudales de aire y agua (sistema
bifásico) circulando al interior de la columna en contracorriente.
Evaluar la caída de presión de una columna de relleno frente a las variaciones de
flujo de aire y agua.
Determinar los parámetros de flujo y localizar los puntos de carga e inundación
de la columna

II. FUNDAMENTO TEÓRICO

2. COLUMNAS Y TORRES

Son dispositivos cilíndricos de funcionamiento vertical en cuyo interior

se disponen materiales de contacto. Las columnas tienen un cierto diámetro y
son de considerable altura de acuerdo a la capacidad de procesamiento y al
grado de separación que se espera en la columna, por la cual circula el sistema
bifásico, generalmente en contracorriente, tal que la fase pesada se inyecta por la
parte superior y la fase ligera por la parte inferior.

Las columnas o torres de separación son los equipos más ampliamente

utilizadas en la transferencia de masa, ya sea en operaciones como: destilación
batch con rectificación, destilación continua, destilación reactiva, absorción de
gases, desorción de gases, humidificadores adiabáticos, extracción líquida, etc.
Son dispositivos que facilitan el contacto continuo entre las fases que circulan
por su interior, incrementando el área de transferencia de masa interfacial
gracias a los dispositivos de contacto interno instalados y promoviendo la
separación de componentes.

Las columnas vacías son usadas para estudios de investigación de

transferencia de masa de pared mojada; sin embargo, a nivel industrial las
columnas no son vacías, por el contrario, se dispone en su interior diversos
dispositivos de contacto mecánico (no confundir con catalizadores o lechos de
interacción), cuya única finalidad es de incrementar la superficie interfacial de
contacto bifásico y bajo algún grado de turbulencia provocar una mayor
velocidad de transferencia de masa. (Arias, 2011, pág. 62)

Figura. 1

Torre de fraccionamiento de destilación

[Link] DE COLUMNA DE RELLENO

Columnas de relleno están conformadas por la carga aleatoria, randómica

o al azar de ciertos materiales particulados inertes denominados como relleno o
empaque dispuestos por una determinada altura al interior de la columna. Los
rellenos tienen la función principal de dispersar al máximo la corriente gas o
vapor con la corriente líquida que fluyen al interior de la columna, vale decir
aumentan la superficie interfacial de contacto de fases para incrementar a la vez
la tasa de transferencia de masa del componente deseado, contacto que sucede
continuamente a lo largo de toda la columna, razón por la cual se denomina
como contacto continuo diferencial; tal que la transferencia de masa entre el gas
y el líquido a cualquier altura de la columna siempre está sujeta a la fuerza
motriz del gradiente de potencial químico entre las fases, definida por el
coeficiente de transferencia de masa interfacial vigente en las condiciones
hidrodinámicas y térmicas de los fluidos en contacto.
Los rellenos o materiales particulados pueden ser de materiales metálicos,
polímeros, cerámicos o vidrio; poseen diversos tamaños de acuerdo al diámetro
de la columna y especialmente son de diversa geometría o forma como una
respuesta a la eficiencia de contacto o incremento del área interfacial de contacto
que puedan proveer éstos al contacto gas/líquido. Son muy conocidos los anillos
Rashig y las sillas o monturas Berl, y otros más eficientes como los anillos Pall,
las monturas Intalox y los anillos Bialecki entre otros.
En la selección y carga de un tipo de relleno se debe maximizar la eficiencia de
contacto que provocaría está a las fases circulantes, minimizar la caída de
presión por altura unitaria en la columna, ofrecer resistencia mecánica y térmica
en la operación unitaria y los aspectos de costo. (Arias, 2011, págs. 63-64)

Figura. 2
Columnas de relleno y sus diferentes formas

3. TORRES EMPACADAS PARA ABSORCION Y DESTILACIÓN

Las torres empacadas se usan para el contacto continuo a contracorriente

de un gas y un líquido y en absorción y también para el contacto de un vapor y
un líquido en la destilación. La torre de la figura siguiente consiste en una
columna cilíndrica que contiene una entrada de gas y un espacio de distribución
en el fondo, una entrada de líquido y un dispositivo de distribución en la parte
superior, la salida del gas en la parte superior, una salida de líquido en el fondo y
el empaque un relleno de la torre. El gas entra en el espacio de distribución que
está debajo de la reacción empacada y se va elevando a través de las aberturas o
intersticios de relleno, así se pone en contacto con el líquido descendiente que
fluye a través de las mismas aberturas. El empaque proporciona una extensa área
de contacto íntima entre el gas y el líquido.
 La característica de un buen empaque es la de tener una gran proporción
de espacios vacíos entre el orden del 60 y el 90%. El relleno permite que
volúmenes relativamente grandes de líquido pasen a contracorriente con
respecto al gas que fluye a través de las aberturas, con caídas de presión del gas
relativamente bajas. (Christie, 1998, pág. 681)

Figura. 3

Columna de relleno

III. MATERIALES Y EQUIPOS

MATERIALES
Columna de Relleno
Compresora
Termocupla

IV. PROCEDIMIENTO EXPERIMENTAL

 Instalar correctamente los dispositivos y conexiones hacia la columna de relleno.

 Carga del tanque de agua, operación de la bomba y del rotámetro.
 Aire comprimido y suministro hacia la columna.
 Nivelación y acondicionamiento de los piezómetros de agua (columna) y
manómetro de mercurio (inyección de aire).
 Determinar la temperatura del agua, el aire.
 Para iniciar la operación y evaluación de la columna, fijar primeramente el
caudal del agua hacia la columna, siendo la primera opción q agua = 0 (relleno
seco). Variar la inyección del aire hacia la columna desde un valor mínimo hasta
el máximo valor que admita el equipo.
  Efectuar la lectura de los piezómetros en cada variación de flujo de aire.
 Fijar en un nuevo caudal el flujo de agua y volver a variar el flujo de aire desde
un mínimo hasta el máximo permisible del equipo, (relleno mojado) reportar
nuevamente las variaciones de altura en los piezómetros. repetir la operación
hasta un máximo caudal de agua.
 El suministro de agua es a través de un rotámetro estándar en l/min.
 Mientras que la inyección del aire se realiza a través de un medidor de orificio,
cuya variación de presión se registra a través de un manómetro de mercurio.

V. CÁLCULOS

Para E. Rotámetro 0 L/min

Cálculos para la obtención de resultados

Calculamos área 1

𝐷1 = 16.2𝑚𝑚 = 0.00162𝑚

𝜋 2 𝜋
𝐴1 = 𝐷 = (0.0162𝑚)2 = 2.01062 × 10−4 𝑚2
4 1 4

Hallamos 𝜷

𝐷1 = 16.2𝑚𝑚 = 0.00162𝑚

𝐷2 = 7.5𝑚𝑚 = 7.5 ∗ 10−3 𝑚

0.0162𝑚
𝛽= = 2.16𝑚
7.5 ∗ 10−3 𝑚

Hallamos densidad de aire

𝑀𝑎𝑖𝑟𝑒 ∗ 𝑃
𝜌= = ⋯ … (1)
𝑅∗𝑇

𝑎𝑡𝑚∗𝑚3
Valor de 𝑅 = 0.082057 𝑘𝑚𝑜𝑙∗𝐾

Calculamos el peso molecular del aire

𝑘𝑔 𝑘𝑔
𝑀𝑎𝑖𝑟𝑒 = 𝑋𝑂2 ∗ 𝑃𝑀𝑂2 + 𝑋𝑁2 ∗ 𝑃𝑀𝑁2 = 0.21 ∗ 32 + 0.79 ∗ 28
𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙
𝑘𝑔 𝑘𝑔
= 28.84 = 29
𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙
Reemplazando en la ecuación (1)
𝑘𝑔 1𝑎𝑡𝑚
𝑀𝑎𝑖𝑟𝑒 ∗𝑃 29 ∗548𝑚𝑚𝐻𝑔∗
𝑘𝑚𝑜𝑙 670𝑚𝑚𝐻𝑔
𝜌𝑎𝑖𝑟𝑒 = 𝑅∗𝑇
= 𝑎𝑡𝑚∗𝑚3
= 0.8660 𝑘𝑔/𝑚3
0.082057 ∗294.15𝐾
𝑘𝑚𝑜𝑙∗𝐾

Densidad del mercurio en tablas

𝑘𝑔
𝜌𝐻𝑔 = 13600
𝑚3

Hallamos caudal

𝐶𝑜 𝐴1 𝜌𝐻𝑔 − 𝜌𝑎𝑖𝑟𝑒
𝑄𝑎𝑖𝑟𝑒 = √ 2𝑔 ∗ (∆ℎ𝐻𝑔 )[ ]
√𝛽4 − 1 𝜌𝑎𝑖𝑟𝑒

𝑘𝑔
0.61 ∗ 2.0106 × 10−4 𝑚2 √ 𝑚 13600 3 − 0.8660 𝑘𝑔/𝑚3
= 29.81 2 ∗ (0.008𝑚)[ 𝑚 ]
4
√2.16 − 1 𝑠 0.8660 𝑘𝑔/𝑚3

𝑚3 3600𝑠 𝑚3
= 1.3361 ∗ 10−3 × = 4.8101
𝑠 1ℎ ℎ

𝐿𝑛(𝑄𝑎𝑖𝑟𝑒 ) = 𝐿𝑛(4.8101) = 1.5707

Hallamos caída de presión en la columna

∆𝑃𝑐 = (∆ℎ𝐼 − ∆ℎ𝑠)𝑔𝜌𝑎𝑔𝑢𝑎 … … . . (2)

Datos:

𝑚
𝑔 = 9.81
𝑠2

∆ℎ𝐼 = ℎ4 − ℎ3 = 0.041𝑚

∆ℎ𝑠 = ℎ6 − ℎ5 = 0.039𝑚

Hallamos densidad del agua de tablas a 20.3 °C por interpolación

𝜌𝑎𝑔𝑢𝑎 = 998.1402 𝑘𝑔/𝑚3

Reemplazando en la ecuación 2

𝑚 𝑘𝑔
∆𝑃𝑐 = (∆ℎ𝐼 − ∆ℎ𝑠)𝑔𝜌𝑎𝑔𝑢𝑎 = (0.041𝑚 − 0.039𝑚)9.81 2
∗ 998.1402 3
𝑠 𝑚
= 19.5835 𝑃𝑎

Hallamos 𝑳𝒏∆𝑷𝒄

𝐿𝑛∆𝑃𝑐 = ln(19.5835 ) = 2.9747𝑃𝑎

Tabla 1
Resultados de la Columna de Relleno Seco
Datos del Co= 0.61 Taire(°C)= 21.1 Tagua(°C)= 20.3
Med. Orificio D1= 16.2mm=0.16m 𝜌𝑎𝑖𝑟𝑒 = 0.866 kg/m3 𝜌𝑎𝑔𝑢𝑎 = 998.14kg/m3
D2= 7.5mm=0.0075m 𝜌𝐻𝑔 = 13600 kg/m3
A1= 0.000201062m2
𝛽= 2.16
g= 9.81𝑚/𝑠 2
Qagua(L/h)
Med. Orif. (mmHg) Caudal del Aire Tope Inferior
h1 h2 Dh(mm) Dh(m) Q(m3/s) Q(m3/h) lnQ h3 h4 Dh (mm) Dh(m)
295 303 8 0.008 0.001336138 4.8100955 1.57071694 220 261 41 0.041
290 309 19 0.019 0.002059126 7.41285501 2.00321566 188 295 107 0.107
285 314 29 0.029 0.00254393 9.15814899 2.21464408 168 313 145 0.145
280 317 37 0.037 0.002873473 10.3445015 2.33645512 133 345 212 0.212
275 323 48 0.048 0.003272855 11.7822796 2.46659667 90 387 297 0.297
Tope
Superior Dpcol
h5 h6 Dh(mm) Dh(m) ∆𝑃𝑐 𝐿𝑛∆𝑃𝑐
195 234 39 0.039 19.5835107 2.97468792
165 264 99 0.099 78.3340429 4.36098228
148 279 131 0.131 137.084575 4.92059807
115 310 195 0.195 166.459841 5.11475409
70 349 279 0.279 176.251597 5.1719125
Realizar la gráfica logarítmica de la caída de presión de la columna frente al
caudal de aire:

Figura. 4.

Grafica 𝐿𝑛∆𝑃𝑐 VS 𝐿𝑛(𝑄𝑎 )

6
y = 2.5645x - 0.9239
5 R² = 0.964

4
Ln(Pcol)

0
0 0.5 1 1.5 2 2.5 3
LnQa

De la ecuación de la gráfica calculamos n

𝑦 = 2.5645𝑥 − 0.9239

𝑛 = 2.5645

Calculamos el porcentaje de error

Valor de n teórico

𝑛 = 1.8

1.8 − 2.56
%𝐸𝑟𝑟𝑜𝑟 = | | ∗ 100 = 42.2
1.8
 RELLENO MOJADO

Para E. Rotámetro 4 L/min

Cálculos para la obtención de resultados

Calculamos área 1

𝐷1 = 16.2𝑚𝑚 = 0.00162𝑚

𝜋 2 𝜋
𝐴1 = 𝐷 = (0.0162𝑚)2 = 2.01062 × 10−4 𝑚2
4 1 4

Hallamos 𝜷

𝐷1 = 16.2𝑚𝑚 = 0.00162𝑚

𝐷2 = 7.5𝑚𝑚 = 7.5 ∗ 10−3 𝑚

0.0162𝑚
𝛽= = 2.16𝑚
7.5 ∗ 10−3 𝑚

Hallamos densidad de aire

𝑀𝑎𝑖𝑟𝑒 ∗ 𝑃
𝜌= = ⋯ … (1)
𝑅∗𝑇

𝑎𝑡𝑚∗𝑚3
Valor de 𝑅 = 0.082057 𝑘𝑚𝑜𝑙∗𝐾

Calculamos el peso molecular del aire

𝑘𝑔 𝑘𝑔
𝑀𝑎𝑖𝑟𝑒 = 𝑋𝑂2 ∗ 𝑃𝑀𝑂2 + 𝑋𝑁2 ∗ 𝑃𝑀𝑁2 = 0.21 ∗ 32 + 0.79 ∗ 28
𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙
𝑘𝑔 𝑘𝑔
= 28.84 = 29
𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙

Reemplazando en la ecuación (1)

𝑘𝑔 1𝑎𝑡𝑚
𝑀𝑎𝑖𝑟𝑒 ∗𝑃 29 ∗548𝑚𝑚𝐻𝑔∗
𝑘𝑚𝑜𝑙 670𝑚𝑚𝐻𝑔
𝜌𝑎𝑖𝑟𝑒 = = 𝑎𝑡𝑚∗𝑚3
= 0.8660 𝑘𝑔/𝑚3
𝑅∗𝑇 0.082057 ∗294.15𝐾
𝑘𝑚𝑜𝑙∗𝐾

Densidad del mercurio en tablas

𝑘𝑔
𝜌𝐻𝑔 = 13600
𝑚3

Hallamos caudal
 𝐶𝑜 𝐴1 𝜌𝐻𝑔 − 𝜌𝑎𝑖𝑟𝑒
𝑄𝑎𝑖𝑟𝑒 = √ 2𝑔 ∗ (∆ℎ𝐻𝑔 )[ ]
√𝛽4 − 1 𝜌𝑎𝑖𝑟𝑒

𝑘𝑔
0.61 ∗ 2.0106 × 10−4 𝑚2 √ 𝑚 13600 3 − 0.8660 𝑘𝑔/𝑚3
= 29.81 2 ∗ (0.008𝑚)[ 𝑚 ]
4
√2.16 − 1 𝑠 0.8660 𝑘𝑔/𝑚3

𝑚3 3600𝑠 𝑚3
= 1.3361 ∗ 10−3 × = 4.8101
𝑠 1ℎ ℎ

𝐿𝑛(𝑄𝑎𝑖𝑟𝑒 ) = 𝐿𝑛(4.8101) = 1.5707

Hallamos caída de presión en la columna

∆𝑃𝑐 = (∆ℎ𝐼 − ∆ℎ𝑠)𝑔𝜌𝑎𝑔𝑢𝑎 … … . . (2)

Datos:

𝑚
𝑔 = 9.81
𝑠2

∆ℎ𝐼 = ℎ4 − ℎ3 = 0.043𝑚

∆ℎ𝑠 = ℎ6 − ℎ5 = 0.033𝑚

Hallamos densidad del agua de tablas a 20.3 °C

𝜌𝑎𝑔𝑢𝑎 = 998.1402 𝑘𝑔/𝑚3

Reemplazando en la ecuación 2

𝑚 𝑘𝑔
∆𝑃𝑐 = (∆ℎ𝐼 − ∆ℎ𝑠)𝑔𝜌𝑎𝑔𝑢𝑎 = (0.043𝑚 − 0.033𝑚)9.81 2
∗ 998.1402 3
𝑠 𝑚
= 97.9176 𝑃𝑎

Hallamos 𝐿𝑛∆𝑃𝑐

𝐿𝑛∆𝑃𝑐 = ln(97.9176 ) = 4.5841𝑃𝑎

TABLA 2
Resultados de la Evaluación de columna de relleno en húmedo
Datos del Co= 0.61 Taire(°C)= 21.1 Tagua(°C)= 20.3
998.1402
Med. Orificio D1= 16.2mm=0.16m 0.8660 kg/m3 𝜌𝑎𝑔𝑢𝑎 = kg/m3
D2= 7.5mm=0.0075m 𝜌𝐻𝑔 13600 kg/m3
A1= 0.000201062m2
B= 2.16
g= 9.81m/s2
Qagua(L/h)
Caudal del Aire
Med. Orif. (mmHg) Tope Inferior
h1 h2 Dh(mm) Dh(m) Q(m3/s) Q(m3/h) lnQ h3 h4 Dh (mm) Dh(m)
295 303 8 0.008 0.00133614 4.8100955 1.57071694 218 261 43 0.043
290 309 19 0.019 0.00205913 7.41285501 2.00321566 172 309 137 0.137
285 314 29 0.029 0.00254393 9.15814899 2.21464408 139 342 203 0.203
280 317 37 0.037 0.00287347 10.3445015 2.33645512 118 362 244 0.244
275 323 48 0.048 0.00327286 11.7822796 2.46659667 49 429 380 0.38
Tope Superior Dpcol
h5 h6 Dh(mm) Dh(m) Pa lnpa
198 231 33 0.033 97.9175536 4.58412583
159 268 109 0.109 274.16915 5.61374525
133 295 162 0.162 401.46197 5.99511281
116 310 194 0.194 489.587768 6.19356375
59 365 306 0.306 724.589897 6.58560584
Realizar la gráfica logarítmica de la caída de presión de la columna frente al
caudal de aire:

Figura. 5.

Grafica 𝐿𝑛∆𝑃𝑐 VS 𝐿𝑛(𝑄𝑎 )

5
Ln(Pcol)

4
y = 2.1689x + 1.2001
3 R² = 0.9948
2

0
0 0.5 1 1.5 2 2.5 3
ln(Qa)

De la ecuación de la gráfica calculamos n

𝑦 = 2.1689𝑥 + 1.2001

𝑛 = 2.1689

Hallamos el porcentaje de error

Valor de n teórico

𝑛 = 1.8

1.8 − 2.17
%𝐸𝑟𝑟𝑜𝑟 = | | ∗ 100 = 20.56
1.8
 Para E. Rotámetro 6 L/min

Cálculos para la obtención de resultados

Calculamos área 1

𝐷1 = 16.2𝑚𝑚 = 0.00162𝑚

𝜋 2 𝜋
𝐴1 = 𝐷 = (0.0162𝑚)2 = 2.01062 × 10−4 𝑚2
4 1 4

Hallamos 𝜷

𝐷1 = 16.2𝑚𝑚 = 0.00162𝑚

𝐷2 = 7.5𝑚𝑚 = 7.5 ∗ 10−3 𝑚

0.0162𝑚
𝛽= = 2.16𝑚
7.5 ∗ 10−3 𝑚

Hallamos densidad de aire

𝑀𝑎𝑖𝑟𝑒 ∗ 𝑃
𝜌= = ⋯ … (1)
𝑅∗𝑇

𝑎𝑡𝑚∗𝑚3
Valor de 𝑅 = 0.082057 𝑘𝑚𝑜𝑙∗𝐾

Calculamos el peso molecular del aire

𝑘𝑔 𝑘𝑔
𝑀𝑎𝑖𝑟𝑒 = 𝑋𝑂2 ∗ 𝑃𝑀𝑂2 + 𝑋𝑁2 ∗ 𝑃𝑀𝑁2 = 0.21 ∗ 32 + 0.79 ∗ 28
𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙
𝑘𝑔 𝑘𝑔
= 28.84 = 29
𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙

Reemplazando en la ecuación (1)

𝑘𝑔 1𝑎𝑡𝑚
𝑀𝑎𝑖𝑟𝑒 ∗𝑃 29 ∗548𝑚𝑚𝐻𝑔∗
𝑘𝑚𝑜𝑙 670𝑚𝑚𝐻𝑔
𝜌𝑎𝑖𝑟𝑒 = = 𝑎𝑡𝑚∗𝑚3
= 0.8660 𝑘𝑔/𝑚3
𝑅∗𝑇 0.082057 ∗294.15𝐾
𝑘𝑚𝑜𝑙∗𝐾

Densidad del mercurio en tablas

𝑘𝑔
𝜌𝐻𝑔 = 13600
𝑚3
Hallamos caudal

𝐶𝑜 𝐴1 𝜌𝐻𝑔 − 𝜌𝑎𝑖𝑟𝑒
𝑄𝑎𝑖𝑟𝑒 = √ 2𝑔 ∗ (∆ℎ𝐻𝑔 )[ ]
√𝛽4 − 1 𝜌𝑎𝑖𝑟𝑒

𝑘𝑔
0.61 ∗ 2.0106 × 10−4 𝑚2 √ 𝑚 13600 3 − 0.8660 𝑘𝑔/𝑚3
= 29.81 2 ∗ (0.008𝑚)[ 𝑚 ]
√2.164 − 1 𝑠 0.8660 𝑘𝑔/𝑚3

−3
𝑚3 3600𝑠 𝑚3
= 1.3361 ∗ 10 × = 4.8101
𝑠 1ℎ ℎ

𝐿𝑛(𝑄𝑎𝑖𝑟𝑒 ) = 𝐿𝑛(4.8101) = 1.5707

Hallamos caída de presión en la columna

∆𝑃𝑐 = (∆ℎ𝐼 − ∆ℎ𝑠)𝑔𝜌𝑎𝑔𝑢𝑎 … … . . (2)

Datos:

𝑚
𝑔 = 9.81
𝑠2

∆ℎ𝐼 = ℎ4 − ℎ3 = 0.035𝑚

∆ℎ𝑠 = ℎ6 − ℎ5 = 0.027𝑚

Hallamos densidad del agua de tablas a 20.3 °C interpolando

𝜌𝑎𝑔𝑢𝑎 = 998.1402 𝑘𝑔/𝑚3

Reemplazando en la ecuación 2

𝑚 𝑘𝑔
∆𝑃𝑐 = (∆ℎ𝐼 − ∆ℎ𝑠)𝑔𝜌𝑎𝑔𝑢𝑎 = (0.035𝑚 − 0.027𝑚)9.81 2
∗ 998.1402 3
𝑠 𝑚
= 78.3340 𝑃𝑎

Hallamos 𝑳𝒏∆𝑷𝒄

𝐿𝑛∆𝑃𝑐 = ln(97.9176 ) = 4.3610𝑃𝑎

TABLA 3
Resultados de la Evaluación de columna de relleno en húmedo
Datos del Co= 0.61 Taire(°C)= 21.1 Tagua(°C)= 20.3
Med. Orificio D1= 16.2mm=0.16m 𝜌𝑎𝑖𝑟𝑒 0.8660 𝑘𝑔/𝑚3 𝜌𝑎𝑔𝑢𝑎 = 998.1402𝑘𝑔/𝑚3

D2= 7.5mm=0.0075m 𝜌𝐻𝑔 13600𝑘𝑔/𝑚3

A1= 0.00020106m2
𝛽 2.16
g= 9.81m/s2
Qagua(L/h)
Med. Orif.
(mmHg) Caudal del Aire Tope Inferior
h1 h2 Dh(mm) Dh(m) Q(m3/s) Q(m3/h) lnQ h3 h4 Dh (mm) Dh(m)
295 303 8 0.008 0.00133614 4.8100955 1.57071694 223 258 35 0.035
290 309 19 0.019 0.00205913 7.41285501 2.00321566 167 314 147 0.147
285 314 29 0.029 0.00254393 9.15814899 2.21464408 153 327 174 0.174
280 317 37 0.037 0.00287347 10.3445015 2.33645512 114 365 251 0.251
275 323 48 0.048 0.00327286 11.7822796 2.46659667 33 444 411 0.411
Tope Superior
Dpcol
h5 h6 Dh(mm) Dh(m) Pa lnpa
200 227 27 0.027 78.3340429 4.36098228
158 269 111 0.111 352.503193 5.86505968
149 278 129 0.129 440.628991 6.08820323
118 307 189 0.189 607.088832 6.40867513
57 367 310 0.31 988.967292 6.89666126
 Realizar la gráfica logarítmica de la caída de presión de la columna frente al
caudal de aire:

Figura 6.

Grafica 𝐿𝑛∆𝑃𝑐𝑜𝑙 VS 𝐿𝑛(𝑄𝑎 )

8
7 y = 2.693x + 0.2193
R² = 0.9749
6
Ln(Pcol)

5
4
3
2
1
0
0 0.5 1 1.5 2 2.5 3
LnQa

De la ecuación de la gráfica calculamos n

𝒚 = 𝟐. 𝟔𝟗𝟑𝒙 + 𝟎. 𝟐𝟏𝟗𝟑

𝑛 = 2.693

Hallamos el porcentaje de error

Valor de n teórico

𝑛 = 1.8

1.8 − 2.693
%𝐸𝑟𝑟𝑜𝑟 = | | ∗ 100 = 49.61
1.8
 Para E. Rotámetro 8 L/min

Cálculos para la obtención de resultados

Hallamos Área:

 𝐷1 = 16.2 𝑚𝑚 = 0.0162 𝑚

𝜋
𝐴= 0.01622 𝑚 → 𝐴 = 2.0612 ∗ 10−4 𝑚2
4

Hallamos 𝜷:

𝐷1 = 16.2 𝑚𝑚 = 0.0162 𝑚

𝐷2 = 7.5 𝑚𝑚 = 0.0075 𝑚

𝐷1 0.0162 𝑚
𝛽= = = 2.16
𝐷2 0.0075 𝑚

Hallamos la densidad del aire:

 𝑀𝑎𝑖𝑟𝑒 = 28.96 𝑘𝑔/𝑚3 ≅ 29 𝑘𝑔/𝑚3

 𝑃 = 548 𝑚𝑚𝐻𝑔 ≅ 0.7211 𝑎𝑡𝑚
 𝑇 = 21.1 °𝐶 ≅ 294.15 𝐾
 𝑅 = 0.082057 𝑎𝑡𝑚. 𝑚3 /𝐾. 𝐾𝑚𝑜𝑙

𝑀∗𝑃 29 𝑘𝑔/𝑚3 ∗ 0.7211 𝑎𝑡𝑚

𝜌𝑎𝑖𝑟𝑒 = = = 0.8664 𝑘𝑔/𝑚3
𝑅 ∗ 𝑇 0.082057 𝑎𝑡𝑚. 𝑚3 /𝐾. 𝐾𝑚𝑜𝑙 ∗ 294.15 𝐾

Hallamos Q:

Densidad de mercurio de tablas

 𝜌𝑚𝑒𝑟𝑐𝑢𝑟𝑖𝑜 = 13600 𝑘𝑔/𝑚3

  Hg   aire 
Q
C0 A1
2 g h 
Hg  
 4 1   aire 

0.61 ∗ 2.0612 ∗ 10−4 𝑚2 𝑚 13600 𝑘𝑔/𝑚3 − 0.8664 𝑘𝑔/𝑚3

= √2 ∗ 9.81 ∗ 0.008 𝑚 ∗ [ ]
√2.164 − 1 𝑠2 0.8664 𝑘𝑔/𝑚3
 𝑄 = 1.3695 ∗ 10−3 𝑚3 /𝑠

CAÍDA DE PRESIÓN EN LA COLUMNA

PC  PI  PS  hI  hS  g agua

 𝑔 = 9.81 𝑚/𝑠 2
 𝑇 = 20.3 °𝐶

 De tablas 𝜌𝑎𝑔𝑢𝑎 = 998.1402 𝑘𝑔/𝑚3

𝑚
∆𝑃𝐶 = (0.259 − 0.221)𝑚 − (0.227 − 0.201)) ∗ 9.81 ∗ 998.1402 𝑘𝑔/𝑚3
𝑠2

∆𝑃𝐶 = 117.5011 𝑘𝑔/𝑚. 𝑠 2 ≅ 117.50117𝑃𝑎

TABLA 4
Resultados de la Evaluación de Relleno húmedo

Datos del Co= 0.61 Taire(°C)= 21.1 Tagua(°C)= 20.3

Med. Orificio D1= 16.2mm=0.16m 𝜌𝑎𝑖𝑟𝑒 = 0.866 kg/m3 𝜌𝑎𝑔𝑢𝑎 = 998.1402kg/m3
D2= 7.5mm=0.0075m 𝜌𝐻𝑔 = 13600 kg/m3
A1= 0.00020612m2 2.0612 ∗ 10−4
𝛽= 2.16
g= 9.81𝑚/𝑠 2
Qagua(L/h)
Med. Orif. Caudal del
(mmHg) Aire Tope Inferior
h1 h2 Dh(mm) Dh(m) Q(m3/s) Q(m3/h) lnQ h3 h4 Dh (mm) Dh(m)
295 303 8 0.008 0.0013695 4.930073 1.5954 221 259 38 0.038
290 309 19 0.019 0.0021105 7.597753 2.0279 161 320 159 0.159
285 314 29 0.029 0.0026074 9.386580 2.2393 135 344 209 0.209
280 317 37 0.037 0.0029451 10.602524 2.3611 90 389 299 0.299
275 323 48 0.048 0.0033545 12.076164 2.4912 48 427 379 0.379
Tope
Superior Dpcol
h5 h6 Dh(mm) Dh(m) ∆𝑃𝑐 𝐿𝑛∆𝑃𝑐
201 227 26 0.026 117.5011 4.7664
157 269 112 0.112 460.2125 6.1317
140 286 146 0.146 616.8806 6.4247
108 318 210 0.21 871.4662 6.7702
83 343 260 0.26 1165.2189 7.0607
Realizar la gráfica logarítmica de la caída de presión de la columna frente al caudal
de aire:

Figura 7.

Grafica 𝐿𝑛∆𝑃𝑐𝑜𝑙 VS 𝐿𝑛(𝑄𝑎 )

Ln Q versus Ln ΔPc
8.0000
y = 2.5213x + 0.8277
7.0000
R² = 0.9848
6.0000

5.0000

4.0000

3.0000

2.0000

1.0000

0.0000
0.0000 0.5000 1.0000 1.5000 2.0000 2.5000 3.0000

De la ecuación de la gráfica calculamos n

𝑦 = 2.5213𝑥 + 0.8277

𝑛 = 2.5213

Calculamos el porcentaje de error

Valor de n teórico

𝑛 = 1.8

1.8 − 2.5213
%𝐸𝑟𝑟𝑜𝑟 = | | ∗ 100 = 40.0722
1.8
 Para E. Rotámetro 10 L/min

Cálculos para la obtención de resultados

Hallamos Área:

 𝐷1 = 16.2 𝑚𝑚 = 0.0162 𝑚

𝜋
𝐴= 0.01622 𝑚 → 𝐴 = 2.0612 ∗ 10−4 𝑚2
4

Hallamos 𝜷:

𝐷1 = 16.2 𝑚𝑚 = 0.0162 𝑚

𝐷2 = 7.5 𝑚𝑚 = 0.0075 𝑚

𝐷1 0.0162 𝑚
𝛽= = = 2.16
𝐷2 0.0075 𝑚

Hallamos la densidad del aire:

 𝑀𝑎𝑖𝑟𝑒 = 28.96 𝑘𝑔/𝑚3 ≅ 29 𝑘𝑔/𝑚3

 𝑃 = 548 𝑚𝑚𝐻𝑔 ≅ 0.7211 𝑎𝑡𝑚
 𝑇 = 21.1 °𝐶 ≅ 294.15 𝐾
 𝑅 = 0.082057 𝑎𝑡𝑚. 𝑚3 /𝐾. 𝐾𝑚𝑜𝑙

𝑀∗𝑃 29 𝑘𝑔/𝑚3 ∗ 0.7211 𝑎𝑡𝑚

𝜌𝑎𝑖𝑟𝑒 = = = 0.8664 𝑘𝑔/𝑚3
𝑅 ∗ 𝑇 0.082057 𝑎𝑡𝑚. 𝑚3 /𝐾. 𝐾𝑚𝑜𝑙 ∗ 294.15 𝐾

Hallamos Q:

 𝜌𝑚𝑒𝑟𝑐𝑢𝑟𝑖𝑜 = 13600 𝑘𝑔/𝑚3

  Hg   aire 
Q
C0 A1
2 g h 
Hg  
 4 1   aire 

0.61 ∗ 2.0612 ∗ 10−4 𝑚2 𝑚 13600 𝑘𝑔/𝑚3 − 0.8664 𝑘𝑔/𝑚3

= √2 ∗ 9.81 ∗ 0.008 𝑚 ∗ [ ]
√2.164 − 1 𝑠2 0.8664 𝑘𝑔/𝑚3
 𝑄 = 1.3695 ∗ 10−3 𝑚3 /𝑠

CAÍDA DE PRESIÓN EN LA COLUMNA

PC  PI  PS  hI  hS  g agua

 𝑔 = 9.81 𝑚/𝑠 2
 𝑇 = 20.3 °𝐶

 De tablas 𝜌𝑎𝑔𝑢𝑎 = 998.1402 𝑘𝑔/𝑚3

𝑚
∆𝑃𝐶 = (0.264 − 0.210)𝑚 − (0.228 − 0.198)) ∗ 9.81 2
∗ 998.1402 𝑘𝑔/𝑚3
𝑠

∆𝑃𝐶 = 235.0021 𝑘𝑔/𝑚. 𝑠 2 ≅ 235.0021𝑃𝑎

TABLA 5
Resultados de la Evaluación de Relleno húmedo
Datos del Co= 0.61 Taire(°C)= 21.1 Tagua(°C)= 20.3
Med. Orificio D1= 16.2mm=0.16m 𝜌𝑎𝑖𝑟𝑒 = 0.866 kg/m3 𝜌𝑎𝑔𝑢𝑎 = 998.14kg/m3
D2= 7.5mm=0.0075m 𝜌𝐻𝑔 = 13600 kg/m3
A1= 0.00020612m2
𝛽= 2.16
g= 9.81𝑚/𝑠 2
Qagua(L/h)
Med. Orif. Caudal del
(mmHg) Aire Tope Inferior
h1 h2 Dh(mm) Dh(m) Q(m3/s) Q(m3/h) lnQ h3 h4 Dh (mm) Dh(m)
295 303 8 0.008 0.0013695 4.930073 1.5954 210 264 54 0.054
290 309 19 0.019 0.0021105 7.597753 2.0279 170 309 139 0.139
285 314 29 0.029 0.0026074 9.386580 2.2393 120 360 240 0.240

Tope
Superior Dpcol
h5 h6 Dh(mm) Dh(m) ∆𝑃𝑐 𝐿𝑛∆𝑃𝑐
198 228 30 0.030 235.0021 5.4596
172 254 82 0.082 558.1301 6.3246
114 282 168 0.168 705.0064 6.5582
Realizar la gráfica logarítmica de la caída de presión de la columna frente al caudal
de aire:

Figura 8.

Grafica 𝐿𝑛∆𝑃𝑐𝑜𝑙 VS 𝐿𝑛(𝑄𝑎 )

Ln Q versus Ln ΔPc
7.0000
y = 1.7496x + 2.6952
6.0000 R² = 0.9845
5.0000

4.0000

3.0000

2.0000

1.0000

0.0000
0.0000 0.5000 1.0000 1.5000 2.0000 2.5000

De la ecuación de la gráfica calculamos n

𝑦 = 1.7496𝑥 + 2.6952

𝑛 = 1.7496

Calculamos el porcentaje de error

Valor de n teórico

𝑛 = 1.8

1.8 − 1.7496
%𝐸𝑟𝑟𝑜𝑟 = | | ∗ 100 = 2.8
1.8
CALCULO DEL CAUDAL ÓPTIMO PARA CADA CASO

 Para el caudal de la Figura 5:

7
punto de inundacion
6
punto de carga
5
Ln(Pcol)

0
0 0.5 1 1.5 2valor optimo 2.5 3
ln(Qa)

Valor en el punto de carga: LnQa=2.003215657

Valor en el punto de inundación: LnQa=2.336455124

Con estos dos puntos podemos hallar el caudal óptimo en LnQa haciendo un
promedio que resulta 2.1698
𝐿𝑛𝑄𝑎 = 2.1698
𝑄𝑎 = 𝑒 2.1698 = 8.7565
 Para el caudal de la Figura 6:

8
7
6
Ln(Pcol)

5
4
3
2
1
0
0 0.5 1 1.5 2 2.5 3
LnQa

Valor en el punto de carga: LnQa=2.003215657

Valor en el punto de inundación: LnQa=2.214644082

Con estos dos puntos podemos hallar el caudal óptimo en LnQa haciendo un
promedio que resulta 2.1089
𝐿𝑛𝑄𝑎 = 2.1089
𝑄𝑎 = 𝑒 2.1698 = 8.2392
 Para el caudal de la Figura 7:

Ln Q versus Ln ΔPc
8.0000

7.0000

6.0000

5.0000

4.0000

3.0000

2.0000

1.0000

0.0000
0.0000 0.5000 1.0000 1.5000 2.0000 2.5000 3.0000

Valor en el punto de carga: LnQa=2.0279

Valor en el punto de inundación: LnQa=2.2393

Con estos dos puntos podemos hallar el caudal óptimo en LnQa haciendo un
promedio que resulta 2.1336
𝐿𝑛𝑄𝑎 = 2.1336
𝑄𝑎 = 𝑒 2.1336 = 8.4452
GRAFICO COMPARATIVO
Este grafico comparativo mide las tendencias de cada una de los caudales de aire como para
Q=0, 4,6,8 y 10 L/m respectivamente.

Figura.9

Grafica de Ln(△pcol) vs LnQa con respecto s cada caudal

Ln (ΔPcol) versus Ln Qa
3

2.5

2 caudal 0 L/h
caudal 240 L/h
1.5
caudal 360 L/h

1 caudal 480 L/h

caudal 600 L/h
0.5

0
0 2 4 6 8

VI. DISCUSIONES
 Las columnas o torres de separación son los equipos más ampliamente utilizadas en
la transferencia de masa, ya sea en operaciones como: destilación batch con
rectificación, destilación continua, destilación reactiva, absorción de gases,
desorción de gases, humidificadores adiabáticos, extracción líquida, etc. Son
dispositivos que facilitan el contacto continuo entre las fases que circulan por su
interior, incrementando el área de transferencia de masa interfacial gracias a los
dispositivos de contacto interno instalados y promoviendo la separación de
componentes.
  Al realizar los cálculos y graficar 𝐿𝑛∆𝑃𝑐 VS 𝐿𝑛(𝑄𝑎𝑖𝑟𝑒 ) en la columna de relleno
seco la ecuación de la gráfica nos da el valor de la pendiente que es igual al valor de
n=2.5645 que es mayor al valor teórico de n=1.8 lo que nos quiere decir que le falta
relleno en nuestra columna, sacando un porcentaje de error del 42.2 %.

VII. CONCLUSIONES
 Se estudió el fenómeno producido por los caudales agua y aire (Sistema bifásico) en
el interior de la columna en contracorriente, para dicho experimento se dejó fluir
agua con caudal conocida y el aire con caudal variable.
 En las columnas de relleno en húmedo al igual que la columna de relleno en seco
los valores de n son mayores al valor teórico, excepto el experimento 5 que el n se
acerca al valor teórico.
 Como los porcentajes de error son altos interpretamos que a la columna se debe
aumentar más relleno.
 En la FIG.8 para el diámetro tenemos demasiada agua por lo que sufre una
inundación
 Determinamos los parámetros de flujo y localizamos el punto de carga e inundación
de la columna.

VIII. CUESTIONARIO
1. Explorar de que factores depende el dimensionamiento de una columna.
Los factores que dependen el dimensionamiento de una columna son:
 Uno de los factores son la concentración o recuperación de soluto que uno
quiere obtener
 Los tipos de rellenos
 El área del relleno dependiendo del caudal
IX. BIBLIOGRAFÍA
1. Christie J. Geankoplis, (1998). Procesos de Transporte y Operaciones Unitarias.
México: Compañia Editorial Continental, S.A. de C.V. México.

2. Alfredo Arias Jara (2011). Fundamentos y Aplicación de Transferencia de Masa.

Ayacucho: PRAKTICO L&C S.R.L.

ANEXO 1

PESO MOLECULAR DEL OXIGENO Y NITROGENO

𝑘𝑔
𝑃𝑀𝑂2 = 32
𝑘𝑚𝑜𝑙

𝑘𝑔
𝑃𝑀𝑁2 = 28
𝑘𝑚𝑜𝑙

Fuente: Arias Jara, Alfredo; Fundamentos y aplicaciones de transferencia

de masa (pág 411)

ANEXO 2

FUENTE: Arias Jara, Alfredo; Fundamentos y aplicaciones de transferencia de masa (pág

406)