UNIVERSIDAD AUTÓNOMA DE NUEVO LEÓN

FACULTAD DE CIENCIAS QUÍMICAS

INGENIERÍA QUÍMICA

Laboratorio de Reactores Químicos

Reporte de Práctica

Saponificación

Equipo 1 Grupo 2

Samantha Virginia Cañamar Garza 1633367

Julián Eliud Martínez Martínez 1624716

Mauricio Mellado Ayala 1720536

David Alberto Palomino Santoyo 1720570

Dr. Gerardo Escamilla Flores

San Nicolás de Los Garza, N.L. 7 de mayo de 2019

INTRODUCCIÓN
Una ecuación cinética caracteriza la tasa de reacción, la determinación de esta se
realiza de manera experimental que consiste en encontrar la dependencia de la tasa
de reacción a la concentración. Los equipos para obtener los datos se pueden
clasificar en reactores intermitentes (o por lotes) y en reactores continuos. El reactor
intermitente es un recipiente donde se mantienen las sustancias mientras
reaccionan, se mide la extensión de la reacción en diversos tiempos, lo que se logra
de distintas maneras:

1. Siguiendo la concentración de un determinado componente.

2. Siguiendo la variación de alguna propiedad física del fluido (por ejemplo, la
conductividad eléctrica o el índice de refracción).
3. Siguiendo la variación de la presión total de un sistema a volumen constante.
4. Siguiendo la variación del volumen de un sistema a presión constante. [1]

La reacción que será sometida al estudio cinético es saponificación, consiste en

mezclar acetato de etilo con una base fuerte que es el hidróxido de sodio que se
convierten en acetato de sodio y alcohol etílico:

𝐶𝐻3 𝐶𝑂𝑂 − 𝐶𝐻2 𝐶𝐻3 + 𝑁𝑎𝑂𝐻 ↔ 𝐶𝐻3 𝐶𝑂𝑂 − 𝑁𝑎 + 𝐶𝐻3 𝐶𝐻2 − 𝑂𝐻

Debido a que la reacción es ligeramente exotérmica en los rangos de 20°C a 40

°C y se estará trabajando en estas condiciones que favorecen el equilibrio hacia
los reactivos [4]] y para su estudio cinético, la reacción se considerará de la forma:

𝐴 + 𝐵 → 𝑃𝑟𝑜𝑑𝑢𝑐𝑡𝑜𝑠

−𝑟𝐴 = 𝑘𝐶𝐴 𝛼 𝐶𝐵 𝛽

Donde α y β son los órdenes de reacción.

Debido a la literatura, se conoce que al tener acetato de etilo con una base fuerte
como el hidróxido de sodio, sin algún uso de un catalizador y trabajando en rangos
de temperatura de 20°C a 50°C se consigue consumir todo el reactivo para producir
el producto de acetato de sodio. Esta reacción se ha reportado como una reacción
de 2do orden global y con 1er orden para cada reactivo. [4] [5]

La cinética de la reacción puede seguirse por conductimetría debido a las diferentes

especies iónicas en la mezcla reaccionante; la conductividad del sistema se debe
principalmente a los iones OH, cuya conductividad es muy superior a la de los iones
acetato formados como producto. Por tanto, la conductividad Ω de la mezcla es una
medida fiable de la concentración, e irá disminuyendo con el tiempo debido a la
desaparición del ión conductor OH (es decir, al consumo del reactivo B). [2]

2
El método de la conductividad eléctrica es útil para el estudio de reacciones en que
intervienen iones que presentan conductividades iónicas relativamente altas,
particularmente los iones H+ y OH. En disolución diluida, la sustitución de un ión por
otro de diferente conductividad iónica será proporcional a la velocidad de variación
de la concentración del ión reactivo. [3]

Con esta información, se propone la Ec.1 para calcular la concentración en términos

de la conductividad iónica:

𝐶𝐴 (𝑡) 𝜅(𝑡)−𝜅∞
= [1]
𝐶𝐴0 𝜅0 −𝜅∞

Dónde:

𝐶𝐴 (𝑡) Concentración del compuesto A en el tiempo t.

𝐶𝐴0 Concentración inicial del compuesto A.

𝜅(𝑡) Conductividad iónica al tiempo t.

𝜅0 Conductividad iónica inicial.

𝜅∞ Conductividad iónica residual.

OBJETIVO
Determinar el orden de la reacción y la energía de activación de saponificación
mediante un ajuste de datos experimentales para obtener la expresión cinética
más adecuada.

METODOLOGÍA
A continuación para la metodología experimental se hizo un diagrama de flujo que
aplica en las tres partes de esta práctica (A, B, C) solo se cambia un paso
dependiendo que sesión de la práctica ser este realizando y se anexa la tabla de
las condiciones iniciales de cada sesión de la práctica.

3
 Se preparó una Se instaló el programa
Se preparó una "HORIBA LAQUA" en
solución de NaOH con solución de acetato de
etilo con la pc, para tomar los
la concentración que datos de conductividad
decidió el equipo. concentración que
propuso el equipo. con el tiempo.

Al momento de vertir la
solución de acetato de Se midió la Se colocó el NaOH
etilo, correr el conductividad ionica sobre una plancha con
programa para tomar del NaOH inicial, un agitador magnético,
los datos de mediante el equipo para mejorar el
conductividad. mencionado mezclado.
anteriormente.

Para parte B: Para parte C: Se repiten

Esperamos hasta que Se repiten los pasos los pasos anteriores pero
los valores de anteriores con las soluciones se
conductividad en el diferentes calientan o enfrian antes
monitor no cambien. concentraciones. de mezclarlas.

Tabla 1. Condiciones iniciales para sesiones A, B y C

Concentración Concentración Temperatura
Sesión Relación
Acetato [M] NaOH [M] [°C]
Exceso
.1 .01 21
A Acetato
Equimolar .1 .1 21
Exceso
0.15 0.1 19.3
Acetato
Exceso
B 0.1 0.005 19.3
Acetato
Exceso
NaOH
15
19
C Equimolar 0.1 0.1 30
35
40

4
Fig. 1.- Sistema simulando un reactor por lotes.

Fig. 2.- Software para captura de datos.

5
METODOLOGÍA DE ANÁLISIS
Sesión A

i) Alimentación de acetato de Etilo en exceso

Primeramente, utilizando la ecuación 1 y las conductividades a través del tiempo

medidas con el potenciómetro, se determinó la concentración de hidróxido de
sodio a través del tiempo.

Tabla 2 – Resultados para corrida con Ac. De Etilo en exceso

Cond [mS/cm] Tiempo [s] Ca [mol/L]

80.30 0 0.02

49.40 3 0.010689

43.50 6 0.008911

40.90 9 0.008127

37.80 12 0.007193

35.90 15 0.00662

33.90 18 0.006018

32.3 21 0.005536

31.2 24 0.005204

29.5 27 0.004692

28.5 30 0.004391

… … …

13.94 471 3.01E-06

13.94 474 3.01E-06

13.95 477 6.03E-06

13.94 480 3.01E-06

13.94 483 3.01E-06

13.94 486 3.01E-06

13.94 489 3.01E-06

13.94 492 3.01E-06

13.94 495 3.01E-06

13.94 498 3.01E-06

6
Debido a que se llevó el experimento con una alimentación en exceso de acetato
de etilo, se escribió la tasa de reacción para el reactivo limitante de la siguiente
manera
𝑑𝐶𝐴 𝛽
−𝑟𝐴 = − = 𝑘𝐶𝐴𝛼 𝐶𝐵 [2]
𝑑𝑡

Considerando que la concentración de acetato de etilo no cambia significativamente

a través del tiempo se puede considerar que

𝐾 = 𝑘𝐶𝐵 𝛽 [3]

Y la ecuación 2 se puede simplificar en

𝑑𝐶𝐴
−𝑟𝐴 = − = 𝐾𝐶𝐴𝛼 [4]
𝑑𝑡

A) Método Integral

Reescribiendo la ecuación 4, agrupando variables obtenemos

𝐶𝐴 𝑡
𝑑𝐶𝐴
−∫ 𝛼 = ∫ 𝐾𝑑𝑡 [5]
𝐶𝐴0 𝐶𝐴 0

Resolviendo la ecuación para 𝛼 = 1, 𝛼 = 2 𝑦 𝛼 = 𝑛 se obtienen las ecuaciones:

𝐶𝐴
ln ( ) = −𝐾𝑡 𝛼 = 1 [6]
𝐶𝐴0

1 1
− = −𝐾𝑡 𝛼 = 2 [7]
𝐶𝐴0 𝐶𝐴

