ECUACIONES DE ESTADO CÚBICAS

Carolay J Cabrales Castro; Anderson J Romero Prieto; Miguel Ángel Terán

Romero

Estudiante de Ingeniería Mecánica, Código: 200127502

Estudiante de Ingeniería Mecánica, Código: 200127234

Estudiante de Ingeniería Mecánica, Código:

Profesor: Johnnys A Bustillo Maury

Laboratorio de Termodinámica II Universidad del Norte

Fecha 10-03-2020

Barranquilla-Colombia

__________________________________________________________________

Resumen

El propósito de este laboratorio fue interpretar los diferentes fenómenos

termodinámicos modelados por el programa aspen HYSYS. Mediante la aplicación
de las Ecuaciones de Estado Cúbicas se realizaron cálculos de las propiedades
del Etileno, de esta manera se obtuvieron propiedades como la entalpia, entropía,
presión, temperatura, densidad, volumen específico y relaciones entre ellos. Datos
que posteriormente fueron comparados con los teóricos brindados por la literatura.
Se realizó la tabulación y se graficaron los resultados y las respectivas
comparaciones para la fase vapor y la fase líquida del compuesto puro,
empleando programas computacionales como Excel y AspenHysys V10. La
ecuación empleada para el cálculo fue Peng-Robinson (PR), haciendo pruebas con
dos casos específicos, el caso Costald y el caso Rackett, que sólo presentaron diferencias
significativas en la densidad y por consiguiente en el volumen específico del Etileno.

Palabras claves: Etileno, Domo, Liquido saturado, Vapor saturado, Isóbara, Isoterma,
Peng-Robinson, propiedades.

Abstract:

By the application of the Cubic State Equations they were made the calculations of the
properties for the Ethylene, in this way properties such as enthalpy, entropy, pressure,
temperature, density, specific volume and relationships between them were obtained.
Later they were compared with the teoric information obtained from the literature. It were
made the tabulation and plotted these results and the appropriate comparisons for the
vapor and liquid phases of the pure compound using some specialized computational
programs like Excel y AspenHysys V10. The cubic state equation used for the
calculations was Peng-Robinson (PR), testing two specific cases, the Costald case and
the Rackett case, which had only significant differences in density and therefore the
specific volume of Ethylene.
Keywords: Ethylene, Dome, Saturated liquid, Saturated steam, Isobar, Isotherm, Peng-
Robinson.
________________________________________________________________________
1. Introducción diferentes valores y/o grados de
precisión.
Los procesos termodinámicos son
más comunes de lo que puede llegar 2. Metodología o Detalles
a pensarse, estos conceptos se Experimentales
encuentran en la cotidianidad, un
mejor entendimiento de dichos
2.1 Equipos
fenómenos, nos permitirán una mayor
comprensión de los eventos -Programa Aspen HYSYS
pequeños, los cuales sirven de base -Excel
para procesos mucho más grandes y 2.2 Procedimiento
significativos. 1. Iniciar Aspen HYSYS
Con la ayuda de las EdE o 2. Por medio del programa
Ecuaciones de Estado se presentó en se modelo un sistema
forma de diagramas la relación entre
presión, temperatura y volumen termodinámico para el
específico. A pesar de que este Etileno puro (C2H4) con
proceso se puede dar en diferentes las siguientes
tipos de gases, el presente informe se características las cuales
enfocó en el Etileno (C2H4), a partir fueron ofrecidas por la
del cual se realizaron todas las guía: 1000 kPa, 300 K, 1
simulaciones y análisis presentados.
Kg/s, fracción de vapor
Las EdE permiten analizar el
comportamiento de una sustancia en (0). Además se tomaron
sus fases vapor y líquido. Para los datos del punto crítico
resolver el problema se empleó una que daba HYSYS y los
de las EEC, la cual es PR (Peng- datos del punto triple de la
Robinson) con la variable de estados web, con el fin de variar la
de densidad de Costald y Rackett; presión y la temperatura
Con esta ecuación se calcula el
con respecto al volumen
volumen específico de la sustancia en
cada estado termodinámico definido específico en el estado de
(Temperatura, Presión) mediante densidad Costald.
métodos iterativos, habiendo ya 3. Se repitió el paso anterior
determinado las condiciones mínimas. con fracción de vapor (1) y
Al mismo tiempo, se hace un análisis la misma variación de
detallado del domo del Etileno, presión entre el punto
especificando el comportamiento del
crítico y punto triple.
mismo y de cierta isobara y cierta
isoterma. Esta ecuación cuenta con
 4. Se halló las isobaras e
