Universidad de Costa Rica

Ciudad Universitaria Rodrigo Facio

Facultad de Ingeniería
Escuela de Ingeniería Química

IQ0534. Laboratorio de Operaciones de Separación por métodos

difusionales

Datos experimentales, intermedios y cálculos de pared

húmeda

Profesora: Ing. Dianella Garro

Asistente: Paúl Solano

Daniel Matamoros Arnáez (B54204 )

Walter Nolasco (B55069 )
Sebastian Zúñiga Rojas (B59104 )

San José, Costa Rica, 2019.

MARCO TEÓRICO

En la ingeniería química, es muy importante el entendimiento de diferentes

procesos químicos y termodinámicos para un mejor aprovechamiento de recursos y
para lograr mayor eficiencia en los resultados requeridos. Por esta razón, el entender
y conocer acerca de la transferencia de masa en equipos como las torres de pared
húmeda es esencial en los profesionales en esta área.

2.1. Modelos para el cálculo de coeficientes de transferencia

de masa

Existen diferentes modelos para expresar el cálculo de los coeficientes de trans-

ferencia de masa entre los que están (McCabe, 2007):

Modelo de película
La teoría de la película permite aproximar la resistencia a la difusión como
una delgada capa de determinado espesor. Frecuentemente este modelo se usa
para modelar problemas complejos de transferencia de masa en reacciones
químicas. Este expresa el coeficiente de transferencia de masa de la siguiente
manera (McCabe, 2007):
DAB
kc = (1)
δ

donde,

kc = coeficiente de transferencia de masa, m/s

m2
DAB = difusividad, s

δ = espesor de película, m

No obstante, la ecuación (1) aplica solamente bajo las suposiciones de que en

la ley de Fick el cambio en la concentración es poco, la fracción molar de la
sustancia no es cercana a uno y que el largo de la película es mucho mayor que
su espesor (Izquierdo y Torres, 2004).

1
2

Modelo de pared mojada

En este modelo el soluto no debe penetrar demasiado la película de líquido.
El coeficiente de transferencia de masa se aproxima como (Cengel, 2007):

r
4 ∗ DAB ∗ v∞
kc = (2)
π∗L
donde,

kc = coeficiente de transferencia de masa, m/s

m2
DAB = difusividad, s

L = longitud, m

v∞ = velocidad en el seno del fluido, m/s

Modelo de capa límite

Este modelo es la mejor forma de aproximar un soluto a otro que se encuen-
tra bajo las mismas condiciones hidrodinámicas, sin embargo, si el régimen
de flujo se vuelve turbulento o se da una separación evidente, esta funciona
solamente como guía para la formulación de correlaciones empíricas. Las rela-
ciones usadas para transferencia de calor pueden usarse para transferencia de
masa simplemente reemplazando el número de Nusselt por el número de Sher-
wood y el número de Prandtl por el número de Schmidt gracias a la analogía
entre transferencia de masa y transferencia de calor. El modelo es el siguiente
(Cengel, 2007):

DAB p
kc = 0,664 ∗ ∗ Rel ∗ Sc1/3 (3)
L
donde,

kc = coeficiente de transferencia de masa, m/s

m2
DAB = difusividad, s

L = longitud, m

Rel = número de Reynolds del líquido, adim.

 3

Sc = número de Schmidt, adim.

La capa límite abarca desde la superficie sobre la que se mueve el flujo hasta
donde se alcanza el 99 % de la velocidad del flujo como si estuviera libre.
Su espesor no es constante ya que existe una transición de flujo laminar a
turbulento conforme se aumente la distancia X horizontal desde donde inicia el
flujo, aunque siempre existe una pequeña subcapa laminar donde el gradiente
de velocidad es muy alto. Esto se puede observar en la Figura 2.1 (Welty,
2008).

Figura 2.1. Diagrama que muestra la variación del espesor de la capa límite (Welty,
2008)

Modelo de la doble resistencia

Existen varios casos donde la materia debe difundirse entre fases distintas y
la velocidad global de transferencia de masa se ve afectada por las velocidades
locales de difusión en ambas fases. La teoría supone que en la interfase existe
un equilibrio donde la resistencia a la transferencia de materia en ambas fases
se convierten en una sola. El modelo es el siguiente (McCabe, 2007):

1 1 m
= + (4)
Ky ky kx

Ky = coeficiente global de transferencia de masa, m/s

ky = coeficiente local de transferencia de masa de fase gas, m/s

4

kx = coeficiente local de transferencia de masa del líquido, m/s

m = pendiente de la recta de equilibrio, adimensional

2.2. Variables que influyen en la transferencia de masa

La transferencia de masa se puede ver influenciada por diferentes variables ya

sean termodinámicas o propias de las sustancias implicadas en este fenómeno. Ade-
más dependiendo de las fases en las que se encuentren las sustancias que participan
en la transferencia de masa, las variables termodinámicas pueden afectar de diferente
manera.
La Ley de Fick, la cual es por la que se rige la transferencia de masa en la
mayoría de los casos, en la siguiente ecuación (Welty, 2008):

JA = −DAB ∇CA (5)

donde,

JA = Flux de molar de A de transferencia de masa, mol/m2 s.

DAB = Coeficiente de difusividad de la sustancia A en la sustancia B, m2 /s.

CA = Concentración de A, mol/m3 .

