UNIVERSIDAD DE CASTILLA- LA MANCHA

FACULTAD DE CIENCIAS Y TECNOLOGÍAS

QUÍMICAS
GRADO EN INGENIERÍA QUÍMICA

LABORATORIO INTEGRADO DE PROCESOS Y DE PRODUCTOS

REACCIÓN

GRUPO 2 TURNO 1:
Francisco Javier Benito Jurado
Paula Bravo García – Calvo
Lucía López Rivilla FEBRERO 2023
Índice

1. OBJETIVO .......................................................................................................................... 1
2. FUNDAMENTO TEÓRICO .............................................................................................. 1
3. MATERIALES .................................................................................................................... 1
4. PROCEDIMIENTO EXPERIMENTAL .......................................................................... 3
5. PRESENTACIÓN DE RESULTADOS............................................................................. 4
5.1. Estudio de la cinética en un reactor de laboratorio. ................................................... 4
5.2. Validación de la ecuación cinética y optimización en un reactor discontinuo de
planta piloto. ........................................................................................................................... 20
5.3. Diseño de reactores continuos. .................................................................................. 26
6. CONCLUSIONES ............................................................................................................. 27
7. BIBLIOGRAFÍA ............................................................................................................... 27
 1. OBJETIVO
El principal objetivo es la obtención de la ecuación cinética completa a escala de laboratorio de
la reacción de saponificación de acetato de etilo con sosa y posteriormente validarla en un reactor
discontinuo de la planta piloto de 40 L de volumen. También se diseñarán reactores discontinuos
y continuos a partir de la ecuación cinética obtenida.

2. FUNDAMENTO TEÓRICO
El ion hidroxilo del hidróxido de sodio posee una fuerte naturaleza. Por este motivo, la evolución
de la reacción puede seguirse midiendo la conductividad iónica en cada momento porque la
cantidad de sosa va a ser menor a medida que transcurre la reacción.
La ecuación cinética puede obtenerse a partir de consideraciones teóricas o mediante el resultado
de un procedimiento empírico de ajuste de curvas, aunque es importante saber que actualmente
el valor de las constantes de velocidad de dicha ecuación, tan solo se pueden determinar por vía
experimental porque los modelos teóricos no resultan válidos.
Los reactores discontinuos suelen trabajar a Tª y volumen constante, debido a esto son los
reactores preferentes para la obtención de los datos cinéticos. Para analizar dichos datos se puede
emplear tanto el método integral como el diferencial, en este caso se emplea el integral.1
El modelo de Arrhenius y el de Eyring relacionan la constante cinética k con la Tª. Ambos indican
que siempre a medida que aumenta la Tª, aumenta k. Los dos modelos serán utilizados en la
práctica para realizar los cálculos.
En cuanto al diseño de reactores, se tendrá en cuenta tanto el reactor mezcla perfecta como el
flujo pistón. El reactor mezcla perfecta indica que la composición se puede considerar uniforme
en todo el reactor; en cambio, en el reactor flujo pistón se trabaja en estado estacionario, las
propiedades no varían con el tiempo, pero sí con la posición en el reactor. Cada uno tiene sus
respectivas ecuaciones que se utilizarán posteriormente en el diseño de reactores.

3. MATERIALES
Se realizará la práctica a dos escalas: laboratorio y planta piloto.
En el laboratorio, la unidad está formada por lo siguiente:
- Matraz de reacción de 2 L termostatizado, colocado sobre un agitador magnético que
varía la velocidad de agitación y, además, está conectado por la parte superior con un
serpentín de refrigeración.
- Sistema de control de la temperatura de tipo ON/OFF basado en un termopar que está
conectado a un selector digital de temperatura y a la placa calefactora.
- Sistema de análisis de datos que está formado por un conductímetro

1
En la siguiente Ilustración 1 se muestra la instalación experimental a escala de laboratorio:

Ilustración 1. Instalación experimental a escala de laboratorio

En la planta piloto se reproduce a una escala mayor cómo funciona el reactor discontinuo de la
unidad del laboratorio, contiene las siguientes partes:
- Reactor termostatizado de 40 L con un motor de agitación.
- Sistema de control de la temperatura.
- Circuitos de alimentación de reactivos con unos depósitos.

La siguiente Ilustración 2 muestra la instalación experimental en la planta piloto:

Ilustración 2. Instalación experimental a escala de planta piloto

2
 4. PROCEDIMIENTO EXPERIMENTAL
Se van a realizar 3 experimentos con la finalidad de determinar la constante cinética de la reacción
a 20ºC, 35ºC y 45ºC. Se siguen los siguientes pasos:
1) Se conecta el paso de agua por el serpentín y se mantiene abierto toda la práctica, además
de introducir la sonda de conductividad en el reactor.
2) Se prepara una disolución de sosa 0,02 N en un matraz de 1 L y se añade al reactor.
3) Se establece la temperatura que se requiere en cada momento y cuando se consiga, se
pone en marcha el sistema de control de temperatura.
4) Se pone en marcha el agitador.
5) Cuando se mezcla la disolución y se estabiliza la temperatura, se debe anotar el dato de
conductividad inicial que marca el conductímetro y también el dato de temperatura.
6) Se añade el volumen necesario de acetato de etilo para que la concentración en el reactor
sea de 0,01M.
7) Cuando se añade el acetato se apunta el dato de la conductividad en ese momento. Se
pone en marcha el cronómetro.
8) Mientras se está produciendo la reacción se deben anotar los datos de conductividad cada
10 segundos durante los primeros minutos y cada 30 segundos cuando la variación de la
conductividad es casi nula.
9) El dato de conductividad a tiempo infinito se tomará transcurrido el tiempo suficiente
para que la reacción se estabilice.
10) Parar todo y limpiar.
11) Repetir a partir del paso 2 con las temperaturas que se indica en el guion.

Además, se realiza un experimento en planta piloto a una temperatura de 37ºC para comprobar la
influencia del tiempo de reacción tanto en la conversión experimental como en la teórica.
El procedimiento que se siguió es el siguiente:
1) Se conectó el paso de agua por el refrigerante.
2) A continuación, se cargó el reactor con agua previamente caliente, hasta la marca de 40
litros.
3) Se estableció la temperatura de consigna a 30ºC, aunque posteriormente se trabajó con
una temperatura de 37ºC, esto es debido a que la manta calefactora se encuentra caliente,
aunque no se siga calentando, haciendo que siga aumentando la temperatura.
4) Se añadió una disolución de NaOH con una concentración de 0,02 M. Esta se disolvió
previamente.
5) Se conectó el agitador. Una vez homogeneizada la solución, se tomó una muestra del
tanque con el objetivo de conocer la verdadera concentración de NaOH.
6) Cuando la temperatura llegó a 37 ºC se añadió acetato de etilo. La cantidad suficiente
para que la concentración de acetato de etilo en de la disolución fuera 0,01 M.
7) Una vez se hubo añadido se tomó tiempo cero, se esperó unos minutos y se comenzó a
tomar muestras de 10 ml.
8) Las muestras se valoraron con HCl 0,05 M empleando fenolftaleína como indicador para
determinar la concentración de NaOH.
9) Se finalizó una vez se consiguieron tres muestras con un mismo volumen de HCl.
10) Se cerró el paso de agua de refrigerante.

3
 5. PRESENTACIÓN DE RESULTADOS

5.1. Estudio de la cinética en un reactor de laboratorio.

Se presenta a continuación la Tabla 1 con los datos para el cálculo de la constante de velocidad
para una temperatura = 20 ºC
EXPERIMENTO A 20 ºC

T0 (ºC) 20
Tabla 1. Datos para el cálculo de la constante de velocidad a 20 °C

(1/(CA0·(M-
Tiempo 1))·ln((M-
µ (mS/cm) T (ᵒC) XA CA(mol/L) CB(mol/L)
(s) XA)/(M·(1-XA)))
(L/mol)
- 10,88 20 - - - -
0 10,85 20 0 0,01 0,02 0
10 10,71 - 0,07143 0,00929 0,01929 1,88702
20 10,66 - 0,09694 0,00903 0,01903 2,61408
30 10,59 - 0,13265 0,00867 0,01867 3,68439
40 10,54 - 0,15816 0,00842 0,01842 4,48927
50 10,48 - 0,18878 0,00811 0,01811 5,50332
60 10,43 - 0,21429 0,00786 0,01786 6,39167
70 10,38 - 0,23980 0,00760 0,01760 7,32255
80 10,33 - 0,26531 0,00735 0,01735 8,29926
90 10,29 - 0,28571 0,00714 0,01714 9,11608
100 10,24 - 0,31122 0,00689 0,01689 10,18482
110 10,19 - 0,33673 0,00663 0,01663 11,31079
120 10,15 - 0,35714 0,00643 0,01643 12,25612
130 10,12 - 0,37245 0,00628 0,01628 12,99298
140 10,08 - 0,39286 0,00607 0,01607 14,01510
150 10,04 - 0,41327 0,00587 0,01587 15,08568
160 10,01 - 0,42857 0,00571 0,01571 15,92269
170 9,97 - 0,44898 0,00551 0,01551 17,08746
180 9,94 - 0,46429 0,00536 0,01536 18,00014
190 9,91 - 0,47959 0,00520 0,01520 18,94867
200 9,88 - 0,49490 0,00505 0,01505 19,93542
210 9,85 - 0,51020 0,00490 0,01490 20,96292
220 9,82 21 0,52551 0,00474 0,01474 22,03400
230 9,8 - 0,53571 0,00464 0,01464 22,77378
240 9,77 - 0,55102 0,00449 0,01449 23,92451
250 9,75 - 0,56122 0,00439 0,01439 24,72063
260 9,72 - 0,57653 0,00423 0,01423 25,96120
270 9,7 - 0,58673 0,00413 0,01413 26,82106
280 9,68 - 0,59694 0,00403 0,01403 27,70880
290 9,65 - 0,61224 0,00388 0,01388 29,09608
300 9,63 - 0,62245 0,00378 0,01378 30,06048
310 9,61 - 0,63265 0,00367 0,01367 31,05868
320 9,59 - 0,64286 0,00357 0,01357 32,09269

4
 330 9,58 - 0,64796 0,00352 0,01352 32,62381
340 9,56 - 0,65816 0,00342 0,01342 33,71571
350 9,54 - 0,66837 0,00332 0,01332 34,84930
360 9,52 - 0,67857 0,00321 0,01321 36,02731
370 9,51 - 0,68367 0,00316 0,01316 36,63390
380 9,49 - 0,69388 0,00306 0,01306 37,88429
390 9,48 - 0,69898 0,00301 0,01301 38,52895
400 9,46 - 0,70918 0,00291 0,01291 39,85955
410 9,45 - 0,71429 0,00286 0,01286 40,54651
420 9,44 - 0,71939 0,00281 0,01281 41,24863
450 9,4 - 0,73980 0,00260 0,01260 44,22078
480 9,36 - 0,76020 0,00240 0,01240 47,48833
510 9,33 - 0,77551 0,00224 0,01224 50,16511
540 9,31 - 0,78571 0,00214 0,01214 52,07269
570 9,28 - 0,80102 0,00199 0,01199 55,14383
600 9,26 - 0,81122 0,00189 0,01189 57,34867
630 9,24 - 0,82143 0,00179 0,01179 59,69612
660 9,22 - 0,83163 0,00168 0,01168 62,20336
690 9,2 - 0,84184 0,00158 0,01158 64,89078
720 9,17 21,6 0,85714 0,00143 0,01143 69,31472
750 9,16 - 0,86224 0,00138 0,01138 70,90939
780 9,15 - 0,86735 0,00133 0,01133 72,57168
810 9,14 - 0,87245 0,00128 0,01128 74,30698
840 9,13 - 0,87755 0,00122 0,01122 76,12133
870 9,12 - 0,88265 0,00117 0,01117 78,02152
900 9,11 - 0,88776 0,00112 0,01112 80,01527
930 9,11 - 0,88776 0,00112 0,01112 80,01527
960 9,1 - 0,89286 0,00107 0,01107 82,11139
990 9,1 - 0,89286 0,00107 0,01107 82,11139
1020 9,09 21,7 0,89796 0,00102 0,01102 84,31995
1050 9,09 - 0,89796 0,00102 0,01102 84,31995
1080 9,08 - 0,90306 0,00097 0,01097 86,65259
1110 9,08 - 0,90306 0,00097 0,01097 86,65259
1140 9,07 - 0,90816 0,00092 0,01092 89,12285
1170 9,07 - 0,90816 0,00092 0,01092 89,12285
2370 8,89 19,8 1 0 0,01 -

En el caso de este experimento, los gramos de NaOH que se precisan son de 0,8 gramos, pero
debido a un error de la báscula, el peso que marcaba era un peso de 0,86 gramos, el cual es un
error aceptable. Sin embargo, ese no era el peso real, y esto se puede contemplar a la hora de
medir la conductividad, debido a que son valores demasiado altos para una supuesta
concentración de 0,02 M.
A continuación, se representan los datos de conductividad medidos experimentalmente con
respecto al tiempo, incluido el dato de conductividad a tiempo infinito, el cual se alcanza cuando
la reacción se ha finalizado. Estos datos están recogidos en el Gráfico 1.

5
 10,8

10,3
Conductividad (mS/cm)

9,8

9,3

8,8
0 500 1000 1500 2000
tiempo (s)

Gráfico 1. Datos experimentales de conductividad frente al tiempo a 20 °C

Se puede observar que los datos presentan una tendencia lineal en los primeros valores de tiempo,
y transcurridos aproximadamente 350 segundos, comienza a estabilizarse hasta llegar al punto de
tiempo infinito, donde la reacción ha finalizado. Estos datos se compararán posteriormente con
los obtenidos a distintas temperaturas.
Para el cálculo de la constante de equilibrio, hay que hacer unos cálculos previos, de los cuales se
mostrará un ejemplo de cálculo para cada uno de los parámetros para el experimento a 45 ºC.
En el Gráfico 2 se muestra la representación de los datos mencionados anteriormente, donde
únicamente se ha representado la zona lineal para, con la pendiente, obtener el valor de k
(L/mol·s). Posteriormente se recogerán los valores de k a distintas temperaturas en la Tabla 4.

45
((1/(Cao·(M-1))·ln((M-Xa)/(M·(1-Xa)))

40

35

30
(L/mol)

25

20

15 y = 0,0997x
R² = 0,9998
10

0
0 50 100 150 200 250 300 350 400
tiempo (s)

Gráfico 2. Obtención de constante de reacción a 20ºC

6
EXPERIMENTO A 35 ºC
Para los datos de temperaturas de 35 ºC, se hacen los mismos cálculos, los cuales están indicados
posteriormente para el ejemplo de 45 ºC.
Se muestran, al igual que para 20 ºC, los datos recogidos y calculados en la Tabla 2.

T0 (ºC) 35

Tabla 2. Datos para el cálculo de la constante de velocidad a 35 ºC

(1/(CA0·(M-
1))·ln((M-
Tiempo CB
µ (mS/cm) T (ᵒC) XA CA(mol/L) XA)/(M·(1-
(s) (mol/L)
XA)))
(L/mol)
- 4,33 35 - - - -
0 4,27 - 0 0,01 0,02 0
10 4,12 - 0,1000 0,0090 0,0190 5,40672
20 4,09 - 0,12000 0,00880 0,01880 6,59580
30 4,04 - 0,15333 0,00847 0,01847 8,66832
40 3,99 - 0,18667 0,00813 0,01813 10,86338
50 3,96 - 0,20667 0,00793 0,01793 12,24407
60 3,93 - 0,22667 0,00773 0,01773 13,67589
70 3,89 - 0,25333 0,00747 0,01747 15,66985
80 3,86 - 0,27333 0,00727 0,01727 17,23330
90 3,83 - 0,29333 0,00707 0,01707 18,85912
100 3,8 - 0,31333 0,00687 0,01687 20,55133
110 3,77 - 0,33333 0,00667 0,01667 22,31436
120 3,74 - 0,35333 0,00647 0,01647 24,15302
130 3,71 - 0,37333 0,00627 0,01627 26,07263
140 3,68 - 0,39333 0,00607 0,01607 28,07902
150 3,66 - 0,40667 0,00593 0,01593 29,46800
160 3,63 - 0,42667 0,00573 0,01573 31,63373
170 3,61 - 0,44000 0,00560 0,01560 33,13571
180 3,59 - 0,45333 0,00547 0,01547 34,68709
190 3,57 36,4 0,46667 0,00533 0,01533 36,29055
200 3,55 - 0,48000 0,00520 0,01520 37,94896
210 3,53 - 0,49333 0,00507 0,01507 39,66545
220 3,51 - 0,50667 0,00493 0,01493 41,44338
230 3,49 - 0,52000 0,00480 0,01480 43,28641
240 3,47 - 0,53333 0,00467 0,01467 45,19851
250 3,45 - 0,54667 0,00453 0,01453 47,18402
260 3,44 - 0,55333 0,00447 0,01447 48,20576
270 3,42 - 0,56667 0,00433 0,01433 50,31036
280 3,4 - 0,58000 0,00420 0,01420 52,50103
290 3,39 - 0,58667 0,00413 0,01413 53,63047
300 3,37 - 0,60000 0,00400 0,01400 55,96158
310 3,36 36,3 0,60667 0,00393 0,01393 57,16496
320 3,34 - 0,62000 0,00380 0,01380 59,65203
330 3,33 - 0,62667 0,00373 0,01373 60,93773

7
 340 3,32 - 0,63333 0,00367 0,01367 62,25296
350 3,31 - 0,64000 0,00360 0,01360 63,59888
360 3,3 - 0,64667 0,00353 0,01353 64,97669
370 3,28 - 0,66000 0,00340 0,01340 67,83321
380 3,27 - 0,66667 0,00333 0,01333 69,31472
390 3,26 - 0,67333 0,00327 0,01327 70,83373
400 3,25 - 0,68000 0,00320 0,01320 72,39188
410 3,24 - 0,68667 0,00313 0,01313 73,99089
420 3,23 - 0,69333 0,00307 0,01307 75,63261
450 3,2 36,2 0,71333 0,00287 0,01287 80,83429
480 3,17 - 0,73333 0,00267 0,01267 86,49974
510 3,15 - 0,74667 0,00253 0,01253 90,57086
540 3,13 36,1 0,76000 0,00240 0,01240 94,90806
570 3,11 - 0,77333 0,00227 0,01227 99,54281
600 3,09 36,0 0,78667 0,00213 0,01213 104,51236
630 3,06 - 0,80667 0,00193 0,01193 112,69428
660 3,05 35,9 0,81333 0,00187 0,01187 115,64319
690 3,04 - 0,82000 0,00180 0,01180 118,71657
720 3,03 - 0,82667 0,00173 0,01173 121,92403
750 3,01 35,8 0,84000 0,00160 0,01160 128,78543
780 3 - 0,84667 0,00153 0,01153 132,46502
810 2,99 - 0,85333 0,00147 0,01147 136,33048
840 2,98 - 0,86000 0,00140 0,01140 140,39939
870 2,98 - 0,86000 0,00140 0,01140 140,39939
900 2,97 35,7 0,86667 0,00133 0,01133 144,69190
930 2,96 - 0,87333 0,00127 0,01127 149,23126
1950 2,77 25 1 0 0,01 -

En este caso, no hubo errores de la báscula, puesto a que los valores de conductividad se
encuentran en un rango adecuado. Todos los cálculos se realizan del mismo modo que en los
casos de distintas temperaturas. Se adjunta el Gráfico 3, donde se representan los datos de
conductividad con respecto al tiempo.

4,3

4,1

3,9
Conductividad (mS/cm)

3,7

3,5

3,3

3,1

2,9

2,7
0 200 400 600 800 1000 1200 1400 1600 1800 2000
tiempo (s)

Gráfico 3. Datos experimentales de conductividad frente al tiempo a 35 ºC

8
En este caso, al aumentar la temperatura, la velocidad de reacción es más rápida, por lo que el
tiempo que tarda en estabilizarse es menor. Esto se ve reflejado en que, la zona lineal tarda menos
en estabilizarse, ya que transcurridos aproximadamente 250 segundos, empieza a estabilizarse.
Esto indica que la pendiente es más pronunciada, y, por ende, la reacción es más rápida.
En el Gráfico 4, se hace la representación para poder calcular la k (L/mol·s) a 35 ºC, y, al igual
que en el caso anterior, se ha representado solo la zona lineal.

80
((1/(Cao·(M-1))·ln((M-Xa)/(M·(1-Xa)))

70

60

50
(L/mol)

40

30
y = 0,1861x
20 R² = 0,9984

10

0
0 50 100 150 200 250 300 350 400
tiempo (s)

Gráfico 4. Obtención de constante de reacción a 35 ºC

En este caso, se puede comprobar que, como se menciona anteriormente, la constante de velocidad
de reacción es superior al caso de 20 ºC.

9
EXPERIMENTO A 45 ºC
En este experimento el procedimiento de cálculo es el mismo y, a continuación, se muestran todos
los datos obtenidos de conductividad y los datos de conversión y concentración calculados. En
este caso, habrá un ejemplo de cálculo para cada uno de los parámetros en el final de este
experimento. Se adjuntan todos estos datos en la Tabla 3.

T0 (ºC) 45

Tabla 3. Datos para el cálculo de la constante de velocidad a 45 ºC

(1/(CAo·(M
-1))·ln((M-
Tiempo
µ (mS/cm) T (ᵒC) XA CA(mol/L) CB (mol/L) XA)/(M·(1-
(s)
XA)))
(L/mol)
- 4,06 45 - - - -
0 3,99 - 0 0,01 0,02 0
10 3,92 - 0,04930 0,00951 0,01951 2,55955
20 3,86 - 0,09155 0,00908 0,01908 4,91592
30 3,8 - 0,13380 0,00866 0,01866 7,43983
40 3,74 - 0,17606 0,00824 0,01824 10,15069
50 3,68 - 0,21831 0,00782 0,01782 13,07121
60 3,64 - 0,24648 0,00754 0,01754 15,14769
70 3,54 - 0,31690 0,00683 0,01683 20,86054
80 3,5 - 0,34507 0,00655 0,01655 23,38388
90 3,46 - 0,37324 0,00627 0,01627 26,06342
100 3,41 - 0,40845 0,00592 0,01592 29,65710
110 3,4 - 0,41549 0,00585 0,01585 30,41126
120 3,36 - 0,44366 0,00556 0,01556 33,55677
130 3,33 - 0,46479 0,00535 0,01535 36,06145
140 3,31 - 0,47887 0,00521 0,01521 37,80661
150 3,28 - 0,50000 0,00500 0,01500 40,54651
160 3,25 - 0,52113 0,00479 0,01479 43,44526
170 3,23 - 0,53521 0,00465 0,01465 45,47362
180 3,21 - 0,54930 0,00451 0,01451 47,58459
190 3,19 - 0,56338 0,00437 0,01437 49,78384
200 3,17 - 0,57746 0,00423 0,01423 52,07760
210 3,15 - 0,59155 0,00408 0,01408 54,47272
220 3,13 - 0,60563 0,00394 0,01394 56,97682
230 3,11 - 0,61972 0,00380 0,01380 59,59834
240 3,09 - 0,63380 0,00366 0,01366 62,34673
250 3,07 - 0,64789 0,00352 0,01352 65,23252
260 3,06 - 0,65493 0,00345 0,01345 66,73059
270 3,05 - 0,66197 0,00338 0,01338 68,26759
280 3,03 - 0,67606 0,00324 0,01324 71,46534
290 3,02 - 0,68310 0,00317 0,01317 73,12989
300 3 - 0,69718 0,00303 0,01303 76,60085
310 2,99 - 0,70423 0,00296 0,01296 78,41190
320 2,98 - 0,71127 0,00289 0,01289 80,27669
330 2,97 - 0,71831 0,00282 0,01282 82,19801

10
 340 2,96 - 0,72535 0,00275 0,01275 84,17882
350 2,95 - 0,73239 0,00268 0,01268 86,22235
360 2,94 - 0,73944 0,00261 0,01261 88,33207
390 2,91 - 0,76056 0,00239 0,01239 95,09763
420 2,88 - 0,78169 0,00218 0,01218 102,61572
450 2,86 - 0,79577 0,00204 0,01204 108,12205
480 2,84 - 0,80986 0,00190 0,01190 114,09147
510 2,82 - 0,82394 0,00176 0,01176 120,59708
540 2,81 - 0,83099 0,00169 0,01169 124,07868
570 2,79 - 0,84507 0,00155 0,01155 131,56768
600 2,78 - 0,85211 0,00148 0,01148 135,60806
630 2,77 - 0,85915 0,00141 0,01141 139,87169
660 2,76 - 0,86620 0,00134 0,01134 144,38182
690 2,74 - 0,88028 0,00120 0,01120 154,25437
720 2,74 - 0,88028 0,00120 0,01120 154,25437
750 2,73 - 0,88732 0,00113 0,01113 159,68591
780 2,72 - 0,89437 0,00106 0,01106 165,50484
810 2,71 - 0,90141 0,00099 0,01099 171,76515
840 2,71 - 0,90141 0,00099 0,01099 171,76515
870 2,7 - 0,90845 0,00092 0,01092 178,53286
900 2,7 - 0,90845 0,00092 0,01092 178,53286
1800 2,57 34,6 1 0 0,01 -

En el Gráfico 5, se adjuntan los datos de conductividad frente al tiempo para este caso.

4,1

3,9

3,7
Conductividad (mS/cm)

3,5

3,3

3,1

2,9

2,7

2,5
0 200 400 600 800 1000 1200 1400 1600 1800
tiempo (s)

Gráfico 5. Datos experimentales de conductividad frente al tiempo a 45 °C

Como se ve en la representación, la pendiente es más pronunciada, lo que indica que la reacción

a esta temperatura es más rápida, y como se ve, empieza a estabilizarse a un tiempo de 140
segundos aproximadamente, esto se ve en el valor de k a esta temperatura, que será un valor mayor
que en los casos anteriores.

11
En el Gráfico 6, se muestra la representación para el cálculo de dicha constante, representando
solo la zona lineal de los datos, puesto que la constante es igual a la pendiente de esta recta.
((1/(Cao·(M-1))·ln((M-Xa)/(M·(1- 120

100

80
Xa))) (L/mol)

60

40 y = 0,2536x
R² = 0,9986
20

0
0 50 100 150 200 250 300 350 400 450
tiempo (s)

Gráfico 6. Obtención de constante de reacción a 45 ºC

Se puede contemplar que esta constante es mucho mayor que a temperaturas más bajas.
Comparando y representado la conductividad frente al tiempo a distintas temperaturas, como en
el Gráfico 7, se ve como varía la velocidad de caída de conductividad a cada temperatura.

12

10
Conductividad (mS/cm)

20 ºC
4
35 ºC

2 45 ºC

0
0 200 400 600 800 1000 1200
tiempo (s)

Gráfico 7. Comparativa de caída de conductividad a distintas Tas

Teniendo en cuenta el error mencionado para el experimento de 20 ºC con respecto a la

concentración real, se ve que, para una temperatura mayor, la conductividad empieza a
estabilizarse mucho antes.
El procedimiento realizado para calcular la constante de velocidad ha sido el mismo para cada
temperatura únicamente cambian los datos de conductividad tomados a cada una de estas.
Se realiza un ejemplo de cálculo de cada una de los parámetros para la temperatura de 45 ºC:

12
13
14
Los distintos valores de k se recogen en la Tabla 4 para cada una de las temperaturas:
Tabla 4. Datos de la constante de velocidad para distinta temperatura y para Arrhenius

Temperatura (ºC) 20 35 45
k (L/mol·s) 0,0997 0,1861 0,2536
ln(k) -2,30559 -1,68147 -1,37200
1/T 0,00341 0,00325 0,00314

Como se contempla en esta tabla, los valores de la constante aumentan con la temperatura.
Para relacionar este cambio en la constante de velocidad con la temperatura, hay dos
metodologías, una de ellas es mediante Arrhenius, y otra mediante Eyring.

ARRHENIUS

15
La representación se muestra en el Gráfico 8:
-1,2

-1,4

-1,6
ln (k)

-1,8 y = -3509,6x + 9,6778

R² = 0,997
-2,0

-2,2

-2,4
0,00310 0,00315 0,00320 0,00325 0,00330 0,00335 0,00340 0,00345
1/T (K-1)

Gráfico 8. Método de Arrhenius

A continuación, se muestran los cálculos necesarios para la energía de activación y para el factor
pre-exponencial:

Tabla 5. Datos de Ea y A

EA(J/mol·K) 29178,8144
A(L/mol·s) 15959,34771
R2 0,997

Si se compara este valor de EA obtenida con el extraído de la bibliografía, se observa que hay una
diferencia considerable, pues el valor de EA para una reacción de saponificación obtenido de la
bibliografía es de 48367,04 J/mol·K. Esto se puede deber al error cometido en el experimento a
20ºC.

16
EYRING

En este método, se precisa representar el ln(k/T) frente a 1/T, para ellos se precisa de los estos
datos, mostrados en la Tabla 5.
Tabla 6. Datos para método de Eyring

Temperatura (ºC) 20 35 45
k (L/mol·s) 0,0997 0,1861 0,2536
ln(k/T) -7,98627 -7,41206 -7,13452
1/T 0,00341 0,00325 0,00314

17
En el Gráfico 9, se hace la representación mencionada anteriormente:

-7,0
-7,1
-7,2
-7,3
-7,4
ln (k)

-7,5
y = -3204,7x + 2,9574
-7,6 R² = 0,9962
-7,7
-7,8
-7,9
-8,0
-8,1
0,00310 0,00315 0,00320 0,00325 0,00330 0,00335 0,00340 0,00345
1/T (K-1)

Gráfico 9. Método de Eyring

Los cálculos para la entalpía y entropía en este método se muestran a continuación:

18
El resumen de los datos obtenidos se muestra en la Tabla 7.
Tabla 7. Datos de ΔH≠ y ΔS≠

kB (J/K) 1,38E-23
h(J·s) 6,626E-34
ΔH≠(J/mol) 26643,8758

ΔS≠(J/(mol·K)) 584,1948428

R2 0,9962

En ambos casos, con los dos métodos, el coeficiente de correlación es bastante bueno ya que es
un valor superior a 0.99, por lo que se ajusta de buen modo a una recta.
La ecuación cinética completa de la reacción de saponificación en función de la temperatura
mediante Arrhenius es la siguiente:

19
 5.2. Validación de la ecuación cinética y optimización en un reactor discontinuo de
planta piloto.

Se ha realizado un experimento en planta piloto a una temperatura de 37 ºC con el objetivo de

observar la influencia del tiempo de reacción tanto en la conversión experimental como en la
conversión teórica.
El procedimiento realizado es semejante al realizado a escala de laboratorio.

20
En la Tabla 8 se adjuntan el volumen de HCl obtenido en las valoraciones realizadas a tres
muestras para conocer la concentración real inicial de NaOH.
Tabla 8. Datos de volumen de HCl recogidos en cada valoración.

Muestra ml HCL 0,05 M

1 3,6
2 3,3
3 3,4
Media 3,43

Una vez se añade el acetato de etilo se toman muestras cuyo volumen = 10 ml para valorarlas y
así obtener la concentración de NaOH a diferentes tiempos.
Para la validación de la ecuación cinética obtenida, se han empleado los datos de NaOH obtenidos
de las valoraciones, a partir de los cuales se puede calcular la concentración de acetato y la
conversión de este a lo largo del tiempo.
Los datos obtenidos se recogen en la siguiente tabla:
Tabla 9. Datos obtenidos en reactor de planta piloto régimen isotermo.

CA
t(s) ml HCl CB (mol/L) XA exp
(mol/L)
0 3,433 0,017167 0,00767 0
120 2,9 0,0145 0,005 0,3478
240 2,5 0,0125 0,003 0,6087
319,8 2,3 0,0115 0,002 0,7391
372 2,3 0,0115 0,002 0,7391
438 2,2 0,0110 0,0015 0,8043
494 2,1 0,0105 0,001 0,8696
551 2,05 0,0103 0,00075 0,9022
610 2 0,0100 0,0005 0,9348
665 2 0,0100 0,0005 0,9348
728 1,9 0,0095 0 1
780 1,9 0,0095 0 1
840 1,9 0,0095 0 1

Los cálculos realizados para completar la Tabla 9 se presentan a continuación:

21
22
Para el cálculo de la conversión teórica se ha empleado la ecuación de diseño de un reactor
discontinuo que se obtendría con el tiempo en el que ha transcurrido la reacción en planta piloto.
Los datos se recogen en la Tabla 10.
Tabla 10. Comparativa de conversión experimental y teórica.

t(s) XA exp XA teórica

0 0 0
120 0,3478 0,3095
240 0,6087 0,5019
319,8 0,7391 0,5926
372 0,7391 0,6409
438 0,8043 0,6923
494 0,8696 0,7291
551 0,9022 0,7614
610 0,9348 0,7902
665 0,9348 0,8135
728 1 0,8367
780 1 0,8534
840 1 0,8704

23
Los cálculos realizados para completar la Tabla 10 son los siguientes:

Si se representan ambas conversiones frente al tiempo se obtiene el siguiente gráfico:

Gráfico 10. Conversión experimental y teórica vs tiempo.

24
Como se puede observar en el Gráfico 10 la conversión teórica como experimental presentan la
misma tendencia. Sin embargo, la experimental presenta una velocidad mayor que la teórica. Esto
puede ser debido al error en el primer proceso experimental a escala laboratorio a la temperatura
de 20 ºC. Aunque como se ha podido comprobar el error no es demasiado grande a la hora de
ajustar la constante cinética, sí que puede influir de una manera más acusada en el cálculo de la
conversión teórica.

A continuación, se analiza el tiempo de no producción para ello se compara la conversión y

tiempos de operación óptimos para dos tiempos de carga y descarga distintos de un reactor
discontinuo.
Tiempo de no producción = 20 min.

Gráfico 11. Representación del método gráfico para la optimización de la producción de un reactor
discontinuo para un tnp=20 min.

Tiempo de no producción = 40 min.

Gráfico 12. Representación del método gráfico para la optimización de la producción de un reactor
discontinuo para un tnp=40 min.

25
Como se puede observar, para un valor de tiempo de no producción, el valor del tiempo óptimo y
de la conversión óptima es mayor que para un tiempo de no producción menor. El motivo se debe
a que cuando el tiempo de carga y descarga es mayor, normalmente vale la pena emplear más
tiempo de operación y conseguir mayores conversiones que cuando los tiempos de carga y
descarga son pequeños. Pues cuando son pequeños, se puede operar más veces consiguiendo
conversiones pequeñas a tiempos menores ya que el esfuerzo es menor que cuando los tiempos
de carga y descarga son altos y por tanto, es preferible operar menos veces y conseguir la
conversión deseada en un solo proceso.

5.3. Diseño de reactores continuos.

En este caso, se seleccionaría el reactor de Flujo Pistón, puesto que su tiempo espacial es menor
que el de Mezcla Perfecta.
La razón de que esto sea así es debido a la definición de tiempo espacial, que es el tiempo que se
tarda en tratar un volumen igual al volumen del reactor. Como en este caso queremos una
conversión determinada, los reactores de Mezcla Perfecta tienen un mayor volumen que los de
Flujo Pistón a una misma conversión, por lo que, si el volumen es mayor, también lo será su
tiempo espacial.
Además de eso, si se trata un mayor caudal o se desea una mayor conversión, esta diferencia será
más pronunciada.

26
 6. CONCLUSIONES
Como se puede contemplar en esta práctica, la variación de la velocidad de una reacción con
respecto a la temperatura se confirma gracias a la ley de Arrhenius. Esto se ha demostrado
siguiendo el procedimiento experimental mostrado anteriormente.
Exceptuando los errores cometidos en el experimento a 20 ºC, los cuales afectan al cálculo de la
conversión teórica en el reactor de 40 litros, puesto a que un error en el cálculo de una de las
constantes genera una ecuación de Arrhenius que no se corresponde del todo con la realidad,
aunque sí que se ajusta a una recta, puesto que, teóricamente la conversión sería más baja de lo
que se ha podido comprobar experimentalmente. De todos modos, la tendencia es la misma tanto
teóricamente como experimentalmente, por lo que, por esa parte, se confirma que se ajusta de un
modo correcto a los datos experimentales. Esto permite concluir que, si no fuera por los errores
cometidos, la ecuación de velocidad sí que sería reproducible y se puede validar para un reactor
discontinuo de planta piloto.
En cuanto a la relación de los tiempos de no producción con la producción global, se puede afirmar
que, por norma general, si estos tiempos de no producción son muy elevados, la conversión y el
tiempo óptimos serán más grandes que si estos tiempos de producción son menores. Esto se debe
a que en ese tiempo muerto no se está produciendo nada, por lo que se opta por dejar que la
reacción transcurra durante más tiempo, ya que su productividad se verá incrementada.
En el caso de que los tiempos de no producción fuera bajos, interesaría detener antes la reacción
pese a que la conversión sea menor que en el caso anterior. Esto significaría que de manera global
se puede hacer más cantidad de producto en experimentos más cortos, aunque tenga menos
pureza.
Siempre hay que encontrar un óptimo entre estos dos casos, porque cada caso depende de la
producción requerida.
Otra de las conclusiones es que, a una conversión dada, los reactores continuos de flujo pistón
tendrán un tiempo espacial menor que los de mezcla perfecta, puesto a que estos primeros precisan
de un volumen menor, y, sabiendo la definición de tiempo espacial, se llega a esta conclusión

7. BIBLIOGRAFÍA
Fogler, H.A. (1991) Elements of Chemical Reaction Engineering. Prentice Hall, Englewood
Cliffs, Nueva York.
Levenspiel, O. (1999) Chemical Reactor Engineering. John Wiley, Nueva York.

27