Diseño de Reactores Químicos: Flujos y Reacciones
Diseño de Reactores Químicos: Flujos y Reacciones
DISEÑO DE REACTORES II
REPASO- ALGUNOS PUNTOS NECESARIOS
BALANCES MOLARES
• Sistemas de Flujo: Pv = zFTRT
Povo zo FToRTo
Pv = zFTRT
Povo = zo FToRTo v = vo (Po/p) (T/To) (FT/FTo)
P cte !!!
v = vo (1 + e xA)
T cte !!!
• Ejercicios:
üEjemplo3-5
üEjemplo 3-6
EJEMPLO: Datos para el ejercicio:
MICRORREACTORES
Medidas típicas:
DISEÑO DE REACTORES II
REACCIONES MULTIPLES
Semestre I/2024
Serie (o consecutivas) A à B àC
Complejas (combinación A + B à C +D
de Serie y Paralelo) A+CàE
Independientes (ni
serie ni paralelo) AàB+C
DàE+F
DEFINICIONES
AàD Se debe definir claramente cual
AàU es el producto deseado!!!!
A à D àU
REACCIONES EN PARALELO:
Fogler 4 Ed.:
Problema 6-9
inciso (a)
inciso (b)
inciso (c)
ALGUNAS RELACIONES UTILES
En un reactor TAC
la cantidad (FR) se puede calcular
con el área del rectángulo
CONT…
Definiciones de
Selectividad instantánea
SB/XY
0
0,1
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
0
0,02
0,04
0,06
0,08
0,1
0,12
0,14
0,16
0,18
0,2
0,22
0,24
0,26
0,28
0,3
0,32
0,34
0,36
0,38
0,4
0,42
CA
CA=0.0397
Cb=0.1628
Cx=0.107
Cy=0.091
Cómo es la curva de rendimiento?
YB
Cómo podemos usar? PFR
TAC
TAC
ß CA
Ejemplo 2:
Considerar el siguiente sistema de reacciones en fase acuosa:
Calcular el valor máximo de Cs que puede obtenerse en una operación a temperatura constante
con CAo de 2.
a) En un reactor CSTR
b) En un reactor flujo pistón
c) Proponga un esquema de dos reactores que maximice la concentración de S en
la salida.
Ejemplo 3:
Considerar las siguientes reacciones en fase acuosa:
SEMESTRE I-2024
En un reactor tubular
flujo pistón (PFR)
Reactor tubular flujo pistón (PFR)
τ’= W/vo
A à B àC
Balance estequiométrico:
Ecuaciones de velocidad:
FA + FB + FC = FAo
(-rA)= k1CA
Si v=vo:
rB,NETA= k1CA – k2CB CA + CB + CC = CAo
CA=CAo 𝑒 !"!τ
Reemplazando CA:
Esta ecuación tiene la forma:
dCB/dτ + k2CB= k1 CAo 𝑒 !"!τ y’ + P(x) y = Q(x)
CB
τ
Ø En la figura se muestra la solución a las ecuaciones diferenciales.
Ø Se obtienen las mismas funciones en reactores batch a volumen constante
(usando t en lugar de τ ).
Ø Comportamiento típico: la existencia de un máximo del producto intermedio.
Ø Cc se puede hallar con la ecuación diferencial 3 o del balance
estequiométrico.
Ø Se puede conocer el máximo derivando la función CB=CB(τ)
Rendimiento óptimo:
dCB/dτ = 0
Para C: -FC + rC V =0
CB,max
CB
τop
Ejercicios:
SEMESTRE I-2024
rE,NETA = …… rE,NETA = ……
Ecuaciones de velocidad:
r1A = k1A f(CA,CB,….)
aA + bB à cC + dD
(−𝑟𝐴) (−𝑟𝐵) 𝑟𝐶 𝑟𝐷
= = =
𝑎 𝑏 𝑐 𝑑
donde: FT= FA + FB + FC + ….
CTo = Po/RTo
%!
Para T=To y P=Po : Cj= Cto(% )
$
• Cambio de volumen:
PV = zNTRT
PV = zNTRT PoVo zo NToRTo
P cte !!!
v = vo (1 + e xA)
T cte !!!
Balance molar (EC. de Diseño):
En un reactor tubular
flujo pistón (PFR)
Reactor tubular flujo pistón (PFR)
τ’= W/vo
Reactor Tanque Agitado Continuo (CSTR)
Para C: -FC + rC V =0
REACCIONES MULTIPLES – PARTE 4
τ’= W/vo
W = 𝜌𝑏 V
Balances molares- Ecuaciones de diseño
- P=Po y T=To
Calculamos algunos
valores necesarios
Formulación de las EDO
CH
CM
Cx
Reactor Tanque Agitado Continuo (CSTR)
Para X: -FX + rX V =0
τ = 0.1 τ= 0.2 τ= 0.5 τ= 1
Comentarios:
Con cuál reactor se alcanza una mayor selectividad global SX/T (para τ= 0.5) ?
Ejercicio reacciones complejas 2
Serie- paralelo
Calcular a 673 oC:
Balances molares
Estequiometría y datos:
Ecuaciones diferenciales
a resolver:
Balances molares
Dada la reacción en fase líquida A à B que se lleva a cabo en un reactor tubular flujo pistón (PFR).
Calcular el volumen para alcanzar el 70% de conversión si la reacción es exotérmica y adiabática.
Cinética
(k en función de T)
Estequiometría
Ecuación de diseño con los datos
BALANCE DE ENERGIA
Velocidad de reacción es función de la composición y la temperatura:
xA Equilibrio xA Equilibrio
T T
Ejemplo : Reacción irreversible de primer orden
AàB
−𝑟𝐴 = k1CA
𝑘1 = 𝑘1𝑜 𝑒 "$!⁄%&
−𝑟𝐴 = 𝑘1𝑜 𝑒 "$!⁄%& CAo (1-xA)
CA= CAo (1-xA) a volumen constante
Ejemplo: Reacción reversible primer orden
AàB
−𝑟𝐴 = k1CA – k2CB. ; CBo=0
𝑘1 𝐶𝐵𝑒𝑞
= = 𝐾𝑐
𝑘2 𝐶𝐴𝑒𝑞
𝐾𝑐
𝑥Ae =
𝐾𝑐 + 1
EXOTERMICA ENDOTERMICA
Equilibrio en reacciones químicas
(dGt)T,P= 0
Gt
eeq
Kp, Kc…
ØCómo afecta la temperatura?
xA xA
no
ico iso
é rm té
rm
t
iso ico
n o
T T
EXOTERMICA ENDOTERMICA
CANTIDAD DE CALOR
BALANCE DE ENERGIA
PERFIL DE TEMPERATURA
Q= DH
SI Q=0 (PROCESO ADIABATICO) DH = 0
Hent Hsal
DH = Hsal-Hent
Ejercicios :
2.- Construir las curvas cinéticas del Ejemplo 9-3 del Levenspiel. La de equilibrio y únicamente
dos curvas cinéticas con los datos necesarios para ver la forma de las curvas.
DISEÑO DE REACTORES NO ISOTERMICOS - 2
En estado estacionario
Q= DH Sin trabajo de eje
DH = Hsal-Hent
Entalpía en función de la temperatura:
en general Cpi=Cpi(T)
Ejercicio 2:
H2 - 6.52
HReactivos DH HProductos
a Tent a Tsal
DHoR DHoP
H2 0 29.30
Ejercicio 4:
Se quema por completo CO a 225 oC y 1 at con 90% de exceso de aire seco que está a
525 oC. Los productos de la combustión salen de la cámara de combustión a 1025 oC.
Calcule el calor generado por kmol de CO que se alimenta. (CO +1/2 O2 à CO2)
Ejercicio 5:
Calcule la máxima temperatura que se puede alcanzar durante la combustión de
CO a 100 oC (iniciales) con 100% de exceso de aire. Suponga que la compustión es
completa.
DHf,298(kJ/kmol) Cp (kJ/kmol C)
o
CO -110600 30.35
O2 0 36.00
N2 0 29.30
OPERACIÓN ADIABATICA
Q=0
(8-29) xA
Balance de energía
To T
Ecuación Lineal
(8-30)
Balance de energía – Proceso adiabático
Forma útil para el diseño junto
al balance molar
&'!'('
= (−𝑟𝐴) Balance molar
')
PROCEDIMIENTO DE SOLUCION PARA UN REACTOR TUBULAR FLUJO PISTON
Concentraciones:
&(
Para T=To y P=Po : Cj= Cto(& )
)
PROCEDIMIENTO POR INTEGRACION
𝐹𝐴𝑜
(−𝑟𝐴)
xA
DISEÑO DE REACTORES NO ISOTERMICOS – 3
𝐹𝐴𝑜
(−𝑟𝐴)
xA
OTRO PROCEDIMIENTO
Ejercicio:
Se va llevar a cabo la conversión del n-butano a iso-butano en fase líquida a alta presión
usando un catalizador líquido que da una velocidad específica de reacción de 31.1 h-1 a 360 oK.
n butano à i butano
Calcule el volumen del reactor tubular flujo pistón (PFR) necesario para procesar
100000gal/día (163 kmol/h) con una conversión del 70%. La alimentación está compuesta por
90% de n-butano y 10% mol de un inerte a 330 oK.
El proceso es adiabático.
Datos adicionales:
DHoRx= -6900 J/mol n butano, E=65.7 kJ/mol, Kc= 3.03 a 60 oC, CAo= 9.3 kmol/m3
CpnB= 141 J/mol.K, CpiB= 141 J/mol.K, CpI (inerte) = 161 J/mol.K
SOLUCIÓN AàB
Ecuación de diseño
Ecuación de velocidad
Constantes en función
de la temperatura
Constante de equilibrio
Fase líquida
(no hay cambio de densidad)
Balance de energía
PROCESO ADIABATICO Q=0
T=T(XA)
REEMPLAZAR DATOS
Y
EVALUACION DE PARAMETROS
Solución por Ecuaciones Diferenciales Ordinarias
(Solución en Polymath)
Conversión en función del Volumen
xA= 0.70 para V=2.5 m3
xA= 0.70 para V=2.5 m3
Velocidad de desaparición de A en función del Volumen
Temperatura en función del Volumen
Conversión de equilibrio de A en función del Volumen
Cálculo para un reactor Tanque
Agitado Continuo (TAC, CSTR)
Cálculos:
𝐹𝐴𝑜
(−𝑟𝐴)
𝑭𝑨𝒐 𝑭𝑨 𝒐
vs. xA
−𝒓𝑨 #𝒓𝑨
xA
xA
Reactor Tubular Flujo Pistón con cambiador de calor
Reemplazando Fi en función de
la conversión
Se lleva a cabo la desintegración catalítica en fase vapor de acetona a ceteno y metano con una cinética de
primer orden con respecto a la acetona con una velocidad especìfica de reacción dada por:
ln k = 34.34 -34 222 /T (k en s-1 y T en oK) .
CH3COCH3 à CH2CO + CH4
Se desea alimentar 7850 kg/h de acetona a un reactor compuesto por 1000 tubos de una pulgada cédula 40.
La temperatura y presión son 1035 oK y 162 kPa respectivamente.
Ecuación de velocidad
Variación de volumen
k en función de T
CASO ADIABATICO
Balance de energía
V
(-rA)
V
CON INTERCAMBIADOR
DE CALOR
Balance de energía
Datos y cálculo de parámetros
x
V
(-rA)
V
CON INTERCAMBIADOR
DE CALOR
Balance de energía
Balance al fluido
refrigerante
V
Ta
V
Fogler p. 504
DISEÑO DE REACTORES NO ISOTERMICOS – 3.1
Se lleva a cabo la desintegración catalítica en fase vapor de acetona a ceteno y metano con una cinética de
primer orden con respecto a la acetona con una velocidad especìfica de reacción dada por:
ln k = 34.34 -34 222 /T (k en s-1 y T en oK) .
CH3COCH3 à CH2CO + CH4
Se desea alimentar 7850 kg/h de acetona a un reactor compuesto por 1000 tubos de una pulgada cédula 40.
La temperatura y presión son 1035 oK y 162 kPa respectivamente.
Ecuación de velocidad
Variación de volumen
k en función de T
CASO ADIABATICO
Balance de energía
V
(-rA)
V
CON INTERCAMBIADOR
DE CALOR
Balance de energía
Datos y cálculo de parámetros
x
V
(-rA)
V
CON INTERCAMBIADOR
DE CALOR
Balance de energía
Balance al fluido
refrigerante
V
Ta
V
Fogler p. 504
REACTOR TANQUE AGITADO CONTINUO CON EFECTOS DE CALOR
Ec. 8-26
Reemplazando el balance
molar : Ec. 8-42
Del balance de energía CASOS:
Balance molar
Caso 3:
T
- Ejemplo 8.8
- Ejemplo 8.9
REACCIONES MULTIPLES EN UN REACTOR TUBULAR FLUJO PISTON
REACCIONES MULTIPLES EN UN REACTOR TANQUE AGITADO CONTINUO
EJEMPLOS
8-10
8.11
CATALISIS Y REACTORES CATALITICOS
Semestre I/2024
!"!
- = (-r’A) [mol/kg .s]
!#
DEFINICIONES
Fogler
Catálisis heterogénea:
Incluye más de una fase generalmente
una de ellas es el catalizador en fase
sólida. Transferencia de masa.
Fogler
PROPIEDADES
v Catalizador típico sílica-alúmina: área superficial 300 m2/g, radio del poro 4
nanómetros, volumen del poro 0.6 cm3/g.
• Ejemplos: Niquel Raney, Platino sobre alúmina, hierro
v Tamices moleculares: poros más pequeños, admiten sólo algunas moléculas.
• De origen natural como ciertas arcillas o artificiales –aluminosilicatos cristalinos.
Pueden ser altamente selectivos.
v Monolitos: Estructuras tipo cerámicas. Panal.
Ø Catalizadores soportados.- Ejemplo: platino sobre alúmina.
Ø Catalizadores no soportados.– Ejemplo: hierro (síntesis de amoníaco)
PASOS DE UNA REACCION CATALITICA
Bird
Fogler
PASOS DE UNA REACCION CATALITICA
1 7
2 6
4 3
5
CAb
𝛿
CAs
Si es la etapa más lenta: Velocidad= (-rA)’=kc (CAb – CAs)
Si no hubiera resistencia a la
transferencia de masa
S : Sitio activo
A.S : Especie A adsorbida A + S à A.S
(átomo, molécula, combinación de átomos)
Ct : Concentración molar total de sitios activos por unidad de masa (mol/g catalizador) Se asume
(número total de sitios dividido por el número de Avogadro) constante
Cv : Concentración molar de sitios vacíos por unidad de masa (mol/g catalizador)
(número total de sitios vacíos dividido por el número de Avogadro)
Pi : Presión parcial de la especie i en fase gaseosa (atm o kPa)
pA pA pA pA
Ø Toda la superficie tiene la misma actividad.- enérgicamente uniforme
Adsorción de un compuesto A:
S
A+S ⇔ A.S
rAD = kA PA Cv
ISOTERMA
pA
A
Adsorción de un compuesto A: S
A+S ⇔ A.S
rAD = kA PA Cv
pA
pA
Adsorción disociativa de un compuesto A2:
A2
A2 + 2S ⇔ 2A.S
A
A
rAD = kA PA2 Cv2 S S
Se puede linealizar
ISOTERMA
Ejercicio:
a) Adsorción molecular
b) Adsorción disociativa
Fogler 10.2.3
Reacción superficial:
Sitio único
Sitio dual
Otro ejemplo de sitio
dual
MECANISMO PROPUESTO
Eliminar las
concentraciones Ci.s
Ecuación de velocidad
RELACION DE CONSTANTES
ASUMIR QUE OTRA ETAPA ES
LA CONTROLANTE: ejemplo
la etapa de Adsorción
Pc
CATALISIS Y REACTORES CATALITICOS
I/2024
Ejemplo:
Una reacción en fase gaseosa catalizada por un sólido tiene la siguiente forma: A+B à C
Grafique las curvas de velocidad inicial ( a la conversión cero) en función de la presión total para
los siguientes casos:
En todos los casos asumir que las presiones iniciales de A y B son iguales.
Mecanismo para incisos
a, b y c
0.8
0.6
0.4
0.2
(-rA)’o
p
Ecuación de velocidad:
(d)
Pc
Tarea:
5) Para la reacción T + H à B + M
Encontrar una expresión de velocidad asumiendo el siguiente
mecanismo:
T + S ↔ T.S
T.S + H ↔ B.S + M (etapa controlante)
B.S ↔ B + S
Análisis cuantitavo de la ecuación de velocidad
- Opción 1:
Tratar de acomodar a una forma lineal.
- Opción 2:
Regresión no lineal
I/2024
Ø Cómo se caracterizan?
Ø No todas las moléculas, grupos de molélulas, burbujas, etc., pasan el mismo tiempo en
el reactor.
Ø Canalizaciones dentro del reactor.
Ø Zonas estancadas (muertas), cortocircuitos en el reactor.
Fogler
Se debe adoptar otra aproximación más práctica…
Fogler
Ø Tres conceptos para describir reactores no ideales:
üDistribución de tiempos de residencia.
üCalidad de mezclado.
üModelo empleado para describir el sistema.
Ø Ordenar conceptos:
ü Primero considerar flujo ideal: PFR o CSTR
ü Siguiente nivel: DTR (macromezcla)
ü Finalmente: Información de micromezclado para calcular la conversión.
FUNCION DE DISTRIBUCION DE TIEMPOS DE RESIDENCIA (DTR)
Fogler
Función de entrada: Pulso
C C
Ejemplo:
Función de entrada Escalón y
respuesta en un Reactor Flujo
Pistón 𝜏 t
0 0
t
DISTRIBUCION DE TIEMPOS DE RESIDENCIA (DTR) A PARTIR
DE MEDICIONES DE RESPUESTA
FUNCION ESCALON
C/Co
1
Función de
Entrada
0 t
C/Co
Función de 1
Salida
0 t
Co
0 t +∆𝑡
Curva C
0
Función de salida
E(t)
0 𝑡
Tiempo medio de residencia
- A densidad constante
- 𝜀 =0
Ejercicio:
I/2024
Función escalón:
REACTOR TUBULAR FLUJO PISTON:
Función pulso:
REACTOR TANQUE AGITADO CONTINUO:
Función escalón:
REACTOR TANQUE AGITADO CONTINUO:
Función pulso:
Partiendo de la curva F ya
Obtenida:
REACTORES REALES
DIAGNOSTICO DE FUNCIONAMIENTO
I/2024
Cola larga
IDEAL
BY PASS
(DERIVACION)
ZONAS
MUERTAS
REACTOR TANQUE AGITADO CONTINUO:
IDEAL
BY PASS
(DERIVACION)
ZONAS
MUERTAS
REACTORES REALES
CALCULO DE LA CONVERSION
I/2024
CONVERSION
(k=0.1 min-1)
Ejemplo13-5
FOGLER
MODELO DE DISPERSION:
Función escalón:
[Link]
Fogler 4ª Ed.
Sistema cerrado-cerrado(condiciones de borde):
MODELO DE TANQUES EN SERIE:
Función pulso:
J.M. Smith
Función pulso:
Fogler
REACTOR CON RECIRCULACION: