2 Transferencia de Masa Mecanismos de Transferencia 2020

21/03/2021
Curso de Transferencia de Masa
Prof. Francisco Javier Garfias Vásquez

Facultad de Química, UNAM
2. Transferencia de masa molecular
1
21/03/2021
2.1 CONCENTRACIÓN

2.1 Concentración
• En una mezcla multicomponente, la concentración másica (o densidad) de la especie

A, A [kg/m3], corresponde a la masa del componente A por unidad de volumen de
la mezcla.
• La concentración másica total es,
𝜌 = ෍ 𝜌𝑖 (1)
𝑖=1
donde “𝑛” corresponde al número de especies en la mezcla.
2
21/03/2021
Ejemplo 1:
Determinar la concentración másica total de una mezcla de gases contenidos en un
volumen de un metro cúbico.
La composición másica es 6.703 kg de oxígeno gas, 21.864 kg de nitrógeno gas, 37.36

kg de argón gas, 1.382 kg de dióxido de carbono gas y 0.037 kg de neón gas.
Solución:
𝑛
kg kg kg kg kg 𝐤𝐠
𝜌 = ෍ 𝜌𝑖 = 6.703 + 21.864 3 + 37.36 3 + 1.382 3 + 0.037 3 = 𝟔𝟕. 𝟑𝟒𝟔 𝟑
m3 m m m m 𝐦
𝑖=1
13

2.1 Concentración
• La fracción másica es:
𝜌𝐴 𝜌𝐴
𝜔𝐴 = = (2)
𝜌=σ𝑛
𝑖=1 𝜌𝑖 𝜌
• La suma de las fracciones debe ser igual a la unidad,
෍ 𝜔𝑖 = 1 (3)
𝑖=1
14
3
21/03/2021
Ejemplo 2:
• Determinar la fracción másica de los gases contenidos en el volumen de un metro
cúbico del ejemplo anterior.
Solución:
kg
𝜌𝐴 𝜌𝐴 6.703
m3
𝜔𝐴 = = = = 0.0995
𝜌=σ𝑛
𝑖=1 𝜌𝑖 𝜌 kg
67.346 3
m
Gas Concentración másica (kg/m3) Fracción másica
Oxígeno 6.703 0.0995

Nitrógeno 21.864 0.3247
Argón 37.36 0.5548
Dióxido de carbono 1.382 0.0205
Neón 0.037 0.0005
67.346 1.0000
22

2.1 Concentración
La concentración molar de la especie A, es:
𝜌𝐴
𝑐𝐴 = (4)
𝑀𝐴
Donde MA es la masa molecular.

Para un gas ideal,
𝑛𝐴 𝑝𝐴
𝑐𝐴 = = (5)
𝑉 𝑅𝑇
Donde pA corresponde a la presión parcial de la especie A en la mezcla, T es la

temperatura en la escala absoluta y R es la constante de los gases (R = 8.314 J/mol K)
23
4
21/03/2021
Ejemplo 3:
Determinar la concentración molar de cada uno de los gases contenidos en el
volumen de un metro cúbico del ejemplo anterior.
Solución:
kg
𝜌𝐴 6.703 kmol
m3
𝑐𝐴 = = kg = 0.2095
𝑀𝐴 32 m3
kmol
Gas Concentración másica (kg/m3) Fracción másica Concentración molar

(kmol/m3)
Oxígeno 6.703 0.0995 0.2095
Nitrógeno 21.864 0.3247 0.7809
Argón 37.36 0.5548 0.9340
Dióxido de carbono 1.382 0.0205 0.0314
Neón 0.037 0.0005 0.0018
67.346 1.0000
29
2.1 Concentración
La concentración molar total es,
𝑐 = ෍ 𝑐𝑖 (6)
𝑖=1
Para una mezcla gaseosa que obedece la ley del gas ideal,
𝑃
𝑐𝐴 = (7)
𝑅𝑇
donde P es la presión total de la mezcla gaseosa.
30
5
21/03/2021
2.1 Concentración
La fracción molar para mezclas de líquidos o sólidos xA y para mezclas gaseosas yA, son
la concentración de las especies entre la concentración molar total:
𝑐𝐴
𝑥𝐴 = 𝑙í𝑞𝑢𝑖𝑑𝑜𝑠 𝑦 𝑠ó𝑙𝑖𝑑𝑜𝑠 (8𝑎)
𝑐
𝑐𝐴
𝑦𝐴 = 𝑔𝑎𝑠𝑒𝑠 (8𝑏)
𝑐
31
2.1 Concentración
Para gases que obedecen la ley del gas ideal,
𝑐𝐴 𝑝𝐴
𝑦𝐴 = = 𝐿𝑒𝑦 𝑑𝑒 𝐷𝑎𝑙𝑡𝑜𝑛 (9)
𝑐 𝑃
La suma de las fracciones molares es por definición igual a 1,
n n
 y = x
i =1
i
i =1
i =1 (10)(10)
32
6
21/03/2021
33
Ejemplo 4:
Determinar la concentración molar total y la fracción mol de
cada uno de los gases del ejemplo anterior.
Solución:
kmol
𝑐𝐴 0.2095 m3
𝑦𝐴 = = = 0.1070
𝑐 kmol
1.9576
m3
Gas Concentración másica (kg/m3) Fracción Concentración molar Fracción mol
másica (kmol/m3)
Oxígeno 6.703 0.0995 0.2095 0.1070
Nitrógeno 21.864 0.3247 0.7809 0.3989
Argón 37.36 0.5548 0.9340 0.4771
Dióxido de carbono 1.382 0.0205 0.0314 0.0160
Neón 0.037 0.0005 0.0018 0.0010
67.346 1.0000 1.9576 1.0000
39
7
21/03/2021
Problema 1: Concentración de la alimentación en un absorbedor de gas

Un gas que contiene 88% (en volumen) de metano CH4, 4% de etano C2H6, 5% de n-propano n-C3H8 y 3% de
n-butano n-C4H10 a 300 K y 500 kPa será absorbido por contacto con un aceite no volátil en un absorbedor de
gases. El objetivo en el proceso es recuperar, en el efluente líquido, lo más posible de hidrocarburo pesado.
Calcular: a) La concentración molar total del gas en la alimentación
b) La densidad del gas a la entrada
c) La composición del gas a la entrada, expresado en fracción masa
Aceite no volátil Gas de salida
Producto líquido Gas alimentado

88% CH4, 4% C2H6, 5% n-C3H8,
3% n-C4H10, 300 K, 500 kPa
40

Solución:
a) La concentración molar total del gas en la alimentación
P 500 kPa 𝐤𝐦𝐨𝐥

c= = = 𝟎. 𝟐
RT 8.314 kPa m3 Τkmol K 300 K 𝐦𝟑
45
8
21/03/2021

Solución:
b) La densidad del gas a la entrada. Considerando una base de 100 kmol de la mezcla de gases.
1926.61 kg 𝐤𝐠 kmol kg 𝐤𝐠
Mprom = = 𝟏𝟗. 𝟐𝟕 ρ = c Mprom = 0.2 19.27 = 𝟑. 𝟖𝟓 𝟑
100 kmol 𝐤𝐦𝐨𝐥 m3 kmol 𝐦
Gas Base: 100 kmol de Masa molecular Masa (kg)

mezcla gaseosa (kg/kmol)
CH4 88 16.04 1411.52
C2H6 4 30.07 120.28
n-C3H8 5 44.09 220.45
n-C4H10 3 58.12 174.36
Total 100 1926.61
53

Solución:
c) La composición del gas a la entrada, expresado en fracción masa.
1411. 52 kg
ωA = = 𝟎. 𝟕𝟑𝟑
1926.61 kg
Gas Base: 100 kmol Masa molecular Masa (kg) Fracción masa
de mezcla (kg/kmol)
gaseosa
CH4 88 16.04 1411.52 0.733
C2H6 4 30.07 120.28 0.062
n-C3H8 5 44.09 220.45 0.114
n-C4H10 3 58.12 174.36 0.091
Total 100 1926.61 1.000
57
9
21/03/2021
58
Problema 2: Concentración de Nitrato de Potasio en una solución de lavado
• En la fabricación de nitrato de potasio, el cloruro de potasio reacciona con una solución caliente de
nitrato de sodio:
KCl + NaNO3 → KNO3 + NaCl
• La mezcla reaccionante se enfría a 293 K para formar cristales puros de KNO3. La solubilidad en el
equilibrio del KNO3 en agua a 293 K es de 24% (en peso), y la densidad de la solución saturada es de
1162 kg/m3. Calcular:
a) La densidad molar de la solución de lavado
b) La composición de la solución expresada en fracción mol
KNO3 + NaCl Cristales lavados
Solución de lavado Solución fresca

24% de KNO3
1162 kg/m3
59
10
21/03/2021
nitrato de sodio:
Usar Base: 100 kg de solución de lavado
Componente Masa (kg) Masa molecular Mol (kmol)

(kg/kmol)
KNO3 24 101.10 0.237
H20 76 18.02 4.218
Total 100 4.455
76 kg
𝑛B = = 𝟒. 𝟐𝟏𝟖 𝐤𝐦𝐨𝐥
18.02 kgΤkmol
72
nitrato de sodio:
Componente Masa (kg) Masa molecular Mol (kmol)
(kg/kmol)
KNO3 24 101.10 0.237
H20 76 18.02 4.218
Total 100 4.455
mtotal 100 kg
Mprom = = = 𝟐𝟐. 𝟒𝟓 𝐤𝐠Τ𝐤𝐦𝐨𝐥
ntotal 4.455 kmol
ρ 1162 kgΤm3
c= = = 𝟓𝟏. 𝟕𝟕 𝐤𝐦𝐨𝐥Τ𝐦𝟑
Mprom 22.45 kgΤkmol
79
11
21/03/2021
nitrato de sodio:
b) La composición de la solución expresada en fracción mol
Componente Mol (kmol) Fracción mol
KNO3 0.237 (0.237/4.455) = 0.053

H20 4.218 (4.218/4.455) = 0.947
Total 4.455 1.000
83
2.2 Velocidad y flux

• La relación empírica básica para estimar la rapidez de difusión molecular, postulada por Fick en 1855 y
referida como la Primera Ley de Fick, cuantifica la difusión del componente o especie A en un sistema
isotérmico e isobárico.
• De acuerdo con la ley de Fick, las especies tienen una velocidad relativa a la velocidad promedio másica o
molar (llamada velocidad de difusión), sólo si existe un gradiente de concentración.
Difusión
84
12
21/03/2021
2.2 VELOCIDAD Y FLUX
85
• En un sistema multicomponente las diversas especies se moverán a
diferentes velocidades con respecto a la velocidad molar o másica
promedio.
• Por ésta razón, una evaluación de la velocidad característica para una
mezcla de gases requerirá del promedio de las velocidades de cada
una de las especies presentes.
86
13
21/03/2021

• Ejemplo de la velocidad relativa a una especie en particular con relación al
movimiento local o velocidad de difusión.
87
• La velocidad másica promedio para una mezcla multicomponente es

definida en términos de su densidad,
n n
 i vi  v i i n
v= i =1
= i =1
=  i vi (11)
(11)
n


i =1
i
i =1
• Donde vi es la velocidad absoluta de la especie 𝑖 relativa al eje de

coordenadas estacionarias. La velocidad v es la velocidad medida con un
tubo pitot.
88
14
21/03/2021
Tubo pitot
𝑝2
𝑝1
∆𝑝 =
89
Problema 3: Velocidad másica promedio

Del problema de un gas que contiene 88% (en volumen) de metano CH4, 4% de etano C2H6, 5% de n-propano
n-C3H8 y 3% de n-butano n-C4H10 a 300 K y 500 kPa, calcular la velocidad másica promedio.
Solución:
Gas Base: 100 Masa molecular Masa Fracción Velocidad
kmol de (kg/kmol) (kg) masa cuadrática
mezcla media
gaseosa ucm (m/s)
CH4 88 16.04 1411.52 0.733 683
C2H6 4 30.07 120.28 0.062 499
n-C3H8 5 44.09 220.45 0.114 412
n-C4H10 3 58.12 174.36 0.091 359
Total 100 1926.61 1.000
Primero, calculamos la velocidad cuadrática media,
J Kg m2 Τs2 Kg m2 Τs2
3𝑅𝑇 3 8.314 300 K 3 8.314 300 K 3 8.314 300 K 𝐦
mol K mol K mol K
ucm = = = = = 𝟔𝟖𝟑
𝑀 0.01604 kg/mol 0.01604 kg/mol 0.01604 kg/mol 𝐬
98
15
21/03/2021

