INTRODUCCION

Hasta el momento se han visto diferentes relaciones termodinámicas que permiten calcular el estado
termodinámico de alguna sustancia, para este ejercicio se centra en una sola especie química, el
agua, en donde se estimará la presión a partir de diferentes métodos:
• Grafico:
o Superficie P-v-T o graficas del factor del factor de compresibilidad.
• Tabular:
o Tablas de vapor.
• Analítico:
o Ecuaciones de estado.
Anteriormente ya se han manejado tablas de vapor y ecuación de Van Der Waals, estas relaciones
resultan muy adecuadas para efectuar operaciones matemáticas para el cálculo de propiedades
termodinámicas, en esta ocasión se presenta la relación p-v-T de Redlich-Kwong.
Redlich-Kwong es una ecuación de dos constantes, la cual es considerada como una de las mejores
ecuaciones de estado de dos constantes:
𝑅𝑇 𝑎
𝑃= −
𝑣̅ − 𝑏 𝑣̅ (𝑣̅ + 𝑏)𝑇 1/2
Análisis dimensional
1
𝑏𝑎𝑟 ∙ 𝑚3 𝑏𝑎𝑟 ∙ 𝑚6 𝐾 2
𝐾
𝑃 [=] 𝑘𝑚𝑜𝑙 ∙𝐾 − 𝑘𝑚𝑜𝑙 2 [ ]
𝑚3 𝑚3 𝑚3 𝑚3 𝑚3 1 = 𝑏𝑎𝑟
( ) 2
𝑘𝑚𝑜𝑙 − 𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙 + 𝑘𝑚𝑜𝑙 𝐾
Esta ecuación fue propuesta en 1949, la cual consta de una naturaleza empírica sin una justificación
rigurosa en términos moleculares, es explicita en la presión, pero no en el volumen o temperatura
específicos, y al igual que Van Der Waals es cubica en su expresión de volumen especifico, y
aunque es más difícil de manipular matemáticamente, es más precisa para las presiones altas.
Tiene su limite a medida que aumenta la presión o densidad. Para una mayor precisión se requiere
de ecuaciones con un mayor número de constantes y Redlich-Kwong también se ha ajustado a ello.
El cálculo de las constantes pueden ser obtenidas para determinadas sustancias ajustando la
ecuación para los datos de p-v-T en términos de presión y temperatura critica a partir del punto
crítico, como en Van Der Waals.
5
𝑅2 𝑇𝑐2
𝑎 = 0.42748
𝑃𝑐
𝑏𝑎𝑟 2 ∙ 𝑚6 2 1/2
𝑚𝑜𝑙 2 ∙ 𝐾2 ∗ 𝐾 𝐾 𝑏𝑎𝑟 ∙ 𝑚6 𝐾 1/2
𝑎[=] [=]
𝑏𝑎𝑟 𝑚𝑜𝑙 2
𝑅𝑇𝑐
𝑏 = 0.08664
𝑃𝑐
𝑚3
𝑏[=]
𝑚𝑜𝑙
 𝑏𝑎𝑟 ∙ 𝑚3
𝑅 = 0.0831447
𝑘𝑚𝑜𝑙 ∙ 𝐾
METODOLOGIA Y PROCEDIMIENTO
11.2 La presión de un depósito de 23,3 m3 no debe sobrepasar los 105 bar. Si se llena con 1000 kg
de vapor de agua a 360°C, estime la presión, en bar, utilizando:
(a) Tablas de vapor
(b) Carta de compresibilidad
(c) Redlich–Kwong.
(d) Van der Waals.
(e) ideal gas equation of state.
Volumen especifico (360°C, P = ¿?):
23,3m3 m3
𝑣̅𝑚𝑎𝑠𝑎 = = 0.0233
1000 𝑘𝑔 𝑘𝑔
m3 18𝑘𝑔 m3
𝑣̅𝑚𝑜𝑙𝑎𝑟 = 0.0233 ( ) = 0.4194
𝑘𝑔 𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙

Tablas de vapor
m3
A 360°C y un 𝑣̅ = 0.0233 𝑘𝑔

El volumen se encuentra en sobrecalentamiento

𝑣𝑓 < 𝑣𝑔 < 𝑣

m3 m3 m3
0.001895 < 0.006950 < 0.0233
𝑘𝑔 𝑘𝑔 𝑘𝑔
Por lo tanto, de la tabla A.4
𝑃 = 100 𝑏𝑎𝑟

Redlich-Kwong
𝑅𝑇 𝑎
𝑃= −
𝑣̅ − 𝑏 𝑣̅ (𝑣̅ + 𝑏)𝑇 1/2
Tabla A-24
𝑏𝑎𝑟 ∙ 𝑚6 𝐾 1/2
𝑎 = 142.59
𝑘𝑚𝑜𝑙 2
𝑚3
𝑏 = 0.02111
𝑘𝑚𝑜𝑙
 0.0831447 (633) 142.59
𝑃= −
0.4194 − 0.02111 0.4194(0.4194 + 0.02111)6331/2
𝑃 = 101.47 𝑏𝑎𝑟

Van Der Waals

𝑅𝑇 𝑎
𝑃= − 2
(𝑉̅ − 𝑏) 𝑉̅
Tabla A-24
𝑏𝑎𝑟 ∙ 𝑚6
𝑎 = 5.531
𝑘𝑚𝑜𝑙 2
𝑚3
𝑏 = 0.0305
𝑘𝑚𝑜𝑙
0.0831447(633) 5.531
𝑃= −
(0.4194 − 0.0305) (0.4194)2
𝑃 = 103.88 𝑏𝑎𝑟

Gas ideal
𝑅𝑇
𝑃=
𝑉̅
0.0831447(633)
𝑃=
(0.4194)
𝑃 = 125.5 𝑏𝑎𝑟

Carta de compresibilidad
De la tabla A-1
𝑇𝑐 = 647.3 𝐾
𝑃𝑐 = 220.9 𝑏𝑎𝑟
𝑇 633 𝐾
𝑇𝑟 = = = 0.978
𝑇𝑐 647.3 𝐾
 Gráficamente z=0.83
𝑃𝑉 𝑍𝑅𝑇
𝑍= ; 𝑃=
𝑅𝑇 𝑉
0.83(0.0831447)(633)
𝑃= = 102.9 𝑏𝑎𝑟
0.4194
0.83(0.0831447)(633)
𝑃= = 104.1 𝑏𝑎𝑟
0.4194

Apoyándonos de la ecuación generalizada de Redlich-Kwong

𝑣𝑅′ 0.42748
𝑍= ′ −
𝑣𝑅 − 0.08664 (𝑣𝑅 + 0.08664)𝑇 3/2
′
𝑅

m3
𝑣𝑃𝑐 0.4194 (220.9 𝑏𝑎𝑟)
𝑣𝑅′ = = 𝑘𝑚𝑜𝑙 = 1.72
𝑅𝑇𝑐 𝑏𝑎𝑟 ∙ 𝑚3
0.0831447 𝑘𝑚𝑜𝑙 ∙ 𝐾 (647.3 𝐾 )

1.72 0.42748
𝑍= − 3 = 0.82
1.72 − 0.08664 (
1.72 + 0.08664)0.9782
𝑃𝑉 𝑍𝑅𝑇
𝑍= ; 𝑃=
𝑅𝑇 𝑉
 0.82(0.0831447)(633)
𝑃= = 102.9 𝑏𝑎𝑟
0.4194

RESULTADOS

Tablas de Redlich- Van Der Carta de

Gas ideal
vapor Kwong Waals compresibilidad
PRESIÓN
100 101.47 103.88 125.5 104.1
(bar)

CONCLUSIONES
Los cálculos de las presiones quedaron muy cercanas, quedando debajo del limite de 105 bar, a
excepción del gas ideal, el cual al tratarse de una presión alta no puede predecir de manera
aproximada el comportamiento del gas, como referencia podemos tomar la fuente de tablas de vapor,
pues sus datos corresponden a fuentes experimentales, y se puede observar el por qué de la
importancia de la ecuación de estado de Redlich-Kwong, pues predice muy aproximadamente el valor
de la presión del sistema, con un error de 1.47%.
REFERENCIAS
Moran, M. J., & Shapiro, H. N. (2014). Appendices to accompany Fundamentals of Engineering
Thermodynamics, Eighth Edition pp (678,720) (8.a ed.). Wiley.
Moran, M. J., & Shapiro, H. N. (2004). Fundamentos de termodinámica técnica, pp (571,637) (2.a
ed.). Reverté.
Primero podemos empezar por calcular la fracción de n-butano:
𝑦𝑛 = 1 − 0.6037 = 0.3963
Seguimos calculando los factores ‘a’ y ‘b’ de la mezcla líquido vapor, de ambos
componentes:
Am → Factor A de metano
An → Factor A de n-butano
Bm → Factor B de metano
Bn → Factor B de n-butano

1
2
0.4278 1
𝐴𝑚 = ( 2.5 ) = 0.0113 𝑝𝑠𝑖𝑎−2
649.67
(673.1) (
343.9 )
1
2
0.4278 −
1
𝐴𝑛 = ( 2.5 ) = 0.0342 𝑝𝑠𝑖𝑎
2
649.67
(550.7) (
765.3 )

0.0867
𝐵𝑚 = = 1.51 × 10−4 𝑝𝑠𝑖𝑎−1
649.67
(673.1) (
765.3 )

0.0867
𝐵𝑛 = = 1.85 × 10−4 𝑝𝑠𝑖𝑎−1
649.67
(550.7) (
765.3 )

Aplicando las reglas de mezclas, obtenemos que:

 1
𝐴 = 0.6037(0.0113) + 0.3963(0.0342) = 0.02037 𝑝𝑠𝑖𝑎−2
𝐵 = 0.6037(1.51 × 10−4 ) + 0.3963(1.85 × 10−4 ) = 1.64 × 10−4
𝐵𝑃 = 1.64 × 10−4 (600) = 0.0984
Ahora podemos sustituir estos valores en la ecuación para obtener Z:
𝑍 3 − 𝑍 2 + 0.1408 𝑍 − 0.0244 = 0
Resolviendo esta ecuación para Z, obtenemos solo una raíz:
𝑍 = 0.8704
Ahora podemos calcular el volumen específico:
𝑍𝑅𝑇 (0.8704)(10.731)(649.67) 𝑓𝑡 3
𝑣= = = 10.11
𝑃 600 𝑙𝑏 𝑚𝑜𝑙
Problemario semana 3

Se calienta monóxido de carbono en un recipiente rígido, de

14.7 psia y 70 °F, hasta 800 °F. Calcule la presión final del
monóxido de carbono, considerándolo como a) un gas ideal, b)
usando la ecuación de Soave y c) usando la ecuación de
Peng-Robinson.

Monóxido Monóxido
de de
carbono carbono

V = cte. V = cte.
P = 14.7 psia P =?
T = 70°F T = 800°F

Lo primero es cambiar las unidades

1 𝑎𝑡𝑚
14.7 𝑝𝑠𝑖𝑎 ( ) = 1.0002 𝑎𝑡𝑚
14.6959 𝑝𝑠𝑖𝑎
5
70°𝐹 − 32 ( ) = 21.11°𝐶
9
5
800°𝐹 − 32 ( ) = 426.66°𝐶
9
a) Usando la ecuación del gas ideal
estado 1
𝑅𝑇
𝑉=
𝑃
𝑎𝑡𝑚 𝐿
(0.082 ) (294.26 𝐾) 𝐿
𝑚𝑜𝑙 𝐾
𝑉= = 24.12
1.0002 𝑎𝑡𝑚 𝑚𝑜𝑙
estado 2
𝑎𝑡𝑚 𝐿
𝑅𝑇 (0.082 𝑚𝑜𝑙 𝐾 ) (699.81 𝐾)
𝑃= = 𝐿
= 2.37 𝑎𝑡𝑚
𝑉 24.12
𝑚𝑜𝑙

b) Usando la ecuación del Soave

estado 1
usando la ecuación de Soave en su forma polinómica
𝑅𝑇 2 1 𝑎𝛼𝑏
𝑉3 − 𝑉 + (𝑎𝛼 + 𝑏𝑅𝑇 − 𝑃𝑏 2 )𝑉 − =0
𝑃 𝑃 𝑃
se calculan los parámetros a, b, ω y
𝑃𝑖 𝑠𝑎𝑡
𝜔 = [− log ( ) ]−1
𝑃𝑐 𝑇 =0.7
𝑟
para Tr
𝑇
𝑇𝑟 =
𝑇𝑐
𝑇
0.7 =
133 𝐾
𝑇 = 93.1 𝐾 = −180.05°𝐶
para la Psat se usa la ecuación de Antoine
𝐵
ln 𝑃𝑖 𝑠𝑎𝑡 = 𝐴 − ; 𝑇[=]°𝐶 𝑦 𝑃[=]𝑏𝑎𝑟
𝐶+𝑇
291.743
ln 𝑃𝑖 𝑠𝑎𝑡 = 3.81912 − = 1.65 𝑏𝑎𝑟
267.996 − 180.05°𝐶
se cambian las unidades a atmosferas
1 𝑎𝑡𝑚
1.65𝑏𝑎𝑟 = ( ) = 1.63 𝑎𝑡𝑚
1.01325 𝑏𝑎𝑟
sustituyendo en la ecuación del factor acéntrico
1.63 𝑎𝑡𝑚
𝜔 = [− log ( ) ] − 1 = 0.32
34.5423 𝑎𝑡𝑚 𝑇𝑟 =0.7
después
𝑎𝑡𝑚 𝐿 2
2 2
𝑅 𝑇𝑐 (0.082 ) (133 𝐾)2 𝑎𝑡𝑚 𝐿2
𝑎 = 0.42747 = 0.42747 𝑚𝑜𝑙 𝐾 = 1.47
𝑃𝑐 34.5423 𝑎𝑡𝑚 𝑚𝑜𝑙 2
𝑎𝑡𝑚 𝐿
𝑅𝑇𝑐 (0.082 ) (133 𝐾) 𝐿
𝑏 = 0.08664 = 0.08664 𝑚𝑜𝑙 𝐾 = 0.027
𝑃𝑐 34.5423 𝑎𝑡𝑚 𝑚𝑜𝑙
1 2
𝛼 = [1 + (0.48508 + 1.5517𝜔 − 0.1561𝜔2 ) (1 − 𝑇𝑟 2 )]
Donde la Tr se obtiene de
𝑇 (21.11 + 273.15) 𝐾
𝑇𝑟 = = = 2.21
𝑇𝑐 133 𝐾
sustituyendo
1 2
𝛼 = [1 + (0.48508 + 1.5517(0.32) − 0.1561(0.32)2 ) (1 − 2.212 )] = 1.46
ahora solo queda calcular V
 𝑎𝑡𝑚 𝐿
