3.8.

- Transporte de masa de fluidos en solidos porosos

El estudio de los fenómenos de transferencia de masa en solidos porosos

La ecuacion de balance de materia del transporte de un fluido en el seno de un sólido, está dado
por:

 Entrada de materia   Salida de materia  Generacion de   Acumulacion  (3.32)

 por difusion molecular  −  por difusion molecular  −  de materia  =  de materia 
       

En este item, estudiaremos el transporte de masa en régimen estacionario y sin generación de

materia (ausencia de reacción química)
3.8.1- Transporte de masa en un sólido poroso de pared plana.-
El transporte de masa en el interior de un sólido poroso, donde no se produce reacción química,
solo se produce por difusión molecular proceso espontáneo); por tanto, es aplicable la ley de Fick.

Figura 3.9.- Transporte de masa unidimensional a través de un elemento de volumen en una pared
plana grande

Para plantear el balance de materia, tomaremos un volumen de control WLx (Figura 3x), luego
rebatimos una porción para señalar la concentración en cada cara así como la dirección de la
transferencia de materia por difusión molecular.

J Ax x WL − J Ax x +x
WL = 0 (3.33)

Llevando al limite

J Ax WL − J Ax x WL
lim x +x
=0 (3.34)
x →0 WLx
Esto permite obtener:

dJ Ax
=0 (3.35)
dx
La ecuacion (3.35) es la ecuacion diferencial de transporte de masa unidimensional en
coordenadas rectangulares

La ley de Fick para difusión molecular

 dCA
J Ax = − DAe (3.36)
dx
Donde DAe es el coeficiente de difusión del componente que se difunde a través del material
poroso.

Reemplazando (3.36) en (3.35)

dC A
d (− DAe )
dx = 0 (3.37)
dx
Para el caso de difusividad efectiva constante, la ecuacion (3.37) se transforma en:

d 2C A
=0 (3.38)
dx 2
La ecuacion (3.38) se puede hacer extensiva para los tres ejes

 2C A  2C A  2C A
+ + =  2C A = 0 (3.39)
x 2
y 2
z 2

Problema 3.8

Transporte de materia difusional en un sólido poroso.

Considerando la ecuacion diferencial de transferencia de masa unidimensional en coordenadas

rectangulares y conociendo las condiciones de frontera de primera clase (condiciones de
frontera de Dirichlet), encuentre:

a) El perfil de concentración
b) La densidad de flujo molar
c) La concentración media

Solución

La ecuacion de transferencia de masa unidireccional en coordenadas rectangulares, en régimen

estacionario y sin generación de materia, esta dada por:

d 2C A
=0 (1)
dx 2
La integración de la ecuacion (1) conduce a:

C A = C1 x + C2 (2)

Con las condiciones de frontera

x = x1 C A = C A1 y x = x2 C A = C A2

Reemplazando las condiciones de frontera resulta:

C2 = C A1 y C1 = (C A2 − C A1 ) / ( x2 − x1 )
Reemplazando en (2) tenemos:

x
C A = C A1 − (C A1 − C A 2 )  (3)
( x2 − x1 )

La ecuacion (3) es la ecuacion de distribución de concentración. Esta ecuacion se puede expresar

de forma adimensional del siguiente modo:

C A − C A1 x
= (4)
(C A2 − C A1 ) ( x2 − x1 )

La ecuacion (4) escrita de forma mas simplificada es la siguiente

CA* = x* (5)

Quiere decir que la concentración adimensional es igual a la longitud adimensional.

Las ecuaciones (4) y (5) pueden ser escritas en términos de la presion parcial y las fracciones
molares:

En términos de la presion parcial

p A − p A1 x
= p*A = x* (6)
( p A 2 − p A1 ) ( x2 − x1 )

En términos de la fracción molar

y A − y A1 x
= y*A = x* (7)
( y A 2 − y A1 ) ( x2 − x1 )

La densidad de flujo molar se obtiene a partir de (3.40), teniendo en consideración que dicho
valor es constante a través de toda la sección transversal de la pared plana. Es decir:

J Ax x = x1
= J Ax x= x 2
(8)

En consecuencia solo basta conocer uno de los términos (8), pues dicho valor es constante.

dC A (C A1 − C A2 )
J Ax = − DAe = − DAeC1 = DAe (9)
x = x1
dx x = x1 ( x2 − x1 )

El flujo molar es calculado multiplicando la densidad de flujo molar por el área de transferencia
de masa. El termino DAe / ( x2 − x1 ) , se conoce como coeficiente de transferencia de masa km ,
de tal forma que (9) puede ser escrita como:

J Ax x = x1
= km (C A1 − C A2 ) (20)

Luego el flujo molar en términos del coeficiente de transferencia de masa será:

QAx x = x1
= J Ax x = x1
WL = kmWL(C A1 − C A2 ) (11)

También se debe tener en cuenta que los flujos molares son constantes, es decir:

QAx x = x1
= QAx x= x 2
(12)
Puesto que el balance de materia esta fundamentada en esta ecuacion (24)

La concentración promedio puede ser encontrada a partir de:

 x 

x2
C A1 − (C A1 − C A2 )  ( x − x )  dx
   LW  C A dx
L W x2 x2
C A dzdydx x1
 1  1
 C A = 0 0 x1
= x1
= 2
= (C A1 + C A 2 ) (13)
LW ( x2 − x1 ) ( x2 − x1 )
  
L W x2
dzdydx 2
0 0 x1

Hasta aquí esta mas o menos arreglada ecuaciones, tablas, figuras

3.8.2- Transporte de masa en un sólido poroso de pared cilíndrica

Consideremos un cilindro de pared porosa en la parte lateral e impermeable en las tapas
superior e inferior, tal como se muestra en la Figura 3.10. En el interior del cilindro, la presion
parcial es del gas es PA1 y en el exterior PA 2 . Conociendo los radios interno y externo del cilindro
desarrollaremos el balance de materia.

Para representar de mejor manera las variables que intervienen en este proceso de
transferencia de masa, se representa el mismo esquema visto desde la parte superior del
cilindro.

Tomando como base la ecuacion (3.32); en régimen estacionario y sin generación de materia, se
desarrolla el balance de materia respectivo.
J Ar r 2 rL − J Ar r +r
2 rL = 0 (3.40)

Llevando al limite

Jr Ar +r
2 rL − J Ar 2 rL
lim r
=0 (3.41)
r →0 2 Lr
Simplificando la ecuacion (3.41) se tiene

d (rJ Ar )
=0 (3.42)
dr
La ley de Fick esta dada por:

dCA
J Ar = − DAe (3.43)
dr
Reemplazando (3.43) en (3.42)

dC A
d (−rDAe )
dr = 0 (3.44)
dr
Dado que el coeficiente de difusión es constante, (3.44) se transforma en:

d dCA
(r )=0 (3.45)
dr dr
La ecuacion (3.45) se puede integrar en términos de la concentración molar; sin embargo, dado
que el sistema es gaseoso, la ecuacion (3.45) se expresa en términos de la presion parcial y luego
se integra. Sabemos que PA = C A RT → dPA = RTC A , luego

d r dpA
( )=0 (3.46)
dr RT dr
Dado que el proceso se lleva a cabo en condiciones isotérmicas, (3.46) se transforma en:

d dpA
(r )=0 (3.47)
dr dr
Las ecuaciones (3.45) y (3.47), representan la ecuacion diferencial de transferencia de masa en
una pared porosa de geometría cilíndrica
Problema 3.9
Transferencia de masa por difusion unidireccional en una pared porosa de pared cilíndrica
Considere que se tiene un cilindro poroso y cuya transferencia de masa se da solo por la pared
lateral (las paredes de las bases son impermeables); en la cual se difunde el gas A; se conoce las
dimensiones de la pared porosa, así como las propiedades físicas y las condiciones a loas que se
encuentra sometida; en base a esta información, encuentre

a) La ecuacion de distribución de presiones parciales

b) La densidad de flujo molar
c) La presion media del gas
Solución

La ecuacion diferencial de transferencia de masa para este sistema expresado en términos de

presiones parciales viene dada por la ecuacion (3.51), es decir:

d dpA
(r )=0 (1)
dr dr
Cuya integración conduce a:

PA = C1 ln r + C2 (2)

Las condiciones de frontera son:

r = R1 PA = PA1 y r = R2 PA = PA2

( PA1 − PA 2 )
Con estas condiciones las constantes C1 y C2 son: C1 = y C2 = PA1 − C1 ln R1 con
ln( R1 / R2 )
lo cual la ecuacion (2) se transforma en:

( PA2 − PA1 )
PA = ln(r / R1 ) + PA1 (3)
ln( R2 / R1 )

a) La ecuacion (3) es la ecuacion de distribución de presiones parciales en la pared

cilíndrica. Por otro lado, la ecuacion (3) puede ser modificada para expresarla en forma
adimensional

PA − PA1 ln(r / R1 )
= (4)
PA2 − PA1 ln( R2 / R1 )

Ahora podemos usar los siguientes cambios de variable para ser reemplazadas en la ecuacion (4):
PA* = ( PA − PA1 ) / ( PA2 − PA1 ) , r * = r / R1 y A = 1/ ln( R2 / R1 ) (4) se transforma en

PA* = A ln r * (5)

b) La densidad de flujo molar en la pared interna se obtiene mediante:

DAe dPA DAe C1 DAe 1 ( PA1 − PA2 )

J Ar r = R1
=− =−  =  = km ( PA1 − PA2 ) (6)
RT dr r = R1 RT R1 RT R1 ln( R2 / R1 )

1 DAe
Con km = , siendo km el coeficiente de transferencia de masa
RTR1 ln( R2 / R1 )

Del mismo modo, la densidad de flujo molar en la pared externa esta dada por:

DAe dPA DAe C1 DAe 1 ( PA1 − PA2 )

J Ar r = R2
=− =−  =  = km ( PA1 − PA2 ) (7)
RT dr r =R2 RT R2 RT R2 ln( R2 / R1 )

1 DAe
Con km = , siendo km el coeficiente de transferencia de masa
RTR2 ln( R2 / R1 )

Como era de esperarse, las densidades de flujo molar interno y externo, no son iguales, es decir:
J Ar r = R1
 J Ar r = R2
(8)

Los flujos molares son iguales, pues en la igualdad del flujo de masa radica el balance de materia
en régimen estacionario y sin reacción química.