n-C3H8 y 3% de n-butano n-C4H10 a 300 K y 500 kPa, calcular la velocidad másica promedio.
Solución:
mezcla media
gaseosa ucm (m/s)
CH4 88 16.04 1411.52 0.733 683
C2H6 4 30.07 120.28 0.062 499
n-C3H8 5 44.09 220.45 0.114 412
n-C4H10 3 58.12 174.36 0.091 359
Total 100 1926.61 1.000
v = ෍ ωi vi = 0.733x683 + 0.062x499 + 0.114x412 + 0.091x359 = 𝟔𝟏𝟏. 𝟐 𝐦/𝐬

i=1
101

• La velocidad molar promedio,
n n
c v c v i i i i n
V= i =1
n
= i =1
=  xi vi (12)
(12)
c
c i =1
i
i =1
• v𝑖 − v, la velocidad de difusión de la especie i con respecto a la velocidad

másica promedio
• v𝑖 − V, la velocidad de difusión de la especie i con respecto a la velocidad

molar promedio
102
16
21/03/2021

n-C3H8 y 3% de n-butano n-C4H10 a 300 K y 500 kPa, calcular la velocidad molar promedio.
Solución:
mezcla media
gaseosa ucm (m/s)
CH4 88 16.04 1411.52 0.733 683
C2H6 4 30.07 120.28 0.062 499
n-C3H8 5 44.09 220.45 0.114 412
n-C4H10 3 58.12 174.36 0.091 359
Total 100 1926.61 1.000
V = ෍ yi vi = 0.88x683 + 0.04x499 + 0.05x412 + 0.03x359 = 𝟔𝟓𝟐. 𝟒 𝐦/𝐬

i=1
108
2.2 Velocidad y flux (densidad de flujo)
El flux másico o molar de una especie es una cantidad vectorial, en unidades de masa o
moles, que atraviesa por un área unitaria normal al vector por unidad de tiempo.
Se define con referencia a:
✓ coordenadas fijas en el espacio,
✓ coordenadas que se mueven con la velocidad másica o molar promedio.
El flux másico de i con respecto a las coordenadas fijas en el espacio:

ni = ρi vi (13)
Si sumamos el flux de los componentes, se obtiene el flux másico total,
n=ρv (14)
109
17
21/03/2021
Problema 5: Flux másico total

n-C3H8 y 3% de n-butano n-C4H10 a 300 K y 500 kPa, calcular el flux másico total.
Solución:
mezcla media
gaseosa ucm (m/s)
CH4 88 16.04 1411.52 0.733 683
C2H6 4 30.07 120.28 0.062 499
n-C3H8 5 44.09 220.45 0.114 412
n-C4H10 3 58.12 174.36 0.091 359
Total 100 1926.61 1.000
kg m 𝐤𝐠
n = ρ v = 3.85 611.2 = 𝟐𝟑𝟓𝟐 𝟐
m3 s 𝐦 𝐬
116
Efecto de las
corrientes de
convección en la
transferencia de
masa
117
18
21/03/2021
• El flux molar de la especie 𝑖 es,

Ni = ci vi (15)
• El flux molar total,
n
N = ෍ ci vi o N=cV (16)
i=1
• El flux másico de difusión de la especie 𝑖, con respecto a la velocidad

másica promedio:
n
ji = ρi vi − v y ෍ ji = 0 (17)
i=1
118
• El flux molar de difusión de la especie i, con respecto a la velocidad molar promedio:

n
J i = ci ( vi − V ) y J
i =1
i =0 (18)
(18)
• El flux másico ni esta relacionado con el flux de difusión másico como,
n i = ji + i v = ji + i n (19)
(19)
• El flux molar Ni esta relacionado con el flux de difusión molar como,
N i = J i + ci V = J i + yi N (20)
(20)
119
19
21/03/2021
120
Problema 5: Flux másico total

n-C3H8 y 3% de n-butano n-C4H10 a 300 K y 500 kPa, calcular la velocidad molar promedio.
Solución:
mezcla media
gaseosa ucm (m/s)
CH4 88 16.04 1411.52 0.733 683
C2H6 4 30.07 120.28 0.062 499
n-C3H8 5 44.09 220.45 0.114 412
n-C4H10 3 58.12 174.36 0.091 359
Total 100 1926.61 1.000
kg m 𝐤𝐠 kmol m 𝐤𝐦𝐨𝐥
n = ρ v = 3.85 611.2 = 𝟐𝟑𝟓𝟐 𝟐 N = c V = 0.2 652.4 = 𝟏𝟑𝟎. 𝟒𝟖 𝟐
m3 s 𝐦 𝐬 m3 s 𝐦 𝐬
124
20
21/03/2021
N i = J i + ci V = J i + yi N (20)
(20)
• Es importante hacer notar que el flux molar Ni , descrito por la ecuación (20) es el
resultado de dos cantidades vectoriales:
Ji El flux molar de difusión resultado de la contribución por

gradiente de concentración
yi N = ci V El flux molar resultado del arrastre del componente i por el seno del
fluido; este término se refiere a la contribución por el movimiento del fluido
• Una o ambas cantidades tienen un efecto importante sobre el total del flux molar Ni.
125
El flux molar resultado del

arrastre del componente i
por el movimiento del
fluido
126
21
21/03/2021
2.3 RELACIONES DE MAXWELL-STEFAN
127
2.3 Relaciones de Maxwell – Stefan
• Maxwell – Stefan (1870)
James Clerk Maxwell Josef Stefan
128
22
21/03/2021
• Objetivos:
– Escribir las ecuaciones de Maxwell-Stefan (MS) para un sistema binario

y para sistemas multicomponentes
– Definir los conceptos de difusividad y factor termodinámico de MS
– Expresar la fuerza motriz para la transferencia de masa en términos de

gradiente de fracción molar para sistemas ideales y no ideales
129

• Considerar una mezcla binaria de gases ideales 1 y 2 a temperatura y
presión constante. De un balance de momentum o cantidad de movimiento,
describiendo las colisiones entre moléculas de especies 1 y de especies 2, se
obtiene (Taylor and Krishna, 1993):
p1 = − f12 y1 y2 ( v1 − v 2 ) (21)
• Donde f12 es un parámetro análogo al factor de fricción ó coeficiente de

arrastre.
130
23
21/03/2021
• Por conveniencia, definimos un coeficiente de arrastre inverso 𝑫12 = 𝑃Τ𝑓12

y rescribiendo la ecuación (21),
𝛻𝑝1 𝑦1 𝑦2 v1 −v2
𝑑1 = =− (22)
𝑃 𝑫12
• Donde , d1 = (1 P ) p1 es la fuerza motriz para la difusión de la especie 1 en

una mezcla de gas ideal a temperatura y presión constante.
131
• La ecuación (22) es la Ecuación de Maxwell – Stefan para la difusión de la

especie 1 en una mezcla binaria ideal de gases.
• El símbolo 𝑫12 es la difusividad de Maxwell – Stefan (MS).
• De un análisis similar para la especie 2, se obtiene

𝛻𝑝2 𝑦1 𝑦2 v2 −v1
𝑑2 = =− (23)
𝑃 𝑫21
• Para todas las aplicaciones consideradas, la presión del sistema es

constante a través del camino de difusión.
132
24
21/03/2021
• Así, las ecuaciones (22) y (23) se simplifican a,
𝑦1 𝑦2 v1 −v2 𝑦1 𝑦2 v2 −v1
𝛻𝑦1 = − y 𝛻𝑦2 = − (24)
𝑫12 𝑫21
• Para mezclas multicomponentes, la ecuación (22) se generaliza,
𝑛
𝑦𝑖 𝑦𝑗 v𝑖 − v𝑗
d𝑖 = − ෍ i = 1, 2, …, n−1 (25)
𝑫𝑖𝑗
𝑗=1
133
• La ecuación (25) también se escribe en términos de flux molar, N i = ci vi ,

𝑛
𝑦𝑖 N𝑗 − 𝑦𝑗 N𝑖
d𝑖 = − ෍ i =1, 2, …, n−1 (26)
𝑐 𝑫𝑖𝑗
𝑗=1
• En términos de flux de difusión Ji,

𝑛
𝑦𝑖 J𝑗 − 𝑦𝑗 J𝑖
d𝑖 = − ෍ i =1, 2, …, n−1 (27)
𝑐 𝑫𝑖𝑗
𝑗=1
• Estas son las ecuaciones de difusión de Maxwell – Stefan (1870) para

sistemas multicomponentes.
134
25
21/03/2021
• También, para una mezcla de gases ideales multicomponente, se cumple

que,
𝑫𝑖𝑗 = 𝑫𝑗𝑖 (28)
• Para una mezcla binaria de gases ideales, donde la fuerza motriz para la
difusión es el gradiente de fracción mol, la ecuación (27) se reduce a,
J1 = − 𝑐𝑫12 d1 = − 𝑐𝑫12 𝛻𝑦1 (29)

•
135
• Para fluidos no ideales la fuerza motriz para la difusión debe estar

definida en términos del gradiente de potencial químico,
x
d𝑖 = RTi ∇T,P μi (30)
• Los subíndices T, P son para enfatizar que el gradiente debe ser calculado
bajo condiciones de temperatura y presión constante. Podemos expresar
la ecuación (30) en términos de los coeficientes de actividad, i, y el
gradiente de fracción mol como,
𝑛−1
di = − ෍ Γ𝑖𝑗 ∇xj = 𝚪∇xi (31)

𝑗=1
136
26
21/03/2021
• Donde la matriz del factor termodinámico  esta dado por,
 ln γ i
i j = δi j + xi (32)
x j
T , P ,
• Donde ij es la delta de Kronecker definida como,

δi j =1 si i = j
δi j =0 si i  j (33)
• El símbolo  es usado en la ecuación (32) para indicar que la diferenciación con respecto
a la fracción mol será llevada a cabo de tal manera que la suma de la fracción mol
siempre será igual a 1.0.
137
• Combinando las ecuaciones (26), (27) y (31),

𝑛 𝑛
𝑦𝑖 N𝑗 − 𝑦𝑗 N𝑖 𝑦𝑖 J𝑗 − 𝑦𝑗 J𝑖
d𝑖 = − ෍ i =1, 2, …, n−1 (26) d𝑖 = − ෍ i =1, 2, …, n−1 (27)
𝑐 𝑫𝑖𝑗 𝑐 𝑫𝑖𝑗
𝑗=1 𝑗=1
𝑛−1
di = − ෍ Γ𝑖𝑗 ∇xj = 𝚪∇xi (31)

𝑗=1
• Se obtiene,
𝑛 𝑛
𝑥𝑖 J𝑗 − 𝑥𝑗 J𝑖 𝑥𝑖 N𝑗 − 𝑥𝑗 N𝑖
Γ∇𝑥𝑖 = ෍ =෍ i =1, 2, …, n−1 (34)
𝑗=1 𝑗=1
138
27
21/03/2021

• Para una mezcla multicomponente de gases ideales, la ecuación (34) es
𝑛 𝑛
𝑦𝑖 J𝑗 − 𝑦𝑗 J𝑖 𝑦𝑖 N𝑗 − 𝑦𝑗 N𝑖
Γ∇𝑦𝑖 = − ෍ =෍ i =1, 2, …, n−1 (35)
𝑗=1 𝑗=1
• Para una mezcla binaria, la ecuación (34) se reduce a,
J1 = − 𝑐𝑫12 Γ ∇𝑥1 (36)
• Donde el factor termodinámico  esta dado por,
 ln γ1
 = 1 + x1 (37)
x1
139
• El factor termodinámico es evaluado para mezclas de líquidos a partir de

modelos de coeficiente de actividad. Para una solución regular, por
ejemplo,
ln γ1 = Ax22 = A (1 − x1 )
2
(38)
• Así, la ecuación (37) llegaría a ser,

 = 1 − 2 Ax1 x2 (39)
