24/10/2022

Curso de Transferencia de Masa

Prof. Francisco Javier Garfias Vásquez
Facultad de Química, UNAM

Concentración Perfil de
de la capa concentración
límite

Especie A

1
 24/10/2022

2.1 INTRODUCCION

Transporte de masa que implica el transporte de materia entre,

Transferencia de masa
 una superficie límite y un fluido móvil
por convección  dos fluidos relativamente inmiscibles en movimiento

Difusión molecular Transferencia de masa a través de,

es responsable de  un disolvente en reposo o
 en un régimen laminar fluyendo en capas

2.1 INTRODUCCION
Flujo turbulento

Espesor muy pequeño de la película y difícil de predecir

• Los remolinos contribuyen considerablemente al proceso de transferencia de masa debido a

su rápido movimiento.
• La transferencia de masa en el régimen turbulento es mucho más rápido que en condiciones
de flujo en régimen laminar.
• Esta situación no puede ser modelada en términos de la ley de Fick.
• Es por esto que la transferencia de masa es determinada a partir de los coeficientes de
transferencia de masa, una aproximación que simplifica el análisis de un problema muy
complejo a resolver.

2
 24/10/2022

2.2 Coeficiente de transferencia de masa

Objetivo:
•Explicar el concepto de coeficiente de transferencia de masa para difusión en flujo turbulento por analogía
con la difusión molecular.
•Definir y usar los coeficientes de transferencia de masa para la difusión del componente A, a través del
componente B en reposo, y para contra difusión equimolar.

Reemplazando, 𝐹 =
𝑐 DAB en la ecuación (1 – 74)
𝑧
𝑐
ΨA − A2ൗ𝑐
𝑁A = Ψ 𝐹 ln 𝑐 (1)
ΨA − A1ൗ𝑐

Coeficiente de transferencia de masa

Donde, 𝒄𝐀ൗ𝒄 corresponde a la fracción mol xA para líquidos, y la fracción mol yA para gases.

Flux= (coeficiente) (diferencia de concentración)

2.2.1 Difusión de A a través de B estancado ( NB = 0, A = 1 )

La ecuación (1) se escribe

1 − 𝑦𝐴2
𝑁𝐴 = 𝐹𝐺 ln 𝑃𝑎𝑟𝑎 𝑔𝑎𝑠𝑒𝑠 (2)
1 − 𝑦𝐴1
1 − 𝑥𝐴2
𝑁𝐴 = 𝐹𝐿 ln 𝑃𝑎𝑟𝑎 𝑙í𝑞𝑢𝑖𝑑𝑜𝑠 (3)
1 − 𝑥𝐴1

Manipulando las ecuaciones (2) y (3),

𝐹𝐺
𝑁𝐴 = 𝑦 − 𝑦𝐴2 = 𝑘𝑦 𝑦𝐴1 − 𝑦𝐴2 (4)
𝑦𝐵,𝑀 𝐴1
𝐹𝐿
𝑁𝐴 = 𝑥 − 𝑥𝐴2 = 𝑘𝑥 𝑥𝐴1 − 𝑥𝐴2 (5)
𝑥𝐵,𝑀 𝐴1
Donde,
𝑦𝐵2 − 𝑦𝐵1 𝑥𝐵2 − 𝑥𝐵1 𝐹𝐺 𝐹𝐿 𝐹𝐺 = 𝑘𝑦 𝑦𝐵,𝑀
𝑦𝐵,𝑀 = 𝑦𝐵2 𝑥𝐵,𝑀 = 𝑥𝐵2 𝑘𝑦 = 𝑘𝑥 =
ln ln 𝑦𝐵,𝑀 𝑥𝐵,𝑀 𝐹𝐿 = 𝑘𝑥 𝑥𝐵,𝑀
𝑦𝐵1 𝑥𝐵1

3
 24/10/2022

Otras definiciones de coeficientes de transferencia de masa son,

𝑁𝐴 = 𝑘𝐺 𝑝𝐴1 − 𝑝𝐴2 = 𝑘𝑐 𝑐𝐴1 − 𝑐𝐴2 (7)

𝑁𝐴 = 𝑘𝐿 𝑐𝐴1 − 𝑐𝐴2 (8)

Donde es fácil mostrar que,

𝑝𝐵,𝑀 𝑝𝐵,𝑀
𝐹𝐺 = 𝒌𝑮 𝑝𝐵,𝑀 = 𝑘𝑐 = 𝒌𝒚 𝑑𝑜𝑛𝑑𝑒 𝒌𝒚 = 𝒌𝑮 𝑷 (9)
𝑅𝑇 𝑷

𝐹𝐿 = 𝑘𝐿 𝑥𝐵,𝑀 𝑐 = 𝑘𝑥 𝑥𝐵,𝑀 (10)

𝑅𝑇 𝑃
De lo anterior, 𝑘𝑐 = 𝑘𝑦 𝑜 𝑘𝑦 = 𝑘𝑐 y 𝑘𝑥 = 𝑘𝐿 𝑐
𝑃 𝑅𝑇

Para soluciones diluidas, 𝑦𝐵,𝑀 ≡ 𝑥𝐵,𝑀 ≅ 1.0 𝑦 𝑝𝐵,𝑀 ≅ 𝑃

Nota importante:
Para soluciones no diluidas usar las ecuaciones 2 y 3
1 − 𝑦𝐴2
𝑁𝐴 = 𝐹𝐺 ln 𝑃𝑎𝑟𝑎 𝑔𝑎𝑠𝑒𝑠 (2)
1 − 𝑦𝐴1
1 − 𝑥𝐴2
𝑁𝐴 = 𝐹𝐿 ln 𝑃𝑎𝑟𝑎 𝑙í𝑞𝑢𝑖𝑑𝑜𝑠 (3)
1 − 𝑥𝐴1

Para soluciones diluidas, usar las ecuaciones 4, 5, 7 y 8

