Modulo 1 Calores Sensible y Latente

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 22

CALORES SENSIBLE Y LATENTE

1. CALOR SENSIBLE

El calor sensible es la energía intercambiada por un cuerpo sin cambiar su fase,


provocando un cambio en su temperatura.

La energía interna depende del estado termodinámico del cuerpo, descrito por
variables como presión y temperatura, o temperatura y volumen. Esta energía es una
función de estado, lo que significa que su valor solo depende del estado inicial y final
del sistema.

𝑈 = 𝑈(𝑃, 𝑇) = 𝑈(𝑇, 𝑉)
𝜕𝑈 𝜕𝑈
𝑑𝑈 = ( ) 𝑑𝑇 + ( ) 𝑑𝑉
𝜕𝑇 𝑉 𝜕𝑉 𝑇

La energía interna se describe con una diferencial total que incluye dos términos: uno
para el cambio con la temperatura a volumen constante (𝐶𝑉 ) y otro para el cambio con
el volumen a temperatura constante. 𝐶𝑉 , la capacidad calorífica a volumen constante
se sustituye directamente en la ecuación. El segundo término puede ser igual a cero
en dos casos:

• Para cualquier proceso a volumen constante sin importar la sustancia.


• Siempre que la energía interna sea independiente del volumen, sin considerar
el proceso. Esto es verdadero para gases ideales y fluidos incompresibles , y
aproximadamente cierto para gases a baja presión.

𝜕𝑈
𝑑𝑈 = 𝐶𝑉 𝑑𝑇 + ( ) 𝑑𝑉
𝜕𝑉 𝑇
Por lo cual esta ecuación se reduciría a este termino

𝑑𝑈 = 𝐶𝑉 𝑑𝑇

Para un sistema cerrado a volumen constante y adiabático se tiene que Q = ∆𝑈 e


integrando esta ecuación quedará escrita de la siguiente manera:
𝑇2
𝑄 = ∆𝑈 = ∫ 𝐶𝑉 𝑑𝑇
𝑇1

Así, el cambio de energía interna a volumen constante se expresa usando 𝐶𝑉 . En casos


donde 𝐶𝑉 es una función de la temperatura, se sustituye por un polinomio para facilitar
la integración. Más adelante, se explorará cómo manejar esta integral para distintos
procesos.
1
Por medio de un proceso similar, se puede calcular el cambio de entalpía de un
sistema termodinámico en función de sus variables independientes.

𝐻 = 𝐻(𝑃, 𝑇)

Se sabe que la entalpía también es una función de estado y se puede desarrollar el


diferencial total de la siguiente manera:
𝜕𝐻 𝜕𝐻
𝑑𝐻 = ( ) 𝑑𝑇 + ( ) 𝑑𝑃
𝜕𝑇 𝑃 𝜕𝑃 𝑇

El primer término representa el cambio de entalpía con respecto a la temperatura a


presión constante, y corresponde a 𝐶𝑃 , la capacidad calorífica a presión constante. El
segundo término describe el cambio de entalpía con respecto a la presión a
temperatura constante. Una vez más, dos circunstancias permiten que el termino final
sea igual a cero:

• Para cualquier proceso a presión constante sin importar la sustancia.


• Cada vez que la entalpía de la sustancia sea independiente de la presión, sin
considerar el proceso. Esto es precisamente cierto para gases ideales, y
aproximadamente válido para gases a baja presión.

𝜕𝐻
𝑑𝐻 = 𝐶𝑃 𝑑𝑇 + ( ) 𝑑𝑃
𝜕𝑃 𝑇

Por lo cual esta ecuación se reduciría a estos términos:

𝑑𝐻 = 𝐶𝑃 𝑑𝑇

Al integrar esta función, se obtiene la siguiente función:


𝑇2
∆𝐻 = ∫ 𝐶𝑃 𝑑𝑇
𝑇1

El cambio de entalpia es igual a la integral desde una temperatura inicial hasta una
temperatura final del 𝐶𝑃 por dT.

2. A continuación, se verá cómo se puede determinar 𝑪𝑷 en función de la


temperatura.

Se debe tener en cuenta que 𝐶𝑃 en el tema de calorimetría y usualmente para todos los
gases es aproximadamente dependiente de la temperatura, no se puede considerar
constante.