𝐶𝐴1−𝑛 − 𝐶𝐴 01−𝑛 = (𝑛 − 1)𝐾𝑡 𝛼 = 𝑛 [8]

Las ecuaciones 6,7 y 8 tienen la forma y = mx. Se realizó una regresión por mínimos
cuadrados para encontrar el coeficiente cinético y posteriormente calcular un error
residual.
2
𝑆𝑟 = ∑(𝐶𝐴𝐸𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙 − 𝐶𝐴𝑇𝑒ó𝑟𝑖𝑐𝑜 ) [9]

B) Método diferencial

Para este método, se procedió a linealizar la ecuación 4 de la siguiente manera:

𝑑𝐶𝐴
𝐼𝑛(− ) = ln(𝐾) + 𝛼 ln(𝐶𝐴 ) [10]
𝑑𝑡

7
Donde la ecuación 10 tiene forma de 𝑦 = 𝑏 + 𝑚𝑥 y de acuerdo con una regresión
por mínimos cuadrados:

𝑛 ∑(𝑥𝑦) − (∑ 𝑥)(∑ 𝑦)
𝑚= [11]
𝑛 ∑(𝑥 2 ) − (∑(𝑥))2

∑(𝑦) ∑(𝑥 2 ) − (∑ 𝑥)(∑ 𝑥𝑦)

𝑏= [12]
𝑛 ∑(𝑥 2 ) − (∑(𝑥))2

Donde los valores de las derivadas se encontraron mediante diferencias finitas a

intervalos de tiempo igualmente espaciados usando las siguientes ecuaciones:

Punto inicial:

𝑑𝐶𝐴 −3𝐶𝐴𝑖 + 4𝐶𝐴𝑖+1 − 𝐶𝐴𝑖+2

= [13]
𝑑𝑡 2∆𝑡
Puntos intermedios:

𝑑𝐶𝐴 𝐶𝐴𝑖+1 − 𝐶𝐴𝑖−1

= [14]
𝑑𝑡 2∆𝑡
Punto final:

𝑑𝐶𝐴 𝐶𝐴𝑖−2 − 4𝐶𝐴𝑖−1 + 3𝐶𝐴𝑖

= [15]
𝑑𝑡 2∆𝑡

ii) Alimentación Equimolar

Primeramente, utilizando la ecuación 1 y las conductividades a través del tiempo

medidas con el potenciómetro, se determinó la concentración de hidróxido de sodio
a través del tiempo.

Tabla 3 – Resultados para corrida en condiciones equimolar

Tiempo
Cond [mS/cm] Ca [mol/L]
[s]
374.00 0 0.1000
233.00 3 0.0538
214.00 6 0.0476
213.00 9 0.0473
211.00 12 0.0466
204.00 15 0.0443

8
 196.00 18 0.0417
209.00 21 0.0460
204.00 24 0.0443
195.00 27 0.0414
190.00 30 0.0398
… … …
68.8 2040 6.55E-05
68.8 2043 6.55E-05
68.7 2046 3.27E-05
68.7 2049 3.27E-05
68.7 2052 3.27E-05
68.7 2055 3.27E-05
68.7 2058 3.27E-05
68.7 2061 3.27E-05
68.7 2064 3.27E-05
68.7 2067 3.27E-05
68.7 2070 3.27E-05

Debido a que se llevó el experimento con una alimentación equimolar de acetato de

etilo e hidróxido de sodio, se escribió la tasa de reacción para el hidróxido de sodio
de la siguiente manera
𝑑𝐶𝐴 𝛽
−𝑟𝐴 = − = 𝑘𝐶𝐴𝛼 𝐶𝐵 [16]
𝑑𝑡

Procediendo a determinar una ecuación que relacionara la concentración de b con

a

𝐶𝐵 = 𝐶𝐵0 − (𝐶𝐴 0 − 𝐶𝐴 ) [17]

Y debido a que hay una alimentación equimolar

𝐶𝐴 0 = 𝐶𝐵 0 [18]

Entonces
𝛽
𝐶𝐵 𝛽 = 𝐶𝐴 [19]

Combinando la ecuación 16 y 19 obtenemos

𝑑𝐶𝐴 𝛼+𝛽 𝛾
− = 𝑘 𝐶𝐴 = 𝑘 𝐶𝐴 [20]
𝑑𝑡
Donde 𝛾 es el orden global de la reacción.

