TALLER DE TERMODINÁMICA N°1

TABLAS TERMODINÁMICAS Y ECUACIONES DE ESTADO

KAYRETH DIAZ PEREZ

ANDREA GONZALEZ ASCENCIO
LUISA MENDEZ ROMERO
MARIA PERDOMO SALGADO
ERIKA RAMIREZ OCAMPO

CORPORACION UNIVERSITARIA DEL CARIBE

CECAR

FACULTAD DE CIENCIAS BASICAS, INGENIERIAS Y ARQUITECTURA

PROGRAMA DE INGENIERIA INDUSTRIAL

2021
 TALLER DE TERMODINÁMICA N°1
TABLAS TERMODINÁMICAS Y ECUACIONES DE ESTADO

1. Complétese la siguiente tabla de propiedades del agua. Indíquese el análisis.

ft 3 Btu
P ,( psia) T (ºF) V ( )
l bm
U(
l bm
) X Descripción de fase

1 11,538 200 30854 1000 0,918 MLV

2 80 312,02 0,1843489688 282 0,00015853 MLV
3 140 700 4,86 1251,4 NA Vapor Sobrecalentado
4 0,099386 100 350 1043,284 NA Vapor Sobrecalentado

Para determinar los valores faltantes en la tabla, se realizará un estudio de las

propiedades que ahí se indican, lo cual permitirá identificar las propiedades en blanco.
Cabe resaltar que, los valores se buscaran en las tablas del sistema Ingles
Para el caso 1, se identifica la descripción de la fase a partir de:
Btu
T =200 ºF y U=1000
l bm

Se busca en la tabla agua saturada. Tabla de temperatura, la temperatura indicada y se

Btu
consideran los valores de U (energía interna en ): U f =168,10 ; U g=1074
l bm

Btu
Teniendo en cuenta, las condiciones dadas (U =1000 ), se tiene que la descripción
l bm
de la fase es Mezcla Liquido-Vapor (MLV). Posteriormente, se calcula la calidad y a
partir de ello, la propiedad de volumen especifico (V) y la presión encontrada de
11,538psia:
Btu Btu
1000 −168,10
U −U f l bm lb m
X = sistema = =0,918
U fg Btu Btu
1074 −168,10
lb m l bm

ft 3 ft 3 ft 3 ft 3
V sistema =V f + X V fg =0,01663
l bm (
+ ( 0,918 ) 33,613
l bm
−0,01663 )
l bm
=30,858
l bm

Para el caso 2, se identifica la descripción de la fase a partir de:

Btu
P=800 psia y U =282
l bm
Se busca en la tabla agua saturada. Tabla de presión, la presión indicada y se consideran
Btu
los valores de U (energía interna en ): U f =281,87 ; U g=1102,3
l bm

Btu
Teniendo en cuenta, las condiciones dadas (U =282 ), se tiene que la descripción de
l bm
la fase es Mezcla Liquido-Vapor (MLV). Posteriormente, se calcula la calidad y a partir
de ello, la propiedad de volumen especifico (V) y la temperatura encontrada de 312,02
Btu Btu
282 −281,87
U −U f l bm l bm
X = sistema = =0,00015853
U fg Btu Btu
1102,3 −281,87
lb m l bm

ft 3 ft 3 ft 3 ft 3
V sistema =V f + X V fg =0,01757
l bm (
+ ( 0,00015853 ) 5,1689
lb m
−0,01757 )
l bm
=0,1843489688
l bm

Para el caso 3, se identifica la descripción de la fase a partir de:

ft 3
P=140 psia y V =4,86
l bm

Se busca en la tabla agua saturada. Tabla de presión, la presión indicada y se consideran

ft 3
los valores de V (volumen especifico en ): V f =0,01802 ; V g =3,2202
lb m

ft 3
Teniendo en cuenta, las condiciones dadas (V =4,86 ), se tiene que la descripción de
l bm
la fase es Vapor Sobrecalentado, en este caso, la propiedad Calidad no aplica.
Posteriormente, se busca en la tabla vapor de agua sobrecalentado, la presión de 140psia
ft 3
y un volumen especifico de 4,86 . de esta manera, se encuentra a T=700ºF y U=1251,4
lb m
Btu
l bm

Para el caso 4, se identifica la descripción de la fase a partir de:

ft 3
T =100 ºF y V =350
l bm

Se busca en la tabla agua saturada. Tabla de temperatura, la temperatura indicada y se

ft 3
consideran los valores de V (volumen especifico en ): V f =0,01613 ; V g =349,83
lb m
 ft 3