Para cuando se tiene transferencia de masa únicamente en la dirección z aplica:

∂CA
JA = −DAB (6)
∂z
donde,

JA = Flux de A molar, mol/m2 s.

DAB = Coeficiente de difusividad de la sustancia A en la sustancia B, m2 /s.

CA = Concentración de A, mol/m3 .

z= distancia en la dirección de la transferencia de masa, m.

 5

Se observa que entre mayor sea la diferencia de concentración entre las dos fases
que intervengan en la transferencia de masa, mayor será la transferencia de masa.
Además entre mayor sea el coeficiente de difusividad la transferencia de masa se ve
favorecida. A este coeficiente lo afectan diversos factores dependiendo de si se da la
transferencia de masa en líquido o en gas. Para gases se observa en el modelo de
Hirschfelder descrito en la siguiente ecuación:

h i1/2
1 1
0,001858T 3/2 MA
+ MB
DAB = 2
(7)
P σAB ΩD
donde,

DAB = Coeficiente de difusividad de la sustancia A en la sustancia B, m2 /s.

T = Temperatura, K.

P = Presión, kPa.

MA = Peso molecular de la sustancia A, g/mol.

MB = Peso molecular de la sustancia B, g/mol.

σAB = Diámetro de colisión, .

ΩD = Integral de colisión, adim.

Al observar esta ecuación se puede determinar que variables termodinámicas

como la temperatura y la presión influyen en el coeficiente de difusividad y por ende
en la transferencia de masa. Al aumentar la temperatura se favorece el fenómeno en
estudio mientras que al aumentar la presión se ve perjudicado. Además, los pesos
moleculares, así como el diámetro de colisión afectan a la transferencia de masa, los
cuales son independientes de algún control pero se puede ver que al utilizar distintas
sustancias, la transferencia de masa es distinta.

2.3. Teoría de película

Esta teoría fue originalmente propuesta en 1904 por Walther Nernst, utilizada
para explicar la transferencia de masa. Esta dice que cercano a donde un fluido
6

interactúa con una frontera, hay una pequeña capa con un regimen laminar que
tiene una resistencia a esta transferencia con un valor igual al que tendría todo el
fluido si se estuviera moviendo. Esto es igual a suponer que la resistencia del proceso
está en esta capa donde ocurre difusión. En el caso de un flujo turbulento, la película
laminar comprendería la subcapa laminar y parte de la capa de amortiguamiento en
un flujo con un núcleo turbulento (Suastegui, 1997).

Dentro de este modelo se encuentra también el de la película estancada, el cual

considera que el fluido tiene dos zonas principales, siendo una la capa laminar y la
otra una película con agitación perfecta donde no hay gradientes de concentración.
Además, se considera que la transferencia es realizada en estado estacionario. La
velocidad local sobre el límite de la frontera del fluido se aproxima a cero, por lo que
la transferencia se realiza por difusión molecular. La Figura 2.2 muestra el cambio
de la concentración de una sustancia de acuerdo al espesor de la película (Izquierdo
y Torres, 2004).

Figura 2.2. Diagrama que muestra la variación de la concentración de acuerdo al espe-

sor de la película laminar (Izquierdo y Torres, 2004)

Esta teoría no es muy utilizada por el hecho de que el espesor de la película no se

puede medir ya que es ficticio (Izquierdo y Torres, 2004) y se realiza la implicación de
que el coeficiente de transferencia de masa es proporcional a la difusividad mediante
una relación lineal, lo cual no siempre se cumple (Suastegui, 1997).
 7

2.4. Teoría de penetración

La teoría de penetración fue propuesta por Higbie en 1935. Este sugirió que en
ciertos fenómenos como la absorción de gases ocurren contactos entre las fases, a
lo que luego ocurre un proceso de mezcla entre estas. La Figura 2.3 muestra estos
procesos, donde un gas está contacto con un líquido y luego ocurre la mezcla.

Figura 2.3. Contacto de gas con un líquido y mezcla de fases (Suastegui, 1997)

De acuerdo a la situación descrita anteriormente, se establece que el tiempo

juega un papel importante en el coeficiente de transferencia de masa en la fase
líquida, ya que pueden afectar los remolinos formados. La ecuación (8) expresa esto
(UAH, sf):

r
2 D
KL = √ (8)
π θ
donde,

KL = Coeficiente de transferencia de masa fase líquida, m/s.

D= Coeficiente de difusividad de A en película, m2 /s.

θ= Tiempo de contacto, s.

Higbie aplicó esta ecuación con la que demostró que si el tiempo de contacto
entre las moléculas que se mueven por difusión es muy corto, no lograrán moverse lo
suficiente. El espesor correspondiente a la penetración es en el que la concentración
de la sustancia varíe en un 1 % (Suastegui, 1997).
8

2.5. Correlaciones de torre de pared húmeda

Para la transferencia de masa en el interior de tuberías se ha encontrado que

