Informe No.

1-Ecuaciones de estado cúbicas

Laboratorio-Termodinámica II.

Informe final- Capitulo 1- Ecuaciones de estado cúbicas

Romero, A-1. Narváez, A-2. Hurtado, J-2.
1-200156179
Resumen análisis termodinámico de sustancias
puras.
En este estudio, se empleó Aspen HYSYS
V14® para llevar a cabo una serie de Palabras clave
simulaciones termodinámicas, con el fin
de comprender y evaluar cómo distintas EdE; Ecuación de estado, Costald, PR;
ecuaciones de estado que influyen en la Peng-Robinson, Vander Waals; VDW,
predicción de propiedades de sustancias Rackett, Isobara e Isoterma.
puras bajo varias condiciones.
Ciclobutano y el Isobutano, se examinó el Abstract
comportamiento de estas sustancias
utilizando las ecuaciones de estado de In the present report, various simulations
Rackett, Costald, y Peng-Robinson. Estas were conducted using the software known
ecuaciones, al aplicarse al volumen as Aspen HYSYS V14®. The simulations
específico, exhiben una capacidad were carried out with the purpose of
destacada para predecir con precisión las understanding how the different state
propiedades del líquido saturado, el vapor equations used, reflect on the behavior of
saturado y el vapor sobrecalentado. Los a substance under different conditions.
resultados demostraron que las ecuaciones The substances used during these
de estado son cruciales para predecir con simulations were Cyclobutane and
precisión propiedades termodinámicas y Isobutane. The previously mentioned
que Aspen HYSYS V14® simplifica
substances were evaluated using different
notablemente estos cálculos. Se notaron
state equations. The purpose of comparing
diferencias importantes al utilizar diversas
correlaciones para el volumen de líquido the state equations is to understand how
saturado. Adicionalmente, se ha adquirido the precision and accuracy of state
conocimiento sobre la utilidad de equation varies depending on the type of
representar estas ecuaciones en diagramas substance analyzed, purpose, or
de presión-volumen (P-v) y temperatura- thermodynamic state it is in. The state
volumen (T-v). Estas representaciones equations used were Rackett and Costald.
gráficas proveen una visión más clara y These equations were compared with
comprensible del comportamiento scientific literature to identify any
termodinámico de la sustancia, facilitando differences between them, showing that
la interpretación de resultados y not all state equations are meant to be used
permitiendo una mejor comprensión de los under the same conditions, making them
fenómenos estudiados. En consecuencia,
imprecise when modeling states such as
la combinación de ecuaciones de estado
cúbicas y representaciones gráficas saturated liquid on different occasions.
constituye una metodología robusta para el
Informe No. 1-Ecuaciones de estado cúbicas
Laboratorio-Termodinámica II.

2. Diseño experimental-
metodología
1. Introducción
En esta experiencia de laboratorio se
En esta sección, se expondrán los
utilizó el programa Aspen HYSYS V14®,
antecedentes fundamentales que respaldan
en el cual se modeló el comportamiento de
el experimento llevado a cabo mediante el
las sustancias a diferentes condiciones y a
programa Aspen HYSYS V14®. La
diferentes ecuaciones de estado. Es por
elección de este programa se basa en su
esto que se puede resumir las simulaciones
capacidad para realizar cálculos utilizando
en cuatro distintas partes. Primero la
la ecuación de Peng-Robinson y otras
simulación utilizando la ecuación de
ecuaciones cúbicas, con la incorporación
estado de los gases, segundo utilizando la
de factores de corrección como Rackett y
ecuación PR-Costald, tercero utilizando la
COSTALD, especialmente en relación con
ecuación Rackett y finalmente las isóbaras
el volumen específico del líquido saturado.
e isotermas de cada sustancia. Cuando se
Estos elementos son cruciales para
realiza cada simulación es importante
comprender y predecir el comportamiento
tener en cuenta que para cada sustancia
termodinámico de sustancias en estados
implementada se debe establecer
críticos.
parámetros que cumplan con las
El objetivo principal de esta experiencia de características de estas mismas.
laboratorio es analizar y contrastar los
Durante este proceso, los materiales que se
resultados obtenidos a través de diferentes
utilizaran son; Aspen HYSYS V14®
modelos termodinámicos, centrándose
especialmente en la ecuación de Peng- Para el primer caso, se utilizó la ecuación
Robinson y Rackett Para lograr esto, se de estado Costald a partir de la sustancia
llevó a cabo la determinación de llamada Isobuteno. Cuando se trabaja con
parámetros y la solución para líquido esta ecuación de estado es importante
saturado, vapor saturado e isotermas. Este definir las siguientes condiciones.
enfoque permitirá identificar
discrepancias entre los modelos, lo que a • Fracción de Vapor: 0.0 y
su vez facilitará mejoras en la precisión de 1.0
los cálculos de propiedades • Presión: 1000 kPa y 1000
termodinámicas. kPa
• Flujo másico 1 y 1
• Presión Crítica: 3670 kPa
• Presión del punto Triple:
1.948 × 10−5 kPa
Informe No. 1-Ecuaciones de estado cúbicas
Laboratorio-Termodinámica II.

