Tarea1 - Gas 1
Tarea1 - Gas 1
Tarea1 - Gas 1
XAVIER DE CHUQUISACA
ASIGNATURA
Condiciones de entrada
P=850 psia
T=75 °F
Componente Composición Mi Mi*Yi Pri Tri
C1 0.8655 16.04 13.88262 1.26577566 1.55762396
C2 0.0697 30.07 2.095879 1.20093136 0.97262243
C3 0.0324 44.1 1.42884 1.37895232 0.80345926
IC4 0.0098 58.12 0.569576 1.59253569 0.72857221
nC4 0.0044 58.12 0.255728 1.55865494 0.69792972
iC5 0.0015 72.15 0.108225 1.75465396 0.64419625
nC5 0.0017 72.15 0.122655 1.74420871 0.63199764
C6 0.0006 86.18 0.051708 1.94057724 0.58483735
C7 0.0002 100.2 0.02004 2.14202913 0.54976659
CO2 0.0011 44.01 0.048411 0.79303695 0.9764592
N2 0.0131 28.01 0.366931 1.72877371 2.35371544
𝑀𝑔𝑎𝑠 = ∑ 𝑋𝑖 𝑀𝑖 = 18.951
𝑖=1
𝑀𝑔𝑎𝑠
𝑆𝐺𝑔 = = 0.654
28.97
2
XH=0 𝑍 = 1 − (∑ 𝑋𝑖 ∗ √𝑏𝑖 ) = 0.9973
𝐻𝑜
𝐻𝑟 = = 1144.839
𝑍
𝑷 850
𝑷𝒓𝑪𝟏 = = = 1.266
𝑷𝒄𝑪𝟏 671.525
𝑻𝒐𝒑𝒆𝒓𝒂𝒄𝒊𝒐𝒏 534.67
𝑻𝒓𝑪𝟏 = = = 1.5576
𝑻𝒄𝑪𝟏 343.26
Componente w mi A B
C1 0.011 0.3915722 0.81517654 0.19445871 0.0632228
C2 0.099 0.52467825 1.01451641 0.58888829 0.09606242
C3 0.153 0.60428722 1.12918039 1.10288205 0.13352574
IC4 0.183 0.64783423 1.19873153 1.64441003 0.17005765
nC4 0.199 0.67086064 1.23300811 1.80400087 0.17374723
iC5 0.227 0.71082431 1.30029244 2.51385577 0.2119107
nC5 0.251 0.74474203 1.32868256 2.652974 0.2147151
C6 0.299 0.81164462 1.41834324 3.67947872 0.25815196
C7 0.351 0.88271885 1.50851443 4.88835084 0.3031284
CO2 0.239 0.72782204 1.01730982 0.38688569 0.06318572
N2 0.039 0.43437759 0.58976742 0.08415021 0.0571431
Pri
𝑨𝑪𝟏 = 0.45724 ∗ 0.81517654 ∗ = 0.790493
Tri2
PrC1
𝑩𝑪𝟏 = 0.0778 ∗ = 0.23930736
Tr C1
𝑵 𝑵
𝑨 = ∑ ∑ 𝑿𝒊 𝑿𝒋 √𝑨𝒊 𝑨𝒋 = 0.2489
𝒊=𝟏 𝒋=𝟏
𝑩 = ∑ 𝑿𝒊 𝑩𝒊 = 0.0699
𝒊=𝟏
𝒁𝟑 − (𝟏 − 𝑩)𝒁𝟐 + (𝑨 − 𝟐𝑩 − 𝟑𝑩𝟐 )𝒁 − (𝑨𝑩 − 𝑩𝟐 − 𝑩𝟑 ) = 𝟎
Z= 0.834
R=10.731 psi*ft3/lbmol*R
𝑃𝑀𝑔𝑎𝑠
𝜌𝑔𝑎𝑠 = = 3.366 𝑙𝑏/𝑓𝑡 3
𝑍𝑅𝑇
Viscosidad por graficas
P=850 T=75 °F SGg=0.654 De graficas 𝜇𝑔𝑎𝑠 = 0.0122
Ppc 667.6565341
Ppr 1.273109685
Tpc 374.53482
Tpr 1.427557523
Lee Gonzales
986
𝑋 = 3.5 + + 0.01 ∗ 𝑀𝑔 = 5.5336
𝑇
𝑌 = 2.4 − 0.2 ∗ 𝑋 = 1.2933
10−4 ∗ (9.4 + 0.02 ∗ 𝑀𝑔𝑎𝑠) ∗ 𝑇 1.5
𝐾= = 0.01095
209 + 19 ∗ 𝑀𝑔 + 𝑇
𝜌𝑔𝑎𝑠 𝑌
𝜇𝑔𝑎𝑠 = 𝐾 ∗ 𝐸𝑋𝑃 (𝑋 ∗ ( ) ) = 0.0124
62.4
𝑵 𝑵 𝑻 𝑻
𝒎𝒊 √ 𝑻𝒄𝒊 𝒎𝒋 √ 𝑻𝒄𝒋
𝑫 = ∑ ∑ 𝑿𝒊 𝑿𝒋 √𝑨𝒊 𝑨𝒋 ∗ ∗ + ∗ = 𝟎. 𝟏𝟑𝟑𝟕
𝟐 𝜶𝒊 𝟐 𝜶𝒋
𝒊=𝟏 𝒋=𝟏
[ ]
Q=3E7 ft3/dia=3.47E2 ft3/s
Condición estándar P=14.7 y T=60°F =>Zo=0.997
𝑃𝑀𝑔𝑎𝑠 Lbm
𝜌𝑜 = = 0.05009 3
𝑍𝑅𝑇 𝑓𝑡
𝑚𝑎𝑠𝑎 𝑓𝑡 3 Lbm
𝜌𝑔𝑎𝑠 = => 𝑚𝑎𝑠𝑎 = 𝜌𝑜 ∗ 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 = 3.47𝐸2 ∗ 0.05009
𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑠 𝑓𝑡 3
𝐹 = 𝑚𝑎𝑠𝑎 = 1.74E1 lbm /s
𝑚
𝑛= = 9.18𝐸 − 1 𝑙𝑏𝑚𝑜𝑙
𝑀𝑔𝑎𝑠
COMPONENTE Yi A B C D H0i
C1 0.8655 19.25 0.05213 1.197E-05 -1.132E-08 -34.1635398
C2 0.0697 5.409 0.1781 -6.938E-05 8.713E-09 -4.06550329
C3 0.0324 -4.224 0.3063 -1.586E-04 3.215E-08 -2.65429784
IC4 0.0098 -1.39 0.3847 -1.846E-04 2.895E-08 -1.06176396
nC4 0.0044 9.487 0.3313 -1.108E-04 -2.822E-09 -0.48005434
iC5 0.0015 -9.525 0.5066 -2.729E-04 5.723E-08 -0.19762661
nC5 0.0017 -3.626 0.4873 -2.580E-04 5.305E-08 -0.22655486
C6 0.0006 -4.413 0.582 -3.119E-04 6.494E-08 -0.09525057
C7+ 0.0002 -5.146 0.6762 -3.651E-04 7.658E-08 -0.03683811
CO2 0.0011 19.8 0.07344 -5.602E-05 1.715E-08 -0.04540074
N2 0.0131 31.15 -0.01357 2.680E-05 -1.168E-08 -0.42469158
𝑇
𝐻0 = ∑ (𝑋𝑖 ∗ ∫ 𝐶𝑝𝑖 ∗ 𝑑𝑇)
𝑇𝑜
𝐵𝑖 𝐶𝑖 𝐷𝑖
𝐻0 = ∑ (𝑋𝑖 ∗ (𝐴𝑖 ∗ (𝑇 − 𝑇𝑜 ) + ∗ (𝑇 2 − 𝑇𝑜 2 ) + ∗ (𝑇 3 − 𝑇𝑜 3 ) + ∗ (𝑇 4 − 𝑇𝑜 4 )))
2 3 4
𝐽 𝐵𝑇𝑈
𝐻0 = −43.4515 = −18.6808
𝑚𝑜𝑙 𝑙𝑏𝑚𝑜𝑙
𝐴 𝐷 𝑍 + 𝐵 ∗ (1 + √2)
(𝐻0 − 𝐻) = 𝑅𝑇 ∗ [1 − 𝑍 + (1 + ) ∗ 𝐿𝑛 ( )] = 627.192
√8𝐵 𝐴 𝑍 + 𝐵 ∗ (1 − √2)
Calculo a la salida
Tsalida=40 °F
Psalida=830 psia
Componente Composición Pri Tri A B
C1 0.8655 1.2359927 1.45566043 0.22532057 0.06605952
C2 0.0697 1.17267415 0.90895365 0.68112028 0.10037261
C3 0.0324 1.34650638 0.75086406 1.27528235 0.13951686
IC4 0.0098 1.55506426 0.68087919 1.90087862 0.17768791
nC4 0.0044 1.5219807 0.65224259 2.08531479 0.18154304
iC5 0.0015 1.71336799 0.60202655 2.90502629 0.22141885
nC5 0.0017 1.70316851 0.59062648 3.06846452 0.22434909
C6 0.0006 1.8949166 0.54655335 4.25783739 0.2697349
C7 0.0002 2.09162845 0.51377835 5.66105862 0.31672937
CO2 0.0011 0.77437726 0.91253927 0.45333678 0.06602078
N2 0.0131 1.68809668 2.19963902 0.09959822 0.05970703
Amezcla=0.2883
Bmezcla=0.07302
Dmezcla=0.1471
Z=0.792
COMPONENTE Yi H0i
C1 0.8655 -622.301363
C2 0.0697 -73.2624993
C3 0.0324 -47.6658732
IC4 0.0098 -19.0851648
nC4 0.0044 -8.64700587
iC5 0.0015 -3.54747202
nC5 0.0017 -4.07213415
C6 0.0006 -1.71213254
C7+ 0.0002 -0.66219139
CO2 0.0011 -0.82997738
N2 0.0131 -7.85926252
𝐴 𝐷 𝑍 + 𝐵 ∗ (1 + √2)
(𝐻0 − 𝐻) = 𝑅𝑇 ∗ [1 − 𝑍 + (1 + ) ∗ 𝐿𝑛 ( )] = 708.605
√8𝐵 𝐴 𝑍 + 𝐵 ∗ (1 − √2)
𝐵𝑖 𝐶𝑖 𝐷𝑖
𝐻0 = ∑ (𝑋𝑖 ∗ (𝐴𝑖 ∗ (𝑇 − 𝑇𝑜 ) + ∗ (𝑇 2 − 𝑇𝑜 2 ) + ∗ (𝑇 3 − 𝑇𝑜 3 ) + ∗ (𝑇 4 − 𝑇𝑜 4 )))
2 3 4
𝐽 𝐵𝑇𝑈
𝐻0 = −789.645 = −339.486
𝑚𝑜𝑙 𝑙𝑏𝑚𝑜𝑙
𝐻𝑠𝑎𝑙𝑖𝑑𝑎 = 𝐻0 − (𝐻0 − 𝐻) = −1048.092
𝑄 = −(𝐻𝑒𝑛𝑡 − 𝐻𝑠𝑎𝑙 ) ∗ 𝑛 = −3.69𝐸2
Sistema Compresor
Calculo de entropía
Amezcla=0.2883
Bmezcla=0.07302
Dmezcla=0.1471
Z=0.792
𝑃𝑜 𝐷 𝑍 + 𝐵 ∗ (1 − √2)
(𝑆 0 − 𝑆) = −𝑅 ∗ [𝐿𝑁(𝑍 − 𝐵) + 𝐿𝑁 ( )+ ∗ 𝐿𝑛 ( )] = 9.0103
𝑃 √8𝐵 𝑍 + 𝐵 ∗ (1 + √2)
COMPONENTE Yi S0i
C1 0.8655 -1.12274575
C2 0.0697 1.28906578
C3 0.0324 0.75831119
IC4 0.0098 0.31058402
nC4 0.0044 0.16846905
iC5 0.0015 0.06877139
nC5 0.0017 0.07599204
C6 0.0006 0.0310612
C7+ 0.0002 0.01186281
CO2 0.0011 0.05942235
N2 0.0131 0.44487062
𝑇
𝐶𝑝𝑖
𝑆 0 = ∑ (𝑋𝑖 ∗ (∫ ∗ 𝑑𝑇 − 𝑅𝐿𝑁(𝑋𝑖)))
𝑇𝑜 𝑇
𝐶𝑖 𝐷𝑖
𝑆 0 = ∑ (𝑋𝑖 ∗ (𝐴𝑖 ∗ (𝑙𝑛𝑇 − 𝑙𝑛𝑇𝑜 ) + 𝐵𝑖 ∗ (𝑇 − 𝑇𝑜 ) + ∗ (𝑇 2 − 𝑇𝑜 2 ) + ∗ (𝑇 3 − 𝑇𝑜 3 )))
2 3
𝐽 𝐵𝑇𝑈
𝑆 0 = 2.0957 = 0.5005
𝑚𝑜𝑙𝐾 𝑙𝑏𝑚𝑜𝑙𝑅
𝑆 = 𝑆 0 − (𝑆 0 − 𝑆) = −8.5097
Algoritmo
-Suponer T salida
Calcular S0
Calcular S salida
Sent=Ssal
Si Fin
Condición de salida
Tsalida =-48.2F
Psalida=400 psia
Amezcla=0.2247
Bmezcla=0.04274
Dmezcla=0.0993
Z=0.8026
𝑃𝑜 𝐷 𝑍 + 𝐵 ∗ (1 − √2)
(𝑆 0 − 𝑆) = −𝑅 ∗ [𝐿𝑁(𝑍 − 𝐵) + 𝐿𝑁 ( )+ ∗ 𝐿𝑛 ( )] = 7.344
𝑃 √8𝐵 𝑍 + 𝐵 ∗ (1 + √2)
COMPONENTE Yi S0i
C1 0.8655 -6.6669241
C2 0.0697 0.66564109
C3 0.0324 0.35905641
IC4 0.0098 0.15003889
nC4 0.0044 0.09501896
iC5 0.0015 0.0391174
nC5 0.0017 0.04175251
C6 0.0006 0.01666279
C7+ 0.0002 0.00629328
CO2 0.0011 0.05193881
N2 0.0131 0.37045797
𝑇
𝐶𝑝𝑖
𝑆 0 = ∑ (𝑋𝑖 ∗ (∫ ∗ 𝑑𝑇 − 𝑅𝐿𝑁(𝑋𝑖)))
𝑇𝑜 𝑇
𝐶𝑖 𝐷𝑖
𝑆 0 = ∑ (𝑋𝑖 ∗ (𝐴𝑖 ∗ (𝑙𝑛𝑇 − 𝑙𝑛𝑇𝑜 ) + 𝐵𝑖 ∗ (𝑇 − 𝑇𝑜 ) + ∗ (𝑇 2 − 𝑇𝑜 2 ) + ∗ (𝑇 3 − 𝑇𝑜 3 )))
2 3
𝐽 𝐵𝑇𝑈
𝑆 0 = −4.871 = −1.1634
𝑚𝑜𝑙𝐾 𝑙𝑏𝑚𝑜𝑙𝑅
𝑆 = 𝑆 0 − (𝑆 0 − 𝑆) = −8.508
COMPONENTE Yi H0i
C1 0.8655 -2022.99062
C2 0.0697 -230.9454
C3 0.0324 -148.690895
IC4 0.0098 -59.7040151
nC4 0.0044 -27.2255483
iC5 0.0015 -11.0513423
nC5 0.0017 -12.7350656
C6 0.0006 -5.35505554
C7+ 0.0002 -2.07132807
CO2 0.0011 -2.72015519
N2 0.0131 -26.6325635
𝐴 𝐷 𝑍 + 𝐵 ∗ (1 + √2)
(𝐻0 − 𝐻) = 𝑅𝑇 ∗ [1 − 𝑍 + (1 + ) ∗ 𝐿𝑛 ( )] = 475.311
√8𝐵 𝐴 𝑍 + 𝐵 ∗ (1 − √2)
𝐵𝑖 𝐶𝑖 𝐷𝑖
𝐻0 = ∑ (𝑋𝑖 ∗ (𝐴𝑖 ∗ (𝑇 − 𝑇𝑜 ) + ∗ (𝑇 2 − 𝑇𝑜 2 ) + ∗ (𝑇 3 − 𝑇𝑜 3 ) + ∗ (𝑇 4 − 𝑇𝑜 4 )))
2 3 4
𝐽 𝐵𝑇𝑈
𝐻0 = −2550.122 = −1096.355
𝑚𝑜𝑙 𝑙𝑏𝑚𝑜𝑙
𝐻𝑠𝑎𝑙𝑖𝑑𝑎 = 𝐻0 − (𝐻0 − 𝐻) = −1571.666
Hentrada=−1048.092
Tablas.
R=1.987
𝑾𝒎𝒆𝒛𝒄𝒍𝒂 = ∑ 𝒘𝒊 ∗ 𝒚𝒊 = 𝟎. 𝟎𝟐𝟓𝟖
Pr 1.273109685
Tr 1.427557523
De tablas
(𝟎)
𝑯𝟎 − 𝑯
( ) = 𝟎. 𝟕𝟑𝟑
𝑹𝑻𝒄
(𝟏)
𝑯𝟎 − 𝑯
( ) = 𝟎. 𝟎𝟗𝟑𝟓
𝑹𝑻𝒄
𝐽 𝐵𝑇𝑈
𝐻0 = −43.4515 = −18.6808
𝑚𝑜𝑙 𝑙𝑏𝑚𝑜𝑙
(𝟎) (𝟏)
0
𝑯𝟎 − 𝑯 𝑯𝟎 − 𝑯
(𝐻 − 𝐻) = 𝑅𝑇𝑐 ∗ [( ) +𝒘∗( ) ] = 547.295
𝑹𝑻𝒄 𝑹𝑻𝒄
Condición de salida.
Pr 1.243154163
Tr 1.334108268
(𝟎)
𝑯𝟎 − 𝑯
( ) = 𝟎. 𝟖𝟒𝟔
𝑹𝑻𝒄
(𝟏)
𝑯𝟎 − 𝑯
( ) = 𝟎. 𝟐𝟎𝟐
𝑹𝑻𝒄
𝐽 𝐵𝑇𝑈
𝐻0 = −789.645 = −339.486
𝑚𝑜𝑙 𝑙𝑏𝑚𝑜𝑙
(𝟎) (𝟏)
0
𝑯𝟎 − 𝑯 𝑯𝟎 − 𝑯
(𝐻 − 𝐻) = 𝑅𝑇𝑐 ∗ [( ) +𝒘∗( ) ] = 633.472
𝑹𝑻𝒄 𝑹𝑻𝒄
Sistema Compresor.
(𝟎)
𝑺𝟎 − 𝑺
( ) = 𝟎. 𝟒𝟔𝟏
𝑹
(𝟏)
𝑺𝟎 − 𝑺
( ) = 𝟎. 𝟐𝟑𝟐
𝑹
(𝟎) (𝟏)
𝑆0 − 𝑆 𝑃0 𝑺𝟎 − 𝑺 𝑺𝟎 − 𝑺
= −𝐿𝑁 +( ) +𝒘∗( )
𝑅 𝑃 𝑹 𝑹
(𝟎) (𝟏)
𝑆0 − 𝑆 1 𝑏𝑎𝑟 𝑺𝟎 − 𝑺 𝑺𝟎 − 𝑺
= −𝐿𝑁 +( ) +𝒘∗( ) = 𝟒. 𝟓𝟏𝟒
𝑅 𝑃 𝑹 𝑹
𝑆0 − 𝑆
𝑺𝒆𝒏𝒕𝒓𝒂𝒅𝒂 = 𝑆 0 − 𝑅 ∗ [ ] = −8.47
𝑅
Condición de salida
Algoritmo
-Suponer T salida
Calcular S0
Calcular S salida
Sent=Ssal
Si Fin
(𝟎)
𝑺𝟎 − 𝑺
( ) = 𝟎. 𝟑𝟕𝟏
𝑹
(𝟏)
𝑺𝟎 − 𝑺
( ) = 𝟎. 𝟑𝟓𝟎
𝑹
(𝟎) (𝟏)
𝑆0 − 𝑆 𝑃0 𝑺𝟎 − 𝑺 𝑺𝟎 − 𝑺
= −𝐿𝑁 +( ) +𝒘∗( )
𝑅 𝑃 𝑹 𝑹
(𝟎) (𝟏)
𝑆0 − 𝑆 1 𝑏𝑎𝑟 𝑺𝟎 − 𝑺 𝑺𝟎 − 𝑺
= −𝐿𝑁 +( ) +𝒘∗( ) = 𝟑. 𝟔𝟗𝟕
𝑅 𝑃 𝑹 𝑹
𝑆0 − 𝑆
𝑺𝒔𝒂𝒍𝒊𝒅𝒂 = 𝑆 0 − 𝑅 ∗ [ ] = −8.47
𝑅
(𝟎)
𝑯𝟎 − 𝑯
( ) = 𝟎. 𝟓𝟖𝟏
𝑹𝑻𝒄
(𝟏)
𝑯𝟎 − 𝑯
( ) = 𝟎. 𝟑𝟖𝟏
𝑹𝑻𝒄
𝐵𝑇𝑈
𝐻0 = −1080.015
𝑙𝑏𝑚𝑜𝑙
(𝟎) (𝟏)
0
𝑯𝟎 − 𝑯 𝑯𝟎 − 𝑯
(𝐻 − 𝐻) = 𝑅𝑇𝑐 ∗ [( ) +𝒘∗( ) ] = 439.696
𝑹𝑻𝒄 𝑹𝑻𝒄