DAe 1 ( PA1 − PA2 ) D (P − P )

QAr r = R1
= 2 R1 L  = 2 L Ae  A1 A2 (9)
RT R1 ln( R2 / R1 ) RT ln( R2 / R1 )

DAe 1 ( PA1 − PA2 ) D (P − P )

QAr r = R2
= 2 R2 L  = 2 L Ae  A1 A2 (10)
RT R2 ln( R2 / R1 ) RT ln( R2 / R1 )

QAr r = R1
= QAr r = R2
(11)

Nota: Los coeficientes de transferencia de masa dependen de tres factores: Las condiciones de
operación (en este caso la temperatura), la geometría (en este caso los radios interno y externo)
y la propiedad física (en este caso el coeficiente de difusión efectiva del gas que difunde)

La ecuacion (9) se puede volver a escribir como:

DAe ( PA1 − PA2 ) ( R2 − R1 )

QAr = 2 L   (12)
r = R1
RT ln( R2 / R1 ) ( R2 − R1 )

Esta ecuacion se puede transformar en:

DAe ( PA1 − PA2 ) 2 L( R2 − R1 ) DAe ( PA1 − PA2 ) A2 − A1

QAr =   =   (13)
r = R1
RT ln(2 LR2 / 2 LR1 ) ( R2 − R1 ) RT ln( A2 / A1 ) ( R2 − R1 )

DAe ( A2 − A1 )
QAr =  ( PA1 − PA2 ) = km  ( PA1 − PA2 ) (14)
r = R1
RT ( R2 − R1 ) ln( A2 / A1 )

también se puede expresar en términos de concentración molar

DAe ( A2 − A1 )
QAr =  (C A1 − C A2 ) = km  (C A1 − C A2 ) (15)
r = R1
( R2 − R1 ) ln( A2 / A1 )

Del mismo modo se aplica a la ecuacion (15) y se obtiene:

QAr r = R2
= km  ( PA1 − PA2 ) = km(C A1 − C A2 ) (16)

c) La concentración media se obtiene mediante:

2
  rP dr =
R
A 2  ( PA2 − PA1 ) 

R
 pA = ln(r / R1 ) + PA1  rdr
0 R1
  (17)
− R12
  rdr R
2 R2 2
2
R1
 ln( R2 / R1 ) 
0 R1

La ecuacion (17) es la ecuacion que permite calcular la presion media.

2   R  ( PA2 − PA1 )  

 pA = R  ln(r / R1 )  rdr + PA1rdr  (18)
R22 − R12 
 1  ln( R2 / R1 )  

 2  ( PA 2 − PA1 )  2
 R1 rdr
R R2
 p A =  ln(r / R1 )  rdr + 2 P (19)
R2 − R12 R2 − R12
2 A1
R1
 ln( R2 / R1 ) 

2 ( PA2 − PA1 ) R2 2 R2

R22 − R12 ln( R2 / R1 ) R1

 p A = r ln( r / R ) dr + PA1 R rdr (20)
R22 − R12
1
1

R2
2 ( PA2 − PA1 ) r 2 r2 r2
 pA = 2 ( ln r − − ln R1 ) + PA1 (21)
R2 − R12 ln( R2 / R1 ) 2 4 2 R 1

 R2 1 
 p A =  2 2 2 −  ( PA 2 − PA1 ) + PA1 (22)
 R2 − R1 2 ln( R2 / R1 ) 

Problema 3.10

Aplicación numérica del cálculo del perfil de presiones, densidad de flujo molar, flujo molar y
presion media en una pared porosa de geometría cilíndrica.

En el interior de un recipiente cilíndrico de pared lateral poroso de radio interno R1 = 20 cm,

radio externo R2 = 21 cm y L = 80 cm de altura, se encuentra encerrado gas metano; la
presion del gas en el interior del cilindro es 760 mmHg y en la parte externa 10 mmHg. Tomando
en consideración que el coeficiente de difusión efectivo es DAe igual a 110−6 m2/s, encuentre
los siguientes parámetros

a) El perfil de presión parcial

b) La densidad de flujo molar interno y externo
c) El flujo molar interno y externo (que deben resultar iguales)
d) La presion promedio

Solución

El programa desarrollado es el siguiente:

d(Y)/d(r) = 1
Y(0) = 0
PA = ((PA2 - PA1) / ln(R2 / R1)) * ln(r / R1) + PA1
Jar1 = ((De / (R * T)) * (1 / R1) * (PA1 - PA2) / ln(R2 / R1)) * 1000
Jar2 = ((De / (R * T)) * (1 / R2) * (PA1 - PA2) / ln(R2 / R1)) * 1000
Qar1 = Jar1 * 2 * pi * R1 * L
Qar2 = Jar2 * 2 * pi * R2 * L
PAm = ((R2 ^ 2 / (R2 ^ 2 - R1 ^ 2)) - (1 / (2 * ln(R2 / R1)))) * (PA2 - PA1) + PA1
PA1 = 760
PA2 = 76
R2 = 0.21
R1 = 0.20
R = 62.4
T = 298
L = 0.8
pi = 3.1416
De = 1E-6
r(0) = 0.20
r(f) = 0.21

a) El perfil de presion parcial es el siguiente

 b) Las densidades de flujo molar interno y externo expresados en mol/m3s, son 0.0037696
y 0.0035901, respectivamente
c) El flujo molar resulta: 0.0037896mol/s
d) La presion promedio igual a 412.4388 mmHg

3.8.3- Transporte de masa en un sólido poroso de pared esférica

Consideremos un recipiente esférico cuya pared es porosa en las tapas superior e inferior, tal como
se muestra en la figura xxx. En el interior de la esfera la presion parcial es del gas es PA1 y en el
exterior PA 2 . Conociendo los radios interno y externo de la esfera desarrollaremos el balance de
materia.
Tomando como base la ecuacion (3.32); en régimen estacionario y sin generación de materia, se
desarrolla el balance de materia respectivo.

J Ar r 4 r 2 − J Ar r +r
4 r 2 = 0 (3.48)

Llevando al limite

Jr Ar +r
4 r 2 − J Ar 4 r 4
lim r
=0 (3.49)
r →0 4r
Simplificando la ecuacion (3.49) se tiene

d (r 2 J Ar )
=0 (3.50)
dr
La ley de Fick esta dada por:

dCA
J Ar = − DAe (3.51)
dr
Reemplazando (3.51) en (3.50)

dC A
d (−r 2 DAe )
dr = 0 (3.52)
dr
Dado que el coeficiente de difusión es constante, (3.52) se transforma en:

d 2 dCA
(r )=0 (3.53)
dr dr
La ecuacion (3.53) se puede integrar en términos de la concentración molar; sin embargo, dado
que el sistema es gaseoso, la ecuacion (3.53) se expresa en términos de la presion parcial y luego
se integra. Sabemos que PA = C A RT → dPA = RTC A , luego
d r 2 dpA
( )=0 (3.54)
dr RT dr
Dado que el proceso se lleva a cabo en condiciones isotérmicas, (3.54) se transforma en:

d 2 dpA
(r )=0 (3.55)
dr dr
La ecuacion (3.55), es la ecuacion diferencial de transferencia de masa en paredes cilíndricas
porosas.

Problema 3.11

Aplicación numérica del cálculo del perfil de presiones, densidad de flujo molar, flujo molar y
presion media en una pared porosa de geometría esférica.

En el interior de un recipiente esférico de pared lateral porosa de radio interno R1 = 20 cm y

radio externo R2 = 21 cm de altura, se encuentra encerrado gas metano; la presion del gas en
el interior del cilindro es 760 mmHg y en la parte externa 10 mmHg. Tomando en consideración
que el coeficiente de difusión efectivo es DAe igual a 110−6 m2/s, encuentre los siguientes
parámetros

e) El perfil de concentración
f) La densidad de flujo molar interno y externo
g) El flujo molar interno y externo (que deben resultar iguales)
h) La presion promedio

Solución

La ecuacion de transferencia de masa para este caso, está dado por la ecuacion (3.59), es decir

d 2 dpA
(r )=0 (1)
dr dr
Cuya integración conduce a:

C1
PA = − + C2 (2)
r
Las condiciones de frontera son:

r = R1 PA = PA1 y r = R2 PA = PA2

( PA1 − PA 2 )
Con estas condiciones las constantes C1 y C2 son: C1 = − y C2 = PA1 + C1 / R1
(1/ R1 − 1/ R2 )
con lo cual la ecuacion (3) se transforma en:

(1/ R1 − 1/ r )
PA = ( PA1 − PA2 ) + PA1 (4)
(1/ R2 − 1/ R1 )
La ecuacion (4) es la ecuacion de distribución de presiones parciales en la pared cilíndrica. Por
otro lado, la ecuacion (4) puede ser modificada para expresarla en forma adimensional

PA − PA1 (1/ R1 − 1/ r )
=− (5)
PA2 − PA1 (1/ R1 − 1/ R2 )

Ahora podemos usar los siguientes cambios de variable para ser reemplazadas en la ecuacion (5):
PA* = ( PA − PA1 ) / ( PA2 − PA1 ) , r* = (1/ r ) / (1/ R1 −1/ R2 ) y A = −1/ R1 / (1/ R1 − 1/ R2 ) (5) se
transforma en

PA* = A + B / r * (6)

Nota: La ecuacion (2) difiere ligeramente de (6) debido a que la primera esta expresada en
términos absolutos y la segunda en forma adimensional.

La densidad de flujo molar en la pared interna se obtiene mediante:

DAe dPA DAe C1 DAe 1 ( PA1 − PA2 )

J Ar =− =−  =  = km ( PA1 − PA2 ) (7)
r = R1
RT dr r = R1 RT R12 RT R12 (1/ R2 − 1/ R1 )

1 DAe
Con km = , siendo km el coeficiente de transferencia de masa
RTR1 (1/ R2 − 1/ R1 )
2

Del mismo modo, la densidad de flujo molar en la pared externa esta dada por:

DAe dPA DAe C1 DAe 1 ( PA1 − PA2 )

J Ar =− =−  =  = km ( PA1 − PA2 ) (8)
r = R2
RT dr r = R2 RT R22 RT R22 (1/ R1 − 1/ R2 )

1 DAe
Con km = , siendo km el coeficiente de transferencia de masa
RTR2 (1/ R1 − 1/ R2 )
2

