RQ Parte3
RQ Parte3
RQ Parte3
Groucho Marx
88
Daniel Rojas Bravo Reactores Químicos
Con este objeto, primero enfocaremos nuestra atención al análisis de los datos cinéticos
de un proceso batch y posteriormente los aplicaremos a un CSTR:
Operación Batch
Ya conocemos que en el modelamiento de un reactor batch, no existen entradas ni
salidas de flujos durante la reacción, y el mezclado es perfecto. Para muchas reacciones en
fase líquida, el cambio de la densidad con la reacción es muy pequeña, por lo que puede
despreciarse.
Consecuentemente, el balance de moles para un reactor batch es:
1 ⎛dN A ⎞
⎜ ⎟= r
V ⎝ dt ⎠ A
que puede representarse en función de la concentración
⎛⎜dC A ⎞⎟ = r
⎝ dt ⎠ A
Generalmente, para análisis experimentales del laboratorio, es mejor procesar estos
datos en términos de la variable medida. Por ejemplo, la concentración es una variable medida
muy común en reacciones en fase líquida, y generalmente la ecuación de balance molar se
aplica directamente a este tipo de reacciones porque no existe un cambio significante en el
volumen. En este caso, la ecuación queda:
dC A
− = −rA
dt
Por ejemplo, para una reacción química Inicialmente, CA=CAo para t=0. Si la reacción
que obedece a una cinética de segundo se realiza isotérmicamente, se puede
orden en la cual A esta siendo consumida, integrar para dar la concentración de
la ley de velocidad es: reactante en cualquier tiempo:
−rA = kCA2 C t
∫ ∫ dt
1 A dC A
la cual podemos combinar con la ecuación − =
kC CA 0
de balance de moles para obtener Ao
dC A 1⎛ 1 1 ⎞
− = kCA2 ⎜ − ⎟=t
dt k ⎝C A C Ao ⎠
dC Este valor es el tiempo necesario para
− A2 = dt
kCA reducir la concentración del reactante desde
89
Daniel Rojas Bravo Reactores Químicos
tt = t f + t e + t c + t
En la mayoría de los casos el tiempo de reacción representa una fracción pequeña del
tiempo del ciclo total. Por ejemplo, en la tabla 4-1, se muestran los diferentes tiempos para un
ciclo de operación en un reactor batch de un proceso de polimerización. “Los tiempos de
polimerización pueden variar entre 5 a 60 horas*". Es claro, que el tiempo de reacción con una
reacción de 60 horas es un problema crítico. Entonces si se reducen los tiempos de reacción
toma importancia el uso de líneas grandes y de bombas para alcanzar una transferencia rápida
y disminuir el tiempo del ciclo.
TABLA 4-1: TIEMPOS DE CICLO TIPICOS PARA UN PROCESO DE POLIMERIZACIÓN BATCH.
ACTIVIDAD TIEMPO (h)
1.-Alimentación de la carga al reactor y agitado 1.5-3.0
2.- Calentamiento hasta la temperatura de reacción 1.0-2.0
3.- llevar a cabo la reacción (varía)
4.- Vaciar y limpiar el reactor 0.5-1.0
___________
Tiempo total excluyendo la reacción 3.0-6.0
DATOS TOMADOS DEL LIBRO “Elements of Chemical Reactor Engineering”, pagina 110, Scott Fogler, Prentice Hall, (1992).
90
Daniel Rojas Bravo Reactores Químicos
Solución:
∫ ∫ kdt
V dt dC A
− =
CA
C Ao 0
2.- La ley de velocidad es
−rA = kCA Aquí estamos usando las condiciones de que
La expresión de velocidad toma esta CA=CAO cuando t=0. La concentración inicial del
forma, porque como el agua está en exceso para componente A después de mezclar los dos
cualquier tiempo, esto hace que la velocidad de volúmenes es 1.0 kmol/m3 (1 mol/litro).
reacción sea independiente de la concentración
del agua (CB=CBo) 5.- Integrando, genera:
C Ao
ln = kt
3.- Estequiometría: En una reacción en fase CA
líquida los cambios de volumen se pueden Por lo tanto, la concentración de óxido de
despreciar (V=Vo): etileno a cualquier tiempo es:
CA =
NA NA
=
C A = C Ao e−kt
V Vo
Con estos cálculos, la concentración del etilenglicol a un tiempo t se obtiene mediante el uso de la
estequiometría de la reacción.
A+ B⎯
⎯→ C
TABLA E-4.1: Tabla estequiométrica
Especie Símbolo Cantidad Cambio Remanente Concentración
Inicial
CH2CH2O A NAo −N X Ao
N A = N Ao (1 − X) C A = C Ao (1 − X)
H 2O B Θ B N Ao − N Ao X N A = N Ao (Θ B − X) C B = C Ao (Θ B − X)
C B = C AoΘ B = C Bo
(CH2OH)2 C 0 N AoX N A = N Ao X C C = C Ao X
NTo NT=NTo-NAo X = C Ao - C A
91
Daniel Rojas Bravo Reactores Químicos
C]
C o- CA ) / A
-1
C C = C Ao − C A = C Ao (1− e−kt )
ln[(A
rearreglando y tomando logaritmos de ambos
lados, produce: -2
C − CC
ln Ao = −kt
C Ao
Esta ecuación se puede graficar con las -3
92
Daniel Rojas Bravo Reactores Químicos
⎛ F X ⎞
V = v o ⎜ Ao ⎟
⎝(−rA ) sale ⎠
⎛C − C A ⎞
V = v o ⎜ Ao ⎟
⎝ −rA ⎠
−rA = kCA
que podemos combinar esta ley de velocidad y la ecuación de balance de moles para dar:
C Ao − C A
τ=
kCA
Resolviendo para la concentración del efluente de A, CA,
C Ao
CA =
1 + τk
Para el caso que estamos tratando, no existe un cambio de volumen durante el curso de
la reacción:
C Ao − C A
X=
C Ao
τk
X=
1 + τk
El producto τk se conoce como el Número de Damkohler.
93
Daniel Rojas Bravo Reactores Químicos
Una reacción de primer orden que serie para una reacción de primer orden se
no sufre cambio de volumen (v=vo) se esta muestra en la figura 4-3, con varios valores
realizando en una serie de dos reactores del número de Damkohler.
CSTR, como se muestra en la siguiente
figura. Observe lo que se había dicho
CAo anteriormente se cumple, cuando el
producto del tiempo espacial y la velocidad
de reacción son relativamente grandes, se
-rA2 dice que Da≥1, y aproximadamente el 90%
-rA1 V2 de conversión es alcanzada en los dos o
V1 tres primeros reactores. Así el costo de
otros reactores adicionales no puede ser
CA1 justificado.
X1 CA2
X2 Con valores pequeños del producto,
τ k , es decir, con Da≤0.1, la conversión
FIGURA 42. Dos reactores CSTR en serie.
La concentración del efluente de A desde el continúa incrementándose con cada reactor
reactor 1 es adicionado.
C Ao La velocidad de desaparición de A
CA 1 = en el n-ésimo reactor es:
1+ τ 1 k1
C Ao
Considerando la ecuación de balance sobre −rA n = kCA n = k
el reactor 2,
( τ k)n
1+
FA1 − FA 2 v o (C A 1 − C A 2 ) 1.1
V2 = = 1
−rA 2 k2 C A 2
0.9
resolviendo para la CA2, la concentración
0.8
que sale del segundo reactor es:
Conversión, X
C A1 CA o 0.7
CA 2 = =
1+ τ 2 k2 (1+ τ 2 k2 )(1+ τ 1 k1) 0.6
0.5
Si en lugar de dos reactores conectados Da=0.1
0.4
en serie, tenemos n reactores conectados Da=0.5
Da=1.0
en serie (τ 1 = τ 2 = τ 3 = .......= τ n ) operando a la 0.3
Da=10
0.2
misma temperatura, la concentración
alcanzada en el último reactor será: 0.1
C Ao C Ao 0
CA n = = 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
(1+ τ k ) (1+ Da ) n
n Número de tanques, n
la conversión para el n-ésimo tanque es: FIGURA 43. Conversión como una función del número de
1 tanques en serie para una reacción de primer orden
X = 1−
(1+ τ k ) n
Una gráfica de la conversión como
una función del número de reactores en
94
Daniel Rojas Bravo Reactores Químicos
95
Daniel Rojas Bravo Reactores Químicos
vo X pies 3 0.8
V= = 15.34 197.3pies 3
k(1− X) min (0.311/ min)(1− 0.8)
= 1480gal = 5.6m 3
Un tanque de 5 pies de diámetro y de aproximadamente 10 pies de altura, es el necesario para
alcanzar el 80% de conversión.
AHORA USAREMOS LOS ARREGLOS DE REACTORES, LLEVANDO AL MISMO TIEMPO LOS DOS
ARREGLOS:
Para arreglo de reactores CSTRs en paralelo. Para arreglo de reactores CSTRs en serie.
Tomando dos reactores CSTRs de 800 gal arreglados en Si ahora, los dos reactores de 800gal son arreglados en
paralelo con 7.67 pies3/min (vo/2), la conversión que se serie. la conversión en el primer reactor es:
alcanza se evalúa con: τ1 k
t k X1 =
X= 1+ τ 1 k
1 +t k en donde
de donde V1 ⎛ 1pie 3 ⎞ 1
V ⎛ 1 pie ⎞⎛
3
1 ⎞ τ1 = = ⎜800gal ⎟x = 6.97 min
τ= ⎜800 galx ⎟⎜ ⎟ = 13.94 min v o1 ⎝ 7.48 gal ⎠ 15.34 pie3 / min
vo / 2 ⎝ 7.48 gal ⎠⎝7.67 pie 3 / min⎠
El número de Damkohler es,
De aquí podemos evaluar el número de Damkohler 0.311
0.311 Da = τ 1k = 6.97min = 2.167
Da = τk = 13.94 min = 4.34 min
min De aquí que la conversión alcanzada en el primer reactor es,
Sustituyendo en la ecuación de la conversión, se obtiene;
2.167
4.34 X1 = = 0.684
X= = 0.81 3.167
1+ 4.34
Para calcular le conversión de salida en el
que representa la conversión alcanzada
segundo reactor, recordaremos que V1=V2=V y que
v o1 = v o2 = v o , entonces τ 1 = τ 2 = τ .
Combinando la ecuación de balance de moles, en
el segundo reactor, con la ley de velocidad , se obtiene:
7 .6 7 pie3 /m in
F (X − X1 ) CAo v o (X 2 − X1 ) v o ⎛X 2 − X1 ⎞
V = Ao 2 = = ⎜ ⎟
−rA 2 kCAo (1 − X 2 ) k ⎝ 1− X2 ⎠
Como ya conocemos todos los valores de los
diferentes parámetros, podemos resolver para la conversión
15 .3 4 pie3 /m in
X2 (que es la conversión alcanzada en el segundo reactor):
X1 + τk 0.684 + 2.167
X2 = = = 0.90
X=0.81 1+ τ k 1+ 2.167
o bien, se obtiene el mismo resultado con la ecuación que
habíamos obtenido anteriormente para un sistema
reaccionante con una cinética de primer orden:
1 1
X = 1− = 1− = 0.90
(1+ τ k) n
(1+ 2.167) 2
7 .6 7 pie3 /m in
vo=15.34 pies3/min
(7.67 dm3/s)
X=0.8 1
X1 =0.68 X2=0.90
96
Daniel Rojas Bravo Reactores Químicos
Evaluando la producción para comprobar la veracidad de Evaluando la producción para comprobar la veracidad
nuestros cálculos: de nuestros cálculos:
lb lb
= 199681148.16 = 221867941
año año
Por lo tanto, 200 millones de libras de etilenglicol por año se pueden producir usando dos reactores
arreglados en serie, de 800 gal cada uno.
NOTA: Como siempre, dos reactores conectados en serie darán más alta conversión que dos reactores
conectados en paralelo, del mismo tamaño y cuando el orden de la reacción es mayor que cero.
97