PROBLEMA 1

Se comprime gas propano desde 5 bar y 105 ºC hasta su estado final a 25 bar y 190
ºC. Determine empleando la ecuación cúbica de SRK:
a) El cambio de entalpia para el proceso en kJ/kmol
b) El requerimiento energético en kJ, si la tasa de flujo másico es 500 kg/h
SOLUCIÓN
Paso 1
Sustancia= propano
T 1=105 ºC

P1=¯¿

T 2=190 ºC

P2=25 ¯¿

Paso 2 Se hace la conversión de la temperatura de entrada a grados kelvin utilizando

la tabla de conversiones. (M. Koretsky, 2013, p.2)
T 1=105 ºC+273.15=378.15 K

Paso 2 Se hace la conversión de la temperatura de salida a grados kelvin utilizando la

tabla de conversiones. (M. Koretsky, 2013, p.2)

T 2=190 ºC+273.15=463.15 K

Paso 3 Se representa el diagrama de flujo del proceso

Paso 3

dm
=∑ ṁ¿−¿ ∑ ṁout +¿ ∑ ṁgen−¿ ∑ ṁcont ¿ ¿ ¿
dt
0=∑ ṁ¿ −∑ ṁout +0−0

∑ ṁ¿=∑ ṁout
ṁ 1=ṁ 2

Paso 3
dE
=∑ Ė ¿ ∑ Ėout + Q̇+ Ẇ
dt
 dE
=∑ ( Ḣ + Ė k + Ė p )1−∑ ( Ḣ + Ėk + Ė p ) 2+ ¿ Q̇+ Ẇ ¿
dt
0=∑ ( Ḣ +0+0 )1−∑ ( Ḣ +0+0 )2+ ¿ Q̇+ 0¿

Q̇= Ḣ 2− Ḣ ¿
^ ṁ H
ṁ Q= ^ 2−ṁ H
^1

^ ṁ( ^
ṁ Q= ^1)
H 2− H
^ =( ^
Q H1)
H 2− ^

Paso 5

RT a
P= −
V −b V (V +b)
Paso Expresado en J m3/mol2
2 2
R TC
[ ]
2
1+ Ω ( 1−T r )
0.5
a=0.42748
PC

Paso Expresado en m3/mol

RT c
b=0.08664
Pc

Paso
2
Ω=0.480+1.574 ω−0.176 ω

Paso

RT a
P= −
V −b V (V +b)
Paso

RT [ V (V + b ) ] −a (V −b )
P=
( V −b ) ( V +b )( V )

RT ( V 2 +Vb )−a ( V −b )
P=
( V 2−b2 ) V
2
RT V + RTVb−aV +ab
P= 3 2
(V −V b )
 P ( V 3−V b2 ) =RT V 2+ RTVb−aV +ab

RT V 2 RTVb aV ab
V 3−V b 2= + − +
P P P P

RT V 2 RTVb aV ab
V 3−V b 2− − + − =0
P P P P

( )
2
RT V
3 a 2 RTb ab
V − + −b − V − =0
P P P P

PV =ZRT
ZRT
V=
P

( ) ( ) RTP +( Pa −b − RTb
P )( P ) P
3 2
ZRT ZRT ZRT ab 2
− − =0
P P

( ) [( ) RTP +( aP −b − RTb ) ] ( P ) P
3 2
ZRT ZRT ZRT ab 2
− −
P P P
=0
(P)
3
ZRT

Z3 −Z 2 +
[ a
P
−b2−
RTb
P
Z ] −
ab
P
=0
( ) ( )
2
RT RT 3
P P

{[ ][ ][ ]} ( )
RTb ab
2
a b P P
Z3 −Z 2 + − − Z− =0
( ) ( ) ( )
2 2 2 3
RT RT RT RT
P
P P P P

{[ ] [ ] [ ]}
2 2 2
3 2 a b P bP ab P
Z −Z + − − Z− =0
( RT )
2
( RT )
2
RT ( RT )
3

3
Z −Z +
2
{[ ] [ ] [ ]}
a
( RT )
2
−
b 2 P2
( RT )
2
−
bP
RT
Z−
aP bP
( RT ) RT
2
=0

aP
A=
( RT )2
 bP
B=
RT

Z3 −Z 2 + { A −B 2−B } Z− A B =0
Paso
2 2
R TC
[ ]
2
1+ Ω ( 1−T r )
0.5
a=0.42748
PC