Como era de esperarse, las densidades de flujo molar interno y externo, no son iguales, es decir:

J Ar r = R1
 J Ar r = R2
(9)

Los flujos molares son iguales, pues en la igualdad del flujo de masa radica el balance de materia
en régimen estacionario y sin reacción química.

4 R12 DAe 1 ( PA1 − PA2 ) 4 DAe ( PA1 − PA2 )

QAr = 4 R J Ar2
=  2 =  (10)
r = R1 r = R1
R1 (1/ R1 − 1/ R2 ) (1/ R1 − 1/ R2 )
1
RT RT

4 R22 DAe 1 ( PA1 − PA2 ) 4 DAe ( PA1 − PA2 )

QAr = 4 R22 J Ar =  2 =  (11)
r = R2 r = R2
RT R2 (1/ R1 − 1/ R2 ) RT (1/ R1 − 1/ R2 )

QAr r = R1
= QAr r = R2
(12)

Nota: Los coeficientes de transferencia de masa dependen de tres factores: Las condiciones de
operación (en este caso la temperatura), la geometría (en este caso los radios interno y externo)
y la propiedad física (en este caso el coeficiente de difusión efectiva del gas que difunde)

La ecuacion (10) se puede volver a escribir como:

 4 DAe ( PA1 − PA 2 ) 4 R2 R1 DAe ( PA1 − PA 2 ) DAe A1 A2
QAr =  =  = ( PA1 − PA 2 ) (13)
r = R1
RT (1/ R1 − 1/ R2 ) RT ( R2 − R1 ) RT ( R2 − R1 )

Escribiendo (13) en términos del área media geométrica.

DAe Amg
QAr = ( PA1 − PA 2 ) = km ( PA1 − PA2 ) = km (C A1 − C A 2 ) (14)
r = R1
RT ( R2 − R1 )

Del mismo modo

QAr r = R2
= km  ( PA1 − PA2 ) = km(C A1 − C A2 ) (15)

La concentración media se obtiene mediante:

2 
   r P sen drd d =
R
2
A 3  C1  2
3 R
R
 pA = (16)
− r + C2  r dr
0 0 R1
 
− R1 1
   r sen drd d R
2 R2 3
2
2
0 0 R1

R2
3  C1 2 C2 3  3  C1 2 PA1 + C1 / R1 3 
 p A = 3 3 
− r + r  = 3  − 2 ( R2 − R1 ) +
2
( R2 − R13 )  (17)
R2 − R1  2 3  R1 R2 − R13 3 

3  C1 2 C ( R 3 − R13 ) 
 p A = 3 
− ( R2 − R12 ) + 1 ( R23 − R13 ) + 2 PA1  (18)
R − R1  2
3
2 3R1 3 

 3C ( R 2 − R12 ) C1  1 1 3 ( R22 − R12 )

 p A =  − 1 23 + + PA1  = ( − )( PA1 − PA 2 ) + PA1 (19)
 2 R2 − R1  (1/ R2 − 1/ R1 ) R1 2 R2 − R1
3 3 3
R1

La ecuacion (19) es la ecuacion que permite calcular la presion media.

El programa desarrollado es el siguiente:
d(Y)/d(r) = 1
Y(0) = 0
PA = (PA2 - PA1) * (1 / R1 - 1 / r) / (1 / R1 - 1 / R2) + PA1
Jar1 = ((De / (R * T)) * (1 / R1 ^ 2) * (PA1 - PA2) / (1 / R1 - 1 / R2)) * 1000 # Ecuacion (7)
Jar2 = ((De / (R * T)) * (1 / R2 ^ 2) * (PA1 - PA2) / (1 / R1 - 1 / R2)) * 1000 # Ecuacion (8)
Qar1 = Jar1 * 4 * pi * R1 ^ 2 # Ecuacion (10)
Qar2 = Jar2 * 4 * pi * R2 ^ 2 # Ecuacion (11)

QAr1A = kmA * (PA1 - PA2) # Ecuacion (13)

kmA = ((De/(R * T)) * Amg / (R2 - R1))*1000

Amg = (A1 * A2) ^ 0.5
A1 = 4 * pi * R1 ^ 2
A2 = 4 * pi * R2 ^ 2

PAm = (1/R1 - (3 / 2) * (R2 ^ 2 - R1 ^ 2) / (R2 ^ 3 - R1 ^ 3)) * (PA1 - PA2) / (1 / R2 - 1 / R1) + PA1 # Ecuacion (19)
PA1 = 760
PA2 = 76
R2 = 0.21
R1 = 0.20
R = 62.4
T = 298
pi = 3.1416
De = 1E-6 # m2/s
r(0) = 0.20
r(f) = 0.21
Resultados
a) El perfil de presion parcial es el siguiente

b) Las densidades de flujo molar interno y externo expresadas en mol/m 3s, son 0.0038623 y
0.0035032, respectivamente.
c) El flujo molar es: 0.0019414 mol/s
d) La presion promedio es 406.8803 mmHg
3.8.4 Difusión en una pared cónica
Consideremos una pared porosa cónica donde el gas se difunde en el sentido divergente-
convergente, tal como se muestra en la figura adjunta.

Para desarrollar el balance de materia en este sistema conico, se requiere esquematizarlo para
mostrar las dimensiones principales necesarias y sus respectivas relaciones geométricas, pues a
pesar de tratarse de un cono, el balance de materia se realizara en un sistema de coordenadas
rectangulares.
Hacemos el balance de materia según lo establecido en la ecuacion (3.32) es decir el balance de
materia difusional sin generación de materia y en régimen estacionario

J Ax A( x) x − J Ax A( x) x +x =0 (3.56)

Llevando al limite

J Ax A( x) x +x − J Ax A( x) x
lim =0 (3.57)
x →0 x
Por definición, la ecuacion (3.57) se transforma en:

d
( J Ax A( x)) = 0 (3.58)
dx
Dado que varia la densidad de flujo de masa como el area de la sección transversal, se aplica la
derivada del producto.

dJ Ax dA( x)
A( x) + J Ax =0 (3.59)
dx dx
El área de transferencia de masa puede ser expresada del siguiente modo

A( x) =  y 2 =  (mx + b) 2 (3.60)

Puesto que y y x están relacionados por la ecuacion de la recta. Por otro lado en base a la ley
de Fourier se tiene:

dCA
J Ax = − DAe (3.61)
dx
Los términos correspondientes de (3.59) son obtenidos a partir de (3.60) y (3.61)

Para el caso de difusividad constante, se tiene:

dJ Ax d 2C A
= − DAe (3.62)
dx dx 2
dA( x)
= 2 m(mx + b) (3.63)
dx
Reemplazando (3.62) y (3.63) en (3.59)
  d 2C A   dC A 
 (mx + b)   − DAe
2
2 
+  − DAe  2 m(mx + b) = 0 (3.64)
 dx   dx 

d 2C A 2m dCA
+ =0 (3.65)
dx 2
(mx + b) dx
La ecuacion (3.65), es la ecuacion de diferencial del transporte de masa en una pared cónica
divergente convergente.

Problema 3.12

Flujo molar en una pared porosa cónica divergente convergente

Encuentre la ecuacion para el cálculo del flujo molar del componente A que se desplaza en la
dirección axial en un cono divergente convergente; se conoce las dimensiones geométricas del
cono, las propiedades físicas del gas y las condiciones de trabajo.

Solución

La ecuacion de transferencia de masa correspondiente, es la ecuacion (3.65)

d 2C A 2m dCA
+ =0 (1)
dx 2
(mx + b) dx

dCA
Para resolver (1) hacemos el siguiente cambio de variables: = z , luego (1) se transforma
dx
en:

dz 2m
+ z=0 (2)
dx (mx + b)
La ecuacion (2) es fácilmente integrable pues se trata de una ecuacion diferencial de variables
separables.

ln z = −2ln(mx + b) + ln C1 (3)

A partir de (3) se tiene:

dCA
= z = C1 (mx + b)−2 (4)
dx
Integrando (4)

C1
CA = − + C2 (5)
m(mx + b)
Condiciones de frontera

x = 0 C A = C A0 y x = L C A = C AL

C A0 − C AL
C1 = (6)
1 1
−
m(mL + b) mb
 C1
C2 = CA0 + (7)
mb
Reemplazando (7) en (5)

C1 C  1 1 
CA = − + C A0 + 1 = C1  −  + C A0 (18)
m(mx + b) mb  mb m(mx + b) 
Ahora se reemplaza (6) en (8)

C A0 − C AL  1 1  C A0 − C AL 1 1  (19)
CA =  mb − m(mx + b)  + C A0 =  b − (mx + b)  + C A0
1 1   1 1  
− −
m(mL + b) mb (mL + b) b

Reordenando

1 1
−
C A − C A0 (mx + b) b  mL + b  x (20)
= = ( )
C AL − C A0 1 1  mx + b  L
−
(mL + b) b

La ecuacion (20), representa la ecuacion de distribución de concentración.

Respecto a la densidad de flujo de materia, se debe tener en cuenta que su valor varia en cada
x por tanto, haciendo uso de la ecuacion (3.65) se tiene:
(21)
dC A C A0 − C AL
J Ax x =0
= − DAe = − DAe (mx + b)−2
dx 1 1
x =0 −
m(mL + b) mb x =0

Evaluando x = 0

(mL + b)
J Ax x =0
= DAe (CA0 − CAL )
Lb (22)

 R12 (mL + b)
QAx x =0
= J Ax x =0
 R12 = DAe (CA0 − CAL ) (23)
Lb
Del mismo modo cuando x = L
dC A C A0 − C AL
J Ax x=L
= − DAe = − DAe (mL + b)−2 (24)
dx 1 1
x=L −
m(mL + b) mb

dC  R22 m(C A0 − C AL )
QAx = J Ax  R − DAe A
2
= − DAe (25)
x=L
(mL + b)2
x=L 2
dx x=L
(mL + b) −
b
La ecuacion (25) es la ecuacion que permite el cálculo del flujo molar de salida del cono, que es
el mismo al de ingreso.
Problema 3.13

Difusión de un gas a través de una pared cónica (I)

