Problemas gases reales

1. Los valores notificados para los coeficientes viriales del vapor de isopropanol a 200ºC
son:
B=-388 cm3mol-1
C=-26000 cm6mol-2
Calcule V y Z para el vapor de isopropanol a 200ºC y 10 bar con los métodos siguientes:
a) La ecuación del gas ideal: PV=RT
b) Ecuación: Z=P*V/(R*T)=1+B*P/(R*T)
c) Ecuación: Z=P*V/(R*T)=1+B/V+C/V2

Donde:
V: volumen molar del gas o vapor

T=200+273.15=473.15 K
R=83.14 cm3*bar*mol-1*K-1

Solución:

a) V=R*T/P
V=83.14 cm3*bar*mol-1*K-1*473.15 K/10 bar=3934 cm3/mol, este es el volumen
ideal del gas

b) Z=P*V/(R*T)=1+B*P/(R*T)
P*V/(R*T)=1+B*P/(R*T)

Simplificando obtenemos:
V=R*T/P+B, este es el volumen real del gas
𝑇 bar 473.15𝐾 𝑐𝑚3 𝑐𝑚3
𝑉 = 𝑅 ∗ + 𝐵 = 83.14 𝑐𝑚3 ∗ ∗ + (−388 ) = 3546
𝑃 mol ∗ K 10𝑏𝑎𝑟 𝑚𝑜𝑙 𝑚𝑜𝑙
Ahora calculamos el factor de comprensibilidad (Z):
𝑉
𝑉 𝑉𝑟𝑒𝑎𝑙 3546
𝑍=𝑃∗ = ∗𝑇 =
𝑅 = = 0.9014
𝑅∗𝑇 𝑃 𝑉𝑖𝑑𝑒𝑎𝑙 3934
c) Para facilitar el proceso de iteración, se escribe la ecuación de la siguiente
forma:
Vi+1=(R*T/P)*(1+B/Vi+C/Vi2)
Si i=0, la primera iteración:
V1=(R*T/P)*(1+B/V0+C/V02)
Donde V0 es una estimación inicial del volumen molar. Para obtener éste se hace
uso del valor que corresponde al gas ideal:
P= 10 bar
T= 200
R= 83.14 cm3*bar*mol-1*K-1
B= -388 cm3mol-1
C= -26000 cm6mol-2

1
 i Vi
0 3933.7691 cm3mol-1
1 3539.15966 cm3mol-2
2 3494.34237 cm3mol-3
3 3488.60035 cm3mol-4
4 3487.85383 cm3mol-5
5 3487.75658 cm3mol-6

𝑉𝑟𝑒𝑎𝑙 3487.75658
𝑍= = = 0.88661
𝑉𝑖𝑑𝑒𝑎𝑙 3934

3933.7691−3487.75658
% de diferencia= ∗ 100 = 12.79%
3487.75658

2) Calcule el volumen molar del líquido saturado y del vapor saturado con la
ecuación de Van der waals para cada una de las siguientes sustancias.

a) Propano a 40ºC donde Psat=13.71 bar

b) Propano a 50ºC donde Psat=17.16 bar.
c) Propano a 60ºC donde Psat=21.22 bar.
d) N-butano a 100ºC donde Psat=15.41 bar.
e) N-butano a 110ºC donde Psat=18.66 bar.
𝑅𝑇 𝑎
𝑃= − 2
𝑉−𝑏 𝑉
Donde:
𝑇𝑐2
𝑎 = 27 ∗ 𝑅 2 ∗
64𝑃𝑐
𝑅𝑇𝑐
𝑏=
8𝑃𝑐
Donde:
Pc=Presión critica
Tc=Temperatura critica
A continuación mostramos la curva de equilibrio liquido-vapor para el agua
en un gráfico T-V.

2
 Resolveremos la ecuación por el método de Newton-Raphson

La ecuación de Van der Waals se puede simplificar de la siguiente forma:

𝑎
(𝑃 + 𝑉 2 ) ∗ (𝑉 − 𝑏) = 𝑅𝑇

Multiplicando y simplificando:

𝑃𝑉 3 − (𝑃𝑏 + 𝑅𝑇)𝑉 2 + 𝑎𝑉 − 𝑎𝑏 = 0

𝑓(𝑉) = 𝑃𝑉 3 − (𝑃𝑏 + 𝑅𝑇)𝑉 2 + 𝑎𝑉 − 𝑎𝑏 = 0

La ecuación presenta 3 soluciones, la solución menor corresponde al volumen del

líquido saturado, y la solución mayor corresponde al volumen del vapor saturado.

Para aplicar el método de Newton-Raphson es necesario determinar la derivada

de f(V) , la cual se muestra a continuación:

′
𝑓(𝑉) = 3𝑃𝑉 2 − 2(𝑃𝑏 + 𝑅𝑇)𝑉 + 𝑎

Tc= 369.8 K Pc=42.48 bar

a= 9387568.87
b= 90.4695504

3
Para la fase vapor:
i Vi
0 1899.000073
1 1667.933747
2 1590.404768
3 1581.327893
4 1581.207313
5 1581.207292

Para la fase lìquida:

i Vi
0 0
1 90.46955038
2 134.9658936
3 151.3411031
4 154.121146
5 154.2041288

4
5
6
7
8
Metodología para calcular el valor de “Z”

3) Estime lo siguiente:
a) El volumen ocupado por 18 kg de etileno a 55ºC y 35 bar.
b) La masa de etileno contenida en un cilindro de 0.25 m3 a 50ºC y 15
bar.

Realice los cálculos usando la ecuación del gas ideal, la ecuación de

Van der Waals, y los gráficos generalizados para Z.

9
Calcule el Z y el V para el Hexafluoruro de Azufre a 75°C y 15 bar mediante las
ecuaciones siguientes.
𝑇° = 75°𝐶 = 348.15𝐾 ∧ 𝑃 = 15𝑏𝑎𝑟 ∧ 𝑅=
83.14𝑐𝑚3 𝑏𝑎𝑟 𝑚𝑜𝑙 −1 𝐾 −1
a) La ecuación virial truncada con los valores experimentales siguientes de los
coeficientes viriales:
𝐵 = −194𝑐𝑚3 𝑚𝑜𝑙 −1 𝐶 = 15300𝑐𝑚6 𝑚𝑜𝑙 −2

Determinación del V:
 Calculo de V por medio de la ecuación de gases ideales:
𝑅𝑇 (83.14𝑐𝑚3 𝑏𝑎𝑟 𝑚𝑜𝑙 −1 𝐾 −1 )(348.15𝐾)
𝑉= =
𝑃 15𝑏𝑎𝑟
= 1929.68𝑐𝑚3 𝑚𝑜𝑙 −1
 Calculo de V usando la ecuación virial truncada:
𝑅𝑇 𝐵 𝐶
𝑉= (1 + + 2 )
𝑃 𝑉 𝑉

Para obtener el valor V usaremos un proceso iterativo, para lo cual la

ecuación se escribe de la siguiente manera.
𝑅𝑇 𝐵 𝐶
𝑉𝑖+1 = (1 + + 2 )
𝑃 𝑉𝑖 𝑉𝑖
Donde i denota el número de iteración y 𝑉0 = 1929.68𝑐𝑚3 𝑚𝑜𝑙 −1
Para 𝑖 = 0
𝑅𝑇 𝐵 𝐶
𝑉1 = (1 + + 2 )
𝑃 𝑉0 𝑉0
83.14 × 348.15 −194 15300
𝑉1 = (1 + + )
15 1929.68 1929.682
𝑉1 = 1743.61 𝑐𝑚3 𝑚𝑜𝑙 −1

Se realizan 6 iteraciones y se muestran los resultados se muestran a

continuación:
Numero
de Volumen
iteración
0 1929.68
1 1743.61
2 1724.69
3 1722.55
4 1722.30
5 1722.27
6 1722.27

Por lo tanto 𝑉 = 1722.27 𝑐𝑚3 𝑚𝑜𝑙 −1

