3. Se verifica una transferencia de masa de una esfera de naftaleno con radio igual a 10 mm.

La
esfera está en un gran volumen de aire en reposo a 52.6°C y 1 atm de presión. La presión de vapor
del naftaleno a 52.6°C es 1.0 mmHg. La difusividad del naftaleno en el aire a 0°C es
−6 2
5.16 ×10 m /s . Calcule la velocidad de evaporación del naftaleno de la superficie en
2
kgmol /s ∙ m .
SOLUCIÓN:

Para el calculo de la velocidad de evaporación del naftaleno de la superficie se utilizará la difusión

en estado estacionario en fluidos sin movimiento y flujo laminar utilizando la siguiente ecuación:

D AB ∙ P
N A= ( p − pA 2)
R ∙ T ∙ r ∙ pB , M A 1

p B 2− p B 1
pB , M =
ln
( )
pB 2
pB 1

Se sabe:

r =10 mm=0.01 m
T aire =52.6° C=325.75 K

paire =1 atm=101325 pa

presión de vapor del Naftaleno=1 mmHg¿ p A 1 =133.322 pa

p A 2 =0

Pa∙ m3
R=8 314
kmol ∙ K

Considerando la presión atmosférica de 1 atm (101 325 pa)

pB 1 =PT − p A 1=( 101 325 pa−133.322 pa )=101 191.678 pa

pB 2=PT − p A 2=( 101 325 pa−0 pa )=101 325 pa

Cálculo de pB , M
 p B 2− p B 1
pB , M =
ln
( )
pB 2
pB 1

Sustituyendo los valores:

(101325−101 191.678) pa
pB , M =
ln( 101 325 pa
101 191.678 pa )
pB , M =101258.3244 pa

Sustituyendo los valores conocidos:

(5.16 ×10−6 m2 /s )(101325 pa)

N A= ( 133.322−0 ) pa
( )
3
Pa ∙ m
8 314 ( 325.75 K )( 0.01 m ) (101258.3244 pa)
kmol ∙ K

−8 2
N A =2.5418 ×10 kgmol /s ∙ m

4. En un tubo de 0.15 m de longitud que contiene N 2 y acetona gaseosos a presión total de 2 atm,
se verifica una contradifusión equimolar en estado estacionario. La presión parcial N 2 es 1 090
mmHg en un extremo y 1 300 mmHg en el otro.

a) Calcule el flujo específico en kgmol /s ∙ m2 a 135°C para el N 2 .

b) Repita a 250°C. ¿Hay un aumento del flujo específico?

c) Repita a 135°C, pero a una presión total de 5 atm. La presión parcial de N 2 permanece a 1 090 y
1 300 mmHg, como en el inciso a). ¿Hay algún cambio del flujo específico?

SOLUCIÓN:

a) Calculo del flujo específico a 135°C para el N2.

Haremos uso de la ecuación de difusión molecular en gases para contradifusión equimolar en

estado estacionario.

D AB
N A= ( p − pA 2)
R ∙ T ∙ z A1
Donde:
 3
Pa∙ m
R=8 314
kmol ∙ K
T =135 ° C=408.15 K
z=0.15 m
p A 1 =1090 mmHg=145321.4 pa

P A 2 =1300 mmHg=173 319.1 pa

Cálculo de D AB a 135 °C (408.15 K) y 2 atm

Ecuación de Chapman Enskog para el coeficiente de difusión:

(
1.8583 ×10−7 T 3/ 2
)
1 /2
1 1
D AB= +
2
P ∙ σ AB ∙ Ω D , AB P M A P MB

Donde:
2
D AB : difusividad en m / s
T :temperatura en K =135° C=408.15 K
P M A : peso molecular de A=28.0 2

P M B : peso molecular de B=58.08

P : presión absoluta en atm=2 atm

σ AB :diametro promedio de colisión en

Ω D , AB : integral de colisión basada en el potencial de Lennard−Jones para difusividad

Cálculo del diámetro promedio de colisión:

El diámetro de las sustancias se calcula con:

( )
1
TC 3
σ =2.44
PC

De tablas

Para el Nitrógeno (A)