En un sistema que se encuentra a la presion atmosférica y 25oC están separados dos ambientes
por un tronco de cono poroso, en el cual el gas difunde axialmente . En el lado izquierdo del
cono divergente convergente se tiene una concentración de 0.15 mol/dm3 y su radio de 15 cm;
en el lado derecho, la concentración es de 0.03 mol/dm3, teniendo un radio de 5 cm. La longitud
del tronco de cono es 20 cm. Considere la difusividad del gas igual a 0.164 cm2/s

a) Obtenga el perfil de concentración

b) La densidad de flujo de materia en cada extremo del cono
c) El flujo molar en cada extremo del cono

Resultados

El programa desarrollado es el siguiente

d(P)/d(X) = 1
P(0) = 1
CA = (CAL - CA0) * ((m * L + b) / (m * X + b)) * (X / L) + CA0
Jax0 = (De * ((m * L + b) / (L * b)) * (CA0 - CAL))1000
JaxL =( -De * m * (CA0 - CAL) / ((m * L + b) - (m * L + b) ^ 2 / b))/1000
Qax0 = Jax0 * pi * R1 ^ 2
QaxL = JaxL * pi * R2 ^ 2
m = (R2 - R1) / L
b = 15
R1 = 15
R2 = 5
CAL = 0.03E-3
CA0 = 0.15E-3
De = 0.164
L = 20
pi = 3.1416
X(0) = 0
X(f) = 20

Resultado

a) El grafico correspondiente es el siguiente

b) Las densidades de flujo molar a la entrada y salida del cono expresado en mol/cm3s son
3.28E-10 y 2.952E-09, respectivamente.
 c) El flujo molar que resulta ser igual en ambos extremos del cono es 2.319E-07 mol/s
Problema 3.14

Transporte de masa en una pared cónica (II)

Usando el procedimiento de métodos numéricos encuentre los perfiles de concentración y

densidad de flujo molar, así como las densidades de flujo molar y los flujos molares a la entrada
y salida. Utilice los siguientes datos

CA0=0.15E-3 mol/cm3, CAL=0.03E-3 mol/cm3, R1=15 cm, R2=5 cm, L=20 cm, DAE=0.164cm2/s

Solución

La ecuacion de transporte de masa se obtiene a partir denla ecuacion (3.58)

d
( J Ax A( x)) = 0 (1)
dx
Derivando la ecuacion (1)

dJ Ax dA( x)
A( x) + J Ax =0 (2)
dx dx
Ahora trataremos de encontrar cada uno de los términos de (2). El área de la sección transversal
en cualquier punto x viene dado por:

A( x) =  (mx + b) 2 (3)

Derivando la ecuacion (3)

dA( x)
= 2 m(mx + b) (3)
dx
Reemplazando (2) y (3) en (2)

dJ Ax 2m
+ J Ax = 0 (4)
dx mx + b
La ley de Fick

dCA
J Ax = − DAe (5)
dx
Se resuelve simultáneamente (4) y (5)

