INSTITUTO POLITÉCNICO NACIONAL

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

 

Departamento de Ingeniería Química Petrolera

 
Termodinamica de equilibrio de fases
 
Nestor Lorenzo Diaz Ramirez
 
Examen 2do. departamental
 
Sabido Torres Mauricio Miguel
 
Grupo: 2PM32
 
Fecha entrega
10 – Junio - 2020
IONAL
DUSTRIAS EXTRACTIVAS

rolera

s
 Using the data in Illustration 6.4-1, and the same reference state, compute the vapor pressure of oxygen ove
to the critical temperature, and also compute the specific volume, enthalpy, and entropy along the vapor-liqu
Add these results to Figs. 6.4-3, 6.4-4, and 6.4-5.

T°C TK Pvap Pliq

-200 73.15 0.11 0.07658932 presion de
-190 83.15 0.47 0.36122868
-180 93.15 1.39 1.42355647
-175 98.15 2.19 2.36861148
-170 103.15 3.31 4.04244633
-160 113.15 6.76 9.25289948

Ln(Pvap)
-150 123.15 12.3 18.3943657 1
-140 133.15 20.54 35.8370935 0
1.4 1.3 1.2 1.1 1
-130 143.15 32.15 55.7520746 -1
-125 148.15 39.44 75.4532265
-2
-120 153.15 47.85 95.1543785
-3

1/Tx1

FORMULAS PARA LOS CALCULOS

