UNIVERSIDAD NACIONAL

JOSÉ FAUSTINO SÁNCHEZ CARRIÓN

ESCUELA ACADÉMICO PROFESIONAL

INGENIERÍA QUÍMICA

TRABAJO DE INVESTIGACIÓN

EJERCICIOS DE REACTORES TANQUE AGITADO (CSTR), FLUJO PISTÓN

(PFR) E INTERMITENTE ISOTÉRMICO

CURSO:

CINÉTICA Y DISEÑO DE REACTORES QUÍMICOS

PRESENTADO POR:

GONZÁLEZ MUÑOZ, CAMILA PAOLA

DOCENTE:

Ing. COCA RAMIREZ VICTOR RAUL

HUACHO – 2023
Solución.
Ec. Diseño para un Reactor Intermitente Isotérmico. Para un reactivo en fase líquida:
d XA
(−r A )V =N Ao
dt
Ordenamos los términos e integramos.
t XA
d XA
∫ dt =N Ao∫ (−r
0 0 A )V

Consideramos que su densidad permanece constante y ε =0, obtenemos la siguiente

ecuación que nos indica el tiempo necesario para que se alcance la conversión de X A :
XA XA
N Ao d X A d XA
τ= ∫
V 0 (−r A )
=C Ao∫ =13 min .
0 (−r A )

Para un reactor de flujo pistón (PFR):

XA
d XA
τ =C Ao ∫
0 (−r A )

Observamos que, debido a las condiciones del reactivo en fase líquida, su densidad
permanece constante y ε =0, la ecuación de diseño para determinar el tiempo necesario
para alcanzar la conversión es igual que la del reactor intermitente isotérmico.
Hallando el tiempo espacial τ :
XA
d XA
τ =C Ao ∫ =13 min .
0 (−r A )

Hallando la velocidad espacial s:

1 1 −1
s= = =0.07692≈ 0.08 min
τ 13 min .

Para un reactor de tanque agitado (CSTR):

V C Ao X A
τ= =
V o (−r A )

Hallando el tiempo espacial τ :

Debido a que desconocemos la cinética de la reacción, no podemos calcular τ
Hallando la velocidad espacial s:
Debido a que desconocemos la cinética de la reacción, no podemos calcular s

Solución:
Ec. Diseño para un Reactor CSTR:
F Ao X A
V CSTR =
(−r A )
Primero, para el volumen actual (V1) del reactor, tenemos:
F Ao X A 1
V 1= … … … …( I )
(−r A )
Y, sabemos que:
F Ao=ν o C Ao

C A=C Ao (1− X A )

Sustituimos en (I), despejando las variables obtenemos:

ν o C Ao X A 1
V 1=
¿¿
ν o C Ao X A 1
V 1=
( k C ) ( 1−X
3 3

A 1)
2 2
Ao

k V1 C Ao X A 1
= … … … (II )
νo
(C ) ( 1−X
3 3

Ao
2
A1 ) 2

Ahora, la ecuación de diseño para el segundo reactor (V2), consideramos la siguiente

relación:
V 2=2 V 1
Del enunciado sabemos que se usará la misma concentración inicial y flujo volumétrico
cuando se escale el reactor. Por lo tanto, la ecuación de diseño para el segundo reactor
es la siguiente:
ν o C Ao X A 2
2 V 1=
( k C ) ( 1−X
3 3

A2)
2 2
Ao

Despejamos:
k 2V 1 C Ao X A 2
= … … …(III )
νo
( C )( 1−X
3 3

Ao
2
A2 ) 2

De acuerdo a la relación, igualamos (II) y (III) para tener el siguiente resultado:

C Ao X A 1
kV1
( )
3 3
νo C Ao2 ( 1−X A 1 ) 2
=
k2V1 C Ao X A 2
νo
( )
3 3
2
C Ao ( 1−X A 2 ) 2

Reemplazamos los valores que el enunciado nos da: C Ao=10 mol /L, y X A 1=0.7 , luego
simplificamos la ecuación:
mol
10 ∙0.7
C Ao X A 1 L

( )
3
kV1 mol 2
3

( )
3 3
10 ( 1−0.7 ) 2
νo C Ao2 ( 1−X A 1 ) 2 1 L
= → =
k2V1 C Ao X A 2 2 mol
10 ∙ X A2
L
νo
( )
3 3
C Ao2 ( 1−X A 2 ) 2 mol 32
( )
3
10
L
( 1−X A 2 ) 2

( )
3 3
mol 2
1.3472 10 ( 1−X A 2 ) 2
1 L
=
2 10 ∙ X A 2

X A 2=0.7944 ≈ 0.79

