UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO

FACULTAD DE ESTUDIOS SUPERIORES CUAUTITLÁN

TAREA 3
Carrera: Ingeniería Química

Asignatura: Equilibrio Químico

Alumno: Vázquez Sánchez José Antonio

Prof: Dr. Ricardo Baltazar Ayala

Grupo: 1201

Fecha: 25 de diciembre del 2020

 TAREAS 3° PARCIAL
A partir de los siguientes datos experimentales del equilibrio líquido vapor

• Propano: (CH3CH2CH3)
• Hexafluoropropano: (CF3CH2CF3)

T= 281.13 K
X1) P* Propano 637 KPa
X2) P* Hexafluoropropano 161 KPa
1° Elaboración del Diagrama del equilibrio líquido-vapor
A partir de los resultados elaboramos el diagrama del equilibrio líquido-vapor, donde la presión está en función
de la composición de propano, agregando los datos de los componentes puros, en donde la asignación de los
ejes del plano cartesiano será, para el eje X serán fracciones molares de cada fase del compuesto de propano
y el eje Y presiones.

P/ Kpa X CH3CH2CH3, liq Y CH3CH2CH3, vap

161.0 0 0
215.6 0.0286 0.2603
261.4 0.0544 0.3943
309.8 0.0861 0.4951
392.9 0.1539 0.6124
475.6 0.2554 0.6960
534.1 0.3662 0.7430
578.3 0.4922 0.7815
598.5 0.5711 0.8003
627.2 0.7222 0.8369
629.1 0.7249 0.8372
633.8 0.7643 0.8484
641.3 0.8168 0.8661
648.5 0.8944 0.9014
649.2 0.8992 0.9045
648.9 0.9213 0.919
647.1 0.9301 0.9256
645.1 0.9664 0.9581
637 1 1

Nomenclatura del grafico

Y CH3CH2CH3

X CH3CH2CH3
2°Calculo de fracciones molares en la fase líquido (X CF3CH2CF3) y vapor (Y CF3CH2CF3) de Hexafluoropropano
La suma de la fracción molar de Propano más la fracción molar del Hexafluoropropano es igual a 1 en cualquier
fase, empleando ese concepto utilizamos las siguientes formulaciones.
𝑋CH3 CH2 CH3 + 𝑋CF3 CH2 CF3 = 1; 𝑋CF3 CH2 CF3 = 1 − 𝑋CH3 CH2 CH3

𝑌CH3 CH2 CH3 + 𝑌CF3 CH2 CF3 = 1; 𝑌CF3 CH2 CF3 = 1 − 𝑌CH3 CH2 CH3

Ejemplo: Para 𝑃 = 215.6 𝑘𝑃𝑎

Compuesto P (kPa) X CH3CH2CH3 Y CH3CH2CH3

X1) CH3CH2CH3 215.6 0.0286 0.2603
𝑋CF3 CH2 CF3 = 1 − 𝑋CH3 CH2 CH3 = 1 − 0.0286 = 0.9714, 𝑋CF3 CH2 CF3 = 0.9714

𝑌CF3 CH2 CF3 = 1 − 𝑌CH3 CH2 CH3 = 1 − 0.2603 = 0.7397, 𝑌CF3 CH2 CF3 = 0.7397

Empleado las formulaciones anteriores le restamos a 1 las fracciones molares de Propano, realizándolo para
cada una de las presiones, como también para ambas fases, información recopilada en la siguiente tabla de
resultados.
Tabla de las fracciones molares en la fase líquido (X CF3CH2CF3) y vapor (Y CF3CH2CF3) de Hexafluoropropano

P/ Kpa X CH3CH2CH3, liq Y CH3CH2CH3, vap X CF3CH2CF3, liq Y CF3CH2CF3, vap

161.0 0 0 1 1
215.6 0.0286 0.2603 0.9714 0.7397
261.4 0.0544 0.3943 0.9456 0.6057
309.8 0.0861 0.4951 0.9139 0.5049
392.9 0.1539 0.6124 0.8461 0.3876
475.6 0.2554 0.6960 0.7446 0.3040
534.1 0.3662 0.7430 0.6338 0.2570
578.3 0.4922 0.7815 0.5078 0.2185
598.5 0.5711 0.8003 0.4289 0.1997
627.2 0.7222 0.8369 0.2778 0.1631
629.1 0.7249 0.8372 0.2751 0.1628
633.8 0.7643 0.8484 0.2357 0.1516
641.3 0.8168 0.8661 0.1832 0.1339
648.5 0.8944 0.9014 0.1056 0.0986
649.2 0.8992 0.9045 0.1008 0.0955
648.9 0.9213 0.919 0.0787 0.0810
647.1 0.9301 0.9256 0.0699 0.0744
645.1 0.9664 0.9581 0.0336 0.0419
637 1 1 0 0
3° Calculo de los coeficientes de actividad del Propano y el Hexafluoropropano
𝑌𝑖 𝑃
𝛾𝑖 =
𝑋𝑖 𝑃𝑖 ∗
Por lo tanto:
𝑌CH3 CH2 CH3 𝑃 𝑌CF3 CH2 CF3 𝑃
𝛾CH3 CH2 CH3 = , 𝛾CF3 CH2 CF3 =
𝑋CH3 CH2 CH3 𝑃CH3 CH2 CH3 ∗ 𝑋CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗

Ejemplo Para 𝑃 = 215.6 𝑘𝑃𝑎

Compuesto P (kPa) Xi Yi P* (kPa)

X1) CH3CH2CH3 59.4 0.0286 0.2603 637
X2) CF3CH2CF3 59.4 0.9714 0.7397 161
Hacemos la sustitución de datos y calculamos el coeficiente de actividad para las presiones mencionadas, con
las formulaciones previamente dichas.
𝑌CH3 CH2 CH3 𝑃 0.2603 ∙ 215.6 𝑘𝑃𝑎
𝛾CH3 CH2 CH3 = ∗ = = 3.0805
𝑋CH3 CH2 CH3 𝑃CH3 CH2 CH3 0.0286 ∙ 637𝑘𝑃𝑎
𝑌CF3 CH2 CF3 𝑃 0.7397 ∙ 215.6 𝑘𝑃𝑎
𝛾CF3 CH2 CF3 = ∗ = = 1.0197
𝑋CF3 CH2 CF3 𝑃CF3 CH2 CF3 0.9714 ∙ 161 𝑘𝑃𝑎

Resolvemos con la anterior metodología donde utilizamos las formulaciones anteriores para el cálculo de los
coeficientes de actividad, para los dos componentes a cada presión correspondiente, información recopilada en
las siguientes tablas de resultados.
Tabla de los coeficientes de actividad del Propano
𝑃CH3 CH2 CH3 ∗= 637 𝑘𝑃𝑎

P/ Kpa X CH3CH2CH3, liq Y CH3CH2CH3, vap 𝛾CH3 CH2 CH3

161.0 0 0 1
215.6 0.0286 0.2603 3.0805
261.4 0.0544 0.3943 2.9744
309.8 0.0861 0.4951 2.7966
392.9 0.1539 0.6124 2.4544
475.6 0.2554 0.6960 2.0347
534.1 0.3662 0.7430 1.7012
578.3 0.4922 0.7815 1.4415
598.5 0.5711 0.8003 1.3166
627.2 0.7222 0.8369 1.1410
629.1 0.7249 0.8372 1.1406
633.8 0.7643 0.8484 1.1045
641.3 0.8168 0.8661 1.0675
648.5 0.8944 0.9014 1.0260
649.2 0.8992 0.9045 1.0252
648.9 0.9213 0.919 1.0161
647.1 0.9301 0.9256 1.0109
645.1 0.9664 0.9581 1.0040
637 1 1 1
Tabla de los coeficientes de actividad del Hexafluoropropano
𝑃CF3 CH2 CF3 ∗= 161 kPa

P/ Kpa X CF3CH2CF3, liq Y CF3CH2CF3, vap 𝛾CF3 CH2 CF3

161.0 1 1 1
215.6 0.9714 0.7397 1.0197
261.4 0.9456 0.6057 1.0400
309.8 0.9139 0.5049 1.0631
392.9 0.8461 0.3876 1.1179
475.6 0.7446 0.3040 1.2061
534.1 0.6338 0.2570 1.3452
578.3 0.5078 0.2185 1.5456
598.5 0.4289 0.1997 1.7309
627.2 0.2778 0.1631 2.2872
629.1 0.2751 0.1628 2.3124
633.8 0.2357 0.1516 2.5320
641.3 0.1832 0.1339 2.9113
648.5 0.1056 0.0986 3.7609
 649.2 0.1008 0.0955 3.8203
648.9 0.0787 0.0810 4.1482
647.1 0.0699 0.0744 4.2780
645.1 0.0336 0.0419 4.9966
637 0 0 1
Nota: La obtención de los coeficientes de actividad de los dos compuestos involucrados, es requerido para los
posteriores desarrollos matemáticos.
4° Calculo del ∆𝐺𝑚𝑒𝑧 para cada una de las presiones requeridas
Empleado la siguiente formulación

∆𝐺𝑚𝑒𝑧 = ∑ 𝑋𝑖 𝑅𝑇 ln(𝛾𝑖 𝑋𝑖 )

∆𝐺𝑚𝑒𝑧 = 𝑅𝑇 ∑ 𝑋𝑖 ln(𝛾𝑖 𝑋𝑖 )
Por lo tanto:

∆𝐺𝑚𝑒𝑧 = 𝑅𝑇[𝑋CH3 CH2 CH3 ln(𝛾CH3 CH2 CH3 𝑋CH3 CH2 CH3 ) + 𝑋CF3 CH2 CF3 ln(𝛾CF3 CH2 CF3 𝑋CF3 CH2 CF3 )]

Sustituimos los valores correspondientes a la formulación anterior, para cada una de las presiones
correspondientes.
𝑅 = 8.314 𝐽/𝑚𝑜𝑙𝐾
𝑇 = 281.13 𝐾
Ejemplo Para 𝑃 = 215.6 𝑘𝑃𝑎

X CH3CH2CH3, liq X CF3CH2CF3, liq 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
0.0286 0.9714 3.0805 1.0197
∆𝐺𝑚𝑒𝑧 = 𝑅𝑇[𝑋CH3 CH2 CH3 ln(𝛾CH3 CH2 CH3 𝑋CH3 CH2 CH3 ) + 𝑋CF3 CH2 CF3 ln(𝛾CF3 CH2 CF3 𝑋CF3 CH2 CF3 )]

∆𝐺𝑚𝑒𝑧 = (8.314 𝐽/𝑚𝑜𝑙𝐾) ∙ (281.13 𝐾 )[0.0286 ∙ ln(3.0805 ∙ 0.0286) + 0.9714 ∙ ln(1.0197 ∙ 0.9714)]
∆𝐺𝑚𝑒𝑧 = (2337.3148𝐽/𝑚𝑜𝑙)[0.0286 ∙ ln(0.0881) + 0.9714 ∙ ln(0.9905)]
∆𝐺𝑚𝑒𝑧 = −183.9366 𝐽/𝑚𝑜𝑙
Resolvemos con la anterior metodología donde utilizamos las formulaciones anteriores para el cálculo del
∆𝐺𝑚𝑒𝑧 , información obtenida en la siguiente tabla de resultados.
Tabla de resultados de los cálculos correspondientes a ∆𝐺𝑚𝑒𝑧

X CH3CH2CH3, liq X CF3CH2CF3, liq 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3 ∆Gmez (J/ mol)
0 1 1 1 0
0.0286 0.9714 3.0805 1.0197 -183.9366
0.0544 0.9456 2.9744 1.0400 -268.5466
0.0861 0.9139 2.7966 1.0631 -348.2106
0.1539 0.8461 2.4544 1.1179 -460.2259
0.2554 0.7446 2.0347 1.2061 -577.9470
0.3662 0.6338 1.7012 1.3452 -641.3452
0.4922 0.5078 1.4415 1.5456 -682.4063
0.5711 0.4289 1.3166 1.7309 -679.2327
0.7222 0.2778 1.1410 2.2872 -621.2025
0.7249 0.2751 1.1406 2.3124 -613.0686
0.7643 0.2357 1.1045 2.5320 -587.0545
0.8168 0.1832 1.0675 2.9113 -530.7520
 0.8944 0.1056 1.0260 3.7609 -407.5229
0.8992 0.1008 1.0252 3.8203 -395.9158
0.9213 0.0787 1.0161 4.1482 -347.9523
0.9301 0.0699 1.0109 4.2780 -331.1059
0.9664 0.0336 1.0040 4.9966 -208.2827
1 0 1 1 0
5° Elaboración del Diagrama de ∆𝐺𝑚𝑒𝑧 con respecto a la composición de Propano CH3CH2CH3
A partir de los resultados del punto anterior, elaboramos el diagrama del ∆𝐺𝑚𝑒𝑧 vs 𝑋CH3 CH2 CH3 , ∆Gibbs de
mezclado en función de la composición de propano, agregando los datos de los componentes puros, en donde
la asignación de los ejes del plano cartesiano será, para el eje X serán fracción molar compuesto de Propano
en la fase líquido y el eje Y ∆Gibbs de mezclado para cada una de las presiones asignadas.

P/ Kpa X CH3CH2CH3, liq ∆Gmez (J/ mol)

161.0 0 0
215.6 0.0286 -183.9366
261.4 0.0544 -268.5466
309.8 0.0861 -348.2106
392.9 0.1539 -460.2259
475.6 0.2554 -577.9470
534.1 0.3662 -641.3452
578.3 0.4922 -682.4063
598.5 0.5711 -679.2327
627.2 0.7222 -621.2025
629.1 0.7249 -613.0686
633.8 0.7643 -587.0545
641.3 0.8168 -530.7520
648.5 0.8944 -407.5229
649.2 0.8992 -395.9158
648.9 0.9213 -347.9523
647.1 0.9301 -331.1059
645.1 0.9664 -208.2827
637 1 0
 ̅̅̅̅
6° Calculo del 𝐺 𝐸 /𝑅𝑇 para cada una de las temperaturas obtenidas experimentalmente

Empleado la siguiente formulación:

̅̅̅̅ 𝐺𝐸 𝑛𝑖 ̅̅̅̅
𝐺𝐸 = = 𝑅𝑇 ∑ ln(𝛾𝑖 ) → 𝐺 𝐸 = 𝑅𝑇 ∑ 𝑋 ln(𝛾 )
𝑖 𝑖
𝑛𝑇 𝑛𝑇
̅̅̅̅ ̅̅̅̅
𝐺𝐸 = 𝑋1 𝑅𝑇 ln(𝛾1 ) + 𝑋2 𝑅𝑇 ln(𝛾2 ) → 𝐺 𝐸 = 𝑅𝑇[𝑋 ln(𝛾 ) + 𝑋 ln(𝛾 )]
1 1 2 2

̅̅̅̅
𝐺𝐸
= 𝑋1 ln(𝛾1 ) + 𝑋2 ln(𝛾2 )
𝑅𝑇
Por lo tanto, reescribimos la ecuación sustituyendo los compuestos empleados durante la experimentación.
Si 𝑋1 = 𝑋CH3 CH2 CH3 , 𝑋2 = 𝑋CF3 CH2 CF3
̅̅̅̅
𝐺𝐸
= 𝑋1 ln(𝛾1 ) + 𝑋2 ln(𝛾2 )
𝑅𝑇
̅̅̅̅
𝐺𝐸
= 𝑋CH3 CH2 CH3 ln(𝛾CH3 CH2 CH3 ) + 𝑋CF3 CH2 CF3 ln(𝛾CF3 CH2 CF3 )
𝑅𝑇
Sustituimos los valores correspondientes a la formulación anterior, y realizamos el cálculo para cada presión, a
excepción de los componentes puros dado que sabemos por teoría que les corresponde un valor de Gibbs de
exceso igual con 0.
𝑅 = 8.314 𝐽/𝑚𝑜𝑙𝐾 𝑇 = 281.13 𝐾
Ejemplo Para 𝑃 = 215.6 𝑘𝑃𝑎

X CH3CH2CH3, liq X CF3CH2CF3, liq 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
0.0286 0.9714 3.0805 1.0197
̅̅̅̅
𝐺𝐸
= 𝑋CH3 CH2 CH3 ln(𝛾CH3 CH2 CH3 ) + 𝑋CF3 CH2 CF3 ln(𝛾CF3 CH2 CF3 )
𝑅𝑇
̅̅̅̅
𝐺𝐸
= 0.0286 ∙ ln(3.0805) + 0.9714 ∙ ln(3.0805) = 0.0511
𝑅𝑇
̅̅̅
𝐺𝐸
Resolvemos con la anterior metodología donde utilizamos las anteriores fórmulas para el cálculo de y
𝑅𝑇
llenamos con la información obtenida en la siguiente tabla de resultados.
̅̅̅̅
𝐺𝐸
Tabla de resultados de los cálculos correspondientes a 𝑅𝑇 del ejercicio

X CH3CH2CH3, liq X CF3CH2CF3, liq 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3 GE/RT
0 1 1 1 0
0.0286 0.9714 3.0805 1.0197 0.0511
0.0544 0.9456 2.9744 1.0400 0.0964
0.0861 0.9139 2.7966 1.0631 0.1444
0.1539 0.8461 2.4544 1.1179 0.2325
0.2554 0.7446 2.0347 1.2061 0.3209
0.3662 0.6338 1.7012 1.3452 0.3825
0.4922 0.5078 1.4415 1.5456 0.4011
0.5711 0.4289 1.3166 1.7309 0.3924
0.7222 0.2778 1.1410 2.2872 0.3251
0.7249 0.2751 1.1406 2.3124 0.3260
0.7643 0.2357 1.1045 2.5320 0.2949
0.8168 0.1832 1.0675 2.9113 0.2491
 0.8944 0.1056 1.0260 3.7609 0.1629
0.8992 0.1008 1.0252 3.8203 0.1574
0.9213 0.0787 1.0161 4.1482 0.1267
0.9301 0.0699 1.0109 4.2780 0.1117
0.9664 0.0336 1.0040 4.9966 0.0579
1 0 1 1 0
Nota: La información de esta tabla será requerida para realizar el siguiente gráfico.

7° Elaboración del Diagrama de ̅̅̅̅

𝐺𝐸 /𝑅𝑇 con respecto a la composición de propano 𝑋CH3 CH2 CH3

P/ Kpa X CH3CH2CH3, liq GE/RT

161.0 0 0
215.6 0.0286 0.0511
261.4 0.0544 0.0964
309.8 0.0861 0.1444
392.9 0.1539 0.2325
475.6 0.2554 0.3209
534.1 0.3662 0.3825
578.3 0.4922 0.4011
598.5 0.5711 0.3924
627.2 0.7222 0.3251
629.1 0.7249 0.3260
633.8 0.7643 0.2949
641.3 0.8168 0.2491
648.5 0.8944 0.1629
649.2 0.8992 0.1574
648.9 0.9213 0.1267
647.1 0.9301 0.1117
645.1 0.9664 0.0579
637 1 0
̅̅̅̅
𝐸
Elaboramos el diagrama del 𝐺 /𝑅𝑇 vs 𝑋CH3 CH2 CH3 , Gibbs de exceso en función de la composición de propano,
con los datos puros, en donde la asignación de los ejes del plano será, para el eje X serán fracción molar de
propano en la fase líquido y el eje Y, la energía de Gibbs de exceso para cada una de las presiones asignadas.
8° Calculo del producto de las fracciones molares del Propano y el Hexafluoropropano
Realizaremos la correspondiente multiplicación de las fracciones molares en fase liquida porque nos ayudaran
a facilitar los posteriores cálculos, ya que será utilizados para las fórmulas que corresponden a los modelos de
las soluciones reales.
Ejemplo Para 𝑃 = 215.6 𝑘𝑃𝑎

P (kPa) X CH3CH2CH3 X CF3CH2CF3

215.6 0.0286 0.9714
(𝑋1 ∙ 𝑋2 ) = 𝑋1 𝑋2

Por lo tanto, reescribimos la ecuación sustituyendo los compuestos empleados durante la experimentación
con sus correspondientes valores numéricos.
Si 𝑋1 = 𝑋CH3 CH2 CH3 , 𝑋2 = 𝑋CF3 CH2 CF3

𝑋1 𝑋2 = 𝑋CH3 CH2 CH3 ∙ 𝑋CF3 CH2 CF3

𝑋CH3 CH2 CH3 ∙ 𝑋CF3 CH2 CF3 = 0.0286 ∙ 0.9714 = 0.0278