√(∝)=1+𝜑(1−√(𝑇/𝑇𝑐) 𝑑𝑎/𝑑𝑇=−0.4572
𝑦=𝑦_1+(𝑥−𝑥_1)/(𝑥_2−𝑥_1 ) ) 𝑇_𝑐^2)/𝑃_𝑐 𝜑√(
(𝑦_2−𝑦_1)

𝑎(𝑇)=0.45724 (𝑅^2 〖𝑇𝑐〗 ▁𝐻 (𝑇,𝑃)−▁(𝐻^𝐼𝐺 ) (𝑇,𝑃)=

𝜑=0.37464+1.54226𝜔−0.26992𝜔^2 𝐿𝑛[(𝑍+(1+√2)𝐵)/(𝑍+(1−√
^2)/𝑃𝑐∝(𝑇)

𝐴=𝑎𝑃/(𝑅^2 𝐵=𝑃𝑏/𝑅𝑇
𝑇^2 )
▁𝑆 (𝑇,𝑃)−▁(𝑆^𝐼𝐺 ) (𝑇,𝑃)=𝑅
𝑏=0.07780 𝐿𝑛[(𝑍+(1+√2)𝐵)/(𝑍+(1−√2
𝑅𝑇𝑐/𝑃𝑐

𝑍^3−(1−𝐵) 𝑍^2+(𝐴−3𝐵^2−2𝐵)𝑍−(𝐴𝐵−𝐵^2−𝐵^3 )=0

𝑍=(𝑃▁𝑉)/𝑅𝑇 ∴ 1/Tx10^2 TK Ln(Pvap)

▁𝑉=𝑍𝑅𝑇/𝑃
1.4 71.4285714 -3
1.3 76.9230769 -2
1.2 83.3333333 -1
1.1 90.9090909 0
1 100 1
0.9 111.111111 2
0.8 125 3
0.7 142.857143 4
0.6 166.666667 5
por pressure of oxygen over the temperature range of −200◦C
ntropy along the vapor-liquid equilibrium phase envelope.

presion de vapor
5
4
3
2
1
0
1.2 1.1 1 0.9 0.8 0.7 0.6
-1

1/Tx10^2

𝑑𝑎/𝑑𝑇=−0.45724 (𝑅^2
𝑇_𝑐^2)/𝑃_𝑐 𝜑√(𝛼/(𝑇𝑇_𝑐 ))

𝐻 (𝑇,𝑃)−▁(𝐻^𝐼𝐺 ) (𝑇,𝑃)=𝑅𝑇(𝑍−1)+(𝑇(𝑑𝑎/𝑑𝑇)−𝑎)/(2√2𝑏)
𝑛[(𝑍+(1+√2)𝐵)/(𝑍+(1−√2)𝐵)]

(𝑇,𝑃)−▁(𝑆^𝐼𝐺 ) (𝑇,𝑃)=𝑅𝐿𝑛(𝑍−𝐵)+((𝑑𝑎/𝑑𝑇))/(2√2𝑏)
(𝑍+(1+√2)𝐵)/(𝑍+(1−√2)𝐵)]

Pvap bar
0.04978707
0.13533528
0.36787944
1
2.71828183
7.3890561
20.0855369
54.59815
148.413159
 Propiedades criticas Oxigeno
Tc K Pc bar Vc m^3/kmol Zc ꙍ R ϕ
154.600000 50.460000 0.073200 0.288000 0.021000 0.083145 0.406908
√(∝)=1+𝜑(1−√(𝑇/𝑇𝑐)
) 𝑏=0.07780 𝑎(𝑇)=0.45724 (𝑅^2
𝜑=0.37464+1.54226𝜔−0.26992𝜔^2 𝑅𝑇𝑐/𝑃𝑐 ^2)/𝑃𝑐∝(𝑇)
𝐵=𝑃𝑏/𝑅𝑇 𝐴=𝑎𝑃/(𝑅^2
𝑇^2 ) 𝑍^3−(1−𝐵) 𝑍^2+(𝐴−3𝐵^2−2𝐵)𝑍−(𝐴𝐵−𝐵^2−𝐵^3 )=0

TK Pvap Pliq α^1/2 α a

1 73.150000 0.110000 0.076589 1.127011 1.061608 1.589470
2 83.150000 0.470000 0.361229 1.108492 1.052849 1.576357
3 93.150000 1.390000 1.423556 1.091057 1.044537 1.563910
4 98.150000 2.190000 2.368611 1.082690 1.040524 1.557903
5 103.150000 3.310000 4.042446 1.074535 1.036598 1.552024
6 113.150000 6.760000 9.252899 1.058796 1.028978 1.540616
7 123.150000 12.300000 18.394366 1.043739 1.021636 1.529622
8 133.150000 20.540000 35.837094 1.029282 1.014535 1.518992
9 143.150000 32.150000 55.752075 1.015358 1.007650 1.508682
10 148.150000 39.440000 75.453227 1.008579 1.004280 1.503637
11 153.150000 47.850000 95.154378 1.001913 1.000956 1.498660
 b
0.019819

𝑎(𝑇)=0.45724 (𝑅^2 〖𝑇𝑐〗

^2)/𝑃𝑐∝(𝑇)

2−𝐵^3 )=0
SOLVE
Avap Bvap Zvap
0.004727 0.00024957 0.9955
0.015501 0.00134734 0.9856
0.036240 0.00355692 0.9664
0.051231 0.00531858 0.9524
0.069842 0.00764893 0.9347
0.117668 0.01424079 0.8878
0.179452 0.02380743 0.8236
0.254566 0.03677063 0.7380
0.342392 0.05353422 0.5892
0.390845 0.06345664 0.5413
0.442261 0.07447436 0.4181
Aliq Bliq Zliq
0.003291 0.000250 0.00249
0.011914 0.001036 0.00132
0.037115 0.003643 0.00489
0.055409 0.005752 0.00793
0.085297 0.009342 0.01328
0.161061 0.019492 0.02968
0.268366 0.035603 0.03897
0.444154 0.064155 0.05883
0.593750 0.092835 0.11553
0.747731 0.121400 0.18304
0.879479 0.148100 0.24217
0.29818
 R ϕ Tc K Pc bar τ*Raiz
0.083145 0.40690843 154.6 50.46 -0.60923509
𝜏=𝑑𝑎/𝑑𝑇=−0.45724 (𝑅^2 𝑇_𝑐^2)/𝑃_𝑐
𝑍=(𝑃▁𝑉)/𝑅𝑇 ∴ 𝜑√(𝛼/(𝑇𝑇_𝑐 ))
▁𝑉=𝑍𝑅𝑇/𝑃

T Pvap Pliq Zvap Zliq Vvap Vliq

73.15 0.11 0.07658932 0.9955 0.00249 55.04261 0.19773
83.15 0.47 0.36122868 0.9856 0.001322 14.49777 0.02530
93.15 1.39 1.42355647 0.9664 0.004886 5.38470 0.02658
98.15 2.19 2.36861148 0.9524 0.0079328 3.54896 0.02733
103.15 3.31 4.04244633 0.9347 0.013279 2.42186 0.02817
113.15 6.76 9.25289948 0.8878 0.02968 1.23555 0.03018
123.15 12.3 18.3943657 0.8236 0.03897 0.68562 0.02169
133.15 20.54 35.8370935 0.738 0.05883 0.39777 0.01817
143.15 32.15 55.7520746 0.5892 0.11553 0.21813 0.02466
148.15 39.44 75.4532265 0.5413 0.183044 0.16906 0.02988
153.15 47.85 95.1543785 0.4181 0.242169 0.11126 0.03241
 α τ
1.0616076825 -0.00590276
1.0528493558 -0.00551356
1.0445365563 -0.0051886
1.0405241227 -0.005045
1.0365977042 -0.00491191
1.028978254 -0.00467257
1.0216355081 -0.00446284
1.0145352993 -0.00427704
1.0076497794 -0.00411092
1.0042801678 -0.00403419
1.0009558937 -0.00396122
 TK Zvap Zliq τ a Bvap Bliq
73.15 0.9955 0.00249 -0.00590276 1.58946979 0.00024957 0.00024957
83.15 0.9856 0.001322 -0.00551356 1.57635657 0.00134734 0.00103553
93.15 0.9664 0.004886 -0.0051886 1.5639104 0.00355692 0.00364279
98.15 0.9524 0.0079328 -0.005045 1.55790287 0.00531858 0.00575235
103.15 0.9347 0.013279 -0.00491191 1.55202412 0.00764893 0.00934151
113.15 0.8878 0.02968 -0.00467257 1.54061606 0.01424079 0.01949239
123.15 0.8236 0.03897 -0.00446284 1.52962229 0.02380743 0.03560346
133.15 0.738 0.05883 -0.00427704 1.51899165 0.03677063 0.06415543
143.15 0.5892 0.11553 -0.00411092 1.50868245 0.05353422 0.09283495
148.15 0.5413 0.183044 -0.00403419 1.50363737 0.06345664 0.12139981
153.15 0.4181 0.242169 -0.00396122 1.49866017 0.07447436 0.14809951

▁𝐻 (𝑇,𝑃)−▁(𝐻^𝐼𝐺 ) (𝑇,𝑃)=𝑅𝑇(𝑍−1)+(𝑇(𝑑𝑎/𝑑𝑇)−𝑎)/(2√2𝑏)
b R 𝐿𝑛[(𝑍+(1+√2)𝐵)/(𝑍+(1−√2)𝐵)]
0.01981883 0.083145

▁𝑆 (𝑇,𝑃)−▁(𝑆^𝐼𝐺 ) (𝑇,𝑃)=𝑅𝐿𝑛(𝑍−𝐵)+((𝑑𝑎/𝑑𝑇))/(2√2𝑏)
𝐿𝑛[(𝑍+(1+√2)𝐵)/(𝑍+(1−√2)𝐵)]
 Hvap Hliq Sliq Svap
-0.03096781 -7.38218501 Err:502 -0.00039672
-0.11928639 -14.3339124 Err:502 -0.00132425
-0.31355822 -14.899687 Err:502 -0.00315994
-0.46943844 -15.1547243 Err:502 -0.00453808
-0.67874175 -15.3755377 Err:502 -0.00632336
-1.28765779 -15.7153557 Err:502 -0.01128938
-2.22136895 -18.4061728 Err:502 -0.01866752
-3.60567377 -20.3374604 Err:502 -0.02968124
-6.13598775 -18.3384954 -0.40812364 -0.05228299
-7.22808189 -16.8030382 -0.25126249 -0.06192905
-9.70622221 -15.989046 -0.20780446 -0.08986784

𝑎/𝑑𝑇)−𝑎)/(2√2𝑏)
-0.00590276
-0.00551356
-0.0051886
𝑑𝑇))/(2√2𝑏) -0.005045
-0.00491191
-0.00467257
-0.00446284
-0.00427704
-0.00411092
-0.00403419
-0.00396122