Manual Lab. Tef

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 59

INSTITUTO POLITÉCNICO NACIONAL

ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS

DEPARTAMENTO DE INGENIERIA QUIMICA INDUSTRIAL

LABORATORIO DE
TERMODINÁMICA DEL EQUILIBRIO DE FASES

PRÁCTICA 1

EQUILIBRIO DE FASES PARA


SUSTANCIAS PURAS

EQUILIBRIO DE FASES PARA SUSTANCIAS PURAS

OBJETIVOS

• Determinar las condiciones de equilibrio para sustancias puras.


• Calcular los calores de vaporización y sublimación de sustancias puras.
• Estimar el punto triple de una sustancia pura.
• Verificar la validez de la ecuación de Claussius - Clapeyron.
• Determinar la entalpía de vaporización del agua y compararlo con el valor bibliogarfico.

1
Práctica 1: Equilibrio de fases para sustancias puras

HIPÓTESIS
“La presión de vapor de un líquido es constante a una temperatura dada, pero aumenta si aumenta la
temperatura”.

EQUILIBRIO DE FASES PARA SUSTANCIAS PURAS

¿Como se alcanza el equilibrio físico entre fases?

Se alcanza el equilibrio material cuando los potenciales químicos de la sustancia son iguales en las

dos fases:

En un sistema cerrado en equilibrio termodinámico, el potencial químico de un componente


dado es el mismo en todas las fases en las que el componente está presente.

Vamos a centrarnos en el estudio del equilibrio de fases en sistemas formados por un solo
componente (sustancia pura). Para especificar el estado termodinámico de un sistema formado por
una sustancia pura, el número variables intensivas independientes que se deben conocer (grados
de libertad) es:

• Si hay presente una fase, L = 1componente - 1fase + 2 = 2 variables, es necesario especificar


por ejemplo la P y la T
• Si hay presente dos fases, L = 1componente - 2fases + 2 = 1 variable, es necesario
especificar sólo P o T
• Si hay presente tres fases, L = 1componente - 3fases + 2 = 0 variables

Por tanto podemos representar cualquier estado de equilibrio del sistema formado por una sustancia
pura mediante un punto en un diagrama bidimensional de presión-temperatura. Este diagrama se
denomina diagrama de fases. Un ejemplo del diagrama de fases de una sustancia pura se muestra
en la figura.

2
Práctica 1: Equilibrio de fases para sustancias puras

En el diagrama de la figura las líneas AB, BD y BC corresponden a valores (P,T) en las que coexisten
dos fases:

• En AB coexisten en equilibrio sólido y gas. La curva AB es la curva de presión de vapor


del sólido
• En BD coexisten en equilibrio sólido y líquido.
• En BC coexisten en equilibrio líquido y gas.

El punto B marca los valores de P y T en los que coexisten tres fases, sólido, líquido y gas, y se
denomina Punto Triple. Este punto, que indica la temperatura mínima a la que el líquido puede
existir, es característico de cada sustancia, y puede emplearse como referencia para calibrar
termómetros.
El punto C indica el valor máximo (PC, TC) en el que pueden coexistir en equilibrio dos fases, y se
denomina Punto Crítico. Representa la temperatura máxima a la cual se puede licuar el gas
simplemente aumentando la presión. Fluidos con T y P mayores que T C y PC se denominan fluidos
supercríticos.

Tabla 1.1 correlaciones para la estimación de la presión de vapor de un componente puro

Efecto de la presión y de la temperatura. Ecuación de Clapeyron.

La ecuación de Clapeyron permite calcular la pendiente de una línea de equilibrio entre dos fases en
el diagrama de fases P-T de un sistema de un componente.

Consideremos un punto cualquiera sobre una línea de equilibrio entre dos las fases, que llamaremos

α y β. La condición para que exista equilibrio de fases es que: , pero para una sustancia

pura por tanto en un punto sobre la curva de equilibrio de dos fases ,y


cualquier variación infinitesimal que suponga un desplazamiento sobre la curva de equilibrio implica

que . O lo que es lo mismo:

3
Práctica 1: Equilibrio de fases para sustancias puras

y reagrupando términos . Por otra parte si se considera que en un cambio de fase

reversible a T y P constantes se tiene que:

Ecuación de Clapeyron

Consideraciones sobre la ecuación de Clapeyron:

• En un cambio de fase líquido-vapor, tanto ΔH como ΔV son positivos, por tanto la pendiente
de la línea de equilibrio líquido-vapor es positiva. Lo mismo sucede con la línea sólido-vapor.
• En un cambio de fase sólido-líquido, ΔH es positivo y en general ΔV también, por lo tanto la
pendiente de esta línea también será positiva. Existen sin embargo algunas excepciones
como el H2O, Ga o Bi debido a una disminución de volumen que sufren estos componentes
al fundirse, en estos casos la pendiente de la línea de equilibrio sólido-líquido será negativa.
• En el cambio de fase sólido-líquido ΔV es mucho menor que en los cambios de fase sólido-
gas o líquido-gas. Por esta razón la pendiente en el primer caso es mucho mayor que en los
últimos.

Aplicación de la ecuación de Clapeyron a distintos cambios de fase

• Equilibrio líquido-vapor y sólido-vapor

En estos dos casos el V molar del gas es mucho mayor que el del líquido o que el del sólido
por lo que puede hacerse la aproximación

Si además se hace la suposición de que el gas se comporta como gas ideal, la ecuación de

Clapeyron se transforma en:

Esta ecuación se suele expresar como llamada ecuación de Clausius Clapeyron.

Si el rango de temperatura analizado es pequeño, se puede suponer que ΔH es constante a lo largo

de la línea de equilibrio, y por tanto:

ecuación de Clausius-Clapeyron integrada

La ecuación de Clasius-Clapeyron sólo es aplicable para obtener la presión de vapor de un sólido


o un líquido a una cierta temperatura, conocido otro punto de equilibrio entre las fases. Esto es
así porque para llegar a esta expresión desde la ec. de Clapeyron se hace la aproximación de

4
Práctica 1: Equilibrio de fases para sustancias puras

despreciar el volumen molar del sólido o del líquido frente al del gas, que además se supone
de comportamiento ideal.

Rudolf Clausius sugirió que la integración de la ecuación de CLAPEYRON permitía la determinación


de la entalpía para los cambios de fase a partir de datos experimentales de presión de vapor a
diferentes temperaturas.

− H vap 1
LnP =  +C
R T 

Esta ecuación es semejante a la ecuación de una recta: y = mx + b


− H
Donde: y = LnP , m= , x= 1 , b=C
R T

Si graficamos LnP contra 1/T y calculamos el valor de la pendiente podemos calcular el valor de
ΔHvap.

DESARROLLO EXPERIMENTAL

Material y equipo
Llave
• Refrigerante de paso
• Termómetro de 0°C – 100°C
• Hervidor (matráz de fondo
plano)
• Parrilla de calentamiento con
agitación magnética condensador

• Soporte universal agua de


• Reóstato recirculación

• Pinzas de nuez
• Pinzas de tres dedos termómetro
• Regla de 60 cm con
aproximación a mm Trampa de vacío
• Manómetro de mercurio Manómetro matráz
• Agua destilada de mercurio concavidad
mezcla

Procedimiento experimental

1. Montar el equipo como se muestra en el esquema.


2. Revisar perfectamente las conexiones para evitar fugas o entradas de aire.
3. Poner el líquido en el hervidor.
4. Arrancar la bomba de vacío hasta su máxima capacidad de succión y observar la diferencia
de niveles creada en las columnas del manómetro.
5. Calentar el líquido regulando el termostato del equipo de tal forma que no se caliente
excesivamente. (con agitación constante).
6. Esperar a que el líquido llegue a ebullición y tanto la presión como la temperatura
permanezcan constantes. (sólo hasta este momento comenzar a tomar los datos de
temperatura y presión alcanzados).

5
Práctica 1: Equilibrio de fases para sustancias puras

7. Cambiar la presión del sistema mediante la inyección de aire, utilizando la válvula


reguladora.
8. Repetir los pasos 6 y 7 hasta alcanzar la presión atmosférica. (10 lecturas)

TABLA DE DATOS EXPERIMENTALES

Temperatura Presión
Experiment
ln P 1/T
o Manométrica Absoluta(P)
T (°C) T (K)
(Δh)mm Hg =Patm-Δh

10

11

6
Práctica 1: Equilibrio de fases para sustancias puras

Cálculos

1.- Grafique los datos experimentales Psat vs T (temperatura) para obtener la curva de presión de
vapor experimental.

7
Práctica 1: Equilibrio de fases para sustancias puras

Curva de presión de vapor del agua

500

450

400
Presión (mm Hg)

350

300

250

200

150

100
330 335 340 345 350 355 360 365

Temperatura (K)

2.- Con los datos experimentales determinar el calor de vaporización (ΔHvap) ya sea por mínimos
cuadrados o calculando la pendiente de la recta, graficando lnP vs 1/T por regresión lineal, aplicando
la ecuación de Clasius-Clapeyron.

8
Práctica 1: Equilibrio de fases para sustancias puras

− H vap  1  Ec. que define la curva de


LnP =  +C ← vaporización
R T 
y= m x +b
Esta ecuación es semejante a la ecuación de una recta:
Donde:
− H vap
y = LnP x=
1
 m= = b=C=
T R

Linealizacion de la ecuacion de Clasius-Clapeyron