• Temperatura crítica: Finalmente, solo harían falta las isobaras e

407.7K isotermas. Cuando se realizan estas curvas,
• Temperatura del punto es importante saber que el proceso de
triple: 113.74 K ejecución no es el mismo expuesto
anteriormente. No obstante, el proceso a
Las anteriores especificaciones fueron seguir también se encuentra en [2] (página
dadas para poder completar las curvas de 28). Es importante tener en cuenta que,
liquido saturado y de vapor saturado, cuando se sigue la metodología expuesta
siguiendo el paso a paso mostrado en [2] en [2] es necesario tener en cuenta las
(paginas 25-26). Cabe recalcar que una vez especificaciones de cada sustancia, dado
se complete estas curvas, solo será que estas especificaciones son las que
necesario cambiar la ecuación de estado delimitan las isotermas y las isóbaras.
utilizada. Es por esto que para realizar la Finalmente es importante recordar que
simulación con Rackett se utilizan las cuando se trabaja con isotermas, se debe
mismas especificaciones, cambiando solo escoger una temperatura constante que
la ecuación de estado. este dentro del domo de liquido y vapor
Para el caso 2 analizado, también se saturado y una presión que es la que
empieza analizando la sustancia, la cual es variara. De esta misma manera cuando se
Ciclobutano. La ecuación de estado con la trabaja con una isóbara se deberá escoger
que se empieza a analizar será la Costald. una presión que este dentro del domo de
Las especificaciones que se tendrán en liquido y vapor saturado y esta será
cuenta para esta sustancia son: constante mientras que la temperatura que
se escoja será la que varíe. Para escoger
• Fracción de Vapor: 0.0 y 1.0 estos valores se tienen en cuenta los
• Presión: 1000 kPa y 1000 kPa valores de presiones y temperatura críticas
• Flujo másico 1 y 1 y de punto triple. Dentro de ese intervalo
• Presión Crítica: 4985 kPa debe estar el valor escogido para las
• Presión del punto Triple: 0.1837 isotermas e isobaras.
kPa
Una vez se completen, todos los pasos
• Temperatura crítica: 459.9 K
mencionados anteriormente se habrá
• Temperatura del punto triple: logrado hacer las simulaciones para ambas
182.4 k sustancias de trabajo, dando como
Al igual que para el primer caso, las resultado distintas graficas.
especificaciones que se presentaron Finalmente, como se solicita en la guía,
anteriormente fueron las que se utilizaron todas las gráficas obtenidas bajo las
para lograr obtener las curvas de liquido y
distintas ecuaciones de estado serán
vapor saturado. Tal como se mencionó en
comparas con los resultados expuestos por
el caso 1, para obtener las curvas
la literatura científica, para así poder
utilizando la ecuación Rackett, se utilizará definir cuales son las ventajas y
las mismas especificaciones, solo desventajas de las distintas ecuaciones de
cambiando la ecuación de estado utilizada. estado.
Informe No. 1-Ecuaciones de estado cúbicas
Laboratorio-Termodinámica II.