Discretizando (34

J i − J i −1 2m
+ Ji = 0 (5)
x mix + b
Expandiendo (4)

2mx
i = 1 J1 − J 0 + J1 = 0
1 mx + b
 2mx
i = 2 J 2 − J1 + J2 = 0
2  mx + b
2mx
i = 3 J3 − J 2 + J3 = 0
3  mx + b
En este caso, debido a que la longitud de 20 cm es relativamente grande, para poder ajustar los
resultados de manera mas coherente con el valor analítico, se procede a tomar mas puntos de
discretización; en este caso discretizamos hasta i=40

2mx
i = 40 J 40 − J 39 + J 40 = 0
40  mx + b
Discretizando la ecuacion (5)

Ci − Ci −1
J i + DAe =0 (6)
x
Expandiendo la ecuacion (6)

C1 − C0
i = 1 J1 + DAe = 0 C0 = CA0
x
C2 − C1
i = 2 J 2 + DAe =0
x
C3 − C2
i = 3 J 3 + DAe =0
x
Y así sucesivamente

C40 − C39
i = 40 J 40 + DAe = 0 C40 = CAL
x
El programa respectivo es el siguiente
f(J0) = J1 - J0 + ((2 * m * DX) / (1 * m * DX + b)) * J1
f(J1) = J2 - J1 + ((2 * m * DX) / (2 * m * DX + b)) * J2
f(J2) = J3 - J2 + ((2 * m * DX) / (3 * m * DX + b)) * J3
f(J3) = J4 - J3 + ((2 * m * DX) / (4 * m * DX + b)) * J4
f(J4) = J5 - J4 + ((2 * m * DX) / (5 * m * DX + b)) * J5
f(J5) = J6 - J5 + ((2 * m * DX) / (6 * m * DX + b)) * J6
f(J6) = J7 - J6 + ((2 * m * DX) / (7 * m * DX + b)) * J7
f(J7) = J8 - J7 + ((2 * m * DX) / (8 * m * DX + b)) * J8
f(J8) = J9 - J8 + ((2 * m * DX) / (9 * m * DX + b)) * J9
f(J9) = J10 - J9 + ((2 * m * DX) / (10 * m * DX + b)) * J10
f(J10) = J11 - J10 + ((2 * m * DX) / (11 * m * DX + b)) * J11
f(J11) = J12 - J11 + ((2 * m * DX) / (12 * m * DX + b)) * J12
f(J12) = J13 - J12 + ((2 * m * DX) / (13 * m * DX + b)) * J13
f(J13) = J14 - J13 + ((2 * m * DX) / (14 * m * DX + b)) * J14
f(J14) = J15 - J14 + ((2 * m * DX) / (15 * m * DX + b)) * J15
f(J15) = J16 - J15 + ((2 * m * DX) / (16 * m * DX + b)) * J16
f(J16) = J17 - J16 + ((2 * m * DX) / (17 * m * DX + b)) * J17
f(J17) = J18 - J17 + ((2 * m * DX) / (18 * m * DX + b)) * J18
f(J18) = J19 - J18 + ((2 * m * DX) / (19 * m * DX + b)) * J19
f(J19) = J20 - J19 + ((2 * m * DX) / (20 * m * DX + b)) * J20
f(J20) = J21 - J20 + ((2 * m * DX) / (21 * m * DX + b)) * J21
f(J21) = J22 - J21 + ((2 * m * DX) / (22 * m * DX + b)) * J22
f(J22) = J23 - J22 + ((2 * m * DX) / (23 * m * DX + b)) * J23
f(J23) = J24 - J23 + ((2 * m * DX) / (24 * m * DX + b)) * J24
f(J24) = J25 - J24 + ((2 * m * DX) / (25 * m * DX + b)) * J25
f(J25) = J26 - J25 + ((2 * m * DX) / (26 * m * DX + b)) * J26
f(J26) = J27 - J26 + ((2 * m * DX) / (27 * m * DX + b)) * J27
f(J27) = J28 - J27 + ((2 * m * DX) / (28 * m * DX + b)) * J28
f(J28) = J29 - J28 + ((2 * m * DX) / (29 * m * DX + b)) * J29
f(J29) = J30 - J29 + ((2 * m * DX) / (30 * m * DX + b)) * J30
f(J30) = J31 - J30 + ((2 * m * DX) / (31 * m * DX + b)) * J31
f(J31) = J32 - J31 + ((2 * m * DX) / (32 * m * DX + b)) * J32
f(J32) = J33 - J32 + ((2 * m * DX) / (33 * m * DX + b)) * J33
f(J33) = J34 - J33 + ((2 * m * DX) / (34 * m * DX + b)) * J34
f(J34) = J35 - J34 + ((2 * m * DX) / (35 * m * DX + b)) * J35
f(J35) = J36 - J35 + ((2 * m * DX) / (36 * m * DX + b)) * J36
f(J36) = J37 - J36 + ((2 * m * DX) / (37 * m * DX + b)) * J37
f(J37) = J38 - J37 + ((2 * m * DX) / (38 * m * DX + b)) * J38
f(J38) = J39 - J38 + ((2 * m * DX) / (39 * m * DX + b)) * J39
f(J39) = J40 - J39 + ((2 * m * DX) / (40 * m * DX + b)) * J40
f(C1) = J1 + De * (C1 - C0) / DX
f(C2) = J2 + De * (C2 - C1) / DX
f(C3) = J3 + De * (C3 - C2) / DX
f(C4) = J4 + De * (C4 - C3) / DX
f(C5) = J5 + De * (C5 - C4) / DX
f(C6) = J6 + De * (C6 - C5) / DX
f(C7) = J7 + De * (C7 - C6) / DX
f(C8) = J8 + De * (C8 - C7) / DX
f(C9) = J9 + De * (C9 - C8) / DX
f(C10) = J10 + De * (C10 - C9) / DX
f(C11) = J11 + De * (C11 - C10) / DX
f(C12) = J12 + De * (C12 - C11) / DX
f(C13) = J13 + De * (C13 - C12) / DX
f(C14) = J14 + De * (C14 - C13) / DX
f(C15) = J15 + De * (C15 - C14) / DX
f(C16) = J16 + De * (C16 - C15) / DX
f(C17) = J17 + De * (C17 - C16) / DX
f(C18) = J18 + De * (C18 - C17) / DX
f(C19) = J19 + De * (C19 - C18) / DX
f(C20) = J20 + De * (C20 - C19) / DX
f(C21) = J21 + De * (C21 - C20) / DX
f(C22) = J22 + De * (C22 - C21) / DX
f(C23) = J23 + De * (C23 - C22) / DX
f(C24) = J24 + De * (C24 - C23) / DX
f(C25) = J25 + De * (C25 - C24) / DX
f(C26) = J26 + De * (C26 - C25) / DX
f(C27) = J27 + De * (C27 - C26) / DX
f(C28) = J28 + De * (C28 - C27) / DX
f(C29) = J29 + De * (C29 - C28) / DX
f(C30) = J30 + De * (C30 - C29) / DX
f(C31) = J31 + De * (C31 - C30) / DX
f(C32) = J32 + De * (C32 - C31) / DX
f(C33) = J33 + De * (C33 - C32) / DX
f(C34) = J34 + De * (C34 - C33) / DX
f(C35) = J35 + De * (C35 - C34) / DX
f(C36) = J36 + De * (C36 - C35) / DX
f(C37) = J37 + De * (C37 - C36) / DX
f(C38) = J38 + De * (C38 - C37) / DX
f(C39) = J39 + De * (C39 - C38) / DX

# J40 = -De * (C40 - C39) / DX

J40 = -De * (C40 - C38) / (2 * DX)
Q00 = pi * R1 ^ 2 * J00
Q40L = pi * R2 ^ 2 * J40L
R = Q00 / Q40L
m = (R2 - R1) / L
R2 = 10
R1 = 15
L = 20
De = 0.164
C0 = 0.15E-3
C40 = 0.03E-3
DX = L / 40
b = 15
pi = 3.1416
J00 = J0 / 1000
J40L = J40 / 1000

J0(0) = 0.00001
J1(0) = 0.00001
J2(0) = 0.00001
J3(0) = 0.00001
J4(0) = 0.00001
J5(0) = 0.00001
J6(0) = 0.00001
J7(0) = 0.00001
J8(0) = 0.00001
J9(0) = 0.00001
J10(0) = 0.00001
J11(0) = 0.00001
J12(0) = 0.00001
J13(0) = 0.00001
J14(0) = 0.00001
J15(0) = 0.00001
J16(0) = 0.00001
J17(0) = 0.00001
J18(0) = 0.00001
J19(0) = 0.00001
J20(0) = 0.00001

J21(0) = 0.00001
J22(0) = 0.00001
J23(0) = 0.00001
J24(0) = 0.00001
J25(0) = 0.00001
J26(0) = 0.00001
J27(0) = 0.00001
J28(0) = 0.00001
J29(0) = 0.00001
J30(0) = 0.00001
J31(0) = 0.00001
J32(0) = 0.00001
J33(0) = 0.00001
J34(0) = 0.00001
J35(0) = 0.00001
J36(0) = 0.00001
J37(0) = 0.00001
J38(0) = 0.00001
J39(0) = 0.00001

C1(0) = 0.001
C2(0) = 0.001
C3(0) = 0.001
C4(0) = 0.001
C5(0) = 0.001
C6(0) = 0.001
C7(0) = 0.001
C8(0) = 0.001
C9(0) = 0.001
C10(0) = 0.001
C11(0) = 0.001
C12(0) = 0.001
C13(0) = 0.001
C14(0) = 0.001
C15(0) = 0.001
C16(0) = 0.001
C17(0) = 0.001
C18(0) = 0.001
C19(0) = 0.001
C20(0) = 0.001
C21(0) = 0.001
C22(0) = 0.001
C23(0) = 0.001
C24(0) = 0.001
C25(0) = 0.001
C26(0) = 0.001
C27(0) = 0.001
C28(0) = 0.001
C29(0) = 0.001
C30(0) = 0.001
C31(0) = 0.001
C32(0) = 0.001
C33(0) = 0.001
C34(0) = 0.001
C35(0) = 0.001
C36(0) = 0.001
C37(0) = 0.001
C38(0) = 0.001
C39(0) = 0.001

Resultado

Hacer los gráficos de CA vs X y JA Vs X

Se ha desarrollado el programa considerando los mismos datos del problema 3.13; el resultado
final para la relación de los flujos molares a la entrada y salida resulta 0.9; es decir una
aproximación de 90%; el resultado final de dicha relación debería ser igual a la unidad. Se analizo
la sensibilidad y se encontró que cuando la diferencia de radios es pequeña, se mejora el valor
de la relación.

S ha ejecutado el programa considerando todos los datos del problema 3.13 excepto que se
cambio el radio final de 5 cm a 10 cm., con lo cual se obtuvo una relación de 0.987; es decir un
resultado mas próximo al verdadero, con lo cual, las densidades de flujo molar a la entrada y
salida expresada en mol/cm2s son 6.443732E-07 y 1.468401E-09, respectivamente. Los flujos
molares a la entrada y salida en el cono expresados en mol/s son 4.554817E-07 y 4.61313E-07,
respectivamente; estos dos valores deberían ser exactamente iguales. Estos errores podrían
atribuirse al método numérico dado que por definición son aproximados

3.9.- Deducción de las ecuaciones de difusión molecular y convectiva de

mezclas binarias y multicomponente
3.10.1 Mezclas binarias

Supongamos que se tiene un fluido estacionario (sin movimiento), entonces la velocidad de

difusión molecular de un componente A desde un punto de origen (referencial) hacia la derecha
será VAD , luego, la densidad de flujo molar estará dada por:

J A* = VADCA (3.66)
Consideremos ahora lo que sucede cuando la totalidad del fluido se mueve con un flujo
convectivo a la derecha y con una velocidad VM . El componente A sigue difundiéndose, pero
ahora su velocidad de difusión se mide respecto al fluido en movimiento.

Si observamos el sistema desde un punto de vista de un eje de referencia fijo, la velocidad de

A, V A , será la suma de las velocidades de difusión y la velocidad media de la fase total

VA = VAD + VM (3.67)

Esquematicamente, puede ser representada como:

Multiplicando la ecuacion (3.67) por C A se tiene

VAC A = VADC A + VM C A (3.68)

A partir de (3.68) se tiene

N A = J A* + VM CA (3.69)

Para un sistema de dos componentes la densidad de flujo molar total será:

N = VM (C A + CB ) = VM C = N A + N B (3.70)

Luego:

N A + NB
VM = (3.71)
C
Reemplazando (3.71) en (3.69)

N A + NB
N A = J A* + CA (3.72)
C
A partir de (3.72) se tiene

dyA
N A = −CDAB + yA ( N A + N B ) (3.73)
dz
La ecuacion (3.73) es la Ley de Fick para mezclas binarias cuando se produce simultáneamente
difusión y convección

Dentro de la difusión de mezclas binarias, se pueden distinguir claramente dos casos principales:

✓ La difusión de un fluido (gas o liquido) en otro fluido o su mezcla estancado (sin

movimiento), en condiciones isotérmicas o no isotérmicas
✓ La contradifusion (flujo de dos fluidos o solutos ) en direcciones opuestas y que pueden
tener una relación equimolar o no equimolar y en condiciones isotérmicas o no
 isotérmicas y a su vez con reacción instantánea o lenta cuando la superficie solida no
porosa es reactiva.

En ambos casos podemos enfocar los problemas de transferencia de masa en los diversos
sistemas de coordenadas considerando que los materiales que se transfieren son líquidos o
gases.

SISTEMA GASEOSO

Sistema gaseoso isotérmico

Transferencia de masa de un líquido volátil en una película de gas estancada

Consideremos el sistema que se muestra en la figura adjunta.

En el fondo del recipiente esta contenido el líquido volátil A , cuyo nivel se mantiene constante;
además se considera que el gas B es insoluble en el líquido B y la cantidad de líquido volátil en
la interfase liquido gas es proporcional a su presion de vapor a una determinada temperatura. En
la parte superior circula una mezcla de A y B, cuyas concentraciones en términos de la fracción
molar son conocidas.
Cuando se alcanza el estado estacionario, existe un movimiento ascendente neto de A, mientras
que B permanece fijo. En este caso se trata de encontrar la ecuacion de distribución de
concentración de A expresado en términos de la fracción molar y simultáneamente la densidad de
flujo molar de A y el flujo molar de A.
Balance de materia

N Az S z − N Az S z +z =0 (3.74)

Llevando al limite

N Az S − N Az S
z +z dN Az
lim z
= 0 (3.75)
z →0 S z dz
Ahora utilizamos la Ley de Fick que considera tanto el transporte difusivo como convectivo

dyA
N A = −CDAB + yA ( N A + NB ) (3.76)
dz
N B = 0 (Debido a que el gas B es estancado, es decir no fluye); de esta manera, la ecuacion
(3.77) se simplifica a:

CDAB dy A
NA = − (3.77)
1 − y A dz

Las ecuaciones (3.75) y (3.77) constituyen las bases para obtener el perfil de concentración del
flujo a través de un gas estancado.

Problema 3.15

Flujo axial a través de un gas estancado en un sistema en un tubo vertical

Encuentre el perfil de concentración y la densidad de flujo molar de la difusión de un gas A a

través de un gas estancado B en un tubo de sección constante

Solución

La ecuacion de conservacion de materia para un sistema no reaccionante y en la dirección axial

viene dada por la ecuacion (3.79)

dN Az
0 (1)
dz
La ley de Fick para el caso de flujo axial a través de un gas estancado viene dado por la ecuacion
(3.81)

CDAB dy A
NA = − (2)
1 − y A dz

Reemplazando (2) en (1)

d  CDAB dy A 
− =0 (3)
dz  1 − y A dz 

Integrando dos veces se obtiene

− ln(1 − y A ) = C1 z + C2 (4)

Condiciones de frontera

z = z1 y A = y A1 y z = z2 y A = y A2

Reemplazando las condiciones de frontera

1 1 − y A2
C1 = ln (5)
z1 − z2 1 − y A1

C2 = − ln(1 − y A1 ) − C1 z1 (6)

1 − yA z − z 1 − y A2
ln = 1 ln (7)
1 − y A1 z1 − z2 1 − y A1
La ecuacion (7) representa el perfil de composición molar a lo largo del tubo. Esta ecuacion
puede volver a escribirse del siguiente modo
z1 − z
1 − y A 2  z1 − z2
y A = 1 − (1 − y A1 )   (8)
 1 − y A1 
La densidad de flujo molar puede ser calculada a partir de la ecuacion (2), teniendo en
consideración que N Az z = z = N Az z = z , es decir es constante debido a que el área de
1 2

transferencia de masa también es constante.

CDAB dy A CDAB C1 (1 − y A1 )
N Az =− =− = −CDAB C1 (9)
z = z1
1 − y A1 dz z = z1 (1 − y A1 )

Ahora reemplazando el valor de la constante

CDAB 1 − y A2
N Az = ln (10)
z = z1
z2 − z1 1 − y A1

La ecuacion (10) puede ser escrita en términos del coeficiente de transferencia de masa,
considerando que: y A1 + yB1 = 1 y y A2 + yB2 = 1

CDAB  1 − y A2  ( y A1 − y A2 )
N Az = ln  (11)
z = z1
z2 − z1  1 − y A1  ( yB2 − yB1 )

CDAB ( y A1 − y A2 ) CDAB ( y A1 − y A2 )
N Az =  =  = km ( y A1 − y A2 ) (12)
z = z1
z2 − z1 ( yB2 − yB1 ) z2 − z1 ( yB )ln
ln( yB 2 / yB1 )

( yB )ln : Media logarítmica de composición molar

Otro modo de simplificar la ecuacion (10) es usando la aproximación: ln(1 − x)  − x . Con lo cual
se tiene:

CDAB 1 − y A2 CDAB
N Az = ln = ( y A1 − y A 2 ) (13)
z = z1
z2 − z1 1 − y A1 z2 − z1

Para que ambas ecuaciones (12) y (13) sean iguales, se requiere que ( yB )ln sea 1; sin embargo,
este valor fluctúa entre 0.8 y 0.9, con lo cual (12) resulta ser mas preciso.

Por otro lado, la ecuacion (12) también puede ser escrita en términos de las presiones parciales

PDAB P( y A1 − y A2 ) PDAB ( p A − pA2 )

N Az =  =  1 = km ( pA1 − p A2 ) (14)
z = z1
RT ( z2 − z1 ) P( yB2 − yB1 ) RT ( z2 − z1 ) ( pB )ln
ln( PyB 2 / PyB1 )

El flujo molar, se obtiene multiplicando la densidad de flujo molar por el área de transferencia
de masa (Area normal a la densidad de flujo que seria el área de la sección del tubo)
Problema 3.16

Flujo de un líquido volátil en un gas estancado

Considerando las ecuaciones anteriores, obtenga el perfil de composición molar y la densidad

de flujo molar dado en el esquema de la figura dada anteriormente.

Datos: Ojo los datos tipear en mathtype

Longitud del tubo desde la superficie libre de la acetona: 10 cm

Coeficiente de difusión de la acetona en el aire: 0.164 cm2/s

Temperatura del medio: 25 oC

Presion de vapor de la acetona: 231.06 mmHg

Presion total 760 mmHg

Presion del aire en la parte superior del tubo: 720 mmHg.

Solución:

Se desarrolla el programa respectivo tomando en cuenta las ecuaciones deducidas en el

problema 3.15.

El programa desarrollado es el siguiente

YA = 1 - (1 - YA1) * ((1 - YA2) / (1 - YA1)) ^ E
YB = 1 - YA
E = (Z1 - Z) / (Z1 - Z2)
PVAP = 231.06
P2 = 720
PT = 760
YA1 = PVAP / PT
YA2 = (PT - P2) / PT
Z1 = 0
Z2 = 10
NAZ = ((C * DAB) / (Z2 - Z1)) * (YA1 - YA2) / MLY
MLY = (YB2 - YB1) / ln(YB2 / YB1)
C = (PT/(R * T)) / 1000
R = 62.4
T = 298
DAB = 0.164
YB1 = 1 - YA1
YB2 = 1 - YA2
d(P)/d(Z) = 1
P(0) = 0
Z(0) = 0
Z(f) = 10

Resultado. El perfil de composición molar de los componentes A y B son mostrados en la Figura

xx
 −7
La densidad de flujo molar es N AZ = 2.067 10 mol/cm2 s

Problema 3.17
Evaporación de metanol en un tubo de Stefan tal como se muestra en la Figura 15
La evaporación de metanol en una corriente de aire seco se estudia en un tubo de Stefan con una
longitud de 0.238 m; las medidas se hacen a 328.5 K, donde la presion de vapor de metanol es de
68.4 Kpa y la presion total es de 99.4 Kpa. El coeficiente de difusión del metanol en el aire es de
1.991x10-5 m2/s, cuyo valor fue estimado por Taylor y Krishna.
a) Calcule la densidad de flujo dentro del tubo de Stefan usando una técnica numérica,
considerando que a la salida del tubo la fracción molar de metanol es igual a cero.
b) Reporte gráficamente la fracción molar de metanol desde la superficie liquida de metanol
hasta la corriente de aire que fluye
c) Compare sus valores con el obtenido por la ecuacion (12) del problema 3.15

