1) Esquema del reactor:

𝑲
2) Reacción involucrada: 𝑪𝑵𝑩𝒓 + 𝑪𝑯𝟑 → 𝑪𝑯𝟑 𝑩𝒓 + 𝑵𝑪𝑵𝑯𝟐
3) Fase líquida Volumen Constante.
Donde 𝐶𝑁𝐵𝑟 : 𝐴 ; 𝐶𝐻3 : 𝐵 ; 𝐶𝐻3 𝐵𝑟 : 𝐶 (Producto Deseado); y 𝑁𝐶𝑁𝐻2 : 𝐷
4) Balance molar:
 Para A: 𝑬𝑨 − 𝑺𝑨 + 𝑮𝑨 − 𝑪𝑨 = 𝑨𝑨 −𝑪𝑨 = 𝑨𝑨 (1)
𝒅𝑵𝑨
Donde: 𝑪𝑨 = (−𝒓𝑨 )𝑽 y 𝑨𝑨 = sustituimos en (1):
𝒅𝒕
𝒅𝑵𝑨
−(−𝒓𝑨)𝑽 = Sabemos que: 𝑁𝐴 = 𝐶𝐴 𝑉 ; sustituimos:
𝒅𝒕

𝒅(𝑪𝑨 𝑽)
−(−𝒓𝑨 )𝑽 =
𝒅𝒕
𝒅(𝑪𝑨 ) 𝒅(𝑽)
−(−𝒓𝑨 )𝑽 = 𝑽 + 𝑪𝑨 (𝑰)
𝒅𝒕 𝒅𝒕
𝒅(𝑽)
Debemos hallar la expresión , haciendo un balance de masa global
𝒅𝒕
𝒅𝒎
𝑬−𝑺+𝑮−𝑪=𝑨 𝒎𝒆 = 𝒅𝒕

Sabemos que: 𝒎 = 𝝆𝑽 sustituimos:

𝒅(𝝆𝑽)
𝝆𝟎 𝒗𝟎 =
𝒅𝒕
𝒅(𝑽) 𝒅(𝝆)
𝝆𝟎 𝒗𝟎 = 𝝆 +𝑽 Como estamos en fase líquida la densidad es constante, por lo
𝒅𝒕 𝒅𝒕
tanto:
𝒅(𝑽)
=𝟎
𝒅𝒕
Por otro lado:
𝒅(𝑽) 𝒅(𝑽)
𝝆𝟎 𝒗𝟎 = 𝝆 ; 𝝆𝟎 = 𝝆 Entonces: 𝒗𝟎 = (II)
𝒅𝒕 𝒅𝒕
𝒕 𝑽
𝒗𝟎 𝒅𝒕 = 𝒅𝑽 ⇒ 𝒗𝟎 ∫ 𝒅𝒕 = ∫ 𝒅𝑽 ⇒ 𝒗𝟎 𝒕 = 𝑽 − 𝑽𝟎 ⇒ 𝑽 = 𝒗𝟎 𝒕 + 𝑽𝟎 (∗)
𝟎 𝑽𝟎

Sustituimos (II) en (I):

𝒅(𝑪𝑨 )
−(−𝒓𝑨 )𝑽 = 𝑽 + 𝑪𝑨 𝒗𝟎
𝒅𝒕
𝒅(𝑪𝑨 )
−(−𝒓𝑨 )𝑽 − 𝑪𝑨 𝒗𝟎 = 𝑽
𝒅𝒕
𝒗𝟎 𝒅(𝑪𝑨 )
−(−𝒓𝑨 ) − 𝑪𝑨 = (𝑰𝑰𝑰)
𝑽 𝒅𝒕

 Para B: 𝑬𝑩 − 𝑺𝑩 + 𝑮𝑩 − 𝑪𝑩 = 𝑨𝑩 ⇒ 𝑬𝑩 − 𝑪𝑩 = 𝑨𝑩 (2)
