UNIVERSIDAD CENTROAMERICANA

JOSÉ SIMEÓN CAÑAS

TAREA 3

REACTORES QUÍMICOS

INTEGRANTES:
JOSÉ ANDRÉS DÍAZ RAMÍREZ 00018117
JOSÉ LUIS MAYORGA DONIS 00353917
ANDREA ISABEL PÉREZ PORTILLO 00027817
ALEJANDRO ARTURO VEGA RECINOS 00029517

CATEDRÁTICA:
MARÍA DOLORES ROVIRA

FECHA DE ENTREGA:

21 DE MAYO DEL 2022

EJERCICIO 1

La antraquinona (AQ) se produce por oxidación del antraceno (AN), resultando como producto
intermedio la antranona (ANT). Las reacciones en fase liquida son de primer orden en serie y pueden
representarse de la siguiente forma:

Las constantes cinéticas dependen de la temperatura de la siguiente forma donde las k están en
segundos (𝑠 −1 ) y la temperatura T en K:

Si la reacción ocurre en un CSTR de volumen 1 m3 a 95°C al que se alimenta un caudal de 1.5m3 /h

de antraceno a una concentración de 0.025 mol.dm-3. Determine las concentraciones de
antraquinona, antraceno y antranona a la salida del reactor.

 Paso 1

Primero enumeramos las reacciones, para este caso solo tenemos una reacción, por lo que
procedemos al paso 2.

 Paso 2

Conocemos que el reactor es un CSTR, por lo que podemos plantear las ecuaciones del balance de
materia, donde obtenemos:
𝐹𝐴𝑜 − 𝐹𝐴 + 𝑟𝐴 𝑉 = 0

−𝐹𝐵 + 𝑟𝐵 𝑉 = 0

−𝐹𝐶 + 𝑟𝐶 𝑉 = 0

 Paso 3

Ya que tenemos el balance de materia, procedemos a encontrar las ecuaciones de velocidad neta
para cada una:
𝑟𝐴 = −𝑘1 𝐶𝐴

𝑟𝐵 = 𝑘1 𝐶𝐴 − 𝑘2 𝐶𝐵

𝑟𝐶 = −𝑘2 𝐶𝐵

 Paso 4

Se sustituye los valores de la velocidad en las ecuaciones del paso 2.

𝑣𝑜 𝐶𝐴𝑜 − 𝑣𝑜 𝐶𝐴 − 𝑘1 𝐶𝐴 𝑉 = 0
 −𝑣𝑜 𝐶𝐵 + (𝑘1 𝐶𝐴 − 𝑘2 𝐶𝐵 )𝑉 = 0

−𝑣𝑜 𝐶𝐶 − 𝑘2 𝐶𝐵 𝑉 = 0

Planteamos nuestras incógnitas: CA, CB, CC, K1 y K2. Los valores de K los podemos conocer, debido a
que tenemos ecuaciones dependientes de la temperatura, por lo que obtenemos:
−6345
𝑘1 = 2.24 ∗ 104 ∗ exp ( ) = 0.000792
95
−9535
𝑘2 = 3.98 ∗ 107 ∗ exp ( ) = 0.000369
95

Al sustituir los valores de K, podemos resolver el sistema de ecuaciones:

𝑚𝑜𝑙
𝐶𝐴 = 8.62283
𝑚3
𝑚𝑜𝑙
𝐶𝐵 = 8.68865
𝑚3
𝑚𝑜𝑙
𝐶𝐶 = 7.68852
𝑚3

EJERCICIO 2

El ácido tereftálico (ATP) es muy empleado en la manufactura de fibras sintéticas (por ejemplo
dacrón) y como intermedio en películas de poliéster (como el mylar). Se estudió la formación de
tereftalato de potasio (S) a partir de benzoato de potasio (A) en un reactor.

Se encontró que los intermediarios (R) que se forman por la disociación de benzoato de potasio
sobre catalizador de CdCl2, reaccionan con tereftalado de potasio en un paso de reacción
autocatalítico.

Se carga A puro al reactor a presión de 110 kPa. Las velocidades específicas de reacción a 410°C son:
𝐸1 = 42.6 𝑘𝑐𝑎𝑙 ⁄𝑚𝑜𝑙 = 42600 𝑐𝑎𝑙 ⁄𝑚𝑜𝑙

𝐸2 = 48.6 𝑘𝑐𝑎𝑙 ⁄𝑚𝑜𝑙 = 48600 𝑐𝑎𝑙 ⁄𝑚𝑜𝑙

𝐸3 = 32 𝑘𝑐𝑎𝑙 ⁄𝑚𝑜𝑙 = 32000 𝑐𝑎𝑙 ⁄𝑚𝑜𝑙