6.5

6.3

6.1

5.9

5.7
ln P

5.5

5.3

5.1

4.9

4.7

4.5
0.0026 0.0027 0.0028 0.0029 0.003 0.0031
1/T K

H vap = − m * R  H vap =

Donde R=8.314 J/mol K

3.- Con la ecuación que define a la curva de vaporización suponer valores de temperatura superiores
al punto triple hasta el primer valor de temperatura experimental. Graficar los valores supuestos.

9
Práctica 1: Equilibrio de fases para sustancias puras

1 Ec. que define la curva de


LnP =  + ← vaporización
8.314 T 

Curva de presión de vapor con Temperaturas Supuestas

200

180

160

140
Presión (mm Hg)

120

100

80

60

40

20

0
270 280 290 300 310 320 330

Temperatura (K)

Tsup. (K)

Pcalc. (mmHg)
4.- Calcular el calor latente de sublimación molar ΔHsub, utilizando el calor de fusión bibliográfico y la
ΔHvap experimental ya calculada.

10
Práctica 1: Equilibrio de fases para sustancias puras

𝛥𝐻 𝑠𝑢𝑏 = 𝛥𝐻 𝑓𝑢𝑠 + 𝛥𝐻𝑣𝑎𝑝

H sub =

5.- De la ecuación de Clausius Clapeyron, encontrar la ecuación que defina a la curva de


sublimación. Es recomendable calcular la constante de integración “c” de la ecuación de Clausius
Clapeyron con los datos de temperatura y presión del punto triple.

− H sub 1 ← Ec. que define la curva de


LnP =  +C
R T  sublimación
Datos necesarios para cálculos (teóricos)
Presión en el Punto Triple del agua
PPT =

Temperatura en el Punto Triple del agua


TPT =

Constante universal de los gases ideales


R=

 fusión 
 
Agua  sub lim ación 
 
 vaporización 

Considerando los datos de bibliografía, determinar la constante de integración “c” para la


expresión que define la curva de sublimación.

− H sub  1 
LnP =  +C
R T 
− H sub  1 
C = ln P −  
R T 

−  
C = ln ( )− 1
 
8.314  

C=

6.- Calcular las presiones a diferentes temperaturas supuestas por debajo de la temperatura del
punto triple y grafique los datos de Psat vs T (temperatura) para obtener la curva de sublimación.

11
Práctica 1: Equilibrio de fases para sustancias puras

− 1 Ec. que define la curva de


LnP =  + ←
sublimación
8.314 T 

Curva de Sublimación supuesta

4.5

3.5
Presión (mm Hg)

2.5

1.5

0.5

0
220 230 240 250 260 270

Temperatura (K)

Tsup. (K)

Pcalc. (mmHg)

7.- Con los gráficos obtenidos representar el diagrama de fases de una sustancia pura (agua).

12
Práctica 1: Equilibrio de fases para sustancias puras

Diagrama de Fases del agua

500

450

400

350
Presión (mm Hg)

300

250

200

150

100

50

0
220 230 240 250 260 270 280 290 300 310 320 330 340 350 360

Temperatura (K)

Curva de sublimación sup. Curva vaporización sup. Curva vaporización exp.


T(K)supuesta P (mm Hg) Tsup(K) P (mm Hg) Tsup(K) P (mm Hg)

Conceptos

13
Práctica 1: Equilibrio de fases para sustancias puras

Entregar de tarea las definiciones de lo siguiente ( derecho a hacer la práctica )

1. Equilibrio termodinámico:

2. Sustancia pura:

3. Presión de vapor:

4. Punto de ebullición:

5. Sistema y alrededores:

6. Variable dependiente e independiente:

7. Presión atmosférica:

8. Presión manométrica:

9. Presión absoluta:

10. Diagrama de fases:

11. Sublimación

12. Vaporización

13. Fusión

Reporte.

14
Práctica 1: Equilibrio de fases para sustancias puras

En equipo: 50%
1. Gráfica de P (mmHg) vs T (°C)
2. Gráfica de LnPabs vs 1/T (°K)
3. Cálculos de la pendiente y ΔHvap
4. Calcular ΔHsub = ΔHvap + ΔHfus [=] J/mol
5. Calcular datos para construir diagrama de fases
6. Gráfica donde se muestre el diagrama de fases
7. Observaciones

Individual: 50%
Conclusiones

Conclusiones.

Bibliografía.

15
INSTITUTO POLITÉCNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS

DEPARTAMENTO DE INGENIERIA QUIMICA INDUSTRIAL

LABORATORIO DE
TERMODINÁMICA DEL EQUILIBRIO DE FASES

PRÁCTICA 2

PROPIEDADES MOLARES
PARCIALES
OBJETIVOS

• Determinar el volumen real de mezclas de soluciones.


• Aplicar el principio de Arquímedes para determinar la densidad de cada una de las soluciones.
• Determinar para cada solución, el cambio de volumen de mezclado.
• Ajustar los datos experimentales a una ecuación (Expansión de Redlich-Kister ).
• Determinar los volúmenes molares parciales de agua y alcohol en distintas soluciones
binarias, determinando las densidades de las mismas.

1
Práctica 2: Propiedades molares parciales

HIPÓTESIS
Si tomamos una propiedad extensiva cualquiera, por ejemplo el volumen de una mezcla binaria, y
suponemos que “al mezclar dos volúmenes cualquiera de dos sustancias puras diferentes, el
volumen resultante no es la suma de los volúmenes iniciales, sino que puede presentar diferencias
significativas”

PROPIEDADES MOLARES PARCIALES

Nos permite establecer las diferencias que presentan en su comportamiento fisicoquímico las
sustancias cuando están en estado puro y cuando forman una mezcla.

A la temperatura T y la Presión P constantes, para una sustancia dada, sus propiedades: volumen,
energía interna, entropía, energía de Helmholtz, entalpía, energía de Gibbs, densidad, curva de
presión de vapor, etc. adquieren valores fijos característicos.

A la misma T y P, para cualquier otra sustancia, el valor numérico de sus propiedades serán
diferentes a los de la anterior, pudiera suceder que en alguna hubiera coincidencia, pero no puede
ser posible que todas coincidan pues entonces se trataría de la misma sustancia.

Cuando dos o más sustancias se mezclan, a T y P constantes, para formar una solución, las
propiedades de la solución son función de las propiedades de los componentes que la forman. La
funcionalidad que guardan las propiedades de la solución respecto a las propiedades de los
componentes es compleja y sólo se determina por medio del experimento.

Soluciones ideales y no ideales.

En raras ocasiones, las propiedades de las soluciones resultan ser aditivas respecto a las
propiedades de sus componentes:
N
M =  xi M i
1
N
M t =  ni M i
1
Donde:

M y Mt = Propiedad molar y propiedad total de la solución.


xi y ni = Frac. mol y número de moles de cada componente.
Mi = Propiedad molar de cada componente en estado puro.
M = Propiedad extensiva: Volumen, Energía interna, Entalpía, Cp y otras, excepto Entropía y las
derivadas de ella.
Las soluciones que obedecen estas expresiones se les llaman SOLUCIONES IDEALES.

Las SOLUCIONES NO IDEALES obedecen ecuaciones análogas:

N
M =  xi M i
1

N
M t =  ni M i
1

Laboratorio de Termodinámica del Equilibrio de Fases 2


Práctica 2: Propiedades molares parciales

donde M i se llama: PROPIEDAD MOLAR PARCIAL y es el valor de la propiedad que tiene un


componente EN SOLUCION, y que difiere del valor de la propiedad como componente puro debido
a las interacciones con las demás sustancias de la solución.

En particular para un sistema de dos componentes:

Una propiedad molar parcial se obtiene mediante la expresión:

(n M )
Mi =
ni T , P , n j i
Aplicaciones:

Las propiedades molares parciales juegan un papel importantísimo en la termodinámica, ejemplos:

Volúmen molar parcial: Permite hacer cálculos exactos en balances de masa.


Entalpía molar parcial: Permite hacer cálculos exactos de balance de energía.
Energía libre de Gibbs molar parcial o Potencial químico: Es el fundamento teórico que permite la
deducción de las ecuaciones básicas para resolver todos los problemas relativos a equilibrio entre
fases y en reacciones químicas.

Estos tres tipos de problemas son parte importante de la razón de ser del Ingeniero Químico.

DESARROLLO EXPERIMENTAL

El método utilizado en la práctica se denomina “MÉTODO DINÁMICO DE CÁLCULO DE


DENSIDADES POR MEDIO DEL FLOTADOR” el cual se explica a continuación:

Material y equipo

• Tubos de Nessler cuidadosamente seleccionados


• Balanza analítica.
• Flotador, Probeta de 40 ml
• Soporte de madera
• Gancho de alambre
• 2 vasos de precipitado
• 2 bureta de 50 ml cada una
• Agua destilada
• Metanol

Procedimiento experimental

1. Montar el equipo como se muestra en el esquema.


2. Revisar que los tubos estén perfectamente limpios y secos
3. Preparar un juego de soluciones Metanol-Agua cuya concentración varíe de 0 a 1 con
intervalos de 0.1 en fracción mol con respecto al etanol. Los volúmenes de cada solución
pueden calcularse con la ayuda de las siguientes fórmulas, tomando en cuenta que los
subíndices 1 son relacionadas al metanol y el 2 serán de agua