2
En termodinámica, una bomba calorimétrica mide el calor intercambiado durante
reacciones químicas o procesos físicos para calcular la entalpía y otras propiedades
energéticas. Permite determinar la energía absorbida o liberada, evaluar el valor
calorífico de combustibles mediante la calorimetría de combustión y analizar diversos
procesos termodinámicos, como la disolución. A través de experimentos realizados
en una bomba calorimétrica se midió 𝐶𝑃 de diversos gases en función de la
temperatura. Conociendo el cambio de temperatura y la cantidad de calor recibido, se
determinó 𝐶𝑃 . Los datos obtenidos se ajustaron mediante un polinomio, que se
analizará a continuación.
𝐶𝑝
= 𝐴 + 𝐵𝑇 + 𝐶𝑇 2 + 𝐷𝑇 −2
𝑅
Este polinomio analiza la dependencia de la capacidad calorífica con la temperatura,
Los coeficientes que acompañan este polinomio se pueden consultar en la bibliografía
del libro de Smith o también en el pdf que se encuentra en este módulo llamado
capacidades caloríficas.
𝐶𝑝
El cociente , donde R es la constante universal de los gases, se ajusta a un polinomio
𝑅
de 4 términos con constantes A, B, C, y D, y la T como variable independiente. Este
polinomio se utiliza en la integral. Para algunas sustancias, las constantes C o D
𝐶𝑝
pueden ser cero. Se usa porque es adimensional.
𝑅

A continuación, se procederá a analizar la tabla C1 del libro de Smith.

En la parte superior, se encuentra el polinomio mencionado anteriormente. También


se indica que la temperatura se mide en Kelvin, con un rango de valores que va desde
298 K hasta Tmáx , que según la tabla es 1500 K. En la tabla también se reportan el
nombre de las especies químicas, su fórmula, la temperatura máxima para cada
especie (que define el rango de aplicación del polinomio), y la capacidad calorífica a

3
298 K. Este último valor es una constante que debe coincidir al sustituir 298 K en el
polinomio. Además, se incluyen los valores de las constantes A, B, C, y D.

La constante A del metano es 1,702. Para las demás constantes, se puede notar que
están afectadas por un exponente. Por ejemplo, se reporta 103 B, 106 C y 10-5 D. Estos
datos indican que la constante está multiplicada por ese exponente. Para calcular la
constante B del metano, que tiene un valor de

103 𝐵 = 9,081
Para encontrar el valor directamente de B se debe despejar y quedaría de la siguiente
manera:
9,081
𝐵 = = 9,081𝑥10−3
103
Se cambia el signo al exponente porque lo que indica esto es que reportan en la tabla
el valor de B multiplicado por 10 3 . Ocurre lo mismo para los valores de C y D, D en este
caso no posee valor, pero si lo tuviera quedaría 1𝑥105 positivo.

Esta sería la forma en la que se debe leer las tablas, no todas las especies tienen una
constante D o C, Si no se encuentran valores reportados directamente se toma como
cero y se elimina ese término del polinomio.

Después de calcular el polinomio de cp se reemplaza en la integral con la formula del


cambio de entalpia.
𝑇2
∆𝐻 = ∫ 𝐶𝑃 𝑑𝑇
𝑇1

𝐶𝑝 𝐶𝑝
Se observa que ese polinomio es Entonces en la integral para que sea
𝑅 𝑅
homogéneamente dimensional se tendría que multiplicar la integral por R.

Para calcular la energía interna esta depende del 𝐶𝑉 y no del 𝐶𝑃 , no se reporta un


polinomio de 𝐶𝑉 se debe tener en cuenta que para un gas ideal la relación entre 𝐶𝑉 y 𝐶𝑃
tiene una ecuación:

𝐶𝑃 − 𝐶𝑉 = 𝑅

Con esta ecuación se puede encontrar una relación directa entre 𝐶𝑉 y 𝐶𝑃 , se despeja
𝐶𝑉 Y si esta ecuación se divide entre R se obtiene una ecuación la cual puede
relacionar el 𝐶𝑃 entre R y el 𝐶𝑉 entre R, así

4
𝐶𝑉 𝐶𝑃
= −1
𝑅 𝑅
Por esta razón no se reportan tablas para 𝐶𝑉 , ya que por medio de 𝐶𝑃 es posible
calcular muy fácilmente esta relación.

Mediante estas dos ecuaciones ya se calcula tanto entalpía utilizando 𝐶𝑃 como


energía interna utilizando 𝐶𝑉 .
𝐶𝑝
Entalpia: = 𝐴 + 𝐵𝑇 + 𝐶𝑇 2 + 𝐷𝑇 −2
𝑅

𝐶𝑉 𝐶𝑃
Energía interna : = −1
𝑅 𝑅

2.1 Para calcular el 𝑪𝑷 de una mezcla de gases


En una mezcla de gases A, B y C, el 𝐶𝑃 de la mezcla se calcula como un promedio
ponderado, de la siguiente manera:

𝐶𝑃𝑚𝑒𝑧𝑐𝑙𝑎 = 𝑦𝐴 𝐶𝑃𝐴 + 𝑦𝐵 𝐶𝑃𝐵 + 𝑦𝐶 𝐶𝑃𝐶

En donde 𝐶𝑃𝐴 , 𝐶𝑃𝐵 y 𝐶𝑃𝐶 son las capacidades caloríficas molares de A, B y C puras en
estado de gas ideal y 𝑦𝐴 , 𝑦𝐵 , y 𝑦𝐶 representan las fracciones molares.

2.1.1 Cálculo de la integral de 𝑪𝑷