9
 A) Método Integral

Reescribiendo la ecuación 20, agrupando variables obtenemos

𝐶𝐴 𝑡
𝑑𝐶𝐴
−∫ 𝛾 = ∫ 𝑘𝑑𝑡 [21]
𝐶𝐴0 𝐶𝐴 0

Resolviendo la ecuación para 𝛾 = 1, 𝛾 = 2 𝑦 𝛾 = 𝑛 se obtienen las ecuaciones:

𝐶𝐴
ln ( ) = −𝑘𝑡 𝛾 = 1 [22]
𝐶𝐴0

1 1
− = −𝑘𝑡 𝛾 = 2 [23]
𝐶𝐴0 𝐶𝐴

𝐶𝐴1−𝑛 − 𝐶𝐴 01−𝑛 = (𝑛 − 1)𝑘𝑡 𝛾 = 𝑛 [24]

Las ecuaciones 22, 23 y 24 tienen la forma y = mx. Se realizó una regresión por
mínimos cuadrados para encontrar el coeficiente cinético y posteriormente calcular
un error residual.

Sesión B

0.1M Acetato de Etilo y 0.05M de NaOH

En primer lugar, usando los valores de las conductividades a través del tiempo medidas con
el potenciómetro, se determinó la concentración de hidróxido de sodio a través del tiempo
hasta llegar al estado estable. Obteniendo la siguiente tabla de datos.

Tabla 4 - Concentraciones experimentales para corrida de 0.1M de Acetato de

Etilo y 0.05M NaOH

Tiempo Conductividad Ca
(s) (S/m) (mol/L)
0 0.847 0.0500
51 0.755 0.0409
102 0.669 0.0324
153 0.611 0.0266
204 0.564 0.0220
255 0.527 0.0183
306 0.498 0.0154
357 0.476 0.0133
408 0.456 0.0113
459 0.441 0.0098

10
 510 0.427 0.0084
549 0.418 0.0075

Método Diferencial

En el caso B no se puede considerar un reactivo en exceso, por lo que se tiene que

considerar que ambos reactivos tienen efecto en la tasa de reacción. Se escribe la tasa de
reacción y se le aplican logaritmos naturales para llegar a una ecuación lineal de 3
parámetros.

−𝑟𝐴 = 𝑘𝐶𝐴 𝛼 𝐶𝐵 𝛽
𝑑𝐶𝐴
ln −𝑟𝐴 = ln (−
𝑑𝑡
) = ln 𝑘 +α ln 𝐶𝐴 + 𝛽 ln 𝐶𝐵 [25]

Calculando la derivada numérica mediante diferencias finitas, y con las concentraciones

experimentales se realizó una regresión lineal múltiple para el cálculo de los órdenes de
reacción y el coeficiente cinético de reacción.

Con los valores de los órdenes de reacción y el coeficiente cinético se puede obtener una
concentración de hidróxido de sodio y calcular un error relativo total.
𝐶𝐴 𝑒𝑥𝑝 −𝐶𝐴 𝑐𝑎𝑙𝑐
𝑒𝑟𝑟𝑜𝑟 𝑡𝑜𝑡𝑎𝑙 = ∑ ( 𝐶𝐴 𝑒𝑥𝑝
) ∗ 100 [26]

0.05M Acetato de Etilo y 0.10M de NaOH

Para esta corrida se usó el pH

0.15M Acetato de Etilo y 0.10M de NaOH

Se repite el mismo procedimiento que en la primera corrida de esta sesión, a continuación

se muestra la tabla de datos obtenidos.

Tabla 5 - Concentraciones experimentales para corrida de 0.15 M de Acetato de

Etilo y 0.10 M NaOH

Tiempo Conductividad Ca
(s) (S/m) (mol/L)
0 3.65 0.1000
51 1.307 0.0231
102 0.976 0.0123
153 0.816 0.0070
204 0.728 0.0041

11
 255 0.677 0.0025
306 0.648 0.0015
354 0.632 0.0010

Sesión C

Obteniendo la concentración inicial (0.05 M), y el cambio de la conductividad en

función del tiempo, se puede calcular la concentración de NaOH.

Sabiendo que el orden de reacción global “n” es 1.63 y la reacción es equimolar,

