Para el sistema metil eter butilico(1)/diclorometano(2) en ELV aT=308ªK, determinar:
a)Los parametros de margules
b)Los coeficientes de actividad para x1=0.35
c)La presión de burbuja para el mismo x1 mediante la formulación GAMA/PHI. Conociendo:
B11= -2060 cc/mol
B22= -860 cc/mol
B12= -790 cc/mol
Solución:
Mediante los datos experimentales dados DATOS EXPERIMENTALES Por la ecuación de ELV
Se utilizan las siguientes ecuaciones: P(Kpa) x1 y1
Psat2 85.265 0 0 𝛾1=(𝑦1∗𝑃1)/(𝑥1∗𝑃𝑠𝑎
83.402 0.033 0.0141
82.202 0.0579 0.0253
80.481 0.0924 𝛾2=((1−𝑦1)∗𝑃1)/((1−𝑥1
0.0416
76.719 0.1665 0.0804
𝑥1∗𝛾1𝑃𝑠𝑎𝑡1=𝑦1∗𝑃 72.422 0.2482 0.1314
68.005 0.3322 0.1975
65.096 0.388 0.2457 𝑔^𝐸/𝑅𝑇=x1
59.651 0.5036 0.3686 2
56.833 0.5749 0.4564
53.689 0.6736 0.5882
51.62 0.7676 0.7176
50.455 0.8476 0.8238
49.926 0.9093 0.9002
49.72 0.9529 0.9502
Psat1 49.624 1 1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
-0.1
-0.2
𝐴21=𝐿𝑛γ├ 2^∞ ](x1=1)
LN(GANMA)
GANMA1
-0.3 GANMA2
Logarithmic (GANMA
Logarithmic (GANMA
Logarithmic (GANMA
-0.4
-0.2
LN(GANMA)
GANMA1
-0.3 GANMA2
Logarithmic (GANMA
Logarithmic (GANMA
Logarithmic (GANMA
-0.4
-0.5
𝐴12=𝐿𝑛γ├ 1^∞ ](
-0.6
X1
or la ecuación de ELV, CALCULAMOS: Datos calculados de forma practica con la ecuación de ELV
Ganma1 Ganma2 Ln(Ganma1) Ln(Ganma2)
(𝑦1∗𝑃1)/(𝑥1∗𝑃𝑠𝑎𝑡1𝑙) 1.0000 0.00000
0.71811 0.9973 -0.33114 -0.00274
0.72382 0.9974 -0.32321 -0.00257
−𝑦1)∗𝑃1)/((1−𝑥1)∗𝑃𝑠𝑎𝑡2) 0.73017 0.9967 -0.31448 -0.00328
0.74654 0.9927 -0.29231 -0.00731
0.77263 0.9813 -0.25795 -0.01884
0.81474 0.9584 -0.20489 -0.04244
𝑔^𝐸/𝑅𝑇=x1∗Lnγ1+x2∗Lnγ 0.83068 0.9410 -0.18551 -0.06084
2 0.87982 0.8899 -0.12803 -0.11670
0.90921 0.8524 -0.09518 -0.15976
0.94475 0.7944 -0.05684 -0.23014
0.97246 0.7357 -0.02792 -0.30699
0.98820 0.6842 -0.01187 -0.37957
0.99602 0.6443 -0.00399 -0.43961
0.99910 0.6166 -0.00090 -0.48361
1.0000 0.00000
POR LA GRAFICA Ln(Ganma) Vs. X1
los parametros de Margules son:
A12: -0.357
A21: -0.557
GANMA1
GANMA2
Logarithmic (GANMA2)
Logarithmic (GANMA2)
Logarithmic (GANMA2)
GANMA1
GANMA2
Logarithmic (GANMA2)
Logarithmic (GANMA2)
Logarithmic (GANMA2)
𝐴12=𝐿𝑛γ├ 1^∞ ](x1=0)
cuación de ELV Datos calculados por el modelo de Margules: Datos calculados con el m
GE/RT Ln(Ganma1)
-0.3570
-0.01357
Lnγ2=𝑥_1^2 (𝐴21+2(𝐴12−𝐴21)𝑥2) -0.3462
-0.02113 -0.3374
-0.03204 -0.3245
-0.05476 Lnγ1=𝑥_2^2 (𝐴12+2(𝐴21−𝐴12)𝑥1) -0.2943
-0.07819 -0.2579
-0.09641 -0.2185
-0.10921 𝑔^𝐸/𝑅𝑇=(A21∗x1+A12∗x2)x1∗x2 -0.1918
-0.12241 -0.1376
-0.12263 -0.1061
-0.11340 -0.0667
-0.09278 A12: -0.357 -0.0359
-0.06791 A21: -0.557 -0.0162
-0.04350 -0.0059
-0.02364 -0.0016
0.0000
COEFICIENTE DE ACTIVIDAD PARA X1=0.35
X1 GANMA1 GANMA2
0.3322 0.8038 0.9543
0.35 0.8107 0.9501
0.388 0.8254 0.9413
Diferencia cuadrada de GE/RT
os calculados con el modelo Margules
Ln(Ganma2) GE/RT Gamma 1 Gamma2 DELTA^2(GE/RT)
0.0000 0.0000 0.6998 1.0000 0.00000E+00
-0.0006 -0.0116 0.7074 0.9994 3.87964E-06
-0.0018 -0.0201 0.7136 0.9982 1.05294E-06
-0.0044 -0.0315 0.7229 0.9956 2.99893E-07
-0.0136 -0.0542 0.7451 0.9865 3.55897E-07
-0.0282 -0.0759 0.7727 0.9722 5.34238E-06
-0.0468 -0.0939 0.8038 0.9543 6.09632E-06
-0.0605 -0.1032 0.8254 0.9413 3.61649E-05
-0.0902 -0.1144 0.8714 0.9138 6.37047E-05
-0.1081 -0.1153 0.8994 0.8975 5.31010E-05
-0.1305 -0.1081 0.9354 0.8777 2.80038E-05
-0.1473 -0.0911 0.9648 0.8631 2.90903E-06
-0.1566 -0.0680 0.9840 0.8551 1.03963E-08
-0.1598 -0.0444 0.9941 0.8523 8.83794E-07
-0.1597 -0.0246 0.9984 0.8524 8.75761E-07
-0.1570 0.0000 1.0000 0.8547 0.00000E+00
2.0268E-04
SUMA:
ncia cuadrada de GE/RT
^2(GE/RT)
Para el sistema metil eter butilico(1)/diclorometano(2) en ELV aT=308.15ªK, deter
c)La presión de burbuja para el mismo x1 mediante la formulación GAMA/PHI. Co
B11= -2060 cc/mol
B22= -860 cc/mol
B12= -790 cc/mol 𝑦1∗∅1∗𝑃𝑏𝑢𝑟𝑏=𝑦1∗𝛾1∗𝑃𝑠𝑎𝑡1
De los datos anteriores:
𝑦2∗∅2∗𝑃𝑏𝑢𝑟𝑏=𝑦2∗𝛾2∗𝑃𝑠𝑎𝑡2
X1= 0.35
X2= 0.65
T= 308.15 ªK
Gamma1 0.81066922
Gamma2 0.95013671
Psat1 49.624 0.49624 bar
Psat2 85.265 0.85265 bar
R= 83.14 bar*cc/mol*K
T=308.15ªK, determinar:
ón GAMA/PHI. Conociendo:
Metodo Gamma/Phi 𝛿12=2𝐵12−𝐵11−𝐵22
𝑠𝑎𝑡1 δ12= 1340 cc/mol 𝐿𝑛∅1=𝑃𝑏𝑢𝑟𝑏/𝑅𝑇(𝐵11+𝑦_2^2 𝛿1
y1= 0.212890448824
y2= 0.787109551176 𝐿𝑛∅2=𝑃𝑏𝑢𝑟𝑏/𝑅𝑇 (𝐵22+𝑦_1^2 𝛿12
∗𝛾2∗𝑃𝑠𝑎𝑡2 Pburb= 0.68343133273 bar 〖𝐸𝐶 1:𝑦 〗 _1 ∅_1 𝑃𝑏𝑢𝑟𝑏−𝑥_1 𝛾
φ1= 0.967725621704 𝑃_𝑠𝑎𝑡1=0
φ2= 0.97890432137
Ec1 -6.34294618E-08 〖𝐸𝐶 2:𝑦 〗 _2 ∅_2 𝑃𝑏𝑢𝑟𝑏−𝑥_2 𝛾
Ec2 7.412883374E-08 𝑃_𝑠𝑎𝑡2=0
Ec3 0
〖𝐸𝐶 3:𝑦 〗 _1+𝑦_2−1=0
F(y1,y2,P) 1.069937194E-08
la presión de burbuja es: 68.3431 Kpa
2−𝐵11−𝐵22
𝑏𝑢𝑟𝑏/𝑅𝑇(𝐵11+𝑦_2^2 𝛿12)
𝑢𝑟𝑏/𝑅𝑇 (𝐵22+𝑦_1^2 𝛿12)
〗 _1 ∅_1 𝑃𝑏𝑢𝑟𝑏−𝑥_1 𝛾_1
0
〗 _2 ∅_2 𝑃𝑏𝑢𝑟𝑏−𝑥_2 𝛾_2
0
〗 _1+𝑦_2−1=0
Para el sistema binario 1-Propanol(1)/agua(2), basado en la ecuación de Margules
Dado t=70°C y {xi}; determine la presión y composición de burbuja
1. Datos obtenidos por tablas para los parametros de Antoine 4.Calculamos la Psat para cada component
1-Propanol(1) AGUA(2) Psat(Kpa)
A 16.1154 16.3872 Etanol(1) 32.6085597
B 3483.670 3885.700 AGUA(2) 31.2557127
C 205.807 230.170
Calculos previos para la tabla
2.Datos dados por el problema: Para el calculo del coeficiente de actividad
T(ºC) 70
T(ºK) 343
R (cal/mol-K) 1.987
3.Parametros de Wilson:
Etanol(1) AGUA(2)
vi (cc/mol) 58.2788 18.07 Para el calculo de la Presión de burbuja
a12 382.3 𝒙𝟐=𝟏−𝒙𝟏
a21 955.45
𝑷𝒃𝒖𝒓𝒃𝒖𝒋𝒂=𝒙𝟏∗𝜸𝟏∗𝑷𝒔𝒂𝒕𝟏+(𝟏−𝒙
Calculamos los parametros de interacción:
Para calcular la composición de vapor:
Λ12 0.176944 𝒚𝟏=(𝒙𝟏∗𝜸𝟏∗𝑷𝒔𝒂𝒕𝟏)/𝑷𝒃𝒖𝒓𝒃𝒖𝒋𝒂
Λ21 0.793810
a Psat para cada componente
5.Realizamos la tabla respectiva para realizar la grafica P-x-y
𝑃_(S𝑎𝑡 𝑖)=𝑒𝑥𝑝(𝐴𝑖−𝐵_𝑖/(𝑇(∘𝐶)
+𝐶𝑖)) X1 Gama(1) gama(2)
0 6.9456 1.00000
0.1 3.4376 1.03421
os para la tabla 0.2 2.2849 1.10970
del coeficiente de actividad (gama) para cada componete aplicamos 0.3 1.7476 1.21223
0.4 1.4500 1.33952
0.5 1.2688 1.49330
0.6 1.1531 1.67752
0.7 1.0787 1.89800
0.8 1.0326 2.16260
0.9 1.0077 2.48172
de la Presión de burbuja 1 1.0000 2.86901
𝒖𝒋𝒂=𝒙𝟏∗𝜸𝟏∗𝑷𝒔𝒂𝒕𝟏+(𝟏−𝒙𝟏)∗𝜸𝟐∗𝑷𝒔𝒂𝒕𝟐
a composición de vapor:
𝑷𝒔𝒂𝒕𝟏)/𝑷𝒃𝒖𝒓𝒃𝒖𝒋𝒂
X1= 0.3
Pburbuja: 64.543
tiva para realizar la grafica P-x-y
Pburbuja y1 77
31.256 0.0000
73
40.302 0.2781
69
42.649 0.3494
65
43.619 0.3920
44.033 0.4295 61
44.023 0.4699 57
PRESIÓN
43.534 0.5182 53
42.419 0.5805 49
40.456 0.6658 45
37.331 0.7922 41
32.609 1.0000
37
33
29
25
0 0.2 0.4 0.6 0.8 1
COMPOSICIÓN X-Y
Kpa
x1
y1
0.8 1 1.2
N X-Y
