12.

3 Lo siguiente es el conjunto de informacion de EVL para el sistema de Acetona(1)/Metanol(2

P/kPa X1 Y1 X2 Y2
72.278 0.0278 0.0647 0.9722 0.9353
75.279 0.057 0.1295 0.943 0.8705
77.524 0.0858 0.1848 0.9142 0.8152
78.951 0.1046 0.219 0.8954 0.781
82.528 0.1452 0.2694 0.8548 0.7306
86.762 0.2173 0.3633 0.7827 0.6367
90.088 0.2787 0.4184 0.7213 0.5816
93.206 0.3579 0.4779 0.6421 0.5221
95.017 0.405 0.5135 0.595 0.4865
93.365 0.448 0.5512 0.552 0.4488
97.646 0.5052 0.5844 0.4948 0.4156
98.462 0.5432 0.6174 0.4568 0.3826
99.811 0.6332 0.6772 0.3668 0.3228
99.95 0.6605 0.6926 0.3395 0.3074
100.278 0.6945 0.7124 0.3055 0.2876
100.467 0.7327 0.7383 0.2673 0.2617
100.999 0.7752 0.7729 0.2248 0.2271
101.059 0.7922 0.7876 0.2078 0.2124
99.877 0.908 0.8959 0.092 0.1041
99.799 0.9448 0.9336 0.0552 0.0664

Presiones de vapor a partir de los datos de equilibrio:

Psat1= 96.885 kPa

Psat2= 68.728 kPa

Calcule los valores experimentales de los coeficientes d

i γ1 γ2 ln⁡γ1 ln⁡γ2 GE/RT

1 1.736237243 1.011737246 0.551720268 0.011668899 0.026682327
2 1.765274349 1.011107235 0.568306117 0.011046003 0.042809829 0.6
3 1.723432618 1.005831936 0.54431801 0.005814997 0.052018555
4 1.706135509 1.001977278 0.534230877 0.001975326 0.057649256
5 1.580431772 1.026319945 0.457698083 0.025979535 0.088665068 0.5

6 1.497196084 1.02691705 0.403594082 0.026561159 0.108490413

7 1.395934719 1.056919 0.333564241 0.055358072 0.132894131
0.4
8 1.284584445 1.102709682 0.250435277 0.097770498 0.152409222
lnY1, lnY2, GERT

9 1.24345535 1.130403483 0.217894077 0.122574634 0.161179008

0.3

0.2
 0.4

lnY1, lnY2, GERT

10 1.185656135 1.104496047 0.170296322 0.099389165 0.131155571 0.3
11 1.165855643 1.193346598 0.153455275 0.176761627 0.164987258
12 1.155098376 1.199924215 0.144185515 0.1822584 0.161577209
0.2
13 1.101787666 1.278053105 0.096934012 0.245337908 0.151368561
14 1.081772458 1.316779865 0.078600861 0.27518926 0.145342622
15 1.061697466 1.373566293 0.05986901 0.317410491 0.138547932 0.1
16 1.044897205 1.431180769 0.043918512 0.358499816 0.128006095
17 1.039369751 1.484582015 0.038614521 0.395133262 0.118759934
18 1.037025229 1.502969819 0.036356258 0.40744303 0.113468089 0
1 2
19 1.017144449 1.6443516 0.016999142 0.497346142 0.061191066
20 1.017865992 1.746712729 0.017708271 0.557735581 0.047517779

a) Ajustar los datos GE / RT a la ecuación de Margules por mínimos cuadrados lineales:

VXi= X1i Vyi= GERT/X1i*X2i

Bajada: -0.018
Intercepto: 0.708

A12= Intercepto A21= Bajada + A12

A12= 0.708 A21= 0.690

De las siguientes ecuaciones obtenemos los valores calculados:

γ1 (X1, X2)=exp [X2^2*{A12+2(A21-A12)*X1}]

γ2 (X1, X2)=exp [X1^2*{A21+2(A12-A21)*X2}]

j= 1…101 X1j= 0.01j - 0.1 X2j= 1 - X1j

Pcalc j= X1j*γ1(X1j, X2j)*Psat1 + X2j*γ2(X1j, X2j)*Psat2

Y1 calc j= (𝑋1𝑗∗γ1(𝑋1𝑗,
𝑋2𝑗)∗𝑃𝑠𝑎𝑡1)/(𝑃𝑐𝑎𝑙𝑐 𝑗)

105
Diagrama: ecuación de margules ajustada a datos GE / RT:
100

Pcalc= X1i*γ1(X1i, X2i)*Psat1 + X2i*γ2(X1i, X2i)*Psat2 95

90
Y1calc=(𝑋1𝑖∗γ1(𝑋1𝑖,
𝑋2𝑖)∗𝑃𝑠𝑎𝑡1)/𝑃𝑐𝑎𝑙𝑐
85
Pi (kPa)

80

75
 95

90

85

Pi (kPa)
80
Desviacion de RMS en P:
√(Σ 〖(𝑃𝑖−𝑃𝑐𝑎𝑙𝑐) 〗^2/𝑛) 75
RMS=
70

65
RMS= 0.851 kPa
60
0 0.1 0.2 0

b) Ajustando los datos de GE / RT a la ecuación de Van Laar por mínimos cuadrados lineales:

Vxi= X1 VYi=(𝑋1𝑖∗𝑋2𝑖)/
𝐺𝐸𝑅𝑇𝑖

Bajada= bajada(VX, VY) Intercepto= intercepto(VX, VY)

Bajada= 0.015 Intercepto= 1.442

A12=1/𝐼𝑛𝑡𝑒𝑟𝑐𝑒𝑝𝑡𝑜 A21= 1/((𝐵𝑎𝑗𝑎𝑑𝑎+𝐼𝑛𝑡𝑒𝑟𝑐𝑒𝑝𝑡𝑜))

A12= 0.693 A21= 0.686

exp⁡[𝐴12∗(1+(𝐴12∗𝑋1)/(𝐴21∗𝑋2))^(−2)]
γ1(X1, X2)=

105
exp⁡[𝐴21∗(1+(𝐴21∗𝑋2)/(𝐴12∗𝑋1))^(−2)]
γ2(X1, X2)=
100

95
𝑃_𝑐𝑎𝑙𝑐=𝑋1∗γ1(X1, X2)∗ 〖
𝑃𝑠𝑎𝑡 〗_1+𝑋2∗γ2(X1, X2)∗ 〖
𝑃𝑠𝑎𝑡 〗_2
90

𝑌1𝑐𝑎𝑙𝑐=(𝑋1∗γ1(𝑋1, 𝑋2)∗ 〖
𝑃𝑠𝑎𝑡 〗_1)/𝑃𝑐𝑎𝑙𝑐 85
Pi (kPa)

80

Diagrama: ecuación de van Laar ajustada a datos GE / RT: 75

70
Desviacion RMS en P:
√(Σ 〖(𝑃𝑖−𝑃𝑐𝑎𝑙𝑐) 〗^2/𝑛) 65
RMS=
60
0 0.1 0.2 0
 70

√(Σ 〖(𝑃𝑖−𝑃𝑐𝑎𝑙𝑐) 〗^2/𝑛) 65

60
0 0.1 0.2 0

RMS= 0.701 kPa

C) Ajustando los datos GE / RT a la ecuación de Wilson por mínimos cuadrados no lineales. Minim

Suponiendo:

A12= 0.5
A21= 1.0

𝑆𝑆𝐸(𝐴12,𝐴21)=∑▒ 〖
〖[𝐺𝐸𝑅𝑇+𝑋1∗ln⁡(𝑋1+𝑋2∗𝐴12)+𝑋2∗ln⁡(𝑋2+𝑋1∗𝐴21)] 〗^2 〗

A12= 0.71
A21= 0.681

γ1(X1, X2)=(𝑒𝑥𝑝[𝑋2∗(𝐴12/(𝑋1+(𝑋2∗𝐴12)
)−𝐴21/(𝑋2+(𝑋1∗𝐴21) ))])/(𝑋1+𝑋2∗𝐴12)

γ2(X1, X2)=(𝑒𝑥𝑝[−𝑋1∗(𝐴12/(𝑋1+
(𝑋2∗𝐴12) )−𝐴21/(𝑋2+(𝑋1∗𝐴21)
))])/(𝑋2+𝑋1∗𝐴21)

𝑃_𝑐𝑎𝑙𝑐𝑗=𝑋1𝑗∗γ1(X1,j X2𝑗)∗ 〖
𝑃𝑠𝑎𝑡 〗_1+𝑋2𝑗∗γ2(X1j, X2𝑗)∗ 〖
𝑃𝑠𝑎𝑡 〗_2

𝑃_𝑐𝑎𝑙𝑐𝑖=𝑋1𝑖∗γ1(X1i, X2𝑖)∗ 〖
𝑃𝑠𝑎𝑡 〗_1+𝑋2𝑖∗γ2(X1,i X2i)∗ 〖
𝑃𝑠𝑎𝑡 〗_2

𝑌1𝑐𝑎𝑙𝑐𝑗=(𝑋1𝑗∗γ1(𝑋1𝑗, 𝑋2𝑗)∗ 〖
𝑃𝑠𝑎𝑡 〗_1)/𝑃𝑐𝑎𝑙𝑐𝑗 𝑌1𝑐𝑎𝑙𝑐𝑖=(𝑋1𝑖∗γ1(𝑋1𝑖, 𝑋2𝑖)∗ 〖
𝑃𝑠𝑎𝑡 〗_1)/𝑃𝑐𝑎𝑙𝑐𝑖