existe una analogía con la transferencia de calor debido a que la forma de las ecua-
ciones básicas para la difusión y para la conducción de calor son muy similares. Al
encontrarse en un régimen de flujo turbulento, el número de Sherwood se puede
relacionar analogamente al número de Nusselt, obteniendo un valor límite de 3.66
cuando se tiene una concentración de pared constante y un tercio dependiente de la
velocidad de flujo para tubos cortos. Cuando se tiene una variación de la concentra-
ción de la pared en la dirección axial, el número de Sherwood tiende a ser un poco
mayor que que para una concentración constante(McCabe, 2007).
Una correlación utilizada para la transferencia de materia en régimen turbulento
en las paredes interiores de una tubería, la cual es análoga a la transferencia de calor
reemplazando Nusselt y Prandtl por Scherwood y Schmidt, es la siguiente(McCabe,
2007):

 0,14
µ
(9)
0,8 1/3
Sh = 0,023Re Sc
µW
donde,

Sh= Número de Sherwood, adim.

Re= Número de Reynolds, adim.

Sc= Número de Schmidt, adim.

µ= Viscosidad dinámica, Pa · s.

Cuando se tiene una película de líquido formando una pared húmeda alrededor
de la pared interna de una tubería y se da la transferencia de masa mediante difusión
en un régimen de flujo laminar se puede utilizar la expresión(McCabe, 2007):

Sh = 0,023Re0,81 Sc0,44 (10)

donde,

Sh= Número de Sherwood, adim.

 9

Re= Número de Reynolds, adim.

Sc= Número de Schmidt, adim.

Además se ha encontrado una correlación para el caso en que la transferencia

de masa se dé en ámbitos elevados del número de Schmitd (430 a 100000). Esto se
determinó mediante la medición la velocidad de disolución de tubos de ácido benzoi-
co en agua y líquidos viscosos, la cual se representa de la siguiente manera(McCabe,
2007):

Sh = 0,0096Re0,913 Sc0,346 (11)

donde,

Sh= Número de Sherwood, adim.

Re= Número de Reynolds, adim.

Sc= Número de Schmidt, adim.

Por otro lado, también es importante recalcar que se puede dar la formación
de ondulaciones o pequeños vórtices debido a un gradiente de tensión interfacial.
A este fenómeno se le conoce como efecto Marangoni. Y la correción que se utiliza
para este caso es la siguiente (Welty, 2008):

 0,15
4Γ
Sh = 0,00814Re 0,83
Sc 0,44
(12)
µ
donde,

Sh= Número de Sherwood, adim.

Re= Número de Reynolds, adim.

Sc= Número de Schmidt, adim.

µ= Viscosidad, Pa · s.

Γ= Caudal por unidad de perímetro mojado, m2 .

10

2.6. Surfactantes

Los surfactantes o tensoactivos son sustancias químicas que modifican la inter-

acción entre dos superficies. Estos agentes modifican la superficie de los líquidos ye
en las interfases entre dos líquidos ó un líquido y un gas. Los surfactantes tienen la
característica de poseer un grupo hidrófilo unido a una larga cadena hidrocarbonada
que es liposoluble. Lo anterior permite a los surfactantes unir la superficie del agua
con otros líquidos en una emulsión o con sólidos en caso de una suspensión. Los
surfactantes son usados como detergentes y humectantes (Fuentes de Piedrahita,
2002).

2.7. Presencia de solutos en la transferencia de masa de las

torres de pared húmeda

La presencia de diferentes sustancias dentro de una corrida en una torre de pared

húmeda puede afectar de diferentes maneras la transferencia de masa que ocurre,
dependiendo de la sustancia que esté disuelta en el líquido, ya que es posible que
ocurran interacciones con el gas presente. Ching et al. (2014) estudiaron la absorción
de dióxido de carbono por parte de una disolución de amoniaco para la reducción de
las concentraciones de este gas en el aire. Variaron la concentración entre 1 % y un
15 % v/v y descubrieron que al aumentar la concentración la absorción de dióxido de
carbono aumentaba notablemente, esto debido a que la disolución formada tiene una
gran capacidad de absorber el gas. Además, debido a estos resultados descubrieron
que el coeficiente global de transferencia de masa del gas aumentó.
Servia et al. (2014) estudiaron el efecto que las disoluciones de piperazina en
la absorción de dióxido de carbono. Al igual que con la disolución de amoniaco, al
aumentar la concentración de piperazina el flujo de dióxido de carbono absorbido
aumentó, esto debido a que ocurre una reacción entre el soluto y el gas.
Jyant (2017) realizó un recuento de procesos de absorción de gases en torres
de pared húmeda. Las disoluciones de amoniaco también pueden ser utilizadas para
absorber dióxido de azufre en este equipo. Además, soluciones de glicinato de sodio
también pueden ser utilizadas para absorber dióxido de carbono. Existe una relación
proporcional entre la cantidad de soluto utilizado y la cantidad de gas absorbido,
lo que indica que los coeficientes de transferencia de masa aumentan debido a la
presencia del soluto en el líquido.

Cuadro 2.1. Datos para ecuación de Antoine con unidades de mmHg y Kelvin
para etanol y agua (Ríos, 2004)
Sustancia Parámetro A Parámetro B Parámetro C
Etanol 18.9119 3803.98 -41.68
Agua 18.3036 3816.44 -46.13

11
 NOMENCLATURA
Mayúsculas

A Constante de Antoine para ◦ C y kPa

Atrans área de transferencia de masa, m2

B Constante de Antoine para ◦ C y kPa

C Constante de Antoine para ◦ C y kPa

CA Concentración de A, mol/m3

D Coeficiente de difusividad, m2 /s.

G flujo volumétrico de aire, m3 /s