(0.082 ) (294.26 𝐾)
3
𝑉 − 𝑚𝑜𝑙 𝐾 𝑉2
1.0002 𝑎𝑡𝑚
1 𝑎𝑡𝑚 𝐿2
+ ((1.47 ) (1.46)
1.0002 𝑎𝑡𝑚 𝑚𝑜𝑙 2
𝐿 𝑎𝑡𝑚 𝐿
+ (0.027 ) (0.082 ) (294.26 𝐾)
𝑚𝑜𝑙 𝑚𝑜𝑙 𝐾
𝑎𝑡𝑚 𝐿2 𝐿
𝐿 2 (1.47 ) (1.46) (0.027 )
𝑚𝑜𝑙 2 𝑚𝑜𝑙
− (1.0002 𝑎𝑡𝑚) (0.027 ) )𝑉 −
𝑚𝑜𝑙 (1.0002 𝑎𝑡𝑚)
=0
𝐿 𝐿2 𝐿3
𝑉 3 − (24.12 ) 𝑉 2 + (2.708 ) 𝑉 − 0.057 =0
𝑚𝑜𝑙 𝑚𝑜𝑙 2 𝑚𝑜𝑙 3
𝐿
𝑉1 = 24.07
𝑚𝑜𝑙
𝐿
𝑉2 = 0.084
𝑚𝑜𝑙
𝐿
𝑉3 = 0.028
𝑚𝑜𝑙
usamos el volumen de gas (el más grande) porque a las condiciones del
problema el monóxido de carbono es gaseoso

estado 2
usando la forma estándar
𝑅𝑇 𝑎𝛼
𝑃= −
𝑣 − 𝑏 (𝑣 + 𝑏)𝑣
𝑎𝑡𝑚 𝐿 𝑎𝑡𝑚 𝐿2
(0.082 ) (699.81 𝐾) (1.47 ) (1.46)
𝑚𝑜𝑙 𝐾 𝑚𝑜𝑙 2
𝑃= −
𝐿 𝐿 𝐿 𝐿 𝐿
(24.07 − 0.027 ) (24.07 − 0.027 ) (24.07 )
𝑚𝑜𝑙 𝑚𝑜𝑙 𝑚𝑜𝑙 𝑚𝑜𝑙 𝑚𝑜𝑙
𝑃 = 2.38 𝑎𝑡𝑚

c) Usando la ecuación de Peng-Robinson

estado 1
usando la ecuación de Peng-Robinson en forma polinómica

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

se calculan los parámetros a, b y

𝑎𝑡𝑚 𝐿 2
2
𝑅 𝑇𝑐 2 (0.082 ) (133 𝐾)2 𝑎𝑡𝑚 𝐿2
𝑎 = 0.45724 = 0.45724 𝑚𝑜𝑙 𝐾 = 1.57
𝑃𝑐 34.5423 𝑎𝑡𝑚 𝑚𝑜𝑙 2
 𝑎𝑡𝑚 𝐿
𝑅𝑇𝑐 (0.082 ) (133 𝐾) 𝐿
𝑏 = 0.0778 = 0.0778 𝑚𝑜𝑙 𝐾 = 0.024
𝑃𝑐 34.5423 𝑎𝑡𝑚 𝑚𝑜𝑙
se usan las mismas Tr y ω del inciso pasado
1 2
2)
𝛼 = [1 + (0.37464 + 1.54226𝜔 − 0.26992𝜔 (1 − 𝑇𝑟 2 )]
1 2
𝛼 = [1 + (0.37464 + 1.54226(0.32) − 0.26992(0.32)2 ) (1 − 2.212 )]
𝛼 = 0.349

se sustituye en la ecuación cúbica

𝑅𝑇 1 𝑎𝛼𝑏 𝑅𝑇𝑏 2
𝑉 3 − ( − 𝑏) 𝑉 2 + (𝑎𝛼 − 2𝑏𝑅𝑇 − 3𝑃𝑏 2 )𝑉 − ( − − 𝑏3 ) = 0
𝑃 𝑃 𝑃 𝑃
𝑎𝑡𝑚 𝐿
(0.082 ) (294.26 𝐾) 𝐿
3
𝑉 −( 𝑚𝑜𝑙 𝐾 − 0.024 ) 𝑉2
1.0002 𝑎𝑡𝑚 𝑚𝑜𝑙

1 𝑎𝑡𝑚 𝐿2
+ ((1.57 ) (0.349)
1.0002 𝑎𝑡𝑚 𝑚𝑜𝑙 2
𝐿 𝑎𝑡𝑚 𝐿
− 2 (0.024 ) (0.082 ) (294.26 𝐾)
𝑚𝑜𝑙 𝑚𝑜𝑙 𝐾
𝐿 2
− 3(1.0002 𝑎𝑡𝑚) (0.024 ) )𝑉
𝑚𝑜𝑙
𝑎𝑡𝑚 𝐿2 𝐿
(1.57 ) (0.349) (0.024 )
𝑚𝑜𝑙 2 𝑚𝑜𝑙
−(
1.0002 𝑎𝑡𝑚

𝑎𝑡𝑚 𝐿 𝐿 2
(0.082 ) (294.26 𝐾) (0.024 ) 3
− 𝑚𝑜𝑙 𝐾 𝑚𝑜𝑙 − (0.024 𝐿 ) ) = 0
1.0002 𝑎𝑡𝑚 𝑚𝑜𝑙

𝐿 𝐿2 𝐿3
𝑉 3 − (24.1 ) 𝑉 2 − (0.61 ) 𝑉 + 7.618 =0
𝑚𝑜𝑙 𝑚𝑜𝑙 2 𝑚𝑜𝑙 3
resolviendo la cúbica
𝐿
𝑉1 = 24.11
𝑚𝑜𝑙
𝐿
𝑉2 = 0.556
𝑚𝑜𝑙
𝐿
𝑉3 = −0.55
𝑚𝑜𝑙
estado 2
usamos la forma estándar de la ecuación de Peng-Robinson y usamos el
volumen de gas (el más grande) porque a las condiciones del problema el
monóxido de carbono es un gas

𝑅𝑇 𝑎𝛼
𝑃= − 2
𝑉 − 𝑏 𝑉 + 2𝑏𝑉 − 𝑏 2
𝑎𝑡𝑚 𝐿
(0.082 ) (699.81 𝐾)
𝑃= 𝑚𝑜𝑙 𝐾
𝐿 𝐿
(24.11 − 0.024 )
𝑚𝑜𝑙 𝑚𝑜𝑙
𝑎𝑡𝑚 𝐿2
(1.57 ) (0.349)
𝑚𝑜𝑙 2
−
𝐿 2 𝐿 𝐿 𝐿 2
(24.11 ) + 2 (0.024 ) (24.11 ) − (0.024 )
𝑚𝑜𝑙 𝑚𝑜𝑙 𝑚𝑜𝑙 𝑚𝑜𝑙
𝑃 = 2.38 𝑎𝑡𝑚
1. Se desea construir un tanque de espera esférico cuya presión máxima de trabajo
será de 150 atm, para almacenar temporalmente 125 kg de CO2. La temperatura
máxima que puede alcanzar el CO2 es de 200°C. Estímese el diámetro interior
del tanque, utilizando: a) La ecuación del gas ideal .b) La gráfica generalizada
del factor de compresibilidad. d) La ecuación de Redlich-Kwong.
a) Gas ideal:
𝑃𝑉 = 𝑛𝑅𝑇
Para n:
𝑚
𝑛=
𝑃𝑀
Para el CO2, PM= 44 g/mol
125000𝑔
𝑛= 𝑔 = 2841 𝑚𝑜𝑙
44
𝑚𝑜𝑙
𝑎𝑡𝑚𝐿
𝑛𝑅𝑇 (2841𝑚𝑜𝑙) (0.08206 𝑚𝑜𝑙𝐾 ) (473𝐾)
𝑉= = = 735.144 𝐿
𝑃 150 𝑎𝑡𝑚
1𝑚3
735.144𝐿 ( ) = 0.735 𝑚3
1000𝐿

Para una esfera:

4
𝑉 = 𝜋𝑟 3 = 0.735 𝑚3
3
3(0.735𝑚3 )
3
𝑟 = = 0.175 𝑚3
4𝜋
3
𝑟 = √0.175 𝑚3 = 0.559 𝑚
𝐷𝑖á𝑚𝑒𝑡𝑟𝑜 = 𝐷 = 2𝑟 = 2(0.559𝑚) = 1.118 𝑚
b) Gráfica generalizada del factor de compresibilidad
Para el dióxido de carbono:
Pc=72.04 atm
Tc=32°C=305K
𝑃 150𝑎𝑡𝑚
𝑃𝑟 = = = 2.082
𝑃𝑐 72.04𝑎𝑡𝑚
𝑇 473 𝐾
𝑇𝑟 = = = 1.55
𝑇𝑐 305 𝐾
Tomando Z=0.85
𝑃𝑉
𝑍=
𝑛𝑅𝑇
𝑎𝑡𝑚𝐿
𝑍𝑛𝑅𝑇 (0.85)(2841 𝑚𝑜𝑙) (0.08206 𝑚𝑜𝑙𝐾 ) (473𝐾)
𝑉= = = 624.87 𝐿
𝑃 150 𝑎𝑡𝑚
1𝑚3
624.87𝐿 ( ) = 0.624 𝑚3
1000𝐿

Para una esfera:

4
𝑉 = 𝜋𝑟 3 = 0.624 𝑚3
3
3(0.624𝑚3 )
𝑟3 = = 0.149 𝑚3
4𝜋
3
𝑟 = √0.149 𝑚3 = 0.530 𝑚
𝐷𝑖á𝑚𝑒𝑡𝑟𝑜 = 𝐷 = 2𝑟 = 2(0.530𝑚) = 1.06 𝑚
c) La ecuación de R-K
La ecuación cúbica de R-K:
RT 2 1 a 2 ab
v3 − ( ) v + ( 0.5 − bRT − Pb )v − 0.5 = 0
P P T PT
Donde
0.08206𝑎𝑡𝑚𝐿 2
2
𝑅 𝑇𝑐 2.5 ( ) (305𝐾)2.5
𝑎 = 0.42748 = 0.42748 𝑚𝑜𝑙𝐾 = 64.916
𝑃𝑐 72.04 𝑎𝑡𝑚
𝑎𝑡𝑚𝐿
𝑅𝑇𝑐 (0.08206 ) (305𝐾)
𝑏 = 0.08664 = 0.08664 𝑚𝑜𝑙𝐾 = 0.0301
𝑃𝑐 72.04 𝑎𝑡𝑚
(0.08206)(473K) 2 1 64.916
v3 − ( )v + (
150 atm 150atm (473K)0.5
2
− 0.0301(0.08206)(473K) − (150atm)(0.0301) )v
(64.916)(0.0301)
− 0.5 =0
(150)(473K)
v3 − 0.258v2 + 0.0112v − 0.283 = 0
Resolviendo la cúbica con programa para resolver ec. Cúbicas:

𝐿
𝑉 = (0.3187 ) (2841 𝑚𝑜𝑙) = 905.426 𝐿
𝑚𝑜𝑙
1𝑚3
905.426𝐿 ( ) = 0.905 𝑚3
1000𝐿

4
𝑉 = 𝜋𝑟 3 = 0.905 𝑚3
3

3
3(0.905𝑚3 )
𝑟 = = 2.132 𝑚3
4𝜋
 3
𝑟 = √2.132 𝑚3 = 1.28 𝑚
𝐷𝑖á𝑚𝑒𝑡𝑟𝑜 = 𝐷 = 2𝑟 = 2(1.28𝑚) = 2.56 𝑚
 EJERCICIO 4.10

Utilícese la ecuación de estado R-K para predecir el volumen molar del líquido
de una mezcla equimolar de etano y n-decano a 100°𝐹 y 1000 𝑝𝑠𝑖𝑎. Compare
su resultado con el valor experimental de 2.13 𝑓𝑡 3 /𝑙𝑏𝑚𝑜𝑙 [J. Chem. Eng. Data.
7, 486 (1962)].

Solución
etano → E
n-decano → D
𝑃 = 1000 𝑝𝑠𝑖𝑎
𝑇 = 100°𝐹
𝑇 = (100 + 459.67)𝑅 = 559.67 𝑅
𝑣𝑙𝑖𝑞 = ?
𝑝𝑠𝑖𝑎∙𝑓𝑡 3
𝑅 = 10.73
𝑙𝑏𝑚𝑜𝑙∙𝑅

Al mencionar que se trata de una mezcla equimolar, se supondrá que ambos

componentes de la mezcla tienen la misma fracción molar en líquido, por lo tanto

𝑥𝐸 = 0.5
𝑥𝐷 = 0.5

Cuando la ecuación R-K se aplica a mezclas, se puede implementar la siguiente

forma equivalente del factor de compresibilidad

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

donde las reglas de mezclas son

𝐴 = ∑𝐶𝑖=1 𝐴𝑖 𝑦𝑖 o bien ∑𝐶𝑖=1 𝐴𝑖 𝑥𝑖 (2)

𝐵 = ∑𝐶𝑖=1 𝐵𝑖 𝑦𝑖 o bien ∑𝐶𝑖=1 𝐵𝑖 𝑥𝑖 (3)
siendo
1/2
𝑎𝑖 1/2 0.4278
𝐴𝑖 = ( 2 2.5 ) = ( 5/2
) (4)
𝑅 𝑇 𝑃𝑐 𝑇𝑟
𝑏𝑖 0.0867
𝐵𝑖 = = (5)
𝑅𝑇 𝑃𝑐 𝑇𝑟
Las propiedades críticas de los componentes mencionados son las siguientes
etano n-decano
𝑇𝑐,𝐸 = 550 𝑅 𝑇𝑐,𝐷 = 1114 𝑅
𝑃𝑐,𝐸 = 709.8 𝑝𝑠𝑖𝑎 𝑃𝑐,𝐷 = 306 𝑝𝑠𝑖𝑎

Calculando la temperatura reducida para cada componente, obtenemos que

559.67 𝑅
𝑇𝑟,𝐸 =
550 𝑅
𝑇𝑟,𝐸 = 1.0176

559.67 𝑅
𝑇𝑟,𝐷 =
1114 𝑅
𝑇𝑟,𝐷 = 0.5024
Calculando los 𝐴𝑖
1/2
0.4278
𝐴𝐸 = [ 5]
(709.8 𝑝𝑠𝑖𝑎)(1.0176)2
𝐴𝐸 = 24.02 × 10−3 𝑝𝑠𝑖𝑎−1/2
1/2
0.4278
𝐴𝐷 = [ 5]
(306 𝑝𝑠𝑖𝑎 )(0.5024)2
𝐴𝐷 = 88.39 × 10−3 𝑝𝑠𝑖𝑎−1/2
Calculando los 𝐵𝑖
0.0867
𝐵𝐸 =
(709.8 𝑝𝑠𝑖𝑎 )(1.0176)
𝐵𝐸 = 1.2 × 10−4 𝑝𝑠𝑖𝑎−1
0.0867
𝐵𝐷 =
(306 𝑝𝑠𝑖𝑎 )(0.5024)
𝐵𝐷 = 5.64 × 10−4 𝑝𝑠𝑖𝑎−1
Calculando 𝐴
𝐴 = 𝐴𝐸 𝑥𝐸 + 𝐴𝐷 𝑥𝐷
𝐴 = (24.02 × 10 𝑝𝑠𝑖𝑎−1/2 )0.5 + (88.39 × 10−3 𝑝𝑠𝑖𝑎−1/2 )0.5
−3

𝐴 = 56.205 × 10−3 𝑝𝑠𝑖𝑎−1/2