𝐹𝐺
𝑁𝐴 = 𝑦 − 𝑦𝐴2 = 𝑘𝑦 𝑦𝐴1 − 𝑦𝐴2 (4)
𝑦𝐵,𝑀 𝐴1
𝐹𝐿
𝑁𝐴 = 𝑥 − 𝑥𝐴2 = 𝑘𝑥 𝑥𝐴1 − 𝑥𝐴2 (5)
𝑥𝐵,𝑀 𝐴1

𝑁𝐴 = 𝑘𝐺 𝑝𝐴1 − 𝑝𝐴2 = 𝑘𝑐 𝑐𝐴1 − 𝑐𝐴2 (7)

𝑁𝐴 = 𝑘𝐿 𝑐𝐴1 − 𝑐𝐴2 (8)

4
 24/10/2022

Problema 1:
Calcular el coeficiente de transferencia de masa kx, basado en la fracción mol, para la absorción de oxígeno
en agua a 310 K donde kL = 3.3x10-5 m/s.

Solución:
T = 310 K
m
kL = 3.3 x 10−5
s
kg
ρB (310 K) = 993 3
m
siendo 𝑘𝑥 = 𝑘𝐿 𝑐 si 𝑥𝐵,𝑀 ≅ 1.0

𝜌 𝜌𝐵 993 kmol
La concentración, 𝐶= = = = 55.2 3
𝑀𝑝𝑟𝑜𝑚 𝑀𝐵 18 m

m kmol 𝐤𝐦𝐨𝐥
Así, 𝑘𝑥 = 3.3 x 10−5 55.2 = 𝟏. 𝟖𝟐𝐱 𝟏𝟎−𝟑 𝟐
s m3 𝒎 𝒔

Problema 2:
Un absorbedor de gas es utilizado para remover amoniaco del aire a contracorriente con agua a 300 K y 1 atm.
En cierto punto del absorbedor, la fracción mol del amoniaco en la fase gas es 0.8, mientras que la concentración
en la interfase gas es 0.732. El flux de amoniaco medido en ese punto es de 4.3x10-4 kmol/m2s. Despreciando la
evaporación del agua, calcular el coeficiente de transferencia de masa en la fase gas.
Solución:

Difusión de amoniaco (A) a través de aire (B) no difundente, así, 𝑦𝐴1 = 0.8 y 𝑦𝐴2 = 0.732

Despejando la ecuación (2),

𝑁𝐴
1 − 𝑦𝐴2 𝐹𝐺 =
𝑁𝐴 = 𝐹𝐺 ln (2) 1 − 𝑦𝐴2
1 − 𝑦𝐴1 ln
1 − 𝑦𝐴1

kmol
𝑁𝐴 4.3 x 10−4
𝐹𝐺 = = m2 s = 𝟏. 𝟒𝟕 𝐱 𝟏𝟎−𝟑 𝐤𝐦𝐨𝐥
1 − 𝑦𝐴2 1 − 0.732 𝐦𝟐 𝐬
ln ln
1 − 𝑦𝐴1 1 − 0.8

5
 24/10/2022

2.2.2 Contra difusión equimolar ( NB = -NA , A = indefinido )

• Partiendo de la ecuación (1 – 82),

𝑁𝐴 = 𝑘𝐺´ 𝑃𝐴1 − 𝑃𝐴2 = 𝑘𝐶´ 𝐶𝐴1 − 𝐶𝐴2 = 𝑘𝑦´ 𝑦𝐴1 − 𝑦𝐴2 = 𝐹𝐺 𝑦𝐴1 − 𝑦𝐴2 (12)

𝑁𝐴 = 𝑘𝐿´ 𝐶𝐴1 − 𝐶𝐴2 = 𝑘𝑥´ 𝑥𝐴1 − 𝑥𝐴2 = 𝐹𝐿 𝑥𝐴1 − 𝑥𝐴2 (13)

• Para soluciones diluidas, donde

𝑃
𝐹𝐺 = 𝑘𝐺´ 𝑃 = 𝑘𝐶´ 𝑅𝑇 = 𝑘𝑦´ (14)

𝐹𝐿 = 𝑘𝐿´ 𝑐 = 𝑘𝑥´ (15)

Problema 3:
Una columna de destilación empacada, es utilizada para separar adiabáticamente una mezcla metanol-agua a
1 atm. El metanol (que es el más volátil), se difunde de la fase líquida hacia la fase vapor, mientras el agua se
difunde en sentido contrario.
Asumiendo que el calor latente de vaporización es casi igual para ambos componentes, el proceso es una
contra-difusión equimolar. En un punto de la columna, el coeficiente de transferencia de masa es 1.62x10-5
kmol/m2 s KPa. La fracción mol de metanol en la fase gas a la interfase es 0.707 y en el seno del gas es 0.656.
Estimar el flux de metanol en ese punto.

Solución:

Difusión de metanol (A) en (B) agua, así, 𝑦𝐴1 = 0.707 y 𝑦𝐴2 = 0.656

kmol kmol
𝑘𝑦´ = 𝑘𝐺´ 𝑃 = 1.62x10−5 101.3 kPa = 1.64x10−3
m2 s kPa m2 s

𝐤𝐦𝐨𝐥
𝑁𝐴 = 𝑘𝑦´ 𝑦𝐴1 − 𝑦𝐴2 = 1.64x10−3 0.707 − 0.656 = 𝟖. 𝟑𝟔𝐱𝟏𝟎−𝟓
𝐦𝟐 𝐬

6
 24/10/2022

2.4 Analogías de transferencia de masa y calor

Existen mas datos experimentales para transferencia de calor que para
transferencia de masa.
Las analogías permiten resolver problemas de transferencia de masa a partir
de información de transferencia de calor.
Para circunstancias análogas,
• los perfiles de temperatura y concentración en forma adimensional y
• los coeficientes de transferencias de masa y calor en forma adimensional,
están dados por los mismos grupos.

2.4 Analogías de transferencia de masa y calor

Para convertir datos de ecuaciones o correlaciones de

transferencia de calor en correlaciones de transferencia de

masa, los grupos adimensionales de transferencia de calor

son reemplazados por los grupos de transferencia de masa

correspondientes.

7
 24/10/2022

