Universidad del Istmo campus Tehuantepec. Ingeniería Química, Diseño de reactores. Santo Domingo Tehuantepec Oax a 03 de julio de 2014. M.C. José Luis Sánchez López
Modelado Matemático de un reactor tipo PBR trifásico para un sistema de reacción de HDS del Dibenzotiofeno con un catalizador de Co-Mo/γ-Al 2O3. Carrasco Carballo Alan Universidad del Istmo Campus Tehuantepec Ciudad Universitaria S/N, Barrio Santa Cruz, 4a. Sección Sto. Domingo Tehuantepec, Oax., México C.P. 70760, tel: 9711206801, correo electrónico:
[email protected]
Resumen En este este trab trabaj ajoo se pres presen enta ta la obte obtenc nció iónn de un mode modelo lo mate matemá máti tico co para para un reac reacto torr tipo tipo PBR PBR de tres tres fase fasess en dond dondee se llev llevaa a cabo cabo una una reac reacci ción ón de hidr hidrod odes esul ulfu fura raci ción ón del del dibe dibenz nzot otio iofe feno no con con un cata catali liza zado dorr come comerc rcia iall de Co-M Co-Mo/ o/γγ-Al Al2O3, donde el modelo comprende el seguimiento de la concentración de los reactivos y productos para esta reacción y la puesta en marcha del mismo modelo con un sistema típico de dibenzotiofeno disuelto en tetralina y con una corriente de H 2 gase gaseos osoo puro puro,, para para lo cual cual se obtu obtuvi vier eron on los los pará paráme metr tros os físi físico cos, s, quím químic icos os y ciné cinéti tico coss nece necesa sari rios os para para el mode modela lado do medi median ante te dive divers rsas as refe refere renc ncia ias, s, obte obteni nien endo do así así perf perfil iles es de conc concen entr trac ació iónn en func funció iónn de la longitud del reactor y permitiendo analizar el efecto de los fenómenos de trasporte sobre la conversión del reactor. Palabras clave: PBR, hidrodesulfuración, hidrodesulfuración, catalizador, dibenzotiofeno, reactor, hidrogeno.
Abstract In this this pap paper we obta obtain in a mat mathema hemattical ical mode modell for for a PBR PBR type ype thre threee phas hase reac reacttor wher wheree a reac reacttion ion take takess plac placee hydr hydrod odes esul ulfu furi riza zati tion on of dibe dibenz nzot othi hiop ophe hene ne with with a comm commer erci cial al cata cataly lyst st Co-M Co-Mo/ o/γγ-Al Al2O 2O33 wher wheree the mode modell incl ncludes udes moni monito tori rinng occu occurs rs the conc concen entr trat atio ionn of the reac reacttant ants and and the the reac reacti tion on products and the implementation of a typical model system dibenzothiophene dibenzothiophene dissolved in tetralin and and a strea tream m of pure pure H2 gas for for which hich the the chem chemiical cal para parame metters ers and and obtai btainned neces ecesssary ary kinet inetiic for for mod modeli eling usin sing diff differ eren entt refe refere rennces, ces, obtai btainning ing con concent centra rati tion on pro profil files bas based on the the len length gth of the reac reacto torr and and allow llowin ingg anal analys ysis is of the the effe effect ct of the phen henomen omenaa of tran transp spor orttati ation of the the conve onvers rsio ionn reactor.
Keywords: PBR, hydrodesulfurization catalyst, dibenzothiophene, reactor, hydrogen.
1. Introducción
mismo lo que permite realizar análisis de un
El consumo de diésel y otros combustibles de uso automotriz representa la mayor demanda de productos derivados del petróleo para el transporte de personas, bienes y servicios en Europa. Las normas de emisión de vehículos diésel inician con el establecimiento de niveles máximos permisibles a los automóviles según la directriz 91/441/EEC. Ésta da origen a la norma Euro I para el período 1992–1995; posteriormente se emite la Directriz 94/12/EEC, que promulga la norma Euro II, vigente para los vehículos de inyección indirecta desde 1996, las cuales son muy similares a las aplicadas en Estados Unidos y otros países, por lo que se ponen como limites valores entre 50 a 20 ppm de contenido de azufre en el diésel y en las gasolinas homólogamente, inclusive en algunos países como Japón se manejan valores aceptables máximos de 10 ppm, por lo que constantemente se ha desarrollado investigación en métodos para eliminar este contaminante. La hidrodesulfuración es un proceso destinado a eliminar el azufre (que es una impureza contaminante) que se encuentra en el combustible al finalizar todos los tratamientos anteriores, tales como destilación fraccionada, destilación por presión reducida, reforming o cracking [8]. Una opción muy utilizada en la actualidad para la investigación de este proceso es el modelado matemático [5] del
sistema reactivo (reactor) sin la necesidad de llevar acabo experimentación directa en un
reactor lo cual involucra gastos, además de evitar riegos innecesarios en la operación de dicho reactor, además la simulación permite realizar pruebas en condiciones extremas y optimizar de forma más económica, claro esto en función de que tan robusto sea el modelo, el principal problema es la obtención de parámetros físicos, químicos, termodinámicos y cinéticos para dicho modelos, para lo cual se puede obtener de experimentación o por medio de correlaciones ya establecidas, o bien se pueden buscar sustancias patrón que sirvan para simular el comportamiento de una carga normal con la ventaja de que para muchos patrones ya existen parámetros reportados en la literatura que facilitan el modelado matemático del reactor teniendo en cuenta que el patrón debe ser adecuado y de ser posible también la carga real debe contenerlo en una composición considerable y en algunos casos se busca que sea el más difícil de aplicar el sistema de reacción que se analice.
otros derivados azufrados[1].
El dibenzotiofeno es el compuesto orgánico de azufre que consta de dos anillos de benceno fusionados a un anillo central de tiofeno. Es un sólido incoloro que es químicamente algo similar al antraceno. Este heterociclo tricíclico (de tres anillos o ciclos fusionados entre sí), y especialmente sus derivados sustituidos con grupos alquilo, se producen ampliamente en las fracciones más pesadas de petróleo. Este compuesto y sus alquil derivados son utilizados como compuestos modelo en estudios hidrodesulfuración debido a la dificultad de su reacción ya que en caso de encontrar las condiciones para la HDS de estos compuestos normalmente es más fácil realizar la HDS en sustancias de menor peso molecular como los tioles, mercaptanos, tiofeno, benzotiofeno, y
Fig. 1 Sistema del reactor a analizar
2. Fundamentos teóricos.
Para el sistema de estudio de un reactor del tipo PBR planteado es del tipo piloto para investigación en un arreglo como el que se muestra en la fig.1 y las condiciones de operación son las que comúnmente operan los reactores de este tipo en la realidad.
Para el sistema a estudiar se plantean las ecuaciones para cada especie en las diversas fases del sistema las cuales provienen de balances de materia para cada especie, incluyendo en estos la influencia de los fenómenos de transporte, en este caso particular los fenómenos de trasferencia de masa para el efecto del modelo difusión-reacción para el catalizador de estudio y así mismo los fenómenos de difusión entre las distintas fases del sistema. 2.1 Hipótesis y consideraciones. Para el análisis del sistema se realizó basado en las siguientes hipótesis: I.
Se trabajara en estado estacionario.
II. III. IV.
V.
El Proceso será isobárico. El proceso será isotérmico. Se despreciaran los efectos de dispersión angular debido a que el flujo se alimenta en una sola dirección y con una concentración angular homogénea. Se desprecian los efectos de dispersión radial, debido a que se conservara una relación L/D lo suficientemente grande para despreciar este efecto.
3.0 Algoritmo del método de solución
Para modelar este proceso se realizara de acuerdo a la siguiente serie de pasos: 1. Plantear los modelos matemáticos de las ecuaciones de balance de materia para las especies a seguir en las fases que se encuentren presente, aplicando las hipótesis que se vean involucradas en esta etapa. 2. Plantear las condiciones de operación
para el caso de estudio a simular. 3. Analizar y observar los modelos obtenidos, para así determinar los VI.
VII. VIII. IX.
X.
Las fracciones de cada fase se mantendrán constantes a lo largo del reactor El factor de eficiencia se mantendrá constante en todo el reactor No existe desactivación del catalizador. El gas de hidrogeno entra en un sistema de flujo pistón por lo que no hay dispersión en el gas. Los productos no se vuelven a adsorber en el catalizador sino se van al líquido.
4.
5.
6.
7.
parámetros que serán necesarios para la simulación del proceso. Realizar Ñel análisis del modelo difusión-reacción para las especies involucradas para obtener los factores de efectividad de cada especie. Obtener todos los parámetros que se necesiten para el modelo ya sea consultando la bibliografía o por el uso de correlaciones. Simular el proceso y obtener perfiles de concentración para las especies de interés en función de la longitud del reactor. Analizar la coherencia de los datos y ver los efectos presentes en el sistema.
4.0 Método de Solución
Para esa reacción se propone un mecanismo de reacción según la literatura se proponen las dos rutas de reacción cada una con su propia cinética como se muestra en la fig. 2.
En la fase liquida se obtiene el modelo: (1) 2 l
l
l
∂C A l ∂ C A ls l u εl ∂∂C t A = εl D AT − 2 ∂ z − k A a s (C A − C AS ) l ∂ z
Como se trabaja en estado estacionario (hipótesis I): (2) 2 l
l
l d C A A εl D AT − ul dC dz dz 2
− k Als a s (C Al − C AS ) = 0
En el catalizador se propone el modelo: (3) s
ls l ′′ A ε s (1 − φ) ∂C ∂t = k A a s (C A − C AS ) − v A εs ρcat S a ηA r A
Como se trabaja en estado estacionario y rearreglando la ecuación 3 obtenemos: (4) k A a s (C A − C AS ) = v A εs ρcat S a ηA r A′′ ls
l
Sustituyendo (4) en (2) podemos simplificar el modelo a: (5) 2 l
l
l d C A A εl D AT − ul dC dz dz 2
− v A εs ρcat S a ηAr A′′ = 0
4.1 Planteamiento de los modelos. Una vez analizado el problema y realizado las consideraciones pertinentes se ha establecido que el modelo será unidireccional (hipótesis IV y V) y se procede a plantear los balances de materia para cada especie de estudio y a su vez se realizaran las simplificaciones adecuadas para llevar lograr así establecer el modelo a analizar. ➢
Modelo del Balance de materia para el DBT (A).
➢
Modelo del balance de materia para el H2
En la fase gas se considera flujo pistón, además de que en la fase gas se considera prácticamente que solo está compuesta por H 2 lo que da sustento al caso de que no existe dispersión en esta fase (hipótesis IX). (6)
εl
∂C g H ∂t
=− u g
2
∂C g H ∂ z
2
el BP (B)
l * − C l ) a (C H − k H H 2 2 l 2
En base a la hipótesis de que los productos solo se van al exterior una vez formados y no se vuelven a introducir en la fase liquida (hipótesis X) se obtiene el modelo: (13)
Trabajando en estado estacionario: (7) g
− u g
dC H
2
dz
l * − C l ) = 0 − k H a (C H H 2 2 l 2
2 l
l
l
∂C B l ∂ C B sl l B εl ∂C ∂t = εl D BT ∂ z 2 − ul ∂ z + k B as (C BS − C B )
En la fase liquida: (8) l ∂C H εl ∂t 2
=ε
l ∂2 C H 2 ∂ z 2
l ∂C H − ul ∂ z 2
ls l * − C l ) as (C Al − C AS ) + k H a (C H − k H H 2 2 l Como 2 se2 trabaja en estado estacionario (hipótesis I): Como se trabaja en estado estacionario (14) (hipótesis I): l 2 l (9) dC B l d C B u εl D BT − + k B sl as (C BS − C Bl ) = 0 l dz dz 2
D l l H 2 T
2 l
l
d C H l 2 D l H 2 T dz 2
ε
− ul
dC H
2
dz
ls l l * − C l ) = 0 − k H − as (C H C H S ) + k H al (C H H 2 2 2 2 2 2
En el catalizador se propone el modelo: (15)
En el catalizador se propone el modelo: (10) ε s (1 − φ)
s C H 2
∂
∂t
s
sl l ′′ B ε s (1 − φ) ∂C ∂t =− k B as (C BS − C B ) + v B εs ρcat S a ηB r B
ls l ′′ = k H − C H 2 S ) − v H 2 εs ρcat S a ηH 2 r H as (C H 2 2
2
Como se trabaja en estado estacionario y rearreglando la ecuación 15 obtenemos: (16)
Como se trabaja en estado estacionario y rearreglando la ecuación 10 obtenemos: (11) ls
l
k H as (C H 2
2
k B as (C BS − C B ) = v B εs ρc S a ηB r B′′ sl
′′ − C H 2 S ) = v H 2 εs ρcat S a ηH 2 r H 2
Sustituyendo (16) en (14) podemos simplificar el modelo a: (17)
Sustituyendo (11) en (9) podemos simplificar el modelo a: (12) 2 l
d C H 2 D l l H 2 T dz 2
ε
2 l
l
− ul
dC H
∂ z
2
l
− v H 2 εs ρcat S a ηH 2
r ′′
H 2
+
➢ ➢
Modelo del balance de materia para
l
l d C B B ε−l D − ul dC l 2 0 BT )dz = dz C H 2 H 2
(C *
l k H a 2 l
+ v B εs ρcat S a ηB r B′′ = 0
Modelo del balance de materia para el CHB (C)
(24) En base a la hipótesis de que los productos solo se van al exterior una vez formados y no se vuelven a introducir en la fase liquida se obtiene el modelo: (18) 2 l
l
r B
= k 1 C A (25)
r C = k 2 C A
l
∂C C l ∂ C C sl l C εl ∂C ∂t = εl DCT ∂ z 2 − ul ∂ z + k C as (C CS − C C )
Como se trabaja en estado estacionario (hipótesis I): (19) 2 l
l
dC C l d C C εl DCT − u 2 l dz dz
(26) − r D = (k1 + 3k 2 )C A Resumiendo los modelos y introduciendo las cinéticas propuestas de la literatura consultada se obtienes los modelos:
l + k sl ) =0 a (C CS − C C C s
(27) 2 l
En el catalizador se propone el modelo: (20) s ε s (1 − ϕ ) ∂C ∂t C
l =− k sl ) + v C εs ρcat S a ηC r ′′C a (C CS − C C C s
g
− u g
(C CS −
) = v C εs ρcat S a ηC r C
dC H
2
dz
l * − C l ) = 0 a (C H − k H H 2 2 l 2
(29) 2 l
εl D Hl 2 T
′′
l C C
+ ε s ηA kC A = 0 (28)
Como se trabaja en estado estacionario y rearreglando la ecuación 20 obtenemos: (21) sl k C as
l
l d C A A εl D AT − ul dC dz dz 2
d C H dz 2
2
l
− ul
dC H
2
∂ z
l a (C + ε s ηH 2 (k1 + 3k 2 )C A + k H 2 l
(30) Sustituyendo (21) en (19) podemos simplificar el modelo a: (22) 2 l
l
l d C C C εC DCT − ul dC dz dz 2
+ v C εs ρcat S a ηC r ′′C = 0
Tomando las cinéticas de la literatura estudiada para esta parte[3]: (23) − r A = k C A
2 l
l
dC B l d C B εl D BT − + ε s ηB k1 C A = 0 u 2 l dz dz
(31) 2 l
l
l d C C C εC DCT − ul dC + ε s ηC k2 C A = 0 dz dz 2
4.2 Condiciones de operación consideraciones pertinentes.
y
Para las pruebas iniciales se escogerá un WHSV de 5.0 g de alimentación/g cat. h, con base en la literatura se ha observado que
se busca una relación de 2000 a 1 (scf / Bbl), [6]Para el análisis inicial se tomara una longitud de 50 cm con lo que se calcularan los demás valores para las condiciones de operación y alimentación, obteniendo los siguientes valores: V R
= 8835.73 cm3
wcat = 2335.32 g −3
v l = 2.78 * 10 cm/ s v g = 0.0618 cm/ s F l = 0.1473 cm3 / s F g = 2.1860 cm3 / s
ε g = 0.2
Se establece una concentración inicial patrón del DBT pequeña debido a las consideraciones tomadas y a que en la realidad se presenta con estos valores: = 0.01
mol cm3
La concentración del hidrogeno en la fase gas inicial dado las condiciones de operación se estima un valor de: g
C H 0 2
4.3 Análisis del sistema del sistema difusión-reacción del catalizador. Se toma el modelo de la literatura para una partícula esférica general en estado estacionario: (32) 2
d C i dr 2
= 0.01 mol cm3
Estos son los valores de condiciones iniciales, y para los productos se considerara que no hay nada en la alimentación y el mismo caso para la concentración del
′′ 2 dC i ± ri S a ρcat
+
r dr
i
Def
=0
Modelo el cual se hará adimensional para el análisis de sus módulos de Thiele por medio de las ecuaciones: (33) γi =
εl = 0.3
L C A0
hidrogeno en el líquido.
C i C iS
λ = Rr
y
Para el caso de cada sustancia tanto reactivo y producto consultando la literatura se establece que cada uno tiene su propio módulo de Thiele y viene dado por la ecuación general: (34) φi = R
√
k i S a ρcat i
Def
Denotando que las constantes tienen la forma presentada en las ecuaciones 23, 24, 25 y 26. Con los módulos de Thiele se procede a calcular el factor de efectividad según la ecuación ya demostrada: (35) ηi = φ32 (φi coth φi − 1) i
4.3 Obtención de químicos y cinéticos.
parámetros
físicos,
✓
Difusividades en la fase liquida.
ecuación[2]: (39)
Debido a las bajas concentraciones de DBT y de los productos de reacción se puede despreciar a estos en el efecto de la Difusividad por lo que se hace la hipótesis de que la difusión es prácticamente binaria para la especia analizada en la tetralina. Se utiliza la siguiente expresión para el cálculo de la Difusividad binaria para el sistema analizado[2]: (36) DiT
✓
=
v0.267 T T 8.98 * 10−8 v0.433 μ T i
Calculo de la Difusividad efectiva.
Se estimaron los valores de la Difusividad efectiva por medio de la ecuación[2]: (37) D e,i
= τϵ (
1 D iT
1 + D1
✓
Datos cinéticos y característicos del Catalizador.
Estos datos se tomaran de la literatura consultada debido a la similitud de condiciones de operación. k = 18.04 k 1
= 15.76 1/h
k 2 S a
1/h
= 2.28 1/h
= 310 m2 / g
ρc = 33.0 lb/ f t 3 El modelo se resolverá utilizando el software Polymath 6.0 para sistemas de ecuaciones ordinarias acopladas en condiciones iniciales por reducción de orden.
)
k .i
Para el caso de estudio se supone que los poros del catalizador son totalmente cilíndricos por lo que la tortosidad del catalizador es 1.0. La Difusividad knusen se calcula mediante la ecuación[2]: (38) D k ,i
r g = S g2ρϵ p
0.5
T = 9700r g ( MW ) i
Donde rg se calcula mediante la
5.0 Resultados y discusión.
5.1 Obtención de parámetros Cada uno de los parámetros se obtuvo mediante las ecuaciones seleccionadas primero se obtuvieron los datos de Difusividades binarias los cuales se presentan en la tabla 1. Tabla 1. Difusividad binaria para las especies de interés. DiT por Volumen ecuación de 2 Especie molar DiT (cm /s) Wilke y (m3/kmol) Chang (cm2/s)
Tetralin 0.0226** ----a BDT 0.1555 5.518*10-4 3.108*10-4 H2 0.0038 2.752*10-3 2.881*10-3 BP 0.1503 6.599*10-4 3.172*10-4 CHB 0.1500* 5.604*10-4 3.175*10-4 *estimado respecto a su peso molecular a partir del BP y del Naftaleno ** Se estimó su valor con un valor cercano a la decalina.
Se realizó además con la correlación de Wilke y Chang y se observó que tienen los mismos órdenes de magnitud dando un mayor confiabilidad en los datos se tomaron los calculados con la ecuación propuesta ya que con esta se han realizado cálculos para este tipo de sistemas. Para el cálculo de las Difusividades efectivas y la knusen primero se calculó el radio de poro y para el sistema seleccionado tiene un valor de: r g= 7.079*10-9 m, con esto se obtuvo los datos de Difusividad necesarios los cuales se muestran en la tabla 2. Tabla 2. Difusividad efectiva y knusen para cada especie de interés.
Especie DBT H2 BP CHB
Difusividad Knusen (cm2/s) 0.01294 0.12409 0.01414 0.01387
Difusividad efectiva (cm2/s) 3.064*10-4 1.562*10 -3 3.675*10 -4 3.124*10-4
Especie DBT H2 BP CHB
Módulo de Factor de Thiele Efectividad 0.142 0.9987 0.070 0.9997 0.121 0.9990 0.05 0.9998
Se puede observar que los factores de efectividad tiene un valor alto esto se debe principalmente a que se está despreciando el efecto del factor de constricción y la tortuosidad del catalizador, es decir se está idealizando demasiado el catalizador lo que aleja de la realidad en modelo pero como una aproximación inicial permito modelar el sistema para un caso donde el catalizador tenga una efectividad alta, la tortuosidad para catalizadores comerciales anda en valores como de 4.0 [2] lo que nos modifica el módulo de Thiele y el factor de efectividad para el caso del DBT a un valor de 0.97 para el factor de efectividad y por lo que se puede observar que se deben contar con dicha información para aumentar la robustez del modelo. Los últimos parámetros que se calcularon fue el coeficiente de trasferencia de masa y el valor de la constante de Henry para dicho sistema obteniendo los valores de: k Li a L
= 0.05128
cm2 s
H = 792.89 atm
Si se observan los datos tiene valores dentro de los órdenes de magnitud coherentes a datos consultados en la literatura, lo cual permitió calcular los datos del módulo de Thiele para cada especie como se observa en la tabla 3.
En el caso del coeficiente se calculó con la correlación antes descrita y para la constante de Henry se hizo mediante extrapolación de datos para sistemas homólogos en un simulador de propiedades, Hysys Aspen.
Tabla 3. Módulo de Thiele y factor de efectividad para cada una de las especies reaccionantes.
Ya con los parámetros obtenidos se procedió a la resolución del sistema de
ecuaciones diferenciales con el auxilio del software Polymath 6.0, obteniendo para el sistema planteado un perfil de concentración en función de la longitud del reactor como el que se muestra en la fig. 3, observando así un comportamiento lógico esperado, vemos como la concentración del reactivo disminuye según se avance en el reactor y se observa que la concentración de H2 se mantiene prácticamente constante, aunque disminuye en una pequeña cantidad por lo que es razonable hacer la hipótesis de que se mantiene constante en el reactor, además si se alimenta en una mayor proporción se observa que el cambio es aún más insignificante, así mismo se observa como la concentración en el líquido del H2 es variable ya que esta está en intercambio directo de una fase a otra variando en función de la reacción y de los fenómenos de transporte presentes en estas dos fases.
operación y más específico se está desperdiciando una gran cantidad de reactor, lo que indica que para este sistema reactivo se puede aumentar las condiciones de flujo para así no modificar los otros parámetros ya calculados por lo que se decidió hacer un análisis para obtener el perfil de concentración del DBT en función de la longitud del reactor a diversas velocidades de la fase liquida para observar si al disminuir el tiempo de contacto es decir disminuir el WHSV y así como modifica la conversión del contaminante como se puede observar en la fig. 4
Fig 4. Perfil de concentración para el DBT a distintas velocidades superficiales en cm/s.
Fig 3. Perfil de concentración en función de la longitud del reactor
Pero al analizar estos valores sobretodo la distancia al que se ha alcanzado prácticamente una conversión total es muy pequeña lo que indica que la reacción es altamente efectiva para dichas condiciones de
En esta grafica se puede analizar que efectivamente al aumentar la velocidad de flujo de la alimentación esto disminuye la conversión del DBT a los productos desulfurados lo que nos dice que incluso para el caso particular se podría aumentar 100 veces el flujo de alimentación y aun así un reactor de 100 cm sería suficiente para lograr conversiones aceptables según las normas consultadas, sin descuidar el efecto que estos cambios tiene en la transferencia de masa del
sistema, pero aun así el modelo permitió analizar estos cambios sin la necesidad de una experimentación compleja y costosa. 6.0 Conclusiones.
Batch reactor stuhy to determine effectiveness factors of commercial HDS catalyst. [3] Y.S. Al-Zeghayer and B. Y. Jibril, Kinetics of Hydrodesulfurization of Dibenzothiophene on Sulfided Commercial Co-Mo/γAl 2O3 catalyst. [4] Rob. J. Berger, E. Hugh Stitt, Guy B. Marin,
En este trabajo se realizó el cálculo de un modelo matemático para un reactor tipo PBR de tres fases en el cual se llevara a cabo una DHS del DBT. Se piensa utilizar un catalizador comercial Co-Mo/γ-Al2O3, este modelo tiene la cualidad de apreciar el seguimiento de la concentración de los reactivos y productos para esta reacción, Se pretende reproducirlo en un simulador para un futuro desarrollo de un modelo planta piloto, hasta llegar a la puesta en marcha de un sistema con combustible y con una corriente de H2 gaseoso impuro, para lo que se adecuaran los parámetros , a partir del uso de este modelo matemático que permitirá analizar el efecto de los fenómenos de trasporte sobre la conversión del reactor. 7.0 Agradecimientos
M.C. José Luis Sánchez López, por sus conocimientos brindados durante todos los cursos del área. Dra. María Guadalupe Hernández Linares, por sus explicaciones de la química detrás de una reacción de HDS. Ing. Gabriel Guerrero Luna, por sus conocimientos brindados en la importancia de esta reacción para la ingeniería química. 8.0 Referencias
[1] L. H. Klemm, Joseph J. Karchesy "Dibenzothiophene from biphenyl and derivatives" Journal of Heterocyclic Chemistry, 1978, Volume 15 Issue 4, Pages 561 563. doi 10.1002/jhet.5570150407 [2] G. Marroquin, J. Ancheyta, C. Esteban, A
Freek Kapteijin, Jacob A. Moulijin, Chemical
reaction kinetics in practice. [5] Andres Mahecha-Botero, Said S.E.H Elnashaie, John R. Grace y C.Jim Lim, A comprehesive Approach to Reaction Engineering [6] F.S. Mederos, J. Ancheyta/Applied Catalysis A: General. [7] H. Scott Fogler, Elementos de Ingenieria de las reacciones quimicas.. [8]Benner,S. Evaluación a escala laboratorio de procesos de eliminación de hidrocarburos aromáticos.