UNIVERSIDAD DE ANTIOQUIA

FACULTAD DE INGENIERIA
DEPARTAMENTO DE INGENIERIA QUÍMICA
CURSO: IQU-458 INGENIERÍA DE LAS REACCIONES QUÍMICAS
Tarea 2

Barrientos Pérez Andrea, Bravo Caballero Yulieth y García Álvarez Tatiana.

1. La alquilación de tolueno con acetileno en presencia de ácido sulfúrico

(Reacción 1) se lleva a cabo en un reactor semicontinuo. Se cargan 6910,5 kg de
tolueno en cada lote, junto con la cantidad requerida de ácido sulfúrico, y el
acetileno se alimenta continuamente al reactor a presión. Bajo intensa agitación,
se puede suponer que el líquido siempre está saturado con acetileno y que el
tolueno se consume en una reacción de pseudo primer orden respecto al tolueno
con una constante de velocidad de 1.1x10-4 s-1 a 8 ° C.

2 C 6 H 5 CH 3+C 2 H 2 →(C 6 H 4 CH 3)2 CHCH 3 Reacción 1

El reactor se apaga durante un período de 15 minutos entre lotes. El volumen de

la mezcla inicial es de 4 m 3, el flujo volumétrico es de 0,5 L / s y el flujo molar de
acetileno es de 10 mol / s. El reactor tiene una capacidad de 10 m 3 (incluido un
20% de espacio libre).

a. ¿Cuál es la conversión de tolueno lograda en un reactor semibatch?

Grafique la conversión, el volumen y las moles de todas las especies en
función del tiempo.

Consideraciones:
 Semi-lotes con respecto al componente A (Tolueno)
 Se agrega controladamente el componente B (Acetileno)
 Reacción en fase líquida (ρ=cte)
 Proceso isotérmico (8°C)
 Intensa agitación
 Reacción de pseudo primer orden respecto al componente A

1
  El volumen es variable

La reacción de alquilación de tolueno con acetileno en presencia de ácido sulfúrico

(Reacción 1) se puede simplificar de la siguente forma:
2 A + B→ C Reacción 2
Donde:
A: Tolueno
B: Acetileno
C: 1-2-ditoliletano (DTE)
A continuación, se divide la reacción 2 por el coeficiente estequiométrico de A (γ A
= 2)
1 1
A+ B → C Reacción 3
2 2
Se realiza el balance de moles para cada una de las especies:
Balance de moles para A
d NA
F A 0−F A +r A∗V = Ec .1
dt
Donde: FAo = Velocidad de flujo molar de A a la entrada [mol/s]
FA = Velocidad de flujo molar de A a la salida [mol/s]
rA = Velocidad de formación de A [mol/L*s]
V = Volumen del reactor (variable) [L]
dNA/dt = Velocidad de acumulación de A [mol/s]

Debido a que no hay entrada ni salida de A se obtiene la ecuación 2.

d NA
r A∗V = Ec .2
dt
Balance de moles para B
d NB
F B 0−F B + r B∗V = Ec . 3
dt
Donde: FBo = Velocidad de flujo molar de B a la entrada [mol/s]

2
FB = Velocidad de flujo molar de B a la salida [mol/s]
rB = Velocidad de formación de B [mol/L*s]
dNB/dt = Velocidad de acumulación de B [mol/s]

Se agrega controladamente B, pero no hay salida de esta especie (FB = 0), por lo
tanto, se obtiene la ecuación 4.
d NB
F B 0 +r B∗V = Ec . 4
dt
Balance de moles para C
d NC
F C 0−FC +r C∗V = Ec . 5
dt

Donde: FCo = Velocidad de flujo molar de C a la entrada [mol/s]

FC = Velocidad de flujo molar de C a la salida [mol/s]
rC = Velocidad de formación de C [mol/L*s]
dNC/dt = Velocidad de acumulación de C [mol/s]

Debido a que no hay entrada ni salida de C se obtiene la ecuación 6.

d NC
r C ∗V = Ec . 6
dt
Sabiendo que el tolueno se consume en una reacción de pseudo primer orden
respecto al tolueno con una constante de velocidad de k A = 1.1x10-4 s-1 se puede
obtener la ecuación para la velocidad de reacción r A.
r A =−k A∗C A E c .7

A partir de la relación de velocidades se determina r B y rC

r A rB rC
= =
γ A γ B γC

Donde: 𝛾 es el coeficiente estequiométrico cada especie.

Así:
−1∗r A
r B=
−2

3
 rA
r B= Ec .8
2
1∗r A
rC =
−2
−r A
rC = Ec . 9
2
Por definición se tiene:
NA
C A= Ec .10
V
NB
C B= Ec . 11
V
NC
C C= Ec . 12
V
Donde: CA = Concentración de A [mol/L]
CB = Concentración de B [mol/L]
CC = Concentración de C [mol/L]
NA = Moles de A [mol]
NB = Moles de A [mol]
NC = Moles de A [mol]

Ahora se debe encontrar una expresión para el volumen debido a que está
variando.
Balance de masa global
dm
ṁ o− ṁ+0= Ec . 13
dt
Donde: ṁ o= Flujo másico que entra [g/s]
ṁ= Flujo másico que sale [g/s]
dm/dt = Variación de la masa respecto al tiempo [g/s]

mi
Teniendo en cuenta la definición de densidad: ρi=
Vi

4
 d ( ρV )
ρ0∗v o −0+0= Ec . 14
dt
Donde: vo = Flujo volumétrico [L/s]
Considerando ρ = ρo se obtiene la ecuación 15.
dV
=v o Ec . 15
dt
Separando variables e integrando desde V o a V y desde to = 0 s hasta t = t se
obtiene la ecuación 16.
V t

∫ dV =v o∗∫ dt
Vo to

V =V o +v o∗t E c .16

Finalmente, la conversión la evaluamos para el compuesto que está dentro del

reactor.
N Ao −N A
X= Ec .17
N A0

Planteando un algoritmo en software POLYMATH se pueden resolver las

ecuaciones diferenciales y determinar la conversión de tolueno lograda en el
reactor semibatch. El algoritmo se muestra en la ilustración 1.
Se debe tener en cuenta que los datos y condiciones iniciales son las siguientes:
mA
N Ao= Ec . 18
PM A

Donde: NAo = Moles inicialea de A [mol]

mA = Masa de A [g]
PMA = Peso molecular de A [g/mol]
6910500 g
N Ao= =75000 mol
g
92 .14
mol
Vo = 4000 L
FBo = 10 mol /s
vo = 0.5 L/s

5
El volumen del reactor es de 10000 L teniendo en cuanta la claridad del 20%, por
lo que la capacidad neta del reactor es 8000 L, entonces cuando se alcance dicha
capacidad t = 8000 s.

Ilustración2.
Ilustración 1.Código
Algoritmo empleado en POLYMATH
POLYMATH

6
 Ilustración 3 Reporte POLYMATH

La conversión de tolueno lograda en el reactor semi-batch es X = 0,5852

Ilustración 3. Conversión de Tolueno en función del tiempo (s)

7
Ilustración 4. Volumen (L) en función del tiempo (s)
 Ilustración 5 Moles de todas las especies (mol) en función del Tiempo (s)

En la Ilustración 5 se evidencia como se esperaba que las moles de A disminuyen

mientras las moles de B y C aumentan conforme transcurre el tiempo.

b. Si necesita producir 17000 kg de 1-2-ditoliletano (DTE), cuántos lotes y cuánto

tiempo toma el proceso si solo tiene un reactor, asumiendo pérdidas del 5% en las
unidades de separación de productos posteriores. ¿Qué volumen debe tener el
reactor para lograr una conversión del 90%?
Consideraciones:
 Semi-lotes con respecto al componente A (tolueno)
 Se agrega controladamente el componente B (acetileno)
 Reacción en fase líquida
 Proceso isotérmico (8°C)
 Intensa agitación
 Reacción de pseudo primer orden respecto al componente A
 El volumen es variable
 Pérdidas del 5% en las unidades de separación de productos posteriores

8
Dado que se requiere producir 17000 kg de 1-2-ditoliletano (DTE) (C) y teniendo
en cuenta pérdidas del 5% en las unidades de separación de productos
posteriores, en cada lote se producirán un poco más kg de C para al final obtener
la cantidad solicitada.
Con una regla de tres se puede determinar la cantidad adicional que se debe
producir.
y = (17000 kg C * 0.05)/ 100 = 850 kg
mC = 17000 kg + 850 kg = 17850 kg C = 17850000 g C
mC: Masa que se debe producir de C [kg]
Con el peso molecular de C se pueden hallar las moles que se producen:
mC
N C= Ec .19
PM C

Donde: NC = Moles de C [mol]

PMC = Peso molecular de C [g/mol]
17850000 g
N Ao= =85000 mol
g
210
mol
Conociendo las moles de C que se deben producir, se puede partir del balance de
moles para C (Ecuación 6) y operando obtener una expresión para el tiempo de
reacción tR.
dNC
r C ∗V =
dt
Separando variables y reemplazando la ecuación 9
d NC
V ∗dt= Ec .20
−r A
2
Reemplazando la ecuación 7 en la ecuacón 20
2∗d N C
V ∗dt = Ec . 21
k A∗C A

Reemplazando la ecuación 10 en la ecuacón 21

9
 2∗d N C
V ∗dt= Ec . 22
k A∗C A
V
Por definición:
N A =N Ao∗( 1−X ) Ec .23

Se reemplaza la ecuación 23 en la ecuación 22 para obtener la ecuación 24

2∗d N C
V ∗dt = Ec . 24
k A∗[ N Ao∗( 1− X ) ]

Integrando la ecuación 24 desde las moles iniciales de C, N Co = 0 (mol) porque

inicialmente no hay nada del compuesto C hasta N C = 85000 (mol) y desde t o = 0 s
a t = t (s)
t NC
2
∫ dt= ∗∫ d N C
k A∗[ N Ao∗(1−X ) ] 0
0

2∗N C
t= Ec .25
k A∗[ N Ao∗(1− X) ]

Reemplazando los valores conocidos y teniendo en cuenta que la conversión X =

0.5852, obtenida anteriormente se calcula t R
2∗85000 mol
t= − 4 −1
=49677 . 0989 s=t R
1. 1 x 10 s ∗[ 75000 mol∗(1−0 .5852)]

Conociendo el tiempo de reacción, se reemplaza el valor obtenido en la expresion

del volumen (Ecuación 16) para conocer el volumen de la reacción.
L
V =4000 L+ 0,5 ∗49677 .0989 s=28838. 5495 L
s
Sabiendo que el reactor tiene una capacidad de 10000 L incluida la claridad del
20%, la capacidad neta del reactor para llevar a cabo la reacción es V R = 8000 L;
para hallar la cantidad de lotes que tienen en el proceso, se divide el volumen de
la reacción por el volumen neto del reactor:
V
Número de lotes= Ec .26
VR

28838. 5495 L
Número de lotes= =3 .6048
8000 L

10
Aproximando, Número de lotes = 4
Continuando, el tiempo que toma producir un lote está comprendido por el tiempo
de apagado entre lotes y el tiempo que toma llegar a la capacidad máxima del
reactor expresado en la ecuación 27.
t 1lote =t apagado +t Capacidad máxima Ec . 27

t 1lote =900 s+ 8000 s=8900 s

Entonces, para conocer el tiempo que toma todo el proceso, se multiplica el

tiempo que dura producir un lote por el número de lotes, así:
t Total =t 1 lote∗Número de lotes Ec . 28

t Total =8900 s∗4=35600 s=9.8889 h

Si se quiere conocer el volumen que debería tener el reactor para alcanzar una
conversión del 90%, se puede determinar siguiendo el procedimiento descrito a
continuación.
Partiendo de la ecuación 2 y la ecuación 23, se expresa el balance de moles para
A, en términos de conversión
−r A∗V dX
= Ec . 29
N Ao dt

Separando variables:
N Ao
Vdt = dX Ec . 30
−r A

Reemplazando la ecuación 7 en la ecuación 30

N Ao
Vdt =
k A∗¿ C A
dX Ec .31 ¿

Reemplazando la ecuación 10 en la ecuación 31

N Ao
Vdt = dX Ec .32
k A∗N A
V
Teniendo en cuenta la ecuacion 23
N Ao
Vdt = dX Ec . 33
k A∗N Ao∗(1− X )
V

11
 1 dX
dt = Ec .34
k A (1− X)

Integrando (por sustitución) la ecuación 34 desde t o = 0 s a t = t (s) y desde X 0 = 0

a X = 0.9 (conversión que se quiere lograr)
t 0,9

∫ dt= k1 ∗∫ (1−X
dX
)
Ec . 35
0 A 0

1
t= ∗ [ ln ( 1−X ) ] ] 0,9 Ec .36
kA 0

Resolviendo la ecuación 36 se obtiene:

2,3036
t= =20932 . 7273 s
1. 1 x 10− 4 s−1
Reemplazando t = 20932.7273 s en la ecuación 16 se obtiene el volumen que
debería tener el reactor para alcanzar una conversión del 90%.
L
V =4000 L+ 0 .5 ∗20932. 7273 s=14466 . 3637 L=14 . 4664 m3
s

La temperatura de reacción se aumenta de 8 °C a 20 °C y la reacción se lleva a

cabo de forma isotérmica. A esta temperatura, la reacción es reversible con una
constante de equilibrio de 10 en las unidades correspondientes. La alimentación
se suspende cuando se alcanza la capacidad máxima, dejando que el sistema
reaccione hasta el equilibrio. La reacción inversa es de segundo orden en tolueno
y de primer orden en DTE. La constante de velocidad de reacción es 7x10-6 L /
(mol.s) para la reacción directa.

c. Grafique la conversión de tolueno y moles de todas las especies en función del