300Gráfica 2. Diagrama T-v (Colstad)
isotermas por medio de
Aspen HYSYS. 250
5. Con la ayuda de Excel se
graficaron los diagramas 200
P-v, T-v, h-v, s-v. 150
6. Se repitieron los pasos del
2 al 5 con el estado de 100
0 0.1 1 10 100
densidad Rackett.
3. Resultados Grafica 2. Domo T-v para el
Una vez seguido el procedimiento estado de densidad Costald
se obtuvieron las siguientes
gráficas: En esta gráfica se observa lo
mismo que la anterior, que la
Para Costald: línea de líquido saturado no es
muy precisa hallada desde el
estado de densidad Costald;
6000 Gráfica 1. Diagrama P-v (Colstad) Además en este diagrama se
5000 observa una isóbara
4000 perfectamente alineada.
3000
2000
1000 Gráfica 3. Diagrama s-v
0 (Colstad)
0 0.01 0.1 1 10 100 1000 10000
8
Grafica 1. Domo P-v para el
estado de densidad Costald 6
4
En esta gráfica podemos 2
observar que la las gráficas de 0
líquido y vapor saturado no 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14
cierran completamente debido a Grafica 3. Diagrama entropía-
que la línea de líquido saturado volumen específico para el
no es muy precisa hallada desde estado de densidad Costald
el estado de densidad Costald; y
que la isoterma mostrada se Este grafico fue hallado con los
encuentra abajo debido a que se datos obtenidos para la isóbara
escogió una baja temperatura. de costald.
 valor de la presión donde iniciaba
Gráfica 4. Diagrama h-v vapor saturado. De la misma
(Colstad) manera, el paso de la iteración
para vapor saturado también era
2500
muy grande es decir (2). Esto nos
2000
permite identificar que los datos
1500 para la línea de líquido saturado
1000 son más precisas. Por otro lado,
500 la isoterma mostrada se
0 encuentra abajo debido a que se
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 escogió una baja temperatura.
Grafica 4. Diagrama Entalpía-
volumen específico para el
estado de densidad Costald
400 Gráfica 6. Diagrama T-v (Racket)
Este grafico al igual que al 350
anterior, fue hallado con los datos 300
obtenidos para la isóbara de 250
costald. 200
150
100
Para Racket: 0 0.1 1 10 100

Grafica 6. Domo T-v para el

Gráfica 5. Diagrama P-v (Racket) estado de densidad Racket.

En esta gráfica se observa lo

mismo que la anterior, que las
líneas de saturación no se
encuentran debido a los grandes
intervalos con los que se realizó
0
0 0.1 1 10 100 la iteración; Además en este
diagrama se observa una isóbara
Grafica 5. Domo P-v para el perfectamente alineada.
estado de densidad Racket.

En esta gráfica podemos

observar que la las gráficas de Gráfica 7. Gráfica h-v (Racket)
líquido y vapor saturado no 3000
cierran completamente debido a
que cuando el programa realizo 2000
las iteraciones, se inició con una 1000
presión de 0,12 para liquido
saturado, lo cual hizo que la 0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14
última iteración no llegara hasta el
 Grafica 7. Diagrama entalpía - este orden de ideas, al momento
volumen específico para el de hallar datos con la ecuación PR
estado de densidad Racket. se recomienda utilizar los estados
por separado. Haciendo uso del
que tiene mayor precisión para
vapor saturado y de igual manera
Gráfica 8. Gráfica s-v (Racket) para liquido saturado. (Costald y
8 Rackett respectivamente)