𝑘1 = 1.08 𝑥 10−3 𝑠 −1 𝑬𝒄. 𝟏

𝑘2 = 1.19 𝑥 10−3 𝑠 −1 𝑬𝒄. 𝟐

𝑘3 = 1.59 𝑥 10−3 𝑑𝑚3⁄𝑚𝑜𝑙 𝑬𝒄. 𝟑

Literal A. Grafique la concentración de A, R y S en función del tiempo en un reactor intermitente
a 410°C observando cuando se produce el máximo de R

Siguiendo los pasos del algoritmo para resolución de reacciones múltiples:

1. Enumerando las reacciones.

𝑘1
𝐴 → 𝑅 (1)
𝑘2
𝑅 → 𝑆 (2)
𝑘3
𝑅 + 𝑆 → 2𝑆 (3)

2. Escribiendo un balance de moles para cada especie presente, utilizando la ecuación de diseño
de un reactor Batch:
𝑑(𝑁𝐴 /𝑉)
= 𝑟𝐴
𝑑𝑡
𝑑(𝐶𝐴 )
= 𝑟𝐴
𝑑𝑡
Sustituyendo para cada especie:
𝑑𝐶𝐴
𝐴→ = 𝑟𝐴
𝑑𝑡
𝑑𝐶𝑅
𝑅→ = 𝑟𝑅
𝑑𝑡
𝑑𝐶𝑆
𝑆→ = 𝑟𝑆
𝑑𝑡
3. Encontrando las velocidades de reacción netas.
𝑟𝐴 = 𝑟𝐴1 = −𝑘1 𝐶𝐴

𝑟𝑅 = 𝑟𝑅1 − 𝑟𝐵2 − 𝑟𝐵3 = 𝑘1 𝐶𝐴 − 𝑘2 𝐶𝑅 − 𝑘3 𝐶𝑅 𝐶𝑆

𝑟𝑆 = 𝑟𝑆2 − 𝑟𝑆3 = 𝑘2 𝐶𝑅 − 𝑘3 𝐶𝑅 𝐶𝑆

Sustituyendo las velocidades de reacción en los balances de masa:

𝑑𝐶𝐴
= −𝑘1 𝐶𝐴 𝑬𝒄. 𝟒
𝑑𝑡
𝑑𝐶𝑅
= 𝑘1 𝐶𝐴 − 𝑘2 𝐶𝑅 − 𝑘3 𝐶𝑅 𝐶𝑆 𝑬𝒄. 𝟓
𝑑𝑡
𝑑𝐶𝑆
= 𝑘2 𝐶𝑅 − 𝑘3 𝐶𝑅 𝐶𝑆 𝑬𝒄. 𝟔
𝑑𝑡
Encontrando el tiempo óptimo para maximizar el producto de R:
1 𝑘1
𝑡𝑜𝑝𝑡 = 𝑙𝑛 ( )
𝑘1 − 𝑘2 𝑘2
 1 1.08𝑥10−3 𝑠 −1 )
𝑡𝑜𝑝𝑡 = 𝑙𝑛 ( )
1.08𝑥10−3 𝑠 −1 −3
− 1.19𝑥10 𝑠 −1 1.19𝑥10−3 𝑠 −1 )

𝑡𝑜𝑝𝑡 = 881.75 𝑠

Además, se debe encontrar la concentración inicial de A para establecer los intervalos al momento de resolver
las ecuaciones diferenciales:
𝑃 110𝑘𝑃𝑎
𝐶𝐴𝑂 = = = 0.019 𝑚𝑜𝑙 ⁄𝑑𝑚3
𝑅𝑇 (8.314 𝑘𝑃𝑎 𝑑𝑚3 ⁄𝑚𝑜𝑙 𝐾 )(683 𝐾)

En el programa Polymath se introdujeron las ecuaciones 1, 2, 3, 4, 5 y 6 planteadas anteriormente

y además se definieron las concentraciones iniciales y un intervalo de tiempo para poder resolver el
sistema de ecuaciones diferenciales, en el caso del tiempo en el punto final se eligió arbitrariamente
pero siendo mayor que el tiempo optimo:

Los valores resultantes son los siguientes:

El gráfico de las concentraciones en función del tiempo resultante es el siguiente:

Literal B. Repita el inciso A, variando la temperatura a 430°C y 390°C

 430°C (703 K)

Utilizando la siguiente ecuación para encontrar una nueva k para cada reacción, debido a que la
constante varía con la temperatura:

𝐸 1 1
𝑘(𝑇) = 𝑘 ′ (𝑇0 )exp ( ( − ))
𝑅 𝑇0 𝑇

Encontrando k1