2.4 Analogía de transferencia de masa y calor

La siguiente tabla muestra los grupos adimensionales mas comunes:

Transferencia de masa Transferencia de calor

Numero de Reynolds Numero de Reynolds
l𝑣𝜌 l𝑣𝜌
l = longitud característica Re = Re =
𝜇 𝜇
Numero de Schmidtl 𝜇 Numero de Prandtl 𝐶𝑝 𝜇
Sc = Pr =
𝜌 DAB 𝑘
Numero de Sherwood Numero de Nusselt
𝐹 l 𝑘𝐺 𝑝𝐵,𝑀 𝑅𝑇 l 𝑘𝑐 𝑝𝐵,𝑀 l
Sh = , , , ℎl
𝑐 DAB 𝑃 DAB 𝑃 DAB Nu =
𝑘
𝑘𝑐′ l 𝑘𝑦′ 𝑅𝑇 l
, , etc.
DAB 𝑃 DAB

2.4 Analogía de transferencia de masa y calor

La siguiente tabla muestra los grupos adimensionales mas comunes:

Transferencia de masa Transferencia de calor

Numero de Grashof Numero de Grashof (β= coeficiente de volumen de expansión)
g l 3 𝜌 ∆𝜌 g l 3 𝜌2 𝛽 ∆𝑇
Gr𝐷 = Gr𝐻 =
𝜇2 𝜇2
Numero de Peclet Numero de Peclet
l𝑣 𝐶𝑝 l 𝑣 𝜌
Pe𝐷 = Re Sc = Pe𝐻 = Re Pr =
DAB 𝑘
Numero de Stanton Numero de Stanton
Sh Sh 𝐹 Nu Nu ℎ
St 𝐷 = = = St 𝐻 = = =
Re Sc Pe𝐷 𝑐 𝑣 Re Pr Pe𝐻 𝐶𝑝 𝑣 𝜌

Factor j de Chilton- Colburn Factor j de Chilton- Colburn

𝑗𝐷 = St 𝐷 Sc 2Τ3 𝑗𝐻 = St 𝐻 Pr 2Τ3

8
 24/10/2022

 Para las analogías de transferencia de masa y calor, las limitaciones son:

a) Las condiciones de flujo y geometría deben ser las mismas
b) La mayoría de los datos de transferencia de calor están basadas en situaciones que no
consideran transferencia de masa.
El uso de las analogías produce coeficientes de transferencia de masa más cercanos a
k’G , k’C o k’y (=F).
El número de Sherwood, Sh, se escribe comúnmente en términos de cualquiera de los
coeficientes, pero cuando es derivado del reemplazo del número de Nusselt, Nu, para
el caso en que N ≠ 0, el número de Sh debe ser tomado como Sh = F l /c DAB y F
utilizado con la ecuación (1), 𝐶A2
ΨA − ൗ
𝑁A = Ψ 𝐹 ln 𝐶 (1)
𝐶A1
ΨA − ൗ
𝐶

Problema:
Para un flujo que se desplaza en un ángulo recto a un cilindro circular, el coeficiente de
transferencia de calor promedio alrededor del cilindro esta dado por la ecuación de Eckert
y Drake (1959),
𝑁𝑢𝑝𝑟𝑜𝑚 = 0.43 + 0.532 𝑅𝑒 0.5 𝑃𝑟 0.31 𝑝𝑎𝑟𝑎 1 ≤ 𝑅𝑒 ≤ 4000

Donde el número de Nusselt promedio, Nuprom, y el número de Reynolds, Re, son

utilizados considerando l = d (diámetro del cilindro) , y las propiedades del fluido son
evaluadas a condiciones promedio de una película de fluido estancado alrededor de un
solido (temperatura y concentración promedio).

Estimar la rapidez de sublimación de un cilindro de UF6 (M=352), de 1 cm diámetro y

10 cm de longitud expuesto a una corriente de aire que fluye normal al eje del cilindro a
una velocidad de 1 m/s. La temperatura de la superficie del solido es de 303 K, a la cual la
presión de vapor del UF6 es 27 kPa. En el seno del aire la presión y la temperatura son 1
atm y 325 K, respectivamente.
En la siguiente diapositiva se muestra un esquema del problema.

9
 24/10/2022

Esquema del
problema Película de
gas estancado
Corriente de aire

v = 10 m/s
P = 1 atm

Sólido UF6 1 NA 2
303 K 325 K
pA= 27 kPa pA= 0 kPa

Solución: Película de
gas
estancado Corriente de aire
Reemplazando 𝑵𝒖𝒑𝒓𝒐𝒎 por 𝑺𝒉𝒑𝒓𝒐𝒎 𝑦 𝑷𝒓 por 𝑺𝒄 en la correlación, la analogía para v = 10 m/s
transferencia de masa es, 1 2 P = 1 atm
Sólido NA

𝑆ℎ𝑝𝑟𝑜𝑚 = 0.43 + 0.532 𝑅𝑒 0.5 𝑆𝑐 0.31 1 ≤ 𝑅𝑒 ≤ 4 000 UF6 303 K

325 K
pA= 0 kPa
pA= 27 kPa

303 + 325
La temperatura promedio es: 𝑇𝑝𝑟𝑜𝑚 = = 314 K
2

PUNTO 1: Gas saturado en UF6

PUNTO 2: Gas libre de UF6

𝑃 𝑀𝑝𝑟𝑜𝑚 101.3 x 72 kg
27 + 0 𝜌𝑝𝑟𝑜𝑚 = = = 2.8 3
𝑝𝐴 𝑝𝑟𝑜𝑚 = kPa = 13.5 kPa 𝑅 𝑇𝑝𝑟𝑜𝑚 8.314 x 314 m
2
13.5 kPa
𝑦𝐴 𝑝𝑟𝑜𝑚 = = 0.133 𝜇𝑝𝑟𝑜𝑚 = ? 𝐌é𝐭𝐨𝐝𝐨 𝐝𝐞 𝐋𝐮𝐜𝐚𝐬 (𝟏𝟗𝟖𝟎)
101.3 kPa