Laboratorio de Termodinámica del Equilibrio de Fases 3


Práctica 2: Propiedades molares parciales

x1 PM 1 VT x 2 PM 2 VT
1 2
V1 = V1 =
x1 PM 1 x 2 PM 2 x1 PM 1 x 2 PM 2
+ +
1 2 1 2

4. Una vez calculados los volúmenes, agregar al tubo de Nessler primero el agua y después
muy lentamente el metanol.
Bureta
5. Después de agregar el metanol verificar como
quedan los volúmenes y que llama la atención del
experimento en ese paso.
Muestras
6. Enseguida agitar cada tubo para mezclar
perfectamente los componentes puros y observar los
cambios que ocurren en los volúmenes (tomar nota).

7. Enfriar las soluciones

8. Apoyándose en el principio de Arquímedes


determinar la densidad real de los componentes
puros y de las soluciones mediante los siguientes
pasos:

a. Pesar el flotador en el aire WF: ponga la


balanza en ceros, cuelgue el flotador en el gancho y tome el peso en gramos
aproximado hasta la segunda cifra decimal.
b. Pesar el flotador sumergido en el agua W2 y en el metanol W1: Ponga en la probeta
aproximadamente 35 ml del líquido. Coloque la probeta en el soporte de madera,
cuelgue el flotador en el gancho de alambre y sumerja el flotador en la probeta de
manera que éste no toque sus paredes ni el fondo de la misma, espere que se
estabilice el peso en la balanza y enseguida tome la lectura del peso del flotador en
cada uno de los líquidos.
c. Pesar el flotador sumergido en cada solución. Proceda como en el paso anterior.

Datos necesarios para cálculos


Hilo

Peso
Sustancia Densidad (  )
Soporte
Molecular metálico
(PM) Flotador
Metanol (1)

Agua (2) probeta

Soporte
Volumen Total. _______________ ml de madera

Balanza
analítica

Laboratorio de Termodinámica del Equilibrio de Fases 4


Práctica 2: Propiedades molares parciales

Cálculos
1.- Con los datos obtenidos calcular:

o El empuje del flotador en el agua.

• E F = W F − W1
o E volumen del flotador
EF
• VF =
 agua
o El empuje del flotador en cada solución

• E Fi = WF − Wi

o La densidad de cada solución


E Fi
• i =
VF

o El peso molecular de cada solución

• PM i = x1 PM1 + x2 PM 2

o El volumen molar real para cada solución

PM i
• Vi real =
i

o El volumen ideal de cada solución

x1 PM 1 x 2 PM 2
• Vi ideal = +
1 2

o El cambio de volumen de mezclado de cada solución

• V1 = Vi real − Vi ideal

Laboratorio de Termodinámica del Equilibrio de Fases 5


Práctica 2: Propiedades molares parciales

Metanol Agua
X1 X2 V1 (ml) V2 (ml) Wi(gr) EFi ρi (g/ml) PMi (gr) Vreal (ml) Videal (ml) ΔVexp (ml)

0 1

0.1 0.9

0.2 0.8

0.3 0.7

0.4 0.6

0.5 0.5

0.6 0.4

0.7 0.3

0.8 0.2

0.9 0.1

1.0 0

Peso del flotador en el aire. WF =

Peso del flotador en el agua. W1 =

Peso del flotador en metanol. W2 =

Laboratorio de Termodinámica del Equilibrio de Fases 7


Práctica 2: Propiedades molares parciales

2.- Graficar.

Composicion vs cambio de volumen experimental


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
0

-0.2

-0.4
ΔVexp

-0.6

-0.8

-1

-1.2

-1.4

X1
3.- Ajuste los datos experimentales a la expansión de Redlich-Kister por el método gráfico.


V = x1 x2 A + B( x1 − x2 ) + C ( x1 − x2 ) 2 
 
V exceso = V = x1 x2 A + B( x1 − x2 ) + C ( x1 − x2 ) 2 = V real − V ideal

a. Realizar un cambio de variable a una ecuación polinómica de segundo orden, para poder
determinar A,B y C mediante este cambio de variable.

Si: 
V = x1 x2 A + B( x1 − x2 ) + C ( x1 − x2 ) 2 
V

= A + B( x1 − x2 ) + C ( x1 − x2 ) 2  Donde:
x1 x2 V
y=
y =  A+ B x + C x 2
 x1 x 2
x = ( x1 − x 2 )

Laboratorio de Termodinámica del Equilibrio de Fases 8


Práctica 2: Propiedades molares parciales

Ajuste por cambio de variable


(x1 - x2)
ΔV/(x1*x2)

b. Graficar X y Y para determinar una ecuación polinómica y obtener las constantes A,B y C.

Ajuste por cambio de variable

-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1


0

-0.5

-1

-1.5
ΔVexp

-2

-2.5

-3

-3.5

-4

-4.5

-5
X1
De gráfico:

y =  + x + x2 
‫؞‬ A=
B=
C=

Laboratorio de Termodinámica del Equilibrio de Fases 9


Práctica 2: Propiedades molares parciales

c. Recalcular ΔVajustada con ecuación de Redlich Kister para cada concentración.


V = x1 x2 A + B( x1 − x2 ) + C ( x1 − x2 ) 2 
X1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
ΔVajus

4.- Graficar a ΔVajustada y compararla con la experimental.

Composicion vs cambio de volumen experimental


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
0

-0.2

-0.4
ΔVajustada

-0.6

-0.8

-1

-1.2

-1.4
X1

5.- Calcular la propiedad molar parcial.

Si
V1 =V10 + V 1 y V2 =V20 + V 2

 nV1   nV2 
V 1 =   V 2 =  
 n1   n2 
T , P , n2 T , P , n1

Laboratorio de Termodinámica del Equilibrio de Fases 10


Práctica 2: Propiedades molares parciales

d. Poner en función del número de moles a la ecuación Riedlich Kister.

Para ΔV1 :

Entonces:

V1 = x1 x2 A + B( x1 − x2 ) + C ( x1 − x2 ) 2 
x2 = 1 − x1

Para ΔV2 :

Entonces:

V2 = x1 x2 A + B( x1 − x2 ) + C ( x1 − x2 ) 2 
x1 = 1 − x2

e. Derivar parcialmente a ΔV1 y ΔV2 para obtener la propiedad molar parcial de las
mismas V 1 y V 2

Laboratorio de Termodinámica del Equilibrio de Fases 11


Práctica 2: Propiedades molares parciales

f. Calcular la propiedad molar parcial V1 y V2 para cada concentración.


0 0
Se sabe que V y V son los componentes puros en la mezcla.
1 2

5.- Calcular Vreal para cada concentración.

real
Vajustada = x1V1 + x 2 V 2

6.- Calcular Videal para cada concentración.

real
Vajustada = x1V10 + x 2V20

Tabla de resultados
V10 = V20 =
Propiedades Molares Parciales
V 1 V 2
real ideal
X1 V1 V2 Vajustada Vajustada
0.0

0.1
0.2

0.3

0.4
0.5

0.6

0.7
0.8

0.9

1.0

Laboratorio de Termodinámica del Equilibrio de Fases 12


Práctica 2: Propiedades molares parciales

ideal real
7.- Graficar Tabla de resultados de Propiedades Molares Parciales. X1 vs Vajustada , X1 vs Vajustada
X1 vs V 1 y X2 vs V 2

Propiedades Molares Parciales.

45

40

35

30

25

20

15

10

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

X1

Laboratorio de Termodinámica del Equilibrio de Fases 13


Práctica 2: Propiedades molares parciales

Conceptos
Entregar de tarea las definiciones de lo siguiente ( para derecho a hacer la práctica )

1. Propiedades molares:

2. Propiedades molares parciales:

3. Mezclas:

4. Principio de Arquímides:

5. Densidad:

6. Volumen real:

7. Volumen ideal:

8. Peso molecular:

9. Cambio de volumen de mezclado:

10. Empuje

11. Picnometría

12. Líquido puro

13. Tubo Nessler

Laboratorio de Termodinámica del Equilibrio de Fases 14


Práctica 2: Propiedades molares parciales

Reporte.

En equipo: 50%
g. Tabla de resultados
h. Gráfica Videal, Vreal vs X1 (Metanol (1) – Agua (2))
i. Ajustes necesarios con la ecuación de Riedlich – Kister
j. Gráficos V1, V2, Vrealajustado vs X1
k. Observaciones

Individual: 50%
Conclusiones

Conclusiones.

Bibliografía.

Laboratorio de Termodinámica del Equilibrio de Fases 15


INSTITUTO POLITÉCNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS

DEPARTAMENTO DE INGENIERIA QUIMICA INDUSTRIAL

LABORATORIO DE
TERMODINÁMICA DEL EQUILIBRIO DE FASES

PRÁCTICAS 3

EQUILIBRIO LÍQUIDO-VAPOR
EQUILIBRIO IDEAL - “LEY DE RAOULT”

1
Práctica 3 y 4: Equilibrio Ideal y No Ideal

EQUILIBRIO IDEAL - LEY DE RAOULT


OBJETIVOS

• Observar el comportamiento de dos sistemas binarios.


• Construir los diagramas de temperatura y composición del líquido y del vapor.
• Comparar con los datos de temperatura y composición obtenidas con la Ley de Raoult.
• Determine la temperatura de burbuja y composición de la fase gaseosa teóricas de cada
solución, usando el algoritmo apropiado que se deriva de la ley de Raoult.

PROPÓSITO
- Determinar las curvas de equilibrio líquido-vapor para mezclas binarias de líquidos
totalmente miscibles.
- Determinar si las mezclas obedecen a la ley de Raoult.
- Realización del test de consistencia termodinámica de los datos obtenidos.

HIPÓTESIS
A una presión constante las mezclas de líquidos miscibles de diferente composición hervirán
a diferentes temperaturas y su vapor será más rico en el componente más volátil.

EQUILIBRIO LÍQUIDO VAPOR

I. INTRODUCCION

Una fase de acuerdo a la definición clásica, se define como materia uniforme en todas sus
partes, tanto en composición química como física, ejemplos de ésta son: agua líquida, hielo,
vapor de agua, solución de sacarosa en agua, la mezcla de gases en una muestra de aire etc.
Dos o más fases se presentan cuando hay una superficie macroscópica de separación entre
las fases homogéneas, ejemplos son: un sistema formado por tres fases (hielo-agua-aire),
sistemas formados por dos fases (líquidos inmiscibles) y mezclas de sólidos o formas
cristalinas diferentes del mismo sólido. El número de fases identificables depende de cómo
se observe la muestra, por ejemplo, una solución de DNA es homogénea a simple vista, y a
nivel microscópico se observa no homogénea.
Existen diferentes tipos de equilibrios entre fases: sólido–líquido (S-L), líquido-líquido (L-
L), líquido-vapor (L-V), etc., y también, pueden existir equilibrios de un solo componente,
de dos componentes llamado binario y de tres componentes denominado ternario El estado
de una muestra de materia se define por dos tipos de variables: intensivas y extensivas. Las
variables intensivas son aquellas que no dependen de la cantidad de materia presente como
la temperatura, presión y concentración; y las variables extensivas son aquellas que si
presentan dependencia de la cantidad de materia presente como por ejemplo la masa y
volumen.

Laboratorio de Termodinámica del Equilibrio de Fases 2


Práctica 3 y 4: Equilibrio Ideal y No Ideal

En un sistema de una o más fases, generalmente se requiere conocer al menos una variable
extensiva por fase, para determinar la masa de cada una de ellas, y el resto de las propiedades
del sistema se pueden determinar a partir de propiedades intensivas.
La regla de las fases de Gibbs define el número mínimo de variables intensivas necesarias
para describir un sistema:
f = c - p + 2 (1)
donde f son los grados de libertad del sistema, y se define como el número mínimo de
variables intensivas necesarias para describir un sistema, c es el número de componentes
químicos y p es el número de fases del sistema.

Ecuación de equilibrio líquido-vapor.

Cuando un sistema multicomponente se encuentra en equilibrio en dos o más fases, se


cumple:

a).- La temperatura de todo el sistema es uniforme y constante (equilibrio térmico).


b).- La presión en todo el sistema es uniforme y constante (equilibrio mecánico).
c).- El potencial químico y por lo tanto la fugacidad parcial de cada componente es el mismo
en cada una de las fases (equilibrio termodinámico).

Para el caso en que las fases que forman el sistema en equilibrio sean las fases líquida y
vapor, la fugacidad parcial de cada componente en esas fases está dado por las expresiones:

Fase líquida:
fˆi =  i xi  i P1 F .P.i
l sat sat

Fase vapor:
fˆi = ˆi yi P
vap

Igualando ambas expresiones, se obtiene la ecuación general de equilibrio líquido-vapor


(método Gamma-phi).
ˆi yi P =  i xi  i sat Pi sat F .P.i

Si el sistema se encuentra a presiones bajas (


F .P.i ˆi , i = 1), y si los componentes del
sat

sistema son de constitución química semejante i = 1, la ecuación general de equilibrio


líquido- vapor adquiere su forma más simple posible conocida como Ley de Raoult.
y i P = xi Pi
sat

Los sistemas líquidos ideales obedecen a la Ley de Raoult, de modo que la presión parcial de
cada uno de los componentes en el vapor, a una temperatura dada, es igual al producto de su
fracción molar en el líquido por la presión de vapor del componente puro a la misma
temperatura. Si el vapor también se comporta idealmente, la presión parcial de un
componente cualquiera en la fase vapor será igual al producto de su fracción molar en el
vapor por la presión total.

Laboratorio de Termodinámica del Equilibrio de Fases 3


Práctica 3 y 4: Equilibrio Ideal y No Ideal

Esto, para una mezcla binaria:

…… (1)
donde y1 e y2 son las fracciones molares de los componentes 1 y 2 en la fase vapor, x1 y x2
son las mismas en la fase líquida, P es la presión total y y las presiones de vapor de los
componentes puros respectivos a la temperatura del sistema.

Modelos de solución en mezclas multicomponentes.

Como consecuencia del análisis de datos de ELV experimentales existentes para sustancias
de constitución química semejante en solución con otra de referencia, se ha podido observar
el efecto que cada grupo químico constituyente de la sustancia tiene sobre el comportamiento
de la solución.

A la luz de esta idea, una solución se considera como una mezcla de ''grupos'' y no como
mezcla de sustancias. Esto ha permitido predecir los valores para los coeficientes de actividad
que una sustancia tiene en solución para la que no se tienen datos experimentales, basándose
únicamente en el conocimiento de la constitución química de la sustancia.

Como ejemplo de métodos de solución predictivos usados en problemas de ELV para


mezclas multicomponentes se menciona el UNIFAC.

Desafortunadamente los métodos predictivos aún no son confiables, y para usos prácticos, se
deben validar en mayor o menor grado por el experimento.

DESARROLLO EXPERIMENTAL

Laboratorio de Termodinámica del Equilibrio de Fases 4


Material y equipo • Reóstato
• Pinzas de tres dedos
• Hervidor • Pinzas de nuez
• Refrigerante • Pipeta muestreadora
• Termómetro • Tubos de ensaye
• Parrilla de calentamiento con • Refractómetro
agitación magnética • Metanol
• Soporte universal • Isopropanol

Procedimiento experimental

1. Montar el equipo como se muestra en el esquema.


2. Preparar mezclas binarias que contengan 10, 20, 40, 60, 80 y 90% de uno de los
componentes:
i. Metanol-Isopropanol

3. Los volúmenes de cada solución pueden calcularse con la ayuda de las siguientes
fórmulas, tomando en cuenta que los subíndices 1 son relacionadas al metanol y el 2
serán de isopropanol.

x1 PM 1 VT x 2 PM 2 VT
1 2
V1 = V1 =
x1 PM 1 x 2 PM 2 x1 PM 1 x 2 PM 2
+ +
1 2 1 2

4. Medir los índices de refracción de cada mezcla a una temperatura determinada (15 a
20o), empleando un refractómetro de Abbe, así como los índices de refracción de los
componentes puros.
5. Graficar la relación entre composición (x) y el índice de refracción (y). (El índice de
refracción permite conocer la composición del líquido y del vapor en equilibrio de
manera indirecta).
6. Poner la solución en el hervidor, suministrar calor con agitación hasta que la solución
ebulla, cuando la temperatura sea constante y haya condensado en la concavidad,
tomar el dato de temperatura.
7. Determinar los puntos de ebullición de cada una de las mezclas en equilibrio del vapor
y del líquido. Lo anterior se debe realizar por medio de destilación simple midiendo
2 ml. de mezcla. (nota: el bulbo del termómetro debe estar inmerso en el vapor
generado).
8. Tomar una muestra del condensado y del destilado y enfriarlas a temperatura
ambiente.
9. Determinar los índices de refracción del destilado (vapor) y del residuo que quedó en
el matraz de destilación (líquido) con el refractómetro de Abbe.
10. Repetir los pasos 6 al 9 para cada una de las mezclas.

1
Práctica 3 y 4: Equilibrio Ideal y No Ideal

condensador

agua de
recirculación

termómetro

matráz
concavidad
mezcla

Cálculos

Con los datos obtenidos calcular:


• Graficar índice de refracción contra composición
• Utilizando la gráfica determinar la composición del condensado (Y).
• Si las gráficas resultan ser de comportamiento ideal. Se tendrá que hacer lo siguiente:
• Aplicar la Ley de Raoult para calcular la composición del líquido y del vapor
aplicando las ecuaciones:
P − P2s P1s
x1 = s y y = x
P1 − P2s 1 1
P
• Calcular las presiones de saturación aplicando la ecuación de Antoine a las
temperaturas seleccionadas entre la temperatura de ebullición del componente 1 y la
temperatura de ebullición del componente 2. (buscar las constantes de Antoine en la
literatura).
• Graficar T vs x1, y1, y y1 vs x1
• Si las gráficas resultan ser de comportamiento no ideal. Se tendrá que hacer lo
siguiente:
• Calcular el coeficiente de actividad (γ) con la ecuación:
y1 P
1 =
x1 P1s
• Para el cálculo de las presiones aplicar la ecuación de Antoine a las temperaturas de
ebullición experimentales (buscar valores de las constantes en la literatura).
• Sacar el logaritmo natural de los coeficientes de actividad experimentales y graficar
Ln γ1 y γ2 vs x1 y x2 respectivamente en la misma gráfica. Esta gráfica servirá para
determinar los valores de las constantes A12 y A21 necesarias para aplicar las
ecuaciones de Van Laar y de Margules.
• Calcular las Ln γ1 y γ2 con las ecuaciones de Van Laar
• Calcular las Ln γ1 y γ2 con las ecuaciones de Margules

Laboratorio de Termodinámica del Equilibrio de Fases 2


Práctica 3 y 4: Equilibrio Ideal y No Ideal

• Comparar los resultados de Ln γ1 y γ2 con los experimentales y decir cual de las dos
ecuaciones se ajustan a los resultados experimentales.
• Identificar los datos del punto azeotrópico.
TABLA DE DATOS EXPERIMENTALES

X1 X2 V1 V2 ηliq T ηvap Y

Conceptos

Entregar de tarea lo siguiente ( derecho a hacer la práctica )

Equilibrio Líquido- Vapor

Comportamiento ideal

Presión de saturación

Ley de Raoult

Laboratorio de Termodinámica del Equilibrio de Fases 3


Práctica 3 y 4: Equilibrio Ideal y No Ideal

Ecuación de Antoine

Reporte.

En equipo: 50%
1. Cálculos de acuerdo a como se verifica en la práctica.
2. Reportar los resultados en forma tabular y gráfica.
3. Realizar las siguientes gráficas:
Índice de refracción Vs Concentración (ambos sistemas)
x1 vs y1 (para ambos sistemas, experimental y corregida)
x1, y1 vs T (para ambos sistemas, experimental y corregida)
4. Tablas de resultados experimentales, Análisis generales
5. Observaciones

Individual: 50%

Conclusiones

Laboratorio de Termodinámica del Equilibrio de Fases 4


Práctica 3 y 4: Equilibrio Ideal y No Ideal

INSTITUTO POLITÉCNICO NACIONAL


ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS

DEPARTAMENTO DE INGENIERIA QUIMICA INDUSTRIAL

LABORATORIO DE
TERMODINÁMICA DEL EQUILIBRIO DE FASES

PRÁCTICAS 4

EQUILIBRIO LÍQUIDO-VAPOR
EQUILIBRIO NO IDEAL - AZEÓTROPOS

Laboratorio de Termodinámica del Equilibrio de Fases 5


Práctica 3 y 4: Equilibrio Ideal y No Ideal

EQUILIDRIO NO IDEAL - AZEÓTROPOS

OBJETIVOS

• Observar el comportamiento de dos sistemas binarios.


• Construir los diagramas de temperatura y composición del líquido y del vapor.
• Comparar con los datos de temperatura y composición obtenidas con la Ley de Raoult.
• Comparar con los datos de temperatura y composición obtenidas con las ecuaciones de Van
Laar y Margules.
• Ajustar los datos experimentales utilizando dos modelos de solución y compararlos con los
datos experimentales.

PROPÓSITO
- Determinación experimental de los coeficientes de actividad de los componentes de
una mezcla binaria.
- Determinar las curvas de equilibrio líquido-vapor para mezclas binarias de líquidos
totalmente miscibles.
- Determinar si las mezclas obedecen a la ley de Raoult o presentan desviación
positiva o negativa.
- Realización del test de consistencia termodinámica de los datos obtenidos.
- Verificación de las ecuaciones de Margules y Van Laar para la determinación de
los coeficientes de actividad

HIPÓTESIS
A una presión constante las mezclas de líquidos miscibles de diferente composición hervirán
a diferentes temperaturas y su vapor será más rico en el componente más volátil.

EQUILIBRIO LÍQUIDO VAPOR

I. INTRODUCCION

La mayoría de las soluciones líquidas de interés práctico no cumplen con la Ley de Raoult.
Para poder tratar estas soluciones llamadas reales, en forma análoga a las soluciones perfectas
se introduce el coeficiente de actividad. Sin embargo, para las soluciones en fase vapor en el
dominio de presiones menores de 2 atm, puede aceptarse que cumplen con la Ley de Dalton.

En consecuencia, para las condiciones de trabajo de la presente experiencia, las ecuaciones a


aplicar son:
….(2)
donde y son los coeficientes de actividad para los componentes (1) y (2)
respectivamente.

Laboratorio de Termodinámica del Equilibrio de Fases 6


Práctica 3 y 4: Equilibrio Ideal y No Ideal

A partir de estas expresiones se calcularán los coeficientes de actividad determinando en el


equilibrio la presión, la temperatura y la composición de la fase vapor para diferentes
composiciones de la fase líquida.

En la práctica, esto significa construir el diagrama de equilibrio de fases T - x a presión


constante o el P - x a temperatura constante. Se debe notar el hecho de que estos coeficientes
de actividad son función de la composición del sistema como así también de la presión y la
temperatura. Trabajando a presión y temperatura constante, la relación de variación de los
coeficientes de actividad con la composición debe cumplir con la Ecuación de Gibbs-
Dühem, que para una mezcla binaria puede plantearse:

…… (3)

Equilibrio líquido-vapor en mezclas azeotrópicas

El equilibrio liquido-vapor de sistemas ideales se pueden describir de acuerdo a las leyes de


Raoult y Dalton. fracciones mol. De acuerdo a la Ley de Dalton la presión total del sistema
queda definida como:
PT = PA + PB (4)
En una mezcla homogénea de dos líquidos volátiles A y B, de acuerdo a la ley de Raoult se
pueden presentar dos casos:

1. Si las fuerzas intermoleculares entre las


moléculas A y B son más débiles que las
que existen entre las moléculas de A y
entre las de B, se dice que existe una
mayor tendencia para considerarse una
mezcla real que en el caso de una solución
ideal. En consecuencia, la presión de
vapor de la solución es mayor que la
predicha por la ley de Raoult para la
misma concentración. Esta conducta da
origen a la desviación positiva. El
mezclado es un proceso endotérmico.

2. Si las moléculas A y B se atraen con


mayor fuerza que entre las del mismo tipo,
la presión de vapor de la solución es menor
que la suma de las presiones parciales,
según la ley de Raoult. Dicha conducta se
le conoce como desviación negativa. El
mezclado es un proceso exotérmico.

Laboratorio de Termodinámica del Equilibrio de Fases 7


Práctica 3 y 4: Equilibrio Ideal y No Ideal

Variación de la composición de líquido y vapor durante la destilación.

Si bien no se conoce la relación analítica que vincula la concentración con el coeficiente de


actividad se han propuesto ecuaciones semiempíricas, que son soluciones de la ecuación (3).
Dos de las más comunes son las siguientes:

MARGULES

ln  1 = x22 A + 2 x1 ( B − A)  ; 
ln  2 = x12 B + 2 x2 ( A − B) 

VAN LAAR

A B
ln  1 = 2
; ln  2 = 2
 Ax   Bx2 
1 + 1  1 + 
 Bx2   Ax1 

Evidentemente que para el uso de tales ecuaciones es necesario conocer el valor de las
constantes A y B, las que se pueden determinar a partir de datos experimentales.
En una mezcla de gases ideales la fugacidad de sus componentes es igual a su presión parcial.
Generalmente los gases puros a bajas presiones se comportan como gases ideales y sus
mezclas también, por lo tanto la fugacidad de los componentes de la mezcla es igual a su
presión parcial. La presión atmosférica a la cual se trabaja en el laboratorio y en la mayoría
de las columnas de destilación se considera a baja presión, entonces el vapor en equilibrio
con el líquido se considera gas ideal.
En cambio los líquidos no se pueden suponer como gases ideales, pero a la presión der
saturación (baja), el coeficiente de fugacidad es considerado igual a 1 y la fugacidad del
líquido saturado o aún comprimido se toma igual a su presión de vapor.

Procedimiento experimental

11. Montar el equipo como se muestra en el esquema.


12. Preparar mezclas binarias que contengan 10, 20, 40, 60, 80 y 90% de uno de los
componentes:
i. Metanol- Cloroformo

13. Los volúmenes de cada solución pueden calcularse con la ayuda de las siguientes
fórmulas, tomando en cuenta que los subíndices 1 son relacionadas al metanol y el 2
serán de isopropanol.

Laboratorio de Termodinámica del Equilibrio de Fases 8


Práctica 3 y 4: Equilibrio Ideal y No Ideal

x1 PM 1 VT x 2 PM 2 VT
1 2
V1 = V1 =
x1 PM 1 x 2 PM 2 x1 PM 1 x 2 PM 2
+ +
1 2 1 2

14. Medir los índices de refracción de cada mezcla a una temperatura determinada (15 a
20o), empleando un refractómetro de Abbe, así como los índices de refracción de los
componentes puros.
15. Graficar la relación entre composición (x) y el índice de refracción (y). (El índice de
refracción permite conocer la composición del líquido y del vapor en equilibrio de
manera indirecta).
16. Poner la solución en el hervidor, suministrar calor con agitación hasta que la solución
ebulla, cuando la temperatura sea constante y haya condensado en la concavidad,
tomar el dato de temperatura.
17. Determinar los puntos de ebullición de cada una de las mezclas en equilibrio del vapor
y del líquido. Lo anterior se debe realizar por medio de destilación simple midiendo
2 ml. de mezcla. (nota: el bulbo del termómetro debe estar inmerso en el vapor
generado).
18. Tomar una muestra del condensado y del destilado y enfriarlas a temperatura
ambiente.
19. Determinar los índices de refracción del destilado (vapor) y del residuo que quedó en
el matraz de destilación (líquido) con el refractómetro de Abbe.
20. Repetir los pasos 6 al 9 para cada una de las mezclas.

condensador