Diagrama P-x,y: Ecuación de Wilson para datos GERT: 105

√(Σ 〖(𝑃𝑖−𝑃𝑐𝑎𝑙𝑐) 〗^2/𝑛) 100

RMS= 95

90

RMS= 0.361 kPa 85

Pi (kPa)

80

75
 95

90

85

Pi (kPa)
80

75

70

65

60
0 0.1 0.2 0.3 0.4
de Acetona(1)/Metanol(2) a 55°C (extraido de D.C. Freshwater y K.A. Pike, J. Chem. Eng Data, vol. 12, pp. 179-183, 1967):

ales de los coeficientes de actividad y el exceso de energía de Gibbs:

Chart Title
0.6

0.5

0.4
lnY1, lnY2, GERT

0.3

0.2
 0.4
lnY1, lnY2, GERT

0.3

0.2

0.1

0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20

X1i
 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

X1i, y1i

os lineales:

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

X1i, y1i
 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
X1i, y1i

ados no lineales. Minimice la suma de los errores cuadrados utilizando la función de minimización de Mathcad
.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

X1i, y1i
83, 1967):
17 18 19 20
ad
 12.6
eter metil ter-butilico/diclorometano T=308.15 K
P/KPa Xi yi Tenemos la ecuacion de margules
85.265 0 0 𝐺^𝐸/𝑅𝑇 = x1 ln γ1 + x2 ln γ2
83.402 0.033 0.0141 Realizamos una extencion a la ecuacion solo aplica para molecula
82.202 0.0579 0.0253 𝐺^𝐸/𝑅𝑇= 〖(𝐴 〗_21 𝑋_1+ 𝐴_12 𝑋_2-C𝑋_1
80.481 0.0924 0.0416 𝑋_2)𝑋_1 𝑋_2
76.719 0.1665 0.0804 𝑋_2^2
ln(𝐴_12+2
γ1= 〖(𝐴 〗_21− 𝐴_12-C 〖)𝑋 〗_1+3𝐶𝑋_1^2)𝑋_1^2
72.422 0.2482 0.1314 𝑋_1^2
ln (𝐴_21+2
γ2= 〖(𝐴 〗_12− 𝐴_21-C 〖)𝑋 〗_2+3𝐶𝑋_1^2)𝑋_2^2
68.005 0.3322 0.1975 𝑥_2= 〖1−𝑥 〗_2 𝑦= 〖1−𝑦 〗_2
65.096 0.388 0.2457 𝑃_𝑠𝑎𝑡1=𝑦_(1(𝑃_𝑡))/𝑥_1 =35.6354KPa 𝑃_𝑠𝑎𝑡2=𝑦_(2(𝑃_𝑡))/𝑥_2 =85.0320
59.651 0.5036 0.3686
56.833 0.5749 0.4564 Psat1= 35.6354 Psat2=
53.689 0.6736 0.5882 Ύ_1 (𝑦_1 Ύ_2 (𝑦_2 (𝑃))/(𝑥_2 (𝑃_𝑠𝑎𝑡2
(𝑃))/(𝑥_1 𝐺^𝐸/𝑅𝑇= x1 ln γ1 + x2 ln
51.62 0.7676 0.7176 (𝑃_𝑠𝑎𝑡1))=1
γ2 =0
50.455 0.8476 0.8238
49.926 0.9093 0.9002
𝐴_12=-0.3 𝐴_21=−0.5
49.72 0.9529 0.9502
49.624 1 1

1.2

0.8

0.6

0.4

0.2

0
45 50 55 60 65 70 75 80 85 90
n de margules x2 y2 Ύ1
x2 ln γ2 1 1 0
o aplica para moleculas homogeneas 0.967 0.9859 1
_2-C𝑋_1 0.9421 0.9747 1.0079586208
0.9076 0.9584 1.0167942588
𝐶𝑋_1^2)𝑋_1^2 0.8335 0.9196 1.0395924247
𝐶𝑋_1^2)𝑋_2^2 0.7518 0.8686 1.0759261527
𝑦= 〖1−𝑦 〗_2 0.6678 0.8025 1.1345579381
(2(𝑃_𝑡))/𝑥_2 =85.0320KPa 0.612 0.7543 1.1567675465
0.4964 0.6314 1.2251963662
85.032 0.4251 0.5436 1.2661124242
𝑦_2 (𝑃))/(𝑥_2 (𝑃_𝑠𝑎𝑡2))=1 0.3264 0.4118 1.3156083421
0.2324 0.2824 1.3542031702
0.1524 0.1762 1.3761108278
0.0907 0.0998 1.3870015675
C=0.2 0.0471 0.0498 1.3912884542
0 0 1.3925478597

0.4

0.3

0.2

0.1

0
0 0.2 0.4 0.6
-0.1

-0.2

-0.3

-0.4

-0.5

-0.6
 Ύ2 Xi lnΎ1 lnΎ2 𝐺^𝐸/𝑅𝑇
1.002740145 0 0 0.002736397534 0
1.000001068 0.033 0 1.067788935E-06 1.0325519E-06
1.000170296 0.0579 0.007927118 0.000170281724 0.000619402545
0.999455098 0.0924 0.016654795 -0.000545050054 0.001044215585
0.995437272 0.1665 0.038828737 -0.004573169454 0.002653247976
0.984023848 0.2482 0.073181828 -0.016105146624 0.006055880504
0.961074551 0.3322 0.126243093 -0.039703296852 0.015424093935
0.943549295 0.388 0.145629517 -0.05810666849 0.02094297167
0.892294342 0.5036 0.20310113 -0.113959220575 0.045712372034
0.854685883 0.5749 0.235951122 -0.157021266383 0.068898559977
0.796597819 0.6736 0.274299176 -0.227405345909 0.110542820184
0.737673487 0.7676 0.303213215 -0.304253981299 0.16203783865
0.68602944 0.8476 0.31926128 -0.376834736065 0.213176246909
0.646052234 0.9093 0.327144271 -0.436874920889 0.257847730717
0.618240087 0.9529 0.330230263 -0.480878406581 0.292027044883
0 1 0.331135062 0 0.331135061861

0.2 0.4 0.6 0.8 1 1.2

12.12 Prepare un diagrama P-x-y a T=60 ºC, para uno de los sistemas binarios listados en la tabla 12.5, con base en
la ecuación (10.5) y en la ecuacion de Wilson

𝑃=𝑋1ɣ1𝑃1𝑠𝑎𝑡+𝑋2ɣ2P2sat 𝑌𝑖𝑃=𝑋𝑖ɣ𝑖𝑃𝑖𝑠𝑎𝑡

X1
DATOS 0
T= 60 ʌ12= 0.1572 0.05
a12= 291.27 ʌ21= 0.4598 0.1
a21= 1448.01 b12= 500.4 P= 70 0.15
V1= 74.05 b21 1636.57 0.2
V2= 18.07 α= 0.5081 0.25
ESPECIE A B C Psat 0.3
1 14.3916 2795.82 230 115.68157568 0.35
2 16.262 3799.89 226.35 19.923719331 0.4
0.45
0.5
0.55
0.6
0.65
0.7
0.75
0.8
0.85
0.9
0.95
1
X2 ɣ1 ɣ2 Peq(kpa) Y1 Y2
1 0 0 0 0 0
0.95 8.108464935 1.220444289 70 0.67 0.33
0.9 4.477808994 1.014982287 70 0.74 0.26
0.85 3.106227861 0.950684767 70 0.77 0.23
0.8 2.390181828 0.922267559 70 0.79 0.21
0.75 1.936349835 0.936906727 70 0.8 0.2
0.7 1.633795173 0.953637204 70 0.81 0.19
0.65 1.414226933 0.983752063 70 0.818 0.182
0.6 9.566666667 1.054020068 455.2748296 0.82 0.18
0.55 1.11609053 1.08596007 70 0.83 0.17
0.5 1.011742789 1.152395274 70 0.836 0.164
0.45 0.924166967 1.249208969 70 0.84 0.16
0.4 0.857238208 1.317525085 70 0.85 0.15
0.35 0.799675244 1.415398377 70 0.859 0.141
0.3 0.743420026 1.639586772 70 0.86 0.14
0.25 0.70999494 1.686432108 70 0.88 0.12
0.2 0.673184122 1.932370124 70 0.89 0.11
0.15 0.647822923 2.108040136 70 0.91 0.09
0.1 0 2.459380158 4.9 0.93 0.07
0.05 0.636957183 0 70 1 0
0 0 0 0
12.14 Elabore un diagrama P-x-y a t=60°C, para uno de los sistemas binarios listados en la tabla 12.5, con base
en la ecuación (10.5) y en la ecuación NRTL.

DATOS ln⁡〖
𝑃𝑠𝑎𝑡 /𝑘𝑃𝑎=𝐴−𝐵/
1-Propanol-Agua (𝑇+𝐶) 〗

Valores de la Tabla 12.5 Ec. NRTL Psat1= 15547.80412063 kPa

