Por: Ing.

Edy Gonzalo Aspi Quispe Incotegas

Ejemplo de Smith 5-2, resuelto paso a paso

Ejemplo 1.8.
Fuente: Ejemplo 5-2, resuelta en Smith. Ingeniería de la cinética química.
Se planea diseñar una planta piloto para la producción de cloruro de alilo. Los reactantes constan
de 4 moles de propileno/mol de cloro y entran al reactor a 392 ºF. La unidad será un tubo vertical
de 2 plg de DI. Si la velocidad de la alimentación combinada es 0.85 lbmol/h.
Los reactantes se precalentarán a 392 °F por separado y se mezclarán a la entrada del reactor. A
esta temperatura baja, los problemas de explosiones durante el mezclado no son serios. El reactor
estará enchaquetado con etilenglicol a ebullición, para que la temperatura en el interior de la pared
sea constante e igual a 392 °F. El coeficiente de transferencia de calor en el interior es
( ).
Datos adicionales y notas. El desarrollo básico del proceso para el cloruro de alilo ha sido reportado
por Groll y Hearne y Fairbairn, Cheney y Cherniavsky. Se encontró que las tres reacciones principales
son la formación del cloruro de alilo.
)

la reacción de adición paralela que produce 1,2-dicloropropano,

)

y la cloración consecutiva del cloruro de alilo para dar 1,3-dicloro-1-propeno

)

Para simplificar el tratamiento cinético del problema, se considerarán únicamente las dos primeras
reacciones. En la Tabla siguiente se muestran los calores de reacción.
( ) 298 K
Reacción 1 -48000
Reacción 2 -79200
Se puede suponer que las capacidades caloríficas molares son constantes e iguales a los valores
de la Tabla siguiente:

Componentes ( )
Propileno (g) 25,3
Cloro (g) 8,6
Cloruro de hidrogeno (g) 1,2
Cloruro de alilo (g) 28,0
1,2-dicloropropano (g) 30,7
No se ha publicado información sobre la velocidad de la reacción, aun cuando se sabe que la
segunda reacción tiene lugar a temperaturas tan bajas como 212 °F mientras que la reacción 1 tiene
una velocidad insignificante a menos de 392 °F. A medida que la temperatura aumenta a más de 392
°F la velocidad de la reacción de la ecuación química 1 aumenta rápidamente, hasta que a 932 °F es
varias veces más veloz que la de la reacción 2. A partir de esta información general y de algunos
Por: Ing. Edy Gonzalo Aspi Quispe Incotegas

datos publicados sobre el efecto del tiempo de residencia y la temperatura sobre la conversión, se
han propuesto las siguientes ecuaciones de velocidad.
( )
( ) ( )

( )
( ) ( )

Donde y están en libra moles de convertido por hora por pie cúbico, T está en grados
Rankine y la presión parcial p en atmósferas.
Determine la conversión a cloruro de alilo en función de la longitud del tubo. Se puede suponer que
la presión es constante e igual a 29.4 , considerando los siguientes:
a) Conversión a cloruro de alilo en función de la longitud del tubo, considerando solo la 1ra
reacción.
b) Conversión a cloruro de alilo en función de la longitud del tubo, considerando solo la 1ra y
la 2da reacción.

Solución.
Datos:

  ( )
  ( )
  ( ) ( )

 ( ) ( )

 ( ) ( )
 [( )( )( )]  ( ) ( )
 ( ) ( )

Formulación de las ecuaciones.

Para la solución de la misma, en este caso nos pide estimar la conversión de cloro y la temperatura
del sistema a lo largo del reactor PFR de 2 pulgadas, para lo cual, las ecuaciones de balance de
masa y energía se debe formular, en función de la longitud del reactor, además es importante ver
que según Smith realiza el estudio hasta una longitud de 20 ft.
 Balance de masa: conversión en función de longitud.
( )
( )

( )
( )

Llevando las ecuaciones anteriores y las velocidades de la reacción en función de la longitud del
reactor, sabiendo que z es la longitud del reactor y d su diámetro, tenemos:

Reemplazando la ecuación C en A y B, tenemos:

Por: Ing. Edy Gonzalo Aspi Quispe Incotegas

( )
( ) ( )

( )
( ) ( )

 Balance de energía: temperatura en función de longitud.

∑ ∫ ( )
̇ ∑ ∫ ( )
̇
| | | |

∑ ∫ ( ) ( )

También sabemos que:

( ) ( )

Reemplazando las ecuaciones G, C y H en F y reordenando, T en función de z y sabiendo que Cp

son constantes, tenemos:

∑ ∫ ( ) ( ) ( ) ( )

∑ ∫ [ ( ) ( ) ( ) ]

Reordenando y en términos de las dos reacciones, tenemos:

( ( ) ( )) ( )
∑

 Tratamiento de la velocidad de reacción en términos de conversión:

Sabiendo que entran al reactor 4 moles de propileno por 1 mol de cloro, la fracción molar de cada
componente es: 4 mol de propileno = 0,8 y 1 mol de cloro = 0,2, por lo tanto la fracción molar, en
términos de conversión que se forman, tanto para la reacción 1 y 2, tenemos:
( ) ( )
( ) ( )
( )
( )
( )
( )
Calculo de la presión parcial de los componentes, para la velocidad de reacción ( ), asumiendo el
comportamiento del gas ideal a una presión total del sistema de .
( ) ( ) ( ) ( )
[ ] [ ]
( ) ( )
( ) ( ) ( ) ( )
[ ] [ ]
( ) ( )
Reemplazando las ecuaciones J y K en las ecuaciones de velocidad tenemos:
Por: Ing. Edy Gonzalo Aspi Quispe Incotegas

( )
( )( )

( ) ( ) ( ) ( ) ( )
[ ] [ ] ( )( )
( ) ( )

( ) [ ( ) ( ) ][ ( ) ( ) ]
{ } ( )( )
[ ( ) ]
( )
( )( )

( ) ( ) ( ) ( ) ( )
[ ] [ ] ( )( )
( ) ( )

( ) [ ( ) ( ) ][ ( ) ( ) ]
{ } ( )( )
[ ( ) ]

Procedimiento de solución para el inciso a):

Para este caso solo consideramos que en el reactor ocurre la reacción 1 ( ),
para lo cual, consideraremos solo las ecuaciones D y reordenando las ecuaciones I y L, témenos:
( )

( ) ( )
[ ( ) ] [ ( ) ]
( )
( )
( )
( )

( ) [ ( ) ][ ( ) ] ( )
{ } { ( ) [( ) ] } ( )

Para la solución mediante Runge Kutta, hasta una longitud de 20 ft, tenemos las siguientes
ecuaciones, seleccionando una . Los pasos a seguir pueden ser las siguientes:

Para conversión Para temperatura

[ ( )] [ ( )]

[ ( )
] [ ( )
]

[ ( )
] [ ( )
]

[ ( )] [ ( )]

( ) ( )

Paso 1. Cálculo del incremente de conversión y temperatura hasta 1 ft por Runge Kutta.
 Sub-paso 1, cálculo para , por Runge Kutta, con base al reactivo limitante (cloro):
Por: Ing. Edy Gonzalo Aspi Quispe Incotegas

 Calculo de incremento de conversión en el reactor bajo, diferencial de longitud:

[ ( )] [ ( )]

( )
{ [ ] } ( )

( ) ( )

[ ( )]

 Calculo de incremento de temperatura en el reactor, bajo diferencial de longitud:

( )
( ) ( )
( ) ( )
[ ( )]

 Sub-paso 2, cálculo para , por Runge Kutta, con base al reactivo limitante (cloro):
 Calculo de incremento de conversión en el reactor bajo, diferencial de longitud:

[ ( )
] [ ( )
]

[ ( )]

( )
{ [ ] } ( )

( ) ( )

[ ( )
]

 Calculo de incremento de temperatura en el reactor, bajo diferencial de longitud:

( )
( ) ( )
( ) ( )
[ ( )]

 Sub-paso 3, cálculo para , por Runge Kutta, con base al reactivo limitante (cloro):
 Calculo de incremento de conversión en el reactor bajo, diferencial de longitud:

[ ] [ ] [ ( )]
( ) ( )

( )
{ [ ] } ( )

( ) ( )

[ ( )]

 Calculo de incremento de temperatura en el reactor, bajo diferencial de longitud:

Por: Ing. Edy Gonzalo Aspi Quispe Incotegas

( )
( ) ( )
( ) ( )
[ ( )]

 Sub-paso 4, cálculo para , por Runge Kutta, con base al reactivo limitante (cloro):
 Calculo de incremento de conversión en el reactor bajo, diferencial de longitud:
[ ( )] [ ( )]
[ ( )]

( )
{ [ ] } ( )

( ) ( )

[ ( )]

 Calculo de incremento de temperatura en el reactor, bajo diferencial de longitud:

( )
( ) ( )
( ) ( )
[ ( )]

 Sub-paso 5, cálculo de conversión para y T, a una longitud de 1 ft, por Runge Kutta:

( ) ( )

( ) ( )

Resumen de resultados para una longitud de 1 pie.

Para 1 ft
Subíndice T X r K m X1 T1
0 851,67 0 0,01361343 0,000349411 0,76373899
1 852,02753 0,00017471 0,01369696 0,000351555 0,71505902
2 852,02753 0,00017578 0,01369058 0,000351391 0,7181168
3 852,388117 0,0001757 0,01378403 0,000353789 0,672773 0,000351515 852,006386

Paso 2. Cálculo del incremente de conversión y temperatura hasta 2 ft por Runge Kutta.
Los sub-pasos descritos del 1 al 5, se repiten, ahora hasta 2 ft, para lo cual las condiciones iniciales
son:

