Informe del Laboratorio de Fisicoquímica. Semestre B-2024, Junio-Noviembre, 2024.
Escuela de Ingeniería Química, Facultad de Ingeniería, Universidad de los Andes (ULA)
EQUILIBRIO LÍQUIDO-LÍQUIDO SISTEMAS BINARIOS
CON MISCIBILIDAD PARCIAL
Estudio de la miscibilidad del sistema 1-Butanol/Agua a través del modelo NRTL
Coello, Mariana¹; Maldonado, Alexandra², Rodríguez, Jesus3
1,2,3
Núcleo Mérida, Universidad de Los Andes,
Mérida, Estado Mérida, Venezuela.
Resumen
La solubilidad de las sustancias es una propiedad que depende de sus propiedades, por lo que algunas veces se puede
presentar casos de miscibilidad parcial entre compuestos. De estas se caracterizan las zonas de dos fases, donde los
delimita las soluciones conjugadas, en las cuales los componentes están en equilibrio químico por lo que se puede utilizar
la energía libre de Gibbs de exceso para calcular sus composiciones. Un método para estimar el coeficiente de actividad,
es por el modelo NRTL que se fundamenta en la teoría de dos líquidos de Scott. Se busca mediante un método experimental
comprender los conceptos de la miscibilidad parcial por medio de la determinación práctica de las composiciones de las
soluciones conjugadas de una mezcla de 1-Butanol/Agua. Además, estimar la posición de equilibrio termodinámico por
medio del modelo NRTL y estudiar su aplicabilidad en la descripción de equilibrios líquido-líquido.
Por medio de la experimentación, se logró determinar que el sistema 1-Butanol/Agua es parcialmente miscible y que el
butanol a temperatura ambiente solo se diluye hasta una concentración del 8%, composición de la solución conjugada. Se
hallaron las composiciones molares de los compuestos. Luego se utilizó el modelo NRTL para determinar sus coeficientes
de actividad, obteniendo valores positivos para ambos componentes, por lo que se tiene una desviación positiva de la ley de
Raoult
Palabras claves: Solubilidad, miscibilidad parcial, 1-butanol, modelo NRTL.
solución saturada. En medio se tienen dos fases (Levine,
[Link]ón 2004).
La solubilidad es una propiedad de los compuestos que
depende de sus propiedades, destacando la polaridad de las
sustancias involucradas. En la mayoría de los casos la
polaridad determina la solubilidad de un compuesto en otro,
siendo que lo semejante disuelve lo semejante. Al ocurrir
esto, se refiere a sustancias miscibles entre sí. Sin embargo,
existen ocasiones en donde dos sustancias pueden mezclarse
y disolverse mutuamente hasta alcanzar un límite máximo.
A esto se le llama miscibilidad parcial.
En los casos más comunes, cuando la presión se
encuentra a 1 atm se obtiene la forma del gráfico como se
muestra en la Imagen 1, donde al moverse por la línea
horizontal de T1, se puede describir las etapas que se Imagen 1. Diagrama de fases liquido-líquido para miscibilidad
forman en una mezcla con miscibilidad parcial. Entre F-G y parcial. Extraído de (Levine, 2004).
E-H se tiene una sola fase, y una solución con pequeñas
cantidades B y luego de A respectivamente. En los puntos G
y E se tienen las soluciones conjugadas, donde en G se
alcanza la solubilidad máxima de la solución B en la Al ocurrir un aumento de temperatura, va
solución C a la temperatura señalada; mientras que en E se disminuyendo progresivamente la zona de inmiscibilidad,
alcanza la cantidad añadida de la sustancia B suficiente para hasta que, al alcanzar Tc, los compuestos son
que todo C pueda disolverse en ella, resultando en una completamente miscibles entre sí.
Informe de Fisicoquímica, B-2024, junio-noviembre, 2024
2
Coello, Mariana; Maldonado, Alexandra, Rodríguez, Jesús
En la presente experimentación se busca entender estos
Una característica de las soluciones conjugadas, es que conceptos de miscibilidad parcial, por medio de la
los componentes se encuentran en equilibrio determinación experimental del equilibrio termodinámico
termodinámico, por lo que sus potenciales químicos pueden de soluciones conjugadas de 1-Butanol/Agua. Además,
igualarse para cada una de ellas. De esta manera y al utilizar estimar esta posición de equilibrio con el modelo NRTL y
el modelo de la Energía Libre de Gibbs de Exceso, se estudiar su aplicabilidad para describir un equilibrio
obtienen Ec. 1 y Ec. 2, donde XA1 y XA2 corresponden a líquido-líquido de dicho sistema.
la concentración de un componente en otro (Calderón S. M.,
Belandria, Cárdenas, Araque, & De Lima, 2022).
[Link] Experimental
L1 L2
x A 1∗γ =x A 2∗γ
A A Ec. 1
Primera experiencia: determinación de la solubilidad
L1 L2 mutua del sistema 1-Butanol-Agua a temperatura ambiente
x B 1∗γ =x B 2∗γ
B B Ec. 2
1- Tomar 2 erlenmeyer de cuello esmerilado con tapón,
Para la determinación del coeficiente de actividad se estos deben estar limpios y secos.
encuentra el modelo Non Random Two Liquids (NRTL).
Este modelo explica que al tener una mezcla de moléculas 2- Preparar las soluciones de 1-Butanol/Agua descritas en la
donde las interacciones son no ideales, estas se agrupan de tabla 2.
manera no aleatoria, y esto resulta en composiciones
localizadas que difieren de la composición global. 3- Tomar nota de la temperatura del laboratorio.
4- Pesar en una balanza analítica el reactivo a usar en cada
muestra.
5- Preparar las soluciones con ayuda de una pipeta acorde a
las muestras.
6- Agitar las soluciones y dejarlas en reposo en un lugar
oscuro (Gaveta)por una semana para asegurar el estado de
equilibrio
Imagen 2. Tipos de células según la teoría de dos líquidos en 7- Transcurrido el tiempo, tomar una inyectadora limpia y
mezclas binarias de Scott y Leland. Extraído de seca.
(Buentello, 2007)
8- Tomar una muestra de 0,3-0,5 ml de la fase orgánica y la
Sus fundamentos están en la teoría de los dos líquidos fase acuosa de la solución con la inyectadora, tomando las
de Scott, donde se explica que existen dos tipos de células precauciones necesarias para no perturbar el sistema.
como se muestra en la figura 2 y la figura 3. En el caso de
las células con molécula 1 que está en el centro, la energía 9- Verter el contenido de la inyectadora en el depósito de
de Gibbs residual es la sumatoria de todas las energías muestras del refractómetro.
residuales de las interacciones efectuadas. Análogamente
ocurre para las células donde la molécula 2 está en el centro 10- Medir el índice de refracción de la muestra y reportar
(Buentello, 2007). Se obtienen entonces expresiones como: los datos de la temperatura.
[ ( )( )]
2
G 2, 1 τ A , B ∙ G2 ,1
ln γ A =X 2B τ B , A +
X A + X B ∙ G2 ,1 (X B+ X A ∙ G1 ,2)2
[ ( )( )]
2
2 G 1, 2 τ B , A ∙ G2 ,1
ln γ B =X A τ A ,B +
X B + X A ∙ G1 ,2 ( X A + X B ∙G2 , 1)
2
Segunda experiencia: construcción de la curva de
calibración para la determinación de la
Informe de Fisicoquímica, B-2024, junio-noviembre, 2024
Equilibrio Líquido-Vapor: Mezclas Azeotrópicas 3
composición de la solución de 1-Butanol y Agua a través Soluciones para estimar la solubilidad mutua butanol-
del índice de refracción aguas
Muestra Masa Butanol Masa Agua
1- Tomar seis contenedores de 5 ml de vidrio con tapón,
(gr) (gr)
para preparar las soluciones de 1-Butanol/Agua descritas en
la tabla – 1 14 10
2 7 15
2- Rotular cada contenedor correctamente y tomar nota de
la temperatura del laboratorio. Tabla 3. Soluciones estándar
Soluciones estándar para la construcción de la
3- Pesar la sustancia a utilizar en la balanza analítica con
curva de calibración
cuidado de no sobrepasar las medidas ya que el 1-Butanol/
Agua son parcialmente miscibles y se podrían formar Estánd X mol Masa Masa Masa
reacciones conjugadas. ar 1- de 1- de total
Butan Butan Agua (g)
4- Utilizar una pipeta para manipular los líquidos a utilizar. ol ol (g) (g)
1 0 0 * *
5- Preparar las soluciones descritas en la tabla.
2 0,5 1,0 0,2432 1,2432
6- Medir el índice de refracción de cada una de las 3 0.5703 1,5 0.2749 1.7749
soluciones estándar. 4 0.6659 1,75 0.2135 1.9635
5 0.8040 2,5 0.1482 2.6482
7- Reportar los datos obtenidos. 6 1 * 0 1
Tabla 4. Construcción de la curva de
3. Discusión y Resultados calibración para la determinación de la
composición de la solución de 1-
3.1 Datos experimentales Butanol y Agua a través del índice de
refracción.
Í
Tabla 1. Determinación de la n
solubilidad mutua del sistema 1- M M di
Butanol (1)-Agua (2) a temperatura as a c
ambiente a s e
E
d a d
s T
Soluciones conjugadas a Temperatura ambiente t
e
Masa
e
Fracci
(°C) 1- t re Fracció
á de ( ón
b o fr n molar
Sistema Fase Índice de T (°C) n Agua ° molar
ut t a (X1)
d (g) C (X2)
refracción a a c
a )
(ƞ) n l ci
r
ol ( ó
1 Acuosa 1,3402 17,5 (g g n
Orgánica 1,3922 ) ) (
2 Acuosa 1,3406 17,5 ƞ
)
Orgánica 1,3908
1,33 17
1 0 0,3367 0 1
34 ,5
1,39 17 0,4672347 0,532765
2 1,007 0,2793
22 ,5 93 21
1,509 1,39 17 0,5739542 0,426045
3 0,2725
2 21 ,5 71 73
1,753 1,39 17 0,5842865 0,415713
4 0,3034
1 51 ,6 98 4
Tabla 2. Soluciones para estimar la 2,500 1,38 17 0,8185509 0,181449
5 0,1348
solubilidad mutua butanol-aguas 0 55 ,7 1 09
Revista Ciencia e Ingeniería. Vol. 37, No. 3, agosto-noviembre, 2016
4
Coello, Mariana; Maldonado, Alexandra, Rodríguez, Jesús
6
0,307
0
1,39 17
1 0 agua es un compuesto polar, por lo que la diferencia de
9 22 ,7 polaridad puede afectar su capacidad para mezclarse. La
polaridad de agua es mucho mayor a la del 1-butanol debido
a la distribución desigual de cargas en su molécula.
3.2 Resultados A pesar de sus diferencias de polaridad comparten
algunas características, ya que ambos son líquidos a
temperatura ambiente, son solventes, pueden formar
puentes de hidrógeno debido a los grupos funcionales OH -
en sus estructuras.
Debido a la similitud de algunas de sus propiedades, se
puede deducir que el 1-butanol y el agua son compuestas
parcialmente miscibles entre sí. Además, se ha comprobado
experimentalmente que el 1-butanol solo se puede
solubilizar en una proporción del 8% en una solución de
agua antes de saturar el sistema y crear dos fases diferentes.
Con el fin de describir el equilibrio termodinámico se
Fig 1. Curva de calibración a partir de las soluciones estándar 1- analizaron los valores obtenidos experimentalmente del
Butanol/Agua. índice de refracción, evaluando una serie de soluciones
estándar (tabla 2) en el refractómetro. Según la información
obtenida (tabla 2), se construyó la figura 1 en la cual se
puede observar el comportamiento de la curva de
Tabla 5. Resultados de la fracción de 1- calibración a partir de las soluciones estándar de 1-
butanol mediante la curva de Butanol/Agua. Mediante esto se hallaron los respectivos
calibración de las soluciones estándar valores de la composición molar, determinándose también
Sistema Fase X 1-butanol Fase X 1-butanol los valores de ( X 1 ) y ( X 2 ), a partir de los datos expresados
1 Orgánica 0,55732 Acuosa 0,0349 en la tabla.
2 Orgánica 0,52503 Acuosa 0,0374
Para poder proporcionar una descripción precisa
del equilibrio liquido – líquido y predecir coeficientes de
actividad, con el fin de saber si existe desviación de la
Tabla 6. Resultados de los coeficientes idealidad, se empleó el modelo NRTL (Non Random Two
de actividad por el modelo NRTL Liquids).
Estimación de los coeficientes de actividad por el modelo
NRTL En los resultados de los coeficientes de actividad por el
modelo NRTL (tabla 5), se pueden observar valores
Sistema XA γA XB γB
positivos (mayores a 1) en ambos componentes. Sabiendo
1 0,2536 1,2958 0,7464 1,035 que mientras el valor de γ se acerque 1, la sustancia se
2 0,1018 1,485 0,8982 1,0059 comportará más como si fuera ideal. Sin embargo,
para γ >1 y γ <1 , la sustancia muestra una desviación
positiva o negativa de la ley de Raoult. En el sistema 1-
Butanol/Agua estudiado, se puede observo una desviación
3.3 Análisis de resultados positiva respecto a la ley de Raoult, que se puede
argumentar por las interacciones intermoleculares entre los
Al realizar un estudio riguroso sobre los sistemas compuestos, lo que afecta de manera importante la
binarios con miscibilidad parcial, necesario considerar todos solubilidad de los mismos.
los factores fisicoquímicos que pueden intervenir en el
sistema. En el caso del sistema 1-Butanol/Agua con una De esta forma, se puede corroborar que el modelo
miscibilidad parcial, se deben tomar en cuenta los factores estudiado ayuda a predecir de manera efectiva y exacta las
que la afectan, tales como la naturaleza de la molécula, el 1- desviaciones de la idealidad del sistema parcialmente
butanol es un alcohol con una cadena hidrocarbonada larga miscible 1-Butanol/Agua, al analizar los coeficientes de
lo que le confiere propiedades lipofílicas, mientras que el actividad de los componentes.
Informe de Fisicoquímica, B-2024, junio-noviembre, 2024
Equilibrio Líquido-Vapor: Mezclas Azeotrópicas 5
McGraw-Hill.
[Link]
[Link]éndice
Los objetivos fueron alcanzados, se determinó
Cálculo de la fracción molar de 1-butanol de las
experimentalmente la posición de equilibrio termodinámico
soluciones preparadas en la tabla 2, para determinar la
de una serie de soluciones dadas; además, se pudo estimar curva de calibración.
la posición de equilibrio termodinámico de una serie de
soluciones conjugadas del sistema 1-Butanol-Agua a
PM 1-butanol= 74,12 g/mol
presión atmosférica usando el modelo Non Random Two
Liquids (NRTL). Sin embargo se entiende que, la PM H2O = 18 g/mol
temperatura es un factor físico clave para que dos sustancias
g
parcialmente miscibles tienden a disolverse entre sí mol=
PM
llegando a un punto medio de estabilidad, por arriba y mol A
X mol A=
debajo de dicho punto en el plano los solventes tenderán a moles totales
sedimentarse y separarse en fases.
1,007 g
mol 1−butanol = =0,0136 mol
mol
74 ,12
g
0 , 2793 g
mol H O= =0,0155 mol
2
mol
18
g
0,0136 mol
mol 1−butanol = =0,467234793
( 0,0136+ 0,0155 ) mol
[Link]
Buentello, G. (Diciembre de 2007). Análisis de modelos
termodinámicos para su aplicación en la
recuperación de proteínas mediante sistemas de
dos fases acuosas. Recuperado el 14 de julio de
2024, de Cálculo de la fracción de 1-butanol a través de la
[Link] curva de calibración, basado en los datos de la tabla 1:
68093/DocsTec_5718.pdf
Calderón, S. M., Belandria, A. L., Cárdenas, L., Araque, Y.,
& De Lima, A. (enero de 2022). Manual de
Trabajo Experimental en Fisicoquimica para
Y =Índice de refracción
Ingenieros Químicos. 2. Mérida, Mérida,
Venezuela: Universidad de Los Andes. X =Fracciónde 1−butanol
Recuperado el 22 de junio de 2024
Calderón, S. M., Belandria, J. I., Cárdenas, A. L., Márquez
Di Santis, L., & Araque, Y. (2017). Manual de
2
Trabajo Experimental en Físicoquimica para Y =−0,1146∗X +0,1674∗X +1,3345
Ingerieros Qúimicos: Experimentación y
Modelado de Fenómenos de Trandormación de
Fase, Distribución entre Fases y Reacciones
Químicas (1st ed.). Mérida, Venezuela: −b ± √ ❑
CreateSpace Independent Publishing Platform. X=
❑
Levine, I. N. (2004). Fisicoquímica (Vol. 1). Madrid:
Revista Ciencia e Ingeniería. Vol. 37, No. 3, agosto-noviembre, 2016
6
Coello, Mariana; Maldonado, Alexandra, Rodríguez, Jesús
−0,1674 ± √ ❑ Para el sistema (1) a temperatura ambiente 17,5°C:
X=
❑
1 mol
moles de C 4 H 10 O=14 g ∙ =0,1889 mol
74 ,12 g
1mol
Se sustituyen los valores de los índices de refracción moles de H 2 O=10 g ∙ =0,555 mol H 2 O
de la tabla 1 para obtener las fracciones experimentales de 18 g
1-butanol en la fase orgánica y en la fase acuosa.
0,1889 mol
X A ( C4 H 10 O )= =0,2536
(1)=¿0,55732 (2)=¿0,5250 ( 0,1889+0,555 ) mol
Para la fase acuosa se procede de la misma forma: X B ( H 2 O ) =1−X A=1−0,2536=0,7464
−0,1674 ± √ ❑ g A , B−g B , B b A ,B
X= −97 ,58 ° K
❑ τ A , B= = = =−0,04
RT RT J
8,314 ∙ 290 , 65 ° K
molK
(1)= 0,0349 (2)= 0,037
g B , A−g A , A bB , A 1312 , 2 ° K
τ B, A= = = =0,5430
RT RT J
8,314 ∙ 290 , 65° K
molK
G1 ,2=exp exp (−a A , B ∙ τ A , B )=exp exp (−0,2860 ∙−0,04038 ) =
Estimación de los coeficientes de actividad por el G2 ,1=exp exp (−a B , A ∙ τ B , A )=exp exp (−0,2860 ∙0,5430 )=0,8
modelo NRTL:
[ ( )( )]
2
G 2, 1 τ A , B ∙ G2 ,1
Tabla 7. Parámetros del modelo NRTL ln γ A =X 2B τ B , A +
X A + X B ∙ G2 ,1 (X B+ X A ∙ G1 ,2)2
Parámetro Valor
[ ( ) (
2
b A , B (° K) -97,58 2 0,8562 −0,
ln γ A =0,7464 0,5430 +
0,2536+0,7464 ∙ 0,8562 ( 0,7464
b B , A (° K) 1312,2
a A, B 0,2860
ln γ A =1,2958
Informe de Fisicoquímica, B-2024, junio-noviembre, 2024
Equilibrio Líquido-Vapor: Mezclas Azeotrópicas 7
ln γ A =1,485
[ ( )( )]
2
G 1, 2 τ B , A ∙ G2 ,1
ln γ B =X 2A τ A ,B +
X B + X A ∙ G1 ,2 (X A + X B ∙ G2 ,1)2
[ ( )( )]
2
2 G 1, 2 τ B , A ∙ G2 ,1
ln γ B =X A τ A ,B +
X B + X A ∙ G1 ,2 (X A + X B ∙ G2 ,1)2
[ ( )( )]
2
2 1,0116 0,5430∙ 0,8562
ln γ B =0,2536 −0,04038 +
0,7464+0,2536 ∙ 1,0116 ( 0,2536+ 0,7464 ∙ 0,8562 )2
[ ( )(
2
2 1,0116
ln γ B =0,1018 −0,04038 +
0,8982+ 0,1018∙ 1,0116 ( 0,10
ln γ B =1,035
ln γ B =1,0059
Para el sistema (2) a temperatura ambiente 17,5°C:
mol ∙ K
τ A , B=−0,04038
J
mol ∙ K
τ B , A =0,5430
J
G1 ,2=1,0116
G2 ,1=0,8562
[ ( )( )]
2
2 G 2, 1 τ A , B ∙ G2 ,1
ln γ A =X τ B , A
B +
X A + X B ∙ G2 ,1 (X B+ X A ∙ G1 ,2)2
[ ( )( )]
2
2 0,8562 −0,04038 ∙ 0,8562
ln γ A =0,8982 0,5430 +
0,1018+0,8982 ∙ 0,8562 ( 0,8982+0,1018 ∙ 1,0116 )2
Revista Ciencia e Ingeniería. Vol. 37, No. 3, agosto-noviembre, 2016