b12 500.4 cal mol^-1 Psat2= 13218.12707084 kPa
b21 1636.57 cal mol^-1
a 0.5343
𝜏12=𝑏12/𝑅𝑇 𝜏21=𝑏21/𝑅𝑇
Tabla B.2 Ctes. De Antoine
1.Propanol Agua 0.1806622702 0.590860214879
A 16.1154 16.3872
B 3483.67 3885.7 𝐺12=exp⁡(−∗τ12) 0.907984609
C 205.807 230.17 𝐺21=exp⁡(−∗τ21) 0.729280671

T= 333.15 K EC.12.21a
60 𝑙𝑛𝛾1= 〖𝑥 2 〗 ^2
R= 8.314 [𝜏21(𝐺21/(𝑥1+𝑥2𝐺21))^2+𝐺12𝜏12/(𝑥2+𝑥1𝐺12)^2 ]

EC.12.21b

𝑙𝑛𝛾2= 〖𝑥 1 〗 ^2
[𝜏12(𝐺12/(𝑥2+𝑥1𝐺12))^2+𝐺21𝜏21/(𝑥1+𝑥2𝐺21)^2 ]

EC.10.5

𝑦1𝑃=𝑥1𝛾1𝑃𝑖𝑠𝑎𝑡

P= x1*P1sat*exp(g1)+x2*P2sat*exp(g2)

Y= XiPisatExp(lng)/R
abla 12.5, con base

Suponemos valores
x1 x2 g1 g2 P y1
0 1 0.754898776 0 13218.12707 0
0.1 0.9 0.580299976 0.009049594 14782.19048 0.187910628
0.2 0.8 0.436645542 0.034266371 15755.19974 0.305427328
0.3 0.7 0.319412904 0.073213275 16375.14686 0.392034577
0.4 0.6 0.22491682 0.123965793 16765.27019 0.464513809
0.5 0.5 0.150146115 0.185009116 16985.52503 0.531823738
0.6 0.4 0.09263633 0.255158586 17058.22591 0.599953921
0.7 0.3 0.050369688 0.333497777 16980.82537 0.67403676
0.8 0.2 0.021696003 0.41933003 16731.85372 0.759691807
0.9 0.1 0.00526979 0.512140326 16272.87804 0.864441963
1 0 0 0.611565204 15547.80412 1

18000

17000

16000
12/(𝑥2+𝑥1𝐺12)^2 ]
15000

14000
P

13000
21/(𝑥1+𝑥2𝐺21)^2 ]
12000

11000

10000
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
X-Y

x1 yeq P 35
0 0 20.007
0.1 0.373 30.48 30

0.2 0.386 30.876

25
0.3 0.395 31.048
0.4 0.414 31.172
20
0.5 0.442 31.085
15

10
 25

20

0.6 0.479 30.657 15

0.7 0.531 29.74
0.8 0.606 28.095 10
0.9 0.732 25.256
1 1 20.275 5

0
0 0.2 0.4 0.6 0.8
0.4 0.6 0.8 1 1.2
 12.16 Para uno de los sistemas binarios listados en la tabla 12.5, con base en la ecuación (10.5) y en la ecuación
Wilson, realice los siguientes cálculos:
a) BUBL P: t = 60°C, Xl = 0.3.
b) DEW P: t = 60°C, YI = 0.3.
c) P, T -vaporización instantánea: t = 60°C, P = (Pburbuja + Procío), Zj = 0.3.

INCIS
DATOS: FORMULAS:
1-PROPANOL AGUA
A(K)= 16.1154 16.3872 𝐴_21=𝑉_1/𝑉_2 𝐸𝑋𝑃 (−𝑎_21)/𝑅𝑇
𝐴_12=𝑉_2/𝑉_1 𝐸𝑋𝑃 (−𝑎_12)/𝑅𝑇
B(K)= 3483.67 3885.7
C(K)= 205.807 230.17
𝛾_1=𝑒𝑥𝑝[𝑥_2 (𝐴_12/(𝑥_1+𝑥_2 𝐴_12 )−𝐴_21/(𝑥_2+𝑥_1 𝐴_21 ))]/(𝑥_1+𝑥_2 𝐴_1
V(cmᶟ/mol)= 75.14 18.07
a₁₂(cal/mol)= 775.48
a₂₁(cal/mol)= 1351.9
𝛾_2=𝑒𝑥𝑝[−𝑥_1 (𝐴_12/(𝑥_1+𝑥_2 𝐴_12 )−𝐴_21/(𝑥_2+𝑥_1 𝐴_21 ))]/(𝑥_2+𝑥_1 𝐴

T(⁰C)= 60
T(K)= 333.15
x₁= 0.3 𝑃_𝐵𝑈𝐵𝐿=∑▒ 〖
𝑥 _𝑖 𝛾_𝑖 〖𝑃 _𝑖 〗 ^𝑆 〗
x₂= 0.7 𝑦_𝑖 𝑃=𝑥_𝑖 𝛾_1 〖𝑃 _𝑖 〗 ^𝑆 (EC.10.5)
R(cal/molK)= 1.9872

INCISO (B)
DATOS: FORMULAS:
y₁= 0.3
𝑃_𝐷𝐸𝑊=1/(∑▒𝑦_𝑖/(𝛾_𝑖 〖𝑃 _𝑖 〗 ^𝑆 ))
y₂= 0.7
T(⁰C)= 60
T(K)= 333.15
PS₁(kPa)= 20.274967368617 𝑦_𝑖 𝑃=𝑥_𝑖 𝛾_1 〖𝑃 _𝑖 〗 ^𝑆
(EC.10.5)
PS₂(kPa)= 20.007044710293
ϒ₁= 2.1244839754152
ϒ₂= 1.279246773451

INCISO (C)
DATOS: FORMULAS: RES
T(⁰C)= 60 x₁=
T(K)= 333.15 𝑍_𝑖=𝑥_𝑖 (1−𝑉)+𝑦_𝑖 𝑉 y₁=
 Z₁= 0.3 V=
x₁(supuesta)= 0.1 𝑃=(𝑃_𝐵𝑈𝐵𝐿+𝑃_𝐷𝐸𝑊)/
y₁(supuesta)= 0.5 2 Z₁=
V(supuesto)= 0.95
PBUBL(kPa)= 30.837916157361
PDEW(kPa)= 29.141774669746
 a ecuación (10.5) y en la ecuación de

INCISO(A)
LAS: RESULTADOS:
lnPs₁(kPa)= 3.009386991 Ps₁(kPa)= 20.274967369
_1/𝑉_2 𝐸𝑋𝑃 (−𝑎_21)/𝑅𝑇 lnPs₂(kPa)= 2.99608444705 Ps₂(kPa)= 20.00704471

𝑥_2+𝑥_1 𝐴_21 ))]/(𝑥_1+𝑥_2 𝐴_12 ) A₁₂= 0.07453721544

A₂₁= 0.53959575618

ln₁= -0.2900951463 exp₁= 0.7481923765

(𝑥_2+𝑥_1 𝐴_21 ))]/(𝑥_2+𝑥_1 𝐴_21 )
ln₂= 0.09763074016 exp₂= 1.1025555804

PBUBL(kPa)= 30.8379161574 y₁= 0.4190345715

y₂= 0.5809654285
=𝑥_𝑖 𝛾_1 〖𝑃 _𝑖 〗 ^𝑆 (EC.10.5)

RESULTADOS:
PDEW(kPa)= 29.1417746697

x₁= 0.20296615616
x₂= 0.79703384384

RESULTADOS
0.099720724 x₂= 0.90027927639
0.367343724 y₂= 0.63265627631
0.748362451

0.299999728 P(kPa)= 29.9898454136

OS:

ϒ₁= 2.124483975
ϒ₂= 1.279246773
 12.16 Para uno de lo sistemas binarios listados en la tabla 12.5, con base en la ecuacion
(10.5) y en la ecuacion de Wilson, realice los siguientes calculos:
a)BUBL P: t=60°C, Xi=0.3
b)DEW P: t=60°C, Yi=0.3
c)P,T-vaporizacion instantanea: t=60°C, P=(Pburbuja+Procio), Zi=0.3
d)Si existe un azeótropo a t=60°C, encuentre Paz y X1az=Y1az

12.17 Desarrolle el problema anterior para la ecuacion NRTL.

T= 60 °C R= 8.314
333.15 K

Calcular P1sat y P2sat ln P/kpa= A-B/°C+C

P1sat= 3.0093869905608 20.27496737 Kpa

P2sat= 2.9960844470483 20.00704471 Kpa

OBTENEMOS ϒ1 Y ϒ2 A PARITR DE ECUACIÓN DE NRTL

Teniendo X1=0.3 entonces X2= 0.7 Para obtener G21 y G12

Ecuacion de NRTL G12= exp (-aτ12)

G21= exp (-aτ21)

a 0.1 a 0.9

min. para G12 1.09 max para G21 0.74

a (-aτ12) G12=exp(-aτ) (-aτ21) G21=exp (-aτ21)

0.1 0.080517856 1.0838481995607 -0.054155356772 0.947284927941217
0.2 -0.01948214 0.9807064064334 -0.108310713543 0.897348734704596
0.3 -0.11948214 0.8873798526485 -0.162466070315 0.850044931492785
0.4 -0.21948214 0.8029344946876 -0.216621427087 0.80523475167594
0.5 -0.31948214 0.7265251750252 -0.270776783858 0.762786743717106
 0.6 -0.41948214 0.6573871635079 -0.32493214063 0.722576385556574
