1.

Hidrógeno gaseoso se almacena a presión elevada en un recipiente rectangular que tiene

paredes de acero de 10 mm de espesor. La concentración molar de hidrógeno en la superficie
interior del acero es de 1 Kmol/m 3, mientras que la concentración del hidrógeno en la
superficie exterior del acero es insignificante. El coeficiente de difusión binaria para el
hidrógeno en el acero es 0,26x 10-12 m2 / s. ¿Cuál es el flujo de difusión molar de hidrógeno a
través del acero?

Solución:

Partiendo de la definición de Flux total, tenemos que:

N=N A + N B

N B =0

El flux total relativo de hidrogeno a un plano esta dado por:

dCA
N A =C A U 0−D v
dZ

Puesto que no hay cantidad de A transportada por el flujo del fluido es decir no hay flujo
convectivo:

dCA
N A =−D v
dZ

Separando variables e integrando:

z CA

N A ∫ dZ=−D v ∫ d C A
0 CA 0

N A (Z−0)=−D v (C A−C A ) 0

D v (C A −C A )
N A= 0

Reemplazando valores:

(0,26 ×10−12 m2 /s)(1 kmol /m 3)

N A=
0,01 m
N A =2,6 ×10−11 kmol/m2 . s

2. Amoníaco gaseoso se difunde a través de N 2 bajo condiciones de estado estable con N 2 sin
difundirse ya que es insoluble en un límite. La presión total es de 1.013 x 10 5 Pa, y la
temperatura es de 298 K. La presión parcial de amoníaco en un punto es 1.333 x 10 4 Pa y en
el otro punto a 20 mm de distancia es 6,666 x 10 3 Pa. La difusividad de la mezcla a 1,013 x 10 5
Pa y 298 K es 2,3 x 10-5 m2/s.

a) Calcular el flujo de amoníaco en Kg-mol / m2s.

b) Hacer lo mismo que (a), pero suponiendo que también se N 2 se difunde, es decir, ambos
límites son permeables a los gases y el flujo es una contra-difusión equimolar. En cuál caso, es
el mayor flujo?

Solución:

Consideraciones:

-Estado estacionario (S.S).

- No hay difusión de B.

Datos:

P=1,013 ×10 5 Pa

T =298 K

P A 1=1,333× 104 Pa

P A =6,666 ×103 Pa

−5 m2
D v =2,3× 10
s

Z=z 2−z 1=B T =0,02m

a) Difusión en 1-vía

Partiendo de la definición de Flux total, tenemos que:

 N=N A + N B

N B =0

El flux molar de A se lleva acabo por el flujo convectivo del fluido y por
la difusión molecular , de estamanera :

dCA
N A =U 0 C A −D v
dz

Expresando laecuación en fracciones molares:

NA
U0= C =ρ M y A d C A =ρM d y A
ρM A

d yA
N A =N A y A −D v ρM
dz

d yA
N A −N A y A =−D v ρM
dz

d yA
N A (1− y A )=−D v ρ M
dz
z2 yA 2
dy A
N A ∫ dz=−Dv ρM ∫
z1 yA 1 (1− y A )

Si u=1− y A du=−dy A

Reemplazando e integrando para cada uno de los límites

z2
1− y A 2
N A ∫ dz=−D v ρM ln
z1
( 1− y A 1 )
PA1 PA2 P
y A1= y A 2= ρ M=
P P RT
 PA 2

( )
1−
P P
N A ( z 2−z 1 )=D v ln
RT P
1− A 1
P

P P−P A 2
N A ( z 2−z 1 )=D v
RT
ln(P−P A 1 )
Dv P P−P A 2
N A=
RTZ
ln (
P−P A 1 )
Reemplazando cada uno de los valores

m2−5 5
(2,3 ×10 )(1,013× 10 Pa)
s ( 1,013× 105 Pa−6,666 ×103 Pa )
N A=
(8,314
Pa. m 3
)(298 K )(0,02 m)
ln
(
( 1,013× 105 Pa−1,333 × 104 Pa ) )
mol . K

