LABORATORIO # 4 - TERMODINAMICA QUIMICA II
Laboratorio # 4: Aplicación del cálculo de propiedades en el análisis
de procesos Termodinámicos y determinación de coeficientes de
actividad en mezclas binarias”
Temas a desarrollar
Cálculo de propiedades termodinámicas en un proceso de
estrangulamiento
Cálculo de coeficientes de actividad.
Base teórica
Los procesos termodinámicos, en su diseño y operación requieren el cálculo de
propiedades termodinámicas, el concepto de propiedad residual, así como el
cálculo de propiedades termodinámicas para gas ideal resultan ser
fundamentales para modelar el comportamiento real de dichos procesos.
La utilización de ecuaciones de estado, en este caso RK, siguen prestando
utilidad para la determinación de las propiedades residuales, así como para la
obtención del valor del factor Z tanto para liquido como para vapor saturado.
Base práctica
Ejemplo 1
Gas propano a 20 bar y 400 K se estrangula en un proceso de flujo en estado
estacionario hasta 1 bar. Estime la temperatura final del propano y su cambio
en la entropía. Las propiedades del propano se obtienen a partir de
correlaciones generalizadas apropiadas.
Desarrollo: En el proceso de estrangulación pueden considerarse las
condiciones 1 a la entrada y las condiciones 2 a la salida del proceso.
Basados en la primera ley de la termodinámica:
∆𝐸𝐶 + ∆𝐸𝑃 + ∆𝐻 = 𝑄 + 𝑊
LABORATORIO # 4 - TERMODINAMICA QUIMICA II
Asumiendo que las pérdidas por energía potencial y cinéticas son
despreciables (ΔEC=ΔEP≅0) y que se trata de un proceso adiabático, donde
no se produce trabajo de eje, se tiene:
∆𝐻 = 0
Por lo que el cambio de entalpia vendría dado solamente por:
Y el cambio de entropía por:
Las variables conocidas son T1, P1 y P2; además de los valores críticos para el
propano y la ecuación para determinar su valor de Cp, tomados del DATA
BANK:
Formula del Cp:
Cp=CPVAP A +(CPVAP B)T + (CPVAP C)T2 +(CPVAP D)T3
CPVAP A=-4.22 Tc=369.8 K
CPVAP B=3.063E-1 Pc= 42.5 bar= 41.94 atm
CPVAP C=-1.586E-4 w = 0.153
CPVPA D=3.215E-8
Para la solución del algoritmo se utilizara la correlación generalizada del
segundo coeficiente del virial:
LABORATORIO # 4 - TERMODINAMICA QUIMICA II
Ejecución del programa:
Paso 1: Declaración de constantes características de la sustancia y funciones
del programa que se utilizan para su ejecución.
LABORATORIO # 4 - TERMODINAMICA QUIMICA II
Paso 2: conocer los valores de entalpia residual 1 y entropía residual 1 en las
condiciones 1.
Paso 3:
Determinación de T2 por el cambio de entalpia total. Cálculo iterativo
Note que la condición del ciclo while AH<=0, es válida por que se ha asumido
que el cambio en la entalpia será cero para las condiciones de trabajo. Sin
embargo, tratándose de un proceso con comportamiento real el valor de AH a
la T2 será lo más cercano a cero.
LABORATORIO # 4 - TERMODINAMICA QUIMICA II
Teniendo el valor de la T2 puede calcularse el valor de H, S. Esto se ejecuta
como:
LABORATORIO # 4 - TERMODINAMICA QUIMICA II
Ejercicio 1
Desarrolle el programa anterior y reporte los resultados obtenidos de la
temperatura a la salida del proceso de estrangulamiento, así como los valores
de ΔS, ΔH. Además, encuentre el valor de las propiedades residuales
respectivas utilizando la correlación del virial en lugar de la ecuación de RW.
¿Cuál sería el valor de ΔH y ΔS si el comportamiento del gas a la salida fuera el
de un gas ideal?
Ejemplo 2: Determinación de coeficientes de actividad.
El cálculo de coeficientes de actividad es de gran importancia en la ingeniería
química, puesto que permite estimar las condiciones del equilibrio líquido vapor
las cuales son requeridas en diferentes procesos de separación.
Hay que recordar que para un sistema multicomponente la energía de Gibbs
(exceso) viene dada por la formula siguiente:
LABORATORIO # 4 - TERMODINAMICA QUIMICA II
Esta ecuación tiene una gran importancia puesto puede utilizarse como base
para encontrar otras propiedades termodinámicas. Para la estimación del
coeficiente de actividad se han desarrollado diferentes métodos siendo uno de
los mas utilizados el modelo de Wilson.
En este modelo se toman en cuenta los efectos de las diferencias de tamaño y
fuerzas de atracción de las moléculas de las distintas especies mediante un
tratamiento basado en el concepto de composición local. Las ecuaciones que
describen dicho modelo son:
GE x x
1 ln x1 12 x2 2 ln x2 21 x1
x1 x2 RT x1 x2 x1 x2
12 21
ln 1 ln x1 12 x2 x2
x1 12 x2 x2 21 x1
12 21
ln 2 ln x2 21 x1 x1
x1 12 x2 x2 21 x1
ln 1 ln 12 1 21
ln 2 ln 21 1 12
Los parámetros ai, j asi como Vi y Vj se muestran en el anexo.
Ejecución del programa
Programe el modelo de Wilson para la determinación de coeficientes de
actividad. Compruebe sus resultados con la resolución del siguiente ejercicio:
Para el sistema acetona (1), metanol (2) , con x1=0.1, encuentre los coeficientes
de actividad de ambas sustancias a la temperatura de 60°c.
LABORATORIO # 4 - TERMODINAMICA QUIMICA II
Ejercicio 2: Encuentre los coeficientes actividad utilizando el modelo de Wilson
para el sistema metanol (1) Benceno (2), con X1=0.45 a 75°C.
Ejercicio 3: Resolver el problema 14, literal C, de la Guía de Discusión 5 utilizando
Excel.
LABORATORIO # 4 - TERMODINAMICA QUIMICA II
Anexos:
LABORATORIO # 4 - TERMODINAMICA QUIMICA II
Referencias Bibliográficas:
1. Hougen, O.A., Watson, K.M. y Ragatz, R.A. (1980). Principios de los
Procesos Químicos, Parte II, Termodiná[Link]ña: 1° Edición en
Español, Editorial Reverte.
2. Reid, Prausnitz and Poling (1987). The properties of gases and Liquids. USA:
4° edición, Editorial McGraw – Hill, Inc.
3. Smith, J.M., Van Ness, H.C., Abbott M.M. (2008). Introducción a la
Termodinámica en Ingeniería Química. México: 7° Edición (5° edición en
español), Mc. Graw-Hill, Co.
4. Faires D. y Burden R. (2002). NumericalMethods. Tercera Edición. Ohio,
Estados Unidos: Brooks Cole.
