UNIVERSIDAD NACIONAL DEL CALLAO
FACULTAD DE INGENIERA QUMICA
ESCUELA PROFESIONAL DE INGENIERIA QUIMICA
Curso: Ingeniera de las Reacciones Qumicas II
Profesor: Ing. Leonardo Flix Machaca Gonzales
Tema:
Diseo de reactores para el tratamiento de Dioxinas
Integrantes:
1)
2)
3)
4)
5)
6)
Aquino Inofuente Julissa
Chuquihuamani Collantes, Cesar
Chavez Melgarejo Jorge
Hancco Ccolque Freddy
Prieto Fernandez Lizbeth
Soras Taype Crithian
CALLAO PERU
2015
INDICE
Pagina
RESUMEN
3
I.
INTRODUCCION
II.
4
FUNDAMENTOS DE DISEO
5
II.1
Tipos de reactores
5
II.2 Modelos de diseo
6
II.3 Factores de efectividad
15
III.
PROCEDIMIENTO DE DISEO
19
III.1 Bases de Diseo
19
III.2 Calculo de capacidad
32
III.3 Diseo de detalles
32
III.4 Especificaciones
34
III.5 Datos de construccin
34
III.6
Elaboracin del plano
36
IV.
DISCUSIN Y RESULTADOS
V.
40
REFERENCIA BIBLIOGRAFICA
41
RESUMEN
El presente trabajo se refiere al diseo de dos reactores, un
reactor homogneo no cataltico as como tambin del diseo
de un reactor heterogneo cataltico de lecho fijo, ambos para
el tratamiento de dioxinas las cuales son consideradas como
micro-contaminantes
persistentes,
que
se
forman
como
subproductos en diferentes procesos qumicos industriales de
combustin que implican al cloro y estando de esa manera es
como contaminan el medio ambiente.
Nuestro trabajo tiene como objetivo hallar la capacidad o
volumen de los reactores, luego dimensionar para cada caso un
reactor en donde se llevar acabo dicho tratamiento a partir de
las condiciones de flujo dadas. As como tambin teniendo
como base los fundamentos de diseo de reactores, los
principios de transferencia de masa, la cintica de reaccin y la
termodinmica.
Cabe recalcar que para el caso del reactor homogneo no
cataltico se dimensionar un reactor tubular con los detalles
correspondientes
para
el
tratamiento
de
la
dioxina.
Paralelamente en el caso del reactor heterogneo cataltico se
dispondr del lecho de partculas de un catalizador como lo es
el xido de titanio activado con xido de wolframio, este
catalizador ser el responsable de disminuir la energa de
activacin y as aumentar la velocidad de reaccin del proceso.
Todo el dimensionamiento, las especificaciones y datos de
construccin de ambos reactores lo hacemos con el uso del
programa Excel. Finalmente la representacin de nuestros
diseos de reactor es mediante el uso de un software o
programa de modelamiento en 2D.
I.
INTRODUCCION
Las dioxinas son sustancias qumicas micro-contaminantes
ms significativas en la actualidad que pertenecen a los
Bifenilos Policlorados (PCB), que son una serie de compuestos
aromticos
organoclorados
conocidos
tambin
como
Contaminantes Orgnicos Persistentes (COP). Los cuales son
sintetizados por el hombre en diferentes industrias como en la
fabricacin de aceites para transformadores, aceites hidrulicos
y
lubricantes,
pinturas,
pegamentos,
ceras,
etc.
Estas
sustancias toxicas se degradan muy lentamente por lo que
permanecen muchos aos en el ambiente (agua, suelo y aire)
causando grandes impactos en la salud pues muchos de ellos se
acumulan
en
la
piel
tejido
adiposo
por
tener
alta
liposolubilidad.
Un informe de OMS de 1986 nos dice que la mayor parte de
la ingesta diaria de dioxinas procede de los vertidos de residuos
clorados, gases de escape de automviles y ciertas industrias
del metal y papel. Adems cabe mencionar que el Per al haber
firmado un convenio como el de Estocolmo se ha comprometido
a establecer medidas necesarias para eliminar la produccin,
uso
de
estas
sustancias
qumicas
las
cuales
se
deben
identificar, etiquetar y retirar las concentraciones mayores a 50
partes por milln y volmenes superiores a 50ml a ms tardar
al 2025. Por todo lo mencionado anteriormente es que parte un
inters o importancia en disear equipos para el tratamiento o
eliminacin de dioxinas.
II.
FUNDAMENTO DEL DISEO DEL REACTOR.
II.1
Tipos
de
Reactores:
Existen
varias
formas
de
clasificarlos.
a) Segn el modo de operacin
Reactores discontinuos: Son aquellos que trabajan por cargas,
es decir se introduce una alimentacin, y se espera un tiempo
dado, que viene determinado por la cintica de la reaccin, tras
el cual se saca el producto.
Reactores continuos: Son todos aquellos que trabajan de forma
continua.
b) Segn las fases que albergan
Reactores homogneos: Tienen una nica fase, lquida o gas.
- Tipos de reactores Homogneos: Reactor batch, Semibatch,
CSTR, PFR.
Representacin de
un Reactor Batch
Representacin
de un reactor
CSTR
Reactores
heterogneos: Tienen varias fases, gas-slido, lquido-slido,
gas-lquido, gas-lquido-slido.
- Tipos de reactores Heterogneos:
a) Reactores
Heterogneos
Catalticos:
Reactor de lecho fijo, reactor de lecho
fluidizado y reactor de lecho empacado.
b) Reactores Heterogneos No Catalticos:
Reactor de lecho fijo, reactor de lecho
mvil,
reactor
de
lecho
fluidizado,
reactor de tanque agitado y reactor de
torre o columna.
II.2
Modelos de Diseo:
Diagrama bsico de un
Reactor lecho fluidizado
Modelos de Diseo para un Reactor Cataltico de Lecho
Fijo.
Para disear reactores de lecho fijo, existen diversos modelos
de diseo y estos modelos se desarrollan sobre la base de la
hiptesis que se plantea sobre el clculo y estos son:
a) Tipo de flujo:
-
Flujo pistn en una sola
direccin.
Flujo
pistn
dispersin
radial.
axial
con
y/o
Representacin de flujo
pistn con dispersin axial.
b) Cintica de las reacciones del proceso:
- Cintica pseudo homognea
- Cintica heterognea
Donde se define al flujo pistn como el paquete de la lnea de
corriente cuya distribucin no vara en el sentido radial,
hacindole solo axialmente.
As mismo, la cintica pseudo homognea se dice cuando el
factor de efectividad es igual a uno (n=1) y el sistema est
constituida de una sola fase de tal modo que las propiedades
del solido-fluido se presentan en conjunto. Por ejemplo: la
densidad, presin, velocidad de reaccin, etc.
En consecuencia, combinado un tipo de flujo con cualquiera de
las cinticas, obtenemos un determinado modelo de diseo as
tenemos:
1. Flujo pistn con cintica Pseudo Homogneo.
2. Flujo pistn con cintica heterognea.
3. Flujo pistn con dispersin axial y/o radial y con cintica pseudo
homognea.
4. Flujo pistn con dispersin axial y/o radial con cintica
heterognea.
A continuacin se detallan los modelos a utilizar en nuestro
diseo.
1. Modelo de flujo pistn con cintica pseudo homognea.
Recordamos que, el flujo pistn es definido como el flujo que se
caracteriza debido a que sus propiedades del sistema son
concentracin, presin y temperatura no varan en el sentido
radial, haciendo solamente en forma axial con la condicin de
que
en
esta
ltima
direccin
no
existe
gradientes
de
temperatura y la concentracin del transporte molecular es
despreciable.
La cintica pseudo homognea expresa la velocidad de
desaparicin de un reactante dado, pero sin discriminar las
etapas de transporte de masa ni calor y as mismo de la
reaccin qumica propiamente dicha. Es decir es la expresin
global, ya que tiene incluida los efectos de transporte de masa
y/o calor.
Lgicamente estas expresiones tienen origen experimental y
deben cuidarse varios aspectos al utilizarla. Como regla general
es necesario asegurar que las condiciones en que se van a
aplicar estas expresiones sean similares a las condiciones en
que se la obtuvo. Al respecto hay que verificar que los tamaos
de partcula sean iguales, al igual que el nmero de Reynolds de
la partcula, el tipo de flujo , el rango de temperaturas y
presiones de trabajo, la composicin de la alimentacin de
reactivos y las conversiones a alcanzar. Debe tenerse presente
que la modificacin de alguna de estas condiciones puede
traducirse en una modificacin de los parmetros de la ecuacin
cintica e incluso puede llegarse al caso de la inaplicabilidad.
En este caso las ecuaciones de balance de masa y energa son
las
usadas
en el
diseno
de
reactores
para
el
sistema
homogneo son:
Balance de Masa:
Entrada-
Salida=
Acumulacin
+Generacin
Qumica.
F A 0F A=0+ (r A ) dW
F A 0F Ad F A=0+ (r A ) dW
d F A =F A 0 dx A
(1)
(2)
(3)
de
Reaccin
(F A 0 dx A )=( r A ) dW
d xA
=r A
dW
FA 0
(4)
(5)
Balance de energa:
dT
( F Cp ) dZ
j
r j A T H
=Q v At
Vj
(6)
En la forma anloga:
dT
( F j Cp j) dZ
j
r j H Qv
=
Vj
B
(7)
Siendo:
Q=
4U
(T t )
D
(8)
La ecuacin (7) se transforma en :
dT
F j Cp j
(
) dZ
j
r j H 4 U (T t)
=
Vj
D B
(9)
Casos:
Caso I: Reactor de lecho fijo isotrmico (T f(Z))
A los fines de dimensionamiento del reactor se debe resolver la
ecuacin de Balance de Materia solamente y la resolucin del
balance de energa permitir encontrar las condiciones que
deben cumplirse para lograr la isotermia.
Caso II: Reactor adiabtico U=0
Balance de Masa:
10
FA 0
d xA
=r A
dW
(10)
Balance de Energa:
dT
F j Cp j
(
) dW
j
T =T 0 + x j
r j H
=0
Vj
(11)
(H ) F jo
Vj
( F Cp )
j
(12)
Caso III: Reactor no isotrmico no adiabtico (NINA)
La ecuacin de balance de materia:
dx r j
=
dw F jo
(13)
La ecuacin de balance de energa:
(H )F jo
dx
d w =rj V j F j Cp j
4 U ( Tt)
D B
( F Cp )
j
Con las condiciones iniciales: W=0 y
problema
consiste
en
(14)
la
X j=X jo
resolucin
simultanea
T=To, el
de
las
ecuaciones (13) y (14) puede resolverse por algn mtodo
adecuado.
Otros casos algo ms complejo, como lo son los auto trminos,
pueden ser resueltos aplicando las ecuaciones de masa y
energa.
Perdida de presin en los reactores de lecho fijo:
Para lechos fijos:
11
Para flujo laminar la ecuacin resultante es:
G
a
)
P 5G 2 a
=
=
3
L gc E
(15)
a= superficie del lecho
G= velocidad mxima (kg/hm2)
gc=constante de gravitacin
=porosidad
=densidad
=viscosidad
Para lechos turbulentos:
=
G
100
a
(16)
Re100
2
0.1
P 0.4 G a
G
=
=( )
3
L
a
gc E
(17)
Para nmero de Reynolds elevados:
P
=2 B
L
(18)
Se debe usar:
P 0.4 G2 a
G 0.1
=
=(
)
L
a
g c v E3
Para nmero de Reynolds elevados:
(19)
12
3n
( 1E)
P
2 fcG2
=
=(
3n
L gc Dp
E3
(20)
Dnde:
f= factor de friccin modificada
=factor de forma de la partcula (que est definida por la
relacin de la superficie de la partcula sobre el volumen de la
misma partcula)
n=1 laminar
Para Re<10 grado de turbulencia es despreciable
n=2 turbulento
La ecuacin Ergun:
2
P 150(1E) v 1.75(1E) Gv
=
+
L
Dp
gc E3
Dp 2
gc E 3
1.50 (1E)
Dp
2
(1 E)
P
Gv
=
L gc . . Dp E3
+1.75G)
(21)
(22)
2. Modelo de flujo pistn con cintica heterognea
Recordemos que el modelo con cintica heterognea discrimina
las fases existentes, considerando que en la fase fluida solo
existe transporte de materia y/o calor, mientras que en la fase
slida se produce la reaccin qumica. Considerando que las
partculas del catalizador, es conocido para que se lleve a cabo
la reaccin, es necesario que los reactantes superen la
resistencia de transferencia de masa externa e interna para
poder
alcanzar
la
superficie
cataltica.
Segn
sean
las
13
magnitudes de estas resistencias, habr diferencias en las
ecuaciones de balance de masa y calor. A los fines de darle una
mayor claridad del desarrollo del tema, se consideran dos casos
diferentes. As tenemos:
a) Flujo
Pistn
con
gradiente
externo
de
concentracin
temperatura
b) Flujo Pistn con gradiente externo e interno de concentracin y
temperatura
En el primer caso consideremos una partcula en el cual puede
existir
solamente
gradiente
eterno
de
concentracin
temperatura, mientras que para el segundo caso gradientes,
tanto en el interior de la partcula como as tambin en la
pelcula del fluido que lo rodea.
Consideremos ahora un reactor tubular relleno de partculas
esfricas. La concentracin del fluido y su temperatura se
modificaran en el sentido axial, como consecuencia a ello la
reaccin que se produce sobre el catalizador a travs del cual el
fluido circula. A su vez, para una posicin axial en el reactor, Z,
existirn partculas baadas por fluido con concentracin Cj (del
reactante j) y una temperatura T.
Alrededor de estas partculas puede definirse una pelcula
donde se establecern las resistencias a la transferencia
externas de masa y calor, de modo que las concentraciones de
j, y la temperatura en la superficie externa de las partculas
sern
C js
gradiente
T s . A su vez dentro de las partculas existir un
radial
denominara con
de
s
C js
temperatura
Ts
concentracin,
a dichas variables.
Caso I: Con gradiente externa de masa y temperatura.
se
14
En este caso el perfil de concentracin y temperatura dentro del
catalizador es despreciable por lo que se considera solamente el
perfil de concentracin y temperatura en la pelcula que le
rodea al catalizador.
Para plantear las ecuaciones de balance de masa, debe tenerse
presente que en la fase fluida existe transporte de masa por
circulacin debido a la corriente fluida que se mueve en
direccin axial. Para cada componente la variacin del flujo de
masa transportado por la fase fluida debe ser igual al flujo
transferido a travs de la pelcula en estado estacionario. A su
vez siempre bajo la hiptesis de estado estacionario, el flujo de
materia a travs de la pelcula debe igualar al flujo de materia
que reacciona sobre la superficie del catalizador.
Consideraciones similares pueden hacerse respecto al balance
de energa. En este caso, sin embargo, debe tenerse en cuenta
que la variacin entlpica del fluido se debe a un flujo de calor
intercambiado con el medio externo.
Todo lo dicho puede expresarse de la siguiente manera:
Cambio de j en dZ=Transporte de j a travs de la pelcula
Balance de Materia y Energa para la fase fluida:
Balance de masa:
Fj
dx j k c a
=
CjC sjs ) .(24)
(
dw p
Balance de calor:
dT
F jC j
(
) dw =
j
hc a s
( T T ) .. ( 25 )
p s
Para la fase slida:
15
Balance de masa
(C SJS , T SS ) .(26)
kc a
S
C j CJS
(
)=r j f
p
Balance de calor
( C SJS ,T SS ) .(27)
hc a s
Hr
T sT ) =
r f
(
p
v j
Donde
kc
(cm/seg)
y hc
en
(cal/C,
seg/cm 2), son
los
coeficientes peliculares de transferencia de masa y de calor
entre el reactor y el medio del intercambio trmico cuya
temperatura es T.
Estos sistemas de ecuaciones (24, 25, 26 y 27) estn sometidos
a las siguientes condiciones.
Z=0
w=0
Cj=Cjo . (27 a )
T = .(27 b)
En las ecuaciones presentes puede reconocerse a xj (o Cj), C jss,
T, Tss como variables dependientes de la posicin Z o de la
masa de catalizador w, lo cual puede considerarse como
variable independiente, de tal modo si se fija el valor de las
variables independientes (longitud del reactor o masa del
catalizador), o bien el valor de alguna de las variables
dependientes (tal como la conversin) en la posicin de salida
del reactor, el conjunto de ecuaciones (23) a (26) permiten
resolver el problema.
16
Debe indicarse que el conjunto de ecuaciones planteadas son
generales y contemplan adems el tratamiento de los distintos
regmenes
tcnicos
en
que
puede
funcionar
el
reactor
(isotrmico, adiabtico, etc.), para lo cual solo es suficiente fijar
las condiciones restrictivas correspondientes en cada caso, osea
U=0 para el adiabtico, T f(w).
Debe quedar en claro que frente a un problema cualquiera, esto
es, cualquiera sea la expresin cintica de la reaccin, an ms
compleja
del
tipo
Hougen-Watson,
la
resolucin
de
las
ecuaciones de (1) a (4) permite hallar la solucin del problema,
sin embargo, el tratamiento puede simplificar en casos donde la
expresin cintica sea del tipo ley de potencia.
r ij =k C nj . ( 28 )
Para la ecuacin (28) se puede encontrar expresiones del factor
de efectividad externo en condiciones de la fase gas.
rj ( T ss ,C sj ) =rj ( T 1 , C j ) ( T 1 , C j ) .(29)
Donde se puede obtener resolviendo la ecuacin (25) y (26)
para cintica de la ley de potencia.
n=( 1nDa ) n eo
nDa
( 30)
1+ nDa
Factor de efectividad externo: caso no isotrmico
kf ( T ) C n1
j
Da=
..(31)
kcf ( T ) a/ p
o=
E
(32 )
RT
17
( H ) . Kc . f ( T ) . a . Cj
.(33)
vj . hc . f ( T ) . a . T
En la ecuacin (29), (30), (31), (32) y (33) se observa que estn
en funcin de la temperatura y la concentracin del flujo que
reemplaza a la ecuacin (29), en (25) y en (26), se obtiene:
Kc . a
(C jC SJS )=rj. f ( T ,Cj ) .n . f ..(34)
p
hc . a s
Hr
( T sT )= v r . f ( T ,Cj ) . n . f (T ,Cj ) ..(35)
b
j
Reemplazando la ecuacin (34),(35) en (23)y (24)se obtiene :
Balance de masa:
F j
dx j
=r j f (T , C ) nf ( T ,C j ) .(36)
dw
Balance de energa:
dT
4U
( F C ) dw =r f ( T , C ) nf (T , C ) D ( T t ) ..(37)
j
hca
H
S
S
=( T S T ) =
K r A f ( T S ) ( 38)
b
VA
II.3
-
Factores de Efectividad:
Factor de efectividad interno isotrmico, definido por .
Donde:
T g=T ss=T s =T
Para una partcula esfrica
Una reaccin irreversible de 1er. Orden:
AR
18
N A =DeA
dC A
dr
De cintica de reaccin:
r A=K s CA n
r A
moles transformados o reaccionados
m2. s
Rd
N A=
Ks=
moles transferidos
m2. s
m 3 transferidos
m2. s
Aplicando el balance de masa para el reactante A en una
cascara como se aprecia en la figura de espesor dr ser:
d (densidad de flujo . area ( S T ) )
dV
=r A
St: rea de seccin transversal de la partcula esfrica.
S T =4 r 2
D eA
dCA
4 r2
dr
( )
4
V = 2
3
D eA
r: radio de partcula (no del poro)
dCA
4 r2
dr
( )
dV=4 r
dr
19
d C 2 2C A
K
( 2A + +
)= ev C nA
r dr
dr
DA
Ks . p . Sg
Kv=
p : densidad de particula
d C 2A 2 2 C A
p Sg K v n
( dr 2 + r + dr = De C A
A
Las condiciones lmites para la ecuacin diferencial:
R0
C A =C SAS
r=
dCA
dr = 0
r=0
Aplicando cambio de variables:
r
= R
Concentracin adimensional:
d C 2A 2 df Ro2 p Sg K s C SAS ' '
+
=
f
d2 dp
D eA
d C A 2 df
+
2 f ' '
2
dp =
d
f=
CA
C SAS
n1
Despejando :
=L
n1
K v C SAS
DeA
Llamado mdulo de Thiele para una Rx. irreversible.
L= longitud caracterstica de la partcula.
Para una partcula esfrica:
3
1
1
s =
S tanh s s
s =
3
( coth s1)
S s
Despejando
20
n1
R 0 K v .C SAS
=
e
3
DA
-
Factor
de
efectividad
definido por
n =
exotrmico
isotrmico,
est
velocidad observadao global
veloc . de la Rx . a las condiciones del fluido
La velocidad de la Rx. Qca n=1
S
r A =K v C AS
La velocidad de transferencia de masa/volumen utilizando la
tcnica de la pelcula.
C AgC
N A =K c a ( ASS )
Kc: coeficiente TM en funcin de a.
a: rea seccin interfacial de toda la partcula.
Sg
a= V particula
C AgC
K v C =K c a ( AS S )
S
AS
Definimos que NDa =
velocidad de reaccion Q
velocidad de transferencia de masa
Para una reaccin de orden n:
21
n1
N Da =
K v C Ag
k C
= v Ag
Kc a
K c a C Ag
Para orden n=1
N Da =
Kv
K ca
r obs =K v C SAS =
Kv C Ag
1+ N Da
K v C Ag
r
1+ N Da
n = obs =
r Aj K v. C Ag
Por lo tanto para una reaccin de 1er orden:
n =
1
1+ N Da
Factor de Efectividad Global Isotrmico:
Se plantea:
C AgC
S
K c a( ASS )=nr =RoK v C AS
Para n
r obs =K v n=
1
C Ag
nKv C Ag
=
K
1+n N Da
1+n v
Kc a
22
Como
nG =
III.
3.1
nKv C Ag
r Aobs 1+n N Da
nG =
=
r Ao
Kv C Ag
n
1+ nN DA
PROCEDIMIENTO DEL DISEO
Bases del diseo:
a) Escala de produccin: Nivel planta piloto
b) Tipo de Proceso: Oxidacin cataltica y no cataltica de la
c)
d)
e)
f)
Dioxina.
Condiciones de operacin:
Temperatura: 1000o C
Presin: 1 atm
Conversin: 0.998
Alimentacin: Qt= 30.57 m3/h = 0.5095 m3/s
= 0.0335 min = 2.01s
Propiedades del Fluido:
Composicin de la dioxina: 0.22%
Peso molecular: 322g/mol
Densidad del fluido: 0.28 kg/m3
Viscosidad del fluido: 5.10967x10-5 kg/m.s
Propiedades del catalizador:
Catalizador empleado: TiO2/WO2
rea especfica del catalizador (Sg): 101 m2/g
Dimetro de la partcula: 0.0001 m
Porosidad del lecho: 0.61
Densidad de la partcula (TiO2): 4.23 g/cm3
Cintica de las reacciones Qumicas:
Reaccin Qumica:
C12H4Cl4O2 + 11 O2 12CO2 + 4HC
Para el Sistema Homogneo:
Reactor Homogneo:
Ecuacin de diseo
23
XA
XA
VR
d XA
=
F A 0 r A
VR
d XA
=
v CA 0 0 r A
Como la cintica es de 1er orden.
r A=K CA
XA
VR
d
d XA
= XA =
v CA 0 0 K C A
K n A 0 (1X A )
V 0 (1+ A y a x A )
1X A
1+4.3,384 .106 X A d XA
x
1+ A Y A X A d XA
1
=
K 0
(1 X A )
A
XA
VR
CA 0
=
v CA O K 0
De los coeficientes estequiometricos de la reaccin
A=
( 12+4 )(1+ 11)
=4
1
n A 9.677107
=
=3.384106
nT
0.286
o
XA
1 (1+1,3536X A )
0.0335=
d XA
K 0
(1X A )
24
XA
XA
XA
1
5
0.0335 K=
d XA +1.353610
d XA
0 (1X A )
0 (1X A )
x A=0.998
Para
0.0335 K=6.2146+7.0612 . 105
K=185.5126 min1
VR
m3
V R =v 0 A =0,0335 min30.57
=1.024 m3
v0 A
min
V R =1.024 m 1.3=1.331 m
Para el Sistema Heterognea:
Para la reaccin:
A+B
R+ S
Se sigue el modelo de Langmuir-Hinshelwood-HougenWatson (LHHW)
Mecanismo: Doble Centro
A + X K 1 AX ( 1 )
B+ X K 2 BX . ( 2 )
I + X K I XI .(3)
AX + BX K 4 RX + SX .. ( 4 ) segunda etapa reaccion superficial
25
RX K 5 R + X (5)
SX K 6 S + X .(6)
tercera etapa desorcion
constantes de equilibrio:
De (1):
K 1=
CA
C =K 1 P A C X =K A P A C X (6)
P A C X A
De (2):
K 2=
CB
C =K 2 PB C X =K B P B C X (7)
PB C X B
De (3):
KI =
CI
C = K I P I C X (8)
P I C X I
De (5):
K 5=
PR C X
P CX
C R= R
=K R P R C X ( 9)
CR
K5
De (6):
K 6=
PS C X
P C X
C S= S
=K S PS C X (10)
CS
K6
Etapa controlante: Reaccin superficial
De (4):
A=K 4 C A C B
A=K 4 K A K B P A P B C 2X (10)
-
Balance de sitios:
[ X ] T =C A +C B + C R +C S +C I + C X
[ X ] T =K A P A C X + K B PB C X + K R P R C X + K S P S C X + K I P I C X +C X
[ X ] T =C X [ K A P A + K B PB + K R P R + K S P S + K I P I +1 ]
2
X
C =
[ X ]T
2
2
[ K A P A + K B P B + K R P R + K S PS + K I PI +1 ]
( 11 )
26
(11) EN (10)
A=
K 4 K A K B P A PB [ X ]T
[ K A P A + K B P B + K R PR + K S P S + K I P I + 1 ]
(12)
Haciendo:
K=K 4 K A K B [ X ] T
K R P R + K S PS + K I PI son mucho menor que 1+ K A P A + K B P B
entonces la ecuacion(12)
A=
K PA PB
[ K A P A + K B P B + K I PI +1 ]
CONSIDERACION:
1. Las constantes de adsorcin para nuestro estudio las
consideraremos despreciables para el proceso. Nuestra
adsorvabilidad es muy dbil.
PB , presin del oxgeno ,
2. La
por su concentracin la tomaremos
constante.
Nuestra ecuacin cintica se transforma en:
A=K P A
La constante cintica viene determinada por:
K=K 4 K A K B [ X ] T
K=K 4 [ X ] T
Constante Cintica Heterognea:
27
K4
AX+BX RX+SX
A=K C4
CB
K het =K 4 K A K B [ X ]T
K het =K hom K A K B [ X ] T
. (i)
Considerando:
Ka 1
KB 1
XT = 1.5 %
Khom =3.115 s-1
Khet =7.0088 s-1
Reemplazando en (i):
Modulo de Thiele:
=
Procedemos hallar
R 2KS
3 D eA
(ii)
D eA :
Donde a partir de datos de la Hoja Tcnica (ANEXO1)
Propiedades Fsicas
Dimetro de la partcula
100x10^-6m
rea especifica
101 m2/gr
Densidad del TiO2 (
4.23 gr/cm3
p
)
Volumen del poro
0.37 cm3/gr
28
Factor de contriccin (
)
Tortuosidad ( )
6.5
Calculamos la porosidad:
vp
=
v p+
1
partcula
(iii)
v p =volumen del poro
partcula =densidad del TiO2
=0.61
Hacemos uso del Mtodo de Dusty Gas para hallar la
Difusividad efectiva
1
1
= e +
e
D A D AK
Donde:
e
D AB
e
D AB=D AB x
D AB=0.63228 m /s
DeAB=0.63228
m2 0.61 x 1
x
s
6.5
( [ ])
1 Y A 1
D
e
AB
MA
MB
.(iv)
29
DeAB=0.05935
m2
s
Adems:
YA
Y A=
1
2 d 2A N A
Y A =3.49 x 1015
DeA =1.3126 x 1 07 m2 /s
Reemplazando en (vi):
Finalmente reemplazando en (ii):
=1.72 x 1 01
Factor de efectividad interno:
Para una partcula esfrica:
s =
3
1
1
S tanh s s
s =0.998
Factor de Efectividad Global:
nG =
Dnde:
N DA =
KV
KC a
n
1+ n N DA
(vi)
(v)
30
a=
6 x (1 b )
p x d p
K V =constante cintica
K C =constante de transferencia de masa
Reemplazando en (vi):
nG =0.939551
Numero de Dankeler:
(vii)
Kv: constante cintica
Kc: constante de transferencia de masa
Constante de Transferencia de Masa: La hallamos a
partir de la Correlacin de Kremsen.
K C=
Sh D AB (1 b )
d p b
1
Sc 3
1
S h=(Re 2
Re=
d p v
(1 )
v =velocidad del gas
=viscosidad cinemtica
v =
Q
A ST
31
A ST =
Sc=
2
D
4
D AB
Caida de presin (P): (interaccin fluido partcula)
Nos permite calcular el dimetro del reactor (D), utilizando el
balance de cantidad de movimiento, est planteada por:
P f U 2
+
F =0
X dP
La cada de presin es igual a la fuerza cintica
x = profundidad del lecho
f
= factor de friccin
dp = dimetro la partcula efectiva
U
= velocidad superficial del fluido
F = densidad del fluido
= porosidad
Integrando la ecuacin (1)
2
P
U
(
) = f * d P * F
L
P = P ingresa - P salida
dp = G *
= -(P salida-P ingresa)
Vparticula
S externa
Para una partcula esfrica: dp = dp
32
Dimetro equivalente: dp =
3d Pc
dPc
2+
LP c
Para el factor de correlacin se utiliza la correlacin de ERGUN:
f = 1.75+
dP UF
f
1 B
150
( 1 B ) (
)
B3
dPG
f
Velocidad de flujo msico:
m
G = A ST
4m
D 2
De la ecuacin de ERGUN tenemos que calcular la viscosidad de
la mezcla de los gases.
P
=
L
G
2
(1 )
g cd p
*(
150( 1 )
+1.75G
dp
(viii)
Utilizamos la ecuacin de ecuacin de Chapman Enskog para cada
gas:
Dnde:
=viscosidad
T = temperatura (k)
M= masa molecular
33
= dimetro de colisin caracterstico de lada molcula en angstrom
(Parmetro de lennard jones)
=funcin integral de colisin
Datos de los gases para calcular para cada uno de ellos su
viscosidad a 1000c
gas
M(g/mol)
(
A
)
O2
32
N2
C12 H 4 Cl4 O2
28
322
0.80256
84
0.75630
54
1.32294
3.46
7
3.79
8
Para el oxigeno
De tablas de parmetros Lennard Jonnes.
=106.7
O2
=3.467
T=1273.15 K
KT
O 2
O 2
1273.15
106.7
= 11.932
(con este valor se busca el parmetro
en la tabla)
=0.8025684
Calculamos su viscosidad del oxgeno a 1000c
34
= 2.6693*
10
321273.15
2
( 3.467 ) 0.8025684
= 5.58499* 10
poise (g/cm*s)
Para el nitrgeno
(De tablas de parmetros lennard jonnes)
=71.4
N2
=3.798
T=1273.15 K
KT
parmetro
N 2
1273.15
71.4
N 2
= 17.83123
(con este valor se busca el
en la tabla)
= 0.75630544
Calculamos su viscosidad del nitrgeno a 1000c
= 2.6693*
105
= 4.61794* 10
Para la
281273.15
2
( 3.7987 ) 0.75630544
poise (g/cm-s)
C12 H 4 Cl4 O2
(De tablas de parmetros lennard jonnes)
=1.21* T b
T b=446.5 C
35
=1.21*(446.5+273.15)
=870.7765
1
dioxina =1.18V 3
dioxina
234.282 31
(
)
=1.18*
1000
=0.7274363
T=1273.15 K
KT
1273.15
= 870.7765
parmetro
N 2
dioxina
= 1.462085
(con este valor se busca el
en la tabla)
= 1.3294
Calculamos su viscosidad de la dioxina a 1000 c
5
= 2.6693*
10
3221273.15
2
( 0.7274363 ) 1.3294
=0.02429 poise (g/cm*s)
Calculando la viscosidad de la mezcla de gases (Ecuacin
de Wilke (error 2%)
36
Especie
O2
(1)
N2
(2)
Fraccin
molar
Peso
molecular
Viscosidad
(poise)
0.209999929
32
5.58499* 10
0.78999734
28
4
4.61794* 10
3.3612*
322
2.429* 10
(3)C 12 H 4 Cl 4 O2
i
1
j
1
2
106
Mi/Mj
1
1.142857
14
0.099378
88
i /j
1
1.20896403
0.875
0.82715447
0.086956
52
0.0190147
0.533209
18
10.0625
43.5007829
11.5
52.5908818
1.172330
43
2.438429
66
1
0.02298809
ij
1
1.028342
22
0.543896
22
0.912209
73
1
xj*ij
1.010646
27
0.969465
32
2.122384
86
Reemplazando todos los parmetros en la ecuacin:
Por lo tanto m = 5.11x10-5 Kg/m.s
Para el Re
dP UF
f
37
F =
MmezclaP
RT
dP
Q=A*U
B = 0.61
Luego el factor de cada:
f = 1.75+
G = A ST
1 B
150
( 1 B ) (
)
3
4m
2
D
Dimetro del reactor:
Integrando y ordenando los parmetros usados para hallar la cada
de presin (P) en la Ec. de cantidad de movimiento, obtenemos la
siguiente ecuacin para hallar el dimetro del reactor:
D 6 D2 =0
...(ix)
Donde:
=
150 f (1EB )
4(1.75) dp' m
w
64 m
3
'
f dp (P) B
Adems:
f =viscosidad de la mezcla
1.75 (1E B )
E B3
38
EB = porosidad dellecho
dp' =dimetro equivalente de lapartcula
m=flujo
msico
w=cantidad de catalizador
f =densidad del fluido
P=cada de presin
B =densidad dellecho
Ecuaciones de diseo:
w=
Q
1
ln
K ' nG
1X A
L=
4w
2
D aire
(x)
...(xi)
Resolviendo el siguiente sistema de Ecuaciones:
nG =
P
=
L
G
2
(1 )
g cd p
n
1+ n N DA
*(
(vi)
150( 1 )
+1.75G
dp
D 6 D2 =0
w=
Q
1
ln
K ' nG
1X A
...(ix)
(x)
(viii)
39
L=
4w
D2 aire
...(xi)
De incognitas: L, W P, nG, D obtenemos:
D= 0.6m
nG = 0.9395
W= 111.35 gr
L= 1.4 m
P =5.57atm
g) Termodinmica de la Reaccin:
86
Keq= 3.65 x10
H (TC) = -4377.041 Kj/mol EXOTRMICO
S (TC) = -5578.274 J/mol
Calculo de Capacidad:
3.2
G =-5578.274
REACTOR HOMOGNEO NO CATALITICO
VR
m3
V R =v 0 A =0,0335 min30,57
=1,0241 m3
v0 A
min
V R =1.024 m3(1.3) 1.300 L
REACTOR HETEROGENEO CATALITICO
D R 2
V R=
L
4
V R=
0.62
0.3977 m3
4
40
V R =0.3977 m 1.3 520 L
3.3
t=
Diseo de Detalles:
REACTOR HOMOGNEO NO CATALITICO:
PD 1
+
2SE 8
t: espesor
P: presin total
D:dimetro del reactor
S:esfuerzo permisible a la tensin=68947.57kpa
E: factor de eficiencia de la junta
t=
101.3 kPa0.6 m
3
+3.17510 m
268947.57 kpa0.7
t=3.8 mm
REACTOR HETEROGENEO CATALITICO:
t=
PD 1
+
2SE 8
t=
101.33 kPa0.6 1
+
268947.570.7 8
t=3.8 mm
41
3.4
Especificaciones:
REACTOR NO CATALTICO HOMOGNEO
Materia de construccin del reactor: acero inoxidable
ASTM316
Tipo de fondo: bridado estndar
Capacidad: 1.300 m3
Acondicionador de temperatura: precalentador elctrico
Sensor de temperatura: termocupla
REACTOR CATALTICO HETEROGNEO
Materia de construccin del reactor: acero inoxidable ASTM316
Tipo de fondo: bridado estndar
Capacidad: 0.517 m3
Acondicionador de temperatura: pre calentador elctrico
Sensor de temperatura: termocupla
Catalizador empleado: TiO2 activado por WO3
3.5 Datos de construccin:
REACTOR HOMOGNEO NO CATALITICO
Dimetro interior= 6m
Espesor de la plancha=
3.8 mm
1
Plg
8
42
Longitud del tubo=4.7m
Longitud del lecho=3.62m
7 uniones universales de plg
01 rotametro para flujos bajos
vlvulas de globo de plg
1 tampn macho con rosca de de plg
1 tubo de de plg de 20 cm de longitud
bridas de fierro de 5 plg de dimetro
3 pirmetros digitales
3 contactadores
3 alarmas de incendido o apagado
20 niples de de plg
1 tablero de control
2 discos de soporte del reactor de 6plg de dimetro
REACTOR HETEROGENEO CATALITICO
Dimetro interior= 6m
Espesor de la plancha=
3.8 mm
1
Plg
8
Longitud del tubo=4.73m
Longitub del lecho=1.9m
7 uniones universales de plg
01 rotametro para flujos bajos
vlvulas de globo de plg
1 tampn macho con rosca de de plg
1 tubo de de plg de 20 cm de longitud
bridas de fierro de 5 plg de dimetro
3 pirmetros digitales
3 contactores
3 alarmas de incendido o apagado
20 niples de de plg
1 tablero de control
2 discos de soporte del reactor de 6plg de dimetro.
3.6 Elaboracin del Plano:
43
44
Vista Lateral Reactor Heterogneo:
Vista
Frontal
Vista
Superior
Vista
Superior
45
Vista Lateral Reactor Homogneo
Vista
Frontal
46
Vista
Superior
IV.
DI
S
SI
N
Y
RESULTADOS:
Las dioxinas son un grupo de
compuestos
qumicos
(las policlorodibenzofuranos (PCD
F)
las policlorodibenzodioxinas (PCDD)
y
que
devienen
de
contaminantes ambientales persistentes. Las dioxinas se en el
medio ambiente por todo el mundo y debido a su persistencia
se van acumulando a lo largo de la cadena alimentaria,
principalmente en el tejido adiposo de los animales por su
solubilidad en las grasas. Los productos de origen animal son
los mayores contribuyentes a la ingesta de dioxinas por los
humanos lo cual tiene una elevada toxicidad lo cual puede
provocar problemas de reproduccin y desarrollo, afectar el
sistema inmunitario, interferir con hormonas y, de ese modo,
causar cncer. Con relacin al catalizador usado
Basado en
informacin secundaria donde Se compara el comportamiento
de diferentes catalizadores para la eliminacin conjunta de NO y
1,2-diclorobenceno (o-DCB) de los gases de incineracin de
residuos slidos urbanos (RSU) donde se forma dioxina. Donde
se han estudiado
distintos catalizadores de VOx/TiO2, VOx-
TiO2/ZSM5, Cu/ZSM5, Fe/ZSM5 y MnOx y se han probado en una
47
corriente
gaseosa
que
simula
la
composicin
de
una
incineradora de RSU. La reduccin selectiva de NO tuvo lugar a
menor
temperatura
detectndose
un
que
la
oxidacin
compromiso
para
la
total
del
eleccin
o-DCB,
de
la
temperatura de operacin entre el aumento de la eliminacin
de o-DCB y la disminucin de la conversin de NO. La prdida
de conversin de o-DCB y NO en la reaccin simultnea
respecto
las
reacciones
independientes
sugiere
la
competencia de ambas reacciones sobre el mismo tipo de
centro activo. Asimismo, existe un compromiso entre actividad y
selectividad, si el catalizador ms activo a baja temperatura
pero presenta alta selectividad a N2O. VOx/TiO2 es ms activo
que Cu/ZSM5 pero la conversin de NO es significativamente
menor con la formulacin convencional (VOx/TiO2).De manera
que usamos el Ti02 activado por WO3 como nuestro catalizador
Reduciendo
significativamente
las
emisiones
de
dioxinas
emitidas por la quema de residuos slidos que provienen de la
incineracin no controlada de desechos que contengan cloro o
V.
derivados de benceno.
REFERENCIA BIBLIOGRFICA
1. LEVENSPIEL,. OCTAVE,
Ingeniera de las Reacciones
Qumicas. Edit. Revert. Barcelona, 2da Edicin 1981.
2. MACHACA; G., L., F., Apuntes de clases del curso de
Ingeniera de las Reacciones Qumicas II (catalticas y no
catalticas), 2011.
3. MACHACA; G., L., F., Anlisis, Diseo y construccin del
Reactor Cataltico de Lecho fijo a nivel banco para la
obtencin
de
derivados
de
etanol
Trabajo
de
Investigacin presentado al VRI_UNAC. Callao-Per, 2005.
4. SMITH, J. M. Ingeniera de la Cintica Qumica, Ed.
Cecsa, Nueva Edicin, Mxico 1987.
5. Perry, H. J. Chemical Engineers, Handbook MacGraw Hill
book Co. New York 4a Edit. 1963
48
6. TREYBAL, Operaciones con transferencia de masa.
7. [Link] Pag. 490
8. [Link]
eos_issuu
9. [Link]
ANEXOS
49
ANEXO 1 :
HOJA
TCNICA