Resolvemos con la anterior metodología y con los datos llenaremos la siguiente tabla resultados
Tabla de resultados correspondiente al producto de las fracciones molares del Propano y Hexafluoropropano
X CH3CH2CH3 X CH3CH2CH3 X CH3CH2CH3 X CH3CH2CH3
0 1 0
0.0286 0.9714 0.0278
0.0544 0.9456 0.0514
0.0861 0.9139 0.0787
0.1539 0.8461 0.1302
0.2554 0.7446 0.1902
0.3662 0.6338 0.2321
0.4922 0.5078 0.2499
0.5711 0.4289 0.2449
0.7222 0.2778 0.2006
0.7249 0.2751 0.1994
0.7643 0.2357 0.1801
0.8168 0.1832 0.1496
0.8944 0.1056 0.0944
0.8992 0.1008 0.0906
0.9213 0.0787 0.0725
0.9301 0.0699 0.0650
0.9664 0.0336 0.0325
1 0 0
Nota: La información de esta tabla será requerida para realizar gráficos y cálculos de las formulaciones
correspondientes a soluciones reales.
9° Determinación del parámetro del modelo de Solución real simple
Para soluciones reales simples tenemos que Si 𝑋1 = 𝑋CH3 CH2 CH3 , 𝑋2 = 𝑋CF3 CH2 CF3
̅̅̅̅
𝐺 𝐸 ̅̅̅̅
𝐺 𝐸
= 𝑋1 ∙ 𝑋2 𝐴0 ∴ = 𝑋CH3 CH2 CH3 ∙ 𝑋CF3 CH2 CF3 𝐴0
𝑅𝑇 𝑅𝑇
ln(𝛾1 ) = 𝐴0 𝑋2 2 ∴ ln(𝛾CH3 CH2 CH3 ) = 𝐴0 𝑋CF3 CH2 CF3 2

ln(𝛾2 ) = 𝐴0 𝑋1 2 ∴ ln(𝛾CF3 CH2 CF3 ) = 𝐴0 𝑋CH3 CH2 CH3 2

 ̅̅̅̅
GE
Para calcular A0 , realizaremos un gráfico de RT en función el producto de las fracciones molares del Propano y
el Hexafluoropropano, con los valores dados y con los datos de los componentes puros, en donde la asignación
de los ejes del plano cartesiano será, para el eje X la multiplicación de las fracciones mol y el eje Y la energía
de Gibbs de exceso entre la temperatura y la constante universal de los gases, realizar la serie de puntos para
cada presión.
Una vez realizado el grafico con los puntos experimentales, agregamos una línea de tendencia con el método
de regresión lineal, donde obtendremos como resultado la pendiente de esta línea recta y su ordenada al origen
que debe ser igual a 0 para este caso.

En el modelo de solución real simple tenemos que la línea de tendencia, tiene el siguiente significado físico
𝑦 = 𝑚𝑥 + 𝑏
Donde la pendiente es igual a 𝑚 = A0 y la ordenada al origen es igual a 𝑏 = 0
La ecuación de la línea recta resultante es: 𝑦 = 1.6504𝑥, si la pendiente tiene el valor de 𝑚 = 1.6504, entonces
el valor del parámetro de este modelo de Solución real simple es igual a A0 = 1.6504
10° Elaboración del grafico equilibrio liquido- vapor a través del parámetro de Solución real simple
A partir del parámetro calculado en el modelo de solución real simple se realizará el diagrama liquido-vapor,
para comprobar si el modelo teórico se ajusta los datos experimentales.
Primero debemos proponer un incremento de valor en el intervalo de 0 a 1, en este incremento de pasos el valor
asignado fue de 0.025, así hasta llegar al número 1, esta serie datos hace referencia a la fracción mol del
propano 𝑋CH3 CH2 CH3 .

La suma de la fracción molar de Propano más la fracción molar del Hexafluoropropano es igual a 1 en cualquier
fase, empleando ese concepto utilizamos la siguiente formulación para calcular la fracción mol del
Hexafluoropropano
𝑋CH3 CH2 CH3 + 𝑋CF3 CH2 CF3 = 1; 𝑋CF3 CH2 CF3 = 1 − 𝑋CH3 CH2 CH3

Ejemplo
𝑋CF3 CH2 CF3 = 1 − 𝑋CH3 CH2 CH3
 𝑋CF3 CH2 CF3 = 1 − 0.025 = 0.975

Empleado esta metodología tenemos la siguiente tabla de resultados de fracciones molares:

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3

0 1
0.0250 0.9750
0.0500 0.9500
0.0750 0.9250
0.1000 0.9000
0.1250 0.8750
0.1500 0.8500
0.1750 0.8250
0.2000 0.8000
0.2250 0.7750
0.2500 0.7500
0.2750 0.7250
0.3000 0.7000
0.3250 0.6750
0.3500 0.6500
0.3750 0.6250
0.4000 0.6000
0.4250 0.5750
0.4500 0.5500
0.4750 0.5250
0.5000 0.5000
0.5250 0.4750
0.5500 0.4500
0.5750 0.4250
0.6000 0.4000
0.6250 0.3750
0.6500 0.3500
0.6750 0.3250
0.7000 0.3000
0.7250 0.2750
0.7500 0.2500
0.7750 0.2250
0.8000 0.2000
0.8250 0.1750
0.8500 0.1500
0.8750 0.1250
0.9000 0.1000
0.9250 0.0750
0.9500 0.0500
0.9750 0.0250
1 0
Con el parámetro del modelo de solución real simple previamente calculado, se determinarán los coeficientes
de actividad para cada uno de los componentes involucrados, mediante la siguiente formulación teórica aplicada
a los datos que nos da el ejercicio.
Parámetro de solución real simple:
A0 = 1.6504
Si 𝑋1 = 𝑋CH3 CH2 CH3 , 𝑋2 = 𝑋CF3 CH2 CF3
 𝑙𝑛𝛾1 = A0 𝑋2 2

𝛾1 = exp(A0 𝑋2 2 )

𝛾CH3 CH2 CH3 = exp(A0 𝑋CF3 CH2 CF3 2 )

𝑙𝑛𝛾2 = A0 𝑋1 2

𝛾2 = exp(A0 𝑋1 2 )

𝛾CF3 CH2 CF3 = exp(A0 𝑋CH3 CH2 CH3 2 )

Ejemplo con las siguientes composiciones molares A0 = 1.6504

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3

0.0250 0.9750
𝛾CH3 CH2 CH3 = exp(A0 𝑋CF3 CH2 CF3 2 )

𝛾CH3 CH2 CH3 = exp( 1.6504 ∙ 0.97502 ) = 4.8014

𝛾CF3 CH2 CF3 = exp(A0 𝑋CH3 CH2 CH3 2 )

𝛾CF3 CH2 CF3 = exp(1.6504 ∙ 0.02502 ) = 1.0010

Seguimos la misma metodología para cada una de las fracciones molares previamente propuestas, la
información esta recopilada en la siguiente tabla de resultados. En donde se determinaron los coeficientes de
actividad para cada una de las fracciones molares de los componentes que intervienen en la solución real.

Tabla de resultados de los coeficientes de actividad para del Propano y Hexafluoropropano

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
0 1 5.2091 1
0.0250 0.9750 4.8014 1.0010
0.0500 0.9500 4.4348 1.0041
0.0750 0.9250 4.1047 1.0093
0.1000 0.9000 3.8069 1.0166
0.1250 0.8750 3.5381 1.0261
0.1500 0.8500 3.2950 1.0378
0.1750 0.8250 3.0750 1.0518
0.2000 0.8000 2.8756 1.0682
0.2250 0.7750 2.6947 1.0871
0.2500 0.7500 2.5303 1.1087
0.2750 0.7250 2.3809 1.1329
0.3000 0.7000 2.2450 1.1601
0.3250 0.6750 2.1212 1.1904
0.3500 0.6500 2.0083 1.2241
0.3750 0.6250 1.9054 1.2612
0.4000 0.6000 1.8115 1.3022
0.4250 0.5750 1.7258 1.3473
0.4500 0.5500 1.6475 1.3968
0.4750 0.5250 1.5760 1.4512
0.5000 0.5000 1.5107 1.5107
0.5250 0.4750 1.4512 1.5760
0.5500 0.4500 1.3968 1.6475
0.5750 0.4250 1.3473 1.7258
 0.6000 0.4000 1.3022 1.8115
0.6250 0.3750 1.2612 1.9054
0.6500 0.3500 1.2241 2.0083
0.6750 0.3250 1.1904 2.1212
0.7000 0.3000 1.1601 2.2450
0.7250 0.2750 1.1329 2.3809
0.7500 0.2500 1.1087 2.5303
0.7750 0.2250 1.0871 2.6947
0.8000 0.2000 1.0682 2.8756
0.8250 0.1750 1.0518 3.0750
0.8500 0.1500 1.0378 3.2950
0.8750 0.1250 1.0261 3.5381
0.9000 0.1000 1.0166 3.8069
0.9250 0.0750 1.0093 4.1047
0.9500 0.0500 1.0041 4.4348
0.9750 0.0250 1.0010 4.8014
1 0 1 5.2091
Una vez obtenidos los coeficientes de actividad para el Propano y Hexafluoropropano, se procederá a calcular
la presión mediante la siguiente expresión teórica:
Si 𝑋1 = 𝑋CH3 CH2 CH3 , 𝑋2 = 𝑋CF3 CH2 CF3

𝑃 = 𝑋1 𝛾1 𝑃1 ∗ + 𝑋2 𝛾2 𝑃2 ∗ → 𝑃 = 𝑋CH3 CH2 CH3 𝛾CH3 CH2 CH3 𝑃CH3 CH2 CH3 ∗ + 𝑋CF3 CH2 CF3 𝛾CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗

Recordando que

X1) P* Propano 637 KPa

X2) P* Hexafluoropropano 161 KPa

Entonces esto quiere decir que, si tenemos los componentes puros, deben de coincidir con sus presiones
parciales, cuando estos componentes están puros eso sería una forma de comprobar los cálculos realizados.
Ejemplo con las siguientes composiciones molares y coeficientes de actividad correspondientes a la segunda
fila de la tabla:
𝑃CH3 CH2 CH3 ∗ = 637 KPa

𝑃CF3 CH2 CF3 ∗ = 161 KPa

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
0.0250 0.9750 4.8014 1.0010
𝑃 = 𝑋CH3 CH2 CH3 𝛾CH3 CH2 CH3 𝑃CH3 CH2 CH3 + 𝑋CF3 CH2 CF3 𝛾CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗
∗

𝑃 = 0.0250 ∙ 4.8014 ∙ 637 kPa + 0.9750 ∙ 1.0010 ∙ 161 kPa = 233.5996 kPa
Empleando la misma metodología para cada una de las fracciones molares previamente propuestas haciendo
el mismo calculo esta información esta recopilada en la siguiente tabla de resultados. En donde determinamos
las presiones para cada una de las fracciones molares de los componentes que intervienen en la solución real.

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3 𝑃/ 𝐾𝑝𝑎
0 1 5.2091 1 161.0000
0.0250 0.9750 4.8014 1.0010 233.5996
0.0500 0.9500 4.4348 1.0041 294.8312
0.0750 0.9250 4.1047 1.0093 346.4142
0.1000 0.9000 3.8069 1.0166 389.8129
0.1250 0.8750 3.5381 1.0261 426.2756
 0.1500 0.8500 3.2950 1.0378 456.8670
0.1750 0.8250 3.0750 1.0518 482.4961
0.2000 0.8000 2.8756 1.0682 503.9393
0.2250 0.7750 2.6947 1.0871 521.8600
0.2500 0.7500 2.5303 1.1087 536.8256
0.2750 0.7250 2.3809 1.1329 549.3214
0.3000 0.7000 2.2450 1.1601 559.7626
0.3250 0.6750 2.1212 1.1904 568.5045
0.3500 0.6500 2.0083 1.2241 575.8509
0.3750 0.6250 1.9054 1.2612 582.0617
0.4000 0.6000 1.8115 1.3022 587.3587
0.4250 0.5750 1.7258 1.3473 591.9309
0.4500 0.5500 1.6475 1.3968 595.9391
0.4750 0.5250 1.5760 1.4512 599.5192
0.5000 0.5000 1.5107 1.5107 602.7855
0.5250 0.4750 1.4512 1.5760 605.8330
0.5500 0.4500 1.3968 1.6475 608.7395
0.5750 0.4250 1.3473 1.7258 611.5674
0.6000 0.4000 1.3022 1.8115 614.3645
0.6250 0.3750 1.2612 1.9054 617.1652
0.6500 0.3500 1.2241 2.0083 619.9908
0.6750 0.3250 1.1904 2.1212 622.8494
0.7000 0.3000 1.1601 2.2450 625.7365
0.7250 0.2750 1.1329 2.3809 628.6334
0.7500 0.2500 1.1087 2.5303 631.5070
0.7750 0.2250 1.0871 2.6947 634.3084
0.8000 0.2000 1.0682 2.8756 636.9709
0.8250 0.1750 1.0518 3.0750 639.4076
0.8500 0.1500 1.0378 3.2950 641.5091
0.8750 0.1250 1.0261 3.5381 643.1393
0.9000 0.1000 1.0166 3.8069 644.1319
0.9250 0.0750 1.0093 4.1047 644.2843
0.9500 0.0500 1.0041 4.4348 643.3523
0.9750 0.0250 1.0010 4.8014 641.0417
1 0 1 5.2091 637.0000
Una vez obtenida la presión para cada una de las fracciones mol podemos calcular la fracción mol en la fase
vapor para el Propano y Hexafluoropropano mediante la siguiente expresión
𝑋𝑖 𝛾𝑖 𝑃𝑖 ∗
𝑌𝑖 =
𝑃
𝑋CH3 CH2 CH3 𝛾CH3 CH2 CH3 𝑃CH3 CH2 CH3 ∗ 𝑋CF3 CH2 CF3 𝛾CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗
𝑌CH3 CH2 CH3 = 𝑌CF3 CH2 CF3 =
𝑃 𝑃
Ejemplo con las siguientes composiciones molares y coeficientes de actividad correspondientes a la segunda
fila de la tabla:
𝑃CH3 CH2 CH3 ∗ = 637 KPa

𝑃CF3 CH2 CF3 ∗ = 161 KPa

𝑃/ 𝐾𝑝𝑎 𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
233.5996 0.0250 0.9750 4.8014 1.0010
𝑋CH3 CH2 CH3 𝛾CH3 CH2 CH3 𝑃CH3 CH2 CH3 ∗
𝑌CH3 CH2 CH3 =
𝑃
 0.0250 ∙ 4.8014 ∙ 637 kPa
𝑌CH3 CH2 CH3 = = 0.3273
233.5996 𝑘𝑃𝑎
𝑋CF3 CH2 CF3 𝛾CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗
𝑌CF3 CH2 CF3 =
𝑃
0.9750 ∙ 1.0010 ∙ 161 kPa
𝑌CF3 CH2 CF3 = = 0.6727
233.5996 𝑘𝑃𝑎
Empleando el mismo cálculo para todas las filas, resulta la información que esta recopilada en la siguiente tabla
de resultados. En donde determinamos las fracciones molares de los componentes que intervienen en la
solución real en la fase vapor.

𝑃/ 𝐾𝑝𝑎 𝑋CH3 CH2 CH3 𝑌CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝑌CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
161 0 0 1 1.0000 5.2091 1
233.5996 0.0250 0.3273 0.9750 0.6727 4.8014 1.0010
294.8312 0.0500 0.4791 0.9500 0.5209 4.4348 1.0041
346.4142 0.0750 0.5661 0.9250 0.4339 4.1047 1.0093
389.8129 0.1000 0.6221 0.9000 0.3779 3.8069 1.0166
426.2756 0.1250 0.6609 0.8750 0.3391 3.5381 1.0261
456.8670 0.1500 0.6891 0.8500 0.3109 3.2950 1.0378
482.4961 0.1750 0.7104 0.8250 0.2896 3.0750 1.0518
503.9393 0.2000 0.7270 0.8000 0.2730 2.8756 1.0682
521.8600 0.2250 0.7401 0.7750 0.2599 2.6947 1.0871
536.8256 0.2500 0.7506 0.7500 0.2494 2.5303 1.1087
549.3214 0.2750 0.7593 0.7250 0.2407 2.3809 1.1329
559.7626 0.3000 0.7664 0.7000 0.2336 2.2450 1.1601
568.5045 0.3250 0.7724 0.6750 0.2276 2.1212 1.1904
575.8509 0.3500 0.7776 0.6500 0.2224 2.0083 1.2241
582.0617 0.3750 0.7820 0.6250 0.2180 1.9054 1.2612
587.3587 0.4000 0.7858 0.6000 0.2142 1.8115 1.3022
591.9309 0.4250 0.7893 0.5750 0.2107 1.7258 1.3473
595.9391 0.4500 0.7924 0.5500 0.2076 1.6475 1.3968
599.5192 0.4750 0.7954 0.5250 0.2046 1.5760 1.4512
602.7855 0.5000 0.7982 0.5000 0.2018 1.5107 1.5107
605.8330 0.5250 0.8011 0.4750 0.1989 1.4512 1.5760
608.7395 0.5500 0.8039 0.4500 0.1961 1.3968 1.6475
611.5674 0.5750 0.8069 0.4250 0.1931 1.3473 1.7258
614.3645 0.6000 0.8101 0.4000 0.1899 1.3022 1.8115
617.1652 0.6250 0.8136 0.3750 0.1864 1.2612 1.9054
619.9908 0.6500 0.8175 0.3500 0.1825 1.2241 2.0083
622.8494 0.6750 0.8218 0.3250 0.1782 1.1904 2.1212
625.7365 0.7000 0.8267 0.3000 0.1733 1.1601 2.2450
628.6334 0.7250 0.8323 0.2750 0.1677 1.1329 2.3809
631.5070 0.7500 0.8387 0.2500 0.1613 1.1087 2.5303
634.3084 0.7750 0.8461 0.2250 0.1539 1.0871 2.6947
636.9709 0.8000 0.8546 0.2000 0.1454 1.0682 2.8756
639.4076 0.8250 0.8645 0.1750 0.1355 1.0518 3.0750
641.5091 0.8500 0.8760 0.1500 0.1240 1.0378 3.2950
643.1393 0.8750 0.8893 0.1250 0.1107 1.0261 3.5381
644.1319 0.9000 0.9048 0.1000 0.0952 1.0166 3.8069
644.2843 0.9250 0.9231 0.0750 0.0769 1.0093 4.1047
643.3523 0.9500 0.9445 0.0500 0.0555 1.0041 4.4348
641.0417 0.9750 0.9699 0.0250 0.0301 1.0010 4.8014
637 1 1 0 0.0000 1 5.2091
 Finalmente elaboramos un diagrama del equilibrio liquido-vapor a partir del modelo de solución real simple en
donde, a partir de los resultados de la tabla anterior pondremos la presión en función de la composición de
propano, agregando los datos de los componentes puros, en donde la asignación de los ejes del plano
cartesiano será, para el eje X serán fracciones molares de cada fase del compuesto de propano y el eje Y
presiones.

Nomenclatura del grafico

Y CH3CH2CH3

X CH3CH2CH3

Nota: Los puntos hacen referencia a los datos experimentales, mientras que la línea continua que está en el
grafico viene siendo la construcción del diagrama del equilibrio a partir del parámetro calculado con el modelo
de la Solución real simple.
11° Determinación de los parámetros del modelo de Margules
Margules propone que
̅̅̅̅
𝐺𝐸
= 𝑋1 𝑋2 (𝐴12 𝑋2 + 𝐴21 𝑋1 )
𝑅𝑇
Que se puede escribir como:
̅̅̅̅
𝐺 𝐸
= 𝐴12 𝑋2 + 𝐴21 𝑋1
𝑅𝑇𝑋1 𝑋2
Si 𝑋1 = 𝑋CH3 CH2 CH3 , 𝑋2 = 𝑋CF3 CH2 CF3
̅̅̅̅
𝐺 𝐸
= 𝐴12 𝑋CF3 CH2 CF3 + 𝐴21 𝑋CH3 CH2 CH3
𝑅𝑇𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3

Donde tendremos que graficar ̅̅̅̅𝐺𝐸 /𝑅𝑇𝑋1 𝑋2 en función de la fracción molar compuesto de Propano en la fase
líquido; realizamos los cálculos correspondientes a para poder obtener los datos requeridos para la elaboración
del grafico mencionado con anterioridad.
̅̅̅̅
𝐺𝐸
Primero tenemos que calcular la expresión dada por , información contenida en la siguiente tabla, en
𝑅𝑇𝑋1 𝑋2
donde se realizó la división de Gibbs de exceso entre, el producto de las fracciones molares del Propano y
Hexafluoropropano realizado para las presiones correspondientes a excepción de los componentes que están
puros.
Ejemplo Para 𝑃 = 215.6 𝑘𝑃𝑎

𝑋CH3 CH2 CH3 ∙ 𝑋CF3 CH2 CF3 ̅̅̅̅

𝐺 𝐸

𝑅𝑇
0.0278 0.0511
Si 𝑋1 = 𝑋CH3 CH2 CH3 , 𝑋2 = 𝑋CF3 CH2 CF3
̅̅̅̅
𝐺𝐸 ̅̅̅̅
𝐺𝐸 0.0511
= = = 1.8410
𝑅𝑇𝑋1 𝑋2 𝑅𝑇𝑋CH3CH2 CH3 𝑋CF3 CH2CF3 0.0278
̅̅̅̅
𝐺𝐸
Tabla de resultados correspondiente al cociente de 𝑅𝑇
entre el producto de las fracciones molares de los
componentes involucrados

𝑋CH3 CH2 CH3 ∙ 𝑋CF3 CH2 CF3 ̅̅̅̅

𝐺 𝐸 ̅̅̅̅
𝐺 𝐸

𝑅𝑇 𝑅𝑇𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3

0 0 -
0.0278 0.0511 1.8410
0.0514 0.0964 1.8736
0.0787 0.1444 1.8357
0.1302 0.2325 1.7856
0.1902 0.3209 1.6875
0.2321 0.3825 1.6480
0.2499 0.4011 1.6046
0.2449 0.3924 1.6020
0.2006 0.3251 1.6204
0.1994 0.3260 1.6346
0.1801 0.2949 1.6370
0.1496 0.2491 1.6649
0.0944 0.1629 1.7243
0.0906 0.1574 1.7371
0.0725 0.1267 1.7476
0.0650 0.1117 1.7184
0.0325 0.0579 1.7840
0 0 -