42600 𝑐𝑎𝑙/𝑚𝑜𝑙 1 1
𝑘1 (703 𝐾) = (1.08𝑥10−3 𝑠 −1 )exp ( ( − ))
1.987 𝑐𝑎𝑙 ⁄𝑚𝑜𝑙𝐾 683 𝐾 703 𝐾

𝑘1 (703 𝐾) = 2.64𝑥10−3 𝑠 −1

Encontrando k 2
 48600 𝑐𝑎𝑙/𝑚𝑜𝑙 1 1
𝑘2 (703 𝐾) = (1.19𝑥10−3 𝑠 −1 )exp ( ( − ))
⁄
1.987 𝑐𝑎𝑙 𝑚𝑜𝑙𝐾 683 𝐾 703 𝐾

𝑘2 (703 𝐾) = 3.3𝑥10−3 𝑠 −1

Encontrando k 3

32000 𝑐𝑎𝑙/𝑚𝑜𝑙 1 1
𝑘3 (703 𝐾) = (1.59𝑥10−3 𝑑𝑚3⁄𝑚𝑜𝑙 )exp ( ( − ))
1.987 𝑐𝑎𝑙 ⁄𝑚𝑜𝑙𝐾 683 𝐾 703 𝐾

𝑘3 (703 𝐾) = 3.1𝑥10−3 𝑑𝑚3⁄𝑚𝑜𝑙

Encontrando el tiempo óptimo para maximizar el producto de R:

1 𝑘1
𝑡𝑜𝑝𝑡 = 𝑙𝑛 ( )
𝑘1 − 𝑘2 𝑘2

1 2.64𝑥10−3 𝑠 −1
𝑡𝑜𝑝𝑡 = 𝑙𝑛 ( )
2.64𝑥10−3 𝑠 −1 −3
− 3.3𝑥10 𝑠 −1 3.3𝑥10−3 𝑠 −1

𝑡𝑜𝑝𝑡 = 338 𝑠

Encontrando la concentración inicial de A:

𝑃 110𝑘𝑃𝑎
𝐶𝐴𝑂 = = = 0.0188 𝑚𝑜𝑙 ⁄𝑑𝑚3
𝑅𝑇 (8.314 𝑘𝑃𝑎 𝑑𝑚3 ⁄𝑚𝑜𝑙 𝐾 )(703 𝐾)

Utilizando Polymath para solucionar el sistema de ecuaciones:

Los resultados de las concentraciones son los siguientes:

El gráfico resultante es el siguiente:

 390°C (663 K)

Utilizando la siguiente ecuación para encontrar una nueva k para cada reacción, debido a que la
constante varía con la temperatura:

𝐸 1 1
𝑘(𝑇) = 𝑘 ′ (𝑇0 )exp ( ( − ))
𝑅 𝑇0 𝑇

Encontrando k1
 42600 𝑐𝑎𝑙/𝑚𝑜𝑙 1 1
𝑘1 (663 𝐾) = (1.08𝑥10−3 𝑠 −1 )exp ( ( − ))
⁄
1.987 𝑐𝑎𝑙 𝑚𝑜𝑙𝐾 683 𝐾 663 𝐾

𝑘1 (663 𝐾) = 4.19𝑥10−4 𝑠 −1

Encontrando k 2

48600 𝑐𝑎𝑙/𝑚𝑜𝑙 1 1
𝑘2 (663 𝐾) = (1.19𝑥10−3 𝑠 −1 )exp ( ( − ))
1.987 𝑐𝑎𝑙 ⁄𝑚𝑜𝑙𝐾 683 𝐾 663 𝐾

𝑘2 (663 𝐾) = 4.04𝑥10−4 𝑠 −1

Encontrando k 3

32000 𝑐𝑎𝑙/𝑚𝑜𝑙 1 1
𝑘3 (663 𝐾) = (1.59𝑥10−3 𝑑𝑚3⁄𝑚𝑜𝑙 )exp ( ( − ))
1.987 𝑐𝑎𝑙 ⁄𝑚𝑜𝑙𝐾 683 𝐾 663 𝐾

𝑘3 (663 𝐾) = 7.81𝑥10−4 𝑑𝑚3⁄𝑚𝑜𝑙

Encontrando el tiempo óptimo para maximizar el producto de R:

1 𝑘1
𝑡𝑜𝑝𝑡 = 𝑙𝑛 ( )
𝑘1 − 𝑘2 𝑘2

1 4.19𝑥10−4 𝑠 −1
𝑡𝑜𝑝𝑡 = 𝑙𝑛 ( )
4.19𝑥10−4 𝑠 −1 −4
− 4.04𝑥10 𝑠 −1 4.04𝑥10−4 𝑠 −1

𝑡𝑜𝑝𝑡 = 2430 𝑠