tiempo. ¿Cuánto tiempo toma alcanzar la conversión máxima y cuál es el valor de
la conversión?
Consideraciones:
 Semi-lotes con respecto al componente A (tolueno)
 Se agrega controladamente el componente B (acetileno)
 Reacción en fase líquida
 Proceso isotérmico (20°C)
 Intensa agitación
 Reacción reversible

12
  Reacción de segundo orden respecto al tolueno y de primer orden respecto
al DTE
 El volumen es variable
 k1 =7x10-6 L/mol*s para la reacción directa, Kc = 10 L/mol
 La alimentación se suspende cuando se alcanza la capacidad máxima

Dado que la reacción es reversible 4, se tiene la reacción directa 5 y la reacción

inversa 6.
2 A + B↔ C Reacción 4
2 A + B k 1 C Reacción 5
→

C k−1 2 A+ B Reacción 6
→

Las velocidades de reacción para las reacciones 5 y 6 son respectivamente:

r 1=k 1∗C 2A Ec . 37
r 2=k−1∗C c Ec .38

La velocidad de reacción será entonces

r =r 1−r 2 Ec .39

r =k 1∗C 2A −k −1∗C c Ec . 40

Teniendo en cuenta que por definición:

k1
Kc= Ec . 41
k−1

Cc
(
r =k 1∗ C 2A−
Kc )
Ec . 42

En el equilibrio r = 0
Cc
(
0=k 1∗ C 2A−
Kc )
Ec . 43

2 Cc
C A= Ec . 44
Kc
Para conocer la velocidad de reacción rA se multiplica r por el coeficiente
estequiométrico de A

13
 Cc
(
r A =−2∗k 1∗ C 2A−
Kc )
Ec . 45

A continuación, se divide la reacción 6 por el coeficiente estequiométrico de A

1 1
A+ B → C Reacción 7
2 2
Y a partir de la relación 1 de velocidades se determina r B y rC
rA
r B= Ec . 46
2
−r A
rC = Ec . 47
2
EL proceso se divide en dos etapas; etapa 1: Se alimenta B hasta que la se
alcanza la capacidad máxima del reactor (8000 L) (proceso semi-batch), etapa 2:
Luego de que se suspende la alimentación de B, se deja reaccionar el sistema
hasta que alcance el equilibrio (proceso batch).
Etapa 1
Se realiza el balance de moles para cada una de las especies.
Balance de moles para A
d NA
F A 0−F A +r A∗V = E c . 48
dt

Debido a que no hay entrada ni salida de A se obtiene la ecuación 49.

d NA
r A∗V = Ec . 49
dt
Balance de moles para B
dNB
F B 0−F B + r B∗V = Ecuación 50
dt

Se agrega controladamente B, pero no hay salida de esta especie (FB = 0), por lo
tanto, se obtiene la ecuación 51.
d NB
F B 0 +r B∗V = Ec . 51
dt
Balance de moles para C

14
 d NC
F C 0−FC +r C∗V = E c .52
dt
Donde: FCo = Velocidad de flujo molar de C a la entrada [mol/s]
FC = Velocidad de flujo molar de C a la salida [mol/s]
rC = Velocidad de formación de C [mol/L*s]
V = Volumen del reactor (variable) [L]
dNC/dt = Velocidad de acumulación de C [mol/s]

Debido a que no hay entrada ni salida de C se obtiene la ecuación 6.

d NC
r C ∗V = Ec .53
dt
Y empleando nuevamente las ecuaciones 10 a 12, 16 y 17, se obtienen las
ecuaciones para trabajar en POLYMATH. Recordando que la capacidad máxima
se alcanza cuando t = 8000 s y teniendo en cuenta los datos iniciales
(presentados en la Ilustración 6), se resuelve el sistema de ecuaciones y se
determina la conversión (X) alcanzada hasta ese momento y las moles de N A, NB y
NC que servirán como condiciones iniciales para la etapa 2.

Ilustración 6 Algoritmo empleado en POLYMATH Etapa 1

15
 Ilustración 7 Reporte POLYMATH Etapa 1

Ilustración 8 Conversión de Tolueno en función del Tiempo (s) Etapa 1

Ilustración 9 Moles de todas las especies (mol) en función del Tiempo (s) Etapa 1
16
Etapa 2

Se realiza el balance de moles para cada una de las especies, en esta

etapa el balance de moles para A y C es el mismo (Ecuación 49 y Ecuación
53 respectivamente), para el balance de moles de B se tiene en cuenta que
ya no se está alimentando, por lo cual el balance de moles de B queda
expresado en la ecuación 54.
d NB
r B∗V = Ec . 54
dt

Siguiendo un algoritmo análogo al de la etapa 1 en POLYMATH y teniendo en

cuenta las nuevas condiciones iniciales (presentadas en el algoritmo) para la
etapa 2, se logra determinar la conversión máxima y el tiempo en que se alcanza.

Ilustración 10 Algoritmo empleado en POLYMATH

Etapa 2

17
Ilustración 11 Reporte POLYMATH Etapa 2

La conversión máxima de tolueno lograda en el proceso batch es X = 0.9295.

Ilustración 12 Tabla de resultados POLYMATH Etapa 2 18

 La conversión máxima se alcanza en t = 7.586x105 s, ya que a partir de este
momento la conversión se vuelve constante.

Ilustración 13 Conversión del Tolueno en función del tiempo (s) Etapa 2

Ilustración 14 Moles de todas las especies (mol) en función del tiempo (s)

19
Para tener una mejor apreciación del comportamiento en todo el proceso de las
moles de cada una de las especies y la conversión en función del tiempo, se
exportan los resultados obtenidos en POLYMATH a un archivo en Excel para
unificar la etapa 1 y la etapa 2 y graficarlas.
1
0.9
0.8
0.7
Conversión Tolueno

0.6
0.5
0.4
0.3
0.2
0.1
0
0 100000 200000 300000 400000 500000 600000 700000 800000 900000
Tiempo (s)

Ilustración 15 Conversión del Tolueno en función del Tiempo (s) Proceso completo

8.00E+04

7.00E+04

6.00E+04

5.00E+04
Moles (mol)

4.00E+04
Na (mol)
3.00E+04 Nb (mol)
Nc (mol)
2.00E+04

1.00E+04

0.00E+00
0 200000 400000 600000 800000
Tiempo (s)

Ilustración 16 Moles de todas las especies (mol) en función del Tiempo (s) Proceso completo

d. Repetir el inciso c en el caso de destilación reactiva, si el DTE se evapora y no

acumula; considere αDTE igual a 0.0407 m3 / kmol. Analizar el efecto de la tasa
de evaporación en la conversión.

20
Consideraciones:
 Semi-lotes con respecto al componente A (tolueno)
 Se agrega controladamente el componente B (acetileno)
 Reacción en fase líquida
 Proceso isotérmico (20°C)
 Intensa agitación
 Reacción reversible
 Reacción de segundo orden respecto al tolueno y de primer orden respecto