̅̅̅̅
𝐺 𝐸
Posteriormente nos ponemos a graficar 𝑅𝑇𝑋CH3CH2 CH3 𝑋CF3CH2CF3
en función de la fracción molar compuesto de
Propano en la fase líquido 𝑋CH3 CH2 CH3 como el Modelo de Margules propone con ayuda de la siguiente tabla de
datos.

𝑋CH3 CH2 CH3 ̅̅̅̅

𝐺𝐸
𝑅𝑇𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3
0 -
0.0286 1.8410
0.0544 1.8736
0.0861 1.8357
0.1539 1.7856
0.2554 1.6875
0.3662 1.6480
0.4922 1.6046
 0.5711 1.6020
0.7222 1.6204
0.7249 1.6346
0.7643 1.6370
0.8168 1.6649
0.8944 1.7243
0.8992 1.7371
0.9213 1.7476
0.9301 1.7184
0.9664 1.7840
1 -

Una vez realizado el grafico con los puntos dados, agregamos una línea de tendencia con el método de
regresión lineal, donde obtendremos como resultado la pendiente de esta línea recta y su ordenada al origen
que debe ser diferente de cero.
En el modelo de Margules tenemos que la línea de tendencia, tiene el siguiente significado físico, en donde es
una función que si se evalúa con la composición de propano obtendremos dichos parámetros
Recodando que 𝑋1 = 𝑋CH3 CH2 CH3

Cuando 𝑋1 = 0 evaluada en la función lineal es igual a 𝐴12 lo que corresponde a la ordenada al origen que es
el valor de b.
Cuando 𝑋1 = 1 evaluada en la función lineal es igual a 𝐴21 lo que corresponde al valor de la ordenada cuando
el valor de propano cuando esta puro
Ejemplo con el siguiente grafico
Para Margules:
En el modelo de Margules tenemos que la ecuación de la línea de tendencia, que tiene siguiente como resultado:
𝑦 = 𝑚𝑥 + 𝑏
𝑦 = −0.1121𝑥 + 1.7781
Por lo tanto, el valor de los parámetros de Margules son:
𝐴12 (𝑋1 = 0) → −0.1121 𝑋1 + 1.7781 = −0.1121 ∙ 0 + 1.7781 = 0 + 1.7781
𝐴12 = 1.7781
𝐴21 = (𝑋1 = 1) → −0.1121 𝑋1 + 1.7781 = −0.1121 ∙ 1 + 1.7781 = −0.1121 + 1.7781
𝐴21 = 1.666
12° Elaboración del grafico equilibrio liquido- vapor a través de los parámetros de Margules
A partir de los parámetros calculados en el modelo Margules se realizará el diagrama liquido-vapor, para
comprobar si el modelo teórico se ajusta los datos experimentales.
Primero debemos proponer un incremento de valor en el intervalo de 0 a 1, en este incremento de pasos el valor
asignado fue de 0.025, así hasta llegar al número 1, esta serie datos hace referencia a la fracción mol del
propano 𝑋CH3 CH2 CH3 .

La suma de la fracción molar de Propano más la fracción molar del Hexafluoropropano es igual a 1 en cualquier
fase, empleando ese concepto utilizamos la siguiente formulación para calcular la fracción mol del
Hexafluoropropano
𝑋CH3 CH2 CH3 + 𝑋CF3 CH2 CF3 = 1; 𝑋CF3 CH2 CF3 = 1 − 𝑋CH3 CH2 CH3

Ejemplo
𝑋CF3 CH2 CF3 = 1 − 𝑋CH3 CH2 CH3

𝑋CF3 CH2 CF3 = 1 − 0.025 = 0.975

Empleado esta metodología tenemos la siguiente tabla de resultados de fracciones molares:

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3

0 1
0.0250 0.9750
0.0500 0.9500
0.0750 0.9250
0.1000 0.9000
0.1250 0.8750
0.1500 0.8500
0.1750 0.8250
0.2000 0.8000
0.2250 0.7750
0.2500 0.7500
0.2750 0.7250
0.3000 0.7000
0.3250 0.6750
0.3500 0.6500
0.3750 0.6250
0.4000 0.6000
0.4250 0.5750
0.4500 0.5500
 0.4750 0.5250
0.5000 0.5000
0.5250 0.4750
0.5500 0.4500
0.5750 0.4250
0.6000 0.4000
0.6250 0.3750
0.6500 0.3500
0.6750 0.3250
0.7000 0.3000
0.7250 0.2750
0.7500 0.2500
0.7750 0.2250
0.8000 0.2000
0.8250 0.1750
0.8500 0.1500
0.8750 0.1250
0.9000 0.1000
0.9250 0.0750
0.9500 0.0500
0.9750 0.0250
1 0
Con los parámetros del modelo de Margules previamente calculados, se determinarán los coeficientes de
actividad para cada uno de los componentes involucrados, mediante la siguiente formulación teórica aplicada a
los datos que nos da el ejercicio.
Parámetros del modelo de Margules:
𝐴12 = 1.7781
𝐴21 = 1.666
Si 𝑋1 = 𝑋CH3 CH2 CH3 , 𝑋2 = 𝑋CF3 CH2 CF3

𝑙𝑛𝛾1 = 𝑋2 2 [𝐴12 + 2(𝐴21 − 𝐴12 )𝑋1 ]

𝛾1 = exp{𝑋2 2 [𝐴12 + 2(𝐴21 − 𝐴12 )𝑋1 ]}

𝛾CH3 CH2 CH3 = exp{𝑋CF3 CH2 CF3 2 [𝐴12 + 2(𝐴21 − 𝐴12 )𝑋CH3 CH2 CH3 ]}

𝑙𝑛𝛾2 = 𝑋1 2 [𝐴21 + 2(𝐴12 − 𝐴21 )𝑋2 ]

𝛾2 = exp{𝑋1 2 [𝐴21 + 2(𝐴12 − 𝐴21 )𝑋2 ]}

𝛾CF3 CH2 CF3 = exp{𝑋CH3 CH2 CH3 2 [𝐴21 + 2(𝐴12 − 𝐴21 )𝑋CF3 CH2 CF3 ]}

Ejemplo con las siguientes composiciones molares

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3

0.0250 0.9750
𝐴12 = 1.7781
𝐴21 = 1.666

𝛾CH3 CH2 CH3 = exp{𝑋CF3 CH2 CF3 2 [𝐴12 + 2(𝐴21 − 𝐴12 )𝑋CH3 CH2 CH3 ]}

𝛾CH3 CH2 CH3 = exp{0.97502 [1.7781 + 2(1.666 − 1.7781)0.0250]} = 5.3923

 𝛾CF3 CH2 CF3 = exp{𝑋CH3 CH2 CH3 2 [𝐴21 + 2(𝐴12 − 𝐴21 )𝑋CF3 CH2 CF3 ]}

𝛾CF3 CH2 CF3 = exp{0.02502 [1.666 + 2(1.7781 − 1.666)0.9750]} = 1.0012

Seguimos la misma metodología para cada una de las fracciones molares previamente propuestas, la
información esta recopilada en la siguiente tabla de resultados. En donde se determinaron los coeficientes de
actividad para cada una de las fracciones molares de los componentes que intervienen en la solución real.
Tabla de resultados de los coeficientes de actividad para del Propano y Hexafluoropropano

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
0 1 5.9186 1
0.0250 0.9750 5.3923 1.0012
0.0500 0.9500 4.9264 1.0047
0.0750 0.9250 4.5132 1.0106
0.1000 0.9000 4.1458 1.0189
0.1250 0.8750 3.8187 1.0295
0.1500 0.8500 3.5268 1.0427
0.1750 0.8250 3.2658 1.0583
0.2000 0.8000 3.0322 1.0766
0.2250 0.7750 2.8226 1.0976
0.2500 0.7500 2.6344 1.1215
0.2750 0.7250 2.4650 1.1483
0.3000 0.7000 2.3125 1.1783
0.3250 0.6750 2.1748 1.2116
0.3500 0.6500 2.0505 1.2485
0.3750 0.6250 1.9381 1.2892
0.4000 0.6000 1.8364 1.3339
0.4250 0.5750 1.7443 1.3829
0.4500 0.5500 1.6609 1.4367
0.4750 0.5250 1.5852 1.4955
0.5000 0.5000 1.5166 1.5597
0.5250 0.4750 1.4544 1.6299
0.5500 0.4500 1.3981 1.7066
0.5750 0.4250 1.3470 1.7902
0.6000 0.4000 1.3008 1.8815
0.6250 0.3750 1.2590 1.9810
0.6500 0.3500 1.2214 2.0897
0.6750 0.3250 1.1875 2.2084
0.7000 0.3000 1.1571 2.3380
0.7250 0.2750 1.1300 2.4796
0.7500 0.2500 1.1059 2.6344
0.7750 0.2250 1.0846 2.8037
0.8000 0.2000 1.0660 2.9890
0.8250 0.1750 1.0500 3.1919
0.8500 0.1500 1.0364 3.4143
0.8750 0.1250 1.0250 3.6583
0.9000 0.1000 1.0159 3.9260
0.9250 0.0750 1.0089 4.2201
0.9500 0.0500 1.0039 4.5434
0.9750 0.0250 1.0010 4.8992
1 0 1. 5.2910
Una vez obtenidos los coeficientes de actividad para el Propano y Hexafluoropropano, se procederá a calcular
la presión mediante la siguiente expresión teórica:
Si 𝑋1 = 𝑋CH3 CH2 CH3 , 𝑋2 = 𝑋CF3 CH2 CF3

𝑃 = 𝑋1 𝛾1 𝑃1 ∗ + 𝑋2 𝛾2 𝑃2 ∗ → 𝑃 = 𝑋CH3 CH2 CH3 𝛾CH3 CH2 CH3 𝑃CH3 CH2 CH3 ∗ + 𝑋CF3 CH2 CF3 𝛾CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗

Recordando que

X1) P* Propano 637 KPa

X2) P* Hexafluoropropano 161 KPa
Entonces esto quiere decir que, si tenemos los componentes puros, deben de coincidir con sus presiones
parciales, cuando estos componentes están puros eso sería una forma de comprobar los cálculos realizados.
Ejemplo con las siguientes composiciones molares y coeficientes de actividad correspondientes a la segunda
fila de la tabla:
𝑃CH3 CH2 CH3 ∗ = 637 KPa

𝑃CF3 CH2 CF3 ∗ = 161 KPa

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
0.0250 0.9750 5.3923 1.0012
𝑃 = 𝑋CH3 CH2 CH3 𝛾CH3 CH2 CH3 𝑃CH3 CH2 CH3 ∗ + 𝑋CF3 CH2 CF3 𝛾CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗

𝑃 = 0.0250 ∙ 5.3923 ∙ 637 kPa + 0.9750 ∙ 1.0012 ∙ 161 kPa = 243.0329 kPa
Empleando la misma metodología para cada una de las fracciones molares previamente propuestas haciendo
el mismo calculo esta información esta recopilada en la siguiente tabla de resultados. En donde determinamos
las presiones para cada una de las fracciones molares de los componentes que intervienen en la solución real.

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3 𝑃/ 𝐾𝑝𝑎
0 1 5.9186 1 161
0.0250 0.9750 5.3923 1.0012 243.0329
0.0500 0.9500 4.9264 1.0047 310.5775
0.0750 0.9250 4.5132 1.0106 366.1193
0.1000 0.9000 4.1458 1.0189 411.7206
0.1250 0.8750 3.8187 1.0295 449.0955
0.1500 0.8500 3.5268 1.0427 479.6696
0.1750 0.8250 3.2658 1.0583 504.6306
0.2000 0.8000 3.0322 1.0766 524.9685
0.2250 0.7750 2.8226 1.0976 541.5096
0.2500 0.7500 2.6344 1.1215 554.9439
0.2750 0.7250 2.4650 1.1483 565.8481
0.3000 0.7000 2.3125 1.1783 574.7043
0.3250 0.6750 2.1748 1.2116 581.9154
0.3500 0.6500 2.0505 1.2485 587.8183
0.3750 0.6250 1.9381 1.2892 592.6942
0.4000 0.6000 1.8364 1.3339 596.7777
0.4250 0.5750 1.7443 1.3829 600.2636
0.4500 0.5500 1.6609 1.4367 603.3133
0.4750 0.5250 1.5852 1.4955 606.0597
0.5000 0.5000 1.5166 1.5597 608.6109
0.5250 0.4750 1.4544 1.6299 611.0536
0.5500 0.4500 1.3981 1.7066 613.4558
0.5750 0.4250 1.3470 1.7902 615.8686
0.6000 0.4000 1.3008 1.8815 618.3281
0.6250 0.3750 1.2590 1.9810 620.8561
0.6500 0.3500 1.2214 2.0897 623.4610
0.6750 0.3250 1.1875 2.2084 626.1382
 0.7000 0.3000 1.1571 2.3380 628.8700
0.7250 0.2750 1.1300 2.4796 631.6251
0.7500 0.2500 1.1059 2.6344 634.3582
0.7750 0.2250 1.0846 2.8037 637.0086
0.8000 0.2000 1.0660 2.9890 639.4990
0.8250 0.1750 1.0500 3.1919 641.7336
0.8500 0.1500 1.0364 3.4143 643.5960
0.8750 0.1250 1.0250 3.6583 644.9461
0.9000 0.1000 1.0159 3.9260 645.6172
0.9250 0.0750 1.0089 4.2201 645.4117
0.9500 0.0500 1.0039 4.5434 644.0971
0.9750 0.0250 1.0010 4.8992 641.3999
1 0 1. 5.2910 637

Una vez obtenida la presión para cada una de las fracciones mol podemos calcular la fracción mol en la fase
vapor para el Propano y Hexafluoropropano mediante la siguiente expresión
𝑋𝑖 𝛾𝑖 𝑃𝑖 ∗
𝑌𝑖 =
𝑃
𝑋CH3 CH2 CH3 𝛾CH3 CH2 CH3 𝑃CH3 CH2 CH3 ∗ 𝑋CF3 CH2 CF3 𝛾CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗
𝑌CH3 CH2 CH3 = 𝑌CF3 CH2 CF3 =
𝑃 𝑃
Ejemplo con las siguientes composiciones molares y coeficientes de actividad correspondientes a la segunda
fila de la tabla:
𝑃CH3 CH2 CH3 ∗ = 637 KPa

𝑃CF3 CH2 CF3 ∗ = 161 KPa

𝑃/ 𝐾𝑝𝑎 𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
243.0329 0.0250 0.9750 5.3923 1.0012
𝑋CH3 CH2 CH3 𝛾CH3 CH2 CH3 𝑃CH3 CH2 CH3 ∗
𝑌CH3 CH2 CH3 =
𝑃
0.0250 ∙ 5.3923 ∙ 637 kPa
𝑌CH3 CH2 CH3 = = 0.3533
243.0329𝑘𝑃𝑎
𝑋CF3 CH2 CF3 𝛾CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗
𝑌CF3 CH2 CF3 =
𝑃
0.9750 ∙ 1.0012 ∙ 161 kPa
𝑌CF3 CH2 CF3 = = 0.6467
243.0329 𝑘𝑃𝑎
Empleando el mismo cálculo para todas las filas, resulta la información que esta recopilada en la siguiente tabla
de resultados. En donde determinamos las fracciones molares de los componentes que intervienen en la
solución real en la fase vapor.

𝑃/ 𝐾𝑝𝑎 𝑋CH3 CH2 CH3 𝑌CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝑌CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
161.0000 0.0000 0.0000 1.0000 1.0000 5.9186 1.0000
243.0329 0.0250 0.3533 0.9750 0.6467 5.3923 1.0012
310.5775 0.0500 0.5052 0.9500 0.4948 4.9264 1.0047
366.1193 0.0750 0.5889 0.9250 0.4111 4.5132 1.0106
411.7206 0.1000 0.6414 0.9000 0.3586 4.1458 1.0189
449.0955 0.1250 0.6771 0.8750 0.3229 3.8187 1.0295
479.6696 0.1500 0.7025 0.8500 0.2975 3.5268 1.0427
504.6306 0.1750 0.7214 0.8250 0.2786 3.2658 1.0583
 524.9685 0.2000 0.7359 0.8000 0.2641 3.0322 1.0766
541.5096 0.2250 0.7471 0.7750 0.2529 2.8226 1.0976
554.9439 0.2500 0.7560 0.7500 0.2440 2.6344 1.1215
565.8481 0.2750 0.7631 0.7250 0.2369 2.4650 1.1483
574.7043 0.3000 0.7689 0.7000 0.2311 2.3125 1.1783
581.9154 0.3250 0.7737 0.6750 0.2263 2.1748 1.2116
587.8183 0.3500 0.7777 0.6500 0.2223 2.0505 1.2485
592.6942 0.3750 0.7811 0.6250 0.2189 1.9381 1.2892
596.7777 0.4000 0.7841 0.6000 0.2159 1.8364 1.3339
600.2636 0.4250 0.7867 0.5750 0.2133 1.7443 1.3829
603.3133 0.4500 0.7891 0.5500 0.2109 1.6609 1.4367
606.0597 0.4750 0.7914 0.5250 0.2086 1.5852 1.4955
608.6109 0.5000 0.7937 0.5000 0.2063 1.5166 1.5597
611.0536 0.5250 0.7960 0.4750 0.2040 1.4544 1.6299
613.4558 0.5500 0.7985 0.4500 0.2015 1.3981 1.7066
615.8686 0.5750 0.8011 0.4250 0.1989 1.3470 1.7902
618.3281 0.6000 0.8040 0.4000 0.1960 1.3008 1.8815
620.8561 0.6250 0.8074 0.3750 0.1926 1.2590 1.9810
623.4610 0.6500 0.8111 0.3500 0.1889 1.2214 2.0897
626.1382 0.6750 0.8154 0.3250 0.1846 1.1875 2.2084
628.8700 0.7000 0.8204 0.3000 0.1796 1.1571 2.3380
631.6251 0.7250 0.8262 0.2750 0.1738 1.1300 2.4796
634.3582 0.7500 0.8328 0.2500 0.1672 1.1059 2.6344
637.0086 0.7750 0.8406 0.2250 0.1594 1.0846 2.8037
639.4990 0.8000 0.8495 0.2000 0.1505 1.0660 2.9890
641.7336 0.8250 0.8599 0.1750 0.1401 1.0500 3.1919
643.5960 0.8500 0.8719 0.1500 0.1281 1.0364 3.4143
644.9461 0.8750 0.8858 0.1250 0.1142 1.0250 3.6583
645.6172 0.9000 0.9021 0.1000 0.0979 1.0159 3.9260
645.4117 0.9250 0.9210 0.0750 0.0790 1.0089 4.2201
644.0971 0.9500 0.9432 0.0500 0.0568 1.0039 4.5434
641.3999 0.9750 0.9693 0.0250 0.0307 1.0010 4.8992
637.0000 1.0000 1.0000 0.0000 0.0000 1.0000 5.2910
Finalmente elaboramos un diagrama del equilibrio liquido-vapor a partir del modelo de Margules en donde a
partir de los resultados de la tabla anterior pondremos la presión en función de la composición de propano,
agregando los datos de los componentes puros, en donde la asignación de los ejes del plano cartesiano será,
para el eje X serán fracciones molares de cada fase del compuesto de propano y el eje Y presiones.

Nomenclatura del grafico

Y CH3CH2CH3

X CH3CH2CH3
Nota: Los puntos hacen referencia a los datos experimentales, mientras que la línea continua que está en el
grafico viene siendo la construcción del diagrama del equilibrio a partir de los parámetros calculados con el
modelo de Margules.
13° Determinación de los parámetros del modelo de Van Laar
Van Laar Propone que
𝑋1 𝑋2 𝑋1 𝑋2
= +
̅̅̅̅
𝐺 /𝑅𝑇 𝐴′21 𝐴′12
𝐸

𝑋1 𝑋2
Primero tenemos que calcular la expresión dada por ̅̅̅̅ , información contenida en la siguiente tabla, en
𝐺 𝐸 /𝑅𝑇
donde se realizó el reciproco de Gibbs de exceso entre RT, ya que con anterioridad contábamos con ese dato
solo obtuvimos el reciproco de esa expresión mencionada sin contar los componentes puros.
Si 𝑋1 = 𝑋CH3 CH2 CH3 , 𝑋2 = 𝑋CF3 CH2 CF3

̅̅̅̅ −1 ̅̅̅̅ −1
𝐺𝐸 𝑅𝑇 𝑋1 𝑋2 𝐺𝐸 𝑅𝑇 𝑋CH3 CH2CH3 𝑋CF3 CH2CF3
(𝑅𝑇𝑋 ) = ̅̅̅̅ (𝑅𝑇𝑋 ) = ̅̅̅̅
1 𝑋2 𝐺𝐸 CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝐺𝐸

̅̅̅̅
𝐺𝐸
Por lo tanto, solo tendremos que sacar el reciproco del valor de 𝑅𝑇𝑋 que ya fue previamente
CH3 CH2 CH3 𝑋CF3 CH2 CF3
obtenido en una tabla de resultados.
Ejemplo Para 𝑃 = 215.6 𝑘𝑃𝑎
̅̅̅̅
𝐺 𝐸 𝑅𝑇 𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3
= 1.8410, = 1.8410 −1 = 0.5432
𝑅𝑇𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 ̅̅̅̅
𝐺 𝐸