Solución
La ecuacion usada para este caso, es la misma del problema anterior, es decir:
Ecuacion de balance de materia

dN Az
0 (1)
dz
Ley de Fick para difusión a través de una pelicua estancada

CDAB dy A
NA = − (2)
1 − y A dz

A partir de la ecuacion (2) se obtiene

dy A N (1 − y A )
=− A (3)
dz CDAB

Utilizaremos la técnica denominada método del disparo; las condiciones de frontera que se
tienen y las requeridas son:

pA1 68.4
z = 0 N A ?? y A = = = 0.688
P 99.4
Para resolver el problema, se requiere conocer la densidad de flujo molar en la interfase liquido
gas; sin embargo no se conoce; por tanto, se procederá de manera iterativa, es decir se asume
dicho valor hasta cumplir la condición cuando z = 0.238m y A = 0 .

A partir de (15) se tiene

CDAB ( y A1 − y A2 )
N Az =  (4)
z = z1
z2 − z1 ( yB )ln

El programa desarrollado es el siguiente

d(YA)/d(z) = -(1 - YA) * NA / (DAB * C)
d(NA)/d(z) = 0
z(0) = 0
z(f) = 0.238
YA(0) = 0.688
DAB = 1.991E-5
C = P / (R * T)
R = 8.31434
T = 328.5
P = 99.4
z1 = 0
z2 = 0.238
YB = 1 - YA
# Este valor se obtiene por tanteo es decir por el método del disparo
NA(0) = 3.5461E-6
# Se considera valor correcto cuando en Z=0.238 se 0btiene YA=0
# Luego se obtiene la solución analítica
NA_CALC = DAB * C * (YA1 - YA2) / ((z2 - z1) * YBLN)
YBLN = (YA1 - 0) / ln((1 - YA2) / (1 - YA1))
YA1 = 0.688 # PA1 / P
PA1 = 68.4
YA2 = 0

Resultado
−7
a) Usando la técnica numérica por el método del disparo, se obtiene N AZ = 3.546 10
mol/m2s
b) La grafica correspondiente es:
 −7
c) La ecuacion (12) corresponde al método analítico y permite obtener N AZ = 3.546 10
mol/m2s, cuyo valor es el mismo que el reportado por el método del disparo.

Problema 3.17
Coeficiente de difusión de clorobenceno en el aire
Para determinar el coeficiente de difusión del clorobenceno en el aire, se evapora el mismo
isotérmicamente a 332K en un tubo abierto de 3 mm de diámetro parcialmente lleno de este
líquido; durante el experimento se obtuvieron los siguientes datos
Nivel de clorobenceno en el tubo antes del experimento: 6.5 cm
Nivel de clorobenceno en el tubo después del experimento: 7.6 cm
Duración del experimento 260000 s.
Determine el coeficiente de difusión del clorobenceno en el aire.

Datos adicionales
Densidad del clorobenceno: 1100 Kg/m3
Presion total: 101 330 N/m2
Presion de vapor del benceno a 332 K: 8373 N/m 2
Solución
La densidad de flujo molar del benceno se puede obtener multiplicando la variación del numero
de moles respecto al tiempo y por unidad de área, es decir:

1 dnA
N Az = (1)
S dt
Con las relaciones auxiliares: nA = mA / M A , mA =  Av y v = S  z

v : Volumen del liquido

S : Área de la sección transversal del capilar
z : Altura del tubo capilar
Los demás términos tienen la nomenclatura clásica conocida.

Reemplazando (1) en términos de las ecuaciones auxiliares

1 1 dmA 1 1 d (  Av)
N Az = = (2)
S M A dt S M A dt

Ahora la ecuacion (22) se puede escribir en términos de la variación de altura

1 S d (  A z )  A dz
N Az = = (3)
S M A dt M A dt

La densidad de flujo molar también puede ser obtenida a partir de la ecuacion (14) del problema
3.15.

PDAB ( p A1 − p A2 )
N Az z = z1
=  (4)
RTz ( pB )ln

Igualando (23) y (24)

 A dz PDAB ( p A1 − p A2 )
=  (5)
M A dt RTz ( pB )ln

Reordenando

PDAB M A ( p A1 − p A2 )
L2 t

 zdz =
L1
RT

 A ( pB )ln 0
dt (6)

Al integrar la ecuacion (7) se tiene

L22 − L12 PDAB M A ( p A1 − p A2 )

=  t (7)
2 RT  A ( pB )ln

De la ecuacion (7) se debe obtener el coeficiente de difusión

0.0762 − 0.06522 101330  DAB 112.5 (8373 − 0)

=   26 104
2 8314  332 1100 101330 − 92957
 −6
Operando se obtiene: DAB = 9.2110 m / s
2

Transferencia de masa de un líquido volátil en una película de gas estancada

El problema es similar al ítem anterior (Transferencia de masa radial en un sistema cilíndrico de

un gas A a través de otro B estancado)

Balance de materia

N Ar r 2 rL − N Ar r +r
2 rL = 0 (3.78)

Llevando al limite

N Ar Ar +r
2 rL − N Ar 2 rL
lim r
=0 (3.79)
r →0 2 Lr
Simplificando la ecuacion (3.79) se tiene

d (rN Ar )
=0 (3.80)
dr
La ley de Fick esta dada por:

dy A
N Ar = −CDAe + y A ( N Ar + N Br ) (3.81)
dr
N Br = 0 (Por tratarse de un gas estancado). La ecuacion (3.85) queda del siguiente modo:

CDAe dy A
N Ar = − (3.82)
1 − y A dr

Las ecuaciones (3.80) y (3.82) son la base para obtener el perfil de composición molar, así como
la densidad de flujo molar en cualquiera de los límites del sistema.

Problema 3.18
Difusión de amoniaco a través de un gas estancado contenido en una pared cilíndrica
Amoniaco gaseoso se difunde a través de una capa estacionaria de aire a 298K y 1 atm de presion
entre dos cilindros de paredes porosas, cuya resistencia a la difusión del gas es despreciable. En
el interior del cilindro más pequeño se mantiene una presion de 76 mmHg y en la parte exterior
del cilindro más grande, la presion es nula. Obtenga la ecuacion de distribución de presion en la
capa estacionaria de aire así como el flujo molar por unidad de longitud de tubo.
Datos
Coeficiente de difusión de amoniaco en aire a 298K y 1 atm: 0.217 cm2/s

Radios de los cilindros: R1 = 0.5 cm y R2 = 2.5 cm

Solución

La ecuacion de transferencia de masa para este caso viene dada por la ecuacion (3.84)

d (rN Ar )
=0 (1)
dr
La ley de Fick para el flujo a través de un gas estancado viene dado por la ecuacion (3.86)

CDAe dy A
N Ar = − (2)
1 − y A dr

Reemplazando (2) en (1)

d  CDAe r dy A 
− =0 (3)
dr  1 − y A dr 

Integrando (3)

r dy A
= C1 (4)
1 − y A dr

Volviendo a integrar (4)

− ln(1 − y A ) = C1 ln r + C2 (5)

Condiciones de frontera

r = 0.005m yA1 = ( pA1 / P) = 76 / 760 = 0.1 r = 0.025m yA2 = ( pA2 / P) = 0 / 760 = 0

1 1 − y A2 1 1− 0
C1 = ln = ln = −0.06546 (6)
ln( R1 / R2 ) 1 − y A1 ln(0.005 / 0.025) 1 − 0.1
C2 = − ln(1 − yA1 ) − C1 ln R1 = − ln(1 − 0.1) − (−0.06546)  ln(0.005) = −0.2415 (7)

Luego, estas constantes calculadas se sustituyen en (7)

ln(1 − y A ) = −0.06546ln r + 0.2415 (8)

La ecuacion (8) es equivalente a:

yA = 1 − 1.273 r 0.6456 (9)

La ecuacion (9) puede ser escrito en términos de la presion parcial

pA = yA P = P(1 − 1.273 r 0.6456 ) (10)

Dado que la presion parcial es 760 mmHg, las ecuaciones (9) y (10) pueden ser representadas
en un gráfico.

La densidad de flujo molar depende del radio, por tanto: N Ar r = R1

 N Ar r = R12
; sin embargo, el
flujo molar es constante tal como señala la ecuacion de conservación de materia.

La densidad de flujo molar puede ser calculada a partir de las ecuacion (2) y (4)

CDAe dy A CDAe C1 C
N Ar =− =− (1 − y A1 ) = −CDAe 1
r = R1
1 − y A1 dr r = R1 (1 − y A1 ) R1 R1

1 CDAe 1 − y A2 P DAe 1 − y A2
N Ar = ln = ln
r = R1
R1 ln( R2 / R1 ) 1 − y A1 RTR1 ln( R2 / R1 ) 1 − y A1

Reemplazando datos

P DAe 1 − y A2 760 10−3 0.217 1− 0

N Ar = ln = ln = 1.16110−6 mol / cm 2 s
r = R1
RTR1 ln( R2 / R1 ) 1 − y A1 62.4  298  0.5 ln(2.5 / 0.5) 1 − 0.1

El flujo molar se obtiene mediante:

QA / L = N Ar r = R1
2 R1 = 2 1.16110−6 mol / cm2 s  0.5cm = 3.648 10−6 mol / cm s

El resultado grafico se muestra en la Figura P_3.18

Figura P_3.18: Perfil de composición molar del flujo de un gas A en otro estancado en una pared
cilíndrica.
Sistema gaseoso no isotérmico

Película gaseosa de forma esférica

Se analiza la transferencia de masa de un gas A a través de una película de gas B, estancado,
donde dicha película esférica tiene una temperatura variable cuya forma es la siguiente:

(T / T1 ) = (r / R1 )n (3.83)

Como consecuencia de ello, también se puede considerar que el coeficiente de difusión es

variable, cuyas formas provienen de expresiones empíricas, por tanto pueden tener diversas
expresiones de la difusividad, tal como la mostrada

( DAB / DAB1 ) = (T / T1 )1/2 (3.84)

donde T1 la temperatura de la superficie interna y DAB1 es el coeficiente de difusión a dicha

temperatura

Figura 3.15: Difusión molecular no isotérmica en una película esferica estancada

El balance de materia para el sistema de coordenadas esféricas está dada por:

N Ar r 4 r 2 − N Ar r +r
4 r 2 = 0 (3.85)

Llevando al limite
 Nr Ar +r
4 r 2 − N Ar 4 r 2
lim r
=0 (3.86)
r →0 4r
Simplificando la ecuacion (3.86) se tiene

d (r 2 N Ar )
=0 (3.87)
dr
La ecuacion (3.87) es la ecuacion de transferencia de masa en régimen estacionaro, sin generacion
La ley de Fick esta dada por:

dyA
N Ar = −CDAB + y A ( N Ar + N Br ) (3.88)
dr
N Br = 0 (Gas estancado)

CDAB dy A
N Ar = − (3.89)
1 − y A dr

Problema 3.19

Transferencia de masa en una película estancada de forma esferica

En la parte externa de una esfera porosa (cuya resistencia a la difusión) tal como que se muestra
en la figura xxxx hay un gas estancabo B, a través del cual un gas A difunde en la dirección
radial. Obtenga la ecuacion de distribución de composición molar, así como la densidad de flujo
molar y el flujo molar, considerando que la temperatura varia en la película en función al radio
de acuerdo con la ecuacion: (T / T1 ) = (r / R1 )
n
y el coeficiente de difusión varia con la
temperatura de acuerdo con: ( DAB / DAB1 ) = (T / T1 ) ,
1/2

Solución
Se utiliza la ecuacion de transferencia de masa dada en (3.87)

d (r 2 N Ar )
=0 (1)
dr
La ley de Fick dado en (3.89)

CDAB dy A
N Ar = − (2)
1 − y A dr

Reemplazando (2) en (1)

d  2 CDAB dy A 
 −r =0 (6)
dr  1 − y A dr 

Por otro lado de las relaciones de coeficiente de difusión y temperatura

La relación entre la temperatura y el radio viene dada por

(T / T1 ) = (r / R1 )n (7)

La relación entre coeficiente de difusión y temperatura

( DAB / DAB1 ) = (T / T1 )1/2 (8)

Relacionando (7) y (8) se obtiene

( DAB / DAB1 ) = (r / R1 ) n /2 (9)

Reemplazando (7) en (6)

d  2 CDAB1 dy 
 −r (r / R1 ) n / 2 A  = 0 (10)
dr  1 − yA dr 

Reemplazando (9) en (6) y eliminando todas las constantes de (10) se obtiene

d  r 2+ n /2 dy A 
 =0 (11)
dr  1 − y A dr 

Integrando la ecuacion (11)

r 2+ n /2 dy A
= C1 (10)
1 − y A dr

Integrando nuevamente

dy A dr
 1− y A
= C1  
r 2 + n /2
(12)

La ecuacion (12) permite obtener

r −1−n /2
− ln(1 − y A ) = −C1 + C2 (13)
1+ n / 2
Condiciones de frontera

r = R1 y A = y A1 y r = R2 y A = y A2

(1 + n / 2) (1 − y A2 )
C1 = −1− n /2
−1− n /2
ln (14)
R 2 − R1 (1 − y A1 )

C1R1−1−n /2
C2 = − ln(1 − y A1 ) + (15)
1+ n / 2
Reemplazando (15) en (13) y reordeando

(1 − y A )
ln =
C1
(1 − y A1 ) 1 + n / 2
( r −1− n /2 − R1−1− n /2 ) (16)
Abora remplazamos el valor de C1 en (16), ademas hacemos el cambio a = −1 − n / 2 con
fines de simplificar la escritura.

(1 − y A ) ( r − R1 ) (1 − y A2 )
a a

ln = a ln (17)
(1 − y A1 ) ( R2 − R1a ) (1 − y A1 )

La ecuacion (17) se puede escribir del siguiente modo

p
 (1 − y A2 ) 
y A = 1 − (1 − y A1 )   (18)
 (1 − y A1 ) 

p=
(r a
− R1a )
(19)
( R2a − R1a )

La densidad de flujo molar se obtiene mediante

CDAB1 dy A
N Ar =− (20)
r = R1
1 − y A1 dr r = R1

Nota: Se debe tener en cuenta que las densidades de flujo molar no son iguales en cada punto
referido al radio; en consecuencia: N Ar r = R1
 N Ar r =R2

La derivada se calcula a partir de la ecuacion (10)

dy A (1 − y A1 )
= C1 2+ n /2 (21)
dr r = R1 R1

Reemplazando (21) en (20)

CDAB1 (1 − y A1 ) CDAB (1 + n / 2) (1 − y A2 )
N Ar =− C1 = − 2+n/21 −1−n /2 −1− n /2
ln (22)
(1 − y A1 ) R1
2+n/2
r = R1
R1 R2 − R1 (1 − y A1 )

Ahora (22) puede ser escrita en términos de la media logarítmica de composición molar

CDAB1 (1 + n / 2) (1 − y A1 ) ( y A1 − y A2 )
N Ar = 2+ n/2 −1− n /2 −1− n /2
ln  (23)
r = R1
R1 R2 − R1 (1 − y A2 ) ( yB 2 − yB1 )

A partir de (23) se puede obtener la expresión

CDAB1 (1 + n / 2) ( y − y A2 )
N Ar =−  A1 (24)
r = R1
R1
2+ n/2
R −1− n /2
2 − R1−1− n /2
( yB 2 − yB1 )
ln( yB 2 / yB1 )

Finalmente, se obtiene

CDAB1 (1 + n / 2) ( y − y A2 )
N Ar =− 2+ n /2 −1− n /2
−1− n / 2
 A1 (25)
r = R1
R1 R 2 − R1 ( yB )ln

El flujo molar, que es constante en cada punto R1 y R2 se obtiene mediante:

 4 CDAB1 (1 + n / 2) ( y − y A2 )
QAr = QAr = 4 R12 N Ar =− −1− n /2 −1− n /2
 A1 (26)
r = R1 r = R2 r = R1
− R1
n /2
R1 R2 ( yB )ln

Problema 3.20
Difusión no isotérmica de un gas a través de una película gaseosa estancada, contenido en un
sistema esférico
Considerando los siguientes datos y las ecuaciones deducidas anteriormente, obtenga:
a) Los perfiles de composición molar de A y B
b) La variación del coeficiente de difusión con el radio
c) Variación de la temperatura con el radio
d) La densidad de flujo molar en la superficie interna (correspondiente al radio más
pequeño)
e) El flujo molar

Datos: P = 760mmHg , T1 = 298K , DAB1 = 0.252cm / s , R1 = 20cm , R2 = 22cm ,

2

y A1 = 0.2 , y A 2 = 0.08
Solución

El programa desarrollado es el siguiente

YA = 1 - (1 - YA1) * ((1 - YA2) / (1 - YA1)) ^ b # Ecuacion (18)