140
28
21/03/2021
2.4 PRIMERA LEY DE FICK PARA MEZCLAS

BINARIAS
141
2.4 Primera ley de Fick para mezclas binarias

• En la misma época que Maxwell y Stefan desarrollaban sus ideas sobre
difusión en mezclas multicomponentes, Adolf Fick estaba por descubrir las
ecuaciones básicas para la difusión a través de estudios experimentales de
mezclas binarias (Fick, 1855).
• El resultado del trabajo de Fick fue la ley que lleva su nombre para mezclas
binarias en un sistema isotérmico e isobárico, esta es,
J1 = − 𝑐D12 ∇𝑥1 (40)
• donde D12 es la difusividad o coeficiente de difusión.
142
29
21/03/2021
• Comparando las ecuaciones (36) y (40),
J1 = − 𝑐𝑫12 Γ ∇𝑥1 (36)
J1 = − 𝑐D12 ∇𝑥1 (40)
• vemos que, para un sistema binario, la difusividad de Fick D12 y la

difusividad de Maxwell – Stefan (MS) D12 están relacionadas por,
D12 = 𝑫12 Γ (41)
143

• Para sistemas ideales,  es igual a 1.0 y las difusividades son idénticas,
D12 = 𝑫12 para sistemas ideales (42)
 ln γ1
• Donde el factor termodinámico  esta dado por,  = 1 + x1 (37)
x1
• La difusividad de Fick incorpora dos aspectos (está implícito):
1. El significado de la inversa del coeficiente de arrastre 𝑫12 , y
2. La no idealidad termodinámica (  ).
• Consecuentemente, la interpretación física del coeficiente de difusión de Fick es menos

transparente que la difusividad de Maxwell – Stefan (MS).
144
30
21/03/2021
• Para un sistema binario, la difusividad de Fick D y la difusividad de

Maxwell-Stefan D, están relacionadas por,
D𝐴𝐵 = 𝐷𝐴𝐵 Γ
• Donde el factor termodinámico  esta dado por,
 ln γ1
 = 1 + x1 (37)
x1
• Para sistemas ideales, 𝚪 es igual a la unidad y las difusividades son idénticas,
D𝐴𝐵 = 𝐷𝐴𝐵
145
• La difusividad de Fick incorpora dos aspectos:
➢ El significado de un coeficiente de arrastre inverso (D), y
➢ La no-idealidad termodinámica (𝚪)
• En consecuencia, la interpretación física del coeficiente de difusión de Fick es
menos transparente que para el caso de la difusividad de Maxwell-Stefan.
146
31
21/03/2021
3. El coeficiente de difusión
147
• El factor de proporcionalidad de la primera ley de Fick se le conoce como coeficiente

de difusión o difusividad (D𝐴𝐵 ). Sus dimensiones son,
𝐽𝐴 M 1 𝐋𝟐
D𝐴𝐵 =− = 2 =
𝑐 ∇𝑥𝐴 Lt MΤL3 x 1ΤL 𝐭
• Son idénticas a las dimensiones fundamentales de otras propiedades de transporte:

❖ viscosidad cinemática 𝜈, y
𝜅
❖ difusividad térmica 𝛼 , o su razón equivalente
𝜌𝐶𝑝
148
32
21/03/2021
El coeficiente de difusión depende de la presión, la temperatura y la
composición del sistema, siendo generalmente los valores,
• mayores para gases entre 0.5 y 1x10-5 m2/s,
• menores para líquidos entre 10-10 a 10-9 m2/s y
• mucho menores para sólidos entre 10-14 a 10-10 m2/s
149
3.1 Difusividad en gases (a baja presión)
• La viscosidad (transferencia de momentum) y sólo dependen de la presión

y la temperatura
• la conductividad calorífica (transferencia de calor)
• La difusividad depende de la presión, de la temperatura y de la composición del sistema
• En mezclas binarias de gases a baja presión, la difusividad 𝑫𝑨𝑩 es

• inversamente proporcional a la presión (𝑫𝑨𝑩  𝟏/𝑷) ,
• proporcional a la temperatura (𝑫𝑨𝑩  𝑻𝟏.𝟓 ), y
• casi independiente de la composición
150
33
21/03/2021
151
152
34
21/03/2021
153
3.1 DIFUSIVIDAD PARA SISTEMAS DE GASES

IDEALES BINARIOS
154
35
21/03/2021
3.1 Difusividad para sistemas de gases ideales binarios (a baja presión)
• Al combinar la teoría cinética y los conceptos de los estados correspondientes, se obtiene

la D𝐴𝐵 (cm2/s), a baja presión:
𝑏
𝑃 D𝐴𝐵 𝑇
1Τ2
=𝑎
1Τ3 5Τ12 1 1 𝑇𝑐𝐴 𝑇𝑐𝐵
𝑝𝑐𝐴 𝑝𝑐𝐵 𝑇𝑐𝐴 𝑇𝑐𝐵 𝑀𝐴 + 𝑀𝐵
Mezclas binarias de gases no polares Agua con gas no polar

(excluir He e H2)
𝑎 = 2.745x10−4 𝑎 = 3.640x10−4
𝑏 = 1.823 𝑏 = 2.334
155
Problema 6: Difusividad en gases (a baja presión)

Estimar el coeficiente de difusión a 1 atm y 293 K para el sistema argón – oxígeno. Los valores de temperatura,
presión crítica y masa para el argón y el oxígeno son:
𝑻𝒄 (K) 𝑷𝒄 (atm) M, g/gmol
Ar 151.2 48 40
O2 154.4 49.7 32
Solución:
1Τ2
1Τ3 5Τ12 1 1
𝑝𝑐𝐴 𝑝𝑐𝐵 = 13.362 𝑇𝑐𝐴 𝑇𝑐𝐵 = 66.082 + = 0.237
𝑀𝐴 𝑀𝐵
𝑏 1.823
𝑇 𝑇
𝑎 = 2.745x10−4 = 8.996x10−4
𝑇𝑐𝐴 𝑇𝑐𝐵 𝑇𝑐𝐴 𝑇𝑐𝐵
8.996𝑥10−4 13.362 66.082 0.237 𝟐
D𝐴𝐵 = = 𝟎. 𝟏𝟖𝟖 𝐜𝐦 Τ𝐬
1
DAB𝑒𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙 − DAB𝑡𝑒ó𝑟𝑖𝑐𝑎 0.200 − 0.188
%Error = x 100% = x 100% = 𝟔 %
DAB𝑒𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙 0.200
156
36
21/03/2021
3.1 Difusividad para sistemas de gases ideales binarios
Una de las expresiones teóricas para el coeficiente de difusión en mezcla de gases de baja
densidad, se deriva de la teoría cinética propuesta por Sutherland, Jeans, Chapman y
Cowling.
• Es el modelo dinámico de gases más simple, considera que
• las moléculas son esferas rígidas,
• no hay fuerzas intermoleculares,
• las colisiones son completamente elásticas y
• diámetros diferentes para una mezcla binaria de gases
157
3.1 Difusividad para sistemas de gases ideales binarios
• El modelo para determinar la difusividad es:
1 1 1ൗ2
2𝑘 3/2𝑁𝐴 𝑇 3/2 (2𝑀𝐴 + 2𝑀𝐵 )
1/2 Sutherland, Jeans,
D𝐴𝐵 = Chapman y Cowling
3𝜋 3/2 𝜎 +𝜎 2
𝑃 𝐴2 𝐵
𝒌 Constante de Boltzmann (1.3806𝑥10−16 𝑒𝑟𝑔𝑖𝑜𝑠ൗK)

𝝈𝒊 Diámetro de Lennard-Jones de moléculas esféricas
𝑵𝑨 Número de Avogadro
P en atm y T en Kelvin
En el caso de los coeficientes de difusión en fase gas, se considera 𝐃𝑨𝑩 = 𝐃𝑩𝑨 ,

no es el mismo caso en fase líquida.
158
37
21/03/2021
3.1 Difusividad para sistemas de gases ideales binarios (bajas a moderadas

presiones)
Versiones más modernas se basan en las fuerzas de atracción y repulsión entre las
moléculas. Hirschfelder y col. (1949), usando el potencial de Lennard-Jones, presentan
una ecuación para el coeficiente de difusión para pares de gases no polares y moléculas
no reactivas:
1Τ2
1 1
0.001858 𝑇1.5 𝑀 + 𝑀 Hirschfelder y col. (1949)
𝐴 𝐵
𝐷𝐴𝐵 = 2
𝑃 𝜎𝐴𝐵 Ω𝐷
𝜴𝑫 Integral de colisión, que es una función de 𝑘𝑇Τ𝜀𝐴𝐵

𝝈𝑨𝑩 Es el diámetro de colisión, parámetro de Lennard-Jonnes en Å
𝜺𝑨𝑩 Es la energía de interacción molecular para el sistema binario A y B, un parámetro de Lennard-Jones en ergs
DAB es el coeficiente de difusión, en cm2/s
MA y MB es la masa molecular de A y B
T es la temperatura en la escala absoluta Kelvin
P es la presión en atm
159

presiones)
𝑎 𝑐 𝑒 𝑔
Ω𝐷 = + + +
𝑇∗ 𝑏 exp 𝑑𝑇 ∗ exp 𝑓𝑇 ∗ exp ℎ𝑇 ∗
𝑎 = 1.06036 𝑏 = 0.15610 𝑐 = 0.193 𝑑 = 0.47635

𝑒 = 1.03587 𝑓 = 1.52996 𝑔 = 1.76474 ℎ = 3.89411
𝑘𝑇
Ω𝐷 , es la Integral de colisión para pares de gases no polares, donde 𝑇∗ = 𝜀
𝐴𝐵
Para un sistema binario, compuesto por pares de moléculas no polares, los parámetros de
Lennard-Jones de componentes puros pueden ser combinados empíricamente por las
siguientes relaciones:
𝜎𝐴 + 𝜎𝐵
𝜎𝐴𝐵 = 𝜀𝐴𝐵 = 𝜀𝐴 𝜀𝐵
2
160
38
21/03/2021
Integral de colision difusional en función de la temperatura adimensional
161

presiones)
En ausencia de datos experimentales, los valores para componentes puros pueden ser
estimados por las siguientes relaciones empíricas:
1ൗ
1ൗ 1ൗ 𝑇𝑐 3
𝜎= 1.18𝑉𝑏 3 𝜎= 0.841𝑉𝑐 3 𝜎 = 2.44
𝑃𝑐
𝜀𝐴 𝜀𝐴
ൗ𝑘 = 0.77𝑇𝑐 ൗ𝑘 = 1.15𝑇𝑏 𝑉𝑏 = 0.285𝑉𝑐1.048
Método de Tyn
y Calus (1975)
3
𝑉𝑏 es el volumen molar al punto normal de ebullición 𝑐𝑚 ൗ𝑔 𝑚𝑜𝑙
3
𝑉𝑐 es el volumen molar crítico 𝑐𝑚 ൗ𝑔 𝑚𝑜𝑙
𝑇𝑏 es la temperatura normal de ebullición en K
𝑃𝑐 es la presión crítica en bar
162
39
21/03/2021

presiones)
• La ecuación de Hirschfelder es usada normalmente para extrapolar datos para rangos
moderados de presión hasta 25 bar.
• Se puede predecir también el coeficiente de difusión a cualquier temperatura en

cualquier presión por debajo de 25 bar, a partir de un valor experimental conocido
aplicando la ecuación:
3ൗ
𝑃1 𝑇2 2 Ω𝐷|𝑇1
DAB𝑇2,𝑃2 = D𝐴𝐵𝑇1,𝑃1
𝑃2 𝑇1 Ω𝐷|𝑇2
163
164
40
21/03/2021
Problema 6: Difusividad de gases ideales

Calcular la difusividad del dióxido de carbono en aire a 20°C y 1 atmósfera usando la ecuación de Hirschfelder.
Comparar con el valor experimental DAB = 0.136 cm2/s a 0°C.
Solución:
Los valores de los parámetros de Lennard-Jones son:
𝝈 𝒆𝒏 Å 𝜺/𝒌 𝒆𝒏 𝑲 M, g/gmol
CO2 3.996 190 44
Aire 3.617 97 29
𝜎𝐴 + 𝜎𝐵 3.996 + 3.617 Diámetro de colisión, parámetro de Lennard-