𝑃𝑉 𝑃𝑉 𝑚 𝑃𝑀 𝑚 𝜇𝑟𝑀 = 𝜉𝑀 𝜇 = 𝑓 𝑇𝑟𝑀 (21)

𝑃𝑉 = 𝑛𝑅𝑇 =𝑛 =𝑛= = = 𝜌
𝑅𝑇 𝑅𝑇 𝑀 𝑅𝑇 𝑉
1ൗ
𝑇𝑐𝑀 6
kg
𝑀𝑝𝑟𝑜𝑚 = 352 0.133 + 29 1 − 0.133 = 72 𝜉𝑀 = 0.176 (22)
kmol 𝑀𝑝𝑟𝑜𝑚 3 𝑃𝑐𝑀 3

10
 24/10/2022

Ecuación

Resultados
experimentales

11
 24/10/2022

Resultados
experimentales

Problema (usando la analogia de Chilton- Colburn):

• En 1934, Chilton y Colburn, basados en datos experimentales, definieron el factor ‘’j’’ para transferencia de
masa que establece la analogía con las transferencia de calor.
𝑗𝐷 = 𝑆𝑡𝐷 𝑆𝑐 2/3 = Ψ 𝑅𝑒 = 𝑗𝐻 = 𝑆𝑡𝐻 𝑃𝑟 2/3 (27)
• Esta analogía es muy útil para relacionar datos de transferencia de calor y masa cuando el efecto del número
de Prandtl, Pr, y del número de Scmidt, Sc, no están dados explícitamente.

• Se desea estimar la rapidez a la cual el benceno se evapora de una superficie mojada cuando se hace pasar
nitrógeno gas a 1 bar y 300 K a una velocidad de 10 m/s . No existen datos de transferencia de masa para esta
situación, pero las mediciones de transferencia de calor indican que, para CO 2 a 300 K y 1 bar, el coeficiente
de transferencia de calor, h, entre el gas y la superficie esta dado por la expresión empírica.

ℎ = 20 𝐺𝑦 0.5 𝑊ൗ 𝐺𝑦 = 𝜌𝜐
𝑚2 𝐾

Donde Gy es la velocidad másica del gas superficial en kg/m2s . Estimar el coeficiente de transferencia de masa.

12
 24/10/2022

Convección

Objetivo:
• Estimar los coeficientes de transferencia de masa por convección en las
siguientes situaciones:
a) Flujo paralelo a una superficie plana
b) Flujo a través de una sola esfera
c) Flujo normal a un solo cilindro
d) Flujo turbulento en tubos circulares
e) Flujo a través de lechos empacados y lechos fluidizados
f) Flujo a través de módulos de membrana de fibra hueca

13
 24/10/2022

Consideraciones…
• Usar los coeficientes correspondientes para resolver problemas típicos de
transferencia de masa.

• Se han considerado las correlaciones de transferencia de masa

desarrollados a partir de analogías de transferencia de calor.

• En esta sección veremos varias correlaciones desarrolladas en la literatura a

partir de datos experimentales de transferencia de masa.

Consideraciones…
• Los datos experimentales fueron obtenidos inyectando gases sobre varias
superficies húmedas de diferente geometría, evaporando líquidos o
haciendo pasar líquidos a través de solidos los cuales se disuelven.

• En la mayoría de los casos, los datos son presentados en términos de

coeficientes tipo k aplicable a sistemas binarios considerando NB=0, siendo
que la mayoría de las concentraciones experimentales son usualmente
diluidas, permitiendo que F sea calculado tomando pB,M/P, xB,M, etc., igual a
la unidad.

14
 24/10/2022

4.3.1 Coeficientes de transferencia de masa para placas planas

Utilizados en la medición de la evaporación de una superficie líquida o la sublimación de un sólido, en placas
de largo L bajo una corriente de aire controlada en condiciones laminares o turbulentas.

Coeficientes de transferencia de masa promedio a lo largo de la superficie:

𝑆ℎ𝐿 = 0.664𝑅𝑒𝐿0.5 𝑆𝑐1/3 𝐿𝑎𝑚𝑖𝑛𝑎𝑟 𝑅𝑒𝐿 < 3x105

𝑆ℎ𝐿 = 0.036𝑅𝑒𝐿0.8 𝑆𝑐1/3 𝑇𝑢𝑟𝑏𝑢𝑙𝑒𝑛𝑡𝑜 𝑅𝑒𝐿 ≥ 3x105

Las ecuaciones son válidas si el número de Schmidt se encuentra en un rango de:

0.6 < 𝑆𝑐 < 2500

También se encuentran expresiones similares bajo las mismas condiciones en términos del factor 𝒋𝑫 :

𝑗𝐷 = 0.664𝑅𝑒𝐿−0.5 𝐿𝑎𝑚𝑖𝑛𝑎𝑟 𝑅𝑒𝐿 < 3x105

𝑗𝐷 = 0.036𝑅𝑒𝐿−0.2 𝑇𝑢𝑟𝑏𝑢𝑙𝑒𝑛𝑡𝑜 𝑅𝑒𝐿 ≥ 3x105

Evaporación de Benceno a lo largo de una placa plana vertical.docx

Efecto del Re sobre el Sh para superficies planas para diferentes Sc

Efecto del Re sobre el Sh para superficies planas para diferentes Sc

Sc=10 Sc=100 Sc=2000
7.00

6.00

5.00

4.00
log (Sh)

3.00
Régimen laminar Régimen Turbulento
2.00

1.00

0.00
0.00 1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00
log (Re)

15
 24/10/2022

Efecto del Re sobre el Sh para superficies planas para diferentes Sc

Régimen laminar Régimen Turbulento

Régimen laminar Régimen Turbulento

Problema
Benceno líquido fluye en una película delgada sobre una placa vertical de 1.5 m de
ancho y 3 m de longitud. La temperatura del líquido es 300 K. Una corriente de
nitrógeno gas (libre de benceno) a 300 K y 1 bar fluye a través del ancho de la placa
paralelo a la superficie a una velocidad de 5 m/s. Calcular la rapidez a la cual el líquido
deba ser alimentado en la punta alta de la placa plana para que la evaporación
enriquezca justo en la parte baja de la placa (yA1= 0.139).
La densidad del benceno líquido a 300 K es 0.88 g/cm3.

