ESCUELA POLITÉCNICA NACIONAL

TRANSFERENCIA DE MASA

NOMBRE: Tania Guerrero R

FECHA: 09 – 04 – 2012

DIFUSIÓN SIMULTÁNEA Y REACCIÓN REVERSIBLE EN

UNA CAPA CATALÍTICA
PLANTEAMIENTO DEL PROBLEMA:

La reacción catalítica en fase gaseosa entre los componentes Ay B es reversible en una capa de
catalizador en un determinado punto de un reactor monolítico.

2A↔B
[ 1]
La expresión de velocidad de la reacción catalítica para el reactivo A viene dada por:

CB
(
r ´ A =−k C A −
2
KC )
[ 2]
Donde la velocidad viene expresada en g-mol/cm 3*s, la constante de velocidad calculada, k,
tiene un valor de 8 x 10 4 cm3/g-mol y la constante de equilibrio tiene un valor de Kc = 6x 10 5
cm3/g-mol. Lacapa de catalizador tiene un espesor L = 0,2 cm y la difusividad efectiva de A en B
para esta capa es De= 0,01 cm2/s. La mezcla de reactivos sólo contiene los gases A y B, por lo
que sólo se necesita considerar la difusión de gases binarios. La concentración total de A y B es
CT= 4x 10-5 g-mol/cm3.

El cálculo del factor de efectividad requiere que se considerela difusión simultánea de los dos
componentes junto con el balance de materia de la expresión de velocidad de la reacción
reversible.

a) Utilice una técnica de diferencia finita implícita para calcular el factor de efectividad
para la reacción dada y represente C A y NA en una tabla de 10 intervalo igualmente
espaciados dentro de la capa de catalizador. Las concentraciones de reactivos en la
superficie de la capa tienen los calores conocidos de C As = 3x10-5 g-mol/cm3 y CBs = 1x10-5
g-mol/cm3.
b) Repita el apartado a) utilizando la técnica del dispara para resolver el problema, e
introduzca los valores calculados en la misma tabla.
c) Compare las dos soluciones de los apartados a) y b). Comente que técnica de
resolución prefiere.

Balances de materia A y B en una placa porosa

El balance de materia del componente A en la dirección z desde la parte superior de la

placa porosa se representa mediante una ecuación diferencial ordinaria de flujo N A:
 d NA CB
dz
2
(
=−k C A −
KC )
[ 3]
Con la condición frontera de que no hay flujo de transferencia de masa en la parte inferior
de la placa porosa, expresado como:

N A|z= L=0
[4]
Para el componente B se podría derivar una ecuación diferencial similar, pero es más fácil
utilizar una relación algebraica entre los flujos molares de A y B, en base a la estequiometria de
la reacción.
−1
NB= N
2 A
[ 5]
Obsérvese que esta ecuación también satisface la condición frontera en el fondo de la capa
porosa donde:
N B|z= L =0
[ 6]
Ley de Fick para la difusión binaria:

La difusión de A en este sistema binario se describe mediante:

d CA
N A =−D e + X A ( N A+ N B)
dz
[ 7]
Donde la difusividad efectiva De de la capa de catalizador viene dada por:

D AB ε p σ
De = ~τ
[ 8]
Donde DAB es la difusividad binaria, εp es la porosidad del catalizador, σ es el factor de
constricción y es la tortuosidad.

Se puede reordenar la Ecuación y la definición de X A en términos de concentración, para dar

lugar a:
CA NA
d C A CT 2
=
( )
−N A

dz De
[ 9]
Con la condición inicial:

C A|z=0=C A S

[ 1 0]
Se podría utilizar una ecuación diferencial similar para el componente B; sin embargo, es más
fácil utilizar el balance total de masa para este sistema en fase gaseosa para calcular C B. Por lo
tanto:
C B=C T −C A
[ 11 ]
Por lo que la Ecuación se puede expresar:

d NA ( C T −C A )
dz [
=−k C A2−
KC ]
[ 12 ]

Figura 1: Difusión con reacción química reversible en una capa de catalizador

RESOLUCIÓN

Este problema requiere resolución simultanea de la ecuación [ 9 ] , con la condición inicial en la

Ecuación diferencial [ 12 ] con la condición de frontera dada por la Ecuación [ 6 ] . Por lo que las
condiciones de frontera para las dos ecuaciones diferenciales están separadas, siendo una la
condición inicial y otra la condición final.

a) Resolución mediante una Diferencia Finita Implícita