agua de
recirculación

termómetro

matráz
concavidad
mezcla

Cálculos

Con los datos obtenidos calcular:


• Graficar índice de refracción contra composición
• Utilizando la gráfica determinar la composición del condensado (Y).
• Si las gráficas resultan ser de comportamiento ideal. Se tendrá que hacer lo siguiente:

Laboratorio de Termodinámica del Equilibrio de Fases 9


Práctica 3 y 4: Equilibrio Ideal y No Ideal

• Aplicar la Ley de Raoult para calcular la composición del líquido y del vapor
aplicando las ecuaciones:
P − P2s P1s
x1 = s y y1 = x1
P1 − P2s P
• Calcular las presiones de saturación aplicando la ecuación de Antoine a las
temperaturas seleccionadas entre la temperatura de ebullición del componente 1 y la
temperatura de ebullición del componente 2. (buscar las constantes de Antoine en la
literatura).
• Graficar T vs x1, y1, y y1 vs x1
• Si las gráficas resultan ser de comportamiento no ideal. Se tendrá que hacer lo
siguiente:
• Calcular el coeficiente de actividad (γ) con la ecuación:
y1 P
1 =
x1 P1s
• Para el cálculo de las presiones aplicar la ecuación de Antoine a las temperaturas de
ebullición experimentales (buscar valores de las constantes en la literatura).
• Sacar el logaritmo natural de los coeficientes de actividad experimentales y graficar
Ln γ1 y γ2 vs x1 y x2 respectivamente en la misma gráfica. Esta gráfica servirá para
determinar los valores de las constantes A12 y A21 necesarias para aplicar las
ecuaciones de Van Laar y de Margules.
• Calcular las Ln γ1 y γ2 con las ecuaciones de Van Laar
• Calcular las Ln γ1 y γ2 con las ecuaciones de Margules
• Comparar los resultados de Ln γ1 y γ2 con los experimentales y decir cual de las dos
ecuaciones se ajustan a los resultados experimentales.
• Identificar los datos del punto azeotrópico.

Cálculo de la concentración de etanol utilizando índices de refracción


El refractómetro que se utiliza en el laboratorio está graduado para medir el índice de
refracción, por ejemplo:
Lectura con el refractómetro (índice de refracción) = 1.35162. De la gráfica que correlaciona
índice de refracción y % w/w de etanol, obtenemos la concentración de etanol en la mezcla
= 26.06%

Conceptos

Laboratorio de Termodinámica del Equilibrio de Fases 10


Práctica 3 y 4: Equilibrio Ideal y No Ideal

Entregar de tarea lo siguiente ( derecho a hacer la práctica )


Comportamiento no ideal

Mezcla y punto azeotrópica (o)

Presión de saturación

Ecuación de Van Laar

Ecuación de Margules

Actividad

Coeficiente de actividad

Reporte.

En equipo: 50%
1. Cálculos de acuerdo a como se verifica en la práctica.
2. Reportar los resultados en forma tabular y gráfica.
3. Realizar las siguientes gráficas:
Índice de refracción Vs Concentración (ambos sistemas)
x1 vs y1 (para ambos sistemas, experimental y corregida)
x1, y1 vs T (para ambos sistemas, experimental y corregida)
4. De acuerdo a los resultados analizar si la mezcla presenta un punto
azeotrópico y además, determinar qué tipo de comportamiento presenta de
acuerdo a la ley de Raoult.
5. Tablas de resultados experimentales, Van Laar y Margules
6. Análisis generales
7. Observaciones

Laboratorio de Termodinámica del Equilibrio de Fases 11


Práctica 3 y 4: Equilibrio Ideal y No Ideal

Individual: 50%
Conclusiones

Laboratorio de Termodinámica del Equilibrio de Fases 12


INSTITUTO POLITÉCNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERÍA QUÍMICA E INDUSTRIAS EXTRACTIVAS

DEPARTAMENTO DE INGENIERIA QUIMICA INDUSTRIAL

LABORATORIO DE
TERMODINÁMICA DEL EQUILIBRIO DE FASES

PRÁCTICA 5

EQUILIBRIO LÍQUIDO-LÍQUIDO
SISTEMAS TERNARIOS

1
Práctica 5: Equilibrio Líquido-Líquido: Sistemas Ternarios

SISTEMAS TERNARIOS

OBJETIVOS

• Estudiar las relaciones de solubilidad del sistema de tres componentes.


• Trazar el diagrama de equilibrio para un sistema estudiado.
• Construir la curva de solubilidad (binodal) para una temperatura dada y determinar
las líneas de reparto.
• Construir los diagramas de temperatura y composición del líquido.
• Ajustar los datos experimentales utilizando dos modelos de solución y compararlos
con los datos experimentales.

PROPÓSITO

Mediante la preparación de mezclas agua – cloroformo – ácido acético, obtener datos


experimentales a fin de trazar la curva binodal o de equilibrio para el sistema en estudio.
Se determinará la curva de solubilidad y las líneas de unión de un sistema integrado por agua,
metanol y un líquido orgánico que puede ser: tolueno, benceno, tricloroetileno, cloroformo o
1-butanol. En este caso cloroformo. Mediante otra serie de experimentos trazar dos líneas de
unión

HIPÓTESIS

Se podrá determinar la curva de solubilidad y las líneas de unión de un sistema ternario


integrado por agua, metanol y cloroformo, para un equilibrio líquido-líquido. Todas las
mezclas cuya composición está por debajo de la curva de solubilidad se separaran en dos
capas. Las que están por encima de ella, formarán una sola fase.

EQUILIBRIO LÍQUIDO LÍQUIDO

I. INTRODUCCION

Es posible tratar todos los equilibrios heterogéneos desde un punto de vista unificado por
medio del principio conocido como Regla de las Fases, con el cual él numero de variables a
que se encuentra sometido un equilibrio heterogéneo cualquiera, queda definido bajo ciertas
condiciones experimentales definidas. La cantidad de variables que pueden fijarse
arbitrariamente en un sistema termodinámico puede calcularse mediante: L = C − F + 2 (1).

Dra. Carmen Reza Laboratorio Termodinámica Química I 2


Práctica 5: Equilibrio Líquido-Líquido: Sistemas Ternarios

Donde L número de grados de libertad, C número de componentes químicos, F número de


fases.

Para simplificar la representación gráfica de las condiciones de equilibrio se suponen


temperatura y presión constantes y, además, generalmente se considera que el sistema en
estudio es condensado, es decir que se desprecia la fase vapor. En estas condiciones, las
variables independientes son las concentraciones de dos componentes, siendo la
concentración del tercero una función de ellas.

De acuerdo a la regla de las fases de Gibbs, cuando observa una sola fase en un sistema de
tres componentes como el que estudiaremos, los grados del libertad son 4. Por lo tanto, para
describirlo completamente habrá que fijar 4 de las 5 variables del sistema (T, P y
concentraciones de cada uno de los tres componentes). Como los gráficos de tantas variables
son muy difíciles de interpretar, se elige mantener algunas de ellas constantes y graficar las
restantes una contra otras. En nuestro caso se trabajará en condiciones de presión y
temperatura constantes y se graficará el número de fases del sistema respecto de las
concentraciones de sus tres componentes en un diagrama triangular, en unidades de
porcentaje en masa.

Para un sistema dado de tres componentes líquidos, existirán composiciones para las cuales
la solubilidad es completa, resultando la mezcla en una sola fase. Entonces, a P y T cte., serán
2 los grados de libertad, debiendo establecerse dos de las tres concentraciones para describir
completamente la situación del sistema. Pero pueden darse composiciones en las cuales se
supera la solubilidad y aparecen dos fases inmiscibles, cada una con los tres componentes en
cierta proporción. Los diagramas de líquidos ternarios son de considerable valor en
problemas de separación y extracción con solventes y gracias a ellos es posible deducir si la
separación buscada puede realizarse y en cuales condiciones de operación para lograr
resultados óptimos. Los diagramas de líquidos ternarios son de considerable valor en
problemas de separación y extracción con solventes y gracias a ellos es posible deducir si la
separación buscada puede realizarse y en cuales condiciones de operación para lograr
resultados óptimos.

La mezcla de tres componentes líquidos, en nuestro caso cloroformo, ácido acético y agua,
en distintas proporciones puede dar lugar a la formación de dos fases. En el caso que nos
ocupa, el agua y el cloroformo son inmiscibles (es decir, la extensión de su mezcla es
prácticamente nula) y sin embargo, el ácido acético es totalmente soluble en cualquiera de
los dos productos. Al añadir ácido acético a una mezcla agua-cloroformo se observa que éste
se distribuye entre las dos fases al mismo tiempo que aumenta la solubilidad mutua entre las
mismas. Este hecho se traduce en que adiciones sucesivas de ácido acético hacen variar la
composición de las fases acuosa y orgánica haciéndolas cada vez más semejantes y
apreciándose, como característica particular de este proceso, un aumento de volumen de una
de las fases y la consiguiente disminución del volumen de la otra hasta que llega a desaparecer
una de ellas.

Nuestro objetivo es construir y aprender a manejar un diagrama ternario, determinando la


curva de solubilidad del sistema ternario por titulación hasta la aparición ó desaparición de
dos fases. Esta curva límite separa la zona de composiciones que dan un sistema monofásico

Dra. Carmen Reza Laboratorio Termodinámica Química I 3


