Reporte Equilibrio Quimico y Reactor Por Lotes Equipo - 2
Reporte Equilibrio Quimico y Reactor Por Lotes Equipo - 2
Reporte Equilibrio Quimico y Reactor Por Lotes Equipo - 2
Ingeniería Química
P á g i n a 1 | 23
Índice
3
Introducción
Objetivo
4
Problema experimental
Generalidades 6
Resultados
Análisis de resultados 9
Conclusión
Memoria de Calculo 10
Anexo
Bibliografía 21
P á g i n a 2 | 23
Introducción
La expresión matemática que representa al Equilibrio Químico se conoce como Ley de Acción
de Masas y se enuncia como: “La relación del producto de las actividades (actividad igual a
concentración en soluciones diluidas) elevadas los coeficientes estequiométricos en la
reacción de productos y reactivos permanece constante al equilibrio.”
Para conocer la constante de equilibrio de una reacción, es necesario conocer las condiciones
a las cuales se está expuesta dicha reacción, de igual forma, considerar la cinética química de
una reacción.
P á g i n a 3 | 23
Objetivo
Establecer y realizar el estudio de la reacción en condiciones controladas de volumen y
temperatura a fin de determinar los parámetros de la ecuación cinética correspondiente.
Definir las variables involucradas en la operación del Reactor por Lotes, así como establecer
las condiciones óptimas de operación del reactor para un tiempo establecido.
Problema Experimental
Determinar experimentalmente la distribución de reactivos y productos para la reacción de
saponificación de acetato de etilo y de butilo, corroborando la dependencia de la misma con la
temperatura y las composiciones al equilibrio, así como los parámetros cinéticos (orden de
reacción, constante de velocidad, energía de activación y factor de frecuencia).
Determinar el cambio de la concentración de uno de los componentes en función del tiempo
durante la saponificación del acetato de etilo utilizando como propiedades la conductividad del
efluente del reactor.
Generalidades
Todas las reacciones tienden al estado de equilibrio químico, el punto en el cual tanto el
proceso hacia adelante como el proceso inverso suceden a la misma velocidad. Dado que las
velocidades hacia adelante y hacia atrás son iguales, las concentraciones de los reactivos y
productos son constantes en el equilibrio. De acuerdo con las concentraciones de todas las
especies de reacción en equilibrio, podemos definir una cantidad llamada la constante de
equilibrio Kc que a veces también se escribe como Keq o K. La c en el subíndice significa
concentración, dado que la constante de equilibrio describe las concentraciones molares, en
moles/L en el equilibrio para una temperatura específica. La constante de equilibrio nos ayuda
a entender si una reacción tiende a mostrar una concentración más alta de productos o de
reactivos en el equilibrio. También podemos usar Kc para determinar si la reacción se
encuentra o no en el equilibrio.
Hay algunas cosas importantes que recordar cuando se calcula de kc:
Los principios básicos para la elección de las condiciones para llevar a cabo una reacción son
el equilibrio la cinética química que es el área de la química que estudia la velocidad o rapidez
con la que ocurre una reacción. La importancia de la cinética química abarca dos aspectos:
P á g i n a 4 | 23
1. Predecir la velocidad que tendrá una reacción en unas condiciones determinadas de
presión, temperatura, concentración, catalizador, etc.
2. Determinar y comprender el mecanismo por el que tiene lugar una reacción.
La velocidad de una reacción química está afectada por varios factores, de los cuales los más
sobresalientes son:
1. Concentración de los reactivos.
2. Temperatura.
3. Catalizador.
4. Estado físico de los reactivos.
La ecuación de velocidad nos indica cómo cambia la velocidad de una reacción a una
temperatura específica conforme modificamos las concentraciones de los reactivos. La
ecuación de velocidad para cualquier reacción química debe determinarse de forma
experimental.
La velocidad de casi todas las reacciones químicas aumenta al aumentar la temperatura.
Puesto que la temperatura no afecta mucho la concentración, la constante de velocidad de
reacción debe aumentar al elevarse la temperatura. La constante de velocidad y por tanto la
velocidad aumenta con la temperatura. Como regla general un aumento de 10 grados en la
temperatura duplica la velocidad de reacción.
Arrhenius encontró que la mayor parte de los datos de velocidades de reacción obedecen la
siguiente ecuación:
P á g i n a 5 | 23
Metodología de trabajo
Resultados
P á g i n a 6 | 23
Tiempo(s) Conductividad (µs)
0 6.3
30 5.84
60 5.67
90 5.5
120 5.38
180 5.26
240 5.15
300 5.05
360 4.95
420 4.88
480 4.8
540 4.73
600 4.67
660 4.6
720 4.55
780 4.5
840 4.47
900 4.43
960 4.37
1020 4.33
1080 4.29
1140 4.26
1200 4.22
1260 4.15
1320 4.16
1380 4.12
1440 4.05
1500 4.07
1560 4.04
1620 4.02
1680 3.98
P á g i n a 7 | 23
1740 3.95
1800 3.94
1860 3.93
Tabla 1. Resultados Experimentales @ 24 ºC
t (seg) Conductivida
d (μS)
0 9.05
30 8.94
60 8.83
90 8.72
120 8.61
150 8.46
180 8.19
210 8.11
240 8.03
270 7.96
300 7.9
330 7.83
360 7.78
390 7.72
420 7.65
450 7.58
480 7.52
510 7.44
540 7.38
570 7.32
600 7.27
630 7.24
680 7.21
690 7.18
720 7.15
750 7.11
780 7.08
810 7.07
840 7.06
870 7.03
900 7
930 6.97
960 6.95
990 6.93
1020 6.91
1050 6.9
P á g i n a 8 | 23
1080 6.86
1110 6.85
1140 6.83
1170 6.81
1200 6.8
1230 6.78
1260 6.77
1290 6.76
1320 6.76
1350 6.75
1380 6.74
1410 6.72
1440 6.71
1470 6.7
1500 6.68
1530 6.67
1560 6.66
1590 6.65
1620 6.64
1650 6.63
1680 6.62
1710 6.62
Tabla 2. Resultados Experimentales @ 15 ºC
Curva de calibración:
Conductivida CA / [OH]
d (μS)
11.6554 0.05
3.95 0.01
Análisis de resultados
El comportamiento de nuestra reacción fue el ideal, mediante el seguimiento de la reacción
por conductimetría, podemos ver que siguió la cinética que esperábamos. La reacción
efectuada a 24 °C fue la que mayor grado de avance efectuó, mientras que a 15 °C avanzo
con mayor rapidez ya que a los 28 min se mantuvo constante el valor de la conductancia para
ambos casos la concentración es similar.
Mediante el método integral fue como se determino el orden de reacción, en donde los valores
de correlación no eran precisamente cercanos a R2=1 sin embargo con esta dificultad
presentada se decidió en tomar valores los más altos para R 2, en ambos casos fue de orden
dos, posteriormente mediante el método diferencial pudimos comprobar el orden y aunque no
es tan viable este método debido una desviación muy significativa para la reacción de 25 °C
decidimos descartar mejor este método que fue precisamente el de las diferencias finitas y
aunque realizando un ajuste polinómico nos daba un orden aceptable, decidimos quedarnos
P á g i n a 9 | 23
con el método integral ya que sirvió para analizar mejor nuestros datos que tenemos sin
descartar algunos errores en la toma de las muestras. Posteriormente con estos resultados
trabajamos y determinamos las constantes de velocidad cinética para ambos experimentos, al
igual que poder obtener el grado de avance en la reacción efectuada.
Para la constante de equilibrio resulto comportarse como se esperaba ya que esta variable
tiene una gran dependencia a la temperatura, es de esperar de una reacción ligeramente
exotérmica, apreciablemente el valor de la constante de equilibrio aumenta conforme también
aumenta la temperatura.
La energía de activación necesaria para transformar los reactivos a productos, aplicando la
ecuación de Arrhenius, es de -27790.2171 J/mol, siendo el signo negativo por tratarse de una
reacción exotérmica; y concurriendo el valor de la constante de Arrhenius de 2.4295 x10-6
J/mol·ºC para 15 ºC y 4.905 x10-6 J/mol·ºC para 24 ºC.
Conclusión
Podemos decir que la temperatura a la que se sometió la reacción tiene mucha influencia en
los datos de concentración, las constantes de equilibrio y reacción, sobre todo, la rapidez con
la que los reactivos son transformados en productos. Nos dimos cuenta de que a una menor
temperatura la rapidez o velocidad de reacción aumenta significativamente, la conversión del
reactivo problema disminuye en gran medida y la constante de reacción aumenta. Por ello, se
propone que la reacción sea realizada a temperatura ambiente y a una temperatura que esté
por debajo de ésta. Si se buscan valores exactos como el orden y la constante de velocidad
cinética, es recomendable utilizar el método diferencial, propondría mejor un ajuste polinómico
en vez de diferencias finitas, aunque constaría de más datos de conductividad, lo cual se
pueden obtener midiendo el tiempo al inicio de la reacción en un intervalo de rango de 30 s y
después en uno de 1 min, ya que la conductividad disminuye mucho al comienzo y después
bajo muy poco antes del equilibrio de la reacción. Sin embargo, en un estudio con mayor
precisión si es recomendable seguir ese modelo, pero como ya teníamos conocimiento de la
reacción y su tendencia al orden dos, el método integral es el indicado en este experimento.
No se cuentan con muchos datos teóricos en la bibliografía consultada para poder comparar,
aun así, nuestros datos no se desvían mucho de los obtenidos en otras experimentaciones,
los valores se pueden modificar en el valor de la energía de activación ya que ahí depende el
tipo de reacción que es. Como referencias tenemos la de algunos artículos en internet de la
saponificación de esta reacción describiendo la idealidad.
Memoria de Calculo
NaOH
n- butanol
Ac. de butilo Ac. de sodio
Inicio CA0 CA0
Rxn XACA0 XACA0
P á g i n a 10 | 23
eq CA0(1-XA) CA0(1-XA) XACA0 XACA0
Si XA es el grado de avance, la fracción que no reacciona está dada por Ę= 1-XA
El valor de XACA0 corresponde a la diferencia de concentraciones entre el último valor de la
tabla de concentraciones molares y la concentración inicial, es decir el primer valor registrado
en la tabla de concentraciones molares. CA0(1-XA) corresponde al último valor en la tabla de
concentraciones molares:
[ Acetato de sodio ][ Butanol ]
Keq=
[ Acetato de butilo ][ Hidroxido de Sodio ]
2
( X A C A 0)
Keq= 2
[ C A 0 (1−X A )]
Sustituyendo a 24 °C
2
(0.030705−0.01718)
Keq= 2
=0.619766
[ 0.01718 ]
Sustituyendo a 15 °C
2
(0.04418−0.03232)
Keq= 2
=0.13465
[ 0.03232 ]
Mediante una curva de calibración se obtiene los valores de la concentración a cada tiempo.
Tabla 3. Datos de concentración a 24 °C
t (seg) CA / [OH]
0 0.03075843
30 0.02607143
60 0.02768259
90 0.0268526
120 0.02626672
180 0.02568085
240 0.0251438
300 0.02465557
360 0.02416734
420 0.02382558
480 0.02343499
540 0.02309323
P á g i n a 11 | 23
600 0.0228003
660 0.02245854
720 0.02221442
780 0.02197031
840 0.02182384
900 0.02162855
960 0.02133561
1020 0.02114032
1080 0.02094503
1140 0.02079856
1200 0.02060327
1260 0.0202615
1320 0.02021268
1360 0.02011504
1440 0.01987092
1500 0.01977328
1560 0.01972445
1620 0.01962681
1680 0.01943152
1740 0.01928505
1800 0.01923622
1860 0.0191874
1920 0.01899211
1980 0.01884564
2040 0.01874799
2100 0.01874799
2160 0.01860153
2220 0.0185527
P á g i n a 12 | 23
2280 0.01845506
2340 0.01835741
2400 0.01825977
2430 0.01821094
2460 0.01816212
2520 0.01806447
2580 0.01801565
2640 0.01791801
2700 0.01782036
2760 0.01782036
2820 0.01767389
2880 0.01762507
2940 0.01757625
3000 0.01752742
3060 0.0174786
3120 0.01738095
3180 0.01733213
3240 0.01728331
3300 0.01723448
3360 0.01718566
P á g i n a 13 | 23
Orden 0
0.035
0.03
0.025
f(x) = − 2.79874164951319E-06 x + 0.0250839149941088
R² = 0.849497692072105
0.02
[O]
0.015
0.01
0.005
0
0 500 1000 1500 2000 2500 3000 3500 4000
tiempo/s
Orden 1
-3
0 500 1000 1500 2000 2500 3000 3500 4000
-3.2
-3.4
Ln[OH]
-3.6
-4
-4.2
tiempo/s
Orden 2
70
60
f(x) = 0.00627815721700513 x + 39.4661423330355
R² = 0.939371348323189
50
40
1/[OH]
30
20
10
0
0 500 1000 1500 2000 2500 3000 3500 4000
tiempo/s
Debido a factor de correlación, el valor mas cercano a 1 se toma, por lo cual la reacción a 25
°C es de orden 2, con un valor de la constate de rapidez de: k = 0.0063 L mol-1 seg-1
Si lo convertimos a minutos el valor de la constante es de k = 0.3767 L mol-1 min-1
P á g i n a 14 | 23
Orden 2 (24 °C)
60
f(x) = 0.376689433020307 x + 39.4661423330355
R² = 0.939371348323189
55
50
1/[OH]
45
40
35
30
0 10 20 30 40 50 60
tiempo/min
Ahora mediante un método diferencial vamos a comprobar que el orden es de 2, haciendo uso
de las ecuaciones necesarias.
Dif. Finitas
-3
-9.2 -9 -8.8 -8.6 -8.4 -8.2 -8 -7.8 -7.6 -7.4
-3.2
-3.4
LnCA
-3.6
f(x) = 0.290527789677672 x − 1.437122635132
R² = 0.79346062607376
-3.8
-4
-4.2
LnRa
Como se puede ver en la gráfica, este método no es tan confiable, ya que no son valores
exactos por lo cual vamos a despreciarlo.
P á g i n a 15 | 23
Tabla 4. Datos de concentración a 15 °C
t (seg) CA / [OH]
0 0.04418473
30 0.04364768
60 0.04311062
90 0.04257357
120 0.04203652
150 0.04130418
180 0.03998596
210 0.03959537
240 0.03920479
270 0.03886303
300 0.03857009
330 0.03822833
360 0.03798422
390 0.03769128
420 0.03734952
450 0.03700776
480 0.03671482
510 0.03632424
540 0.0360313
570 0.03573837
600 0.03549425
630 0.03534778
680 0.03520131
690 0.03505484
720 0.03490838
750 0.03471308
P á g i n a 16 | 23
780 0.03456662
810 0.03451779
840 0.03446897
870 0.0343225
900 0.03417603
930 0.03402956
960 0.03393192
990 0.03383427
1020 0.03373663
1050 0.0336878
1080 0.03349251
1110 0.03344369
1140 0.03334604
1170 0.0332484
1200 0.03319957
1230 0.03310193
1260 0.03305311
1290 0.03300428
1320 0.03300428
1350 0.03295546
1380 0.03290664
1410 0.03280899
1440 0.03276017
1470 0.03271135
1500 0.0326137
1530 0.03256488
1560 0.03251605
1590 0.03246723
P á g i n a 17 | 23
1620 0.03241841
1650 0.03236958
1680 0.03232076
1710 0.03232076
Orden 0
0.046
0.044
0.042
0.036
0.034
0.032
0.03
0 200 400 600 800 1000 1200 1400 1600 1800
tiempo/s
Orden 1
-2.9
0 200 400 600 800 1000 1200 1400 1600 1800
-3
-3.1
[O]
-3.2
f(x) = − 0.000160763079223672 x − 3.20373207451562
R² = 0.866229077669127
-3.3
-3.4
-3.5
tiempo/s
Orden 2
32
f(x) = 0.00441850106141075 x + 24.5780302008588
R² = 0.891193085209453
30
28
26
[O]
24
22
20
0 200 400 600 800 1000 1200 1400 1600 1800
tiempo/s
P á g i n a 18 | 23
Orden 2 (15 °C)
35
20
[O] 15
10
0
0 5 10 15 20 25 30
tiempo/s
-3
-3.1
R² = 0.978068717579932 -3.2
-3.3
-3.4
-3.5
LnRa
m=n n 2.6315
lnK=b K 1711704.13 L mol^-1.6315 min^-1
En este caso nos da un orden 2.6315 por lo que el método diferencial se ajusta al orden que
nos proporciona el método integral y sin embargo tomaremos el valor de la constante de este
último.
Mediante la ecuación de Arrhenius se establece que se forma primero una especie de alta
energía (complejo activado), el complejo activado requiere para formarse una energía de
activación.
P á g i n a 19 | 23
Energía de activación para un proceso a dos temperaturas se expresa como:
ln
( ) (
K 1 Ea 1
K2
= −
1
R T2 T1 )
ln ( 0.3767
0.2651 8.314 J /mol ∙ K 297.15 K 288.15 K )
)=
Ea
( 1
−
1
J
Ea=−27790.2171
mol
De la ecuación de Arrhenius despejamos A para obtener el factor de frecuencia.
Factor de frecuencia (A)
Ea 1
lnk=lnA− ∙
RT T
− ( RTEa )
k=A∙e
Despejando el factor de frecuencia no da de la siguiente forma.
A 15 °C
( )
J
−27790.2171
mol
A=0.2651∙ exp −
(8.314 J /mol ∙ K )(288.15 K )
−6 L
A=2.4295 ×10
mol ∙ min
A 24 °C
( )
J
−27790.2171
mol
A=0.3767 ∙ exp
( 8.314 J /mol ∙ K )(297.15 K )
−6 L
A=4.905 ×10
mol ∙ min
Para ambos casos (24 °C y 15 ° C) el coeficiente de correlación fue más cercano a uno fue el
de orden dos, por este motivo la expresión que representa la velocidad de reacción es la
siguiente:
P á g i n a 20 | 23
d [ OH ] 2
=k [ OH ]
dt
Separando variables
Cat t
d [ OH ]
∫ 2
=k ∫ t
Ca0 [ OH ] 0
d [ OH ]
2
=¿
[ OH ]
Integrando
1 1
=kt +
[ OH ] t [ OH ] 0
Conversión experimental
Para 15 °C
0.0441847−0.03232
X A= =0.2685
0.0441847
P á g i n a 21 | 23
Para 24 °C
0.0307584−0.0171856
X A= =0.4412
0.0307584
Como podemos observar, la gráfica que tiene mayor tendencia lineal es la de orden dos, por
lo que nuestra saponificación de acetato de butilo es de segundo orden. Para determinar la
rapidez o velocidad de reacción, partimos de la ecuación generalizada para reactor por lotes
de orden mayor a cero:
2
Ra=k C A
Ra= ( dCdt )
A
Anexos
Graficando el inverso de rapidez contra el grado de avance de la reacción, obtenemos la
siguiente gráfica a cada temperatura:
5000
4000
3000
2000
1000
0
0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45
XA % conv. Teórica
P á g i n a 22 | 23
880 XA vs 1/Ra (15 °C)
860
840
820
800
1/Ra
780
760
740
720
0 0.05 0.1 0.15 0.2 0.25 0.3
XA
% conv. Teórica
Bibliografía
Levenspiel, O., (2010), Ingeniería de las reacciones químicas, México DF, México, 3ra ed., Limusa.
Recuperado: https://www.ugr.es/~focana/dfar/aplica/cinetiHidroli/cinetiHidroliAcMet.pdf
Perry, R.H. (1992). Manual del Ingeniero Químico. México. Mc Graw Hill
Universidad Autónoma de México. (Enero 2015). Datos Tecnicos de Acetato de Etilo. Obtenido de
http://www.quimica.unam.mx/IMG/pdf/8acetatodeetilo.pdf
P á g i n a 23 | 23