Articulo de referencia

UNIFAC

Azeótropo en silla de montar calculado con UNIFAC a 1 atm. Las líneas rojas representan las composiciones de vapor y las líneas azules, las composiciones de líquido. La imagen g...

Azeótropo en silla de montar calculado con UNIFAC a 1 atm. Las líneas rojas representan las composiciones de vapor y las líneas azules, las composiciones de líquido. La imagen gira para mostrar con mayor claridad la forma de silla de montar de los equilibrios vapor-líquido.

En termodinámica estadística , el método UNIFAC ( coeficientes de actividad de grupos funcionales UNI QUAC ) [ 1 ] es un sistema semiempírico para la predicción de la actividad de no electrolitos en mezclas no ideales . UNIFAC utiliza los grupos funcionales presentes en las moléculas que componen la mezcla líquida para calcular los coeficientes de actividad . Mediante el uso de las interacciones de cada uno de los grupos funcionales presentes en las moléculas, así como algunos coeficientes de interacción binaria, se puede calcular la actividad de cada una de las soluciones. Esta información se puede utilizar para obtener información sobre equilibrios líquidos, lo cual es útil en muchos cálculos termodinámicos, como el diseño de reactores químicos y los cálculos de destilación .

El modelo UNIFAC fue publicado por primera vez en 1975 por Fredenslund, Jones y John Prausnitz , un grupo de investigadores de ingeniería química de la Universidad de California . Posteriormente, ellos y otros autores han publicado numerosos artículos sobre UNIFAC, ampliando las capacidades del modelo mediante el desarrollo de nuevos parámetros o la revisión de los existentes. UNIFAC representa un intento de estos investigadores por proporcionar un modelo flexible de equilibrio líquido-líquido para un uso más amplio en química y en las disciplinas de ingeniería química y de procesos .

Introducción

Un problema particular en el campo de la termodinámica del estado líquido es la obtención de constantes termodinámicas fiables. Estas constantes son necesarias para predecir con éxito el estado de energía libre del sistema; sin esta información, es imposible modelar las fases de equilibrio del mismo.

Obtener estos datos de energía libre no es tarea fácil y requiere experimentos minuciosos, como la calorimetría , para medir con precisión la energía del sistema. Incluso cuando se realiza este trabajo, resulta inviable intentar llevarlo a cabo para cada clase posible de compuestos químicos y sus mezclas binarias o de mayor complejidad. Para paliar este problema, se emplean modelos de predicción de energía libre, como UNIFAC, que predicen la energía del sistema a partir de unas pocas constantes medidas previamente.

Es posible calcular algunos de estos parámetros mediante métodos ab initio como COSMO-RS , pero los resultados deben interpretarse con precaución, ya que las predicciones ab initio pueden ser inexactas. De igual modo, UNIFAC puede ser inexacto, y para ambos métodos es recomendable validar experimentalmente las energías obtenidas a partir de estos cálculos.

Correlación UNIFAC

La correlación UNIFAC intenta simplificar la predicción de interacciones entre moléculas describiendo dichas interacciones en función de los grupos funcionales unidos a la molécula. Esto se hace para reducir la cantidad de interacciones binarias que sería necesario medir para predecir el estado del sistema.

Actividad química

El coeficiente de actividad de los componentes de un sistema es un factor de corrección que tiene en cuenta las desviaciones de los sistemas reales respecto a una solución ideal , que puede medirse experimentalmente o estimarse a partir de modelos químicos (como UNIFAC). Al añadir un factor de corrección, conocido como la actividad (ai{\displaystyle a_{i}}Al relacionar la actividad del componente i con la fracción de fase líquida de una mezcla líquida, se pueden explicar algunos de los efectos de la solución real. La actividad de una sustancia química real es función del estado termodinámico del sistema, es decir, la temperatura y la presión.

Conociendo los coeficientes de actividad y sus componentes y proporciones relativas, se pueden calcular fenómenos como la separación de fases y los equilibrios vapor-líquido . UNIFAC pretende ser un modelo general para la predicción precisa de los coeficientes de actividad.

Parámetros del modelo

El modelo UNIFAC divide el coeficiente de actividad para cada especie en el sistema en dos componentes; un coeficiente combinatorioγdo{\displaystyle \gamma ^{c}}y un componente residualγr{\displaystyle \gamma ^{r}}. Para eli{\displaystyle i}-molécula, los coeficientes de actividad se desglosan según la siguiente ecuación:

lnγi=lnγido+lnγir.{\displaystyle \ln \gamma _{i}=\ln \gamma _{i}^{c}+\ln \gamma _{i}^{r}.}

En el modelo UNIFAC, se requieren tres parámetros principales para determinar la actividad de cada molécula en el sistema. En primer lugar, está el área de superficie del grupo.Q{\displaystyle Q}y contribuciones de volumenR{\displaystyle R}obtenidos a partir del área superficial y los volúmenes de Van der Waals . Estos parámetros dependen exclusivamente de los grupos funcionales individuales en las moléculas anfitrionas. Finalmente, está el parámetro de interacción binaria.τij{\displaystyle \tau _{ij}}, que está relacionado con la energía de interacciónUi{\displaystyle U_{i}}de pares moleculares (ecuación en la sección "residual"). Estos parámetros deben obtenerse mediante experimentos, ajuste de datos o simulación molecular.

Combinacional

El componente combinatorio de la actividad se ve influenciado por varios términos en su ecuación (a continuación), y es el mismo que para el modelo UNIQUAC .

lnγido=lnϕiincógnitai+z2qilnθiϕi+Liϕiincógnitaij=1norteincógnitajLj,{\displaystyle \ln \gamma _{i}^{c}=\ln {\frac {\phi _{i}}{x_{i}}}+{\frac {z}{2}}q_{i}\ln {\frac {\theta _{i}}{\phi _{i}}}+L_{i}-{\frac {\phi _{i}}{x_{i}}}\displaystyle \sum _{j=1}^{n}x_{j}L_{j},}

dóndeθi{\displaystyle \theta _{i}}yϕi{\displaystyle \phi _{i}}son los componentes fraccionarios de área y segmento ponderado molar para eli{\displaystyle i}-ésima molécula en el sistema total y se definen mediante la siguiente ecuación;Li{\displaystyle L_{i}}es un parámetro compuesto der{\displaystyle r},z{\displaystyle z}yq{\displaystyle q}.z{\displaystyle z}es el número de coordinación del sistema, pero se ha descubierto que el modelo es relativamente insensible a su valor y con frecuencia se cita como una constante con un valor de 10.

θi=incógnitaiqij=1norteincógnitajqj,ϕi=incógnitairij=1norteincógnitajrj,Li=z2(riqi)(ri1),z=10,{\displaystyle \theta _{i}={\frac {x_{i}q_{i}}{\displaystyle \sum _{j=1}^{n}x_{j}q_{j}}},\quad \phi _{i}={\frac {x_{i}r_{i}}{\displaystyle \sum _{j=1}^{n}x_{j}r_{j}}},\quad L_{i}={\frac {z}{2}}(r_{i}-q_{i})-(r_{i}-1),\quad z=10,}

qi{\displaystyle q_{i}}yri{\displaystyle r_{i}}se calculan a partir de las contribuciones de área superficial y volumen del grupoQ{\displaystyle Q}yR{\displaystyle R}(Generalmente se obtiene a través de valores tabulados), así como el número de ocurrencias del grupo funcional en cada molécula.νk{\displaystyle \nu _{k}}de tal manera que:

ri=k=1norteνkRk,qi=k=1norteνkQk.{\displaystyle r_{i}=\displaystyle \sum _{k=1}^{n}\nu _{k}R_{k},\quad q_{i}=\displaystyle \sum _{k=1}^{n}\nu _{k}Q_{k}.}

Residual

El componente residual de la actividadγr{\displaystyle \gamma ^{r}}se debe a las interacciones entre los grupos presentes en el sistema, y ​​el artículo original hace referencia al concepto de "solución de grupos". El componente residual de la actividad para eli{\displaystyle i}-ésima molécula que contienenorte{\displaystyle n}Los grupos funcionales únicos se pueden escribir de la siguiente manera:

lnγir=knorteνk(i)[lnΓklnΓk(i)],{\displaystyle \ln \gamma _{i}^{r}=\displaystyle \sum _{k}^{n}\nu _{k}^{(i)}\left[\ln \Gamma _{k}-\ln \Gamma _{k}^{(i)}\right],}

dóndeΓk(i){\displaystyle \Gamma _{k}^{(i)}}es la actividad de un grupo aislado en una solución que consiste únicamente en moléculas de tipoi{\displaystyle i}. La formulación de la actividad residual asegura que la condición para el caso límite de una sola molécula en una solución de componente puro, la actividad es igual a 1; como por definición deΓk(i){\displaystyle \Gamma _{k}^{(i)}}, uno encuentra quelnΓklnΓk(i){\displaystyle \ln \Gamma _ {k}-\ln \Gamma _ {k}^{(i)}}será cero. La siguiente fórmula se utiliza para ambosΓk{\displaystyle \Gamma _{k}}yΓk(i){\displaystyle \Gamma _{k}^{(i)}}

lnΓk=Qk[1lnmetroΘmetroΨmetrokmetroΘmetroΨkmetronorteΘnorteΨnortemetro].{\displaystyle \ln \Gamma _{k}=Q_{k}\left[1-\ln \displaystyle \sum _{m}\Theta _{m}\Psi _{mk}-\displaystyle \sum _{m}{\frac {\Theta _{m}\Psi _{km}}{\displaystyle \sum _{n}\Theta _{n}\Psi _{nm}}}\right].}

En esta fórmulaΘmetro{\displaystyle \Theta _{m}}es la suma de la fracción de área del grupometro{\displaystyle m}, sobre todos los diferentes grupos y es algo similar en forma, pero no es lo mismo queθi{\displaystyle \theta _{i}}.Ψmetronorte{\displaystyle \Psi _{mn}}es el parámetro de interacción de grupo y es una medida de la energía de interacción entre grupos. Se calcula mediante una ecuación de Arrhenius (aunque con una pseudoconstante de valor 1).incógnitanorte{\displaystyle X_{n}}es la fracción molar del grupo , que es el número de gruposnorte{\displaystyle n}en la solución dividida por el número total de grupos.

Θmetro=QmetroincógnitametronorteQnorteincógnitanorte,{\displaystyle \Theta _{m}={\frac {Q_{m}X_{m}}{\displaystyle \sum _{n}Q_{n}X_{n}}},}
Ψmetronorte=exp[UmetronorteUnortemetroRT],incógnitametro=jνmetrojincógnitajjnorteνnortejincógnitaj,{\displaystyle \Psi _{mn}=\exp \left[-{\frac {U_{mn}-U_{nm}}{RT}}\right],\quad X_{m}={\frac {\displaystyle \sum _{j}\nu _{m}^{j}x_{j}}{\displaystyle \sum _{j}\displaystyle \sum _{n}\nu _{n}^{j}x_{j}}},}

Umetronorte{\displaystyle U_{mn}}es la energía de interacción entre los grupos m y n , con unidades del SI de julios por mol y R es la constante de los gases ideales . Nótese que no es el caso que Umetronorte=Unortemetro{\displaystyle U_{mn}=U_{nm}}lo que da lugar a un parámetro no reflexivo. La ecuación para el parámetro de interacción grupal se puede simplificar a lo siguiente:

Ψmetronorte=expametronorteT.{\displaystyle \Psi _{mn}=\exp {\frac {-a_{mn}}{T}}.}

De este modoametronorte{\displaystyle a_{mn}}aún representa la energía neta de interacción entre gruposmetro{\displaystyle m}ynorte{\displaystyle n}pero utiliza unidades de temperatura absoluta algo inusuales ( kelvin del SI ). Estos valores de energía de interacción se obtienen a partir de datos experimentales y suelen estar tabulados.

Véase también

Referencias

  1. Aage Fredenslund, Russell L. Jones y John M. Prausnitz, "Estimación de la contribución de grupos de los coeficientes de actividad en mezclas líquidas no ideales", AIChE Journal , vol. 21 (1975), pág. 1086

Lecturas adicionales

  • Grupos estructurales y parámetros de UNIFAC
  • AIOMFAC, modelo en línea basado en UNIFAC, es un modelo de contribución de grupos para el cálculo de coeficientes de actividad en mezclas orgánicas - inorgánicas.
Obtenido de " https://en.wikipedia.org/w/index.php?title=UNIFAC&oldid=1326486222 "