Práctica 5: Equilibrio Líquido-Líquido: Sistemas Ternarios

de las que dan un sistema bifásico. Esta práctica abarcará las relaciones de solubilidad del
sistema de tres componentes como lo son el cloroformo-ácido acético-agua, se trazará el
diagrama de equilibrio y se construirá la curva de solubilidad (binodal) para una temperatura
dada y por último se determinarán las composiciones de los puntos críticos y máximos.

DESARROLLO EXPERIMENTAL

Material y reactivos

• Probetas 100 ml
• Pipetas de 10 ml
• Buretas de 50 ml
• Balanza
• Matraces Erlenmeyer
• Vasos de precipitado
• Soporte Universal
• Pinzas para buretas
• Agua destilada
• Cloroformo
• Ácido acético
• Solución de NaOH 0.5N
• Fenolftaleína

El procedimiento experimental se divide en dos partes:

1. Construcción de la curva de solubilidad (binodal)


2. Construcción de las líneas de unión

Determinación de la curva de solubilidad

1.- Colocar 10 ml de cloroformo en un matraz erlenmeyer.


2.- Titular lentamente con agua, agitar intensamente después de agregar cada gota
hasta la aparición de turbidez (agitar intensamente después de la adición de cada gota,
ya que para este punto solo se necesitan unas cuantas gotas (4 gotas)) tras por lo
menos un minuto de intensa agitación.
3.- Transformar los volúmenes empleados en masas mediante las densidades, y
número de moles mediante pesos moleculares, de las sustancias a la temperatura de
trabajo (T ambiente), mediante la siguiente ecuación:

m
m=V n=
PM

4.- Calcular la composición porcentual en masa del punto y graficarla en un diagrama


triangular.

Dra. Carmen Reza Laboratorio Termodinámica Química I 4


Práctica 5: Equilibrio Líquido-Líquido: Sistemas Ternarios

grT= grCH3COOH + grCHCl3 + grH2O

% p mx
= * 100
px grT

5.- Agregar al sistema 2 ml de ácido acético, con lo cual desaparecerá la turbidez, y


volver a titular con agua, hasta nueva aparición de turbidez.
6.- Transformar nuevamente los volúmenes en masas (ahora teniendo en cuenta
también el ácido acético) y volver a calcular y graficar la composición porcentual en
masa.
7.- Continuar procediendo de esta forma hasta llegar al máximo de la curva de
solubilidad, tratando de obtener puntos equiespaciados, para lo cual puede ser
necesario adicionar sólo de 2 a 3 ml de ácido acético a medida que se avanza,
(asegurarse que en cada caso la mezcla pierda la turbidez, hasta alcanzar un volumen
máximo de ácido adicionado de 15 ml aproximadamente).
8.- Una vez alcanzado el máximo de la curva, repetir todo el procedimiento pero
partiendo de 10 ml de agua y titulando con cloroformo.
9.- Titular lentamente con agua (agitando intensamente después del agregado de cada
gota) hasta la aparición de turbidez (lechoso) tras por lo menos un minuto de intensa
agitación.
NOTA: Cuando la sustancia titulante es el cloroformo un indicador de que se esta
próximo al equilibrio es cuando la solución se torna turbia, pero se observa que al
dejar de agitar esta se vuelve transparente, por tanto cuando se llega a este punto se
continua titulando gota a gota agitando intensamente después de la adición de cada
gota, dejar de agitar y ver si en el fondo del matraz existe una pequeña gota de
cloroformo sin disolver en ese momento se detiene la titulación anotando el volumen
gastado de cloroformo.

Determinación de las líneas de unión:

1.- Colocar 10 ml de agua e igual volumen de cloroformo en 3 vasos de precipitado


previamente pesados.
2.- Agregar 2, 5 y 10 ml de ácido acético en cada uno de los tres vasos de precipitado,
etiquetar como solución 1, 2 y 3 respectivamente, pesar cada solución anotando el
valor en la tabla correspondiente.
3.- Transferir la solución 1 a un embudo de separación, agitar vigorosamente y esperar
a que las fases se separen.
4.- Mientras se separan las fases, calcular la composición porcentual en masa del
sistema y marcar el punto correspondiente en el diagrama triangular (pasos 3 y 4 en
determinación de la curva de solubilidad).
5.- Medir (ml) y pesar (gr) cada una de las fases, anotando los datos en la tabla
correspondiente.

Dra. Carmen Reza Laboratorio Termodinámica Química I 5


Práctica 5: Equilibrio Líquido-Líquido: Sistemas Ternarios

6.- Transferir aproximadamente 2 ml de cada una de las fases a matraces erlenmeyer,


diluir con 20 ml de agua y titular con solución de NaOH (0,5 M), usando fenolftaleína
como indicador.
7.- Calcular el porcentaje en masa de ácido acético de las muestras extraídas, marcar
los puntos correspondientes en el diagrama y trazar la línea de unión.
(Estos dos puntos deben estar alineados con el punto correspondiente a la
composición total del sistema).

a) Para LA FASE SUPERIOR ACUOSA, por ejemplo para un volumen gastado de


NaOH de 0.8 ml, y con los pesos y volúmenes de cada fase, se sabe por el principio
de equivalencia que:
VNaOH M NaOH = Va.acético M a.acético
por lo que la concentración del ácido acético presente en la alícuota de 2 ml queda
como:
V M 0.8mL * O.5M
M a.acético = NaOH NaOH = = 0 .2 M
Va.acético− alícuta 2mL

Con el valor de la molaridad (gmol/l) se calcula la cantidad en gramos de ácido acético


en la fase superior acuosa como sigue:

Wa.acético fase sup . acuosa


 gmol 
= M a.acético   V fase sup acuosa (L ) PM  g 

 L   gmol 

 gmol   g 
Wa.acético fase sup . acuosa = 0.2   0.026 (L ) 60   = 0.312 g
 L   gmol 

Como la fase superior acuosa esta integrada por ácido acético y agua, por el peso que
se tomó de la solución se tiene el valor de esta fase de:
W fase sup acuosa = 28.27 g

Para obtener el porciento de ácido acético en esta fase:

Wa.acético 0.312 g
= 100 = 100 = 1.1%
fase sup . acuosa
% a.acético fase sup . acuosa
W fase sup acuosa 28.27

b) Para LA FASE INFERIOR ORGANICA, se cuenta como datos el volumen


gastado de NaoH, el peso y el volumen de esta fase, se procede de la misma forma
que en inciso anterior.

8.- Repetir los pasos 3, 4, 5, 6 y 7 para las soluciones 2 y 3.

Dra. Carmen Reza Laboratorio Termodinámica Química I 6


CURVA DE SOLUBILIDAD

Volumen Volumen Volumen Peso Peso Peso Peso Porcent. Porcent. Porcent.
Clorofor Agua Ac.Acético Clorofor Agua Ac.Acét. Solución Clorofor Agua Ac.Acét.
(ml) (ml) (ml) (gr) (gr) (gr) (gr) (%) (%) (%)

10 0
10 3
10 (2)5
10 (2)7
10 (2)9
10 (2)11
10 (2)13
10 (2)15
10 0
10 (2)3
10 (2)5
10 (2)7
10 (2)9
10 (2)11
10 (2)13
10 (2)15

7
Práctica 5: Equilibrio Líquido-Líquido: Sistemas Ternarios

LÍNEAS DE UNIÓN
DATOS DE LAS SOLUCIONES

Volumen Volumen Volumen Peso Peso Peso Porcentaje Porcentaje Porcentaje


Cloroformo Agua Ac.Acético Cloroformo Agua Ac.Acético Cloroformo Agua Ac.Acético
(ml) (ml) (ml) (gr) (gr) (gr) (%) (%) (%)
10 10 2
10 10 5
10 10 10

FASE ACUOSA

Volumen Concent. Volumen Peso Concent. Peso Volumen Peso Porcentaje


NaOH NaOH Ac.Acético Solución Ac.Acético F-Acuosa F-Acuosa F-Acuosa Ac.Acético
(ml) (M) (ml) (gr) (M) (gr) (ml) (gr) (%)
2
2
2

FASE ORGÁNICA

Volumen Concent. Volumen Peso Concent. Peso Volumen Peso Porcentaje


NaOH NaOH Ac.Acético Solución Ac.Acético F-Orgánica F-Orgánica F-Orgánica Ac.Acético
(ml) (M) (ml) (gr) (M) (gr) (ml) (gr) (%)
10 2
10 2
10 2

Dra. Carmen Reza Laboratorio Termodinámica Química I 8


Práctica 5: Equilibrio Líquido-Líquido: Sistemas Ternarios

Dra. Carmen Reza Laboratorio Termodinámica Química I 9


Recomendaciones

• Hay que ser muy cuidadoso al momento de utilizar el ácido acético, debido a que este es muy
volátil, así que para la obtención de unos resultados óptimos es necesario que antes de verter
la cantidad del ácido en la bureta para la valoración, la mezcla de acetato y de agua esta lista
para titularse.
• Seria favorable para aumentar nuestros conocimiento de los comportamientos de sistemas de
tres líquidos con un par parcialmente miscibles, trabajar a diferentes temperaturas, y así
observar la variación de la curva de solubilidad a una temperatura alta con respecto a una
menor.
• Realizar experimentos con sistemas de tres líquidos con dos pares parcialmente solubles, dos
líquidos parcialmente soluble y un sólido y tres pares parcialmente solubles.