Resultados
Sustancia: Isobutano
Con los datos obtenidos y usando el
graficador de Excel, se obtuvo las tres Figura 3. Comparación de la curva de
curvas de líquido y vapor saturados para el líquido y vapor saturados T-v para Costald
Isobutano. y EOS

Figura 1. Comparación de la curva de

líquido y vapor saturados P-v para Costald
y EOS. Figura 4. Comparación de la curva de
líquido y vapor saturados T-v para Rackett
y EOS.
Luego se realizaron las isóbaras e
isotermas, para las sustancias. Cabe
recalcar que cuando se realizaron las
isobaras e isotermas, se escogió dos
temperaturas y dos presiones para poder
analizar el comportamiento de estas.

Figura 2. Comparación de la curva de

líquido y vapor saturados P-v para Rackett
y EOS.

Esta misma comparación se realizó con la

temperatura y volumen específico.
Informe No. 1-Ecuaciones de estado cúbicas
Laboratorio-Termodinámica II.

Figura 5. Curva de presiones constantes a

1000 kPa y 2000kPa

Figura 7. Comparación de la curva de

líquido y vapor saturados P-v para Costald
y EOS

Figura 6. Curva de temperaturas

constantes a 300 K y 400 K
Se presentó un desafío con respecto a la
sustancia ciclobutano, ya que la literatura
científica consultada carecía de
información sobre las propiedades
específicas de dicha sustancia. Ante esta
limitación, se llevó a cabo una búsqueda
en otras fuentes científicas con el objetivo Figura 8. Comparación de la curva de
de comprender el comportamiento del líquido y vapor saturados P-v para Rackett
ciclobutano y generar gráficos de curvas y EOS
de saturación. Sin embargo, esta búsqueda
resultó infructuosa, ya que no se logró
obtener información sobre el
comportamiento de la sustancia en
cuestión. Por lo tanto, la única
comparación posible se centrará en las
curvas de líquido saturado y vapor
saturado para las dos configuraciones:
Costald y Rackett.

Figura 9. Comparación de la curva de

líquido y vapor saturados T-v para Costald
y EOS
Informe No. 1-Ecuaciones de estado cúbicas
Laboratorio-Termodinámica II.

Análisis de resultados
La contribución fundamental al desarrollo
de este laboratorio se basó en el empleo del
software Aspen HYSYS V14®. Este
programa facilitó el cálculo preciso de
volúmenes mediante ecuaciones de estado
(EdE) PR, incluyendo correcciones como
Rackett y Costald. Utilizando la
Figura 10. Comparación de la curva de herramienta, se establecieron variables y
líquido y vapor saturados T-v para Rackett parámetros para las sustancias bajo
y EOS estudio, como Isobutano y Ciclobutano.
Estas variables, que abarcan volumen
Una vez se obtienen las curvas de líquido
específico, densidad, entalpía, entropía,
y vapor saturado, se busca hacer las
presión, y puntos críticos de presión y
gráficas de isotermas e isóbaras.
temperatura, se compararon entre las
sustancias y sus respectivos casos
variantes.
Al analizar las gráficas del laboratorio, con
énfasis en el volumen y el comportamiento
de las sustancias con las correcciones
durante la simulación, se destacan
diferencias notables. Para el monóxido de
Isobutano, las gráficas (Figura 1 y Figura
3) indican que la corrección Costald
permite simular de manera más precisa un
Figura 11. Curva de presiones constantes a gas ideal, obteniendo datos más precisos y
1000 kPa y 1500 kPa mejores resultados. En el caso del
ciclobutano (C4H8), como se observa en
las gráficas (Figura 7 y Figura 9), se infiere
que la corrección Costald sigue siendo la
más cercana para simular la sustancia
como un gas ideal.
Luego de obtener tanto los resultados para
el Ciclobutano como para el Isobutano, se
decidió calcular los volúmenes específicos
para estas sustancias teniendo en cuenta
Figura 12. Curva de temperaturas dos temperaturas distintas. Esto se realizó
constantes a 300 K y 400 K utilizando el proceso expuesto en el libro
guía, esto con el propósito de lograr
Informe No. 1-Ecuaciones de estado cúbicas
Laboratorio-Termodinámica II.

