GUÍA 2: REACCIONES HETEROGÉNEAS

ANÁLISIS DE REACCIONES Y REACTORES

PRESENTADO POR:

JESÚS ALBERTO MIZGER ORTEGA

FECHA:

26 de junio de 2024

FACULTAD DE INGENIERÍA QUÍMICA

UNIVERSIDAD NACIONAL DE LITORAL

SANTA FE, ARGENTINA

2024
Problema 1:

𝐴+𝐵 →𝐶 𝑟 = 𝑘𝑣 𝐶𝐴
Datos

𝑇 = 700°𝐶 𝑃 = 1 𝑎𝑡𝑚
𝐷𝑐𝑎𝑡
𝐷𝑐𝑎𝑡 = 3𝑚𝑚 → 𝑟𝑐𝑎𝑡 = = 1,5𝑚𝑚 = 0,15 𝑐𝑚
2
𝐷𝑝
𝐷𝑝 = 100 Å = 1𝑥10−6 𝑐𝑚 → 𝑟𝑝 = ⁄ = 5𝑥10−7 𝑐𝑚
2
𝜀𝑠 = 0,5 𝑘𝑣 = 3,752 𝑠 −1

¿Cuál es el valor de factor de efectividad que podría esperarse?

𝜂 =?

1 𝛪1 (2𝜙) 𝑣𝑝 𝑘𝑣 𝜀𝑠 𝒟𝐴𝐵
𝜂= 𝜙= √ 𝔇𝑒 =
𝜙 𝛪0 (2𝜙) 𝑠𝑥 𝔇𝑒 𝜏

La difusividad molecular 𝒟𝐴𝐵 se puede calcular a partir de la fórmula de Chapman – Enskog

(predomina para poros mayores a 10.000 Å). Por otra parte, la difusividad Knudsen 𝐷𝑘 tiene en
cuenta las interacciones de la molécula con las paredes del poro, predomina cuando los poros
son menores a 1000 Å (las probabilidades de colisión con las paredes del poro son mayores) y
tiene la siguiente expresión:

𝑇
𝒟𝐴𝐾 = 9700𝑟𝑝 √
𝑃𝑀

𝑔𝑟
⁄𝑚𝑜𝑙 𝑐𝑚 973 𝐾
𝒟𝐴𝐾 = 9700 √ ∙ 5𝑥10−7 𝑐𝑚√
𝐾 𝑠 𝑔𝑟
106,167 ⁄𝑚𝑜𝑙

2
𝒟𝐴𝐾 = 0,014 𝑐𝑚 ⁄𝑠
Se procede a calcular la difusividad efectiva:

𝜀𝑠 𝒟𝐴𝐾
𝔇𝑒 =
𝜏
Los valores reportados en la literatura para la tortuosidad sugieren que este parámetro puede
encontrarse entre 3-7 según el libro de Hill, y entre 1,4-10 según Fogler. En ausencia de datos
experimentales se puede tomar un valor de 3 como promedio, por tanto:
2
0,5 ∙ 0,014 𝑐𝑚 ⁄𝑠 2
𝔇𝑒 = = 2,33𝑥10−3 𝑐𝑚 ⁄𝑠
3
Con el valor de la difusividad efectiva, se calcula el módulo Thiele:

0,15 𝑐𝑚 3,752 𝑠 −1
𝜙= √ 2 = 3,0096
2 2,33𝑥10−3 𝑐𝑚 ⁄𝑠

Con todos estos datos se calcula el factor de efectividad:

1 𝛪1 (2 ∙ 3,0096)
𝜂=
3,0096 𝛪0 (2 ∙ 3,0096)
Lo anterior puede resolverse haciendo uso de las funciones de Bessel, por tanto, nos queda.
1 62,4465
𝜂= = 0,3033
3,0096 68,4227
¿Cuál debería ser el tamaño de los cilindros para que en ensayos de laboratorio se pueda
obtener la cinética intrínseca?

𝑣𝑝 𝑘𝑣 𝑟 𝑘𝑣
𝜙= √ → 𝜙= √
𝑠𝑥 𝔇𝑒 2 𝔇𝑒

Despejando el radio:
2𝜙
𝑟=
𝑘𝑣
√
𝔇𝑒

Para cualquier tipo de geometría, si 𝜙 > 5 hay fuerte control difusional interno y el factor de
efectividad puede ser aproximado a 𝜂 ∼ 1/𝜙. En cambio, si 𝜙 < 1/3 las resistencias internas
al transporte de masa son despreciables indicando que 𝜂 ∼ 1. En el rango 1⁄3 ≤ 𝜙 ≤ 5 el
proceso sufre un control mixto, es decir que la velocidad de reacción intrínseca o cinética junto
con la velocidad de transferencia de masa controlan el proceso.

2 ∙ (1⁄3)
𝑟= = 0,0166 𝑐𝑚 = 0,166 𝑚𝑚
3,752 𝑠 −1
√ 2
2,33𝑥10−3 𝑐𝑚 ⁄𝑠

Para obtener un tamaño de cilindros con los cuales se pueda tener una cinética intrínseca, su
radio debe ser de 0,166 mm.

Problema 2:

𝐴→𝑃 𝑟 = 𝑘𝐶𝐴𝑚
Parte 1.
 𝑅 (𝑚 + 1) 𝑘𝑣 𝑚−1
𝜙= √ 𝐶
3 2 𝔇𝑒 𝑠

Partiendo de la ecuación del módulo de Thiele generalizado con respecto a la cinética de la

reacción:

𝑠 −1⁄2
𝑣𝑝 𝑟𝐴 (𝐶𝑠𝑠 ) 𝐶𝑠
𝜙= [∫ 𝔇𝑒 (𝐶′) 𝑟𝑣 (𝐶′) 𝑑𝑐]
𝑠𝑥 √2 𝐶𝑒𝑞

Para una esfera:

𝑠 −1⁄2
4𝜋𝑅 3 𝑟𝐴 (𝐶𝑠𝑠 ) 𝐶𝑠
𝜙= [∫ 𝔇𝑒 (𝐶′) 𝑟𝑣 (𝐶′) 𝑑𝑐]
3 ∙ 4𝜋𝑅 2 √2 𝐶𝑒𝑞

Simplificando:

𝑠 −1⁄2
𝑅 𝑟𝐴 (𝐶𝑠𝑠 ) 𝐶𝑠
𝜙= [∫ 𝔇𝑒 (𝐶′) 𝑟𝑣 (𝐶′) 𝑑𝑐]
3 √2 𝐶𝑒𝑞

Aplicando las condiciones de contorno:

𝐶𝑒𝑞 = 0 𝑖𝑟𝑟𝑒𝑣𝑒𝑟𝑠𝑖𝑏𝑙𝑒 → 𝑒𝑞𝑢𝑖𝑙𝑖𝑏𝑟𝑖𝑜

𝐶𝑠𝑠 = 𝐶𝑠
−1⁄2
𝑅 𝑘𝑣 𝐶𝑠𝑚 𝐶𝑠
𝜙= [∫ 𝔇𝑒 𝑘𝑣 𝐶𝑠𝑚 𝑑𝐶𝑠 ]
3 √2 0

𝔇𝑒 independiente de 𝐶 → 𝔇𝑒 = 𝑐𝑡𝑒
𝑅 𝑘𝑣 𝐶𝑠𝑚
𝜙= 1⁄
3 1⁄ 1⁄ 𝐶 2
√2𝔇𝑒 2 𝑘𝑣 2 [∫0 𝑠 𝐶𝑠𝑚 𝑑𝐶𝑠 ]

Integrando y aplicando los límites de integración, se tiene:

1⁄ 1⁄ 1
𝑅 𝑘𝑣 2 𝐶𝑠𝑚 𝑅 𝑘𝑣 2 𝐶𝑠𝑚 ∙ (𝑚 + 1) ⁄2
𝜙= 1⁄ ⇒ 𝜙= 1⁄ 𝑚⁄ 1⁄
3 1⁄ 𝐶 𝑚+1 2 3
√2𝔇𝑒 2 𝐶𝑠 2 ∙ 𝐶𝑠 2
√2𝔇𝑒 2 [𝑚𝑠 + 1 − 0]

Reorganizando:
1⁄
𝑚−1
𝑅 𝑘𝑣 2 1
𝜙= 1⁄ 𝐶𝑠
2
∙ (𝑚 + 1) ⁄2
3
√2𝔇𝑒 2

Finalmente se tiene que la expresión del módulo de Thiele generalizado para orden m es:
 𝑅 (𝑚 + 1) 𝑘𝑣 𝑚−1 𝑅 2 𝑚 + 1 𝑘𝑣 𝐶𝑠𝑚
𝜙= √ 𝐶 ⇒ 𝜙2 = ( ) ( )
3 2 𝔇𝑒 𝑠 3 2 𝔇𝑒 𝐶𝑠

El termino 𝑘𝑣 𝐶𝑠𝑚 representa la velocidad máxima posible y 𝔇𝑒 𝐶𝑠 el flujo máximo posible

difusivo dentro de la partícula.

Parte 2.

Aplicando el balance de masa:

1 𝑑 2 𝑑𝐶𝐴
𝔇𝑒 (𝑟 ) = 𝑘𝑣 𝐶𝐴𝑚 𝐸𝑐. 1
𝑟 2 𝑑𝑟 𝑑𝑟
Condiciones de contorno

𝑟=𝑅 𝐶𝐴 = 𝐶𝐴𝑠
𝑑𝐶𝐴
𝑟=0 =0 𝑝𝑜𝑟 𝑠𝑖𝑚𝑒𝑡𝑟í𝑎
𝑑𝑟
Adimensionalizando la ecuación 1, se tiene:
𝑟 𝐶𝐴
𝑟∗ = 𝐶𝐴∗ =
𝑅 𝐶𝐴𝑠

Despejando:

𝑟 = 𝑅 ∙ 𝑟∗ 𝐶𝐴 = 𝐶𝐴∗ 𝐶𝐴𝑠
Remplazando las variables adimensionales en la ecuación 1, se obtiene:

1 𝑑 𝑑(𝐶𝐴∗ 𝐶𝐴𝑠 )
𝔇𝑒 ((𝑅 ∙ 𝑟 ∗ )2
) = 𝑘𝑣 (𝐶𝐴∗ 𝐶𝑠𝑠 )𝑚
(𝑅 ∙ 𝑟 ∗ )2 𝑑(𝑅 ∙ 𝑟 ∗ ) 𝑑(𝑅 ∙ 𝑟 ∗ )