Teniendo en cuenta, las condiciones dadas (V =350 ), se tiene que la descripción de
l bm
la fase es Vapor Sobrecalentado, en este caso, la propiedad Calidad no aplica.
Posteriormente, se busca en la tabla vapor de agua sobrecalentado, la temperatura de
ft 3
100ºF para todas las presiones de la tabla, se busca un volumen especifico de 350 .
l bm
Btu
de esta manera, se encuentra a P= y U= . (NOTA) Debido a que los valores no se
l bm
encuentran en la tabla, se procede a realizar una extrapolación para hallar los valores de
presión (P) y energía interna (U) de acuerdo a los datos de saturación

Presión (psia) Temperatura (ºF)

X 100
1 101,69
5 162,18
10 193,16
X=0,099386

Presión (psia) Energía interna (

Btu
¿
l bm
0,099386 Y
1 1043,7
5 1062,7
10 1072
Y=1043,2843

Complétese la siguiente tabla de propiedades del refrigerante 134ª (Tetrafluoroetano).

Indíquese el análisis

m3 Kj
P,¿ T (ºC) V ( )
Kg
U(
Kg
) X Descripción de fase
1 3,379 4 0,05096958952 204,8 0,8427 MLV
2 6 21.55 0,0341 240,86514 0,9941 MLV
3 2.8 20 0,07997 246,13 NA Vapor Sobrecalentado
4 4 68,62 0,06610314 284,75 NA Vapor Sobrecalentado

Para el caso 1, se identifica la descripción de la fase a partir de:

 Kj
T =4 ºC y U =204,8
Kg
Se busca en la tabla Refrigerante 134a saturado. Tabla de temperatura, la temperatura
Kj
indicada y se consideran los valores de U (energía interna en ):
Kg
U f =56,99 ; U g=232,38

Kj
Teniendo en cuenta, las condiciones dadas (U =204,8 ), se tiene que la descripción
Kg
de la fase es Mezcla Liquido-Vapor (MLV). Posteriormente, se calcula la calidad y a
partir de ello, la propiedad de volumen especifico (V) y la presión encontrada de
337,90kPa o 3,379bar:
kPa → ¯¿
337,9 kPa∗¿
Kj Kj
204,8 −56,99
U −U f Kg Kg
X = sistema = =0,8427
U fg Kj Kj
232,38 −56,99
Kg Kg

m3 ( m3 m3 m3
V sistema =0.0007804+ X V fg =
Kg (
+ 0,8427 0.060338
)
Kg
−0.0007804
Kg )
=0,05096958952
Kg

Para el caso 2, se identifica la descripción de la fase a partir de:

m3
P=6 ¯¿ 600 kPa y V =0,0341
Kg
Se busca en la tabla Refrigerante 134a saturado. Tabla de temperatura, la temperatura
m3
indicada y se consideran los valores de V (volumen especifico en ):
Kg
V f =0,0008199 ; V g =0,034295

m3
Teniendo en cuenta, las condiciones dadas (U =0,0341 se tiene que la descripción de
Kg
la fase es Mezcla Liquido-Vapor (MLV). Posteriormente, se calcula la calidad y a partir
de ello, la propiedad de energía interna (V) y la presión encontrada de 600kPa o 6bar:

m3 m3
0,0341 −0,0008199
V −V f Kg Kg
X = sistema = 3
=0,9941
V fg m m3
0,034295 −0,0008199
Kg Kg
 Kj Kj Kj
U sistema =U f + X U fg =81.02
Kg (
+ ( 0,9941 ) 160.81 )
Kg
=240,86514
Kg

Para el caso 3, se identifica la descripción de la fase a partir de:

P=2,8 ¯¿ 280 kPa=0,28 MPa y T =20 ºC
Se busca en la tabla Refrigerante 134a saturado. Tabla de presión, la presión dada y se
considera la temperatura de saturación y la temperatura del sistema. Teniendo que:
T sistema >T saturacion; es decir, 20 ºC>−1,23ºC. Por tanto, la descripción de la fase
corresponde a vapor sobrecalentado. Para este caso, la propiedad Calidad no aplica.
Luego, nos ubicamos en la tabla Refrigerante 134a sobrecalentado, se busca la presión
de 0,28MPa y la temperatura de 20ºC para hallar las propiedades de volumen
especifico(V) y energía interna (U), identificando:

m3 Kj
V =0,07997 y U=246,13
Kg Kg

Para el caso 4, se identifica la descripción de la fase a partir de:

Kj
P=4 ¯¿ 400 kPa=0,4 MPa y U =284,75
Kg
Se busca en la tabla Refrigerante 134a saturado. Tabla de temperatura, la temperatura
Kj
indicada y se consideran los valores de U (energía interna en ):
Kg
U f =63,62 ; U g =235,07

Kj
Teniendo en cuenta, las condiciones dadas (U =284,75 ), se tiene que la descripción
Kg
de la fase es Vapor Sobrecalentado. En este caso, la propiedad Calidad no aplica.
Luego, nos ubicamos en la tabla Refrigerante 134a sobrecalentado, se busca la presión
Kj
de 0,4MPa y buscamos a que temperatura la energía interna es igual a 284,75
Kg
Debido a que, el valor de la energía interna no se encuentra explícitamente en la tabla,
se procede a realizar una interpolación (por medio de la calculadora) como se muestra
Kj
en la siguiente tabla, para la energía interna de 284,75 , se tendrá una temperatura
Kg
(X) de 68,62ºC
Kj
TºC U( ¿
Kg
50 268,9
60 277,32
 X 284,75
70 285,86
80 294,53

Así mismo, se realiza una interpolación para hallar el valor del volumen especifico a la
m3
temperatura encontrada. Teniendo así, un volumen especifico de 0,06610314
Kg

m3
TºC V( ¿
Kg
50 0,061724
60 0,064104
68,62 Y
70 0,066443
80 0,068747

2. Calcule el volumen específico de vapor de agua sobrecalentado a 3.5 MPa (35bar) y

450°C, de acuerdo con:
a. La ecuación de Gas Ideal.
b. La carta de compresibilidad generalizada (método analítico y método gráfico)
c. Las tablas de vapor.
Determine el error que se comete en los dos primeros casos.

a. Para la ecuación de Gas Ideal

Ŕ T
PV =RT →V =
P

3,5 MPa∗ ( 1000 kPa

1 MPa )
=3500 kPa

T =450 ºC +273,15 K =723,15 K

KPa∗m 3 (
V=
( 0,4615
Kg∗K )
723,15 K )
=0,09533
m3
3500 KPa Kg

b. Para la carta de compresibilidad:

Método Grafico:
z=0,961
T 723,15 K P 3,5 MPa
T r= = =1,117 P r= = =0,158
T r 647,1 K Pr 22,06 MPa
 m3
V (volumen especifico)=Z∗V ideal →V =1 ( 0,961 )( 0,09533 ) =0,09161
Kg
El error puede definirse como:
V calculado −V real 0,09196−0,09161
%error= = ∗100=0,4 %
V real 0,09196

Método Analítico:

Z=Z 0 +W Z 1
P ¯¿
W =0,345 Pr = =35 ¿
Pc T 723,15 K
220,55 ¯¿ =0,1586 Tr= = =1,117 ¿
T c 647,1 K

Tr
Pr 1,05 1,10 1,117 1,15 1,20
0,05 0,9855 0,9874 0,9891 0,9904
0,1 0,9707 0,9747 0,9780 0,9808
0,1586 0,9518 0,9587 0,9600071 0,9643 0,9689
0,2 0,9401 0,9485 0,9554 0,9611
0,4 0,8743 0,8930 0,9081 0,9205
0
A partir de la doble interpolación, se tiene que el valor de Z =0,9600071

Tr
Pr 1.05 1,10 1,117 1,15 1,20
0,05 -0,0015 0 0,0011 0,0019
0,1 -0,0029 0,0001 0,0023 0.0039
0,1586 -0,0041 0,00082 0,00136 0,0043 0,0068
0,2 -0,0054 0,0007 0,0052 0,0084
0,4 -0,0092 0,0038 0,0127 0,0190
1
A partir de la doble interpolación, se tiene que el valor de Z =0 , 00136

Z=Z 0 +W Z 1=( 0,9600071 )+ ( 0,345 ) ( 0,00136 )=0,9604763

m3 m3
V ( Volumen especifico )=Z∗V gi =0,9604763∗0,09533 =0,091562
Kg Kg

c. Para las tablas de vapor:

 m3
P=3,5 MPa T =450 ºC V =0,9196
Kg
Ru T c A
P=
V́ 2
∗ 1−
( ∇T 3)∗( V́ +B )− 2
V́

(8,314)(175) 42∗10 4 132,339

P=
(0,916)2
∗ 1−
(
( 39 ) ( 175 ) 3 )
∗( 0 , ,916 +0,5084 ) −
0,9162
=1502 kPa

3. Un cilindro de 5ft3 contiene 50Lb de metano y se encuentran en un sitio donde recibe

directamente el calor solar, el manómetro indica que la presión es de 665 Lb/pulg2
manométricas. ¿Cuál es la temperatura del metano en el cilindro? Resolver por Ec. De
estado de gases ideales y Van der Waals.

V =5 f t 3 →0.141583 m 3
m=50 Lb →22.679851 Kg
Lb
P=665 manometricas → 4584.785 Kpa
pulg2

Kpa∗m 3
Ŕ=0.5182
Kg∗K

PV
Ecuacion de gases ideles → PV =m Ŕ T →T =
m Ŕ

1 m3
V →5 f t ∗ 3
(
35.315 f t 3
=0.141583 m3)
1 Kg
m→ 50 Lb∗( 2.2046 )=22.679851 Kg
2 2 2
Lb 4.448 N 12 pulg 3 ft 1 yd
∗( )∗( ) ∗( ) ∗(
m)
P →665 2
pulg 1 Lb 1 ft 1 yd 0.9144

Lb 4.448 N 144 pul g 2 9 f t2 1 y d2

¿ 665
pulg 2
∗
1 Lb (∗
1 f t2 )( ∗
1 y d2
∗ )( )(
0.83612736 m2 )
N Pa∗1 Kpa
¿ 4584785.17 →=4584785.17 =4584.785 Kpa
m 2
1000 Pa
 PV 4584.785 Kpa∗0.141583 m3
T= = =55.232217 K
m Ŕ Kpa∗m3
22.679851 Kg∗0.5182
Kg∗K

a
Ecuacion de Van Der Waals → P+
( V2 )
( V −b )=RT

P=4584.785 Kpa

V =0.141583 m3

Kpa∗m 3
R=0.5182
Kg∗K
T =?
T c =191,1 K

Pc =4,64 Mpa → 4640 Kpa

2
Kpa∗m3
a=
2
27 R T c 2

=
(
27 0.5182
Kg∗K )
∗(191,1 K )2
=0.8916245
m6 Kpa
64 P c 64(4640 Kpa) Kg2

Kpa∗m3
b=
RTc
=
(
0.5182
Kg∗K )
∗(191,1 K )
=0.0026677
m3
8 Pc 8(4640 Kpa) Kg

m6 Kpa
a=0.8916245
Kg 2

m3
b=0.0026677
Kg

( P+ Va ) (V −b )=RT
2

Kpa∗m3
(
643.07567= 0.5182
Kg∗K
T)
 643.07567
=T
Kpa∗m3
0.5182
Kg∗K
T =1240.97968 K

4. Un tanque de 3.27 m3 contiene 100Kg de nitrógeno a 175K. Determine la presión del

tanque, usando:
a) la ecuación de un gas ideal
b) la ecuación de Van der Waals
c) la ecuación Bridgeman
d) la ecuación de Benedict.
Compare los resultados con un valor real de 1505 KPa.