YB = 1 - YA
b = (r ^ a - R1 ^ a) / (R2 ^ a - R1 ^ a)
C = (P/(R * T1)) / 1000
P = 760
R = 62.4
T1 = 298
DAB1 = 0.252
n = 0.5
R1 = 20
R2 = 22
YA1 = 0.2
YA2 = 0.08
YB2 = 1 - YA2
YB1 = 1 - YA1
MLYB = (YB2 - YB1) / ln(YB2 / YB1)
NAR1 = -C * DAB1 * (1 + n / 2) * (YA1 - YA2) / (R1 ^ d * D * MLYB)
Q1 = NAR1 * 4 * 3.1416 * R1 ^ 2
Q2 = -4 * 3.1416 * C * DAB1 * (1 + n / 2) * (YA1 - YA2) / (R1 ^ (n / 2) * D * MLYB) # Ecuacion (26)
D = R1 ^ (n / 2) * (R2 ^ a - R1 ^ a)
a = -1 - n / 2
d=2+n/2
DAB = DAB1 * (r / R1) ^ (n / 2)
T = T1 * (r / R1) ^ n
r(0) = 20
r(f) = 22
d(S)/d(r) = 0
S(0) = 0

Result ado

a) Se muestran los perfiles de composición molar de A y B.

Figura P_3.20A Perfiles de composición molar de A y B en función al radio de la película esferica

 b) Curva de variación del coeficiente de difusión

Figura P_3.20B : Variación del coeficiente de difusion de A en funcion al radio de la película esferica

c) Perfil de temperatura en función al radio de la película esférica.

 d) La densidad de flujo molar en la superficie interna es: 4.087 10−5 mol / cm2 s
e) El flujo molar es 1.904 10−3 mol / s

Película gaseosa de forma cilíndrica

Para deducir las ecuaciones del proceso de transferencia de masa en una película estancada no
isotérmica de forma cilíndrica, se debe aplicar el balance de materia, juntamente con la Ley de
Fick y la ecuacion que relaciona la difusividad con la temperatura.
El balance de materia para el sistema de coordenadas cilíndricas está dada por:

N Ar r 2 rL − N Ar r +r
2 rL = 0 (3.90)

Llevando al limite

Nr Ar +r
2 rL − N Ar 2 rL
lim r
=0 (3.91)
r →0 2 Lr
Simplificando la ecuacion (3.91) se tiene

d (rN Ar )
=0 (3.92)
dr
La ley de Fick esta dada por:

dyA
N Ar = −CDAB + y A ( N Ar + N Br ) (3.93)
dr
N Br = 0 (Gas estancado)

CDAB dy A
N Ar = − (3.94)
1 − y A dr

Problema 3.21
Difusión no isotérmica de un gas a través de una película estancada no isotérmica en un sistema
cilíndrico
Consideremos que en el espacio anular comprendido entre dos cilindros verticales de paredes
porosas (cuya resistencia al la transferencia de masa es despreciable), hay un gas estancado (aire).
En el interior del cilindro mas pequeño se encuentra un gas que se difunde hacia el exterior. Las
propiedades físicas, la geometría y las condiciones de operación son conocidas.
Encuentre la ecuacion de distribución de fracción molar, así como la distribución de temperatura
y de la densidad de flujo molar, conociendo que el coeficiente de difusión varia con la temperatura
según la relación: ( DAB / DABA ) = (T / TA )
1/2
Solución
Este problema es relativamente complejo debido al numero de variables involucradas (fracción
molar, temperatura, coeficiente de difusión, radio). Comparado con el problema anterior, ahora
no se conoce la relación entre la temperatura y el coeficiente de difusión vs el radio, por lo que
tendrá que efectuarse además del balance de materia, el balance de energía respectivo.
La ecuacion de balance de masa para este caso viene dada por (3.92)

d (rN Ar )
=0 (1)
dr
La ley de Fick para este caso esta dado por (3.94)

CDAB dy A
N Ar = − (2)
1 − y A dr

Reemplazando (2) en (1)

d  CDAB dy A 
 −r =0 (3)
dr  1 − y A dr 

Por otro lado de las relaciones de coeficiente de difusión y temperatura

DAB = DABA (T / TA )1/2 (4)

Reemplazando (4) en (3)

d  CDAB0 (T / T0 ) dy A 
1/2

 −r =0 (5)
dr  1 − yA dr 

En la ecuacion (5) se puede simplificar todas las constantes y obtener

d  r T 1/2 dy A 
 =0 (6)
dr  1 − y A dr 

La ecuacion (6) es una ecuacion diferencial ordinaria con dos variables dependientes; para ser
resuelta, se requiere de la ecuacion de conservación de energía
En esta etapa se puede proceder de dos modos
a) Realizar el balance de energía, obtener el perfil de temperatura (en función del radio) y
reemplazar en (6) previamente habiendo obtenido sus derivadas.
Derivando (6), para lo cual se puede separar en dos factores

d  r 1/2 dy A   r   1/2 dy A   r   1/2 dy A 

 T =
   T − T  =0 (7)
dr  1 − y A dr   1 − y A   dr   1 − y A   dr 

Ahora se deriva cada una de sus partes

r (1 − y A ) − (1 − y A )r 1/2 dy A  r   1/2  dy A 1/2  dy A 


T +   ( T ) dr + T   =0 (8)
(1 − y A ) 2 dr  1 − y A    dr  
 
Aplicando las derivadas respectivas

 1 r dy A  1/2 dy A  r   1 −1/2 dT dy A 1/2 d y A 

2

 +  +   + =0
dr  1 − y A   2 dr 2 
T T T (9)
 (1 − y A ) (1 − y A ) dr 
2
dr dr

Multiplicando por (1 − y A )

 r dy A  1/2 dy A 1 dT dy A d 2 yA 
1 + T + r  T −1/2 + T 1/2 =0
dr 2 
(10)
 (1 − y A ) dr  dr 2 dr dr

Dividiendo por T 1/ 2

 r dy A  dy A  1 −1 dT dy A d 2 y A 
1 +  + r  2 T dr dr + dr 2  = 0 (11)
 (1 − y A ) dr  dr  

Dividiendo por r

1 1 dy A  dy A  1 −1 dT dy A d 2 y A 
 +  + T + =0
dr 2 
(12)
 r (1 − y A ) dr  dr  2 dr dr

Ahora reordenando
2
d 2 y A 1 dy A 1  dy A  1 dT dy A
+ +   + =0 (13)
dr 2
r dr (1 − y A )  dr  2T dr dr

Ahora, se requiere hacer uso del balance de energía para obtener el valor de la temperatura

d (rqr )
=0 (14)
dr
Integrando la ecuacion (14) se obtiene

T = C1 ln r + C2 (15)

Las condiciones de frontera son:

r = RA T = TA y r = RB T = TB

Usando estas condiciones se obtiene:

(TA − TB ) (TA − TB )
C1 = , C2 = TA − C1 ln RA = TA − ln RA
ln( RA / RB ) ln( RA / RB )

Reemplazando en (15) se obtiene la ecuacion de distribución de temperatura

T = C1 ln r + TA − C1 ln RA = C1 ln(r / RA ) + TA (16)

Finalmente se tiene

(TA − TB )
T= ln(r / RA ) + TA (17)
ln( RA / RB )

La ecuacion (15) y su derivada deben ser reemplazadas en (13)

2
d 2 y A 1 dy A 1  dy A  C1 dy A
+ +   + = 0 (18)
dr 2
r dr (1 − y A )  dr  2r (C1 ln r + C2 ) dr

Ahora discretizamos la ecuacion (18) omitiendo el subíndice A

yi +1 − 2 yi + yi −1 1 yi − yi −1 1  yi − yi −1  yi − yi −1
2
C1
+ + + = 0 (19)
r 2 ir r (1 − yi )  r  2ir (C1 ln ir + C2 ) r

Reordenando

1 ( y − yi −1 ) 2 C1 ( yi − yi −1 )
yi +1 − 2 yi + yi −1 + ( yi − yi −1 ) + i + =0 (20)
i (1 − yi ) 2i (C1 ln ir + C2 )

Expandiendo la ecuacion (20)

1 ( y − y0 ) 2 C1 ( y1 − y0 )
i = 1 y2 − 2 y1 + y0 + ( y1 − y0 ) + 1 + =0
1 (1 − y1 ) 2 1 (C1 ln(1 r ) + C2 )

r = RA y0 = y A (Condición de frontera)

1 ( y − y1 ) 2 C1 ( y2 − y1 )
i=2 y3 − 2 y2 + y1 + ( y2 − y1 ) + 2 + =0
2 (1 − y2 ) 2  2  (C1 ln(2  r ) + C2 )

1 ( y − y2 ) 2 C1 ( y3 − y2 )
i=3 y4 − 2 y3 + y2 + ( y3 − y2 ) + 3 + =0
3 (1 − y3 ) 2  3  (C1 ln(3  r ) + C2 )

Y así sucesivamente

1 ( y − y8 ) 2 C1 ( y9 − y8 )
i=9 y10 − 2 y9 + y8 + ( y9 − y8 ) + 9 + =0
9 (1 − y9 ) 2  9  (C1 ln(9  r ) + C2 )

r = RB y10 = yB (Condición de frontera)

 b) Resolver en forma simultánea el balance de materia y el perfil de temperatura obtenido a
partir de la integración del balance de energía

A partir de la ecuacion (13) se tiene

2
d 2 y A 1 dy A 1  dy A  1 dT dy A
+ +   + =0 (21)
dr 2
r dr (1 − y A )  dr  2T dr dr

Discretizando (21)

yi +1 − 2 yi + yi −1 1 yi − yi −1 1  yi − yi −1  (Ti − Ti −1 ) yi − yi −1
2

+ + + = 0 (22)
r 2
ir r (1 − yi )  r  2Ti r

Reordenando

1 ( y − yi −1 ) 2 r (Ti − Ti −1 )
yi +1 − 2 yi + yi −1 + ( yi − yi −1 ) + i + ( yi − yi −1 ) = 0 (23)
i (1 − yi ) 2Ti

Se expande (23) con las mismas condiciones dadas en el item (a), considerando T0 = TA y
T10 = TB

Del Balance de energía

T = C1 ln r + C2 (24)

La temperatura se puede expresar para cada punto

Ti = C1 ln(ir + RA ) + C2 (25)

Se expande (25) de i=1…9.