Ecuaciones de Variación de

Propiedades

GASES REALES
 GASES REALES (Usando ecuaciones de estado)
La variación de propiedades para cualquier gas están dadas por las
ecuaciones: 1, 2, 3 y 4.
 

Para un gas ideal, muchos términos de estas ecuaciones se anulan.

Sin embargo, estas simplificaciones no son validas para gases reales.
1.- Los coeficientes de dP y dv, raras veces se anula en las expresiones
de dh y du
2.- Cp y Cv son funciones de P,v y T en los gases reales, de modo que la
evaluación de la integral en las ecuaciones es complicada.
3.- Generalmente no se conoce la trayectoria a lo largo del cual se
desarrolla el proceso.
Se desarrolla un método alternativo aprovechando que las propiedades
u, h y s son funciones de estado que no dependen de la trayectoria.
 f
2 4
3
i

htotal=hf-hi
htotal = h1+ h2+ h3+ h4
El método consiste en descomponer la trayectoria del proceso en
tramos o etapas.
1.- La variación total de la propiedad será la suma de la variación de dicha
propiedad de cada tramo.
2.- Puesto que Cp y Cv solo se conocen para un gas ideal a presiones
bajas, con frecuencia se elige el tramo donde ocurre variación de
temperatura como el de bajas presiones.
Así un proceso cualesquiera se puede descomponer en 3 etapas.

2
. P2
Entalpia:  

P1 .
1   (1)

Entropía:  
P*
a b
  (2)
T1 T2
Tramo 1: (1-a) proceso isotérmico a T1.
Ecuaciones: Integrando:
.
2
P2

   

 
P1
.
1

Tramo 2: (a-b) proceso isobárico a P*. P*

a b
    T1 T2

   

Tramo 3: (b-2) proceso isotérmico a T2.

 
 
La variación total de la propiedad
  (3)

  (4)
  
 
    (5)

Funciones de desviación desviación de la entalpia

similarmente
 
   
(6)

Función de desviación de la entropía

Las funciones de desviación pueden ser obtenidas a partir de la relación

de los datos PVT.
 Funciones de desviación
Por ejemplo: Desviación de la entalpia

A temperatura constante.

 
(7)

Las ecuaciones de estado pueden ser explícitos o implícitos en V.

Por ejemplo: explícitos en V:

(
𝑃𝑣 = 𝑅𝑇 + 𝑏 −
𝑎
𝑇)𝑃

Hallar la derivada del volumen respecto a T es simple:

𝑅𝑇 𝑎
𝑣= +𝑏 −
𝑃 𝑇

( )
𝜕𝑣 𝑅 𝑎
= + 2
𝜕𝑇 𝑃 𝑃 𝑇
 Por lo general, en una ecuación de estado, el volumen esta implícito, por lo
tanto hallar la derivada del volumen respecto a la temperatura resulta
complicado. Por ejemplo la ecuación de Van der Waals.

𝑅𝑇 𝑎
𝑃= − 2
( 𝑣 −𝑏) 𝑣

( )
𝜕𝑣
𝜕𝑇 𝑃
=?

Sin embargo, la termodinámica provee otras funciones para hallar la

solución. Por ejemplo: la energía libre de Helmholtz (ecuación 4).

𝑑𝐴=− 𝑠𝑑𝑇 − 𝑃𝑑𝑣

A temperatura constante.
𝑑𝐴=− 𝑃𝑑𝑣
Integrando la expresión desde el comportamiento ideal hasta una presión P.
𝑅𝑇 𝑅𝑇
𝑑𝐴=− 𝑃𝑑𝑣 + 𝑑𝑣 − 𝑑𝑣
𝑣 𝑣

 
(8)
Otras funciones de desviación:
 
 
 
 
  (9)

Desviación de la entropía: a volumen constante  

  (10)

Integrando la ecuación (8)

 
(11)

  (12)

Las funcione de desviación pueden ser obtenidas usando las relaciones

P,V,T: Ecuación Virial, Ecuaciones Cubicas y Correlaciones generalizadas.
 RT a
Por Ejemplo: Para la ecuación de R-K: P  2
v  b v  bv
Reemplazando en (8)
 

Integrando

𝐴 − 𝐴∗ =− 𝑅𝑇𝑙𝑛 ( 𝑣 − 𝑏 ) −
𝑎
𝑏 (
𝑙𝑛
𝑣)
𝑣 +𝑏
( )
𝑣
− 𝑅𝑇𝑙𝑛𝑣 − 𝑅𝑇𝑙𝑛 ∗
𝑣

∗
𝐴− 𝐴 =− 𝑅𝑇𝑙𝑛 ( 𝑣−𝑏
𝑣
𝑎
− 𝑙𝑛
𝑏 )
𝑣+𝑏
𝑣 ( 𝑣
− 𝑅𝑇𝑙𝑛 ∗
𝑣) ( ) (a)

