“Universidad Autónoma del Estado de México”

Facultad de Química

Ingeniería Química
Transferencia de calor
“Serie de problemas 2”
Nombre del profesor(a): Cesar Perez Alonso

Integrantes:

 Gómez García Ricardo

 Iturbe Tabaco Víctor Hugo
 Pérez Ramírez José Manuel
 López Thompson Lorrie Anne
 Torres Llerena Aldo Enrique
No. equipo: 4 Fecha de entrega:
Grupo: 65 04/04/2024

Sexto Semestre
 ESCALA
CRITERIOS
COMPETENTE (10.0 – 8.0) SUFICIENTE (7.9 – 6.0) INSUFICIENTE (5.9 – 0.0)
El proyecto delimita el tema, así El proyecto delimita el tema, así como El proyecto delimita el tema, pero
Identifica y define
como plantea el problema de plantea el problema de forma clara y no plantea el problema de forma
el problema de
forma clara y precisa y lo define precisa, pero no lo define en términos clara y precisa y ni lo define en
diseño preliminar
en términos objetivos. objetivos. términos objetivos.
Realiza una recopilación, análisis No realiza una recopilación,
Realiza parcialmente una recopilación,
y evaluación de la información análisis y evaluación de la
Recopila análisis y evaluación de la información
disponible para el diseño información disponible para el
información disponible para el diseño preliminar
preliminar del equipo de diseño preliminar del equipo de
necesaria del equipo de transferencia de masa.
transferencia de masa. Realiza transferencia de masa. No realiza
Realiza una investigación bibliográfica
una investigación bibliográfica una investigación bibliográfica
Establece correctamente las No establece correctamente las
Establece correctamente las
necesidades, requerimientos y necesidades, requerimientos y
Busca soluciones necesidades, requerimientos y
especificaciones del equipo a especificaciones del equipo a
creativas especificaciones del equipo a diseñar,
diseñar y busca soluciones diseñar, ni busca soluciones
pero no busca soluciones creativas.
creativas. creativas.
Propone al menos dos soluciones
realizables que cumplen con las
especificaciones de diseño del Propone sólo una solución realizable Propone una solución irrealizable
Propone y evalúa equipo considerando que cumple con las especificaciones y que no satisface las
alternativas requerimientos y restricciones. de diseño del equipo considerando especificaciones de diseño del
Define claramente uno o varios requerimientos y restricciones equipo.
criterios para seleccionar la
solución más apropiada.
Describe de manera correcta y
Describe sólo las principales Describe de manera imprecisa las
Describe el diseño completa las especificaciones
especificaciones técnicas y de especificaciones técnicas y de
preliminar del técnicas y de funcionamiento del
funcionamiento del equipo, y lo funcionamiento del equipo y lo
equipo de equipo, y lo representa mediante
representa mediante planos, representa mediante dibujos
transferencia de planos, diagramas y dibujos
diagramas y dibujos utilizando una incompletos sin respetar la
calor utilizando una simbología
simbología apropiada. simbología apropiada.
apropiada.
Presenta un documento No presenta un documento
Presenta un documento ordenado con
ordenado y completo, con toda la ordenado y completo, con toda la
Integra un información parcial utilizada para el
información de apoyo utilizada información de apoyo utilizada
documento del diseño del equipo, y cita
para el diseño del equipo, y cita para el diseño del equipo, ni cita
diseño preliminar correctamente las fuentes de consulta
correctamente las fuentes de correctamente las fuentes de
de calidad reconocida.
consulta de calidad reconocida. consulta de calidad reconocida.
26.10

En una cámara de combustión caliente se difunde oxígeno a través de aire hasta la

superficie del carbón en donde reacciona para formar CO y/o CO 2. La fracción mol de
oxígeno a z = 0 es 0.21. La reacción en la superficie es instantánea. En la película
gaseosa no se produce ninguna reacción.

Determinar la rapidez de difusión del oxígeno por hora a través de un metro cuadrado
de área si

a) Solo se produce monóxido de carbono en la superficie del carbón;

2 C+O2 →2 CO
Solución

Esquema de interpretación del sistema:

Por conveniencia el punto de referencia z = 0 se cambió a la superficie de la película

estancada.

z=0 yA = yA0

z=δ yA = yAδ