0.7 -0.51948214 0.5948285036785 -0.379087497402 0.684485719323984
0.8 -0.61948214 0.5382230874426 -0.433242854174 0.648403005306612
0.9 -0.71948214 0.4870043887689 -0.487398210945 0.614222394158743

T= 60 °C 333.15 K
X1= 0.3 (1-X1) 0.7

a) Presion de burbuja b) Presion de rocío

P= Σ XiPisat P= 1/ΣYi/Pisat

P= X1P1sat+ X2P2sat P= 1/ Y1/P1sat+ 1/Y2/ P2sat

P= 26.78953492 P= 96.16471701

d) Si existe azeótropo a T=60°C. Encuentre Paz y Xaz=Yaz

X1=Y1=0.3

α12= Y1/x1/Y1/X1=1 entonces hay un azeótropo

α12= 1

Volatilidad relativa

α12= ϒ1 P1sat
ϒ2 P2sat

ϒ1 = P1sat 20.27 1.0133914159839

ϒ2 P2sat 20

entonces ln (1.01339142)
0.0133025435125
Ln ϒ1/ϒ2= A( 1- 2X1 )
Ln ( 1.01339142)= A(1-2X1)
0.01330254= A( 1-2X1) 0.0052507030562
0.0133054/A = ( 1-2X1)
0.0052507-1/-2 = X1
 X1= 0.505250703

Calculamos coefciente de actividad del azeótropo

ϒaz= exp A( X1)^2 Paz= ϒazP1sat puesto que X1=Y1=1

ϒaz= 1.909310202
Paz= 38.71120205 Kpa

Calculo de A

ln ϒ1 - ln ϒ2= AX1^2- AX2^2 P1sat/P2sat 1.013391416

Ln ϒ1/ϒ2= A( X2^2- X1^2)

si ϒ1/ϒ2= P1sat/P2sat

Ln P1sat/P2sat= A( X2^2- X1^2)

despejamos

A= ln P1sat/P2sat/( X2^2- X1^2)

A= 2.53347854
 Valores para ecuación de Antoine
a b c
1-propanol 16.1154 3483.67 205.807
agua 16.3872 3885.7 230.17

NRTL

a obtener G21 y G12 Estimaciones de parámetros

b12 b21 τ12 τ21
τ12= b12/RT 150 1000 0.054155357 0.361035712
τ21= b21/RT 200 1500 0.072207142 0.541553568
300 2000 0.108310714 0.722071424
350 2500 0.126362499 0.90258928
400 3000 0.144414285 1.083107135
450 3500 0.16246607 1.263624991
500 4000 0.180517856 1.444142847
550 4500 0.198569641 1.624660703

Sustituyendo parámetros para encontrar coef. de actividad

ln ϒ1 0.206778218487 ϒ1= 1.229709815

ln ϒ2 0.111264896587 ϒ2= 1.11769094
 c) P y T de vaporización instantánea Z1=0.3

P= PB+ PR/2 Z1= 0.3

P= 61.477125963105 Kpa

K1= Ps/P 0.3297969293617

K2= Ps/P 0.3254388424452

z= Y1V + X1(1-V)

V Y1 X1 Z
0.1 0.8 0.85 0.845
0.2 0.7 0.6 0.62
0.3 0.6 0.2 0.32
0.4 0.5 0.4 0.44
0.5 0.4 0.3 0.35
0.6 0.3 0.5 0.38
0.7 0.2 0.7 0.35
0.8 0.15 0.9 0.3

P= X1ϒ1P1sar/ Y1 Tsat=( B/A-lnP)-C

Pvaporización= 149.59395818724 Kpa T1sat 107.8259273

T2sat 111.301657

T= X1P1sat + X2P2sat
T= 108.173500268539
 12.23 Desarrolle el problema anterior para la ecuación NRTL

PARÁMETROS
Parámetro b [cal/mol K] Parámetros a

Componentes 1 = Acetona 2 = Metano 3 = Agua 1 = Acetona 2 = Metanol

1 = Acetona 0 184.7 631.05 0.000 0.308

2 = Metanol 222.64 0 -253.88 0.308 0.000
3 = Agua 1197.41 845.21 0 0.534 0.299
R [cal/mol K]= 1.98721

DATOS P (kPa) = 101.33 x1 = 0.300000

x2 = 0.400000
x3 = 0.300000

FORMULAS

P (kPa) 101.33 Estimació CÁLCULOS PARA LA FASE LIQUIDA

n Inicial
T (K) 300.000 x1 0.300000
P1°(kPa) 33.18102 x2 0.400000
P2°(kPa) 18.51843 x3 0.300000
P3°(kPa) 3.55701 Ln (Gama 1) 0.422911
t 12 0.309815 Ln (Gama 2) 0.047698
t 13 1.058519 Ln (Gama 3) 0.502728
t 21 0.373455 Gama 1 1.52639813
t 23 -0.425857 Gama 2 1.04885346
t 31 2.008528 Gama 3 1.65322487
t 32 1.417750
G 12 0.908876
G 13 0.568218
G 21 0.891212
G 23 1.135987
G 31 0.341928
G 32 0.654114
or para la ecuación NRTL

ROS
ámetros a Constantes de Antoine

3 = Agua A B C

0.534 14.31450 2756.22 228.060

0.299 16.57850 3638.27 239.500
0.000 16.38720 3885.70 230.170

A LA FASE LIQUIDA CÁLCULO DE LA FUNCIÓN ERROR

K 1 = (P1° g1) / P 0.499827
K 2 = (P2° g2) / P 0.191682
K 3 = (P3° g3) / P 0.058034
y1=K1 x1 0.149948
y2=K2x2 0.076673
y3=K3 x3 0.017410
Sy = y1 + y2 + y3 0.244031
E = ln ( Sy) -1.410461
 P calculada 24.727644
P calc / P 0.24403083
E = ln (Pcalc/P) -1.410E+00

Composición del vapor

y1 0.614463
y2 0.314193
y3 0.071344
Total 1.000000
Temperatura de burbuja
Tburb(°K) 300.000000
 12.23 (b)
PARÁMETROS

Parámetro b[cal/mol K] Parámetros a

Componentes 1 = Acetona 2 = Metanol 3 = Agua 1 = Acetona

1 = Acetona 0 184.7 631.05 0.000

2 = Metanol 222.64 0 -253.88 0.308
3 = Agua 1197.41 845.21 0 0.534
R [cal/mol K]= 1.98721

DATOS P (kPa) = 101.33 y1 = 0.300000

y2 = 0.400000
y3 = 0.300000

Iteración 1
Estimación Inicial

P (kPa) 101.33 101.33

T (K) 300.000 301.000 RESULTADO
P1°(bar) 34.61300 34.61300
P2°(bar) 19.48918 19.48918
P3°(bar) 3.77166 3.77166
t 12 0.309815 0.308785
t 13 1.058519 1.055003
t 21 0.373455 0.372214
t 23 -0.425857 -0.424442
t 31 2.008528 2.001855
t 32 1.417750 1.413040
G 12 0.908876 0.909164
G 13 0.568038 0.569106
G 21 0.891212 0.891553
G 23 1.135987 1.135505
G 31 0.341928 0.343149
G 32 0.654114 0.655037
 12.23 (b)
PARÁMETROS

Parámetros a Constantes de Antoine

2 = Metanol 3 = Agua A B C

0.308 0.534 14.31450 2756.22 228.060

0.000 0.299 16.57850 3638.27 239.500
0.299 0.000 16.38720 3885.70 230.170

Iteración 1 Estimacion inicial

CÁLCULOS PARA LA FASE LIQUIDA
x1 0.079712 0.079712
x2 0.188759 0.188759
x3 0.731528 0.731528
Ln (Gama 1) 1.267406 1.265583
Ln (Gama 2) 0.276764 0.276359
Ln (Gama 3) 0.114844 0.114463
Gama 1 3.55162815 3.54515838
Gama 2 1.31885540 1.31832129
Gama 3 1.12169851 1.12127153
CÁLCULOS DE LAS RAZONES DE EQUILIBRIO
K 1 = (P1° g1) / P 1.213189 1.210979
K 2 = (P2° g2) / P 0.253660 0.253558
K 3 = (P3° g3) / P 0.041751 0.041735
x1= y1/K1 0.247282 0.247733
x2= y2/K2 1.576911 1.577550
x3= y3/K3 7.185398 7.188134
Sx = x1 + x2 + x3 9.009591 9.013418
E = ln ( Sx ) 2.198290 2.198714
P calculada 11.246903 11.242128
 P calc / P 0.110993 0.110946
E = ln (Pcalc/P) -2.198290 -2.198714
Nueva T 16.4835
x1 normalizada 0.027447
x2 normalizada 0.175026
x3 normalizada 0.797528
Total= 1.000000 RESULTADO
 12.23 C)
PARÁMETROS Parámetros b [cal/mol K] Parámetros a

Componentes 1 = Acetona 2 = Metanol 3 = Agua 1 = Acetona

1 = Acetona 0 184.7 631.05 0.000

2 = Metanol 222.64 0 -253.88 0.308
3 = Agua 1197.41 845.21 0 0.534
R [cal/mol K]= 1.98721