16
 24/10/2022

4.3.2 Coeficientes de transferencia de masa para una esfera

La transferencia de masa de esferas individuales es correlacionada con el número de Sherwood mediante

la adición directa de términos que representan la transferencia por difusión puramente molecular y la
transferencia por convección forzada.
• Para valores de Re muy bajo y no hay efecto de convección natural (𝑅𝑒 ≥ 0.4𝐺𝑟𝐷0.5 𝑆𝑐 −1/6 )

𝑆ℎ = 2.0 + 𝑪 𝑅𝑒 𝒎 𝑆𝑐 1Τ3
donde C y m son constantes de la correlación

• Ecuación de Brian and Hales (1969) para transferencia de masa en un fluido líquido:

2/3 1/2
𝑆ℎ = 4 + 1.21𝑃𝑒𝐷 𝑃𝑒𝐷 < 10000

• Para PeD mayor a 10000, Levich (1962) recomienda:

1/3
𝑆ℎ = 1.01𝑃𝑒𝐷 𝑃𝑒𝐷 > 10000

4.3.2 Coeficientes de transferencia de masa para una esfera

• La ecuación de Froessling (1939) y Evnochides y Thodos (1959) correlaciona los datos para
transferencia de masa en gases en el rango de 𝟐 ≤ 𝑹𝒆 ≤ 𝟏𝟐𝟎𝟎𝟎 𝑦 𝟎. 𝟔 < 𝑺𝒄 ≤ 𝟐. 𝟕:

𝑆ℎ = 2.0 + 0.552 𝑅𝑒 1Τ2 𝑆𝑐 1Τ3

Cuando la transferencia ocurre en presencia de convección natural:

• Se utilizan las ecuaciones de Steinberger y Treybal (1960) según:
0.62
𝑆ℎ = 𝑆ℎ𝑐𝑛 + 0.347 𝑅𝑒𝑆𝑐 1Τ2
1 ≤ 𝑅𝑒 ≤ 3x104
0.6 ≤ 𝑆𝑐 ≤ 3200
0.25
𝑆ℎ𝑐𝑛 = 2.0 + 0.569 𝐺𝑟𝐷 𝑆𝑐 𝐺𝑟𝐷 𝑆𝑐 < 108
1/3
𝑆ℎ𝑐𝑛 = 2.0 + 0.0254 𝐺𝑟𝐷 𝑆𝑐 𝑆𝑐 0.244 𝐺𝑟𝐷 𝑆𝑐 > 108
Evaporación de una gota de agua cayendo a través de aire.docx

17
 24/10/2022

Problema: (evaporación de una gota de agua cayendo en aire)

Estimar la distancia a la cual, una gota de agua esférica inicialmente de 1 mm de diámetro

debe caer en un medio de aire seco y en reposo a 323 K y 1 atm para alcanzar un volumen
de 50% al evaporarse. Asumir que la velocidad de la gota es igual a la velocidad terminal a
partir de su diámetro promedio durante el proceso, y que la T del agua permanece a 293K.

2.5.3 Coeficiente de transferencia de masa en un cilindro

Varios investigadores han estudiado la sublimación de un cilindro sólido en aire que fluye
normal al eje. Resultados adicionales en la disolución de un cilindro sólido en una corriente
turbulenta de agua han sido presentados. Bedingfield y Drew (1950) correlacionaron los
datos disponibles en la expresión,
𝑘𝐺 𝑃 𝑆𝐶0.56 −0.4
= 0.281 𝑅𝑒 (45)
𝐺𝑀

400 < 𝑅𝑒 < 25000

0.6 < 𝑆𝑐 < 2.6

Donde,
Re es utilizado con el diámetro del cilindro como longitud característica
GM = v c es la velocidad molar másica del gas
P = presión total

18
 24/10/2022

Problema
Repetir el problema del UF6 utilizando correlaciones de transferencia de masa.
Estimar la rapidez de sublimación de un cilindro de UF6 (M=352), de 1 cm diámetro y
10 cm de longitud expuesto a una corriente de aire que fluye normal al eje del cilindro a
una velocidad de 1 m/s. La temperatura de la superficie del solido es de 303 K, a la cual la
presión de vapor del UF6 es 27 kPa. En el seno del aire la presión y la temperatura son 1
atm y 325 K, respectivamente.
En la siguiente diapositiva se muestra un esquema del problema.

Esquema del
problema Película de
gas estancado
Corriente de aire

v = 1 m/s
P = 1 atm

Sólido UF6 1 NA 2
303 K 325 K
pA= 27 kPa pA= 0 kPa

19
 24/10/2022

4.4 Flujo turbulento en tubos circulares

Se ha estudiado extensivamente la transferencia
de masa de la pared interior de un tubo a un
fluido en movimiento. La mayoría de los datos se
han obtenido de torres de pared húmeda, donde
el líquido volátil fluye dentro de la pared del tubo
circular hacia abajo, mientras que un gas fluye
hacia arriba.

• La medición de la evaporación permite calcular

el coeficiente de transferencia de masa para la
fase gas.

• Usar diferentes gases y líquidos genera

variaciones del Sc.

4.4 Flujo turbulento en tubos circulares

Gilliland y Sherwood (1934) estudiaron la vaporización de nueve diferentes líquidos en
aire, obteniendo la siguiente correlación:
𝑆ℎ = 0.023𝑅𝑒 0.83 𝑆𝑐 0.44

2000 < 𝑅𝑒 < 35000 0.6 < Sc < 2.5

Donde,
 La longitud característica en Sh y Re es el diámetro del tubo y
 Las propiedades del gas son evaluadas en el seno de la corriente de gas a dichas condiciones.

Linton y Sherwood (1950) extendieron el rango del número de Sc,

𝑆ℎ = 0.023𝑅𝑒 0.83 𝑆𝑐 1/3