Resumen de resultados.
Subíndice T X r K m X1 T1
0 852,387144 0,00035152 0,01376862 0,000353394 0,67222732
1 852,701786 0,00052821 0,0138408 0,000355246 0,62928369
2 852,701786 0,00052914 0,01383512 0,000355101 0,63197974
3 853,019124 0,00052907 0,01391807 0,00035723 0,59209245 0,000706735 853,018285
Por: Ing. Edy Gonzalo Aspi Quispe Incotegas

Paso 3, Cálculo del incremente de conversión y temperatura hasta 3 ft por Runge Kutta.
Para una longitud de 3 ft, para lo cual las condiciones iniciales son:

Resumen de resultados.
Subíndice T X r K m X1 T1
0 853,018285 0,00070673 0,01390236 0,000356826 0,59149064
1 853,295084 0,00088515 0,01396433 0,000358417 0,55359864
2 853,295084 0,00088594 0,01395928 0,000358287 0,55597671
3 853,574262 0,00088588 0,01403279 0,000360174 0,52089633 0,001065136 853,573541

Paso 4, Cálculo del incremente de conversión y temperatura hasta 4 ft por Runge Kutta.
Para una longitud de 4 ft, para lo cual las condiciones iniciales son:

Resumen de resultados.
Subíndice T X r K m X1 T1
0 853,573541 0,00106514 0,01401682 0,000359764 0,5202446
1 853,816943 0,00124502 0,01406962 0,000361119 0,48680367
2 853,816943 0,0012457 0,01406514 0,000361004 0,48890199
3 854,062443 0,00124564 0,01413017 0,000362673 0,45805826 0,00142625 854,061827

Paso 5, Cálculo del incremente de conversión y temperatura hasta 8 ft por Runge Kutta.
Para una longitud de 8 ft, para lo cual las condiciones iniciales son:

Resumen de resultados.
Subíndice T X r K m X1 T1
0 855,198905 0,00252188 0,014319 0,00036752 0,3095828
1 855,343541 0,00270564 0,01434391 0,000368159 0,28927266
2 855,343541 0,00270596 0,01434115 0,000368089 0,29054779
3 855,489452 0,00270593 0,01438039 0,000369096 0,27220752 0,002890068 855,489143

Paso 6, Cálculo del incremente de conversión y temperatura hasta 12 ft por Runge Kutta.
Para una longitud de 12 ft, para lo cual las condiciones iniciales son:

Resumen de resultados.
Subíndice T X r K m X1 T1
0 856,160676 0,00399981 0,01444309 0,000370705 0,18156734
1 856,245289 0,00418516 0,01445057 0,000370897 0,16922621
2 856,245289 0,00418526 0,0144489 0,000370854 0,17000273
3 856,330679 0,00418523 0,01447197 0,000371446 0,15925141 0,00437075 856,330555

Paso 7, Cálculo del incremente de conversión y temperatura hasta 16 ft por Runge Kutta.
Para una longitud de 16 ft, para lo cual las condiciones iniciales son:
Por: Ing. Edy Gonzalo Aspi Quispe Incotegas

Resumen de resultados.
Subíndice T X r K m X1 T1
0 856,719147 0,00548436 0,01445698 0,000371062 0,10377594
1 856,767289 0,00566989 0,01445381 0,00037098 0,09628468
2 856,767289 0,00566985 0,01445281 0,000370954 0,09675768
3 856,815904 0,00566984 0,01446593 0,000371291 0,0906215 0,005855396 856,815894

Paso 8, Cálculo del incremente de conversión y temperatura hasta 20 ft por Runge Kutta.
Para una longitud de 20 ft, para lo cual las condiciones iniciales son:

Resumen de resultados.
Subíndice T X r K m X1 T1
0 857,032661 0,00696662 0,01440407 0,000369704 0,05656509
1 857,058675 0,00715147 0,01439451 0,000369458 0,05202824
2 857,058675 0,00715135 0,01439391 0,000369443 0,05231604
3 857,084977 0,00715134 0,01440098 0,000369624 0,04898658 0,007336138 857,085034

Resumen total de los resultados

z ft X T ºR
0 0 851,67
1 0,00035152 852,387144
2 0,00070673 853,018285
3 0,00106514 853,573541
4 0,00142625 854,061827
5 0,00178966 854,491
6 0,00215498 854,867988
7 0,00252188 855,198905
8 0,00289007 855,489143
9 0,00325927 855,743464
10 0,00362925 855,966072
11 0,00399981 856,160676
12 0,00437075 856,330555
13 0,00474192 856,478605
14 0,00511316 856,607382
15 0,00548436 856,719147
16 0,0058554 856,815894
17 0,00622617 856,899387
18 0,00659661 856,971184
19 0,00696662 857,032661
20 0,00733614 857,085034
Por: Ing. Edy Gonzalo Aspi Quispe Incotegas

Procedimiento de solución para el inciso b):