mol Kg mol
N A =3,4334 ×10−3 2
=3,4334 × 10−6
m s m2 s

b) Contra difusión equimolar

Partiendo de la definición de Flux total de la especie A, tenemos que:

N A =U 0 C A + J A

dCA
N A =U 0 C A −D v
dz

Enla contra difusión equimolar no hay flujo convectivo del fluido , sinoque el flux molar
de A se debe a ladifusión molecular , de esta manera :

dCA
N A =−D v
dz
 z2 C A2

N A ∫ dz=−Dv ∫ d C A
z1 C A1

C A= ρ M y A

dC A =ρM dy A

BT yA

N A ∫ dz=−D v ρM ∫ d y A
0 yA 1

N A . B T =−D v ρM ( y A − y A 1 )

Teniendo en cuenta que para un gas ideal se cumple que:

P
ρ M =C A=
RT

PA
y A=
P

−Dv . P P A P A 1
N A=
RT BT P (−
P )
m2 −5 5
(−2,3 ×10 )(1,013× 10 Pa)
s 6,666 × 103 Pa 1,333× 104 Pa
N A=
(8,314
Pa. m 3
)(298 K )(0,02 m)
( −
1,013 ×10 5 Pa 1,013 ×105 Pa )
mol . K

mol
N A =3,09318 ×10−3
m2 s

Kgmol
N A =3,09318 ×10−6
m2 s

En conclusión : El flux molar de A en una sola vía es mayor que cuando se presenta
contradifusión equimolar .
3. Helio y nitrógeno gaseosos se encuentran contenidos en un conducto de 5 mm de diámetro
y 0,1 m de largo a 298 K y a una presión constante uniforme de 1 atm. La presión parcial de
Helio en un extremo del tubo es de 0,06 atm y en el otro extremo es 0,02 atm. Usando un valor
de difusividad de 0.687 x 10 -4 m2/s. Calcular para estado estacionario en contra-difusión
equimolar:

a) El flujo de Helio en Kg-mol/ m2s

b) El Flujo de Nitrógeno en Kg-mol/ m2s

c) La presión parcial del Helio en un punto a 0,05 m desde cualquier extremo

Solución:

Condiciones del Sistema

A=Helio( He)

B=Nitrógeno( N 2 )

T =298 K
P=1 atm

P A 1=0,06 atm

P A 2 =0,02 atm

−4 m2
D v =0,687 ×10
s

Z=0,1 m

Consideraciones

- Estado Estacionario (S.S).

- Asumiendo Contradifusión Equimolar.

a) Partiendo de la definición de Flux total, tenemos que:

N=N A + N B

N A =U 0 C A + J A

Como se presenta solocontra difusión equimolar no hay flujo convectivo del fluido

N A =J A

dCA
N A =−D v
dz
C A= ρ M y A

dC A =ρ M dy A

BT yA

N A ∫ dz=−Dv ρM ∫ d y A
0 yA 1

N A . B T =−D v ρM ( y A − y A 1 )

Teniendo en cuenta que para un gas ideal se cumple que:

 P
ρ M =C A=
RT
PA
y A=
P

−Dv . P P A 2 P A 1
N A=
RT BT P ( −
P )
Reemplazando cada uno de los valores

−4 m2
(−0,687× 10 )(1 atm)
s 0,02 atm 0,06 atm
N A=
(8,205 ×10−5
atm m 3
)(298 K)(0,1m)
( 1 atm
−
1 atm )
K mol

mol
N A =1,1238 ×10−3
m2 s

Kgmol
N A =1,1238 ×10−6
m2 s

b)

Partiendo de la definición de Flux total de la especie B, tenemos que:

N B =U 0 C B +J B

No hay flujo convectivo del componente B por el movimiento del fluido , por tanto :

N B =J B

d CB
N B =−D v
dz

C B= ρM y B dC B =ρM dy B
 Z y B2

N B ∫ dz =−Dv ρM ∫ d y B
0 y B1

