Producción de Estireno

Diseño de Plantas II

Sección 1
Reactor R-201

Figura 1.1 Diagrama del proceso

Tabla 1. Parámetros obtenidos del pfd

Corriente 10 11
Descripció Salida del horno F- Salida del Reactor R-
Parámetros n 201 201
Unidade
Flujo s ton/h 34,456 34,456

Estado Físico S/L/G/V V V

Densidad kg/m^3 0,645 0,321

Temperatura ºC 730,000 700,000

Presión barg 0,01 0,01

Etilbenceno 0,442 0,005

Estireno 0,027 0,279

Hidrogeno --- 0,005

Metano --- 184 ppm

Etileno --- 0,055

Fracción Másica

Benceno 1 ppb 0,132

Tolueno 1ppm 973 ppm

Agua 0,523 0,523

Hidróxido de Sodio --- ---

Ácido Sulfúrico --- ---

Aire --- ---

Fuel Oil --- ---

Entalpia kcal/kg -1272,966 -1160,143

 Producción de Estireno
Diseño de Plantas II
En el proceso de producción del estireno se emplea como materia prima etilbenceno
que originalmente se encuentra alrededor de 35˚C como líquido, la cual se mezcla con una
corriente de vapor de agua (relación vapor/etilbenceno =8 mol/mol), la presencia de esta
gran cantidad de vapor de agua sobrecalentado se encarga de suministrar calor a la reacción
(reacción endotérmica), además que consigue disminuir la presión parcial de los reactantes,
evitando la formación de depósitos carbonosos y ayudando al avance de la reacción. Esta
corriente se debe acondicionar a una temperatura de 730 °C para ser alimentada a un
reactor multitubular en fase gas, que trabaja con un catalizador sólido de óxido de hierro
(Fe2O3) en forma de hematita promovidos por los óxidos de potasio, cerio y cromo. La
corriente a la salida del reactor contiene alta composición de estireno, así como también de
etilbenceno sin reaccionar, agua, tolueno y benceno que debe ser retirados en procesos
posteriores para la obtención del producto con las especificaciones requeridas.

Descripción del Reactor:

El reactor catalítico es el equipo donde una reacción química catalítica tiene lugar
de manera controlada, para nuestro proceso se utiliza un reactor catalítico heterogéneo
donde podemos encontrar que los reactantes (etilbenceno) se encuentra en fase gaseosa y el
catalizador en fase sólida (óxido de hierro). Actualmente, el 90% de todos los procesos
catalíticos se basa en reacciones catalizadas por un sólido. Ello se explica por la fácil
separación del sólido del medio de reacción. Ventaja que se suma a las que normalmente se
tienen cuando se utiliza catalizador como:

1) El catalizador aumenta la velocidad de reacción respecto a la reacción sin

catalizador, por lo que puede llevarse a cabo a menor presión y temperatura. En
consecuencia, de forma más segura y fácil de controlar.

2) Aumenta la selectividad al acelerar fundamentalmente una de las reacciones

del sistema entre las varias que pueden tener lugar.

Además, de la característica de ser un reactor catalítico heterogéneo, se usa un

reactor de lecho fijo, el cual consiste en uno o más tubos empacados con partícula de
catalizador que se opera en posición vertical. Las partículas catalíticas pueden variar de
forma: granulares, gránulos compactados, cilíndricas, esféricas, etc. En algunos casos,
especialmente con catalizadores metálicos como el platino, no se emplean partículas de
 Producción de Estireno
Diseño de Plantas II
metal, sino que éste se presenta en forma de mallas de alambre. El lecho catalizador
consiste de un conjunto de capas de este material. Debido a la necesidad de mantener el
calor en todo el reactor se sobrecalienta toda la corriente antes de entrar al reactor como se
ilustra en la Figura 2

Figura 2. Reactor catalítico con precalentamiento [ CITATION Sta90 \l 8202 ]

En la actualidad, para la producción del estireno se utiliza un reactor de lecho fijo

multibular permitiendo el contacto del etilbenceno en fase gas, con el óxido de
hierro(catalizador) en fase sólida. Pese a que la reacción tiene un alto rendimiento y
selectividad, se producen simultáneamente dos reacciones secundarias, las tres reacciones
mencionadas son endotérmicas, por lo que es de suma importancia controlar la temperatura
para no afectar la velocidad de reacción. Debido a la necesidad de introducir calor (reacción
endotérmica), no es posible usar un solo tubo de gran diámetro empacado con el
catalizador. En su lugar, el reactor deberá construirse con cierto número de tubos
engarzados en un solo cuerpo como se ilustra en la Figura 3. El intercambio energético con
los alrededores se obtiene circulando, o algunas veces hirviendo un fluido en el espacio que
hay entre los tubos. Si el efecto calorífico no es considerable, (como en nuestro reactor,
debido al flujo de vapor de agua sobrecalentado que entra con el reactante), no es necesario
este parámetro de diseño.

