Examen Final A
Examen Final A
Examen Final A
𝒌𝟏
𝑨 → 𝟐𝑩 𝒌𝟏 = 𝟎. 𝟐𝟗𝒉−𝟏
𝒌𝟐
𝑩 →𝑪 𝒌𝟐 = 𝟎. 𝟏𝟕𝒉−𝟏
Ecuaciones cinéticas:
𝑟1 = 𝑘1 𝐴
𝑟2 = 𝑘2 𝐵
𝑑𝐴
− = 𝑟𝐶𝐴
𝑑𝑡
𝑑𝐵
− = 𝑟𝐶𝐵
𝑑𝑡
𝑑𝐶
− = 𝑟𝐶𝐶
𝑑𝑡
En donde:
𝑟𝐶𝐴 = 𝑟1
𝑟𝐶𝐵 = 𝑟2 − 2𝑟1
𝑟𝐶𝐶 = −𝑟2
1. BATCH
Si se desea maximizar la concentración de la especie “B” en un reactor por lotes a 20 ◦C con
concentración inicial de “A” de 0.5 M, determine:
Para el componente A
𝑑𝐴
− = 𝑘1 𝐴
𝑑𝑡
𝐴 𝑡
∫ 𝑑𝐴 = −𝑘1 ∫ 𝑑𝑡
𝐴0 0
𝐴
𝑙𝑛 ( ) = −𝑘1 𝑡
𝐴0
𝑨 = 𝑨𝟎 𝒆−𝒌𝟏 𝒕
Para el componente B
𝑑𝐵
− = 𝑘2 𝐵 − 2𝑘1 𝐴
𝑑𝑡
Sustituyendo la concentración de A:
𝑑𝐵
− = 𝑘2 𝐵 − 2𝑘1 𝐴0 𝑒 −𝑘1 𝑡
𝑑𝑡
En donde:
𝑀(𝑡, 𝐵) = 𝑘2 𝐵 − 2𝑘1 𝐴0 𝑒 −𝑘1 𝑡
𝜕𝑀(𝑡, 𝐵)
𝑀𝐵 = ( ) = 𝑘2
𝜕𝐵 𝑡
y
𝑁(𝑡, 𝐵) = 1
𝜕𝑁(𝑡, 𝐵)
𝑁𝑡 = ( ) =0
𝜕𝑡 𝐵
𝜇 = 𝜇(𝑡)
𝑀𝐵 − 𝑁𝑡
= 𝑘2
𝑁
El cociente es función de una sola variable (t), por lo tanto el factor integrante será:
𝑀𝐵 −𝑁𝑡
𝑑𝑡
𝜇(𝑡) = 𝑒 ∫ 𝑁 = 𝑒 𝑘2 ∫ 𝑑𝑡 = 𝑒 𝑘2 𝑡
𝜕𝑀(𝑡, 𝐵)
𝑀𝐵 = ( ) = 𝑘2 𝑒 𝑘 2 𝑡
𝜕𝐵 𝑡
𝑁(𝑡, 𝐵) = 𝑒 𝑘2 𝑡
𝜕𝑁(𝑡, 𝐵)
𝑁𝑡 = ( ) = 𝑘2 𝑒 𝑘 2 𝑡
𝜕𝑡 𝐵
𝑀𝐵 = 𝑁𝑡
𝜕𝛹
𝑎) ( ) = 𝑀(𝑡, 𝐵)
𝜕𝑡 𝐵
𝛹(𝑡, 𝐵) = 𝛼 en donde { 𝜕𝛹
𝑏) ( ) = 𝑁(𝑡, 𝐵)
𝜕𝐵 𝑡
𝑘1
𝛹(𝑡, 𝐵) = 𝐵𝑒 𝑘2 𝑡 − 2𝐴0 ( ) 𝑒 (𝑘2 −𝑘1 )𝑡 + 𝛼(𝐵)
𝑘2 − 𝑘1
𝜕𝛹(𝑡, 𝐵)
( ) = 𝑒 𝑘2 𝑡 + 𝛼´(𝐵)
𝜕𝐵 𝑡
𝑒 𝑘2 𝑡 + 𝛼´(𝐵) = 𝑒 𝑘2 𝑡
∴ 𝛼´(𝐵) = 0
𝛼(𝐵) = ∫ 𝛼´(𝐵)𝑑𝐵
𝛼(𝐵) = 0 + 𝛼(𝐵)
𝑘1
𝛹(𝑡, 𝐵) = 𝐵𝑒 𝑘2 𝑡 − 2𝐴0 ( ) 𝑒 (𝑘2 −𝑘1 )𝑡 = −𝛼(𝐵)
𝑘2 − 𝑘1
Para obtener el valor de la contante de integración 𝛼(𝐵), se aplica como condición de frontera que,
a un tiempo cero se tiene una concentración inicial de B igual a cero 𝐵(𝑡 = 0) = 𝐵0 = 0.
𝑘1
𝛼(𝐵) = 2𝐴0 ( )
𝑘2 − 𝑘1
𝑘1 𝑘1
𝛹(𝑡, 𝐵) = 𝐵𝑒 𝑘2 𝑡 − 2𝐴0 ( ) 𝑒 (𝑘2 −𝑘1 )𝑡 + 2𝐴0 ( )
𝑘2 − 𝑘1 𝑘2 − 𝑘1
𝑘1
𝛹(𝑡, 𝐵) = 𝐵𝑒 𝑘2 𝑡 + 2𝐴0 ( ) [1 − 𝑒 (𝑘2 −𝑘1 )𝑡 ]
𝑘2 − 𝑘1
𝒌𝟏
𝑩 = 𝟐𝑨𝟎 ( ) (𝒆−𝒌𝟏𝒕 − 𝒆−𝒌𝟐𝒕 )
𝒌𝟐 − 𝒌𝟏
Para obtener el tiempo requerido para obtener la máxima concentración de B usamos la siguiente
ecuación:
𝑘1
𝐵 = 2𝐴0 ( ) (𝑒 −𝑘1 𝑡 − 𝑒 −𝑘2 𝑡 )
𝑘2 − 𝑘1
En solver se utiliza como celda objetivo la ecuación obtenida con un valor de cero y variando t:
Objetivo -3.73E-09
t (h) 4.451
𝑑2 𝐵 𝑘1
2
= 2𝐴0 ( ) (𝑘1 2 𝑒 −𝑘1 𝑡 − 𝑘2 2 𝑒 −𝑘2 𝑡 )
𝑑𝑡 𝑘2 − 𝑘1
𝑑2 𝐵 0.29ℎ−1 −1 −1
2
= (2)(0.5𝑀) ( −1 ) [(0.29ℎ−1 )2 𝑒 −(0.29ℎ )(4.451ℎ) − (0.17ℎ−1 )2 𝑒 −(0.17ℎ )(4.451ℎ) ]
𝑑𝑡 0.17ℎ
𝑑2 𝐵
= −0.023
𝑑𝑡 2
Ya que el valor de la segunda derivada es negativo (𝐵´´ < 0) podemos concluir que se trata de un
máximo. Por lo tanto el tiempo que se requiere para maximizar la concentración de la especie B es:
𝒕 = 𝟒. 𝟒𝟓𝒉
𝑘1
𝐵 = 2𝐴0 ( ) (𝑒 −𝑘1 𝑡 − 𝑒 −𝑘2 𝑡 )
𝑘2 − 𝑘1
0.29ℎ−1 −1 −1
𝐵 = (2)(0.5𝑀) ( ) [𝑒 −(0.29ℎ )(4.451ℎ) − 𝑒 −(0.17ℎ )(4.451ℎ) ]
0.17ℎ−1
𝑄(𝐴0 − 𝐴) = 𝑉𝑟𝐶𝐴
𝑉
(𝐴0 − 𝐴) = 𝑘 𝐴
𝑄 1
Despejando A se obtiene:
𝐴0
𝐴=
1 + 𝜏𝑘1
𝟎. 𝟓𝑴
𝑨= = 𝟎. 𝟏𝟓𝟏𝑴
𝟏+ (𝟖𝒉)(𝟎. 𝟐𝟗𝒉−𝟏 )
𝑄(𝐵0 − 𝐵) = 𝑉𝑟𝐶𝐵
𝑉
(𝐵0 − 𝐵) = (𝑘 𝐵 − 2𝐾1 𝐴)
𝑄 2
Sustituyendo la concentración de A:
𝐴0
(𝐵0 − 𝐵) = 𝜏 [𝑘2 𝐵 − 2𝐾1 ( )]
1 + 𝜏𝑘1
𝐴0
−𝐵 = 𝜏 [𝑘2 𝐵 − 2𝐾1 ( )]
1 + 𝜏𝑘1
Despejando B:
𝐴0
2𝜏𝐾1 ( ) 2𝜏𝑘1 𝐴0
1 + 𝜏𝑘1
𝐵= =
1 + 𝜏𝑘2 (1 + 𝜏𝑘2 )(1 + 𝜏𝑘1 )
𝑄(𝐶0 − 𝐶) = 𝑉𝑟𝐶𝐶
𝑉
(𝐶0 − 𝐶) = − 𝑘 𝐵
𝑄 2
Sustituyendo la concentración de B:
𝑉 2𝜏𝑘1 𝐴0
(𝐶0 − 𝐶) = − 𝑘2 [ ]
𝑄 (1 + 𝜏𝑘2 )(1 + 𝜏𝑘1 )
Sabemos que la concentración inicial de C es cero (𝐶0 = 0), por lo tanto la concentración de C estará
dada por:
2𝜏 2 𝑘1 𝑘2 𝐴0
𝐶=
(1 + 𝜏𝑘2 )(1 + 𝜏𝑘1 )
𝑑𝐴 1
− = 𝑟
𝑑𝐿 𝜈 𝐶𝐴
𝑑𝐴 1
− = 𝑘 𝐴
𝑑𝐿 𝜈 1
𝐴 𝐿
𝑑𝐴 1
∫ = − 𝑘1 ∫ 𝑑𝐿
𝐴0 𝐴 𝜈 0
𝐴 1
𝑙𝑛 ( ) = − 𝑘1 𝐿
𝐴0 𝜈
Despejando A:
1
𝐴 = 𝐴0 𝑒 −𝜈𝑘1 𝐿
𝑑𝐵 1
− = 𝑟
𝑑𝐿 𝜈 𝐶𝐵
𝑑𝐵 1
− = (𝑘 𝐵 − 2𝑘1 𝐴)
𝑑𝐿 𝜈 2
Sustituyendo la concentración de A:
𝑑𝐵 1 1
− = (𝑘2 𝐵 − 2𝑘1 𝐴0 𝑒 −𝜈𝑘1 𝐿 )
𝑑𝐿 𝜈
En donde:
1 1
𝑀(𝐿, 𝐵) = (𝑘2 𝐵 − 2𝑘1 𝐴0 𝑒 −𝜈𝑘1 𝐿 )
𝜈
𝜕𝑀(𝐿, 𝐵) 𝑘2
𝑀𝐵 = ( ) =
𝜕𝐵 𝐿
𝜈
y
𝑁(𝐿, 𝐵) = 1
𝜕𝑁(𝐿, 𝐵)
𝑁𝐿 = ( ) =0
𝜕𝐿 𝐵
𝜇 = 𝜇(𝐿)
𝑀𝐵 − 𝑁𝐿 𝑘2
=
𝑁 𝜈
El cociente es función de una sola variable (L), por lo tanto el factor integrante será:
𝑀𝐵 −𝑁𝐿 𝑘2 𝑘2
𝑑𝐿 ∫ 𝑑𝐿
𝜇(𝐿) = 𝑒 ∫ 𝑁 =𝑒𝜈 =𝑒𝜈𝐿
𝑘2 𝑘 2 𝐿 𝑘1 𝑘 −𝑘
( 2 1 )𝐿
𝑘2
( 𝐵𝑒 𝜈 − 2 𝐴0 𝑒 𝜈 ) 𝑑𝐿+𝑒 𝜈 𝐿 𝑑𝐵 = 0
𝜈 𝜈
𝑘2 𝑘2 𝐿 𝑘1 𝑘 −𝑘
( 2 1 )𝐿
𝑀(𝐿, 𝐵) = 𝐵𝑒 𝜈 − 2 𝐴0 𝑒 𝜈
𝜈 𝜈
𝜕𝑀(𝐿, 𝐵) 𝑘2 𝑘2
𝑀𝐵 = ( ) = 𝑒𝜈𝐿
𝜕𝐵 𝐿
𝜈
𝑘2
𝑁(𝐿, 𝐵) = 𝑒 𝜈 𝐿
𝜕𝑁(𝐿, 𝐵) 𝑘2
𝑁𝐿 = ( ) = 𝑒 𝑘2 𝐿
𝜕𝐿 𝐵
𝜈
𝑀𝐵 = 𝑁𝐿
𝜕𝛹
𝑎) ( ) = 𝑀(𝐿, 𝐵)
𝜕𝐿 𝐵
𝛹(𝐿, 𝐵) = 𝛼 en donde { 𝜕𝛹
𝑏) ( ) = 𝑁(𝐿, 𝐵)
𝜕𝐵 𝐿
1 𝑘2 𝑘 −𝑘
( 2 1 )𝐿
𝛹(𝐿, 𝐵) = ∫ 𝑀(𝐿, 𝐵)𝑑𝐿 = ∫ [ (𝑘2 𝐵𝑒 𝜈 𝐿 − 2𝑘1 𝐴0 𝑒 𝜈 )] 𝑑𝐿
𝜈
𝑘2 𝑘1 𝑘 −𝑘
( 2 1 )𝐿
𝛹(𝐿, 𝐵) = 𝐵𝑒 𝜈 𝐿 − 2𝐴0 ( )𝑒 𝜈 + 𝛼(𝐵)
𝑘2 − 𝑘1
𝜕𝛹(𝐿, 𝐵) 𝑘2
( ) = 𝑒 𝜈 𝐿 + 𝛼´(𝐵)
𝜕𝐵 𝐿
𝜕𝛹(𝐿, 𝐵)
( ) = 𝑁(𝐿, 𝐵)
𝜕𝐵 𝐿
𝑘2 𝑘2
𝑒 𝜈 𝐿 + 𝛼´(𝐵) = 𝑒 𝜈 𝑡
∴ 𝛼´(𝐵) = 0
𝛼(𝐵) = ∫ 𝛼´(𝐵)𝑑𝐵
𝛼 = 0 + 𝛼(𝐵)
Para obtener el valor de la contante de integración 𝛼(𝐵), se aplica como condición de frontera que,
a una longitud de cero se tiene una concentración inicial de B igual a cero.
𝑘1
𝛼(𝐵) = 2𝐴0 ( )
𝑘2 − 𝑘1
𝑘2 𝑘1 𝑘 −𝑘
( 2 1 )𝐿 𝑘1
𝛹(𝐿, 𝐵) = 𝐵𝑒 𝜈 𝐿 − 2𝐴0 ( )𝑒 𝜈 + 2𝐴0 ( )
𝑘2 − 𝑘1 𝑘2 − 𝑘1
𝑘2 𝑘1 𝑘 −𝑘
( 2 1 )𝐿
𝛹(𝐿, 𝐵) = 𝐵𝑒 𝜈 𝐿 + 2𝐴0 ( ) [1 − 𝑒 𝜈 ]
𝑘2 − 𝑘1
𝑘1 𝑘1 𝑘2
𝐵 = 2𝐴0 ( ) (𝑒 − 𝜈 𝐿 − 𝑒 − 𝜈 𝐿 )
𝑘2 − 𝑘1
𝑄 = 𝜈𝐴 𝑇
𝑃𝐵 = 𝑄𝐵
4. Ajuste de CSTR
Con el fin de optimizar el proceso, se realizan pruebas en un reactor continuo perfectamente agitado
que opera en régimen permanente a 30 ◦C con concentración de entrada de “A” de 0.5 M, y se
miden las concentraciones variando el tiempo de residencia. Los resultados obtenidos se muestran
a continuación:
𝑄(𝐴0 − 𝐴) = 𝑉𝑘1 𝐴
Despejando 𝑘1 :
𝐴0 − 𝐴
𝑘1 =
𝐴𝜏
𝐸𝑎 1
𝑙𝑛(𝑘120°𝐶 ) = −𝑙𝑛(𝐴) − ( )
𝑅 𝑇1
𝐸𝑎 1
𝑙𝑛(𝑘130°𝐶 ) = −𝑙𝑛(𝐴) − ( )
𝑅 𝑇2
𝑘120°𝐶 𝐸𝑎 1 1
𝑙𝑛 ( )=− ( − )
𝑘130°𝐶 𝑅 𝑇2 𝑇1
Despejando 𝐸𝑎 :
𝑘120°𝐶
𝑙𝑛 ( )𝑅
𝑘130°𝐶
𝐸𝑎 =
1 1
(𝑇 − 𝑇 )
1 2
𝟎. 𝟐𝟗𝒉−𝟏 𝒄𝒂𝒍
𝒍𝒏 ( ) (𝟏. 𝟗𝟖𝟕𝟐 ) 𝒄𝒂𝒍
𝟎. 𝟓𝟕𝟕𝒉−𝟏 𝒎𝒐𝒍 𝑲
𝑬𝒂 = = −𝟏𝟐𝟏𝟒𝟗. 𝟑𝟓
𝟏 𝟏 𝒎𝒐𝒍
( − )
𝟐𝟗𝟑. 𝟏𝟓𝑲 𝟑𝟎𝟑. 𝟏𝟓𝑲
(𝐶0 − 𝐶) = −𝜏𝑘2 𝐵
Seleccionamos el balance en términos de C, ya que es una ecuación más sencilla y tenemos datos
de las concentraciones y del tiempo de residencia.
La concentración inicial de C es cero (𝐶0 = 0). Por lo tanto despejando 𝑘2 se obtiene:
𝐶
𝑘2 =
𝜏𝐵
Mediante los datos proporcionados se obtienen los valores de 𝑘2 para cada tiempo de residencia y
posteriormente se obtiene un promedio de todos los valores:
Ƭ B C k2 (C)
0 - - -
0.5 0.2 0.029 0.2900
1 0.279 0.081 0.2903
1.5 0.331 0.144 0.2900
2 0.334 0.193 0.2889
2.5 0.35 0.253 0.2891
3 0.333 0.289 0.2893
Promedio 0.2896
𝐸𝑎 1
𝑙𝑛(𝑘220°𝐶 ) = −𝑙𝑛(𝐴) − ( )
𝑅 𝑇1
𝐸𝑎 1
𝑙𝑛(𝑘230°𝐶 ) = −𝑙𝑛(𝐴) − ( )
𝑅 𝑇2
𝑘220°𝐶 𝐸𝑎 1 1
𝑙𝑛 ( )=− ( − )
𝑘230°𝐶 𝑅 𝑇2 𝑇1
Despejando 𝐸𝑎 :
𝑘220°𝐶
𝑙𝑛 ( )𝑅
𝑘230°𝐶
𝐸𝑎 =
1 1
(𝑇 − 𝑇 )
1 2
𝟎. 𝟏𝟕𝒉−𝟏 𝒄𝒂𝒍
𝒍𝒏 ( −𝟏 ) (𝟏. 𝟗𝟖𝟕𝟐 ) 𝒄𝒂𝒍
𝟎. 𝟐𝟖𝟗𝟔𝒉 𝒎𝒐𝒍 𝑲
𝑬𝒂 = = −𝟗𝟒𝟎𝟕. 𝟒𝟗
𝟏 𝟏 𝒎𝒐𝒍
( − )
𝟐𝟗𝟑. 𝟏𝟓𝑲 𝟑𝟎𝟑. 𝟏𝟓𝑲