Análisis Del Factor de Compresibilidad de Los Gases

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 3

Análisis del factor de compresibilidad de los gases

Debido a que en el reactor se trabajará a una presión muy alta, debe realizarse un análisis de la
no idealidad de los gases, ya que no son condiciones favorables para modelizar a través de la
ecuación de estado de gases ideales.

Para ello, calcularemos el factor de compresibilidad (Z) de la mezcla, cuyo valor nos indicará
cuán lejos de un modelo ideal se encuentra la mezcla. Recordemos que un valor de Z=1
corresponde al gas ideal, y que valores significativamente distantes a 1 representan no idealidad.

Para el cálculo, utilizaremos la ecuación de Redlich-Kwong-Soave , que es válida en a altas


presiones y temperaturas.

La expresión es la siguiente:
𝑅∗𝑇 𝑎
𝑃= − (1)
𝑉−𝑏 𝑉∗ (𝑉 + 𝑏)
Donde a y b están definidos por las siguientes expresiones:

0.42747 ∗ 𝑅 2 ∗ 𝑇𝑐 2
𝑎= ∗𝛼
𝑃𝑐
0.08664 ∗ 𝑅 ∗ 𝑇𝑐
𝑏=
𝑃𝑐
𝛼 = [1 + (0.48508 + 1.5517 ∗ 𝑤 + 0.15613 ∗ 𝑤 2 ) ∗ (1 − 𝑇𝑟 0.5 )]2
Siendo Tc la temperatura crítica, Pc la presión crítica y w el factor acéntrico, todos los términos
referidos a la mezcla. Tr es la temperatura reducida, siendo el cociente entre la temperatura T
y la temperatura crítica.

Para el cálculo de los parámetros mencionados en el párrafo anterior, utilizaremos la regla


sencilla de mezclado, debido a que raramente podemos conocer los mismos de bibliografía
(este procedimiento es propuesto por “Introducción a la termodinámica en Ingeniería Química,
Smith, Van Ness):

𝑇𝑐(𝑚𝑒𝑧𝑐𝑙𝑎) = ∑ 𝑦𝑖 ∗ 𝑇𝑐𝑖

𝑃𝑐(𝑚𝑒𝑧𝑐𝑙𝑎) = ∑ 𝑦𝑖 ∗ 𝑃𝑐𝑖

𝑤(𝑚𝑒𝑧𝑐𝑙𝑎) = ∑ 𝑦𝑖 ∗ 𝑤𝑖

Consideraremos despreciables las fracciones molares de inertes, ya que son muy pequeñas.
Tomando alimentación estequiométrica al reactor: yH2=2/3 ; yN2=1/3

De bibliografía (“Manual del Ingeniero Químico, Perry) obtenemos los valores de Tc, Pc y w
para el nitrógeno y el hidrógeno:

Tc (H2)= 33.25 K Pc(H2)= 12.8 atm w(H2)= -0.22

Tc (N2)= 126.05 K Pc(N2)= 33.5 atm w(N2)= 0.04


La temperatura de entrada al reactor será 773 K (la más desfavorable, ya que al tratarse de
una reacción exotérmica irá aumentando y favorecerá la idealidad) y la presión 300 atm (el
reactor operará en un rango 250-300 atm, analizaremos el caso más desfavorable).

Sabiendo que el factor de compresibilidad (Z) está definido de la siguiente forma:


𝑃∗𝑉
𝑍=
𝑅∗𝑇
Acomodando la expresión (1) nos queda:

𝑍 3 − 𝑍 2 + (𝐴∗ − 𝐵∗ − (𝐵∗ )2 ) ∗ 𝑍 − 𝐴∗ ∗ 𝐵∗ = 0
Siendo:
𝑎∗𝑃
𝐴∗ =
𝑅2 ∗ 𝑇 2
𝑏∗𝑃
𝐵∗ =
𝑅∗𝑇
De esta forma, queda la ecuación cúbica expresada en función de Z, que será útil para analizar
a idealidad de la mezcla.

Realizando los cálculos correspondientes:

𝑇𝑐(𝑚𝑒𝑧𝑐𝑙𝑎) = 64.2 𝐾
𝑃𝑐(𝑚𝑒𝑧𝑐𝑙𝑎) = 19.7 𝑎𝑡𝑚
𝑤(𝑚𝑒𝑧𝑐𝑙𝑎) = −0.133
𝑇𝑟 = 12.04
𝛼 = 0.0929
𝑎 = 0.05587
𝑏 = 0.02315
𝐴∗ = 0.004172
𝐵∗ = 0.1096

𝑍 3 − 𝑍 2 − 0.1174 ∗ 𝑍 − 0.0001885 = 0
Resolviendo la ecuación cúbica llegamos a tres soluciones, de la cual una sola tendrá sentido
físico:

𝑍1 = 1.106
𝑍2 = −0.105
𝑍3 = −0.001628
El coeficiente de compresibilidad calculado para la mezcla es 1.106. A su vez, a través del
software REFPROP (utilizado en la industria petrolera), simulamos la mezcla y el programa
calculó el valor del coeficiente. El mismo dio como resultado 1.0939, muy cercano al calculado
mediante la ecuación de estado.
Como conclusión de este análisis, se puede observar una desviación de la idealidad de la mezcla,
debida a las altas presiones de trabajo. En el desarrollo del trabajo práctico modelaremos
nuestro sistema de reacción a partir de gases ideales, ya que simplifica la operatoria matemática,
pero poniendo atención a que esta suposición no es tan cierta.

También podría gustarte