Figura 3. Diagrama de un reactor multibular de lecho fijo

 Producción de Estireno
Diseño de Plantas II
Parámetros de Diseño que definen el Reactor R-301
Aunque la intención es el diseño detallado del reactor, esto llevaría una serie de
pasos complejos y de experiencia industrial que se adquiere con estudios experimentales y
con conocimientos más detallado de diseño de rectores de los cuales no se cuenta a este
nivel, por ende, lo expresado a continuación es para indicar una base de diseño adecuada,
de manera de tener un estimado de cuáles deberían ser los parámetros adecuados para la
construcción del reactor. Mediante la simulación en Matlab se logró modelar el
comportamiento en estado estacionario del proceso químico estudiado, determinando
presiones, temperaturas y velocidades de flujo, además de las dimensiones características.

Tabla 2. Datos de la Industria Química de Reactores [ CITATION Sta90 \l 8202 ]

Tipo de Fase del Condiciones de operación Tiempo de residencia y

Producto Catalizador
reactor Reactor P[atm] y T [°C] velocidad espacial
Estireno a 1 atm 0,3 s
Reactor Óxidos de
partir de gas
multibular metales 600-650 °C 700 por hora
etilbenceno

Reacciones

En el reactor R-301, ocurren las siguientes reacciones:

Reacción Principal

La reacción de deshidrogenación de etilbenceno es reversible, ocurre con bajo

esfuerzo, además de ser altamente endotérmica y la conversión se limita al equilibrio
existente, teniendo un ∆H de reacción de 124,9 Kj/mol a una temperatura de 600 °C.
[ CITATION Ull98 \l 8202 ]

Reacciones Segundarias
 Producción de Estireno
Diseño de Plantas II
Además de la reacción principal, hay otras reacciones que no están en el equilibrio y
que están presentantes en el proceso, estas reacciones están limitadas por las condiciones
con las que opera el reactor, estas reacciones son:

Cinética Química

Para determinar la rapidez de la reacción, es necesario conocer diversos factores,

como la temperatura de la reacción, tipo de reacción, concentración de reactivos, presencia
de catalizador.

 Temperatura de Reacción

La temperatura se debe controlar, ya que, si se comienza a disminuir, se desfavorece

la reacción de deshidrogenación, aumentando las reacciones segundarias de formación de
benceno, tolueno, etileno, hidrogeno y metano, si se utiliza una temperatura muy alta se
corre el riesgo de la formación de depósitos carbonosos (coque) lo cual disminuiría la
efectividad del reactor y su tiempo útil, ocasionando paradas continuas del proceso. Para
este fin el reactor se trabajará isotérmicamente a 730°C, utilizando un reactor multitubular,
como si fuera un intercambiador de calor de tubo y coraza, en donde los reactantes pasaran
por los tubos y por la coraza se tendrá un fluido que mantenga la temperatura, es decir que
suministre calor al sistema debido a que es una reacción endotérmica.

 Presión
 Producción de Estireno
Diseño de Plantas II
En los reactores de lecho fijo, la caída de presión es muy pequeña, hasta el punto de
despreciar sus efectos (Smith, 1986). En nuestro caso trabajaremos con la ecuación de
Ergun para simular la caída de presión que se da por el roce de los reactantes con el lecho
fijo de catalizador, además si contemplamos que se utiliza una serie de tubos, estos
disminuyen la caída por los mismos. Lo que hace que nuestro proceso logre trabajar a
presiones bajas, favoreciendo así el avance de la reacción y selectividad de la misma a la
formación de estireno.

 Catalizador

Los óxidos inorgánicos con promotores son los catalizadores más empleados en la
industria para la producción de estireno-monómero, particularmente los óxidos de hierro en
forma de hematita (Fe2O3) promovidos por los óxidos de potasio, cerio y cromo. Si bien
estos catalizadores presentan alta actividad, selectividad y bajo costo por ser abundantes, su
desactivación rápida, su área superficial especifica baja, su alta toxicidad y su mayor
polución ambiental debido a la presencia de cromo constituyen notorias desventajas.
[ CITATION Sun06 \l 8202 ]

Figura 4. Forma del Óxido de hierro

La composición más común de este catalizador es 90% Fe203, 4 % Cr203 y 6 %
K2C03. A la temperatura de reacción el carbonato potásico es considerablemente volátil y
esto hace que sea lentamente perdido. Por esto es que el carbonato debe ser repuesto por
adición a la alimentación. Comercialmente se suelen utilizar partículas de catalizador
extraídas de 4-6 mm de diámetro. Para mejorar la selectividad del catalizador se utilizan
tamaños más pequeños, trayendo como consecuencia un aumento en la presión. El tiempo
de vida del catalizador es de 1 a 2 años.
 Producción de Estireno
Diseño de Plantas II
Las propiedades físicas del catalizador utilizadas para nuestro dimensionamiento se
enumeran en Tabla 3. Las mismas se miden en el laboratorio o viene dado por el fabricante
del catalizador.

Tabla 3. Datos industriales para la geometría y las propiedades del reactor[ CITATION Eln93 \l 8202 ]

Densidad global del

Porosidad del Catalizador Diámetro de Partícula
catalizador
kg Cat
2146.27 0,35 0.0047 m
m3

Consideraciones de Diseño para el dimensionamiento del Reactor R-301:

Como se señaló anteriormente, el reactor R-301 consiste en un intercambiador con

una coraza, que contiene en su interior una cantidad de tubos. Para estimar el diámetro de la
coraza se debe conocer entonces los parámetros de diseño de los tubos empleados y el
número de tubos, de manera que se pueda calcular el área ocupada por el set de tubos,
como se muestra a continuación:

π∗d 2
A tubo = Ec. 1
4
Donde Atubo = área de los tubos y d= diámetro de los tubos

En principio, para obtener una base de cálculo se define el diámetro del tubo como
el más común utilizado para intercambiadores de tubo y coraza de 4 plg de (0,1016 m) con
un espesor de 0,049 plg. En el cual el diámetro interno de los tubos es 3,951 plg (0,1003
m).

π∗(0,1003)2
A tubo = =0.00791 m2
4
Para calcular el número de tubos que necesito en mi proceso divido el flujo másico
entre el promedio de la densidad de entrada y salida de mi proceso, con esto obtengo el
flujo volumétrico entrante a mi reactor, si a este valor le dividimos la velocidad del fluido
por el área efectiva de cada tubo, obtendremos un estimado de los tubos que necesitamos.
 Producción de Estireno
Diseño de Plantas II
Fm
ϑ0 ρ( p , T )
N ° tubos = =
v∗A tubo∗θ v∗Atubo∗θ

Donde Fm = flujo másico de entrada, ρ=densidad del fluido, v = velocidad del

fluido, y θ es el área hueca que deja el catalizador.

Todos los datos de propiedades se tomaron de la simulación de Aspen 11.0 y del pfd
realizado para la primera entrega.

Para el calcula de la velocidad máxima permitidita para los gases en los tubos, se
tiene la ecuación:

1800 1800 ft
v= = =65,51
√ Presion∗PM s
√ 26,25951 psi∗28,75
lb
lbmol

Entonces tenemos que el número de tubos utilizados es:

Kg
34456
h
Kg
0,483 3
ϑ0 m
N ° tubos = = =358,5 ≈ 359 tubos
v∗A tubo∗θ ft

( )
∗0,3048 m
s
∗3600 s
1 ft
65,51 ∗0.00791 m2∗0,35
1h

Teniendo el número de tubos y el área de un tubo podemos calcular el área de flujo

del intercambiador y con esta el diámetro de flujo para introducirlo en los cálculos del
simulador.

A tubos=N ° tubos∗A tubo=359∗0.00791 m2=2,83969 m2

d tubos=
√ 4∗Atubos
π √
=
4∗2,83969 m2
π
=1,9015 m

Seguidamente se puede calcular entonces el diámetro de la coraza del reactor,

considerando la siguiente expresión, que corresponde al cálculo del número de tubos de un
intercambiador de uno o dos pasos:
 Producción de Estireno
Diseño de Plantas II

N ° tubos =
0,875∗( CTP
CL )
∗Dc 2

Pr 2 ¿ do2

Donde:

CTP = constante que relaciona el espaciamiento entre los tubos y la coraza con los
pasos a través de los tubos; cuyo valor es de 0,93 para un solo paso.

CL= constante de configuración y arreglo; el cual vale 0,87 para arreglos

triangulares, considerando que la configuración triangular es menos costosa y más eficiente
para la transferencia de calor.

PR = es la relación entre el espaciamiento entre tubo (pich) y el diámetro del tubo;

para tubos de 4 plg el pich corresponde a un valor (1,25*4) plg = 5 plg. Entonces PR = 5/4
= 1,25

√ √
CL
Dc=
Pr2 ¿ do2∗N ° tubos∗ ( CTP )= ( 0,87
1,252 ¿ 0,1016 m2∗359∗
0,93 )
=2,488 m
0,875 0,875

Longitud y Volumen del Reactor R-301:

Para un reactor de lecho se tiene la siguiente relación [ CITATION Lev04 \l 8202 ]:

L
esta entre 2 y 5
D

Asumiendo una relación de 3.5, como criterio de diseño podemos calcular la

longitud de los tubos utilizados (este valor se corroborará con el avance de la reacción y si
la misma se favorece dejando por más tiempo los tubos).

L=3.5∗1.9015 m=6.65525 m=21,8348 pie ≈ 22 pie=6.7056 m

Asumiendo que el reactor es cilíndrico, se tiene que su volumen es:

2 3
V = A tubos∗L=2,83969 m ∗6.7056 m=19.0418 m

Peso del catalizador

Considerando que el peso del catalizador sólido debe estimarse utilizando la
ecuación descrita por [ CITATION Fog011 \l 8202 ] tenemos que:
 Producción de Estireno
Diseño de Plantas II
W =L∗Atubos∗ρc∗(1−θ)

Donde L es la longitud de los tubos con lecho (m), Atubos es el área catalizada (m2),
ρc la densidad del catalizador (Kg/m3) y Φ es la porosidad del catalizador. Todas estas
propiedades y valores antes calculados, entonces podemos estimar el peso del catalizador.

2 kg Cat
W =2,83969m ∗6.7056 m∗2146.27 ∗( 1−0.35 )=26564,5 Kgde Cat
m3

Tiempo de residencia
El tiempo de residencia (τ) se define como la relación entre el volumen de reactor y
el caudal volumétrico (υ):
3
19.0418 m ∗0.35
∗3600 s
Kg
34456
h
Kg
0,483 3
V Vreactor∗θ m
τ= = = =0.33633 segundos
υ0 υ0 1h

Como se puede observar el tiempo de residencia es baja, debido a la cinética de las

reacciones que son bastante espontaneas haciendo que su velocidad de reacción sea elevada
y en consecuencia un tiempo de residencia bajo, por lo que es ventajoso ya que el costo del
reactor es bajo en comparación a otros procesos de cinéticas mucho más exigentes, además
si lo comparamos con el tiempo de residencia utilizado en la industria (ver Tabla 2) la
discrepancia es mínima por lo que se asume que se está trabajando con parámetros
utilizados en la industria.

Fundamentos Teóricos para el Cálculo de los Flujos

Recapitulando, las reacciones involucradas son:

Reacción 1. Reacción principal reversible

C6H5-C2H5 → C6H5-CH=CH2 + H2

A  B + C

Reacciones Secundarias.

Reacción Secundaria 1.
 Producción de Estireno
Diseño de Plantas II
C6H5-C2H5 → C6H6 + CH2=CH2

A  D + E

Reacción Secundaria 2.

C6H5-C2H5 + H2 → C6H5-CH3 + CH4

A + C  F + G

Para el diseño de reactores catalíticos las ecuaciones se deducen de manera análoga

a los equipos tubulares. No obstante, la principal diferencia entre reactores en los que
interviene reacciones homogéneas, y aquellos en los que ocurren reacciones heterogéneas
gas-sólido; es que en este último la velocidad de reacción se basa en la masa del catalizador
sólido (W) y no con respecto al volumen del reactor (V). Para ello se sustituye la
coordenada de la masa del catalizador y se supone en que el PBR no existen gradientes
radiales de concentración, temperatura ni velocidad de reacción. De manera que al plantear
un balance molar en un diferencial ΔW se obtiene:

Moles Entran−MolesSalen + MolesGeneran =Moles Acumulan

Fi ( Wo )−Fi ( Wo+ ∆ W ) +ri∗∆ W =Moles Acumulan

Dicho esto, la ecuación de diseño para un reactor flujo pistón catalizado será:

dX −ri
=
dW F ¿

Finalmente, para definir todas las variables de la ecuación de diseño, se expone a

continuación la cinética para la reacción principal y la secundaria con mayor avance, de la
siguiente manera:

Ley de la velocidad

R 1=8.314 ; constarnte de los gasesideales [J /mol∗k ]

k 1=(2.2215 e 16/3600)∗exp(−272230/(T∗R 1)) ;(Kmol / .̄ m3. s)

¯¿
Ke=exp(16.12−(15.350/T )); ¿
 Producción de Estireno
Diseño de Plantas II
3.
ra=( k 1.∗(Pa−( Pb∗Pc)/Ke)); (kmol / m s)

k 2=(2.4217 e 20 /3600)∗exp(−352790/(T∗R 1)) ;(Kmol / .̄ m3. s )

3.
rb=k 2.∗Pa ;(kmol/m s )

k 3=(3.8224e17 /3600)∗exp(−313060/(T∗R 1));( Kmol / .̄ m3.s )

3. s
rc=(k 3.∗Pa.∗Pc); (kmol /m )

Donde:

T = temperatura del flujo gaseoso (K)

Pi = presión parcial del compuesto i (bar)

Las ecuaciones de las constantes de reacción tienen la forma de Arrhenius y sus

valores fueron tomas de la literatura [ CITATION Eln93 \l 8202 ] y la constante de
equilibrio para la primera reacción fue tomada del [ CITATION Ull98 \l 8202 ]. Los valores

j
para la energía de activación están en unidades de y el factor pre-exponencial está
mol

Kmol
en unidades de .
¯¿ m3∗h

Empleando la estequiometria se estima la presión parcial de cada componente,

mediante la siguiente ecuación:

Fi
∗¿
Ft
Pi= ∗Pt
T

Sabiendo que tenemos un reactor isotérmico y para facilitar los cálculos podemos
tomar las presiones parciales de los compuestos involucrados en las reacciones del reactor
R-301 como:

Cto=Po/(R*To); % concentración inicial [Kmol/m^3]

Ci=Cto*(Fi/Ft) * (To/T); % cambio de la concentración del compuesto i [kmol/m3]

Pi=Ci*R*T; % variación de la presión del compuesto i [bar]

 Producción de Estireno
Diseño de Plantas II
Por consiguiente, la variación de la presión para cada componente involucrado en
función del peso del catalizador se plantea de la expresión para la caída de presión que
viene dada por Ergun [ CITATION Fog011 \l 8202 ]:

−α
∗T
2∗y
∗F T
dP TO
= ∗Po
dW F¿

Alfa queda definida por la siguiente ecuación:

2∗β
∗1
Atubos∗( 1−θ )∗ρc
α=
Po

β=
( G∗(1−θ)
ρ o∗gc∗Dp∗θ 3
∗
)(
( ( 150 )∗( 1−θ )∗μ )
Dp
∗1.75∗G )
Donde:
2
G = Flujo Másico por unidad de Área Kg/(m ∗s)

ρo = Densidad del Fluido a la entrada del reactor Kg/m

3

μ = Viscosidad promedio de la mezcla reaccionante P(1 P=1 Kg /(m∗s))

gc = Factor de corrección Kg∗m/ s2∗N

Dp = diámetro de partículas del catalizador en el lecho fijo m

θ = porosidad

Para estimar el perfil de temperatura se plantea un balance de energía obteniendo la

ecuación de diseño siguiente:

Ua
∗(Ta−T ) + ∑ (−ri ) + [ −∆ H rxij (T ) ]
dT ρo
=
dW ∑ Fj∗Cpj
 Producción de Estireno
Diseño de Plantas II
Empleando el simulador de Matlab se realizó el código que diseña un reactor
tubular debido a que todos los componentes están en fase gaseosa, el cual funciona
isotérmicamente, trabajando con una temperatura de 730°C y a 1 atm de presión.

El código se encuentra en los anexos, en donde podemos graficar el comportamiento

del flujo másico en función del volumen del reactor recorrido, así como la caída de presión
por el lecho y la conversión del etilbenceno a estireno del mismo.

Para el comportamiento de la temperatura en el lecho no fue necesaria graficarla ni

calcularla, ya que, el flujo de vapor que entra acompañado del etilbenceno es lo
suficientemente grande para suministrar el calor que necesita la reacción, y así mantenerla
isotérmicamente (con una misma temperatura en todo el reactor), es decir la temperatura no
varía significativamente entre la entrada y salida del rector R-301.

Figura 5. Flujo másico de los compuestos en función del volumen del reactor
Producción de Estireno
Diseño de Plantas II

Figura 6. Comportamiento de la caída de presión a medida que el fluido pasa por el reactor

Tabla 4. Dimensiones Características y Especificaciones del Reactor

Código del Equipo R-301

Diámetro de los tubos [plg] 4

N° de tubos utilizados 359

Velocidad de los reactivos por los tubos [ pie/s] 65.51

Tiempo de residencia [segundos] 0.33633

Área de Sección Trasversal de los tubos [m2] 2,83969
Longitud de los tubos [pie] 22
Arreglo de los tubos triangular
Volumen del reactor [m3] 19.0418
Diámetro de la coraza [m] 2,488
Material del catalizador Óxido de hierro
Peso del catalizador [Kg] 26564,5
 Producción de Estireno
Diseño de Plantas II

Diagrama del equipo

Bibliografía

Elnashaie, S., & Elshishini, S. (1993). MODELLING, SIMULATION AND

OPTIMIZATION OF INDUSTRIAL FIXED BED CATALYTIC REACTORS.
Amsterdam: Gordon and Breach Science Publishers S.A.

Fogler, H. .. (2001). Elemento de Ingeniería de las Reacciones Químicas. Mexico: Pearson.

Kern, D. (1999). Procesos de transferencia de calor. Mexico: CONTINENTAL, S.A.

Lee, S. (2006). Encyclopedia of Chemical Processing. Columbia, Missouri: Taylor &

Francis Group.

Levenspiel, O. (2004). Ingenieria de las Reacciones Quimicas. Mexico: Limusa-Wiley.

Ullmann’s Encyclopedia of Industrial Chemistry. (1985). Germany: VCH: Weinheim.

Walas, S. M. (1990). Chemical Process Equipment Selection and Design. Washington:

Library Cataloguing.
 Producción de Estireno
Diseño de Plantas II

Código principal para la simulación de la producción de estireno por

deshidrogenación catalítica del etilbenceno en un Reactor multitubular isotérmico con
catalizador de oído de hierro.
clear all
close all
Fao=Fto*yao; % flujo inicial molar del etilbenceno [kmol/s]
clc
Fbo=Fto*ybo; % flujo inicial molar del estireno [kmol/s]
format long
Fco=Fto*yco; % flujo inicial molar del hidrogeno [kmol/s]
global To Po Fto Fio Cto R FMT Rhoi Gc G Visc Ua Ar Ta Dp
Fdo=Fto*ydo; % flujo inicial molar del benceno [kmol/s]
rhoc phi
Feo=Fto*yeo; % flujo inicial molar del etileno [kmol/s]
%% Declaro las variables de operación
Ffo=Fto*yfo; % flujo inicial molar del tolueno [kmol/s]
To=720+273.15; %temperatura de operación [°K]
Fgo=Fto*ygo; % flujo inicial molar del metano [kmol/s]
T=To; %temperatura en todo el lecho del reactor [°K]
Fio=Fto*yio; % flujo inicial molar del Agua [kmol/s]
Po=1.01325; % presión inicial [bar]
R=0.083144721345; % constante de los gases ideales
[m^3*bar/Kmol*k]
%% Declaro las variables y las propiedades del gas a Cto=Po/(R*To);%concentración inicial [Kmol/m^3]
reaccionar Cao=yao*Cto; %concentración del reactivo
FMT= 34456/3600; %Flujo Másico total [=]Kg/s limitante(etilbenceno) [Kmol/m^3]
Rhoi=0.645; %Densidad del Fluido [=] Kg/m^3
Fto=1162.4668/3600; %flujo molar total que entra al reactor
%% Declaro las dimensiones del reactor
[kmol/s]
dtubo=0.1003554 ; %Diámetro de un tubo
yao=0.12350909;
Atubo=(pi*(dtubo)^2)/4; %área transversal por tubo[m^2]
ybo=0.00754332;
Nt=359; % número de tubos utilizados
yco=1.6748e-09;
Ar=Nt*Atubo; %área transversal del lecho[m^2]
ydo=0.00013935;
dr=(Ar*4/pi)^(1/2); %diámetro del reactor [m]
yeo=0.00794051;
L1=3.5*dr; %longitud del reactor [m]
yfo=5.3221e-07;
L= 6.7056; %longitud del reactor [m]
ygo=2.7032e-05;
Vf=Ar*L; %Volumen del Reactor [m^3]
yio=0.86084015;
V=[0 Vf];

PMA=106; %Peso Molecular del etilbenceno [=]Kg/Kmol

%% Declaro las especificaciones del catalizador
PMB=104; %Peso Molecular del estireno [=]Kg/Kmol
Dp=0.0047; %diámetro de partícula [m]
PMC=2; %Peso Molecular del hidrogeno [=]Kg/Kmol
rhoc=2146.27; %densidad del catalizador [Kg/m^3]
PMD=78; %Peso Molecular del benceno [=]Kg/Kmol
phi=0.35; % fracción hueca%
PME=28; %Peso Molecular del etileno [=]Kg/Kmol
w=L*Ar*(1-phi)*rhoc; %masa del catalizador [Kg]
PMF=92; %Peso Molecular del tolueno [=]Kg/Kmol
PMG=16; %Peso Molecular del metano [=]Kg/Kmol
PMI=18; %Peso Molecular del Agua %% Declaro otras variables a considerar
[=]Kg/Kmol tres=(Vf*phi*3600)/71337.4741201; % tiempo de residencia
[s]
 Producción de Estireno
Diseño de Plantas II
Gc=1; % [Kg*m/s^2*N] plot(W,P,'g','linewidth',2)
G=FMT/Ar; %Flujo Másico por unidad de Área [Kg/ grid on
(m^2*s)] xlabel('Volumen del reactor (m)')
Visc=0.341044e-4; %Viscosidad promedio de la mezcla ylabel('Caída de Presión (bar)')
reaccionante [=]P (1 P = 1 Kg/(m*s))
Ua= 1; % [Kcal/kg. °K]
X1=Fb./Fao; %Conversión de la reacción de deshidrogenación
Ta=1073.15; %temperatura de chaqueta [°K]
del etilbenceno

%% Calculo del reactor multitubular operando

figure(3)
isotérmicamente.
plot(W,X1,'b','linewidth',2)
grid on
Masa=[0 w]; % Masa de Catalizador utilizada xlabel('Volumen del reactor (m)')
Io=[Fao Fbo Fco Fdo Feo Ffo Fgo To 1]; ylabel('Conversion del etilbenceno a estireno')
legend('Reacción 1')
%RESOLUCION DEL SISTEMA DE ECUACIONES
DIFERENCIALES función utilizada para resolver las
% Resolviendo las Edo por el metodo del ode45 con
[W,x1]=ode15s('compuestoMultipleIso',V,Io); ecuaciones diferenciales:
function solucion1=compuestoMultipleIso (W, Io)
%% FLUJOS MOLARES POR TUBO
format long
global To Po Fto Fio Cto R Rhoi Gc G Visc Ua Ar Ta Dp rhoc
Fa=x1(:,1); %Flujo Molar de etilbenceno a lo largo del phi
reactor [=]Kmol/s
Fb=x1(:,2); %Flujo Molar de estireno a lo largo del reactor
%% valores iniciales de convergencia
[=]Kmol/s
Fc=x1(:,3); %Flujo Molar de hidrogeno a lo largo del
reactor [=]Kmol/s Fa=Io(1);
Fd=x1(:,4); %Flujo Molar de benceno a lo largo del reactor Fb=Io(2);
[=]Kmol/s Fc=Io(3);
Fe=x1(:,5); %Flujo Molar de etileno a lo largo del reactor Fd=Io(4);
[=]Kmol/s Fe=Io(5);
Ff=x1(:,6); %Flujo Molar de tolueno a lo largo del reactor Ff=Io(6);
[=]Kmol/s Fg=Io(7);
Fg=x1(:,7); %Flujo Molar de metano a lo largo del reactor T=Io(8);
[=]Kmol/s Y=Io(9);
Fi=Fio; % flujo molar del inerte [kmol/s]
Ft=Fa+Fb+Fc+Fd+Fe+Ff+Fg+Fio; % Flujo total
%% FLUJOS MASICOS TOTALES

FAMTot=Fa*PMA*3600*(1/1000); %Flujo Másico de

etilbenceno a lo largo del reactor [=]Ton/h %% relaciones estequiométricas
FBMTot=Fb*PMB*3600*(1/1000); %Flujo Másico de estireno
a lo largo del reactor [=]Ton/h
Ca=Cto*(Fa/Ft)*(To/T); % cambio de la
FCMTot=Fc*PMC*3600*(1/1000); %Flujo Másico de
concentración del etilbenceno [kmol/m3]
hidroheno a lo largo del reactor [=]Ton/h
Cb=Cto*(Fb/Ft)*(To/T); % cambio de la
FDMTot=Fd*PMD*3600*(1/1000); %Flujo Másico de
concentración del estireno [kmol/m3]
benceno a lo largo del reactor [=]Ton/h
Cc=Cto*(Fc/Ft)*(To/T); % cambio de la
FEMTot=Fe*PME*3600*(1/1000); %Flujo Másico de etileno
concentración del hidrogeno [kmol/m3]
a lo largo del reactor [=]Ton/h
Cd=Cto*(Fd/Ft)*(To/T); % cambio de la
FFMTot=Ff*PMF*3600*(1/1000); %Flujo Másico de tolueno
concentración del benceno [kmol/m3]
a lo largo del reactor [=]Ton/h
Ce=Cto*(Fe/Ft)*(To/T); % cambio de la
FGMTot=Fg*PMG*3600*(1/1000); %Flujo Másico de metano
concentracion del etileno [kmol/m3]
a lo largo del reactor [=]Ton/h
Cf=Cto*(Ff/Ft) *(To/T); % cambio de la
FWMTot=Fio*PMI*3600*(1/1000);%Flujo Másico de agua a
concentración del tolueno [kmol/m3]
lo largo del reactor [=]Ton/h
Cg=Cto*(Fg/Ft)*(To/T); % cambio de la
figure
concentración del metano [kmol/m3]
plot(W,FAMTot,'g',W,FBMTot,'y',W,FCMTot,'r',W,FDMTot,'c',
W,FEMTot,'m',W,FFMTot,'k',W,FGMTot,'b','linewidth',2)
grid on Pa=Ca*R*T; % variación de la presión del etilbenceno
legend('etilbenceno','estireno','hidrogeno','benceno','etileno','tol [bar]
ueno','metano') Pb=Cb*R*T; % variación de la presión del estireno
xlabel('Volumen del reactor (m)') [bar]
ylabel('Flujo Masico (t/h)') Pc=Cc*R*T; % variación de la presión del hidrogeno
[bar]
Pd=Cd*R*T; % variación de la presión del benceno
P=Po*x1(:,9); %Presión a lo largo del reactor [=]bar
[bar]
Pe=Ce*R*T; % variación de la presión del etileno [bar]
figure(2) Pf=Cf*R*T; % variación de la presión del tolueno [bar]
 Producción de Estireno
Diseño de Plantas II
Pg=Cg*R*T; % variación de la presión del metano
[bar]

%% Ley de la velocidad

R1=8.314; % constante de
los gases ideales [J/mol*k]
k1=(2.2215e16/3600)*exp(-272230/(T*R1)) ;
% (Kmol/bar.m^3.s)
Ke=exp(16.12-(15.350/T)) ; % (bar)
ra=(k1.*(Pa-(Pb*Pc)/Ke)); %
(Kmol/m^3.s)
k2=(2.4217e20/3600)*exp(-352790/(T*R1)) ;
% (Kmol/bar.m^3.s)
rb=k2.*Pa; % (Kmol/m^3.s)
k3=(3.8224e17/3600)*exp(-313060/(T*R1)) ;
% (Kmol/bar.m^3.s)
rc=(k3.*Pa.*Pc); %
(kmol/m^3.s)

%% balance Molar

%NETOS
rAn=-ra-rb-rc;
rBn=ra;
rCn=ra-rc;
rDn=rb;
rEn=rb;
rFn=rc;
rGn=rc;

%% Estimación de Caída de Presión a lo largo de Reactor

Betha=(G/(Rhoi*Gc*Dp))*((1-phi)/(phi^3))*(((150*(1-
phi)*Visc)/(Dp))+1.75*(G)); % Factor Beta [=]Kg/(s^2*m^2)
Alfa=(2*Betha/(Ar*rhoc*(1-phi)*Po*1*10^5));% Factor Alfa
[=]1/Kg

%% Procedimiento a resolver
dFa=rAn ; % cambio del flujo molar del etilbenceno
[kmol/s]
dFb=rBn ; % cambio del flujo molar del estireno
[kmol/s]
dFc=rCn ; % cambio del flujo molar del hidrogeno
[kmol/s]
dFd=rDn ; % cambio del flujo molar del benceno
[kmol/s]
dFe=rEn ; % cambio del flujo molar del etileno
[kmol/s]
dFf=rFn ; % cambio del flujo molar del tolueno
[kmol/s]
dFg=rGn ; % cambio del flujo molar del metano
[kmol/s]
dT=0*((Ua*(Ta-T))+(rA1)*(dHr1)+(rA2*(dHr2)))/
(Fa*CpCO+Fb*CpH2+Fc*CpM+Cpw*Fd+Fe*CpCO2+CpN2
*Fi); %variación de la temperatura del reactor con respecto a la
masa del catalizador
dYdW=(-Alfa/(2.*Y)).*(T./To).*(Ft./Fto);
solucion1=[dFa;dFb;dFc;dFd;dFe;dFf;dFg;dT;dYdW];
end