comparar la precisión y exactitud del termodinámico de sustancias puras,

programa utilizado. Lastimosamente, los considerando diferentes tipos de
cálculos realizados al ser comparados con ecuaciones (como la del gas ideal y las
los datos utilizados dan un error ecuaciones cúbicas) y la importancia de
demasiado alto, mostrando que los factores de corrección en estos
probablemente hay algún error con la parte cálculos. La EdE de Peng-Robinson se ha
práctica realizada. Estos errores se pueden convertido en la ecuación de estado más
atribuir básicamente a tres aspectos; error popular para sistemas de gas natural en la
humano, equivocaciones al momento de industria petrolera durante la década de
calcular algún dato, error de cálculo, al 1970. el rendimiento del PR y el SRK,
utilizar una calculadora científica para están bastante cerca de un empate; Están
cada proceso realizado, los valores de "codo a codo", excepto por un
respuesta siempre cargan cierto error, que comportamiento ligeramente mejor del PR
al realizar múltiples pasos empieza a ser en el punto crítico. Un rendimiento
apreciable y finalmente el error de ligeramente mejor en condiciones críticas
búsqueda científica, dado que si se utilizó hace que el PR sea algo más adecuado para
algún dato erróneo este creara un valor sistemas de gas/condensado. En el marco
totalmente erróneo. de este experimento de laboratorio, se
Ciclobutano- Volumen específico (m3/kg)
implementó el software Aspen HYSYS
Temperatura V. Experimental V. Analítico Error% V14® para modelar de manera detallada el
Liq.Saturado 0,2125 0,001117 189,241719 comportamiento de diversas sustancias
bajo variadas condiciones y mediante el
100 K Vap.Saturado 0,2158 0,07935 171,959672
uso de distintas ecuaciones de estado. Este
Liq.Saturado 0,3135 0,00127 245,850394
proceso permitió desglosar la simulación
200 K Vap.Saturado 0,4521 0,0744 507,66129 en varias fases secuenciales. Inicialmente,
se efectuó la modelización empleando la
Tabla1. Volumen específico calculado
ecuación de estado de los gases ideales; en
Ciclobutano.
una segunda fase, se procedió con la
Ciclobutano- Volumen específico (m3/kg) aplicación de la ecuación de estado de
Temperatura V. Experimental V. Analítico Error%
Peng-Robinson-Costald; seguidamente, se
Liq.Saturado 0,1568 0,01607 87,57 utilizó la ecuación de estado de Rackett
100 K Vap.Saturado 0,1467 0,0866 69,39953811 para una tercera fase de modelado; y, por
último, se realizaron análisis de las
Liq.Saturado 0,1471 0,0704 108,9488636
isóbaras e isotermas correspondientes a
200 K Vap.Saturado 0,1235 0,7317 83,12149788 cada sustancia estudiada.
Tabla1. Volumen específico calculado
Isobutano.
Conclusión
Cuestionario
El propósito primordial de este
En el capítulo dado se explicó cómo se experimento de laboratorio consistió en
utilizan las ecuaciones de estado para llevar a cabo un análisis y comparación
predecir el comportamiento meticulosos de los datos obtenidos
Informe No. 1-Ecuaciones de estado cúbicas
Laboratorio-Termodinámica II.

mediante la aplicación de distintos precisión para ambos compuestos, según