𝒅𝑵𝑩
Donde: 𝑬𝑩 = 𝑭𝑩 = 𝑪𝑩 𝒗𝟎 ; 𝑪𝑩 = (−𝒓𝑩 )𝑽 y 𝑨𝑩 = sustituimos en (2):
𝒅𝒕

𝒅𝑵𝑩 𝒅(𝑪𝑩 𝑽)
𝑪𝑩𝟎 𝒗𝟎 − (−𝒓𝑩 )𝑽 = ⇒ 𝑪𝑩𝟎 𝒗𝟎 − (−𝒓𝑩 ) =
𝒅𝒕 𝒅𝒕
𝒅(𝑪𝑩 ) 𝒅(𝑽)
⇒ 𝑪𝑩𝟎 𝒗𝟎 − (−𝒓𝑩 )𝑽 = 𝑽 + 𝑪𝑩
𝒅𝒕 𝒅𝒕
𝒅(𝑪𝑩 ) 𝒅(𝑪𝑩 )
⇒ 𝑪𝑩𝟎 𝒗𝟎 − (−𝒓𝑩 )𝑽 = 𝑽 + 𝑪𝑩 𝒗𝟎 ⇒ 𝑪𝑩𝟎 𝒗𝟎 − (−𝒓𝑩 )𝑽 − 𝑪𝑩 𝒗𝟎 = 𝑽
𝒅𝒕 𝒅𝒕

𝒗𝟎 𝒗𝟎 𝒅(𝑪𝑩 ) 𝒅(𝑪𝑩 ) 𝒗𝟎
⇒ 𝑪𝑩𝟎 − (−𝒓𝑩 ) − 𝑪𝑩 = ⇒ = (𝑪𝑩𝟎 − 𝑪𝑩 ) − (−𝒓𝑩 ) (𝑰𝑽)
𝑽 𝑽 𝒅𝒕 𝒅𝒕 𝑽

 Para C: 𝑬𝑪 − 𝑺𝑪 + 𝑮𝑪 − 𝑪𝑪 = 𝑨𝑪 ⇒ 𝑮𝑪 = 𝑨𝑪 (3)
𝒅𝑵𝑪
Donde 𝑮𝑪 = (𝒓𝑪 )𝑽 y 𝑨𝑪 = Sustituimos en (3):
𝒅𝒕

𝒅𝑵𝑪 𝒅(𝑪𝑪 𝑽) 𝒅(𝑪𝑪 ) 𝒅(𝑽)

(𝒓𝑪 )𝑽= ⇒ (𝒓𝑪 )𝑽 = ⇒ (𝒓𝑪 )𝑽 = 𝑽 + 𝑪𝑪
𝒅𝒕 𝒅𝒕 𝒅𝒕 𝒅𝒕

𝒅(𝑪𝑪 ) 𝒅(𝑪𝑪 )
⇒ (𝒓𝑪 )𝑽 = 𝑽 + 𝑪𝑪 𝒗𝟎 ⇒ (𝒓𝑪 )𝑽−𝑪𝑪 𝒗𝟎 = 𝑽
𝒅𝒕 𝒅𝒕

𝒅(𝑪𝑪 ) 𝒗𝟎
⇒ = (𝒓𝑪 )−𝑪𝑪 (𝑽)
𝒅𝒕 𝑽
  Para D: 𝑬𝑫 − 𝑺𝑫 + 𝑮𝑫 − 𝑪𝑫 = 𝑨𝑫 ⇒ 𝑮𝑫 = 𝑨𝑫 (4)
𝒅𝑵𝑫
Donde 𝑮𝑫 = (𝒓𝑫 )𝑽 y 𝑨𝑫 = Sustituimos en (4):
𝒅𝒕

𝒅𝑵𝑫 𝒅(𝑪𝑫 𝑽) 𝒅(𝑪𝑫 ) 𝒅(𝑽)

(𝒓𝑫 )𝑽 = ⇒ (𝒓𝑫 )𝑽 = ⇒ (𝒓𝑫 )𝑽 = 𝑽 + 𝑪𝑫
𝒅𝒕 𝒅𝒕 𝒅𝒕 𝒅𝒕

𝒅(𝑪𝑫 ) 𝒅(𝑪𝑫 )
⇒ (𝒓𝑫 )𝑽 = 𝑽 + 𝑪𝑫 𝒗𝟎 ⇒ (𝒓𝑫 )𝑽−𝑪𝑫 𝒗𝟎 = 𝑽
𝒅𝒕 𝒅𝒕
𝒅(𝑪𝑫 ) 𝒗𝟎
⇒ = (𝒓𝑫 )−𝑪𝑫 (𝑽𝑰)
𝒅𝒕 𝑽
5) Las leyes de velocidad están dadas por:
(−𝒓𝑨 ) = 𝒌𝑪𝑨 𝑪𝑩 (VII); (−𝒓𝑩 ) = 𝒌𝑪𝑨 𝑪𝑩 (VIII); (𝒓𝑪 ) = 𝒌𝑪𝑨 𝑪𝑩 (IX) y
(𝒓𝑫 ) = 𝒌𝑪𝑨 𝑪𝑩 (X)
6) Sustituimos (VII) en (III):

𝒅(𝑪𝑨 ) 𝒗𝟎
= −𝒌𝑪𝑨 𝑪𝑩 − 𝑪𝑨 (𝑿𝑰)
𝒅𝒕 𝑽

Sustituimos (VIII) en (IV):

𝒅(𝑪𝑩 ) 𝒗𝟎
= (𝑪𝑩𝟎 − 𝑪𝑩 ) − 𝒌𝑪𝑨 𝑪𝑩 (𝑿𝑰𝑰)
𝒅𝒕 𝑽

Sustituimos (IX) en (V):

𝒅(𝑪𝑪 ) 𝒗𝟎
= 𝒌𝑪𝑨 𝑪𝑩 −𝑪𝑪 (𝑿𝑰𝑰𝑰)
𝒅𝒕 𝑽

Sustituimos (X) en (VI):

𝒅(𝑪𝑫 ) 𝒗𝟎
= 𝒌𝑪𝑨 𝑪𝑩 −𝑪𝑫 (𝑿𝑰𝑽)
𝒅𝒕 𝑽

Partiendo de las ecuaciones (XI), (XII) y (XIII) planteamos un sistema de ecuaciones

diferenciales, que al introducir los valores iniciales del problema como t=0 ; 𝐶𝐴0 =0.05
mol/dm3 𝐶𝐶 = 𝐶𝐵 = 𝐶𝐷 = 0 Y k=0,22 mol/dm3 podemos obtener una solución empleando el
Software POLYMATH.
El volumen requerido del reactor para obtener el máximo beneficio, es cuando l
Concentracion de la especie C sea máxima, y mediante el reporte proporcionado por el
software podemos determinar el tiempo en que se alcanza dicha concentracion, que
corresponde al valor de 220,29 segundos. Así sustituimos este valor en (*) y nos queda:

𝒅𝒎𝟑
𝑽 = 𝒗𝟎 𝒕 + 𝑽𝟎 ⇒ 𝑽 = 𝟎, 𝟎𝟓 ∗ 𝟐𝟐𝟎, 𝟐𝟗𝒔 + 𝟓𝒅𝒎𝟑 ⇒ 𝑽 = 𝟏𝟔, 𝟎𝟏𝟒𝟓𝒅𝒎𝟑
𝒔
 1) CSTR: se trata de un reactor donde se asume que no hay gradientes de temperatura
ni concentración, la temperatura y composición de la corriente de salida es igual a la
que se tiene dentro del reactor y se trabaja con régimen estacionario.

2) Se consideran las siguientes reacciones:

k1
A+B → C
k2
B+B → D
3) Fase líquida → Volumen constante
4) Balance molar para cada componente:
a) Para A:

EA -SA -CA+GA =AA

Pero, A no se genera ni se acumula, por lo tanto:
EA -SA-CA =0 (1)
Sabemos que EA=FA0=CA0 V0 (2), SA=FA=CAV0 (3) y CA=(-rA )V (4), sustituyendo
(2), (3) y (4) en (1)
CA0V0-CA V0 -(-rA )V=0
(-rA)V=CA0V0 -CAV0
V CA0-CA
= (5)
V0 (-rA )
V
Por otra parte se tiene el tiempo espacial τ= V (6), sustituyendo (6) en (5)
0
CA0 -CA
τ= (7)
(-rA)
b) Para B:
EB -SB -CB+GB=AB
Pero B no se genera ni se acumula, por lo tanto:
EB -SB -CB=0 (8)
Sabemos que EB=FB0 =CB0V0 (9), SB=FB =CBV0 (10) y CB =(-rB)V (11), sustituyendo
(9), (10) y (11) en (8)
CB0 V0 -CBV0-(-rB)V=0
(-rB)V=CB0 V0 -CBV0
 V CB0 -CB
= (12)
V0 (-rB )
Sustituyendo (6) en (12)
CB0-CB
τ= (13)
(-rB)
c) Para C:
EC-SC -CC+GC=AC
Pero C no entra al reactor, no se consume y no se acumula, por lo tanto:
GC=SC (14)
Sabemos que GC =rC V (15) y SC =FC =CC V0 (16)
rCV=CCV0
V CC
= (17)
V0 rC
Sustituyendo (6) en (17)
CC
τ= (18)
rC
d) Para D:
ED-SD-CD+GD=AD
Pero D no entra al reactor, no se consume y no se acumula, por lo tanto:
GD=SD (19)
Sabemos que GD=rDV (20) y SD=FD=CDV0 (21)
rDV=CD V0
V CD
= (22)
V0 rD
Sustituyendo (6) en (22)
CD
τ= (23)
rD

5) Cinética de la reacción, asumiendo reacción elemental

(-rA)=k1 CACB=rC (24)

(-rB)1 rC
=
1 1
(-rB)1=rC =k1 CACB (25)
(-rB)2 rD
=
2 1
(-rB)2=2rD (26)
rD=k2 CB 2 (27)
Sustituyendo (27) en (26)
(-rB)2=2k 2CB 2 (28)
(-rB)=(-rB)1 +(-rB)2 (29)
Sustituyendo (25) y (28) en (29)
(-rB)=k 1CA CB +2k 2CB2 (30)
 6) Involucrando los pasos 4 y 5

Sustituyendo (24) en (7)

CA0-CA
τ= (31)
k 1CA CB
Sustituyendo (30) en (13)
CB0-CB
τ= (32)
k 1 CA CB +2k 2CB2
Sustituyendo (24) en (18)
CC
τ= (33)
k 1CA CB
Sustituyendo (27) en (23)
CD
τ= (34)
k 2CB 2
Todos los τ son iguales, puesto que el volumen de reactor es el mismo y la
alimentación también, así, se procede a igualar las ecuaciones necesarias. Igualando (31)
con (32)
CA0-CA CB0 -CB
=
k 1CACB k 1CACB +2k2 CB2
CB0-CB k 1CACB +2k2CB2
=
CA0-CA k 1CACB
CB0 -CB k 1 CA+2k2 CB
=
CA0 -CA k 1CA
(k 1CA+2k2 CB )(CA0 -CA)=k1 CA (CB0-CB )
k 1CA0CA+2k2 CA0 CB -k1CA2 -2k2 CACB=k1 CB0CA -k 1CACB
CB (2k2 CA0 + k1 CA-2k2CA)=k 1CA2
k 1CA2
CB= (35)
2k2 CA0 + k1 CA -2k2 CA
Igualando (31) y (33)
CA0-CA CC
=
k 1CACB k 1CACB
CC=CA0 -CA (36)
Igualando (32) y (33)
CB0-CB CC
=
k 1CACB +2k 2CB2 k 1 CACB
CB0-CB k 1CACB+2k2 CB 2
=
CC k 1 CACB
k 1CA(CA0 -CB)
CC= (37)
k 1CA +2k 2CB
Sustituyo (35) en (37)
 k C 2
k 1CA (CA0- [2k C + 1k AC -2k C ])
2 A0 1 A 2 A
CC=
k C 2
k 1CA +2k2 [2k C + 1k AC -2k C ]
2 A0 1 A 2 A
k 1 CA 2
(CA0- [
2k2 CA0 + k1 CA-2k 2CA])
CC=
2k C
1+ [2k C + k2 CA -2k C ]
2 A0 1 A 2 A
CA0(2k2 CA0+ k 1CA -2k2 CA)-k1CA2
2k 2CA0+ k1CA-2k 2CA
CC=
2k 2CA0+ k 1CA-2k2 CA +2k 2CA
2k 2CA0+ k1CA-2k 2CA
2k2 CA02 +k 1CA0CA-2k2 CA0CA -k1CA
CC =
2k 2CA0+ k 1CA
2k 2CA0(CA0-CA)+k 1CA (CA0-1)
CC= (38)
2k2 CA0 + k1 CA
F
Pero, CC = VC (39), sustituyendo (39) en (38)
0
FC 2k2 CA0 (CA0 -CA)+k1 CA(CA0 -1)
=
V0 2k 2CA0+ k 1CA
V0(2k2 CA0 (CA0 -CA)+k1 CA(CA0 -1))
FC= (40)
2k2 CA0 + k1 CA
7) Costo del proceso
$ $ $
Costo total=FA0 *100 +FB*100 +2 *V (41)
mol mol L*h R
8) Algoritmo para el tanteo
a) Supongo 0<CA <0,1y encuentro V0 en (36)
b) Sustituyo CA en (35) y encuentro CB
c) Con V0 y CB encuentro FB
d) Con V0 y CA obtengo FA
e) Sustituyo CA en (36) y encuentro CC
f) Sustituyo CA, CB y CC en (33) para obtener τ
g) Con V0 y τ encuentro el volumen del reactor
h) Tanteo hasta que el costo sea el mínimo, es decir, hasta que el volumen del
reactor sea lo menor posible.
 1. Dibujar el esquema del reactor a trabajar, y escribir la ecuación de reacción
que ocurre en dicho reactor:

2.

𝒌𝟏 𝒌𝟐
𝑨→𝑩→𝑪

3. Realizamos el balance general molar del compuesto A:

𝑬−𝑺−𝑪+𝑮=𝑨 ;
𝒅𝑵𝑨
Donde: C: −𝜸𝑨 ∗ 𝑽 y A:
𝒅𝒕

Así:
𝒅𝑵𝑨
−(−𝜸𝑨 ) ∗ 𝑽 =
𝒅𝒕
𝒅𝑪𝑨
⟹ (−𝜸𝑨 ) = − ; (1)
𝒅𝒕

Realizamos el balance general molar del compuesto B:

𝑬−𝑺−𝑪+𝑮=𝑨 ;
𝒅𝑵𝑩
Donde: C: −𝜸𝑩 ∗ 𝑽 ; G: 𝜸𝑩 ∗ 𝑽 y A:
𝒅𝒕
Así:
𝒅𝑵𝑩
−(−𝜸𝑩 ) ∗ 𝑽 + 𝜸𝑩 ∗ 𝑽 =
𝒅𝒕
𝒅𝑵𝑩
⟹ 𝑽 ∗ (−(−𝜸𝑩 ) + 𝜸𝑩 ) =
𝒅𝒕
𝒅𝑪𝑩
⟹ −(−𝜸𝑩 ) + 𝜸𝑩 = ; (2)
𝒅𝒕
Realizamos el balance general molar del compuesto C:
𝑬−𝑺−𝑪+𝑮=𝑨 ;
𝒅𝑵𝑪
Donde: G: 𝜸𝑪 ∗ 𝑽 y A:
𝒅𝒕
Así:
 𝒅𝑵𝑪
𝜸𝑪 ∗ 𝑽 =
𝒅𝒕
𝒅𝑵𝑪
𝜸𝑪 =
𝒅𝒕𝑽
𝒅𝑪𝑪
𝜸𝑪 = ; (3)
𝒅𝒕

4. Determinamos la ley de velocidad en términos de la concentración del

reactivo A, asumiendo un orden de reacción 1:
(−𝜸𝑨 ) = 𝒌𝟏 ∗ 𝑪𝑨 𝒏 ; Donde n: orden de reacción
⟹ (−𝜸𝑨 ) = 𝒌𝟏 ∗ 𝑪𝑨 ; (4)
Determinamos la ley de velocidad en términos de la concentración del compuesto B,
asumiendo un orden de reacción 1:
𝜸𝑩 = 𝒌𝟏 ∗ 𝑪𝑨 𝒏
(−𝜸𝑩 ) = 𝒌𝟐 ∗ 𝑪𝑩 𝒏 ; Donde n: orden de reacción
⟹ 𝜸𝑩 = 𝒌𝟏 ∗ 𝑪𝑨 ; (5)
⟹ (−𝜸𝑩 ) = 𝒌𝟐 ∗ 𝑪𝑩 ; (6)
Determinamos la ley de velocidad en términos de la concentración del compuesto C,
asumiendo un orden de reacción 1:
(𝜸𝑪 ) = 𝒌𝟐 ∗ 𝑪𝑩 ; (7)

5. Consideramos fase líquida a V ctte.

Sustituimos (4) en (1):
𝒅𝑪𝑨
𝒌𝟏 ∗ 𝑪𝑨 = −
𝒅𝒕
𝒕𝑹 𝑪𝑨
⟹ 𝒌𝟏 ∗ ∫ 𝒅𝒕 = − ∫ 𝒅𝑪𝑨
𝟎 𝑪𝑨𝒐
𝑪𝑨𝒐
⟹ 𝒌 ∗ 𝒕𝑹 = 𝑳𝒏 ( )
𝑪𝑨
∴ 𝑪𝑨 = 𝑪𝑨𝒐 𝒆−𝒌𝟏 𝒕𝑹 ; (8)
Luego, sustituimos (5) y (6) en (2):
𝒅𝑪𝑩
𝒌𝟏 ∗ 𝑪𝑨 − 𝒌𝟐 ∗ 𝑪𝑩 =
𝒅𝒕
Sustituimos (8) en la expresión anterior:
𝒅𝑪𝑩
⟹ 𝒌𝟏 ∗ 𝑪𝑨𝒐 𝒆−𝒌𝟏 𝒕𝑹 − 𝒌𝟐 ∗ 𝑪𝑩 =
𝒅𝒕
 𝒅𝑪𝑩
⟹ + 𝒌𝟐 ∗ 𝑪𝑩 = 𝒌𝟏 ∗ 𝑪𝑨𝒐 𝒆−𝒌𝟏 𝒕𝑹
𝒅𝒕
Aplicando factor integrante:

𝒆∫ 𝑷(𝒕𝑹) = 𝒆∫ 𝒌𝟐∗𝒅𝒕𝑹 = 𝒆𝒌𝟐 ∗𝒕𝑹

𝒅𝑪𝑩 𝒌 ∗𝒕
⟹ 𝒆 𝟐 𝑹 + 𝒌𝟐 ∗ 𝑪𝑩 ∗ 𝒆𝒌𝟐∗𝒕𝑹 = 𝒌𝟏 ∗ 𝑪𝑨𝒐 𝒆−𝒌𝟏𝒕𝑹 𝒆𝒌𝟐∗𝒕𝑹
𝒅𝒕
𝒅𝑪𝑩 𝒌 ∗𝒕
⟹ 𝒆 𝟐 𝑹 + 𝒌𝟐 ∗ 𝑪𝑩 ∗ 𝒆𝒌𝟐∗𝒕𝑹 = 𝒌𝟏 ∗ 𝑪𝑨𝒐 𝒆(𝒌𝟐 −𝒌𝟏)𝒕𝑹
𝒅𝒕
𝒕𝑹
𝒅 𝒌 ∗𝒕
⟹∫ (𝒆 𝟐 𝑹 ∗ 𝑪𝑩 ) = 𝒌𝟏 ∗ 𝑪𝑨𝒐 ∫ 𝒆(𝒌𝟐 −𝒌𝟏 )𝒕𝑹 𝒅𝒕
𝒅𝒕 𝟎

𝒌𝟐 ∗𝒕𝑹
𝒆(𝒌𝟐 −𝒌𝟏 )𝒕𝑹 − 𝟏
⟹𝒆 ∗ 𝑪𝑩 = 𝒌𝟏 ∗ 𝑪𝑨𝒐 ∗ ( )
(𝒌𝟐 − 𝒌𝟏 )

𝒆−𝒌𝟐 ∗𝒕𝑹 (𝒆(𝒌𝟐 −𝒌𝟏 )𝒕𝑹 − 𝟏)

⟹ 𝑪𝑩 = 𝒌𝟏 ∗ 𝑪𝑨𝒐 ∗ ( )
(𝒌𝟐 − 𝒌𝟏 )

𝒆−𝒌𝟏 ∗𝒕𝑹 − 𝒆−𝒌𝟐 ∗𝒕𝑹

∴ 𝑪𝑩 = 𝒌𝟏 ∗ 𝑪𝑨𝒐 ∗ ( ) ; (9)
(𝒌𝟐 − 𝒌𝟏 )
Ahora, para hallar el 𝑡𝑜𝑝 , derivamos (9) e igualamos la expresión obtenida a 0:

𝒅𝑪𝑩 𝒆−𝒌𝟏 ∗𝒕𝑹 − 𝒆−𝒌𝟐 ∗𝒕𝑹

= 𝒌𝟏 ∗ 𝑪𝑨𝒐 ∗ ( )=𝟎
𝒅𝒕 (𝒌𝟐 − 𝒌𝟏 )

𝒌𝟐 𝒆−𝒌𝟐 ∗𝒕𝑹 − 𝒌𝟏 𝒆−𝒌𝟏 ∗𝒕𝑹 = 𝟎

𝒌𝟐 𝒆−𝒌𝟐 ∗𝒕𝑹 = 𝒌𝟏 𝒆−𝒌𝟏 ∗𝒕𝑹
Aplicando Logaritmo neperiano ambos lados:
𝐥𝐧(𝒌𝟐 ) −𝒌𝟐 ∗ 𝒕𝑹 = 𝐥𝐧(𝒌𝟏 ) −𝒌𝟏 ∗ 𝒕𝑹
𝐥𝐧(𝒌𝟐 ) − 𝐥𝐧(𝒌𝟏 )
𝒕𝒐𝒑 =
𝒌𝟐 − 𝒌𝟏
Sustituyendo los valores:
𝐥𝐧(𝟎. 𝟎𝟏) − 𝐥𝐧(𝟎. 𝟒)
𝒕𝒐𝒑 = = 𝟗. 𝟒𝟔𝒉
𝟎. 𝟎𝟏 − 𝟎. 𝟒
Ahora, encontramos el valor de 𝐶𝐵𝑚𝑎𝑥 , para ello evaluamos 𝑡𝑜𝑝 en la ecuación (7):
 𝒆−𝒌𝟏 ∗𝒕𝒐𝒑 − 𝒆−𝒌𝟐 ∗𝒕𝒐𝒑
𝑪𝑩𝒎𝒂𝒙 = 𝒌𝟏 ∗ 𝑪𝑨𝒐 ∗ ( )
(𝒌𝟐 − 𝒌𝟏 )
−𝟏 −𝟏 ∗𝟗.𝟒𝟔𝒉
−𝟏 )
𝒎𝒐𝒍 𝒆−𝟎.𝟒𝒉 ∗𝟗.𝟒𝟔𝒉 − 𝒆−𝟎.𝟎𝟏𝒉 𝒎𝒐𝒍
𝑪𝑩𝒎𝒂𝒙 = (𝟎. 𝟒𝒉 ∗ (𝟏. 𝟔 )∗( ) = 𝟏. 𝟒𝟔
𝒅𝒎𝟑 (𝟎. 𝟎𝟏 − 𝟎. 𝟒) 𝒅𝒎𝟑

Involucramos (7) en (3), nos queda:

𝒅𝑪𝑪
𝒌𝟐 ∗ 𝑪𝑩 = ; (10)
𝒅𝒕
Luego, de la ecuación (10) sustituimos 𝐶𝐵 por 𝐶𝐵𝑚𝑎𝑥 :
𝒅𝑪𝑪
𝒌𝟐 ∗ 𝑪𝑩𝒎𝒂𝒙 =
𝒅𝒕
𝒌𝟐 ∗ 𝑪𝑩𝒎𝒂𝒙 ∗ 𝒅𝒕 = 𝒅𝑪𝑪
Integrando ambos lados
𝒕𝒐𝒑 𝑪𝑪
𝒌𝟐 𝑪𝑩𝒎𝒂𝒙 ∫ 𝒅𝒕 = ∫ 𝒅𝑪𝑪
𝟎 𝑪𝑪𝟎

𝒕𝒐𝒑 𝑪𝑪
𝒌𝟐 𝑪𝑩𝒎𝒂𝒙 ∫ 𝒅𝒕 = ∫ 𝒅𝑪𝑪
𝟎 𝟎

𝒌𝟐 𝑪𝑩𝒎𝒂𝒙 𝒕𝒐𝒑 = 𝑪𝑪
𝒎𝒐𝒍 𝒎𝒐𝒍
𝑪𝑪 = 𝟎. 𝟎𝟏𝒉−𝟏 ∗ 𝟏. 𝟒𝟔 𝟑
∗ 𝟗. 𝟒𝟔𝒉 = 𝟎. 𝟏𝟑𝟖
𝒅𝒎 𝒅𝒎𝟑
Finalmente determinamos los costos
𝒎𝒐𝒍 $
 Costo A = 𝟏. 𝟔 𝒅𝒎𝟑 − 𝟓𝟎𝟎𝒅𝒎𝟑 ∗ 𝟏𝟎 𝒎𝒐𝒍 = 𝟖𝟎𝟎𝟎$
𝒎𝒐𝒍 $
 Costo B= 𝟏. 𝟒𝟓𝟔 𝒅𝒎𝟑 − 𝟓𝟎𝟎𝒅𝒎𝟑 ∗ 𝟓𝟎 𝒎𝒐𝒍 = 𝟑𝟔𝟒𝟎𝟎$
𝒎𝒐𝒍
 Costo C= 𝟎. 𝟏𝟑𝟖 𝒅𝒎𝟑 − 𝟓𝟎𝟎𝒅𝒎𝟑 ∗ 𝟑𝟎$[𝒆𝟎.𝟓∗𝟎.𝟏𝟑𝟖 − 𝟏] = 𝟏𝟒𝟕. 𝟖𝟕$

Así, los beneficios están dados por

Costo B – Costo A - Costo C= (36400-8000-147.87)$
Beneficio= 28252.13$
A continuación se presenta una gráfica obtenida a través del uso del Software POLYMATH
donde se graficaron las concentraciones CA,CB y Cc con respecto al tiempo.