1.

Determine el perfil de reacción del siguiente sistema en un reactor por lotes en

fase líquida y densidad constante después de 4 horas, cuando la concentración
inicial de A y B es de 10 mol/m3 y las constantes de velocidad k1 = 0.1 mol/m3·h,
k2 = 0.1 mol/m3·h.
𝐴+𝐵⟶𝐶
𝐵+𝐶⟶𝐷
𝑑𝐶𝐴
= −𝑘1 ∗ [𝐴] ∗ [𝐵]
𝑑𝑡
𝑑𝐶𝐵
= −𝑘1 ∗ [𝐴] ∗ [𝐵] − 𝑘2 ∗ [𝐵] ∗ [𝐶]
𝑑𝑡
𝑑𝐶𝐶
= 𝑘1 ∗ [𝐴] ∗ [𝐵] − 𝑘2 ∗ [𝐵] ∗ [𝐶]
𝑑𝑡
𝑑𝐶𝐷
= 𝑘2 ∗ [𝐵] ∗ [𝐶]
𝑑𝑡

k1 = 0.1 mol/m3·h
k2 = 0.1 mol/m3·h
t_total = 4 horas
𝐶𝐴0 = 10 mol/m3
𝐶𝐵0 = 10 mol/m3
𝐶𝐶0 = 0 mol/m3 (concentración inicial de C)
𝐶𝐷0 = 0 mol/m3 (concentración inicial de D)
t = 0 tiempo inicial
𝐶𝐴 = 𝐶𝐴0 concentración inicial de A
𝐶𝐵 = 𝐶𝐵0 concentración inicial de B
𝐶𝐶 = 𝐶𝐶0 concentración inicial de C
𝐶𝐷 = 𝐶𝐷0 concentración inicial de D
r1 = k1 * 𝐶𝐴 * 𝐶𝐵
r2 = k2 * 𝐶𝐵 * 𝐶𝐶
𝑑𝐶𝐴 = -r1
𝑑𝐶𝐵 = -r1 - r2
 𝑑𝐶𝐶 = r1 - r2
𝑑𝐶𝑐 = r2
𝐶𝐴 += 𝑑𝐶𝐴 * tincremento
𝐶𝐵 += 𝑑𝐶𝐵 * tincremento
𝐶𝐶 += 𝑑𝐶𝑐 * tincremento
𝐶𝐷 += 𝑑𝐶𝐷 * tincremento

𝐶𝐴 = 3.2939
𝐶𝐶 =3.7583
𝐶𝐵 = 0.3464
𝐶𝐷 = 2.9476

2. La reacción en fase líquida (A y B equimolar) mostrada se lleva isotérmicamente

en un reactor PFR y en un reactor CSTR a 27°C, con flujo de 2 dm3 /s y CA0=0.1
Kmol/m3 .
A+B=C
𝑘 = 0.01 𝑑𝑚3 /𝑚𝑜𝑙 ∙ 𝑠 𝑎 300 𝐾
𝐸 = 10,000 𝑐𝑎𝑙 /𝑚𝑜𝑙
𝐻𝐵 ° 273 = −15 𝑘𝑐𝑎𝑙 /𝑚𝑜𝑙
𝐻𝐴 , ° 273 = −20 𝑘𝑐𝑎𝑙 /𝑚𝑜𝑙 ,
𝐻𝐶 ° 273 = −41 𝑘𝑐𝑎𝑙 /𝑚𝑜𝑙
𝐶𝑃𝐴 = 𝐶𝑃𝐵 = 15 𝑐𝑎𝑙/ 𝑚𝑜𝑙 ∙ 𝐾 𝐶𝑃𝐶 = 30 𝑐𝑎𝑙 /𝑚𝑜𝑙 ∙ k
(a) Calcule el volúmen de un PFR que debe tenerse para lograr una conversión del
85%.
(b) Grafica la en función del volumen del reactor.
a formulacion de diseño del PFR se puede expresar como:

V = V0 * X / (-rA)
Donde: V es el volumen del reactor PFR V0 es el volumen de entrada del flujo X es
la conversión de A rA es la velocidad de reacción de A

La velocidad de reacción de A se puede determinar utilizando la ecuación de

Arrhenius y la ley de acción de masas:

rA = k * CA * CB

Donde: k es la constante de velocidad CA y CB son las concentraciones molares de

A y B, respectivamente

500 dm -> X=0.977

250dm-> X_1=0.967 and X_2=0.992

CSTR:
X=0.977.

3.Una reacción elemental en serie en fase líquida se efectúa en un reactor