Luego, tabulamos:
V2 XA2
0.0625 0.19
0.0125 0.31
 0.25 0.44
0.5 0.58
1 0.7
2 0.79
4 0.86
8 0.91
16 0.94
32 0.96

En conjunto a Polymath, se muestra la figura 1, que muestra cómo cambia la conversión

con el cambio del volumen del reactor, manteniendo la información brindada en el
problema (flujos volumétricos de entrada, concentración de entrada, V 1 y X A 1).

Fig 1. Gráfica que explica la variación de la conversión con respecto al factor de

escalamiento del reactor
Cómo se ve en la gráfica, entre mayor sea el volumen del reactor, mayor será la
conversión. Esto se debe a que, al mantener el flujo y concentraciones de entrada
iguales que el reactor a escalar, el tiempo de residencia t aumenta, favoreciendo la
reacción.
Solución.
Sistema de densidad constante porque es líquido.
k 1 Cℜ X Ae 0.04
K= = = = =4 → X Ae =0.8
k 2 C Ae 1−X Ae 0.01

k 1 (1−X Ae )
−r A =k 1 C Ao ( 1− X A ) −k 2 C Ao X A =k 1 C Ao ( 1− X A ) − C Ao X A
X Ae

{ [
−r A =k 1 C Ao 1−X A 1+
1−X Ae
X Ae ]}
=k 1 C Ao ( 1−1.25 X A )

XA XA
V d XA d XA
τ = =C Ao ∫ =C Ao∫
0 (k 1 C Ao ( 1−1.25 X A ) )
Vo 0 (−r A )

k1 τ
X Ae (
=−ln 1−
XA
X Ae )
0.04 ( 2000
100 )
0.8
=−ln 1−(X
0.8
A
)=1
Hallamos la conversión para el reactivo A:

X A =0.8 ( 1−e−1 )=0.506

Solución.
Sistema de densidad constante porque es líquido
k 1 (1−X Ae )
−r A =k 1 C Ao ( 1− X A ) −k 2 C Ao X A =k 1 C Ao ( 1− X A ) − C Ao X A
X Ae

{ [
−r A =k 1 C Ao 1−X A 1+
1−X Ae
X Ae ]}
X Ae =0.8

V 2000 C Ao X A C Ao X A XA
τ= = = = =

{ [ ]}
V o 100 −r A 1−X Ae 0.04 (1−1.25 X A )
k 1 C Ao 1− X A 1+
X Ae

Hallamos la conversión para el reactivo A:

0.8−X A =X A → X A =0.4

Solución.
Sistema de densidad variable porque la alimentación F total, varía, por lo tanto, el flujo
volumétrico también presentará cambios.
 1
k τ=( 1+ ε A X A ) ln +ε X
1− X A A A

mol
100
F min L
ν o= Ao = =50
C Ao mol min
2
L
V X −V X 2.5−1
ε A= A=1 A=0
= ( 1 )=1.5
VX A=0
1

10 ( 2250 )=2.5 ln 1−X1 −1.5 X A

A

4.4=2.5 ln (1− X A )−1.5 X A=f ( X A )

Tabulamos los datos obtenidos:

XA 0.7 0.8 0.75
fXA 4.05 5.22 4.59

Con ayuda de Polymath obtenemos la siguiente gráfica:

X A =0.73
SOLUCIÓN

Sistema de densidad constante porque es un líquido

C A 0=0 ,1 mol /L

C B 0=0 , 4 mol / L

V 0=400 L/min
X A =0 , 9

CA 0 X A
τ m=
−r A

C A=C A 0 ( 1− X A )

C B=C A 0 ( M B −b /a X A )

C B 0 200
M B= = =2
C A 0 100
b=a=1
2
−r A =200 ( 0 , 1 ) (1−X A )(2−X A )

0 ,1(0 ,9)
τ m= 2
=49 , 9 min
200 ( 0 ,1 ) (1−X A )(2− X A )
3
V m =τ m . V 0=49 ,9 ( 400 )=19960 L ≈ 20 m

SOLUCIÓN:

Sistema de densidad variable porque es gaseoso y como varía F total, el flujo volumétrico
también variará
 C A 0 ( 1−X A ) C A 0 ( 1− X A ) C A 0 ( 1−X A )
C Af = = =