Ahora, reemplazando en (10):

[ ]
∗
∗ 𝜕 ( 𝐴− 𝐴 )
𝑠 −𝑠 =−
𝜕𝑇 𝑣

Luego,
∗
𝑠 −𝑠 =𝑅𝑙𝑛 ( 𝑣 −𝑏
𝑣
− ) 𝑎
2 𝑏𝑇
𝑙𝑛
𝑣 +𝑏
𝑣 ( 𝑣
+ 𝑅𝑙𝑛 ∗
𝑣 ) ( ) (b)
 Luego, h=𝑢+ 𝑃𝑣
A
h= 𝐴+ 𝑇𝑠+ 𝑃𝑣
h − h ∗ = ( 𝐴 − 𝐴∗ ) +𝑇 ( 𝑠 − 𝑠∗ ) + 𝑃 (𝑣 −𝑣 ∗ ) (c)

Reemplazando (a) y (b) en ( c)

∗
h −h =− 𝑅𝑇𝑙𝑛 ( 𝑣 −𝑏
𝑣 ) − 𝑙𝑛 (
𝑎
𝑏
𝑣+ 𝑏
𝑣 ) 𝑣
𝑣 ∗
( 𝑣 2 𝑏𝑇 𝑣
− 𝑅𝑇𝑙𝑛 ( )+𝑇 𝑅𝑙𝑛 (
𝑣 −𝑏
) −
𝑎
𝑙𝑛 (
𝑣 +𝑏
) 𝑣 )
+ 𝑅𝑙𝑛 ( ) + 𝑃 ( 𝑣 − 𝑣 )
𝑣
∗
∗

h −h =−
∗ 3𝑎
2𝑏
𝑙𝑛 (
𝑣+𝑏
𝑣 )
+ ( 𝑍 − 1 ) 𝑅𝑇 (d)

O también

h −h =−
∗ 3𝑎
2𝑏
𝑙𝑛 (
𝑣+𝑏
𝑣
−
𝑎
)
+
𝑅𝑇𝑏
(𝑣+ 𝑏) (𝑣 − 𝑏)
Asimismo:
𝑢 −𝑢∗ =( 𝐴 − 𝐴 ∗ )+ 𝑇 ( 𝑠 − 𝑠∗ )
Reemplazando (a) y (b)

𝑢 −𝑢∗ =−
3𝑎
2𝑏
𝑙𝑛 (
𝑣 +𝑏
𝑣 ) (e)

Además, 𝐺 − 𝐺 =( 𝐴 − 𝐴 ) + 𝑅𝑇 ( 𝑍 − 1) (f)
∗ ∗
 RT a
Para la ecuación de Soave-R-K: P  2
v  b v  bv
( h − h∗) =( 𝐴 − 𝐴∗ ) +𝑇 ( 𝑠 − 𝑠 ∗ ) + 𝑅𝑇 ( 𝑍 −1)
( h − h∗ ) = 𝑎 𝑙𝑛
𝑏 ( 𝑍
𝑍+𝐵
∗ )
−
𝑇 𝜕𝑎
𝑏 𝜕𝑇
𝑙𝑛
𝑍
(
𝑍+ 𝐵
∗ )
+ 𝑅𝑇 ( 𝑍 − 1) (g)

( )
0.5
𝜕𝑎 0.42748 𝑎 𝑃 𝑐
=− 𝑅 𝑓 𝑤
𝜕𝑇 𝑇 𝑇𝑐

𝑓 𝑤 =0.48+ 1.574 𝑤 −0.176 𝑤 2

(h)
 Propiedades residuales
Se definen como la diferencia entre el valor de la propiedad de un gas real y el valor
de la propiedad de un gas ideal, la misma presión, temperatura y composición.

Asi, cualquier propiedad termodinámica (V,H,U,S,G, etc) a T y P del sistema.

𝑀 𝑅 =𝑀 − 𝑀 𝑔𝑖

h 𝑅 =h − h𝑔𝑖 𝐺 𝑅 =𝐺 − 𝐺𝑔𝑖
𝑣 𝑅 =𝑣 − 𝑣 𝑔𝑖

𝑠 𝑅=𝑠 − 𝑠 𝑔𝑖

Aplicable a sustancias puras y mezclas.

Ejercicio 1:
Derivar las funciones de desviación para una sustancia pura usando la
ecuación de Redlich–Kwong.
Ejercicio 2:
Derivar las funciones de desviación para una sustancia pura usando la
ecuación de Van der Waals.
Ejercicio 3:
Hallar el calor de vaporización del isobutano a 15.72 bar y 360.93 K
usando las funciones de desviación para la ecuación de Redlich-Kwong.
El valor experimental es 13653 J/mol.
Ejercicio 4:
Hallar el calor de vaporización del isobutano a 15.72 bar y 360.93 K
usando las funciones de desviación para la ecuación de Van der Waals.
El valor experimental es 13653 J/mol.
Ejercicio 5:
Resolver el ejercicio 4, pero usando las funciones de desviación para la
ecuación de Soave.
Ejercicio 6:
Encuentre los valores para la entalpia residual y la entropía residual para
el gas n-butano a 500 K y 50 bar usando la ecuación de R-K.
 Solución 3:
Hallar el calor de vaporización del isobutano a 15.72 bar y 360.93 K usando las
funciones de desviación para la ecuación de Redlich-Kwong. El valor
experimental es 13653 J/mol.

La desviación de la entalpia para R-K

es:
h −h ∗ =−
3𝑎
2𝑏
𝑙𝑛( 𝑣+𝑏
𝑣 )+ ( 𝑍 − 1 ) 𝑅𝑇

Como: ∆ h 𝑣 =h𝑔 −h 𝑓

h𝑔 − h∗ =−
3𝑎
2𝑏
𝑣 +𝑏
𝑙𝑛 𝑔
𝑣𝑔 ( )
+ ( 𝑍 𝑔 −1 ) 𝑅𝑇

∗
h 𝑓 −h = −
3𝑎
2𝑏
𝑙𝑛
𝑣𝑓(
𝑣𝑓 +𝑏
)
+ ( 𝑍 𝑓 −1 ) 𝑅𝑇

Luego,
Solución 6:
Encuentre los valores para la entalpia residual y la entropía residual para
el gas n-butano a 500 K y 50 bar usando la ecuación de R-K.
h 𝑅 =h − h∗
La desviación de la entalpia para R-K
es: ∗
h −h =−
3𝑎
2𝑏
𝑙𝑛 (
𝑣 𝑔+ 𝑏
𝑣𝑔 ) + ( 𝑍 𝑔 − 1 ) 𝑅𝑇

Similarmente,
𝑠 𝑅=𝑠 − 𝑠 ∗

𝑠 −𝑠 ∗ = 𝑅𝑙𝑛
( 𝑣𝑔 − 𝑏
𝑣𝑔
−
)
𝑎
2𝑏𝑇
𝑣 +𝑏
𝑙𝑛 𝑔
𝑣𝑔 ( 𝑣
+ 𝑅𝑙𝑛 ∗𝑔
𝑣 ) ( )
𝑣 ∗ Corresponde a P = 1 bar
Isobutano
T 360.93 Tc 408.1
P 15.72 Pc 36.5
R 83.14 w 0.176
videal 1908.88805
fw 0.75086822
Tr 0.88441558 Z3-Z2+(A-B-B2)Z-AB
a 14715713.6
b 80.538251 0.21293094
A* 0.25690222 0.01083901
B* 0.04219119
Z 0.075 0.728
F(obj) -7.2305E-05
Vg 1389.6705
Vf 143.166604

Term1 71256.1557
Term3 0.38998083
Term2 11045.7381
hv 95158.823
9515.8823 J/mol
Funciones de desviación : Correlaciones Generalizadas

( ) ( )
Sea 𝜕𝑣 𝑍𝑅 𝑅𝑇 𝜕 𝑍
𝑃𝑣 𝑍𝑅𝑇 = +
𝑍= ; 𝑣= ;
𝑅𝑇 𝑃 𝜕 𝑇 𝑃 𝑃 𝑃 𝜕𝑇 𝑃

Desviación de la entalpia: reemplazando en la ecuación de la entalpia

[ ( )] ( ) 𝑑𝑃
2
𝜕𝑣 𝑅𝑇 𝜕 𝑍
𝑑h = 𝑣 − 𝑇 𝑑𝑃=−
𝜕𝑇 𝑃 𝑃 𝜕𝑇 𝑃

𝑃 ( 𝜕𝑇 )
2
𝑅𝑇 𝜕 𝑍
h −h =−∫
Integrando ∗
𝑑𝑃 (13)
𝑃
∗ 𝑃

Introduciendo las propiedades reducidas

𝑃 2❑

( )
∗
h −h 𝑇𝑟 𝜕 𝑍
=∫ 𝑑 𝑃𝑟 (14)
𝑅 𝑇 𝑐 𝑃 𝑃 𝑟 𝜕𝑇 𝑟
∗ 𝑃𝑟
 Desviación de la entropía: similarmente

𝑑𝑠=− ( )
𝜕𝑣
𝜕𝑇
𝑑𝑃
como ( )
𝜕𝑣
=
𝜕 𝑇 𝑃 𝑃 𝑃 𝜕𝑇( )
𝑍𝑅 𝑅𝑇 𝜕 𝑍
+
𝑃
𝑃

