1
UNIVERSIDAD NACIONAL DE TRUJILLO
FACULTAD DE INGENIERÍA QUÍMICA
PROGRAMA DE INGENIERÍA QUÍMICA
Proceso de Isomerización del n-Pentano
Anthony F. Canepa,
Mario S. Polo,
Jhorch J. Sarmiento,
Luis A. Reyes, Fernández
D. Valderrama,
Javier J. Villacorta,
Pierre A. Vidales (coordinador)
Termodinámica II
Ms. Rosa E. Nomberto
Trujillo – Perú
2022
2
Índice
Resumen ................................................................................................................................ 3
INTRODUCCIÓN: ............................................................................................................... 5
Proceso de Isomerización del n-Pentano............................................................................... 7
Descripción del proceso .................................................................................................... 7
Condiciones de operación ................................................................................................. 8
Tipo de Reactor ..................................................................................................................... 8
Balance equilibrio químico y equilibrio de fases. ................................................................. 8
Conclusiones ....................................................................................................................... 14
Referencias .......................................................................................................................... 15
Apéndice ………………………………………………………………………………… 16
3
Resumen
Se define isomerización como el proceso químico mediante el cual una molécula es
transformada en otra que posee los mismos átomos pero dispuestos de forma distinta. De este
modo, se dice que la primera molécula es un isómero de la segunda, y viceversa. En algunos
casos y para algunas moléculas, la isomerización puede suceder espontáneamente. De hecho,
algunos isómeros poseen aproximadamente la misma energía de enlace, lo que conduce a que
se presenten en cantidades más o menos iguales que se Inter convierten entre sí. La diferencia
de energía existente entre dos isómeros se denomina energía de isomerización.
El n-pentano es el componente principal de la gasolina Light Straight Run, y la
isomerización del n-pentano da como resultado un aumento significativo del índice de octano
(RON) de 61,7 a 92,3. Por lo tanto, es de gran importancia industrial estudiar los efectos del
tamaño de partícula de zeolita y los límites internos del grano sobre la isomerización de n-
pentano sobre Pt / Beta.
En la presente investigación tenemos como objetivo describir los parámetros
termodinámicos para hallar el equilibrio de fases en el proceso de isomerización del n-
pentano. Para ello determinamos los coeficientes de actividad utilizando los datos de
VL(cm3/mol) y δ (J/cm3)1/2. Para luego relacionar el valor de Ka con las composiciones de
ambas fases, mediante las fugacidades, la energía libre de Gibbs y temperaturas obteniendo la
relación entre el Ka y phi de la mezcla y los Y en ambas fases de la reacción.
Finalmente obtuvimos el equilibrio de fases en el proceso de isomerización del n-
pentano. Y partiendo de la regla de fases se pudo determinar un solo grado de libertad, de
modo que, una variable fija el estado del sistema.
Palabras clave: isomerización, zeolita, molécula
4
Abstract
Isomerization is defined as the chemical process by which one molecule is transformed
into another that has the same atoms but arranged differently. Thus, the first molecule is said
to be an isomer of the second, and vice versa. In some cases and for some molecules,
isomerization can happen spontaneously. In fact, some isomers have approximately the same
bond energy, which leads to their being present in more or less equal amounts that they
interconvert with each other. The difference in energy between two isomers is called the
isomerization energy.
n-pentane is the main component of Light Straight Run gasoline, and isomerization of
n-pentane results in a significant increase in octane number (RON) from 61.7 to 92.3.
Therefore, it is of great industrial importance to study the effects of zeolite particle size and
internal grain boundaries on the isomerization of n-pentane on Pt/Beta.
In the present investigation we aim to describe the thermodynamic parameters to find
the phase equilibrium in the isomerization process of n-pentane. For this we determined the
activity coefficients using the data of VL(cm3/mol) and δ (J/cm3)1/2. To then relate the value
of Ka with the compositions of both phases, through the fugacities, the Gibbs free energy and
temperatures, obtaining the relationship between the Ka and phi of the mixture and the Y in
both phases of the reaction.
Finally we obtained the phase equilibrium in the isomerization process of n-pentane.
And starting from the phase rule, it was possible to determine a single degree of freedom, so
that a variable fixes the state of the system.
Keywords: isomerization, zeolite, molecule
5
Introducción:
En el presente informe se presentará sobre la “Isomerización del n-pentano” en
como investigadores como Hu et al., han visto aplicaciones a gran escala de catalizadores está
actualmente restringida en la isomerización del n-pentano debido a su rápida desactivación y
alto costo aun así nos habla como los ácidos sólidos como la H-mordenita ha sido estudiado
en la isomerización de alcanos incluyendo al n-pentano, Ye et al., en el 2018, nos hablan que
el Pt/Beta es más activo y selectivo que otros para catalizar por lo que podría aumentar el
interés potencial del catalizador en uso industrial para la isomerización. En el caso Zhao et al.
en el 2020, nos habla como el n-pentano es uno de los mejores procedimientos para obtener
gas limpio de alta calidad y como el desarrollo de catalizadores para obtener este producto es
importante para la producción de combustibles limpios. Por parte de Ejtemaei et al, en el
2016, Se investigó adicionalmente de pentano a isopentano sobre el catalizador preparado en
condiciones atmosféricas. Una evaluación de las variables de reacción, incluida la
temperatura de reacción, la relación molar de hidrógeno a hidrocarburo y la velocidad
espacial, en la alimentación del norte se realizó la isomerización de pentano.
En el caso de Zhao, H et al. (2020) La isomerización de alcanos, como el n-pentano,
es uno de los mejores procedimientos para lograr una línea de gas limpio de alta calidad. Por
lo tanto, el desarrollo de catalizadores de alta actividad para la isomerización de alcanos es
muy importante para la producción de combustibles limpios. Así mismo se investigaron los
efectos de La y Ni en la estructura y el rendimiento de isomerización del catalizador, donde
los catalizadores LA-Ni-SZA-I y La-Ni-SZA-W fueron probados para la isomerización de n-
pentano a una presión de 2.0 MPa, consecuentemente obtuvieron mayores rendimientos para
dichos catalizadores a una temperatura de 220 °C y 280 °C, respectivamente. Se describen los
parámetros termodinámicos para hallar el equilibrio de fases en el proceso de isomerización
del n-pentano. Y partiendo de la regla de fases se pudo determinar un solo grado de libertad,
6
de modo que, una variable fija el estado del sistema, es decir la temperatura. Para su
modelado, Parputs et al., en 2014, de un proceso de n-pentano por destilación catalítica, se
nos es necesario es necesario comprender los parámetros cinéticos del proceso basado en un
catalizador ZnO2/SO4. Se realizó una serie de experimentos con un reactor de flujo clásico
para determinar las constantes de velocidad de reacción para la isomerización de n-pentano.
7
Proceso de Isomerización del n-Pentano
Descripción del proceso
En el presente trabajo, determinamos los coeficientes de actividad utilizando los datos
de VL(cm3/mol) y δ (J/cm3)1/2. Luego se relacionó el valor de Ka con las composiciones de
ambas fases, mediante las fugacidades, la energía libre de Gibbs y temperaturas
Se obtuvo la relación entre el Ka, y phi de mezcla y los Y en ambas fases de la
reacción.
Luego, mediante el equilibrio de fases se encontró la relación entre las fugacidades en
fase gaseosa y fase líquida, también una relación entre los valores de Yi, y los valores
correspondientes de 𝑃 𝑠𝑎𝑡 , ∅𝑠𝑎𝑡 y los coeficientes de actividad para ambas especies y
̂ 𝑠𝑎𝑡 utilizando la ecuación del virial, utilizando los datos de 𝐵0 , 𝐵1 y B,
posteriormente el ∅
Luego, se determinó 𝑃 𝑠𝑎𝑡 , utilizando datos de las constantes para la ecuación de
Antoine tanto para el n – pentano como para el isopentano.
Se calcularon los valores de XA y XB, asumiendo que 𝛾𝐴 y 𝛾𝐵 son iguales a 1. Con los
valores encontrados de XA y XB y las ecuaciones de determinación de los coeficientes de
actividad por método de soluciones regulares y datos de tabla, se determinaron los nuevos
valores de 𝛾𝐴 y 𝛾𝐵 . Posteriormente, con los datos de ∅𝐴𝑠𝑎𝑡 , ∅𝐵𝑠𝑎𝑡 , 𝑃𝐴𝑠𝑎𝑡 , 𝑃𝐵𝑠𝑎𝑡 𝛾𝐴 , 𝛾𝐵 , 𝑋𝐴 y 𝑋𝐵 , se
calculó el valor de Ka, esperando hallar un resultado igual al primero. Dado que no fue el
caso, se repitió el proceso anterior, hasta obtener el valor deseado. Con los datos de la última
iteración se evaluaron valores de XA y XB. Y con estos últimos datos se terminaron los
̂𝐴 y ∅
valores de YA y YB, asumiendo que ∅ ̂ 𝐵 son igual a 1.
Luego, se calculó el valor de la presión. Con el dato de presión y composición Y que
̂𝐴 y ∅
se obtuvo, se determinaron los valores de ∅ ̂ 𝐴 de mezcla, utilizando datos de tabla.
8
Condiciones de operación
Temperatura a 320 K
Tipo de Reactor
La isomerización del n-pentano se realiza en un equipo de flujo continuo con un
reactor tubular de lecho fijo en vidrio Pyrex, a presión atmosférica. Se utilizó H2 (en la
reacción de hidroisomerización) o una mezcla 1.5H2: N2 (en la reacción de
deshidroisomerización), como gas de arrastre manteniendo una determinada relación gas de
arrastre/ hidrocarburo. Para alcanzar esta relación molar, se hizo pasar el gas de arrastre a
través de un recipiente de vidrio (saturador) que contiene n-pentano sumergido en un baño de
hielo a 0ºC. Al pasar por el saturador el gas de arrastre se satura con el hidrocarburo
alcanzando una presión parcial igual a la presión de n-pentano como componente puro a 0ºC,
logrando la relación molar especificada para cada reacción. Este método es particularmente
conveniente para líquidos volátiles, ya que se logra saturar la corriente de gas de arrastre
fácilmente. (Lopez et al, 2011, p. 157)
Balance equilibrio químico y equilibrio de fases.
1. Regla de las fases de Gibbs
Se considera la regla de fases para el sistema de isomerización del n-pentano
La regla de las fases de Gibbs describe el número de grados de libertad (F) en un
sistema cerrado en equilibrio, en términos del número de fases separadas (P), el número de
componentes químicos (C) del sistema y 2 el número de variables no composicionales (por
ejemplo; presión o temperatura). (Josiah Willard., 1875).
9
Entonces se tiene:
F= C−M−P+2
Siendo:
F: grados de libertad
P: número de fases
C: componentes del sistema
M: reacciones independientes
Despejando en la ecuación de grados de libertad para la isomerización del n-
pentano, se tiene 2 fases, 2 componentes y 1 reacción de isomerización, se tiene:
F= 2−2−1+2
F=1
Por lo tanto, al tener solo un grado de libertad el n-pentano ejerce una presión
de descomposición fija a una temperatura fija.
2. Relación del Ka con las composiciones de ambas fases
𝑛 − 𝑝𝑒𝑛𝑡𝑎𝑛𝑜 ⇄ 𝐼𝑠𝑜𝑝𝑒𝑛𝑡𝑎𝑛𝑜
𝐴⇄𝐵
En una reacción:
𝑓̂𝑖
𝐾𝑎 = 𝜋(𝑎̂𝑖 )𝜈 ; = 𝑎̂𝑖
𝑓𝑜
𝑆𝑖 𝑓𝑖 𝑜 = 1 ⇒ 𝑎̂𝑖 = 𝑓̂𝑖
𝑓̂𝐵𝑉 𝑓̂𝐵𝐿
𝐾𝑎 = = … (2)
𝑓̂ 𝑉
𝐴 𝑓̂ 𝐿
𝐴
10
3. Se evaluó Ka a T del sistema de Rx; Ka a T= 320 K
°
∆𝐺(𝑇)
ln 𝐾𝑎 = − … (3)
𝑅𝑇
Usar la ecuación:
∆𝐺 ° = 𝐴 + 𝐵𝑇 + 𝐶𝑇 2 … (4)(T: °K)
∆𝐺 ° = ∆𝐴 + ∆𝐵𝑇 + ∆𝐶𝑇 2 … (5)
Se reemplazó en (5) los datos presentados en la tabla 4, adjunta al Apéndice
2, con una temperatura de 320 °K del sistema:
KJ KJ ∗ K KJ ∗ K 2
∆𝐺 = −8.304 + 0.00651(320𝐾) + 0.000001006(320K)2
gmol gmol gmol
KJ
∆𝐺 = 6.324
gmol
Reemplazar en la ecuación (3):
KJ
6.324
gmol
ln 𝐾𝑎 = −
𝐾𝐽
0.008314 ∗ 320𝐾
𝐾𝑔𝑚𝑜𝑙
𝐾𝑎 = 10.773
4. Se relacionó las fases líquida y vapor
Se reemplazó la ecuación (15) en (14) y se obtuvo para la constante de equilibrio:
𝛾𝐵 (1 − 𝑥𝐴 ) ∅𝐵𝑠𝑎𝑡 𝑃𝐵𝑠𝑎𝑡
𝐾𝑎 = … (14)
𝛾𝐴 𝑥𝐴 ∅𝐴𝑠𝑎𝑡 𝑃𝐴𝑠𝑎𝑡
11
̂ 𝐬𝐚𝐭 para ambas especies
5. Determinación de ∅
Usando datos de la tabla 1, se pudo determinar los valores de ∅̂ 𝒔𝒂𝒕 para ambas
especies que fueron expresando en la tabla 5, adjunta al Apéndice 2. Para ello se
utilizó las ecuaciones (16) a la (19) presentadas en el Apéndice 2.
Se determina 𝑃𝑠𝑎𝑡 para ambas especies con los datos de la tabla 2 adjunta al
Apéndice 1, y la ecuación de Antoine
Para el n-pentano:
1064.840
𝑃𝐴𝑠𝑎𝑡 = 10(3.97786−46.85+233.460)
𝑃𝐴𝑠𝑎𝑡 = 1.44 𝑏𝑎𝑟
Para el isopentano:
1022.880
𝑃𝐴𝑠𝑎𝑡 = 10(3.92023−46.85+232.014)
𝑃𝐴𝑠𝑎𝑡 = 1.87 𝑏𝑎𝑟
Determinar el valor de XA y XB con la ecuación (14) y (15) presentadas en el
Apéndice 2, suponiendo los 𝛾𝐴 y 𝛾𝐵 iguales a 1.
𝑥𝐴 =0.10650744
𝑥𝐵 =0.89349256
6. Determinación de los Coeficientes de Actividad para las composiciones de ambas
fases
Se determinó los nuevos valores de 𝛾𝐴 y 𝛾𝐵 con los valores de XA y XB
obtenidos en el paso 5, también con las ecuaciones de determinación de los
12
coeficientes de actividad por método de soluciones regulares y los valores de
volúmenes de líquido molar y solubilidad dadas en la tabla 3.
𝛾𝐴 = 1.0012569
𝛾𝐵 = 1.0000176
El despeje de estas ecuaciones se presentó ordenadamente en el Apéndice 3.
7. Determinación de Ka, las composiciones de ambas fases y la presión del sistema.
Se determina el nuevo valor de Ka con los datos de XA, XB, 𝛾𝐴 y 𝛾𝐵 obtenidos en
los pasos anteriores, además de los datos saturados en la ecuación (14):
𝛾𝐵 (1 − 𝑥𝐴 ) ∅𝐵𝑠𝑎𝑡 𝑃𝐵𝑠𝑎𝑡
𝐾𝑎 =
𝛾𝐴 𝑥𝐴 ∅𝐴𝑠𝑎𝑡 𝑃𝐴𝑠𝑎𝑡
(1.00002) (0.89375)(0.936793)(1.86693)
𝐾𝑎 =
(1.00126)(0.10625)(0.94624)(1.44334)
𝐾𝑎 = 10.7582
Se continuó realizando 2 iteraciones más, repitiendo el procedimiento para la
obtención de XA, XB, 𝛾𝐴 y 𝛾𝐵. Estas interacciones se pueden observar en la Tabla 6 en
el Apéndice 4. Con la última iteración se evaluó XA y XB, obteniendo:
XA = 0.10614
XB = 0.89386
Se determinó los valores de YA y YB con el despeje de la ecuación (10)
presentada en el apéndice 2, los valores obtenidos en la última iteración y los valores
de ∅̂𝐴 y ∅̂𝐵 que serán igual a 1, se obtuvo:
13
𝑦𝐴 = 0.08495
𝑦𝐵 = 0.91505
Luego, se determinó el valor de la presión despejando y reemplazando los datos
obtenidos en la ecuación (11).
𝛾𝐴 𝑥𝐴 ∅𝐴𝑠𝑎𝑡 𝑃𝐴𝑠𝑎𝑡 𝛾𝐴 𝑥𝐴 ∅𝐴𝑠𝑎𝑡 𝑃𝐴𝑠𝑎𝑡
𝑦𝐴 = → 𝑃=
̂𝐴 𝑃
∅ ̂ 𝐴 𝑦𝐴
∅
(1.00126)(0.10613518)(0.94624)(1.44334)
𝑃=
(1) (0.08495)
𝑃 = 1.70846 𝑏𝑎𝑟
Con este dato de presión y composición de Y, se determinaron los valores de ∅̂𝐴
y ∅̂𝐵 de mezcla, utilizando los datos de la Tabla 1 y las ecuaciones (21) hasta la (30).
Con los datos determinados de ∅̂ 𝐴 y ∅̂ 𝐵 presentados en el apéndice 3, se
vuelve a determinar YA,YB y la presión, realizando sucesivamente 5 iteraciones, las
cuales pueden observarse en la tabla 9, en el Apéndice 6.
En la quinta iteración, recién se asumió los valores de YA, YB, y Presión del
sistema, los resultados fueron:
YA=0.08532
YB=0.91468
Psistema=1.82237 bar
14
Conclusión
Como final de este trabajo de investigación formativa, donde hemos aplicado nuestros
conocimientos fundamentales y científico, paralelamente los principios termodinámicos que
nos ha permitido determinar las composiciones de equilibrio de las fases líquido y vapor, para
la isomerización del n-pentano teniendo la condición de operación (temperatura) a 320K,
donde hemos evaluado el número de variables independientes necesarias para poder
encontrar una solución, recurriendo a la regla de fase de Gibbs, por lo que se obtuvo 1 grado
de libertad, que esto quiere decir que la isomerización del n-pentano ejerce una presión de
descomposición fija a una temperatura fija.
Aplicando nuestros conocimientos hemos podido leer las propiedades termodinámicas
requeridas para este sistema. Es decir, esto nos ha permitido evaluar los coeficientes de
fugacidad, conociendo la temperatura critica, presión critica, factor acéntrico entre otro. (para
la fase gas) y también los parámetros para la fase liquida.
En conclusión, hemos podido encontrar el valor del coeficiente de equilibrio Ka que es
de 10.77, y juntamente también las composiciones como YA que tiene un valor de 0.08532 y
YB que tiene un valor de 0.91468, conociendo que la presión del sistema final es de 1.82237
bar.
Esta operación (isomerización n-pentano) se lleva a cabo en un equipo de flujo
continuo con un reactor tubular de lecho fijo. Donde en este presente informe damos a
conocer la descripción de este proceso, las condiciones de operación, y además el balance del
que lleva a realizar tanto el equilibrio químico como el equilibrio de fases.
15
Referencias
Ejtemaei, M., Charchi Aghdam, N., Babaluo, A., Tavakoli, A., & Bayati, B. (2017). n-
pentane isomerization over Pt-Al promoted sulfated zirconia nanocatalyst. Scientia
Iranica, 24(3), 1264–1271. https://doi.org/10.24200/sci.2017.4110
Hu, P., Zheng, J., Jiang, W., Zhong, L., & Zhou, S. (2020). Isomerization of n-pentane
catalyzed by amide-AlCl3-based ionic liquid analogs with various additives. Chinese
Journal of Chemical Engineering, 28(1), 152–157.
https://doi.org/10.1016/j.cjche.2019.04.018
López, C., Ortiz, C., Sazo, V., & García, L. (2011). Transformación de n-pentano sobre
catalizadores de platino soportado en tamices moleculares mesoporosos sba-15 con
aluminio. Revista de la Facultad de Ingeniería U.C.V, 26(3), 155–162.
https://www.researchgate.net/publication/255761477_Transformacion_de_n-
pentano_sobre_catalizadores_soportados_en_tamices_moleculares_mesoposorsos_
SBA-15_con_aluminio
Prausnitz, J. Poling, B & O’Connell J. (2001). The Properties of Gases and Liquids (Fifth
ed.).McGRAW-HILL. doi:10.1036/0070116822
Ye, G., Sun, Y., Guo, Z., Zhu, K., Liu, H., Zhou, X., & Coppens, M. O. (2018). Effects of
zeolite particle size and internal grain boundaries on Pt/Beta catalyzed isomerization
of n-pentane. Journal of Catalysis, 360, 152–159.
https://doi.org/10.1016/j.jcat.2018.01.033
Zhao, H., Song, H., Zhao, L., & Li, F. (2020). Influence of sulfating method on La–Ni–
S2O82–/ZrO2–Al2O3 solid superacid catalyst for n-pentane isomerization. Progress
in Reaction Kinetics and Mechanism, 45, 146867831989772.
https://doi.org/10.1177/1468678319897724
16
Apéndice
Apéndice 1: Tablas con datos y parámetros esenciales
Tabla 1
Constantes básicas para la determinación de coeficiente de fugacidad
Vc,
Especie Tc, K Pc, bar Zc w
cm3/mol
n-pentano 469.70 33.70 311.00 0.268 0.252
Isopentano 460.39 33.81 308.30 0.272 0.229
Nota. Datos extraídos de Poling, B. E., Prausnitz, J. M., & O'Connell, J. P. (2001)
Tabla 2
Constantes para la ecuación de Antoine
Especie A B C
n-pentano 3.97786 1064.840 232.014
Isopentano 3.92023 1022.880 233.460
Nota. Datos extraídos de Poling, B. E., Prausnitz, J. M., & O'Connell, J. P. (2001)
Tabla 3
Datos de volúmenes de líquido molar y solubilidad para la determinación de coeficientes
de actividad de mezclas binarias.
Especie VL(cm3/mol) δ (J/cm3)1/2
n-pentano 116 14.5
Isopentano 117 13.9
Nota. Datos extraídos de Poling, B. E., Prausnitz, J. M., & O'Connell, J. P. (2001, p.307).
17
̂ 𝒔𝒂𝒕
Apéndice 2: Datos y ecuaciones para la determinación Ka y ∅
Tabla 4
Datos de coeficientes para la determinación de la energía libre de Gibbs
Especie µ A B C
n-pentano -1 -149.141 0.45748 0.000044417
Isopentano 1 -157.445 0.46399 0.000043411
∆ 8.304 -0.00651 0.000001006
Nota. Datos extraídos de Yaws, C.L. & Chiang, P.Y.
Para determinar Ka primero relacionar ambas fases líquida y vapor, para el equilibrio de
fases:
𝑓̂𝑖𝐿 = 𝑓̂𝑖𝑉 … (6)
Despejando la ecuación (2) para la relación Ka de n-pentano (líquido) e isopentano
(gas), teniendo en cuenta las 2 fases en que se lleva a cabo la reacción:
𝑓̂𝐵𝑉
𝐾𝑎 = … (7)
𝑓̂𝐴𝐿
Del equilibrio de fases se conoce que:
𝑓̂𝑖𝐿 = 𝛾𝑖 𝑥𝑖 ∅𝑖𝑠𝑎𝑡 𝑃𝑖𝑠𝑎𝑡 … (8)
̂ 𝑖 𝑦𝑖 𝑃 … (9)
𝑓̂𝑖𝑉 = ∅
De la ecuación (6), se tiene:
̂ 𝑖 𝑦𝑖 𝑃 … (10)
𝛾𝑖 𝑥𝑖 ∅𝑖𝑠𝑎𝑡 𝑃𝑖𝑠𝑎𝑡 = ∅
De la ecuación (10), despejando “y” para el componente A y B.
18
𝛾𝐴 𝑥𝐴 ∅𝐴𝑠𝑎𝑡 𝑃𝐴𝑠𝑎𝑡
𝑦𝐴 = … (11)
̂𝐴 𝑃
∅
𝛾𝐵 𝑥𝐵 ∅𝐵𝑠𝑎𝑡 𝑃𝐵𝑠𝑎𝑡
𝑦𝐵 = … (12)
̂𝐵 𝑃
∅
Reemplazando la ecuación (9), en (7):
̂ 𝐵 𝑦𝐵 𝑃
∅
𝐾𝑎 =
̂ 𝐴 𝑦𝐴 𝑃
∅
̂ 𝐵 𝑦𝐵
∅
𝐾𝑎 = … (13)
̂ 𝐴 𝑦𝐴
∅
Además, reemplazando la ecuación (11) y (12), en (13):
̂ 𝐵 𝛾𝐵 𝑥𝐵 ∅𝐵𝑠𝑎𝑡 𝑃𝐵𝑠𝑎𝑡 (∅
∅ ̂ 𝐵 𝑃)
𝐾𝑎 =
̂ 𝐴 𝛾𝐴 𝑥𝐴 ∅𝐴𝑠𝑎𝑡 𝑃𝐴𝑠𝑎𝑡 (∅
∅ ̂ 𝐴 𝑃)
𝛾𝐵 𝑥𝐵 ∅𝐵𝑠𝑎𝑡 𝑃𝐵𝑠𝑎𝑡
𝐾𝑎 = … (14)
𝛾𝐴 𝑥𝐴 ∅𝐴𝑠𝑎𝑡 𝑃𝐴𝑠𝑎𝑡
Fase líquida, para nuestro sistema de 2 componentes:
∑ 𝑥𝑖 = 1 = 𝑥𝐴 + 𝑥𝐵
𝑥𝐵 = 1 − 𝑥𝐴 … (15)
̂ 𝑠𝑎𝑡 para ambas especies, utilizando las siguientes ecuaciones:
Se determinó ∅
0.422
𝐵0 = 0.083 − … (16)
𝑇𝑟 1.6
0.172
𝐵1 = 0.139 − … (17)
𝑇𝑟 4.2
𝑅 𝑇𝑐 0
𝐵= (𝐵 + 𝑤 𝐵1 ) … (18)
𝑃𝑐
19
𝐵 𝑃 𝑠𝑎𝑡
̂ 𝑠𝑎𝑡 =
𝑙𝑛∅ … (19)
𝑅 𝑇 𝑠𝑎𝑡
Tabla 5
Datos obtenidos para el cálculo coeficiente de fugacidad saturado
Especie 𝐁𝟎 … (𝟏𝟔) 𝑩𝟏 … (𝟏𝟕) 𝑩 … (𝟏𝟖) ̂ 𝒔𝒂𝒕 … (𝟏𝟗)
𝒍𝒏∅ ̂ 𝒔𝒂𝒕
∅
n-pentano -0.62562 -0.531534 -880.1417 -0.080406 0.922742
Isopentano -0.63168 -0.546686 -856.8357 -0.080285 0.922853
Suponer que 𝛾𝐴 y 𝛾𝐵 son iguales a 1 y determinar el valor de XA y XB con la ecuación
(14) y (15).
𝛾𝐵 𝑥𝐵 ∅𝐵𝑠𝑎𝑡 𝑃𝐵𝑠𝑎𝑡
𝐾𝑎 = … (14)
𝛾𝐴 𝑥𝐴 ∅𝐴𝑠𝑎𝑡 𝑃𝐴𝑠𝑎𝑡
Despejando para 𝑥𝐴 :
1
𝑥𝐴 =
𝑘 𝛾 ∅𝐴𝑠𝑎𝑡 𝑃𝐴𝑠𝑎𝑡
1 + ( 𝑎 𝐴 𝑠𝑎𝑡 )
𝛾𝐵 ∅𝐵 𝑃𝐵𝑠𝑎𝑡
𝑥𝐵 = 1 − 𝑥𝐴 … (15)
Apéndice 3: Ecuaciones para la Determinación de los Coeficientes de Actividad para las
composiciones de ambas fases
Primero reemplazamos en la ecuación:
𝑥𝑗 𝑉𝑗𝐿
𝛷𝑗 = … (𝑖)
∑𝑖 𝑥𝑖 𝑉𝑖𝐿
20
𝑥𝐴 𝑉𝐴𝐿
𝛷𝐴 =
(𝑥𝐴 𝑉𝐴𝐿 ) + (𝛿𝑗̅ 𝑉𝐵𝐿 )
𝑥𝐴 𝑉𝐴𝐿
𝛷𝐵 =
(𝑥𝐴 𝑉𝐴𝐿 ) + (𝛿𝑗̅ 𝑉𝐵𝐿 )
𝛷𝐴 = 0.1054404
𝛷𝐵 = 0.8945596
Reemplazar en la ecuación:
𝛿 ̅ = ∑ 𝛷𝑗 𝛿𝑗 … (𝑖𝑖𝑖)
𝑗
𝛿 ̅ = 𝛷𝐴 𝛿𝐴 + 𝛷𝐵 𝛿𝐵
𝛿 ̅ = 13.9632
Reemplazar en la ecuación:
(𝛿1 − 𝛿 ̅)2
ln 𝛾1 (𝑥1 ) = … (𝑖𝑖)
𝑅𝑇
𝑉𝐴𝐿 (𝛿𝐴 − 𝛿 ̅)2
ln 𝛾𝐴 =
𝑅𝑇
𝑉𝐵𝐿 (𝛿𝐵 − 𝛿 ̅)2
𝑙𝑛𝛾𝐵 =
𝑅𝑇
21
Apéndice 4: Interacciones para la determinación de nuestras variables
Tabla 6
Datos obtenidos en las iteraciones para valores de XA, XB, 𝛾𝐴 , 𝛾𝐵 y Ka
Iteración 1 Iteración 2
𝜸𝑨 1.0012572 1.00125725
𝜸𝑩 1.0000176 1.00001756
𝑿𝑨 0.1061352 0.10613518
𝑿𝑩 0.8938648 0.89386482
Ka 10.771496 10.771500
Apéndice 5: Ecuaciones para la determinación de los coeficientes de fugacidad de la
mezcla.
𝑃
̂𝐴 =
ln ∅ (𝐵 + 𝑦𝑗2 𝛿𝑖𝑗 ) … . (20)
𝑅𝑇 𝑖𝑖
𝑃
̂𝐵 =
ln ∅ (𝐵𝑗𝑗 + 𝑦𝑖2 𝛿𝑖𝑗 ) … . (21)
𝑅𝑇
𝛿𝑖𝑗 = 2𝐵𝑖𝑗 − 𝐵𝑖𝑖 − 𝐵𝑗𝑗 … (22)
𝑅𝑇𝑐𝑖𝑗 0
𝐵𝑖𝑗 = (𝐵 + 𝑤𝑖𝑗 𝐵1 ) …. (23)
𝑃𝑐𝑖𝑗
0.422
𝐵 0 = 0.083 − … (24)
𝑇𝑟 1.6
0.172
𝐵1 = 0.139 − … . (25)
𝑇𝑟 4.2
22
𝑤𝑖 + 𝑤𝑗
𝑤𝑖𝑗 = … (26)
2
𝑍𝑐𝑖 + 𝑍𝑐𝑗
𝑍𝑐𝑖𝑗 = … . (27)
2
1
𝑇𝑐𝑖𝑗 = (𝑇𝑐𝑖 × 𝑇𝑐𝑗 )2 (1 − 𝐾𝑖𝑗 ) … . (28)
1 1 3
𝑉𝑐𝑖 3 + 𝑉𝑐𝑗 3
𝑉𝑐𝑖𝑗 = [ ] … . (29)
2
𝑍𝑐𝑖𝑗 × 𝑅𝑇𝑐𝑖𝑗
𝑃𝑐𝑖𝑗 = … (30)
𝑉𝑐𝑖𝑗
Apéndice 6: Obtención de datos del trabajo
Tabla 7
Datos cálculos para la mezcla de n-pentano con isopentano.
𝒘𝒊𝒋 . . . (𝟐𝟔) 𝒁𝒄𝒊𝒋 … (𝟐𝟕) 𝑻𝒄𝒊𝒋 … (𝟐𝟖) 𝑽𝒄𝒊𝒋 … (𝟐𝟗) 𝑷𝒄𝒊𝒋 … (𝟑𝟎) 𝑩𝟎 … (𝟐𝟒) 𝑩𝟏 … (𝟐𝟓)
0.2405 0.27 465.0217 309.64804 33.710395 -0.684414 -
0.687582
Tabla 8
Datos cálculos para la mezcla y determinación de los coeficientes de fugacidad.
𝑩𝒊𝒋 … (𝟐𝟑) 𝜹𝒊𝒋 … (𝟐𝟐) ̂ 𝑨 … (𝟐𝟎)
𝐥𝐧 ∅ ̂ 𝑩 … (𝟐𝟏)
𝐥𝐧 ∅ ̂𝑨
∅ ̂𝑩
∅
-974.5626 -212.1477 -0.0679288 -0.055123 0.934327 0.94636866
23
Tabla 9
̂𝐴, ∅
Datos obtenidos de las iteraciones para YA, YB, ∅ ̂𝐵 y P
Iteración Iteración Iteración Iteración Iteración
1 2 3 4 5
̂𝑨
∅ 0.93813 0.93373 0.93351 0.93349 0.93349
̂𝑩
∅ 0.94074 0.93801 0.93788 0.93787 0.93787
YA 0.08517 0.08531 0.08531 0.08532 0.08532
YB 0.91483 0.91469 0.91469 0.91468 0.91468
P 1.81652 1.82207 1.82236 1.82237 1.82237
Apéndice 7: Enlace al Vídeo de exposición
https://drive.google.com/file/d/1RDqc2IomyRGmWHG0naqFTlbY9QJyOqJS/view?usp
=sharing