Consideraciones para resolver el sistema

  El sistema opera en estado estacionario
 La concentración yA0 > yAδ
 La temperatura y presión permanecen constantes por lo que el coeficiente de
difusividad DAB permanece constante
 El sistema está concentrado, por lo tanto, el mecanismo de transporte que
predomina en el sistema es convectivo.
 El flux molar es unidireccional en la coordenada Z
 El área superficial del sistema es constante.
 Catalizador no poroso

1. Definir el sistema coordenado “rectangulares”

2. Tomar un elemento diferencial de la misma geometría del sistema
∆z

∆y

∆x

ΔV = Δz∗Δy∗Δx

ΔV = Δz∗A s
3. Balance general para la especie A

entra por difusión en el plano z ) ( por difusiónen el plano z +∆ z ) ( reacción química en el volumen A ∆ z ) (
( Flujo molar de la especie A que − Flujo molar de A que sale + Rapidez de generación de A por = Ra
s

El sistema al estar en estado estacionario el término de la rapidez de acumulación y

rapidez de generación son igual a 0, por tanto, el balance queda de la siguiente
manera:
N A A s|z −N A A s|z +∆ z =0 … … … ec ( 1 )

Ahora divido la ecuación entre el volumen diferencial

 N A|z− N A|z +∆ z
=0
Δz
Aplicando la definición de limite

− lim
Δ Z →0
( N A|Z + Δ Z − N A|Z=Z
Δz )
=¿ … … … ec (2)¿

Evaluar el límite cuando z tiende a 0 obtenemos

−ⅆ N A
=0 … … … ec(3)
ⅆz
Resolviendo la ecuación diferencial por variables separables:
N A =C 1 … … … ec ( 4 )
Aplicando la primera ley de Fick
d yA
N A =−c D AB + y A (N A + N B )… … … ec (5 )
dz
De la estequiometria de la reacción se sabe que:

N B=−2 N A … … … ec ( 6 )
Sustituyendo la ecuación

d yA
N A =−c D AB + y A (N A −2 N A )… … … ec (7 )
dz

d yA
N A =−c D AB −yA NA
dz
d yA
N A + y A N A=−c D AB
dz
d yA
N A (1+ y A )=−c D AB
dz
d yA
−c D AB
dz
N A=
(1+ y ¿¿ A )¿
−c D AB d y A
N A= … … … ec ( 8 )
(1+ y ¿¿ A)dz ¿
 Resolviendo la ecuación diferencial por variables separables, además de utilizar
sustitución simple, y finalmente aplicando los respectivos limites:

z=δ y A = y Aδ
d yA
∫ N A dz=−c D AB ∫ (1+ y ¿¿ A )
… … … ec ( 9 ) ¿
z=0 y A= y A 0

u=(1+ y ¿¿ A)¿
du=dy A

y A = y Aδ
du
N A δ=−c D AB ∫ u
y A = y A0

N A δ=−c D AB ln ¿ ¿

1+ y Aδ
N A δ=−c D AB ln ⁡( )
1+ y A 0

1+ y Aδ
N A δ=−c D AB ln ⁡( )
1+ y A 0

−c D AB 1+ y Aδ
N A= ln ⁡( )… … … ec ( 10 )
δ 1+ y A 0

b) Solo se produce dióxido de carbono en la superficie del carbón;

C+ O2 → C O2
Solución

Esquema de interpretación del sistema:

Por conveniencia el punto de referencia z = 0 se cambió a la superficie de la película

estancada.

z=0 yA = yA0

z=δ yA = yAδ
Consideraciones para resolver el sistema

 El sistema opera en estado estacionario

 La concentración yA0 > yAδ
 La temperatura y presión permanecen constantes por lo que el coeficiente de
difusividad DAB permanece constante
 El sistema está concentrado, por lo tanto, el mecanismo de transporte que
predomina en el sistema es convectivo.
 El flux molar es unidireccional en la coordenada Z
 El área superficial del sistema es constante.
 Catalizador no poroso

1. Definir el sistema coordenado “rectangulares”

2. Tomar un elemento diferencial de la misma geometría del sistema
∆z

∆y

∆x

ΔV = Δz∗Δy∗Δx

ΔV = Δz∗A s
3. Balance general para la especie A

entra por difusión en el plano z ) ( por difusiónen el plano z +∆ z ) ( reacción química en el volumen A ∆ z ) (
( Flujo molar de la especie A que − Flujo molar de A que sale + Rapidez de generación de A por = Ra
s

El sistema al estar en estado estacionario el término de la rapidez de acumulación y

rapidez de generación son igual a 0, por tanto, el balance queda de la siguiente
manera:
N A A s|z −N A A s|z +∆ z =0 … … … ec ( 1 )
Ahora divido la ecuación entre el volumen diferencial
N A|z− N A|z +∆ z
=0
Δz
Aplicando la definición de limite

− lim
Δ Z →0
( N A|Z + Δ Z − N A|Z=Z
Δz )=¿ … … … ec (2)¿

Evaluar el límite cuando z tiende a 0 obtenemos

−ⅆ N A
=0 … … … ec(3)
ⅆz
Resolviendo la ecuación diferencial por variables separables:
N A =C 1 … … … ec ( 4 )
Aplicando la primera ley de Fick
d yA
N A =−c D AB + y A (N A + N B )… … … ec (5 )
dz
De la estequiometria de la reacción se sabe que:

N B=−N A … … … ec ( 6 )

Sustituyendo la ecuación 6 en la ecuación 5

d yA
N A =−c D AB + y A (N A −N A )
dz
d yA
N A =−c D AB … … … ec ( 7 )
dz
Resolviendo la ecuación diferencial por variables separables y aplicando límites de
integración con respecto al diagrama:
z=δ y A = y Aδ

∫ N A dz=−c D AB ∫ d y A … … … ec ( 8 )
z=0 y A= y A 0

N A δ=−c D AB ( y Aδ − y A 0 )
 c D AB
N A= ( y A 0 − y Aδ )… … … ec ( 10 )
δ

c) Se produce la siguiente reacción instantánea en la superficie del carbón

4 C+3 O2 →2 CO+2 C O2
Solución

Esquema de interpretación del sistema:

Por conveniencia el punto de referencia z = 0 se cambió a la superficie de la película

estancada.

z=0 yA = yA0

z=δ yA = yAδ

Consideraciones para resolver el sistema

 El sistema opera en estado estacionario

 La concentración yA0 > yAδ
 La temperatura y presión permanecen constantes por lo que el coeficiente de
difusividad DAB permanece constante
 El sistema está concentrado, por lo tanto, el mecanismo de transporte que
predomina en el sistema es convectivo.
 El flux molar es unidireccional en la coordenada Z
 El área superficial del sistema es constante.
 Catalizador no poroso

1. Definir el sistema coordenado “rectangulares”

2. Tomar un elemento diferencial de la misma geometría del sistema
∆z

∆y

∆x
 ΔV = Δz∗Δy∗Δx

ΔV = Δz∗A s
3. Balance general para la especie A

entra por difusión en el plano z ) ( por difusiónen el plano z +∆ z ) ( reacción química en el volumen A ∆ z ) (
( Flujo molar de la especie A que − Flujo molar de A que sale + Rapidez de generación de A por = Ra
s

El sistema al estar en estado estacionario el término de la rapidez de acumulación y

rapidez de generación son igual a 0, por tanto, el balance queda de la siguiente
manera:
N A A s|z −N A A s|z +∆ z =0 … … … ec ( 1 )

Ahora divido la ecuación entre el volumen diferencial

N A|z− N A|z +∆ z
=0
Δz
Aplicando la definición de limite

− lim
Δ Z →0
( N A|Z + Δ Z − N A|Z=Z
Δz )
=¿ … … … ec (2)¿

Evaluar el límite cuando z tiende a 0 obtenemos

−ⅆ N A
=0 … … … ec(3)
ⅆz
Resolviendo la ecuación diferencial por variables separables:
N A =C 1 … … … ec ( 4 )
Aplicando la primera ley de Fick
d yA
N A =−c D AB + y A (N A + N B + N C )… … … ec ( 5 )
dz
De la estequiometria de la reacción se sabe que:
 3 N B =−2 N A
−2
N B= N … … … ec ( 6 )
3 A
3 N C =−2 N A
−2
N C= N … … … ec (7 )
3 A

Sustituyendo las ecuaciones 6 y 7 en la ecuación 5

d yA 2 2
N A =−c D AB + y A (N A − N A − N A ) … … … ec ( 6 )
dz 3 3

d yA 1
N A =−c D AB − y N
dz 3 A A
1 d yA
N A+ y A N A=−c D AB
3 dz
1 d yA
N A (1+ y )=−c D AB
3 A dz
d yA
−c D AB
dz
N A=
1
(1+ y ¿ ¿ A)¿
3
−c D AB d y A
N A=
1
(1+ y ¿ ¿ A)dz … … … ec ( 7 ) ¿
3
Resolviendo la ecuación diferencial por variables separables, además de utilizar
sustitución simple, y finalmente aplicando los respectivos limites:

z=δ y A = y Aδ
d yA
∫ N A dz=−c D AB ∫ 1
… … … ec ( 8 ) ¿
z=0 y A= y A 0
(1+ y ¿¿ A )
3

1
u=(1+ y ¿ ¿ A)¿
3
 1
du= dy
3 A
3 du=dy A

y A = y Aδ
du
N A δ=−c D AB 3 ∫ u
y A = y A0

N A δ=−3 c D AB ln ¿ ¿

1
1+ y
3 Aδ
N A δ=−3 c D AB ln ⁡( )
1
1+ y
3 A0

1
1+ y
−3 c D AB 3 Aδ
N A= ln ⁡( )… … … ec ( 9 )
δ 1
1+ y
3 A0

26.12

El siguiente diagrama ilustra la difusión en fase gaseosa en la vecindad de una

superficie catalítica. El componente A se difunde a través de una película estancada
hasta la superficie catalítica donde instantáneamente se convierte en B. El producto B
se difunde alejándose de la superficie catalítica, otra vez hacia la película estancada.
 a) Si la reacción sobre la superficie es A → B, determinar una ecuación que
predecirá el perfil de concentraciones; o sea, yA a varios niveles de z.
b) Si la reacción sobre la superficie es A → 3B determinar una ecuación que
predecirá el perfil de concentraciones.

a.) Solución

Esquema de interpretación del sistema:

Por conveniencia el punto de referencia z = 0 se cambio a la superficie de la película

estancada.

z=0 yA = yA0

z=δ yA = yAδ

Consideraciones para resolver el sistema

 El sistema opera en estado estacionario

 La concentración yA0 > yAδ
 La temperatura y presión permanecen constantes por lo que el coeficiente de
difusividad DAB permanece constante
 El sistema está concentrado, por lo tanto, el mecanismo de transporte que
predomina en el sistema es convectivo.
 El flux molar es unidireccional en la coordenada Z
 El área superficial del sistema es constante.
 Catalizador no poroso
4. Definir el sistema coordenado “rectangulares”
5. Tomar un elemento diferencial de la misma geometría del sistema
∆z

∆y

∆x

ΔV = Δz∗Δy∗Δx

ΔV = Δz∗A s
6. Balance general para la especie A

entra por difusión en el plano z ) ( por difusiónen el plano z +∆ z ) ( reacción química en el volumen A ∆ z ) (
( Flujo molar de la especie A que − Flujo molar de A que sale + Rapidez de generación de A por = Ra
s

El sistema al estar en estado estacionario el término de la rapidez de acumulación y

rapidez de generación son igual a 0, por tanto, el balance queda de la siguiente
manera:
N A A s|z −N A A s|z +∆ z =0 … … … ec ( 1 )

Ahora divido la ecuación entre el volumen diferencial

N A|z− N A|z +∆ z
=0
Δz
Aplicando la definición de limite

− lim
Δ Z →0
( N A|Z + Δ Z − N A|Z=Z
Δz )
=¿ … … … ec (2)¿

Evaluar el límite cuando z tiende a 0 obtenemos

−ⅆ N A
=0 … … … ec(3)
ⅆz
Resolviendo la ecuación diferencial por variables separables:
 N A =C 1 … … … ec ( 4 )
Aplicando la primera ley de Fick
d yA
N A =−c D AB + y A (N A + N B )… … … ec (5 )
dz
De la estequiometria de la reacción se sabe que:

N B=−N A … … … ec ( 6 )

Sustituyendo la ecuación 6 en la ecuación 5

d yA
N A =−c D AB + y A (N A −N A )
dz
d yA
N A =−c D AB … … … ec ( 7 )
dz
Resolviendo la ecuación diferencial por variables separables y aplicando límites de
integración con respecto al diagrama:
z=δ y A = y Aδ

∫ N A dz=−c D AB ∫ d y A … … … ec ( 8 )
z=0 y A= y A 0

N A δ=−c D AB ( y Aδ − y A 0 )… … … ec ( 9 )

c D AB
N A= ( y A 0 − y Aδ )… … … ec ( 10 )
δ
Ahora si integro desde 0 a una distancia z y de una concentración yA0 a una
concentración yA
z=z y A= y A

∫ N A dz=−c D AB ∫ d y A … … … ec ( 11 )
z=0 y A = y A0

N A z=−c D AB ( y A − y A 0)… … … ec ( 12 )
Si divido la ecuación 12 entre la ecuación 9 obtengo el perfil de concentraciones
N A z −c D AB ( y A − y A 0)
=
N A δ −c D AB ( y Aδ − y A 0)
 z ( y A− y A 0 )
=
δ ( y Aδ − y A 0 )
z
( y Aδ − y A 0 ) δ = y A − y A 0

z
y A =( y Aδ − y A 0 ) + y A 0 … … … ec (13 )
δ
b.) Solución
El procedimiento que da solución a este inciso es análogo al inciso anterior por lo
que hasta la ecuación 5 del inciso a, el procedimiento se obviara. Además, a modo
de generalizar la resolución de este problema la reacción química tendrá los
siguientes coeficientes estequiométricos:
nA → mB

n N B=−m N A
−m
N B= N A … … … ec ( 14 )
n

Sustituyendo la ecuación 9 en la ecuación 5

d yA m
N A =−c D AB + y A (N A − N A )… … … ec (15 )
dz n

d yA m
N A =−c D AB + y A N A (1− )
dz n
m d yA
N A − y A N A (1− )=−c D AB
n dz
n−m d yA
N A (1− y ¿ ¿ A )=−c D AB ¿
n dz
d yA
−c D AB
dz
N A=
n−m
(1− y ¿ ¿ A )¿
n
 −c D AB d y A
N A=
n−m
(1− y ¿ ¿ A )dz ¿
n
z=δ y A = y Aδ
d yA
∫ N A dz=−c D AB ∫ n−m
… … … ec ( 16 ) ¿
z=0 y A= y A 0
(1− y ¿¿ A )
n
Resolviendo la ecuación diferencial por variables separables, además de utilizar
sustitución simple, y finalmente aplicando los respectivos limites:

n−m
u=(1− y ¿ ¿ A )¿
n
n−m
du=−dy A
n
−n
du=dy A
n−m

y A = y Aδ
n du
n−m y ∫
N A δ=−c D AB−
=y
A
u A0

n
N A δ=c D AB ln ¿ ¿
n−m

n−m
1− y Aδ
n n
N A δ=c D AB ln ⁡( )
n−m n−m
1− y A 0
n

n 1− y Aδ
N A δ=c D AB ln ⁡( )… … … ec ( 17 )
n−m 1− y A 0

Proponiendo la siguiente expresión

n−m
β= … … … ec ( 18 )
n

1 n
= … … … ( 19 )
β n−m

Sustituyendo la ecuación 19 en la ecuación17 se obtiene

 c D AB 1− y Aδ
N A δ= ln ⁡( )… … … ec ( 20 )
β 1− y A 0

Ahora si integro desde 0 a una distancia z y de una concentración yA0 a una

concentración yA
z=z y A= y A
d yA
∫ N A dz=−c D AB ∫ n−m
… … … ec (21 ) ¿
z=0 y A = y A0
(1− y ¿¿ A )
n

n 1− y A
N A z=c D AB ln ⁡( )… … … ec ( 22 )
n−m 1− y A 0
Si divido la ecuación 22 entre la ecuación 17 obtengo el perfil de concentraciones
n 1− y A
c D AB
ln ⁡( )
NA z n−m 1− y A 0
=
N Aδ n 1− y Aδ
c D AB ln ⁡( )
n−m 1− y A 0
1− y A
ln ⁡( )
z 1− y A 0
=
δ 1− y Aδ
ln ⁡( )
1− y A 0
1− y Aδ z 1− y A
ln ⁡( ) =ln ⁡( )
1− y A 0 δ 1− y A 0

( )
z
1− y Aδ δ 1− y A
=
1− y A 0 1− y A 0

y A =1−¿

26.23

Se utilizará Hexafluoruro de Uranio, UF6, como combustible en un reactor experimental

de energía nuclear. Para el diseño del reactor se requiere la rapidez de producción de
UF6. El UF6 se prepara por exposición de pelets de uranio puro de forma esférica al gas
flúor a 1000 ºK y 1 atm de presión.

F2(g)

UF6(g)

Se piensa que la etapa que controla el proceso es la difusión molecular del Fluor hasta
la superficie de los pelets. Si la reacción

U +3 F → U F6

Ocurre en forma irreversible e instantánea sobre la superficie del pelet de uranio y si la

difusividad del flúor a través del hexafluoruro de uranio es de 0.273 cm 2/s a 1000 K y 1
atm, determinar los moles de UF 6 que se producen cuando el diámetro del pelet es de
0.4 cm.

Datos Suposiciones
 DFUF6 = 0.273 cm2/s.  Estado pseudo - estacionario.
 T = 1000 ºK  YA0 > YA
 P = 1 atm.  Gas concentrado.
 P, T = ctes  DAB = cte.
 JA unidireccional (eje r).
 As cte.
 Contra difusión.
1) Ssitemas de coordenadas
Rectangulares
2) Elemento diferencial:
 r