Esta resolución implica utilizar fórmulas de diferencias finitas para las ecuaciones diferenciales
y las condiciones frontera. El sistema resultante de ecuaciones algebraicas se puede resolver
para las variables en varios puntos nodales utilizados en el problema. En la Figura 1 se
representan los 10 intervalos con 11 puntos nodales para este problema. La Ecuación (9) se
puede rescribir utilizando una fórmula de diferencia central de segundo orden para la primera
derivada como:
 CA NA
d C A C A −C A
= n
= n+ 1 n−1
( )( )
CT
n

2
−N A n

n
para ( 2 ≤n ≤ 10 )
dz 2∆ z De
[ 13 ]
La condición inicial de la Ecuación[ 10 ] da lugar a CAI = CAs.
Para el nodo 11 en el fondo de la capa de catalizador, se puede utilizar, para la primera
derivada, la aproximación de diferencia posterior de segundo orden, dando lugar a:
CA NA
dCA 3 C A −4 C A +C A ( )( )
11

CT 2
−N A 11

|
11
n 11 10 9
= =
dz n=11 2∆ z De
[ 15 ]

La condición frontera para NA, Ecuación [ 4 ] , del miembro derecho de la Ecuación [ 14 ] se

transforma en cero; por lo tanto, la simplificación resultante da lugar a la siguiente ecuación
algebraica en el nodo 11:

4 C A −C A
CA =11 ( 3
10 9

)
[ 16 ]
De forma similar, para el flujo NA, la Ecuación diferencial [ 12 ] se transforma en

d N A N A −N A C −C A )
dz
= n

2∆ z
2 ( T
=−k C A −
n+ 1

KC
para ( 2 ≤n ≤ 10 )
n −1

[ n
n

]
[ 17 ]

La condición final de la Ecuación [ 4 ] , da lugar a flujo cero en el fondo de la capa de catalizador.

N A =0
11

[ 1 8]
Para el nodo 1, en la parte superior de la capa de catalizador, se puede utilizar una
aproximación por diferencia anterior de segundo orden para la primera derivada, dando lugar
a:

−3 N A + 4 N A −N A C −C A )
dNA
dz | n=1
=
2∆ z
1 22 ( T
=−k C A −
KC
3

[ 1
1

]
[ 1 9]

Cálculo del factor de la efectividad

El factor de efectividad se puede obtener mediante un balance de materia en la superficie de

la capa de catalizador. En este caso, el flujo superficial de A que actúa sobre un área arbitraria
debe de ser igual a la velocidad de reacción media de la capa de catalizador por debajo del
área arbitraria. Por lo tanto, en términos de nodos diferenciales finitos para un área arbitraria
Ac.
N A AC =(−r ´ A ) AVG A C L
1

[ 20 ]
Por lo tanto, el factor de efectividad se puede expresar como:

NA 1

(−r ´ A )|AVG L
Ƞ= =
( −r ´ A )|SURFACE ( C T −C A
[k C A 2−
1
KC
1
)
]
[ 21 ]

Se resolvió POLYMATH usando el Solucionador de Ecuaciones No Lineales Simultaneas, pero se

lo resolvió de manera simultánea no lineal e implícita para evitar las referencias cíclicas de las
variables. Detallándose el programa en el siguiente ejemplo de cálculo

k = 8 * 10 ^ 4
Kc = 6 * 10 ^ 5
z = 0.2 / 10
De = 0.01
Ct = 4 * 10 ^ -5
Ca1 = 3 * 10 ^ -5
Na11 = 0
Na10 = 1.361 * 10 ^ -7

f(Ca2) = Ca2 - Ct * De * (Ca3 - Ca1) / (z * Na2) + 2 * Ct

Ca2(0) = 0
Ca2(min) = 0
Ca2(max) = 1

f(Na1) = Na1 - (2 * z * k * (Ca1 ^ 2 - (Ct - Ca1) / Kc) + 4 * Na2 - Na3) / 3

Na1(0) = 0
Na1(min) = 0
Na1(max) = 1

b) Resolución mediante una Diferencia Finita Implícita Método del Disparo

Los perfiles de flujo y la concentración se pueden determinar durante la resolución numérica

de la Ecuación diferencial [ 9 ] con una condición inicial conocida y la Ecuación diferencial [ 12 ],
que requiere la determinación de una condición inicial que satisfaga la condición frontera de
flujo cero cuando z= L, tal como se expresa en la Ecuación [ 4 ] .