modelos termodinámicos. Se pone un lo evidenciado por las gráficas de
énfasis particular en la evaluación de la comportamiento volumétrico y de estado
ecuación de estado de Peng-Robinson y el de las sustancias.
modelo de Rackett, con el fin de
diferenciar su precisión Esta experiencia Referencias
de laboratorio permitió entender de
manera más profunda cómo estos modelos [1] Constantes y coeficientes de
pueden predecir el comportamiento de propiedades físicas.
diferentes sustancias bajo variadas https://termoapunefm.files.wordpress.com
condiciones de presión y temperatura y se /2011/10/41-1-
representó el comportamiento
termodinámico de una sustancia de manera
más clara y comprensible mediante las
[2] González A., Acuña F., Yepes Y., Bula
ecuaciones relacionadas mediante gráficos
A. (2022). Casos de estudio prácticos en
presión-volumen (P-v) y temperatura-
termodinámica.
volumen (T-v). Estos diagramas nos
permiten apreciar de forma gráfica cómo [3] Isobaric Properties for Ciclobutane.
cambian la presión y la temperatura de una (n.d.).
sustancia en función de su volumen, https://webbook.nist.gov/cgi/cbook.cgi?I
facilitando así la interpretación y el D=C75285&Units=SI&Mask=4#Thermo-
análisis de sus propiedades Phase
termodinámicas. A su vez en los resultados
notaron diferencias importantes al usar [4] Isobaric Properties for Isobutane.
diversas correlaciones para el volumen de (n.d.).
líquido saturado en cuanto al modelo https://webbook.nist.gov/cgi/cbook.cgi?I
utilizado. Por otra parte, y no menos D=C75285&Units=SI&Mask=4#Thermo-
importante se destacó el uso del software Phase
Aspen HYSYS para el avance en el
laboratorio, permitiendo la estimación [5] Y. A. Cengel, M. A. Boles, and M.
exacta de volúmenes con la aplicación de Kanoglu, Termodinámica, 9 ed. Ciudad de
la ecuación de estado Peng-Robinson México: McGraw-Hill, 2019
(PR), y mejoras mediante las correcciones
de Rackett y Costald. Se definieron y Anexos
analizaron variables críticas como
volumen específico, densidad, entalpía, Ciclobutano
entropía, presión, y puntos críticos para
T = 100K
Isobutano y Ciclobutano. Las
comparaciones de estas variables, bajo P = 1000 kPa
diferentes escenarios, revelaron
diferencias significativas en el W = 0.2575
comportamiento de las sustancias,
M = 56.107 g/mol
particularmente en términos de simulación
de gas ideal. Los resultados mostraron que Pc = 3499 kPa
la corrección de Costald es superior en
Tc = 460 K
Informe No. 1-Ecuaciones de estado cúbicas
Laboratorio-Termodinámica II.

Ru = 8.1416 kJ/mol K 19.32 ∗ 0.1481 ∗ 100 𝐾 𝑚3