Determinación del Z.
10
 𝑃𝑉 15𝑏𝑎𝑟 × 1722.27 𝑐𝑚3 𝑚𝑜𝑙 −1
𝑍= =
𝑅𝑇 83.14𝑐𝑚3 𝑏𝑎𝑟 𝑚𝑜𝑙 −1 𝐾 −1 × 348.15𝐾
𝑍 = 0.8925
b) Ecuación de Van der Waals
𝑅𝑇 𝑎
𝑃= − 2
𝑉−𝑏 𝑉
𝑎 = 7.857 ∧ 𝑏 = 0.0879
Se despeja V para realizar un proceso iterativo de la siguiente manera:
𝑅𝑇
𝑉𝑖+1 = 𝑎 +𝑏
𝑃+ 2
𝑉𝑖

Después de 2 iteraciones se obtiene lo siguiente:

Numero
de Volumen
iteración
0 1929.68
1 1929.77
2 1929.77

Por lo tanto 𝑉 = 1929.77 𝑐𝑚3 𝑚𝑜𝑙 −1

Determinación del Z.
𝑃𝑉 15𝑏𝑎𝑟 × 1929.77 𝑐𝑚3 𝑚𝑜𝑙 −1
𝑍= =
𝑅𝑇 83.14𝑐𝑚3 𝑏𝑎𝑟 𝑚𝑜𝑙 −1 𝐾 −1 × 348.15𝐾
𝑍 = 1.00
c) Grafico Z
𝑃 15 𝑏𝑎𝑟
Z = ¿? Pr = 𝑃𝑐 → = 0.40 𝑏𝑎𝑟
37.59 𝑏𝑎𝑟

𝑇 348.15 𝐾
Tr = 𝑇𝑐 → 318.5 𝐾 = 1.09 𝐾
Z = 0.90 según gràfico.
PV = ZRT

𝑎𝑡𝑚.𝑙
14.85 atm × 𝑉 = 0.90 × 0.082 × 348.15 𝐾
𝑚𝑜𝑙.𝐾

𝑎𝑡𝑚. 𝑙
0.90 × 0.082 × 348.15𝐾
𝑉= 𝑚𝑜𝑙. 𝐾
14.85𝑎𝑡𝑚

𝑉 = 1.73 litros

11
Calcule Z y V para el etano a 50° C y 15 Bar con las ecuaciones siguientes:

a) La ecuación virial truncada [ecuación, con los valores experimentales siguientes

de los coeficientes viriales:

𝑩 = −156.7 𝑐𝑚3 𝑚𝑜𝑙 −1

𝑪 = 9650𝑐𝑚6 𝑚𝑜𝑙 −2
𝑻 = 50 + 273.15 = 323.15𝐾
𝑏𝑎𝑟. 𝑐𝑚3
𝑹 = 83.14
𝑚𝑜𝑙. 𝐾

𝑅. 𝑇 83.14𝑏𝑎𝑟 ∗ 𝑐𝑚3 ∗ 𝑚𝑜𝑙 −1 ∗ 𝐾 −1 ∗ 323.15𝐾 𝑐𝑚3

𝑽= = = 𝟏𝟕𝟗𝟏. 𝟏𝟏𝟐𝟕𝟑𝟑
𝑃 15𝑏𝑎𝑟 𝑚𝑜𝑙

𝒄𝒎𝟑
𝑽 = 1791.112733 → 𝑽𝒐𝒍𝒖𝒎𝒆𝒏 𝒅𝒆𝒍 𝑮𝒂𝒔 𝑰𝒅𝒆𝒂𝒍
𝒎𝒐𝒍

Ecuación Virial 2do coeficiente

𝑷∗𝑽 𝑩∗𝑷
𝒁= =𝟏+
𝑹∗𝑻 𝑹∗𝑻
𝑻
𝑺𝒊𝒎𝒑𝒍𝒊𝒇𝒊𝒄𝒂𝒏𝒅𝒐 ⇒ 𝑽 = 𝑹 ∗ +𝑩
𝑷

𝑏𝑎𝑟. 𝑐𝑚3 323.15𝐾 3 −1 )

