Tarea 1 Ing de Reacciones, Angela Teran
Tarea 1 Ing de Reacciones, Angela Teran
Tarea 1 Ing de Reacciones, Angela Teran
FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA QUÍMICA
INGENIERÍA DE LAS REACCIONES QUÍMICAS
740045M
HOMEWORK 1
Para la solución de este problema usaremos el método del semiperiodo de reacción para una
reacción de orden n :
dCA n
=−k A C A
dt
Cuya integración da:
C1−n 1−n
A −C A 0
=−k A t , para n ≠ 1
1−n
Que puede reacomodarse como:
( )
1−n
CA n−1
=1+C A 0 ( n−1 ) k A t
CA 0
1
Sustituyendo C A= C A 0 y t=t 1 / 2:
2
n−1
2 −1
t 1 / 2=
( n−1 ) C n−1
A0 k A
Por lo que una representación de ln t 1 / 2 frente a ln C A 0 da una línea recta de pendiente 1−n.
Con los datos proporcionados por el problema, se calculan los tiempos de vida media para puntos
seleccionados de concentración inicial, los resultados se muestran en la tabla a continuación,
donde los puntos se calculan por medio de interpolación lineal.
t(h) CA(mol/L) CA0 C(t_1/2) t_0 t_f t_1/2 CA0 ln(t_1/2) LN(CA0]
0 2 2 1 0 3.70588235 3.70588235 2 1.30992138 0.69314718
3 1.08 1.75 0.875 0.81521739 4.80882353 3.99360614 1.75 1.38469462 0.55961579
6 0.74 1.5 0.75 1.63043478 5.91176471 4.28132992 1.5 1.45426369 0.40546511
9 0.56 1.25 0.625 2.44565217 7.91666667 5.47101449 1.25 1.69946406 0.22314355
12 0.46 1 0.5 3.70588235 10.8 7.09411765 1 1.95926594 0
ln t_1/2 vs ln C_A0
2
1.9
f(x) = − 0.953832724674465 x + 1.92042470408957
1.8 R² = 0.963114267207734
1.7
Ln t_1/2
1.6
1.5
1.4
1.3
1.2
-0.3 0.2 0.7
ln C_A0
−( n−1 )=−0.9538
Entonces
n=1.9538 ≈ 2
2−1
2 −1
ln =1.9204
( 2−1 ) k A
Entonces:
k A =0.1465 L / mol∗h
Con base en este resultado y con el fin de evitar errores debidos a la interpolación lineal, se
procede a ajustar los datos a un modelo de segundo orden. Partiendo de la ecuación para una
cinética de orden n y sustituyendo para n=2:
( )
−1
CA
=1+C 1A 0 ( 1 ) k A t
CA 0
De donde:
CA 0
=1+ k A C A 0 t
CA
2.5
2
1.5
1
0.5
0
0 2 4 6 8 10 12 14
t(h)
C A 0 k A =0.2805
Entonces:
Análisis estadístico
Estadísticas de la regresión
0.9998350
Coeficiente de correlación múltiple 9
0.9996702
Coeficiente de determinación R^2 1
0.9995602
R^2 ajustado 8
0.0279059
Error típico 7
Observaciones 5
El R2 ajustado indica que la regresión de la ecuación linealizada indica que el modelo es un muy
buen ajuste de los datos cinéticos.
ANÁLISIS DE VARIANZA
Grados de libertad
Suma de cuadrados
Promedio de los cuadrados F Valor crítico de F
Regresión 1 7.08160773 7.08160773 9093.63779 2.5421E-06
Residuos 3 0.00233623 0.00077874
Total 4 7.08394396
Respecto al análisis de varianza podemos notar que el valor del estadístico F es bastante mayor
que el valor de F crítico, bajo la prueba de hipótesis:
H 1 : en caso contrario
Por lo que es una conclusión fuerte el hecho de que los coeficientes son distintos de cero.
Nótese que se tiene un error en este análisis debido a que de antemano sabemos que el
intercepto es 1, así que lo más conveniente en este caso, es realizar un ajuste C A 0 / C A−1 vs t y
hacer una regresión para el intercepto igual a cero:
ANÁLISIS DE VARIANZA
Grados de Suma de Promedio de los Valor crítico
libertad cuadrados cuadrados F de F
Regresió 33438.464
n 1 21.4424675 21.4424675 4 3.60624E-07
Residuos 4 0.00256501 0.00064125
Total 5 21.4450325
Nótese que ahora podemos concluir fuertemente que la regresión ajusta de manera muy
adecuada los datos cinéticos. Respecto al intervalo de confianza para la constante cinética k A ,
tenemos:
Coeficientes Error típico Estadístico t Probabilidad Inferior 95% Superior 95% Inferior 95.0%Superior 95.0%
Intercepción 0 #N/D #N/D #N/D #N/D #N/D #N/D #N/D
Variable X 1 0.28180941 0.00154111 182.861873 5.365E-09 0.27753062 0.28608821 0.27753062 0.28608821
0.2775 / 2≤ k A ≤ 0.2861 / 2
0.1388 ≤ k A ≤ 0.1431
Es claro que este análisis se pudo hacer también para el orden de reacción n tomando como base
el ajuste lineal a la cinética de orden general, aunque tomando las recomendaciones de la
literatura, se aproxima al valor más cercano (si el error relativo no es muy alto), en este caso n=2.