𝑑𝑠=−(
𝜕𝑣
𝜕𝑇 ) 𝑑𝑃=−
[𝑃𝑃
+ (
𝑍𝑅 𝑅𝑇 𝜕 𝑍
𝑃 𝜕𝑇 ) ] 𝑑𝑃
𝑃

Integrando e introduciendo las propiedades reducidas

[ ( ) ] 𝑑𝑙𝑛 𝑃
∗ 𝑃
𝑠 −𝑠 𝜕𝑍
=∫ 𝑍 +𝑇 𝑟 𝑟
(15)
𝑅 𝑃
∗ 𝜕 𝑇𝑟 𝑃𝑟

Las funciones de desviación:

h∗ − h y 𝑠∗ − 𝑠
𝑅𝑇𝑐 𝑅

Se pueden representar en graficas o tablas en función de las propiedades

reducidas Pr y Tr
Método de Lee-Kesler
Propiedad residual 𝑀 𝑅 =𝑀 − 𝑀 𝑔𝑖
Ejemplo: h 𝑅 =h − h∗
𝑠 𝑅=𝑠 − 𝑠 ∗
𝑣 𝑅 =𝑣 − 𝑣 ∗
Desviación de la entalpia residual:

( ) ( )
𝑜 1
h𝑅 h𝑅 h𝑅 (16)
= +𝑤
𝑅 𝑇𝑐 𝑅 𝑇 𝑐 𝑅𝑇𝑐
están dadas en graficas y tablas en función de Pr y Tr

Desviación de la entropía residual:

( ) ( ) ( )
𝑅 𝑅 𝑜 𝑅 1 ∗
𝑠 𝑠 𝑠 𝑃 (17)
= +𝑤 + ln
𝑅 𝑅 𝑅 𝑃
están dadas en graficas y tablas en función de Pr y Tr
 Ejercicio 7
Estime la desviación de la entalpia y entropía para propileno a 398,15
K y 100 bar. Bier et al reportaron los siguientes valores experimentales:

h∗ − h=244,58 J/g 𝑠 ∗ − 𝑠=1,4172 J/g K

DATOS:
M = 42,081 Tc = 364,9 K Pc= 46 bar Zc = 0,274 w = 0,144

Usando:
a) La ecuación de R-K y de VDW
b) Las graficas generalizadas
c) El método de Lee-Kesler
Solución Desviación de la entalpia

( )
Propileno ∗ 3𝑎 𝑣+𝑏 𝑎 𝑅𝑇𝑏
T 398.15
h −h =− 𝑙𝑛 − +
2𝑏 𝑣 (𝑣+ 𝑏) (𝑣 − 𝑏)
P 100
Tc 364.9 91392.33412 bar cc/mol
Pc 46
R 83.14 217,182 J/g
M 42.081
217,182 J/g
a 8188216.324
b 57.14056041 Desviación de la entropía
A* 0.74726704
B* 0.172618666
Ecuacion cubica
Z3-Z2+(A*-B*-B*2)-
∗
𝑠 −𝑠 =𝑅𝑙𝑛 ( 𝑣 −𝑏
𝑣
−) 𝑎
2 𝑏𝑇
𝑙𝑛
𝑣 +𝑏
𝑣 ( ) 𝑣
+ 𝑅𝑙𝑛 ∗
𝑣 ( )
A*B*=0 Term1 -42.6696456
Z 0.43 Term2 -60.7352876
F(obj) -9.92359E-05 Term3 -453.041521
v ideal 331.02191 s-s* -556.446454
Vgas 142.3394213
v* = 33102.191
𝑠 −𝑠 ∗ =−556.446454 𝑏𝑎𝑟 𝑐𝑐 /𝑚𝑜𝑙𝐾
Term 1 -72545.26427
Term2 -41047.80968 J/g K
Term3 22200.73983
h-h* -91392.33412 J/g K
b) Usando gráficos Desviación de la entalpia
Tr 1.09112086
Pr 2.17391304 cal/mol

cal/mol k
cal/mol
J/g
Desviación de la entropía
cal/mol K grafico

Utilizando la ecuación (17)

cal/mol K

J/g K
c) Método de Lee-Kesler Desviación de la entalpia

( ) ( )
𝑅 𝑅 𝑜 𝑅 1
Tr 1.09112086 h h h
Pr 2.17391304 = +𝑤
𝑅 𝑇𝑐 𝑅 𝑇 𝑐 𝑅𝑇𝑐
Interpolando de tablas
( )
1

( ) h𝑅
𝑜
h𝑅 =−1.4186
=−3.0325
𝑅 𝑇𝑐 𝑅 𝑇𝑐
Pr (h/RTc)º (h/RTc)1
2 -2.965 -1.261
2.174 h𝑅
=− 3.0325+0.144 ( − 1.4186 ) =−3.2368
3 -3.353 -2.167 𝑅 𝑇𝑐
1 -0.388 -0.906
0.174
(h/RTc)º -3.032512
h 𝑅 =−3.2368 𝑅 𝑇 𝑐
(h/RTc)1 -1.418644
J/g
Interpolando de tablas
Pr (s/R)º (s/R)1 J/g
2 -2.081 -1.241
2.174 Desviación de la entropía