Calculando B
𝐵 = 𝐵𝐸 𝑥𝐸 + 𝐵𝐷 𝑥𝐷
𝐵 = (1.2 × 10 𝑝𝑠𝑖𝑎−1 )0.5 + (5.64 × 10−4 𝑝𝑠𝑖𝑎−1 )0.5
−4

𝐵 = 3.42 × 10−4 𝑝𝑠𝑖𝑎−1

Sustituyendo los valores obtenidos en la ecuación (1)
(56.205 × 10−3 )2 (56.205 × 10−3 )2
𝑍 3 − 𝑍 2 + (3.42 × 10−4 )(410.3) ( −4
− (3.42 × 10−4 )(410.3) − 1) 𝑍 − (3.42 × 10−4 )2 (410)2 = 0
3.42 × 10 3.42 × 10−4

Resolviendo en una calculadora CASIO fx-991ES, obtenemos

𝑍 = 0.4403
la cual es la única solución real.
De la ecuación generalizada
𝑃𝑣 = 𝑍𝑅𝑇 (6)

Despejando el volumen específico, 𝑣, obtenemos

𝑍𝑅𝑇
𝑣=
𝑃
Sustituyendo los valores
𝑝𝑠𝑖𝑎 ∙ 𝑓𝑡 3
(0.4403) (10.73 ) (559.67 𝑅)
𝑙𝑏𝑚𝑜𝑙 ∙ 𝑅
𝑣=
1000 𝑝𝑠𝑖𝑎
𝑣 = 2.64 𝑓𝑡 3 /𝑙𝑏𝑚𝑜𝑙

Calculando el porcentaje de error para un mejor análisis de resultados se tiene

que
|𝑣𝑒𝑥𝑝 − 𝑣𝑡𝑒𝑜 |
%𝐸 = × 100
𝑣𝑒𝑥𝑝
|2.13 − 2.64|
%𝐸 = × 100
2.13
%𝐸 ≅ 24%

Por lo tanto, con un error del 24%, el volumen molar del líquido de una mezcla
equimolar de etano y n-decano a 100°𝐹 y 1000 𝑝𝑠𝑖𝑎 es de 2.64 𝑓𝑡 3 /𝑙𝑏𝑚𝑜𝑙.
INTRODUCCION
La ecuación de Redlich-Kwong es una ecuación de estado utilizada en la termodinámica
para describir el comportamiento de los gases. Fue propuesta en 1949 por los químicos
Joseph Redlich y Myron Kwong, y es una alternativa a la ecuación de estado de Van der
Waals.

La ecuación de Redlich-Kwong tiene en cuenta las interacciones entre las moléculas del
gas, lo que la hace más precisa en comparación con la ecuación de Van der Waals, que no
considera estas interacciones. La ecuación de Redlich-Kwong se utiliza comúnmente para
calcular las propiedades termodinámicas de los gases, como la presión, el volumen y la
temperatura.

La ecuación de Redlich-Kwong tiene la forma de una ecuación cúbica de estado, y se

puede utilizar para describir gases reales en condiciones de alta presión y baja
temperatura. Aunque la ecuación de Redlich-Kwong es una mejora en la precisión en
comparación con la ecuación de Van der Waals, todavía tiene limitaciones y se han
desarrollado ecuaciones de estado más precisas para describir el comportamiento de los
gases.
OBJETIVO
Realizar un ejercicio relacionado con el tema visto en clase el cual en esta ocasión fue la
ecuación de Redlich - Kwong

EJERCICIO
Un recipiente de 30 000 cm3 contiene una masa de 500 g de amoniaco gaseoso.
El recipiente está sumergido en un baño que tiene una temperatura constante
de 65’C. Calcule la presión del gas mediante la ecuación de Redlich – Kwong.

𝑣 = 30000𝑐𝑚3 = 30𝐿
𝑇 = 65°𝐶 = 338.15 𝐾
𝑚 = 500𝑔
𝑇𝑐 = 405.5 𝐾
𝑃𝑐 = 11.28 𝑀𝑝𝑎 = 111.32 𝑎𝑡𝑚
𝑎𝑡𝑚 ∙ 𝐿
𝑅 = 0.08205
𝑚𝑜𝑙 ∙ 𝐾
𝑔
𝑃𝑀 = 17.03
𝑚𝑜𝑙
 500𝑔
𝑛= 𝑔 = 29.36 𝑚𝑜𝑙
17.03
𝑚𝑜𝑙
5 𝑎𝑡𝑚 ∙ 𝐿 2 5 1
0.42748𝑅 2 𝑇𝑐2 (0.42748) (0.08205 ) (405.5𝐾)2 𝑎𝑡𝑚 ∙ 𝐿2 ∙ 𝐾 2
𝑎= = 𝑚𝑜𝑙 ∙ 𝐾 = 85.6
𝑃𝑐 (111.32 𝑎𝑡𝑚) 𝑚𝑜𝑙 2

𝑎𝑡𝑚 ∙ 𝐿
0.08664𝑅𝑇𝑐 (0.08664) (0.08205 𝑚𝑜𝑙 ∙ 𝐾 ) (405.5 𝐾) 𝐿
𝑏= = = 0.128
𝑃𝑐 (111.32𝑎𝑡𝑚) 𝑚𝑜𝑙

𝑛𝑅𝑇 𝑛2 𝑎
𝑃= −
𝑣 − 𝑛𝑏 [(√𝑇)(𝑣 − 𝑏𝑛)][𝑣]
1
𝑎𝑡𝑚 ∙ 𝐿2 ∙ 𝐾 2
𝑎𝑡𝑚 ∙ 𝐿 (29.36 𝑚𝑜𝑙)2 (85.6 )
(29.36 𝑚𝑜𝑙) (0.08205 ) (338.15 𝐾) 𝑚𝑜𝑙2
𝑃= 𝑚𝑜𝑙 ∙ 𝐾 − = 25.94 𝑎𝑡𝑚
𝐿 𝐿
(30𝐿 − (29.36 𝑚𝑜𝑙) (0.128 𝑚𝑜𝑙)) [(√338.15𝐾) (30𝐿 − (0.128 𝑚𝑜𝑙) (29.36 𝑚𝑜𝑙))] [30𝐿]

CONCLUSION

Se puede concluir que la ecuación de Redlich-Kwong proporciona una descripción

razonablemente precisa del comportamiento de gases a altas presiones y bajas
temperaturas, pero no es adecuada para describir el comportamiento de gases a altas
temperaturas y bajas presiones.
Además, se puede utilizar la ecuación de Redlich-Kwong para predecir el comportamiento
de mezclas de gases, lo que es importante en la industria química. La ecuación de Redlich-
Kwong también es útil para la modelización de procesos termodinámicos, como la
compresión y expansión de gases, y para la simulación de procesos industriales, como la
producción de gases a alta presión.
Introducción
Las ecuaciones de estado son muy importantes en los estudios termodinámicos de
diferentes sistemas. Estas ecuaciones permiten estimar propiedades
termodinámicas de compuestos y sustancias, así como, evaluar los estados de
sustancias puras y de mezclas, ya que son modelos matemáticos que no requieren
el uso de diagramas termodinámicos, eso sí, requieren un cálculo numérico preciso.

Problema

Primero calculamos las fracciones molares de la mezcla:

0.2733
𝑋𝑁 = 14.0067 = 0.4469
0.2733 0.726
+
14.0067 30.07
𝑋𝐶2𝐻6 = 1 − 0.4469 = 0.5530

Usando Redlich-Kwong calculamos los parámetros A y B:

1
0.4278 2
𝐴=( )
𝑃𝑐 𝑇𝑟2.5
0.0867
𝐵=( )
𝑃𝑐 𝑇𝑟
Calculamos los parámetros tanto del etano como del nitrógeno:
1
2
0.4278
𝐴𝑁 = ( ) = 5.4019𝑥10−3 𝑝𝑠𝑖𝑎−1
859.67 °𝑅 2.5
492.9 𝑝𝑠𝑖𝑎 ∗ ( 221.3 °𝑅 )
1
2
0.4278
𝐴𝐶2𝐻6 = ( ) = 0.0140 𝑝𝑠𝑖𝑎−1
859.67 °𝑅 2.5
709.8 𝑝𝑠𝑖𝑎 ∗ ( )
550°𝑅

0.0867
𝐵𝑁 = ( ) = 4.5280𝑥10−5 𝑝𝑠𝑖𝑎−1
859.67 °𝑅
492.9 𝑝𝑠𝑖𝑎 ∗ ( 221.3 °𝑅 )
 0.0867
𝐵𝐶2𝐻6 = ( ) = 7.8147𝑥10−5 𝑝𝑠𝑖𝑎−1
859.67 °𝑅
709.8 𝑝𝑠𝑖𝑎 ∗ ( )
550 °𝑅
Después usamos la sumatoria de parámetros usando la regla de las mezclas:
𝑐

𝐴 = ∑ 𝐴𝑖 𝑦1 = 0.4469 ∗ 5.4019𝑥10−3 𝑝𝑠𝑖𝑎−1 + 0.5530 ∗ 0.0140 𝑝𝑠𝑖𝑎−1

𝑖=1
= 0.0101𝑝𝑠𝑖𝑎−1
𝑐

𝐵 = ∑ 𝐵𝑖 𝑦1 = 0.4469 ∗ 4.5280𝑥10−5 𝑝𝑠𝑖𝑎−1 + 0.5530 ∗ 7.8147𝑥10−5 𝑝𝑠𝑖𝑎−1

𝑖=1

= 6.3450𝑥10−5 𝑝𝑠𝑖𝑎−1
𝐵𝑃 = 6.3450𝑥10−5 𝑝𝑠𝑖𝑎−1 ∗ 2000𝑝𝑠𝑖𝑎 = 0,1269
Ahora usando la fórmula del factor de compresibilidad Z de R-K donde:
𝐴2 𝐴2
𝑍 3 − 𝑍 2 + 𝐵𝑃 ( − 𝐵𝑃 − 1) 𝑍 − (𝐵𝑃)2 = 0
𝐵 𝐵

𝑍 3 − 𝑍 2 + 0.0610𝑍 + 0.0258 = 0
Resolviendo la ecuación cúbica nos da:
𝑍1 = 0.9004
𝑍2 = 0.2262
𝑍3 = −0.1266

Usando formula del gas ideal con el factor de compresibilidad tenemos;

𝑝𝑠𝑖𝑎 ∗ 𝑓𝑡 3
𝑍𝑅𝑇 0.9004 ∗ 10.731 ∗ 859.67°𝑅 𝑓𝑡 3
𝑣= = 𝑙𝑏𝑚𝑜𝑙 ∗ °𝑅 = 4.15
𝑃 2000 𝑝𝑠𝑖𝑎 𝑙𝑏𝑚𝑜𝑙
La masa molar es dada como:
𝑐
𝑙𝑏 𝑙𝑏 𝑙𝑏
∑ 𝑀𝑖 𝑦𝑖 = 0.4469 ∗ 14.0067 + 0.5530 ∗ 30.07 = 22.88
𝑙𝑏𝑚𝑜𝑙 𝑙𝑏𝑚𝑜𝑙 𝑙𝑏𝑚𝑜𝑙
𝑖=1

Sacando el volumen especifico:

 𝑓𝑡 3
𝑣 4.15 𝑓𝑡 3
= 𝑙𝑏𝑚𝑜𝑙 = 0.1815
𝑀 22.88 𝑙𝑏 𝑙𝑏
𝑙𝑏𝑚𝑜𝑙
OBJETIVO

Utilizar y aplicar diferentes ecuaciones de estado en un fluido térmico para conocer

sus propiedades en las diferentes condiciones que esté presente.

INTRODUCCIÓN

La ley de gas ideal trata a las moléculas de un gas, como partículas puntuales con
colisiones perfectamente elásticas. Esto funciona bien en muchas circunstancias
experimentales, con gases diluidos. Pero las moléculas de gas no son masas
puntuales, y hay circunstancias donde las propiedades de las moléculas tienen un
efecto medible experimentalmente. Existen varias ecuaciones, llamadas ecuaciones
de estado que incluyen estas variaciones de las moléculas dentro de sus variantes
para poder determinar con mayor precisión las propiedades termodinámicas de la
sustancia que se trabaje.

Redlich—Kwong

La ecuación de estado de Redlich—Kwong es una ecuación empírica algebraica que

relaciona la temperatura, la presión y el volumen de gases. Generalmente es más
precisa que la van der Waals y las ecuaciones de estado de gas ideales a temperaturas
por encima de la temperatura crítica. Fue formulado por Otto Redlich y Joseph Neng
Shun Kwong en 1949, quienes demostraron que una simple ecuación de estado de
dos parámetros bien podría reflejar la realidad en muchas situaciones. Redlich y
Kwong revisaron la Ecuación de Estado de van der Waals y propusieron las siguientes
expresiones:

Donde…
Donde a es la interacción entre las moléculas y b el volumen ocupado por estas en
esta ecuación y en las demás que se describirán a continuación.

La ecuación cúbica de Redlich—Kwong es:

𝑅𝑇 1 𝑎 𝑎𝑏
𝑣3 − 𝑣2 ( )+ ( − 𝑃𝑏 2 − 𝑅𝑇𝑏) 𝑣 − =0
𝑃 𝑃 √𝑇 𝑃√𝑇

Soave

En 1972, Soave propuso una modificación importante a la ecuación de Redlich –

Kwong. Todas las modificaciones a la ecuación de Van Der Waals se habían centrado
en la dependencia de la temperatura del parámetro atractivo. Soave amplió esto al
proponer una dependencia de dos variables para “a”. Fue la primera vez que “a” se
expresó no sólo en función de la temperatura, sino también en función de la forma
esfericidad de las moléculas.

Como recordamos, el factor acéntrico de Pitzer es una medida de la configuración y

esfericidad de la molécula. También puede verse como una medida de la deformidad
de la molécula. Dicha ecuación esta dada por:

Donde…

2
𝛼 = [1 + (0.48508 + 1.55171𝜔 − 0.15613𝜔2 )(1 − √𝑇𝑟)]

La ecuación cúbica de Soave es:

 𝑅𝑇 1 𝑎𝛼𝑏
𝑣3 − 𝑣2 ( ) + (𝑎𝛼 − 𝑏𝑅𝑇 − 𝑃𝑏 2 )𝑣 − =0
𝑃 𝑃 𝑃

Y en términos de la compresibilidad queda como:

Donde…

Peng- Robinson

Esta ecuación es un modelo semiempírico que modifica el término de presión de

atracción de la ecuación de van der Waals con el fin de determinar con mayor precisión
el valor de la densidad de la fase líquida y el valor de la presión de vapor de sustancias
puras y mezclas multicomponentes (Peng & Robinson, 1975).

El modelo, al igual que el planteado por Soave, propone que el término de

interacciones intermoleculares es función de las propiedades críticas, del factor
acéntrico y de la temperatura. El factor de ajuste del modelo inicialmente se estimó
como el valor numérico necesario para que la ecuación propuesta pudiera calcular la
presión de vapor y las condiciones de equilibrio de las sustancias reportadas en la
literatura. La ecuación queda expresada como:

𝑅𝑇 𝑎𝛼
𝑃= − 2
𝑉 − 𝑏 𝑣 + 2𝑏𝑣−𝑏 2

