E j emplo 4-1 D eter ter mi nación de la k a par par tir ti r de dato datoss inte i nterr mitentes mitentes
Se desea diseñar un CSTR para producir 200 millones de libras de etilenglicol al año, hidrolizando óxido de etileno. Sin embargo, antes de realizar el diseño, es necesario llevar a cabo y analizar un experimento en un reactor intermitente para determinar la constante específica de velocidad de la reacción, k. Como la reacción será efectuada isotérmicamente, la velocidad de reacción específica deberá determinarse sólo a la temperatura de reacción del CSTR. A altas temperaturas, hay significativa formación de subproductos, mientras que a temperaturas inferiores a 40°C la reacción no' se realiza a velocidad significativa; en consecuencia, se eligió una temperatura de 55°C. Como el agua suele estar presente en exceso, su concentración puede considerarse constante en el curso de la reacción. La reacción es de primer orden para el óxido de etileno
En el experimento de laboratorio se mezclaron 500 ml de una solución 2 M (2 kmol/m3) de óxido de etileno en agua con 500 ml de agua, que contenía 0,9% en peso de ácido sulfúrico como catalizador. La temperatura se mantuvo a 55°C. La concentración de etilenglicol se registró en función del tiempo (tabla E4-1.1). Usando los datos de la tabla E4-1.1, determine la velocidad de reacción específica a 55°C TABLA E4 1.1. DATOS DE CONCENTRACIÓN CONTRA TIEMPO Tiempo (min) 0,0 0,5 1,0 1,5 2,0 3,0 4,0 6,0 10,0
Concentración de etilenglicol (kmol/m3)a 0,000 0,145 0,270 0,376 0,467 0,610 0,715 0,848 0,957
l kmol/m3 = 1 mol/dm3 = 1 mol/L.
a
Análisis
En este ejemplo usaremos el algoritmo para resolución de problemas (de la A a la G) dado en el CDROM y en el sitio www.engin.umich .edu/~problemsolving . Quizás usted desee seguir este algoritmo para resolver otros ejemplos ejemplos de este capítulo capítulo y problemas que se dan al al final del mismo. mismo. Sin embargo, embargo, para ahorrar espacio no lo repetiremos repetiremos para otros problemas de ejemplo. A. Formulación del problema. problema. Determine la velocidad de reacción especifica, B. Diagrama
k A .
C. Identificar C.1. Teorías pertinentes Balance molar: molar:
dN A
Ley de velocidad: velocidad: C.2. Variables
dt
r A
r AV
k AC A
Conocidos: concentración de etilenglicol en función del tiempo Desconocidos: Desconocidos:
1. Concentración de óxido de etileno en función del tiempo C A ¿? 2. Velocidad de reacción específica, k ¿? 3. Volumen del reactor, V ¿? C.4. Entradas y salidas: reactivos alimentados en una sola etapa al reactor intermitente C.5. Información faltante: Ninguna; aparentemente no habrá que consultar otras fuentes.
A
D. Suposiciones y aproximaciones
Suposiciones 1. Bien mezclado 2. Todos los reactivos entran de manera simultánea 3. No hay reacciones secundarias secundarias 4. Tiempo de llenado despreciable 5. Operación isotérmica Aproximaciones Aproximaciones
1. Agua en exceso, de manera que su concentración sea prácticamente constante (es decir C C ) B
B0
E. Especificación. Especificación. El problema no está sobreespecificado ni subespecificado F. Material relacionado. relacionado. En este problema se usan balances molares desarrollados en el capítulo
1 para reactor intermitente, así como las leyes de estequiometría y velocidad desarrolladas en el capítulo 3 G. Usar un algoritmo. Para una reacción isotérmica, usar el algoritmo de ingeniería de reacciones químicas de las figuras 4-1 y 4-2. Solución
C. Identificar C.1. Teorías pertinentes Balance molar: molar:
dN A
Ley de velocidad: velocidad: C.2. Variables
dt
r A
r AV
k AC A
Conocidos: concentración de etilenglicol en función del tiempo Desconocidos: Desconocidos:
1. Concentración de óxido de etileno en función del tiempo C A ¿? 2. Velocidad de reacción específica, k ¿? 3. Volumen del reactor, V ¿? C.4. Entradas y salidas: reactivos alimentados en una sola etapa al reactor intermitente C.5. Información faltante: Ninguna; aparentemente no habrá que consultar otras fuentes.
A
D. Suposiciones y aproximaciones
Suposiciones 1. Bien mezclado 2. Todos los reactivos entran de manera simultánea 3. No hay reacciones secundarias secundarias 4. Tiempo de llenado despreciable 5. Operación isotérmica Aproximaciones Aproximaciones
1. Agua en exceso, de manera que su concentración sea prácticamente constante (es decir C C ) B
B0
E. Especificación. Especificación. El problema no está sobreespecificado ni subespecificado F. Material relacionado. relacionado. En este problema se usan balances molares desarrollados en el capítulo
1 para reactor intermitente, así como las leyes de estequiometría y velocidad desarrolladas en el capítulo 3 G. Usar un algoritmo. Para una reacción isotérmica, usar el algoritmo de ingeniería de reacciones químicas de las figuras 4-1 y 4-2. Solución
1. El balance molar para un reactor intermitente bien mezclado es 1 dN A V
dt
r A
(E4-1.1)
2. La ley de velocidad es r kC
A
(E4-1.2)
A
Como hay agua en exceso, la concentración de agua en cualquier tiempo, t, es prácticamente la misma que la concentración inicial y la ley de velocidad es independiente de la concentración de H 2O C C 0 3. Estequiometría. Fase líquida, sin cambio ca mbio de volumen, V V (tabla E4-1.2) B
B
Especies CH2CH2O H 2O
0
TABLA E4-1.2 TABLA ESTEQUIOMÉTRICA Símbolo Inicial Cambio Remanente A N X N N N 1 X B N N X N N X A 0
A 0
A0
B
A0
A
A 0
A0
B
B
Concentración C C 0 1 X C C X A
A
A0
B
B
C B C A0 B C B0
CH OH 2
C
2
0
N A 0 X
N T 0
Recuerde que C A
B
N A0 X
NT
NT 0
N A
V
C A0 X
B
N B0 N A0
)
N A V 0
dt
CC
N A0X
es el número inicial de moles de A a B (es decir,
1 dN A
V0
NC
d NA
/V 0
dt
dC A
dt
4. Combinando la ley de velocidad y el balance molar, tenemos
dC A
r A
dt
kC A
(E4-1.3)
5. Evaluación. Para operación isotérmica, k es constante, de modo que podemos integrar esta ecuación (E4-1.3), reordenando:
dC A C A
kdt
usando la condición inicial que dice que t 0 , sea C C . La concentración inicial de A tras mezclar los dos volúmenes juntos es de 1.0 kmol/m3 (1 mol/L) . Integrando, se obtiene A
C A
dC A
C A 0
C A
ln
t
0
kdt
k
t
A0
dt
0
C A0 kt C A
La concentración de óxido de etileno en cualquier tiempo t es C A C A0 e kt
(E4-1.4) (E4-1.5)
La concentración de etilenglicol en cualquier tiempo t puede obtenerse de la estequiometria de reacción: A
B
C
NC N A0 X N A0 N A
Para reacciones en fase líquida V CC
NC V
V 0
N C C A0 C A C A 0 1 e V 0
kt
(E4-1.6)
Reacomodando y calculando el logaritmo en ambos lados, se obtiene ln
C A0 C C kt C A0
(E4-1.7)
Vemos que la gráfica de ln C 0 C /C 0 , en función de t, será una línea recta con pendiente -k. Usando la tabla E4-1.1, podemos construir la tabla E4-1.3 y emplear Excel para graficar ln C 0 C / C * 0 en función de t. A
A
C
C
A
A
TABLA E4-1.3. DATOS PROCESADOS C A0
C C
t
C C
(min)
(kmol/m3)
C A0
0,0 0,5 1,0 1,5 2,0 3,0 4,0 6,0 10,0
0,000 0,145 0,270 0,376 0,467 0,610 0,715 0,848 0,957
1,000 0,855 0,730 0,624 0,533 0,390 0,285 0,152 0,043
C A0 C C C A0 0,0000 -0,1567 -0,3147 -0,4716 -0,6292 -0,9416 -1,2553 -1,8839 -3,1466
ln
Con la pendiente de la gráfica de ln C 0 C /C 0 contra como se muestra en la figura de Excel E4-4.1. A
C
A
t
, podemos encontrar
0.0000 0.0
2.0
4.0
6.0
8.0
10.0
-0.5000 -1.0000 ) 0 A
C / ) c -1.5000 C 0 A
C -2.0000 ( ( n l
-2.5000
y = -0.3144x R² = 1
-3.0000 -3.5000
t (min)
Figura E4-1.1 Gráfica de los datos en Excel.
12.0
k,
Pendiente k
k
0,311min
0,311min
1
1
La ley de velocidad se transforma en: r A
kC A
r A 0,311min
1
C A
Ahora podemos usar la ley de velocidad para diseñar un CSTR industrial. E jemplo 4-2 Producción de 200 millones de libras al año en un CSTR
Se produjeron cerca de 12,2 mil millones de toneladas métricas de etilenglicol (EG) en el 2000, el cual fue el 26° producto químico más fabricado en la nación en ese año, con base en libras globales. Aproximadamente la mitad del etilenglicol se emplea como anticongelante y la otra mitad en la manufactura de poliésteres. Dentro de la categoría de poliésteres, el 88% se usó para fibras y el 12% para manufactura de envases y películas. El precio de venta en el 2004 del etilenglicol fue de 0,28 centavos de dólar por libra. Se desean producir 200 millones de libras al año de etilenglicol. El reactor operará isotérmicamente. Una solución de 1 lb mol/pies 3 de óxido de etileno (OE) en agua se alimenta al reactor (de la figura E4-2.1) junto con una solución volumétricamente igual de agua que contiene 0,9% en peso de catalizador H2SO4. La constante específica de velocidad de reacción es de 0,311 min -1), según se determinó en el ejemplo 4-1. Mukesh y Wasteel dieron recientemente directrices prácticas para aumentar la escala del reactor. (a) Para lograr una conversión del 80%, determine el volumen necesario del CSTR. (b) Si dos reactores de 800 galones estuvieran ordenados en paralelo, ¿cuál sería la conversión correspondiente? (e) Si dos reactores de 800 galones estuvieran ordenados en serie, ¿cuál sería la conversión correspondiente? Solución Suposición : El único producto de reacción que se forma es el etilenglicol (EG).
Figura E4-2.1 CSTR sencillo
La tasa especificada de producción de etilenglicol (EO) en lb mol/min es F C 2 108
1año 1día 1h 1lb mol lb lb mol 6,137 año 365días 24h 60min 62lb min
Por la estequiometría de la reacción FC
FA 0 X
encontramos el flujo molar que se requiere para el óxido de etileno, para una conversión del 80%: F A0
F A0
F C X
6,137
lb mol
min 0,8
7,67
lb mol min
g mol lb mol 453,59237 g 1 min 57,997 min 1lb 60s s
7,67
(a) Ahora calcularemos el volumen del
CSTR para lograr una conversión del 80% usando el
algoritmo CRE. 1. Ecuación de diseño:
Para un CSTR: V
F A0 X
(E4-2.1)
kC A
(E4-2.2)
r A
2. Ley de velocidad r A
3. Estequiometría. Fase líquida C A
F A
0
F A
: 0 1 X 0
C A 0 1 X
0
(E4-2.3)
4. Combinando:
Introduciendo (E4-2.3) en (E4-2.2):
1 X
r A kC A kC A0
reemplazando en (E4-2.1): V
F A0 X r A
FA0 X
kC A0 1
X
además sabemos de (E4-2.3) que 0 V
0 X k 1 X
F A 0
/ C A 0 , entonces:
(E4-2.4)
5. Evaluando
El flujo volumétrico de entrada de la corriente A, con es
A 0
F A 0 C A 01
7 ,67 lb mol / min 1lbmol/pie 3
Por el enunciado del problema
B0
7,67
C A01
1lbmol/pie 3 antes de la mezcla
pie 3 min
A 0
pie3 1lbmol lbmol lb FB0 B0C B 01 7,67 26,6 62, 4 min pie3 18lb min
El flujo volumétrico global de entrada del líquido es
0
A 0
B0
15,34
pies 3 (7,24 dm3/s) min
Sustituyendo en la ecuación (E4-2.4), recordando que
k
0,311min
1
, da
pies3 15,34 0,80 min 0 X V 197,3pies 3 1480gal 5,6m 3 1 k 1 X 0,311min 1 0,80
Se requiere un tanque de 5 pies de diámetro y aproximadamente 10 pies de alto para alcanzar una conversión del 80%. (b) CSTR en paralelo. Para dos CSTR en paralelo de 800 galones (como se muestra en la figura E42.2), alimentando 7,67 pies3/min 0 / 2 a cada reactor, la conversión lograda puede calcularse
reacomodando la ecuación (E4-2.4) Es necesario partir de una ecuación parecida: V
V 0 V
0
X
k 1 X
k
(E4-2.4)
X k 1 X
k
X
0
1
X
k 1 X X k kX X k X
kX
k X 1 k
para obtener: X
donde
k 1 k
(E4-2.5)
1pies3 1 800 gal 13, 94 min 7,48gal 7,62 pies3 min 0 2 V
El número de Damkohler es Da k 13, 94 min
Sustituyendo la ecuación (E4-2.5)
0,311 min
4, 34
Figura E4·2.2 CSTR en paralelo.
k Da X
Da 1 Da
4, 34 1 4, 34
0,81
La conversión que se logra a la salida de cada CSTR en paralelo es del 81%. (c) CSTR en serie. Si se ordenan los reactores de 800 galones en serie, la conversión en el
primer reactor [véase ecuación (E4-2.5)] es X 1
1 k 1 1 k
(E4-2.6)
Donde Entra todo, reactor por reactor: 1pies3 1 1 800 gal 6, 97 min 7,48gal 15,34 pies3 min 01 V 1
El número de Damkohler es Da1 X 1
1 k 6, 97 min 2, 167
1 2, 167
0,311
2, 167 3, 167
min
2, 167
0,684
Para calcular la conversión en la salida de segundo reactor, recordemos que 01 02 ; por lo tanto, 1
2
V1
V2
V y
El balance molar para el segundo reactor es Entradas-Salidas+Generación=0 F A1 FA2 rA2V 0
Basando la conversión en el número global de moles que reaccionan hasta un punto, por mol de A alimentado al primer reactor, F 1 F 0 1 X 1 y F 2 F 0 1 X 2 Reacomodando A
A
A
A
V
F A 1
F A 2
r
F A 0 1 X1
F A 0
A0
2
r
A 2
V
F 1 X A2
FA 0 X
FA 0
FA 0 X2
r A 2
X2
F A0
X 2
X1
r A 2
además:
r A 2
kC A 2
k
F A 2
k
F A0 1
0
kC A 0 1
X2
0
Combinando el balance molar para el segundo reactor [véase ecuación (2-24)] con la ley de velocidad, obtenemos V
F A0 X2
X1
r A2
1
C A00 X2
kC A0
X1
X2
X k1
0
2
X1
X2
(E4-2.7)
Despejando la conversión que sale del segundo reactor, se obtiene X 2
X1
Da
1 Da
X1
k
1 k
0,684 2,167 1 2, 167
0,90
Podríamos haber obtenido el mismo resultado empleando la ecuación (4-11): X 1
X 2
1
1 Da
1
n
1
1
1 2,167
2
1
1 k
n
(4-11)
0,90
Es posible producir 200 millones de libras de EG al año usando dos reactores de 800 galones (3,0 m3) en serie.
Figura E4-2.3 CSTR en serie.
Dos CSTR de igual tamaño colocado en serie (como se muestra en la figura E4-2.3) darán una conversión superior a dos CSTR en paralelo del mismo tamaño, si el orden de reacción es mayor de cero. E jemplo 4-3 Producción de 300 millones de libras al año de etileno en un PF R: diseño de un PF R a escala industrial
El etileno ocupa el cuarto lugar en Estados Unidos en lo que respecta al total de producto químico sintetizado al año y es el producto químico orgánico que ocupa el primer lugar de producción anual. Se produjeron más de 50 mil millones de libras en el 2000, que se vendieron a 0,27 centavos de dólar por libra. El 67% del etileno producido se emplea en la fabricación de plásticos, el 20% para óxido de etileno, el 16% para bicloruro de etileno y etilenglicol, el 5% para fibras y el 5% para disolventes. Determine el volumen del PFR necesario para producir 300 millones de libras de etileno al año por desintegración catalítica de una corriente de alimentación de etano puro. La reacción es irreversible
y sigue una ley de velocidad elemental. Se desea lograr una conversión del 80% de etano haciendo que el reactor funcione isotérmicamente a 1100 K Y una presión de 6 atm. Solución: C2 H 6
C2H 4
H 2
Sean A C2 H6 ; B C2 H4 y C H 2 . En símbolos,
A B C
Como deseamos que el lector se familiarice tanto con unidades métricas como con unidades del sistema inglés, trabajaremos algunos ejemplos con unidades inglesas. El flujo molar del etileno que sale del reactor es F B
300 10
lb
6
año
1año
365días
1días
24h
1h
3600s
lbmol
28lb
A continuación se calcula la velocidad de alimentación molar de etano, mol/s de etileno cuando se logra una conversión del 80%. FB
0, 340
F A 0
lbmol s
, para producir 0,34 lb
FA0 X
0,34
F A0
0,8
0,425
lbmol s
1. Ecuación de diseño del PFR: dX
F A0
dV
(2-15)
r A
Reacomodando e integrando para el caso en que no hay caída de presión y la operación es isotérmica, se obtiene V
F A0
X
0
dX
r A
(E4-3.1)
2. Ley de velocidad: con k 0,072s -1 en 1000K La energía de activación es de 82 kcal/g mol. r A
kC A
(E4-3.2)
3. Estequiometria. Para operación isotérmica y caída de presión despreciable, la concentración de etano se calcula como sigue: Fase gaseosa, T y P constantes: 0
C A C C
F T 0 1 X F T 0
F A
F A0 1 X
1 X C A0 0 1 X 1 X
C A0 X
1 X
(E4-3.3)
(E4-3.4)
4. A continuación combinaremos las ecuaciones (E4-3.1) a (E4-3.3) para obtener V F A0
X
0
X dX FA0 0 kC A
X 1 X dX dX FA 0 0 kC 1 X A0 1 X kC A0 1 X
V
F A0 C A0
X
0
1 X dX k 1 X
(E4-3.5)
5. Evaluar. Como la reacción se realiza isotérmicamente, podemos sacar k de la integral y recurrir al apéndice A.1 para efectuar la integración:
X
1 x dx 1x
0
1 ln
1 1x
x
(A-5)
entonces, integrando (E4-3.5): V
F A 0
kC A 0
X
1 X dX
0
1X
F A 0 kC A 0
1 ln X 1 1X
(E4-3.6)
6. Evaluación de parámetros: 6 atm 1 0,73pies3 atm lb mol °R × 1980°R RT 0 lbmol C A0 0,00415 dm 3 y A 0 1 1 1 1 1 C A0 y A0 C T 0
y A0 P0
La constante de velocidad k se da a 1000 K, por lo que es necesario calcular k en las condiciones de reacción, que son de 1100 K. E 1
1
E T2 T 1 T T 2 1 1 2 82000cal gmol 1100 1000 K 0,072 1 exp k T2 3,07 s s 1,987 cal g molK 1000K 1100K k T1 exp k T1 exp R T T R
k T2
(E4-3.7)
Sustituyendo la ecuación (E4-3.6) se obtiene V
V
F A 0 kC A 0
1 1 ln 1 X X 0,425lbmol s
1 1 1 ln 1 X 1 X 3,07 s 0,00415lbmol pies3
(E4-3.8)
1 V 33,36pies 3 2 ln X 1 X
Para
X
0, 8 :
V
1 0,8 80,7 pies 3 33,36pies3 2 ln 1 0,8
Se decidió usar un banco de tubos cédula 80 de 2 pulgadas en paralelo y de 40 pies de largo. Para el tubo cédula 80, el área de la sección transversal, A , es 0,0205 pies2. El número de tubos que necesitamos es de C
n
80,7pies3
0,0205pies 40pies 2
98,4
(E4-3.9)
Para determinar los perfiles de conversión y las concentraciones a lo largo del reactor, z, dividiremos la ecuación del volumen (E4-3.8) entre el área de la sección transversal, AC .
z
V
AC
(E4-3.10)
La ecuación (E4-3.9), junto con A = 0,0205 pies2, así como las ecuaciones (E4-3.8) y (E4-3.3) se emplearon para obtener la figura E4-3.1. Al usar un banco de 100 tubos, tendremos el volumen del reactor necesario para preparar 300 millones de libras al año de etileno a partir de etano. En la figura E4-3.1 se muestran los perfiles de conversión y concentración a lo largo de cualquiera de los tubos. C
0.0040
0.9
z vs CC,CB
0.0035 ) 3 s e i 0.0030 p / l o m0.0025 b l ( n 0.0020 ó i c a 0.0015 r t n e 0.0010 c n o C 0.0005
0.8
z vs CA
0.7
z vs X
0.6
n
ó 0.5 i s r
e v o C
0.4 n 0.3 0.2 0.1
0.0000 0.000
0 1.802
4.008
6.727 10.117 14.423 20.059 27.795 39.365
Distancia a los largo del reactor z × 10 -2 (pies)
Figura E4·3.1 Perfiles de conversión y concentración.
Los resultados de los cálculos empleados para elaborar la figura E4-3.1 se muestran en la siguiente Tabla:
X
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8
pies3
V
z pies
(E4-3.8)
(E4-3.10)
0,0 3,7 8,2 13,8 20,7 29,6 41,1 57,0 80,7
0,00 180,19 400,81 672,68 1011,67 1442,35 2005,89 2779,50 3936,46
z 10
2
pies
0,000 1,802 4,008 6,727 10,117 14,423 20,059 27,795 39,365
C A
lbmol dm
(E4-3.3) 0,0042 0,0034 0,0028 0,0022 0,0018 0,0014 0,0010 0,0007 0,0005
3
CC
C B
lbmol dm
(E4-3.4) 0,0000 0,0004 0,0007 0,0010 0,0012 0,0014 0,0016 0,0017 0,0018
3
E jemplo 4-4 Cálculo de la caída de presión en un lecho empacado
Grafique la caída de presión para 60 pies de longitud en tubería de 1 ½ pulgadas cédula 40, empacada con partículas de catalizador cuyo diámetro es de ¼ de pulgada, si pasan 104,4 Ib m/h de gas por el lecho. La temperatura es constante a lo largo del tubo; es de 260°C. La fracción de vacío es de 45% y las propiedades del gas son similares a las del aire a esa temperatura. La presión de entrada es de 10 atmósferas. Solución
Al final del reactor,
z
L,
y la ecuación (4-34) se transforma en
2 z P y 1 0 P0 P0
1/2
(4-34)
1 2
2 L 1 0 P0
P P0
(E4-4.1)
donde es una constante que depende únicamente de las propiedades del lecho empacado y las condiciones de entrada. 0
0
G 1 150 1 1,75 G gc 0 DP 3 DP m
G
Ac
(4-25) (E4-4.2)
Para tubos de 1 ½ pulgadas cédula 40, A
0,01414 pies2 : 104,4lbm h lbm 7383,3 G 2 0,01414 pies h pies2 c
Para aire a 260°C y 10 atmósferas 0,0673lbm piesh
0 0
0,413lbm pies 3 104,4lbm h m 0,413lbm pies3 0
252,8pies3 h
Del enunciado del problema, 1 pulg 0,0208pies, 0,45 4 lbm pies gc 4,17 108 lbf h 2 DP
Sustituyendo estos valores en la ecuación (4-25), se tiene G 1 150 1 1,75G 3 gc 0 DP DP 7383 lbm pies2 h 1 0,45 0 4,17 108 lbm pies lbf h 2 0,413lbm pies3 0,0208pies 0,453 0
150 1 0 , 45 0 , 0673 lb m pies h lbm 1,75 7383,3 0,0208pies pies 2 h
(E4-4.3)
Término1 Término2 lb f h lbm lb f 266,9 +12 920,8 0 0,01244 164,1 2 pies h pieslbm pies 3
(E4-4.4)
Observamos que el término de flujo turbulento, Término 2, predomina. 1pies 2 1atm lbf 0 164,1 3 2 pies 144pulg 14,7 lbf /pulg2 0
y
0,0775
atm pies
25,8
2 L P 1 0 P0 P0
12
kPa m
(E4-4.5)
0,155 2×0,0775atm pies×60pies 1 10atm
P 0,265P0 2,65atm
268kPa
P P0 P 10 2,65 7,35
atm
12
(E4-4.7)
Ahora usaremos los datos para graficar los perfiles de presión y el flujo volumétrico. Recordando la ecuación (3-41), para el caso en que 0 y T T F P T 0 T 0 FT 0 P T 0
0
P0 P
0
y
(E4-4.6)
0
(3-41) (E4-4.8)
Las ecuaciones (4-34) y (E4-4.8) se emplearon para construir la tabla E4-4.1 TABLA E4-4.1 PERFILES DE P Y v z
(pies) 0 10 20 30 40 g50 60
y
(4-34) 1,000 0,919 0,831 0,732 0,617 0,475 0,265
P P0 y (atm) 10,00 9,19 8,31 7,32 6,17 4,75 2,65
pies
3
h
(E4-4.8) 252,78 274,98 304,29 345,54 409,94 532,56 952,99
Para c 120lbm pies3
2 0,0775 atm pies 2 0 0,00165lbm -1 3 2 c 1 Ac P0 120 lbm pies 1-0,45 0,1414pies 10atm
Los valores de la tabla E4-4.1 se emplearon para obtener la figura E4-4 .1.
10
1000
9
900
8
800
7
700 600
6
P 5 (atm) 4
(pies3/h)400
3
300
2
200
1
100
0
0
υ
0
20
40
500
0
60
20
40
60
z (pies)
z (pies) Figura E4-4.1 Perfiles P y υ.
E jemplo 4-5 E fecto de la caída de presión sobre el perfil de conversión
Reconsidere el lecho empacado del ejemplo 4-4 para el caso reacción de segundo orden 2 A B C
en 20 metros de tubo de 1 ½ pulgadas cédula 40 empacado con catalizador. Las condiciones de flujo y lecho empacado del ejemplo siguen siendo las mismas, excepto que se transforman a unidades SI, o sea, P0 10 atm = 1013 kPa , y Flujo volumétrico de entrada: 0 7, 15 m3 /h 252 pies/h Tamaño de las partículas del catalizador: D p 0, 006 m (aproximadamente ¼ de pulgada)]
Densidad del catalizador sólido:
c
1923 kg/m3 120 Ibm /pies3
Área de la sección transversal del tubo de 1 ½ pulgadas cédula 40: Parámetro de la caída de presión: a 0 25 ,8 kPa m Longitud del reactor: L 20 m
Ac
2
0, 0013 m
Modificaremos el tamaño de partícula para observar qué ocurre con el perfil de conversión. Sin embargo, asumiremos que la velocidad de reacción específica, k, no se ve afectada por el tamaño de partícula, suposición que sabemos, según se explica en el capítulo 12, sólo es válida para partículas pequeñas. (a) Primero, calcule la conversión en ausencia de caída de presión. (b) A continuación, calcule la conversión tomando en cuenta la caída de presión. (e) Por último, determine cómo variaría la respuesta aparente al inciso (b) si el diámetro de partícula del catalizador se duplicara. La concentración de entrada de A es de 0,1 kmol/m 3 y la velocidad de reacción específica es 12 m 6 k kmol.kg cat.h Solución
Usando la ecuación (4-38) kC A0 W W 1 2 0 X kC A 0 W W 1 1 0 2
(4-38)
Para la densidad volumétrica del catalizador b
c 1
1923 1 0, 45 1058 kg m3
El peso del catalizador en el tubo de 20 metros de 1 ½ pulgadas y cédula 40 es de kg Ac b L 0, 0013m 2 1058 3 20m 27, 5 kg m 6 27,5kg kC A 0 W 12m kmol 0, 1 4, 6 3 kmolkg cal h 0 m 7,15 m 3 h
W
(a)
Primero calcularemos la conversión para
P 0
(es decir,
0)
kC A0 W 0
X 1 (b)
kC A 0 W
4, 6 0,82 1 4, 6
(E4-5.1)
0
A continuación, calcularemos la conversión con caída de presión. Recordando la ecuación (4-29) y sustituyendo la densidad volumétrica 1 1058 kg m3 b
2 0
P0 Ac c 1
c
(4-29)
kPa 2 25,8 2 0 m 0,037kg -1 kg P0 Ac b 1013kPa 0,0013m 2 1058 3 m
(E4-5.2)
Entonces
0 , 037 27 , 5 W 1 1 0, 49 2 2 kC A0 W W 1 2 4, 6 0, 49 v0 2,36 X 0,693 kC A 0 W W 1 4, 6 0, 49 3, 26 1 1 2 v0
(E4-5.3)
(E4-5.4)
Observamos que la conversión predicha descendió de 82,2% a 69,3% por la caída de presión. No sólo resultaría embarazoso, sino también todo un desastre económico despreciar la caída de presión y obtener, como resultado, una conversión real significativamente menor. (c) Roberto el preocupón se pregunta: ¿Qué pasaría si aumentaremos el tamaño del catalizador
en un factor de 2? Vemos por la ecuación (E4-4.5) que el segundo término de la ecuación de Ergun predomina; o sea, 1,75
150 1 DP
(E4-5.5)
Por lo tanto, de la ecuación (4-25) 0
G 1 150 1 G 1,75 gc 0 DP 3 DP
(4-25)
tenemos que 0
1,75G2 1 0 gc DP 3
(E4-5.6)
Observamos para las condiciones dadas por la ecuación (E4-4.4) Término1 Término2 lb f h lbm lb f 0 0,01244 266,9 +12 920,8 164,1 pies 2 h pieslbm pies 3
(E4-4.4)
que el parámetro de caída de presión varía inversamente con el diámetro de la partícula. 1
0
DP
por lo tanto, 1
DP
Para el caso 2,
2
DP 2
1
DP 1 DP 2
2 DP 1
1
0, 037kg-1 2
0, 0185 kg-1
(E4-5.7)
Sustituyendo este nuevo valor de a en la ecuación (E4-5.4)
4, 6 1 X 2
0,0185 27,5 2
0,0185 27, 5 1 4, 6 1 2
3,43 4,43
0,774
Aumentando el tamaño de partículas, se reduce el parámetro de caída de presión; por lo tanto, aumentan la velocidad de reacción y la conversión. Sin embargo, en los capítulos 10 y 12 se explica que cuando los efectos de difusión interpartícula son importantes en las partículas del catalizador, no siempre ocurrirá un incremento de conversión al aumentar el tamaño de partícula. Para partículas de mayor tamaño, se requiere más tiempo para que determinado número de moléculas de reactivo y de producto se difundan entrando y saliendo de las partículas del catalizador donde experimentan la reacción (véase figura 10-6). En consecuencia, la velocidad de reacción específica disminuye al aumentar el tamaño de partícula, k ~ l/Dp [véase ecuación (12-35)], lo cual a la vez reduce la conversión. Para partículas de diámetro pequeño, la constante de velocidad, k, es grande y tiene su valor máximo, pero la caída de presión también es alta, produciendo menor velocidad de reacción. Con partículas de diámetro grande, la caída de presión es baja, pero también lo son la constante de velocidad, k, y la velocidad de reacción, por lo cual la conversión es baja. Por lo tanto, se observa baja conversión con diámetros de partículas grandes o pequeños y se alcanza un óptimo en etapas intermedias. Este óptimo se muestra en la figura E4-5.1. Véase problema P4-23.
Figura E4-5.1 Diámetro óptimo de partícula.
Para minimizar la caída de presión, ¿por qué no empacar el catalizador en un tubo de gran diámetro para reducir la densidad de flujo másico, G, reduciendo en consecuencia P ? Son dos los motivos por los cuales no se incrementa el diámetro del tubo: 1) Hay una mayor probabilidad de que el gas forme un canal sin pasar por la mayor parte del catalizador, ocasionando conversión baja (véase figuras 13-2 y 13-10); 2) la proporción entre el área superficial de transferencia de calor respecto del volumen del reactor (peso del catalizador) se reducirá, dificultando la transferencia de calor en reacciones altamente exotérmicas o endotérmicas. A continuación procederemos (ejemplo 4-6) a combinar la caída de presión con la reacción en un lecho empacado, cuando hay cambio de volumen con la reacción y, por lo tanto, no se puede obtener una solución analítica. E jemplo 4-6 Cálculo de X en un reactor con caída de presión
Se produjeron aproximadamente 7 mil millones de libras de óxido de etileno en Estados Unidos durante 1997. El precio de venta en ese año fue de 0,58 centavos de dólar la libra, lo cual equivale a un valor comercial de 4,0 mil millones de dólares. Más del 60% del óxido de etileno producido se emplea para fabricar etilenglicol. Los principales usos finales del óxido de etileno son como anticongelante (30%), en poliéster (30%), tensoactivos (10%) y disolventes (5%). Se desea calcular el peso del catalizador necesario para lograr una conversión de 60% al preparar óxido de etileno por oxidación catalítica de etileno con aire en fase gaseosa.
El etileno y el oxígeno se alimentan en proporción estequiométrica a un reactor empacado que opera isotérmicamente a 260°C. El etileno se alimenta a razón de 0,30 lb mol/s y presión de 10 atm. Se propone usar 10 bancos de tubería de diámetro de 1½ pulgadas y cédula 40 empacada con catalizador, se tienen 100 tubos por banco. En consecuencia, el flujo molar en cada tubo será de 3×10 -4 lb mol/s. Las propiedades del fluido reaccionante deben considerarse idénticas a las del aire a esta temperatura y esta presión. La densidad de las partículas del catalizador de ¼ de pulgada es de 120 lbm/pies 3 y la fracción de vacío en el lecho es de 0,45. La ley de velocidad es
r A
1 3
2 3
kPA PB
lb mol lbm cal h
con k 0,0141
lbmol a 260°C atmlbm cat h
Solución 1. Balance molar diferencial: F A0
dX dW
r A
(E4-6.1)
2. Ley de velocidad: r A
13
23
kPA PB
13
k CA RT
C
B
23
RT
(E4-6.2) (E4-6.3)
r A kRTCA1 3CB2 3
3. Estequiometría. Fase gaseosa, isotérmica 0 1 X P0 C A0 1 X P CA0 1 X y donde y 1 X 1 X P0 F C X 2 C B B A0 B y 1 X C A
F A
Para alimentación estequiométrica B CB
C A0 2
F B0
: P P0
(E4-6.4) (E4-6.5)
1
F A0
P
2
1 X y 1 X
Para operación isotérmica, la ecuación (4-30) se transforma en dy dW dy dW
2y
2y
1 X
T T0
(4.30)
1 X
(E4-6.6)
4. Combinando la ley de velocidad y las concentraciones: r A kRTCA1 3CB2 3 C 1 X r A kRT0 A0 y X 1 23
13
CA0 1 X y X 2 1
Factorizando 1 /2 y recordando que (E4-6.7) como sigue
P A 0 1 3
2 3
C A 0 RT0 ,
(E4-6.7)
podemos simplificar la ecuación 2 3
C 1 X C A0 1 X kRT0 A 0 y y 1 1 X X 2 3 2 3 C 1 X 1 1 X 1 y kPA 0 y r A kRT0 A 0 1 X 2 1 X 2 1 r A 2
2 3
1 X r A k y 1 X
(E4-6.8)
1 donde k kP A0 2
2 3
0,63kPA0
5. Evaluación de parámetros por tubo (es decir, dividir las velocidades de alimentación entre 1000): Etileno: 4
3 10 lbmol s 1,08lbmol h
1,5 10 lbmol s 0,54lbmol h
F A0
Oxigeno: F B0 I
Inertes
F I
4
N 2 :
1,5 10 4 lbmol s
4
F I 5,64 10
0,79molN2 0,21molO2
lbmol s 2,03lbmol h
Sumando: FT 0 FA 0 FB0 F1 3,65lbmol h y A 0
F A 0 1,08 0,30 F T 0 3,65
y A 0 0, 3 1
1 0,15 2 1
P A0 y A0 P0 3,0atm
1 k kP A0 2
Para calcular
23
0, 0141
lb mol atm lbm cat h
3 atm 0, 63 0, 0266
lb mol h lb m cat
: 2 0
P0 Ac c 1
(4-29)
se requiere la densidad de flujo másico, G. Los flujos másicos de cada especie a la entrada son m A0
1,08
mB0
0, 54
mI 0 2,03
lbmol h
28
32
lbmol h
lbm lbmol lbm lbmol
30, 24 lbm h
17, 28lbm h
lbmol lbm 28 56,84lbm h h lbmol
El flujo másico global es de mT 0
G
104,4
lbm h
mT 0 10,4lbm h lbm 7383,3 2 Ac 0,01414pies h pies 2
La densidad de flujo másico, la temperatura y la presión son iguales que en el ejemplo 4-4. En consecuencia, podemos usar el valor de 0 del ejemplo 4-4, para calcular 0 0,0775
atm pies
2 0 Ac
1
c P0
2 0,0755 atm pies 0,01414pies2 0,55 120 lbm cal pies3 10atm
0,0166 lbm cal
6. Resumen. Combinando las ecuaciones (E4-6.1) (E4-6.8), Y resumiendo dX k 1 X y dW F A 0 1 X 1 X dy dW
F A0
2y
k 0,0266
1,08
lbm cal
(E4-6.10)
lbmol
(E4-6.11)
h lbm cal
lbmol
0,0166
(E4-6.9)
h
.
(E4-6.12)
(E4-6.13)
0,15
(E4-6.14)
Tenemos las condiciones frontera W 0, X 0 y y 1, 0 , Y W f 60 lbm . Aquí estamos suponiendo un límite superior de integración de 60 lbm con la expectativa de lograr un 60% de conversión con este peso de catalizador. Si no se logra la conversión del 60%, supondremos un peso más alto y volveremos a efectuar los cálculos. Gran número de paquetes ordinarios de software para resolver ecuaciones diferenciales (es decir, resolvedores EDO), que en la actualidad son muy amigables para el usuario, se encuentran disponibles de forma comercial. Usaremos el Polymath para resolver los ejemplos del texto impreso. Con el Polymath, simplemente se introducen las ecuaciones (E4-6.9) y (E4-6.10), así como los valores de los parámetros correspondientes (ecuaciones (E4-6.11) a (E4-6.14) a la computadora con las condiciones frontera; además, se resuelven y presentan como se muestra en las figuras E4-6.1 y E46.2. Se introducen las ecuaciones (E4-6.9) y (E4-6.10) como ecuaciones diferenciales, por lo que los valores de los parámetros se fijan usando ecuaciones explícitas. La ley de velocidad puede introducirse como una ecuación explícita, para generar una gráfica de la velocidad de reacción conforme varía a lo largo del reactor, usando la función de Polymath para graficar. El CD-ROM contiene todos los programas de soluciones MatLab y Polymath que se emplean para resolver los problemas de ejemplo y también un ejemplo con Aspen. En consecuencia, se puede cargar el programa Polymath directamente del CD-ROM, que tiene programadas las ecuaciones (E4-6.9) a (E4-6.14), y correr el programa para diferentes valores de los parámetros. También es interesante observar lo que ocurre con el flujo volumétrico a todo lo largo del reactor. Recordando la ecuación (3-45),
0 1 X
P0 T P T0
0 1 X T P P0
T 0
(3-45)
Sea f la proporción entre el flujo volumétrico, , y el flujo volumétrico de entrada, 0 , en cualquier punto a lo largo del reactor. Para operación isotérmica, la ecuación (3-45) se transforma en
f
0
1 X
y
(E4-6.15)
La figura E4-6.2 muestra X, y (es decir, y = P/P o) y fa lo largo del reactor. Vemos que tanto la conversión como el flujo volumétrico aumentan a todo lo largo del reactor, mientras que la presión disminuye. Para reacciones en fase gaseosa, con órdenes mayores de cero, esta reducción de presión provocará que la velocidad de reacción sea menor que en el caso en el que no haya caída de presión. Introducimos los datos y ecuaciones a Polymath: d(X)/d(W)=-raprime/FA0 X(0)=0 W(0) = 0 W(f ) = 60 d(y)/d(W)=-alpha*(1+eps*X)/2/y y(0)=1 eps=-0.15 kprime=0.0266 FA0=1.08 alpha=0.0166 raprime=-kprime*(1-X)/(1+eps*X)*y f =(1+eps*X)/y rate=raprime
Exportamos los resultados a Excel y graficamos: 3 2.5 2
-rA:102 1.5 1 0.5 0 0
10
20
30
40
50
60
70
W
Figura E4-6.1 Perfil de velocidad de reacción
a lo largo del reactor empacado.
4 3.5 3
f X y
2.5 2 1.5 1 0.5 0 0
10
20
30
40
50
60
70
W
Figura E4-6.2 Reproducción gráfica de Excel
Por el perfil de conversión (que se muestra en la figura E4-6.2) o la tabla de resultados Polymath (que no se muestra en el texto, pero está disponible en el CD), encontramos que se logra una conversión del 60% con 44,5 lb m de catalizador en cada tubo. Observamos en la figura E4-6.2 que el peso del catalizador necesario para aumentar la conversión un 1 %, de 65% a 66% (3,5 lbm) , es de 8.5 veces mayor que el requerido (0,41 lb m) para elevar la conversión 1 % en la entrada del reactor. Además, durante el último incremento de 5% de la conversión, la presión disminuye de 3,8 atm a 2,3 atm. Este peso del catalizador es de 44,5 lb m/tubo y corresponde a una caída de presión de aproximadamente 5 atm. Si hubiéramos despreciado erróneamente la caída de presión, el peso del catalizador se hubiese encontrado integrando la ecuación [E4-6.9] con y 1 , para dar
W
W
1 X 1 ln k 1 X
F A 0
1 1 0, 15 ln 0, 15 0, 6 35, 3 lbm 0, 0266 1 0, 6 1, 08
(E4-6.16) (E4-6.17)
Si hubiésemos usado este peso del catalizador en el reactor, no habríamos tenido el suficiente catalizador para lograr la conversión deseada. Para este peso de catalizador (es decir, 35,300 lb m total, 35,3 lbm/tubo) la figura E4-6.2 da una conversión de sólo 53%. E jemplo 4-7 R eacción en fase gaseosa utilizando un microrr eactor en términos de flujo molar
La reacción en fase gaseosa 2 NOCl 2 NO
Cl2
se efectúa a 425°C y 1641 kPa (16,2 atm) . Se alimenta NOCl puro y ésta sigue una ley de velocidad elemental. Se desean producir 20 toneladas de NO al año en un sistema de microrreacción usando un banco de 10 microrreactores en paralelo. Cada microrreactor tiene 100 canales y cada canal mide 0,2 mm cuadrados y tiene 250 mm de largo.
Grafique los flujos molares en función del volumen a lo largo del reactor. El volumen de cada canal es 10 dm . 5
3
Información adicional
Para producir 20 toneladas al año de NO con una conversión del 85% se requeriría un flujo de alimentación de 0,0226 mol/s de NOCl o de 2,26 10 5 mol/s por canal. La constante de velocidad es
k
0, 29
dm
3
mol s
a 500K con E
24
kcal mol
Solución
Para un canal,
Aunque este problema en particular se resolvería por conversión, demostremos cómo puede resolverse usando flujos molares como variable en el balance de moles. Primero escribiremos la reacción en forma simbólica y después dividiremos entre el coeficiente estequiométrico del reactivo frontera, NOCl. 2 NOCl
2 NO
Cl2
2 A 2 B C
A B
1
C
2
1. Balances molares para las especies A, B Y C: dF A dV
r A
(E4-7.1)
r B
(E4-7.2)
r C
(E4-7.3)
dF B dV dF C dV
2. Ley de velocidad:
r A kCA2 , k 0,29
dm3 a 500K mols
(E4-7.4)
3. Estequiometria: Fase gaseosa con
T
T 0
y
P
P0 ,
entonces
F T
0
F T 0
a. Velocidades de reacción r A
r B
1
1
rB
rC
r C 1 2
r A 1 2
r A
b. Concentración Aplicando la ecuación (3-42) a las especies A, B Y C, las concentraciones son C A CT 0
con
FT
FA
F A F F , C B CT 0 B , CC C T C FT FT F T 0
FB
4. Combinando: la ley de
r A
FC
(E4-7.5)
velocidad en términos de los flujos molares tenemos
F kC T 2 0 A F T
2
Combinando todo dF A dV dFB dV dF C dV
F kC T 0 A F T
2
(E4-7.6)
2
F kC T 0 A F T
2
(E4-7.7)
2
F kC T 0 A F T
2
(E4-7.8)
2
5. Evaluando: C T 0
P0 RT 0
1641kPa kPadm 8,314 molK 3
698K 0, 286
mol dm
3
0, 286
mol cm
3
Al usar Polymath u otro resolvedor de EDO, no hay que combinar los balances molares, las leyes de velocidad y la estequiometría en realidad, como se hizo en el paso de combinación de ejemplos previos en este capítulo, pues el revolvedor de EDO hace todo eso. ¡Gracias por contar con el resolvedor de EDO! El programa Polymath y sus resultados se muestran tanto en la tabla E4-7.1 como en la figura E4-7.1 E jemplo 4-8 R eactor de membrana
Según el Departamento de Energía de Estados Unidos, se pudo lograr un ahorro energético de 10 billones de BTU al año, usando reactores catalíticos de membrana, para sustituir los reactores convencionales para reacciones de deshidrogenación como la deshidrogenación de etilbenceno para dar estireno:
y de butano para dar buteno: C 4 H 10
C4H8
H2
La deshidrogenación de propano es otra reacción que ha comprobado tener éxito en un reactor de membrana. C3 H 8
C
3
H6
H2
Todas las reacciones de des hidro gen ación anteriores pueden representarse simbólicamente como sigue A
B C
y ocurren del lado del catalizador en un RMICA. La constante de equilibrio para esta reacción es muy pequeña a 227°C (es decir, K 0,05 mol/dm3 ). La membrana es permeable a B (por ejemplo, H2), pero no a A y C. A gaseoso puro entra al reactor a 8,2 atm y 227°C a razón de 10mol/min. C
Tomaremos la velocidad de difusión de B que sale del reactor, por unidad de volumen del reactor, R , como proporcional a la concentración de B (o sea, R k C ). B
B
C
B
(a) Realizaremos balances molares diferenciales para A, B Y C, con la finalidad de obtener un
conjunto de ecuaciones diferenciales acopladas para resolverlas. (b) Graficaremos los flujos molares de cada especie en función del espacio tiempo. 3 (c) Calcularemos la conversión de V 400 dm . Información adicional : Aunque esta reacción es una reacción catalítica gas-sólido, usaremos la
densidad volumétrica del catalizador para escribir los balances en términos de volumen del reactor, en lugar del peso del catalizador (recuerde que r r ). Para la densidad volumétrica del catalizador b 1,5 g/cm 3 y un tubo de diámetro interno de 2 cm que contenga partículas de catalizador, la velocidad de reacción específica, k, y el coeficiente de transporte, kc, son k 0,7 min -1 y kC 0,2 min-1 , respectivamente. A
A
b
Solución
Elegiremos el volumen del reactor en vez del peso del catalizador como variable independiente en este ejemplo. El peso del catalizador, W, y el volumen del reactor se relacionan fácilmente a través de la densidad volumétrica del catalizador, (o sea, W bV ). Primero, realizaremos balances molares para el elemento del volumen V que se muestra en la figura 4-13(d). b
1. Balances molares: Balance para A en el lecho catalítico:
Flujode Flujode entrada - salida + Generación = Acumulación F A FA rA V 0 V V V
Dividiendo entre
V
y encontrando el límite cuando V 0 , se tiene
dF A dV
r A
(E4-8.1)
Balance para B en el lecho catalítico:
El balance para B se obtiene mediante la ecuación (4-4 1). dF B dV
rB
RB
(E4-8.2)
donde R B es el flujo molar de B a través de la membrana por unidad de volumen del reactor. El balance molar para e se efectúa de manera idéntica al de A y la ecuación resultante es dF C dV
r C
(E4-8.3)
2. Ley de velocidad:
C C r A k C A B C K C
(E4-8.4)
3. Transporte de salida del reactor. Aplicaremos la ecuación (4-42) para el caso en el que la
concentración de B del lado del barrido es cero, C BS = 0, para obtener R K C B
C
B
(E4-8.5)
donde kc es el coeficiente de transporte. En este ejemplo, asumiremos que la resistencia de la especie B que sale de la membrana es una constante y, en consecuencia, kc es una constante. 4. Estequiometría. Recuerde la ecuación (3-42) para el caso de temperatura y presión
constantes, en operación isotérmica y sin caída de presión (T = T 0, P = P0), se tiene lo siguiente Concentraciones: C A
C T 0
CB
C T 0
CC
C T 0
FT
FA
F A F T F B F T F C F T
(E4-8.6)
(E4-8.7)
(E4-8.8)
FB
FC
(E4-8.9)
Velocidades relativas: r A
1
r B 1
rB
r A
rC
r C 1
r A
Combinando y resumiendo:
(E4-8.10)
dF A
dF B
r k A
dV
C
F F B
C T 0
T
dF C dV
r
A
dV
r
A
F C F F r kC F K F F F F F F A
A
T0
B
C
T
C
T
T
T 0
T
A
B
C
6. Evaluación de parámetros: P 0
k
0,7 min , KC
830,6 kPa
C T 0
RT0
3
8,314k Pa dm
1
F A 0
10 mol
F B 0
F CD
mol K 500 K 3
0,2
0,05mol dm , kC 0,2min
mol dm
3
1
min
0
7. Solución numérica. Las ecuaciones (E4-8.1) a (E4-8 .10) fueron resueltas por Polymath y MatLab, otro resolvedor de EDO. A continuación se muestran los perfiles de flujo molar. En la tabla E4-8.l se muestran los programas Polymath y en la figura E4-8.l , los resultados de la solución numérica en las condiciones iniciales (de entrada). V
0:
F 0,
F A
A
FB
0,
FC
0
(c) Por la figura E4-8.1, vemos que el flujo molar de salida de A a 400 dm3 es de 4 mol/min, para lo cual la conversión correspondiente es F 0 F 10 4 0,6 X A
A
F A0
10
E jemplo 4-9 R eactor semicontinuo isotérmico con r eacción de segundo orden
La producción de bromuro de metilo es una reacción irreversible en fase líquida que sigue una ley de velocidad elemental. La reacción CNBr
CH 3 NH 2
CH 3 Br
NCNH 2
se efectúa isotérmicamente en un reactor semicontinuo. Se al imenta una solución acuosa de
metilamina (B) con una concentración de 0 .025 mol/dm3 con flujo volumétrico de 0.05 dm3/s a una solución acuosa de cianuro de bromo (A) contenida en un reactor recubierto de vidrio. Él volumen inicial del líquido en el recipiente debe ser de 5 dm3 con una concentración de cianuro de bromo de 0.05 mol/dm3. La constante de velocidad específica de la reacción es k
3
2,2 dm s mol
Encuentre las concentraciones de cianuro de bromo y bromuro de metilo, así como la velocidad de reacción en función del tiempo. Solución
Simbólicamente, escribiremos la reacción como sigue A B C D
Es una reacción elemental; por lo tanto, la ley de velocidad es r kC C Sustituyendo la ley de velocidad en las ecuaciones (4-56) y (4-58), se obtiene A
dC A dt dC B dt V
A
B
kC ACB
kC ACB
V0
v0 V v0 V
C A
C
B0
(E4-9.2)
CB
(E4-9.3)
t
v0
De manera similar, para C y D tenemos dN C dt dN C dt
rCV
0
(E4-9.4)
rAV
d CcV
dt
V
(E4-9.1)
(E4-9.5)
dCC dt
CC
dV dt
V
dCC
dt
v0 CC
(E4-9.6)
Entonces dCC
dt dC D dt
kC AC B
kC AC B
v0CC V
(E4-9.7)
(E4-9.8)
v0C D V
También podríamos calcular la conversión de A. X
X
N A0
N A
N A0
C A0V0
(E4-9.9)
C V
C A0V 0
A
(E4-9.10) Las condiciones (E4-9.2) a (E4-9.10) se resuelven fácilmente con ayuda de un revolvedor de EDO como Polymath (tabla E4-9.1)
Las concentraciones de cianuro de bromo (A) y metilamina se muestran en función del tiempo en la figura E4-9.1, en tanto que la velocidad se muestra en la figura E4-9.2.
Capítulo 5 Recolección y análisis de datos de velocidad
Determinación de la dependencia de la ley de velocidad del hidrógeno
La tabla E5-5.2, sugiere que la dependencia de r' CH respecto de P H no se puede representar con una ley de potencia. Comparando la prueba 4 con la 5 y la prueba 5 con la 6, veremos que la velocidad de reacción primero aumenta al incrementarse la presión parcial del hidrógeno y luego disminuye al aumentar P H2. Es decir, parece haber una concentración de hidrógeno en la que la velocidad es máxima. Un conjunto de leyes de velocidad, congruente con estas observaciones, es: r CH
l. A concentraciones bajas de H donde velocidad podría tener la forma
4
2
rCH
aumenta conforme
P H
2
se incrementa, la ley de
4
P H 1 2
(E5-5.6)
2. A concentraciones altas de r CH
4
H 2
r CH
donde
4
disminuye conforme
P H
2
aumenta,
1
1 P H
(E5-5.7) Nos gustaría encontrar una ley de velocidad que sea congruente con datos de velocidad de reacción a concentraciones tanto bajas como altas de hidrógeno. La aplicación de material del capítulo 10 sugiere que las ecuaciones (ES-S .6) Y (E5-5.7) pueden combinarse en la forma 2
r CH
4
P H 1 2
1 1 P H
(E5-5.8) En el capítulo 10 veremos que esta combinación y las leyes de velocidad similares que tienen concentraciones de reactivo (o presiones parciales) en el numerador y el denominador son comunes en catálisis heterogénea. 2
Veamos si la ley de velocidad resultante (E5-5.5) es cualitativamente consistente con la velocidad observada. P H
1. Para la condición 1: a rCH
2
bajo,
b P H
2
y la ecuación (E5-5.S) se reduce a
2
1
4
P H 1 2
(E5-5.9)
La ecuación (E5-5.9) es congruente con la tendencia al comparar las pruebas 4 y 5. 2. Para la condición 2: A a
P H
2
alto,
b P H
2
2
y la ecuación (E5-5.S) se reduce a
1
1
r CH
4
P P
1
H 2
2
P
2 1
(E5-5.10) donde 2 · La ecuación (E5-5.1O) es congruente con las tendencias, al comparar las pruebas 5 y 6. Combinando las ecuaciones (E5-5.S) y (E5-5.5) H 2
2
1
H 2
r CH
aPCO P H 1
4
2
1 bP H
2
(E5-5.11) Ahora utilizaremos el programa de regresión de Polymath para encontrar los valores de los parámetros a, y . Los resultados se muestran en la tabla E5-5.3. 2
2
1
La ley de velocidad correspondiente es 0,61
r CH
4
0,025 P H 2
1 2, 49 P H 2
(E5-5.12) Podríamos usar la ley de velocidad dada por la ecuación (E5-5.l2) tal cual, pero sólo hay seis datos, pero nos debería preocupar la posibilidad de extrapolar la ley de velocidad sobre un intervalo más amplio de presiones parciales. Podríamos tomar más datos o llevar a cabo un análisis teórico del tipo que se discute en el capítulo 10 para reacciones heterogéneas. Si asumimos que el hidrógeno experimenta adsorción disociativa sobre la superficie del catalizador, esperaríamos una dependencia del hidrógeno elevada a la ½ , Como 0.61 es cercano a 0.5, vamos a efectuar la regresión de datos de nuevo fijando . Los resultados se muestran en la tabla ES.SA.
La ley de velocidad ahora es 12
r CH
4
0,018 PCO P H 2 1 1,49 P H 2
r CH
donde se da en mol/gcat . s y la presión parcial en atm. Podríamos también haber fijado 1 1 2 y 2 1, 0 Y reordenar la ecuación (E5-5.11) en la forma 4
12
PCO P H
2
rCH
4
1 a
b a
P H
2
(E5-5.13)