( ) ( ) ( )
𝑅 𝑅 𝑜 𝑅 1 ∗
3 -2.202 -2.067
𝑠 𝑠 𝑠 𝑃
1 -0.121 -0.826 = +𝑤 + ln
0.174 𝑅 𝑅 𝑅 𝑃
(s/R)º -2.102054

( ) ( )
𝑜 1
(s/R)1 -1.384724 𝑠𝑅 𝑠𝑅
=− 2.10205 =−1.3847
𝑅 𝑅 𝑅
𝑠 =− 6.90662 𝑅
𝑠𝑅 1
J/mol k =−2.10205+ 0.144(− 1.3847)+ ln ⁡( )
𝑅 100
J/g K
 Procesos Isentrópicos
Requiere 3 etapas

∆ h 1− 2=∆ h1 −𝑎 +∆ h𝑎 −𝑏 +∆ h𝑏−2
∆ h 1− 2=( h − h 𝑃 1) 𝑇 1 +( h2 − h1 ) − ( h −h 𝑃 2) 𝑇 2
∗ ∗ ∗
(18)

Donde: =

Similarmente

∆ 𝑆1 −2=∆ 𝑆 1− 𝑎+∆ 𝑆 𝑎− 𝑏 +∆ 𝑆 𝑏−2

∆ 𝑆1 −2=( 𝑆 −𝑆 𝑃 1 )𝑇 1 +(𝑆 2 − 𝑆1 ) − ( 𝑆 − 𝑆 𝑃 2 )𝑇 2
∗ ∗ ∗ (19)

Donde:
=

Otra alternativa es usar fugacidades en vez de presiones

Se sabe h=𝑢+𝑃𝑣
𝑑h =𝑑𝑢+𝑃𝑑𝑣 +𝑣𝑑𝑃 ; 𝑑h =𝑇𝑑𝑠+𝑣𝑑𝑃
𝜕h 𝑣 (20)
𝑑𝑆= − 𝜕 𝑃
𝑇 𝑇
Integrando desde el estado ideal a un estado real:
∗ 𝑃
( 𝑆 − 𝑆 )= h − h −∫ 𝑣 𝑑𝑃
∗
(21)
𝑇 𝑃
𝑇 ∗

Para resolver el termino integral se recurra la energía libre de Gibbs, donde:

𝑑𝐺=𝑣𝑑𝑃
Para gas ideal 𝑑 𝐺∗=𝑣∗ 𝑑𝑃 =𝑅𝑇𝑑𝑙𝑛𝑃
Para gas real 𝑑𝐺=𝑣𝑑𝑃=𝑅𝑇𝑑𝑙𝑛 𝑓 (22)

Donde: 𝑓 = 𝑓𝑢𝑔𝑎𝑐𝑖𝑑𝑎𝑑
 𝑃

Luego, ∫ 𝑇𝑣 𝑑𝑃=− 𝑅𝑙𝑛 𝑓𝑓∗

𝑃
∗
( ) (23)

Reemplazando en (25)

h − h∗
( )
𝑓 (24)
( 𝑆 − 𝑆 )=
∗
− 𝑅𝑙𝑛 ∗
𝑇 𝑓

Luego para un proceso de tres etapas

∆ 𝑆1 −2=( 𝑆 −𝑆 𝑃 1 )𝑇 1 +(𝑆 2 − 𝑆1 ) − ( 𝑆 − 𝑆 𝑃 2 )𝑇 2
∗ ∗ ∗

∆ 𝑆1 −2=( h∗ − h 𝑃 1 )𝑇 1+(𝑆2 − 𝑆1)∗ − ( h∗ − h 𝑃 2 )𝑇 2

Donde:

= (25)
Ejercicio 8
Se comprime Helio en un compresor isotérmico reversible. El Helio
ingresa a 300 K y 2000 kPa y se comprime hasta 12000 kPa. Calcule
la potencia requerida por el compresor y la cantidad de calor
desprendido
a) Si el gas se comporta como gas ideal
b) Si el gas satisface la ecuación de estado:

Donde:
(
𝑃𝑣 = 𝑅𝑇 + 𝑏 −
𝑎
𝑇)𝑃

∗ 𝑘𝐽 𝑚3 𝐾 𝑚3
𝐶 𝑃 =20,8 𝑎=0,355 𝑏=0,0152
𝑘𝑚𝑜𝑙 𝐾 𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙

Ejercicio 9
Si el compresor del ejercicio anterior opera en forma adiabática
reversible, Calcule la temperatura de salida y la potencia que debe
suministrarse para obtener 12000 kPa de presión de salida.
Ejercicio 10
Ingresa etano a una tobera a 5000 kPa, 600 K y 180 m/s. El gas se
expande isotérmicamente hasta 1000 kPa. Determine la velocidad
máxima de salida del gas.
a) Si el gas se comporta como gas ideal
b) Si el gas satisface la ecuación de estado:
𝑃𝑣=𝑅𝑇 +𝑎𝑣
Donde: kPa

Ejercicio 11
Si la tobera del ejercicio anterior opera en forma adiabática
reversible, Calcule la temperatura de salida y la velocidad máxima
de salida del gas. 𝑘𝐽
∗
𝐶 𝑃 =52,5
𝑘𝑚𝑜𝑙 𝐾

a) Si el gas se comporta como gas ideal

b) Si el gas satisface la ecuación de estado:
𝑃𝑣=𝑅𝑇 +𝑎𝑣
Ejercicio 8 Proceso Isotérmico
Se comprime Helio en un compresor isotérmico reversible. El Helio
ingresa a 300 K y 2000 kPa y se comprime hasta 12000 kPa. Calcule
la potencia requerida por el compresor y la cantidad de calor
desprendido
a) Si el gas se comporta como gas ideal
b) Si el gas satisface la ecuación de estado:

Donde:
𝑃𝑣 = 𝑅𝑇 + 𝑏 − ( 𝑎
𝑇)𝑃

∗ 𝑘𝐽 𝑚3 𝐾 𝑚3
𝐶 𝑃 =20,8 𝑎=0,355 𝑏=0,0152
𝑘𝑚𝑜𝑙 𝐾 𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙
Solución:
𝑇 1 =300 𝐾 𝑃 2=12000 𝑘𝑃𝑎
compresor
𝑃 1=2000 𝑘𝑃𝑎

1.- Balance de energía: Proceso FEES

( ) ( )
2 2
𝑉𝑒 𝑉𝑠 ˙ 𝑊˙ =0
h 𝑒+ +𝑔 𝑧 𝑒 𝑚
˙ 𝑒 − h𝑠+ + 𝑔 𝑧𝑠 𝑚
˙ 𝑠 + 𝑄+
2 2

Luego, por unida de masa ˙ =−𝑞+ ∆ h1 −2

𝑊 (1)
2.- Balance de entropía: Proceso FEES
❑ ˙
❑
𝑞
∑ 𝑠𝑠 𝑚𝑠 − ∑ 𝑠𝑒 𝑚 𝑒= 𝑇
˙
𝑠 𝑠
𝑞 (2)
∆ 𝑠1 − 2=
𝑇
Si se calcula se determina el calor desprendido, luego calculando y usando la
ecuación (1) se obtiene la potencia del compresor.
Como no hay cambio de temperatura, entonces el proceso será de 1 etapa.
Ecuaciones

[ ( )]
𝑃2
𝜕𝑣
𝑑h =𝐶 𝑃 𝜕 𝑇 + 𝑣 − 𝑇 𝜕𝑃
𝜕𝑇 𝑃

𝑃1
𝑇 =300 𝐾
𝑑𝑠=𝐶 𝑃
𝑇 ( ) 𝜕𝑃
𝜕𝑇 𝜕 𝑣
−
𝜕𝑇 𝑃

Para proceso isotérmico 𝜕𝑇 =0 [ 𝜕𝑇 ) ]𝜕 𝑃

𝑑h = 𝑣 − 𝑇 (
𝜕𝑣
𝑃
(3)

𝑑𝑠=−( )
𝜕𝑣
𝜕𝑃 (4)
𝜕𝑇 𝑃
a) Para Gas ideal: 𝑃𝑣=𝑅𝑇

( )
𝜕𝑣
=
𝜕𝑇 𝑃 𝑃
𝑅

Reemplazando en (1)

[
𝑑h = 𝑣 −
𝑅𝑇
𝑃 ]
𝜕 𝑃= 0
∆ h 1− 2=0 (5)

Calculo de
𝑃2
𝑅 ∆ 𝑠1 − 2=− Rln ⁡( ) (6)
𝑑𝑠 =− 𝑑 𝑃 𝑃1
𝑃

∆ 𝑠1 − 2=−14,896

Usando (2)
𝑞=300
˙ ( −14,896 )=− 4468,8
Reemplazando en (1)
kJ/kmol