( ) P)
( ( 1+ε A )
fe T P 0
( 1+ ε A ) T
0

ε A= ( 1−3
3 ) ( 1 )=
2
3
1−X A
330=660 X A =0 , 75
2 
1− X A
3
0 ,75
d X A C A0 660
τ p=C A 0 ∫ = X A= 0 ,75=9 , 17 min
0 −r A −r A 54

( )
FA0 540
V =τ p .V 0=τ p =9 ,17 =7 ,5 L
CA0 660

SOLUCIÓN:

Sistema de densidad variable porque es gaseoso y como varía F total durante el transcurso de la
reacción, el flujo volumétrico varía
τ V C A0 X A
m=¿ = ¿
V0 −r A

ε A= ( 3−1
1 )
( 0 , 5 )=1

1−X A 1−X A
−r A =0 , 6 C A=0 , 6 C A 0 =0 ,6 C A 0
1+ ε A X A 1+ X A
τ V 1000 C A0 X A (1+ X A ) X A (1+ X A )
m=¿ = = = ¿
V 0 180 0 , 6 C A 0(1−X A ) 0 ,6 (1− X A )

2
3 X A +13 X A−10=0

−13 ± √ 169+ 40(3)

X A= =0 , 67
2(3)
SOLUCIÓN:

El sistema es líquido, así que es de densidad constante y ∴ ε A =0

τ=
V
=
2L
V 0 0 , 5 L/min
=4 min
60 s
1 min(=240 s )
2 2 XA
k τ P C A 0=2 ε A ( 1+ ε A ) ln ( 1−X A ) +ε A X A +(ε A +1)
1− X A
XA
k τ P C A 0=
1−X A
k τP C A0 0 ,05 (240)(1)
X A= = =0 , 92
1+k τ P C A 0 1+0 , 05(240)(1)

SOLUCIÓN:
El sistema es de densidad constante porque es líquido
C A 0−C A C A 0−C A (100−C A )V 0
τ m=  −r A = ∴ −r A =
−r A τm V
v (L/min) 1 6 24
CA 4 20 50
(mmol/L)
−r A 96 480 1200
ln(−r A ¿=lnk+ n ln C A
n
−r A =k C A 
 1
ln 1200−ln 96 CA 50 −1
n= =1 ∴ k= = =0,0417 min
ln 50−ln 4 −r A 1200
mmol
−r A =0,0417 C A
Lmin

Problema 5.21

Solución
Como la reacción se encuentra en fase liquida
Se aplicará la siguiente expresión
CA 1

d CA
τ =∫
CA 0
−r A

Esto se puede asemejar a la expresión de

 { }
CA f −1
d CA ∆ CA 1
1

1 1
τ =∫ ≈ + +∑
C −r A
A0
2 −r A 0 −r Af i=1 −r Ai

Además
Como se tiene que hallar el tiempo entre CA = 0,3 hasta CA = 1,3 mol/L
Se tendrá que

{ }
1.3 f −1
d CA ∆ CA 1 1 1
τ =∫ ≈ + +∑
0.3 −r A 2 −r A 0 −r Af i=1 −r Ai

Donde para una mejor resolución de los datos completaremos la escala desde 0.3 a 1.3
Siendo la grafica

Y los datos aproximados faltantes son

Ca (mol/L) 0.8 0.9 1 1.1 1.2 1.3
-ra (mol/min 0.06 0.053 0.05 0.0475 0.046 0.045
Por lo cual si empleamos los datos de las dos tablas en la expresión se tendrá

{ [ ]}
1.3
d C A 1.3−0.3 1 1 1 1 1 1 1 1 1 1 1
τ =∫ ≈ + + + + + + + + + +
0.3 −r A 2 0.5 0.045 0.6 0.5 0.25 0.1 0.06 0.053 0.05 0.0475 0.046

Efectuando
1.3
d CA
τ =∫ ≈ 70.107 min
0.3 −r A

Problema 5.22
Solución
Como la reacción está en fase liquida
Diremos que
C A=C A 0 ( 1− X A )

Siendo así
C A=1.5 ( 1−0.8 )=0.3 mol/ L

Ahora bien
Para determinar el tamaño se empleará en base al tiempo espacial
CA 1

d CA
τ =∫
CA 0
−r A

Siendo así que se tendrá que integrar desde 0.3 hasta 1.5 mol/L
1.5
d CA
τ =∫
0.3 −r A

Pudiendo expresar esto de la siguiente manera

1.5
d C A 1.3 d C A 1.5 d C A
τ =∫ =∫ +∫
0.3 −r A 0.3 −r A 1.3 −r A

Del cual se tiene que

1.3
d CA
τ =∫ ≈ 70.107 min
0.3 −r A

Por lo cual se vende determinar desde 1.3 hasta 1.5 de la siguiente manera
Tomando en cuenta el grafico anterior completamos los datos faltantes
Ca (mol/L) 1.3 1.4 1.5
-ra (mol/min 0.045 0.0445 0.044
Reemplazamos los valores para este, en la siguiente fórmula

{ ( )}
1.5
d C A ∆ CA 1 1 1
τ =∫ ≈ + +2.
1.3 −r A 2 −r A 0 −r Af −r Ai

Reemplazando los valores

{ )}
1.5

(
d C A 1.5−1.3 1 1 1
τ =∫ ≈ + +2.
1.3 −r A 2 0.045 0.0455 0.044

Efectuando
 1.5
dCA
τ =∫ ≈ 28.989 min
1.3 −r A

Por lo tanto
τ total=70.107 min+28.989 min=99.096 min

Problema 5.24

Solución
Como se encuentra en fase gaseosa
Aplicaremos la siguiente expresión

V C A 0∗X A
τ= =
V0 −r A

Despejando ra
V 0∗C A 0∗X A F A 0∗X A F A 0∗X A
−r A = = =
V V 0.1 L
Por otor lado
CA
1−
CA 0 V XA=1−V XA=0 ( 5−1 )
X A= ; donde ε A= = =4
CA V XA =0 1
1+ε A
CA 0

Por lo cual crearemos la siguiente tabla

Fao Ca Xa -ra
300 16
 1000 30
3000 50
5000 60
Ahora
Sea la reacción de velocidad
n
−r A =k . C A

Linealizando se tendría que

ln (−r A ) =ln
⏟ k + n. lnC A
⏟
y b +mx

Así que bien

La nueva tabla y respectiva grafica serian

Del cual se obtiene que

n=1.0113 ≈ 1
ln k =4.5648 → k =96.04
Por lo tanto
−r A =96. C A

Problemas 5.25
Solución
cómo nos dan de dato el tiempo de permanencia, entonces consideramos que
t=τ
Del cual
CA −C A CA −C A
τ= entrada salida
→−r A= entrada salida

−r A τ

Así que dicho esto

Los datos tabulados son

De estos datos completamos los datos faltantes

Ca (mol/L) 0.2 0.3 0.4 0.5 0.6 0.7 0.8
-ra (mol/min 0.0056 0.0011 0.0021 0.0034 0.0042 0.0046 0.0048
Donde el tiempo se calculará mediante la siguiente expresión
Que va desde 0.2 hasta 0.8

{ }
0.8 f −1
d CA ∆ CA 1 1 1
τ =∫ ≈ + +2 ∑
0.2 −r A 2 −r A 0 −r Af i=1 −r Ai

Que empleando los datos de las tablas se obtiene

{ [ ]}
0.8
d C A 0.8−0.2 1 1 1 1 1 1 1
τ =∫ ≈ + +2 + + + +
0.2 −r A 2 0.0056 0.0048 0.0011 0.0021 0.0034 0.0042 0.0046

Resolviendo
 0.8
d CA
τ =∫ ≈ 1879.145 s
0.2 −r A

Problema 5.26

Como será una reacción completa

Tan solo aplicaremos
CA −C A
τ= entrada salida

−r A

Siendo así que

CA −C A 0.8−0.2
τ= entrada salida
= =1071.42 s ≈ 0.298 h
−r A 0.00056

Problema 5.28

Solución
Supongamos que la reacción es
n
−r A =k . C A

Si n = 1 (primer orden)
−d C A
−r A = =k .C A
dt
 dCA
=−k dt
CA

Integrando
CA t
dCA
∫ CA
=−k ∫ dt → ln C A −ln C A 0=−k . t
C A0 0

ln C A−ln C A 0=−k . t → ln
( )
CA
C A0
=−k . t

Del cual al ver que se tiene datos de presión

PA
C A=
RT
Por lo tanto

()
PA
RT
ln =−k . t
PA0
RT

Linealizando
PA
ln
P
⏟ bx
a 0
( )
⏟.t
=−k

Por lo cual con los datos realizaremos una grafica

Del cual se obtuvo que

−1
−k =−0.0096 → k =0.0096 s
Del cual se tiene que
−r A =0.0096 C A

Con esto
Aplicaremos la siguiente expresión para el cálculo del volumen
k∗V
k .τ= =−ln ( 1−X A )depejando
V0
 −V 0
V= . ln ( 1−X A )
k
Donde
F A 0 F A 0 F A 0 .T 0 . R
V 0= = =
C A 0 P AO PA0
RTo

Reacomodando los valores

(100 mol/h)(100+273)( 0.082)
V 0=
1atm (0.8)

L
V 0=3823.25 ∗
1h
h 3600 s (
=1.06
L
s )
Por lo tanto
−V 0 −1.06∗ln ( 1−0.95 )
V= . ln ( 1−X A )=
k 0.0096
V =−(−330.778 )=330.778 L