N B . Z=−D v ρM ( y B− y B 1)

Teniendo en cuenta que para un gas ideal se cumple que:

P
ρ M =C B =
RT

PB
yB =
P

−Dv . P P B 2 P B 1
NB=
RTZ (P
−
P )
P=P A 1+ P B 1 P=P A 2 + PB 2

PB 1 =P−P A 1 PB 2=P−P A 2

PB 1 =1 atm−0,06 atm PB 2=1 atm−0,02 atm

PB 1 =0,04 atm PB 2 =0,08 atm

m2 (
(−0,687 ×10−4
s )
1 atm )
0,08 atm 0,04 atm
NB=
atm m 3 ( 1 atm
−
1 atm )
(8,205 ×10 −5
K mol ) ( 298 K )( 0.1 m )

mol
N B =−1,1238 ×10−3
m2 s

Kmol
N B =−1,1238 ×10−6
m2 s

c)

El flux molar de Helio viene dado por

N A =J A

dCA
N A =−D v
dz

C A= ρM y A dC A=ρ M dy A

z yA

N A ∫ dz=−Dv ρM ∫ d y A
0 yA 1

N A . Z =−Dv ρM ( y A − y A 1 )

Teniendo en cuenta que para un gas ideal se cumple que:

P
ρ M =C A=
RT
PA
y A=
P

−Dv . P P A 2 P A 1
N A=
RTZ P ( −
P )
En este caso Z=0,05 m

P A 2 =Presión en Z=0,05 m

P A 1=Presión en Z=0 m

La presión parcial de He en un punto de 0,05 m desde cualquier extremo es :

NAP
PA2= + PA1
−D v P
( RTZ )
N A RTZ
PA2= + PA1
(−Dv )
 3
mol −5 atm m
(1,1238× 10−3 )(8,205 ×10 )(298 K)( 0,05 m)
m2 s K mol
PA2= 2
+0,06 atm
−4 m
−0,687 × 10
s

P A 2 =0,04001 atm

4. Estimar la difusividad de las siguientes mezclas de gases:

a) Acetona-aire a 1 atm. y 0 ° C.
b) Nitrógeno-dióxido de carbono a 1 atm. y 25 ° C.
c) Cloruro de Hidrógeno-aire a 200 kN/m2 y 25 ° C.
d) Benceno en aire a 2 atm. y 100 ° C.

Solución:
Utilizando la modificación de wilke-lee del método de Hirschfelder-Bird-Spotz para calcular la
difusividad de la mezcla de gases:

1 1 1 1

D AB=
(
10−4 1,084−0,249
√ +
MA MB ) √
T 3 /2 +
MA MB
Pt (r ¿¿ AB)2 f (kT /ε AB)(1) ¿

Donde:

M A =peso molecular de A , Kg/ Kmol

M B =peso molecular de B , Kg/ Kmol
Pt = presion absoluta , N /m2
r A + rB
r AB =separaccion molecular durante el choque= ,nm
2
ε AB=energia de atraccion molecular=√ ε A ε B
T =temperatura , K
k =constante de Boltzmann
kT
f ( )
ε AB
=funcion de choque dada por la figura 2.5

a) Acetona-aire a 1 atm. y 0 ° C
Datos:

A: Acetona (C3H6O)
B: Aire

M A =58,04 Kg/Kmol
M B =29 Kg/Kmol
Pt =101325 N /m2
r B=0,3711nm
T =273,15 K
r A =1,18 v A 1/3

Donde: v=volumen molar del líquido en el punto de ebullición normal

Utilizando los datos de la tabla 2.3 Robert E. Treybal:

m3
v A =3 0,0148 +6 0,0031 +0,0014=0,0644
( ) ( )
Kmol
1 /3
r A =1,18(0,0644 ) =0,4729 nm
 r A +r B 0,4729+0,3711
r AB = = =0,4220 nm
2 2

εB
=78,6
k

εA
=1,21 T b
k

Donde T b es la temperatura de ebullicion del compuesto .

εA
=1,21 ( 329,4 K ) =398,574
k