b) Para Gas real: (

𝑃𝑣 = 𝑅𝑇 + 𝑏 −
𝑎
𝑇
𝑃) 𝑣=
𝑅𝑇
𝑃
+ 𝑏−(𝑎
𝑇 )
( )
𝜕𝑣 𝑅 𝑎
= + 2
𝜕𝑇 𝑃 𝑃 𝑇
Reemplazando en (1) [
𝑑h = 𝑣 −
𝑅𝑇
𝑃
−
𝑎
𝑇
𝜕𝑃
]
𝑑h = [ 𝑅𝑇
𝑃
+ 𝑏−
𝑎
𝑇(−
𝑅𝑇
𝑃
−
𝑎
𝑇
𝜕𝑃 ) ]
[ 2𝑎
] [ ]
𝑃2
𝑑h = 𝑏 − 𝜕𝑃 2𝑎
𝑇 ∆ h 1− 2=∫ 𝑏− 𝜕𝑃
𝑃1
𝑇

∆ h 1− 2= 𝑏 − ( 2𝑎
𝑇 )
( 𝑃2 − 𝑃1) (6)

Reemplazando datos:

(
∆ h 1− 2= 0,0152 −
2(0,355)
300
(12000 −2000) )
∆ h 1− 2=128,33 𝑘𝐽 / 𝑘𝑚𝑜𝑙
Calculo de
𝑑𝑠=− ( ) 𝜕𝑣
𝜕𝑇 𝑃
𝜕𝑃 𝑑𝑠 =−
( 𝑅 𝑎
+
𝑃 𝑇2 ) 𝜕𝑃
𝑃

𝑃2

∆ 𝑠1 − 2=∫ −
𝑃1
( 𝑅 𝑎
+
𝑃 𝑇2 ) 𝜕𝑃
𝑃
Reemplazando datos: (7)

∆ 𝑠1 − 2=−14,936 𝑘𝐽 /𝑘𝑚𝑜𝑙𝐾

Luego,
kJ/kmol

Reemplazando en (2):

kJ/kmol
 Ejercicio 9 Proceso Isentropico
Si el compresor del ejercicio anterior opera en forma adiabática
reversible, Calcule la temperatura de salida y la potencia que debe
suministrarse para obtener 12000 kPa de presión de salida.
a) Si el gas se comporta como gas ideal
b) Si el gas satisface la ecuación de estado:
Donde:
𝑘𝐽
𝑃𝑣 = 𝑅𝑇 + 𝑏 − ( 𝑚 𝐾
𝑎
)
𝑇 3
𝑃
𝑚3
∗
𝐶 𝑃 =20,8 𝑎=0,355 𝑏=0,0152
𝑘𝑚𝑜𝑙 𝐾 𝑘𝑚𝑜𝑙 𝑘𝑚𝑜𝑙
Solución:
𝑇 1 =300 𝐾 compresor
𝑃 2=12000 𝑘𝑃𝑎
𝑃 1=2000 𝑘𝑃𝑎

1.- Balance de energía: Proceso FEES

( ) ( )
2 2
𝑉𝑒 𝑉𝑠 ˙ 𝑊
h 𝑒+ +𝑔 𝑧 𝑒 𝑚
˙ 𝑒 − h𝑠+ + 𝑔 𝑧𝑠 𝑚
˙ 𝑠 + 𝑄+ ˙ =0
2 2

˙ =( h𝑠 − h𝑒) 𝑚
𝑊 ˙

Luego, por unida de masa ˙ =∆ h 1− 2

𝑊 (1)
2.- Balance de entropía: Proceso FEES
0
∆ 𝑠1 − 2=0 (2)

proceso isentropico

a) Gas ideal 𝑃𝑣=𝑅𝑇 ( )

𝜕𝑣
=
𝜕𝑇 𝑃 𝑃
𝑅

𝑑𝑠=𝐶 𝑃
𝜕𝑇 𝜕 𝑣
𝑇
−
𝜕𝑇( ) 𝜕𝑃 𝑃
𝜕𝑇 𝑅

( )
12 8,314 /20,8
𝑑𝑠=𝐶 𝑃 − 𝜕𝑃
𝑇 𝑃 𝑇 2=300  
2
𝜕𝑇 𝑅
𝐶𝑃 = 𝜕𝑃 K
𝑇 𝑃
Integrando Calculo de ∆ h
𝑇2

𝐶 𝑃 𝑙𝑛( )
𝑇2
𝑇1
=𝑅𝑙𝑛
𝑃2
𝑃1 ( ) ∆ h=∫ 𝐶 𝑃 𝑑𝑇 =20,8(𝑇 2 − 𝑇 1 )
𝑇1

kJ/kmol

( )
𝑅/ 𝐶 𝑃
𝑃2
𝑇 2=𝑇 1   Reemplazando en (1)
𝑃1
6537,5 kJ/kmol
 b) Gas real 𝑃𝑣 = 𝑅𝑇 + 𝑏 −( 𝑎
𝑇
𝑃 )
.
2

( )
𝜕𝑣 𝑅 𝑎
= + 2
𝜕𝑇 𝑃 𝑃 𝑇
P2