La diferencial de entalpia es igual a 𝐶𝑃 por diferencial de temperatura y al integrar


queda:
𝑇2
∆𝐻 = ∫ 𝐶𝑃 𝑑𝑇
𝑇1

𝐶𝑝
Teniendo en cuenta que el polinomio planteado anteriormente es el cociente de , en
𝑅
ese caso para sustituir ese polinomio en la integral se va a hacer lo siguiente, al ser R
constante se multiplica por R al inicio de la integral para que dé simplemente 𝐶𝑃 .
𝑇2 𝐶𝑝
∆𝐻 = 𝑅 ∫ 𝑑𝑇
𝑇1 𝑅
Se multiplica el polinomio por R para obtener las mismas unidades de entalpía
seleccionando las unidades de R que mejor convenga.

𝐶𝑝
= 𝐴 + 𝐵𝑇 + 𝐶𝑇 2 + 𝐷𝑇 −2
𝑅

5
Como se vio anteriormente como plantear el polinomio de 𝐶𝑃 en función de la
temperatura, ahora se aplica para el cálculo de entalpía.

La evaluación de la integral se consigue mediante la sustitución de 𝐶𝑃 como función de


la temperatura
𝑇 𝐶𝑝 𝐵 𝐶 𝐷 𝜏− 1
∫ 𝑑𝑇 = 𝐴𝑇0 (𝜏 − 1) + 𝑇02 (𝜏 2 − 1) + 𝑇03 (𝜏 3 − 1) + ( )
𝑇0 𝑅 2 3 𝑇0 𝜏

Donde
𝑇
𝜏=
𝑇0

𝑇 = 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑓𝑖𝑛𝑎𝑙
𝑇0 = 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑖𝑛𝑐𝑖𝑎𝑙

Mediante un procedimiento algebraico esta ecuación se puede reorganizar y quedará


como resultado la capacidad calorífica media, así:

〈𝐶𝑃 〉𝐻 𝐵 𝐶 𝐷
= 𝐴+ 𝑇0 (𝜏 − 1) + 𝑇20 (𝜏2 + 𝜏 + 1) +
𝑅 2 3 𝜏𝑇20

Así quedara la ecuación para hallar la diferencia de entalpia, teniendo en cuenta que
R es la constante de los gases.

∆𝐻 = 〈𝐶𝑃 〉𝐻 (𝑇 − 𝑇0 )
Ejercicios para comprender estas fórmulas anteriormente vistas

Ejercicio 1

Calcule el calor necesario para aumentar la temperatura de 1 mol de metano de 260 a


600 °C en un proceso de flujo estable a una presión lo suficientemente baja para que
el metano se considere un gas ideal.

En la tabla podemos consultar el valor de las constantes.

6
Solución

Se procede a convertir las temperaturas a Kelvin utilizando la siguiente relación:

T (k) =T (°C) + 273.15

Se convierten la temperatura inicial y final del sistema dando como resultado:

𝑇0 = 533,15 𝐾

𝑇 = 873,15 𝐾

Luego se halla el termino 𝜏 con la ecuación vista anteriormente:


𝑇 873,15
𝜏= = = 1,6377
𝑇0 533,15

Luego, se calcula la cantidad de calor necesario. Al tratarse de un proceso a presión


constante, el calor es igual al cambio de entalpía. Así, se sustituye en la ecuación de la
siguiente manera:
873,15 𝐶𝑝
𝑄 = ∆𝐻 = 𝑅 ∫ 𝑑𝑇
533,15 𝑅
Se puede desarrollar esta integral también de forma manual reemplazando el
𝐶𝑝
polinomio de en la integral y resolviendo, así.
𝑅
873,15
2 −2
𝑄 = ∆𝐻 = 𝑅 ∫ 𝐴 + 𝐵𝑇 + 𝐶𝑇 + 𝐷𝑇 𝑑𝑇
533,15

En la tabla C1 del libro de Smith se consultaron los valores de las constantes,

𝑅 = 8,314

𝐴 = 1,702

𝐵 = 9,081𝑋10 −3

𝐶 = −2,164𝑋10−6
Conociendo la solución a la integral anteriormente mencionada.
𝑇𝐶𝑝 𝐵 𝐶 𝐷 𝜏− 1
∫ 𝑑𝑇 = 𝐴𝑇0 (𝜏 − 1) + 𝑇02 (𝜏 2 − 1) + 𝑇03 (𝜏 3 − 1) + ( )
𝑇0 𝑅 2 3 𝑇0 𝜏

Como para el metano no existe valor para la constante D, se toma como cero y se
elimina ese término.

7
La solución a la integral sin el termino D quedara así:
𝑇 𝐶𝑝 𝐵 𝐶
∫ 𝑑𝑇 = 𝐴𝑇0 (𝜏 − 1) + 𝑇02 (𝜏 2 − 1) + 𝑇03 (𝜏 3 − 1)
𝑇0 𝑅 2 3