Información recopilada en la siguiente tabla de resultados

̅̅̅̅
𝐺𝐸 𝑅𝑇 𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3
𝑅𝑇𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 ̅̅̅̅
𝐺𝐸
- -
1.8410 0.5432
1.8736 0.5337
1.8357 0.5448
1.7856 0.5600
1.6875 0.5926
1.6480 0.6068
1.6046 0.6232
1.6020 0.6242
1.6204 0.6171
1.6346 0.6118
1.6370 0.6109
1.6649 0.6006
1.7243 0.5799
1.7371 0.5757
1.7476 0.5722
1.7184 0.5819
1.7840 0.5605
- -
𝑅𝑇 𝑋CH3CH2CH3 𝑋CF3 CH2CF3
Posteriormente nos ponemos a graficar ̅̅̅̅ en función de la fracción molar compuesto de
𝐺𝐸
propano en la fase líquido 𝑋CH3 CH2 CH3 como el Modelo de Van Laar propone con ayuda de la siguiente tabla de
datos, a excepción del componente puro.
 𝑋CH3 CH2 CH3 𝑅𝑇 𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3
̅̅̅̅
𝐺𝐸
0 -
0.0286 0.5432
0.0544 0.5337
0.0861 0.5448
0.1539 0.5600
0.2554 0.5926
0.3662 0.6068
0.4922 0.6232
0.5711 0.6242
0.7222 0.6171
0.7249 0.6118
0.7643 0.6109
0.8168 0.6006
0.8944 0.5799
0.8992 0.5757
0.9213 0.5722
0.9301 0.5819
0.9664 0.5605
1 -

Una vez realizado el grafico con los puntos dados, agregamos una línea de tendencia con el método de
regresión lineal, donde obtendremos como resultado la pendiente de esta línea recta y su ordenada al origen
que debe ser diferente de cero.
En el modelo de Van Laar tenemos que la línea de tendencia, tiene el siguiente significado físico, en donde es
una función que si se evalúa con la composición de propano obtendremos dichos parámetros
Recodando que 𝑋1 = 𝑋CH3 CH2 CH3

Cuando 𝑋1 = 0 evaluada en la función lineal es igual a 1/𝐴′12 lo que corresponde a la ordenada al origen que
es el valor de b.
Cuando 𝑋1 = 1 evaluada en la función lineal es igual a 1/𝐴′21 lo que corresponde al valor de la ordenada cuando
el valor de propano cuando esta puro
Ejemplo con el siguiente grafico
Para Van Laar: 𝑅𝑇 𝑋1 𝑋2
̅̅̅̅
𝐺𝐸

En el modelo de Margules tenemos que la ecuación de la línea de tendencia, que tiene siguiente como resultado:
𝑦 = 𝑚𝑥 + 𝑏
𝑦 = 0.0356𝑥 + 0.5644
Por lo tanto, el valor de los parámetros de Van Laar son los siguientes y tenemos que sacar el reciproco para
despejar el termino deseado:
1
(𝑋 = 0) → 0.0356 𝑋1 + 0.5644 = 0.0356 ∙ 0 + 0.5644 = 0 + 0.5644
𝐴′12 1
1 1
= 0.5644, 𝐴′12 = = 1.7717 → 𝐴′12 = 1.7717
𝐴′12 0.5644
1
= (𝑋1 = 1) → 0.0356 𝑋1 + 0.5644 = 0.0356 ∙ 1 + 0.5644 = 0.0356 + 0.5644
𝐴′21
1 1
= 0.6, 𝐴′21 = = 1.666666667 → 𝐴′21 = 1.6666
𝐴′21 0.6
14° Elaboración del grafico equilibrio liquido- vapor a través de los parámetros de Van Laar
A partir de los parámetros calculados en el modelo Margules se realizará el diagrama liquido-vapor, para
comprobar si el modelo teórico se ajusta los datos experimentales.
Primero debemos proponer un incremento de valor en el intervalo de 0 a 1, en este incremento de pasos el valor
asignado fue de 0.025, así hasta llegar al número 1, esta serie datos hace referencia a la fracción mol del
propano 𝑋CH3 CH2 CH3 .

La suma de la fracción molar de Propano más la fracción molar del Hexafluoropropano es igual a 1 en cualquier
fase, empleando ese concepto utilizamos la siguiente formulación para calcular la fracción mol del
Hexafluoropropano
𝑋CH3 CH2 CH3 + 𝑋CF3 CH2 CF3 = 1; 𝑋CF3 CH2 CF3 = 1 − 𝑋CH3 CH2 CH3

Ejemplo
𝑋CF3 CH2 CF3 = 1 − 𝑋CH3 CH2 CH3

𝑋CF3 CH2 CF3 = 1 − 0.025 = 0.975

Empleado esta metodología tenemos la siguiente tabla de resultados de fracciones molares:

 𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3
0 1
0.0250 0.9750
0.0500 0.9500
0.0750 0.9250
0.1000 0.9000
0.1250 0.8750
0.1500 0.8500
0.1750 0.8250
0.2000 0.8000
0.2250 0.7750
0.2500 0.7500
0.2750 0.7250
0.3000 0.7000
0.3250 0.6750
0.3500 0.6500
0.3750 0.6250
0.4000 0.6000
0.4250 0.5750
0.4500 0.5500
0.4750 0.5250
0.5000 0.5000
0.5250 0.4750
0.5500 0.4500
0.5750 0.4250
0.6000 0.4000
0.6250 0.3750
0.6500 0.3500
0.6750 0.3250
0.7000 0.3000
0.7250 0.2750
0.7500 0.2500
0.7750 0.2250
0.8000 0.2000
0.8250 0.1750
0.8500 0.1500
0.8750 0.1250
0.9000 0.1000
0.9250 0.0750
0.9500 0.0500
0.9750 0.0250
1 0
Con los parámetros del modelo de Van Laar previamente calculados, se determinarán los coeficientes de
actividad para cada uno de los componentes involucrados, mediante la siguiente formulación teórica aplicada a
los datos que nos da el ejercicio.
Parámetros del modelo de Van Laar:
𝐴′12 = 1.7717
𝐴′21 = 1.6666
Si 𝑋1 = 𝑋CH3 CH2 CH3 , 𝑋2 = 𝑋CF3 CH2 CF3

𝐴 𝑋 −2 𝐴 𝑋 −2
𝑙𝑛𝛾1 = 𝐴12 (1 + 𝐴12 𝑋1 ) 𝛾1 = exp {𝐴12 (1 + 𝐴12 𝑋1 ) }
21 2 21 2
 𝐴12 𝑋CH3 CH2 CH3 −2
𝛾CH3 CH2 CH3 = exp {𝐴12 (1 + ) }
𝐴21 𝑋2

𝐴 𝑋 −2 𝐴 𝑋 −2
𝑙𝑛𝛾2 = 𝐴21 (1 + 𝐴21 𝑋2 ) 𝛾2 = exp {𝐴21 (1 + 𝐴21 𝑋2 ) }
12 1 12 1

−2
𝐴21 𝑋CF3 CH2 CF3
𝛾CF3 CH2 CF3 = exp {𝐴21 (1 + ) }
𝐴12 𝑋CH3 CH2 CH3

Ejemplo con las siguientes composiciones molares

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3

0.0250 0.9750
𝐴′12 = 1.7717
𝐴′21 = 1.6666

𝐴12 𝑋CH3 CH2 CH3 −2

𝛾CH3 CH2 CH3 = exp {𝐴12 (1 + ) }
𝐴21 𝑋2

1.7717 ∙ 0.0250 −2
𝛾CH3 CH2 CH3 = exp {1.7717 (1 + ) } = 5.3603
1.6666 ∙ 0.9750
−2
𝐴21 𝑋CF3 CH2 CF3
𝛾CF3 CH2 CF3 = exp {𝐴21 (1 + ) }
𝐴12 𝑋CH3 CH2 CH3

1.6666 ∙ 0.9750 −2
𝛾CF3 CH2 CF3 = exp {1.6666 (1 + ) } = 1.0012
1.7717 ∙ 0.0250

Seguimos la misma metodología para cada una de las fracciones molares previamente propuestas, la
información esta recopilada en la siguiente tabla de resultados. En donde se determinaron los coeficientes de
actividad para cada una de las fracciones molares de los componentes que intervienen en la solución real,
cuando hay componentes puro el coeficiente debe ser por teoría 1 solo que en este caso hay división entre cero
y provoca una indeterminación por la formula, sin embargo, en ese caso colocaremos 1 para el coeficiente este
puro, para poder seguir con los cálculos adecuadamente.
Tabla de resultados de los coeficientes de actividad para del Propano y Hexafluoropropano

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
0 1 5.8814 1
0.0250 0.9750 5.3603 1.0012
0.0500 0.9500 4.8989 1.0047
0.0750 0.9250 4.4895 1.0106
0.1000 0.9000 4.1256 1.0188
0.1250 0.8750 3.8015 1.0294
0.1500 0.8500 3.5122 1.0425
0.1750 0.8250 3.2535 1.0581
0.2000 0.8000 3.0219 1.0762
0.2250 0.7750 2.8140 1.0971
0.2500 0.7500 2.6272 1.1209
0.2750 0.7250 2.4591 1.1475
0.3000 0.7000 2.3077 1.1774
0.3250 0.6750 2.1710 1.2105
0.3500 0.6500 2.0474 1.2472
0.3750 0.6250 1.9357 1.2876
 0.4000 0.6000 1.8346 1.3320
0.4250 0.5750 1.7430 1.3808
0.4500 0.5500 1.6600 1.4343
0.4750 0.5250 1.5846 1.4927
0.5000 0.5000 1.5163 1.5566
0.5250 0.4750 1.4543 1.6265
0.5500 0.4500 1.3981 1.7027
0.5750 0.4250 1.3472 1.7859
0.6000 0.4000 1.3011 1.8767
0.6250 0.3750 1.2594 1.9759
0.6500 0.3500 1.2218 2.0841
0.6750 0.3250 1.1879 2.2024
0.7000 0.3000 1.1575 2.3316
0.7250 0.2750 1.1304 2.4728
0.7500 0.2500 1.1062 2.6273
0.7750 0.2250 1.0849 2.7963
0.8000 0.2000 1.0663 2.9815
0.8250 0.1750 1.0502 3.1845
0.8500 0.1500 1.0366 3.4072
0.8750 0.1250 1.0252 3.6517
0.9000 0.1000 1.0160 3.9203
0.9250 0.0750 1.0089 4.2157
0.9500 0.0500 1.0040 4.5409
0.9750 0.0250 1.0010 4.8993
1 0 1 5.2945
Una vez obtenidos los coeficientes de actividad para el Propano y Hexafluoropropano, se procederá a calcular
la presión mediante la siguiente expresión teórica:
Si 𝑋1 = 𝑋CH3 CH2 CH3 , 𝑋2 = 𝑋CF3 CH2 CF3