P (kPa) 101.330
T (K) 300.50
P1°(bar) 33.89085
P2°(bar) 18.99852
P3°(bar) 19.87352
t 12 0.309299099
t 13 1.056757967
t 21 0.372833522
t 23 -0.425148107
t 31 2.005185893
t 32 1.415390859
G 12 0.9090203
G 13 0.5685729
G 21 0.8913823
G 23 1.1357456
G 31 0.3425386
G 32 0.6545757

xi
0.025440494792
0.816974533875
0.157584971334
1
yi Nombre
0.07721619426 1 = Acetona
0.888018050427 2 = Metanol
0.034765755313 3 = Agua
 Nombre
1 = Acetona
2 = Metanol
3 = Agua

z i K i  1 
V
1 K i  1 
F

z i K i  1
S   V
1  K i  1
F
z i K i  1 
2

2
 V
 1  K i  1 
 F 

z i K i  1 
2
S   V
2
1  K i  1
 F 

V  V  S
     
 F  nueva  F  S 

z i K i  1
S   V
1  K i  1
F

z i K i  1 
2
S   V
2
1  K i  1
 F 

V  V  S
     
 F  nueva  F  S 

z i K i  1
S   V
1  K i  1 
F
 z i K i  1
S   V
1  K i  1 
F

z i K i  1
2
S   V
2
1  K i  1
 F 

V  V  S
     
 F  nueva  F  S 

z i K i  1
S   V
1  K i  1 
F

Zi
xi 
V
1 K i  1 
F
S x   xi

Zi Ki
yi 
V
1  K i  1 
F
S y   yi

xi
xi 
Sx

yi
yi 
Sy
 12.23 C)
Parámetros a Constantes de Antoine

2 = Metanol 3 = Agua A B C zi

0.308 0.534 14.31450 2756.22 228.060 0.300000

0.000 0.299 16.57850 3638.27 239.500 0.200000
0.299 0.000 16.38720 3885.70 230.170 0.500000

Estimación
inicial
FASE LÍQUIDA
x1 0.900000 K 1 = (P1° g1) / P 0.336998
x2 0.090000 K 2 = (P2° g2) / P 0.312195
x3 0.010000 K 3 = (P3° g3) / P 0.856152
Estimació
Suma 1.000000 n inicial
Ln (Gama 1) 0.007561 V/F 0.843400
Ln (Gama 2) 0.509895
Ln (Gama 3) 1.473687 Error Eq 1 -2.9330E-01
Gama 1 1.00758917 Error Eq 2 6.1902E-02
Gama 2 1.66511696 Error Eq 3 8.9144E-01
Gama 3 4.36529894 Error BM 1 -1.5063E-01
Estimació
n inicial Error BM 2 -1.1000E-01
FASE VAPOR Error BM 3 2.6063E-01
y1 0.010000
y2 0.090000
y3 0.900000
Suma 1.000000

Iteración 1

Iteración 1
zi
0.300000
0.200000
0.500000
 Ki
0.33700
0.31220
0.856152

V/F
0.5000

-0.2975328574
-0.20966532342
-0.07749803283
-0.58469621364

-0.2950860041
-0.21979773922
-0.01201189018
-0.52689563351

V/F
-0.60970024509

-0.14164357807
-0.0969179403
-0.0661246481
-0.30468616648

-0.06687634403
-0.04696543576
-0.00874493817
-0.12258671797

V/F
-3.09517478529

-0.06516827729
-0.04396497254
-0.0497663398
-0.15889958963

-0.01415634788
-0.00966459405
-0.00495337715
-0.02877431909

V/F
0.471258360514 RESULTADO

-0.28928666191
-0.20353275419
-0.07715432258
-0.56997373868

Cálculo de nuevas composiciones:

xi
0.4363287580
0.2959165120
0.5363596196
1.2686048896
yi
0.1470420961
0.0923837579
0.4592052970
0.6986311509

Cálculo de nuevas composiciones

xi
0.3439437776
0.2332613680
0.4227948544
1.0000000000
yi
0.2104717144
0.1322353831
0.6572929025
1.0000000000
 12.25 Las posibles ecuaciones de correlacion para lnƔ₁ en un sistema binario liquido se dan a continuación.
Para uno de estos casos determine mediante integración de la ecuación de Gibbs/Duhem [ecuación (11.100)], la
ecuación correspondiente para lnƔ₂. ¿Cuál es la ecuación correspondiente para GE/RT? Observe que por su
definición, Ɣ₁=1 para Xi=1.

a)lnƔ₁=Ax²₂
b)lnƔ₁=x²₂(A+Bx₂)
c)lnƔ₁=x²₂(A+Bx₂+Cx²₂)

De la ecuación 11.100 para un sistema binario y diviendo en por Xi

𝑥_1 (ⅆ ln⁡〖𝛾 _1 〗)/(ⅆ𝑥_1 )+(𝑥_2 ⅆ ln⁡〖𝛾 _2

de 〗)/(ⅆ𝑥_1 )=0
donde (ⅆ ln⁡〖𝛾 _2 〗)/(ⅆ𝑥_1 )=−𝑥_1/𝑥_2
(ⅆ ln⁡〖𝛾 _1 〗)/(ⅆ𝑥_1 )=𝑥_1/𝑥_2 (ⅆ
ln⁡〖𝛾 _1 〗)/(ⅆ𝑥_2 )
Integrar recordando que para lnƔ₂= 1 para Xi=0

ln⁡〖𝛾 _2 〗=ln⁡(1)+∫130_0^(𝑥_1)▒ 〖𝑥 _1/𝑥_2 (ⅆ ln⁡〖𝛾 _1 〗)/(ⅆ𝑥_2 ) ⅆ𝑥_1 〗=∫130_0^(𝑥_1)▒ 〖𝑥 _1/𝑥_2

(ⅆ ln⁡〖𝛾 _1
〗)/(ⅆ𝑥_2 )=2𝐴𝑥_2

Inciso a)
ln⁡〖𝛾 _2 〗=2𝐴∫130_0^(
de donde 𝑥_1)▒ 〖𝑥 _1 ⅆ𝑥_1 〗 o

𝐺^𝐸/𝑅𝑇=𝐴𝑥_1 𝑥_2

Inciso b)

(ⅆ ln⁡〖𝛾 _1 〗)/(ⅆ𝑥_2 )=2𝑥_2 (𝐴+𝐵𝑥_2 )+𝑥_2^2

𝐵=2𝐴𝑥_2+3𝐵𝑥_2^2=2𝐴𝑥_2+3𝐵𝑥_2 (1−𝑥_1 )

ln⁡〖𝛾 _2 〗=2𝐴∫130_0^(𝑥_1)▒ 〖𝑥 _1 ⅆ𝑥_1 〗+3𝐵∫130_0^(𝑥_1)▒ 〖𝑥 _1 ⅆ𝑥_1 〗

− 3𝐵∫130_0^(𝑥_1)▒ 〖𝑥 _1^2 ⅆ
de donde

ln⁡〖𝛾 _2 〗=𝐴𝑥_1^2+3𝐵/2 𝑥_1^2−𝐵_1^3 ln⁡〖𝛾 _2 〗=𝑋_1^2 (𝐴+3𝐵/2−𝐵𝑥_1 )=𝑥

o [𝐴+𝐵/2 (1+2𝑥_2 )]
 ln⁡〖𝛾 _2 〗=𝑋_1^2 (𝐴+3𝐵/2−𝐵𝑥_1 )=𝑥
[𝐴+𝐵/2 (1+2𝑥_2 )]

𝐺^𝐸=𝑥_1 𝑥_2^2 (𝐴+𝐵𝑥_2 )+𝑥_2 𝑥_1^2 (𝐴+3𝐵/2−𝐵𝑥_1 )

Reducción algebraica

𝐺^𝐸/𝑅𝑇=𝑥_1 𝑥_2 [𝐴+𝐵/2

(1+𝑥_2 )]

Inciso c)

(ⅆ ln⁡〖𝛾 _1 〗)/(ⅆ𝑥_2 )=2𝑥_2 (𝐴+𝐵𝑥_2+𝐶𝑥_2^2 )+𝑥_2^2 (𝐵+2𝐶𝑥_2 )=2𝐴𝑥_2+3𝐵𝑥_2^2+4𝐶𝑥_2^3

2𝐴𝑥_2+3𝐵𝑥_2 (1−𝑥_1 )+4𝐶𝑥_2 (1−𝑥_1 )^2

de donde

ln⁡〖𝛾 _2 〗=2𝐴∫130_0^(𝑥_1)▒ 〖𝑥 _1 ⅆ𝑥_1 〗+3𝐵∫130_0^(𝑥_1)▒ 〖𝑥 _1

o
(1−𝑥_1 )ⅆ𝑥_1 〗+4𝐶∫130_0^(𝑥_1)▒ 〖𝑥 _1 (1−𝑋_1 )^2 ⅆ𝑥_1 〗

ln⁡〖𝛾 _2 〗=((2𝐴+3𝐵+4𝐶)/2)
𝑥_1^2−((3𝐵+8𝐶)/3) 𝑥_1^3+𝐶𝑥_1^4