Habs Humedad absoluta, kgH2 O /kgaireseco

Hen humedad absoluta a la entrada, kgH2 O/kgaireseco

Hsal humedad absoluta a la salida, kgH2 O/kgaireseco

HR Humedad relativa, %

J Flux molar de transferencia de masa, mol/m2 s.

K Coeficiente de transferencia de masa, m/s.

KG Coeficiente global de transferencia de masa, m · kgH2 O /kgaireseco · s,

longitud, m

Masa molecular del agua, kg/kmol

Masa molecular del aire, kg/kmol

Presión total, Pa

Flujo de agua, m3 /s

Número de Reynolds, adim.

Número de Schmidt, adim.

13
14

Número de Sherwood, adim.

Volumen, m3

coeficiente de transferencia de masa, m/s

Tiempo, s

velocidad en el seno del fluido, m/s

Fracción molar de agua, adim.

Distancia en la dirección de la transferencia de masa, m

espesor de película, m

Caudal por unidad de perímetro mojado, m2 .

Viscosidad, Pa · s.

Integral de colisión, adim.

Tiempo de contacto, s

Se refiere a la sustancia A en la sustancia B

Se refiere a entrada

Se refiere a líquido

Se refiere a salida

Se refiere a saturación

Se refiere a la pared.
BIBLIOGRAFÍA

Cengel, Y. (2007). Transferencia de Calor y Masa. México: McGraw-Hill.

Ching, HEE and Co, LMP and Tan, SIC and Roces, SA and Dugos, NP and Robles,
J and Uy, MM and others (2014). Carbon dioxide absorption in a fabricated
wetted-wall column using varying concentrations of aqueous ammonia. ASEAN
Journal of Chemical Engineering, 2:9–21.

Fuentes de Piedrahita, C. Román, A. (2002). Los Surfactantes: clases, propiedades y

uso con herbicidas, volumen 1. Cali: Centro Internacional de Agricultura Tropical.

Gutiérrez, , Teixidó, C. M., y López, J. C. (1987). Diseño de columnas de paredes

mojadas. Ingeniería química, (217):161–171.

Himmelblau, D. (1997). Principios Básicos y Cálculos en Ingeniería Química. (6ta.

edición). México: Prentice Hall.

Izquierdo, J. y Torres, J. (2004). Cinética de las Reacciones Químicas. Metodologia

(Universitat de Barcelona). Publicacions i Edicions de la Universitat de Barcelona.

Jyant, S. (2017). Wetted wall column: Review on studies and investigations. Inter-
national Journal of Research and Review, 4(6):24–27.

McCabe, W., S. J. . H. P. (2007). Operaciones Unitarias en Ingeniería Química.

México: McGraw-Hill.

Ríos, L. G. (2004). Modelos termodinámicos para el equilibrio vapor–liquido a bajas

presiones: fase liquida, modelo de wilson. Scientia et technica, 1(24).

Servia, A., Laloue, N., Grandjean, J., Rode, S., y Roizard, C. (2014). Modeling of
the co2 absorption in a wetted wall column by piperazine solutions. Oil & Gas
Science and Technology–Revue d’IFP Energies nouvelles, 69(5):885–902.

Smith, J., Van Ness, H., y Abbott, M. (2007). Introducción a la termodinámica en

ingeniería química. (7ma edición). New York: McGraww Hill.

15
16

Suastegui, A. (1997). Análisis de los modelos en transferencia de masa para los

empaques estructurados. (Tesis de Título en ingeniertía química) Instituto Tecno-
lógico de Toluca, Toluca.

UAH (s.f.). Transferencia de materia gas-líquido. https://ebuah.uah.es/

dspace/bitstream/handle/10017/477/G.%20TRANSFERENCIA%20DE%
20MATERIA%20GAS-L%C3%8DQUIDO.pdf?sequence=6&isAllowed=y. Acce-
sado el 23 de agosto del 2019.

Welty, J. (2008). Fundamentals of Momentum, Heat, and Mass Transfer (5ta. edi-
ción), volumen 1. New Jersey: John Wiley and Sons, Inc.
APÉNDICES

A. DATOS EXPERIMENTALES

Cuadro A.1. Dimensiones de la columna, datos de presión atmosférica y coefi-

ciente de difusividad del agua en aire (Welty, 2008)
Dato Valor Dato Valor
Presión Total, 88110 Coeficiente de 2.634
Pt/(Pa) difusividad agua
en aire a 25°C,
Dab/(m2 P a/s)

Presión Total, 88.11 Coeficiente de 2.98945E-05

Pt/(kPa) difusividad agua
en aire a 25°C,
Dab/(m2 /s)

Diámetro de 0.03 Largo de 1.1

Columna, D/(m) Columna, L/(m)

Área transversal, 0.000706858 Área 0.207345115

At/({m2 }) transferencia,
At/({m2 })

Cuadro A.2. Flujo de aire y datos para obtener el flujo de agua que atraviesa la
torre
Corrida Flujo de aire Volumen de Agua, Tiempo flujo de
rotámetro, V/(mL) agua, t/(s)
Qa/(LPM)
1 61 160 9.93
2 58 160 9.93
3 50 160 9.93
4 40 160 9.93
5 30 160 9.93
6 20 160 9.93

17
18

Cuadro A.3. Temperaturas del agua y del aire a la entrada y a la salida de la