al DTE
 El volumen es variable
 k1 =7x10-6 L/mol*s para la reacción directa
 Kc = 10 L/mol
 La alimentación se suspende cuando se alcanza la capacidad máxima
 Proceso de destilación reactiva desde el comienzo, el DTE (C) se evapora y
no acumula
 αC = 0.0407 L / mol
Repitiendo el procedimiento del inciso anterior (inciso c) teniendo en cuento a la
destilación reactiva desde el comienzo se modela el proceso en dos etapas; etapa
1: Se alimenta B hasta que la se alcanza la capacidad máxima del reactor (8000 L)
(proceso semi-batch), etapa 2: Luego de que se suspende la alimentación de B, se
deja reaccionar el sistema hasta que alcance el equilibrio (proceso batch).
Etapa 1
Se realiza el balance de moles para cada una de las especies, en esta
etapa el balance de moles para A y B se obtuvo previamente expresado en
las ecuaciones 49 y 51 respectivamente, para el balance de moles de C se
tienen en cuenta la destilación reactiva (C se evapora y no se acumula), por
lo tanto, el balance de moles de C queda expresado en la ecuación 56.
d NC
F C 0−FC +r C∗V = Ec . 55
dt
d NC
−F C + r C∗V = Ec . 56
dt
Ahora necesitamos definir el volumen, se realiza un balance de masa global
Balance de masa global
d ( ρV )
F Bo∗M wB −F C∗M wC + 0= Ec . 57
dt
Donde: Mwi = Peso molecular de la especie i [g/mol]

21
Considerando ρ = ρo (densidad [g/L])
F Bo∗M wB F C∗M wC dV
− = Ec .58
ρ ρ dt
Por definición:
M wi
=α i Ec . 59
ρ
Donde: αi [L/mol]
Así, la ecuación 58 queda
dV
F Bo∗α B −FC ∗α C = Ec . 60
dt
Teniendo en cuenta que
F Bo∗α B =v o Ec . 61

La ecuación 60 se puede reescribir así:

dV
v o−F C∗α C = Ec . 62
dt
A continuación, se debe definir el flujo de C, para lo cual se tiene en cuenta que
este se está evaporando y no se acumula
d NC
=0
dt
−F C + r C∗V =0 Ec . 63

r C ∗V =F C Ec . 64

Luego, se define el balance de moles para A (Ecuación 49) en términos de la

conversión y sabiendo que:
N A =N Ao∗( 1−X ) Ec .65

Reemplazando la ecuación 65 en la expresión del balance y operando se obtiene

la ecuación 66
d N Ao∗( 1−X )
r A∗V =
dt
−r A∗V dX
= Ec . 66
N Ao dt

Ahora, ya se tienen todas las ecuaciones para trabajar en POLYMATH y encontrar

el tiempo en el que se llega a la capacidad neta del tanque (8000 L) y así obtener,
22
los valores de X, NA, NB y NC en dicho momento, las cuales servirán como
condiciones iniciales para la etapa siguiente del proceso. Se debe tener en cuenta
que se emplean las ecuaciones 10 a 12 (con C C = 0) y las ecuaciones 45 a 47.
El algoritmo empleado en POLYMATH para la primera etapa del proceso, se
presenta en la Ilustración 17.

Ilustración 17 Algoritmo empleado en POLYMATH Etapa 1 Destilación

reactiva

23
Ilustración 18 Reporte POLYMATH Etapa 1 Destilación
reactiva

La capacidad máxima del reactor se alcanza a t = 10000 s, el cual es un tiempo

mayor al que se obtenía anteriormente para el aforo, esto es debido que además
de que se está ingresando B también sale C, por lo es evidente que dicha
capacidad efectivamente se alcanzaría en un tiempo mayor. También se puede
observar que la conversión del tolueno cuando se alcanza la capacidad máxima
del reactor es X = 0.6552, lo cual es un valor mayor a la conversión alcanzada en
el proceso semi-batch sin destilación reactiva obtenido anteriormente X = 0.5912
(inciso c).

Etapa 2
Realizando el balance de moles para todas las especies se obtienen nuevamente
las ecuaciones 49 y 54 para A y C (no hay acumulación de C), mientras que para
el balance de moles de B se tiene en cuenta que ya no se está alimentando, por lo
cual se obtiene la ecuación 54.
Se debe tener en cuenta que se emplean las ecuaciones 10 a 12 (con C C = 0), las
ecuaciones 45 a 47 y las ecuaciones 62, 64 y 66.
El algoritmo empleado en POLYMATH para la etapa 2 del proceso se presenta en
la Ilustración 19.

24
Ilustración 19 Algoritmo empleado en POLYMATH Etapa 2
Destilación reactiva

25
Ilustración 20 Reporte POLYMATH Etapa 2 Destilación
reactiva

De la ilustración 20 se puede evidencia que para un tiempo alto (t = 2’000.000 s)

se obtiene una conversión elevada (X = 0.9964), haciendo una inspección
variando los valores finales para t, se evidencia que la conversión seguirá
aumentando poco a poco pero no llega a la unidad debido a que la reacción es
reversible. Se toma como tiempo máximo de conversión t = 800000 s para realizar
una comparación de la variación de la conversión de tolueno en función del tiempo
con el proceso anterior (cuando no había destilación reactiva), el resultado se
muestra en la ilustración 21.
Los resultados obtenidos en los algoritmos de POLYMATH se exportan a un
archivo de Excel para unificar las etapas 1 y las etapas 2 de ambos procesos y
graficarlas.

26
 1
0.9
Conversión Tolueno 0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 100000 200000 300000 400000 500000 600000 700000 800000
Tiempo (s)

Conversión sin destilación reactiva Conversión con destilación reactiva

Ilustración 21 Conversión de tolueno en función del Tiempo (s) para los procesos sin y con destilación reactiva

En la ilustración 21 se puede evidenciar que la conversión para un mismo intervalo

de tiempo es mayor cuando en el proceso existe destilación reactiva, esto se debe
a que, al sacar constantemente C, por principio de Le Chatelier, el equilibrio se
desplaza a la derecha de la reacción (Reacción 4) aumentando la conversión.

2. El estireno se obtiene por deshidrogenación catalítica de etilbenceno:

C 6 H 5 ∙CH 2∙ CH 3 ↔C 6 H 5 ∙ CH =CH 2+ H 2 Reacción 7

La ecuación de velocidad para esta reacción se da en la Ecuación 67.

1
r ' =k P Et − ( )
P P Ec .67
K P St H

Donde PEt, PSt y PH son presiones parciales de Etilbenceno, Estireno e Hidrógeno,

respectivamente.
El reactor constará de una serie de tubos de ¼ de diámetro nominal cada uno, con
cédula 40 rellenos de catalizador con una densidad aparente de 1500 kg / m 3. El
etilbenceno se diluye con vapor, las velocidades de alimentación por unidad de
área de sección transversal son etilbenceno 1.6x10-3 kmol / (m 2.s), vapor 2.9x10-2