𝑣= = 0.2158
1000 𝑘𝑃𝑎 𝑘𝑔
𝑘𝐽
𝑅𝑢 8.1416𝑚𝑜𝑙 𝐾 𝑘𝐽
𝑅= = 𝑔 = 0.1451 Para 200k
𝑀 56.108𝑚𝑜𝑙 𝑘𝑔∗𝐾
T = 200K
𝑚 = −0.2699 (0.2575)2 +
1.54226 (0.2575) + 0.37464 = 0.7538 P = 1000 kPa
100 𝐾 W = 0.2575
𝑇𝑟 = = 0.6522
460 𝐾
M = 56.107 g/mol
𝛼 = ( 1 + 0.4480 ∗ ( 1 − (0.6522)0.5 ))2
= 1.18 Pc = 3499 kPa
𝑘𝐽 Tc = 460 K
0.457235 ∗ (0.1481 )2 ∗ ( 460)2
𝑘𝑔 ∗ 𝐾
𝑎=
3499 𝑘𝑃𝑎 Ru = 8.1416 kJ/mol K
𝑚3
= 0.00155 ( )2 𝑘𝐽
𝑘𝑔 𝑅𝑢 8.1416𝑚𝑜𝑙 𝐾 𝑘𝐽
𝑅= = 𝑔 = 0.1451
𝑀 56.108𝑚𝑜𝑙 𝑘𝑔∗𝐾
𝑘𝐽
0.077796 ∗ ( 0.1481 ) ∗ (460 𝐾)2
𝑘𝑔 ∗ 𝐾
𝑏= 𝑚 = −0.2699 (0.2575)2 +
3499 𝑘𝑃𝑎
𝑚3 1.54226 (0.2575) + 0.37464 = 0.7538
= 0.6064 ( )2
𝑘𝑔
200 𝐾
𝑎𝑃 0.6064 ∗ 1000 𝑘𝑃𝑎
𝑇𝑟 = = 0.6522
460 𝐾
𝐴= 2
= = 0.007066
(𝑅𝑇) (0.1481 ∗ 100 𝐾)2
𝛼 = ( 1 + 0.4480 ∗ ( 1 − (0.6522)0.5 ))2
= 1.18
𝑏𝑃 0.6064 ∗ 1000 𝑘𝑃𝑎 𝑘𝐽
𝐵= = = 40.945 0.457235 ∗ (0.1481 )2 ∗ ( 460)2
(𝑅𝑇) (0.1481 ∗ 100 𝐾) 𝑘𝑔 ∗ 𝐾
𝑎=
3499 𝑘𝑃𝑎
𝑉1 𝑚3
= 0.00155 ( )2
3(40.945 − 2(0.007066) − 3(0.0007066) − (−(1 − 0.007066)2 𝑘𝑔
=
9
𝑘𝐽
= 13.97 0.077796 ∗ ( 0.1481 ) ∗ (460 𝐾)2
𝑘𝑔 ∗ 𝐾
𝑏=
𝑉2 = −2.6 3499 𝑘𝑃𝑎
𝑚3
= 0.6064 ( )2
𝑘𝑔
Luego, se encontraron las raíces reales, esto se hizo 𝑎𝑃 0.6064 ∗ 1000 𝑘𝑃𝑎
determinando Z. Al hallar Z, se utiliza el menor 𝐴= 2
= = 0.001766
(𝑅𝑇) (0.1481 ∗ 200 𝐾)2
valor para el líquido Saturado y el valor mayor para
el vapor saturado.

Los valores utilizados son; 𝑏𝑃 0.6064 ∗ 1000 𝑘𝑃𝑎

𝐵= = = 20.47
(𝑅𝑇) (0.1481 ∗ 200 𝐾)
V= 18.4 - 19.32
18.4 ∗ 0.1481 ∗ 100 𝐾 𝑚3
𝑣= = 0.2125
1000 𝑘𝑃𝑎 𝑘𝑔
Informe No. 1-Ecuaciones de estado cúbicas
Laboratorio-Termodinámica II.

𝑉1 𝑘𝐽
0.457235 ∗ (0.1481 )2 ∗ ( 460)2
3(20.47 − 2(0.001766) − 3(0.001766) − (−(1 − 0.001766) 2 𝑘𝑔 ∗ 𝐾
= 𝑎=
9 3700𝑘𝑃𝑎
= 7.15 𝑚3
= 0.00155 ( )2
𝑘𝑔
𝑉2 = −1.2
𝑘𝐽
0.077796 ∗ ( 0.1481 ) ∗ (460 𝐾)2
Luego, se encontraron las raíces reales, esto se hizo 𝑘𝑔 ∗ 𝐾
𝑏=
determinando Z. Al hallar Z, se utiliza el menor 3700 𝑘𝑃𝑎
valor para el líquido Saturado y el valor mayor para 𝑚3
= 0.6064 ( )2
el vapor saturado. 𝑘𝑔

Los valores utilizados son; 𝑎𝑃 0.6064 ∗ 1000 𝑘𝑃𝑎

𝐴= 2
= = 0.001766
(𝑅𝑇) (0.1481 ∗ 200 𝐾)2
V= 6.75 – 7.21
7.21 ∗ 0.1481 ∗ 200 𝐾 𝑚3
𝑣= = 0.3135
1000 𝑘𝑃𝑎 𝑘𝑔 𝑏𝑃 0.6064 ∗ 1000 𝑘𝑃𝑎
𝐵= = = 20.47
(𝑅𝑇) (0.1481 ∗ 200 𝐾)
6.75 ∗ 0.1481 ∗ 200 𝐾 𝑚3
𝑣= = 0.4521
1000 𝑘𝑃𝑎 𝑘𝑔 𝑉1
3(20.47 − 2(0.001766) − 3(0.001766) − (−(1 − 0.001766)2
=
9
= 7.15
Para el Isobutano
𝑉2 = −1.2
T = 100K
Luego, se encontraron las raíces reales, esto se hizo
P = 1000 kPa determinando Z. Al hallar Z, se utiliza el menor
valor para el líquido Saturado y el valor mayor para
W = 0.1825 el vapor saturado.
M = 58.12 g/mol Los valores utilizados son;

Pc =3700 kPa V= 6.75 – 7.21

Tc = 407.7 K 7.21 ∗ 0.1481 ∗ 200 𝐾 𝑚3

𝑣= = 0.1568
1000 𝑘𝑃𝑎 𝑘𝑔
Ru = 8.1416 kJ/mol K
6.75 ∗ 0.1481 ∗ 200 𝐾 𝑚3
𝑘𝐽
𝑣= = 0.1467
𝑅𝑢 8.1416𝑚𝑜𝑙 𝐾 𝑘𝐽 1000 𝑘𝑃𝑎 𝑘𝑔
𝑅= = 𝑔 = 0.1451
𝑀 56.108𝑚𝑜𝑙 𝑘𝑔∗𝐾 Para 200k

𝑚 = −0.2699 (0.1825)2 + T = 200K

1.54226 (0.1825) + 0.37464 = 0.7538
P = 1000 kPa
200 𝐾
𝑇𝑟 = = 0.6522 W = 0.1825
460 𝐾

𝛼 = ( 1 + 0.4480 ∗ ( 1 − (0.6522)0.5 ))2 M = 58.12 g/mol

= 1.18
Pc =3700 kPa
Tc = 407.7 K
Informe No. 1-Ecuaciones de estado cúbicas
Laboratorio-Termodinámica II.

Ru = 8.1416 kJ/mol K
𝑘𝐽
𝑅𝑢 8.1416𝑚𝑜𝑙 𝐾 𝑘𝐽
𝑅= = 𝑔 = 0.1451
𝑀 56.108𝑚𝑜𝑙 𝑘𝑔∗𝐾

𝑚 = −0.2699 (0.1825)2 +
1.54226 (0.1825) + 0.37464 = 0.7538
200 𝐾
𝑇𝑟 = = 0.6522
460 𝐾

𝛼 = ( 1 + 0.4480 ∗ ( 1 − (0.6522)0.5 ))2

= 1.18
𝑘𝐽
0.457235 ∗ (0.1481 )2 ∗ ( 460)2
𝑘𝑔 ∗ 𝐾
𝑎=
3700𝑘𝑃𝑎
𝑚3
= 0.00155 ( )2
𝑘𝑔
𝑘𝐽
0.077796 ∗ ( 0.1481 ) ∗ (460 𝐾)2
𝑘𝑔 ∗ 𝐾
𝑏=
3700 𝑘𝑃𝑎
𝑚3
= 0.6064 ( )2
𝑘𝑔
𝑎𝑃 0.6064 ∗ 1000 𝑘𝑃𝑎
𝐴= 2
= = 0.001766
(𝑅𝑇) (0.1481 ∗ 200 𝐾)2

𝑏𝑃 0.6064 ∗ 1000 𝑘𝑃𝑎

𝐵= = = 20.47
(𝑅𝑇) (0.1481 ∗ 200 𝐾)

𝑉1
3(20.47 − 2(0.001766) − 3(0.001766) − (−(1 − 0.001766)2
=
9
= 7.15

𝑉2 = −1.2

Luego, se encontraron las raíces reales, esto se hizo

determinando Z. Al hallar Z, se utiliza el menor
valor para el líquido Saturado y el valor mayor para
el vapor saturado.

Los valores utilizados son;

V= 6.75 – 7.21
7.21 ∗ 0.1481 ∗ 200 𝐾 𝑚3
𝑣= = 0.1471
1000 𝑘𝑃𝑎 𝑘𝑔
6.75 ∗ 0.1481 ∗ 200 𝐾 𝑚3
𝑣= = 0.1235
1000 𝑘𝑃𝑎 𝑘𝑔