ε AB=√ ( 78,6 )∗( 398,574 )=176,997 k

kT k 273,15
= ≈ 1,543
ε AB 176,997 k

De la figura 2,5: Función de choque para la difusión Robert E. Treybal.

kT
f ( )
ε AB
≈ 0,59

1 1 1 1
√ +
M A MB
=
√
+ =0,227
58,04 29

Reemplazando en la ecuación (1)

3
−4 2
10 ∗( (1,084−0,249 )∗( 0,227 )) ¿ ( 273,15 ) ∗(0,227)
D AB= 2
=9,8892 ×10−6 m2 / s
( 101325 ) ¿ ( 0 , 4220 ) ∗(0,59)
 b) Nitrógeno-dióxido de carbono a 1 atm. y 25 ° C.

Datos:

A: Nitrógeno (N2)

B: Dióxido de carbono (CO2)

M A =28 Kg /Kmol
M B =44 Kg / Kmol
Pt =101325 N /m2
r A =0,3798 nm
r B=0,3941 nm
T =298,15 K

r A +r B 0,3798 nm+0,3941 nm
r AB = = =0,3827 nm
2 2

εB
=195,2 K
k

εA
=1,21 T b=1,21 ( 77,2 K ) =78,41
k

ε AB=√ ( 195,2 )∗( 78,41 )=123,716 k

kT k 298,15
= ≈ 2,4099
ε AB 123,716 k

De la figura 2,5: Función de choque para la difusión Robert E. Treybal.

kT
f ( )
ε AB
≈ 0,457

1 1 1 1
√ +
M A MB
=
√+ =0,2417
28 44

Reemplazando en (1)
 3
10−4∗( (1,084−0,249)∗( 0,2417 ) ) ¿ ( 298,15 ) 2∗(0,2417)
D AB= 2
=1,878 ×10−5 m2 /s
( 101325 )∗( 0 ,3827 ) ∗(0,457)

c) Cloruro de hidrógeno-aire a 200 kN/m2 y 25 ° C.

Datos:

A: Cloruro de Hidrógeno (HCl)

B: Aire

M A =36,6 Kg /Kmol
M B =29 Kg/Kmol
Pt =200000 N /m2
r A =0,3339
r B=0,3711nm
T =298,15 K

r A +r B 0,3339 nm+ 0,3711 nm

r AB = = =0,3525 nm
2 2

εA
=1,21 T b=1,21 ( 188,15 K ) =227,66
k
εB
=78,6 K
k

ε AB=√ ( 227,66 )∗( 78,6 )=133,769 k

kT k 298,15
= ≈ 2,23
ε AB 133,769 k

De la figura 2,5: Función de choque para la difusión Robert E. Treybal.

f ( εkT ) ≈ 0,5
AB
 1 1 1 1
√ +
M A MB
=
√
+ =0,2488
36,46 29

Reemplazando en (1)

D AB=10−4∗¿ ¿

d) Benceno en aire a 2 atm. y 100 ° C.

Datos

A: Benceno
B: Aire

M A =78,11 Kg /Kmol
M B =29 Kg/Kmol
Pt =202650 N /m2
r A =0,5349 nm
r B=0,3711nm

r A +r B 0,5349 nm+ 0,3711 nm

r AB = = =0,453 nm
2 2

εA
=412,3
k
εB
=78,6 K
k

ε AB=√ ( 412,3 )∗( 78,6 )=180,018 k

kT k 373.15
= ≈ 2,1
ε AB 180,018 k

De la figura 2,5: Función de choque para la difusión Robert E. Treybal.

kT
f ( )
ε AB
≈ 0,52

1 1 1 1
√ +
M A MB
=
√
+ =0,2175
78,11 29
 Reemplazando en (1)
3
−4 2
10 ∗( (1,084−0,249)∗( 0,2175 ) )∗( 373,15 ) ∗(0,2175)
D AB= 2
=7,464 ×10−6 m2 /s
( 202650 )∗( 0,453 ) ∗(0,52)