3) Flujo molar
Para determinar el flujo molar es necesario aplicar la siguiente ecuación:

(−W A ) =N A A S =( 4 π r 2 ) N A ec(1)

De la ley de Fick sabemos:

d yA
N A =−c D AB + y A (N A + N B ) ec(2)
dr
De la estequiometria se sabe:
−1
N =N B
3 A
Sustituyendo en (2):
d yA −1
N A =−c D AB + y A (N A N )
dr 3 A
d yA 2
N A =−c D AB + y A( N )
dr 3 A
Despejando N A :

−c D AB d y A
N A= ec(3)
2 dr
(1− y A )
3

Sustituyendo (3) en (1):

c D AB d y A
(−W A ) =( 4 π r )
2
2 dr ec(4)
(1− y A )
3
Resolviendo (4) por separación de variables:

dr d yA
(−W A ) 2
=( 4 π ) c D AB
2
r (1− y A )
3

Aplicando límites de integración:

r=r ∞ y A = y A∞
dr d yA
(−W A ) ∫ r
2
=( 4 π ) c D AB ∫
2
r= R y =y
A
(1− y A )
AR

Resolviendo la ecuación deferencial, donde y A ∞ = 1 y y A R=0:

(−W A ) =6 c D AB πR ln 3 ec(5)
De la estequiometria sabemos que:

1
W =W B
3 A

∴ W B =2 c D AB πR ln 3 ec(6)

Donde:
P
c=
RT
Sustituyendo datos en (6):

(( ) ( )
2
101,325 Pa m
W B=2 0.0000273 ( 0.002m ) πLn 3
)
Pam
3 s
8.314 ( 1000 K )
molK

−6 mol
W B=4.59329 x 10
s
26.24.

Determinar el tiempo que se requiere para reducir el diámetro del pelet de uranio que
se describe en el problema 26.23 de 0.4 a 0.2 cm. El peso específico del uranio es
19.05.

Solución.

Para encontrar el tiempo en que el pelet reducirá su diámetro desde 0.4 a 0.2 cm, es
necesario recordar cómo se define matemáticamente qué tan rápido se consumen los
moles del sólido y la relación con su radio, esto para cualquier instante de tiempo:

4 3
d( π R )
−d mc −d (ρ c V ¿ M c ) − ρc 3 −ρ c dR
= = = 4 π R2 (ec .1)
dt dt Mc dt Mc dt

De acuerdo a la estequiometria de la reacción, por cada mol de uranio que se

consume, se consumen 3 moles de flúor. Esto representa que 3 moles de flúor están en
contacto con el pelet de uranio, siendo así que el pelet se reduce 3 veces más rápido
de lo que se consume el flúor. Esto significa que:

d mc
−3 =¿ ¿
dt

Tomando la expresión del flujo de A obtenida en el problema anterior, y la ec.1, se

obtiene lo siguiente:

ρc 2 dR
−3 4πR =C c D AB 6 π ln ⁡(3) R(ec .2)
Mc dt

Reacomodando, es posible llegar a la siguiente expresión:

ρc
dt=−2 RdR
Cc D AB M c ln ⁡(3)
Se quiere calcular el tiempo que tardará el pelet en reducir su diámetro de 0.4 cm a 0.2
cm, o 0.2 cm a 0.1 cm, o 0.002 m a 0.001 m. Estos últimos dos valores serán nuestros
límites inferior y superior de la integral, respectivamente.

t 0.001 m
ρc
∫ dt =¿−2 ∫ RdR ¿
Cc D AB M c ln ⁡(3) 0.002 m
0

ρc
t=−2 ¿¿
C c D AB M c ln ⁡(3)

Tomando los datos del problema anterior se tiene que:

2 2
cm cm
D AB=0.273 =0.273 ¿
s s

Para obtener la densidad del uranio se utiliza el dato brindado por el problema, el cuál
menciona que el uranio tiene un peso específico de 19.05. Se sabe que el peso
específico es igual a la densidad de la especie química de interés entre la densidad del
agua, la cual es de 1000 kg/m3, sabiendo esto, se calcula la densidad del uranio:

ρc kg
19.05= ρ c=19050 3
kg m
1000 3
m

Por su parte, el peso molecular del uranio es el siguiente:

MC=238.029 g/mol=0.238029 kg/mol

Fuente: Uranium - Element information, properties and uses | Periodic Table. (s. f.).
https://www.rsc.org/periodic-table/element/92/uranium.

Finalmente, para el cálculo de la concentración del gas es necesario usar la ecuación

del gas ideal:

PV =nRT ( ec .4 )

La concentración es igual al número de moles presentes en un volumen dado, por lo

que, reacomodando se obtiene que:
 n P
=Cc = (ec .5)
V RT

La presión que nos da el problema es de 1 atm o 101325 Pa, y la temperatura es de

1000 K. La constante universal de los gases, R, es igual a 8.314 Pam 3/molK.
Sustituyendo valores en la ec. 5 se obtiene la concentración del flúor:

101325 Pa mol
C c= =12.1873
( )
3
Pa m3 m
8.314 ( 1000 K )
molK

Ya se dispone de todos los valores necesarios para calcular el tiempo que tardará el
pelet en disminuir su diámetro de 0.4 cm a 0.2 cm, por lo que se procede a sustituir
valores en la ec. 3:

kg
−2(19050 3
)
m
t= ¿¿
( )
2
mol −5 m kg
12.1873 3 2.73× 10 (0.238029 )ln ⁡(3)
m s mol

s
t=−437.9056 ×1 0
6
2
(−1.5 ×1 0−6 m2 )
m

t=656.8584 s=10.94764 min

Esto quiere decir que, para que el diámetro del pelet se reduzca de 0.4 cm a 0.2 cm, se
requieren de 10.94763 minutos.
26.28
La varilla cilíndrica de grafito (carbono puro) de 20 cm de largo se inserta en una
atmosfera oxidante a 1145 K y 1.013 ∙1 05 Pa de presión. El proceso oxidante está
limitado por la difusión a contraflujo de oxigeno con el monóxido de carbono que se
forma sobre la superficie catalítica. Bajo las condiciones del proceso de combustión,
2
m
puede suponerse que la difusividad del oxígeno en la mezcla gaseosa es de 1 ∙10−5 .
s
a) Determinar los moles de CO que se producen en la superficie de la varilla por
segundo en el momento en que el diámetro de la varilla es 1 cm y la concentración de
oxígeno a 1 cm de distancia radial de la varilla es .4 en mol. Suponer que es un proceso
en estado estacionario.
b) ¿Cuál sería la composición del oxígeno a una distancia radial de 1 cm desde el centro
de la varilla?
O2( g) +2 C( s ) →2 C O(g)

A+2 C → 2 B

r =R 1 r =R 2
y A =0 y A = y A 0

Suposiciones
- El proceso ocurre en estado estacionario
- y A 0> y A
- El sistema es concentrado, por lo tanto, los mecanismos de transporte de materia
son tanto la difusión como la convección libre
- El proceso ocurre a Temperatura y presión constantes, por lo tanto, el coeficiente
de difusión también es constante
- El flux es unidireccional (r)
- El área perpendicular al flujo es variable
- Re. laminar
- En la superficie catalítica la reacción es instantánea

Datos
R1=1 cm

R2=2 cm

L=20 cm
y A 0=.4
2
−5 m
D AB=1 ∙ 10
s
5
P=.013 ∙ 10 Pa
T =1415 K
I) Sistema de coordenadas: Cilíndricas
II) Elemento diferencial
r

V =πLr ∆ r
A s=2 Lπr

III) Balance de Materia

J A 2 π ∆ rL|r=r −J A 2 π ∆ rL|¿r +∆ r =0 … … … ec (1)

d yA
N A =C D AB + y A ( N A + N B ) … … … .ec ( 2)
dr
De la estequiometria de la reacción se sabe que:
N B=−2 N A … … … . ec ( 3 )

Sustituyendo la ec (3) en ec (2)

d yA
N A =C D AB + y A ( N A −2 N A )
dr
 d yA
N A =C D AB + y A (−N A )
dr
C D AB d y A
N A=
( 1+ y A ) dr
Dado que se están evaluando las condiciones en la superficie del catalizador la taza de
consumo de A será específicamente para ese punto.
(−W A )|r= R=N A ∙ A S

(−W A )|r= R=
( C D AB d y A
( 1+ y A ) dr )
( 2 Lπr ) … … … ec (4 )

( )
R2 yA 0
dr dy
(−W A )|r= R∫ =( 2 Lπ C D AB ) ∫ ( 1+ yA )
R r 1 0 A

(−W A )|r= R ln
( ( )) R2
R1
=2 LπC D AB ln (1+ y A 0 )

2 LπC D AB ln (1+ y A 0)
(−W A )|r= R= … … … ec (5)

( )
R
ln 2
R1

IV) Sustituyendo valores en la Ec .5

La concentración de oxigeno no es proporcionada, sin embargo, se puede obtener
mediante un modelo de gases reales y lo datos dados por el problema.
El factor acéntrico (Z) se obtuvo mediante la ecuación de Van der Waals.
5
P (1.013 ∙10 Pa) mol
C= = =10.642
( ) m3
zRT Pa∙ m 3
( 1.0002 ) 8.314 ( 1145 K )
mol ∙ K

( m ) ( s )
2
mol m −5
2 π ( .2 m ) 10.642 1 ∙10
3
∙ ln ( 1+.4 )
(−W A )|r= R =
ln ( 21 )
mol
(−W A )|r= R=0.00065 2
s ∙m
(−W A )|r= R=−2 (−W B )|r =R
 mol
(−W B )|r =R=−0.0013 2
s∙m
Nota: El signo negativo no representa que se esté desplazando una cantidad negativa
de materia, si no que el flujo de la especie B tiene sentido opuesto al flujo de la especie
A.
mol
PCO =0.0013 a¿
s

V) Retomando la ec (3)
Esta vez se desarrollará un perfil general, por lo que los limites superiores serán
valores generales los cuales podrán ser sustituidos para cualquier valor de r entre los
límites del sistema ( R1 y R 2 ¿
r yA
d yA
N A ∫ dr =C D AB∫ … … … ec ( 6 )
R1 0 ( 1+ y A )
R2 yA 0
d yA
N A ∫ dr =C D AB ∫ … … … ec (7 )
R1 0 ( 1+ y A )

NA ( r −R 1 )=C D AB ln(1+ y A )

Dividiendo ambos perfiles se obtiene:

NA ( r −R1 )=C D AB ln (1+ y A )
N A ( R2−R1 ) =C D AB ln (1+ y A 0)

( r −R1 ) ln (1+ y A )
= … … … ec ( 8 )
( R2−R1 ) ln ( 1+ y A 0 )

Sustituyendo valores en la ec ( 7 ) :

( 1−1 ) cm ln(1+ y A )
=
( 2−1 ) cm ln(1+.4 )
ln ( 1+ y A )=0

y A =0 b ¿

Esto solo confirma que el perfil obtenido en la ecuación es correcto dado que al
considerar que la reacción en la superficie catalítica es instantánea, la concentración de
A en dicha superficie es 0, lo cual se consideró desde un inicio para el desarrollo del
perfil general mostrado en la ec ( 8 ) .
Conclusiones
En la ec ( 8 ) , se observa una relación directamente proporcional entre el radio al exterior
del cilindro y la concentración del oxígeno.
Al considerarse la reacción instantánea y esta solo ocurrir en la superficie del
catalizador, el único factor del que depende la generación de monóxido de carbono es
del transporte de materia del oxígeno hacia la superficie catalítica.
El coeficiente estequiométrico del solido no es relevante para los fluxes de las especies
A y B, dado que este no se desplaza.