ANEXO 1

Tratamiento de datos:
Consiste en obtener la curva de solubilidad y las líneas de unión y debe hacerse en forma
paralela a la parte operativa. Ubique en el vértice superior el 100% de ácido acético, en el
vértice izquierdo el 100% de agua y en el derecho el 100% de cloroformo. Asegúrese de
comprender su significado y funcionamiento. Ubique en los lados correspondientes al punto
que indica la concentracion en % m/m de la solucion de partida que acaba de preparar.

Punto de saturación en el sistema ternario:

Pesar nuevamente el erlenmeyer con la solución de ácido acético en agua (o de ácido acético
en cloroformo). Titular la solución de acético en agua con cloroformo (o bien la solución de
ácido acético en cloroformo, con agua) vertiendo gota a gota desde una bureta de 10ml, con
buena agitación hasta un punto de turbidez permanente. Es el punto de saturación del sistema
ternario, a partir del cual se separan dos fases. Pesar nuevamente el erlenmeyer y calcular por
diferencia la masa de cloroformo (o de agua) agregada hasta llegar a ese punto de saturación.
Calcular la nueva composición porcentual del sistema. Registrar todos los datos en su tabla.
Ubicar ese punto en el diagrama ternario.

Recopilación de datos:

En cada repetición del procedimiento descrito, cada alumno o grupo de alumnos habrá
determinado uno de los puntos que pertenecen a la curva de solubilidad del sistema agua-
ácido acético-cloroformo a presión y temperatura constantes. Ubicando todos los datos
experimentales en un gráfico común a todos los grupos, podrá dibujarse la curva más
probable que separa en el diagrama las zonas de composiciones que resultan en una sola fase
de aquellas que resultan en sistemas bifásicos. Cada grupo agregará a su propio diagrama
ternario la curva límite común, indicando el número de fases en cada zona, las unidades de
composición en los ejes, y la temperatura y presión a la que fue realizada la experiencia (esto
es, en la que el diagrama es aplicable).

10
Práctica 5: Equilibrio Líquido-Líquido: Sistemas Ternarios

ANEXO 2: Cálculos típicos de la composición de las mezclas de tres componentes.

Ácido Acético:

ρCH3COOH= 1.050 gr/ ml

PMCH3COOH= 60 gr/ mol

Cloroformo:

ρCHCl3= 1.48 gr/ ml

PMCHCl3= 119.4 gr/ mol

Agua:

ρH2O= 1.00 gr/ ml

PMH2O= 18 gr/ mol

Calculo de las líneas de Unión:

Para la fase orgánica formada por cloroformo y agua

Para la disolución:

nNaOH= [NaOH] * V

si VNaOH = 5.4 ml = 0.0054 L

nNaOH=

Para la fase acuosa formada por agua y cloroformo:

nNaOH= [NaOH] * V

si VNaOH = 4.8 ml = 0.0048 L

nNaOH=

Dra. Carmen Reza Laboratorio Termodinámica Química I 11


Práctica 5: Equilibrio Líquido-Líquido: Sistemas Ternarios

las fracciones molares serán:

Xfor= 0.47059

Xfac= 0.52941

Los porcentajes en peso serán:

Disolución # 1:

Dra. Carmen Reza Laboratorio Termodinámica Química I 12


REPORTES DE PRÁCTICAS DE LABORATORIO
En equipo: 50%
• Carátula bien definida, destacando el título completo de la práctica y nombres de alumnos
que conforman el equipo.
• Objetivos que se persiguen en la práctica y una breve relación de los pasos a seguir para
lograrlos (no tiene que ver con el desarrollo experimental) ejemplo: relación que se estudia,
ecuación que se aplica, etapas de que consta el trabajo.
• Introducción (no copiar de la práctica ni de internet).
• Procedimiento experimental (diagramas con dibujos si son necesarios)
• Gráficas (si se piden o si son necesarias).
• Cálculos y resultados (los necesarios y bien desarrollados)
• Diagramas (si son necesarios)
• Observaciones (deben redactarse solo en los casos en los que se haya producido un
cambio significativo o dificultad en el desarrollo experimental, lo cual deberá tenerse en
cuenta al elaborar las conclusiones).

Individual: 50%
1. Conclusiones

IMPORTANTE
• las conclusiones deberán ser redactadas de manera individual por cada estudiante, lo que permitirá al
mismo realizar un análisis profundo de sus resultados. En las conclusiones deberá tocarse los siguientes
puntos (teniendo en cuenta las características propias de la práctica realizada).
a. Correspondencia de los resultados obtenidos con los esperados o con valores teóricos, si
corresponden entre si, si son lógicos o no según e procedimiento realizado.
b. Información práctica o teórica que brindan estos resultados.
c. Si quedó demostrado lo que se pretendía.
d. Si los resultados no son apropiados, o no son los esperados, deben evaluarse las posibles
causas que provocan estos resultados.
e. Las conclusiones deben estar escritas a mano, en tinta y letra clara, evitando tachaduras y
borrones o cualquier otra manifestación de descuido en su presentación.
f. Las conclusiones deberán ir al final del reporte en equipo y tener nombre y firma del alumno
correspondiente.

NO SE ACEPTAN CONCLUSIONES QUE:


i. describan lo que se hizo en la práctica: ejemplo: primero se arrancó la bomba de
vacío luego se comenzó a calentar el líquido en el hervidor luego… eso
NO es conclusión.
ii. exponga lo que se aprendió en la práctica: ejemplo: en la práctica aprendí a calcular
el calor de vaporización…. Eso NO es conclusión.
iii. que sean párrafos de libros…. Deben ser conclusiones de cada alumno no del
libro.
iv. que sean muy largas o muy cortas. (tres conclusiones como mínimo utilizando
aproximadamente 1 cuartilla por las tres)
v. que sean muy parecidas a las de los demás compañeros (redacción)

1. Los alumnos que presenten conclusiones que sean de ese estilo tendrán calificaciones
reprobatorias y no se les podrá tomar en cuenta otra conclusión de la misma práctica.
2. Aquellos alumnos que presenten muy buenas conclusiones en la primer práctica y
tengan calificaciones mayores del 95% se abstendrán de presentar conclusiones en las
prácticas siguientes, no así tendrán que seguir ayudando a sus compañeros de equipo,
tanto en la práctica y en los reportes de prácticas.
3. Se califica también la redacción y ortografía.

1
Práctica 5: Equilibrio Líquido-Líquido: Sistemas Ternarios

ORIENTACIONES PARA EL TRABAJO DEL LABORATORIO.

Son de estricto cumplimiento las siguientes medidas y orientaciones:


1. Encontrarse en el área de actividad de laboratorio a la hora indicada para su
comienzo.
2. Haber estudiado previamente los aspectos teóricos correspondientes y los pasos a
seguir en la práctica, conociendo el motivo por el cual se realizará cada uno de ellos
(se hará un examen antes del comienzo de la práctica teniendo validez en la
evaluación final) aquella persona que no conteste correctamente no tendrá derecho
a realizar la práctica.
3. Traer las Hojas de Seguridad de los reactivos utilizados en la práctica (obligatorio
haberlo leído con anticipación)
4. Entregar el informe correspondiente de la práctica anterior.
5. No portar anillos y otras prendas que puedan afectarse por contacto con ácidos u
otros reactivos.
6. Pedir permiso para abandonar el área de laboratorio siempre y cuando haya un
motivo justificado e impostergable para ello.
7. Permanecer en su puesto de trabajo durante la sesión de la práctica.
8. Mantener únicamente en la mesa de trabajo los materiales de laboratorio necesarios
para realizar la práctica.
9. Minimizar condiciones de riesgo.
10. Dejar limpio y ordenado el puesto de trabajo una vez finalizada la práctica; y pedir
inspección del profesor en turno.
11. Se prohíbe fumar e ingerir alimentos durante la sesión de práctica.
12. Cumplir con todas las medidas de protección y seguridad en el laboratorio entre las
que se encuentran:
a. Utilizar bata de laboratorio.
b. Utilizar lentes de seguridad.
c. Utilizar mascarilla de seguridad cuando así se requiera.
d. Evitar utilizar pupilentes mientras se realiza la práctica.
e. Evitar que se derramen los reactivos sobre la mesa de trabajo.
f. Evitar utilizar mas reactivo del necesario.
g. Evitar pipetear reactivo de su envase original (evitar contaminaciones, en
tal caso no vaciar mas reactivo del necesario en un vaso de precipitado)
h. No pipetear con la boca.
i. Informar al profesor en turno cuando haya ocurrido un accidente o
incidente durante la práctica.
j. No conectar equipos eléctricos sin la previa autorización del profesor
corroborando siempre que el voltaje del equipo corresponda con el voltaje
de la línea a la cual se conecta.
13. Realizar un buen TRABAJO EN EQUIPO.

PARA LA REALIZACIÓN DEL TRABAJO EN EL LABORATORIO ES NECESARIO:

Papel toalla (1 rollo por equipo).


Jabón (1 bolsa o envase por grupo)
Guantes para lavar material (un par por equipo)
Escobillón (3 diferentes tamaños por grupo (vaso de pp, pipeta, bureta, embudo, etc.))
Franela (1 por grupo)
Guantes de Látex (15 – 20 pares por grupo)

Dra. Carmen Reza Laboratorio Termodinámica Química I 2

También podría gustarte