Termodinámica I
Fugacidad y Coeficiente de
fugacidad de gases puros
MSc. Pablo Díaz Bravo
Fugacidad y coeficiente de fugacidad
La fugacidad de un gas real es la presión parcial efectiva que reemplaza la presión parcial mecánica.
Se sabe de la ecuación (3): 𝑑𝐺=𝑣𝑑𝑃− 𝑠𝑑𝑇
A T constante 𝑑𝐺=𝑣𝑑𝑃
∗ ∗ 𝑑𝑃
Para un gas ideal: 𝑑 𝐺 =𝑣 𝑑𝑃 =𝑅𝑇 = 𝑅𝑇𝑑𝑙𝑛𝑃 (1)
𝑃
(2)
Para un gas real: 𝑑𝐺=𝑣𝑑𝑃=𝑅𝑇𝑑𝑙𝑛𝑓
Donde f es la fugacidad
Haciendo (2)-(1) 𝑑 ( 𝐺 − 𝐺∗ ) =( 𝑣 −𝑣 ∗ ) 𝑑𝑃= 𝑅𝑇𝑑𝑙𝑛 ( )
𝑓
𝑃
(3)
El termino (G-G*) se denomina energía residual de Gibbs, y a la relación (f/P) coeficiente de
fugacidad.
𝐺 𝑅 =( 𝐺 − 𝐺∗ ) (4)
𝜑= ( )
𝑓
𝑃
(5)
En general, las propiedades residuales se aplican para gases y miden la
desviación de la propiedad real respecto a la propiedad del gas ideal.
𝑀 𝑅 =𝑀 − 𝑀 ∗ (6)
Así tenemos:
De la ecuación (3):
(7)
Luego,
Integrando,
(8)
Usando la ecuación Virial
Truncando la ecuación en dos términos:
Reemplazando en (8) e integrando.
𝑙𝑛 𝜑=
𝐵𝑃
𝑅𝑇
=( )
𝐵 𝑃𝑐 𝑃𝑟
𝑅𝑇𝑐 𝑇𝑟
(9)
Usando Ecuaciones Cubicas
Por ejemplo: La ecuación de Van der Waals:
Se sabe,
Derivando,
Reemplazando en (8)
(10)
Resolveremos el termino integral usando la ecuación de Van der Waals.
RT a
P 2
v b v
(11)
Reemplazando en (33)
Usando condiciones de contorno
Para gas ideal,
Reemplazando en (11)
Luego,
(12)
Para la ecuación de Redlich-Kwong
(13)
Para la ecuación de Peng-Robinson
(14)
Usando Correlaciones Generalizadas
Método de 2 parámetros
De (8)
Introduciendo las propiedades reducidas
(15)
El coeficiente de fugacidad es tabulado en cartas generalizadas.
Método de 3 parámetros
De la ecuación Virial truncado en dos términos,
Usando la correlación de Pitzer,
Luego
Reemplazando en la ecuación (15)
Entonces,
(16)
(17)
Donde:
Ejercicio 1
Determine el coeficiente de fugacidad para 1-buteno a 200 ºC y 70 bar. usando:
a) La ecuación virial
b) La ecuación de Van der Waals Van der Waals
c) La ecuación de Redlich-Kwong a 12723267.8
d) La ecuación de Soave R-K b 107.960673
e) Correlaciones Generalizadas con 2 parámetros A* 0.57590969
f) Correlaciones generalizadas con 3 parámetros B* 0.19217304
Solución Z 0.466
Propiedades criticas F(obj) 6.7675E-06
Ecuacion virial v 261.793607
1-BUTANO
lnPhi -0.4745988
T (ºC) 200 -0.26592613 Phi 0.6221346
P (bar) 70 0.03458601
Tc 420 BPc/RTc -0.2593202 𝒁 𝟑 − ( 𝟏+ 𝑩 ∗ ) 𝒁 𝟐+ 𝑨 ∗ 𝒁 − 𝑨 ∗ 𝑩∗=𝟎
Pc 40.43
w 0.191 lnphi -0.39867481 Redlich-Kwong
R 83.14 Phi 0.67120894 a 12148558.67
b 74.82970151
Tr 1.12619048 A* 0.54989589
Pr 1.73138758 B* 0.133198978
Videal 561.788857 Z3-Z2+(A-B-B2)Z-AB=0
Z 0.492
F(obj) 7.17502E-05
v 276.4001177
lnPhi -0.47213622
Phi 0.62366855
c) Correlaciones Generalizadas con 2 parámetros Tr 1.12619048
Pr 1.73138758
Videal 561.788857
De grafico 𝝋=𝟎 ,𝟔𝟐
d) Correlaciones Generalizadas con 3 parámetros
𝜑 𝑜=0.61 𝜑 1=1,093
𝜑=0.61 (1,093 )0.191 =𝟎 , 𝟔𝟐𝟎𝟓
Ejercicio 1
Determine el coeficiente de fugacidad para Isobutano a 200 °C y 50 bar
usando:
a) La ecuación virial
b) La ecuación de Van der Waals
c) La ecuación de Redlich-Kwong
d) La ecuación de Peng-Robinson
e) Correlaciones Generalizadas con 2 parámetros
f) Correlaciones generalizadas con 3 parámetros
Ejercicio 2
Determine el coeficiente de fugacidad para etano a 200 °C y 50 bar usando:
a) La ecuación virial
b) La ecuación de Van der Waals
c) La ecuación de Redlich-Kwong
d) La ecuación de Peng-Robinson
e) Correlaciones Generalizadas con 2 parámetros
f) Correlaciones generalizadas con 3 parámetros