intermitente de 500 dm3. La concentración inicial de A es 1.6 mol/dm3, el producto
deseado es B. Si cada reacción sea irreversible, grafique las concentraciones en
función del tiempo. Compara con un CSTR.
𝐴⟶𝐵⟶𝐶
k1 = 0.4 ℎ −1
𝑘2 = 0.01 ℎ −1 100°c
r1 = -k1 * [A]
r2 = -k2 * [B]
Para el reactor intermitente de 500 dm3, la concentración inicial de A es 1,6
mol/dm3. Para el CSTR, también asumiremos una concentración inicial de A de 1,6
mol/dm3 y un volumen de 500 dm3.
Las leyes de velocidad para estas reacciones son las siguientes
𝑟 1 = 𝐾 1 𝐶𝐴
𝑟 2 = 𝐾 2 𝐶𝐵
Para el reactor intermitente, la ecuación de diseño se puede expresar como:

𝑉 = 𝑉0 + 𝑛 ∗ (𝑟 1 ∗ 𝑡1 + 𝑟 2 ∗ 𝑡2),
Para CSTR
𝑉𝐶𝑆𝑇𝑅 = 𝑛 ∗ (𝐹1 ∗ 𝑡1 + 𝐹2 ∗ 𝑡2),
Ahora, podemos proceder a graficar las concentraciones en función del tiempo para
comparar el reactor intermitente y el reactor CSTR.
import numpy as np
import matplotlib.pyplot as plt
CA0 = 1.6
k1 = 0.4
k2 = 0.01
t1 = np.linspace(0, 10, 100)
t2 = np.linspace(0, 10, 100)
CB_intermittent = (k1 * CA0 * (1 - np.exp(-k1 * t1))) / (k1 + k2)
CC_intermittent = (k1 * k2 * CA0 * (1 - np.exp(-k1 * t1)) * (1 - np.exp(-k2 * t2))) / ((k1
+ k2) * k2)
F1 = CA0 * 500
F2 = CB_intermittent[-1] * 500
CB_CSTR = (F1 * t1) / (F1 + k1)
CC_CSTR = (F1 * k2 * t1 + F2 * t2) / (F1 + k1)
plt.plot(t1, CB_intermittent, label='CB (Intermitente)')
plt.plot(t1, CC_intermittent, label='CC (Intermitente)')
plt.plot(t1, CB_CSTR, label='CB (CSTR)')
plt.plot(t1, CC_CSTR, label='CC (CSTR)')
plt.xlabel('Tiempo (h)')
plt.ylabel('Concentración (mol/dm3)')
plt.title('Concentraciones en función del tiempo')
plt.legend()
plt.grid(True)
plt.show()
De este código usado en Python obtenemos la siguiente gráfica
4. Traza el perfil isotérmico para la reacción mostrada que se realiza en fase gas
con condiciones iniciales de T0 = 300 K, a una presión de 5 bar y V0 = 0,5 m³. El
calor de reacción es –6280 KJ/kmol y la capacidad calorífica de A, R y S son 125.6,
104.7 y 83.7 kJ/kmol K, respectivamente.
𝑘 = 1014𝑒𝑥𝑝 (−10 000/T)h^-1
𝐴⟶𝑅+𝑆
El balance de energia para un reactor en fase gas sin transferencia de calor externo
se puede expresar como:

𝑉0 * 𝐶𝑝𝐴 * (𝑇 − 𝑇0) + Δ𝐻𝑟 * 𝑉 * 𝑟 = 0

Donde: 𝑉0 es el volumen inicial del reactor 𝐶𝑝𝐴 es la capacidad calorífica de A 𝑇0

es la temperatura inicial 𝑇 es la temperatura en el reactor Δ𝐻𝑟 es el calor de reacción
𝑉 es el volumen en el reactor en un momento dado 𝑟 es la velocidad de reacción

La velocidad de reacción se puede expresar utilizando la ecuación de Arrhenius:

𝑟 = 𝑘 * 𝑁𝐴 * 𝐶𝐴 * exp(-Δ𝐻𝑟 / (𝑅 * 𝑇))

Donde: 𝑘 es la constante de velocidad 𝑁𝐴 es el número de Avogadro 𝐶𝐴 es la

concentración de A 𝑅 es la constante de los gases ideales
# Parámetros
T0 = 300 # K
P0 = 5 # bar
V0 = 0.5 # m^3
delta_Hr = -6280 # kJ/kmol
Cp_A = 125.6 # kJ/kmol K
Cp_R = 104.7 # kJ/kmol K
Cp_S = 83.7 # kJ/kmol K
k = lambda T: 1014 * np.exp(-10000 / T) # h^-1
R = 8.314 # J/mol K