𝜎𝐴𝐵 = = = 𝟑. 𝟖𝟎𝟕 Å
2 2 Jonnes
𝜀𝐴𝐵 𝜀𝐴 𝜀𝐵 Energía de interacción molecular para el sistema

= = 190 x 97 = 𝟏𝟑𝟓. 𝟖 𝐊 binario A en B (parámetro de Lennard-Jones en
𝑘 𝑘 𝑘
ergs /k) en Kelvin
𝑘𝑇 293
𝑇∗ = = = 𝟐. 𝟏𝟓𝟖 para 𝐓 = 𝟐𝟗𝟑 𝐊 (𝟐𝟎°𝐂)
𝜀𝐴𝐵 135.8
165

Comparar con el valor experimental DAB = 0.136 cm2/s a 0°C (Welty, J.R., 1984).
Solución:
1.06036 0.193 1.03587 1.76474

Ω𝐷 = + + + = 𝟏. 𝟎𝟒𝟖
2.1580.1561 exp(0.47635 x2.158) exp(1.52996 x2.158) exp(3.899411 x2.158)
1Τ2
1 1
0.001858 𝑇1.5 𝑀 + 𝑀
𝐴 𝐵
𝐷𝐴𝐵 = 2
𝑃 𝜎𝐴𝐵 Ω𝐷
0.001858 293 1.5 0.0572 0.5 𝐜𝐦𝟐 𝐦𝟐

D𝐴𝐵 = 2 = 𝟎. 𝟏𝟒𝟕 = 𝟏. 𝟒𝟕𝐱𝟏𝟎−𝟓
1 3.807 1.048 𝐬 𝐬
3ൗ
𝑃1 𝑇2 2 Ω𝐷|𝑇1
Para comparar, utilicemos la ecuación a diferente temperatura, DAB𝑇2,𝑃2 = D𝐴𝐵𝑇1,𝑃1
𝑃2 𝑇1 Ω𝐷|𝑇2
166
41
21/03/2021

Comparar con el valor experimental DAB = 0.136 cm2/s a 0°C (Welty, J.R., 1984).
Solución:
En el punto 1: T1 = 273 K, P1 = 1 atm, Ω𝐷|𝑇1 = 1.073 y DAB = 0.136 cm2/s

En el punto 2: T2 = 293 K, P2 = 1 atm, Ω𝐷|𝑇2 = 1.047 y 𝐃𝑨𝑩 = ? 𝐜𝐦𝟐/𝐬
3ൗ 3ൗ
𝑃1 𝑇2 2 Ω𝐷|𝑇1 293 2 1.073 𝟐
DAB𝑇2,𝑃1 = D𝐴𝐵𝑇1 ,𝑃1 = 0.136 = 𝟎. 𝟏𝟓𝟓 𝐜𝐦 Τ𝐬 a 293 K (20°C)
𝑃2 𝑇1 Ω𝐷|𝑇2 273 1.048

%Error = x 100% = x 100% = 𝟓. 𝟐 %
167
168
42
21/03/2021
3.1 Difusividad para sistemas de gases ideales binarios (bajas a moderadas presiones)
Brokaw propuso un método para estimar la difusividad para mezclas binarias de gases
que contienen compuestos polares usando la ecuación de Hirschfelder, evaluando para
la integral de colisión con:
2
0.196𝛿𝐴𝐵 1/2
1.94𝑥103 𝜇𝑃2
Ω𝐷 = Ω𝐷0 + 𝛿𝐴𝐵 = 𝛿𝐴 𝛿𝐵 𝛿=
𝑇∗ 𝑉𝑏 𝑇𝑏
𝜇𝑃 es el momento dipolar, en Debyes.

𝑐𝑚3
𝑉𝑏 es el volumen molar del componente líquido a su temperatura normal de ebullición, 𝑔 𝑚𝑜𝑙
El diámetro de colisión y la energía de interacción se evalúa con,
1ൗ
1.585𝑉𝑏 3
1
𝜎𝐴𝐵 = 𝜎𝐴 𝜎𝐵 ൗ2 𝜎=
1 + 1.3𝛿 2
𝜀ൗ = 1.18 1 + 1.3𝛿 2 𝑇
𝑘 𝑏
169
3.1 Difusividad para sistemas de gases ideales binarios (bajas a moderadas presiones)
• Uno de los métodos más ampliamente utilizado para determinar el DAB de un sistema
binario de gases es el propuesto por Wilke y Lee (1955):
0.98
3.03 − 0.5 10−3 𝑇1.5
𝑀𝐴𝐵 Wilke y Lee (1955)
D𝐴𝐵 =
𝑃 𝑀𝐴𝐵0.5 𝜎𝐴𝐵2 Ω𝐷
2 𝑀𝐴 𝑀𝐵
• Donde 𝑀𝐴𝐵 = 𝑀𝐴 +𝑀𝐵
• P es la presión en bars
• Esta correlación retiene la misma forma general de la ecuación de Hirschfelder, con una constante
empírica basada en datos experimentales.
170
43
21/03/2021

Calcular la difusividad del vapor de disulfuro de carbono (CS2) en aire a 273 K y 1 bar usando la ecuación de
Wilke-Lee. Comparar el valor estimado con el valor experimental DAB = 8.94 x 10-6 m2/s.
Solución:
Los valores de los parámetros de Lennard-Jones son:
𝝈 𝒆𝒏 Å 𝜺/𝒌 𝒆𝒏 𝑲 M, g/gmol
CS2 4.483 467 76
Aire 3.620 97 29
𝜎𝐴 + 𝜎𝐵 4.483 + 3.62
𝜎𝐴𝐵 = = = 𝟒. 𝟎𝟓𝟐 Å Diámetro de colisión, parámetro de Lennard-
2 2 Jonnes
𝜀𝐴𝐵 𝜀𝐴 𝜀𝐵 Energía de interacción molecular para el sistema
= = 467 x 97 = 𝟐𝟏𝟐. 𝟖 𝐊 binario A en B (parámetro de Lennard-Jones en
𝑘 𝑘 𝑘
ergs /k) en Kelvin
𝑘𝑇 273
= = 𝟏. 𝟐𝟖𝟑
𝜀𝐴𝐵 212.8
184

Calcular la difusividad del vapor de disulfuro de carbono (CS2) en aire a 273 K y 1 bar usando la ecuación de
Wilke-Lee. Comparar el valor estimado con el valor experimental DAB = 8.94 x 10-6 m2/s
Solución:
1.06036 0.193 1.03587 1.76474

Ω𝐷 = + + + = 𝟏. 𝟐𝟖𝟐
1.2830.1561 exp(0.47635 𝑥 1.283) exp(1.52996 𝑥 1.283) exp(3.899411 𝑥 1.283)
−1
1 1
M𝐴𝐵 =2 + = 𝟒𝟏. 𝟗𝟖𝟏
76 29
0.98
10−3 3.03 − 273 1.5
41.981 𝐜𝐦𝟐 𝐦𝟐
D𝐴𝐵 = = 𝟎. 𝟎𝟗𝟓𝟐 = 𝟗. 𝟓𝟐𝐱𝟏𝟎−𝟔
1 4.052 2 1.282 41.981 𝐬 𝐬
8.94x10−6 − 9.52x10−6
%Error = x 100% = 𝟔. 𝟓%
8.94x10−6
185
44
21/03/2021
186
187
45
21/03/2021

Estimar la difusividad del C3H5Cl en aire a 298 K y 1 bar usando la ecuación de Wilke –
Lee. Comparar el valor estimado con el valor experimental DAB = 0.098 cm2/s.
No hay parámetros de Lennard-Jones. Por lo tanto, deben ser estimados a partir de las ecuaciones:
𝑉𝑏 = 3 14.8 + 5 3.7 + 24.6 = 𝟖𝟕. 𝟓 𝐜𝐦𝟑Τ𝐦𝐨𝐥 𝑜 𝑉𝑏 = 0.285 𝑉𝑐 1.048 = 0.285 234 1.048 = 𝟖𝟔. 𝟕 𝐜𝐦𝟑Τ𝐦𝐨𝐥
1ൗ 1ൗ
3
𝜎𝐴 = 1.18 𝑉𝑏 = 1.18 87.5 3 = 𝟓. 𝟐𝟒 Å
T (°C) 10
La temperatura normal de ebullición es,
T (k) 283.15
Tb = 318.3 K
𝜀𝐴
= 1.15 𝑇𝑏 = 1.15 318.3 = 𝟑𝟔𝟔 𝐊
𝑘
y MA = 76.5,
188
2C 6H 1O
V C2H5OH 14.8 3.7 7.4 59.2
La viscosidadProblema
de la solución con 0.05 mol de
8: Difusividad dealcohol/litro
gases idealesde agua (dilución infinita) es 1.45 cP.
Lee. Comparar el valor estimado con2.26el valor experimental D = 0.098 cm2/s.
AB
Con los valores del aire :

𝜀𝐵 7.96E-06 2
𝑐 /𝑠
𝜎𝐵 = 𝟑. 𝟔𝟐 Å = 𝟗𝟕 𝐊 y MB = 𝟐𝟗
𝑘
𝜎𝐴 + 𝜎𝐵 5.24 + 3.62
= 𝟒. 𝟒𝟑 Å Valor teórico 8.30E-06 2
𝜎𝐴𝐵 = = 𝑐 /𝑠
2 2 Error 4%
𝜀𝐴𝐵
= 366 𝑥 97 = 𝟏𝟖𝟖. 𝟒𝟐 𝐊
𝑘
𝑘𝑇 298
= = 𝟏. 𝟓𝟖𝟐
𝜀𝐴𝐵 188.42
189
46
21/03/2021

Lee. Comparar el valor estimado con el valor experimental DAB = 0.098 cm2/s.
1.06036 0.193 1.03587 1.76474

Ω𝐷 = + + + = 𝟏. 𝟏𝟕𝟓
1.5820.1561 exp(0.47635 x 1.582) exp(1.52996 x 1.582) exp(3.899411 x 1.582)
−1
1 1
M𝐴𝐵 = 2 + = 𝟒𝟐. 𝟎𝟓𝟕
76.5 29
0.98
10−3 3.03 − 298 1.5
42.06 𝐜𝐦𝟐 𝐦𝟐
D𝐴𝐵 = = 𝟎. 𝟎𝟗𝟗𝟐 = 𝟗. 𝟗𝟐𝐱𝟏𝟎−𝟔
1 4.43 2 1.175 42.06 𝐬 𝐬
0.0992 − 0.098
%Error = x 100% = 𝟏. 𝟐%
0.098
190
191
47
21/03/2021
3.2 DIFUSIVIDAD EFECTIVA DE GASES EN

MEZCLA MULTICOMPONENTE
192
3.2 Difusividad efectiva de gases en mezcla multicomponente
• Las ecuaciones de Maxwell-Stefan (MS) para difusión en sistemas

multicomponente hacen uso de la difusividad efectiva o aproximación
pseudo binaria.
• La difusividad efectiva es definida bajo la suposición de que la difusión del

componente 𝒊 solo depende de su propio gradiente de concentración,
donde Di,ef es una característica de la difusividad de la especie 𝑖 en la
mezcla la cual proviene de la difusividad binaria de MS,
J i = -c Di ,ef xi (56)
193
48
21/03/2021

En sistemas gaseosos, la difusividad binaria de MS es independiente de la
composición, así la difusividad efectiva del componente 𝒊 en una mezcla
gaseosa multicomponente puede ser derivada de la ecuación de MS según,
n
N i − yi  N j
j =1
Di ,ef = (57)
( y j Ni − yi N j )
n
1

j =1 D i j
j i
194
Para una reacción química general, podemos escribir,  A = 0

i
i i (58)
Donde, i , es el coeficiente estequiométrico y

Ai , es la especie presente en la fórmula química
La convención de signos para i es,

positiva (+) para productos y
negativa (-) para reactivos.
195
49
21/03/2021
3.3 Difusividad efectiva en mezcla multicomponente

