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

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 … … … ( 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 )=¿ … … …(2)¿

Evaluar el límite cuando z tiende a 0 obtenemos

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

N B=−N A … … … ( 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 … … … (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 … … … ( 8)
z=0 y A= y A 0

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

c D AB
N A= ( y A 0 − y Aδ )… … … ( 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 … … … (11 )
z=0 y A = y A0

N A z=−c D AB ( y A − y A 0)… … … ( 12 )
Si divido la ecuación 12 entre la ecuación 9
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 … … … ( 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 … … … ( 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 )… … … ( 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
… … … ( 1 6) ¿
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 A δ=−c D AB− ∫
n−m y = y u A 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 ⁡( )… … … ( 1 7 )
n−m 1− y A 0
Proponiendo lo siguiente

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

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

Sustituyendo la ecuación 19 en 17 se obtiene

c D AB 1− y Aδ
N A δ= ln ⁡( )… … … ( 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
… … … ( 21 ) ¿
z=0 y A = y A0
(1− y ¿¿ A )
n

n 1− y A
N A z=c D AB ln ⁡( )… … … ( 2 2 )
n−m 1− y A 0
Si divido la ecuación 22 entre la ecuación 17
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−¿