𝜅(𝑡)−𝜅∞
𝐶𝐴 (𝑡) = 𝐶𝐴0
𝜅0 − 𝜅∞
puede graficar la siguiente ecuación para obtener el parámetro de coeficiente
cinético a partir del coeficiente m una vez realizada la regresión lineal.

𝐶𝐴1−𝑛 − 𝐶𝐴01−𝑛 = (𝑛 − 1)𝑘𝑡 [27]

y = mx
Donde

𝑌 = 𝐶𝐴1−𝑛 − 𝐶𝐴01−𝑛 [28]

𝑚 = (𝑛 − 1)𝑘

𝑥=𝑡
Posteriormente, ya conocido el valor del coeficiente cinético a varias temperaturas,
se puede obtener la energía de activación con la siguiente función:
𝐸𝑎 1
ln 𝑘 = ln 𝐴 − ( ) [29]
𝑅 𝑇

𝑦 = 𝑚𝑥 + 𝑏
Donde:

y = ln 𝑘

1
𝑥=
𝑇
Se grafica X vs. Y para obtener m y b, a partir del coeficiente m se obtuvo la energía
de activación

𝑏 = ln 𝐴 [30]

𝐸𝑎
𝑚=−
𝑅

12
RESULTADOS Y DISCUSIÓN
Parte A

i) Alimentación de acetato de Etilo en exceso

A) Método integral

0.12000 Método integral

0.10000

Datos Exp
0.08000
Orden n
Ca [mol/L]

0.06000
Orden 1

0.04000 Orden 2

0.02000

0.00000
0 500 1000 1500 2000
t [s]

Fig A1 - Resultados para el método integral

Una vez hechas las regresiones lineales con diferentes órdenes se obtuvo la
gráfica anterior, que compara las concentraciones de NaOH calculadas con un
orden de reacción 1, 2 y de 1.0641.

Tabla 6 - Comparación de resultados para el método integral

Orden K Sr
1 0.01848 1.11E-03
1.0641 0.09606 1.44E-05
2 227.186 5.77E-04

13
 B) Método diferencial

Método Diferencial
0
-16 -14 -12 -10 -8 -6 -4
-2

y = 0.9272x + 2.8553
R² = 0.8728 -4

-6
Ln(-ra)

-8

-10

Ln(-ra) vs Ln(Ca) -12

-14
Ln(Ca)

Fig A2 - Resultados para el método diferencial

Para el método diferencial se encontró que el orden de reacción es de 0.9272 y la

velocidad de reacción es de 0.3536 [L/mol]-0.0728 s-1

ii) Alimentación equimolar

14
 A) Método integral

0.12000
Método integral

0.10000
Datos Exp
0.08000
Orden n
Ca [mol/L]

0.06000 Orden 1

Orden 2
0.04000

0.02000

0.00000
0 200 400 600 800 1000
t [s] 1200 1400 1600 1800 2000

Fig A3 - Resultados del método integral

Una vez hechas las regresiones lineales con diferentes órdenes se obtuvo la
gráfica anterior, que compara las concentraciones de NaOH calculadas con un
orden de reacción 1, 2 y de 1.5792.
Tabla 7 - Comparativa de resultados para el método integral

Orden k Sr
1 0.00362 0.2101
1.5792 0.1585 7.72E-03
2 2.4144 0.0248

15
 Método Diferencial
-5
-5.5 -5 -4.5 -4 -3.5 -3 -2.5
-6
y = 1.7052x - 2.4964
-7

-8
Ln(-ra)

-9

-10

-11

-12

-13
Ln(Ca)

Fig A1 - Resultados del método diferencial

En la gráfica anterior se muestran los resultados obtenidos con el método diferencial, en

el que se tuvo que hacer una regresión linear, con este método se obtuvo un orden global
de 1.7052 y un coeficiente de velocidad de reacción de 0.08325 [L/mol]0.2468 s-1
Parte B

Utilizando los parámetros siguientes, obtenidos en la parte A. Se graficaron las

concentraciones predichas de acuerdo con las concentraciones iniciales de los 3 casos
estudiados en esta parte y se compararon con las experimentales en las siguientes
graficas.

Alpha 1.0641
Beta 0.5151
K 0.1585 (Mol/L) ^-0.5792/s
Se obtuvieron los siguientes parámetros.

Tabla 8 - Parámetros obtenidos

Alpha 1.015
Beta 0.576
K 0.15572907 (Mol/L) ^-0.591/s