T C =126.2 K

PC =33.5 atm=33.9438 ¯
¿
 ( )
1
TC 3
σ =2.44 =2.44 ¿ ¿
PC

De tablas

Para el Acetona (B)

T C =508.0 K

PC =47.0 atm=47.6228 ¯
¿

σ =2.44
( )
T C 13
PC
=2.44 ¿ ¿

Diámetro de colisión para mezcla:

1 1
σ AB= ( σ A +σ B )= ( 3 .7829+5.3712 ) Å
2 2
σ AB=4.5770 Å

Cálculo de la energía característica:

De tablas de predicción del coeficiente de difusión para el Nitrógeno (A):

ε
=71.4 k
k
De tablas de predicción del coeficiente de difusión para el Acetona(B):

ε
=560.2 k
k
Para la mezcla:

ε AB
k √
ε ε
= A ∙ B =√ (71.4 k ) ( 560.2k )
k k

ε AB
=199.9957 K
k
Cálculo de la función integral de colisión:
 ¿ T 408.15 K
T= = =2.04 07
ε AB /K 199.9957 K

Interpolando, se obtiene:

T¿ ΩD
2.00 1.0754
2.0407 1.0681
2.1 1.0576

La función integral de colisión es 1.0681

Sustituyendo los valores conocidos en la ecuación de Chapman – Enskog:

( )
−7 3/ 2 1 /2
1.8583 ×10 T 1 1
D AB= 2
+
P ∙ σ AB ∙ Ω D , AB P MA P MB

−7
D AB=1.8583 ×10 ¿ ¿
−6 2
D AB=7.8758 ×10 m / s

Sustituyendo los valores conocidos:

D AB
N A= ( p − pA 2)
R ∙ T ∙ z A1
Donde:
3
Pa∙ m
R=8 314
kmol ∙ K
T =135 ° C=408.15 K
z=0.15 m
p A 1 =1090 mmHg=145321.4 pa

P A 2 =1300 mmHg=173 319.1 pa

 −6 2
7.8758 ×10 m / s (145 321.4−173 319.1 ) pa
N A=
( )
3
Pa ∙ m
8 314 ( 408.15 K ) ( 0.15 m)
kmol ∙ K

N A =4.3320 ×10−7 kgmol/ s ∙ m2

b) Repita a 250°C (523.15 K). ¿Hay un aumento del flujo específico?

7.8758 ×10−6 m2 / s ( 145 321.4−173319.1 ) pa

N A=
( 8 314
Pa ∙ m3
kmol ∙ K )
( 523.15 K )( 0.15 m )

−7 2
N A =3.3797 ×10 kgmol /s ∙m

c) Repita a 135°C (408.15K), pero a una presión total de 5 atm. La presión parcial de N 2
permanece a 1 090 y 1 300 mmHg, como en el inciso a). ¿Hay algún cambio del flujo específico?

D AB ∙ PT
N A= ( p − pA2)
R ∙ T ∙ z ∙ pB , M A 1

p B 2− p B 1
pB , M =
ln
( )
pB 2
pB 1

Considerando la presión de 5 atm (506 625 pa)

pB 1 =PT − p A 1=( 506 625 pa−145 321.4 pa )=361 303.6 pa

pB 2=PT − p A 2=( 506 625 pa−173 319.1 pa ) =333 305.9 pa

Cálculo de pB , M
 p B 2− p B 1
pB , M =
ln
( )
pB 2
pB 1

Sustituyendo los valores:

(3 61303.6−3 33 305.9) pa
pB , M =
ln(361 303.6 pa
333 305.9 pa )
pB , M =347 116.5841 pa

Sustituyendo valores:

D AB ∙ PT
N A= ( p − pA2)
R ∙ T ∙ z ∙ pB , M A 1
−6 2
(7.8758 ×10 m / s)(506 625 pa)
N A= (1 73 319.1−145 321.4 ) pa
( 8 314
Pa ∙ m3
kmol ∙ K )
( 408.15 K ) ( 0.15 m) (347 116.5841 pa)

−7 2
N A =6.3227 ×10 kgmol/ s ∙m