6
4
2
0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 5. Solución al Cuestionario
Grafica 8. Diagrama entropía-
volumen específico para el
estado de densidad Racket. a) De acuerdo con el diagrama de la
Figura 1, ¿es posible generar vapor de
agua a 263 K? ¿es posible observar la
4. Análisis de Resultados
fusión de hielo seco a presión
A partir de los resultados
obtenidos a través del programa atmosférica?
Aspen HYSYS durante toda la
práctica, se puede afirmar que la
ecuación de Peng-Robinson tiene
mucha similitud a los valores
reales que se hallaron por
métodos experimentales y se
encuentran tabulados en este
informe. También se puede
realizar una comparación entre los
diferentes estados de densidades
Costald y Rackett, es importante
resaltar y considerar que, Costald
R// Respecto a la isoterma trazada en la
cuenta con una mayor precisión al
gráfica, podemos afirmar que si es
momento de hallar valores que se
encuentran en la línea de vapor posible generar vapor de agua a 263
saturado. Por el contrario Rackett K a muy bajas presiones.
es más preciso en los valores No es posible observar la fusión del
cercanos a la línea de líquido hielo seco a temperatura ambiente
saturado. De esto podemos debido a que esta presión está por
concluir que ambos son métodos debajo del punto triple del CO2, el
opuestos y a la vez cambio de estado que se puede dar
complementarios. Cada uno tiene es sublimación.
una función y uso específico. En
 27 R 2 T e2
b) ¿Bajo qué condiciones los gases reales a=
64 P e
se pueden modelar satisfactoriamente
con la ecuación de gas ideal?, ¿cuál es kJ 2
el factor de compresibilidad de CO2 a
a=
(
27∗ 0.5182
Kg∗K )
∗(190.6)2
7000 kPa y 380 K? 64∗(4599 Kpa)
R// Los gases reales se pueden modelar con
la ecuación de gas ideal cuando el
factor de comprensibilidad (Z) muy m3
a=0.8948
cercano a 1, esta condición se cumple Kg
cuando las presiones reducidas (Pr) son Para b:
muy bajas (Pr << 1) R2T e
b=
independientemente de la 8 Pe
Temperatura; también se puede
presentar el mismo comportamiento si kJ 2
(0.5182 ) (190.6)2
la temperatura reducida (Tr) es mayor Kg∗K
b=
a 2, independientemente de la presión; 8∗(4599 Kpa)
a excepción si Pr >> 1).

380 K m3
Tr= =1.24 ≅ 1.20 b=1.39x10−3
304.2 K Kg

7 Mpa
Pr= =0.94
7.39 Mpa d) ¿Cuáles son los parámetros de la EdE
PR para CO2 y H2O a 1000 kPa y 650 K?

Con estos datos, si los usamos con la Para el CO2:

carta de compresibilidad, podemos
decir que: Para determinar el parámetro a se
Z= 0,81 necesita conocer α , y para conocer α se
necesita conocer m y Tr. M depende del
c) ¿Cuáles son las constantes a y b de la
EdE VDW para una sustancia cuya factor acéntrico de Pitzer, w, que difiere y
masa molar, T c y Pc son está establecida en cada gas. En el caso
respectivamente 16,043 g mol-1, 190,6 del CO2, estos son los datos de w y Tc.
K y Pc 4599 kPa? ¿cuál podría ser
dicha sustancia? ω = 0.228 ; Tc = 304.2K
R// Observando la masa molar y
comparándola con la tabla A1, se debido a que ω ≤ 0,49m se determina
determina que la sustancia en cuestión mediante la siguiente fórmula.
es el metano CH4.
Para a: m=−0,26992ω 2+1,54226 ω+0,37464
 Se aplican las mismas ecuaciones que
m=−0,26992¿ ( 0.228 )2 +1,54226∗(0.228)+0,37464
en el caso del CO2, pero se utilizan los
m=0.712 siguientes valores de w y Tr:

Luego de tener m, procedemos a hallar ω = 0.3436, Tc =647,1K

la temperatura reducida
debido a que ω ≤ 0,49m se determina
mediante la siguiente fórmula.
T
T r=
Tc m=−0,26992ω 2+1,54226 ω+0,37464
650 K 2
m=−0,26992∗( 0.3436 ) +1,54226∗( 0.3436 ) +0,37464
T r= =2.13675
304.2 K
m=0.87269
Teniendo los valores necesarios,
procedemos a hallar α T
T r=
Tc
α 0,5=1+m(1−T 0,5
r )
650 K
0,5 0.5 T r= =1.0045
α =1+0.712∗(1−2.13675 ) 647.1 K

