Libro Nuevo Notas Del Curso de Diseño de Reactores
Libro Nuevo Notas Del Curso de Diseño de Reactores
Libro Nuevo Notas Del Curso de Diseño de Reactores
DE REACTORES QUIMICOS
DISENO
por
Antonio Flores T.
Departamento de Ciencias
Universidad Iberoamericana
Primavera 2001
E-mail: [email protected]
URL: 200.13.98.241/antonio
September 5, 2001
Contenido
I
REACTORES ISOTERMICOS
1 Balance de masa
1.1 Introduccion . . . . . . . . . . . . . . . . . . . .
1.2 Balance general de masa . . . . . . . . . . . . .
1.2.1 Sistema homogeneo . . . . . . . . . . . .
1.2.2 Sistema heterogeneo . . . . . . . . . . .
1.3 Velocidad de reaccion . . . . . . . . . . . . . . .
1.4 Midiendo el progreso de una reaccion . . . . . .
1.4.1 Extension de la reaccion . . . . . . . . .
1.4.2 Conversion . . . . . . . . . . . . . . . . .
1.4.3 Relacion entre la extension de la reaccion
. . . . . . . . . . . . .
. . . . . . . . . . . . .
. . . . . . . . . . . . .
. . . . . . . . . . . . .
. . . . . . . . . . . . .
. . . . . . . . . . . . .
. . . . . . . . . . . . .
. . . . . . . . . . . . .
y grado de conversion
2
2
2
3
5
6
7
7
8
9
2 Reactor Batch
2.1 Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.2 Ecuaciones de dise
no . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.3 Reacciones m
ultiples . . . . . . . . . . . . . . . . . . . . . . . . . . . .
10
10
10
13
.
.
.
.
19
19
20
22
25
4 Reactor Tubular
4.1 Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.2 Reacciones m
ultiples . . . . . . . . . . . . . . . . . . . . . . . . . . . .
29
29
34
5 Reactor SemiBatch
5.1 Introduccion . . . . . . . . .
5.2 Reactor Semibatch . . . . .
5.3 Reacciones m
ultiples . . . .
5.4 Operacion Semibatch-Batch
.
.
.
.
36
36
36
46
49
6 Relaciones Estequiom
etricas
6.1 Reactores de flujo discontinuo . . . . . . . . . . . . . . . . . . . . . . .
6.2 Reactores de flujo continuo . . . . . . . . . . . . . . . . . . . . . . . . .
53
53
59
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
7 Configuraci
on de Reactores
7.1 Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
7.2 Determinado el mejor arreglo de reactores para una cierta conversion .
64
64
68
8 Sensibilidad param
etrica
8.1 Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
8.2 Ecuaciones de sensibilidad parametrica . . . . . . . . . . . . . . . . . .
8.3 Evaluacion de sensibilidad parametrica . . . . . . . . . . . . . . . . . .
73
73
73
74
80
ii
Lista de Figuras
1.1
1.2
1.3
2
4
5
2.1
2.2
2.3
14
15
18
3.1
3.2
29
30
5.1
5.2
5.3
5.4
5.5
5.6
5.7
5.8
Reactor Semibatch. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Reactor Semibatch. . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Variacion del grado de conversion con el tiempo. . . . . . . . . . . . . .
Perfil de concentraciones. . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Perfil de concentracion en un reactor semibatch. . . . . . . . . . . . . .
Perfil de n
umero de moles en un reactor de operacion semibatch-batch.
36
40
43
44
44
45
48
51
6.1
densidad constante
. . . . . . . . . . .
densidad constante
. . . . . . . . . . .
densidad constante
. . . . . . . . . . .
( =
. . .
( =
. . .
( =
. . .
0) y
. . .
0) y
. . .
0) y
. . .
57
58
63
64
65
65
7.4
7.5
7.6
7.7
7.8
7.9
7.10
7.11
7.12
Arreglo tubular-tanque
. . . . . . . . . . . . .
. . . . . . . . . . . . .
. . . . . . . . . . . . .
. . . . . . . . . . . . .
. . . . . . . . . . . . .
. . . . . . . . . . . . .
. . . . . . . . . . . . .
. . . . . . . . . . . . .
agitado.
. . . . .
. . . . .
. . . . .
. . . . .
. . . . .
. . . . .
. . . . .
. . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
66
68
69
70
70
71
71
72
72
8.1
8.2
8.3
76
77
78
iv
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
Lista de Tablas
Parte I
REACTORES ISOTERMICOS
Captulo 1
Balance de masa
1.1
Introducci
on
1.2
Fn0
REACTOR
F1
F2
Fn
F lujo de masa
que entra
al reactor
M asa
+ consumida/
generada
F lujo de masa
que abandona
el reactor
2
M asa
acumulada
=
en el reactor
(1.2.1)
Cap. 1
Balance de Masa
F lujo de masa
que entra
= Fio
al reactor
M asa
consumida/ = Gi
generada
F lujo de masa
que abandona
el reactor
M asa
acumulada
en el reactor
= Fi
dNi
dt
masa
1
i
notese que los anteriores terminos (Fio , Gi , Fi y dN
) poseen unidades de tiempo
.
dt
Substituyendo dichos terminos en la ecuacion 1.2.1 obtenemos una expresion general
del balance de masa en el reactor,
dNi
= Fio + Gi Fi
(1.2.2)
dt
la forma de evaluar el termino Gi depende de si el reactor donde se efectua la reaccion
deseada es de propiedades homogeneas o heterogeneas.
1.2.1
Sistema homog
eneo
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 1
Balance de Masa
(1.2.3)
(1.2.4)
dependiendo del tipo de reactor en cuestion la ecuacion anterior puede tomar formas
diferentes seg
un se muestra a continuacion.
Reactor continuo tanque agitado
dNi
= Fio Fi + ri V
dt
(1.2.5)
dNi
= ri V
dt
(1.2.6)
Reactor batch
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 1
Balance de Masa
V2
V1
Vn
Reactor semibatch
dNi
= Fio + ri V
dt
1.2.2
(1.2.7)
Sistema heterog
eneo
Si las propiedades del sistema cambian de un punto a otro en el interior del reactor,
entonces la cantidad de masa formada y/o consumida no puede evaluarse mediante la
ecuacion 1.2.3. Aquellos sistemas de reaccion cuyas propiedades (p.e. concentracion,
temperatura, presion) cambian de un punto a otro del reactor se denominan sistemas de
reaccion de propiedades heterogeneas. Ejemplo de un sistema de reaccion heterogeneo
es el reactor tubular de flujo tapon.
Para derivar el termino Gi en sistemas heterogeneos suponemos que el volumen original puede dividirse en una serie de n subvolumenes tan peque
nos de manera que en
el interior de dichos subvolumenes las propiedades pueden considerarse homogeneas.
Por ejemplo para el sistema mostrado en la figura 1.3 el cual ha sido dividido en n
subvolumenes:
la cantidad de masa formada/consumida en cada volumen de reaccion Gi puede
evaluarse de la siguiente manera.
primer subvolumen: Gi1 = ri1 V1
segundo subvolumen: Gi2 = ri2 V2
UIA-Santa Fe
Cap. 1
Balance de Masa
n
*
(Gi )j =
j=1
n
*
rij Vj
(1.2.8)
j=1
+ V
(1.2.9)
ri dV
+ V
ri dV =
dNi
dt
(1.2.10)
1.3
Velocidad de reacci
on
Hemos antes mencionado que cuando se alcanzan las condiciones de equilibrio las velocidades de reaccion de los componentes A y B de la reaccion ?? son iguales: es decir
la velocidad neta de formacion es cero. La velocidad neta es simplemente la diferencia
entre las velocidades de reaccion de A y B.
Frecuentemente la velocidad de reaccion de un componente i se denota por ri . Las
moles
unidades de ri son (tiempo)(volumen)
. En el caso de la reaccion ?? tenemos dos velocidades
de reaccion: rA y rB . A la velocidad de reaccion de un reactivo (componente A) se
le antepone el signo negativo - para recordarnos precisamente que la velocidad de
reaccion es la de un reactivo. A la velocidad de reaccion de un producto (B) se le
antepone el signo positivo +.
Fsicamente rA denota las moles del reactivo A que han reaccionado en un cierto
intervalo de tiempo y en un volumen dado de volumen del sistema.
rB denota las moles del producto B que se han obtenido en un cierto intervalo
de tiempo si la reaccion se efectua en un volumen dado del sistema.
en cualquier caso debe notarse que, al menos al nivel de la definicion de velocidad
de reaccion, si uno modifica tanto el tiempo que se permite para que transcurra la
reaccion y/o el volumen del sistema disponible para efectuar la reaccion, estaremos en
la posibilidad de modificar la velocidad de reaccion algo que, ciertamente, es deseable
en muchos sistemas de reaccion.
Normalmente solo necesitamos una expresion de velocidad de reaccion por cada
reaccion qumica. Es decir, por ejemplo en el caso de la reaccion ??, no necesitamos
tener reportadas tanto rA como rB . Si tenemos la expresion matematica de rA
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 1
Balance de Masa
(1.3.11)
masa
(es
En el dise
no de reactores qumicos, nos interesa conocer no la (tiempo)(volumen
masa
2
ltimo termino puede
decir la velocidad de reaccion) sino la tiempo . Desde luego este u
obtenerse facilmente multiplicando la velocidad de reaccion por el volumen del sistema,
(1.3.12)
rA V
Formalmente, para un sistema de volumen constante V
se define matematicamente como,
ri =
, la velocidad de reaccion
1 dNi
V dt
(1.3.13)
1.4
dCi
dt
(1.3.14)
1.4.1
Extensi
on de la reacci
on
Con anterioridad hemos visto que existe una relacion entre las velocidades de reaccion
de los distintos componentes involucrados en una reaccion dada. Por ejemplo para la
reaccion dada por la ecuacion ?? dicha relacion esta expresada por la ecuacion 1.3.11:
rA
rB
=
a
b
masa
esto se debe a que los balances de masa se expresan como tiempo
.
un sistema de reaccion dado se considera de volumen constante si la presi
on, temperatura, y el
n
umero global de moles (sobre todo para reacciones en fase gas) no cambian durante el tiempo que
dura la reacci
on. El n
umero global de moles se refiere simplemente a la diferencia entre las moles
de producto y moles de reactivo, esta diferencia puede obtenerse f
acilmente restando los coeficientes
estequiometricos de lo productos de los coeficientes estequiometricos de los reactivos. Si la diferencia
es cero entonces se dice que no hay variacion en el n
umero global de moles del sistema.
3
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 1
Balance de Masa
(1.4.15)
(1.4.16)
donde el subndice o denota las moles iniciales de los componentes en cuestion. Nota
que la ecuacion anterior expresa un hecho importante: si tuvieramos, digamos, n componentes presentes en la mezcla de reaccion la igualdad anterior se cumpliria para
todos los n componentes. A este cociente se le llama la extension de la reaccion y se
denota por la letra ,
=
NA NAo
Nn Nno
NB NBo
= ... =
=
a
b
n
(1.4.17)
Ni Nio
i
(1.4.18)
(1.4.19)
1.4.2
i d
V dt
(1.4.20)
Conversi
on
Cuando se trabaja con reacciones qumicas resulta muy conveniente el poder tener una
idea precisa de la cantidad de materia que ha reaccionado (es decir de la masa que se
ha transformado). Es decir, podriamos definir un ndice que variara u
nicamente en un
intervalo bien definido (por ejemplo de 0 a 100). Si para la reaccion,
A B
este ndice valiera 0 significara que absolutamente nada del reactivo A ha reaccionado,
mientras que si el ndice es igual a 100 esto significara entonces que toda la masa
disponible del reactivo A se ha convertido en el producto B. El ndice que se ha
propuesto para este proposito se llama la conversi
on, y se denota por Xi donde el
subindice i denota el reactivo al que esta referida la conversion. La conversion se define
de la siguiente manera.
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 1
Balance de Masa
Fio Fi
(1.4.21)
Fio
donde Fio , y Fi representan el flujo masico o molar del componente i en la corriente de alimentacion y de producto, respectivamente.
Sistemas batch.
Xi =
Nio Ni
(1.4.22)
Nio
donde Nio y Ni representan la masa o moles del componente i al comienzo y final
del periodo de tiempo durante el que ocurre la reaccion, respectivamente.
Xi =
como puede observarse facilmente Xi solo puede tomar valores entre 0 y 1 (o bien entre
0 y 100 como porcentaje).
Notese que la conversion, tal como ha sido definida anteriormente, no es aplicable
a productos. Por ejemplo, para la reaccion antes considerada efectuada en un reactor
de flujo continuo, si la conversion estuviese referida al producto B, pero solamente el
componente A estuviera presente en la corriente de alimentacion (es decir FBo = 0)
la definicion de conversion XB estara indefinida. Un argumento similar se cumple
tambien para el reactor batch.
1.4.3
Relaci
on entre la extensi
on de la reacci
on y grado de
conversi
on
Existe una relacion estrecha entre y Xi . Dicha relacion puede obtenerse facilmente
reescribiendo las ecuaciones 1.4.18 y 1.4.22 como,
i = Nio Ni
(1.4.23)
Nio Xi = Nio Ni
(1.4.24)
c
2001,
Antonio Flores T.
Nio
UIA-Santa Fe
(1.4.25)
Captulo 2
Reactor Batch
2.1
Introducci
on
2.2
Ecuaciones de dise
no
Supongamos que cargamos un reactor batch con una cierta cantidad de un reactivo
A y que la reaccion endotermica que se efectua esta dada por: A B. El reactivo
se calienta hasta una cierta temperatura T donde comienza la reaccion anterior. El
reactor se mantiene operando a la temperatura T hasta que se logra un cierto objetivo
de operacion.
Algunas de las preguntas que nos gustara responder, en relacion a la operacion del
reactor, son:
Que tiempo debe permanecer el reactivo A en el reactor hasta que se alcance una
conversion XA ?
10
Cap. 2
Reactor Batch
11
(2.2.2)
sustituyendo:
+ CA
CAi
+ t
dCA
= k
dt
0
(2.2.3)
(2.2.4)
UIA-Santa Fe
Cap. 2
Reactor Batch
12
XA = 1 ekt
1
t = ln(1 XA )
k
1
= ln(1 .9) = 46 min
.05
de la ecuacion anterior podemos observar que a medida que la velocidad de reaccion es mayor
(k se incrementa) se podra lograr la misma conversion en menores tiempos. Aumentos en
k pueden lograrse, por ejemplo, aumentando la temperatura a la que opera el reactor. Esto
significa , en terminos practicos, que la disminucion del tiempo de reaccion podra implicar
aumentar el consumo del medio de calentamiento.
En la mayora de las situaciones practicas resulta difcil, o imposible, el disponer
de la solucion analtica de la ecuacion de dise
no. En casos como este se debe recurrir
a una solucion n
umerica o aproximada del problema en cuestion. Una situacion como
esta se muestra en el siguiente problema.
Ejemplo 2 Se dispone de un reactor batch donde se alimentan ciertas cantidades iniciales
de los reactivos A, B y C. La cinetica de la reaccion esta dada por:
rA = kCA.5 CB CC
donde k = 1x106 . El volumen inicial del reactor es 3000 lts con una concentraccion
del reactivo A de 800 mol/lt. Las moles iniciales de los reactivos B y C son 2160000 y
1920000, respectivamente.
Determinar el perfil del grado de conversion.
dXA
rA V
=
dt
NAo
donde,
rA = kCA.5 CB CC
CAo (1 XA )
CA =
1 + XA
CAo (B XA )
CB =
1 + XA
CAo (C XA )
CC =
1 + XA
CAo (D XA )
CD =
1 + XA
V = Vo (1 + XA )
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 2
Reactor Batch
13
ca=cao*(1-xa)/cden
cb=cao*(thetab-xa)/cden
cc=cao*(thetac-xa)/cden
cd=cao*(thetad+xa)/cden
cden=1+epsilon*xa
v=vo*cden
nao=cao*vo
epsilon=-2*yao
yao=1/(thetaa+thetab+thetac)
k=1e-6
cao=800
thetaa=1
thetab=.9
thetac=.8
thetad=0 vo=3000
t(0)=0
xa(0)=0
t(f)=500
En la figura 2.1 se muestra la variacion del perfil del grado de conversion con respecto
al tiempo, mientras que en la figura 2.2 se muestra la forma como las concentraciones de
los reactivos y productos cambian con el tiempo.
2.3
Reacciones m
ultiples
UIA-Santa Fe
Cap. 2
14
Reactor Batch
0.8
0.7
Grado de Conversion
0.6
0.5
0.4
0.3
0.2
0.1
0
0
50
100
150
200
250
300
Tiempo (min)
350
400
450
500
A 1 B
k
A 2 C
como hemos visto anteriormente las ecuaciones de dise
no de los tres componentes involucrados en la reaccion estan dados por:
dNA
= rA V
dt
dNB
= rB V
dt
dNC
= rC V
dt
dado que el reactivo A aparece en dos reacciones su velocidad de reaccion global (rA )
estara dada por la suma algebraica de 2 contribuciones individuales correspondientes
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 2
15
Reactor Batch
1500
Concentracion (mol/lt)
Cd
1000
Ca
500
Cb
Cc
0
0
50
100
150
200
250
300
Tiempo (min)
350
400
450
500
(2.3.5)
donde,
rA1 = k1 CA
rA2 = k2 CA
(2.3.6)
(2.3.7)
los compuestos B y C solo aparecen en una reaccion por lo que sus velocidades de
reaccion globales estaran dadas por:
rB = k1 CA
rC = k2 CA
(2.3.8)
(2.3.9)
1
2
3
A
B
C
D
las cuales pueden escribirse como un sistema de reacciones simples (de un solo paso)
como se muestra:
k
A 1 B
k1
B A
k
B 2 C
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 2
16
Reactor Batch
k2
C B
k
C 3 D
k
D 3 C
=
=
=
=
k1 CA + k1 CB
k1 CA k1 CB k2 CB + k2 CC
k2 CB k2 CC k3 CC + k3 CD
k3 CC k3 CD
(2.3.14)
(2.3.15)
(2.3.16)
(2.3.17)
notese que las ecuaciones tambien se pueden escribir en terminos de las concentraciones:
dCA
dt
dCB
dt
dCC
dt
dCD
dt
= rA
(2.3.18)
= rB
(2.3.19)
= rC
(2.3.20)
= rD
(2.3.21)
Ejemplo 3 Determinar, durante los primeros 30 minutos de operacion, el perfil de composiciones de los reactivos y productos que intervienen en el siguiente esquema de reaccion:
k
1
2
3
A
B
C
D
que ocurre en un reactor batch el cual se carga inicialmente con solo 10 mol/lt del reactivo
A. Las constantes cineticas estan dadas por:
UIA-Santa Fe
Cap. 2
Reactor Batch
(1)
(2)
(3)
(4)
=
=
=
=
ra;
rb;
rc;
rd;
fx = fx;
% -- End of file batchmr.m -c
2001,
Antonio Flores T.
UIA-Santa Fe
17
Cap. 2
18
Reactor Batch
10
C
A
C
B
CC
CD
Concentracion (mol/lt)
10
15
Tiempo (min)
20
25
30
Figura 2.3: Variacion del perfil de concentraciones en un reactor batch con reacciones
m
ultiples.
c
2001,
Antonio Flores T.
UIA-Santa Fe
Captulo 3
Reactor Continuo de Tanque
Agitado
3.1
Introducci
on
Alimentacion
Calentamiento/
enfriamiento
Producto
19
Cap. 3
3.2
20
Ecuaciones de Dise
no
FA = QCA
(3.2.2)
si suponemos que el volumen (V ) ocupado por la mezcla de reaccion permanece constante durante todo el transcurso de la reaccion,
Q
dCA
= (CAo CA ) + rA
dt
V
(3.2.3)
La ecuacion de dise
no esta dada por,
Q
dCA
= (CAo CA ) kCA
dt
V
(3.2.4)
(3.2.5)
donde,
a=
integrando,
Q
QCAo
, b= +k
V
V
+ CA
CAi
c
2001,
Antonio Flores T.
+ t
dCA
=
dt
a bCA
0
UIA-Santa Fe
(3.2.6)
Cap. 3
21
obtenemos,
,C
ln(a bCA ) ,, A
,
=t
,
b
CAi
bt
QCAo
+
Q + KV
(1 ebt )
(3.2.7)
En este caso,
Q
dCA
= (CAo CA ) kCA2
dt
V
la cual la reescribimos como,
dCA
= a bCA cCA2
dt
(3.2.9)
(3.2.10)
donde,
a=
integrando,
Q
QCAo
, b=
, c=k
V
V
+ CA
CAi
obtenemos,
2 arctanh
+ t
dCA
=
dt
a bCA cCA2
0
2 cCA b
4 ac + b2
.,CA
,
,
,
,
CAi
1
=t
4 ac + b2
(3.2.11)
(3.2.12)
t 4 ac + b2
2 cCAo + b
1
b
+ arctanh
tanh
4 ac + b2
CA =
2c
2
4 ac + b2
4 ac + b2
(3.2.13)
1
las condiciones de operacion conocidas como el estado estacionario de un proceso dado implican
que las propiedades que caracterizan al sistema no cambian con respecto al tiempo. Operacionalmente
esta condicion puede obtenerse, aproximadamente, permitiendo que el equipo opere durante un periodo
largo de tiempo. Esto equivale a tomar el limite de la ecuacion 3.2.7 cuando t .
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 3
22
CA =
Q2 + 4QkV CAo
2kV
(3.2.14)
rA = kCA CB
En este caso,
dCA
Q
= (CAo CA ) kCA CB
dt
V
si expresamos a CB como funcion de CA ,
(3.2.15)
CB = CBo (CAo CA )
(3.2.16)
(3.2.17)
donde,
a=
Q
QCAo
, b=
+ kCBo + kCAo
V
V
/
, c=k
3.3
Q
V
2 + 3
4QCAo k
+ kCBo + kCAo
2k
Q
V
+ kCBo + kCAo
22
(3.2.18)
Dise
no de un reactor tanque agitado
UIA-Santa Fe
Cap. 3
23
Evaluar
CA
V
Q
CAo
Q(CAo CA )
kCA
lt
100 min
(70 10) mol
lt
= 6000 lts
=
1
(.1 min
)(10 mol
)
lt
=
ahora supongamos que la concentracion de la corriente de alimentacion de este reactor, que opera en estado estacionario, experimenta un incremento instantaneo de tipo
escalon el cual modifica CAo hasta 72 mol/lt. Las preguntas que nos gustara responder
son:
Como cambia la concentracion a la salida del reactor con respecto al tiempo ?
Cuanto tiempo debe transcurrir para que el reactor alcanze un nuevo estado
estacionario ?
ambas preguntas se pueden responder facilmente usando la ecuacion de dise
no en forma
dinamica:
CA =
CA =
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 3
24
Solucion analitica
10.3
10.25
Concentracion
10.2
10.15
10.1
10.05
10
0
10
15
20
25
30
Tiempo (min)
35
40
45
50
UIA-Santa Fe
Cap. 3
25
10.3
10.25
Concentracion
10.2
10.15
10.1
10.05
10
0
10
15
20
25
30
Tiempo (min)
35
40
45
50
3.4
Reacciones m
ultiples
1
2
3
A
B
C
D
A 1 B
k
B 1 A
k
B 2 C
k2
C B
k
C 3 D
k3
D C
las ecuaciones de dise
no para cada componente involucrado en el sistema de reaccion
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 3
26
= FAo FA + rA V
(3.4.20)
= FBo FB + rB V
(3.4.21)
= FCo FC + rC V
(3.4.22)
= FDo FD + rD V
(3.4.23)
=
=
=
=
k1 CA + k1 CB
k1 CA k1 CB k2 CB + k2 CC
k2 CB k2 CC k3 CC + k3 CD
k3 CC k3 CD
(3.4.24)
CAo CA
+ rA
CBo CB
=
+ rB
CCo CC
=
+ rC
CDo CD
=
+ rD
(3.4.25)
(3.4.26)
(3.4.27)
(3.4.28)
1
2
3
A
C
D
B
suponiendo que la corriente de alimentacion esta solo compuesta del reactivo A con una
concentracion de 10 mol/lt y que el tiempo de residencia en el reactor es de 30 minutos. El
reactor opera en estado estacionario, emplear los mismos valores numericos de las constantes
de velocidad de reaccion que fueron usados en el ejemplo 2.3.
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 3
27
CA k1 CA + k1 CB = CAo
CB k1 CA k1 CB k2 CB + k2 CC = CBo
1
CC k2 CB k2 CC k3 CC + k3 CD = CCo
o en forma matricial,
CD k3 CC k3 CD = CDo
(1 + k1 )
k1
0
0
(1 + k1 + k2 )
k2
0
k1
(1 + k2 + k3 )
k3
0
k2
0
0
k3
(1 + k3 )
CA
CB
CC
CD
CAo
CBo
CCo
CDo
x=Ab
las concentraciones obtenidas en el CSTR al cabo de 30 minutos de tiempo de residencia
son:
CA = 2.252, CB = 1.281, CC = 1.332, CD = 5.136
a continuacion se muestra el listado Matlab usado para resolver este ejemplo.
%
% Multiple reactions in a steady-state isothermal CSTR
%
k1=.2; k1p=.15; k2=1; k2p=.8; k3=.9; k3p=.2; theta=30;
cao=10; cbo=0; cco=0; cdo=0;
%
% define elements of the A matrix
%
a = zeros(4,4);
a(1,1) = 1+k1*theta;
a(1,2) = -k1p*theta;
a(2,1) = -k1*theta;
a(2,2) = 1+k1p*theta+k2*theta;
a(2,3) = -k2p*theta;
a(3,2) = -k2*theta;
a(3,3) = 1+k2p*theta+k3*theta;
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 3
a(3,4) = -k3p*theta;
a(4,3) = -k3*theta;
a(4,4) = 1+k3p*theta;
%
% Elements of the independent vector
%
b = [cao cbo cco cdo];
%
% solve the system of AE
%
x = inv(a)*b;
fprintf (Ca = %8.3f \n,x(1))
fprintf (Cb = %8.3f \n,x(2))
fprintf (Cc = %8.3f \n,x(3))
fprintf (Cd = %8.3f \n,x(4))
% -- End of file cstrmr.m --
c
2001,
Antonio Flores T.
UIA-Santa Fe
28
Captulo 4
Reactor Tubular
4.1
Introducci
on
Generalmente este tipo de reactor es un tubo (o tal vez una familia de tubos) donde se
llevan a efecto las reacciones involucradas. En este reactor existe una corriente continua
de alimentacion, y una corriente continua de extraccion de los productos (ver figura 4.1).
El patron de flujo en el interior del reactor se supone que obedece el modelo de flujo
Calentamiento/enfriamiento
Producto
Alimentacion
Clase 4
30
Reactor Tubular
casos como estos la perdida de presion a lo largo del reactor puede ser significativa y
debera tomarse en cuenta en el modelamiento de la operacion del reactor.
En sistemas de propiedades heterogeneas, el modelamiento (balances de masa y
energa) se realiza tomando un peque
no elemento diferencial de volumen del sistema
original (ver figura 4.2).
V
Fio(y)
Fi (y+ y)
dNi
dt
(4.1.1)
(4.1.2)
(4.1.4)
entonces,
1 dNi dFi
+
= ri A
(4.1.5)
y dt
dy
expresando tanto al n
umero de Moles Ni como al flujo de masa hacia el reactor Fi en
terminos de la concentracion:
1 d(V Ci ) d(QCi )
+
= ri A
y dt
dy
c
2001,
Antonio Flores T.
UIA-Santa Fe
(4.1.6)
Clase 4
31
Reactor Tubular
(4.1.7)
(4.1.8)
(4.1.9)
FA = FAo (1 XA )
(4.1.10)
(4.1.11)
(4.1.12)
obtenemos,
A+L
1 + CA dCA
=
dL
k CAo CA
Q 0
kA
CA = CAo e Q L
(4.1.13)
(4.1.14)
rA = kCA2
tenemos,
A+L
1 + CA dCA
=
dL
k CAo CA2
Q 0
CA =
c
2001,
Antonio Flores T.
QCAo
Q + kACAo L
UIA-Santa Fe
(4.1.15)
(4.1.16)
Clase 4
32
Reactor Tubular
rA = kCA CB
=
dL
k CAo CA2 + (CBo CAo )CA
Q 0
obtenemos,
donde,
u=
CBo CAo
eu
CBo CAo eu
0
(4.1.18)
kAL
(CBo CAo )
Q
(4.1.19)
CA = CAo
(4.1.17)
Ejemplo 5 Determine el volumen de un reactor tubular necesario para producir 300 millones de libras de etileno/ano a partir del crakeo de una corriente de etano puro. La reaccion
es de primer orden, irreversible y elemental,
C2 H6 C2 H4 + H2
dXA
= rA
dV
+ XA
0
+ XA
dXA
dXA
= FAo
rA
kCA
0
1 XA
1 + XA
FAo + XA (1 + XA )
V =
dXA
(1 XA )
kCAo 0
1
FAo
(1 + )ln
XA
V =
kCAo
1 XA
:
esta ecuacion permite entonces evaluar el volumen V que se requiere para lograr la conversion deseada XA suponiendo que la operacion sea a temperatura, y presion constantes.
A continuacion evaluamos los parametros necesarios.
c
2001,
Antonio Flores T.
UIA-Santa Fe
Clase 4
33
Reactor Tubular
FAo .
dado que la corriente de alimentacion es etano puro FB = FAo XA entonces,
FAo
300x106
FB
=
=
XA
.8
lb
ano
= .425
lbmol
seg
CAo
si la corriente de alimentacion es etano puro yAo = 1 entonces,
6 atm
P yAo
lbmol
=
= .00415
CAo =
3 atm
pie
RT
pie3
(.73 molo R )(1980 o R)
k a 1100 o K.
dado que tenemos el valor de la constante de velocidad de reaccion k disponible a
1000 o K, pero la reaccion ocurre a 1100 o K, debemos usar alguna forma de corregir
el valor de k para poder utilizarla a esta ultima temperatura. Esta correccion por
efectos de temperatura se puede realizar facilmente definiendo T1 = 1000 o K y
T2 = 1100 o K. Escribiendo las ecuacion de Arrhenius a estas dos temperaturas,
E
RT
E
RT
k(T1 ) = Ae
k(T2 ) = Ae
Ae RT2
k(T2 )
=
E
k(T1 )
Ae RT1
de esta ecuacion,
k(T2 ) = k(T1 )e
E
R
entonces,
k(T2 ) = .072
<
1
e
seg
cal
gmol
cal
1.987
gmolo K
82000
1
T1
T1
2
1
1
1100
[ 1000
] o1K
= 3.06
1
seg
lbmol
)
seg
1
V =
(1 + 1)ln
.8 = 80.9 pie3
1
lbmol
1 .8
(3.06 seg )(.00415 pie3 )
:
Normalmente los reactores tubulares no estan formados por un unico tubo sino por un grupo
de tubos de tamano mas pequeno, cuyo volumen total es igual al volumen del reactor cuando
se utiliza un solo tubo. Por ejemplo si decidimos emplear tubos cuya longitud individual
es 40 pies, de 2 de diametro cedula 40 (cuya area de seccion transversal es .0205 pie2 ),
entonces el numero requerido de tubos individuales n (o sea de reactores individuales) estara
dado por,
n=
80.9 pie3
= 98.65 tubos 100 tubos
(40 pies)(.0205 pie2 )
c
2001,
Antonio Flores T.
UIA-Santa Fe
Clase 4
4.2
Reactor Tubular
34
Reacciones m
ultiples
1
2
3
A
B
C
D
A 1 B
k1
B A
k
B 2 C
k
C 2 B
k
C 3 D
D 3 C
= rA
(4.2.20)
= rB
(4.2.21)
= rC
(4.2.22)
= rD
(4.2.23)
=
=
=
=
k1 CA + k1 CB
k1 CA k1 CB k2 CB + k2 CC
k2 CB k2 CC k3 CC + k3 CD
k3 CC k3 CD
UIA-Santa Fe
(4.2.24)
Clase 4
Reactor Tubular
35
= rA
(4.2.25)
= rB
(4.2.26)
= rC
(4.2.27)
= rD
(4.2.28)
c
2001,
Antonio Flores T.
UIA-Santa Fe
Captulo 5
Reactor SemiBatch
5.1
Introducci
on
5.2
Reactor Semibatch
Cap. 5
37
Reactor SemiBatch
d(CA V )
dNA
=
dt
dt
(5.2.1)
notese que, como se dijo antes, el volumen que ocupa la mezcla de reaccion V no
es constante. Por esta razon el volumen permanece dentro de la derivada, con
respecto al tiempo, en la anterior ecuacion. La ecuacion anterior debe entonces
diferenciarse y reescribirse como,
FAo + rA V = V
dV
dCA
+ CA
dt
dt
(5.2.2)
si de esta ecuacion deseamos evaluar, por ejemplo, la forma como la concentracion CA vara con respecto al tiempo t, necesitamos entonces encontrar una
expresion matematica que nos diga la forma como el volumen de la mezcla de
reaccion V vara con respecto al tiempo t. En otras palabras necesitamos formular
una ecuacion de donde podamos evaluar la derivada dV
. Esta expresion puede
dt
determinarse facilmente realizando el balance de masa global sobre el reactor
semibatch,
dM
FT o =
(5.2.3)
dt
donde FT o es el flujo total alimentado al reactor, y M es el holdup o masa de
lquido en el interior del reactor. Si en la ecuacion anterior substituimos M por
la siguiente equivalencia,
M = V
(5.2.4)
obtenemos,
FT o =
d(V )
dt
(5.2.5)
y si la densidad es constante,
dV
FT o
=
= Qo
dt
dCA
= Qo (CAo CA ) + rA V
dt
(5.2.6)
dCA
dt
obtenemos,
(5.2.7)
c
2001,
Antonio Flores T.
dV = Qo
UIA-Santa Fe
+ t
0
dt
(5.2.8)
Cap. 5
38
Reactor SemiBatch
(5.2.10)
NB
V
(5.2.11)
las moles del componente B se evaluan planteando el balance de masa para este
reactivo,
dNB
= rB V
(5.2.12)
dt
debido a la estequiometra de la reaccion,
por lo cual,
rB = r A
(5.2.13)
dNB
= rA V
dt
(5.2.14)
el termino
>
dNB =
(5.2.15)
rA V dt
integrando,
FAo
tenemos,
+ t
0
dt +
FAo t +
despejando,
dNA
dt
rA V dt =
+ NA
NAo
dNA
(5.2.16)
rA V dt = NA NAo
(5.2.17)
rA V dt = NA NAo FAo t
(5.2.18)
sustituyendo esta u
ltima ecuacion en la ecuacion 5.2.15 tenemos,
+ NB
NBo
c
2001,
Antonio Flores T.
(5.2.19)
Cap. 5
39
Reactor SemiBatch
entonces,
(5.2.20)
(5.2.21)
en esta ecuacion las moles del reactivo A (denotadas por NA ) se calculan simplememente de la definicion de concentracion para el reactivo A,
NA = CA V
(5.2.22)
NC
V
(5.2.23)
las moles del componente C se evaluan planteando el balance de masa para este
producto,
dNC
(5.2.24)
= rC V
dt
debido a la estequiometra de la reaccion,
por lo cual,
rA = rC
(5.2.25)
dNC
= rA V
dt
(5.2.26)
el termino
A,
>
dNC =
rA V dt
(5.2.27)
rA V dt = NA NAo FAo t
sustituyendo esta u
ltima ecuacion en la ecuacion 5.2.27 y depejando tenemos,
NC = NCo NA + NAo + FAo t
(5.2.28)
UIA-Santa Fe
Cap. 5
40
Reactor SemiBatch
En este modo de operacion el reactor se carga con una cierta cantidad inicial del
reactivo A mientras que el reactivo B se agrega al reactor como se muestra en
la figura 5.2. Suponiendo que A es el reactivo lmite, balance de masa para el
componente A,
dNA
(5.2.30)
rA V =
dt
del balance global de masa,
V = Vo + Qo t
(5.2.31)
dXA
= rA (Vo + Qo t)
dt
(5.2.32)
NAo (1 XA )
NA
=
V
Vo + Qo t
(5.2.33)
NB
V
(5.2.34)
y,
CB =
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 5
41
Reactor SemiBatch
el n
umero de moles del componente B que existen en cualquier instante NB se
puede evaluar del balance de masa para este componente,
FBo + rB V =
dNB
dt
dNB
dt
+ t
0
dt +
entonces,
FBo t +
el termino
ecuacion,
>t
0
+ t
0
+ t
0
rA V dt =
+ NB
NBo
dNB
rA V dt = NB NBo
(5.2.35)
rA V dt =
+ NA
NAo
dNA
rA V dt = NAo
+ XA
0
dXA = NAo XA
(5.2.36)
sustituyendo ahora las expresiones para CA y CB (dadas por las ecuaciones 5.2.33
y 5.2.36) en la ecuacion 5.2.32,
NAo
c
2001,
Antonio Flores T.
UIA-Santa Fe
(5.2.37)
Cap. 5
Reactor SemiBatch
42
Ejemplo 6 La reaccion elemental A + B C se realiza en un reactor semibatch. Inicialmente se cargan 60 lts del reactivo A con una concentracion de 55 mol/lt. El reactor se
alimenta con un flujo volumetrico continuo del reactivo B equivalente a 50 lt/min con una
concentracion de 30 mol/lt. La constante de velocidad de reaccion es 3.5x103 lt/(molmin). Determinar el tiempo necesario para convertir el 100 % del reactivo A. Graficar
tambien los perfiles de concentracion de reactivos y productos.
Previamente hemos obtenido la ecuacion de diseno para esta situacion,
dXA
k(1 XA )(NBo + FBo t NAo XA )
=
dt
Vo + Qo t
las concentracciones estaran dadas por,
NA
NAo (1 XA )
=
V
Vo + Qo t
NBo + FBo t NAo XA
=
Vo + Qo t
NAo XA
=
Vo + Qo t
CA =
CB
CC
el siguiente listado muestra el programa Polymath empleado para resolver este problema.
d(xa)/d(t)=k*(1-xa)*(nbo+fbo*t-nao*xa)/cden
k=3.5e-3
nbo=0
vo=60
cao=55
cbo=30
q=50
fbo=cbo*q
nao=vo*cao
cden=vo+q*t
ca=nao*(1-xa)/cden
cc=nao*xa/cden
cb=(nbo+fbo*t-nao*xa)/cden
t(0)=0
xa(0)=0
t(f)=50
Como se muestra en la figura 5.3 se requieren aproximadamente 50 minutos para consumir
el reactivo A en su totalidad. En la figura 5.4 se muestran la variacion de los perfiles de
concentracion.
En la mayora de las ocasiones las ecuaciones de dise
no que resultan del modelamiento
de reactores semibatch son complicadas de integrar analiticamente. Por esta razon
la aproximacion numerica de la solucion es un camino viable para resolver problemas
relativos a diferentes aspectos de este tipo de reactores (tales como modelamiento,
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 5
43
Reactor SemiBatch
1
0.9
0.8
Grado de Conversion
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0
10
15
20
25
30
Tiempo (min)
35
40
45
50
simulacion, control etc). Para ilustrar el empleo del procedimiento numerico de solucion
considere la situacion mostrada en la figura 5.5. Se desean analizar los siguientes dos
modos de operacion de un reactor semibatch con objeto de decidir el tipo de operacion
que produce la mayor cantidad de producto. Las ecuaciones de dise
no para cada
situacion son las siguientes:
AB
dNA
= FAo + rA V
dt
dNB
= rB V
dt
rA = kCA CB
(5.2.38)
(5.2.39)
BA
dNA
= rA V
dt
dNB
= FBo + rB V
dt
rB = rA
(5.2.40)
(5.2.41)
en ambos casos la variacion del volumen con respecto al tiempo esta representado por
la siguiente ecuacion diferencial:
dV
=q
dt
c
2001,
Antonio Flores T.
UIA-Santa Fe
(5.2.42)
Cap. 5
44
Reactor SemiBatch
60
Concentracion (mol/lt)
50
40
Ca
30
Cb
20
10
Cc
0
0
10
15
20
25
30
Tiempo (min)
35
40
45
50
(a)
(b)
Figura 5.5:
Ejemplo 7 Usando los datos del ejemplo 6 calcular los perfiles de concentracion de productos y reactivos suponiendo que el reactor semibatch puede ser operado de las 2 formas
anteriormente discutidas.
En las figuras 5.6(a) y 5.6(b) se muestran los perfiles de concentraciones como resultado de
los modos diferentes de operar al reactor batch. A continuacion se muestran los programas
Matlab usados para generar estos resultados.
clear all; clc
%
% Dynamic simulation of a semi-batch reactor. The program allows for
% 2 simulation modes:
% (1) A is fed to a reactor containing reactant B (wcase=1)
% (2) B is fed to a reactor containing reactant A (wcase.ne.1)
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 5
45
Reactor SemiBatch
60
60
50
50
CA
CB
Cc
30
30
20
20
10
10
40
Concentracion (mol/lt)
Concentracion (mol/lt)
40
C
A
C
B
C
10
15
Tiempo (min
20
25
30
(a)
10
15
Tiempo (min
(b)
Figura 5.6:
%
global cao cbo q vo k
wcase = input( Which case (A->B = 1, B->A =/ 1) ?);
cao=55; cbo=30; q=50; vo=60; k=3.5e-03;
time = linspace(0,30);
if wcase == 1
x0 = [0 cbo*vo 0 vo]; % initial conditions
else
x0 = [cao*vo 0 0 vo];
end
[t,x] = ode15s(caseab,time,x0,[],wcase)
na
nb
nc
v
ca
=
=
=
=
=
plot(time,ca,red,time,cb,green,time,cc,blue)
xlabel(Tiempo (min)
legend(C_{A},C_{B},C_{c})
UIA-Santa Fe
20
25
30
Cap. 5
Reactor SemiBatch
46
ca = na/v;
cb = nb/v;
ra
rb
rc
= -k*ca*cb;
= ra;
= -ra;
if wcase == 1
fao = q*cao;
fx (1) = fao+ra*v;
fx (2) = rb*v;
else
fbo = q*cbo;
fx (1) = ra*v;
fx (2) = fbo+rb*v;
end
fx (3) = rc*v;
fx (4) = q;
fx = fx;
% -- End of file caseab.m --
5.3
Reacciones m
ultiples
1
2
3
A
B
C
D
A 1 B
k1
B A
k
B 2 C
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 5
Reactor SemiBatch
47
k2
C B
k
C 3 D
k
D 3 C
= FAo + rA V
(5.3.43)
= rB V
(5.3.44)
= rC V
(5.3.45)
= rD V
(5.3.46)
= q
(5.3.47)
=
=
=
=
k1 CA + k1 CB
k1 CA k1 CB k2 CB + k2 CC
k2 CB k2 CC k3 CC + k3 CD
k3 CC k3 CD
(5.3.48)
FAo
CA q
+ rA
V
V
CB q
= rB
V
CC q
= rC
V
CD q
= rD
V
=
UIA-Santa Fe
(5.3.49)
(5.3.50)
(5.3.51)
(5.3.52)
Cap. 5
48
Reactor SemiBatch
Ejemplo 8 Determinar los perfiles de concentracion de todos los componentes involucrados en el esquema de reaccion dado por:
k
1
2
3
A
B
C
D
CA
CB
CC
CD
3.5
Concentracion (mol/lt)
2.5
1.5
0.5
10
15
Tiempo (min)
20
25
30
UIA-Santa Fe
Cap. 5
Reactor SemiBatch
49
[t,x] = ode15s(sbatchmr,time,x0);
ca=x(:,1); cb=x(:,2); cc=x(:,3); cd=x(:,4);
plot(t,ca,red,t,cb,green,t,cc,blue,t,cd,magenta);
xlabel(Tiempo (min)), ylabel(Concentracion (mol/lt))
legend(C_{A},C_{B},C_{C},C_{D})
% -- End of file runsbatchmr.m --
(1)
(2)
(3)
(4)
(5)
=
=
=
=
=
fao/v+ra-ca*q/v;
rb-cb*q/v;
rc-cc*q/v;
rd-cd*q/v;
q;
fx = fx;
% -- End of file sbatchmr.m --
5.4
Operaci
on Semibatch-Batch
UIA-Santa Fe
Cap. 5
Reactor SemiBatch
50
Q (lt/min)
19
0
19
0
= rA V
= FBo + rB V
= rC V
= q
Batch
dNA
dt
dNB
dt
dNC
dt
dV
dt
c
2001,
Antonio Flores T.
= rA V
= rB V
= rC V
= 0
UIA-Santa Fe
Cap. 5
51
Reactor SemiBatch
en ambos casos,
rA = kCA CB
rB = rA
rC = rA
Los perfiles de variacion del numero de moles obtenidos durante los primeros 400 minutos
de operacion se muestran en la figura 5.8. El programa Matlab usado para efectuar los
calculos se muestra a continuacion.
4
x 10
2.5
Moles
N
a
N
b
Nc
1.5
0.5
50
100
150
200
Tiempo (min)
250
300
350
400
UIA-Santa Fe
Cap. 5
Reactor SemiBatch
= x(4);
if time >0 & time < 10 | time > 130 & time < 140
q = qf; fbo=q*cbo;
else
q = 0; fbo = 0;
end
f
f
f
f
(1)
(2)
(3)
(4)
=
=
=
=
ra*v;
fbo+rb*v;
rc*v;
q;
f = f;
%-- End of the estcont.m file --
c
2001,
Antonio Flores T.
UIA-Santa Fe
52
Captulo 6
Relaciones Estequiom
etricas
6.1
(6.1.2)
moles convertidas de A
= [moles iniciales de A]
moles iniciales de A
= [moles convertidas de A]
de manera semejante el n
umero de moles del reactivo B en cualquier instante se puede
obtener empleando una ecuacion analoga,
b
(6.1.3)
NB = NBo NAo XA
a
que el segundo termino del lado derecho de esta ecuacion representa el n
umero de
moles del reactivo B en cualquier instante puede verificarse facilmente analizando las
unidades de la ecuacion,
?
b
moles convertidas de B
NAo XA =
[moles convertidas de A]
a
moles convertidas de A
= [moles convertidas de B]
de manera semejante se puede obtener una ecuacion similar a la ecuacion 6.1.3 para
cualquier otro reactivo o producto. Por ejemplo para el producto C la ecuacion correspondiente del n
umero de moles en cualquier instante de tiempo estara dada por,
c
NC = NCo + NAo XA
(6.1.4)
a
53
Cap. 6
Relaciones Estequiom
etricas
54
los resultados para la reaccion elemental dada por la ecuacion 6.1.1 se resumen en la
siguiente tabla.
Componente
A
B
C
D
Inertes
Moles
iniciales
NAo
NBo
NCo
NDo
NIo
Moles consumidas/
Moles restantes
formadas al tiempo t
al tiempo t
NAo XA
NA = NAo NAo XA
b
a NAo XA
NB = NBo ab NAo XA
c
N X
NC = NCo + ac NAo XA
a Ao A
d
N X
ND = NDo + ad NAo XA
a Ao A
NIo
donde el n
umero total de moles NT esta dado por,
donde,
NT = NT o + NAo XA
(6.1.5)
(6.1.6)
y,
b
d c
+ 1
(6.1.7)
a a a
usando la informacion contenida en esta tabla podemos expresar la concentracion de
cualquier componente en funcion de la conversion del reactivo lmite (y de algunos
otros terminos).
Por ejemplo para el componente A por definicion sabemos que,
=
CA =
NA
V
NAo (1 XA )
V
(6.1.8)
(6.1.9)
NCo + ac NAo XA
NC
=
CC =
V
V
(6.1.10)
NDo + ad NAo XA
ND
=
V
V
(6.1.11)
Nio
NAo
(6.1.12)
CB =
para el producto C,
para el producto D,
CD =
si definimos el siguiente termino,
i =
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 6
55
Relaciones Estequiom
etricas
NAo
NAo
NAo
NBo
NAo
NCo
NAo
NDo
NAo
ab XA
+ ac XA
V
+ ad XA
=
=
(6.1.13)
(6.1.14)
NAo B ab XA
(
NAo C + ac XA
V
NAo D + ad XA
=
(6.1.15)
V
V
para completar este desarrollo debemos encontrar la forma como el volumen V vara
con respecto a la conversion XA . Esta variacion ocurre principalmente para la gran
mayora de reacciones en fase gas, y para unas pocas en fase lquida (por ejemplo en
reacciones de polimerizacion). Sin embargo, en algunos casos el volumen V ocupado
por la mezcla de reaccion permanece aproximadamente constante durante el curso de
la reaccion lo cual permite simplificar los calculos como se muestra a continuacion.
En algunos casos el volumen ocupado por la mezcla de reaccion vara durante el
transcurso de la reaccion. En reactores batch esto puede ocurrir principalmente para
reacciones en fase gas. Un caso concreto donde el volumen cambia ocurre en reacciones
donde el n
umero total de moles de los reactivos es diferente del n
umero total de moles
de los productos. Por ejemplo considere la reaccion,
A + 3B 2C
(6.1.16)
(6.1.17)
la cual se supondra valida en cualquier punto del reactor batch, y para cualquier
instante de tiempo t. En la ecuacion anterior P es la presion, Z es el factor de compresibilidad, NT es el n
umero total de moles, R es la constante de los gases ideales, y
T es la temperatura.
Es claro que durante el curso de la reaccion estos valores de P, T y NT pueden ir
cambiando. Por ejemplo al inicio de la reaccion las condiciones de operacion estaran
dadas por ciertos valores iniciales de presion, temperatura, y n
umero de moles denotados por Po , To , y NT o , respectivamente. Entonces en el instante t = 0 (inicio de la
reaccion) la ecuacion de estado, dada por la ecuacion 6.1.17, se escribira como,
Po Vo = Zo NT o RTo
(6.1.18)
UIA-Santa Fe
(6.1.19)
Cap. 6
56
Relaciones Estequiom
etricas
Po
P
0/
T
To
0/
Z
Zo
0/
NT
NT o
(6.1.20)
de la tabla estequiometrica 1 ,
NT = NT o + NAo XA
substituyendo esta ecuacion en la ecuacion 6.1.20,
V = Vo
Po
P
T
To
0/
0/
T
To
0/
Z
Zo
0-
NT o + NAo XA
NT o
(6.1.21)
NAo
.
NT o
Po
P
0/
Z
(1 + yAo XA )
Zo
0
(6.1.22)
= yAo
(6.1.23)
?
@?
@
cambio en el numero total de moles
moles de A
=
(6.1.24)
moles de A
numero de moles iniciales
?
@
cambio en el numero total de moles
(6.1.25)
=
numero de moles iniciales
entonces la ecuacion 6.1.22 se reescribe como,
V = Vo
Po
P
0/
T
(1 + XA )
To
0
(6.1.26)
NAo (1 XA )
CAo (1 XA ) P
NA
= 1P 2 1 T 2
=
V
(1 + XA )
Po
Vo Po To (1 + XA )
/
0/
To
T
(6.1.27)
en la siguiente tabla se resumen las ecuaciones para todos los componentes involucrados
en la reaccion en cuestion.
CA
NA
V
CB
NB
V
CC
NC
V
CD
ND
V
=
=
=
NAo (1XA )
)( TT )(1+XA )
Vo (
Po
P
Vo (
Po
P
)( TTo )(1+XA )
NAo [C + ac XA ]
Vo ( PPo )( TT
)(1+XA )
NAo [D + ad XA ]
Vo (
Po
P
CAo (1XA )
(1+XA )
CAo [B ab XA ]
(1+XA )
CAo [C + ac XA ]
(1+XA )
CAo [D + ad XA ]
(1+XA )
NAo [B ab XA ]
)( )(1+XA )
T
To
1
c
aC
P
Po
21
P
Po
P
Po
P
Po
To
T
21
21
21
To
T
To
T
To
T
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 6
57
Relaciones Estequiom
etricas
0.9
=0
0.8
=3
0.7
0.6
0.5
0.4
0.3
0.2
0.1
10
20
30
40
50
60
Tiempo (min)
70
80
90
100
UIA-Santa Fe
Cap. 6
58
Relaciones Estequiom
etricas
0.9
=0
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
10
20
30
40
50
60
Tiempo (min)
70
80
90
100
UIA-Santa Fe
Cap. 6
6.2
59
Relaciones Estequiom
etricas
@?
de manera semejante puede obtenerse el flujo para cualquier otro componente involucrado en la reaccion en terminos de la conversion XA . Este paso es totalmente semejante
a lo que se hizo para el reactor batch.
Para el reactivo B,
b
FB = FBo FAo XA
a
que el segundo termino del lado derecho de esta ecuacion representa el n
umero de
moles del reactivo B en cualquier instante puede verificarse facilmente analizando las
unidades de la ecuacion,
b
FAo XA =
a
@?
expresiones semejantes pueden obtenerse para cualquier otro componente. Por ejemplo
para el producto C el n
umero de moles en un punto del reactor donde la conversion es
XA esta dada por,
c
(6.2.28)
FC = FCo + FAo XA
a
los resultados para la reaccion simple en cuestion se resumen en la siguiente tabla.
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 6
Relaciones Estequiom
etricas
Componente
A
B
C
D
Inertes
60
Moles a
Moles consumidas/
Moles a
la entrada
formadas
la salida
FAo
FAo XA
FA = FAo FAo XA
b
FBo
a FAo XA
FB = FBo ab FAo XA
c
FCo
F X
FC = FCo + ac FAo XA
a Ao A
d
FDo
F X
FD = FDo + ad FAo XA
a Ao A
FIo
FIo
donde,
FT = FT o + FAo XA
(6.2.29)
(6.2.30)
y,
d c
b
(6.2.31)
+ 1
a a a
usando la informacion contenida en esta tabla podemos expresar la concentracion
de cualquier componente involucrado en la reaccion Ci en funcion de la conversion de
A.
Para el reactivo A por definicion la concentracion de la especie A en un cierto punto
del reactor se determina de la ecuacion,
=
CA =
donde FA es el flujo molar de la especie A en
(
(
volumen total
.
tiempo
)
moles de A
.
volumen total
FA
Q
(
moles de A
tiempo
(6.2.32)
)
, y Q es el flujo volumetrico en
para el producto C,
para el producto D,
FBo ab FAo XA
FB
=
CB =
Q
Q
(6.2.34)
FCo + ac FAo XA
FC
CC =
=
Q
Q
(6.2.35)
FDo + ad FAo XA
FD
CD =
=
Q
Q
(6.2.36)
Este comentario solo es aplicable al reactor tubular ya que en el reactor tanque agitado FA y Q
pueden medirse mas facilmente a la entrada del reactor
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 6
61
Relaciones Estequiom
etricas
Fio
FAo
el anterior grupo de ecuaciones se puede reescribir de la siguiente manera,
i =
CB =
CC =
FAo
FAo
FAo
FBo
FAo
FCo
FAo
FDo
FAo
ab XA
+ ac XA
Q
+ ad XA
=
=
(6.2.38)
(6.2.39)
FAo B ab XA
Q
FAo C + ac XA
Q
(6.2.37)
FAo D + ad XA
=
(6.2.40)
Q
Q
a continuacion debemos encontrar la forma como el flujo volumetrico Q vara con
respecto a la conversion XA . De manera semejante al caso del reactor batch, en caso
de que el volumen ocupado por la mezcla de reaccion, en un reactor de flujo continuo,
cambie a lo largo del reactor, el flujo volumetrico Q debera expresarse como una funcion
del grado de conversion XA .
Para reactores de flujo continuo recuerdese de la ecuacion 6.2.32 la definicion de
concentracion, y si suponemos que la ecuacion de estado dada por la ecuacion 6.1.17
sigue siendo valida en cualquier punto del reactor, entonces,
CD =
CT =
P
FT
NT
=
=
V
ZRT
Q
(6.2.41)
en particular esta ecuacion es valida a las condiciones que prevalecen a la entrada del
reactor,
NT o
FT o
Po
CT o =
(6.2.42)
=
=
V
Zo RTo
Qo
donde el subindice o denota que las propiedades en cuestion se miden a la entrada
del reactor. Dividiendo la ecuacion 6.2.41 entre la ecuacion 6.2.42,
FT
Q
FT o
Qo
P
ZRT
1
2
Po
Zo RTo
(6.2.43)
FT
FT o
0/
Po
P
0/
T
To
(6.2.44)
FT o + FAo XA
FT o
./
UIA-Santa Fe
Po
P
0/
T
To
(6.2.45)
Cap. 6
62
Relaciones Estequiom
etricas
FAo
.
FT o
Po
P
0/
T
To
(6.2.46)
Finalmente,
Q = Qo (1 + XA )
Po
P
0/
T
To
(6.2.47)
donde esta definida en exactamente la misma forma como fue usada en la ecuacion
6.1.23.
La dependencia del flujo volumetrico Q con respecto a la conversion XA , dada
por la ecuacion anterior, puede ahora substituirse en la ecuacion de definicion de la
concentracion. Por ejemplo para el componente A,
CA =
CAo (1 XA ) P
FA
FAo (1 XA )
1 21 2
=
=
Q
(1 + XA )
Po
Qo PPo TTo (1 + XA )
/
0/
To
T
(6.2.48)
en la siguiente tabla se resumen las ecuaciones para todos los componentes involucrados
en la reaccion en cuestion.
CA
FA
Q
CB
FB
Q
CC
CD
=
=
FC
Q
FD
Q
FAo (1XA )
Qo ( PPo )( TT )(1+XA )
CAo (1XA )
(1+XA )
CAo [B ab XA ]
(1+XA )
CAo [C + ac XA ]
(1+XA )
CAo [D + ad XA ]
(1+XA )
=
=
FAo [B ab XA ]
Qo (
Po
P
)( )(1+XA )
Qo (
Po
P
)( )(1+XA )
T
To
FAo [C + ac XA ]
T
To
FAo [D + ad XA ]
Qo ( PPo )( TT
)(1+XA )
P
Po
21
P
Po
P
Po
P
Po
To
T
21
21
21
To
T
To
T
To
T
UIA-Santa Fe
Cap. 6
63
Relaciones Estequiom
etricas
1
0.9
=0
0.8
=1
0.7
XA
0.6
0.5
0.4
0.3
0.2
0.1
10
15
20
25
Longitud (m)
30
35
40
45
50
c
2001,
Antonio Flores T.
UIA-Santa Fe
Captulo 7
Configuraci
on de Reactores
7.1
Introducci
on
productos
Tubular-Tubular.
Tanque agitado-Tubular.
Tubular-Tanque agitado.
64
Cap. 7
65
Configuraci
on de Reactores
reactivos
productos
productos
Tanque agitado-Tubular.
Ejemplo 12 Se propone realizar la reaccion elemental A + B C + D en un sistema
de reactores conectados en serie consistente de un reactor tanque agitado y de un reactor tubular. Determinar con cual de los dos posibles arreglos, tanque agitado-tubular o
tubular-tanque agitado, se obtendra mayor conversion de ambos reactivos. La corriente
de alimentacion al primer reactor de la serie esta formada por una mezcla de los reactivos
A y B con concentraciones de 100 y 170 mol/lt, respectivamente. El flujo volumetrico
alimentado al primer reactor es 50 lt/min. El volumen del reactor tanque agitado es 1000
lts. La longitud del reactor tubular es 10 mts, y el radio es 30 cmts. La constante de
velocidad de reaccion es 2.5x104 lt/(mol-min).
Las ecuaciones de diseno para ambos reactores estan dadas por:
Reactor tanque agitado.
f (xA ) = FAo xA + rA V
donde,
rA = kCA CB
2
= kCAo
(1 xA )(B xA )
CBo
B =
CAo
Reactor tubular.
rA A
dxA
=
dL
FAo
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 7
66
Configuraci
on de Reactores
reactivos
productos
Figura 7.4: Arreglo tubular-tanque agitado.
=
=
=
=
fsolve(cstrseq,xg,[],[],caf,cbf);
100*(cbf-cb1)/cbf; cb10=cb1;
ode15s(tubseq,[0 lenghtub],0,[],ca1,cb1);
100*(cb10-cb1)/cb10;
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 7
Configuraci
on de Reactores
ux = length(xa2);
else
[lenght,xa1] = ode15s(tubseq,[0 lenghtub],0,[],caf,cbf);
xb1
= 100*(cbf-cb1)/cbf; cb10=cb1;
xa2
= fsolve(cstrseq,xg,[],[],ca1,cb1);
xb2
= 100*(cb10-cb1)/cb10;
ux = length(xa1);
end
xag = 100*(caf-ca1)/caf;
xbg = 100*(cbf-cb1)/cbf;
disp( ); disp( );
disp ( ----------------------------);
disp ( Reactor
Xa
Xb);
disp ( ----------------------------);
if kinds == 1
fprintf( CSTR
%6.2f %6.2f \n,100*xa1,xb1);
fprintf( Tubular
%6.2f %6.2f \n,100*xa2(ux),xb2);
else
fprintf( Tubular
%6.2f %6.2f \n,100*xa1(ux),xb1);
fprintf( CSTR
%6.2f %6.2f \n,100*xa2,xb2);
end
disp ( ----------------------------);
disp( );
fprintf ( Reactant A Global Conversion = %6.2f \n,xag);
fprintf ( Reactant B Global Conversion = %6.2f \n,xbg);
%-- End of the runseq.m file -function fx = cstrseq (x,caf,cbf)
global qf k vcstr areatub
global ca1 cb1
xa = x(1);
fao = qf*caf;
thetab = cbf/caf;
ca1 = caf*(1-xa);
cb1 = caf*(thetab-xa);
ra = -k*ca1*cb1;
fx = fao*xa+ra*vcstr;
%-- End of the cstrseq.m file -c
2001,
Antonio Flores T.
UIA-Santa Fe
67
Cap. 7
68
Configuraci
on de Reactores
7.2
X
A1
1
V
1
X
A2
V
2
Figura 7.5:
o visceversa ? para que el volumen total ocupado por los dos reactores (V1 + V2 ) sea
el menor posible. Notese que esta decision implica fundamentalmente seleccionar el
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 7
69
Configuraci
on de Reactores
C
Ao
FAo
X
A1
1
V
1
X
A2
V
2
Figura 7.6:
(7.2.1)
(7.2.2)
XA1
1
=
rA1
CAo
(7.2.3)
o bien como:
(7.2.4)
(7.2.5)
(7.2.6)
entonces,
UIA-Santa Fe
Cap. 7
Configuraci
on de Reactores
70
(7.2.7)
1
CAo
1
r
A
X
A1
Figura 7.7:
K
-
1
r
A
Area=
N
Area=
1
CAo
X
A1
2
CAo
X
A2
Figura 7.8:
UIA-Santa Fe
Cap. 7
Configuraci
on de Reactores
K
-
1
r
A
71
M
Area=
X
A1
Area= 1
CAo
2
CAo
X
A2
Figura 7.9:
en cualquier caso debe notarse que el volumen total del sistema de reaccion
(obtenido como la suma de los volumenes V1 y V2 ) sera mnimo cuando el area
definida por el rectangulo KLM N sea maxima.
Existe un procedimiento grafico simple para obtene las coordenas (-1/rA1 ,XA1 )
de tal forma que el area del rectangulo KLM N sea maxima. A continuacion
detallamos este procedimiento. Supongamos que deseamos calcular las coordenas
(y, x) de tal manera que el area ocupada por un rectangulo debajo de la curva
mostrada sea maxima:
M(x,y)
y
Area=xy
Figura 7.10:
(7.2.8)
UIA-Santa Fe
(7.2.9)
Cap. 7
72
Configuraci
on de Reactores
de donde,
y
dy
=
dx
x
(7.2.10)
este resultado significa que el area ocupada por el rectangulo sera maxima cuando
el punto M esta ubicado en un lugar (x, y) donde la pendiente de la curva es igual
a la pendiente de la diagonal N L del rectangulo como se muestra en la siguiente
grafica:
y
M(x,y)
pendiente= -
dy
dx
Diagonal=y/x
Figura 7.11:
Diagonal del
reactangulo
L
Pendiente de la
curva en M
-1
rA
N
X A1
X A2
Figura 7.12:
c
2001,
Antonio Flores T.
UIA-Santa Fe
Captulo 8
Sensibilidad param
etrica
8.1
Introducci
on
8.2
f
f
S+
x
p
(8.2.2)
donde,
S=
ademas,
Sij =
xi
pj
(8.2.3)
(8.2.4)
para poder comparar entre si las magnitudes de los coeficientes de sensibilidad frecuentemente se acostumbra usar alg
un tipo de escalamiento para este proposito. Una
de las posibles formas de escalamiento consiste en multiplicar Sij por la magnitud de
la variable independiente. En terminos de este tipo de escalamiento:
Sij = pj
xi
pj
73
xi
ln(pj )
(8.2.5)
Clase 8
8.3
74
Sensibilidad Param
etrica
Evaluaci
on de sensibilidad param
etrica
Uno de los procedimientos que pueden usarse para evaluar los coeficientes de sensibilidad parametrica consiste en integrar simultaneamente las ecuaciones que representan la
conducta dinamica del proceso en cuestion (8.2.1) y el grupo de ecuaciones diferenciales
lineales que representan a los coeficientes de sensibilidad (8.2.2).
De acuerdo con este comentario el procedimiento para evaluar los Sij consiste en
los siguientes pasos.
Para cada paso de integracion:
Integrar el modelo matematico
x = f (x, t, p), x(to ) = xo
Evaluar las derivadas parciales dadas por,
f (x, t, p)
x
f (x, t, p)
p
f (x, t, p)
f (x, t, p)
S+
x
p
.
2A 1 2B 2 C
esta dado por el siguiente conjunto de ecuaciones:
dCA
= k1 CA2
dt
dCB
= k1 CA2 k2 CB2
dt
dCC
= k2 CB2
dt
con k1 = 1, k2 = 2. Inicialmente el reactor se carga con solo 1 mol/lt del reactivo A.
Determinar la manera como los coeficientes cineticos k1 , k2 afectan a las conversiones
CA , CB y CC durante los primeros 20 minutos de operacion del reactor.
Evaluando el jacobiano,
2k1 CA
0
0
f
= 2k1 CA 2k2 CB 0
x
0
2k2 CB 0
c
2001,
Antonio Flores T.
UIA-Santa Fe
Clase 8
75
Sensibilidad Param
etrica
2
= CA CB2
p
0
CB2
definiendo,
-
CA
CA
, S12 =
S11 =
k1
k2
.
.
CB
CB
S21 =
, S22 =
k1
k2
.
.
CC
CC
S31 =
, S32 =
k1
k2
la matriz de sensibilidad parametrica S esta dada por:
S 11 S 12
2k1 CA
0
0
S11 S12
CA2
0
o bien,
S 11 S 12
2k1 CA S11
2k1 CA S12
CA2
0
2
CB2
S21 S 22 = 2k1 CA S11 2k2 CB S21 2k1 CA S12 2k2 CB S22 + CA
2k
C
S
2k
C
0
C
S
S31 S32
2 B 21
2 B 22
B
,
,
,
,
,
,
,
,
,
CA (0) = 1
CB (0) = 0
CC (0) = 0
S11 (0) = 0
S12 (0) = 0
S21 (0) = 0
S22 (0) = 0
S31 (0) = 0
S32 (0) = 0
UIA-Santa Fe
Clase 8
76
Sensibilidad Param
etrica
1
C
a
C
b
Cc
0.9
0.8
Concentration
0.7
0.6
0.5
0.4
0.3
0.2
0.1
10
Time
12
14
16
18
20
clear all
clc
global k1 k2
k1 = 1; k2=2;
time = linspace(0,20);
x0 = [1 0 0 0 0 0 0 0 0];
[t,x] = ode15s (odessa_demo,time,x0);
ca
cb
cc
s1
s2
s3
s4
s5
s6
=
=
=
=
=
=
=
=
=
x(:,1);
x(:,2);
x(:,3);
k1*x(:,4); % "Scaled" sensitivity coefficients
k2*x(:,5);
k1*x(:,6);
k2*x(:,7);
k1*x(:,8);
k2*x(:,9);
figure(1)
plot(t,ca,t,cb,t,cc), ylabel(Concentration), xlabel(Time)
legend (C_{a},C_{b},C_{c})
c
2001,
Antonio Flores T.
UIA-Santa Fe
Clase 8
77
Sensibilidad Param
etrica
K1
0.3
C
a
C
b
C
0.2
Sensitivity Coefficients
0.1
0.1
0.2
0.3
0.4
10
Time
12
14
16
18
20
figure (2)
plot(t,s1,t,s3,t,s5), ylabel (Sensitivity Coefficients),xlabel(Time)
legend(C_{a},C_{b},C_{c}), title(K_{1})
figure (3)
plot(t,s2,t,s4,t,s6), ylabel (Sensitivity Coefficients), xlabel(Time)
legend(C_{a},C_{b},C_{c}), Title(K_{2})
UIA-Santa Fe
Clase 8
78
Sensibilidad Param
etrica
K2
0.15
C
a
C
b
C
0.1
Sensitivity Coefficients
0.05
0.05
0.1
0.15
0.2
10
Time
12
14
16
18
20
s1 = x(4);
s2 = x(5);
s3 = x(6);
s4 = x(7);
s5 = x(8);
s6 = x(9);
%
% Math model
%
ca2 = ca*ca;
cb2 = cb*cb;
f (1) = -k1*ca2;
f (2) = k1*ca2-k2*cb2;
f (3) = k2*cb2;
%
% Sensivity states
%
t1
= 2*k1*ca;
t2
= 2*k2*cb;
f (4) = -t1*s1-ca2;
f (5) = -t1*s2;
f (6) = t1*s1-t2*s3+ca2;
f (7) = t1*s2-t2*s4-cb2;
f (8) = t2*s3;
f (9) = t2*s4+cb2;
f
= f;
c
2001,
Antonio Flores T.
UIA-Santa Fe
Clase 8
Sensibilidad Param
etrica
c
2001,
Antonio Flores T.
UIA-Santa Fe
79
Captulo 9
Problemas de Reactores Isot
ermicos
Preguntas para Discutir
Explicar como opera un reactor tanque agitado.
Que significa el suponer que el contenido de un reactor esta perfectamente mezclado ?
A que nos referimos cuando decimos que un RCTA es un sistema de propiedades
homogeneas ?
Que es un sistema de propiedades heterogeneas ?
Que es el estado estacionario de un proceso ?
Porque se plantea el modelo matematico de un RCTA en terminos solo de las
corrientes de entrada y salida del reactor ? Porque es no necesario tomar un
elemento diferencial de volumen para derivar el modelo de este reactor ?
Cual es la diferencia entre la operacion dinamica y en estado estacionario de un
proceso ?
Que es el analisis de sensibilidad ? Cual es su utilidad ? Como se realiza ?
80
Cap. 9
81
1. Derivar el modelo matematico (en forma dinamica y bajo condiciones estacionarias de operacion) de un reactor tanque agitado isotermico donde se lleva a cabo
la reaccion A + B C suponiendo los siguientes tipos de cineticas:
rA = kCA2 CB
rA = kCA CB2
2. Realizar el analisis de sensibilidad, en estado estacionario, de los 2 modelos derivados en el ejercicio anterior.
3. Se dispone de un reactor tanque agitado donde se lleva a cabo la reaccion A +
3B C. La corriente de alimentacion esta formada por un flujo de 200 lts/min
con una concentracion de 100 y 90 mol/lt de los reactivos A y B, respectivamente.
El reactor opera a 300 o K y tiene un volumen de 100 lts. El factor preexponencial
es igual 50000 lts/(mol-min), y la energa de activacion es 6500 cal/mol.
Calcular la concentracion de la corriente de productos que abandona el reactor.
Determinar la temperatura a la que debe operar el reactor para que consumir
25 % del reactivo A.
suponer que el reactor opera en estado estacionario.1
4. Derivar el modelo matematico de un reactor tubular isotermico e isobarico en
estado estacionario donde se lleva a cabo la reaccion A + 2B C suponiendo
los siguientes tipos de cineticas:
rA = kCA2 CB
rA = kCA CB2
5.
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 9
82
(mol/lt) T (o K)
100
300
200
600
100
800
seleccionar alg
un punto de operacion y justificar su eleccion. La cinetica de la
reaccion esta dada por: rA = kCA3 , el factor preexponencial A es igual a 105
lt2 /(mol2 -min) y la energa de activacion E es igual a 2500 cal/mol.
9. En la operacion industrial de algunos reactores semibatch es practica com
un
emplear perfiles optimos de flujos volumetricos de alimentacion. Considere el
reactor semibatch planteado en el problema anterior. Proponer alg
un perfil de
alimentacion tal que la cantidad consumida del reactivo A sea maxima usando
la menor cantidad posible del reactivo B en el menor tiempo posible. La u
nica
variable que se puede modificar es el flujo volumetrico alimentado. Reportar los
perfiles de concentracion de todos los compuestos involucrados durante todo el
periodo de operacion del reactor.
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 9
83
A/B
(1)
(2)
(3)
A/B
A/B
(4)
(5)
10. Seleccionar, entre los esquemas mostrados, el mejor esquema de reaccion para
efectuar la reaccion del problema 1.
en la siguiente tabla se muestran las condiciones de alimentacion y de arranque
para cada sistema de reaccion (el subndice f se refiere a las condiciones de la
corriente de alimentacion, mientras que el subndice o se refiere a las condiciones
de proceso en el interior del reactor al arranque de este).
Parametros Unidades (1) (2) (3) (4) (5)
QAf
(lt/min) 10
10
10
QBf
(lt/min)
10
10 10
CAf
(mol/lt) 300
300
300
CBf
(mol/lt)
80
80 80
CAo
(mol/lt)
300 150 150 150
CBo
(mol/lt) 80
40 40 40
Vo
(lts)
30 30 60 60 60
11. Para efectuar las siguientes reacciones,
CA2
A + B 3P1 , rA = k1
CB
k1
k
A 2 P2 , rA = k2 CA
se dispone del siguiente reactor tubular con alimentacion frontal y lateral:
P1 es el producto deseado y P2 el producto indeseado del anterior esquema de
reaccion. La corriente de alimentacion del reactivo A esta formada por 60 lt/min
con concentracion de 800 mol/lt. La alimentacion frontal del reactivo B esta compuesta por 355 lt/min con concentracion de 50 mol/lt. Para cualquiera de las
3 alimentaciones laterales la concentracion del reactivo B es tambien 50 mol/lt.
Las constantes de velocidad de reaccion k1 , k2 son .09 min1 y .03 min1 , respectivamente.
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 9
84
B
Productos
B
A
(1)
(2)
(3)
(4)
Seleccionar valores de los flujos volumetricos Q1 , Q2 y Q3 , asi como de los volumenes V1 , V2 , V3 y V4 de tal forma que la selectividad SP 1P 2 sea maxima.
12. Resolver nuevamente el problema anterior si en vez de usar un reactor tubular se
emplea la siguiente secuencia de reactores tanque agitado colocados en serie,
A
B
(1)
(2)
(3)
(4)
la alimentacion al primer reactor tanque agitado esta formado por las corrientes
de alimentacion frontal descritas en el problema anterior para el reactor tubular.
13. Se desea efectuar la reaccion:
A + 2B C
en un sistema de dos reactores conectados en serie. Los reactores a considerar son
tanque agitado y tubular ambos operando en estado estacionario. Seleccionar la
combinacion de reactores, de las cuatro posibles, que requiere el menor volumen
total para efectuar la reaccion descrita antes. La cinetica de la reaccion esta dada
por:
rA = kCA CB2
la conversion que se debe lograr en el primer reactor de la serie es 60 % y en
el segundo reactor se requiere una conversion global del 90 %. Las condiciones
de la corriente de alimentacion del primer reactor de cualquier combinacion son:
CAo =20 mol/m3 , CBo = 60 mol/m3 , FAo = 100 mol/min, FBo = 300 mol/min,
To =300 o C. El valor de k reportado a 300 o C es .001 m6 /(min-mol2 ).
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 9
85
Kubickova,Z., Kubicek,M., Marek, M., Fed-Batch Operation of Stirred Reactors, Chemical Engineering Science, 42,2,327-333 (1987).
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 9
86
Productos
UIA-Santa Fe
Cap. 9
87
A + B 1 3C , rA = k1 CA CB
k
D + C 2 P ,
rD = k2 CD CC2
ABC
(1)
(2)
A 1 B
k
B + C 2 A + C
k
2B 3 C + B
el modelo dinamico que representa la forma como las concentraciones de los
compuestos que participan en la reaccion cambian con respecto al tiempo esta
dado por:
dCA
= k1 CA + k2 CB CC
dt
c
2001,
Antonio Flores T.
UIA-Santa Fe
Cap. 9
88
= k1 CA k2 CB CC k3 CB2
= k3 CB2
c
2001,
Antonio Flores T.
UIA-Santa Fe