Distribución de Tiempos de Residencia y Modelos para Reactores No Ideales
Distribución de Tiempos de Residencia y Modelos para Reactores No Ideales
Distribución de Tiempos de Residencia y Modelos para Reactores No Ideales
Después de estudiar la DTR, diga ¿qué problemas se están presentando en el reactor durante el
período de operación deficiente (línea delgada)?. El lecho se volvió a empacar y se hizo de nuevo
la prueba con trazador con los resultados mostrados en la figura P-13-5 (línea gruesa).
La variación de HDS ha sido simulada como una reacción de primer orden, de acuerdo a las DTR
presentadas, ¿bajo qué condiciones piensa que se obtendría la mayor conversión?, ¿cuál sería el
caso óptimo?
Respuesta:
La línea de operación deficiente de la DTR de este reactor, presenta un pico máximo de salida
de fluido en un corto tiempo y posteriormente, una salida del resto de fluido en un tiempo
mucho mayor formando una cola larga, lo que es indicativo de canalización y zonas muertas
dentro del lecho del reactor, respectivamente. La canalización genera una salida rápida de
gran parte del fluido en corto tiempo y las zonas muertas hacen que el fluido permanezca
dentro del reactor un tiempo mayor al tiempo de residencia.
El reactor reempacado presenta una DTR más acorde con el tipo de reactor, sin embargo,
también se observa una cola bastante larga, que indica la presencia de zonas muertas dentro
del reactor.
Si tomamos en cuenta que el reactor en operación deficiente presenta una canalización del
fluido lo que causa que parte de la alimentación salga del reactor es un tiempo de residencia
menor al tiempo para el cual fue diseñado, entonces, podemos concluir que el reactor con
lecho reempacado debería obtener una mayor conversión. Sin embargo, si bien el lecho
reempacado presenta una mejor DTR, también presenta zonas muertas. A modo de aclarar un
poco lo antes explicado se calculó el tiempo medio a partir de cada DTR, de la Figura 1
podemos calcular el área bajo la curva para obtener los valores de función de distribución de
tiempos de residencia. Esta evaluación, se realizó mediante la utilización de Polymath 5.0 con
los datos digitalizados por medio del programa Engauge Digitizer 4.0. Posteriormente, con los
valores de t y E(t) se calcula el tiempo medio, a través de:
2
C (t )
E (t ) = ∞
Ec. 1
∫ C(t )dt
0
∞
tm = tE (t )dt
∫ Ec. 2
0
P-13-5 Digitalizada
90 500
80 450
400
70
350
60
300
50
C (t)
250
40
200
30
150
20
100
10 50
0 0
0 10 20 30 40 50 60
t (min)
Operación Deficiente Lecho Reempacado
3
80 500
450
70
400
60
350
50
300
C (t)
40 250
200
30
150
20
100
10
50
0 0
0 5 10 15 20 25 30 35 40 45 50 55 60
t (min)
Operación Deficiente tm (Operación Deficiente)
Lecho Reempacado tm (Lecho Reempacado)
Como se observa en la Figura 2, la operación deficiente presenta una salida de fluido mucho
más adelantada que la salida de fluido de el lecho reempacado (primer pico para ambas), sin
embargo, observamos que en el lecho reempacado el fluido también sale antes del tiempo
medio, lo que indica que sigue existiendo canalización del fluido. También observamos, que en
ambos casos existe una cola larga con un tiempo medio de hasta alrededor de 3.5 veces el
tiempo medio, confirmando la existencia de zonas muertas.
Un caso óptimo sería aquel en el que al reempacar el reactor la DTR presentara un tiempo
medio simétrico, esto puede lograrse si: el reempacado del reactor es tal que deje suficiente
espacio entre las partículas para que el fluido circule (evitar zonas muertas) y si se colocan
distribuidores que eviten que el lecho se canalice.
4
Operación Deficiente
-1
t (min) C(t) ∫C(t)dt E(t) (min ) tE(t) tm
2.71 0.48 632.35 7.67E-04 2.08E-03 15.73
2.77 20.56 3.25E-02 9.01E-02
2.90 36.62 5.79E-02 1.68E-01
3.04 44.12 6.98E-02 2.12E-01
3.08 49.21 7.78E-02 2.40E-01
3.24 63.40 1.00E-01 3.25E-01
3.31 75.18 1.19E-01 3.94E-01
3.53 77.61 1.23E-01 4.33E-01
3.66 77.62 1.23E-01 4.49E-01
4.11 73.91 1.17E-01 4.80E-01
4.35 58.68 9.28E-02 4.04E-01
4.36 50.12 7.93E-02 3.46E-01
4.53 39.96 6.32E-02 2.86E-01
5.03 32.51 5.14E-02 2.59E-01
6.15 27.25 4.31E-02 2.65E-01
8.30 22.34 3.53E-02 2.93E-01
11.47 18.58 2.94E-02 3.37E-01
14.33 17.74 2.81E-02 4.02E-01
17.97 16.16 2.56E-02 4.59E-01
21.21 15.62 2.47E-02 5.24E-01
24.86 13.77 2.18E-02 5.41E-01
28.12 12.16 1.92E-02 5.41E-01
31.01 9.71 1.54E-02 4.76E-01
35.08 5.49 8.68E-03 3.05E-01
39.26 3.15 4.98E-03 1.96E-01
44.46 1.70 2.68E-03 1.19E-01
47.70 0.89 1.40E-03 6.68E-02
50.28 0.14 2.15E-04 1.08E-02
60.00 0.00 0.00E+00 0.00E+00
5
Tabla 2. Lecho Reempacado
Lecho Reempacado
-1
t (min) C(t) ∫C(t)dt E(t) (min ) tE(t) tm
0.00 0.00 5704.59 0.00E+00 0.00E+00 14.21
1.16 3.05 5.35E-04 6.22E-04
1.68 6.06 1.06E-03 1.78E-03
2.32 10.57 1.85E-03 4.30E-03
2.96 21.04 3.69E-03 1.09E-02
3.33 40.43 7.09E-03 2.36E-02
3.70 56.84 9.96E-03 3.69E-02
4.07 77.72 1.36E-02 5.54E-02
4.31 94.87 1.66E-02 7.17E-02
4.55 109.78 1.92E-02 8.76E-02
4.73 129.16 2.26E-02 1.07E-01
4.97 147.80 2.59E-02 1.29E-01
5.33 176.87 3.10E-02 1.65E-01
5.99 229.80 4.03E-02 2.41E-01
6.83 300.61 5.27E-02 3.60E-01
7.56 351.31 6.16E-02 4.65E-01
8.10 393.05 6.89E-02 5.58E-01
8.45 425.85 7.47E-02 6.31E-01
8.90 450.45 7.90E-02 7.03E-01
9.13 465.37 8.16E-02 7.45E-01
9.50 477.00 8.36E-02 7.94E-01
9.76 480.00 8.41E-02 8.21E-01
10.41 480.35 8.42E-02 8.76E-01
10.81 469.94 8.24E-02 8.90E-01
11.21 454.33 7.96E-02 8.93E-01
11.68 432.01 7.57E-02 8.85E-01
13.05 358.35 6.28E-02 8.20E-01
13.95 300.31 5.26E-02 7.34E-01
14.56 266.83 4.68E-02 6.81E-01
15.24 240.05 4.21E-02 6.41E-01
16.18 202.87 3.56E-02 5.75E-01
16.78 182.79 3.20E-02 5.38E-01
18.17 154.57 2.71E-02 4.92E-01
19.55 132.31 2.32E-02 4.53E-01
21.12 111.55 1.96E-02 4.13E-01
22.76 93.03 1.63E-02 3.71E-01
24.39 76.75 1.35E-02 3.28E-01
25.70 64.91 1.14E-02 2.92E-01
27.72 50.14 8.79E-03 2.44E-01
30.00 37.62 6.59E-03 1.98E-01
32.79 25.87 4.54E-03 1.49E-01
35.06 19.31 3.38E-03 1.19E-01
37.27 15.72 2.76E-03 1.03E-01
39.79 12.90 2.26E-03 9.00E-02
41.87 10.05 1.76E-03 7.37E-02
45.82 5.82 1.02E-03 4.68E-02
48.09 3.73 6.54E-04 3.14E-02
49.90 3.10 5.43E-04 2.71E-02
51.26 1.69 2.97E-04 1.52E-02
52.62 0.25 4.42E-05 2.33E-03
56.44 0.20 3.53E-05 1.99E-03
60.00 0.00 0.00E+00 0.00E+00
6
Problema 14-11
Ocurre una reacción irreversible de segundo orden en un CSTR isotérmico, pero no ideal. El
volumen del reactor es de 1000 dm3 y la velocidad de flujo de la corriente de reactivo es de
1 dm3/s. A la temperatura del reactor, k = 0.005 dm 3 mol ⋅ s . La concentración de A en la corriente
de alimentación es de 10 mol/dm3. Se obtiene la DTR mediante una prueba con trazador en este
reactor a la velocidad de alimentación y temperatura de reacción deseadas. A partir de la DTR en
cuestión:
(a) Estime (1) las conversiones de segregación máxima y (2) segregación mínima (mezclado
máximo) que se pueden obtener en este reactor debido a las diferentes condiciones de
micromezclado que puede haber.
(c) ¿Qué conversiones se puede obtener modelando este sistema como dos CSTR en paralelo?
(d) Compare estas conversiones con las de un CSTR perfectamente mezclado y un PFR del mismo
volumen
7
Respuesta:
(a.1) En este caso tomamos la ecuación de un reactor por lotes y obtenemos la expresión de la
conversión en función del tiempo, tomando en cuenta que la reacción es de segundo orden.
dX
N A0 = −r AV Ec. 3
dt
− r A = kC A2 Ec. 4
kC A0 t
X (t ) = Ec. 5
1 + kC A0 t
Con esta expresión podemos entonces calcular la conversión promedio utilizando la función de
distribución de tiempos de residencia
∞
X = ∫ X (t )E (t )dt
o
Ec. 6
Para evaluar la integral se halla el área bajo la curva de la gráfica X(t)E(t) vs t, que se presenta en
la Figura 3. En la Tabla 3, se muestran los valores calculados.
8
2.0E-03
1.5E-03
X(t)E(t)(s )
-1
1.0E-03
5.0E-04
0.0E+00
0 2000 4000 6000 8000 10000 12000 14000 16000 18000 20000
t(s)
9
Tabla 3. Conversión Asumiendo Segregación Máxima
Segregación Máxima
-1 -1
t (s) E(t) (s ) X(t) X(t)E(t)(s ) X promedio
0 3.250E-03 0 0 0.869
5 3.187E-03 0.200 6.374E-04
10 3.124E-03 0.333 1.041E-03
25 2.945E-03 0.556 1.636E-03
40 2.776E-03 0.667 1.851E-03
70 2.468E-03 0.778 1.920E-03
100 2.194E-03 0.833 1.828E-03
175 1.637E-03 0.897 1.469E-03
250 1.224E-03 0.926 1.133E-03
325 9.184E-04 0.942 8.652E-04
400 6.913E-04 0.952 6.584E-04
700 2.366E-04 0.972 2.300E-04
1000 9.755E-05 0.980 9.564E-05
2500 2.691E-05 0.992 2.670E-05
4000 1.839E-05 0.995 1.830E-05
7000 8.689E-06 0.997 8.664E-06
10000 4.104E-06 0.998 4.096E-06
15000 1.176E-06 0.999 1.174E-06
20000 3.369E-07 0.999 3.366E-07
dX E (λ )
= −kC A0 (1 − X )2 + X Ec. 7
dλ 1 − F (λ )
Para resolver el balance se utilizan deltas finitos, lo que convierte la Ec. 7 en:
E (λ )
∆X = ∆λ − kC A0 (1 − X )2 + X Ec. 8
1 − F (λ )
Para el cálculo se asume que cuando λ 0 = 100 entonces X = 0 , los valores intermedios de E(t) se
calculan por interpolación, se toma ∆λ = −10 , se debe acotar que no se toma el tiempo máximo ya
que valores de ∆λ superiores pueden hacer la que la expresión salte a valores de conversión
desproporcionados. En la Tabla 4, se presentan los resultados para mezclado máximo,
observándose que la conversión obtenida mediante este modelo es de 0.766, y como es de
esperarse es menor que la conversión obtenida para el modelo de segregación total para este
orden de reacción.
10
Tabla 4. Conversión Asumiendo Mezclado Máximo
Mezclado Máximo
-1 -1
∆λ (s) λ0 (s) λf (s) E(λ0 ) (s ) F(λ0) E(λ0 )/(1-F(λ0)) (s ) X0 Xf
-10 100 90 2.194E-03 0.269 3.000E-03 0.000 0.500
90 80 2.285E-03 0.246 3.032E-03 0.500 0.610
80 70 2.377E-03 0.223 3.059E-03 0.610 0.667
70 60 2.468E-03 0.199 3.081E-03 0.667 0.702
60 50 2.571E-03 0.174 3.111E-03 0.702 0.725
50 40 2.673E-03 0.147 3.136E-03 0.725 0.740
40 30 2.776E-03 0.120 3.156E-03 0.740 0.750
30 20 2.889E-03 0.092 3.181E-03 0.750 0.758
20 10 3.005E-03 0.063 3.205E-03 0.758 0.763
10 0 3.124E-03 0.032 3.227E-03 0.763 0.766
0 3.250E-03 0.000 3.250E-03 0.766 0.769
(b) Para obtener la conversión de tanques en serie, debemos calcular mediante la varianza, el
número de CSTR en serie del mismo volumen y utilizar la ecuación de diseño, tomando en cuenta
que es una reacción de segundo orden.
∫ (t − tm ) E (t )dt
2
σ2 = Ec. 9
0
∞
tm = tE (t )dt
∫ Ec. 10
0
τ2
n= Ec. 11
σ2
FA0 X = −r AV Ec. 12
C An =
(
− 1 ± n −1 1 + 2 − 1 ± 1 + 4 k τ i C A0 ) (Para n>1) Ec. 13
2kτi
11
− 1 ± 1 + 4 k τ i C A0
C A1 = (Para n=1) Ec. 14
2kτi
τ
τi = Ec. 15
n
C A0 −
(
− 1 ± n −1 1 + 2 − 1 ± 1 + 4 k τ C
A0 )
2 k τ Ec. 16
X = (Para n>1)
C A0
− 1 ± 1 + 4 k τ C A0
C A0 −
2kτ Ec. 17
X = (Para n=1)
C A0
3.50E-01
3.00E-01
2.50E-01
2.00E-01
tE(t)
1.50E-01
1.00E-01
5.00E-02
0.00E+00
0 2000 4000 6000 8000 10000 12000 14000 16000 18000 20000
t(s)
12
3500
3000
2500
2000
(t-tm) E(t)
2
1500
1000
500
0
0 2000 4000 6000 8000 10000 12000 14000 16000 18000 20000
t(s)
Figura 5. Varianza
Tanques en Serie
-1 2 2 2
t (s) E(t) (s ) tE(t) tm (s) (t-tm) (t-tm) E(t) ∫(t-tm) E(t)dt τ (s) n CA1 X
0 3.250E-03 0.00E+00 977 955352 3.10E+03 4.761E+06 1000 0.210 1.318 0.868
5 3.187E-03 1.59E-02 945603 3.01E+03
10 3.124E-03 3.12E-02 935903 2.92E+03
25 2.945E-03 7.36E-02 907106 2.67E+03
40 2.776E-03 1.11E-01 878758 2.44E+03
70 2.468E-03 1.73E-01 823413 2.03E+03
100 2.194E-03 2.19E-01 769868 1.69E+03
175 1.637E-03 2.86E-01 643879 1.05E+03
250 1.224E-03 3.06E-01 529141 6.48E+02
325 9.184E-04 2.98E-01 425653 3.91E+02
400 6.913E-04 2.77E-01 333415 2.30E+02
700 2.366E-04 1.66E-01 76962 1.82E+01
1000 9.755E-05 9.76E-02 510 4.97E-02
2500 2.691E-05 6.73E-02 2318247 6.24E+01
4000 1.839E-05 7.36E-02 9135984 1.68E+02
7000 8.689E-06 6.08E-02 36271458 3.15E+02
10000 4.104E-06 4.10E-02 81406932 3.34E+02
15000 1.176E-06 1.76E-02 196632722 2.31E+02
20000 3.369E-07 6.74E-03 361858512 1.22E+02
13
Como se observa, el n obtenido es menor a uno, debido a que el orden de reacción es distinto a
uno en este caso debemos aproximar al entero superior, por lo tanto, tendríamos un solo tanque
CSTR, con una conversión de X = 0.868.
Dos CSTR en paralelo tendrían cada uno un volumen V1 = V2 = V 2 , con un flujo de inicial
FA01 = FA0 2 = FA0 2 . Siendo que ambos reactores trabajan a las mismas condiciones de
temperatura entonces, r A1 = r A2 , lo que hace que sus conversiones sean iguales entre ellas, e
iguales a la conversión total de un CSTR de un volumen V, por lo tanto, la conversión de este caso
es igual a la anterior X = 0.868.
(d.1) CSTR
La conversión obtenida para este caso es la misma del caso anterior X = 0.868.
(d.2) PFR
dX
FA 0 = −r A Ec. 18
dV
k C A0 τ
X = Ec. 19
1 + k C A0 τ
14