Escogemos un diagrama P-T con 3 etapas debido a que

la temperatura cambia de T1 a T2.
P1
.
1

[ ( )]
Ecuaciones:
𝜕𝑣
𝑑h =𝐶 𝑃 𝜕 𝑇 + 𝑣 − 𝑇 𝜕𝑃
𝜕𝑇 𝑃
P*
a b

( ) 𝜕𝑃
T1 T2
𝜕𝑇 𝜕 𝑣
𝑑𝑠=𝐶 𝑃 −
𝑇 𝜕𝑇 𝑃

Etapa 1: (1-a) proceso isotérmico a T1. 𝑑𝑠=− ( ) 𝜕𝑃

𝜕𝑣
𝜕𝑇 𝑃

( )
𝑃∗

( )
∗
𝑅 𝑎 𝑃 𝑎
∆ 𝑠1 − 𝑎=− ∫ +   𝜕 𝑃=− 𝑅𝑙𝑛 − 2 (
𝑃 ∗ − 𝑃 1)
𝑃1 𝑃 𝑇1 2
𝑃1 𝑇1

Etapa 2: (a-b) proceso isobárico a P*.

( )
𝑇2
𝜕𝑇 𝑇
∆ 𝑠𝑎 − 𝑏=∫ 𝐶 𝑃 =𝐶 𝑃 𝑙𝑛 2
𝑇1 𝑇 𝑇1

Etapa 3: (b-2) proceso isotérmico a T2. ∆ 𝑠𝑏 − 2=− 𝑅𝑙𝑛

( )
𝑃2
𝑃
∗
−
𝑎
𝑇2
2(
𝑃 2 − 𝑃 ∗)
Sumando ∆ 𝑠1 − 2=∆ 𝑠1 − 𝑎+∆ 𝑠 𝑎 −𝑏 +∆ 𝑠 𝑏−2

( ) 𝑇2 𝑃2
( ) ( )
∗
𝑃 𝑎 𝑎
− 2 ( 𝑃 − 𝑃 1) + 𝐶 𝑃 𝑙𝑛 − 𝑅𝑙𝑛 ∗ − 2 ( 𝑃 2 − 𝑃 )
∗ ∗
∆ 𝑠1 − 2=− 𝑅𝑙𝑛
𝑃1 𝑇1 𝑇1 𝑃 𝑇2

0=− 𝑅𝑙𝑛
( ) (
𝑃2
𝑃1
−𝑎
𝑃2
𝑇2
2
−
𝑃1
𝑇1
2
+𝐶 𝑃 𝑙𝑛
𝑇2
𝑇1) ( )
Reemplazando datos,

−8.314 𝑙𝑛
12
2 ( )
− 0,355
12000 2000
𝑇2
2
−
300
2
(
+20,8 𝑙𝑛
𝑇2
300
=0
) ( )
Resolviendo:
K
Calculo de ∆ h
Etapa 1: (1-a) proceso isotérmico a T1.
[ ( )]
𝑑h = 𝑣 − 𝑇
𝜕𝑣
𝜕𝑇 𝑃
𝜕𝑃

𝑑h =
[ 𝑅𝑇
𝑃
+𝑏 −
𝑎 𝑅𝑇 𝑎
𝑇1
−
𝑃
−
𝑇1
𝜕𝑃
]
∆ h 1− 𝑎= 𝑏 −
[ 2𝑎
𝑇1 ] ( 𝑃∗ − 𝑃 1)
 Etapa 2: (a-b) proceso isobárico a P*
𝑑h =𝐶 𝑃 𝜕 𝑇
∆ h 𝑎− 𝑏 =𝐶 𝑃 (𝑇 2 −𝑇 1)

Etapa 3: (b-2) proceso isotérmico a T2.

[
∆ h 𝑏 −2 = 𝑏 −
2𝑎
𝑇2 ]
( 𝑃2− 𝑃 )
∗

Sumando: ∆ h 1− 2=∆ h1 −𝑎 +∆ h𝑎 −𝑏 +∆ h𝑏−2

∆ h 1− 2= 𝑏−
[ 2𝑎
𝑇1
( 𝑃
]
∗
− 𝑃 1) +𝐶 𝑃 (𝑇 2 −𝑇 1 )+ 𝑏 −
2𝑎
𝑇2[ ( ]
𝑃2 − 𝑃 )
∗

∆ h 1− 2=𝑏 ( 𝑃 2 − 𝑃 1 ) −2 𝑎
𝑃2 𝑃1
−
𝑇2 𝑇1 ( ) +𝐶 𝑃 (𝑇 2 − 𝑇 1)
Reemplazando datos,

=0,152 ( 12000 − 2000 ) −2(0,355) (

300 )
12000 2000
∆ h 1− 2 − +20,8(614 −300)
614
kJ/kmol

6676,1 kJ/kmol