𝑃 = 𝑋1 𝛾1 𝑃1 ∗ + 𝑋2 𝛾2 𝑃2 ∗ → 𝑃 = 𝑋CH3 CH2 CH3 𝛾CH3 CH2 CH3 𝑃CH3 CH2 CH3 ∗ + 𝑋CF3 CH2 CF3 𝛾CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗

Recordando que
X1) P* Propano 637 KPa
X2) P* Hexafluoropropano 161 KPa
Entonces esto quiere decir que, si tenemos los componentes puros, deben de coincidir con sus presiones
parciales, cuando estos componentes están puros eso sería una forma de comprobar los cálculos realizados.
Ejemplo con las siguientes composiciones molares y coeficientes de actividad correspondientes a la segunda
fila de la tabla:
𝑃CH3 CH2 CH3 ∗ = 637 KPa

𝑃CF3 CH2 CF3 ∗ = 161 KPa

𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
0.0250 0.9750 5.3603 1.0012
𝑃 = 𝑋CH3 CH2 CH3 𝛾CH3 CH2 CH3 𝑃CH3 CH2 CH3 + 𝑋CF3 CH2 CF3 𝛾CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗
∗

𝑃 = 0.0250 ∙ 5.3603 ∙ 637 kPa + 0.9750 ∙ 1.0012 ∙ 161 kPa = 242.5214 kPa
Empleando la misma metodología para cada una de las fracciones molares previamente propuestas haciendo
el mismo calculo esta información esta recopilada en la siguiente tabla de resultados. En donde determinamos
las presiones para cada una de las fracciones molares de los componentes que intervienen en la solución real.
 𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3 𝑃/ 𝐾𝑝𝑎
0 1 5.8814 1 161.0000
0.0250 0.9750 5.3603 1.0012 242.5214
0.0500 0.9500 4.8989 1.0047 309.6964
0.0750 0.9250 4.4895 1.0106 364.9836
0.1000 0.9000 4.1256 1.0188 410.4227
0.1250 0.8750 3.8015 1.0294 447.7083
0.1500 0.8500 3.5122 1.0425 478.2502
0.1750 0.8250 3.2535 1.0581 503.2224
0.2000 0.8000 3.0219 1.0762 523.6042
0.2250 0.7750 2.8140 1.0971 540.2125
0.2500 0.7500 2.6272 1.1209 553.7301
0.2750 0.7250 2.4591 1.1475 564.7275
0.3000 0.7000 2.3077 1.1774 573.6821
0.3250 0.6750 2.1710 1.2105 580.9931
0.3500 0.6500 2.0474 1.2472 586.9943
0.3750 0.6250 1.9357 1.2876 591.9647
0.4000 0.6000 1.8346 1.3320 596.1372
0.4250 0.5750 1.7430 1.3808 599.7054
0.4500 0.5500 1.6600 1.4343 602.8301
0.4750 0.5250 1.5846 1.4927 605.6435
0.5000 0.5000 1.5163 1.5566 608.2536
0.5250 0.4750 1.4543 1.6265 610.7473
0.5500 0.4500 1.3981 1.7027 613.1928
0.5750 0.4250 1.3472 1.7859 615.6417
0.6000 0.4000 1.3011 1.8767 618.1308
0.6250 0.3750 1.2594 1.9759 620.6826
0.6500 0.3500 1.2218 2.0841 623.3065
0.6750 0.3250 1.1879 2.2024 625.9986
0.7000 0.3000 1.1575 2.3316 628.7423
0.7250 0.2750 1.1304 2.4728 631.5074
0.7500 0.2500 1.1062 2.6273 634.2495
0.7750 0.2250 1.0849 2.7963 636.9089
0.8000 0.2000 1.0663 2.9815 639.4093
0.8250 0.1750 1.0502 3.1845 641.6555
0.8500 0.1500 1.0366 3.4072 643.5314
0.8750 0.1250 1.0252 3.6517 644.8970
0.9000 0.1000 1.0160 3.9203 645.5848
0.9250 0.0750 1.0089 4.2157 645.3959
0.9500 0.0500 1.0040 4.5409 644.0952
0.9750 0.0250 1.0010 4.8993 641.4052
1 0 1 5.2945 637.0000
Una vez obtenida la presión para cada una de las fracciones mol podemos calcular la fracción mol en la fase
vapor para el Propano y Hexafluoropropano mediante la siguiente expresión
𝑋𝑖 𝛾𝑖 𝑃𝑖 ∗
𝑌𝑖 =
𝑃
𝑋CH3 CH2 CH3 𝛾CH3 CH2 CH3 𝑃CH3 CH2 CH3 ∗ 𝑋CF3 CH2 CF3 𝛾CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗
𝑌CH3 CH2 CH3 = 𝑌CF3 CH2 CF3 =
𝑃 𝑃
Ejemplo con las siguientes composiciones molares y coeficientes de actividad correspondientes a la segunda
fila de la tabla:
𝑃CH3 CH2 CH3 ∗ = 637 KPa
𝑃CF3 CH2 CF3 ∗ = 161 KPa

𝑃/ 𝐾𝑝𝑎 𝑋CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
242.5214 0.0250 0.9750 5.3603 1.0012
𝑋CH3 CH2 CH3 𝛾CH3 CH2 CH3 𝑃CH3 CH2 CH3 ∗
𝑌CH3 CH2 CH3 =
𝑃
0.0250 ∙ 5.3603 ∙ 637 kPa
𝑌CH3 CH2 CH3 = = 0.3520
242.5214 𝑘𝑃𝑎
𝑋CF3 CH2 CF3 𝛾CF3 CH2 CF3 𝑃CF3 CH2 CF3 ∗
𝑌CF3 CH2 CF3 =
𝑃
0.9750 ∙ 1.0012 ∙ 161 kPa
𝑌CF3 CH2 CF3 = = 0.6480
242.5214 𝑘𝑃𝑎
Empleando el mismo cálculo para todas las filas, resulta la información que esta recopilada en la siguiente tabla
de resultados. En donde determinamos las fracciones molares de los componentes que intervienen en la
solución real en la fase vapor.

𝑃/ 𝐾𝑝𝑎 𝑋CH3 CH2 CH3 𝑌CH3 CH2 CH3 𝑋CF3 CH2 CF3 𝑌CF3 CH2 CF3 𝛾CH3 CH2 CH3 𝛾CF3 CH2 CF3
161.0000 0.0000 0.0000 1.0000 1.0000 5.8814 1
242.5214 0.0250 0.3520 0.9750 0.6480 5.3603 1.0012
309.6964 0.0500 0.5038 0.9500 0.4962 4.8989 1.0047
364.9836 0.0750 0.5877 0.9250 0.4123 4.4895 1.0106
410.4227 0.1000 0.6403 0.9000 0.3597 4.1256 1.0188
447.7083 0.1250 0.6761 0.8750 0.3239 3.8015 1.0294
478.2502 0.1500 0.7017 0.8500 0.2983 3.5122 1.0425
503.2224 0.1750 0.7207 0.8250 0.2793 3.2535 1.0581
523.6042 0.2000 0.7353 0.8000 0.2647 3.0219 1.0762
540.2125 0.2250 0.7466 0.7750 0.2534 2.8140 1.0971
553.7301 0.2500 0.7556 0.7500 0.2444 2.6272 1.1209
564.7275 0.2750 0.7628 0.7250 0.2372 2.4591 1.1475
573.6821 0.3000 0.7687 0.7000 0.2313 2.3077 1.1774
580.9931 0.3250 0.7736 0.6750 0.2264 2.1710 1.2105
586.9943 0.3500 0.7777 0.6500 0.2223 2.0474 1.2472
591.9647 0.3750 0.7811 0.6250 0.2189 1.9357 1.2876
596.1372 0.4000 0.7842 0.6000 0.2158 1.8346 1.3320
599.7054 0.4250 0.7868 0.5750 0.2132 1.7430 1.3808
602.8301 0.4500 0.7893 0.5500 0.2107 1.6600 1.4343
605.6435 0.4750 0.7917 0.5250 0.2083 1.5846 1.4927
608.2536 0.5000 0.7940 0.5000 0.2060 1.5163 1.5566
610.7473 0.5250 0.7963 0.4750 0.2037 1.4543 1.6265
613.1928 0.5500 0.7988 0.4500 0.2012 1.3981 1.7027
615.6417 0.5750 0.8015 0.4250 0.1985 1.3472 1.7859
618.1308 0.6000 0.8045 0.4000 0.1955 1.3011 1.8767
620.6826 0.6250 0.8078 0.3750 0.1922 1.2594 1.9759
623.3065 0.6500 0.8116 0.3500 0.1884 1.2218 2.0841
625.9986 0.6750 0.8159 0.3250 0.1841 1.1879 2.2024
628.7423 0.7000 0.8209 0.3000 0.1791 1.1575 2.3316
631.5074 0.7250 0.8266 0.2750 0.1734 1.1304 2.4728
634.2495 0.7500 0.8333 0.2500 0.1667 1.1062 2.6273
636.9089 0.7750 0.8410 0.2250 0.1590 1.0849 2.7963
639.4093 0.8000 0.8499 0.2000 0.1501 1.0663 2.9815
641.6555 0.8250 0.8602 0.1750 0.1398 1.0502 3.1845
643.5314 0.8500 0.8721 0.1500 0.1279 1.0366 3.4072
 644.8970 0.8750 0.8860 0.1250 0.1140 1.0252 3.6517
645.5848 0.9000 0.9022 0.1000 0.0978 1.0160 3.9203
645.3959 0.9250 0.9211 0.0750 0.0789 1.0089 4.2157
644.0952 0.9500 0.9432 0.0500 0.0568 1.0040 4.5409
641.4052 0.9750 0.9693 0.0250 0.0307 1.0010 4.8993
637.0000 1.0000 1.0000 0.0000 0.0000 1 5.2945
Finalmente elaboramos un diagrama del equilibrio liquido-vapor a partir del modelo de Van Laar en donde a
partir de los resultados de la tabla anterior pondremos la presión en función de la composición de propano,
agregando los datos de los componentes puros, en donde la asignación de los ejes del plano cartesiano será,
para el eje X serán fracciones molares de cada fase del compuesto de propano y el eje Y presiones.

Nomenclatura del grafico

Y CH3CH2CH3

X CH3CH2CH3

Nota: Los puntos hacen referencia a los datos experimentales, mientras que la línea continua que está en el
grafico viene siendo la construcción del diagrama del equilibrio a partir de los parámetros calculados con el
modelo de Van Laar.