# Conversión de unidades
P0 *= 1e5 # bar to Pa
V0 *= 1e3 # m^3 to L

# Definición de la ecuación diferencial

def reactor_profile(T, V):
Ca = P0 / (R * T) # Concentración de A
r = k(T) * Ca # Velocidad de reacción
dT = -delta_Hr * r / (V0 * Cp_A) # Cambio de temperatura
return dT

# Resolución del perfil isotérmico utilizando el método de Euler

V = np.linspace(0, V0, 100)
T = np.zeros_like(V)
T[0] = T0

for i in range(1, len(V)):

dV = V[i] - V[i-1]
dT = reactor_profile(T[i-1], V[i-1])
T[i] = T[i-1] + dT * dV

# Gráfico del perfil isotérmico

plt.plot(V, T)
plt.xlabel('Volumen (L)')
plt.ylabel('Temperatura (K)')
plt.title('Perfil isotérmico')
plt.grid(True)
plt.show()
5.Cual será el mejor reactor de flujo para obtener la selectividad de B en las
reacciones mostradas, para n=0, 1 y 2.

Para determinar el mejor reactor de flujo para obtener la selectividad de B en las

reacciones mostradas, debemos comparar los resultados de selectividad para
diferentes valores de "n".
Primera reacción
,𝑟 1 = 𝐾 1 𝐶𝐴 , la selectividad de B se puede expresar como
𝑟1
𝑆𝐵 = , donde r2 es la velocidad de reacción de formación de productos no
(𝑟 1 +𝑟 2 )
deseados. En este caso, la selectividad de B no depende de "n".
Segunda reacción
𝑟1
𝑟 1 = 𝐾2 𝐶𝐴𝑛 , la selectividad de B se puede expresar como 𝑆𝐵 = ,. Aquí es
(𝑟 1 +𝑟 2 )
donde el valor de "n" tiene un impacto en la selectividad de B.
Para n = 0, la segunda reacción se convierte en una reacción de orden cero con
1
respecto a A, y la selectividad de B se puede expresar como 𝑆𝐵 = 𝑟 2 , SB = 1 /
(1+ )
𝑟1

(1 + r2/r1). En este caso, la selectividad de B no depende de la concentración inicial

de A.
Para n = 1
la segunda reacción se convierte en una reacción de primer orden con respecto a
1
A, y la selectividad de B se puede expresar como 𝑆𝐵 = 𝑟 2 , SB = 1 / (1 + (r2 /
(1+ )
𝐾 2 𝐶𝐴

(K2 * CA))). En este caso, la selectividad de B dependerá de la concentración inicial

de A.
Para n = 2
La segunda reacción se convierte en una reacción de segundo orden con respecto
1
a A, y la selectividad de B se puede expresar como como 𝑆𝐵 = 𝑟 2 ,.Aquí, la
(1+ )
𝐾2 𝐶2

selectividad de B también dependerá de la concentración inicial de A.

Por lo tanto, para determinar el mejor reactor de flujo en términos de selectividad de
B, debemos analizar las condiciones específicas del problema, como la
concentración inicial de A y los valores de las constantes de velocidad (K1 y K2), y
comparar los resultados de selectividad para los diferentes valores de n. No es
posible determinar de manera general cuál será el mejor reactor de flujo sin tener
en cuenta estos valores.

6.Cual será el mejor reactor de flujo para obtener la selectividad de B en las

reacciones mostradas, siendo k1= 2 L/mol min, k2= 1 min-1, flujo volumétrico de 20
L/min, V de 100 L y una concentración inicial del reactante de 4 mol/L.
A=B r1=K1CA^2
A=C r2=K2CA
Dado que las reacciones son de primer orden en ambos casos, podemos calcular
las concentraciones de A, B y C en función del volumen utilizando las ecuaciones
cinéticas y las condiciones iniciales.
Para el reactor de flujo en serie (PFR, por sus siglas en inglés), podemos utilizar las
ecuaciones de diseño y resolverlas numéricamente para obtener las
concentraciones en el volumen final. Sin embargo, para el reactor de flujo continuo
de tanque agitado (CSTR, por sus siglas en inglés), podemos utilizar el equilibrio de
estado estacionario para obtener directamente la concentración de B.
Para determinar el mejor reactor de flujo para obtener la selectividad de B en las
reacciones dadas, podemos comparar el rendimiento de cada reactor en términos
de la conversión y la selectividad.
PFR
𝑑𝐶𝐴
= −𝑟1
𝑑𝑣
Integrando esta ecuación diferencial, obtenemos la relación de conversión en
función del volumen del reactor:
X = (𝑪𝑨𝟎 -𝑪𝑨 ) /𝑪𝑨𝟎 = K1V / ((𝒗𝟎 * (𝑪𝑨𝟎 )
CSTR
V = 𝒗𝟎 * / (-r1) * (𝑪𝑨𝟎 - 𝑪𝑨 )
Para calcular la selectividad de B, podemos utilizar la siguiente fórmula:
𝑪𝑩
𝑆𝐵 = 𝑪 ,.
(𝑪𝑩 + 𝑩 )
𝑪𝑪

Comparando ambos reactores, podemos calcular la selectividad de B para

diferentes valores de conversión y determinar cuál reactor proporciona la mayor
selectividad.
7. Para estimar el rendimiento de B en un CSTR y un PFR al reaccionar el 90% de
A cuando CAo = 2 mol/L, k2 = 0.2 L/mol min y k1 = 0.1 min^-1, debemos resolver
las ecuaciones cinéticas y calcular las concentraciones de A, B y C en ambos
reactores.
La fórmula cinética para la isomerización de A a B es r1 = k1 * CA, y la fórmula
cinética para la dimerización de A a C es r2 = k2 * CA^2.
Dado que el 90% de A debe reaccionar, podemos establecer la concentración de A
al final de la reacción como 0.1 * CAo. Utilizaremos esta concentración para calcular
las concentraciones finales de B y C.
Reactor CSTR:
En un CSTR, la concentración en el reactor se mantiene constante. Por lo tanto,
cuando se alcanza el equilibrio, la velocidad de producción de B debe ser igual a la
velocidad de consumo de A. Podemos establecer la siguiente ecuación:
𝑟1 = 𝐾 1 𝐶𝐴 − 𝐾 2 𝐶 𝐴2

Al reaccionar el 90% de A, tenemos:

2
0.1 𝐶𝐴0 − 0.2𝐶 𝐴0 = 0.9 𝐶𝐴0

Resolviendo esta ecuación, encontraremos el valor de CAo que cumple con la

condición de reacción del 90% de A. Con este valor, podemos calcular la
concentración de B en el equilibrio.
Reactor PFR:
𝑑𝐶𝐴
= −𝑟1
𝑑𝑣

Podemos integrar esta ecuación desde 0 hasta V para obtener la conversión en

función del volumen. Luego, podemos calcular la concentración de B en el volumen
V correspondiente a una conversión del 90% de A.

roble
8.Escriba las ecuaciones de balance de masa de reactor continuos que deben
resolverse para el
sistema de reacción mostrado
CH2+2 O2= CO2+H2O r1=k1CCH4CO
CO+H2O=CO+H2 r2=kf2CCOCH2O-Kb2CCO2CH2
CH4+H2O=CO+3H2 r3=k3CCH4CH20
CH4+CO2=2CO+2H2 r4=k4CCH4CCO2
Para el compuesto CH4 (metano):
Entrada: FCH4_in - FCH4_out + r1 - r3 - r4 = 0
Acumulación: d(V C_CH4)/dt = F C_CH4_in - F C_CH4_out + V r1 - V r3 - V r4
Para el compuesto CO (monóxido de carbono):
Entrada: FCO_in - FCO_out + r2 + r4 = 0
Acumulación: d(V C_CO)/dt = F C_CO_in - F C_CO_out + V r2 + V r4 - V r3
Para el compuesto O2 (oxígeno):
Entrada: FO2_in - FO2_out + 2r1 = 0
Para el compuesto CO2 (dióxido de carbono):
Entrada: FCO2_in - FCO2_out + r1 + r2 - 2r4 = 0
Acumulación: d(V C_CO2)/dt = F C_CO2_in - F C_CO2_out - V r1 + V r2 + V r4
Para el compuesto H2O (agua):
Entrada: FH2O_in - FH2O_out - r1 + r2 + r3 = 0
Dónde:
F es el flujo molar de cada compuesto.
Los subíndices "in" y "out" representan la entrada y salida del reactor,
respectivamente.
r1, r2, r3 y r4 son las velocidades de reacción correspondientes a cada reacción.
k1, kf2, Kb2, k3 y k4 son las constantes de velocidad de reacción para las reacciones
r1, r2, r3 y r4, respectivamente.
C es la concentración molar de cada compuesto.