Para la ecuación de gas ideal:

Volumen especifico del nitrógeno

3,27 m3
V= =0,0327
100 Kg
RT
P=
V

kPa∗m 3
(0,2968 )(175 K)
Kg∗K
P= =1588,37 kPa
m3
0,0327
Kg

La constante de Van der Waals:

27 R 2 T 2cr
a= =(27)¿ ¿
64 Pcr
R T cr (0,2968)(126,2) m3
b= = =0,00138
8 P cr (8)(3390) Kg
P=1495 kPa

Para la constante de Bridgeman:

m3
V́ =M v =( 28,013 ) ( 0,0327 )=0,916
Kmol
 ( av )=136,2315 (1− 0,02617
A=A 0 1−
0,916 )
=132,339

b −0,00691
B=B ( 1− )=0,05046 ( 1−
0,916 )
0 =0,05084
v

4 m3 K 3
C=4,2∗10
Kmol
Para la ecuación de Benedict:
−γ
RuT C 1 bR T −a aα c γ
P=
V́ (
+ B0 RU T −A 0− 0 2 + U 3
T2 V V́ )
+ 6 + 3 2 1+ −2 e V
V V́ T v ( )
−2

( 8,314 )( 175 ) 8,164 x 10 5 1 ( 0,002328 ) ( 8,314 ) ( 175 )−

P=
( 0,916 ) ( (
+ ( 0,04074 ) ( 8,314 ) ( 175 )−( 106,73 ) −
1752 ))( ( 0,0327 )
2
)(
+
( 0,916 )3

327,73
P=1519,82 kPa