• Para la difusión con reacción química heterogenea, la relación entre flux está
especificada por la estequiometria de la reacción descrita por la siguiente ecuación:
N1 N2 N3 Ni
= = = = (59)
1 2 3 i
• Por ejemplo, si el amoniaco va ser craqueado en un catalizador sólido de acuerdo a la
siguiente reacción,
1 3
NH3 → N 2 + H 2
2 2
bajo circunstancias donde el NH3 (1) se difunde hacia la superficie del catalizador, y N2 (2)
y H2 (3) se difunden en sentido contrario, la estequiometria de la reacción fija las
siguientes relaciones,
1 3
N 2 = − N1 y N3 = − N1
2 2
196
• Si todos los componentes tienen N = 0 excepto N1, la ecuación se escribirá

para esta reacción (mezcla estancada),
1 − y1 1
D1,ef = n
= n ' (60)
yi

yi
i = 2 D 1i
i = 2 D 1i
• Donde yi’ es la fracción mol del componente i libre del componente 1,

que se escribe,
y
yi' = i i = 2, 3, , n (61)
1 − y1
197
50
21/03/2021
Problema 9: Difusividad efectiva en una mezcla gaseosa multicomponente

• Amoniaco va ser craqueado en un catalizador sólido de acuerdo a la siguiente reacción,
1 3
NH3 → N2 + H2
2 2
• En una parte del equipo donde la presión es 1 bar y la temperatura 300 K, el análisis del gas es 40% NH3
(1), 20% N2 (2) y 40% H2 (3) en volumen. Estimar la difusividad efectiva de amoniaco en la mezcla
gaseosa.
Mezcla de gases
Catalizador
198
Primero, calculamos la difusividad binaria D12 del amoniaco en nitrógeno y la difusividad binaria
D13 del amoniaco en hidrógeno a partir de:
Gas , e A /, en K M, g/mol

NH3 2.9 467 17
N2 3.798 71.4 28
H2 2.827 59.7 2
Así, las difusividades son D12 = 0.237 cm2 Τs y D13 = 0.728 cm2 Τs
1 3
N2 = − N1 y N3 = − N1
2 2
1 3 4
N1 + N2 + N3 = N1 − N1 − N1 = N1 − N1 = N1 − 2N1 = −N1
2 2 2
199
51
21/03/2021

Solución: n
N i − yi  N j
j =1
Di ,ef = (57)
( y j Ni − yi N j )
n
1

j =1 D i j
j i
Sustituyendo en la ecuación para mezclas de gases multicomponentes,
N1 − y1 N1 + N2 + N3 N1 − y1 −N1
𝐷1,𝑒𝑓 = =
y2 N1 − y1 N2 y3 N1 − y1 N3 1 3
+ y2 N1 − y1 − N1 y3 N1 − y1 − N1
D12 D13 2 2
+
D12 D13
N1 + y1 N1 N1 1 + y1
𝐷1,𝑒𝑓 = =
1 3 1 3
y2 N1 + y1 N1 y3 N1 + y1 N1 y2 N1 + y1 N1 y3 N1 + y1 N1
2 + 2 2 + 2
D12 D13 D12 D13
205

Solución:
Cancelando N1 arriba y abajo de la ecuación,
1 + y1
𝐷1,𝑒𝑓 =
1 3
y2 + y1 y3 + y1
2 + 2
D12 D13
y sustituyendo datos,
1 + 0.4 𝐜𝐦𝟐
𝐷1,𝑒𝑓 = = 𝟎. 𝟒𝟓𝟕
0.2 + 0.5x0.4 0.4 + 1.5x0.4 𝐬
+
0.237 0.728
210
52
21/03/2021
211
Problema 10: Difusividad efectiva en una mezcla gaseosa multicomponente estancada

Amoniaco va ser absorbido de una mezcla estancada de nitrógeno e hidrógeno por contacto con ácido
sulfúrico 2N. En una parte del equipo donde la presión es 1 bar y la temperatura 300 K, el análisis del gas es
40% NH3 (1), 20% N2 (2) y 40% H2 (3) en volumen. Estimar la difusividad efectiva de amoniaco en la mezcla
gaseosa.
Solución:
La difusividad binaria es la misma que en el ejemplo anterior. Sin embargo, en este caso solo el amoniaco
se difunde a través de la interface líquida siendo qué el nitrógeno gas y el hidrógeno gas son solubles en el
ácido sulfúrico. Siendo así, se aplica la ecuación,
yi y2 0.2 y3 0.4
yi′ = siendo y2′ = = = 0.333 y y3′ = = = 0.667
1 − y1 1 − y1 1 − 0.4 1 − y1 1 − 0.4
Sustituyendo,
1 1 𝐜𝐦𝟐
𝐷1,𝑒𝑓 = = = 𝟎. 𝟒𝟑𝟏
y2′ y3′ 0.333 0.667 𝐬
D12 + D13 0.237 + 0.728
225
53
21/03/2021
3.3 DIFUSIVIDAD EN LÍQUIDOS DILUIDOS
226
3.3 Difusividad en líquidos diluidos
En contraste con los gases:

• La teoría de líquidos y las características de su transporte son aún inadecuados para
permitir un tratamiento riguroso.
• Los valores de la difusividad para líquidos son más pequeños en varios ordenes de
magnitud.
• La difusividad depende de la concentración, la cual causa cambios de viscosidad.
• A medida que la fracción mol de cada componente se acerca a la unidad, el factor

termodinámico  se aproxima también a 1.0, y la difusividad de Fick y la difusividad de
MS son iguales.
• La difusividad obtenida bajo estás condiciones es la difusividad a dilución infinita 𝑫°𝑨𝑩 .
227
54
21/03/2021
• Algunas especies se difunden como moléculas, mientras que otras como electrolitos que se
ionizan en solución.
• La neutralidad eléctrica de la solución indica que los iones se difunden en la misma razón.
• Existen dos teorías, el agujero de Eyring y la teoría hidrodinámica, explican la difusión de
solutos no electrolitos en soluciones de baja concentración.
• En el concepto de Eyring el líquido ideal es tratado

como un espacio cuasi cristalino con agujeros, el
fenómeno de transporte es descrito como saltos
de las moléculas de soluto en los agujeros del
modelo.
• La teoría hidrodinámica establece que la
difusividad líquida se relaciona a la movilidad
molecular del soluto.
228

Las leyes de la hidrodinámica ofrecen relaciones entre la fuerza y la velocidad, de la cual se obtuvo
una ecuación que ha sido desarrollada por Stokes-Einstein:
°
𝑘𝑇
𝐷𝐴𝐵 =
6 𝜋 𝑟𝐴 𝜇𝐵
°
𝐷𝐴𝐵 , es la difusividad de A diluido en B;
𝑟A, es el radio de la partícula de soluto;
𝜇𝐵 , es la viscosidad del disolvente B
Esta ecuación describe con éxito la difusión de partículas coloidales, partículas esféricas grandes o
moléculas grandes a través del disolvente que actúa como medio continuo para la especie que
difunde.
Esta ecuación considera que el líquido se encuentra adherido sobre la superficie de la partícula. Si
el fluido no tendiera a adherirse, la constante de la ecuación de Stokes-Einstein cambiaría su valor
de 6 a 4.
Los valores obtenidos con éste modelo, presentan errores cercanos al 40%, por lo que es
recomendable utilizar alguna de las siguientes correlaciones.
229
55
21/03/2021

Wilke y Chang (1955) propusieron la siguiente correlación utilizada para no electrolitos en una disolución
infinitamente diluida:
° 1ൗ
D𝐴𝐵 𝜇𝐵 7.4x10−8 Φ𝐵 𝑀𝐵 2 °
𝐷𝐴𝐵 , la difusividad a dilución infinita en cm2/s
= 0.6
𝑇 𝑉𝑏𝐴 𝑐𝑚3
𝑉𝑏𝐴 , el volumen molar al punto normal de ebullición 𝑔 𝑚𝑜𝑙
Φ𝐵 , el parámetro de asociación para el disolvente B

Disolventes comunes Φ𝐵 𝑔
MB, la masa molecular del disolvente B en
Agua 2.26 𝑔 𝑚𝑜𝑙
Metanol 1.9 T, la temperatura absoluta en Kelvin

Etanol 1.5 B, la viscosidad dinámica del disolvente B en centipoises
Benceno, éter, heptano y otros 1.0
disolventes desasociados
Cuando los valores de volumen molar no están disponibles
Tyn y Calus recomiendan:
𝑉𝑏𝐴 = 0.285𝑉𝑐1.048
230
231
56
21/03/2021
Problema 11: Difusividad de líquidos en solución acuosa

Estimar la difusividad del etanol (C2H6O) a dilución infinita en agua a 288 K, utilizando la correlación de
Wilke - Chang. Comparar el valor obtenido con el valor experimental 1x10-5 cm2/s.
Solución:
° 1ൗ
D𝐴𝐵 𝜇𝐵 7.4x10−8 Φ𝐵 𝑀𝐵 2
= 0.6
𝑇 𝑉𝑏𝐴
Con el VC = 167.1 cm3/mol, calculamos VbA,
𝑉𝑏𝐴 = 0.285 𝑉𝑐 1.048 = 0.285 167.1 1.048 = 𝟔𝟎. 𝟗 𝐜𝐦𝟑Τ𝐦𝐨𝐥
La viscosidad del agua a 288 K es 1.153 cP y el Φ𝐵 = 2.26, sustituyendo,
1ൗ 1ൗ
°
7.4x10−8 𝑇 Φ𝐵 𝑀𝐵 2 7.4x10−8 288 2.26 x 18 2 𝐜𝐦𝟐
D𝐴𝐵 = 0.6 = 0.6
= 𝟏. 𝟎𝟎𝟐 𝐱 𝟏𝟎−𝟓
𝜇𝐵 𝑉𝑏𝐴 1.153 60.9 𝐬
DAB𝑒𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙 − DAB𝑡𝑒ó𝑟𝑖𝑐𝑎 1 − 1.002

%Error = x 100% = x 100% = 𝟎. 𝟐%
DAB𝑒𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙 1
232
233
57
21/03/2021
Scheibel propuso una relación para que el parámetro de asociación, Φ𝐵 , fuera

eliminado:
𝐷°𝐴𝐵 𝜇𝐵 𝑘
= 1/3
𝑇 𝑉 𝑏𝐴
2
3𝑉𝑏𝐵 3
donde 𝑘 = 8.2𝑥10−8 1+
𝑉𝑏𝐴
Excepto para:
- Benceno como disolvente, si 𝑉𝑏𝐴 < 2𝑉𝑏𝐵 , 𝑢𝑠𝑎𝑟 𝑘 = 18.9𝑥10−8
- Otros disolventes orgánicos, si 𝑉𝑏𝐴 < 2.5𝑉𝑏𝐵 , 𝑢𝑠𝑎𝑟 𝑘 = 17.5𝑥10−8
Reid, Prausnitz y Sherwood recomendaron esta ecuación para solutos difundiéndose en disolventes
orgánicos, aunque puede tener de errores de alrededor del 20%.
234

Hayduk y Laudie propusieron una ecuación más simple para evaluar la difusividad a dilución infinita de no
electrolitos en agua,
° cm2
𝐷𝐴𝐵 = 13.26x10−5𝜇𝐵−1.14 𝑉𝑏𝐴
−0.589
en [ ]
s
Por otra parte, Hayduk y Minhas (1982) presentaron dos correlaciones para no electrolitos,
(a) Para solutos en solución acuosas:
°
𝐷𝐴𝐵 = 1.25x10−8 𝜇𝐵𝜀 (𝑉𝑏𝐴
−0.19
−0.292)𝑇1.52
9.58
donde 𝜀= − 1.12
𝑉𝑏𝐴
(b) Para soluciones no acuosas:
0.27
°
𝑉𝑏𝐵 𝑇 1.29 𝜎𝐵1.25
𝐷𝐴𝐵 = 1.55x10−8
𝑉𝑏𝐴 𝜇𝐵0.92 𝜎𝐴0.105
0.42
Donde…
235
58
21/03/2021
Donde…
 , es la tensión superficial a la temperatura del punto normal de ebullición, en dinas/cm
La letra A es para el soluto yla letra B para el disolvente
La tensión superficial se calcula por el método de Brock y Bird (1955) limitado para líquidos no polares:
  PC  
 Tbr ln  1.013  
 = PC TC ( 0.132 C − 0.278 )(1 − Tbr ) ;  C =0.9076 1 +
