TAREA NÚMERO 9
PRESENTADO POR:
Coordinador: NATALIA GONZALEZ MEJIA
C.C. 1.128.444.379
Auditor: SANTIAGO VASQUEZ AGUDELO
C.C. 80.873.955
Secretario: EDWIN G. SANCHEZ T.
C.C. 1.128.425.574
PROFESORA: CONSUELO MONTES
INGENIERIA DE LAS REACCIONES QUIMICAS
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERIA
INGENIERIA QUIMICA
MEDELLIN
2009
P8-21 Se efectuará la reacción irreversible de primer orden en fase liquida
A→ B
En un CSTR con camisa. Se alimenta A puro al reactor a razón de 0.5 gmol/min. La curva
de generación de calor para esta reacción y sistema de reactor,
− ∆H Rxn
º
G (T ) =
1 + 1/τ κ
Se muestra en la figura P8-21.
Información adicional:
Calor de reacción (constante): -100 cal/g mol de A
Capacidad calorífica de A y B: 2 cal/g mol·ºC
UA: 1 cal/min·ºC
Temperatura ambiente, Ta: 100ºC
a) ¿A qué temperatura se debe precalentar el fluido de entrada para que el reactor opere
con una conversión elevada?
Primero se debe hacer el balance de moles de A en el reactor (CSTR):
XFA0 XC A0
FA0 − FA + rAV = 0 ⇒ V = ⇒ τ =
− rA − rA
Se reemplaza la expresión de ley de velocidad y la expresión de concentración, teniendo
en cuenta que para una reacción en fase liquida el volumen es constante:
donde
− rA = kC A FA
CA =
v0
Se obtiene:
X
τk=
1− X
Cuando se reemplaza esta expresión de en la ecuación de generación de calor
τk
dada, se obtiene que la generación de calor es igual a:
G (T ) = − X ∆H Rxn
º
De este resultado, podemos afirmar que la conversión será alta para temperaturas de
operación mayores que 170ºC. Para temperaturas menores que esta, podemos asegurar
que la conversión es mala. Además sabemos que del balance de energía en un CSTR,
la expresión de conversión obtenida es:
UA
FA 0
( Ta − T ) + ∑ θ i C pi ( T − T0 ) UA
FA 0
( Ta − T ) + C pA ( T − T0 )
X BE = ⇒ X BE =
− ∆H Rxn − ∆H Rxn
Donde pues se alimenta A puro al reactor. De esta expresión, podemos
∑θ C i pi = C pA
asegurar que el término en el numerador representa el calor retirado del sistema por
intercambio de calor. Reorganizando términos, se obtiene que:
UA kTa + T0
Qret = C pA (1 + k )( T − Tc ) donde k= y Tc =
C pA FA0 1+ k
Se dan valores de y se grafican las rectas resultantes en la misma gráfica de G(T).
T0
Se obtiene la siguiente gráfica:
De la grafica podemos ver que para una temperatura de entrada de 205ºC, el reactor
tiene tres temperaturas de estado estable. En este punto no es recomendable trabajar ya
que cualquier cambio leve de temperatura puede hacer que la temperatura de operación
en el reactor caiga rápidamente, y por consiguiente se baja la conversión. Entonces se
desea operar el reactor de manera que se puede asegurar una alta conversión. Esto se
puede afirmar para temperaturas mayores que 210ºC (sin incluir esta temperatura). Para
temperaturas mayores que 210ºC, la recta de calor retirado solo tendrá una intercepción
con la curva de calor generado, y no existe la posibilidad de un bajón abrupto de
temperatura, y la conversión será alta.
b) Determine la temperatura correspondiente del fluido en el CSTR a esa temperatura de
entrada.
La temperatura del fluido dentro del CSTR corresponde a la temperatura de operación
del reactor y se encuentra con la intersección de la curva de generación de calor con la
recta de calor retirado. Según la gráfica construida en la parte (a), se pueden ver dos
intercepciones cuando se usa To = 210ºC. La primera intersección es en T = 160ºC y la
otra es en T = 178ºC. De estas dos temperaturas, es obvio que se prefiere operar a la
temperatura mayor, pues es allí donde se alcanzan altas conversiones.
c) Suponga que ahora el fluido se calienta 5ºC más que la temperatura de la parte a) y
luego se enfría 10ºC, permaneciendo en ese nivel. ¿Cuál será la conversión?
El procedimiento descrito consiste en calentar hasta 215ºC luego enfriar hasta 205ºC. Se
lee el punto de intersección entre las curvas de calor generado y calor retirado para una
temperatura de entrada de 205ºC y se despeja de la expresión: y
X G (T ) = − X ∆H Rxn
º
se encuentra que .
X ≈ 0.88
d) Calcule la temperatura de extinción para este sistema de reacción.
De la grafica construida en la parte (a), se toman los interceptos de la curva de calor
generado con las diferentes rectas de calor retirado para construir la gráfica de TS vs. T0.
Se obtienen los siguientes datos:
T0 TS
170 139
180 144
190 149
198 154 170
205 158 164 175
210 160 178
220 184
La temperatura de extinción de este sistema de reacción corresponde a la primera
línea vertical (T0 = 198ºC). Si el reactor está operando a 170ºC con una temperatura
de alimentación de 198ºC, cualquier cambio leve de temperatura puede hacer caer
bruscamente la temperatura, y no se obtendría una alta conversión.
P8-30 (Múltiples reacciones con efectos térmicos) El estireno se puede producir a partir de
etilbenceno por la siguiente reacción:
Etilbenceno ↔ estireno + H2
[J. Zinder y B. Subramaniam, Chem. Eng. Sci., 49, 5585 (1194)] Se alimenta etilbenceno a
razón de 0.00344 kmol/s a un reactor PFR de 10.0 m3 junto con vapor de agua inerte a una
presión total de 2.4 atm. La proporción molar vapor de agua/etilbenceno es inicialmente [es
decir, partes (a) y (c)] 14.5:1, pero puede variarse. Dados los datos siguientes, determine las
velocidades de flujo molar de estireno, benceno y tolueno a la salida, para las siguientes
temperaturas de entrada cuando el reactor se opera adiabáticamente.
a) T0= 800 K
Primero, se tiene que hacer un balance de masa para cada especie que entra al reactor
o que interviene en una de las reacciones:
rEB=dFEBdV rSt=dFStdV rH2=dFH2dV rB=dFBdV
rEt=dFEtdV rT=dFTdV rMet=dFMetdV
Luego del balance de energía, teniendo en cuenta que el reactor opera adiabáticamente,
se obtiene la expresión:
dTdV=-
r1S∆H1S+r2B∆H2B+r3T∆H3TFEBCpEB+FStCpSt+FH2CpH2+FBCpB+FEtCpEt+FTC
pT+FMetCpMet+FH2OCpH2O
La presión total es dada (P = 2.4 atm) y según la ley Raoult:
Pi = yiP y
yi = Fi/FT
De esta manera se encuentran expresiones de presión para utilizar en las leyes de
velocidad.
De la relación molar de vapor de agua/etilbenceno, se deduce que:
FH2O = 14.5 FEBo
Luego se encuentran las velocidades de reacción netas para cada componente:
rEB=-r1S+r2B+r3T rEt=r2B
rSt=r1S rH2=r1S-r3T
rB=r2B rMet=r3T
rT=r3T
Se resuelven estas ecuaciones simultáneamente en PolyMath, y se obtiene que para T0=
800 K, los flujos de estireno, tolueno y benceno son:
FSt=8.975×10-4 mol/s
FT=3.589×10-5 mol/s
FB=1.078×10-5 mol/s
POLYMATH Results
No Title 02-26-2009, Rev5.1.233
Calculated values of the DEQ variables
Variable initial value minimal value maximal value final value
V 0 0 10, 10,
Feb 0,00344 0,0024958 0,00344 0,0024958
Fst 0 0 8,975E-04 8,975E-04
Fb 0 0 1,078E-05 1,078E-05
Ft 0 0 3,589E-05 3,589E-05
Fh2 0 0 8,616E-04 8,616E-04
Fet 0 0 1,078E-05 1,078E-05
Fmet 0 0 3,589E-05 3,589E-05
T 800, 765,25883 800, 765,25883
P 2,4 2,4 2,4 2,4
Feb0 0,00344 0,00344 0,00344 0,00344
Fh2o 0,04988 0,04988 0,04988 0,04988
Ftt 0,05332 0,05332 0,0542283 0,0542283
Peb 0,1548387 0,1104588 0,1548387 0,1104588
Ph2 0 0 0,0381329 0,0381329
Pst 0 0 0,0397212 0,0397212
b5 1,302E-06 1,302E-06 1,302E-06 1,302E-06
p 2137, 2137, 2137, 2137,
q 0,4 0,4 0,4 0,4
r2B 2,991E-06 5,165E-07 2,991E-06 5,165E-07
r3t 0 0 4,197E-06 4,153E-06
b6 -0,004931 -0,004931 -0,004931 -0,004931
b1 -17,34 -17,34 -17,34 -17,34
b2 -1,302E+04 -1,302E+04 -1,302E+04 -1,302E+04
b3 5,051 5,051 5,051 5,051
b4 -2,314E-10 -2,314E-10 -2,314E-10 -2,314E-10
Kp1 0,0459123 0,0196664 0,0459123 0,0196664
r1S 2,138E-04 2,483E-05 2,138E-04 2,483E-05
reb -2,167E-04 -2,167E-04 -2,95E-05 -2,95E-05
rst 2,138E-04 2,483E-05 2,138E-04 2,483E-05
rb 2,991E-06 5,165E-07 2,991E-06 5,165E-07
rt 0 0 4,197E-06 4,153E-06
ret 2,991E-06 5,165E-07 2,991E-06 5,165E-07
rh2 2,138E-04 2,068E-05 2,138E-04 2,068E-05
rmet 0 0 4,197E-06 4,153E-06
ODE Report (RKF45)
Differential equations as entered by the user
[1] d(Feb)/d(V) = reb
[2] d(Fst)/d(V) = rst
[3] d(Fb)/d(V) = rb
[4] d(Ft)/d(V) = rt
[5] d(Fh2)/d(V) = rh2
[6] d(Fet)/d(V) = ret
[7] d(Fmet)/d(V) = rmet
[8] d(T)/d(V) = -(r1S*118000+r2B*105200+r3t*(-53900))/
(Feb*299+Fst*273+Fh2*30+Fet*90+Ft*249+Fb*201+Fmet*68+Fh2o*40)
Explicit equations as entered by the user
[1] P = 24e-1
[2] Feb0 = 344e-5
[3] Fh2o = 145e-1*Feb0
[4] Ftt = Feb+Fst+Fb+Ft+Fh2+Fet+Fmet+Fh2o
[5] Peb = Feb/Ftt*P
[6] Ph2 = Fh2/Ftt*P
[7] Pst = Fst/Ftt*P
[8] b5 = 1302e-9
[9] p = 2137
[10] q = 4e-1
[11] r2B = p*(1-q)*exp(132392e-4-25000/T)*Peb
[12] r3t = p*(1-q)*exp(2961e-4-11000/T)*Peb*Ph2
[13] b6 = -4931e-6
[14] b1 = -1734e-2
[15] b2 = -13020
[16] b3 = 5051e-3
[17] b4 = -2314e-13
[18] Kp1 = exp(b1+b2/T+b3*ln(T)+((b4*T+b5)*T+b6)*T)
[19] r1S = p*(1-q)*exp(-8539e-5-10925/T)*(Peb-Pst*Ph2/Kp1)
[20] reb = -r1S-r2B-r3t
[21] rst = r1S
[22] rb = r2B
[23] rt = r3t
[24] ret = r2B
[25] rh2 = r1S-r3t
[26] rmet = rt
Independent variable
variable name : V
initial value : 0
final value : 10
b) T0= 930 K
Utilizando el mismo procedimiento que la parte (a), se obtiene que para T0= 930 K, los
flujos de estireno, tolueno y benceno son:
FSt=1.937×10-3 mol/s
FT=2.037×10-4 mol/s
FB=2.167×10-4 mol/s
POLYMATH Results
No Title 02-26-2009, Rev5.1.233
Calculated values of the DEQ variables
Variable initial value minimal value maximal value final value
V 0 0 10, 10,
Feb 0,00344 0,001083 0,00344 0,001083
Fst 0 0 0,0019366 0,0019366
Fb 0 0 2,167E-04 2,167E-04
Ft 0 0 2,037E-04 2,037E-04
Fh2 0 0 0,0017417 0,001733
Fet 0 0 2,167E-04 2,167E-04
Fmet 0 0 2,037E-04 2,037E-04
T 930, 850,57827 930, 850,57827
P 2,4 2,4 2,4 2,4
Feb0 0,00344 0,00344 0,00344 0,00344
Fh2o 0,04988 0,04988 0,04988 0,04988
Ftt 0,05332 0,05332 0,0554734 0,0554734
Peb 0,1548387 0,0468544 0,1548387 0,0468544
Ph2 0 0 0,0753926 0,0749755
Pst 0 0 0,0837866 0,0837866
b5 1,302E-06 1,302E-06 1,302E-06 1,302E-06
p 2137, 2137, 2137, 2137,
q 0,4 0,4 0,4 0,4
r2B 2,36E-04 5,804E-06 2,36E-04 5,804E-06
r3t 0 0 3,105E-05 1,465E-05
b6 -0,004931 -0,004931 -0,004931 -0,004931
b1 -17,34 -17,34 -17,34 -17,34
b2 -1,302E+04 -1,302E+04 -1,302E+04 -1,302E+04
b3 5,051 5,051 5,051 5,051
b4 -2,314E-10 -2,314E-10 -2,314E-10 -2,314E-10
Kp1 0,6304376 0,1397033 0,6304376 0,1397033
r1S 0,001442 5,871E-06 0,001442 5,871E-06
reb -0,001678 -0,001678 -2,632E-05 -2,632E-05
rst 0,001442 5,871E-06 0,001442 5,871E-06
rb 2,36E-04 5,804E-06 2,36E-04 5,804E-06
rt 0 0 3,105E-05 1,465E-05
ret 2,36E-04 5,804E-06 2,36E-04 5,804E-06
rh2 0,001442 -8,776E-06 0,001442 -8,776E-06
rmet 0 0 3,105E-05 1,465E-05
ODE Report (RKF45)
Differential equations as entered by the user
[1] d(Feb)/d(V) = reb
[2] d(Fst)/d(V) = rst
[3] d(Fb)/d(V) = rb
[4] d(Ft)/d(V) = rt
[5] d(Fh2)/d(V) = rh2
[6] d(Fet)/d(V) = ret
[7] d(Fmet)/d(V) = rmet
[8] d(T)/d(V) = -(r1S*118000+r2B*105200+r3t*(-53900))/
(Feb*299+Fst*273+Fh2*30+Fet*90+Ft*249+Fb*201+Fmet*68+Fh2o*40)
Explicit equations as entered by the user
[1] P = 24e-1
[2] Feb0 = 344e-5
[3] Fh2o = 145e-1*Feb0
[4] Ftt = Feb+Fst+Fb+Ft+Fh2+Fet+Fmet+Fh2o
[5] Peb = Feb/Ftt*P
[6] Ph2 = Fh2/Ftt*P
[7] Pst = Fst/Ftt*P
[8] b5 = 1302e-9
[9] p = 2137
[10] q = 4e-1
[11] r2B = p*(1-q)*exp(132392e-4-25000/T)*Peb
[12] r3t = p*(1-q)*exp(2961e-4-11000/T)*Peb*Ph2
[13] b6 = -4931e-6
[14] b1 = -1734e-2
[15] b2 = -13020
[16] b3 = 5051e-3
[17] b4 = -2314e-13
[18] Kp1 = exp(b1+b2/T+b3*ln(T)+((b4*T+b5)*T+b6)*T)
[19] r1S = p*(1-q)*exp(-8539e-5-10925/T)*(Peb-Pst*Ph2/Kp1)
[20] reb = -r1S-r2B-r3t
[21] rst = r1S
[22] rb = r2B
[23] rt = r3t
[24] ret = r2B
[25] rh2 = r1S-r3t
[26] rmet = rt
Independent variable
variable name : V
initial value : 0
final value : 10
c) T0= 1100 K
Utilizando el mismo procedimiento que la parte (a), se obtiene que para T0= 1100 K, los
flujos de estireno, tolueno y benceno son:
FSt=1.654×10-3 mol/s
FT=1.274×10-4 mol/s
FB=1.608×10-3 mol/s
POLYMATH Results
No Title 02-26-2009, Rev5.1.233
Calculated values of the DEQ variables
Variable initial value minimal value maximal value final value
V 0 0 10, 10,
Feb 0,00344 5,123E-05 0,00344 5,123E-05
Fst 0 0 0,0017492 0,0016538
Fb 0 0 0,0016076 0,0016076
Ft 0 0 1,274E-04 1,274E-04
Fh2 0 0 0,0016486 0,0015265
Fet 0 0 0,0016076 0,0016076
Fmet 0 0 1,274E-04 1,274E-04
T 1100, 980,70596 1100, 981,65258
P 2,4 2,4 2,4 2,4
Feb0 0,00344 0,00344 0,00344 0,00344
Fh2o 0,04988 0,04988 0,04988 0,04988
Ftt 0,05332 0,05332 0,0565892 0,0565814
Peb 0,1548387 0,002173 0,1548387 0,002173
Ph2 0 0 0,0699442 0,0647476
Pst 0 0 0,0742029 0,0701497
b5 1,302E-06 1,302E-06 1,302E-06 1,302E-06
p 2137, 2137, 2137, 2137,
q 0,4 0,4 0,4 0,4
r2B 0,015039 1,363E-05 0,015039 1,363E-05
r3t 0 0 1,068E-04 3,298E-06
b6 -0,004931 -0,004931 -0,004931 -0,004931
b1 -17,34 -17,34 -17,34 -17,34
b2 -1,302E+04 -1,302E+04 -1,302E+04 -1,302E+04
b3 5,051 5,051 5,051 5,051
b4 -2,314E-10 -2,314E-10 -2,314E-10 -2,314E-10
Kp1 7,6863544 1,4544589 7,6863544 1,4761353
r1S 0,0088597 -1,641E-05 0,0088597 -1,562E-05
reb -0,0238986 -0,0238986 -1,306E-06 -1,306E-06
rst 0,0088597 -1,641E-05 0,0088597 -1,562E-05
rb 0,015039 1,363E-05 0,015039 1,363E-05
rt 0 0 1,068E-04 3,298E-06
ret 0,015039 1,363E-05 0,015039 1,363E-05
rh2 0,0088597 -2,015E-05 0,0088597 -1,892E-05
rmet 0 0 1,068E-04 3,298E-06
ODE Report (RKF45)
Differential equations as entered by the user
[1] d(Feb)/d(V) = reb
[2] d(Fst)/d(V) = rst
[3] d(Fb)/d(V) = rb
[4] d(Ft)/d(V) = rt
[5] d(Fh2)/d(V) = rh2
[6] d(Fet)/d(V) = ret
[7] d(Fmet)/d(V) = rmet
[8] d(T)/d(V) = -(r1S*118000+r2B*105200+r3t*(-53900))/
(Feb*299+Fst*273+Fh2*30+Fet*90+Ft*249+Fb*201+Fmet*68+Fh2o*40)
Explicit equations as entered by the user
[1] P = 24e-1
[2] Feb0 = 344e-5
[3] Fh2o = 145e-1*Feb0
[4] Ftt = Feb+Fst+Fb+Ft+Fh2+Fet+Fmet+Fh2o
[5] Peb = Feb/Ftt*P
[6] Ph2 = Fh2/Ftt*P
[7] Pst = Fst/Ftt*P
[8] b5 = 1302e-9
[9] p = 2137
[10] q = 4e-1
[11] r2B = p*(1-q)*exp(132392e-4-25000/T)*Peb
[12] r3t = p*(1-q)*exp(2961e-4-11000/T)*Peb*Ph2
[13] b6 = -4931e-6
[14] b1 = -1734e-2
[15] b2 = -13020
[16] b3 = 5051e-3
[17] b4 = -2314e-13
[18] Kp1 = exp(b1+b2/T+b3*ln(T)+((b4*T+b5)*T+b6)*T)
[19] r1S = p*(1-q)*exp(-8539e-5-10925/T)*(Peb-Pst*Ph2/Kp1)
[20] reb = -r1S-r2B-r3t
[21] rst = r1S
[22] rb = r2B
[23] rt = r3t
[24] ret = r2B
[25] rh2 = r1S-r3t
[26] rmet = rt
Independent variable
variable name : V
initial value : 0
final value : 10
d) Encuentre la temperatura de entrada ideal para la producción de estireno con una
proporción vapor de agua/etilbenceno de 58:1.
Se toma el mismo procedimiento que en las partes (a)-(c), variando la proporción de
vapor de agua/etilbenceno y la temperatura de alimentación. Para cada variación de
temperatura se calcula el flujo de estireno producido y se obtienen los siguientes datos
para construir una gráfica de FSt vs.T0:
T0
FSt
4,8670E-
800 04
9,4390E-
850 04
1,5354E-
900 03
2,0661E-
950 03
2,2170E-
975 03
2,2453E-
985 03
2,2518E-
990 03
2,2533E-
995 03
2,2495E-
1000 03
1,9628E-
1050 03
1,4619E-
1100 03
Como se puede observar tanto en la tabla como en la gráfica, para una temperatura de
entrada alrededor de 995 K, el flujo de estireno es máximo.
e) Determine la proporción de vapor de agua/etilbenceno ideal para la producción de
estireno a 900 K.
Se toma el mismo programa construido en PolyMath en las partes (a)-(c), se fija un valor
de T0 = 900 K y se varia la proporción de vapor de agua/etilbenceno desde 10 hasta 35,
para ver el efecto sobre la producción de estireno.
H2OEB FSt
1,5420E-
10
03
1,7554E-
15
03
1,8639E-
20
03
1,8932E-
24
03
1,8939E-
25
03
1,8927E-
26
03
1,8731E-
30
03
1,8265E-
35
03
Como se puede observar tanto en la tabla como en la gráfica, para una proporción de
agua/etilbenceno alrededor de 25, el flujo de estireno es máximo.
f) ¿Cual considera que sea el objetivo de este problema?
El objetivo del problema es aprender a maximizar el flujo de un producto deseado, con
respecto a otras variaciones como la temperatura de alimentación y proporción de
inertes alimentados al reactor.
g) Haga otra pregunta o sugiera otro calculo que pueda hacerse en este problema.
Se podría variar la presión total del sistema, y ver su efecto sobre la producción de
estireno. También se podría variar la cantidad absoluta de flujo entrante al reactor, es
decir, un flujo molar de etilbenceno mayor o menor a la cantidad dada en el problema.
Otro parámetro que se puede variar es el volumen de reactor; se podría observar cómo
cambia la producción de estireno con volúmenes de reactor mayores o menores que el
volumen dado.
Información adicional:
Capacidades caloríficas
Metano 68 J/mol·K Etilbenceno 299 J/mol·K
Etileno 90 J/mol·K Hidrógeno 30 J/mol·K
Benceno 201 J/mol·K Vapor de agua 40 J/mol·K
Tolueno 249 J/mol·K
Estireno 273 J/mol·K
ρ = 2137 kg/m3 de pastilla
Φ = 0.4
∆HRxn1 = 118,000 kJ/kmol de etilbenceno
∆HRxn2 = 105,200 kJ/kmol de etilbenceno
∆HRxn3 = -53,900 kJ/kmol de etilbenceno
b
K p1 = exp b1 + 2 + b3 ln(T ) + [ ( b4T + b5 )T + b6 ] T atm
T
b1 = −17.34 b4 = −2.314 × 10 −10
b2 = −1.302 × 10 4 b5 = 1.302 × 10−6
b3 = 5.051 b6 = −4.937 × 10−3
Las leyes de velocidad para la formación de estireno (St), benceno (B) y tolueno (T),
respectivamente, son las siguientes. (EB = etilbenceno.)
10,925 PSt PH 2
r1S = ρ (1 − φ ) exp 0.08539− PEB − (kmol / m 3 ·s )
T K p1
25,000
r2 B = ρ (1 − φ ) exp13.2392 − ( PEB ) (kmol / m ·s )
3
T
r3T = ρ (1 − φ ) exp 0.2961−
11,000
( ) 3
PEB PH 2 ( kmol / m ·s )
T
La temperatura T está en Kelvin.