𝒄𝒎𝟑
𝑽 = 83.14 ∗ + (−156.7𝑐𝑚 𝑚𝑜𝑙 = 𝟏𝟔𝟑𝟒. 𝟒𝟏𝟐𝟕𝟑𝟑
𝑚𝑜𝑙. 𝐾 15𝑏𝑎𝑟 𝒎𝒐𝒍

𝒄𝒎𝟑
𝑽 = 𝟏𝟔𝟑𝟒. 𝟒𝟏𝟐𝟕𝟑𝟑 → 𝑽𝒐𝒍𝒖𝒎𝒆𝒏 𝒅𝒆𝒍 𝑮𝒂𝒔 𝑰𝒅𝒆𝒂𝒍
𝒎𝒐𝒍

𝒄𝒎𝟑
𝑽𝒓𝒆𝒂𝒍 𝟏𝟔𝟑𝟒. 𝟒𝟏𝟐𝟕𝟑𝟑
𝒁= = 𝒎𝒐𝒍 = 𝟎. 𝟗𝟏𝟐𝟓𝟏𝟐
𝑽𝒊𝒅𝒆𝒂𝒍 𝒄𝒎𝟑
𝟏𝟕𝟗𝟏. 𝟏𝟏𝟐𝟕𝟑𝟑
𝒎𝒐𝒍
𝒁 = 𝟎. 𝟗𝟏𝟐𝟓𝟏𝟐

b) Ecuación de Van der Waals

Datos:
𝑻 = 323.15𝐾
𝑷 = 15𝑏𝑎𝑟
𝑻𝒄 = 305.3𝐾
𝑷𝒄 = 48.72𝑏𝑎𝑟
𝑏𝑎𝑟. 𝑐𝑚3
𝑹 = 83.14
𝑚𝑜𝑙. 𝐾
Usando Excel Determinamos los valores:

12
Vi+1=(R*T/P)*(1+B/Vi+C/Vi2) Ecuación del Virial
P= 15 bar
T= 50 ºC 323.15 K
cm3*bar*mol-
1
R= 83.14 *K-1
Tc= 305.3 K Pc= 48.72 bar
Ec. Newton Raphson a= 5578920.33

b= 65.1237736

Para la fase vapor:

i Vi
0 1791.112733
1 1661.726215 𝑓(𝑥𝑖 )
2 1638.920629 𝑥𝑖+1 = 𝑥𝑖 −
3 1638.231917 𝑓′(𝑥𝑖 )
4 1638.231298
5 1638.231298

Para la fase líquida:

i Vi
0 0
1 65.1237736
2 118.2889537
3 -61.75372168
4 27.21320185
5 83.79793431

13
c) Gráficos de factor “Z”

Determinamos Tr y Pr para hallar “Z”

𝑇 323.15𝐾
𝑇𝑟 ⟹ = = 1.058467082 = 1.05
𝑇𝑐 305.3𝐾

𝑝 15𝑏𝑎𝑟
𝑃𝑟 ⟹ = = 0.3078817734 = 0.31
𝑃𝑐 48.72𝑏𝑎𝑟

𝑇𝑟 = 1.06
𝑃𝑟 = 0.31

𝐷𝑒𝑙 𝑑𝑖𝑎𝑔𝑟𝑎𝑚𝑎 𝑑𝑒𝑡𝑒𝑟𝑚𝑖𝑛𝑎𝑚𝑜𝑠 "𝒁":

𝒁 = 𝟎. 𝟗𝟎

𝑏𝑎𝑟. 𝑐𝑚3
𝑍 ∗ 𝑅 ∗ 𝑇 0.90 ∗ 83.14 𝑚𝑜𝑙. 𝐾 ∗ 323.15𝐾 𝑐𝑚3
𝑉= = = 𝟏𝟔𝟏𝟐. 𝟎𝟎𝟏𝟒𝟔
𝑃 15𝑏𝑎𝑟 𝑚𝑜𝑙

𝑐𝑚3
𝑽 = 𝟏𝟔𝟏𝟐. 𝟎𝟎𝟏𝟒𝟔
𝑚𝑜𝑙

14