23 13 11 9    ; T = T T (54
1 − Tbr
br b C
 
 
y   PC  
 Tbr ln  1.013  
 = PC2 3TC1 3 ( 0.132 C − 0.278 )(1 − Tbr )
11 9
;  C =0.9076 1 +    ; T = T T (54)
1 − Tbr
br b C
 
 
Donde PC y TC son la presión y la temperatura crítica en bar y en K, respectivamente

𝑇𝑏 es la temperatura normal en el punto de ebullición
236
• La ecuación anterior para calcular la tensión superficial para líquidos no polares presenta
las siguientes restricciones:
1. El método no debe ser usado para la difusión en disolventes viscosos. Se consideran disolventes
viscosos cuando presentan valores de B por arriba de 20 cP.
2. Si el soluto es agua, el VbA debe ser considerado como un dímero

(VbA = 37.4 cm3/mol).
3. Si el soluto es un ácido orgánico y el disolvente otro que no sea agua, metanol, o butanol, el ácido debe
ser considerado un dímero con el doble valor esperado de VbA.
4. Para solutos no polares difundiéndose en alcoholes monohidróxidos (p.e. metanol, etanol, etc.) , el
valor de VbB debe ser multiplicado por 8B .
237
59
21/03/2021
Problema 11: Difusividad de líquidos en solución acuosa

Repetir el problema anterior de estimar la difusividad del etanol (C2H6O) a dilución infinita en agua a 288 K, pero
utilizando la correlación de Hayduk y Minas. Comparar el valor obtenido con el valor experimental 1x10-5 cm2/s.
Solución:
°
𝐷𝐴𝐵 = 1.25x10−8 𝜇𝐵𝜀 (𝑉𝑏𝐴
−0.19
−0.292)𝑇1.52
Con 𝑉𝑏𝐴 = 𝟔𝟎. 𝟗 𝐜𝐦𝟑Τ𝐦𝐨𝐥
y la viscosidad de 1.153 cP, sustituyendo,
9.58 9.58
𝜀= − 1.12 = − 1.12 = −0.963
𝑉𝑏𝐴 60.9
°
𝐷𝐴𝐵 = 1.25x10−8 1.153 −0.963 60.9−0.19 − 0.292 288 1.52 = 𝟎. 𝟗𝟗𝟎𝟔𝐱𝟏𝟎−𝟓 𝐜𝐦𝟐 Τ𝐬
DAB𝑒𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙 − DAB𝑡𝑒ó𝑟𝑖𝑐𝑎 1 − 0.9906

%Error = x 100% = x 100% = 𝟎. 𝟗𝟒%
DAB𝑒𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙 1
238
Problema 12: Difusividad de líquidos en solución no acuosa
Estimar la difusividad del ácido acético (C2H4O2) en Parámetro Ácido acético Acetona
solución diluida de acetona (C3H6O) a 313 K usando la Tb, K 390.4 329.2
correlación de Wilke - Lee. Comparar el valor estimado TC, K 594.8 508

PC, bars 57.9 47.0
con el valor experimental de 4.04x10-5 cm2/s .
VC, bars 171 209
, cP --- 0.264
Solución: M 60 58
Con el VC = 171 cm3/mol, calculamos 𝑉𝑏𝐴 = 0.285 𝑉𝑐 1.048 = 0.285 171 1.048 = 𝟔𝟐. 𝟑𝟖 𝐜𝐦𝟑Τ𝐦𝐨𝐥
1ൗ 1ൗ
°
7.4x10−8 𝑇 Φ𝐵 𝑀𝐵 2 7.4x10−8 313 1 x 58 2 𝐜𝐦𝟐
El Φ𝐵 = 1.0, sustituyendo, D𝐴𝐵 = 0.6 = 0.6
= 𝟓. 𝟓𝟗𝟓 𝐱 𝟏𝟎−𝟓
𝜇𝐵 𝑉𝑏𝐴 0.264 62.38 𝐬

%Error = x 100% = x 100% = −𝟑𝟖. 𝟓 %
239
60
21/03/2021
Del problema anterior, estimar la difusividad del ácido Parámetro Ácido acético Acetona
acético (C2H4O2) en solución diluida de acetona (C3H6O) a Tb, K 390.4 329.2
313 K usando la correlación de Hayduk y Minas. TC, K 594.8 508
Comparar el valor estimado con el valor experimental de PC, bars 57.9 47.0
VC, bars 171 209
4.04x10-5 cm2/s .
, cP --- 0.264
Solución: M 60 58
De acuerdo con la restricción 3 mencionada antes, el volumen molar del ácido acético a utilizar debe ser el doble,
𝑉𝑏𝐴 = 2x62.38 = 𝟏𝟐𝟒. 𝟕𝟔 𝐜𝐦𝟑Τ𝐦𝐨𝐥, el volumen molar de la acetona se calcula como,
𝑉𝑏𝐵 = 0.285 𝑉𝑐 1.048 = 0.285 209 1.048 = 𝟕𝟔. 𝟗𝟖 𝐜𝐦𝟑Τ𝐦𝐨𝐥
Se calcula la tensión superficial para la acetona (B): 𝑇𝑏 329.2

𝑇𝑏𝑟 = = = 𝟎. 𝟔𝟒𝟖
𝑇𝐶 508
𝛼𝐶 = 𝟕. 𝟑𝟏𝟗
𝒅𝒊𝒏𝒂𝒔
𝜎𝐶 = 𝟐𝟎
𝒄𝒎
240
Del problema anterior, estimar la difusividad del ácido Parámetro Ácido acético Acetona
acético (C2H4O2) en solución diluida de acetona (C3H6O) a Tb, K 390.4 329.2
313 K usando la correlación de Hayduk y Minas. TC, K 594.8 508
Comparar el valor estimado con el valor experimental de PC, bars 57.9 47.0
VC, bars 171 209
4.04x10-5 cm2/s .
, cP --- 0.264
Solución: M 60 58
Se calcula la tensión superficial para la ácido acético (A): 𝑇𝑏 390.4
𝑇𝑏𝑟 = = = 𝟎. 𝟔𝟓𝟔
𝑇𝐶 594.8
𝛼𝐶 = 𝟕. 𝟗𝟏𝟎
𝒅𝒊𝒏𝒂𝒔
𝜎𝐶 = 𝟐𝟔. 𝟐
𝒄𝒎
°
0.27
𝑉𝑏𝐵 𝑇1.29 𝜎𝐵1.25 77 0.27 313 1.29 20 1.25 −𝟓
𝐜𝐦𝟐
𝐷𝐴𝐵 = 1.55x10−8 = 1.55x10−8 = 𝟑. 𝟖𝟒 𝐱 𝟏𝟎
0.42
𝑉𝑏𝐴 𝜇𝐵0.92 𝜎𝐴0.105 124.76 0.42 0.264 0.92 26.2 0.105 𝐬

%Error = x 100% = x 100% = 𝟒. 𝟗𝟓 %
241
61
21/03/2021
3.3 Difusividad en líquidos diluidos (altas temperaturas)
Tyne recomienda la siguiente ecuación para extrapolar a temperaturas más altas:
° 𝑛
𝐷𝐴𝐵 𝑇1 𝑇𝑐 − 𝑇2
° =
𝐷𝐴𝐵 𝑇2 𝑇𝑐 − 𝑇1
Donde 𝑇𝑐 es la temperatura crítica del disolvente B en escala absoluta Kelvin y 𝒏 es un

exponente relativo al calor latente de vaporización del disolvente 𝜟𝑯𝝊 :
Δ𝐻𝜐 (kJ/kmol) 7900-30000 30000-39000 39000-46000 46000-50000 >50000
𝑛 3 4 6 8 10
242
3.3 Difusividad en líquidos (electrolitos) diluidos
La difusividad de una sal univalente en una solución diluida es dada por la ecuación
de Nernst:
°
2𝑅𝑇
𝐷𝐴𝐵 =
1 1
+ ℱ
𝜆0+ 𝜆0−
𝜆0+ 𝑦 𝜆0− son las conductividades iónicas molares (concentración cero) en  cm2/mol
ℱ Es la constante de Faraday, 96487 Coulomb/mol
R, es la constante universal de los gases, 8.314 J / mol K
T, es la temperatura en escala absoluta Kelvin
243
62
21/03/2021
3.3 Difusividad en líquidos (electrolitos) diluidos
La teoría de difusión de sales en un medio acuoso a bajas concentraciones fue

desarrollada aún más para disoluciones infinitas, por la ecuación de Nernst – Haskell:
1 1
𝑅𝑇 +
° 𝑛+ 𝑛−
𝐷𝐴𝐵 =
1 1
+ ℱ2
𝜆0+ 𝜆0−
1 1
La ecuación ha sido extendida para iones polivalentes, donde + 𝑛− ,
𝑛+
es un parámetro para las respectivas cargas del catión y del anión de la sal
244
Valores de
conductividad molar
iónica en agua a
25°C
245
63
21/03/2021
VALORES
EXPERIMENTALES DE
DIFUSIVIDAD DE
ELECTROLITOS EN
DISOLUCIÓN ACUOSA
246
En mezclas líquidas multicomponente, es difícil obtener valores numéricos de difusividad

que relacionen el flux con el gradiente de concentración.
Cuando un soluto diluido se difunde a través de una disolución homogénea de una mezcla
de disolventes, se simplifican las correlaciones empíricas al modelo sugerido por Perkins y
Geankoplis:
n
D1,0ef M0.8 =  x j D 10j  0.8
j (62)
j =2
La ecuación fue probada para ocho sistemas ternarios presentando errores menores al
20%, excepto para el caso del CO2 como soluto.
247
64
21/03/2021
Para el CO2 como soluto diluido difundiéndose en una mezcla de disolventes, Takahashi y
col. (1982) recomiendan utilizar,
13
 j 
13
 M  n
0
DCO2 , ef   =  0
x j D CO2 j
  (63)
 VM  j =1
j  CO2
 Vj 
Donde VM es el volumen molar de la mezcla, y Vj es el volumen molar del componente j,

ambos a la temperatura T, en cm3/mol.
Las pruebas con sistemas ternarios incluyendo el CO 2 como soluto, muestran desviaciones
menores a 4%.
248
Problema 15: Difusividad efectiva de un soluto diluido en una mezcla homogénea

de disolventes
Estimar el coeficiente de difusión del ácido acético (1) a muy baja concentración difundiéndose en una
mezcla de disolventes que contiene 40% en masa de alcohol etílico (3) en agua (2) a temperatura de
298 K usando la ecuación de Hayduk y Minhas (1982). Sustituyendo en la ecuación se obtiene 1.32x10-5
cm2/s. Utilizar los datos de la siguiente tabla:
Parámetro Ácido acético Etanol
Tb, K 390.4 351.4
Tc, K 594.8 513.9
Pc, bar 57.9 61.4
Vc, cm3/mol 171.0 167.0
 a 298 K, cP ---- 1.043
249
65
21/03/2021
3.4 DIFUSIVIDAD EN LÍQUIDOS CONCENTRADOS
250
3.4 Difusividad en líquidos concentrados

El cálculo de la difusividad de MS en líquidos concentrados es obtenida por medio de la fórmula de Vignes
recomendada por Reid y col. (1987):
𝑫12 = 𝑫°12 𝑥2 𝑫°21 𝑥1
Esta ecuación tiene menos efectividad en las mezclas que contienen compuestos de asociación, como el
alcohol.
Ahora se obtiene la difusividad de Fick como,
D12 = 𝑫12 Γ
Una modificación hecha por Leffler y Cullian puede ser más certera en soluciones más concentradas:
𝑫12 𝜇 = 𝑫°12 𝜇2 𝑥2 𝑫°21 𝜇1 𝑥1
251
66
21/03/2021
Problema 14: Difusividad de líquidos en solución concentrada

Estimar el coeficiente de difusión de MS y de Fick para la mezcla acetona (1) – benceno (2) de composición
x1=0.7808 a 298 K. Los coeficientes de difusión a dilución infinita son:
0
m2 0
m2
𝐷12 = 2.75 x 10−9 𝑦 𝐷21 = 4.15 x 10−9
s s
A partir de la ecuación NRTL, para este sistema, el factor de corrección termodinámico es igual a =0.871:
El valor experimental de D12 a dicha concentración es 3.35 x 10-9 m2/s.
Solución:
m2
𝑫12 = 𝑫°12 𝑥2 𝑫°21 𝑥1 = 2.75 x 10−9 0.2192 4.15 x 10−9 0.7808 = 3.792 x 10−9
s
𝐦𝟐
D12 = 𝑫12 Γ = 3.792 x 10−9 x 0.871 = 𝟑. 𝟑 𝐱 𝟏𝟎−𝟗
𝐬

%Error = x 100% = x 100% = 𝟏. 𝟓 %
252
253
67
21/03/2021
4. Difusión molecular en estado estacionario
254
4.1 FLUX MOLAR Y ECUACIÓN DE CONTINUIDAD
255
68
21/03/2021
4.0 Transporte de masa por difusión y convección
• El transporte de masa por difusión es conocido como difusión
molecular y es el resultado del movimiento individual de las
moléculas a través de una sustancia
256

• La teoría cinética explica el movimiento molecular como una consecuencia
de la energía térmica del sistema:
– La teoría cinética se simplifica suponiendo que la velocidad de una molécula es
uniforme, hasta que colisiona con otra.
– Al chocar, la velocidad, la magnitud y dirección cambian.
– La distancia promedio entre choques se llama camino libre medio.
– Las moléculas viajan en zig-zag por las colisiones.
257
69
21/03/2021
• La distancia neta
recorrida por las
moléculas en un
determinado
tiempo es la razón
de la difusión.
258

• Mientras la difusión molecular es un proceso lento, la difusión convectiva se
alcanza gracias a la agitación mecánica o por el movimiento turbulento de un fluido.
• Los mecanismos de transferencia de masa dependen de la dinámica del sistema:
El movimiento aleatorio de las Un fluido en movimiento

moléculas en un fluido quieto
Transferencia de masa molecular Transferencia de masa convectiva
259
70
21/03/2021
• En una solución, donde no hay uniformidad en la concentración, el sistema

buscará por medio de la difusión el equilibrio:
– Sustancias con alta concentración se mueven a un espacio en dirección decreciente de

la fracción molar del soluto, esta razón es conocida como gradiente de concentración.
• La primera ley de Fick para mezclas binarias es definida con el gradiente de

concentración y con la difusividad de un soluto A en una solución B:
𝐽𝐴 = − 𝑐D𝐴𝐵 ∇𝑥𝐴 𝑜 𝐽𝐴 = −D𝐴𝐵 ∇𝑐𝐴
260
✓ Es una relación para el flux molar en un sistema isotérmico e isobárico

✓ Escrito en un solo eje de coordenadas:
𝑑𝑐𝐴
𝐽𝐴,𝑧 = −D𝐴𝐵
𝑑𝑧
✓ El signo negativo enfatiza que la difusión ocurre en sentido de la caída de

concentración.
261
71
21/03/2021
4.1 Flux molar y la ecuación de continuidad

Para abordar el problema de transferencia de masa desde otro punto de vista, se necesita:
✓ Una ecuación diferencial (ED) para describir la

situación física y condiciones límite
✓ Para simplificar la ED de transferencia de masa,

eliminar términos dependiendo de la condición
física del problema
✓ Aplicar un balance de materia en el elemento

de volumen diferencial
✓ Obtiener la distribución de flux másico con un perfil de concentraciones
262
4.1 Flux molar y la ecuación de continuidad

La ecuación de continuidad para el componente A en una mezcla, describe el cambio de concentración de A con respecto
al tiempo, en un punto fijo en el espacio, con reacción química de A y transferencia de masa de A
• En términos vectoriales, la ecuación resultante se escribe como (Bird y col., 1960),
cA
  NA + − RA = 0 (64)
t
Reacción química
Vector de divergencia
de flux molar Acumulación
N A, x N A, y N A, z
– El primer término en coordenadas rectangulares, se escribe como,   N A = + + = 0 (65)
x y z
• Considerando la difusión sin acumulación (en estado estacionario) y sin reacción química para un componente i,
  Ni = 0 i = 1, 2, , n (66)
• En coordenadas rectangulares, con difusión sólo en la dirección z, la ecuación se simplifica a,
dNi , z
= 0 i = 1, 2, , n (67)
dz n
Ni = Ni , z para
siendo el flux molar constante = constante i i = 1, en
cada componente 2, la,mezcla,
n y  N = constante
i =1
i (68)
263
72
21/03/2021
4.2 DIFUSIÓN MOLECULAR DE GASES EN

ESTADO ESTACIONARIO
264
4.2 Difusión molecular de gases en estado estacionario

En un sistema binario (de la especie A en B), el flux molar en la dirección z, relativo a la velocidad molar
promedio de A, se describe según la primera ley de Fick,
𝑑𝑦𝐴
𝐽𝐴,𝑧 = −𝑐𝐷𝐴𝐵
𝑑𝑧
El flux molar de la especie A, que es consecuencia del movimiento global del fluido, se escribe como,
𝑁𝑦𝐴 = 𝑦𝐴 𝑁𝐴 + 𝑁𝐵
Siendo el flux molar 𝑁𝐴 relativo a coordenadas estacionarias,
𝑑𝑦𝐴
𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 𝑁𝐴 + 𝑁𝐵 a partir de N i = J i + ci V = J i + yi N (20)
𝑑𝑧
O de manera general para mezclas multicomponentes:
𝑛
𝑁𝐴 = −𝑐𝐷𝐴𝑀 𝛻𝑦𝐴 + 𝑦𝐴 ෍ 𝑁𝑖
𝑖=1
Contribución del
Contribución del gradiente movimiento de A en
de concentración el seno de la solución
265
73
21/03/2021
Retomando lo anterior, para el caso de una dimensión, la ecuación N i = J i + ci V = J i + yi N (20)se

escribe:
dy dy
N A = −cDAB A + yA N = −cDAB A + yA ( N A + N B ) = 0 (69)
dz dz
𝑑𝑦𝐴
Separando las variables, 𝑁𝐴 − 𝑦𝐴 𝑁𝐴 + 𝑁𝐵 = −𝑐𝐷𝐴𝐵 𝑑𝑧
𝑁𝐴 −𝑦𝐴 𝑁𝐴 + 𝑁𝐵 𝑐𝐷𝐴𝐵
=
−𝑑𝑦𝐴 𝑑𝑧
Para obtener,
−dyA dz
= (70)
N A − yA ( N A + N B ) cDAB
Para una mezcla gaseosa, si la presión y la temperatura son constantes, c y DAB lo son también,
independientemente de la posición y de la composición
266
La ecuación puede ser integrada entre dos condiciones de frontera,
a z = z1 yA = yA1
a z = z2 yA = y A 2
Donde 1 indica el inicio de la trayectoria de difusión (yA alta) y 2 el final de la trayectoria de

difusión (yA baja). Definiendo que z2 – z1 =  , al integrar se obtiene,
1  N A − yA 2 ( N A + N B )  δ
ln  = (71)
( N A + N B )  N A − yA1 ( N A + N B )  cDAB
267
74
21/03/2021
1  N A − yA 2 ( N A + N B )  δ
De la ecuación anterior, ln  = (71)
( N A + N B )  N A − yA1 ( N A + N B )  cDAB
𝑵𝐀 𝑁A − 𝑦A2 𝑁A + 𝑁B δ 𝑵𝐀
Se multiplica de ambos lados por NA, ln =
𝑁A + 𝑁B 𝑁A − 𝑦A1 𝑁A + 𝑁B 𝑐 DAB
Se divide entre 𝑁A + 𝑁B en el numerador y denominador del logaritmo,

𝑁A
𝑁A − 𝑦A2 δ 𝑵𝐀
𝑁𝐴 + 𝑁𝐵
ln =
𝑁A + 𝑁B 𝑁A 𝑐 DAB
𝑁𝐴 + 𝑁𝐵 − 𝑦𝐴1
Se despeja 𝑵𝐀 para finalmente obtener,
𝑁A
𝑁A 𝑐DAB 𝑁𝐴 + 𝑁𝐵 − 𝑦A2
𝑵𝐀 = ln
𝑁A + 𝑁B δ 𝑁A
𝑁𝐴 + 𝑁𝐵 − 𝑦𝐴1
268
• Retomando la ecuación anterior,

 NA 
NA cDAB  ( N + N ) − yA 2 
NA = ln  A B
 (71)
( NA + NB ) δ  NA
− yA1 
 ( NA + NB ) 
 
• Definiendo la fracción de flux molar, A: NA
A = (72)
( NA + NB )
• Sustituyendo: cDAB   A − yA 2 
NA =  A ln   (73)
z   A − yA1 
• Si la transferencia de masa en la fase gas bajo condiciones ideales, se puede sustituir

en la ecuación, c = P RT ,(7)
DAB P   A − yA 2 
NA =  A ln   (74)
RT z   A − yA1 
269
75
21/03/2021
Problema 16: Producción de carbonilo de Níquel , en estado estacionario, flux binario en una dimensión
El carbonilo de níquel (A) se produce pasando monóxido de carbono (B) a 323 K y 1 atm sobre una placa de níquel. La
siguiente reacción toma lugar sobre la superficie del sólido: Ni(s) + 4CO(g) → Ni ( CO )4 (g)
La reacción es muy rápida, por lo que la presión parcial del CO en la superficie del metal es prácticamente cero. Los gases se
difunden a través de una película de 0.625 mm de espesor. En estado estacionario, estimar la rapidez de producción de
carbonilo de níquel, en mol/m2 de superficie sólida. La composición en el seno de la fase gas es 50% en mol de CO. La
difusividad binaria del gas bajo estás condiciones es DAB= 20 mm2/s.
Película gaseosa
Placa de níquel
270
Solución:
De la ecuación,
Ni(s) + 4CO(g) → Ni ( CO )4 (g)
𝑁B 𝑁A Película gaseosa
Siendo,
𝑁B = −4𝑁A , 𝑁A + 𝑁B = −3𝑁A
Siendo,
𝑁A 𝑁A 𝟏
ΨA = = =−
𝑁A + 𝑁B −3𝑁A 𝟑
Placa de níquel
271
76
21/03/2021
Solución:
1
DAB 𝑃 ΨA − 𝑦A2 1 DAB 𝑃 − − 𝑦A2
𝑁A = ΨA
R𝑇𝑧
ln
ΨA − 𝑦A1
=− ln 3
3 R𝑇𝑧 1
− 3 − 𝑦A1
Película gaseosa
1
𝑁A =
DAB 𝑃
ln 3 + 𝑦A1
3R𝑇𝑧 1
3 + 𝑦A2
1
𝑁A =
2x10−5 1.013x105
ln 3 + 1.0 = 𝟎. 𝟏𝟖𝟗 𝐦𝐨𝐥/𝐦𝟐 𝐬
3 8.314 323 0.625𝑥10−3 1
3 + 0.5
Placa de níquel
272
4.2 Difusión molecular de gases en estado estacionario (A a través de B estancado)

• Aplicando la ley de conservación de materia en la especie A a través de B (Caso 1), ésta puede entrar o
salir del sistema en virtud de la difusión o del movimiento global del fluido, además se puede formar o
destruir a causa de reacciones químicas:
NA
NA = constante y NB = 0,  A = =1
NA + NB
DAB P   A − yA 2  DAB P 1 − yA 2 
NA =  A ln   (74) NA = ln   (75)
RT z   A − yA1  RT z  1 − yA1 
cA pA DAB P  P − pA 2 
Sustituyendo en la ecuación yA = = , NA = ln   (76)
c P RT z  P − pA1 
pB2 − pB1 = pA1 − pA2 pA1 − pA2

Si sabemos que y =1
pB2 − pB1
D P
RT pRT
D P
( ( ) )
N A =N A = AB AB pA1 −pAp1A−2 pA2 (78)(78) siendo ( (
pB , M p=B , M p=B2 −pBp2 B−
siendo ) )
pBln1  plnB2 ppBB21 pB1 
zB , M z
B,M p
1
273
77
21/03/2021
4.2 Difusión molecular de gases en estado estacionario (Contra difusión equimolar)

Para el caso de contradifusión equimolar en estado estacionario (caso 2):
𝑁𝐴 𝑁𝐴
NB = - NA = cte, σ 𝑁𝑖 =0, Ψ=𝑁 = = 𝐼𝑛𝑑𝑒𝑡𝑒𝑟 𝑖𝑛𝑎𝑑𝑜
𝐴 + (−𝑁𝐴 ) 0
𝑑𝑦𝐴 𝑑𝑦𝐴
• Partiendo de la ecuación 𝑁𝐴 = −𝑐𝐷𝐴𝐵 + 𝑦𝐴 𝑁𝐴 − 𝑁𝐴 = −𝑐𝐷𝐴𝐵
𝑑𝑧 𝑑𝑧
𝑃 𝑑𝑦𝐴
• En términos de presión, 𝑁𝐴 = −𝐷𝐴𝐵 𝑅𝑇 𝑑𝑧
𝑃 𝑃
• Integrando, 𝑁𝐴 = −𝐷𝐴𝐵 𝑦𝐴2 − 𝑦𝐴1 = 𝐷𝐴𝐵 𝑦𝐴1 − 𝑦𝐴2
𝑅𝑇𝑧 𝑅𝑇𝑧
𝐷𝐴𝐵
• En términos de presión parcial, 𝑁𝐴 = 𝑝𝐴1 − 𝑝𝐴2
𝑅𝑇𝑧
274
Problema 12
Se esta difundiendo oxígeno (A) a través de monóxido de carbono (B) en condiciones de estado estacionario, con el
monóxido de carbono sin difundirse. La presión total es 1x105 N/m2, y la temperatura es 0°C. La presión parcial del
oxígeno en placas separadas por 2 mm es, respectivamente 13000 y 6500 N/m2. Calcular la rapidez de difusión del
oxígeno en kmol/s, a través de cada metro cuadrado de las dos placas.
275
78
21/03/2021
Problema 13
Calcular la rapidez de difusión del oxígeno (A), suponiendo que el gas no se esta difundiendo es una mezcla de metano
(B) e hidrógeno (C) en la relación de volumen 2:1,
𝑚2
𝐷𝑂2 −𝐶𝐻4 = 1.86 x 10−5
𝑠
𝑚2
𝐷𝑂2 −𝐻2 = 6.99 x 10−5
𝑠
276
4.3 DIFUSIÓN MOLECULAR DE LÍQUIDOS EN

ESTADO ESTACIONARIO
277
79
21/03/2021
4.3 Difusión molecular de líquidos en estado estacionario (A a través de B estancado)

A partir de la ecuación,
DAB 𝑃 ΨA − 𝑦A2
𝑁A = ΨA ln
R𝑇𝑧 ΨA − 𝑦A1
AB D ρ ΨA − 𝑥A2
Nos queda para el caso de líquidos 𝑁𝐴 = ΨA 𝑧 M ln
ΨA − 𝑥A1
prom
Siendo NA=constante y NB=0 entonces ΨA=1

Así la ecuación anterior se escribe, °
𝐷𝐴𝐵 ρ
𝑁𝐴 = 𝑥𝐴1 − 𝑥𝐴2
𝑧𝑥𝐵,𝑀 M prom
𝑥𝐵2 − 𝑥𝐵1
Donde: 𝑥𝐵,𝑀 =
𝑥
𝑙𝑛 𝐵2
𝑥𝐵1
278
4.3 Difusión molecular de gases en estado estacionario (Difusión equimolar) ***
¨Contra difusión equimolar en estado estacionario¨

NA=-NB
°
𝐷𝐴𝐵 ρ
𝑁𝐴 = 𝑥𝐴1 − 𝑥𝐴2
𝑧 𝑀 𝑝𝑟𝑜𝑚
279
80
21/03/2021
Problema 14 (Difusión molecular de líquidos en estado estacionario)

• Un cristal de CuSO45H2O se disuelve en un gran tanque de agua pura a 273 K.
Estimar la rapidez a la cual el cristal se disuelve calculando el flux de CuSO4 de la superficie
del cristal hacia el seno de la solución. Asumir que la difusividad molecular ocurre a través
de la película líquida de espesor 0.01 mm alrededor del cristal. En la superficie del cristal,
del lado del líquido la concentración de la solución es saturada en CuSO4 mientras que
lejos de la superficie, del otro lado de la película liquida la solución es agua pura. La
solubilidad del CuSO4 en agua a 273 K es 24.3 g de cristal/100 g de agua y la densidad de la
solución saturada es 1140 kg/m3. La densidad del CuSO4 a 273 K, en solución acuosa y a
dilución infinita es 999.8 kg/m3. Corrección de  por temperatura T/334 Agua .
280
2.8 Analogías entre transferencia de momentum, calor y masa

Números adimensionales
• Al haber un fluido en movimiento existe un gradiente de velocidad y, a su vez, de
transferencia de masa. En algunas ocasiones hay gradiente de temperatura originado por
la transferencia de calor.
• En transferencia de momentum al considerarse el perfil de velocidades en un flujo se
habla de una fuerza a través de una unidad de área paralela a la superficie, también
conocida como esfuerzo cortante:
𝑑𝑢 𝜇 𝑑 𝑢𝜌 𝑑 𝑢𝜌
𝜏 = −𝜇 =− = −𝜈 Ecuación análoga a la primera ley de Fick
𝑑𝑧 𝜌 𝑑𝑧 𝑑𝑧
𝜈; viscosidad cinemática, conocida como la

difusión del momentum.
𝑢𝜌; es interpretada como la concentración
volumétrica de momentum.
𝜏, el esfuerzo cortante también se ve como
un flux viscoso de momentum.
𝜈(m2/s); tiene las mismas unidades que DAB
281
81
21/03/2021
2.8 Analogías entre transferencia de momentum, calor y masa.

Números adimensionales.
• El número de Schmidt es la razón adimensional entre las difusividades de momentum y

masa:
μ 𝑣 (𝐷𝑖𝑓𝑢𝑠𝑖𝑣𝑖𝑑𝑎𝑑 𝑝𝑜𝑟 𝑇. 𝑀𝑜 𝑒𝑛𝑡𝑢 )

𝑆𝑐 = =
ρ𝐷𝐴𝐵 𝐷𝐴𝐵 (𝐷𝑖𝑓𝑢𝑠𝑖𝑣𝑖𝑑𝑎𝑑 𝑝𝑜𝑟 𝑇. 𝑀𝑎𝑠𝑎)
- A altas concentraciones molares el paralelismo de los valores entre la difusividad

molecular y de momentum se rompe.
- Para un gas el Sc se encuentra alrededor de la unidad, mientras que la teoría cinética más
sofisticada presenta valores experimentales entre 0.67 y 0.83.
- En mezclas binarias el Sc aumenta su rango arriba de 5 unidades.
- En el caso de líquidos el Sc es mucho mayor, por ejemplo la autodifusión del agua es 297
a 25°C.
282

Números adimensionales.
Cuando hay un gradiente de temperatura entre el fluido y el plano, la razón de
transferencia de calor en régimen laminar es:
𝑑𝑇 𝑘 𝑑 𝑇𝐶𝑝 𝜌 𝑑 𝑇𝐶𝑝 𝜌
𝑞 = −𝑘 =− = −𝛼
𝑑𝑧 𝐶𝑝 𝜌 𝑑𝑧 𝑑𝑧
donde,
𝑘, es la conductividad térmica del fluido,
𝑇 𝐶𝑝 𝜌, es vista como una concentración térmica volumétrica,
𝑑𝑖𝑠𝑡𝑎𝑛𝑐𝑖𝑎 2
𝛼, es la difusividad térmica con unidades de
𝑡𝑖𝑒𝑚𝑝𝑜
Es una ecuación análoga a la primera ley de Fick y al esfuerzo cortante.
283
82
21/03/2021

Números adimensionales
• El número de Prandtl es la razón entre la difusividad térmica y de momentum:
𝜈 𝐶𝑝 𝜇
𝑃𝑟 = =
𝛼 𝑘
Valores experimentales tienen un rango entre 0.65 y 0.9 para gases de baja
presión.
Los datos cuantitativos más altos son de líquidos, por ejemplo 7.02 para agua.
• El número Lewis, es la razón entre la difusividad másica y térmica:
𝛼 𝑆𝑐
𝐿𝑒 = =
𝐷𝐴𝐵 𝑃𝑟
Es importante para resolver problemas de transferencia simultánea de masa y

calor, como son las operaciones de humidificación.
284
Difusividad
másica en
gases
285
83
21/03/2021
Difusividad
másica en
líquidos
286
Difusividad
másica en
sólidos
287
84
21/03/2021
288
289
85
21/03/2021
290
291
86
21/03/2021
292
293
87
21/03/2021
294
298
88

2 Transferencia de Masa Mecanismos de Transferencia 2020

Cargado por

Copyright:

Formatos disponibles

2 Transferencia de Masa Mecanismos de Transferencia 2020

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

2 Transferencia de Masa Mecanismos de Transferencia 2020

Cargado por

Copyright:

Formatos disponibles

21/03/2021

Curso de Transferencia de Masa

Prof. Francisco Javier Garfias Vásquez

2. Transferencia de masa molecular

2. Transferencia de masa molecular

• En una mezcla multicomponente, la concentración másica (o densidad) de la especie

donde “𝑛” corresponde al número de especies en la mezcla.

La composición másica es 6.703 kg de oxígeno gas, 21.864 kg de nitrógeno gas, 37.36

2. Transferencia de masa molecular

• La suma de las fracciones debe ser igual a la unidad,

Gas Concentración másica (kg/m3) Fracción másica

Oxígeno 6.703 0.0995

2. Transferencia de masa molecular

Donde MA es la masa molecular.

Donde pA corresponde a la presión parcial de la especie A en la mezcla, T es la

Gas Concentración másica (kg/m3) Fracción másica Concentración molar

donde P es la presión total de la mezcla gaseosa.

Para gases que obedecen la ley del gas ideal,

La suma de las fracciones molares es por definición igual a 1,

Problema 1: Concentración de la alimentación en un absorbedor de gas

Aceite no volátil Gas de salida

Producto líquido Gas alimentado

Problema 1: Concentración de la alimentación en un absorbedor de gas

a) La concentración molar total del gas en la alimentación

P 500 kPa 𝐤𝐦𝐨𝐥

Problema 1: Concentración de la alimentación en un absorbedor de gas

Gas Base: 100 kmol de Masa molecular Masa (kg)

Problema 1: Concentración de la alimentación en un absorbedor de gas

Problema 2: Concentración de Nitrato de Potasio en una solución de lavado

KNO3 + NaCl Cristales lavados

Solución de lavado Solución fresca

Problema 2: Concentración de Nitrato de Potasio en una solución de lavado

Componente Masa (kg) Masa molecular Mol (kmol)

Problema 2: Concentración de Nitrato de Potasio en una solución de lavado

Problema 2: Concentración de Nitrato de Potasio en una solución de lavado

Componente Mol (kmol) Fracción mol

KNO3 0.237 (0.237/4.455) = 0.053

2.2 Velocidad y flux

2.2 VELOCIDAD Y FLUX

2.2 Velocidad y flux

• En un sistema multicomponente las diversas especies se moverán a

diferentes velocidades con respecto a la velocidad molar o másica

• Por ésta razón, una evaluación de la velocidad característica para una

mezcla de gases requerirá del promedio de las velocidades de cada

una de las especies presentes.

2.2 Velocidad y flux

2.2 Velocidad y flux

• La velocidad másica promedio para una mezcla multicomponente es

• Donde vi es la velocidad absoluta de la especie 𝑖 relativa al eje de

2.2 Velocidad y flux

Problema 3: Velocidad másica promedio

Primero, calculamos la velocidad cuadrática media,

Problema 3: Velocidad másica promedio

v = ෍ ωi vi = 0.733x683 + 0.062x499 + 0.114x412 + 0.091x359 = 𝟔𝟏𝟏. 𝟐 𝐦/𝐬

2.2 Velocidad y flux

• v𝑖 − v, la velocidad de difusión de la especie i con respecto a la velocidad

• v𝑖 − V, la velocidad de difusión de la especie i con respecto a la velocidad

Problema 4: Velocidad másica promedio

V = ෍ yi vi = 0.88x683 + 0.04x499 + 0.05x412 + 0.03x359 = 𝟔𝟓𝟐. 𝟒 𝐦/𝐬

2.2 Velocidad y flux (densidad de flujo)