4000 < 𝑅𝑒 < 70000 0.6 < Sc < 3000

Humidificación de aire en una columna de pared mojada.docx
Información extra: Transferencia de masa en lecho empacado y fluidizado.docx

20
 24/10/2022

Problema
Entrada del
Una corriente de agua fluye hacia abajo al líquido
interior de una columna de pared mojada Recipiente de
cilíndrica mientras una corriente de aire alimentación

seco fluye al interior de la columna hacia Pared de la

tubería
arriba. El diámetro interior de la columna
es de 25.4 mm y la longitud de la sección Película líquida

mojada es de 1.5 m. El aire seco entra con

una velocidad másica de 10 kg/m2s, a 308 K Sección de pared
mojada
y 1 atm de presión. El agua entra a 295 K.
Estimar la presión parcial del agua en el aire
al dejar la columna. Asumir que la fuerza
del aire y del agua permanecen constantes. Sección calmante
Salida del
líquido Gas

𝑘𝐺 𝑎 𝑃 𝑍
𝑝𝐴 𝑆𝑎𝑙𝑖𝑑𝑎 = 𝑷𝑨 1 − 𝑒𝑥𝑝 − (52)
𝐺𝑀

Presión parcial Presión de vapor

6 1 − 𝜀 𝑆ℎ 𝑐 𝑧 𝐷𝐴𝐵
𝑝𝐴 𝑆𝑎𝑙𝑖𝑑𝑎 = 𝑃𝐴 1 − 𝑒𝑥𝑝 − (54)
𝐺𝑀 𝐷2

21
 24/10/2022

2.5.5 Transferencia de masa en lechos empacados y lechos fluidizados

Uso de lechos Adsorción

empacados y Intercambio iónico
fluidizados Cromatografía
Reacciones gaseosas catalizadas en superficies solidas

AUMENTO DE LA SUPERFICIE DE INTERCAMBIO

Correlación de Sherwood et al (1950) para gases,
𝑗𝐷 = 1.17𝑅𝑒 −0.415 (55)
10< 𝑅𝑒 < 2500
𝑑𝑃 𝐺𝑦
Donde, 𝑅𝑒 =
𝜇
Gy = velocidad másica de gas por sección total de área de la columna
dP = diámetro de la esfera de la misma área superficial por unidad de volumen de una
partícula

Para el caso de transferencia de masa en ambos gas y líquido en lecho fijo y

fluidizado de esferas correlacionado por Gupta y Thodos (1962),

0.863
𝜀 𝑗𝐷 = 0.010 + (56)
𝑅𝑒 0.58 − 0.483
1 < 𝑅𝑒 < 2100

Donde ε es la fracción vacía del lecho

El área de transferencia de masa es:

6(1 − 𝜀)
𝑎= (57)
𝑑𝑝

Si el dp ↘ entonces a ↗ , sin embargo la p a través de la columna empacada es

un factor limitante del dp

22
 24/10/2022

Problema
La columna del problema anterior es empacada con un lecho de esferas de vidrio, de
3.5 mm de diámetro. El agua es alimentada por arriba atravesando el lecho empacado a
una velocidad de flujo suficiente para reemplazar el agua perdida por evaporación, y
manteniendo la superficie empacada mojada. Bajo estas condiciones, asumir que se aplica
la ecuación (55).
a) Estimar la profundidad de empaque requerido, si la presión parcial del agua en el aire a
la salida es de 99% de la presión de vapor para el agua en pA.
b) b) Estimar la caída de presión del gas a través del lecho.

Ecuaciones utilizadas para resolver el problema,

6 1 − 𝜀 𝑆ℎ 𝑐 𝑧 𝐷𝐴𝐵
𝑝𝐴 𝑆𝑎𝑙𝑖𝑑𝑎 = 𝑷𝑨 1 − 𝑒𝑥𝑝 − (58)
𝐺𝑀 𝑑𝑝2

Presión parcial Presión de vapor

𝑑𝑝 𝑑𝑝
𝜺 = 0.406 + 0.571 𝑝𝑎𝑟𝑎 < 0.14 (59)
𝐷 𝐷
∆𝑝 𝜀 3 𝑑𝑝 𝜌 150 1 − 𝜀
= + 1.75 (60)
𝑧 1−𝜀 𝐺 2 𝑅𝑒
𝑦

𝑧 𝐺𝑦 2
Δ𝑝 = 2𝑓 (61)
𝐷 𝜌
𝑓 −1/2 = 4 log 𝑅𝑒 𝑓 1/2 − 0.4

23
 24/10/2022

2.5.6 Transferencia de masa en módulos de membrana de fibra hueca

Los módulos de membrana de fibra hueca, son para la transferencia de masa, el equivalente de los
intercambiadores de carcasa y tubos para transferencia de calor.
Separación de gases
Osmosis inversa
Estos módulos son utilizados en los procesos de
Filtración
Diálisis, etc.
Para líquidos atravesando un módulo de fibras se utiliza la correlación para transferencia de masa de
Costello y colaboradores (1993),
−0.47 30 < 100 < 75
1−
𝑆ℎ = 0.53 1 − 1.1 𝑅𝑒 0.53 𝑆𝑐 0.33 (62)
 1−
20 < 𝑅𝑒 < 350

Donde,
2𝑅𝑘𝐿 2 𝑅 v0 𝜌
• 𝑆ℎ = para un modulo de longitud 𝐿 𝑅𝑒 =
𝐷𝐴𝐵 𝜇
• R es el radio de la fibra V0 = velocidad superficial basado en el modulo vacío

•  es la fracción de fibra empacada

Flujo a través del lumen

Las soluciones teóricas han sido desarrolladas para dos condiciones de frontera:
a) Concentración constante en la pared (la más común), y
b) Flujo constante en la pared

La correlación correspondiente es (Cussler, 1997):

1/3
1/3
𝑑𝑖 1/3
𝑆ℎ = 1.62 𝑃𝑒𝐷 = 1.62 𝐺𝑧 (63)
𝐿