Nuestra ecuación para calcular el calor queda de la siguiente forma:


𝐵 2 2 𝐶
𝑄 = ∆𝐻 = 𝑅 [𝐴𝑇0 (𝜏 − 1) + 𝑇0 (𝜏 − 1) + 𝑇03 ( 𝜏 3 − 1)]
2 3
Ahora si se procede a calcular el calor reemplazando los valores de las constantes que
ya tenemos:

𝑄 = ∆𝐻 = (8,314) [1,702(533,15) (1,6377 − 1)

9,081𝑋10 −3
+ (533,15)2 (1,6377 2 − 1)
2
2,164𝑋10 −6
− (533,15)3 (1,6377 3 − 1)]
3
Esto nos da un valor de calor que quedaría en J/mol. Dado que tenemos 1 mol de
metano, al multiplicar ese valor por 1 mol, el resultado se expresaría simplemente en
Joules. Por lo tanto, el calor calculado es igual a:

𝑸 = 𝟏𝟗𝟕𝟕𝟖 𝑱

Ejercicio 2
¿Cuál será la temperatura final cuando se agrega una cantidad de calor de
0,4𝑋106 (Btu) a 25 (lb mol) de amoniaco inicialmente a 500 (°F) en un proceso de flujo
estable a 1 (atm)?

Primero se puede notar que las unidades están dadas en el sistema inglés, se puede
calcular el 𝐶𝑃 en el sistema inglés o bien convertir todo al sistema internacional.

Se sabe que es un proceso a presión constante, al decir que es un proceso de flujo


estable se determina que es un sistema abierto, entonces a presión constante la
cantidad de calor transferido es el número de moles por el cambio de entalpía.

𝑄 = 𝑛∆𝐻

Se procede a calcular inicialmente el cambio de entalpía.

𝑄 0,4𝑋106
∆𝐻 = = = 16000 (𝐵𝑡𝑢)(𝑙𝑏 𝑚𝑜𝑙 )−1
𝑛 25

8
Esto significa que en este proceso de flujo por cada mol que fluye se transfieren 16000
Btu, luego este factor se convierte al sistema internacional, (los factores de conversión
los pueden encontrar en el pdf de tablas de conversión de unidades).
16000
∆𝐻 = = 37218𝐽𝑚𝑜𝑙 −1
0,4299
También se procede a transformar la temperatura de °F a K

(°𝐹 − 32) ∗ 5
𝑇0 = + 273,17
9
Reemplazando 500 °F

(500 − 32) ∗ 5
𝑇0 = + 273,17 = 533,15 𝐾
9
Con estos valores ya se puede pasar a plantear el polinomio de Cp.
𝑇
𝑄 = ∆𝐻 = ∫ 𝐶𝑝 𝑑𝑇
533,15

Se procede a buscar los valores de las constantes en la tabla C1para el amoniaco.

𝐴 = 3,578

𝐵 = 3,020𝑋10−3

𝐷 = −0,186𝑥105

Luego se reemplazan todos los valores que se tienen en la fórmula, quedando de la


siguiente manera:
𝑇
−2
37218𝐽𝑚𝑜𝑙 −1 = 8,314 𝐽/𝑚𝑜𝑙𝐾 ∫ 3,578 + 3,020𝑋10 −3 𝑇 − 0,186𝑥105 𝑇 𝑑𝑇
533,15𝐾

Se basarán en esta ecuación y se procederá con la resolución algebraica. Luego, se


integrará la ecuación, obteniendo el resultado.

3,020𝑋10−3 2
4476,55 𝐾 = 3,578(𝑇 − 533,15𝐾 ) + (𝑇 − 533,15𝐾 2 )
2
+ 0,186𝑋105 (𝑇 −1 − 533,15𝐾 −1 )

Al resolver estas operaciones algebraicamente se puede obtener una ecuación asi:

4476,55 = 3,578𝑇 − 1907,61 + 0,00151𝑇 2 − 429,2159 + 0,186𝑋105 𝑇 −1 − 34,887

9
En esta ecuación, se puede plantear el polinomio, identificando la temperatura como
la incógnita. Aunque no es posible despejar directamente la temperatura, se agrupan
los términos constantes que no incluyen T, dando como resultado:

3,578𝑇 + 0,00151𝑇 2 + 0,186𝑋105 𝑇 −1 = 6848,2

Para encontrar el valor de T, se puede usar el método de Newton-Raphson, detallado


en el PDF, o la función "Buscar Objetivo" en Excel.

3,578𝑇 + 0,00151𝑇 2 + 0,186𝑋105 𝑇 −1 − 6848,2 = 0

Se iguala la ecuación a cero para aplicar el método gráfico.

Mediante el método grafico se puede resolver haciendo una columna de temperaturas


y otra de la función evaluada en cada temperatura, así:

La temperatura incrementa porque se sabe gracias al título que se agrega una cantidad
de calor, luego lo que se hace es graficar estos datos y buscar en donde esa ecuación
da cero.