Donde…

𝑅2 𝑇𝑐 2
𝑎 = 0.45724
𝑃𝑐
 𝑅𝑇𝑐
𝑏 = 0.07780
𝑃𝑐

2
𝛼 = [1 + (0.37464 + 1.54226𝜔 − 0.26992𝜔2 )(1 − √𝑇𝑟)]

La ecuación cúbica de Peng-Robinsón es:

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

DESARROLLO

1) Calcular el volumen para el líquido a 230°C y 2.795 MPa con la ecuación de

Redlich Kwong, Soave y Peng Robinsón y comparar con los datos obtenidos de
tablas. Vf = 0.001209 m3/Kg, Vg = 0.07158 m3/Kg.

Para resolver el presente ejercicio se necesitan algunos datos, los cuales se

obtuvieron de diferentes fuentes: Tc = 647.13K, Pc = 216.3 atm, w = 0.344,
PM=18.015 g-mol

Ahora, ya que conocemos estos datos es conveniente convertir nuestras unidades

de temperatura y presión:

T = 230 + 273.15 = 503.15K

P = 2.795Mpa = 27.58 atm

Tr = T/Tc = 0.7775

a) Resolución con Redlich Kwonn

𝑅𝑇 1 𝑎 𝑎𝑏
𝑣3 − 𝑣2 ( ) + ( − 𝑃𝑏 2 − 𝑅𝑇𝑏) 𝑣 − =0
𝑃 𝑃 √𝑇 𝑃√𝑇
 Lo primero es obtener los parámetros a y b:
𝑎𝑡𝑚𝐿 2
(0.08206 ) (647.13𝐾)2.5 𝑎𝑡𝑚 𝐿2 𝐾 1/2
𝑎 = 0.4278 𝑚𝑜𝑙𝐾 = 141.64
(216.3 𝑎𝑡𝑚) 𝑛2
𝑎𝑡𝑚𝐿
(0.08206 ) (647.3𝐾) 𝐿
𝑏 = 0.086640 𝑚𝑜𝑙𝐾 = 0.02125
8(216.3 𝑎𝑡𝑚) 𝑛

Después sustituimos estos valores en la ecuación cúbica:

𝑎𝑡𝑚𝐿
(0.08206 ) (503.15𝐾 )
𝑣3 − 𝑣2 [ 𝑚𝑜𝑙𝐾 ]
27.58𝑎𝑡𝑚

1
𝑎𝑡𝑚 𝐿2 𝐾 2
(141.64 )
1 𝑛2 𝐿 2 𝑎𝑡𝑚𝐿 𝐿
+ − (27.58𝑎𝑡𝑚) (0.0212 ) − (0.08206 ) (503.15𝐾) (0.0212 ) 𝑣
27.58𝑎𝑡𝑚 √(503.15𝐾) 𝑛 𝑚𝑜𝑙𝐾 𝑛

[ ]

1
𝑎𝑡𝑚 𝐿2 𝐾 2 𝐿
(141.64 2 ) (0.0212 )
𝑛 𝑛
− =0
(27.58𝑎𝑡𝑚)√(503.15𝐾 )
[ ]

𝐿 2 𝐿2 𝐿3
𝑣 3 − 1.49 𝑣 + 0.197 2 𝑣 − 4.85𝑥10−3 3 = 0
𝑛 𝑛 𝑛

Usando un software para ecuaciones cúbicas obtenemos 3 diferentes

volúmenes:

𝐿 1𝑚3 1𝑛 1000𝑔 𝑚3
𝑣1 = 0.03237 ( )( )( ) = 0.001798
𝑛 1000𝐿 18.015𝑔𝑚𝑜𝑙 1𝐾𝑔 𝐾𝑔

𝐿 𝑚3
𝑣2 = 1.346 = 0.074
𝑛 𝐾𝑔

𝐿 𝑚3
𝑣3 = 0.02227 = 0.00618
𝑛 𝐾𝑔
Los valores que nos interesan son el volumen mayor y menor que corresponden al
volumen del vapor y líquido respectivamente.

Ahora, con estos datos obtenidos comparamos el porcentaje de error con los
obtenidos experimentalmente:

|𝐸𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙 − 𝑇𝑒ó𝑟𝑖𝑐𝑜|
%𝐸𝑟𝑟𝑜𝑟 = ∗ 100
𝐸𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙

𝑚3 𝑚3
|0.07158 − 0.074 |
𝐾𝑔 𝐾𝑔
%𝐸𝑟𝑟𝑜𝑟 𝑣𝑎𝑝𝑜𝑟 = ∗ 100 = 3.38%
𝑚3
0.07158
𝐾𝑔

𝑚3 𝑚3
|0.0001209 − 0.001798 |
𝐾𝑔 𝐾𝑔
%𝐸𝑟𝑟𝑜𝑟 𝑙𝑖𝑞𝑢𝑖𝑑𝑜 = ∗ 100 = 48.71%
𝑚3
0.001209
𝐾𝑔

b) Resolución con Soave

𝑅𝑇 1 𝑎𝛼𝑏
𝑣3 − 𝑣2 ( ) + (𝑎𝛼 − 𝑏𝑅𝑇 − 𝑃𝑏 2 )𝑣 − =0
𝑃 𝑃 𝑃

2
𝛼 = [1 + (0.48508 + 1.55171𝜔 − 0.15613𝜔2 )(1 − √𝑇𝑟)]

Lo primero es obtener los parámetros a, b y alfa:

𝑎𝑡𝑚𝐿 2
(0.08206 ) (647.13𝐾)2 𝑎𝑡𝑚 𝐿2
𝑎 = 0.4278 𝑚𝑜𝑙𝐾 = 5.53
(216.3 𝑎𝑡𝑚) 𝑛2
 𝑎𝑡𝑚𝐿
(0.08206 ) (647.3𝐾) 𝐿
𝑏 = 0.086640 𝑚𝑜𝑙𝐾 = 0. .02125
(216.3 𝑎𝑡𝑚) 𝑛
2
𝛼 = [1 + (0.48508 + 1.55171(0.344) − 0.15613(0.344)2 )(1 − √0.7775)] = 1.251

Después sustituimos estos valores en la ecuación cúbica:

𝑎𝑡𝑚𝐿
(0.08206 ) (503.15𝐾 )
𝑣3 − 𝑣2 [ 𝑚𝑜𝑙𝐾 ]
27.58𝑎𝑡𝑚

2
1 𝑎𝑡𝑚 𝐿 𝑎𝑡𝑚𝐿 𝐿 𝐿 2
+ [(5.53 )(1.251 ) − (0.08206 ) (503.15𝐾) (0.0212 ) − (27.58𝑎𝑡𝑚) (0.0212 ) ]𝑣
27.58𝑎𝑡𝑚 𝑛2 𝑚𝑜𝑙𝐾 𝑛 𝑛

𝑎𝑡𝑚 𝐿2 𝐿
(5.53 2 )(1.251) (0.02125 )
𝑛 𝑛
− =0
(27.58𝑎𝑡𝑚)
[ ]

3
𝐿 2 𝐿2 −3
𝐿3
𝑣 − 1.497 𝑣 + 0.1978 2 𝑣 − 5.53𝑥10 =0
𝑛 𝑛 𝑛3

Usando un software para ecuaciones cúbicas obtenemos 3 diferentes

volúmenes:

𝐿 1𝑚3 1𝑛 1000𝑔 𝑚3
𝑣1 = 0.0371 ( )( )( ) = 0.00205
𝑛 1000𝐿 18.015𝑔𝑚𝑜𝑙 1𝐾𝑔 𝐾𝑔

𝐿 𝑚3
𝑣2 = 1.354 = 0.0752
𝑛 𝐾𝑔