α =0.450

Ahora que tenemos α , encontramos el α 0,5=1+m(1−T 0,5

r )

primero parámetro a:
α 0,5=1+ ( 0.87269 )∗( 1−1.0045 0.5 )
0,457235 R 2 T 2c
a= α α =0.999
Pc
0,457235 R 2 T 2c
0,457235∗( 0.1889 ) ∗304.2 2 2 a= α
a= 0.450 Pc
7390
0,457235∗( 0.4615 )2 ¿ ( 647.1 )2
a=0.09193 a= ( 0.999 )
( 22060 )
Para encontrar el parámetro b se utiliza
a=1.8466
la siguiente ecuación.
0,077796 R T c
0,077796 R T c b=
b= Pc
Pc
0,077796∗( 0.4615 )∗( 647.1)
0,077796∗( 0.1889 )∗( 304.2) b=
b= =6.049 x 10
−4 (22060)
(7390)
b=1.05316 x 10−3
Para el H2O
 e) ¿Cuáles son los volúmenes de Líquido ( 1.05316 x 10−3 )∗(1000)
saturado para CO2 y H2O a 1000 kPa de B=
( 0.4615 )∗(452.87)
acuerdo con la EdE PR?
B=5.03 x 10−3
Para el H2O:
Ahora reemplazamos A y B en la
Primero hallamos los parámetros a y b ecuación de Peng-Robinson para hallar
z.
m=−0,26992ω 2+1,54226 ω+0,37464
Z3 −( 1−B ) Z 2 + ( A−2 B−3 B ) Z−( AB−B2 −B 3) =0
2
m=−0,26992 ( 0.3436 ) +1,54226 ( 0.3436 ) +0,37464
Z3 −( 1−5.03 x 10−3 )∗Z 2+ ( 0.0549−2∗( 5.03 x 10−3 )−3∗( 5.03
m=0.87269
Z=0.964

Teniendo z, reemplazamos en la fórmula

452.87 K
T r= =0.7 ZRT
647.1 K v=
P

kJ
0,5
α =1+m(1−T ) 0,5
r
v=
(
( 0.964 )∗ 0.4615
kgK )
∗(452.87 K)

1000 kPa
α 0,5=1+0.87269∗(1−0.70.5 )
m3
α =1.3 v=0.2014
kg
0,457235 R 2 T 2c
a= α
Pc
f) ¿Cuáles son los volúmenes de Líquido
2 2
0,457235¿ ( 0.4615 ) ¿ ( 647.1 ) saturado para CO2 y H2O a 1000 kPa de
a= ( 1.3 )
( 22060 ) acuerdo con las correlaciones de
Rackett y Yamada-Gunn?
a=2.4
 Para el CO2
Con estos parámetros hallamos A y B  Por Rackett

aP 1000
A= Pr = =0,1353
( RT )2 7390

( 2.4 )∗(1000) 650

T r= =2,136
A= 304,2
¿¿¿
A=0.0549 Zc=0,96

bP
B=
RT
 R T c (1 +(1−T ) v Lsat =1.798 x 10−6
2
7
Lsat )
v = Z r

Pc c
 Por Yamada-Gunn
8,31∗( 304,2 )
2
7

v Lsat = ¿ ( 0,96 )1+( 1−2,136) =0,31 453,2

7390 T 'r= =0,7003
647,1
 Por Yamada-Gunn 2

T 'r=
233,9
=0,7689
' 3 '
v exp=0,0009 m /kg v scr =v exp ⁡ −( 1−T [ ' 7
r) ln ⁡( 0.29056−0.0
304,2 2

2
[ 7
v scr=¿ −( 1−0,7003 ) ∗ln ⁡(0.29056−0.08775∗(0,3436)) ]
'
v scr=v exp ⁡ −( 1−T [ ' 7
r ) ln ⁡(0.29056−0.08775 ω) ] v scr=0,0008
2
[
v scr=¿ −( 1−0,7689 ) 7∗ln ⁡(0.29056−0.08775∗(0,228)) 1−T
v Lsat =v scr ( 0,29056−0,08775 ω )( ) ] r
2
7

v scr=0,00077 Lsat
2
(1−0,7) 7
v =0,0008∗( 0,29056−0,08775∗(0,3436) )
2
7
1−T r )
v Lsat =v scr ( 0,29056−0,08775 ω )( v Lsat =3,08∗10−4
2

Lsat ( 1−2,136 ) 7a)

¿Cuál es la correlación para la
v =0,00077 ¿ ( 0,29056−0,08775∗(0,228) )
corrección de volumen de líquido
v Lsat =1,98∗10−4 saturado que está implementada por
defecto en Aspen HYSYS® V10?
R// Por defecto es Costald.
 Para el H2O b) ¿Cuál era el vínculo profesional entre
 Por Rackett los autores de la EdE SRK?
R// Esta ecuación fue desarrollada por Redlich
Pv y Kwong mientras trabajaban para la
Zc=
RT Shell Development Company en
3 Emeryville, California.
m
( 1000 kPa )∗(0.001127 ) Kwong comenzó a trabajar para esa
kg
Zc= compañía en 1944, donde conoció a
kJ
(
0.4615
kgK
∗(452.88 K )) Otto Redlich cuando él se unió al grupo
en 1945. La ecuación surgió de su
trabajo ahí. Ellos buscaban una manera
Zc=5.39 x 10−3
algebraica sencilla para relacionar
R T c (1 +(1−T ) ) 2
7 presiones, volúmenes, y temperaturas
v Lsat = Zc r

Pc de los gases con los que ellos

trabajaban, estos gases eran en su
kJ mayoría no polares. Por ello es

v Lsat =
(
0.4615
kgK
∗647.1 K) ¿ 5.39 x 10−3(1+ ( 1−0.7 ) )
2
importante resaltar que esta ecuación
7

22060 es menos precisa para gases con

 enlaces de hidrógeno. Fue presentada Lopez-Echeverry, J. S., Reif-
conjuntamente en Portland, Oregón en Acherman, S.,
el Symposium on Thermodynamics and & Araujo-Lopez, E. (2017). Peng-
Molecular Structure of Solutions en Robinson equation of state: 40 years
1948, como parte de la 14th junta de la through cubics.
sociedad química americana (American
Aspentech. (2017). Aspen HYSYS
Chemical Society). Reif-Acherman,
V10.
Simón (2008).

Yamada, T., & Gunn, R. D. (1973).

6. Conclusiones y Recomendaciones Saturated liquid molar volumes.
En este informe se realizó una Rackett equation. Journal of
comparación entre los resultados Chemical and Engineering Data,
teóricos y experimentales obtenidos 18(2), 234–236.
y analizados por el grupo mediante
el uso del software Aspen Hysys. Rackett, H. G. (1970). Equation of
Para la obtención de estos datos se state for saturated liquids. Journal of
utilizó el Etileno (C2H4). Se utilizó la Chemical and Engineering Data,
ecuación de Peng-Robinson y se 15(4), 514–517.
agregó la variante de estados de Peng, D.-Y., & Robinson, D. B.
densidad de Costald y Rackett. (1976). A new two-constant equation
Basados en los resultados de la of state. Industrial & Engineering
simulación, se hizo una detallada Chemistry Fundamentals, 15(1), 59–
explicación de cada variante y se 64.
pudo llegar a la conclusión que 8. Anexos
Costald presenta mayor exactitud
para valores cercanos a la línea de
vapor saturado, mientras que
Rackett, presenta mayor exactitud
para valores cercanos a la línea de
liquido saturado. De ahí que esta
variante nos llevó a conocer que
ambos son complementarios y
contrarios. Este laboratorio nos
mostró la importancia de realizar un
buen análisis, especifico y detallado
para los casos de estudio que se
llevaron a cabo, asimismo,
comprendimos como se pueden
afectar los diferentes métodos y
variaciones usados.

7. Referencias bibliográficas