Resolviendo los paréntesis:

1 1 𝑠 𝑑 ∗2
𝑑𝐶𝐴∗ 𝑚
𝔇𝑒 𝐶𝐴 (𝑟 ) = 𝑘𝑣 𝐶𝐴𝑠
𝑟∗ 2 𝑅2 𝑑𝑟 ∗ 𝑑𝑟 ∗

Reorganizando los términos para tener una expresión del módulo de Thiele:

1 𝑑 ∗2
𝑑𝐶𝐴∗ 𝑅 2 𝑘𝑣 𝑠 𝑚−1 𝑚
(𝑟 ) = 𝐶 𝐶
𝑟∗2 𝑑𝑟 ∗ 𝑑𝑟 ∗ 𝔇𝑒 𝐴
𝑅2 𝑘𝑣 𝑠 𝑚−1
El termino 𝔇𝑒 𝐴
𝐶 , se puede expresar en función del módulo de Thiele, por lo que:

𝑘𝑣 𝑚−1 2
𝜙̅ = 𝑅 ∙ √ 𝐶𝐴𝑠 Por tanto: 𝜙̅ = 3 ∙ 𝜙√
𝔇𝑒 𝑚+1
 1 𝑑 ∗2
𝑑𝐶𝐴∗
(𝑟 ) = 𝜙̅ 2 𝐶 𝑚
𝑟∗2 𝑑𝑟 ∗ 𝑑𝑟 ∗

Ahora el factor de efectividad esta dado por la siguiente ecuación y se expresa en

coordenadas radiales:
1 𝑚 2
1 4 3 ∫ 𝑘𝐶 4𝜋𝑟 𝑑𝑟
𝑣 ∫ 𝑟𝐴 𝑑𝑣 𝜋𝑅
𝜂∗ = → 𝜂∗ = 3 𝑚
𝑟𝐴𝑠 𝑘𝐶𝐴𝑠
En variables adimensionales:
1 𝑠𝑚 𝑚 ∗2 ∗2 ∗
4 3 ∫ 𝑘𝐶𝐴 𝐶 4𝜋𝑟 𝑅 𝑅𝑑𝑟
𝜋𝑅 1
𝜂∗ = 3 𝑠 𝑚 = 3 ∫ 𝐶 𝑚 𝑟 ∗ 2 𝑑𝑟 ∗
𝑘𝐶𝐴 0

Llevando el sistema de ecuaciones a Matlab se obtiene la siguiente gráfica:

A valores bajos de modulo Thiele, el perfil se vuelve plano independientemente de la geometría

y valores elevados del módulo, el perfil se vuelve cada vez más lineal.
Parte 3.

Aplicando el balance de masa:

1 𝑑 2 𝑑𝐶𝐴
𝔇𝑒 (𝑟 ) = 𝑘𝑣 𝐶𝐴𝑚
𝑟 2 𝑑𝑟 𝑑𝑟
𝑃𝑎𝑟𝑎 𝑚 = 0
1 𝑑 2 𝑑𝐶𝐴
𝔇𝑒 (𝑟 ) = 𝑘𝑣 𝐸𝑐. 1
𝑟 2 𝑑𝑟 𝑑𝑟
Condiciones de contorno

𝑟=𝑅 𝐶𝐴 = 𝐶𝐴𝑠
𝑑𝐶𝐴
𝑟=0 =0
𝑑𝑟
Adimensionalizando la ecuación 1, se tiene:
𝑟 𝐶𝐴
𝑟∗ = 𝐶𝐴∗ =
𝑅 𝐶𝐴𝑠

Despejando:

𝑟 = 𝑅 ∙ 𝑟∗ 𝐶𝐴 = 𝐶𝐴∗ 𝐶𝐴𝑠
Remplazando las variables adimensionales en la ecuación 1, se obtiene:

1 𝑑 ∗ )2
𝑑(𝐶𝐴∗ 𝐶𝐴𝑠 )
𝔇𝑒 ((𝑅 ∙ 𝑟 ) = 𝑘𝑣
(𝑅 ∙ 𝑟 ∗ )2 𝑑(𝑅 ∙ 𝑟 ∗ ) 𝑑(𝑅 ∙ 𝑟 ∗ )

Resolviendo los paréntesis:

1 1 𝑠 𝑑 ∗2
𝑑𝐶𝐴∗
𝔇𝑒 𝐶𝐴 (𝑟 ) = 𝑘𝑣
𝑟∗2 𝑅2 𝑑𝑟 ∗ 𝑑𝑟 ∗

Reorganizando los términos para tener una expresión del módulo de Thiele:

1 𝑑 ∗2
𝑑𝐶𝐴∗ 𝑅 2 𝑘𝑣
(𝑟 ) =
𝑟∗2 𝑑𝑟 ∗ 𝑑𝑟 ∗ 𝔇𝑒 𝐶𝐴𝑠
𝑅 2 𝑘𝑣
El termino 𝔇 𝑠 , se puede expresar en función del módulo de Thiele:
𝑒 𝐶𝐴

1 𝑑 ∗2
𝑑𝐶𝐴∗
(𝑟 ) = 𝜙̅ 2
𝑟∗2 𝑑𝑟 ∗ 𝑑𝑟 ∗

Como el orden de la reacción es cero, la ecuación diferencial se puede resolver por variables
separables:

∗2
𝑑𝐶𝐴∗ 𝑑𝐶𝐴∗ 𝑟∗3
𝑑 (𝑟 ) = 𝜙̅ 2 𝑟 ∗ 2 𝑑𝑟 ∗ Integrando: 𝑟 ∗2 ̅
=𝜙 2
+ 𝐶1
𝑑𝑟 ∗ 𝑑𝑟 ∗ 3
Se reorganizan términos y se vuelve a integrar:

𝑟 ∗ 𝐶1 𝑟 ∗ 2 𝐶1
𝑑𝐶𝐴∗ = [𝜙̅ 2 + ] 𝑑𝑟 ∗ ⇒ 𝐶𝐴∗ = 𝜙̅ 2 + ∗ + 𝐶2
3 𝑟 ∗2 6 𝑟
Las nuevas condiciones de contorno adimensionales son:
𝑑𝐶𝐴∗
𝑟∗ = 0 =0
𝑑𝑟 ∗
𝑟∗ = 1 𝐶𝐴∗ = 1
Aplicando la primera condición de contorno:
𝑑𝐶𝐴∗ 𝑟∗ 𝐶
= ̅2 + 1
𝜙
𝑑𝑟 ∗ 3 𝑟 ∗2
(0) 𝐶1
0 = 𝜙̅ 2 + ⇒ 𝐶1 = 0
3 (0)2
Aplicando la condición de contorno dos:

𝑟∗2 (1)2 𝜙̅ 2
𝐶𝐴∗ = 𝜙̅ 2 + 𝐶2 ⇒ 1 = 𝜙̅ 2 + 𝐶2 ⇒ 𝐶2 = 1 −
6 6 6
Remplazando el valor de las constantes en la expresión de 𝐶𝐴∗ , se obtiene:

𝑟 ∗2 0 𝜙̅ 2 𝜙̅ 2 ∗ 2
𝐶𝐴∗ = 𝜙̅ 2 + ∗+1− ⇒ 𝐶𝐴∗ = (𝑟 − 1) + 1
6 𝑟 6 6
La anterior ecuación representa el perfil de concentración adimensional cuando la
concentración no es cero. Para validar eso, la concentración en radio cero debe ser mayor o
igual a cero, por tanto:

𝜙̅ 2 ∗ 2
(𝑟 − 1) + 1 ≥ 0 Por lo cual 𝜙̅ es: 𝜙̅ ≤ √6
6
Lo anterior es válido cuando 𝜙̅ ≤ √6, de lo contrario puede incidir en el tamaño de partícula.

Si para 𝜙̅ ≤ √6 es válido lo anterior, entonces debe existir un rango de 𝑟 ∗ donde la

concentración de A no sea cero antes de llegar al centro, por lo cual se plantea que dicho rango
es 𝑟𝑐 ≤ 𝑟 ∗ ≤ 1. Para hallar a 𝑟𝑐 se realiza el balance de masa y se establecen las condiciones
de contorno:

1 𝑑 ∗2
𝑑𝐶𝐴∗
(𝑟 ) = 𝜙̅ 2 𝑟𝑐 ≤ 𝑟 ∗ ≤ 1
𝑟∗2 𝑑𝑟 ∗ 𝑑𝑟 ∗

𝑟 ∗ 2 𝐶1
𝐶𝐴∗ = 𝜙̅ 2 + ∗ + 𝐶2 𝑟𝑐 ≤ 𝑟 ∗ ≤ 1
6 𝑟
Las nuevas condiciones de contorno adimensionales para 𝑟𝑐 son:
 𝑑𝐶𝐴∗
𝑟 ∗ = 𝑟𝑐 =0
𝑑𝑟 ∗
𝑟 ∗ = 𝑟𝑐 𝐶𝐴∗ = 0
Aplicando la primera condición de contorno:
𝑑𝐶𝐴∗ 𝑟∗ 𝐶
= ̅2 + 1
𝜙
𝑑𝑟 ∗ 3 𝑟 ∗2
(𝑟𝑐 ) 𝐶1 (𝑟𝑐 )3
0 = 𝜙̅ 2 + ⟹ 𝐶1 = −𝜙̅ 2
3 (𝑟𝑐 )2 3
Aplicando la condición de contorno dos, se integra y se obtiene el valor de 𝐶2 :

𝑑𝐶𝐴∗ 𝑟 ∗ 𝜙̅ 2 (𝑟𝑐 )3 𝑟 ∗ 2 𝜙̅ 2 (𝑟𝑐 )3

= ̅2 −
𝜙 ⇒ 𝐶𝐴∗ = 𝜙̅ 2 + + 𝐶2
𝑑𝑟 ∗ 3 3𝑟 ∗ 2 6 3𝑟 ∗
(𝑟𝑐 )2 𝜙̅ 2 (𝑟𝑐 )3 (𝑟𝑐 )2 𝜙̅ 2 (𝑟𝑐 )2
0 = 𝜙̅ 2 + + 𝐶2 ⇒ 0 = 𝜙̅ 2 + + 𝐶2
6 3(𝑟𝑐 ) 6 3

1 1 𝜙̅ 2 (𝑟𝑐 )2
0 = 𝜙̅ 2 (𝑟𝑐 )2 ( + ) + 𝐶2 ⇒ 𝐶2 = −
6 3 2
Remplazando el valor de las constantes en la expresión de 𝐶𝐴∗ , se obtiene:

𝑟 ∗ 2 𝜙̅ 2 (𝑟𝑐 )3 𝜙̅ 2 (𝑟𝑐 )2
𝐶𝐴∗ = 𝜙̅ 2 + −
6 3𝑟 ∗ 2
Para encontrar el valor de 𝑟𝑐 , se estable la condición de contorno:

𝐶𝐴∗ = 1 𝑟∗ = 1
Aplicando la condición de contorno:

𝜙̅ 2 𝜙̅ 2 (𝑟𝑐 )3 𝜙̅ 2 (𝑟𝑐 )2
1= + −
6 3 2
2(𝑟𝑐 )3 − 3(𝑟𝑐 )2 − 6𝜙̅ −2 + 1 = 0

Al graficar los valores de la concentración y el radio adimensional cuando 𝜙̅ < √6, 𝜙̅ = √6 y

con el valor calculado de 𝑟𝑐 , se obtiene:
 𝑟 ∗ 2 𝜙̅ 2 (𝑟𝑐 )3 𝜙̅ 2 (𝑟𝑐 )2
𝐶𝐴∗ = 𝜙̅ 2 + −
6 3𝑟 ∗ 2
El valor de 𝑟𝑐 para obtener 𝐶𝐴∗ debe ser de 3/2, esto condiciona que 𝑟 ∗ debe estar entre 0,8844
y 1, si no se cumple lo anterior la concentración será negativa.

Problema 3:

𝐴+𝐵 ↔𝑃
Mecanismo:
𝐶𝐴𝑆1
1. 𝐴 + 𝑆1 ↔ 𝐴𝑆1 𝑘𝐴 =
𝐶𝑆1 𝐶𝐴
𝐶𝐵𝑆2
2. 𝐵 + 𝑆2 ↔ 𝐵𝑆2 𝑘𝐵 =
𝐶𝑆2 𝐶𝐵
𝐶𝑃𝑆1 𝐶𝑆2
3. 𝐴𝑆1 + 𝐵𝑆2 → 𝑃𝑆1 + 𝑆2 𝑘𝑠 → ∞ =
𝐶𝐴𝑆1 𝐶𝐵𝑆2
𝐶𝑃𝑆1
4. 𝑃𝑆1 ↔ 𝑃 + 𝑆1 𝑘𝑃 =
𝐶𝑆1 𝐶𝑃

Balance de sitios:

𝐶𝑇𝑆 = 𝐶𝑆1 + 𝐶𝐴𝑆1 + 𝐶𝑃𝑆1

1

𝐶𝑇𝑆2 = 𝐶𝑆2 + 𝐶𝐵𝑆2

Parte 1.

Si la adsorción de A es la etapa controlante:

𝐴 + 𝑆1 ↔ 𝐴𝑆1 𝑟𝐴1 = 𝑘1 𝐶𝑆1 𝐶𝐴 − 𝑘1−1 𝐶𝐴𝑆1

Se despejas las concentraciones desconocidas y se remplazan en el balance de sitios 𝐶𝑇𝑆 :

1

𝐶𝐵𝑆2
𝑘𝐵 = → 𝐶𝐵𝑆2 = 𝑘𝐵 𝐶𝑆2 𝐶𝐵
𝐶𝑆2 𝐶𝐵
𝐶𝑃𝑆1 𝐶𝑆2 𝐶𝑃𝑆1 𝐶𝑆2
𝑘𝑠 = → 𝐶𝐴𝑆1 = ≈0 𝑘𝑠 → ∞
𝐶𝐴𝑆1 𝐶𝐵𝑆2 𝑘𝑠 𝐶𝐵𝑆2
𝐶𝑃𝑆1
𝑘𝑃 = → 𝐶𝑃𝑆1 = 𝐶𝑆1 𝐶𝑃 𝑘𝑃
𝐶𝑆1 𝐶𝑃

Reemplazamos:

𝐶𝑇𝑆1 = 𝐶𝑆1 + 𝐶𝑃𝑆1 → 𝐶𝑇𝑆1 = 𝐶𝑆1 + 𝐶𝑆1 𝐶𝑃 𝑘𝑃

𝐶𝑇𝑆1
𝐶𝑇𝑆1 = 𝐶𝑆1 (1 + 𝐶𝑃 𝑘𝑃 ) → 𝐶𝑆1 =
1 + 𝐶𝑃 𝑘𝑃
Reemplazado en la velocidad de reacción cuando 𝑘𝑠 → ∞:
𝑘1 𝐶𝐴 𝐶𝑇𝑆1 𝑘1∗ 𝐶𝐴
𝑟𝐴1 = → 𝑟𝐴1 =
1 + 𝐶𝑃 𝑘𝑃 1 + 𝐶𝑃 𝑘𝑃
Si la etapa controlante es la reacción B:

𝐵 + 𝑆2 ↔ 𝐵𝑆2 𝑟𝐴2 = 𝑘2 𝐶𝑆2 𝐶𝐵 − 𝑘2−1 𝐶𝐵𝑆2

Se despejas las concentraciones desconocidas y se remplazan en el balance de sitios 𝐶𝑇𝑆2 :

𝐶𝐴𝑆1
𝑘𝐴 = → 𝐶𝐴𝑆1 = 𝑘𝐴 𝐶𝑆1 𝐶𝐴
𝐶𝑆1 𝐶𝐴
𝐶𝑃𝑆1 𝐶𝑆2 𝐶𝑃𝑆1 𝐶𝑆2
𝑘𝑠 = → 𝐶𝐵𝑆2 = ≈0 𝑘𝑠 → ∞
𝐶𝐴𝑆1 𝐶𝐵𝑆2 𝑘𝑠 𝐶𝐴𝑆1

Reemplazamos:

𝐶𝑇𝑆2 = 𝐶𝑆2 + 𝐶𝐵𝑆2 → 𝐶𝑇𝑆2 = 𝐶𝑆2

Reemplazado en la velocidad de reacción cuando 𝑘𝑠 → ∞:

𝑟𝐴2 = 𝑘2 𝐶𝐵 𝐶𝑇𝑆2 → 𝑟𝐴2 = 𝑘2∗ 𝐶𝐵

Si la etapa controlante es la reacción superficial:

𝐴𝑆1 + 𝐵𝑆2 → 𝑃𝑆1 + 𝑆2 𝑟𝐴3 = 𝑘3 𝐶𝐴𝑆1 𝐶𝐵𝑆2

Se despejas las concentraciones desconocidas y se remplazan en el balance de sitios 𝐶𝑇𝑆1 y
𝐶𝑇𝑆2 :

𝐶𝐴𝑆1
𝑘𝐴 = → 𝐶𝐴𝑆1 = 𝑘𝐴 𝐶𝑆1 𝐶𝐴
𝐶𝑆1 𝐶𝐴
𝐶𝐵𝑆2
𝑘𝐵 = → 𝐶𝐵𝑆2 = 𝑘𝐵 𝐶𝑆2 𝐶𝐵
𝐶𝑆2 𝐶𝐵

Reemplazamos:

Para 𝐶𝑇𝑆 ;
1

𝐶𝑇𝑆1 = 𝐶𝑆1 + 𝐶𝐴𝑆1 + 𝐶𝑃𝑆1 → 𝐶𝑇𝑆1 = 𝐶𝑆1 + 𝑘𝐴 𝐶𝑆1 𝐶𝐴 + 𝐶𝑆1 𝐶𝑃 𝑘𝑃

𝐶𝑇𝑆1
𝐶𝑇𝑆1 = 𝐶𝑆1 (1 + 𝑘𝐴 𝐶𝐴 + 𝐶𝑃 𝑘𝑃 ) → 𝐶𝑆1 =
1 + 𝑘𝐴 𝐶𝐴 + 𝐶𝑃 𝑘𝑃
Para 𝐶𝑇𝑆2 ;

𝐶𝑇𝑆 = 𝐶𝑆2 + 𝐶𝐵𝑆2 → 𝐶𝑇𝑆 = 𝐶𝑆2 + 𝑘𝐵 𝐶𝑆2 𝐶𝐵

2 2

𝐶𝑇𝑆2
𝐶𝑇𝑆2 = 𝐶𝑆2 (1 + 𝑘𝐵 𝐶𝐵 ) → 𝐶𝑆2 =
1 + 𝑘𝐵 𝐶𝐵
Reemplazado en la velocidad de reacción cuando:

𝑟𝐴3 = 𝑘3 𝐶𝐴𝑆1 𝐶𝐵𝑆2 → 𝑟𝐴3 = 𝑘3 𝑘𝐴 𝐶𝑆1 𝐶𝐴 𝑘𝐵 𝐶𝑆2 𝐶𝐵

𝑘3 𝐶𝐴 𝐶𝐵 𝑘𝐴 𝑘𝐵 𝐶𝑇𝑆1 𝐶𝑇𝑆2 𝑘3∗ 𝐶𝐴 𝐶𝐵
𝑟𝐴3 = → 𝑟𝐴3 =
[1 + 𝑘𝐴 𝐶𝐴 + 𝐶𝑃 𝑘𝑃 ][1 + 𝑘𝐵 𝐶𝐵 ] [1 + 𝑘𝐴 𝐶𝐴 + 𝐶𝑃 𝑘𝑃 ][1 + 𝑘𝐵 𝐶𝐵 ]
Si la etapa controlante es la obtención del producto P:

𝑃𝑆1 ↔ 𝑃 + 𝑆1 𝑟𝐴4 = 𝑘4 𝐶𝑃𝑆1 − 𝑘4−1 𝐶𝑆1 𝐶𝑃

De los valores previamente obtenidos cuando se asumieron las demás etapas controlantes,
se dijo que 𝐶𝐴𝑆1 = 0, lo cual conlleva a que 𝐶𝑆1 = 0, por lo que el balance de sitios para 𝐶𝑇𝑆1 :

𝐶𝑃𝑆1 ≈ 𝐶𝑇𝑆1

Para que 𝐶𝑇𝑆1 ≈ 𝐶𝑃𝑆1 , debe cumplirse que 𝐶𝐴𝑆1 = 0 ∧ 𝐶𝑆1 = 0

Lo que significa que:

𝑟𝐴4 = 𝑘4 𝐶𝑇𝑆1 = 𝑘4∗

Parte 2 y 3.

Analizando la suposición de etapa controlante la adsorción de A:

Los datos iniciales plantean que pueden presentarse cambios en la velocidad de reacción al
cambiar la concentración de A, se grafican los datos experimentales, se calculan los
parámetros y se recalculan las velocidades de reacción para hacer el ajuste, lo cual se
evidencia a continuación:

Para lo cual se obtienen los siguientes valores:

𝑘1∗ = 0,2567 𝑘𝑃 = 19,3057 𝐸𝑟𝑟𝑜𝑟𝐶𝑢𝑎𝑑𝑟𝑎 = 3,93𝑥10−6

Analizando la adsorción de B:

El modelo para cuando la adsorción de B es la etapa controlante, propone que la velocidad de

reacción solo dependa de la concentración de B, pero se evidencian cambios de la velocidad
de reacción inicial cuando la concentración de B se mantiene constante, por lo que se puede
despreciar este modelo.

Analizando la reacción superficial:

Los datos iniciales plantean que pueden presentarse cambios en la velocidad de reacción al
cambiar la concentración A, B y P; se grafican los datos experimentales, se calculan los
parámetros y se recalculan las velocidades de reacción para hacer el ajuste, lo cual se
evidencia a continuación:
 𝑘3∗ = 1,1089 𝑘𝐴 = 0,5445 𝑘𝐵 = 2,0573 𝑘𝑃 = 15,1798 𝐸𝑟𝑟𝑜𝑟𝐶𝑢𝑎𝑑𝑟𝑎 = 3,97𝑥10−6

El modelo cuando la etapa controlante es la reacción del producto plantea que la velocidad de
reacción es igual a la constante sin importar la concentración de la especie, cuando se sabe
de los datos que esta velocidad si varia producto de la conversión de A, por lo cual no es tenido
en cuenta dicho modelo.

Analizando los datos del modelo 1 y 3, se puede apreciar que los errores cuadráticos son muy
parecidos, pero los ajusten no tienen la misma precisión

Problema 4:

𝐴→𝐵 𝑟𝐴 = 𝑘 ∗ 𝐶𝐴

0 L e+L e+2L

Condiciones de contorno:

𝐶𝐴(2𝐿+𝑒) = 𝐶𝐴𝑠
𝑑𝐶𝐴
| =0 𝑝𝑜𝑟 𝑠𝑖𝑚𝑒𝑡𝑟𝑖𝑎 𝑒𝑛 𝑒𝑙 𝑐𝑒𝑛𝑡𝑟𝑜 𝑑𝑒 𝑙𝑎 𝑙á𝑚𝑖𝑛𝑎
𝑑𝑧 𝑧=0
𝑑2 𝐶𝑠
𝔇𝑒 − 𝑘𝑣 𝐶𝑠 = 0 → 𝑠𝑒 𝑐𝑢𝑚𝑝𝑙𝑒
𝑑𝑧 2
 𝑘𝑣
𝐶𝑠 (𝑧) cosh √ ⁄𝔇𝑒 𝑧
= ← 𝑠𝑒 𝑐𝑢𝑚𝑝𝑙𝑒 𝑒𝑛 𝐿 ≤ 𝑥 ≤ 𝑒 + 𝐿
𝐶𝑠𝑠 𝑘
cosh √ 𝑣⁄𝔇 𝐿
𝑒

En el material inerte solamente hay difusión por lo que el perfil va a ser lineal, pero en el
material catalítico el perfil es pronunciado.
𝑑𝐶𝐴
0≤𝑥≤𝐿 𝑥=0 ; | =0
𝑑𝑥 𝑥=0
Integrando:

𝑑2 𝐶𝐴
𝔇𝑒 =0 → 𝐶𝐴 = 𝐶1 𝑥 + 𝐶2
𝑑𝑥 2
Aplicando la condición de borde:

𝐶1 = 0 → 𝐶2 = 𝐶𝐴 (1)
Evaluando en L:
𝑑𝐶𝐴
𝑁𝐴 = −𝔇𝑒 | =?
𝑑𝑥 𝑥=𝐿
𝑑𝐶𝐴 𝑑(𝐶𝑧 )
| =0 𝑝𝑜𝑟𝑞𝑢𝑒 =0
𝑑𝑥 𝑥=𝐿 𝑑𝑥
En la parte catalítica:

𝑑2 𝐶𝐴 𝑑2 𝐶𝐴 𝑘𝑣 𝐶𝐴
𝔇𝑒 − 𝑘𝑣 𝐶𝐴 = 0 → − =0
𝑑𝑥 2 𝑑𝑥 2 𝔇𝑒
Integrando:

𝑘𝑣 𝑘𝑣
𝑟2 = → 𝑟 = ±√
𝔇𝑒 𝔇𝑒

Entonces:
𝑘 𝑘
√− 𝑣 𝑥 √ 𝑣𝑥
𝐶𝐴 = 𝐶3 ∙ 𝑒 𝔇𝑒 + 𝐶4 ∙ 𝑒 𝔇𝑒

Lo que equivale a:

𝑘𝑣 𝑘𝑣
𝐶𝐴 = 𝐶3 ∙ cosh (√ 𝑥) + 𝐶4 ∙ sinh (√ 𝑥)
𝔇𝑒 𝔇𝑒

Si x se mueve a (x-L):
 𝑘𝑣 𝑘𝑣
𝐶𝐴 = 𝐶3 ∙ cosh (√ (𝑥 − 𝐿)) + 𝐶4 ∙ sinh (√ (𝑥 − 𝐿))
𝔇𝑒 𝔇𝑒

𝑑𝐶𝐴 𝑘𝑣 𝑘𝑣 𝑘𝑣 𝑘𝑣
= 𝐶3 ∙ sinh (√ (𝑥 − 𝐿)) ∙ √ + √ ∙ 𝐶4 ∙ cosh (√ (𝑥 − 𝐿))
𝑑𝑥 𝔇𝑒 𝔇𝑒 𝔇𝑒 𝔇𝑒

𝑑𝐶𝐴 𝑘𝑣 𝑘𝑣
| = 𝐶4 ∙ √ ∙ cosh(0) = 𝐶4 ∙ √
𝑑𝑥 𝑥=𝐿 𝔇𝑒 𝔇𝑒

Sabiendo que:

𝑑𝐶𝐴 𝑘𝑣
| =0 → 𝐶4 ∙ √ =0 → 𝐶4 = 0
𝑑𝑥 𝑥=𝐿 𝔇𝑒

Por lo que:

𝑘𝑣
𝐶𝐴 = 𝐶3 ∙ cosh (√ (𝑥 − 𝐿)) (2)
𝔇𝑒

Evaluando 𝐶𝐴 𝑒𝑛 𝑥 = 𝐿

𝐶𝐴 = 𝐶3 ∙ cosh(0) → 𝐶𝐴 (𝑥 = 𝐿) = 𝐶3 = 𝐶2
Evaluando en 𝑒 + 𝐿 ≤ 𝑥 ≤ 𝑒 + 2𝐿

𝑑2 𝐶𝐴
𝔇𝑒 =0 → 𝐶𝐴 = 𝐶5 𝑥 + 𝐶6
𝑑𝑥 2
𝐶𝐴 (𝑥 = 2𝐿) = 𝐶𝐴𝑠
𝐶𝐴𝑠 = 𝐶5 (𝑒 + 2𝐿) + 𝐶6 → 𝐶6 = 𝐶𝐴𝑠 − 𝐶5 (𝑒 + 2𝐿)

𝐶𝐴 = 𝐶5 𝑥 + 𝐶𝐴𝑠 − 𝐶5 (𝑒 + 2𝐿) (3) 𝑣á𝑙𝑖𝑑𝑜 𝑝𝑎𝑟𝑎 𝑒𝑙 𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑜 𝑒 + 𝐿 ≤ 𝑥 ≤ 𝑒 + 2𝐿

Como 𝐶𝐴 en 𝐿 ≤ 𝑥 ≤ 𝑒 + 𝐿

𝑘𝑣
𝐶𝐴 = 𝐶3 ∙ cosh (√ (𝑥 − 𝐿))
𝔇𝑒

Evaluamos en 𝑒 + 𝐿

𝑘𝑣
𝐶𝐴 (𝑥=𝑒+𝐿) = 𝐶3 ∙ cosh (√ (𝑒 + 𝐿 − 𝐿)) (2)
𝔇𝑒
 𝑘𝑣
𝐶𝐴 (𝑥=𝑒+𝐿) = 𝐶3 ∙ cosh (√ 𝑒)
𝔇𝑒

𝐶𝐴 (𝑥=𝑒+𝐿) = 𝐶5 (𝑒 + 𝐿) + 𝐶𝐴𝑠 − 𝐶5 (𝑒 + 2𝐿) (3)

𝐶𝐴 (𝑥=𝑒+𝐿) = 𝐶𝐴𝑠 − 𝐶5 (𝐿)

Como el perfil debe se continuo en ambos lados:

𝑘𝑣 𝑘𝑣 1 𝐶𝐴𝑠
𝐶𝐴𝑠 − 𝐶5 (𝐿) = 𝐶3 ∙ cosh (√ 𝑒) ⟹ 𝐶5 = −𝐶3 ∙ cosh (√ 𝑒) + (4)
𝔇𝑒 𝔇𝑒 𝐿 𝐿

Teniendo en cuenta que el flujo es el mismo en ambos lados:

𝐶𝐴 = 𝐶5 𝑥 + 𝐶6
𝑑𝐶𝐴 𝑑𝐶𝐴
= 𝐶5 = |
𝑑𝑥 𝑑𝑥 𝑥=𝑒+𝐿

𝑘𝑣 𝑑𝐶𝐴 𝑘𝑣 𝑘𝑣
𝐶𝐴 = 𝐶3 ∙ cosh (√ (𝑥 − 𝐿)) ⟹ = 𝐶3 ∙ √ (sinh (√ (𝑥 − 𝐿)))
𝔇𝑒 𝑑𝑥 𝔇𝑒 𝔇𝑒

𝑑𝐶𝐴 𝑘𝑣 𝑘𝑣 𝑘𝑣 𝑘𝑣
| = 𝐶3 ∙ √ ∙ sinh (√ 𝑒) ⟹ 𝐶5 = 𝐶3 ∙ √ ∙ sinh (√ 𝑒) (5)
𝑑𝑥 𝑥=𝑒+𝐿 𝔇𝑒 𝔇𝑒 𝔇𝑒 𝔇𝑒

Igualando (4) y (5):

𝑘𝑣 𝑘𝑣 𝑘𝑣 1 𝐶𝐴𝑠
𝐶3 ∙ √ ∙ sinh (√ 𝑒) = −𝐶3 ∙ cosh (√ 𝑒) +
𝔇𝑒 𝔇𝑒 𝔇𝑒 𝐿 𝐿

𝑘𝑣 𝑘𝑣 𝑘𝑣
𝐿𝐶3 ∙ √ ∙ sinh (√ 𝑒) + 𝐶3 ∙ cosh (√ 𝑒) = 𝐶𝐴𝑠
𝔇𝑒 𝔇𝑒 𝔇𝑒

𝑘𝑣 𝑘𝑣 𝑘𝑣
𝐶3 (𝐿 ∙ √ ∙ sinh (√ 𝑒) + cosh (√ 𝑒)) = 𝐶𝐴𝑠
𝔇𝑒 𝔇𝑒 𝔇𝑒

𝐶𝐴𝑠
𝐶3 =
𝑘 𝑘 𝑘
(𝐿 ∙ √𝔇𝑣 ∙ sinh (√𝔇𝑣 𝑒) + cosh (√𝔇𝑣 𝑒))
𝑒 𝑒 𝑒
Problema 5:
𝑟𝑚 ∙ 𝐶𝑠𝑏
⁄(𝐾 + 𝐶 )
𝐾 𝑀 𝑠𝑏
𝑘 = 𝑀⁄𝐶 𝜙𝑒 =
𝑠𝑏 𝑘𝑐 𝑎 ∙ 𝐶𝑠𝑏
El factor de efectividad externo se define como:
(𝑟𝑚 ∙ 𝐶𝑠𝑠 )
𝑟𝑠𝑠 ⁄(𝐾 + 𝐶 )
𝑀 𝑠𝑠
𝜂𝑒𝑥𝑡 = =
𝑟 (𝑟𝑚 ∙ 𝐶 )
𝑠𝑏 ⁄
(𝐾𝑀+ 𝐶𝑠𝑏 )
Reorganizando los términos:
(𝑟𝑚 ∙ 𝐶𝑠𝑠 )(𝐾𝑀+ 𝐶𝑠𝑏 ) 𝐶𝑠𝑠 (𝐾𝑀+ 𝐶𝑠𝑏 )
𝜂𝑒𝑥𝑡 = 𝜂𝑒𝑥𝑡 =
(𝑟𝑚 ∙ 𝐶𝑠𝑏 )(𝐾𝑀 + 𝐶𝑠𝑠 ) 𝐶𝑠𝑏 (𝐾𝑀 + 𝐶𝑠𝑠 )
Se sabe que:
𝐶𝑠𝑠
𝐾𝑀 = 𝑘 ∙ 𝐶𝑠𝑏 𝛽= ⁄𝐶
𝑠𝑏

Se puede definir el factor de efectividad de la siguiente manera:

𝛽 ∙ 𝐶𝑠𝑏 (𝑘 + 1) 𝐶𝑠𝑏 𝛽(𝑘 + 1)
𝜂𝑒𝑥𝑡 = ∙ =
𝑘 ∙ 𝐶𝑠𝑏 + 𝐶𝑠𝑠 𝐶𝑠𝑏 (𝑘 + 𝛽)
Donde:
𝛽 1−𝛽
=
𝑘+𝛽 𝑁𝐷𝑎
Por lo tanto:
(1 − 𝛽)(𝑘 + 1)
𝜂𝑒𝑥𝑡 =
𝑁𝐷𝑎
El número de Damköhler se define como:
𝑟𝑚
𝑁𝐷𝑎 =
𝑘𝑐 𝑎 ∙ 𝐶𝑠𝑏
Relacionándolo con el módulo de Thiele:
𝑟𝑚 ∙ 𝐶𝑠𝑏
⁄(𝐾 + 𝐶 )
𝑀 𝑠𝑏
𝜙𝑒 =
𝑘𝑐 𝑎 ∙ 𝐶𝑠𝑏
𝑟𝑚 𝐶𝑠𝑏 𝑟𝑚 𝐶𝑠𝑏
𝜙𝑒 = ∙ = ∙
𝑘𝑐 𝑎 ∙ 𝐶𝑠𝑏 (𝐾𝑀 + 𝐶𝑠𝑏 ) (𝑘
𝑘𝑐 𝑎 ∙ 𝐶𝑠𝑏 𝐶𝑠𝑏 + 1)
𝑁𝐷𝑎
𝜙𝑒 = → 𝑁𝐷𝐴 = 𝜙𝑒 (𝑘 + 1)
(𝑘 + 1)
Finalmente se llega a la expresión del factor de efectividad externo:
(1 − 𝛽)(𝑘 + 1) (1 − 𝛽)
𝜂𝑒𝑥𝑡 = =
𝜙𝑒 (𝑘 + 1) 𝜙𝑒
Para realizar las gráficas, se parte de:
1−𝛽 𝛽 1−𝛽 𝛽
= ⇒ =
𝑁𝐷𝑎 𝑘 + 𝛽 𝜙𝑒 (𝑘 + 1) 𝑘 + 𝛽

(1 − 𝛽)(𝑘 + 𝛽) = 𝛽(𝜙𝑒 (𝑘 + 1)) ⟹ 𝑘 + 𝛽 − 𝑘𝛽 − 𝛽 2 = 𝛽𝜙𝑒 (𝑘 + 1)

𝑘 + 𝛽(1 − 𝑘) − 𝛽 2 = 𝛽𝜙𝑒 (𝑘 + 1) ⟹ 𝑘 + 𝛽(1 − 𝑘) − 𝛽 2 − 𝛽𝜙𝑒 (𝑘 + 1) = 0

𝛽 2 + [𝜙𝑒 (𝑘 + 1) + (𝑘 − 1)]𝛽 − 𝑘 = 0
Aplicando la ecuación de la cuadrática:

−(𝜙𝑒 (𝑘 + 1) + (𝑘 − 1)) ± √[𝜙𝑒 (𝑘 + 1) + (𝑘 − 1)]2 + 4𝑘 (1 − 𝛽)

𝛽= 𝜂𝑒𝑥𝑡 =
2 𝜙𝑒

−( 𝑁𝐷𝐴 + (𝑘 − 1)) ± √[ 𝑁𝐷𝐴 + (𝑘 − 1)]2 + 4𝑘 (1 − 𝛽)(𝑘 + 1)

𝛽= 𝜂𝑒𝑥𝑡 =
2 𝑁𝐷𝐴
Se grafica el factor de efectividad externo vs. el módulo de Thiele externo, cuando varia beta en
función de k.

En otra grafica se relaciona el factor de efectividad vs. el número de Damköhler, cuando beta
varía en función de k.

A continuación, se presentan graficas en escala normal y escala logarítmica, como las arroja
Matlab.
A valores mayores de 𝑁𝐷𝐴 𝑦 𝜙, la efectividad disminuye, provocando que la velocidad baje en
comparación al modulo de Thiele.

Problema 6:

𝑒𝑥𝑡 𝑖𝑛𝑡
𝑆𝑒𝑥𝑡 = 𝑟𝑀−𝑀 𝑆𝑖𝑛𝑡 = 𝑟𝑀−𝑀
Como:
𝑟𝑟𝑒𝑎𝑙 𝑟𝑖𝑛𝑡 𝑖𝑛𝑡 𝑒𝑥𝑡
𝜂𝑖𝑛𝑡 = ≈ ⟹ 𝑟𝑀−𝑀 = 𝜂𝑖𝑛𝑡 ∙ 𝑟𝑀−𝑀
𝑟𝑖𝑑𝑒𝑎𝑙 𝑟𝑒𝑥𝑡
En ausencia de limitaciones difusivas:
𝑒𝑥𝑡
𝑆𝑒𝑥𝑡 : 𝑟𝑀−𝑀
𝑒𝑥𝑡 𝑖𝑛𝑡
𝑆𝑖𝑛𝑡 : 𝑟𝑀−𝑀 = 𝑟𝑀−𝑀
Para definir el factor de efectividad global
𝑒𝑥𝑡 𝑒𝑥𝑡
𝐺
𝑣𝑒𝑙. 𝑟𝑥𝑛 𝑟𝑒𝑎𝑙 𝑟𝑀−𝑀 ∙ 𝑆𝑒𝑥𝑡 + 𝜂𝑖𝑛𝑡 ∙ 𝑟𝑀−𝑀 ∙ 𝑆𝑖𝑛𝑡
𝜂 = = 𝑒𝑥𝑡
𝑣𝑒𝑙. 𝑟𝑥𝑛 𝑖𝑑𝑒𝑎𝑙 𝑟𝑀−𝑀 ∙ (𝑆𝑒𝑥𝑡 + 𝑆𝑖𝑛𝑡 )
Por tanto, el factor de efectividad global es:
𝑆𝑒𝑥𝑡
⁄𝑆 + 𝜂𝑖𝑛𝑡 𝑎𝑟 + 𝜂𝑖𝑛𝑡
𝑖𝑛𝑡
𝜂𝐺 = ⟹ 𝜂𝐺 =
𝑆𝑒𝑥𝑡 𝑎𝑟 + 1
⁄𝑆 + 1
𝑖𝑛𝑡

Cuando 𝑆𝑒𝑥𝑡 ≪ 𝑆𝑖𝑛𝑡 𝑎𝑟 ≈ 0 → 𝜂 𝐺 ≈ 𝜂𝑖𝑛𝑡

𝜂 𝐺 =?
𝜂𝑖𝑛𝑡 =? → 𝑠𝑒 ℎ𝑎𝑙𝑙𝑎 𝑐𝑜𝑛 𝜙

𝑟𝑚 8𝑥102 𝑚𝑁⁄𝑠
( ) ( ⁄ )
𝑅𝑜 √ ⁄𝐾𝑀 0,2 𝑐𝑚 √ 5𝑥10−2 𝑚𝑁⁄𝑠
𝜙= = 2 = 84,3274
3 𝐷𝑒𝑓𝑓 3 0,01 𝑐𝑚 ⁄𝑠

Como 𝜙 ≫ 3
1 1+𝛽 1 𝐶𝑠𝑠 5 𝑚𝑀
𝜂𝜙 ≫ 3 = ∙ ∙ √2 ∙ (𝛽 − ln(1 + 𝛽)) ⁄2 ⟹ 𝛽= = = 100
𝜙 𝛽 𝐾𝑀 5𝑥10−2 𝑚𝑀
1 1 + 100 1
𝜂𝜙 ≫ 3 = ∙ ∙ √2 ∙ (100 − ln(1 + 100)) ⁄2
84,32 100
0,7 + 0,1654
𝜂𝜙 ≫ 3 = 0,1654 ⟹ 𝜂𝐺 = = 0,509
0,7 + 1

Problema 7:

2
𝐶𝐻𝑆
1. 𝐻2 + 2𝑆 ↔ 2𝐻𝑆 𝑘𝐻2 =
𝑃𝐻2 𝐶𝑆2
𝐶𝐴𝑆
2. 𝐴 + 𝑆 ↔ 𝐴𝑆 𝑘𝐴 =
𝐶𝐴 𝐶𝑆
𝐶𝐴𝐻𝑆 𝐶𝑆
3. 𝐻𝑆 + 𝐴𝑆 ↔ 𝐴𝐻𝑆 + 𝑆 𝑘𝑆1 =
𝐶𝐻𝑆 𝐶𝐴𝑆
 𝐶𝑆 𝐶𝐵𝑆
4. 𝐻𝑆 + 𝐴𝐻𝑆 ↔ 𝐵𝑆 + 𝑆 𝑘𝑆2 =
𝐶𝐻𝑆 𝐶𝐴𝐻𝑆
𝐶𝐵𝑆
5. 𝐵𝑆 ↔ 𝐵 + 𝑆 𝑘𝐵 =
𝐶𝐵 𝐶𝑆
Balance de sitios:

𝐶𝑇𝑆 = 𝐶𝑆 + 𝐶𝐻𝑆 + 𝐶𝐴𝑆 + 𝐶𝐴𝐻𝑆 + 𝐶𝐵𝑆

Asumiendo que la adsorción de hidrógeno es la etapa controlante:

1 2
𝐻2 + 2𝑆 ↔ 2𝐻𝑆 𝑟 = 𝑘𝑑 𝑃𝐻2 𝐶𝑆2 − 𝑘𝑑−1 𝐶𝐻𝑆
2
= 𝑘𝑑 (𝑃𝐻2 𝐶𝑆2 − 𝐶 )
𝑘𝐻𝑆 𝐻𝑆
Se despejas las concentraciones desconocidas y se remplazan en el balance de sitios 𝐶𝑇𝑆 :

𝐶𝐴𝑆
𝑘𝐴 = → 𝐶𝐴𝑆 = 𝑘𝐴 𝐶𝐴 𝐶𝑆 𝐸𝑐. 1
𝐶𝐴 𝐶𝑆
𝐶𝐴𝐻𝑆 𝐶𝑆 𝐶𝐴𝐻𝑆 𝐶𝑆 𝐶𝐴𝐻𝑆 𝐶𝑆 𝐶𝐴𝐻𝑆
𝑘𝑆1 = → 𝐶𝐻𝑆 = = = 𝐸𝑐. 2
𝐶𝐻𝑆 𝐶𝐴𝑆 𝑘𝑆1 𝐶𝐴𝑆 𝑘𝑆1 𝑘𝐴 𝐶𝐴 𝐶𝑆 𝑘𝑆1 𝑘𝐴 𝐶𝐴

𝐶𝑆 𝐶𝐵𝑆 𝐶𝐵𝑆 𝐶𝑆 𝐶𝐵𝑆 𝐶𝑆 𝑘𝑆 𝑘𝐴 𝐶𝐴 𝐶𝐵𝑆 𝐶𝑆

𝑘𝑆2 = → 𝐶𝐴𝐻𝑆 = = → 𝐶𝐴𝐻𝑆 = √ 1 𝐸𝑐. 3
𝐶𝐻𝑆 𝐶𝐴𝐻𝑆 𝑘𝑆2 𝐶𝐻𝑆 𝑘 𝐶𝐴𝐻𝑆 𝑘𝑆2
𝑆2 𝑘 𝑘 𝐶
𝑆1 𝐴 𝐴

𝐶𝐵𝑆
𝑘𝐵 = → 𝐶𝐵𝑆 = 𝑘𝐵 𝐶𝐵 𝐶𝑆 𝐸𝑐. 4
𝐶𝐵 𝐶𝑆
Sustituyendo la ecuación 4 en la ecuación 3, se tiene:

𝑘𝑆 𝑘𝐴 𝐶𝐴 𝑘𝐵 𝐶𝐵 𝐶𝑆 𝐶𝑆 𝑘𝑆 𝑘𝐴 𝑘𝐵
𝐶𝐴𝐻𝑆 = √ 1 =√ 1 𝐶𝐴 𝐶𝐵 ∙ 𝐶𝑆 𝐸𝑐. 5
𝑘𝑆2 𝑘𝑆2

Sustituyendo la ecuación 5 en la ecuación 2, se llega:

1 𝑘𝑆 𝑘𝐴 𝑘𝐵 𝑘𝐵 𝐶𝐵
𝐶𝐻𝑆 = ∙√ 1 𝐶𝐴 𝐶𝐵 ∙ 𝐶𝑆 = √ ∙ 𝐶𝑆 𝐸𝑐. 6
𝑘𝑆1 𝑘𝐴 𝐶𝐴 𝐾𝑆2 𝑘𝑆2 𝑘𝑆1 𝑘𝐴 𝐶𝐴

Reemplazando las ecuaciones 6, 1, 5 y 4 en el balance de sitios totales y despejando 𝐶𝑆 :

𝐶𝑇𝑆 = 𝐶𝑆 + 𝐶𝐻𝑆 + 𝐶𝐴𝑆 + 𝐶𝐴𝐻𝑆 + 𝐶𝐵𝑆

𝑘𝐵 𝐶𝐵 𝑘𝑆 𝑘𝐴 𝑘𝐵
𝐶𝑇𝑆 = 𝐶𝑆 + √ ∙ 𝐶𝑆 + 𝑘𝐴 𝐶𝐴 𝐶𝑆 + √ 1 𝐶𝐴 𝐶𝐵 ∙ 𝐶𝑆 + 𝑘𝐵 𝐶𝐵 𝐶𝑠
𝑘𝑆2 𝑘𝑆1 𝑘𝐴 𝐶𝐴 𝑘𝑆2
 𝐶𝑇𝑆
𝐶𝑆 = 𝐸𝑐. 7
𝑘𝑆 𝑘𝐴 𝑘𝐵 𝑘𝐵 𝐶𝐵
1+√ 1 𝐶𝐴 𝐶𝐵 + √ + 𝑘𝐴 𝐶𝐴 + 𝑘𝐵 𝐶𝐵
𝑘𝑆2 𝑘𝑆2 𝑘𝑆1 𝑘𝐴 𝐶𝐴

Reemplazando la ecuación 6 en la expresión de velocidad de reacción:

2
1 𝑘𝐵 𝐶𝐵 1 𝑘𝐵 𝐶𝐵 2
𝑟 = 𝑘𝑑 (𝑃𝐻2 𝐶𝑆2 − (√ ∙ 𝐶𝑆 ) ) 𝑟 = 𝑘𝑑 (𝑃𝐻2 − )𝐶 𝐸𝑐. 8
𝑘𝐻𝑆 𝑘𝑆2 𝑘𝑆1 𝑘𝐴 𝐶𝐴 𝑘𝐻𝑆 𝑘𝑆2 𝑘𝑆1 𝑘𝐴 𝐶𝐴 𝑆

Reemplazando la ecuación 7 en la 8 se tiene:

2

1 𝑘𝐵 𝐶𝐵 𝐶𝑇𝑆
𝑟 = 𝑘𝑑 (𝑃𝐻2 − )
𝑘𝐻𝑆 𝑘𝑆2 𝑘𝑆1 𝑘𝐴 𝐶𝐴
𝑘𝑆 𝑘𝐴 𝑘𝐵 𝑘𝐵 𝐶𝐵
1+√ 1 𝐶𝐴 𝐶𝐵 + √ + 𝑘𝐴 𝐶𝐴 + 𝑘𝐵 𝐶𝐵
𝑘𝑆2 𝑘𝑆2 𝑘𝑆1 𝑘𝐴 𝐶𝐴
[ ]
De donde:

𝑘𝐻𝑆 𝑘𝑆2 𝑘𝑆1 𝑘𝐴 𝑘𝑆 𝑘𝐴 𝑘𝐵 𝑘𝐵

𝐾1 = 𝐶𝑇2𝑆 𝑘𝑑 𝐾𝑒𝑞 = 𝐾2 = √ 1 𝐾3 = √
𝑘𝐵 𝑘𝑆2 𝑘𝑆2 𝑘𝑆1 𝑘𝐴

Por lo tanto, se llega a la expresión final de la velocidad de reacción:

1 𝐶
𝐾1 (𝑃𝐻2 − 𝐾 𝐶𝐵 )
𝑒𝑞 𝐴
𝑟= 2
𝐶
(1 + 𝐾2 √𝐶𝐴 𝐶𝐵 + 𝐾3 √ 𝐶𝐵 + 𝑘𝐴 𝐶𝐴 + 𝑘𝐵 𝐶𝐵 )
𝐴

Si la reacción se lleva a cabo con superficie saturada, se tiene:

1 𝐶
𝐾1 (𝑃𝐻2 − 𝐾 𝐶𝐵 )
𝑒𝑞 𝐴
𝑟= 2
𝐶
(𝐾2 √𝐶𝐴 𝐶𝐵 + 𝐾3 √ 𝐶𝐵 + 𝑘𝐴 𝐶𝐴 + 𝑘𝐵 𝐶𝐵 )
𝐴
Para lo cual se obtienen los siguientes valores:

𝐾1 = 42,273 𝐾𝑒𝑞 = 4,490 𝐾2 = 276,721 𝐾3 = 250,678

𝑘𝐵 = 744,470 𝑘𝐴 = 209,185 𝐸𝑟𝑐𝑑 = 0,004

Asumiendo que la adsorción de estireno es la etapa controlante:

1
𝐴 + 𝑆 ↔ 𝐴𝑆 𝑟 = 𝑘𝑑 𝐶𝐴 𝐶𝑆 − 𝑘𝑑−1 𝐶𝐴𝑆 = 𝑘𝑑 (𝐶𝐴 𝐶𝑠 − 𝐶 )
𝑘𝐴 𝐴𝑆
Se despejas las concentraciones desconocidas y se remplazan en el balance de sitios 𝐶𝑇𝑆 :
2
𝐶𝐻𝑆
𝑘𝐻2 = → 𝐶𝐻𝑆 = √𝑘𝐻2 𝑃𝐻2 ∙ 𝐶𝑆 𝐸𝑐. 1
𝑃𝐻2 𝐶𝑆2
𝐶𝐴𝐻𝑆 𝐶𝑆 𝐶𝐴𝐻𝑆 𝐶𝑆 𝐶𝐴𝐻𝑆 𝐶𝑆 𝐶𝐴𝐻𝑆
𝑘𝑆1 = → 𝐶𝐴𝑆 = = = 𝐸𝑐. 2
𝐶𝐻𝑆 𝐶𝐴𝑆 𝑘𝑆1 𝐶𝐻𝑆 𝑘𝑆1 √𝑘𝐻2 𝑃𝐻2 ∙ 𝐶𝑆 𝑘𝑆1 √𝑘𝐻2 𝑃𝐻2

𝐶𝑆 𝐶𝐵𝑆 𝐶𝐵𝑆 𝐶𝑆 𝐶𝐵𝑆 𝐶𝑆 𝐶𝐵𝑆

𝑘𝑆2 = → 𝐶𝐴𝐻𝑆 = = → 𝐶𝐴𝐻𝑆 = 𝐸𝑐. 3
𝐶𝐻𝑆 𝐶𝐴𝐻𝑆 𝑘𝑆2 𝐶𝐻𝑆 𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2 ∙ 𝐶𝑆 𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2

𝐶𝐵𝑆
𝑘𝐵 = → 𝐶𝐵𝑆 = 𝑘𝐵 𝐶𝐵 𝐶𝑆 𝐸𝑐. 4
𝐶𝐵 𝐶𝑆
Sustituyendo la ecuación 4 en la ecuación 3, se tiene:
 𝑘𝐵 𝐶𝐵 𝐶𝑆
𝐶𝐴𝐻𝑆 = 𝐸𝑐. 5
𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2

Sustituyendo la ecuación 5 en la ecuación 2, se llega:

𝑘𝐵 𝐶𝐵 𝐶𝑆
𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2 𝑘𝐵 𝐶𝐵
𝐶𝐴𝑆 = = ∙𝐶 𝐸𝑐. 6
𝑘𝑆1 √𝑘𝐻2 𝑃𝐻2 𝑘𝑆2 𝑘𝑆1 𝑘𝐻2 𝑃𝐻2 𝑆

Reemplazando las ecuaciones 1, 6, 5 y 4 en el balance de sitios totales y despejando 𝐶𝑆 :

𝐶𝑇𝑆 = 𝐶𝑆 + 𝐶𝐻𝑆 + 𝐶𝐴𝑆 + 𝐶𝐴𝐻𝑆 + 𝐶𝐵𝑆

𝑘𝐵 𝐶𝐵 𝑘𝐵 𝐶𝐵
𝐶𝑇𝑆 = 𝐶𝑆 + √𝑘𝐻2 𝑃𝐻2 ∙ 𝐶𝑆 + ∙ 𝐶𝑆 + ∙ 𝐶𝑆 + 𝑘𝐵 𝐶𝐵 𝐶𝑠
𝑘𝑆2 𝑘𝑆1 𝑘𝐻2 𝑃𝐻2 𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2

𝐶𝑇𝑆
𝐶𝑆 = 𝐸𝑐. 7
𝑘𝐵 𝐶𝐵 𝑘𝐵 𝐶𝐵
1 + √𝑘𝐻2 𝑃𝐻2 + + + 𝑘𝐵 𝐶𝐵
𝑘𝑆2 𝑘𝑆1 𝑘𝐻2 𝑃𝐻2 𝑘𝑆 √𝑘𝐻 𝑃𝐻
2 2 2

Reemplazando la ecuación 6 en la expresión de velocidad de reacción:

1 𝑘𝐵 𝐶𝐵 1 𝑘𝐵 𝐶𝐵
𝑟 = 𝑘𝑑 (𝐶𝐴 𝐶𝑠 − ∙𝐶 ) 𝑟 = 𝑘𝑑 (𝐶𝐴 − )𝐶 𝐸𝑐. 8
𝑘𝐴 𝑘𝑆2 𝑘𝑆1 𝑘𝐻2 𝑃𝐻2 𝑆 𝑘𝐴 𝑘𝑆2 𝑘𝑆1 𝑘𝐻2 𝑃𝐻2 𝑆

Reemplazando la ecuación 7 en la 8 se tiene:

1 𝑘𝐵 𝐶𝐵 𝐶𝑇𝑆
𝑟 = 𝑘𝑑 (𝐶𝐴 − )∙
𝑘𝐴 𝑘𝑆2 𝑘𝑆1 𝑘𝐻2 𝑃𝐻2 1 + 𝑘 𝑃 + 𝑘 𝐶 𝑘𝐵 𝐶𝐵
√ 𝐻2 𝐻2 𝑘 𝑘 𝐵𝑘 𝐵 𝑃 + + 𝑘𝐵 𝐶𝐵
𝑆2 𝑆1 𝐻2 𝐻2 𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2

De donde:
𝑘𝐻2 𝑘𝑆2 𝑘𝑆1 𝑘𝐴 𝑘𝐵 𝑘𝐵
𝐾1 = 𝐶𝑇𝑆 𝑘𝑑 𝐾𝑒𝑞 = 𝐾2 = + + 𝑘𝐵 𝐾3 = √𝑘𝐻2 𝑃𝐻2
𝑘𝐵 𝑘𝑆2 𝑘𝑆1 𝑘𝐻2 𝑃𝐻2 𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2

Por lo tanto, se llega a la expresión final de la velocidad de reacción:

1 𝐶 1
𝐾1 (𝐶𝐴 − 𝐾 𝑃 𝐵 ) 𝐾1 (𝐶𝐴 − ∗ 𝐶𝐵 )
𝑒𝑞 𝐻2 ∗
𝐾𝑒𝑞
𝑟= 𝐾𝑒𝑞 = 𝐾𝑒𝑞 ∙ 𝑃𝐻2 𝑟=
1 + 𝐾3 + 𝐾2 𝐶𝐵 1 + 𝐾3 + 𝐾2 𝐶𝐵
Si la reacción se lleva a cabo con superficie saturada, se tiene:
1
𝐾1 (𝐶𝐴 − ∗ 𝐶𝐵 )
𝐾𝑒𝑞
𝑟=
𝐾3 + 𝐾2 𝐶𝐵
Para lo cual se obtienen los siguientes valores:
∗
𝐾1 = 41,784 𝐾𝑒𝑞 = 1,49𝑥108 𝐾2 = 14913,080 𝐾3 = 297,492 𝐸𝑟𝑐𝑑 = 0,047

Asumiendo que la inserción del primer hidrógeno es la etapa controlante:

1
𝐻𝑆 + 𝐴𝑆 ↔ 𝐴𝐻𝑆 + 𝑆 𝑟 = 𝑘𝑑 𝐶𝐴𝑆 𝐶𝐻𝑆 − 𝑘𝑑−1 𝐶𝐴𝐻𝑆 𝐶𝑆 = 𝑘𝑑 (𝐶𝐴𝑆 𝐶𝐻𝑆 − 𝐶 𝐶)
𝑘𝑆1 𝐴𝐻𝑆 𝑆

Se despejas las concentraciones desconocidas y se remplazan en el balance de sitios 𝐶𝑇𝑆 :

2
𝐶𝐻𝑆
𝑘𝐻2 = → 𝐶𝐻𝑆 = √𝑘𝐻2 𝑃𝐻2 ∙ 𝐶𝑆 𝐸𝑐. 1
𝑃𝐻2 𝐶𝑆2
𝐶𝐴𝑆
𝑘𝐴 = → 𝐶𝐴𝑆 = 𝑘𝐴 𝐶𝐴 𝐶𝑆 𝐸𝑐. 2
𝐶𝐴 𝐶𝑆
𝐶𝑆 𝐶𝐵𝑆 𝐶𝐵𝑆 𝐶𝑆 𝐶𝐵𝑆 𝐶𝑆 𝐶𝐵𝑆
𝑘𝑆2 = → 𝐶𝐴𝐻𝑆 = = → 𝐶𝐴𝐻𝑆 = 𝐸𝑐. 3
𝐶𝐻𝑆 𝐶𝐴𝐻𝑆 𝑘𝑆2 𝐶𝐻𝑆 𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2 ∙ 𝐶𝑆 𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2

𝐶𝐵𝑆
𝑘𝐵 = → 𝐶𝐵𝑆 = 𝑘𝐵 𝐶𝐵 𝐶𝑆 𝐸𝑐. 4
𝐶𝐵 𝐶𝑆
Sustituyendo la ecuación 4 en la ecuación 3, se tiene:
𝑘𝐵 𝐶𝐵 𝐶𝑆
𝐶𝐴𝐻𝑆 = 𝐸𝑐. 5
𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2
Reemplazando las ecuaciones 1, 2, 5 y 4 en el balance de sitios totales y despejando 𝐶𝑆 :

𝐶𝑇𝑆 = 𝐶𝑆 + 𝐶𝐻𝑆 + 𝐶𝐴𝑆 + 𝐶𝐴𝐻𝑆 + 𝐶𝐵𝑆

𝑘𝐵 𝐶𝐵
𝐶𝑇𝑆 = 𝐶𝑆 + √𝑘𝐻2 𝑃𝐻2 ∙ 𝐶𝑆 + 𝑘𝐴 𝐶𝐴 𝐶𝑆 + ∙ 𝐶𝑆 + 𝑘𝐵 𝐶𝐵 𝐶𝑠
𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2

𝐶𝑇𝑆
𝐶𝑆 = 𝐸𝑐. 6
𝑘𝐵 𝐶𝐵
1 + √𝑘𝐻2 𝑃𝐻2 + 𝑘𝐴 𝐶𝐴 + + 𝑘𝐵 𝐶𝐵
𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2

Reemplazando las ecuaciones 1, 2 y 5 en la expresión de velocidad de reacción:

1 𝑘𝐵 𝐶𝐵 𝐶𝑆
𝑟 = 𝑘𝑑 (√𝑘𝐻2 𝑃𝐻2 𝐶𝑆 ∙ 𝑘𝐴 𝐶𝐴 𝐶𝑆 − 𝐶𝑆 ∙ )
𝑘𝑆1 𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2

1 𝑘𝐵 𝐶𝐵
𝑟 = 𝑘𝑑 (√𝑘𝐻2 𝑃𝐻2 ∙ 𝑘𝐴 𝐶𝐴 − ) 𝐶2 𝐸𝑐. 7
𝑘𝑆1 𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2 𝑆

Reemplazando la ecuación 6 en la 7 se tiene:

2

1 𝑘𝐵 𝐶𝐵 𝐶𝑇𝑆
𝑟 = 𝑘𝑑 √𝑘𝐻2 𝑃𝐻2 𝑘𝐴 (𝐶𝐴 − )
𝑘𝑆1 𝑘𝑆2 𝑘𝐴 𝑘𝐻2 𝑃𝐻2 1 + 𝑘 𝑃 + 𝑘 𝐶 + 𝑘𝐵 𝐶𝐵
√ 𝐻2 𝐻2 𝐴 𝐴 + 𝑘𝐵 𝐶𝐵
[ 𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2 ]
De donde:
𝑘𝐻2 𝑘𝑆2 𝑘𝑆1 𝑘𝐴 𝑘𝐵
𝐾1 = 𝑘𝐴 𝐾3 𝐶𝑇2𝑆 𝑘𝑑 𝐾𝑒𝑞 = 𝐾2 = + 𝑘𝐵 𝐾3 = √𝑘𝐻2 𝑃𝐻2
𝑘𝐵 𝑘𝑆2 √𝑘𝐻2 𝑃𝐻2

Por lo tanto, se llega a la expresión final de la velocidad de reacción:

1 𝐶𝐵 1
𝐾1 (𝐶𝐴 − ) 𝐾1 (𝐶𝐴 −
∗ 𝐶𝐵 )
𝐾𝑒𝑞 𝑃𝐻2 ∗
𝐾𝑒𝑞
𝑟= 𝐾𝑒𝑞 = 𝐾𝑒𝑞 ∙ 𝑃𝐻2 𝑟=
(1 + 𝐾3 + 𝐾2 𝐶𝐵 + 𝑘𝐴 𝐶𝐴 )2 (1 + 𝐾3 + 𝐾2 𝐶𝐵 + 𝑘𝐴 𝐶𝐴 )2
Si la reacción se lleva a cabo con superficie saturada, se tiene:
1
𝐾1 (𝐶𝐴 −∗ 𝐶𝐵 )
𝐾𝑒𝑞
𝑟=
(𝐾3 + 𝐾2 𝐶𝐵 + 𝑘𝐴 𝐶𝐴 )2
Para lo cual se obtienen los siguientes valores:

𝐾1 = 4609,480 𝐾𝑒𝑞 = 1,092𝑥109 𝐾2 = 1596,130 𝐾3 = 2,291

𝑘𝐴 = 841,219 𝐸𝑟𝑐𝑑 = 0,0065

Asumiendo que la inserción del segundo hidrógeno es la etapa controlante:

1
𝐻𝑆 + 𝐴𝐻𝑆 ↔ 𝐵𝑆 + 𝑆 𝑟 = 𝑘𝑑 𝐶𝐴𝐻𝑆 𝐶𝐻𝑆 − 𝑘𝑑−1 𝐶𝐵𝑆 𝐶𝑆 = 𝑘𝑑 (𝐶𝐴𝐻𝑆 𝐶𝐻𝑆 − 𝐶 𝐶)
𝑘𝑆2 𝐵𝑆 𝑆

Se despejas las concentraciones desconocidas y se remplazan en el balance de sitios 𝐶𝑇𝑆 :

2
𝐶𝐻𝑆
𝑘𝐻2 = → 𝐶𝐻𝑆 = √𝑘𝐻2 𝑃𝐻2 ∙ 𝐶𝑆 𝐸𝑐. 1
𝑃𝐻2 𝐶𝑆2
𝐶𝐴𝑆
𝑘𝐴 = → 𝐶𝐴𝑆 = 𝑘𝐴 𝐶𝐴 𝐶𝑆 𝐸𝑐. 2
𝐶𝐴 𝐶𝑆

𝐶𝐴𝐻𝑆 𝐶𝑆 𝑘𝑆1 𝐶𝐻𝑆 𝐶𝐴𝑆 𝑘𝑆1 √𝑘𝐻2 𝑃𝐻2 𝐶𝑆 ∙ 𝑘𝐴 𝐶𝐴 𝐶𝑆

𝑘𝑆1 = → 𝐶𝐴𝐻𝑆 = = = 𝑘𝐴 𝑘𝑆1 √𝑘𝐻2 𝑃𝐻2 𝐶𝐴 𝐶𝑆 𝐸𝑐. 3
𝐶𝐻𝑆 𝐶𝐴𝑆 𝐶𝑆 𝐶𝑆
𝐶𝐵𝑆
𝑘𝐵 = → 𝐶𝐵𝑆 = 𝑘𝐵 𝐶𝐵 𝐶𝑆 𝐸𝑐. 4
𝐶𝐵 𝐶𝑆
Reemplazando las ecuaciones 1, 2, 3 y 4 en el balance de sitios totales y despejando 𝐶𝑆 :
 𝐶𝑇𝑆 = 𝐶𝑆 + 𝐶𝐻𝑆 + 𝐶𝐴𝑆 + 𝐶𝐴𝐻𝑆 + 𝐶𝐵𝑆

𝐶𝑇𝑆 = 𝐶𝑆 + √𝑘𝐻2 𝑃𝐻2 ∙ 𝐶𝑆 + 𝑘𝐴 𝐶𝐴 𝐶𝑆 + 𝑘𝐴 𝑘𝑆1 √𝑘𝐻2 𝑃𝐻2 𝐶𝐴 𝐶𝑆 + 𝑘𝐵 𝐶𝐵 𝐶𝑠

𝐶𝑇𝑆
𝐶𝑆 = 𝐸𝑐. 5
1 + √𝑘𝐻2 𝑃𝐻2 + 𝑘𝐴 𝐶𝐴 + 𝑘𝐴 𝑘𝑆1 √𝑘𝐻2 𝑃𝐻2 𝐶𝐴 + 𝑘𝐵 𝐶𝐵

Reemplazando las ecuaciones 1, 3 y 4 en la expresión de velocidad de reacción:

1
𝑟 = 𝑘𝑑 (√𝑘𝐻2 𝑃𝐻2 𝐶𝑆 ∙ 𝑘𝐴 𝑘𝑆1 √𝑘𝐻2 𝑃𝐻2 𝐶𝐴 𝐶𝑆 − 𝑘 𝐶 𝐶 ∙𝐶 )
𝑘𝑆2 𝐵 𝐵 𝑆 𝑆

1
𝑟 = 𝑘𝑑 (𝑘𝑆1 𝑘𝐻2 𝑃𝐻2 𝑘𝐴 𝐶𝐴 − 𝑘 𝐶 ) 𝐶2 𝐸𝑐. 6
𝑘𝑆2 𝐵 𝐵 𝑆

Reemplazando la ecuación 5 en la 6 se tiene:

2
1 𝑘𝐵 𝐶𝑇𝑆
𝑟 = 𝑘𝑑 𝑘𝑆1 𝑘𝐻2 𝑃𝐻2 𝑘𝐴 (𝐶𝐴 − 𝐶𝐵 ) [ ]
𝑘𝑆2 𝑘𝑆1 𝑘𝐴 𝑘𝐻2 𝑃𝐻2 1 + √𝑘𝐻2 𝑃𝐻2 + 𝑘𝐴 𝐶𝐴 + 𝑘𝐴 𝑘𝑆1 √𝑘𝐻2 𝑃𝐻2 𝐶𝐴 + 𝑘𝐵 𝐶𝐵

De donde:
𝑘𝐻2 𝑘𝑆2 𝑘𝑆1 𝑘𝐴
𝐾1 = 𝑘𝐴 𝑘𝑆1 𝐾32 𝐶𝑇2𝑆 𝑘𝑑 𝐾𝑒𝑞 = 𝐾2 = 𝑘𝐴 𝑘𝑆1 √𝑘𝐻2 𝑃𝐻2 + 𝑘𝐴 𝐾3 = √𝑘𝐻2 𝑃𝐻2
𝑘𝐵
Por lo tanto, se llega a la expresión final de la velocidad de reacción:
1 𝐶 1
𝐾1 (𝐶𝐴 − 𝐾 𝑃 𝐵 ) 𝐾1 (𝐶𝐴 −
∗ 𝐶𝐵 )
𝑒𝑞 𝐻2 ∗
𝐾𝑒𝑞
𝑟= 𝐾𝑒𝑞 = 𝐾𝑒𝑞 ∙ 𝑃𝐻2 𝑟=
(1 + 𝐾3 + 𝐾2 𝐶𝐴 + 𝑘𝐵 𝐶𝐵 )2 (1 + 𝐾3 + 𝐾2 𝐶𝐴 + 𝑘𝐵 𝐶𝐵 )2
Si la reacción se lleva a cabo con superficie saturada, se tiene:
1
𝐾1 (𝐶𝐴 − ∗ 𝐶𝐵 )
𝐾𝑒𝑞
𝑟=
(𝐾3 + 𝐾2 𝐶𝐴 + 𝑘𝐵 𝐶𝐵 )2
Para lo cual se obtienen los siguientes valores:

𝐾1 = 4609,480 𝐾𝑒𝑞 = 1,092𝑥109 𝐾2 = 1596,130 𝐾3 = 2,291

𝑘𝐴 = 841,219 𝐸𝑟𝑐𝑑 = 0,0065