𝐿 𝑚3
𝑣3 = 0.02227 = 0.0056
𝑛 𝐾𝑔

Los valores que nos interesan son el volumen mayor y menor que corresponden al
volumen del vapor y líquido respectivamente.

Ahora, con estos datos obtenidos comparamos el porcentaje de error con los
obtenidos experimentalmente.
 |𝐸𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙 − 𝑇𝑒ó𝑟𝑖𝑐𝑜|
%𝐸𝑟𝑟𝑜𝑟 = ∗ 100
𝐸𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙

𝑚3 𝑚3
|0.07158 − 0.0754 |
𝐾𝑔 𝐾𝑔
%𝐸𝑟𝑟𝑜𝑟 𝑣𝑎𝑝𝑜𝑟 = ∗ 100 = 5.05%
𝑚3
0.07158
𝐾𝑔

𝑚3 𝑚3
|0.0001209 − 0.00205 |
𝐾𝑔 𝐾𝑔
%𝐸𝑟𝑟𝑜𝑟 𝑣𝑎𝑝𝑜𝑟 = ∗ 100 = 69.56%
𝑚3
0.001209
𝐾𝑔

c) Resolución con Peng Robinsón

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

𝑅2 𝑇𝑐 2
𝑎 = 0.45724
𝑃𝑐

𝑅𝑇𝑐
𝑏 = 0.07780
𝑃𝑐

2
𝛼 = [1 + (0.37464 + 1.54226𝜔 − 0.26992𝜔2 )(1 − √𝑇𝑟)]

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

Lo primero es obtener los parámetros a, b y alfa:

𝑎𝑡𝑚𝐿 2
(0.08206 ) (647.13𝐾)2 𝑎𝑡𝑚 𝐿2
𝑎 = 0.45724 𝑚𝑜𝑙𝐾 = 3.6
(216.3 𝑎𝑡𝑚) 𝑛2

𝑎𝑡𝑚𝐿
(0.08206 ) (647.3𝐾) 𝐿
𝑏 = 0.07780 𝑚𝑜𝑙𝐾 = 0. .0148
(216.3 𝑎𝑡𝑚) 𝑛
 2
𝛼 = [1 + (0.37464 + 1.54226(0.344) − 0.26992(0.344)2 )(1 − √0.7775)]
= 1.217

Después sustituimos estos valores en la ecuación cúbica:

𝑎𝑡𝑚𝐿
(0.08206 ) (503.15𝐾 ) 𝐿
𝑣3 − 𝑣2 [ 𝑚𝑜𝑙𝐾 − 0. .0148 ]
27.58𝑎𝑡𝑚 𝑛

(1.217) (3.6
𝑎𝑡𝑚 𝐿2 ) 𝐿 𝑎𝑡𝑚𝐿
𝑛2 (0. .0148 ) (0.08206 ) (503.15𝐾) 𝐿 2
𝑛 𝑚𝑜𝑙𝐾
+ −2 − 3 (0.0148 ) 𝑣
27.58 𝑎𝑡𝑚 27.58𝑎𝑡𝑚 𝑛
[ ]

(1.217) (3.6
𝑎𝑡𝑚 𝐿2 𝐿 𝐿 2
𝑎𝑡𝑚𝐿
) (0. .0148 ) (0.0148 ) (0.08206 ) (503.15𝐾) 3
𝑛2 𝑛 𝑛 𝑚𝑜𝑙𝐾 𝐿
− + + (0.0148 ) = 0
27.58 𝑎𝑡𝑚 27.58𝑎𝑡𝑚 𝑛
[ ]

𝐿 2 𝐿2 𝐿3
𝑣 3 − 1.482 𝑣 + 0.114 2 𝑣 − 2.0176𝑥10−3 3 = 0
𝑛 𝑛 𝑛

Usando un software para ecuaciones cúbicas obtenemos 3 diferentes

volúmenes:

𝐿 1𝑚3 1𝑛 1000𝑔 𝑚3
𝑣1 = 0.027 ( )( )( ) = 0.001499
𝑛 1000𝐿 18.015𝑔𝑚𝑜𝑙 1𝐾𝑔 𝐾𝑔

𝐿 𝑚3
𝑣2 = 1.4017 = 0.0779
𝑛 𝐾𝑔

𝐿 𝑚3
𝑣3 = 0.0533 = 0.002
𝑛 𝐾𝑔

Los valores que nos interesan son el volumen mayor y menor que corresponden al
volumen del vapor y líquido respectivamente.

Ahora, con estos datos obtenidos comparamos el porcentaje de error con los
obtenidos experimentalmente.
 |𝐸𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙 − 𝑇𝑒ó𝑟𝑖𝑐𝑜|
%𝐸𝑟𝑟𝑜𝑟 = ∗ 100
𝐸𝑥𝑝𝑒𝑟𝑖𝑚𝑒𝑛𝑡𝑎𝑙

𝑚3 𝑚3
|0.07158 − 0.0779 |
𝐾𝑔 𝐾𝑔
%𝐸𝑟𝑟𝑜𝑟 𝑣𝑎𝑝𝑜𝑟 = ∗ 100 = 8.82%
𝑚3
0.07158
𝐾𝑔

𝑚3 𝑚3
|0.0001209 − 0.001499 |
𝐾𝑔 𝐾𝑔
%𝐸𝑟𝑟𝑜𝑟 𝑣𝑎𝑝𝑜𝑟 = ∗ 100 = 23.98%
𝑚3
0.001209
𝐾𝑔

RESULTADOS

Los resultados numéricos aparecen subrayados en amarillo.

BIBLIOGRAFÍA

Reviews, C. T. & Cram101 Incorporated. (2013b). Studyguide for Thermodynamics: an

Engineering Approach by Yunus Cengel, ISBN 9780077366742. Cram101

Incorporated.

Cheng, C. H., & Abdel-Hady, S. (2018). Advanced Thermodynamics Engineering. Scitus

Academics LLC.

Libretexts. (2022, 2 noviembre). 16.2: Ecuaciones de Estado de van der Waals y Redlich-

Kwong. LibreTexts Español.

[Link]

3%B3rica/Qu%C3%ADmica_F%C3%ADsica_(LibreTexts)/16:_Las_propiedades_

de_los_gases/16.02:_Ecuaciones_de_Estado_de_van_der_Waals_y_Redlich-

Kwong
PROBLEMARIO
Un recipiente totalmente cerrado contiene butano en fase líquida y gaseosa a
𝟏𝟓𝟎°𝑭. Calcular la densidad de cada fase utilizando la ecuación cúbica de
estado de Soave-Redlich-Kwong.
Solución:
Es necesario conocer los valores críticos de presión y temperatura, así como el
factor acéntrico para el propano, estos de la Tabla 3.1 del capítulo 3, de donde
se obtienen los siguientes valores:

𝑃𝑐 = 616, 𝑇𝑐 = 206.06°𝐹 y 𝜔 = 0.1522

Calculando la temperatura reducida con la ecuación 3.7:

𝑇 610°𝑅
𝑇𝑟 = = = 0.9158
𝑇𝑐 666.06°𝑅

Se determina 𝑚, con la ecuación 6.30:

𝑚 = 0.480 + 1.574𝜔 − 1.76

= 0.480 + 1.574(0.1522) − 1.76 (0.1522)2
= 0.6787

se calcula el valor de 𝛼 con la ecuación 6.29:

2 2
𝛼 = [1 + 𝑚(1 − 𝑇𝑟 1/2 )] = [1 + 0.6787(1 − (0.9158)1/2 )] = 1.0592
Posteriormente se calculan los valores de las constantes 𝑎𝑐 y 𝑏 con las
ecuaciones 6.32 y 6.33: