Reporte Equilibrio Quimico y Reactor Por Lotes Equipo - 2

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 23

Universidad Nacional Autónoma de México

Facultad de Estudios Superiores Cuautitlán


Campo 1

Ingeniería Química

Laboratorio Experimental Multidisciplinario VI


Profesora: Paula Álvarez Fernández
Irving Cerino Carrasco
Reporte Equilibrio Químico Y Reactor Por Lotes
Grupo: 2853
Alumnos:
Cortes Cabañas Jair Raúl
Piñón González José Adrián
Sánchez Escalera Atzel Rafael
Trejo Galván José Raúl
Semestre: 2023-II

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

Equilibrio químico es la denominación que se hace a cualquier reacción reversible cuando se


observa que las cantidades relativas de dos o más sustancias permanecen constantes, es
decir, el equilibrio químico se da cuando la concentración de las especies participantes no
cambia, de igual manera, en estado de equilibrio no se observan cambios físicos a medida
que transcurre el tiempo; siempre es necesario que exista una reacción química para que
exista un equilibrio químico, sin reacción no sería posible.

En un sistema en equilibrio, se dice que se encuentra desplazado hacia la derecha si hay


mayor cantidad de productos presentes que de reactivos y se encuentra desplazado a la
izquierda cuando exista mayor cantidad de reactivos que de productos

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.

En este informe experimental se estudió la reacción de saponificación de acetato de etilo con


el fin de conocer los parámetros cinéticos.

El reactor discontinuo es un recipiente en el que están contenidas las sustancias mientras


reaccionan, lo que se ha de determinar es el avance de la reacción para diversos tiempos.
Algunas de las aplicaciones del reactor por lotes son en la producción farmacéutica y en
procesos de fermentación, producción a pequeña escala, para reacciones en fase líquida que
requieran largos tiempos de reacción, se aplica en el control de contaminación, plantas pilotos
y en calorimetría de reacciones.

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:

 Kc es una constante para una reacción específica a una temperatura específica. Si


cambias la temperatura de la reacción, entonces Kc también cambia.
 Los sólidos y líquidos puros, incluyendo disolventes, no se consideran para la expresión
de equilibrio.
 La reacción debe estar balanceada con los coeficientes escritos como el mínimo valor
entero posible para poder obtener el valor correcto 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:

Esta ecuación se conoce como ecuación de Arrhenius, donde k es la constante de velocidad,


Ea es la energía de activación, R es la constante de los gases, T la temperatura absoluta y A
es el factor de frecuencia.
El reactor por lotes o intermitente es un reactor integral, por lo que es el equipo preferido para
la obtención de datos cinéticos en sistemas homogéneos. En este reactor no hay intercambio
de materia en el exterior y se introduce una mezcla reaccionante a las condiciones de presión
y temperatura tales que la velocidad de reacción cambia proporcionalmente con la
concentración en función del tiempo. Una vez alcanzado determinado el grado de avance de
la misma, el reactor se vacía.
Es un sistema cerrado, opera isotérmicamente, estado no estacionario y a volumen constante,
por lo tanto, la masa total del mismo es constante y su tiempo de residencia de todos los
elementos del fluido es el mismo. Su variable más importante es la conversión y con esta
podemos obtener el tiempo de residencia y el volumen del equipo.

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

f(x) = − 0.000131611367643228 x − 3.68202076281013


-3.8 R² = 0.899972462549347

-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.

Usando el método de las pendientes obtenemos lo siguiente:

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.04 f(x) = − 5.88092879784282E-06 x + 0.0405606691042666


R² = 0.838721824615172
0.038
[O]

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

Obtenemos el orden dos mediante la correlación, que es un valor más cercano a 1.


Orden dos; K = 0.0044 L mol-1 seg-1
Si lo convertimos a minutos el valor de la constante es de k = 0.2651 L mol-1 min-1

P á g i n a 18 | 23
Orden 2 (15 °C)
35

30 f(x) = 0.265110063684645 x + 24.5780302008588


R² = 0.891193085209453
25

20

[O] 15

10

0
0 5 10 15 20 25 30
tiempo/s

Mediante el método diferencial vamos a comprobar que el orden es de 2

Método diferencial (15 °C)


-2.9
-6.78 -6.76 -6.74 -6.72 -6.7 -6.68 -6.66 -6.64 -6.62

-3

-3.1

f(x) = 2.63154890500948 x + 14.3525449766072


LnCA

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.

Tabla 5. Valores de la constante cinética a cada temperatura


T (°C) K (L mol-1 min-1)
24 0.3767
15 0.2651

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

Para conocer la conversión teórica se hace uso de la siguiente ecuación:


tk C A 0
X A=
1+tk C A 0
Para 15 °C
L mol
(28.5 min)(0.2651 )(0.04418 )
mol ∙ min L
X A= =0.25028
L mol
1+(28.5 min)(0.2651 )(0.04418 )
mol ∙ min L
Para 24 °C
L mol
(56 min)(0.3767 )(0.03075 )
mol ∙ min L
X A= =0.3935
L mol
1+(56 min)(0.3767 )(0.03075 )
mol ∙ min L

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:

9000 XA vs 1/Ra (25 °C)


8000
7000
6000
1/Ra

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

Como se describe, la conversión teórica es menor a la experimental.

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

Fogler,H.S.(2008).Elementos de ingenierías de las reacciones químicas. 4ed.Pearson Educación.

Perry, R.H. (1992). Manual del Ingeniero Químico. México. Mc Graw Hill

Gaido, A. Ruiz, A. &Bonami, F. (2002). Saponificación del Acetato de Etilo. Página 23

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

También podría gustarte