Cálculo del factor de efectividad

El factor de efectividad para la capa de catalizador viene dado por la velocidad media de la
capa de catalizador dividida por la velocidad en la superficie de la capa de catalizador. Por lo
tanto:
 L L
CB CB

Ƞ=
∫ −k 1
0
( ( C A 2−
KC
dz
)) ( =
∫
0
C A2 −
KC
dz )
CB CB
−k 1 C A 2 − ( S
KC
L S

) ( C A 2−
S
KC
L S

)
[ 22 ]
que se puede diferenciar respecto a z, dando lugar a:

CB
dȠ
=
( C A2 −
KC )
dz CB
( C A 2−
S
KC
L S

)
[ 23 ]
La condición inicial para la Ecuación [ 22 ] es:

Ƞ|z=0 =0

[ 24 ]
y el factor de efectividad es el valor de r| a z = L. Por lo tanto, el incluir la Ecuación diferencial
[ 22 ] con la condición límite de la Ecuación [ 23 ] completa la formulación del problema.
Condición inicial estimada para NA

La técnica del disparo requiere que sea determinada la condición inicial de N A, la cual da un
valor final próximo a cero. Todas las técnicas requieren una suposición inicial de N A en la
superficie de la capa de catalizador. Con frecuencia esto suele ser difícil, ya que una condición
inicial incorrecta puede dar lugar a resultados numéricos que generan perfiles los cuales no
pueden ser optimizados fácilmente, para mejorar la suposición inicial.

Un balance de materia sencillo en la superficie del catalizador puede proporcionar un

mayorante o cota superior en el valor inicial cuando se supone una resistencia difusional
despreciable. Este balance establece que la velocidad de transferencia de masa molar en la
superficie es igual a la velocidad de reacción, en la capa de catalizador, que está teniendo
lugar a las concentraciones de A y B en la superficie. Por lo tanto:

( C T −C A )
NA |z=0=−r A L=kL
EST [ CA −
S
2
KC
S

]
[ 25 ]
Resolución del valor límite dividido

El Solucionador de Ecuaciones Diferenciales Ordinarias Simultáneas de POLYMATH se puede

utilizar con una "técnica de disparo", para resolver las ecuaciones diferenciales. Con esta
técnica, se puede determinar la condición inicial para N A.
Se resolvió POLYMATH usando el Solucionador de Ecuaciones Diferenciales Ordinarias
Simultáneas, detallándose el programa en el siguiente ejemplo de calculo:

De = 0.01
Ct = 4 * 10 ^ -5
k = 8 * 10 ^ 4
Kc = 6 * 10 ^ 5
err = Na1
z(0) = 0
z(f) = 1
d(Ca1)/d(z) = Ca1 * Na1 / (2 * Ct) - Na1
Ca1(0) = 0
d(Na1)/d(z) = -k * (Ca1 ^ 2 - (Ct - Ca1) / Kc)
Na1(0) = 4.271e-06

RESULTADOS

Diferencia Finita
NODO Método del Disparo
Z Implícita
S
CA NA CA NA
1 0 3E-05 4.271E-06 3E-05 4.250E-06
2 0.02 2.549E-05 3.0645E-06 2.526E-05 3.066E-06
3 0.04 2.165E-05 2.269E-06 2.154E-05 2.239E-06
4 0.06 1.887E-05 1.661E-06 1.866E-05 1.649E-06
5 0.08 1.657E-05 1.242E-06 1.644E-05 1.218E-06
6 0.1 1.493E-05 9.076E-07 1.475E-05 8.954E-07
7 0.12 1.362E-05 6.632E-07 1.349E-05 6.469E-07
8 0.14 1.273E-05 4.548E-07 1.258E-05 4.482E-07
9 0.16 1.209E-05 2.904E-07 1.196E-05 2.823E-07
10 0.18 1.174E-05 1.361E-0.7 1. 161E-05 1.362E-0.7
11 0.2 1.162E-05 0 1.149E-05 -5.150E-10
EFECTIVIDAD 0,3022 0,3007
TABLA 1. Comparación de las soluciones del método del disparo y de diferencia finita impícita

COMPARANDO LOS MÉTODOS DE RESOLUCIÓN

En comparación entre los dos método de cálculo se puede decir que mejor es el de Diferencia
Finita Implícita ya que por una parte cumple con la condiciones de borde establecidas y se
puede tener mayor precisión en los valores obtenidos usando más espacios nodales entre los
puntos de estudio en la capa de catalizador.