Además se debe encontrar la concentración inicial de A para establecer los intervalos al momento de
resolver las ecuaciones diferenciales:
𝑃 110𝑘𝑃𝑎
𝐶𝐴𝑂 = = = 0.0199 𝑚𝑜𝑙 ⁄𝑑𝑚3
𝑅𝑇 (8.314 𝑘𝑃𝑎 𝑑𝑚3 ⁄𝑚𝑜𝑙 𝐾)(663𝐾)

Resolviendo el sistema de ecuaciones en Polymath:

El gráfico resultante es el siguiente:

Literal C. ¿Cuáles serían las concentraciones de salida de un CSTR que operara a 410°C y un tiempo
de residencia de 1200 s?

Para encontrar el sistema de ecuaciones en este caso se realizó una relación entre el tiempo de
residencia de un reactor intermitente y un reactor CSTR. Tomando en cuenta la ecuación de diseño
de un reactor batch:
𝑑𝑁𝐴
= 𝑟𝐴 𝑉
𝑑𝑡
Tomando en cuenta que no hay cambio de moles, se tiene que 𝑣 = 𝑣0 , por lo tanto, al dividir entre
𝑣0 la ecuación anterior y sabiendo que 𝜏 = 𝑉⁄𝑣0 se tiene lo siguiente:
𝑑𝑁𝐴 𝑟𝐴 𝑉
=
𝑣0 𝑑𝑡 𝑣0
𝑑𝑁𝐴
𝜏= 𝐸𝑐. 1
𝑣0 𝑟𝐴 𝑑𝑡

Por otro lado, la ecuación de diseño de un CSTR dividida entre 𝑣0 da el siguiente resultado:
𝐹𝐴𝑂 − 𝐹𝐴
𝑉=
−𝑟𝐴
𝑉 𝐹𝐴𝑂 − 𝐹𝐴
=
𝑣0 −𝑟𝐴 𝑣0
𝐹𝐴𝑂 − 𝐹𝐴
𝜏= 𝐸𝑐. 2
−𝑟𝐴 𝑣0

Igualando la ecuación 1 y 2:
𝑑𝑁𝐴 𝐹𝐴𝑂 − 𝐹𝐴
=
𝑣0 𝑟𝐴 𝑑𝑡 −𝑟𝐴 𝑣0
𝑑𝑁𝐴
− = 𝐹𝐴𝑂 − 𝐹𝐴 𝐸𝑐. 3
𝑑𝑡
La ecuación de diseño para un CSTR es la siguiente, por lo tanto se puede despejar 𝐹𝐴𝑂 − 𝐹𝐴 :

𝐹𝐴𝑂 − 𝐹𝐴
𝑉=
−𝑟𝐴

−𝑟𝐴 𝑉 = 𝐹𝐴𝑂 − 𝐹𝐴

Sustituyendo en la ecuación 3:
𝑑𝑁𝐴
− = −𝑟𝐴 𝑉
𝑑𝑡
𝑑(𝑁𝐴 /𝑉)
= 𝑟𝐴
𝑑𝑡
𝑑(𝐶𝐴 )
= 𝑟𝐴
𝑑𝑡
Como resultado se tiene la misma ecuación que se utilizó en el literal A para plantear las ecuaciones
diferenciales de cada especie, por lo tanto dichas ecuaciones serán las mismas en este caso y
además, siendo la temperatura 410 ºC, se utilizan las mismas constantes dadas por el problema y la
misma concentración inicial de A 𝐶𝐴𝑂 0.019 𝑚𝑜𝑙 ⁄𝑑𝑚3 .
𝑑𝐶𝐴
= −𝑘1 𝐶𝐴
𝑑𝑡
 𝑑𝐶𝑅
= 𝑘1 𝐶𝐴 − 𝑘2 𝐶𝑅 − 𝑘3 𝐶𝑅 𝐶𝑆
𝑑𝑡
𝑑𝐶𝑆
= 𝑘2 𝐶𝑅 − 𝑘3 𝐶𝑅 𝐶𝑆
𝑑𝑡
𝑘1 = 1.08 𝑥 10−3 𝑠 −1

𝑘2 = 1.19 𝑥 10−3 𝑠 −1

𝑘3 = 1.59 𝑥 10−3 𝑑𝑚3⁄𝑚𝑜𝑙

Utilizando Polymath y tomando en cuenta que t=1200 s, se tiene lo siguiente:

Las concentraciones finales para cada componente son las siguientes:

𝐶𝐴 = 0.0052 𝑚𝑜𝑙⁄𝑑𝑚3

𝐶𝐵 = 0.0063 𝑚𝑜𝑙⁄𝑑𝑚3

𝐶𝐶 = 0.0074 𝑚𝑜𝑙⁄𝑑𝑚3

Además, la tendencia del gráfico es la siguiente: