UNIFAC

UNIFAC^[1] (Abkürzung für Universal Quasichemical Functional Group Activity Coefficients) ist ein Verfahren zur Abschätzung von Aktivitätskoeffizienten, das in der Verfahrenstechnik und der Technische Chemie häufig eingesetzt wird.

Prinzipien

UNIFAC ist eine Gruppenbeitragsmethode, die auf dem Prinzip der Mischung aus Strukturgruppen basiert. Dies steht im Gegensatz zu der normalen Betrachtungsweise der Mischung aus Molekülen. Die Gruppen sind dabei üblicherweise funktionelle Gruppen, bspw. Alkohol- oder Carbonylgruppe, aber auch kleinere Molekülfragmente wie einzelne Atome, jedoch werden diese kleinen Fragmente stets unter Berücksichtigung ihrer chemischen Umgebung betrachtet. Einige wenige kleinere Moleküle (etwa Wasser) werden vollständig als eigenständige Gruppe definiert.

Gruppenbeitragsmethodenprinzip

Das Prinzip der Gruppenbeitragsmethoden besteht im Wesentlichen daraus, dass nur die Eigenschaften einiger Dutzend Strukturgruppen bekannt sein müssen, statt die Eigenschaften von einigen Millionen Stoffen für die Berechnung von Stoffgemischeigenschaften zu benötigen. Mit diesen wenigen Strukturgruppen als Bausteinen kann eine sehr große Anzahl von Molekülen konstruiert werden.

UNIFAC-Gruppenbeiträge

UNIFAC benutzt zur Vorhersage zwei Arten von Gruppenbeiträgen:

Additive Beiträge, die sehr kleinen Gruppen (Untergruppen) zugeordnet werden. Dies sind Gruppenvolumina und Gruppenoberflächen.
Wechselwirkungsparameter zwischen größeren Gruppen, die mehrere ähnliche Untergruppen umfassen.

Die Gruppenvolumina- und oberflächen basieren auf angenäherten van der Waalsschen Oberflächen und Volumina^[2] und sind somit Konstanten im Modell, die einen physikalisch-chemischen Hintergrund haben.

Die Wechselwirkungsparameter werden an experimentelle Aktivitätskoeffizienten sowie an Phasengleichgewichtsdaten, aus denen Aktivitätskoeffizienten abgeleitet werden können, mit nichtlinearen Optimierverfahren angepasst. Die Wechselwirkungsparameter werden somit ausschließlich empirisch ermittelt.

Gleichungen

Hier wird auf eine vollständige Aufführung der recht komplexen Gleichungen verzichtet, sondern nur versucht, mit vereinfachten Formulierungen einige Grundprinzipien zu verdeutlichen. Ausführliche Formulierungen können im Originalartikel gefunden werden oder auch in Lehrbüchern^[3].

Die Berechnung des Aktivitätskoeffizienten erfolgt additiv über zwei Terme:

$ln\;\gamma _{i}=ln\;\gamma _{i}^{C}+ln\;\gamma _{i}^{R}$

Kombinatorischer Anteil

γ^C wird als kombinatorischer Anteil bezeichnet (der Gruppe i) und wird aus den van der Waalsschen Oberflächen Q und Volumina R berechnet:

$ln\gamma _{i}^{C}=f\left(\sum \nu _{k}\cdot R_{k},\sum \nu _{k}\cdot Q_{k},x\right)$ mit ν als Gruppenhäufigkeit und x als Molenbrüche der Komponenten, aus denen die Mischung besteht.

Restanteil

γ^R wird als Restanteil bezeichnet (der Gruppe i) und wird aus angepassten Wechselwirkungsparametern berechnet:

$ln\;\gamma ^{R}=f\left(\psi ,x,X\right)$ mit X als Gruppenmolenbruch, x als Komponentenmolenbruch und ψ als

$\psi =exp\left({\frac {-a_{nm}}{T}}\right)$ mit a_nm als Wechselwirkungsparameter.

Geschichte

UNIFAC wurde in den 1970er Jahren entwickelt, wobei der Schwerpunkt zu Beginn alleine auf die Vorhersage von Dampf-Flüssig-Gleichgewichten von Mischungen einfacher, im wesentlichen organischer Substanzen und von Wassergemischen gelegt wurde. Da für die Modellentwicklung Aktivitätskoeffizienten für eine Vielzahl von Gemischen benötigt wurden, wurde mit dem Aufbau einer Faktendatenbank (der Dortmunder Datenbank) begonnen, die noch heute existiert, wenn auch in stark erweiterter und veränderter Form.

Da Aktivitätskoeffizienten es auch erlauben, mit einfachen thermodynamischen Beziehungen Fest-Flüssig-Gleichgewichte (Löslichkeiten von Feststoffen in Flüssigkeiten) und Flüssig-Flüssig-Gleichgewichte (siehe hierzu Mischungslücke) zu berechnen, wurde das UNIFAC-Modell in den achtziger Jahren auch zunehmend für die Berechnung dieser Phasengleichgewichte verwendet. Die dabei zutage tretenden Modellschwächen des Originalmodells führte zum einen zur Entwicklung von speziell parametrisierten Modellen für die Berechnung von Flüssig-Flüssig-Gleichgewichten (UNIFAC-LLE) und später zu einem Modell für die Abschätzung von Octanol-Wasser-Verteilungskoeffizienten^[4]. Zum anderen wurde das Modell selbst erweitert, bspw. um temperaturabhängige Wechselwirkungsparameter:

$\psi =exp\left({{\frac {-a_{nm}}{T}}+{\frac {-b_{nm}\cdot T}{T}}+{\frac {-c_{nm}\cdot T^{2}}{T}}}\right)$

Diese erweiterten Modelle werden als modifizierte UNIFAC-Modelle^[5] bezeichnet.

Bedeutung

UNIFAC (insbesondere die neueren Varianten) ist heute das am häufigsten verwendete Verfahren zur Abschätzung von Stoffgemischdaten. Die Vorhersagen werden bspw. in der Prozessimulation verwendet, dem heute vorherrschenden Verfahren zur Auslegung und Optimierung chemischer Prozesse, Anlagen und ganzer Fabriken. UNIFAC wird außerdem bei der Prozesssynthese verwendet, bei dem, sehr allgemein formuliert, Stoffe mit bestimmten Eigenschaften für bestimmte Aufgaben gesucht werden. Diese Aufgabe kann bspw. die eines Schleppmittels für die Azeotrop- oder Extraktivrektifikation sein.

Aktuelle Entwicklungen

UNIFAC wird in einer Reihe von Forschungsgruppen genutzt und weiterentwickelt. Aktuelle Entwicklungen (Auswahl) zielen auf

die Berechnung von Elektrolyt-haltigen Mischungen^[6]
die Abschätzung von Viskositäten^[7]
die Integration von UNIFAC in Mischungsregeln für Zustandsgleichungen^[8]
die Vorhersage von UNIFAC-Wechselwirkungsparametern^[9]
die Erweiterung des UNIFAC-Modells auf spezielle Stoffgruppen^[10]

Literatur

↑ Å. Fredenslund, J. Gmehling, P. Rasmussen; „Vapor-Liquid Equilibria Using UNIFAC - A Group Contribution Method;“ Elsevier, Amsterdam (1977)
↑ Bondi A., "van der Waals Volumes and Radii", J.Phys.Chem., 68(3), 441-451, 1964
↑ J. Gmehling, B. Kolbe; „Thermodynamik“, 2. Auflage; VCH-Verlag, Weinheim (1992)
↑ Gudrun Wienke, „Messung und Vorausberechnung von n-Octanol/Wasser-Verteilungskoeffizienten", Doktorarbeit, Univ. Oldenburg, 1-172, 1993
↑ Weidlich U., Gmehling J., „A Modified UNIFAC Model. 1. Prediction of VLE, hE, and γ^∞“", Ind.Eng.Chem.Res., 26(7), 1372-1381, 1987
↑ Polka H.-M., „Experimentelle Bestimmung und Berechnung von Dampf-Flüssig-Gleichgewichten für Systeme mit starken Elektrolyten", Doktorarbeit, Univ. Oldenburg, 1-144, 1993
↑ Gaston-Bonhomme Y., Petrino P., Chevalier J.L., „UNIFAC-VISCO Group Contribution Method for Predicting Kinematic Viscosity: Extension and Temperature Dependence“, Chem.Eng.Sci., 49(11), 1799-1806, 1994
↑ Holderbaum T., „Die Vorausberechnung von Dampf-Flüssig-Gleichgewichten mit einer Gruppenbeitragszustandsgleichung“, Fortschrittsber.VDI Reihe 3, 243, 1-154, 1991
↑ H. E. Gonzàlez, J. Abildskov, R. Gani, P. Rousseaux, B. Le Bert, „A method for prediction of UNIFAC group interaction parameters“, AIChE J., 53, 1620-1633, 2007
↑ Zhang Y., Zhang X., Zhang W., Qu H., Wang W., „Prediction of Solid-Liquid Equilibrium to Synthetic Nitro-Musk by using UNIFAC Group-Contribution Method“, Huaxue-gongcheng, 33(5), 69-71, 2005

Weblinks

UNIFAC-Industriekonsortium an der Carl von Ossietzky Universität Oldenburg

[1] Å. Fredenslund, J. Gmehling, P. Rasmussen; „Vapor-Liquid Equilibria Using UNIFAC - A Group Contribution Method;“ Elsevier, Amsterdam (1977)

[2] Bondi A., "van der Waals Volumes and Radii", J.Phys.Chem., 68(3), 441-451, 1964

[3] J. Gmehling, B. Kolbe; „Thermodynamik“, 2. Auflage; VCH-Verlag, Weinheim (1992)

[4] Gudrun Wienke, „Messung und Vorausberechnung von n-Octanol/Wasser-Verteilungskoeffizienten", Doktorarbeit, Univ. Oldenburg, 1-172, 1993

[5] Weidlich U., Gmehling J., „A Modified UNIFAC Model. 1. Prediction of VLE, hE, and γ^∞“", Ind.Eng.Chem.Res., 26(7), 1372-1381, 1987

[6] Polka H.-M., „Experimentelle Bestimmung und Berechnung von Dampf-Flüssig-Gleichgewichten für Systeme mit starken Elektrolyten", Doktorarbeit, Univ. Oldenburg, 1-144, 1993

[7] Gaston-Bonhomme Y., Petrino P., Chevalier J.L., „UNIFAC-VISCO Group Contribution Method for Predicting Kinematic Viscosity: Extension and Temperature Dependence“, Chem.Eng.Sci., 49(11), 1799-1806, 1994

[8] Holderbaum T., „Die Vorausberechnung von Dampf-Flüssig-Gleichgewichten mit einer Gruppenbeitragszustandsgleichung“, Fortschrittsber.VDI Reihe 3, 243, 1-154, 1991

[9] H. E. Gonzàlez, J. Abildskov, R. Gani, P. Rousseaux, B. Le Bert, „A method for prediction of UNIFAC group interaction parameters“, AIChE J., 53, 1620-1633, 2007

[10] Zhang Y., Zhang X., Zhang W., Qu H., Wang W., „Prediction of Solid-Liquid Equilibrium to Synthetic Nitro-Musk by using UNIFAC Group-Contribution Method“, Huaxue-gongcheng, 33(5), 69-71, 2005

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]