ln⁡〖𝛾 _2 〗=𝑥_1^2 [𝐴+3𝐵/2+2𝐶−(𝐵+8𝐶/3)𝑥,

+𝐶𝑥_1^7 ]

ln⁡〖𝛾 _2 〗=𝑥_1^2 [𝐴+𝐵/2 (1+2𝑥_2 )+𝐶/3

(1+2𝑥_2+3𝑥_2^2 )]

Reduciendo la ecuacion
𝐺^𝐸/𝑅𝑇=𝑥_1 𝑥_2 [𝐴+𝐵/2 (1+𝑥_2 )+𝐶/3
(1+𝑥_2+𝑥_2^2 )]
 dan a continuación.
m [ecuación (11.100)], la
Observe que por su

r Xi

(ⅆ𝑥_1 )=−𝑥_1/𝑥_2
(ⅆ𝑥_1 )=𝑥_1/𝑥_2 (ⅆ
_2 )

0_0^(𝑥_1)▒ 〖𝑥 _1/𝑥_2 (ⅆ ln⁡〖𝛾 _1 〗)/(ⅆ𝑥_2 ) ⅆ𝑥_1 〗

ln⁡〖𝛾 _2 〗=𝐴𝑥_1^2

30_0^(𝑥_1)▒ 〖𝑥 _1^2 ⅆ𝑥_1 〗

^2 (𝐴+3𝐵/2−𝐵𝑥_1 )=𝑥_1^2
2 )]
^2 (𝐴+3𝐵/2−𝐵𝑥_1 )=𝑥_1^2
2 )]

ln⁡〖𝛾 _2 〗=(2𝐴+3𝐵+4𝐶) ∫130_0^(𝑥_1)▒ 〖𝑥 _1 ⅆ𝑥_1 〗

− (3𝐵+8𝐶)
∫130_0^(𝑥_1)▒ 〖𝑥 _1^2 ⅆ𝑥_1 〗+4𝐶∫130_0^(𝑥_1)▒ 〖𝑥 _1^3 ⅆ𝑥_1 〗
 12.26 A 25°C y a la presion atmosferica, el cambio de volumen de mezclado de Datos
combinaciones liquidas binarias de las especies 1 y 2 se da por la ecuación: ΔV=
x1x2(45x1+25x2) en donde ΔV está en cm mol. En estas condiciones, V1= 110 y V2=
90 cm mol. Calcule los volumenes molares parciales V1 y V2, en una mezcla que
contiene 40% mol de la especie 1 a las condiciones indicadas.
 T V1(cm^3/mol) V2(cm^3/mol) X1
25 110 90 0.4

Formulas V= ΔV+ΣXiVi
ΔV= X1X2(45X1+25X2) V=
X1X2(45X1+25X2)+X1V1+X2V2
X2= 1-X1 dV/dX1=(d/dX1)
_ ΔV= V-ΣXiVi (X1X2(45X1+25X2)+X1V1+X2V2)
dV/dX1=
_ V1 = V+X2(dV/dX1) 45X1X2+X2(45X1+25X2)+V1
V2 = V-X1(dV/dX1)
Calculos
X2 ΔV(cm^3/mol) V(cm^3/mol)

_ 0.6 _7.92 105.92

V1(cm^3/mol) (dV/dX1) V2(cm^3/mol)
190.28 140.6 49.68
12.27 El cambio de volumen de mezclado (cm3 mol-1) para el sistema etanol(1)/éter metil butílico(2) a 25°C se d

∆𝑉=X1X2[-1.026+0.220(X1-X2)]

Si V1=58.63 y V2=118.46 cm3 mol-1. ¿qué volumen de la meezcla se forma cuando 750 cm3 de la espe
con 1500 cm3 de la especie 2 a 25°C?¿cuál sería el volumen si se formara una solución ideal?

Datos
V1 58.63
V2 118.46
V especie 1 750
V especie 2 1500

moles1 12.79208596 x1 0.502545393

moles 2 12.66250211 x2 0.497454607
moles 25.45458807
 Formulas
1)/éter metil butílico(2) a 25°C se da por la ecuación:
𝑚𝑜𝑙𝑒𝑠1=(𝑉 𝑒𝑠𝑝𝑒𝑐𝑖𝑒1)/𝑣1
X1-X2)]

orma cuando 750 cm3 de la especie pura 1 se mezclan 𝑚𝑜𝑙𝑒𝑠2=(𝑉

ra una solución ideal? 𝑒𝑠𝑝𝑒𝑐𝑖𝑒2)/𝑣2

𝑚𝑜𝑙𝑒𝑠=𝑚𝑜𝑙𝑒𝑠1+𝑚𝑜𝑙𝑒𝑠2

𝑥1=(𝑚𝑜𝑙𝑒𝑠1 )/𝑚𝑜𝑙𝑒𝑠 𝑥2=1−𝑥2

(𝐸𝑞.12.27) 𝑉=∆𝑉+𝑥1∗𝑉1+𝑥2∗𝑉2

𝑉𝑡𝑜𝑡𝑎𝑙=𝑉∗𝑚𝑜𝑙𝑒𝑠

∆𝑉 -0.25621337 (𝐸𝑞.11.81) 𝑉𝑖𝑑𝑒𝑎𝑙=(𝑥1∗𝑉1+𝑥2∗𝑉2)∗𝑚𝑜𝑙𝑒𝑠

V 88.13649578
Vtotal 2243.478194
Videal 2250
12.30 ¿Cual es el efecto termico cuando 20 kg de LiCl(s) se agregan a 125 kg de una
solucion acuosa que contiene 10% en peso de LiCl en un proceso isotermico a 25°C?

Moles de LiCl 0.2948808681 Moles de LiCl agregado 0.471809389

Moles de H2O 6.2447960033 Relacion de moles de la solucion original

Relacion de moles de la solucion final

Δ𝐻_1 741.20732723 Q 756.60010854366

Δ𝐻_2 15.392781316
riginal 21.17735221

inal 8.145135464
12.33 Una solucion al 20% mol de LiCl/H2O, se prepara por seis diferentes procesos de mezclado:
a) Mezcle LiCl(s) con H2O(l)
b) Mezcle H2O(l) con una solucion al 25% mol de LiCl/H2O
c) Mezcle LiCl * H2O con H2O(l)
d) Mezcle LiCl(s) con una solucion al 10% mol de LiCl/H2O
e) Mezcle una solucion al 25% mol de LiCl/H2O con una solucion al 10% mol de LiCl/H2O
f) Mezcle LiCl*H2O(s) con una solucion al 10% mol de LiCl/H2O.

En todos los casos, el mezclado es isotermico a 25°C. Para cada efecto determine el efecto termico en J/MOL de la solucio

T=Cte
P=Cte

a) LiCl + H2O ----------------> LiCl(4H2O)

DATOS: Con la figura 12.14 calcular

X1= 0.2
𝐻=𝑋𝑖𝐻 ̃
𝐻 ̃=𝐻/
𝑋𝑖
𝐻 ̃=-25.5 KJ/mol 𝐻=(0.2 MOL)(-25.5KJ/MOL)
𝐻=-5.1KJ
=-5100J

b) LiCl(3H2O)-------------------> LiCl + 3H2O (1)

LiCl + 4H2O-----------------> LiCl(4H2O) (2)
LiCl(3H2O) + H2O---------> LiCl(4H2O)
DATOS: 𝐻=(0.2MOL)(20.756-25.5 KJ/MOL)
X1=0.2MOL 𝐻=-0.948KJ
𝐻=-948J

𝐻 ̃1=20.756KJ/MOL
𝐻 ̃2=-25.5KJ/MOL

c) LiCl*H2O--------------------> Li + 1/2 Cl2 + H2 +1/2 O2 (1)

H2 + 1/2 O2---------------> H2O (2)
Li +1/2 Cl2----------------> LiCl (3)
LiCl + 4H2O---------------> LiCl(4H2O) (4)
LiCl*H2O + 3H2O---------------> LiCl(4H2O)

DATOS: 𝐻=(0.2MOL)(712.58-285.83-408.61-25.5KJ/MOL)
X1=0.2MOL 𝐻=-1.742KJ
𝐻=-1472J

𝐻 ̃1=712.58KJ/MOL
(𝐻2) ̃=-285.83KJ/MOL
(𝐻3) ̃=-408.61KJ/MOL
(𝐻4) ̃=-25.5KJ/MOL
de mezclado:

Cl/H2O
FIGURA 12.14
fecto termico en J/MOL de la solucion final.

LiCl(4H2O)

ura 12.14 calcular 𝐻

 ̃

𝐻=𝑋𝑖𝐻 ̃

L)(-25.5KJ/MOL)
=-5.1KJ
=-5100J

2O (1)
(2)
2O)
0.756-25.5 KJ/MOL)