Paso Se determina el valor de la constante universal de gases ideales en las unidades

adecuadas para el problema
¯¿
R=0.00008314 m3 ¿
mol K
Paso Se extrae las propiedades del componente puro(propano)

Tabla B.1: Propiedades de especies puras

Fuente: (J. Smith, H. Van Ness, M. Abbott & M. Swihart,2018, p. 650)

T c =369.8 K

PC =42.48 ¯
¿

ω=0.152

Paso 6

T
T r=
Tc

Paso para la temperatura de entrada

378.15 K
T r= =1.0226
369.8 K
Paso 6 para la temperatura de salida

463.15 K
T r= =1.2524
369.8 K
Paso
2
Ω=0.480+1.574(0.152)−0.176(0.152)
Paso

Ω=0.7152
Paso para la entrada
 2 2
R TC
[ ]
2
1+ Ω ( 1−T r )
0.5
a=0.42748
PC

a=0.42748 ¿ ¿ ¿
Paso
−6 6 ¯¿
a=9.3600 x 1 0 m 2
¿
mo l
Paso Se realiza la conversión de unidades a J.m3/mol2, de la tabla A.2

Fuente: (J. Smith, H. Van Ness, M. Abbott & M. Swihart,2018, p. 649)

−6 6 ¯
¿ x 1J
a=9.3600 x 1 0 m ¿
mo l 2 J m3
10 m . ¯−5
¿ =0.9360 3
2
¿
mol

Paso Para la salida

a=0.42748 ¿ ¿ ¿
Paso
−6 6 ¯¿
a=7.9607 x 10 m 2
¿
mol
Paso Se realiza la conversión de unidades a J.m3/mol2, de la tabla A.2

Fuente: (J. Smith, H. Van Ness, M. Abbott & M. Swihart,2018, p. 649)

−6 6 ¯
¿ x 1J
a=7.9607 x 10 m ¿
mol 2 ¯−5 J m3
3
10 m . ¿ =0.7961 ¿
mol 2

Paso

RT c
b=0.08664
Pc

Paso

b=0.08664 ¿ ¿
Paso

−5 m3
b=6.2705 x 1 0
mol
Paso

aP
A=
( RT )2

Paso Para la entrada

J m3 ¯
A =0.9360 x5 ¿ ¿
mol 2 ¿ ¿¿

Paso

3 J
A =4.7348 x 1 0
m ¯
3
¿¿

Paso Se multiplica por el equivalente de 1J =10 -5 bar.m3

3 J
A =4.6574 x 10
¯
m3 ¯¿ x 1 0−5 m3 ¿ ¿ ¿
1J

Paso
−2
A =4.7348 x 1 0

Paso Para la salida

J m3 ¯
A =0.7961 x 25 ¿ ¿
mo l 2 ¿¿ ¿

Paso

J
A =1.3423 x 1 04
m ¯¿ ¿
3

Paso Se multiplica por el equivalente de J =10 -5 bar.m3

 4 J
A =1.3423 x 1 0
3 ¯
m ¯¿ x 1 0 m ¿ =¿ ¿ ¿
3 −5
1J

Paso

A =0.1342

Paso

bP
B=
RT

Paso para la entrada

m3 ¯
B =6.2705 x 10 −5
x5 ¿ ¿
mol ¿¿

Paso

B =9.9724 x 1 0−3

Paso para la salida

3
m ¯
B =6.2705 x 10−5 x 25 ¿ ¿
mol ¿¿

Paso

−2
B =4.0711 x 1 0

Paso Para hallar Z1

Z −Z + { A −B −B } Z− A B =0
3 2 2

Z −Z + [ 4.7348 x 10 −( 9.9724 x 1 0 ) −9.9724 x 1 0 ] Z−( 4.7348 x 10

−3 2
3 2 −2 −3 −2
)( 9.9724 x 1 0−3 ) =0
 Paso Se resuelve y la ecuación se reduce a

Z3 −Z 2 +0.0373 Z −4.7217 x 1 0−4=0

Paso se resuelve la ecuación polinómica con la calculadora Casio (fx-991 LA PLUS)

Z1 =0.9617

Z2 =0.0191i

Z3 =0.0191i

Paso 7

Zl =Mínimo