16
Con estos valores de los órdenes de reacción y el coeficiente cinético se calculó una
concentración de hidróxido de sodio y se calculó un error relativo total.

𝐶𝐴 𝑒𝑥𝑝 − 𝐶𝐴 𝑐𝑎𝑙𝑐
𝑒𝑟𝑟𝑜𝑟 𝑡𝑜𝑡𝑎𝑙 = ∑ ( ) ∗ 100 = 4.68%
𝐶𝐴 𝑒𝑥𝑝

Concentraciones experimentales vs
calculadas de corrida Acetato 0.10 M y
NaOH 0.05 M
0.06000
0.05000
Ca (mol/L)

0.04000
0.03000
Valores Experimentales
0.02000
Valores calculados
0.01000
0.00000
0 200 400 600
tiempo (s)

Concentraciones experimentales vs
calculadas de corrida Acetato 0.15 M y
NaOH 0.10 M
0.12000
0.10000
Ca (mol/L)

0.08000
0.06000
Valores Experimentales
0.04000
Valores calculados
0.02000
0.00000
0 100 200 300 400
tiempo (s)

Fig B1, B2 Predicciones de las concentraciones para las condiciones iniciales en el caso B.

A) 0.1M Acetato de Etilo y 0.05M de NaOH

17
 Concentraciones experimentales vs calculadas
0.06000

0.05000

0.04000
Ca (mol/L)

0.03000
Valores Experimentales
0.02000 Valores calculados

0.01000

0.00000
0 100 200 300 400 500 600
tiempo (s)

Fig A3 - Resultados para método diferencial con 0.1M Acetato y 0.05M de NaOH.

Tabla 9 - Comparación de concentraciones experimentales y calculadas para

corrida de 0.1M de Acetato de Etilo y 0.05M NaOH

Tiempo Ca Ca calculada
(s) experimental (mol/L)
(mol/L)
0 0.0500 0.0500
51 0.0404 0.0406
102 0.0319 0.0330
153 0.0263 0.0271
204 0.0218 0.0223
255 0.0181 0.0186
306 0.0146 0.0153
357 0.0124 0.0127
408 0.0106 0.0105
459 0.0089 0.0086
510 0.0072 0.0071
549 0.0059 0.0062

B) 0.05M Acetato de Etilo y 0.10M de NaOH

Para esta corrida se usó el pH

18
C) 0.15M Acetato de Etilo y 0.10M de NaOH

En primer lugar, usando los valores de las conductividades a través del tiempo medidas con
el potenciómetro, se determinó la concentración de hidróxido de sodio a través del tiempo
hasta llegar al estado estable. Obteniendo la siguiente tabla de datos.

Tabla 10 - Concentraciones experimentales para corrida de 0.15 M de Acetato de

Etilo y 0.10 M NaOH

Tiempo Conductividad Ca
(s) (S/m) (mol/L)
0 3.65 0.1000
51 1.307 0.0231
102 0.976 0.0123
153 0.816 0.0070
204 0.728 0.0041
255 0.677 0.0025
306 0.648 0.0015
354 0.632 0.0010

Método Diferencial

En el caso B no se puede considerar un reactivo en exceso, por lo que se tiene que

considerar que ambos reactivos tienen efecto en la tasa de reacción. Se escribe la tasa de
reacción y se le aplican logaritmos naturales para llegar a una ecuación lineal de 3
parámetros.

−𝑟𝐴 = 𝑘𝐶𝐴 𝛼 𝐶𝐵 𝛽

𝑑𝐶𝐴
ln −𝑟𝐴 = ln (− ) = ln 𝑘 +α ln 𝐶𝐴 + 𝛽 ln 𝐶𝐵
𝑑𝑡

Calculando la derivada numérica mediante diferencias finitas, y con las concentraciones

experimentales se realizó una regresión lineal múltiple para el cálculo de los órdenes de
reacción y el coeficiente cinético de reacción. Se obtuvieron los siguientes parámetros.

Tabla 11 - Parámetros obtenidos

Alpha 1.00448
Beta 0.52963
K 0.15270 (Mol/L) ^-0.5341/s
Con estos valores de los órdenes de reacción y el coeficiente cinético se calculó una concentración
de hidróxido de sodio y se calculó un error relativo total.

