Análisis Del Factor de Compresibilidad de Los Gases
Análisis Del Factor de Compresibilidad de Los Gases
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.
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.
𝑇𝑐(𝑚𝑒𝑧𝑐𝑙𝑎) = ∑ 𝑦𝑖 ∗ 𝑇𝑐𝑖
𝑃𝑐(𝑚𝑒𝑧𝑐𝑙𝑎) = ∑ 𝑦𝑖 ∗ 𝑃𝑐𝑖
𝑤(𝑚𝑒𝑧𝑐𝑙𝑎) = ∑ 𝑦𝑖 ∗ 𝑤𝑖
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:
𝑍 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.
𝑇𝑐(𝑚𝑒𝑧𝑐𝑙𝑎) = 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.