𝐻 ̃1=-25.5K
𝐻 ̃2=4/9(32

+ H2 +1/2 O2 (1)
(2)
(3)
(4)
LiCl(4H2O)

408.61-25.5KJ/MOL) 𝐻 ̃1=5/6(20
𝐻 ̃2=1/6(32
𝐻 ̃3=-25.5K
𝐻 ̃1=5/8(71
𝐻 ̃2=5/8(-2
𝐻 ̃3=3/8(32
𝐻 ̃4=5/8(-4
𝐻 ̃5=-25.5K
A 12.14

d) LiCl + 4H2O------------------> LiCl(4H2O) (1)

4/9[LiCl(9H2O)---------------> LiCl + 9H2O] (2)
5/9LiCl + 4/9LiCl(9H2O)---------------> LiCl(4H2O)
 DATOS: 𝐻 =(0.20MOL)(-25.5 + 14.4KJ/MOL)
X1=0.20MOL 𝐻 =-2.22KJ
𝐻=-2220J

𝐻 ̃1=-25.5KJ/MOL
𝐻 ̃2=4/9(32.4KJ/MOL)

e) 5/6[LiCl(3H2O)---------------------> LiCl + 3H2O] (1)

1/6[LiCl(9H2O)---------------------> LiCl + 9H2O] (2)
LiCl + 4H2O------------------------> LiCl(4H2O)] (3)
5/6 LiCl +4H2O + 1/6 LiCl(9H2O)-----------> LiCl(4H2O)

DATOS: 𝐻=(0.2MOL)(17.296 + 5.4 - 25.5KJ/MOL)

X1=0.20MOL 𝐻=-0.5608KJ
𝐻=-560.8J

𝐻 ̃1=5/6(20.756KJ/MOL)
𝐻 ̃2=1/6(32.4KJ/MOL)
𝐻 ̃3=-25.5KJ/MOL

f) 5/8[LiCl*H2O-------------------> Li + 1/2 Cl2 + H2 +1/2 O2] (1)

5/8[H2 + 1/2 O2--------------> H2O ] (2)
3/8[LiCl(9H2O)---------------->LiCl + 9H2O] (3)
5/8[Li + 1/2 Cl2---------------> LiCl] (4)
LiCl + 4H2O--------------------> LiCl(4H2O) (5)
5/8 LiCl*H2O +3/8 LiCl(9H2O)-------------------> LiCl(4H2O)
DATOS:
X1=0.20MOL
𝐻=(0.2MOL)(445.36 - 178.64 + 12.15 - 255.38 - 25.5KJ/MOL)
𝐻=-0.402KJ
𝐻=-402J
 𝐻=(0.2MOL)(445.36 - 178.64 + 12.15 - 255.38 - 25.5KJ/MOL)
𝐻=-0.402KJ
𝐻=-402J

𝐻 ̃1=5/8(712.58KJ/MOL)
𝐻 ̃2=5/8(-285.83KJ/MOL)
𝐻 ̃3=3/8(32.4KJ/MOL)
𝐻 ̃4=5/8(-408.61KJ/MOL)
𝐻 ̃5=-25.5KJ/MOL
 (1)
2)
(3)
(4)
(5)

J/MOL)
J/MOL)
12.34 Una masa de 12 kg/s de Cu(NO3)2 ∙6H2O junto con 15 kg/s de agua, ambos a 25ºC, alimentan a un tanque donde ocur
resulta pasa a traves de un intercambiador de calor que ajusta su temperatura a 25ºC. ¿Cual es la rapidez de transferencia d
-Para Cu(NO3)2, 〖∆𝐻 ° 〗 _𝑓298=−302.9 𝑘𝐽
-Para Cu(NO3)2∙6𝐻2𝑂, 〖∆𝐻 ° 〗 _𝑓298=−2110.8 𝑘𝐽
-El calor de la solucion de 1 mol de Cu(NO3)2 en agua a 25ºC es -47,84 kJ, en este caso independiente de ñ para l

DATOS DESARROLLO
m de Cu(NO3)2= 12 kg n1= 0.041 kmol/seg
m de H2O= 15 kg n2= 0.833 kmol/seg
PM de Cu(NO3)2= 295.61 kg/mol (6∙𝑛1+𝑛2)/
𝑛1=
PM de H2O= 18.015 kg/mol 26.51
Base de t= 1 seg
T= 298.15 ºK
 alimentan a un tanque donde ocurre el mezclado. La solucion que
ual es la rapidez de transferencia de calor en el intercambiador?
−302.9 𝑘𝐽
8=−2110.8 𝑘𝐽
este caso independiente de ñ para los valores de interes.

REACCIÓN MOLAR
6(𝐻_2+1⁄2 𝑂_2→𝐻_2 𝑂_((𝑙) ) ) (1) VALORES DE LA TA
𝐶𝑢+𝑁_2+3𝑂_2→𝐶𝑢(𝑁𝑂_3 )_(2 ) (2)
𝐶𝑢(𝑁𝑂_3 )_2∙6𝐻_2 𝑂→𝐶𝑢+𝑁_2+6𝑂_2+6𝐻_2 (3)
𝐶𝑢(𝑁𝑂_3 )_2+20,51𝐻_2 𝑂→𝐶𝑢(𝑁𝑂_3 )_2 (20,51 𝐻_2 𝑂) (4) ΔH1=
ΔH2=
Cu(N𝑂_3 )_2 ∙6𝐻_2 𝑂+14,51𝐻_2 𝑂_((𝑙) )→𝐶𝑢(𝑁𝑂_3 )_2 (20,51𝐻_2 𝑂)
ΔH3=
ΔH4=

∆𝐻=∆𝐻1+∆𝐻2+∆𝐻3+∆𝐻4

ΔH= 45.08 kJ

CALCULO PARA 1 MOL DE Cu(NO3)2*6H2O SOBRE LA BASE DE 1 SEGUNDO

𝑄=𝑛1∙∆𝐻/𝑚𝑜𝑙

Q= 1.830 kJ/seg
VALORES DE LA TABLA C.4

-1714.98 kJ
-302.9 kJ
2110.8 kJ
-47.84 kJ
 Ejercicio 12.37

ΔHƭ @25°C
CaCl2 en 10 mol de H2O -862.74 Kj CaCl2
CaCl2 en 15 mol de H2O -867.85 Kj
CaCl2 en 20 mol de H2O -870.06 Kj Procedimiento
CaCl2 en 25 mol de H2O -871.07 Kj Ca + Cl2+ H2O→┴ 𝐶𝑎𝐶𝑙2(𝐻20)
CaCl2 en 50 mol de H2O -872.91 Kj 𝐶𝑎𝐶𝑙2→┴ Ca + Cl2
CaCl2 en 100 mol de H2O -873.82 Kj
CaCl2 en 300 mol de H2O -874.79 Kj 𝐻2𝑂+𝐶𝑎𝐶𝑙2→┴ 𝐶𝑎𝐶𝑙2(𝐻2𝑂)
CaCl2 en 500 mol de H2O -875.13 Kj
CaCl2 en 1000 mol de H2O -875.54 Kj
 ΔHƭ @25°C -60
0 200 400 600 800 1000 1200
-795.8 kJ

-65
l2+ H2O→┴ 𝐶𝑎𝐶𝑙2(𝐻20)
𝐶𝑎𝐶𝑙2→┴ Ca + Cl2
-70

+𝐶𝑎𝐶𝑙2→┴ 𝐶𝑎𝐶𝑙2(𝐻2𝑂) Axis Title

Column P

-75

-80

-85
 ΔH ň
-66.94 10
-72.05 15
-74.26 20
-75.27 25
-77.11 50
-78.02 100
-78.99 300
-79.33 500
-79.74 1000
 Ejercicio 12.38

Procedimiento
Necesitamos datos del ejercicio anterior y
hacer balance de ecuaciones para obetener
los ΔH.

𝐶𝑎𝐶𝑙2 →𝐶𝑎+𝐶𝑙2
2[Ca + Cl2 + 12.5H20→CaCl2(12.5 H2O)]
𝐶𝑎𝐶𝑙2(25𝐻2𝑂)→𝐶𝑎+𝐶𝑙2+25𝐻2𝑂

𝐶𝑎𝐶𝑙2+𝐶𝑎𝐶𝑙2(25𝐻2𝑂) →2[𝐶𝑎𝐶𝑙2(12.5𝐻20)]
 De ahí se obtienen los ΔH de cada compuesto final.

ΔH(CaCl2) 795.8
ΔH(CaCl2(25H2O)) 871.07
ΔH(CaCl2(12.5H2O)) -1730.59

ΔH(total) -63.72
Debido al proceso
Q=ΔH(total) -63.72
12.45. Una masa de 400 (lbm) de solución acuosa al 35% en peso de NaOH a 130 (°F), se mezcla con 175 (lbm) de solu
a) ¿Cuál es el efecto térmico si la temperatura final es 80 (°F)
b) Si el mezclado es adiabático ¿Cuál es la temperatura final?

DATOS 1) Calculamos con la tabla 12.19 las

1)
m1 (lbm) 400
%X1 35
X1 0.35
T (°F) 130
2)
m2 (lbm) 175
%X2 10
X2 0.1
T (°F) 200
3)
T (°F) 80
ENTALPÍAS CALCULADAS
H1 (Btu/lbm) = 100
H2 (Btu/lbm) = 152

a)
Fracción masa final (X3)
𝑋_3=( 〖 %𝑋 〗 _1 𝑚_1+ 〖 %𝑋 〗 _2 𝑚_2)/

X3 (%) = 27.3913043478261

ENTALPÍA FINAL
H3 Btu/lbm = 41
 MASA FINAL
𝑚_3=𝑚_1+𝑚_2

m3 (lbm) = 575

EFECTO TÉRMICO