Donde,
𝑑𝑖 𝑘𝐿
𝑆ℎ = para un modulo de longitud 𝐿 𝐺𝑧 = 𝑃𝑒𝐷 𝑑𝑖ൗ𝐿 numero de Graetz
𝐷𝐴𝐵
𝑑𝑖 v 𝜌
di es el diámetro interno de la fibra 𝑅𝑒 = donde v es la velocidad promedio del fluido
𝜇

24
 24/10/2022

Problema
El agua caliente de proceso debe ser tratada antes para retirar el aire que causa problemas
de corrosión en los refrescadores. Las membranas de fibra de polipropileno micro poroso
pueden ser utilizadas para una rápida remoción de oxígeno disuelto en agua haciendo uso
de un equipo compacto como se muestra en la siguiente figura.
Para un proceso determinado, se necesita
un equipo para tratar 40,000 kg/h. Diseñar
un equipo de membrana molecular para Agua a la salida
este propósito, asumiendo que la unidad es
capaz de remover el 99% del oxígeno
disuelto presente en el agua a 298 K. Paquete de fibra

Agua alimentada

Solución
El agua fluye a través de la zona central a 298 K, paralelo a las fibras, a v0 = 10 cm/s. El N2
puro a 298 K, 1 atm y con una rapidez de 40 L/min, es usado como gas de separación a
contracorriente fluyendo a través del módulo. El diámetro externo de las fibras es 290 mm,
el factor de empaque es de 40% y el área superficial por unidad de volumen es a = 46.84 m-1.
Cálculo del flujo volumétrico: Cálculo del diámetro de la coraza para el flujo de
agua y la velocidad específica v0 = 10 cm/s :
40000 kgΤh 0.5 0.5
𝑄𝐿 = = 0.0111 m3 Τs 4 𝑄𝐿 4 0.0111 m3 Τs
3600 s 1000 kg 𝐷= = = 0.376 m
1h m3 𝜋 v0 𝜋 0.1 mΤs

Las propiedades de la mezcla de oxígeno diluido en agua a 298 K son 𝝆 = 𝟏𝟎𝟎𝟎 𝐤𝐠Τ𝐦𝟑 , 𝛍 = 𝟎. 𝟗 𝐜𝐏 y
𝐃𝑨𝑩 = 𝟏. 𝟗𝟑𝐱𝟏𝟎−𝟓 𝐜𝐦𝟐 Τ𝐬 . Los números adimensionales son Sc = 𝟒𝟔𝟔 𝑦 Re = 𝟑𝟐. 𝟐𝟐

Sustituyendo en la ecuación (62):

−0.47
1 − 0.4 0.53 x 0.33
Sh = 0.53 1 − 1.1x0.4 x x 32.2 467 = 11.74
0.4

25
 24/10/2022

Calculando kL del lado de la coraza:

11.74 D𝐴𝐵 11.74 x 1.93x10−9 m2 Τs

𝑘𝐿 = = = 0.078𝑥10−3 mΤs
𝑑 290x10−6 m

Para éste caso, la resistencia a la transferencia de masa reside del lado de la coraza. Entonces, el
volumen total del módulo para flujo a contracorriente está dado por:

L 𝑥𝑒𝑛𝑡𝑟𝑎𝑑𝑎
𝑉𝑇 = 𝑙𝑛 1−𝐴 +𝐴 (64)
𝑘𝐿 𝑎 𝑐 1 − 𝐴 𝑥𝑠𝑎𝑙𝑖𝑑𝑎

L
Siendo, 𝐴 = m V (65), donde L y V son los flujos molares de líquido y gas y m es la relación de equilibrio de
concentración de gas en la concentración de líquido. Así:
kg 1h 1 kmol kmol
L = 40000 x x = 0.617 ,
h 3600 s 18 kg s

L 1 min 1 m3 101.3 kPa kmol K

y del gas ideal V = 40 x x x x = 2.726x10−5 kmolΤs
min 60 s 1000 L 298 K 8.314 kPa m3

De la solubilidad del oxígeno en agua a 298 K, m= 4.5x104 a partir la correlación de Davis y Cornwell (1998).
Entonces,
0.617
𝐴= = 0.503
4.5x104 2.726x10−5

Para 99% de remoción de oxígeno disuelto,

𝑥𝑒𝑛𝑡𝑟𝑎𝑑𝑎 1
= 1−0.99 = 100
𝑥𝑠𝑎𝑙𝑖𝑑𝑎

Sustituyendo:

0.617 ln 100 1 − 0.503 + 0.503

V𝑇 = = 0.24 m3
0.0078 46.84 55.5 1 − 0.503

Siendo así, la longitud del módulo es,

V𝑇 0.24 4
𝑍= 2 = 2
= 2.16 m
𝜋𝐷 3.14159 0.376
4

26
 24/10/2022

Difusión molecular mas convección y reacción química

Casos particulares de la ecuación general difusional en estado estacionario

1.- Difusión y reacción química en la frontera
En las reacciones catalíticas a menudo, donde A y B se difunden, la relación entre los
flujos NA y NB en estado estacionario es controlado por la estequiometria de la reacción.
Un ejemplo es A → 2B
𝑑𝑥
Siendo NB= -2NA , Así 𝑁𝐴 = −𝑐𝐷𝐴𝐵 𝑑𝑧𝐴 + 𝑥𝐴 𝑁𝐴 + 𝑁𝐵
𝑑𝑥𝐴
Sustituyendo NB= - 2NA 𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑥𝐴 𝑁𝐴 − 2𝑁𝐴
𝑑𝑧
Re arreglando e integrando,
𝑧2 =𝛿 𝑥𝐴2
𝑑𝑥𝐴
𝑁𝐴 න 𝑑𝑧 = −𝑐𝐷𝐴𝐵 න
𝑧1=0 𝑥𝐴1 1 + 𝑥𝐴

𝑐𝐷𝐴𝐵 1 + 𝑥𝐴1
𝑁𝐴 = 𝑙𝑛 (74)
𝛿 1 + 𝑥𝐴2

Superficie catalítica
NA

NB