10
Se puede ver que la curva corta en cero en un valor de 1250 K entonces esa es la
solución T=1250 K

Si no se puede graficar, se puede interpolar utilizando los mismos datos para


determinar el valor de la temperatura.

Como se observan en los datos hay un cambio de positivo a negativo en la función para
los valores de temperatura 1200 y 1300, estos son los valores que se utilizan para la
interpolación.

3. CALOR LATENTE
El calor latente es la energía necesaria para cambiar el estado de una sustancia sin
alterar su temperatura. Se divide en dos tipos principales: el calor latente de fusión,
que es la energía requerida para pasar de sólido a líquido (o viceversa) en el punto de
fusión, y el calor latente de vaporización, que es la energía necesaria para pasar de
líquido a gas (o viceversa) en el punto de ebullición (Smith et al., 2017).

Se está familiarizado con el agua, cuyo estado sólido, el hielo, tiene una estructura
cúbica ordenada. Mientras que para el agua solo existe esta forma, algunos metales y
aleaciones pueden presentar diferentes fases sólidas.

11
Ahora se verá qué ecuaciones existen para poder calcular los calores latentes. Pero
antes se puede observar que existen algunos valores ya predeterminados que se
pueden encontrar en la bibliografía.

Por ejemplo, en la tabla se presentan propiedades de diferentes materiales, incluyendo


el punto de fusión, el punto de ebullición, y los calores latentes de fusión y
vaporización.
𝐽
Para el aluminio, el calor latente de fusión es de 3,97 𝑋 105 𝐾𝑔 . La fusión es el cambio
de fase de sólido a líquido, que requiere adicionar calor. En el caso del aluminio, se
𝐽
necesita agregar 3,97 𝑋 105 para lograr este cambio de estado.
𝐾𝑔

El calor latente de vaporización corresponde al cambio de estado de líquido a vapor, lo


que también requiere la adición de energía. Por ejemplo, para el agua, se necesitan
𝐽
2,26 𝑋 106 𝐾𝑔 para convertirlo en vapor.

Las transiciones de fase opuestas tienen el mismo valor de energía, pero con signo
opuesto. Por ejemplo, la solidificación es la transformación de líquido a sólido y para
𝐽
el aluminio sería −3,97 𝑋 103 . De manera similar, la condensación, el proceso
𝐾𝑔
𝐽
opuesto a la evaporación, tiene un valor negativo, como −2,26 𝑋 106 para el agua.
𝐾𝑔

12
También existen varias ecuaciones para poder estimar estos datos, de las primeras
ecuaciones que se presentaron para la estimación de calores latentes es la siguiente
ecuación que se denomina

4.1 REGLA DE TROUTON

∆𝐻𝑛
~10
𝑅𝑇𝑛

La ecuación muestra que el cociente entre el calor latente de ebullición y el producto


de la constante de los gases ideales y la temperatura de ebullición es
aproximadamente 10. Este valor es adimensional.

Para el argón, el cociente es cerca de 8.0. Como se puede notar casi todos los valores
están cercanos a 10, sin embargo, no es exactamente 10. Por lo que la regla de Trouton
lleva usualmente cierto grado de error.

SUSTANCIA

Ar 8
N2 8,7
O2 9,1
HCl 10,4
C6H6 10,5
H2S 10,6
H2O 13,1

4.2 REGLA DE RIEDEL

Otra regla que se va a analizar es la REGLA DE RIEDEL la cual también es una ecuación
empírica su nombre es debido a su inventor, en este caso parte de la regla de Trouton.

13
∆𝐻𝑛
𝑅𝑇𝑛

He introduce nuevas variables

∆𝐻𝑛 1,092( ln 𝑃𝑐 − 1,013)


=
𝑅𝑇𝑛 0,930 − 𝑇𝑟𝑛
∆𝐻𝑛
La relación adimensional se mantiene, pero se incorporan constantes críticas
𝑅𝑇𝑛
como la presión crítica y la temperatura reducida. Estas propiedades críticas, que se
consultan en tablas ( o en el PDF llamado propiedades de las especies puras), definen
las condiciones en las cuales un gas no se puede comprimir más y se utilizan en
diversas ecuaciones de estado.

A partir de las constantes críticas, se definen las variables reducidas. Por ejemplo, la
temperatura reducida se calcula dividiendo la temperatura del gas entre su
temperatura crítica.
𝑇
𝑇𝑟 =
𝑇𝑐

Las propiedades reducidas se refieren a la temperatura reducida, presión y volumen


reducidos, calculados con respecto a sus valores críticos. Estas variables son
adimensionales.

Siempre se calculan los calores latentes a la temperatura de ebullición, que es la


temperatura en la que ocurre el cambio de fase. Esta expresión empírica se basa en
datos ajustados mediante métodos estadísticos para diferentes gases, y los errores
rara vez superan el 5%. Estas propiedades para especies puras se pueden conseguir
en el apendice B del libro de Smith o en el pdf Proiedades de las especies puras.

Por ejemplo, al verificar para el caso del agua, se observa esta precisión.

∆𝐻𝑛 1,092( ln 𝑃𝑐 − 1,013)


=
𝑅𝑇𝑛 0,930 − 𝑇𝑟𝑛

∆𝐻𝑛 1,092(ln 220,55 − 1,013)


=
𝑅𝑇𝑛 0,930 − 0,577

Se buscan los valores de Presión crítica para el agua de 220,55 bar y con la temperatura
critica se halla el valor de temperatura reducida lo cual da 0,577, al reemplazar en la
formula da como resultado 13,56

14
∆𝐻𝑛
= 13,56
𝑅𝑇𝑛

Luego esto se multiplica por la constante de los gases R y la temperatura de ebullición


la cual es 100 °C, se convierte a kelvin para poder hallar el valor de ∆𝐻𝑛 lo cual daría un
valor de:

∆𝐻𝑛 = 13,56 (8,314)(373,15) = 42065 𝐽𝑚𝑜𝑙 −1


Esto significaría que para que el agua cambie de estado líquido a estado vapor se
requieren adicionar 42065 𝐽 de energía por cada mol de agua.

4.3 LA CORRELACIÓN DE WATSON

Otra ecuación útil para estimar los calores latentes a diferentes temperaturas es la
Correlación de Watson. Esta ecuación permite calcular el calor latente de vaporización
∆𝐻2 a la temperatura 𝑇𝑟2 utilizando el calor latente de vaporización ∆𝐻1 a la
temperatura 𝑇𝑟1 .

∆𝐻2 1 − 𝑇𝑟2 0,38


=( )
∆𝐻1 1 − 𝑇𝑟1
Esta ecuación relaciona los calores latentes de vaporización con las temperaturas
reducidas, las temperaturas reducidas como ya se vio anteriormente se relacionan con
la temperatura a la que se encuentra el compuesto y su temperatura critica.
𝑇
𝑇𝑟 =
𝑇𝑐

A continuación, se verá un ejemplo para entender mejor esta ecuación.

Ejemplo

Dado que el calor latente de vaporización del agua a 100 °C es 2257 J g -1, calcule el
calor latente a 300 °C.

Solución

Se procede a reconocer las variables dadas en el enunciado, las temperaturas deben


convertirse a kelvin debido a que la temperatura critica en la literatura se consigue en
Kelvin.

∆𝐻2 = ?

∆𝐻1 = 2257 J g − 1

15
373,15
𝑇𝑟1 = = 0,577
647,1

573,15
𝑇𝑟2 = = 0,886
647,1

Una vez se obtienen los valores, se despeja de la ecuación el ∆𝐻2

1 − 𝑇𝑟2 0,38
∆𝐻2 = ∆𝐻1 ( )
1 − 𝑇𝑟1
Y se reemplazan los valores

1 − 0,886 0,38
∆𝐻2 = 2257 ( )
1 − 0,577
Lo que dará un valor de calor latente de:

∆𝐻2 = 1371J g − 1

El valor proporcionado en las tablas de vapor es de 1406 J g-1. Como se puede comparar
los valores son muy cercanos, se aproximan mucho entonces esta correlación de
Watson es confiable.

4.4 ECUACIÓN DE CLAPEYRON

Otra ecuación que se puede utilizar para los calores latentes de sustancias puras es la
denominada ECUACIÓN DE CLAPEYRON que se encuentra expresada en su forma
diferencial así:

𝑑𝑃 𝑠𝑎𝑡𝑢𝑟𝑎𝑑𝑜
∆𝐻 = 𝑇∆𝑉
𝑑𝑇
En donde:

∆𝐻 = 𝑐𝑎𝑙𝑜𝑟 𝑙𝑎𝑡𝑒𝑛𝑡𝑒

∆𝑉 = 𝑐𝑎𝑚𝑏𝑖𝑜 𝑑𝑒 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑞𝑢𝑒 𝑎𝑐𝑜𝑚𝑝𝑎ñ𝑎 𝑒𝑙 𝑐𝑎𝑚𝑏𝑖𝑜 𝑑𝑒 𝑓𝑎𝑠𝑒

𝑃 𝑠𝑎𝑡𝑢𝑟𝑎𝑑𝑜 = 𝑝𝑟𝑒𝑠𝑖ó𝑛 𝑑𝑒 𝑠𝑎𝑡𝑢𝑟𝑎𝑐𝑖ó𝑛

La ecuación de Clapeyron relaciona el calor latente de vaporización con la temperatura


y el cambio de volumen durante el cambio de fase. En la vaporización, el volumen del
vapor es mayor que el del líquido, por lo que el cambio de volumen es positivo. La
ecuación también incluye el cambio en la presión de saturación con respecto a la
temperatura. Con esta ecuación, la temperatura ya es conocida: es la temperatura de
ebullición del material en kelvin. El cambio de volumen, que se obtiene de tablas,
también se incluye en el cálculo.