Z g=Máximo

Paso

Z1 =0.9617

Paso Para hallar Z2

Z −Z + { A −B −B } Z− A B =0
3 2 2

Z −Z + [ 0.1342−( 4.0711 x 1 0 ) −4.0711 x 1 0 ] Z− ( 0.1342 ) ( 4.0711 x 10

−2 2
3 2 −2 −2
)=0

Paso Se resuelve y la ecuación se reduce a

Z3 −Z 2 +0.0864 Z+ 5.4232 x 1 0−3=0

Paso se resuelve la ecuación polinómica con la calculadora Casio (fx-991 LA PLUS), obteniendo

Z1 =0.8969

Z2 =0.1448

Z3 =−0.0418

Paso 7

Zl =Mínimo

Z g=Máximo
Paso

Z1 =0.8969

Paso

Zg

Paso
T2

H =H +∫ C giP dT
gi
2
gi
R
T0

T1

H =H +∫ C giP dT
gi
1
gi
R
TO

T2 T1

( H 2 −H 1 )=∫ C P dT −∫ C P dT
gi gi gi gi

TO TO

T2

(H −H ) =∫ C P dT
gi gi gi
2 1
T1

Paso Se escribe la ecuación de la capacidad calorífica del gas ideal

gi
CP 2
= A+ BT + C T
R
Paso Se extrae los valores de A, B, C, D, de la tabla C. 1 de Capacidades térmicas de los gases
en el estado de gas ideal.

Fuente: (J. Smith, H. Van Ness, M. Abbott & M. Swihart,2018, p. 656)

A=1.213

B=2.8785 x 1 0−4
−6
C=−8.8240 x 1 0
kJ
R=8.314
kmol K
Paso

gi
CP −4 −6 2
=1.213+2.8785 x 1 0 T −8.8240 x 1 0 T
R
Paso

[∫ ]
T2
kJ
( H −H )=8.314 kmol
gi
2
gi
1
K
( 1.213+2.8785 x 1 0−2 T −8.8240 x 1 0−6 T 2 ) dT
T1

gi
Paso se calcula el H 1

[ ]
378.15
kJ
gi
H =8.314
1
kmol K
∫ ( 1.213+2.8785 x 1 0−4−8.8240 x 10−6 T 2 ) dT
0

Paso se realiza el cálculo utilizando una calculadora manual Casio

gi kJ
H 1 =19602.2113
kmol
Paso límites de integración desde T1 = 378.15 K a T2=463.15 K
463.15 K

(H gi
2 −H ) =
gi
1 ∫ ( 1.213+2.8785 x 1 0−4 T −8.8240 x 1 0−6 T 2 ) dT
378.15º C

Paso Se distribuye la integral a todos los términos

[∫ ]
190 ºC 190 463.15 K

(H gi
2 −H ) =
gi
1 ( 1.213 ) dT + ∫ ( 2.8785 x 1 0 −4
T ) dT − ∫ ( 8.8240 x 1 0−6 T 2 ) dT
105 ºC 378.15 K 378.15 K

Paso Se resolvió cada integral con una calculadora manual CASIO

( H gi2 −H 1gi)=R [ 103.105+ 3.6089−16.7696 ]

Paso

( H gi2 −H 1gi)=R [ 89.9443 ]

Paso ( H 2 −H 1 ) =[]
gi gi

( H gi2 −H 1gi)=10818.4147

kJ
Paso se reemplaza el valor de R=8.314
kmol K
Paso se calcula la integral de la capacidad calorífica del gas ideal propano, empleando el límite
de integración de T1 = 105 ºC a T2 = 190 ºC
T2

∫ C giP dT
T1

Paso
ρ

( ) dρρ + Z−1
R
H ∂Z
=−T ∫
RT O ∂T ρ

Paso
ρ2
H R2
RT
=−T 2∫
O
∂Z
∂T ( ) dρρ + Z −1
ρ
2

ρ1
H R1
RT
=−T 2∫
O
∂Z
∂T ( ) dρρ + Z −1
ρ
1

Paso

∆ H R=( H R2 −H R1 )

Paso 7

∆ H =( H ig2 −H ig1 ) + ( H R2 −H R1 )

Paso 8

∆ H =−91285.3
kJ
mol (
− −98 565.3
kJ
Kmol )
kJ
∆ H =7 280
kmol