“Año de la unidad, la paz y el desarrollo”

UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS

FACULTAD DE QUÍMICA E INGENIERÍA QUÍMICA
Escuela Académico Profesional de Ingeniería Química

“TUTORIAL 5 – MODELOS DE REACCIONES LINEALES COMPLEJAS”

CURSO: Diseño de Reactores

ALUMNOS:

Caballero Remache, Carlos Gonzalo

Castillón Flores, Jhover Alfredo

Lucana Florez, Joel Fernando

Roldan Sagón, Joselyn Milagros Angie

Tomasto Cristobal, Piero Leonardo

PROFESOR: Rubén Alfredo Palomino Infante

SEMESTRE ACADÉMICO 2023-II

Ciudad Universitaria
Lima-Perú
2023
 RESOLUCIÓN DE PROBLEMAS

Asuma valores para las K ij y obtenga los perfiles correspondientes a las

trayectorias de las concentraciones en función del tiempo.
Considere que inicialmente solo tiene la especie A cuya concentración es unitaria
en unidades arbitrarias.

PROBLEMA 1

SOLUCIÓN:
Derivamos la ecuación de la velocidad:
ⅆ𝐶𝐴
= 𝑘2 𝐶𝐵 − 𝐶𝐴 (𝑘1 + 𝑘9 )
ⅆ𝑡
ⅆ𝐶𝐵
= 𝑘1 𝐶𝐴 + 𝑘4 𝐶𝐶 + 𝑘6 𝐶𝐷 − (𝑘2 + 𝑘3 + 𝑘5 )𝐶𝐵
ⅆ𝑡
ⅆ𝐶𝐶
= 𝑘3 𝐶𝐵 − (𝑘4 + 𝑘8 )𝐶𝐶
ⅆ𝑡
ⅆ𝐶𝐷
= 𝑘9 𝐶𝐴 + 𝑘5 𝐶𝐵 − (𝑘6 + 𝑘7 )𝐶𝐷
ⅆ𝑡
ⅆ𝐶𝐸
= 𝑘8 𝐶𝐶 + 𝑘7 𝐶𝐷
ⅆ𝑡

ⅆ𝐶𝐴
ⅆ𝑡
ⅆ𝐶𝐵
−(𝑘1 + 𝑘9 ) 𝑘2 0 0 0 𝐶𝐴
ⅆ𝑡 𝑘4 𝑘6 𝐶𝐵
ⅆ𝐶𝐶 𝑘1 −(𝑘2 + 𝑘3 + 𝑘5 ) 0
= 0 𝑘3 −(𝑘4 + 𝑘8 ) 0 0 ∙ 𝐶𝐶
ⅆ𝑡 0 −(𝑘6 + 𝑘7 ) 0 𝐶𝐷
𝑘9 𝑘5
ⅆ𝐶𝐷
( 0 0 𝑘8 𝑘 7 0) (𝐶𝐸 )
ⅆ𝑡
ⅆ𝐶𝐸
( ⅆ𝑡 )
Donde A es la matriz de coeficientes cinéticos:
−(𝑘1 + 𝑘9 ) 𝑘2 0 0 0
𝑘1 −(𝑘2 + 𝑘3 + 𝑘5 ) 𝑘4 𝑘6 0
𝐴= 0 𝑘3 −(𝑘4 + 𝑘8 ) 0 0
𝑘9 𝑘5 0 −(𝑘6 + 𝑘7 ) 0
( 0 0 𝑘8 𝑘7 0)
Suponiendo valores de 𝑘𝑖
𝑘1 = 0.00223
𝑘2 = 0.00125
𝑘3 = 0.00128
𝑘4 = 0.00367
𝑘5 = 0.00055
𝑘6 = 0.00423
𝑘7 = 0.00321
𝑘8 = 0.00112
𝑘9 = 0.00187
Utilizando el método de Rungen-Kutta, el perfil de concentraciones está dado por:
 PROBLEMA 2

SOLUCIÓN:
Derivamos la ecuación de la velocidad:
ⅆ𝐶𝐴
= −𝑘1 𝐶𝐴 + 𝑘2 𝐶𝐵
ⅆ𝑡
ⅆ𝐶𝐵
= 𝑘1 𝐶𝐴 − (𝑘2 + 𝑘3 + 𝑘5 )𝐶𝐵 + 𝑘6 𝐶𝐷 + 𝑘4 𝐶𝐶 + 𝑘9 𝐶𝐸
ⅆ𝑡
ⅆ𝐶𝐶
= 𝑘3 𝐶𝐵 − (𝑘4 + 𝑘8 )𝐶𝐶
ⅆ𝑡
ⅆ𝐶𝐷
= −(𝑘6 + 𝑘7 )𝐶𝐷 + 𝑘5 𝐶𝐵
ⅆ𝑡
ⅆ𝐶𝐸
= −𝑘9 𝐶𝐸 − 𝑘7 𝐶𝐷 + 𝑘8 𝐶𝐶
ⅆ𝑡
Matriz de coeficientes cinéticos:
−𝑘1 𝑘2 0 0 0
𝑘1 −(𝑘2 + 𝑘3 + 𝑘5 ) 𝑘4 𝑘6 𝑘9
𝐴= 0 𝑘3 −(𝑘4 + 𝑘8 ) 1 1
0 𝑘5 0 −(𝑘6 + 𝑘7 ) 1
[ 0 0 𝑘6 𝑘7 −𝑘9 ]

Suponemos valores:
𝑘1 = 1, 𝑘2 = 2, 𝑘3 = 3, 𝑘4 = 1, 𝑘5 = 2, 𝑘6 = 1, 𝑘7 = 2, 𝑘8 = 1, 𝑘9 = 2, 𝑘1 0 = 1

−1 2 0 0 0
1 −7 1 1 2
𝐶= 0 3 −2 1 1
0 2 0 −3 1
[0 0 1 2 −2]
Establecemos la matriz identidad y variable

1 0 0 0 0 λ 0 0 0 0
0 1 0 0 0 0 λ 0 0 0
𝐼= 0 0 1 0 0 → Iλ = 0 0 λ 0 0
0 0 0 1 0 0 0 0 λ 0
[0 0 0 0 1] [ 0 0 0 0 λ]

Eigenvalores
(𝐶 − 𝜆𝐼) = 0 … (1)
Reemplazando en (1) tenemos:
−1 − λ 2 0 0 0
1 −7 − λ 1 1 2
[(𝐶 − 𝜆𝐼)] = 0 3 −2 − λ 1 1 =0
0 2 0 −3 − λ 1
[ 0 0 1 2 −2 − λ]

Finalmente, se determinó la determinante de la matriz:

−λ5 − 15λ4 − 69λ3 − 99λ2 − 9λ + 47 = 0
Se encontró una solución utilizando el método de Newton – Rapshon:
Eigenvectores:
𝜆1 = 0.55, 𝜆2 = −1.22, 𝜆3 = 18.92, 𝜆4 = 0, 𝜆5 = 0,
Para comprobar la solución se procedió a hacer la gráfica de la ecuación con la finalidad
de encontrar los valores interceptan con el eje “X” y hacen 0, la ecuación.
 Gráfica 1. Determinante
Del gráfico podemos afirmar con certeza que solo existen 3 soluciones que hacen 0 la
ecuación.
De la siguiente expresión
𝐴𝑋 = 𝜆𝑖 𝑋
• Para 𝜆1 = 0.55
−1 2 0 0 0 𝑥1 𝑥1
1 −7 1 1 2 𝑥2 𝑥2
0 3 −2 1 1 𝑥3 = 0.55 𝑥3
0 2 0 −3 1 𝑥4 𝑥4
[0 0 1 [
2 −2] 5𝑥 ] [ 𝑥5 ]

−𝑥1 2𝑥2 0 0 0 0.55𝑥1

𝑥1 −7𝑥2 𝑥3 𝑥4 2𝑥5 0.55𝑥2
0 3𝑥2 −2𝑥3 𝑥4 𝑥5 = 0.55𝑥3
0 2𝑥2 0 −3𝑥4 𝑥5 0.55𝑥4
[ 0 0 𝑥3 2𝑥4 −2𝑥5 ] [0.55𝑥5 ]
De la matriz obtenemos el siguiente sistema de ecuaciones:
−𝑥1 + 2𝑥2 = 0.55𝑥1
𝑥1 − 7𝑥2 + 𝑥3 + 𝑥4 + 2𝑥5 = 0.55𝑥2
3𝑥2 − 2𝑥3 + 𝑥4 + 𝑥5 = 0.55𝑥3
2𝑥2 − 3𝑥4 + 𝑥5 = 0.55𝑥4
𝑥3 + 2𝑥4 − 2𝑥5 = 0.55𝑥5
Resolviendo el sistema de ecuaciones, tenemos:
𝑥1 = 0, 𝑥2 = 0, 𝑥3 = 0, 𝑥4 = 0, 𝑥5 = 0,
Además, 𝑥𝑖 = 𝐶𝑗 (i= 1, 2, 3, 4, 5; j= A, B, C, D, E).

Análogamente para 𝜆2 = −1.22 𝑦 𝜆3 = 18.92, se obtuvieron los siguientes valores,

respectivamente.
• Para 𝜆2 = −1.22
𝑥1 = 0, 𝑥2 = 0, 𝑥3 = 0, 𝑥4 = 0, 𝑥5 = 0,
• Para 𝜆3 = 18.92
𝑥1 = 0, 𝑥2 = 0, 𝑥3 = 0, 𝑥4 = 0, 𝑥5 = 0,
 PROBLEMA 3

SOLUCIÓN:
Asumiendo valores para 𝑘𝑖𝑗
𝑘1 = 1 𝑘4 = 3 𝑘7 = 1
𝑘2 = 2 𝑘5 = 2 𝑘8 = 1
𝑘3 = 3 𝑘6 = 2 𝑘9 = 3
Entonces
ⅆ𝐶𝐴
= −𝑘1 𝐶𝐴 − 𝑘9 𝐶𝐴 + 𝑘2 𝐶𝐵
ⅆ𝑡

ⅆ𝐶𝐵
= −𝑘2 𝐶𝐵 − 𝑘3 𝐶𝐵 − 𝑘5 𝐶𝐵 + 𝑘1 𝐶𝐴 + 𝑘4 𝐶𝐶 + 𝑘6 𝐶𝐷
ⅆ𝑡

ⅆ𝐶𝐶
= −𝑘4 𝐶𝐶 − 𝑘8 𝐶𝐶 + 𝑘3 𝐶𝐵
ⅆ𝑡
ⅆ𝐶𝐷
= −𝑘6 𝐶𝐷 − 𝑘7 𝐶𝐷 + 𝑘9 𝐶𝐴 + 𝑘5 𝐶𝐵
ⅆ𝑡
ⅆ𝐶𝐵
= 𝑘1 𝐶𝐴 − 𝑘2 𝐶𝐵 − 𝑘3 𝐶𝐵 + 𝑘4 𝐶𝐶 − 𝑘5 𝐶𝐵 + 𝑘6 𝐶𝐷
ⅆ𝑡
ⅆ𝐶𝐶
= 𝑘3 𝐶𝐵 − 𝑘4 𝐶𝐶 − 𝑘8 𝐶𝐶
ⅆ𝑡
ⅆ𝐶𝐷
= 𝑘5 𝐶𝐵 − 𝑘6 𝐶𝐷 − 𝑘7 𝐶𝐷 + 𝑘9 𝐶𝐴
ⅆ𝑡
ⅆ𝐶𝐸
= 𝑘7 𝐶𝐷 + 𝑘8 𝐶𝐶
ⅆ𝑡
Reorganizando
ⅆ𝐶𝐴
= −(𝑘1 + 𝑘9 )𝐶𝐴 + 𝑘2 𝐶𝐵
ⅆ𝑡
 ⅆ𝐶𝐵
= 𝑘1 𝐶𝐴 − (𝑘2 + 𝑘3 + 𝑘5 )𝐶𝐵 + 𝑘4 𝐶𝐶 + 𝑘6 𝐶𝐷
ⅆ𝑡
ⅆ𝐶𝐶
= 𝑘3 𝐶𝐵 − (𝑘4 + 𝑘8 )𝐶𝐶
ⅆ𝑡
ⅆ𝐶𝐷
= 𝑘9 𝐶𝐴 + 𝑘5 𝐶𝐵 − (𝑘6 + 𝑘7 )𝐶𝐷
ⅆ𝑡
ⅆ𝐶𝐸
= 𝑘8 𝐶𝐶 + 𝑘7 𝐶𝐷
ⅆ𝑡

De forma matricial

ⅆ𝐶̅ = [𝑘𝑖𝑗 ]𝐶̅

𝐶𝐴
𝐶𝐵
𝐶̅ = 𝐶𝐶
𝐶𝐷
[𝐶𝐸 ]
−(𝑘1 + 𝑘9 ) 𝑘2 0 0 0
𝑘1 −(𝑘2 + 𝑘3 + 𝑘5 ) 𝑘 4 𝑘 6 0
𝑘𝑖𝑗 = 0 𝑘3 −(𝑘4 + 𝑘8 ) 0 0
𝑘9 𝑘5 0 −(𝑘6 + 𝑘7 ) 0
[ 0 0 𝑘8 𝑘7 0 ]
Con los valores de k iniciales

−4 2 0 0 0 𝐶𝐴
ⅆ𝐶̅ 1 −7 3 2 0 𝐶𝐵
= 0 3 −4 0 0 𝐶𝐶
ⅆ𝑡 3 2 0 −3 0 𝐶𝐷
[ 0 0 1 1 0 ] [𝐶𝐸 ]

Resolvemos la siguiente determinante

|[𝐴 − 𝜆]𝐽| = 0

−4 2 0 0 0 𝜆1 0 0 0 0 1 0 0 0 0
1 −7 3 2 0 0 𝜆2 0 0 0 0 1 0 0 0
| 0 3 −4 0 0 − 0 0 𝜆3 0 0 0 0 1 0 0 || = 0
|
3 2 0 −3 0 0 0 0 𝜆4 0 0 0 0 1 0
[
( 0 0 1 1 0 ] [ 0 0 0 0 𝜆5 ]) (0 0 0 0 1 )
−4 − 𝜆1 2 0 0 0
1 −7 − 𝜆 2 3 2 0
| 0 3 −4 − 𝜆3 0 0 || = 0
|
3 2 0 −3 − 𝜆4 0
0 0 1 1 −𝜆5
Tomando como 𝜆1 = 𝜆2 = 𝜆3 = 𝜆4 = 𝜆5 = 𝜆 para obtener un polinomio
 −𝜆5 − 18𝜆4 − 102𝜆3 − 207𝜆2 − 92𝜆 = 0
De las raíces del polinomio
𝜆1 = 0
𝜆2 = −4
𝜆3 = −0.607
𝜆4 = −4.055
𝜆5 = −9.337
De la siguiente expresión
𝐴𝑋 = 𝜆𝑖 𝑋
• Para 𝜆1 = 0

−4 2 0 0 0 𝑥1 𝑥1
1 −7 3 2 0 2𝑥 𝑥2
0 3 −4 0 0 𝑥3 = 0 𝑥3
3 2 0 −3 0 𝑥4 𝑥4
[ 0 0 1 1 0] [𝑥5 ] [𝑥5 ]
−4𝑥1 + 2𝑥2 0 −4𝑥1 + 2𝑥2 = 0
𝑥1 − 7𝑥2 + 3𝑥3 + 2𝑥4 0 𝑥1 − 7𝑥2 + 3𝑥3 + 2𝑥4 = 0
3𝑥2 − 4𝑥3 = 0 → 3𝑥2 − 4𝑥3 = 0
3𝑥1 + 2𝑥2 − 3𝑥4 0 3𝑥1 + 2𝑥2 − 3𝑥4 = 0
[ 𝑥3 + 𝑥4 ] [0] 𝑥3 + 𝑥4 = 0
Resolviendo el sistema de ecuaciones
𝑥1 =0
𝑥2 =0
𝑥3 =0
𝑥4 =0
𝑥5 =0

• Para 𝜆2 = −4, se sigue el mismo procedimiento

−4 2 0 0 0 𝑥1 𝑥1
1 −7 3 2 0 𝑥2 𝑥2
0 3 −4 0 0 𝑥3 = −4 𝑥3
3 2 0 −3 0 𝑥4 𝑥4
[ 0 0 1 1 0] [𝑥5 ] [𝑥5 ]

−4𝑥1 + 2𝑥2 −4𝑥1 2𝑥2 = 0

𝑥1 − 7𝑥2 + 3𝑥3 + 2𝑥4 −4𝑥2 𝑥1 − 3𝑥2 + 3𝑥3 + 2𝑥4 = 0
3𝑥2 − 4𝑥3 = −4𝑥3 → 3𝑥2 = 0
3𝑥1 + 2𝑥2 − 3𝑥4 −4𝑥4 3𝑥1 + 2𝑥2 + 𝑥4 = 0
[ 𝑥3 + 𝑥4 ] [−4𝑥5 ] 𝑥3 + 𝑥4 + 4𝑥5 = 0
 3𝑥4
𝑥1 = −
9
𝑥2 = 0
5𝑥5
𝑥3 = −
9
𝑥4 = ℝ
𝑥4
𝑥5 = −
9

• Para 𝜆3 = −0.607
−4 2 0 0 0 𝑥1 𝑥1
1 −7 3 2 0 𝑥2 𝑥2
0 3 −4 0 0 𝑥3 = −0.607 𝑥3
3 2 0 −3 0 𝑥4 𝑥4
[ 0 0 1 1 0] [𝑥5 ] [𝑥5 ]

−4𝑥1 + 2𝑥2 −0.607𝑥1 −3.393𝑥1 + 2𝑥2 = 0

𝑥1 − 7𝑥2 + 3𝑥3 + 2𝑥4 −0.607𝑥2 𝑥1 − 6.393𝑥2 + 3𝑥3 + 2𝑥4 = 0
3𝑥2 − 4𝑥3 = −0.607𝑥3 → 3𝑥2 − 3.393𝑥3 = 0
3𝑥1 + 2𝑥2 − 3𝑥4 −0.607𝑥4 3𝑥1 + 2𝑥2 − 2.393𝑥4 = 0
[ 𝑥3 + 𝑥4 ] [−0.607𝑥5 ] 𝑥3 + 𝑥4 + 0.6074𝑥5 = 0

𝑥1 =0
𝑥2 =0
𝑥3 =0
𝑥4 =0
𝑥5 =0

• Para 𝜆4 = −4.055
−4 2 0 0 0 𝑥1 𝑥1
1 −7 3 2 0 𝑥2 𝑥2
0 3 −4 0 0 𝑥3 = −4.055 𝑥3
3 2 0 −3 0 𝑥4 𝑥4
[ 0 0 1 1 0] [𝑥5 ] [𝑥5 ]

−4𝑥1 + 2𝑥2 −4.055𝑥1 0.055𝑥1 + 2𝑥2 = 0

𝑥1 − 7𝑥2 + 3𝑥3 + 2𝑥4 −4.055𝑥2 𝑥1 − 2.945𝑥2 + 3𝑥3 + 2𝑥4 = 0
3𝑥2 − 4𝑥3 = −4.055𝑥3 → 3𝑥2 + 0.055𝑥3 = 0
3𝑥1 + 2𝑥2 − 3𝑥4 −4.055𝑥4 3𝑥1 + 2𝑥2 + 1.055𝑥4 = 0
[ 𝑥3 + 𝑥4 ] [−4.055𝑥5 ] 𝑥3 + 𝑥4 + 4.055𝑥5 = 0

𝑥1 =0
𝑥2 =0
𝑥3 =0
𝑥4 =0
𝑥5 =0

• Para 𝜆5 = −9.337
 −4 2 0 0 0 𝑥1 𝑥1
1 −7 3 2 0 𝑥2 𝑥2
0 3 −4 0 0 𝑥3 = −9.337 𝑥3
3 2 0 −3 0 𝑥4 𝑥4
[ 0 0 1 1 0] [𝑥5 ] [𝑥5 ]

−4𝑥1 + 2𝑥2 −9.337𝑥1 5.337𝑥1 + 2𝑥2 = 0

𝑥1 − 7𝑥2 + 3𝑥3 + 2𝑥4 −9.337𝑥2 𝑥1 + 2.337𝑥2 + 3𝑥3 + 2𝑥4 = 0
3𝑥2 − 4𝑥3 = −9.337𝑥3 → 3𝑥2 + 5.337𝑥3 = 0
3𝑥1 + 2𝑥2 − 3𝑥4 −9.337𝑥4 3𝑥1 + 2𝑥2 − 2.393𝑥4 = 0
[ 𝑥3 + 𝑥4 ] [−9.337𝑥5 ] 𝑥3 + 𝑥4 + 9.337𝑥5 = 0

𝑥1 =0
𝑥2 =0
𝑥3 =0
𝑥4 =0
𝑥5 =0
 PROBLEMA 4

SOLUCIÓN:
Derivamos la ecuación de la velocidad:

−(𝑘1 + 𝑘9 ) 𝑘2 0 0 0
𝑘1 −(𝑘 2 + 𝑘 3 + 𝑘 5 ) 𝑘 4 𝑘 6 0
𝑥⃗ ′ = 0 𝑘3 −(𝑘 4 + 𝑘 8 ) 0 0 𝑥⃗
𝑘9 𝑘5 0 −(𝑘 6 + 𝑘 7 ) 0
( 0 0 𝑘8 𝑘 7 0)

Sea 𝑘1 = 𝑘2 = 𝑘3 = 𝑘4 = 𝑘5 = 𝑘6 = 𝑘7 = 1
 −2 1 0 0 0
1 −3 1 1 0
𝐴= 0 1 −2 0 0
1 1 0 −2 0
[0 0 1 1 0]
(𝐴 − 𝜆𝐼)𝑥 = 0

Eigenvalores, valores 𝜆 que deben cumplir con la siguiente condición:

→ det(𝐴 − 𝜆𝐼)𝑥 = 0
−2 1 0 0 0 1 0 0 0 0
1 −3 1 1 0 0 1 0 0 0
𝐴 − 𝜆𝐼 = 0 1 −2 0 0 −𝜆 0 0 1 0 0 =0
1 1 0 −2 0 0 0 0 1 0
[0 0 1 1 0] [0 0 0 0 1]
−2 − 𝜆 1 0 0 0
1 −3 − 𝜆 1 1 0
0 1 −2 − 𝜆 0 0 =0
1 1 0 −2 − 𝜆 0
[ 0 0 1 1 −𝜆]
Por el método de cofactores, determinamos la determinante:
−2 − 𝜆 1 0 0 0
1 −3 − 𝜆 1 1 0
0 1 −2 − 𝜆 0 0 =0
1 1 0 −2 − 𝜆 0
[ 0 0 1 1 −𝜆]
−3 − 𝜆 1 1 0 1 1 1 0
(−2 − 𝜆) [ 1 −2 − 𝜆 0 0 ] − [ 0 −2 − 𝜆 0 0 ]=0
1 0 −2 − 𝜆 0 1 0 −2 − 𝜆 0
0 1 1 −𝜆 0 1 1 −𝜆
1 1 1
−3 − 𝜆 1 1
(−2 − 𝜆)(−𝜆) [ 0 −2 −𝜆 0 ]=0
1 −2 − 𝜆 0 ]+𝜆[ 1 0 −2 − 𝜆
1 0 −2 − 𝜆

𝜆(2 + 𝜆)[(−3 − 𝜆)(−2 − 𝜆)2 − ((−2 − 𝜆) + (−2 − 𝜆))] + 𝜆[(−2 − 𝜆)2 − (−2 − 𝜆)] = 0
𝜆(2 + 𝜆)[(−3 − 𝜆)(2 + 𝜆)2 + 2(2 + 𝜆)] + 𝜆[(2 + 𝜆)2 + (2 + 𝜆)] = 0

𝜆(2 + 𝜆)[(2 + 𝜆)2 (−3 − 𝜆) + 3(2 + 𝜆) + 1] = 0

𝜆(2 + 𝜆)(𝜆3 + 7𝜆2 + 13𝜆 + 5) = 0

Eigenvalores:
𝜆1 = 0
𝜆2 = −2
𝜆3 = −0.5188
𝜆4 = −2.3111
𝜆3 = −4.1701
𝐴𝑥 = 𝜆𝑥
Determinando los eigenvectores:
Sabiendo que:
(𝐴 − 𝜆𝐼)𝑥 = 0
−2 − 𝜆 1 0 0 0 𝑥1 𝑥1
1 −3 − 𝜆 1 1 0 𝑥2 𝑥2
0 1 −2 − 𝜆 0 0 𝑥3 = 𝑥3
1 1 0 −2 − 𝜆 0 𝑥4 𝑥4
[ 0 0 1 1 −𝜆] [ 5 ] [𝑥5 ]
𝑥

Para 𝜆 = 0
−2 1 0 0 0 𝑥1 0
1 −3 1 1 0 𝑥2 0
0 1 −2 0 0 𝑥3 = 0
1 1 0 −2 0 𝑥4 0
[0 0 1 1 0] [𝑥5 ] [0]

−2𝑥1 𝑥2 0 0 0 0
𝑥1 −3𝑥2 𝑥3 𝑥4 0 0
0 𝑥2 −2𝑥3 0 0 = 0
1 𝑥2 0 −2𝑥4 0 0
[ 0 0 𝑥3 𝑥4 0] [0]
−2𝑥1 + 𝑥2 = 0(1)

𝑥1 − 3𝑥2 + 𝑥3 + 𝑥4 = 0(2)

𝑥2 − 2 𝑥3 = 0(3)
𝑥1 + 𝑥2 − 2𝑥4 = 0(4)
𝑥3 + 𝑥4 = 0(5)

Hacemos un cambio de variable:

𝑥2 = 𝑠𝑥5 = 𝑡

En (1)

−2𝑥1 + 𝑠 = 0

𝑥1 = 0.5𝑠

En (3)

𝑠 − 2𝑥3 = 0

𝑥3 = 0.5𝑠

En (4)
 0.5𝑠 + 𝑠 − 2𝑥4 = 0
𝑥4 = 0.75𝑠

NOVA: No cumple con el sistema de ecuaciones

0.5 𝑠 0.5 0
𝑠 1 0
𝑥1 = 0.5 = 𝑠 0.5 + 𝑡 0
0.75 𝑠 0.75 0
[ 𝑡 ] [ 0 ] [0]

JIVA: Cumple A𝑥1 = 𝜆1 𝑥1

0
0
𝑥1 = 0
0
[0]
Para 𝜆2 = −2
0 1 0 0 0 𝑥1 0
1 −1 1 1 0 𝑥2 0
0 1 0 0 0 𝑥3 = 0
1 1 0 0 0 𝑥4 0
[0 0 1 1 2] [𝑥5 ] [0]
0 𝑥2 0 0 0 0
𝑥1 −𝑥2 𝑥3 𝑥4 0 0
0 𝑥2 0 0 0 = 0
𝑥1 𝑥2 0 0 0 0
[0 0 𝑥3 𝑥4 2𝑥5 ] [0]
𝑥2 = 0(1)𝑥1 − 𝑥2 + 𝑥3 + 𝑥4 = 0(2)𝑥2 = 0(3)𝑥1 + 𝑥2 = 0(4)
𝑥3 + 𝑥4 + 2𝑥5 = 0(5)
Entonces:
𝑥2 = 0
𝑥1 = 0
𝑥3 = −𝑥4
𝑥5 = 0
Queda:
0 𝑥1
0 𝑥2
𝜆2 = 1 = 𝑥3
−1 𝑥4
[ 0 ] [𝑥5 ]
Para 𝜆3 = −0.5188
 −1.4812 1 0 0 0 𝑥1 0
1 −2.4812 1 1 0 𝑥2 0
0 1 −1.4812 0 0 𝑥3 = 0
1 1 0 −1.4812 0 𝑥4 0
[ 0 0 1 1 −0.5188] [𝑥5 ] [0]
−1.4812 𝑥1 𝑥2 0 0 0 0
𝑥1 −2.4812𝑥2 𝑥3 𝑥4 0 0
0 𝑥2 −1.4812𝑥3 0 0 = 0
𝑥1 𝑥2 0 −1.4812 𝑥4 0 0
[ 0 0 𝑥3 𝑥4 −0.5188 𝑥5 ] [0]

−1.4812 𝑥1 + 𝑥2 = 0(1)
𝑥1 − 2.4812𝑥2 + 𝑥3 + 𝑥4 = 0(2)
𝑥2 − 1.4812 𝑥3 = 0(3)
𝑥1 + 𝑥2 − 1.4812 𝑥4 = 0(4)
𝑥3 + 𝑥4 − 0.5188 𝑥5 = 0(5)
Entonces:
𝑥2 = 𝑠
𝑥1 = 0.675128 𝑠
𝑥3 = 0.675128 𝑠
𝑥4 = 1.13093 𝑠
𝑥5 = −3.48122 𝑠
0.675128 𝑠 0.675128 𝑥1 0.675128
𝑠 1 𝑥2 1
𝜆3 = 0.675128 𝑠 = 𝑠 0.675128 → 𝑥3 = 0.675128
1.13093 𝑠 1.13093 𝑥4 1.13093
[−3.48122 𝑠] [−3.48122] [𝑥5 ] [−3.48122]
Para 𝜆4 = −2.31111
0.3111 1 0 0 0 𝑥1 0
1 −0.6889 1 1 0 𝑥2 0
0 1 0.3111 0 0 𝑥3 = 0
1 1 0 0.3111 0 𝑥4 0
[ 0 0 1 1 ]
2.3111 5[𝑥 ] [0]
0.3111 𝑥1 𝑥2 0 0 0 0
𝑥1 −0.6889𝑥 𝑥 3 𝑥4 0 0
2
0 𝑥2 0.3111𝑥3 0 0 = 0
𝑥1 𝑥2 0 0.3111 𝑥4 0 0
[ 0 0 𝑥3 𝑥4 2.3111 𝑥5 ] [0]

0.3111 𝑥1 + 𝑥2 = 0(1)
𝑥1 − 0.6889𝑥2 + 𝑥3 + 𝑥4 = 0(2)
𝑥2 + 0.3111 𝑥3 = 0(3)
𝑥1 + 𝑥2 + 0.3111 𝑥4 = 0(4)
𝑥3 + 𝑥4 + 2.3111 𝑥5 = 0(5)
Entonces:
𝑥2 = 𝑠
−𝑠
𝑥1 =
0.31111
En (3):
−𝑠
𝑥3 =
0.31111
En (4):
−𝑠 0.31111 𝑠
+ + 0.3111 𝑥4 = 0
0.31111 0.31111
0.6889 𝑠
𝑥4 =
0.311112
En (5):
−𝑠 0.6889 𝑠
+ + 2.3111 𝑥5 = 0
0.31111 0.311112
−0.3778889 𝑠
𝑥5 =
0.311113
Queda:
−3.2144 𝑠 −3.2144 𝑥1 −3.2144
𝑠 1 𝑥2 1
𝜆4 = [−3.2144 𝑠] = 𝑠 −3.2144 → 𝑥3 = −3.2144
7.1180 𝑠 7.1180 𝑥4 7.1180
−1.6891 𝑠 [−1.6891] [𝑥5 ] [−1.6891]

Para 𝜆5 = −4.1701
2.1701 1 0 0 0 𝑥1 0
1 1.1701 1 1 0 𝑥2 0
0 1 2.1701 0 0 𝑥3 = 0
1 1 0 2.1701 0 𝑥4 0
[ 0 0 1 1 4.1701 ] [𝑥 5 ] [ 0]
2.1701 𝑥1 𝑥2 0 0 0 0
𝑥1 1.1701𝑥2 𝑥3 𝑥4 0 0
0 𝑥2 2.1701𝑥3 0 0 = 0
𝑥1 𝑥2 0 2.1701𝑥4 0 0
[ 0 0 𝑥3 𝑥4 4.1701 𝑥5 ] [0]
2.1701 𝑥1 + 𝑥2 = 0(1)
𝑥2 + 2.1701 𝑥3 = 0(3)
𝑥1 + 𝑥2 + 2.1701 𝑥4 = 0(4)
𝑥3 + 𝑥4 + 4.1701 𝑥5 = 0(5)
Entonces:
𝑥2 = 𝑠
En (1):
𝑥1 = −0.4608 𝑠
En (3):
𝑥3 = −0.46085 𝑠
En (4):
−0.4608𝑠 + 𝑠 + 2.1701𝑥4 = 0
𝑥4 = −0.2485
En (5):
−0.46085 − 0.2485𝑠 + 4.1701𝑥5 = 0
𝑥5 = 0.1701
Queda:
−0.4608 𝑠 −0.4608 𝑥1 −0.4608
𝑠 1 𝑥2 1
𝜆5 = [−0.4608 𝑠] = 𝑠 −0.4608 → 𝑥3 = −0.4608
−0.2485 𝑠 −0.2485 𝑥4 −0.2485
0.1701 𝑠 [ 0.1701 ] [𝑥5 ] [ 0.1701 ]
Usando:
𝐶𝐴
𝐶𝐵
𝑌(𝑡) = 𝐶𝐶 = ∑ 𝐾𝑖 . 𝑋𝑖 . 𝑒 𝜆𝑖.𝑡 ; 𝑖 = 5
𝐶𝐷
[𝐶𝐸 ]
Para 𝑡 = 0, 𝐶𝐶 = 1 ⋀ 𝐶𝐵 = 𝐶𝐴 = 𝐶𝐷 = 𝐶𝐸 = 0
0 0 0 0.675128 −3.2144 −0.4608
0 0 0 1 1 1
1 = 𝑘1 . 0 + 𝑘2 . 1 + 𝑘3 . 0.675128 + 𝑘4 . −3.2144 + 𝑘5 . −0.4608
0 0 −1 1.13093 7.1180 −0.2485
[0] [1] [1] [−3.48122] [−1.6891] [ 0.1701 ]
Obtenemos el siguiente sistema de ecuaciones:
0 = 0. 𝑘1 + 0. 𝑘2 + 0.6751. 𝑘3 − 3.2144. 𝑘4 − 0.4608. 𝑘5
0 = 0. 𝑘1 + 0. 𝑘2 + 1. 𝑘3 + 1. 𝑘4 + 1. 𝑘5
0 = 0. 𝑘1 + 1. 𝑘2 + 0.6751. 𝑘3 − 3.2144. 𝑘4 − 0.4608. 𝑘5
0 = 0. 𝑘1 − 1. 𝑘2 + 1.13093. 𝑘3 + 7.1180. 𝑘4 − 0.2485. 𝑘5
{0 = 1. 𝑘1 + 1. 𝑘2 − 3.4812. 𝑘3 − 1.6891. 𝑘4 + 0.1701. 𝑘5
Resolviendo el sistema de ecuaciones mediante el método de la matriz inversa:

0 0 0.6751 −3.2144 −0.4608

0 0 1 1 1
𝐴= 0 1 0.6751 −3.2144 −0.4608
0 −1 1.1309 7.1180 −0.2485
( 1 1 −3.4812 −1.6891 0.1701 )
 0
𝑘1
0
𝑘2 = (𝐴−1 ). 1
𝑘3 0
𝑘4 (0)
(𝑘5 )

0 1.00005 1.00003 1.00001 1.00001 0

𝑘1 0 −1 −5.5 𝑥 10−7 1 −5.5 𝑥 10−7 0
𝑘2 = 0 0.37916 0.33521 0.22633 0.22633 . 1
𝑘3 0 −0.20674 −0.02904 −0.09336 −0.09336 0
𝑘4 (1 −0.17241 0.69378 −0.31970 −0.31970 ) (0)
(𝑘5 )
Por lo tanto, obtenemos los siguientes valores para:
𝑘1 = 1.00001; 𝑘2 = 1; 𝑘3 = 0.2263; 𝑘4 = 0.09336; 𝑘5 = −0.31970
Así, obtenemos las siguientes ecuaciones para cada concentración en un 𝑡 = 𝑡

𝐶𝐴 = 0.1528 𝑒 𝜆3 .𝑡 − 0.3001𝑒 𝜆4 .𝑡 + 0.1473 𝑒 𝜆5 .𝑡

𝐶𝐵 = 0.2263 𝑒 𝜆3 .𝑡 − 0.0933𝑒 𝜆4 .𝑡 + 0.173 𝑒 𝜆5 .𝑡

𝐶𝐶 = 𝑒 𝜆2 .𝑡 + 0.1528 𝑒 𝜆3 .𝑡 − 0.3001𝑒 𝜆4 .𝑡 + 0.0797 𝑒 𝜆5 .𝑡

𝐶𝐷 = − 𝑒 𝜆2 .𝑡 + 0.2559 𝑒 𝜆3 .𝑡 − 0.6645𝑒 𝜆4 .𝑡 + 0.043 𝑒 𝜆5 .𝑡

𝐶𝐸 = 𝑒 𝜆1 .𝑡 − 0.7879 𝑒 𝜆3 .𝑡 − 0.1577𝑒 𝜆4 .𝑡 + 0.0294 𝑒 𝜆5 .𝑡

 PROBLEMA 5

SOLUCIÓN:
Derivación de las ecuaciones de velocidad
(ⅆ𝐶𝐴 )
= 𝐶𝐵 𝐾2 + 𝐶𝐷 𝐾1 0 − 𝐶𝐴 𝐾1
ⅆ𝑡
(ⅆ𝐶𝐵 )
= 𝐶𝐴 𝐾1 + 𝐶𝐷 𝐾6 + 𝐶𝐶 𝐾4 + 𝐶𝐸 𝐾9 − 𝐶𝐵 (𝐾3 + 𝐾2 + 𝐾5 )
ⅆ𝑡
(ⅆ𝐶𝐶 )
= 𝐶𝐵 𝐾3 − 𝐶𝐶 (𝐾1 + 𝐾8 )
ⅆ𝑡
(ⅆ𝐶𝐷 )
= 𝐶𝐵 𝐾5 − 𝐶𝐴 (𝐾10 + 𝐾7 + 𝐾6 )
ⅆ𝑡
(ⅆ𝐶𝐸 )
= 𝐶𝐷 𝐾7 + 𝐶𝐶 𝐾8 − 𝐶𝐸 𝐾9
ⅆ𝑡
Matriz de coeficientes estequiométricos

𝑘1 = 1, 𝑘2 = 2, 𝑘3 = 3, 𝑘4 = 1, 𝑘5 = 2, 𝑘6 = 1, 𝑘7 = 2, 𝑘8 = 1, 𝑘9 = 2, 𝑘1 0 = 1
 Eigenvalores (𝐶 − 𝜆𝐼) = 0 … … … … … . (1)
Reemplazando en 1 se tiene
−𝜆5 − 12𝜆4 − 28𝜆3 + 55𝜆2 + 17𝜆 + 80 = 0
𝐸𝑖𝑔𝑒𝑛𝑣𝑎𝑙𝑜𝑟𝑒𝑠 𝜆1 = −4𝜆2 = −0.6𝜆3 = −1.8𝜆4 = 2.3𝜆5 = −7.9
Eigenvectores (𝐶 − 𝜆𝐼) = 0 … … … … . . (2)
𝑃𝑎𝑟𝑎 𝜆1 = −4
𝑋(𝐶 − 𝜆𝐼) = 0 → 𝑋𝐶 − 𝑋𝜆 = 0 → 𝑋𝐶 = 𝑋𝜆
Reemplazando

Tenemos el sistema de ecuaciones:

3𝑋1 + 6𝑋2 + 𝑋3 = −3𝑋1 → 6𝑋1 + 6𝑋2 + 𝑋3 = 0 → 𝑋3 = 0, 𝑋1 = −1𝑋2 = 1
5𝑋1 + 2𝑋2 + 7𝑋3 = −3𝑋2 → 5𝑋1 + 5𝑋2 + 7𝑋3 = 0 → 𝑋3 = 0, 𝑋1 = −1𝑋2 = 1
2𝑋1 + 2𝑋2 + 𝑋3 = −3𝑋3 → 2𝑋1 + 2𝑋2 + 4𝑋3 = 0 → 𝑋3 = 0, 𝑋1 = −1𝑋2 = 1

(9 + √113)
𝑃𝑎𝑟𝑎 𝜆2 =
2
𝑋(𝐶 − 𝜆𝐼) = 0 → 𝑋𝐶 − 𝑋𝜆 = 0 → 𝑋𝐶 = 𝑋𝜆
Reemplazando
 𝑇𝑒𝑛𝑒𝑚𝑜𝑠 𝑒𝑙 𝑠𝑖𝑠𝑡𝑒𝑚𝑎 ⅆ𝑒 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛𝑒𝑠

(9 + √113)
3𝑋1 + 6𝑋2 + 𝑋3 = 𝑋1 → 𝑋3 = 56, 𝑋1 = 3 + 11√113𝑋2 = 95 + 3√113
2
(9 + √113)
5𝑋1 + 2𝑋2 + 7𝑋3 = 𝑋2 → 𝑋3 = 56, 𝑋1 = 3 + 11√113 𝑋2
2
= 95 + 3√113

(9 + √113)
2𝑋1 + 2𝑋2 + 𝑋3 = 𝑋3 → 𝑋3 = 56, 𝑋1 = 3 + 11√113 𝑋2 = 95 + 3√113
2

(9 − √113)
𝑃𝑎𝑟𝑎 𝜆3 =
2
 PROBLEMA 6

SOLUCIÓN:
ⅆ𝐶𝐴
= −(𝑘5 + 𝑘3 )𝐶𝐴 + 𝑘6 𝐶𝐵 + 𝑘1 𝐶𝐷
ⅆ𝑡
ⅆ𝐶𝐵
= −(𝑘6 + 𝑘4 )𝐶𝐵 + 𝑘5 𝐶𝐴 + 𝑘2 𝐶𝐷
ⅆ𝑡
ⅆ𝐶𝐶
= 𝑘3 𝐶𝐴 + 𝑘4 𝐶𝐵
ⅆ𝑡
ⅆ𝐶𝐷
= −𝑘2 𝐶𝐷 − 𝑘1 𝐶𝐷 = −(𝑘2 + 𝑘1 )𝐶𝐷
ⅆ𝑡
De forma matricial, se tiene:
−(𝑘5 + 𝑘3 ) 𝑘6 0 𝑘1
𝑘5 −(𝑘6 + 𝑘4 ) 0 𝑘2
𝐴=[ ]
𝑘3 𝑘4 0 0
0 0 0 −(𝑘2 + 𝑘1 )
Se asume los valores para:

𝑘1 = 1, 𝑘2 = 0.1, 𝑘3 = 2, 𝑘4 = 0.5, 𝑘5 = 2.5, 𝑘6 = 0.2

Reemplazando los valores se tiene:
−4.5 0.2 0 1
𝐴 = [ 2.5 −0.7 0 0.1 ]
3 0.5 0 0
0 0 0 −1.1
Teorema de valor propio:

𝐴 − 𝜆𝐼 = 0
−4.5 0.2 0 1 𝜆 0 0 0
ⅆ𝑒𝑡 ([ 2.5 −0.7 0 0.1 ] − [0 𝜆 0 0]) = 0
3 0.5 0 0 0 0 𝜆 0
0 0 0 −1.1 0 0 0 𝜆
−4.5 − 𝜆 0.2 0 1
ⅆ𝑒𝑡 | 2.5 −0.7 − 𝜆 0 0.1 | = 0
3 0.5 −𝜆 0
0 0 0 −1.1 − 𝜆
Se obtiene la ecuación:

𝜆4 + 6.3𝜆3 + 8.37𝜆2 + 2.915𝜆 = 0

Los valores obtenidos son:

𝜆1 = 0, 𝜆2 = −1.1, 𝜆3 = −4.63, 𝜆4 = −0.57

De la siguiente expresión:

𝐴𝑋 = 𝜆𝑖 𝑋

- Para 𝜆1 = 0, se tiene:
−4.5 0.2 0 1 𝑥1 𝑥1
𝑥 𝑥
0 0.1 ] [ 2 ] = 0 [ 2 ]
[ 2.5 −0.7
𝑥3 𝑥3
3 0.5 0 0
0 0 0 −1.1 4 𝑥 𝑥4
0.2𝑥2 + 2𝑥4 − 4.5𝑥1 0𝑥1
2.5𝑥1 + 0.1𝑥4 − 0.7𝑥2 0𝑥
[ ] = [ 2]
3𝑥1 + 0.5𝑥2 0𝑥3
−1.1𝑥4 0𝑥4
Entonces:

0.2𝑥2 + 2𝑥4 − 4.5𝑥1 = 0

2.5𝑥1 + 0.1𝑥4 − 0.7𝑥2 = 0
3𝑥1 + 0.5𝑥2 = 0
−1.1𝑥4 = 0
Los valores obtenidos son:

𝑥1 =0
𝑥2 =0
𝑥3 =1
𝑥4 =0
Análogamente se sigue con los siguientes valores de 𝜆𝑖

- Para 𝜆2 = −1.1, se tiene:

𝑥1 = 0.204
𝑥2 = −1.527
𝑥3 = 0.137
𝑥4 = 1
- Para 𝜆3 = −4.63, se tiene:

𝑥1 = −1.726
𝑥2 = 1.098
𝑥3 = 1
𝑥4 = 0

- Para 𝜆4 = −0.57, se tiene:

 𝑥1 = −0.045
𝑥2 = −0.877
𝑥3 = 1
𝑥4 = 0
La matriz modal conteniendo los eigenvectores correspondientes es:
𝑋𝜆1 𝑋𝜆2 𝑋𝜆3 𝑋𝜆4

0 0.204 −1.726 −0.045

𝑀=[0 −1.527 1.098 −0.877]
1 0.137 1 1
0 1 0 0
La solución general de la ecuación planteada para el ejercicio 5 sería:

𝐶𝑖 (𝑡) = ∑ 𝛼𝑖 𝑒 𝜆𝑗𝑡 𝑋𝑖

Se puede representar de la siguiente manera:

𝐶1 (𝑡) 𝑋11 𝑋21 𝑋31 𝑋41
𝐶2 (𝑡) 𝑋 𝑋 𝑋 𝑋
= 𝛼1 𝑒 𝜆1 𝑡 [ 12 ] + 𝛼2 𝑒 𝜆2 𝑡 [ 22 ] + 𝛼3 𝑒 𝜆3 𝑡 [ 32 ] + 𝛼4 𝑒 𝜆4 𝑡 [ 42 ]
𝐶3 (𝑡) 𝑋13 𝑋23 𝑋33 𝑋43
[𝐶4 (𝑡)] 𝑋14 𝑋24 𝑋34 𝑋44

Las constantes 𝛼𝑖 se determinan con las condiciones iniciales:

∀𝒕 = 𝟎; 𝐶1 (0) = 1; 𝐶2 (0) = 𝐶3 (0) = 𝐶4 (0) = 0

Sus valores son:
1 0 0.204 −1.726 −0.045
0 0
[ ] = 𝛼1 [ ] + 𝛼2 [−1.527 ] + 𝛼3 [ 1.098 ] + 𝛼4 [−0.877]
0 1 0.137 1 1
0 0 1 0 0
Finalmente se obtienen los siguientes valores:

𝛼1 = 1.26
𝛼2 = 0
𝛼3 = −0.561
𝛼4 = −0.702
 PROBLEMA 7

SOLUCIÓN:
Derivación de las ecuaciones de velocidad
(ⅆ𝐶𝐴 )
= 𝐶𝐵 𝐾6 + 𝐶𝐷 𝐾1 − 𝐶𝐴 𝐾3 − 𝐶𝐴 𝐾5
ⅆ𝑡
(ⅆ𝐶𝐵 )
= 𝐶𝐴 𝐾5 + 𝐶𝐷 𝐾2 − 𝐶𝐵 (𝐾6 + 𝐾4 )
ⅆ𝑡
(ⅆ𝐶𝐶 )
= 𝐶𝐵 𝐾4 + 𝐶𝐴 𝐾3 − 𝐶𝐶 K 7
ⅆ𝑡
(ⅆ𝐶𝐷 )
= −𝐶𝐷 (𝐾2 + 𝐾1 )
ⅆ𝑡
(ⅆ𝐶𝐸 )
= 𝐶𝐶 𝐾7
ⅆ𝑡
Matriz de coeficientes estequiométricos
−𝐾3 − 𝐾5 𝐾6 0 𝐾1 0
𝐾5 (−𝐾6 − 𝐾4 ) 0 𝐾2 0
𝐾3 𝐾4 − 𝐾7 0 0
0 0 0 − (𝐾2 − 𝐾1 ) 0
[0 0 𝐾7 0 0]

Suponiendo los siguientes valores:

𝑘1 = 1, 𝑘2 = 1, 𝑘3 = 2, 𝑘4 = 1, 𝑘5 = 1, 𝑘6 = 1, 𝑘7 = 1
−3 1 0 1 0 1 0 0 0 0
2 −2 0 1 0 0 1 0 0 0
2 1 −1 0 0 −𝜆 0 0 1 0 0 =0
0 0 0 −2 0 0 0 0 1 0
[0 0 1 0 0 ] [0 0 0 0 1]
−3 − 𝜆 1 0 1 0
2 −2 − 𝜆 0 1 0
1 1 −1 − 𝜆 0 0 =0
0 0 0 −2 − 𝜆 0
[ 0 0 1 0 −𝜆]
Dando como resultado:
−𝜆5 − 8𝜆4 − 21𝜆3 − 22𝜆2 − 8𝜆 = 0
De donde se obtiene las siguientes raíces:
𝐸𝑖𝑔𝑒𝑛𝑣𝑎𝑙𝑜𝑟𝑒𝑠: 𝜆1 = 0, 𝜆2 = −1, 𝜆3 = −1, 𝜆4 = −2 𝑦 𝜆5 = −4
Eigenvectores (𝐶 − 𝜆𝐼) = 0 … … … … . . (2)

𝑃𝑎𝑟𝑎 𝜆2 = −1
𝑋(𝐶 − 𝜆𝐼) = 0 → 𝑋𝐶 − 𝑋𝜆 = 0 → 𝑋𝐶 = 𝑋𝜆

Reemplazando:
−3𝑋1 𝑋2 0 𝑋4 0 −𝑋1
2𝑋1 −2𝑋2 0 𝑋3 0 −𝑋2
2𝑋1 𝑋2 −𝑋3 0 0 = −𝑋3
0 0 0 −2𝑋4 0 −𝑋4
[ 0 0 𝑋3 0 0] [ −𝑋5 ]
Obteniendo un sistema de ecuaciones:
−3𝑋1 + 𝑋2 + 𝑋4 = −𝑋1 → 4𝑋1 + 𝑋2 + 𝑋4 = 0
2𝑋1 − 2𝑋2 + 𝑋3 = −𝑋2 → 2𝑋1 − 𝑋2 + 𝑋3 = 0
2𝑋1 + 𝑋2 − 𝑋3 = −𝑋3 → 2𝑋1 + 2𝑋2 + 4𝑋3 = 0
−2𝑋4 = −𝑋4
𝑋3 = −𝑋5 → 𝑋3 + 𝑋5 = 0

Quedando los siguientes Eigenvectores:

𝑋1 = 0, 𝑋2 = 0, 𝑋3 = 3.925 ∗ 10−17 , 𝑋4 = 6.871 ∗ 10−1 𝑦 𝑋5 = −1.48 ∗ 10−1

Análogamente, para las demás raíces:

𝑋1 = 0, 𝑋2 = 0, 𝑋3 = 7.85 ∗ 10−17 , 𝑋4 = −6.871 ∗ 10−1 𝑦 𝑋5 = −4.46 ∗ 10−1
𝑋1 = 0, 𝑋2 = 7.07 ∗ 10−1 , 𝑋3 = −7.07 ∗ 10−1 , 𝑋4 = −2.29 ∗ 10−1 𝑦 𝑋5 = 7.43 ∗ 10−1

𝑋1 = 0, 𝑋2 = 0, 𝑋3 = 0, 𝑋4 = 0 𝑦 𝑋5 = 2.97 ∗ 10−1

𝑋1 = 1, 𝑋2 = 7.07 ∗ 10−1 , 𝑋3 = 7.07 ∗ 10−1 , 𝑋4 = 5.726 ∗ 10−2 𝑦 𝑋5 = −3,71 ∗ 10−1

 PROBLEMA 8

SOLUCIÓN:
Ecuación de la velocidad:
ⅆ𝐶𝐴
= −(𝑘3′ + 𝑘1 + 𝑘4 )𝐶𝐴 + 𝐶𝐵 𝑘4′ + 𝐶𝑆 𝑘3 + 𝐶𝑅 𝑘1′
ⅆ𝑡
ⅆ𝐶𝐵
= −(𝑘4′ + 𝑘5 )𝐶𝐵 + 𝐶𝐴 𝑘4 + 𝐶𝑅 𝑘5 ′
ⅆ𝑡
ⅆ𝐶𝑆
= −(𝑘3 + 𝑘2′ )𝐶𝑆 + 𝐶𝐴 𝑘3′ + 𝐶𝑅 𝑘2
ⅆ𝑡
ⅆ𝐶𝑅
= −(𝑘1′ + 𝑘2 + 𝑘5′ )𝐶𝑅 + 𝐶𝐴 𝑘1 + 𝐶𝑆 𝑘2′ + 𝐶𝐵 𝑘5
ⅆ𝑡
Asumiendo valores para 𝑘𝑖𝑗 :

𝑘1 = 1 𝑘1′ = 2
𝑘4 = 3 𝑘4′ = 1
𝑘2 = 2 𝑘2′ = 1
𝑘5 = 2 𝑘5′ = 3
𝑘3 = 3 𝑘3′ = 2
Matriz de coeficientes estequiométricos:
−(𝑘3′ + 𝑘1 + 𝑘4 ) 𝑘4′ 𝑘3 𝑘1′
𝑘4 −(𝑘5 + 𝑘4′ ) 0 𝑘5′
𝑘𝑖𝑗 = ′
−(𝑘3 + 𝑘2 )
𝑘3′ 0 𝑘2
′
[ 𝑘1 𝑘 5 𝑘 2 −(𝑘1 + 𝑘2 + 𝑘5′ ) ]
′

Con los valores de k iniciales:

−6 1 3 2 𝐶𝐴
ⅆ𝐶̅ 𝐶
= [ 3 −3 0 3 ] [ 𝐵 ]
ⅆ𝑡 2 0 −4 2 𝐶𝑆
1 2 1 −6 𝐶𝑅
Resolviendo:
|[𝐴 − 𝜆]𝐽| = 0

−6 1 3 2 𝜆1 0 0 0 1 0 0 0
3 −3 0 3 𝜆2 0 0 0 1 0 0 )| = 0
| ([ ]−[0 0 ]) (0
2 0 −4 2 0 0 𝜆3 0 1 0
1 2 1 −6 0 0 0 𝜆4 0 0 0 1

−6 − 𝜆1 1 3 2
−3 − 𝜆 0 3
| 3 2
2 |=0
2 0 −4 − 𝜆 3
1 2 1 −6 − 𝜆4

Suponiendo que 𝜆1 = 𝜆2 = 𝜆3 = 𝜆4 = 𝜆 para obtener un polinomio

𝜆4 + 19𝜆3 + 113𝜆2 + 195𝜆 − 42 = 0
De las raíces del polinomio
𝜆1 = 0.19
𝜆2 = −3.57
𝜆3 = −7.51
𝜆4 = −8.11
De la siguiente expresión
𝐴𝑋 = 𝜆𝑖 𝑋
• Para 𝜆1 = 0.19
−6 1 3 2 𝑥1 𝑥1
𝑥 𝑥
[ 3 −3 0 3 ] [𝑥2 ] = 0.19 [𝑥2 ]
2 0 −4 2 3 3
1 2 1 −6 𝑥4 𝑥4

−6𝑥1 + 𝑥2 + 3𝑥3 + 2𝑥4 0.19𝑥1 −6𝑥1 + 𝑥2 + 3𝑥3 + 2𝑥4 = 0.19𝑥1

3𝑥1 − 3𝑥2 + 3𝑥4 0.19𝑥2 3𝑥1 − 3𝑥2 + 3𝑥4 = 0.19𝑥2
[ ]= [ ]→
2𝑥1 − 4𝑥3 − 2𝑥4 0.19𝑥3 2𝑥1 − 4𝑥3 − 2𝑥4 = 0.19𝑥3
𝑥1 + 2𝑥2 + 𝑥3 − 6𝑥4 0.19𝑥4 𝑥1 + 2𝑥2 + 𝑥3 − 6𝑥4 = 0.19𝑥4

La solución general:

0
𝑥 = [0 ]
0
0
De igual forma se calcular para los demás valores
• Para 𝜆2 = −3.57, se sigue el mismo procedimiento

𝑥1 0
𝑥2
𝑥 = [𝑥 ] = [0]
3 0
𝑥4 0
• Para 𝜆3 = −7.51

𝑥1 0
𝑥2
𝑥 = [𝑥 ] = [0]
3 0
𝑥4 0

• Para 𝜆4 = −8.11
𝑥1 0
𝑥2
𝑥 = [𝑥 ] = [0]
3 0
𝑥4 0
 PROBLEMA 9 – REACCIÓN BITRIANGULAR

SOLUCIÓN:
Asumiendo valores para 𝑘𝑖𝑗

𝑘1 = 1 𝑘1′ = 2 𝑘6 = 3 𝑘6′ = 1
𝑘4 = 1
𝑘2 = 2 𝑘2′ = 1 𝑘7 = 2 𝑘7′ = 3
𝑘5 = 2
𝑘3 = 3 𝑘3′ = 2 𝑘8 = 2 𝑘8′ = 4
Entonces
ⅆ𝐶𝐴
= −𝑘3′ 𝐶𝐴 − 𝑘1 𝐶𝐴 + 𝑘3 𝐶𝑆 + 𝑘1′ 𝐶𝑅 − 𝑘4 𝐶𝐴
ⅆ𝑡
ⅆ𝐶𝑆
= −𝑘3 𝐶𝑆 − 𝑘2′ 𝐶𝑠 + 𝑘3′ 𝐶𝐴 + 𝑘2 𝐶𝑅
ⅆ𝑡
ⅆ𝐶𝑅
= −𝑘2 𝐶𝑅 − 𝑘1′ 𝐶𝑅 + 𝑘2′ 𝐶𝑆 + 𝑘1 𝐶𝐴 + 𝑘5 𝐶𝑉
ⅆ𝑡
ⅆ𝐶𝑇
= −𝑘8′ 𝐶𝑇 + 𝑘4 𝐶𝐴 + 𝑘8 𝐶𝑉 − 𝑘6 𝐶𝑇 + 𝑘6′ 𝐶𝑈
ⅆ𝑡
𝑑𝐶𝑉
= −𝑘8 𝐶𝑉 − 𝑘5 𝐶𝑉 − 𝑘7′ 𝐶𝑉 + 𝑘8′ 𝐶𝑇 + 𝑘7 𝐶𝑈
𝑑𝑡

ⅆ𝐶𝑈
= −𝑘6′ 𝐶𝑈 − 𝑘7 𝐶𝑈 + 𝑘6 𝐶𝑇 + 𝑘7′ 𝐶𝑉
ⅆ𝑡
Reorganizando
ⅆ𝐶𝐴
= −(𝑘3′ + 𝑘1 + 𝑘4 )𝐶𝐴 + 𝑘3 𝐶𝑆 + 𝑘1′ 𝐶𝑅
ⅆ𝑡
ⅆ𝐶𝑆
= 𝑘3′ 𝐶𝐴 − (𝑘3 + 𝑘2′ )𝐶𝑆 + 𝑘2 𝐶𝑅
ⅆ𝑡
ⅆ𝐶𝑅
= 𝑘1 𝐶𝐴 + 𝑘2′ 𝐶𝑆 − (𝑘2 + 𝑘1′ )𝐶𝑅 + 𝑘5 𝐶𝑉
ⅆ𝑡
ⅆ𝐶𝑇
= 𝑘4 𝐶𝐴 − (𝑘8′ + 𝑘6 )𝐶𝑇 + 𝑘8 𝐶𝑉 + 𝑘6′ 𝐶𝑈
ⅆ𝑡
ⅆ𝐶𝑉
= −(𝑘8 + 𝑘5 + 𝑘7′ )𝐶𝑉 + 𝑘8′ 𝐶𝑇 + 𝑘7 𝐶𝑈
ⅆ𝑡
 ⅆ𝐶𝑈
= 𝑘6 𝐶𝑇 + 𝑘7′ 𝐶𝑉 − (𝑘6′ − 𝑘7 )𝐶𝑈
ⅆ𝑡
De forma matricial

ⅆ𝐶̅ = [𝑘𝑖𝑗 ]𝐶̅

𝐶𝐴
𝐶𝑆
𝐶
𝐶̅ = 𝑅
𝐶𝑇
𝐶𝑉
[𝐶𝑈 ]
𝑘𝑖𝑗
−(𝑘3′ + 𝑘1 + 𝑘4 ) 𝑘3 𝑘1′ 0 0 0
′
𝑘3 ′
−(𝑘3 + 𝑘2 ) 𝑘2 0 0 0
𝑘1 ′ ′
−(𝑘2 + 𝑘1 ) 0 𝑘5 0
= 𝑘2 ′
𝑘4 0 −(𝑘8 + 𝑘6 ) 𝑘8 𝑘6′
0 ′ ′
0 0 0 𝑘8 −(𝑘8 + 𝑘5 + 𝑘7 ) 𝑘7
′ ′
[ 0 0 0 𝑘6 𝑘7 −(𝑘96 + 𝑘7 )]
Con los valores de k iniciales

−4 3 2 0 0 0 𝐶𝐴
2 −4 2 0 0 0 𝐶𝑆
ⅆ𝐶̅ 1 2 −4 0 2 0 𝐶𝑅
=
ⅆ𝑡 1 0 0 −7 2 1 𝐶𝑇
0 0 0 4 −7 2 𝐶𝑉
[ 0 0 0 3 3 −3] [𝐶 ]
𝑈

Resolvemos la siguiente determinante

|[𝐴 − 𝜆]𝐽| = 0
−4 3 2 0 0 0 𝜆1 0 0 0 0 0 1 0 0 0 0 0
2 −4 2 0 0 0 0 𝜆2 0 0 0 0 0 1 0 0 0 0
| 1 2 −4 0 2 0 0 0 𝜆3 0 0 0 0 0 1 0 0 0 |
− =0
| 1 0 0 −7 2 1 0 0 0 𝜆4 0 0 0 0 0 1 0 0 |
0 0 0 4 −7 2 0 0 0 0 𝜆5 0 0 0 0 0 1 0
([ 0 0 0 3 3 −3] [ 0 0 0 0 0 𝜆6 ]) (0 0 0 0 0 1)
−4 − 𝜆1 3 2 0 0 0
2 −4 − 𝜆2 2 0 0 0
| 1 2 −4 − 𝜆3 0 2 0 |
=0
| 1 0 0 −7 − 𝜆4 2 1 |
0 0 0 4 −7 − 𝜆5 2
0 0 0 3 3 −3 − 𝜆6
Tomando como 𝜆1 = 𝜆2 = 𝜆3 = 𝜆4 = 𝜆5 = 𝜆6 = 𝜆 para obtener un polinomio
𝜆6 + 29𝜆5 + 314𝜆4 + 1538𝜆3 + 3114𝜆2 + 1260𝜆 − 432 = 0
De las raíces del polinomio
 𝜆1 = −6
𝜆2 = 0.215
𝜆3 = −0.883
𝜆4 = −5.628
𝜆5 = −6.774
𝜆6 = −9.930
De la siguiente expresión
𝐴𝑋 = 𝜆𝑖 𝑋
• Para 𝜆1 = −6

−4 3 2 0 0 0 𝑥1 𝑥1
2 −4 2 0 0 0 𝑥2 𝑥2
1 2 −4 0 2 0 𝑥3 𝑥3
𝑥 = −6 𝑥
1 0 0 −7 2 1 4 4
0 0 0 4 −7 2 𝑥5 𝑥5
[ 0 0 0 3 3 −3] [𝑥6 ] [𝑥6 ]
−4𝑥1 + 3𝑥2 + 2𝑥3 −6𝑥1 −4𝑥1 + 3𝑥2 + 2𝑥3 = −6𝑥1
2𝑥1 − 4𝑥2 + 2𝑥3 −6𝑥2 2𝑥1 − 4𝑥2 + 2𝑥3 = −6𝑥2
𝑥1 + 2𝑥2 − 4𝑥3 + 2𝑥5 −6𝑥3 𝑥1 + 2𝑥2 − 4𝑥3 + 2𝑥5 = −6𝑥3
= →
𝑥1 − 7𝑥4 + 𝑥6 + 2𝑥5 −6𝑥4 𝑥1 − 7𝑥4 + 𝑥6 + 2𝑥5 = −6𝑥4
4𝑥4 + 2𝑥6 − 7𝑥5 −6𝑥5 4𝑥4 + 2𝑥6 − 7𝑥5 = −6𝑥5
[ 3𝑥4 − 3𝑥6 + 3𝑥5 ] [−6𝑥6 ] 3𝑥4 − 3𝑥6 + 3𝑥5 = −6𝑥6
Resolviendo el sistema de ecuaciones
Para este caso el sistema de ecuaciones presenta infinitas soluciones
• Para 𝜆2 = 0.215, se sigue el mismo procedimiento
𝑥1 = 0
𝑥2 = 0
𝑥3 = 0
𝑥4 = 0
𝑥5 = 0
𝑥6 = 0
• Para 𝜆3 = −0.883
𝑥1 = 0
𝑥2 = 0
𝑥3 = 0
𝑥4 = 0
𝑥5 = 0
𝑥6 = 0
• Para 𝜆4 = −5.628
 𝑥1 =0
𝑥2 =0
𝑥3 =0
𝑥4 =0
𝑥5 =0
𝑥6 =0
• Para 𝜆5 = −6.774
𝑥1 =0
𝑥2 =0
𝑥3 =0
𝑥4 =0
𝑥5 =0
𝑥6 =0

• Para 𝜆6 = −9.930

𝑥1 =0
𝑥2 =0
𝑥3 =0
𝑥4 =0
𝑥5 =0
𝑥6 =0