16
Para estimar la presión de saturación, se puede utilizar la ecuación de Antoine.

En donde:

P0 = Presión de saturación del sistema

T = Temperatura del sistema

A, B y C = Constantes que dependen de cada sustancia

Esta ecuación solamente se utiliza para un sistema en equilibrio liquido-vapor, las


constantes se pueden leer en la tabla B2 del libro de Smith o también en el pdf titulado tablas
ecuación de Antoine.

Las unidades de la presion para la ecuacion de Antoine son de Kpa y para la


temperatura en °C; esto porque las constantes tambien tienen unas unidades.

La tabla proporciona las especies químicas, las constantes para la ecuación de


Antoine, el intervalo de temperaturas aplicables, el calor latente normal y la
temperatura normal de ebullición a 1 atmósfera. La ecuación de Antoine permite
calcular la presión de saturación y estimar el punto de ebullición al despejar la
temperatura, que es cuando la presión atmosférica iguala la presión de saturación.

Regresando a la ecuación de Clapeyron de esta forma podemos estimar la presión de


saturación.

𝑑𝑃 𝑠𝑎𝑡𝑢𝑟𝑎𝑑𝑜
∆𝐻 = 𝑇∆𝑉
𝑑𝑇

17
Esta ecuación también se puede escribirla de la siguiente manera: (Cengel & Boles, 2019)

𝑑𝑃 ℎ𝑓𝑔
( ) =
𝑑𝑇 𝑠𝑎𝑡 𝑇 ∗ 𝑉𝑓𝑔

En donde

ℎ𝑓𝑔 = ∆𝐻 𝑜 𝑒𝑛𝑡𝑎𝑙𝑝í𝑎 𝑓𝑖𝑛𝑎𝑙 − 𝑒𝑛𝑡𝑎𝑙𝑝í𝑎 𝑖𝑛𝑐𝑖𝑎𝑙

𝑉𝑓𝑔 = ∆𝑉 𝑜 𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑓𝑖𝑛𝑎𝑙 − 𝑣𝑜𝑙𝑢𝑚𝑒𝑛 𝑖𝑛𝑖𝑐𝑖𝑎𝑙


𝑑𝑃
Si se tiene un diagrama de presión-temperatura la pendiente sería el diferencial (𝑑𝑇 )
𝑠𝑎𝑡

A continuación, se va a analizar un ejercicio para comprender como utilizar esta ecuación de


Clapeyron.

Ejercicio
En este ejercicio se va a calcular la entalpía de vaporización de una sustancia a partir de datos
de P-V-T que se pueden encontrar en tablas.

Mediante la ecuación de Clapeyron, estime el valor de la entalpía de vaporización del


refrigerante 134a a 20 °C y compárelo con el valor tabulado.

𝑑𝑃 ℎ𝑓𝑔
( ) =
𝑑𝑇 𝑠𝑎𝑡 𝑇 ∗ 𝑉𝑓𝑔

Se procede a despejar la entalpía de vaporización de la ecuación de Clapeyron

𝑑𝑃
ℎ𝑓𝑔 = 𝑇 ∗ 𝑉𝑓𝑔 ( )
𝑑𝑇 𝑠𝑎𝑡

Para resolver este ejercicio se utilizó la tabla A11 del libro de Cengel, estas tablas las pueden
conseguir también en el PDF llamado Tablas Cengel.

Se sabe que se utilizara el refrigerante 134a a 20 °C el cambio en el volumen también lo puedo


leer de las tablas.

18
𝑚3
𝑉𝑓𝑔 = (𝑉𝑔 − 𝑉𝑓 ) = 0,035969 − 0,0008161 = 0,035153
𝐾𝑔
𝑑𝑃
Se requiere también saber cuánto vale la pendiente (𝑑𝑇 ) a 20 °C
𝑠𝑎𝑡

Para esto se utiliza un delta P y un delta T pequeño alrededor de 20 °C


𝑑𝑃 ∆𝑃 𝑃𝑠𝑎𝑡 24 °𝐶 − 𝑃𝑠𝑎𝑡 16 °𝐶
( ) =( ) =
𝑑𝑇 𝑠𝑎𝑡 20 °𝐶 ∆𝑇 𝑠𝑎𝑡 20 °𝐶 24 °𝐶 − 16 °𝐶

Estos valores de P se obtienen también de las tablas, se procede a reemplazar en la


ecuación:
𝑑𝑃 ∆𝑃 646,18 − 504,58 𝐾𝑃𝑎 𝐾𝑃𝑎
( ) =( ) = = 17,70
𝑑𝑇 𝑠𝑎𝑡 20 °𝐶 ∆𝑇 𝑠𝑎𝑡 20 °𝐶 8 °𝐶 𝐾

Sabiendo que ∆𝑇(°𝐶 ) = ∆𝑇 (𝐾 ) al sustituir se obtiene:

𝑚3 𝐾𝑃𝑎 1 𝐾𝐽
ℎ𝑓𝑔 = (293,15𝐾 ) (0,035153 ) (17,70 )( )
𝐾𝑔 𝐾 1𝐾𝑃𝑎 ∗ 𝑚3
1 𝐾𝐽
Este último factor (1𝐾𝑃𝑎∗𝑚3 ) se coloca ahí para eliminar los metros cúbicos y los
kilopascales.
𝐾𝐽
ℎ𝑓𝑔 = 182,40
𝐾𝑔

Para hacer la comparación se revisa el valor de entalpia de a 20 °C tabulado

19
Esta pequeña diferencia entre los valores se debe a la aproximación utilizada al determinar
la pendiente de la curva de saturación a 20 °C.

4.5 ECUACIÓN DE CLAUSIUS-CLAPEYRON

La ecuación de Clapeyron se puede simplificar para varios cambios de fase liquido-vapor o


solido-vapor con algunas suposiciones.

A bajas presiones el volumen especifico que ocupa el vapor es mucho mayor que el volumen
del liquido 𝑉𝑔 ≫ 𝑉𝑓 por lo que el cambio de volumen 𝑉𝑓𝑔 ≅ 𝑉𝑔 . Si se considera el vapor como
𝑅𝑇
gas ideal entonces 𝑉𝑔 = . Si se sustituyen y reorganizan estas dos aproximaciones en la
𝑃
ecuación original de Clapeyron se puede obtener lo siguiente:

𝑑𝑃 ℎ𝑓𝑔 𝑑𝑃
( ) = ( )
𝑃 𝑠𝑎𝑡 𝑅 𝑇 2 𝑠𝑎𝑡

Al integrar la ecuación de Antoine, se obtiene un logaritmo natural en el lado derecho y


una diferencial de temperatura . La ecuación final se presenta como:

𝑃2 ℎ𝑓𝑔 1 1
ln ( ) = ( − )
𝑃1 𝑠𝑎𝑡 𝑅 𝑇1 𝑇2 𝑠𝑎𝑡

Esta ecuación se conoce como la ecuación de Clausius-Clapeyron. También se aplica


en regiones sólido-vapor, donde el calor latente se sustituye por la entalpía de
sublimación.

A continuación, se presenta un ejercicio para saber cómo utilizar esta ecuación.

Ejercicio

Se ha estudiado la presión de vapor del etanol en función de la temperatura obteniéndose los


siguientes datos:

20
Calcular gráficamente ∆𝐻𝑣𝑎𝑝 y la temperatura de ebullición normal.

Esta ecuación se debe linealizar para que quede como la ecuación de una línea recta, para
tener en cuenta la ecuación de una línea recta es 𝑦 = 𝑏𝑥 + 𝑐.

Entonces al linealizar la ecuación quedara de la siguiente manera:

∆𝐻𝑣𝑎𝑝 1
ln 𝑃 = − ∗ +𝐶
𝑅 𝑇
1
En donde la y seria ln P, la x de la ecuación es el inverso de la temperatura 𝑇 y una constante
c que representa el intercepto.

Entonces partiendo de estos datos se trascriben a una hoja de cálculo en Excel, se requieren
para hallar esos datos la inversa de la temperatura contra el ln P

Se agregan las 2 columnas para calcular la inversa de la temperatura y el logaritmo natural de


la presión y se grafica. Por la linealización que se le hizo a la ecuación de inicio se conoce que
la pendiente debería ser negativa.

Como se puede observar el valor de la pendiente de esta recta es -37,312. Entonces si


∆𝐻𝑣𝑎𝑝
− = −37,312, se puede despejar ∆𝐻𝑣𝑎𝑝 y el valor de R dependerá del sistema de
𝑅
unidades que se esté manejando.

En este caso se está utilizando el sistema internacional, el valor de R seria 8,314 J/mol
𝐽
∆𝐻𝑣𝑎𝑝 = (−37,312)(−8,314) = 310,212
𝑚𝑜𝑙
Como la pendiente es negativa y el R también es negativo daría positiva la entalpía de
vaporización.

REFERENCIAS

1. (Smith et al., 2017)


Smith, J. M., Van Ness, H. C., Abbott, M. M., & Swihart, J. C. (2017). Introducción
a la termodinámica en ingeniería química. McGraw-Hill Education.

21
2. (Cengel & Boles, 2019)
Cengel, Y. A., & Boles, M. A. (2019). THERMODYNAMICS: AN ENGINEERING
APPROACH (9a ed.). McGraw-Hill Education.

Videos Modulo 1
1. https://youtu.be/oyvOqJikX6A?si=QivyDSV5jEMiQu_8
2. https://www.powtoon.com/ws/bJMwV2oj5Sx/1/m
3. https://www.powtoon.com/ws/e4htS7uGsx7/1/m

22

También podría gustarte