FACULTAD NACIONAL DE INGENIERIA
INGENIERIA QUIMICA –ALIMENTOS
FUGACIDAD
Gases:
𝑓𝑖 = 𝑃 ∗ ∅𝑖 (𝑃𝑢𝑟𝑜) 𝑓̂𝑖 = 𝑦𝑖 ∗ ∅̂𝑖 ∗ 𝑃(𝑀𝑒𝑧𝑐𝑙𝑎 𝑜 𝑆𝑜𝑙𝑢𝑐𝑖𝑜𝑛)
Líquidos:
𝑠𝑎𝑡 𝑉𝐿𝑠𝑎𝑡 (𝑃 − 𝑃𝑖𝑠𝑎𝑡 )
𝑓𝑖 = ∅𝑖 ∗ 𝑃𝑖𝑠𝑎𝑡 ∗ exp ( )
𝑅𝑇
*P de saturación es sinónimo de Presión de Vapor.
Determinación del coeficiente de fugacidad:
Experimentalmente:
𝑃
𝑑𝑃
𝐿𝑛 ∅𝑖 = ∫ (𝑍𝑖 − 1) ∗
0 𝑃
Correlación generalizada:
𝑃𝑟 0
𝐿𝑛(∅𝑖 ) = [ (𝐵 + 𝑤𝐵1 )] ; ∅𝑖 = (∅0 )(∅1 )𝑤
𝑇𝑟
Ecuacion Cubica:
𝐿𝑛 ∅𝑖 = 𝑍𝑖 − 1 − 𝐿𝑛(𝑍𝑖 − 𝛽𝑖) − 𝑞𝑖 ∗ 𝐼𝑖
[Link] REA M. PRQ 2221
FACULTAD NACIONAL DE INGENIERIA
INGENIERIA QUIMICA –ALIMENTOS
Ejercicio N.11 Para la siguiente información recopilada de los factores de compresibilidad de
Isobutano a 360 °K y 7 bar
a) Calcule el valor de fugacidad por el método de las correlaciones generalizadas (Lee- Keesler y Pitzer)
b) Compare los valores de fugacidad obtenidos anteriormente con un modelo experimental; para este fin se
tiene los siguientes datos experimentales:
P[Bar] 0,1 0,5 2 4 6 8 10 12 14 15,41
Z 0,997 0,989 0,964 0,936 0,907 0,876 0,841 0,801 0,755 0,717
c) Compare el valor de fugacidad obtenido anteriormente con uso de la ecuación cubica.
a1) Correlación generalizada Lee Keesler fi=P*oi
∅𝑖 = (∅0 )(∅1 )𝑤
Recopilando datos: Isobutano
Tc=408.1 °K Pc=36,48 bar w=0,181
𝑃 7
𝑃𝑟 = = = 0.192
𝑃𝑐 36.48
𝑇 360
𝑇𝑟 = = = 0.882
𝑇𝑐 408.1
[Link] REA M. PRQ 2221
FACULTAD NACIONAL DE INGENIERIA
INGENIERIA QUIMICA –ALIMENTOS
Entrando al apéndice pag 708,709
Interpolando:
𝜙 0 = 0.9077 ; 𝜙 1 = 0.94
Reemplazando:
∅𝑖 = (0,908)(0,94)0,181 = 0,899
𝑓𝑖 = 𝑃 ∗ ∅𝑖 = 7 ∗ 0,9 = 6,3 [𝐵𝑎𝑟]
a2) Correlación generalizada Pitzer
𝑃𝑟 0
𝐿𝑛(∅𝑖 ) = [ (𝐵 + 𝑤𝐵1 )]
𝑇𝑟
0,422 0,172
𝐵0 = 0,083 − = −0.433 ; 𝐵 1
= 0,139 − = −0.152
0.8821,6 0.8824,2
0,192
∅𝑖 = 𝑒𝑥𝑝 [ (−0,433 + 0,181 ∗ (−0,152))] = 0.905
0,882
𝑓𝑖 = 𝑃 ∗ ∅𝑖 = 7 ∗ 0,905 = 6,335𝐵𝑎𝑟
b)
𝑃 𝑑𝑃
𝐿𝑛 ∅𝑖 = ∫0 (𝑍𝑖 − 1) ∗
𝑃
[Link] REA M. PRQ 2221
FACULTAD NACIONAL DE INGENIERIA
INGENIERIA QUIMICA –ALIMENTOS
𝑓𝑖 = 𝑃 ∗ ∅𝑖
𝑍𝑖−1
Regresionando (EXCEL): 𝑣𝑠 𝑃
𝑃
𝑍𝑖−1
=-3E-06x4 + 0,0001x3 - 0,0014x2 + 0,0075x - 0,0283
𝑃
R² =0,9085
𝑃
𝑑𝑃
𝐿𝑛 ∅𝑖 = ∫ (𝑍𝑖 − 1) ∗
0 𝑃
7
𝐿𝑛∅𝑖 = ∫ (−3 ∗ 10−6 P 4 + 0,0001P 3 − 0,0014P 2 + 0,0075P − 0,0283 ) ∗ 𝑑𝑃
0
𝐿𝑛∅𝑖 =-0.124
∅𝑖 = exp(−0,124) = 0.88
𝑓𝑖 = 7 ∗ 0,88 = 6.16 [𝐵𝑎𝑟]
c)
𝐿𝑛 ∅𝑖 = 𝑍𝑖 − 1 − 𝐿𝑛(𝑍𝑖 − 𝛽𝑖) − 𝑞𝑖 ∗ 𝐼𝑖
R-K
P = 7 bar y 360 °K isobutano Redlich Kwong: Z=0.8991 ;β=0.0188; q=5,955 Ii=
[Link] REA M. PRQ 2221
FACULTAD NACIONAL DE INGENIERIA
INGENIERIA QUIMICA –ALIMENTOS
𝑍−𝛽 𝛼(𝑇𝑟) ∗ 𝜓 𝑃𝑟
𝑍 = 1+β−q∗β∗ ;𝑞 = ; 𝛽 =Ω∗
(𝑍 + 𝜖𝛽)(𝑍 + 𝜎𝛽) Ω ∗ 𝑇𝑟 𝑇𝑟
1
Tabla 3.1: RK: 𝜖 = 0 ; 𝜎 = 1 ; 𝛼(𝑇𝑟) = 𝑇𝑟 −2 ; Ω = 0.08664 ; 𝜓 = 0.42748
Ec 6,65 b : Z=0.8991; β=0.0188; q=5,955
1 𝑍 + 𝜎𝛽
𝐼= ln ( ) = 2.069 ∗ 10−2
𝜎−𝜖 𝑍 + 𝜖𝛽
𝐿𝑛 ∅𝑖 = 0.899 − 1 − 𝐿𝑛(0.899 − 0.0188) − 5.955 ∗ 2,069 ∗ 10−2 = −9,66 ∗ 10−2
∅𝑖 = exp(9,66 ∗ 10−2 ) = 0.907
𝑓 = 0.907 ∗ 7 = 6,35 [𝐵𝑎𝑟]
Ejercicio N.12 Calcule la fugacidad del Isobutileno líquido a su temperatura de ebullición y 200 [Bar].
𝑠𝑎𝑡 𝑉𝐿𝑠𝑎𝑡 (𝑃 − 𝑃𝑖𝑠𝑎𝑡 )
𝑓𝑖 = ∅𝑖 ∗ 𝑃𝑖𝑠𝑎𝑡 ∗ exp ( )
𝑅𝑇
DATOS Para hallar volumen liquido saturado (Con Ecuacion de Racket):
2
𝐿𝑠𝑎𝑡 (1−𝑇𝑟)7
TN =266, 3°K (Apendice) 𝑉 = 𝑉𝑐 ∗ 𝑍𝑐
2 2
𝐿𝑠𝑎𝑡 (1−𝑇𝑟 )7 (1−0,637 )7
P=200[Bar] 𝑉 = 𝑉𝑐 ∗ 𝑍𝑐 = 238,9 ∗ 0,275
[Link] REA M. PRQ 2221
FACULTAD NACIONAL DE INGENIERIA
INGENIERIA QUIMICA –ALIMENTOS
𝑆𝐴𝑇 𝑐𝑚3
𝑉𝐿 = 90,88 [ 𝑚𝑜𝑙 ]
Psat=1atm≈1 [Bar]
De Tablas para isobutileno: Tr=Tn/Tc=266,3/417,9=0,637
Vc=238,9 cm3/mol
Zc=0,275
w=0,194
Tc=417,9 K
Pc=40,00 Bar
Calculo coeficiente de fugacidad:
𝑃𝑟 0
∅𝑖 = ∅𝑖𝑠𝑎𝑡 = 𝑒𝑥𝑝 [ (𝐵 + 𝑤𝐵1 )]
𝑇𝑟
𝑃 𝑠𝑎𝑡 1 𝑏𝑎𝑟 𝑇𝑁 266,3
𝑃𝑟 = = = 0,025 𝑇𝑟 = = = 0.637
𝑃𝐶 40 𝐵𝑎𝑟 𝑇𝐶 417.9
0.422 0.172
𝐵0 = 0.083 − 𝐵1 = 0.139 −
𝑇𝑟 1,6 𝑇𝑟 4,2
𝐵0 = −0.785 𝐵1 = −1.004
[Link] REA M. PRQ 2221
FACULTAD NACIONAL DE INGENIERIA
INGENIERIA QUIMICA –ALIMENTOS
0,025
∅𝑖𝑠𝑎𝑡 = 𝑒𝑥𝑝 [ ((−0,785) + 0,194(−1,004))] = 0.962
0,637
Finalmente:
𝑠𝑎𝑡 𝑉𝐿𝑠𝑎𝑡 (𝑃 − 𝑃𝑖𝑠𝑎𝑡 )
𝑓𝑖 = ∅𝑖 ∗ 𝑃𝑖𝑠𝑎𝑡 ∗ exp ( )
𝑅𝑇
𝑐𝑚3
90,88 [ ] ∗ (200 − 1)[𝑏𝑎𝑟]
𝑚𝑜𝑙
𝑓𝑖 = 0,962 ∗ 1[𝐵𝑎𝑟] ∗ exp ( )
𝑐𝑚3 ∗ 𝑏𝑎𝑟
83,14 [ ] ∗ 266,3 [𝐾]
𝑚𝑜𝑙 ∗ 𝐾
𝑓𝑖 = 2.18 [𝐵𝑎𝑟]
[Link] REA M. PRQ 2221