27
kmol /(m2.s). El reactor funciona a una presión promedio de 1.2 bar y la
temperatura se mantiene a 833 K en todo momento.

A 833 K: k = 6.6x10-11 (kmol.m2) / (N.s.kg de catalizador)

KP = 1.0x104 N / m2.

a. Si se requiere una conversión fraccional de etilbenceno de 0,55, calcule la

longitud y el número de tubos necesarios para producir 15 toneladas de estireno
por día. ¿Cuánta masa de catalizador es necesaria para lograr el requisito de
producción? Ignore la caída de presión.
Consideraciones:
 Sistema sin caída de presión
 Sistema isobárico, Isotérmico
 Sistema en fase gaseosa

Partiendo de la ecuación de diseño:

dX −r ' A
= Ec .68
dt F A0

1
(
r A ' =−k P Et − )
P P Ec . 69
K P St H

En el enunciado dan flux molares, entonces se debe multiplicar por el área de

sección transversal para obtener el flujo másico; como se muestra a continuación:

F A 0=FluxA∗Ac Ec . 70

F V 0=FluxV apor∗Ac Ec .71

F T 0=F V 0 + F A 0 Ec .72

Con cédula 40 y D=1/4 in se determina el área de sección transversal

Ac=6.717*10-5 m2

28
Por definición:
P j=C j RT Ec . 73

Donde P j es la presión parcial del componente j (Pa)

C j es la concentración de la especie j (mol/L)

R es la constante de los gases (Pa*L/mol*K)

Evaluando para cada componente:

P A =C A RT Ec . 74

PB =C B RT Ec .75

PC =C C RT Ec .76

C AO (θ j + γ j X)
Por definición: C j= Ec .77
(1+ EX)

E=Y A 0 δ Ec . 7 8

Donde:
Y A 0es la fracción de vapor de A

δ es la suma de los coeficientes estequiométricos de las especies

F T 0 flujo molar total inicial

F A0
Y A 0= Ec .79
FT 0

C A 0=Y A 0∗C T 0 Ec . 80

P0
C T 0= Ec . 81
RT 0
Donde:
C T 0 es la concentración total inicial.

29
Evaluando para cada componente:
C AO (1−X )
C A= Ec . 82
(1+ EX )
C AO X
C B= Ec . 83
( 1+ EX )
C AO X
C C= Ec .84
(1+ EX)

F A=F A 0 ( 1−X ) Ec . 85

F B=F A 0 X Ec . 86

F C =F A 0 X Ec .87

Para determinar la longitud de los tubos del reactor se usa la siguiente ecuación:
w
z= Ec . 88
ρb∗Ac

Donde:
w es el peso del catalizador (Kg)
ρb densidad del catalizador (Kg/m3)
Ac área de sección transversal (m2)
z longitud de los tubos (m)

Para determinar el número de tubos (n) se usa la siguiente ecuación:

15000
n= Ec . 89
FB0

Con
86400∗F A 0 X∗MM B
F B 0= Ec . 9 0
1000
Donde:
MMB peso molecular de B, 104.15g/mol

Se resuelve el sistema de ecuaciones desde la ecuación 68 a la 90 excepto

la 73 y 77 en Scilab y, para una conversión de X=0.55 se obtuvo un peso

30
de catalizador de w=0.2491 Kg; se obtiene una longitud de tubos z=2.4726
metros y un numero de tubos n=28200.7321.
El algoritmo utilizado se muestra en la ilustración 22.

31
32
Ilustración 22 Algoritmo empleado en Scilab

b. Grafique la conversión y los flujos molares de todas las especies en función

de la longitud del reactor.
Con las ecuaciones planteadas para los flujos de cada componente 85, 86 y 87, se
resuelve en Scilab mediante el algoritmo empleado en el literal a. En la ilustración
23 se muestra la variación de los flujos molares vs longitud del reactor, y la
conversión vs la longitud del reactor.

Ilustración 23. Conversión vs longitud del reactor

33
Ilustración 24. Flujos molares vs longitud del reactor.

3. Los ciclohexanoles alquilados como 2-metilciclohexanona son productos

intermedios importantes en la industria de las fragancias y los perfumes. Un
trabajo informado [Ind. Ing. Chem. Re., 28, 693 (1989)] se centró en la
hidrogenación catalizada en fase gaseosa de o-cresol a 2-metilciclohexanona,
Reacción 8, que luego puede hidrogenarse a 2-metilciclohexanol. Se encontró que
la reacción representada con la Reacción 3 en un catalizador de níquel-sílice era
de orden cero en o-cresol y de primer orden en hidrógeno con una velocidad de
reacción específica a 170 ° C de 1,74 mol de o-cresol / (kgcat.min.atm). La mezcla
de reacción entra en el reactor de lecho compacto a una presión total de 5 atm. La
alimentación molar consta de un 67% de H2 y un 33% de o-cresol a un caudal
molar total de 40 mol / min.

Reacción . 8

a. Sin tener en cuenta la caída de presión, grafique la conversión de o-cresol

y las concentraciones de cada especie en función del peso del catalizador.
¿Cuál es la relación entre el peso del catalizador necesario para lograr el
90% de conversión y el peso necesario para lograr el 10% de conversión en
el PBR?

34
 Consideraciones:
 Sistema en fase gaseosa
 Sistema sin caída de presión
 Proceso isotérmico
 Reacción de orden cero en O-cresol y de primer orden en Hidrógeno.

La reaccion a analizar se puede representar de la siguiente manera:

A+2 B→ C Reacción 9

Donde A representa al o-cresol, B al Hidrógeno y la C a 2-metilciclohexanona.

De la reacción 9, se define el reactivo límite como la especie A (o-cresol), porque

la especie B (hidrógeno) se encuentra en exceso para que pueda ocurrir la
reacción.

Balance de Moles para A

F A|w − F A|w+∆ w +r ' A ∆ w=0 Ec .9 1

Donde:
FA es el flujo molar referente al reactivo limite (mol/min)
ΔW se refiere a un diferencial del peso del catalizador
r’A es la velocidad de reaccion del compuesto A
Sacando el límite cuando ∆ w tiende a cero:

F A|w − F A|w+ ∆ w + r ' A ∆ w

lim =0 Ec . 92
∆ w→ 0 ∆w

−d F A '
+r A =0 Ec . 93
dw

Organizando la ecuación 3, se obtiene:

d FA '
=r A Ec . 9 4
dw

Por definición:
35
F A=F A 0∗( 1−X ) Ec . 9 5

Donde FA0 es el flujo molar inicial del compuesto A y X es la conversión.

Derivando F A, se obtiene:

d F A =−F A 0 dX Ec . 9 6

Reemplazando la ecuación 6 en la ecuación 4:

dX
−F A 0 =r ' A Ec . 97
dw

Organizando la ecuación 7, se obtiene la expresión de conversión en términos del

peso del catalizador:

dX −r A
= Ec . 9 8
dw F A 0
Velocidad global de reacción:

r '=K P B
Velocidad de desaparición respecto al reactivo limite (o-cresol).

r A =γ .r ' Ec . 9 9
Donde:
γ =−1 (Coeficiente estequiométrico de A)

−r A =K PB Ec .10 0
Por definición, para gas ideal:

P j=C j∗R∗T Ec .1 0 1
Para la especie B:
PB =C B∗R∗T Ec .1 0 2

Por definición, para el cálculo de concentraciones en un reactor PBR:

36
 C A 0 (θ j+ γ j X )
C j= Ec .1 03
( 1+ εX )

Para A:
C (1−X )
C A= A 0 Ec .1 0 4
( 1+ εX )

Para B:

C A 0 (θ B−2 X )
C B= Ec .1 0 5
( 1+ εX )
Para C:

C A 0(θC + X )
C C= Ec .1 0 6
( 1+ εX )

Por definición:
y A 0 0,33
θA= = =1 Ec .1 0 7
y A 0 0,33
y 0,67
θ B= B 0 = =2,03 Ec .1 0 8
y A 0 0,33

yC 0 0
θC = = =0 Ec .1 0 9
y A 0 0,33

Por definición:

C A0
y A0= Ec .110
CT 0
Despejando C A 0:
C A 0=C T 0∗y A 0 Ec .11 1

Con:

P0
C T 0= ( )
RT0
Ec .112

Reemplazando C T 0 en la ecuación 111, se obtiene:

P0
C A 0= ( )
RT0
∗y A 0 Ec .113

37
Reemplazando la ecuación 105 en la ecuación 102:

C A 0 ( θ B−2 X )∗R∗T
PB = Ec .11 4
( 1+εX )

La ecuación 114 se reemplaza en la ecuación 100:

C A 0 ( θ B−2 X )∗R∗T
−r A =K ' Ec .11 5
( 1+εX )

Con:
ε = y A 0∗δ Ec . 116

δ =(−1−2+ 1 )=−2

y A 0 =0.33

Luego el valor de la ecuación 116 es:

ε =−2∗0.33=−0.66 Ec .117

Finalmente, la ecuación 98 queda expresada como:

dX K∗C A 0 ( θ B−2 X )∗R∗T

= Ec .11 8
dW F A 0 ( 1+ εX )

Las ecuaciones 104, 105, 113, 115, y 118 fueron ingresadas en MATLAB para
obtener las gráficas de conversión de o-cresol y las concentraciones de cada
especie en función del peso del catalizador. El algoritmo se muestra en la
ilustracion 25.:

38
39
Ilustración 24.Algoritmo empleado en MATLAB

Ilustración 25.Conversión en función del peso del catalizador

En la ilustración 26 se puede observar que a medida que aumenta el peso del

catalizador, aumenta la conversión; lo cual concuerda con lo esperado.

40
Ilustración 26.Concentraciones en función del peso del catalizador
A medida que aumenta el peso del catalizador la concentración de los reactivos
(especies A y B) disminuyen, mientras que para el producto aumenta debido a que
éste se va formando (especie C), lo cual es acorde con la ilustración 26 porque a
medida que aumenta la conversión de A va disminuyendo su concentración inicial.
La relación entre el peso del catalizador necesario para lograr el 90% de
conversión y el peso necesario para lograr el 10% de conversión en el PBR es
13.1, como se muestra en el resultado obtenido con MATLAB:

Ilustración 27. Relación del peso del catalizador

b. Teniendo en cuenta la caída de presión en el lecho empacado usando un valor

de α = 0.027 kgcat-1, rehaga la parte a. Grafique la presión frente al peso del
catalizador.

41
Consideraciones:
 Sistema isotérmico
 Sistema en fase gas
 Sistema con caída de presón

La expresión para concentraciones con caida de presión e isotérmico es:

C A 0( θ j + γ j X )
∗P
( 1+εX )
C j= Ec .119
P0

Para A:
C A 0 ( 1− X )
∗P
( 1+ εX )
C A= Ec . 120
P0

Para B:

C A 0 (θB −2 X )
∗P
( 1+ εX )
C B= Ec . 121
P0

Para C:

C A 0 (θC + X)
∗P
( 1+εX )
C C= Ec .123
P0

Con valores de θ B,θC , ε ,δ, y A 0 iguales a la parte a

La ecuación 114 queda:

C A 0 ( θB −2 X )∗R∗T
∗P
( 1+ εX )
PB = Ec . 124
P0

La ecuación 100 se reemplaza en la ecuación 124:

42
 K∗C A 0 ( θ B −2 X )∗R∗T
∗P
( 1+ εX )
−r A = Ec . 125
P0

Finalmente, la ecuación 98 queda expresada como:

K∗C A 0 ( θB −2 X )∗R∗T
∗P
dX F A 0 ( 1+εX )
= Ec . 126
dW P0

A partir de la ecuación de Ergun, se obtiene la ecuación para determinar la caída

de presión en función del peso del catalizador:

dP −∝ P0
dW
=
2 P
P0 ()
( 1+εX ) Ec .127

Las ecuaciones 120, 121, 123, 125, 126 y 127 fueron ingresadas en MATLAB para
obtener las gráficas de conversión de o-cresol, las concentraciones de cada
especie, y la presión en función del peso del catalizador. El algoritmo se muestra
en la ilustración 29.

43
Ilustración 28.Algoritmo empleado en MATLAB

44
Ilustración 309.Conversión en función del peso del catalizador

Ilustración 31. Concentraciones en función del peso del catalizador

45
Ilustración 32. Presión en función del peso del catalizador

La relación entre el peso del catalizador necesario para lograr el 90% de

conversión y el peso necesario para lograr el 10% de conversión en el PBR es
13.231,como se muestra en el resultado calculado con MATLAB ilustración 33.

Ilustración 33. Relación del peso del catalizador con caída de presión

46
 c. Otro ingeniero sugiere que cambiar el tamaño del catalizador sería mejor para
la reacción y recomienda que el diámetro de partícula debe ser la tercera parte, el
ingeniero sugiere que con una cantidad de 4 kg de catalizador la reacción debe
lograr una mejor conversión de la condición inicial. Desde el punto de vista de un
ingeniero, ¿qué puedes decir al respecto?
El ingeniero propone esto :
Dp1
D p 2= Ec . 128
3
La ecuación de Ergum es:

dP −G 1−∅ 150 ( 1−∅ ) µ P T

=
( )[
dZ ρ0 g c D p ∅ 3 Dp ]( )( )
+1,75 G
P0 T0
( 1+ εX ) Ec . 129

Donde:

D p : diámetro de párticula (m)

G : Flujo másico (Kg/m2*s)

∅ : Porosidad de la partícula
gc : Factor de conversión (Kg*m/s2)

µ: Viscosidad del gas (Kg/m*s)

ρ0: Densidad del gas (Kg/m3)

De la ecuación 129 vemos que predomina el primer término, es decir:

3∗150 ( 1−∅ ) µ
≫1.75 G
Dp1

Osea que domina el flujo laminar (1.75 = 0), por tanto, de la ecuación 129, se tiene
que:

β 0=
−G
( )[ 150( 1−∅
1−∅
ρ0 g c D p ∅ 3 D
)µ
p
] Ec .130
De la ecuación 130 observamos que el parámetro de caída de presión varía
inversamente con el diámetro de la partícula al cuadrado:

47
 1
β0
D p2

Por lo tanto:

1
α
D p2

D p1
Para el caso que propone el ingeniero, D p 2=
3

∝2∗D2p 2=∝1∗D2p 1 Ec . 131

Despejando para ∝2:
∝ ∗D2
∝2= 1 2 p 1 Ec .132
Dp2
Con:
2 D 2p 1
D p 2= 2
3

D 2p 1
9= 2
D p2

∝2=9∗∝1 Ec . 133
Donde:
∝ :Factor de caída de presión (Kgcata-1)

El sistema de ecuaciones planteadas, se programa en en MATLAB y los

resultados obtenidos se muestran en la ilustración 34.

48
Ilustración 3410. Algoritmo para ver la variación de conversión y velocidad de
reacción

49
Ilustración 35. Velocidad de reacción con Dp/3

Ilustración 36. Variación de la conversión con Dp/3 y 4 kg de catalizador

50
Al disminuir el diametro de la particula, aumenta el parámetro de caida de presión
α, y por lo tanto la velocidad de reacción y la conversión disminuyen como se
puede ver en las ilustraciones 35 y 36.
Como se puede ver en la ilustración 30, cuando el peso del catalizador es de 4 kg
la conversión era aproximadamente 1, y en la ilustración 36 es aproximadamente
0.92, es decir, es mas baja que la anterior.
El parámetro que hizo que disminuyera la conversión fue el diámetro de la
partícula, es decir, reducir el diametro un tercio de su valor inicial implica que se
debería aumentar la cantidad de catalizador para lograr una conversión
aproximadamente de 1.

d. Suponga que el reactor opera en condiciones no isotérmicas y la temperatura

está relacionada linealmente con la conversión de tal manera que para una
conversión fraccional igual a 1, la temperatura es 270 ° C y la conversión es
cero a 170 ° C. Rehaga las partes b y c. Además, suponga para los cálculos
que la constante de velocidad de reacción no varía significativamente en este
rango de temperatura.
Se utilizan las ecuaciones obtenidas en los literales b y c, y adicionalmente se
tiene una ecuación para la temperatura:
Para T=170ºC X=0
T=270ºC X=21
Con estos datos se obtiene la ecuación de una linea recta:
T −T 0=m∗( X−X 0 ) Ec . 134
T 2−T 1
m= Ec . 135
X 2− X 1

270−170
m= =100 Ec .136
1−0
T −170=100∗( X −0 ) Ec . 137

T =100 X +170 Ec . 138

51
4. W. E. Walles y A. E. Platt [Ind. Ing. Chem., 59 (6), 41 (1967)] informó que la
descomposición térmica de 5-metil-2-oxazolidinona es autocatalítica, Reacción 4.

Reacción. 10

La expresión de velocidad para esta reacción en fase líquida es:

r =k 1 C 2A + k 2 C A C B Ec .1 39

A la temperatura de interés:
k1 = 1.02x10-2 L/(mol*s)
k2 = 7.5x10-2 L/(mol*s)
La reacción se lleva a cabo en un sistema de tres reactores CSTR en serie, la
alimentación consta de una especie A disuelta en un inerte (I) con una
concentración de 0,6 M, el caudal volumétrico de la alimentación es de 180 L/min.
Los volúmenes de los reactores son 70, 100 y 60 L, respectivamente.

52
Determine:
a. ¿Cuál es la conversión total de estado estable de A?

Consideraciones:
 Sistema en estado estacionario.
 Sistema en fase liquida, (ρ=constante y v 0=v )
 Sistema con entrada y salida continua.
 Sistema con buena agitación y propiedades homogéneas dentro de
los reactores.
 El reactivo limite es A.

Teniendo en cuenta las consideraciones anteriores, el balance de moles

para el CSTR es:

F j 0−F j
V CSTR = Ec . 140
−r Aj

La ecuación 2 se aplica a cada uno de los reactores.

Reactor 1:
F A 0−F A 1
V CSTR ,1= Ec . 141
−r A 1
Por definición:

F A 1=F A 0∗( 1−X 1 ) Ec . 142

Reemplazando la ecuación 142 en la 141, se tiene:

F A0 X1
V CSTR ,1= Ec . 143
−r A 1

A partir de la ecuación 139 se plantea la velocidad de desaparición del

reactivo límite A, teniendo en cuenta que:

r A =γ .r ' Ec . 144

53
 r A 1=−2 ( k 1 C2A 1+ k 2 C A 1 C B 1 ) Ec .145

Reemplazando la ecuación 145 en la 143, se obtiene:

F A0 X1
V CSRT ,1= 2
Ec . 146
2(k 1 C A 1+ k 2 C A 1 C B 1)

Por definición: C j=C AO ( θ j −γ j X ) Ec . 1 47

Aplicando las ecuaciones 147 se determinan las concentraciones C A y CB

en términos de la conversión X.

Se tiene que:
FB0
θ B= =0 Ec .148
FA 0

C A 1=C AO ( 1−X 1 ) Ec . 149

C AO X 1
C B 1= Ec .1 50
2

Por definición:
F A 0=C AO v 0 Ec .1 51

Reactor 2:
Partiendo de la ecuación de diseño para el CSTR:

F A 1−F A 2
V CSTR ,2= Ec .1 52
−r A 2

Para calcular −r A 2 y las Cj se realiza el mismo procedimiento que se hizo

para el reactor 1.

r A 2=−2(k 1 C 2A 2+ k 2 C A 2 C B 2) Ec . 153

C A 2=C AO ( 1−X 2 ) Ec . 15 4

C AO X 2
C B 2= Ec .1 55
2

54
 Se tiene que:
F A 1=F A 0 ( 1−X 1 ) Ec . 156

F A 2=F A 0 ( 1−X 2) Ec . 157

Reemplazando las ecuaciones 156 y 157 en la ecuación 152, la ecuación

para el volumen del reactor 2 se puede expresar como sigue:

F A 0 (X 2 −X 1 )
V CSTR ,2= Ec .158
−r A 2

Reactor 3.
Se realiza el mismo procedimiento del reactor 2, y se obtienen las
siguientes ecuaciones:

r A 3=−2 ( k 1 C 2A 3+ k 2 C A 3 C B 3 ) Ec . 159

C A 3=C AO ( 1−X 3 ) Ec . 160

C AO X 3
C B 3= Ec .161
2

F A 0 (X 3 −X 2 )
V CSTR ,3= Ec . 162
−r A 3

Finalmente, se estima la conversión de A mediante un solver en la

herramienta computacional Excel. Para ello, se usan las ecuaciones 141,
145, 149, 150, 152, 153, 154, 155, 158 A 161. Los resultados obtenidos se
muestran en la tabla 1.

Tabla 1.Concentración de las especies y conversiones en cada reactor

Concentración (mol/L)   Conversión X

CA1 0.3846      
CA2 0.1765   X1 0.3590
CA3 0.0985   X2 0.7058
CB1 0.1077   X3 0.8359

55
 CB2 0.2117    
CB3 0.2508      

La conversión de A en estado estacionario corresponde a la conversión

final del sistema, es decir, a la conversión en el último reactor; la cual es
0.8359. Se puede evidenciar que la conversión de A aumenta conforme
aumenta la cantidad de reactores CSTR dispuestos en serie.

b. Determine el tiempo necesario para alcanzar el estado estable (es decir,

cuando CA que sale del tercer reactor es el 99% del valor del estado
estable).

En estado estacionario CA_SS=0.0995 mol/L

Partiendo del balance de moles:

dN j
F j 0−F j 1 +r j V = Ec . 163
dt

Por definición:
F j=C j v 0 Ec . 164

Cj
N j= Ec .165
V

Reemplazando las ecuaciones 164 y 165 en la 163 y teniendo en cuenta

que el volumen es constante:

V∗dC j
C 0 j v 0−C j v 0 +r j V = Ec .166
dt

Al dividir por el volumen V y agrupando términos, se obtiene:

(C ¿ ¿ 0 j−C j ) dC j
v0 + r Aj= Ec . 167 ¿
V dt

Por definición:
Vj
τ j= Ec .168
v0

Reemplazando la ecuación 168 en la 167, resulta:

56
 (C ¿ ¿ 0 j−C Aj ) dC
+r j = j Ec . 169¿
τj dt

Aplicando la ecuación 168 y 169 a cada reactor, se obtiene:

V1
τ1= Ec .170
v0
V
τ 2= 2 Ec .171
v0
V
τ 3 = 3 Ec .172
v0

Para el compuesto A:

(C ¿ ¿ A 0−C A 1 ) dC
+r A 1= A 1 Ec .173 ¿
τ1 dt
(C ¿ ¿ A 1−C A 2) dC
+ r A 2= A 2 Ec . 174 ¿
τ2 dt
(C ¿ ¿ A 2−C A 3) dC
+r A 3 = A 3 Ec .175 ¿
τ3 dt

Para el compuesto B:

(C ¿ ¿ B 0−C B 1) dC
+r B 1 = B 1 Ec .176 ¿
τ1 dt

(C ¿ ¿ B 1−C B 2 ) dC
+ r B 2= B 2 Ec . 177 ¿
τ2 dt

(C ¿ ¿ B 2−C B 3 ) dC
+ r B 3= B 3 Ec .178 ¿
τ3 dt

Para el compuesto C.
(C ¿ ¿ C 0−CC 1) dC
+r C 1= C 1 Ec .198 ¿
τ1 dt

57
(C ¿ ¿ C 1−C C 2) dC
+ r C 2= C 2 Ec . 180 ¿
τ2 dt

(C ¿ ¿ C 2−C C 3) dC
+r C 3= C3 Ec .181 ¿
τ3 dt

Para obtener las velocidades de aparición de B y C se aplica la expresión para

velocidades relativas:
rA rB rC
= = Ec . 182
−1 1/2 1/2

−r A 1
r B 1= Ec . 183
2

−r A 2
r B 2= Ec . 184
2

−r A 3
r B 3= Ec . 185
2

−r A 1
r C 1= Ec .186
2

−r A 2
r C 2= Ec .187
2

−r A 3
r C 3= Ec . 188
2

Haciendo uso de las ecuaciones 170 a la 188 se resuelve en el software

Polymath. El código se muestra en las ilustraciones iiiii.

58
Ilustración. Algoritmo utilizado en Polymath

59
Ilustración 11. Algoritmo utilizado en Polymath

Del código, para una concentración de C A= 0.0995 mol/L se obtiene un tiempo

ts=4.1572 minutos, correspondiente al tiempo en que se alcanza el estado
estacionario.

c. Grafique la concentración de todas las especies que salen de cada tanque

en función del tiempo.

60
Ilustración 12. Concentraciones vs tiempo

d. (6%) Varíe los parámetros del sistema (υ0, V, k) y determine sus efectos en
las condiciones de estado estacionario. Escriba un párrafo que describa las
tendencias que encontró, incluida una discusión del parámetro que tiene la
mayoría de los efectos en los resultados.

Para dar solución a este literal, se hacen variaciones de los parámetros de

υ0, V, y k al algoritmo usado en el literal 4.b. para el estado estacionario, por
lo cual no se considera necesario mostrar el código empleado; los
resultados obtenidos se muestran en las ilustraciones de las
concentraciones del compuesto A vs el tiempo.

 Variación del flujo volumétrico υ0.

61
lustración 13. Concentración de A en cada reactor vs tiempo(min). Para: ʋ0 =180 (L/min)

Ilustración 14. Concentración de A en cada reactor vs tiempo(min). Para: ʋ0 =360 (L/min)

62
Ilustración 15. Concentración de A en cada reactor vs tiempo(min). Para: ʋ0 =90 (L/min)

Al duplicar el flujo volumétrico se observa una mayor cantidad de concentración de

A en cada reactor, esto se debe a que se dispone de más reactivo; pero al
aumentar el flujo volumétrico bajo las mismas condiciones de operación infiere una
menor velocidad de reacción y por ende una menor conversión del producto
deseado. Mientras que, al reducir el flujo volumétrico la velocidad de reacción
aumenta y el proceso se lleva a cabo más rápido, aunque con menos cantidad de
concentraciones, esto se puede evidenciar en las separaciones de las curvas
entre el reactivo y los productos porque están mas distantes.

 Variación del volumen V.

63
Ilustración 16. Concentración de A en cada reactor vs tiempo (min). Para: Vj=2Vj (L)

Ilustración 17. Concentración de A en cada reactor vs tiempo(min). Para: V=V/2 (L)

64
  Variación de K

Ilustración 18. Concentración de A en cada reactor vs tiempo(min). Para: kj =10kj

Ilustración 19. Concentración de A en cada reactor vs tiempo(min). Para: kj =kj/10

65
De la variación de los tres parámetros que intervienen en el proceso el que
más influye es la constante k cuando se aumenta o se reduce en un factor
de 10, esto se puede evidenciar en las ilustraciones 26 y 27 al compararlas
con la ilustración 21 que es la del proceso original sin ninguna modificación.
Cuando el valor de la constante se reduce 10 veces de su valor inicial, se
observa que las concentraciones de los tres reactores tienden a juntarse lo
que infiere una baja conversión y, por consiguiente, un proceso menos
eficiente, mientras que al aumentar 10 veces el valor de k se observa que
las curvas de concentración presentan bastante separaciones razón por la
cual se puede concluir que se esta dando una gran conversión de A.

66