z
Z1 Z2 Difusión del componente A en una
reacción heterogénea en la
superficie
Siendo que la “reacción es instantánea”, xA2= 0 , debido a que A no existe en la superficie
catalítica. Así, la ecuación (74) describe la rapidez global del proceso de difusión más que
la reacción química.

27
 24/10/2022

Reacción lenta en la superficie

Si la reacción heterogénea en la superficie no es instantánea pero es lenta para la reacción
A → 2B, y la reacción es de primer orden,
𝑁𝐴𝑧=𝛿 = 𝑘1′ 𝑐𝐴 = 𝑘1′ 𝑐 𝑥𝐴 (75)
Donde,
𝑘1′ es la constante de la rapidez de la reacción heterogénea de primer orden en m/s.
La ecuación (74) se escribe, para las condiciones que 𝑥𝐴2 a z = δ y a partir de la ecuación (75),
𝑁𝐴𝑍=𝛿 𝑁𝐴
𝑥𝐴 = 𝑥𝐴2 = = ′ (76)
𝑘1′ 𝑐 𝑘𝑐
Para el estado estacionario 𝑁𝐴𝑧=𝛿 = 𝑁𝐴 . Sustituyendo la ecuación (76) en la ecuación (74),
𝑐𝐷𝐴𝐵 1 + 𝑥𝐴1
𝑁𝐴 = 𝑙𝑛 (77)
𝛿 𝑁
1 + ′𝐴
𝑘1 𝑐

La rapidez en la ecuación (77) es menor que en la ecuación (74).

Ecuacion de continuidad para mezclas binarias, ****Revisar P en mis notas ****

𝑑𝑝 𝑎 𝑑𝑛𝐴𝑥 𝑑𝑛𝐴𝑦 𝑑𝑛𝐴𝑧
+ + + = 𝑟𝐴 (66)
𝑑𝑡 𝑑𝑥 𝑑𝑦 𝑑𝑧

En notación vectorial:

𝑑𝑝 𝑎
+ 𝛻 ∗ 𝑛𝐴 = 𝑟𝐴 (67)
𝑑𝑡
Dividiendo de ambos lados de la ecuación (66) por MA

𝑑𝐶𝐴 𝑑𝑁𝐴𝑥 𝑑𝑁𝐴𝑦 𝑑𝑁𝐴𝑧

+ + + = 𝑅𝐴 (68)
𝑑𝑡 𝑑𝑥 𝑑𝑦 𝑑𝑧

Donde RA es los kmol de A generados/ 𝑚2 𝑠. Sustituyendo la definición de 𝑁𝐴 = 𝐶𝐴 𝜐𝐴 y

la ley de Fick, 𝑑𝑥𝐴
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝐶𝐴 𝜐𝑀 (69)
𝑑𝑧

Donde 𝜐𝑀 es la velocidad molar promedio

28
 24/10/2022

Escribiendo la ecuación (68) en las tres direcciones, se escribe:

𝑑𝐶𝐴
+ 𝛻 ∗ 𝐶𝐴 𝜐𝑀 − 𝛻 ∗ 𝑐𝐷𝐴𝐵 𝛻𝑥𝐴 = 𝑅𝐴 (70)
𝑑𝑡
Esta es la ecuación general.

Casos particulares de la ecuación de continuidad.

1.- Ecuación para ρ y𝐷𝐴𝐵 constantes

en la difusión de gases la presión total P es a menudo constante. Entonces siendo
𝑃
𝑐 = , c es constante para una temperatura T, A partir de la ecuación general (70)
𝑅𝑇
𝛻𝐶𝐴
y sustituyendo 𝛻𝑥𝐴 =
𝑐

𝑑𝐶𝐴
+ 𝐶𝐴 𝛻 ∗ 𝜐𝑀 + 𝜐𝑀 ∗ 𝛻𝐶𝐴 − 𝐷𝐴𝐵 𝛻 2 𝐶𝐴 = 𝑅𝐴 (71)
𝑑𝑡

2. Contra difusión equimolar para gases

A p= cte sin reacción, c=cte 𝜈𝑀 =0 , 𝐷𝐴𝐵 = 𝑐𝑡𝑒 y 𝑅𝐴 = 0

𝑑𝐶𝐴 𝑑 2 𝐶𝐴 𝑑 2 𝐶𝐴 𝑑 2 𝐶𝐴
= 𝐷𝐴𝐵 + + (72)
𝑑𝑡 𝑑𝑥 2 𝑑𝑦 2 𝑑𝑧 2
Esta es la ecuación para estado no estacionario en líquidos o solidos.

3. Ecuacion para ρ y 𝑫𝑨𝑩 constantes (líquidos)

En soluciones diluidas , ρ y 𝐷𝐴𝐵 son considerados constantes. A partir de la ecuación (67)
𝛻𝜌
sustituyendo 𝑁𝐴 = −𝜌𝐷𝐴𝐵 𝛻𝜔𝐴 + 𝜌𝐴 𝜐 siendo que para ρ=cte, 𝛻𝜔𝐴 = 𝐴 y que
𝜌
𝛻 ∗ 𝜐 = 0, sustituyendo y dividiendo entre la MA, se obtiene:

𝑑𝐶𝐴
+ 𝜐 ∗ 𝛻𝐶𝐴 − 𝐷𝐴𝐵 𝛻 2 𝐶𝐴 = 𝑅𝐴 (73)
𝑑𝑡

29
 24/10/2022

Problema
Un gas puro A se difunde de un punto 1 a un punto 2 a una presión parcial de 101.32 kPa con
un espesor de 2 mm a la pared. En el punto 2 se lleva a cabo una reacción, en la superficie
catalítica, para A→2B.
El compuesto B se difunde de regreso en estado estacionario. La presión total es P=101.32 kPa
2
La temperatura es 300 K y 𝐷𝐴𝐵 = 0.15𝑥10−4 𝑚 Τ𝑠.
Calcular NA y 𝑥𝐴2 para:

a) Rapidez de reacción instantánea

b) Rapidez de reacción lenta donde 𝑘1´ = 5.63𝑥10−3 𝑚Τ𝑠

30