Q =𝑚_3 𝐻_3− 〖 (𝑚 〗 _1 𝐻_1+ 〖 𝑚〗

Q (Btu) = -43025
se mezcla con 175 (lbm) de solución al 10% en peso a 200 (°F).
ura final es 80 (°F)
temperatura final?

amos con la tabla 12.19 las entalpías con las temperaturas dadas

b)
En un proceso adiabático Q = 0
〖 %𝑋 〗 _2 𝑚_2)/(𝑚_1+𝑚_2 ) ENTALPÍA FINAL
𝐻_3=(𝑚_1 𝐻_1+𝑚_2 𝐻_2)/𝑚_3

H3 (lbm) = 115.826087

Calculamos T final con la tabla 12.19 y H3

T final (°F) = 165
〗 _1 𝐻_1+ 〖 𝑚〗 _1 𝐻_1)
12.46 Un evaporador de simple efecto concentra una solución acuosa de H2SO4 de 20 a 70% en peso. La
relacion y la temperatura de la alimentacion son 25(lbm)(s)-1 y 80(°F), respectivamente. El evaporador se
mantiene a una presion absoluta de 1.5(psia), a la cual el punto de ebullicion del H2SO4 al 70% es 217(°F).
¿Cuál es la rapidez de transferencia de calor en el evaporador?

m1 25
Temperatura 80
Pabs 1.5 m2=
X1 0.2 P.ebullicion 217 m3=
X2 0.7
H1 -24
H2 -55
H3 1157.7
peso. La
orador se
s 217(°F). 𝑀2=(x1 x m1)/𝑋2
𝑀3=𝑚1 −𝑚2 𝑄=𝑚2 𝑥 𝐻2+ 𝑚3 𝑥 𝐻3 −

Rapidez de transferencia
7.142857143 Q= 20880.35714
17.85714286
𝑄=𝑚2 𝑥 𝐻2+ 𝑚3 𝑥 𝐻3 − 𝑚1 𝑥 𝐻1
12.52 Una carga de solucion al 40% en peso, de NaOH, en agua a la presion atmosferica y 80(°F)
se calienta en un tanque aislado por inyeccion a presion de vapor a traves de una valvula
conectada a una linea que contiene vapor saturado a 35 (psia). Se detiene el proceso cuando la
solucion de NaOH alcanza una concentracion de 38% en peso. ¿A que temperatura ocurre esto?

UTILIZANDO LA GRAFICA 12.17 A 80° F Formulas

m1 1 𝑚2=(𝑥1∗𝑚1)/𝑥3−𝑚1
x1 0.4
H1 77
𝑚3=𝑚1+𝑚2
UTILIZANDO TABLA F4 A 35 psia
x3 0.38 𝐻3=(𝑚1∗𝐻1+𝑚2∗𝐻2)/𝑚
H2 1161.1

m2 0.052631579
m3 1.052631579
H3 131.205

UTILIZANDO LA GRAFICA 12.17 A H3 OBTENIDA A 38% LA TEMPERATURA ES

T (°F) 155
 Formulas
=(𝑥1∗𝑚1)/𝑥3−𝑚1

3=𝑚1+𝑚2

3=(𝑚1∗𝐻1+𝑚2∗𝐻2)/𝑚3
12.55 Una solucion liquida que contiene 2 lb/mol de H2SO4 y 15 lb/mol de H2O a 100°F, formando una
solucion mas concentrada de H2SO4. Si el proceso ocurre isotermicamente, determine el calor transferido.

Solucion inicial 0.420591472 PM1

Solucion final 0.538451966 PM2

Datos de la figura 12.17 a 100°F Demezclar la solucion inicial

H agua 68 𝜟𝑯 118.1851032
H1 -75
H acido 9
H2 -101

Reaccionar 1 mol de SO3 (g) con 1 mol de H2O (l) a partir de 1 mol de H2SO4(l).
Ignoramos el efecto de la reaccion, tomando el valor a 100°F igual al valor a 77°F (25°C)

ΔH SO3 -395720 Calculamos el ΔH reaccion

ΔH H2O -285830 ΔH reaccion -132439
ΔH H2SO4 -813989

Finalmente mezclar los sustituyente a partir de la solucion final

ΔH mezclado -137.231334

A partir de los datos obtenidos, calculamos el calor transferido Q

Q -236565.769
100°F, formando una
ne el calor transferido.

98.0768
18.0148

5°C)
 12.58 Los liquidos siguientes, todos a presion atmosferica y 120(°F), se mezclan: 25(lbm) de agua pura, 40(lbm) de ácido
sulfúrico puro y 75(lbm) de ácido sulfúrico al 25% en peso.

a) ¿Cuánto calor se libera si el mezclado es isotérmico a 120(°F)?

b) El proceso de mezclado se realiza en dos etapas: primera, se mezclan ácido sulfúrico puro y la solución al 25% en peso y
se extrae el calor total del inciso a); segunda, el agua pura se agrega adiabaticamente. ¿Cuál es la temperatura de la solución
intermedia formada en la primera etapa?

1.- 25 lbm H20 pura INCISO A)

M4= ?
2.- H2SO4 40 lbm puro X4=? tenemos que:

M1= 40 M2
X1 1 X2
3.- 75 lbm H2SO4 .25

con la información de la tabla 12.17 obtenemos H

H1= 88 BTU/lbm
H2= 14 BT/lbm
H3= -7 BTU/lbm

necesitamos obtener la composición salida del sistema

X4M4 = X1M1+X2M2+ X3M3

X4= X1M1+X2M2+ X3M3 / M4

X4= 0.419642857

calculamos el calor liberado

Q= M4H4 - ΣMH

Q= M4H4 - M1H1+M2H2+MH3
Q= -12515 BTU
a, 40(lbm) de ácido

ción al 25% en peso y

eratura de la solución

INCISO B)

M1 40 M2 75
X1= 1 X2= 0.25

M3= M1 + M2
25 M3 75 M3= 115 lbm
0 X3 0.25

X3M3= X1M1+X2M2
12.17 obtenemos H X3= X1M1+X2M2 /M3

necesitamos la masa salida del sistema X3= 0.510869565

M4= Σ M
M4= M1 + M2 + M3
M4 = 140 lbm calculamos H 3

ición salida del sistema Q= M3H3- ΣMH

Q + ΣMH /M3= H3

H3= -108.521739
H4 según la figura 12.17 = -63 BTU/lbm
si nos fijamos en la tabla 12.17 por H y su composición
la temperatura de la solución se encuentra a 80°F
12.64 20000 lbm/hrs de una solucion de 80% en peso de H2SO4 en agua a 120°F se diluye de
manera continua con agua fria a 40°F para producir una corriente que contiene 50% en peso de
H2SO4 a 140°F.

a) ¿Cuál es la relación de flujo de masa del agua de enfriamiento en lbm/hrs?

b) ¿Cuál es la rapidez de transferencia de calor en Btu/hrs para el proceso de mezclado? ¿El

calor se agrega o se elimina?

c) Si el mezclado ocurre de manera adiabática, ¿Cuál sería la temperatura de la corriente de

productos? Suponga las mismas condiciones de entrada y la misma composición del producto
que en el inciso b).

Formulas
H2SO4 .8(A)+0(B)=.5(C)
H2O .2(A)+1(B)=.5(C)
Calculos
Corriente B Corriente C
(lbm/hrs) (lbm/hrs)

12000 32000

INCISO (B)
DATOS: FORMULAS:
Corriente A= 20000 𝑄_(𝑡𝑟𝑎𝑛𝑠.)= 〖𝐶𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒〗 _𝐶 𝐻_3−( 〖𝐶𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒〗 _𝐴 𝐻_1+ 〖𝐶𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒〗 _𝐵 𝐻_2)
Corriente B= 12000
Corriente C= 32000
H₁(BTU/lb)= -92
H₂(BTU/lb)= 7
H₃(BTU/lb)= -70

INCISO (C)
DATOS: FORMULA:
Corriente A= 20000
Corriente B= 12000 𝐻_(𝑎𝑑𝑖𝑏.)=( 〖𝐶𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒〗 _𝐴 𝐻_1+ 〖𝐶𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒〗 _𝐵 𝐻_2)/ 〖𝐶𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒〗 _𝐶
Corriente C= 32000
H₁(BTU/lb)= -92
H₂(BTU/lb)= 7
 Corriente A
X1(H2SO4) X2(H20) T1(°F)
(lbm/hrs)
Datos 20000 0.8 0.2 120
X1(H2SO4) X2(H20) T2(°F)
0 1 40
X1(H2SO4) X2(H20) T3(°F)
0.5 0.5 140
B=?
1 H2O
0 H2SO4
Balance general de materia
A=20000 lbm/hrs C=? A+B=C

Mezclador

.8 H2SO4 .5 H2SO4
.2 H2O .5 H2O

NCISO (B)
AS: RESULTADO:
_𝐴 𝐻_1+ 〖𝐶𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒〗 _𝐵 𝐻_2) Qtrans.(BTU/h)= -484000

NCISO (C)
LA: RESULTADO:
Hadib.(BTU/h)= -54.875
𝑛𝑡𝑒〗 _𝐵 𝐻_2)/ 〖𝐶𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒〗 _𝐶
T(⁰F)Aprox.= 170
H₁(BTU/lb)

-92
H₂(BTU/lb)
7
H₃(BTU/lb)
-70