Tarea 1

UNIVERSIDAD NACIONAL DE COLOMBIA
FACULTAD DE INGENIERIA
SIMULACION DE PROCESOS QUIMICOS
David Leonardo Barinas Patarroyo Andres David Tellez Rocha
Anlisis de reactor seudohomogneo
INTRODUCCION
Se desea analizar el comportamiento de un reactor seudohomogneo
unidimensional, donde se llevar acabo la oxidacin parcial de o-xileno para
obtener anhidro ftlico como nico producto; solo se considera el caso 1 donde
se obtiene como producto nico el anhidro ftlico
A B
Tenemos:
r A =k 1 P A PO 0
a
b
dP A
=A 1 P A e T
dz
a1
b
dT
=B1 P A e T C 1 (T T c )
dz
1
dTc
=D 1(T T c )
dz
A 1=
d b PM PO 0
Gg
B 1=
d b ( H 1) PO0
G g Cp g
C1 =
4U
G g Cp g d t
D 1=
n d t tn U
W c Cpg
Y las condiciones de frontera son las siguientes:

Arreglo en paralelo:
P A =P A 0
T =T 0
T c =T c 0
Arreglo en contracorriente:
P A =P A 0
z=0
T =T 0
T c =T cL
z=L
Para el desarrollo en la hoja de clculo se atribuye un valor negativo al flujo

msico del refrigerante ( W c <0 ) para el flujo en contracorriente y un valor del
flujo msico del refrigerante positivo para el arreglo en paralelo, y con un flujo
muy grande del refrigerante se analizar el comportamiento del sistema con
temperatura de este constante.
Tambin se tendr en cuenta las sensibilidades paramtricas segn el modelo
de reactor con el fin de verificar que el mtodo propuesto para la resolucin de
las ecuaciones diferenciales sea convergente; obteniendo las siguientes
ecuaciones:
Para llegar a las ecuaciones variacionales se desarroll el siguiente mtodo:
d
13=
P A
T c
dz
P A
a
T A1 a1 P A b T
12=
=
e
2
dz
T
S
d
P A
a
b
P A
=A 1 e T
dz
S
13=
S 1=S11 + S 12+ S13

a
b
d S1
a P S
=A 1 e T S 1+ 1 A2 2
dz
T
1
T
Tc
23=
=C1
dz
S
d
T
a
b
PA
21=
=B1 e T
dz
S
d
T
a
a1 b T
T
22=
=B1 P A 2 e
C 1
dz
T
S
d
S 2=S21 +S 22+ S 23
a
b
d S2
a P S
=B1 e T S 1+ 1 A2 2 C1 (S 2S 3)
dz
T
1
Tc
PA
31=
=0
dz
S
d
Tc
T
32=
= A1
dz
S
d
d
33=
Tc
Tc
=D 1
dz
S
S 3=S31 +S 32+ S 33
d S3
=D1 ( S2S 3)
dz
Teniendo las siguientes condiciones de frontera:
z=0
S 1=S2=0
S 3=1
RESULTADOS
Grafica 1
La primer grfica construida fue la grfica 1.1, se encuentra en la pestaa
Esquemas de enfriamiento dentro el archivo Excel adjunto; el mtodo
numrico empleado para la resolucin del sistema de ecuaciones diferenciales
fue el mtodo de Euler con un tamao de paso de 0,0015. Se verific la
reproducibilidad de los datos presentados en el artculo con los datos hallados
mediante la hoja de clculo, se nota que el arreglo en paralelo se tiene una
nica respuesta al modelo del reactor, mientras que para el arreglo en
contracorriente se presentan tres posibles soluciones y para una condicin
isotrmica del fluido de control se obtiene una lnea recta de pendiente 1; se
analizaron valores de flujo de la sustancia refrigerante de -20, -50, -100,
infinito, 100, 25 y 15, con el fin de comprender el comportamiento de la
temperatura del fluido de reaccin para un reactor de dos metros de longitud al
variar la temperatura del fluido de control.
En la grfica se evidencia que para el flujo en contracorriente se obtienen
multiplicidad de estados estacionarios mediante la forma de la grfica, para el
flujo isotrmico del fluido de control se obtiene una lnea recta de pendiente
uno y para el flujo en paralelo una nica respuesta.
-100
-50
-20
100
15
25
1000000
660
640
620
Tco = To K
600
580
560
540
575
595
615
635
655
675
695
TcL (K)
Grafica 2.
Se procedi a realizar la grfica 1.3 del artculo, donde se presentan los perfiles
de temperatura con arreglo en contracorriente, como ya se mencion
anteriormente con la posibilidad de soluciones mltiples y que se originen tres
estados estacionarios alto, medio y bajo; para hacer esta curva se modificaron
las conversiones, para una mnima conversin (x=0.2) se obtuvo la lnea baja
en la pestaa de soluciones mltiples, para el caso de una moderada
conversin (x=0.8) se obtuvo la lnea media, y para una alta conversin
(x=0.99) se obtuvo la lnea alta
855
825
795
765
735
Temperatura (K)
705
baja
alta
media
675
645
615
585
0.00
0.20
0.40
0.60
0.80
1.00
1.20
1.40
1.60
1.80
diastancia axial (m)
Grafica 3.
Se procedi con la construccin de la grfica 1.5, pestaa sensibilidad del
artculo para comparar los tres arreglos (paralelo, infinito y contracorriente) con
una longitud del reactor de 1.5 m, una conversin de 0.54 y un flujo msico de
refrigerante de 50 kg/s. lo que se puede concluir grficamente es que el punto
caliente se forma a menor longitud a condicin contracorriente, la tendencia
del arreglo en paralelo es una curva ms amortiguada y con un punto caliente
situado hacia la mitad del reactor, lo cual por aspectos de seguridad es el
arreglo ms recomendado.
2.00
670
665
660
655
650
Temperatura de la mezcla (K)
645
640
635
630
625
620
0.2 0.4 0.6 0.8
1.2 1.4 1.6
Longitud del reactor (m)

Paralelo
infinito
contracorriente
Grafica 4.
Se procedi con la construccin de la grfica 1.6, pestaa sensibilidad que
nos relaciona la sensibilidad paramtrica S2 como ya sabemos es referida a la
mezcla de reaccin con la longitud del reactor; y bajo anlisis de las mismas
condiciones del grfico anterior, la cual nos dice la sensibilidad al cambio que
tiene la temperatura de la mezcla reactante. Se concluye que para un arreglo
en contracorriente la temperatura presenta una gran sensibilidad al inicio del
reactor, variando entre 0 y 9.5; nuevamente se aprecia que para el arreglo en
paralelo la sensibilidad est amortiguada en un intervalo entre 0 y 3, y tambin
su pico esta desplazado un poco ms a la derecha respecto al arreglo en
contracorriente, por lo que se concluye que el arreglo ms apropiado para la
operacin de este tipo de reactor es el arreglo en paralelo.
12
10
6
Sensibilidad S2
4
0.2
0.4
0.6
0.8
-2
Longitud del reactor (m)
Paralelo
infinito
contracorriente
1.2
1.4
1.6

Tarea 1

Cargado por

Copyright:

Formatos disponibles

Tarea 1

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tarea 1

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD NACIONAL DE COLOMBIA

Y las condiciones de frontera son las siguientes:

Para el desarrollo en la hoja de clculo se atribuye un valor negativo al flujo

S 1=S11 + S 12+ S13

Teniendo las siguientes condiciones de frontera:

diastancia axial (m)

0.2 0.4 0.6 0.8

1.2 1.4 1.6

Longitud del reactor (m)

También podría gustarte