19
 𝐶𝐴 𝑒𝑥𝑝 − 𝐶𝐴 𝑐𝑎𝑙𝑐
𝑒𝑟𝑟𝑜𝑟 𝑡𝑜𝑡𝑎𝑙 = ∑ ( ) ∗ 100 = 16.79%
𝐶𝐴 𝑒𝑥𝑝

Concentraciones experimentales vs
calculadas
0.12000

0.10000

0.08000
Ca (mol/L)

0.06000
Valores Experimentales
0.04000
Valores calculados
0.02000

0.00000
0 100 200 300 400
-0.02000
tiempo (s)

Fig B4 - Resultados para método diferencial con 0.15 M Acetato y 0.10 M de NaOH.

Tabla 12 - Comparación de concentraciones experimentales y calculadas para

corrida de 0.15 M de Acetato de Etilo y 0.10 M NaOH

Tiempo Ca exp Ca calc

(s) (mol/L) (mol/L)
0 0.1000 0.1000
51 0.0231 0.0293
102 0.0112 0.0152
153 0.0059 0.0075
204 0.0035 0.0044
255 0.0020 0.0017
306 0.0011 0.0011
354 0.0006 0.0004

Parte C

A continuación, se presentan las gráficas de cada temperatura del experimento. Se

obtuvieron los valores del coeficiente cinético para diferentes temperaturas.

20
 1. 15 °C

CA1-n - CA01-n vs t
35

CA1-n - CA01-n [=] (mol/lt)-0.5

30 y = 0.0104x
R² = 0.97
25
20
15
10
5
0
0 500 1000 1500 2000 2500 3000
t [=] segundos

Figura C.1. Regresión lineal a 15°C

2. 19 °C

CA1-n - CA01-n vs t
20
CA1-n - CA01-n [=] (mol/lt)-0.5(seg)-1

18
16 y = 0.1104x
14 R² = 0.9771
12
10
8
6
4
2
0
0 20 40 60 80 100 120 140 160
t [=] segundos

Figura C.2. Regresión lineal a 19°C

3. 30 °C

21
 CA1-n - CA01-n vs t
30

CA1-n - CA01-n [=] (mol/lt)-0.5(seg)-1

25
y = 0.0821x
20 R² = 0.9913

15

10

0
0 50 100 150 200 250 300
-5
t [=] segundos

Figura C.3. Regresión lineal a 30°C

4. 35 °C

CA1-n - CA01-n vs t
45
CA1-n - CA01-n [=] (mol/lt)-0.5(seg)-1

40 y = 0.0361x
35 R² = 0.9896
30
25
20
15
10
5
0
0 200 400 600 800 1000 1200
-5
t [=] segundos

Figura C.4. Regresión lineal a 35°C

5. 40 °C

22
 CA1-n - CA01-n vs t
45

CA1-n - CA01-n [=] (mol/lt)-0.5(seg)-1

40 y = 0.0615x
35 R² = 0.9756
30
25
20
15
10
5
0
0 100 200 300 400 500 600
t [=] segundos

Figura C.5. Regresión lineal a 40°C

De acuerdo con estas regresiones, se obtienen los siguientes valores de

coeficiente cinético a diferentes temperaturas, estos se muestran en la tabla 1.

En la figura C.6 se muestra la tendencia que tienen los datos aplicando la

ecuación 2, se realizó una regresión lineal a esta gráfica para obtener la energía
de activación.

Tabla 13. Coeficiente cinético en función de la temperatura.

Temperatura [=] Coeficiente cinético “k” [=] (mol/lt)-0.5(seg)-1

K
288.15 0.0208
292.15 0.2208
303.15 0.153
308.15 0.0722
313.15 0.123

23
 ln k vs 1/T
0
0.00315 0.0032 0.00325 0.0033 0.00335 0.0034 0.00345 0.0035
-0.5
-1
-1.5
y = -5850.3x + 16.141
-2 R² = 0.5776
ln k

-2.5
-3
-3.5
-4
-4.5
1/T

Figura C.6. Regresión para la obtención de la energía de activación.

Se obtiene una energía de activación de 48.639 kJ/mol.

Discusiones

Parte A

