UNIVERSIDAD NACIONAL DE INGENIERÍA

Facultad de Ingeniería de Petróleo, Gas Natural y Petroquímica

Escuela Profesional de Ingeniería Petroquímica
PRÁCTICA DIRIGIDA N° 02
CURSO: Diseño de Reactores (PQ417A)

Profesor : Ruiz Carlos

Alumno : Código:
Ciclo Académico : 2020-2
Día y hora : 03 de diciembre del 2020 11:00 – 13:00

PREGUNTAS:

1) La cinética de la descomposición del éter en fase gaseosa dada por ( CH3 )2 O(g) CH4(g) +H2(g) +CO(g) se analiza en

un cilindro con pistón móvil a 90ºC en el que se coloca una mezcla de 30% de éter, 20% de metano, 10% de
hidrógeno, 10% de monóxido de carbono y 30% molar de inertes, determinando una presión constante en el
reactor de 10 atm. A la temperatura de operación se determina que la constante de equilibrio K es 1,173 y k1 es
0,0976.

a) Determine la expresión
dx A
dt
( )
= f x A , ,k1 ,K,C A0 para el éter metílico.

b) Mediante el método de Euler determine el aumento de volumen después de 2 minutos t = 1 min .

Solución
Para la reacción

( CH3 )2 O(g) CH4(g) +H2(g) +CO(g)

Se obtiene la ecuación de velocidad

1 dNA N NNN
− = k1CA − k2CBCC CD = k1 A − k2 B C3 D
V dt V V
Y en base al balance de masa:

dx A  C2A (B + x A )( C + x A )(D + x A ) 
a) = k 1 ( 1 − x A ) − 0 
( 1 + x A )
2
dt  K 

2 1
Donde B = ;C = D = .
3 3
Para  :
( CH3 )2 O(g) CH4(g) H2(g) CO(g) Inertes
N x =0 3 2 1 1 3 = 10
N x =1 0 5 4 4 3 = 16

N x =1 − N x = 0 16 − 10
De  = = = 0,6
N x =0 10

Formato 02 FIP-P2
 UNIVERSIDAD NACIONAL DE INGENIERÍA
Facultad de Ingeniería de Petróleo, Gas Natural y Petroquímica
Escuela Profesional de Ingeniería Petroquímica
Para C A 0 :

3
P
PA0
10 0,3  10
C A0 = = = = 0,1008M
RT RT 0,082  363
Reemplazando:

 2  1  1 
 0,10082  + x A  + x A  + x A  
dx A
= 0,0976 (1 − x A ) −  3  3  3 
 (1 + 0,6x A ) 
2
dt 1,173
 
2  1  1 
 + x A  + x A  + x A 
= 0,0976 (1 − x A ) − 8,452  10 −4   3  3 
dx A 3
(1 + 0,6x A )
2
dt

b) Aplicando el método de Euler

 2  1  1 
  + x Ai  + x Ai  + x Ai  
x Ai+1 = x Ai + t
dxA
( )
= x Ai + t 0,0976 1 − x Ai − 8,452  10 −4 
3  3  3 
 ( ) 
2
dt x Ai ,ti 1 + 0,6x Ai
 
 2  1  1 
  + x Ai  + x Ai  + x Ai  
x Ai+1 ( )
= x Ai + 0,0976 1 − x Ai − 8,452  10 −4 
3  3  3 
 ( ) 
2
1 + 0,6x Ai
 

Obteniendo los siguientes resultados:

t xA f'(xA)
0 0.0000 0.0970
1 0.0970 0.0871
2 0.1841
Para el volumen:
V = V0 (1 + x A )
V V
−1 = = x A = 0,6  0,1841 = 0,1104 = 11,04%
V0 V0

Formato 02 FIP-P2
 UNIVERSIDAD NACIONAL DE INGENIERÍA
Facultad de Ingeniería de Petróleo, Gas Natural y Petroquímica
Escuela Profesional de Ingeniería Petroquímica
2) La reacción representada por la ecuación química 2NO(g) + O2(g) → 2NO2(g) tiene una cinética que solo depende del

monóxido de nitrógeno. En laboratorio se lleva a cabo esta reacción utilizando una mezcla de 40% de NO y 60%
de aire “enriquecido” que contiene 50% de oxígeno y 50% de nitrógeno; a la temperatura de operación el
nitrógeno no reacciona por lo que se considera un gas inerte. A una temperatura de 80ºC, el avance de la reacción
se verifica con un manómetro, cuya lectura se indica en la tabla:
t, min 0 5 10 15 20
Pman, atm 1,00 0,89 0,83 0,79 0,76

a) Mediante un balance de masa, determine la relación entre la presión y la presión parcial de monóxido de
nitrógeno.
b) Si se considera que la cinética solo depende del monóxido de nitrógeno, determine el valor de k y el orden de
la reacción por el método integral (compruebe 1er y 2do orden).
Solución
Del balance de masa:
NA = NA0 − NAx
NAx
NB = NB0 − NBx = NB0 −
2
NC = NC0 + NCx = NC0 + NAx

Sumando los términos:

NAx NA0 − NA 0,4N0 − NA RT
N = N0 − = N0 − = N0 − = 0,8N0 + 0,5NA 
2 2 2 V
P = 0,8P0 + 0,5PA
PA = 2P − 1,6P0

a
*Deduzca la fórmula PA = PA0 − (P − P0 )
n
En la tabla:

t, min 0 5 10 15 20
Pman, atm 1,00 0,89 0,83 0,79 0,76
Pabs, atm 2,00 1,89 1,83 1,79 1,76
PA, atm 0,8 0,58 0,46 0,38 0,32

Con la presión parcial de A obtenida se aplica el método integral:

Formato 02 FIP-P2
 UNIVERSIDAD NACIONAL DE INGENIERÍA
Facultad de Ingeniería de Petróleo, Gas Natural y Petroquímica
Escuela Profesional de Ingeniería Petroquímica
Para el orden 1 Para el orden 2
dC A dC A
− = kC A − = kC2A
dt dt
2
dP P dPA P 
− A =k A − =k A 
RTdt RT RTdt  RT 
dPA
− = kPA dPA P2
dt − =k A
dt RT
 PA  1 1  k 
ln  0  = kt − =  t
 PA  PA PA0  RT 

 PA  1 1
t PA ln  0  −
 PA  PA PA0
0 0,80 0,0000 0,00000
5 0,58 0,3216 0,4741
10 0,46 0,5534 0,9239
15 0,38 0,7444 1,3816
20 0,32 0,9163 1,8750

2.00000

1.80000

y = 0.0931x
1.60000
R² = 0.9999
1.40000

1.20000
n=1

1.00000 n=2
Lineal (n = 1)
0.80000
Lineal (n = 2)
0.60000

0.40000
y = 0.0488x
R² = 0.9922
0.20000

0.00000
0 5 10 15 20 25

k
Por lo que el orden de la reacción es 2. Así = 0,0931 y k = 2,695 .
RT

Formato 02 FIP-P2
 UNIVERSIDAD NACIONAL DE INGENIERÍA
Facultad de Ingeniería de Petróleo, Gas Natural y Petroquímica
Escuela Profesional de Ingeniería Petroquímica
3) La reacción en fase gaseosa 2N2O5(g) → 2N2O4(g) + O2(g) es elemental. Cuando se introduce 20 mmol de pentóxido

de dinitrógeno puro a 100°C en un reactor de volumen constante de 1 litro, la presión se eleva un 20% del valor
inicial en 20 minutos. Determine el número de moles obtenido de tetróxido de dinitrógeno (N 2O4) en este tiempo
y compare este valor con el obtenido en los casos:
a) En que el reactor sea de volumen variable (con volumen inicial de 1 litro) y la alimentación contenga 20 mmol
de N2O5.
b) En que el reactor sea de volumen variable (con volumen inicial de 1 litro) y la alimentación contenga 20 mmol
de N2O5 y 30% de inertes.
Solución
a
PA = PA0 − (P − P0 )
n
2
PA = PA0 − (P − P0 )
1
PA = 3P0 − 2P
1 1 k
− = t
PA PA0 RT

Para t = 20 min

P = 1,2P0
PA = 3P0 − 2P = 3P0 − 2 (1,2P0 ) = 0,6P0

1 1 k
− = t
0,6P0 P0 RT
2 k
= t
3P0 RT
2  RT   1  2  V   1  2  1   1  5
k =    =    =   =
3  P0   t  3  N0   t  3  0,02   20  3

 nB = 8,0 mmol
a) A volumen variable

1 dNA N2
− = kC2A = k A2
V dt V
N2A (1 − x A ) (1 − x A )
2 2
dN N2 dx dx
− A = k A  NA0 A = k 0  A = kC A0
dt V dt V0 (1 + x A ) dt (1 + x A )
(1 + x A ) dx = kC dt
(1 − x A )
2 A A 0

xA
 ln (1 − x A ) + (1 +  ) = kC A0 t
1 − xA

Formato 02 FIP-P2
 UNIVERSIDAD NACIONAL DE INGENIERÍA
Facultad de Ingeniería de Petróleo, Gas Natural y Petroquímica
Escuela Profesional de Ingeniería Petroquímica

N x =1 − N x = 0 3−2
con  = = = 0,5
N x =0 2

 x  2
0,5ln (1 − x A ) + (1,5 )  A  =  x A = 0,3755
 1 − xA  3
NB = NA0 x A = 7,51 mmol

N x =1 − N x = 0 ( 7 + 3,5 + 3) − ( 7 + 3) = 0,35
b) con  = =
N x =0 ( 7 + 3)
 x  2
0,35ln (1 − x A ) + (1,35)  A  =  x A = 0,3822
 1 − xA  3
NB = NA0 x A = 7,64 mmol

4) La descomposición del pentóxido de dinitrógeno en fase gaseosa 2N2O5(g) → 2N2O4(g) + O2(g) tiene lugar en un

reactor con pistón móvil a presión constante de 1 atm y una temperatura constante de 120°C. Para la
determinación de la cinética de la reacción se alimenta una mezcla de 70% de pentóxido de dinitrógeno y 30% de
inertes y se da inicio a la reacción obteniéndose los siguientes resultados

t (min) 0 2 4 6 8 10
V (l) 1,000 1,136 1,193 1,225 1,246 1,260

a) En base al balance de masa, deduzca la el valor del coeficiente de expansión volumétrica (ε ).

b) Determine la ecuación cinética por el método integral.

Solución
a) En base a la fórmula
2N2O5(g) 2N2O4(g) O2(g) Inertes
N x =0 7 0 0 3 = 10
N x =1 0 7 3,5 3 = 13,5

N x =1 − N x = 0 13,5 − 10
De  = = = 0,35
N x =0 10

Con el balance de masa:

Formato 02 FIP-P2
 UNIVERSIDAD NACIONAL DE INGENIERÍA
Facultad de Ingeniería de Petróleo, Gas Natural y Petroquímica
Escuela Profesional de Ingeniería Petroquímica
NA = NA0 − NA0 x A
NB = NB0 + NA0 x A
xA
NC = NC0 + NA0
2
NI = NI
_______________
xA   NA  x 
N = N0 + NA0 = N0 1 +  0  A 
2   N0  2 
  NA  x  RT
N = N0 1 +  0  A  
  N0  2  P
  NA  x    7 x 
V = V0 1 +  0  A  = V0 1 +   A  = V0 (1 + 0,35x A )
  N0  2    10  2 

b) Para orden 1
1 dNA N
− = kC A = k A
V dt V
dN dx dx
− A = kNA  NA0 A = kNA0 (1 − x A )  A = k (1 − x A )
dt dt dt
− ln (1 − x A ) = kt

Para orden 2

1 dNA N2
− = kC2A = k A2
V dt V
N2A0 (1 − x A ) (1 − x A )
2 2
dNA N2A dx A dx
− = k  NA0 =k  A = kC A0
dt V dt V0 (1 + x A ) dt (1 + x A )
(1 + x A ) dx = kC dt   ln 1 − x + 1 +   x A  = kC t
( A ) ( ) 
(1 − x A )  1 − xA 
2 A A 0 A 0

V
−1
V0
En ambos casos, de V = V0 (1 + x A ) se obtiene x A =


 x 
t V xA − ln (1 − x A )  ln (1 − x A ) + (1 +  )  A 
 1 − xA 
0 1.000 0.0000 0.0000 0.0000
2 1.136 0.3900 0.4942 0.6900
4 1.193 0.5516 0.8021 1.3800
6 1.225 0.6429 1.0297 2.0700
8 1.246 0.7023 1.2115 2.7600
10 1.260 0.7442 1.3633 3.4500

Formato 02 FIP-P2
 UNIVERSIDAD NACIONAL DE INGENIERÍA
Facultad de Ingeniería de Petróleo, Gas Natural y Petroquímica
Escuela Profesional de Ingeniería Petroquímica
Aplicando la regresión lineal:

4.0000

3.5000
y = 0.345x
R² = 1
3.0000

2.5000
n=1

2.0000 n=2
Lineal (n = 1)
1.5000 Lineal (n = 2)
y = 0.1532x
R² = 0.9788
1.0000

0.5000

0.0000
0 2 4 6 8 10 12

PA0 0,7P0 0,7 (1)

El orden de la reacción es 2 y kC A0 = 0,345 donde CA0 = = = = 0,02172 M . Así k = 15,88 .
RT RT 0,082  393

5) La reacción A ⎯⎯
k
→B ⎯⎯
k
→ C se lleva a cabo en un reactor discontinuo en el cual, la concentración del producto B
asciende hasta llegar a los 5 minutos al valor de 1,2 M; a partir de este tiempo, esta concentración disminuye. En
base a estos datos, determine los valores numéricos de k y la concentración inicial de A.
Solución

Para la reacción en serie A ⎯⎯

k1
→B ⎯⎯
k2
→C
dCA dCB dCC
− = k 1C A − = k2CB − k1CA − = −k2CB
dt dt dt
C A = C A 0 e − k1 t
dCB dC
− = k 2 C B − k 1 C A  − B = k 2 C B − k 1 C A 0 e − k1 t
dt dt
dCB
+ k 2 C B = k 1 C A 0 e − k1 t
dt
Aplicando el factor de integración

FI = e 
k2 dt
= ek2 t
dCB
ek2t + k2ek2tCB = k1CA0 e−k1tek2t
dt

Formato 02 FIP-P2
 UNIVERSIDAD NACIONAL DE INGENIERÍA
Facultad de Ingeniería de Petróleo, Gas Natural y Petroquímica
Escuela Profesional de Ingeniería Petroquímica
d ( ek2 tCB )
= k 1 C A 0 e ( − k1 + k 2 ) t
dt
Para el caso de k1 = k2

d ( ektCB )
= kCA0
dt

d ( ekt CB ) = kC A0 dt

d ( ektCB ) = kC A0  dt
CB ,t t
CB0 ,0 0

ekt CB = kC A0 t  CB = kC A0 e−kt t

*Desarrolle el caso de k1  k2

Graficando CB
Por lo que el tiempo en el que se obtiene la
Cb
concentración máxima de B es:
1.2
1
1
0.8 t=
k
0.6
0.4
Y la concentración máxima es:
0.2
0
1
0 5 10 15 20 −k  
1
CBmax = kC A0 e−kt t = kC A0 e k
  = C A0 e
−1

k 
Se observa que la concentración de este reactivo C A0
CBmax =
e
presenta un valor máximo. Para hallar ese valor máximo
y el tiempo en el que se obtiene: De los datos

dCB 1
=0 k = = 0,2 min−1
dt t
(
d kC A0 e−kt t ) = d(e t) = 0 −kt C A0 = eCBmax = 2,718  1,2 = 3,262 M
dt dt
−ke−kt t + e−kt = 0
e−kt ( −kt + 1) = 0

Formato 02 FIP-P2