torre
Corrida Temperatura Temperatura Temperatura Temperatura
aire entrada, aire salida, agua superior, agua inferior,
Tai/(°C) Tao/(°C) Ths/(°C) Thi/(°C)
1 22.5 23.4 22 21.5
2 22.8 24.2 22 22
3 22.8 22.8 22 21
4 22.6 22.7 22.5 22
5 22.9 24.5 22 21
6 22.8 23.5 22.5 20

Cuadro A.4. Humedades relativas del aire a la entrada y a la salida de la torre

Corrida Humedad aire entrada, Humedad aire salida,
HRi/( %) HRo/( %)

1 2.2 74.5
2 2.7 72.4
3 2.9 75.4
4 2.9 75.1
5 1.9 70.1
6 3.7 72.2

B. DATOS INTERMEDIOS

Cuadro B.5. Datos de flujos y velocidades de agua y aire

Corrida Velocidad de Flujo de aire, Flujo de agua, Velocidad de
aire, V/(m/s) Qa/(m3 /s) Qh/(m3 /s) Agua,
V/(m/s)
1 1.438289115 0.001016667 1.61128E-05 0.022794934
2 1.367553585 0.000966667 1.61128E-05 0.022794934
3 1.178925504 0.000833333 1.61128E-05 0.022794934
4 0.943140404 0.000666667 1.61128E-05 0.022794934
5 0.707355303 0.000500000 1.61128E-05 0.022794934
6 0.471570202 0.000333333 1.61128E-05 0.022794934
 19

Cuadro B.6. Datos de presiones de agua obtenidas en las diferentes partes de la

torre
Corrida Logaritmo de Presión de Presión de Presión de
presión de saturación del agua entrada, agua salida,
saturación, agua, Pen /(kPa) Psal /(kPa)
lnPsat / Psat /(kPa)
(adim.)
1 0.962858939 2.619173839 0.057621824 1.95128451
2 0.978150549 2.659533015 0.071807391 1.896281859
3 0.947536949 2.579348762 0.074801114 1.974857074
4 0.993411869 2.700432294 0.078312537 1.966999553
5 0.947536949 2.579348762 0.049007626 1.836040861
6 0.932184488 2.540051835 0.093981918 1.891043512

Cuadro B.7. Datos de humedad absoluta a la entrada y salida de la torre

Corrida Humedad Humedad salida, Delta Humedad
entrada, Hsal /(kgH2 O /kgaire seco ) Absoluta,
Hen /(kgH2 0 /kgaire seco ) ∆H/(kgH2 0 /kgaire seco )
1 0.000406633 0.014072724 0.013666091
2 0.000506821 0.013667318 0.013160497
3 0.000527969 0.014246628 0.013718659
4 0.000552776 0.014188649 0.013635873
5 0.000345809 0.013223895 0.012878086
6 0.000663498 0.013628735 0.012965237

Cuadro B.8. Fracciones molares del agua en el aire

Corrida Fracción Fracción Fracción Diferencia
molar agua molar agua molar agua logarítmica
entrada salida, saturación, fracciones
yen /(adim) ysa /(adim.) ysat /(adim.) molares,
yalm /(adim.)
1 0.00065398 0.0221460 0.0297262 0.0159882
2 0.00081497 0.0215218 0.0301842 0.0169596
3 0.00084895 0.0224135 0.0292742 0.0151706
4 0.00088880 0.0223244 0.0306484 0.0168254
5 0.00055621 0.0208381 0.0292742 0.0165566
6 0.00106664 0.0214623 0.0288282 0.0153722
20

Cuadro B.9. Flux de agua y coeficiente global de transferencia de masa

Corrida Flux de agua, Coeficiente global de
NA /(m kgagua /s kgaire ) transferencia de masa,
KG /(m kgagua /s kgaire )
1 6.70084E-05 0.004191123
2 6.13557E-05 0.003617759
3 5.51362E-05 0.003634417
4 4.38428E-05 0.002605754
5 3.10547E-05 0.001875667
6 2.08432E-05 0.001355909

Cuadro B.10. Números adimensionales, factor de fricción de fanning y esfuerzo

cortante
Corrida Número de Número de Reynolds Factor de Esfuerzo
Schmidt, Sherwood, del agua, fricción, cortante,
Sc/(adim.) Sh/(adim.) Re/(adim.) f/(adim.) τ /(kg/ms2 )
1 0.5123033 4.20592128 683.8480267 0.007822041 0.019310358
2 0.5140450 3.63053228 683.8480267 0.007822041 0.017423691
3 0.5118294 3.64724949 683.8480267 0.007822041 0.012980793
4 0.5113561 2.61495396 683.8480267 0.007822041 0.008312117
5 0.5146800 1.88228994 683.8480267 0.007822041 0.004658204
6 0.5129358 1.36069659 683.8480267 0.007822041 0.002074355

Cuadro B.11. Caudal por unidad de perímetro mojado

Corrida Espesor, Diámetro Reynolds Caudal por 4Γ/µG
δ/(m) por donde del aire, unidad de /(kg 2 /s5 )
pasa el Re/(adim.) perímetro
aire, mojado,
Dh/(m) Γ/(kgm/s2 )
1 0.0006608 0.02934 2755.3528 9.3731E-07 0.20513769
2 0.0006608 0.02934 2610.9669 9.3772E-07 0.20493151
3 0.0006608 0.02934 2260.5768 9.3869E-07 0.20552116
4 0.0006608 0.02934 1810.1354 9.3971E-07 0.20582551
5 0.0006608 0.02934 1348.8338 9.4050E-07 0.20543269
6 0.0006608 0.02934 902.2803 9.4107E-07 0.20585297
 21

Cuadro B.12. Error en el número de Sherwood obtenido experimentalmente con-

tra el obtenido con la correlación de Polymath sin considerar efec-
tos Mangaroni
Corrida Número de Número de Porcentaje de
Sherwood Sherwood error, %E/(adim.)
calculado, experimental,
Shc/(adim.) She/(adim.)
1 4.1950454 4.205921284 0.259255455
2 3.8487054 3.630532282 5.668740398
3 3.3153928 3.64724949 10.00957383
4 2.664276 2.614953964 1.851235988
5 1.9968715 1.882289935 5.738053986
6 1.3341254 1.36069659 1.99165608

Cuadro B.13. Error en el número de Sherwood obtenido experimentalmente con-

tra el obtenido con la correlación de Polymath considerando efec-
tos Mangaroni
Corrida Número de Número de Porcentaje de
Sherwood Sherwood error, %E/(adim.)
calculado, experimental,
Shc/(adim.) She/(adim.)
1 4.1914792 4.205921284 0.34455817
2 3.8437979 3.630532282 5.548304655
3 3.3120043 3.64724949 10.12212424
4 2.6596101 2.614953964 1.679048224
5 1.9953003 1.882289935 5.66382738
6 1.3343153 1.36069659 1.977140609
22

C. MUESTRA DE CÁLCULO

C.1. Cálculo del flujo de agua

Para el cálculo del flujo de agua se divide el flujo entre el tiempo:

V
Q= (C.1)
t

donde,

Q= Flujo de agua, m3 /s

V = Volumen, m3

t= Tiempo, s

Utilizando los datos del Cuadro A.2 columna 2 y 3, fila 2 se obtiene el dato
tabulado en el Cuadro B.5 columna 4 fila 3:

160mL 1L 1m3
Q= ∗ ∗
s 1000mL 1000L

Q = 1, 611 ∗ 10−5 m3 /s

C.2. Cálculo de la presión de saturación

Para el cálculo de la presión de saturación se utilizará la ecuación de An-

toine (Smith et al., 2007):

B
lnP sat = A − (C.2)
T +C

donde,

P sat = Presión de saturación, kPa.

 23

A = Constante de Antoine para ◦ C y kPa.

B = Constante de Antoine para ◦ C y kPa.

C = Constante de Antoine para ◦ C y kPa.

Utilizando un promedio de los datos del Cuadro A.3 fila 2, columna 4 y 5, así
como las constantes de Antoine en el Cuadro A.14 se puede obtener el dato en el
Cuadro B.6 columna 3 fila 2:

1687, 537
lnP sat = 5,11564 −
21, 75°C + 230,17
Psat = 2, 619kP a

C.3. Cálculo de la presión parcial de agua

Para el cálculo de la presión parcial del agua se utilizará la siguiente ecua-

ción (Himmelblau, 1997):

PA
HR = · 100 (C.3)
PAsat

donde,

HR = Humedad relativa, %

PA = Presión parcial de agua, Pa

Como se observa en la fila 2 y columna 2 del Cuadro A.4 y el Cuadro B.6

columna 3 fila 2 se puede obtener el dato en el Cuadro B.6 columna 4 fila 2:

HR ∗ PAsat
PA =
100

2, 2 % ∗ 2, 619kP a
PA =
100
PA = 0, 0576kP a
24

C.4. Cálculo de la humedad absoluta del aire

Para el cálculo de la humedad absoluta del aire se utiliza la siguiente ecua-

ción (Himmelblau, 1997):

PA MA
Habs = · (C.4)
PT − PA MB

donde,

Habs = Humedad absoluta, kgH2 O /kgaireseco

PA = Presión parcial de agua, Pa

PT = Presión total, Pa

MA = Masa molecular del agua, kg/kmol

MB = Masa molecular del aire, kg/kmol

Utilizando los datos del Cuadro B.6 columna 3 fila 4 y del Cuadro A.1 columna
2 fila 3, se puede obtener el dato en el Cuadro B.7 columna 2 fila 2:

0, 0576kP a 18, 02kg/kmol

Habs = ·
88, 11kP a − 0, 0576kP a 29kg/kmol
Habs = 0, 000406kgH2 O /kgaireseco

C.5. Cálculo de la diferencia en la humedad absoluta

Para el cálculo del delta humedad absoluta del aire se utiliza la siguiente
ecuación:

∆H = Hsal − Hen (C.5)

donde,
 25

Hsal = Humedad absoluta de salida, kgH2 O /kgaireseco

Hen = Humedad absoluta de entrada, kgH2 O /kgaireseco

∆Hen = Diferencia de Humedad absoluta, kgH2 O /kgaireseco

Utilizando los datos del Cuadro B.7 columna 3 y 4, fila 2 se puede obtener el
dato en el Cuadro B.7 columna 4 fila 2:

∆H = (0, 0140 − 0, 000406)kgH2 O /kgaireseco

∆H = 0, 001366kgH2 O /kgaireseco

C.6. Cálculo de la fracción molar de agua

Para el cálculo de la fracción molar de agua se utiliza la siguiente ecuación:

PA
yA = (C.6)
P

donde,

yA = Fracción molar de agua, adim.

PA = Presión de vapor de agua, kPa.

P = Presión total, kPa.

Utilizando los datos del Cuadro B.6 columna 4, fila 2 y el dato en el Cuadro
A.1 columna 2 fila 3, se puede obtener el dato en el Cuadro B.8 columna 2 fila 2:

0, 0576kP a
yA =
88, 11kP a
yA = 0, 000653

C.7. Cálculo del flux

Para el cálculo del flúx de agua se usa esta ecuación (McCabe, 2007):
26

G
NA = ∗ (Hsal − Hen ) (C.7)
Atrans

donde,

G = Flujo volumétrico de aire, m3 /s

Atrans = Área de transferencia de masa, m2

Hen = Humedad absoluta a la entrada, kgH2 O/kgaireseco

Hsal = Humedad absoluta a la salida, kgH2 O/kgaireseco

NA = Flux de agua, m · kgH2 O /kgaireseco · s

Utilizando los datos de el Cuadro B.7 en las columnas 2 y 3 en la fila 2; los

datos de la columna 3 y fila 2 de el Cuadro B.5 se tiene que:

0,0101m3 /s
NA = ∗ (0,01407kgH2 O /kgaireseco − 0,0004066kgH2 O /kgaireseco )
0,2073m2

m · kgH2 O
NA = 0,000067
kgaireseco · s

C.8. Cálculo de la diferencia logarítmica de fracciones molares

Para el cálculo de la diferencia logarítmica de fracciones molares se utiliza

la siguiente ecuación:

yAsat − yAo − (yAsat − yAi )

(yAsat − yA )lm = y sat −y
(C.8)
ln( yAsat −yAo
Ai
)
A

donde,

yA = Fracción molar de agua, adim.

yAo = Fracción molar de agua de salida, adim.

yAi = Fracción molar de agua de entrada, adim.

 27

Utilizando los datos del Cuadro B.8 columna 2, 3 y 4 fila 2, se puede obtener
el dato en el mismo Cuadro B.8 columna 5 fila 2:

0, 0297 − 0, 0221 − (0, 0297 − 0, 00065)

(yAsat − yA )lm = 0,0297−0,0221
ln( 0,0297−0,00065 )

(yAsat − yA )lm = 0, 0159

C.9. Cálculo del coeficiente global de transferencia de masa

Para el cálculo del coeficiente global de transferencia de masa se utiliza la

siguiente ecuación (Welty, 2008):

NA = KG · (yAsat − yA )lm (C.9)

donde,

KG = Coeficiente global de transferencia de masa, m · kgH2 O /kgaireseco · s

Utilizando los datos del Cuadro B.8 columna 5 fila 2 y del Cuadro B.9 fila 2
columna 2, se puede obtener el dato en el mismo Cuadro B.9 columna 3 fila 2:

6, 7 ∗ 10−5
KG =
0, 0159
KG = 0, 00419m · kgH2 O /kgaireseco · s

C.10. Cálculo del número de Reynolds

Para el cálculo del número de Reynolds se utiliza la siguiente ecuación

(Welty, 2008):

ρ·D·v
Re = (C.10)
µ

donde,
28

ρ = Densidad, kg/m3

D = Diámetro de tubería

v = Velocidad de flujo, m/s

µ = Viscosidad, Pa s

Utilizando los datos del Cuadro A.1 columna 2 fila 3, del Cuadro B.5 columna
2 fila 2 y los datos de densidades y viscosidades en el Anexo, se puede obtener el
dato en el Cuadro B.11 columna 3 fila 2:

1, 19kg/m3 · 0, 03m · 1, 43m/s

Re =
1, 82 ∗ 10−5
Re = 2755

C.11. Cálculo del número de Schmidt del aire

Para el cálculo del número de Schmidt se usa esta ecuación (McCabe, 2007):

µ
Sc = (C.11)
ρ ∗ DAB

donde,

µ = Viscosidad dinámica, Pa ∗ s

ρ = Densidad, kg/m3

DAB = Difusividad, m2 /s

Sc = Número de Schmidt, adim.

Según los datos de la columna 4 y fila 3 de el Cuadro A.1 y los datos de las
columnas 2 y 3 en la fila 2 de el Cuadro A.15 se tiene que el número de Schmidt es:

0,000018277Pa ∗ s
Sc =
1,1933kg/m3 ∗ 0,00002989m2 /s

Sc = 0,512
 29

C.12. Cálculo del número de Sherwood del aire

Para el cálculo del número de Sherwood se usa esta ecuación (McCabe, 2007):

k∗D
Sh = (C.12)
DAB
donde,

µ = Viscosidad dinámica, Pa ∗ s

ρ = Densidad, kg/m3

DAB = Difusividad, m2 /s

Sc = Número de Sherwood, adim.

Según los datos de la columna 4, fila 3; columna 2 y fila 4 de el Cuadro A.1 y

la columna 3, fila 2 de el Cuadro B.9 se tiene que:

0,03m
Sh =
0,00002989m2 /s
Sh = 4,205

C.13. Cálculo del factor de fricción

Para el cálculo del factor de fricción respecto al Reynolds del agua se usa esta

ecuación (McCabe, 2007):

f = 0,04 ∗ Re−1/4 (C.13)

donde,

Re= Número de Reynolds, adim.

f = Factor de fricción, adim.

Según los datos de la columna 4, fila 1 de el Cuadro B.10 se tiene que:

f = 0,04 ∗ (683,848)−1/4

f = 0,00782
30

C.14. Cálculo del esfuerzo cortante

Para el cálculo del esfuerzo cortante se usa esta ecuación:

τg = f ρg vg2 (C.14)

donde,

f = Factor de fricción, adim.

ρg = Densidad del gas, kg/m3

vg = Velocidad del gas, m/s

Con el dato de la columna 5, fila 2 de el Cuadro B.10, la columna 2 y fila 2 de

el Cuadro A.15 y la columna 2 y fila 2 de el Cuadro B.5 se tiene que:

τg = 0,00782 ∗ 1,1933sikg/m3 ∗ (1,4382m/s)2

τ = 0,0193kg/(m ∗ s2 )

C.15. Cálculo del espesor de película

Para el cálculo del espesor de película se usa esta ecuación (Gutiérrez et al.,

1987):
3 ∗ νl2 1/3 8/15
δ=( ) ∗ Rel ∗ 400−1/5 (C.15)
g
donde,

g = Aceleración de la gravedad, m/s2

Rel = Reynolds del agua, adim

νl = Viscosidad cinemática del agua, m2 /s

Con el dato de la columna 4 y fila 2 de el Cuadro B.10, mas los datos de las
columnas 1 y 2 en la fila 2 del Cuadro A.16:

3 ∗ (998 ∗ 10−6 P a ∗ s/998kg/m3 )2 1/3

δ=( ) ∗ 683,8488/15 ∗ 400−1/5
9,81m/s2
 31

δ = 0,0006608m

C.16. Cálculo del caudal por unidad de perímetro mojado

Para el cálculo del del caudal por unidad de perímetro mojado se utiliza
la siguiente ecuación (Gutiérrez et al., 1987):

ρgδ 3 τG δ 2
Γ= − (C.16)
3 2

donde,

Γ = Caudal por unidad de perímetro mojado, kgm/s2

δ = Espesor de película de agua, m

τG = Esfuerzo cortante del aire, kg/ms2

Utilizando los datos del Cuadro B.10 columna 6 fila 2, del Cuadro B.11 fila 2
columna 2 y del Cuadro de propiedades en el Anexo, se puede obtener el dato en el
mismo Cuadro B.11 columna 5 fila 2:

1000kg/m3 · 9, 81m/s2 · (0, 00066m)3 0, 00193kg/ms2 (0, 00066m)2

Γ= −
3 2

Γ = 9, 373 ∗ 10−7 kgm/s2

4∗Γ
= 0, 205kg 2 /s5
µG

4 ∗ 9, 373 ∗ 10−7
−5
= 0, 205kg 2 /s5
1, 8277 ∗ 10

C.17. Cálculo del número de Sherwood con la correlación obtenida por

Polymath
32

Se obtuvo la siguiente correlación con el programa Polymath para calcular

el número de Sherwood sin considerar efectos de Mangaroni:

Sh = 0, 006Re0,9924 Sc2,0508 (C.17)

donde,

Sh = Número de Sherwood, adim.

Re = Número de Reynolds, adim.

Sc = Número de Schmidt, adim.

Utilizando los datos del Cuadro B.10 columna 2 fila 2 y del Cuadro B.11 fila 2
columna 4 se obtiene el dato del Cuadro B.12 columna 2 fila 2:

Sh = 0, 006 ∗ (2755, 35)0,9924 (0, 5123)2,0508

Sh = 4, 195

El procedimiento es el mismo para cuando se consideran efectos de Mangaroni

solo que se utiliza como término adicional el caudal por unidad de perímetro mojado.

C.18. Cálculo del error entre el número de Sherwood obtenido experi-

mentalmente y el obtenido con la correlación de Polymath

Para el cálculo del error se considera la siguiente ecuación:

T eorico − Experimental
Error = · 100 (C.18)
T eorico

Utilizando los datos del Cuadro B.12 columna 2 y 3 fila 2 se obtiene el dato del
Cuadro B.12 columna 4 fila 2:
 33

4, 206 − 4, 195
Error = · 100
4, 206
Error = 0, 295 %
ANEXOS

A. Datos obtenidos de referencias bibliográficas

Cuadro A.14. Constantes de Antoine para el agua (Smith et al., 2007)

Sustancia Constante A Constante B Constante C
Agua 5.11564 1687.537 230.17

Cuadro A.15. Propiedades reológicas del aire para las diferentes corridas (Welty,
2008)
Corrida Densidad del aire, Viscosidad del aire,
ρ/(kg/m )3
µ/(Pa s)
1 1.19338 1.8277E-05
2 1.19105 1.8303E-05
3 1.19401 1.8269E-05
4 1.19465 1.8262E-05
5 1.19021 1.8313E-05
6 1.19253 1.8286E-05

Cuadro A.16. Propiedades reológicas del agua (Welty, 2008)

Densidad del agua, ρ/(kg/m3 ) Viscosidad del agua, µ/(Pa s)
998 998*10−6

35