Para la corrida en exceso de acetato de etilo con los resultados podemos ver una
gran diferencia en la exactitud de los métodos, quedando el método integral como
el mejor para predecir el orden de reacción, así como el coeficiente de velocidad de
reacción. Esta diferencia se le podría atribuir a el paso con el que se toman los datos
(cada 3 segundos), ya que entre más grande es éste, menos exactitud se tiene.
Según lo calculado el orden de reacción correcto para esta corrida es de 1.0641 con
un coeficiente de velocidad de 0.09606 [L/mol]0.0641 s-1.

Para la corrida equimolar comparado los resultados del método integral, contra el
diferencial se tiene que nuevamente el método integral tiene mayor exactitud que
el diferencial. Esto se puede deber a los tiempos largos a los que se dejó correr el
experimento y a que los datos se tomaran cada 3 segundos se afecte la exactitud
del método, en esta ocasión el método tiene un mejor ajuste ya que la reacción se
es más lenta que en la corrida pasada. Con los cálculos que se hicieron el orden
global de la reacción para la corrida equimolar queda como 1.5792 con un
coeficiente de velocidad de reacción de 0.1585 [L/mol]0.5792 s-1.

24
 Parte B

De acuerdo con los parámetros obtenidos en la parte A, se es capaz de predecir las

concentraciones a través del tiempo a diferentes condiciones iniciales, sin embargo,
esto no es cierto para todas las condiciones iniciales, tal es el caso de 0.15 M de
Acetato y 0.10 M de Hidróxido.

A pesar de esto los valores de el coeficiente cinético y los ordenes de reacción

encontrados con la regresión lineal múltiple, no son muy diferentes a los
encontrados en la parte A.

Parte C

Se obtuvo la energía de activación a partir de la experimentación de la reacción a

cinco distintas temperaturas. Se puede visualizar de la figura C.1 a la C.5 el
comportamiento de la reacción.

Se puede observar que en la figura C.6 el comportamiento es distinto al que

esperamos, puesto que hay un pico en la temperatura de 30 °C que afecta en el
ajuste. Esto representa que el experimento fue difícil de controlar a la temperatura
deseada, podría mejorarse con distinta agitación y con un sistema que controle
mejor la temperatura.

CONCLUSIÓN

Parte A

Parte B

En el artículo de “Estimation of Parameters of Arrhenius Equation for Ethyl Acetate

Saponification Reaction” se concluye que no es posible encontrar un order de
reacción para cada reactivo que se ajuste a todas las condiciones iniciales. Estos
ordenes y el coeficiente cinetico cambian de acuerdo a la temperatura de reacción,
condiciones de exceso y concentraciones iniciales. Esto va de acuerdo a nuestros
resultados donde los parametrós cineticos no son los mismos para todas las
condiciones, pero presentan una varianza leve.

Parte C

25
En el mismo artículo usado en la parte B, referencia 4, se muestra que la energía
de activación calculada es de 43.094 kJ/mol, mientras que nuestra experimentación
nos arrojó un valor de 48.63 kJ/mol. El error absoluto es de 12.86% Este porcentaje
de error se puede deber a que, los experimentos realizados en el artículo científico
se realizaron con un equipo mejor para el control de las variables como la
temperatura, y con eso disminuye la desviación en los resultados.

REFERENCIAS
1. Levenspiel, O. (2004). Ingeniería de las Reacciones Químicas 3a ed.. México: Limusa
Wiley. pp. 38-41.
2. Ramírez López, R., & Hernández Pérez, I. (2015). Diseño de reactores homogénos (1st
ed.). México: Cengage Learning Editores.
3. Gutiérrez Martín, F., Atanes Sánchez, E., & Cambra Pereira, A. (2016). Estudio cinético en
un reactor discontinuo de tanque agitado. Madrid: Universidad Politécnica de Madrid.
4. Mukhtar, Ahmad, Umar Shafiq, Ali Feroz Khan, Hafiz Abdul Qadir, and Masooma
Qizilbash, ‘Estimation of Parameters of Arrhenius Equation for Ethyl Acetate
Saponification Reaction’, Research Journal of Chemical Sciences Res. J. Chem. Sci, 5
(2015), 2231–2606.
5. Das, Kuheli & Sahoo, Parashuram & Magapu, Sai Baba & Murali, N & Swaminathan,
Pradeep. (2011). Kinetic Studies on Saponification of Ethyl Acetate Using an Innovative
Conductivity-Monitoring Instrument with a Pulsating Sensor. International Journal of
Chemical Kinetics. 19. 648-656. 10.1002/kin.20597.

26