UNIFAC[1] ist ein Verfahren zur Abschätzung von Aktivitätskoeffizienten.

Prinzipien

Bearbeiten

UNIFAC (Abkürzung für Universal Quasichemical Functional Group Activity Coefficients) 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 Fragmente 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

Bearbeiten
 
Beispiel einer Gruppenzuordnung in UNIFAC

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

UNIFAC-Gruppenbeiträge

Bearbeiten

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 (Dampf-Flüssig, Flüssig-Flüssig, Fest-Flüssig, Mischungswärmen), aus denen Aktivitätskoeffizienten abgeleitet werden können, mit nicht-linearen Optimierverfahren angepasst. Die Wechselwirkungsparameter werden somit ausschließlich empirisch ermittelt.

Gleichungen

Bearbeiten

Ausführliche Formulierungen können im Originalartikel gefunden werden oder auch in Lehrbüchern[3][4].

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

 

Kombinatorischer Anteil
Bearbeiten

γC wird als kombinatorischer Anteil bezeichnet und wird aus van der Waalsschen Oberflächen (F) und Volumina (V) sowie Molenbrüchen (x) berechnet.

 

mit

 

 

Neben den Molenbrüchen werden van der Waalsche Oberflächen qi und Volumina ri der Moleküle benötigt. Diese lassen sich aus tabellierten Werten für die Gruppen (Gruppenoberfläche Qk und Gruppenvolumen Rk) bestimmen.

 

 

νk(i) ist die Häufigkeit der Gruppe k im Molekül i.

Restanteil
Bearbeiten

γR wird als Restanteil bezeichnet und wird schlussendlich aus angepassten Wechselwirkungsparametern berechnet.

Der Restanteil wird aus Gruppenaktivitätskoeffizienten Γk berechnet.

 

Die Gruppenaktivitätskoeffizienten in der Mischung Γk und im Reinstoff Γk(i) werden über die Beziehung

 

berechnet.

Θ ist der Oberflächenanteil

 

und X der Gruppenmolenbruch.

 

Der Parameter ψ enthält die anpassbaren Parameter des UNIFAC-Modells.

  mit anm als Wechselwirkungsparameter zwischen den Gruppen.

Geschichte

Bearbeiten

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, wenn auch in stark erweiterter und veränderter Form, noch existiert.

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) 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 zur Abschätzung von Octanol-Wasser-Verteilungskoeffizient[5]. Zum anderen wurde das Modell erweitert, bspw. um temperaturabhängige Wechselwirkungsparameter:

 

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

Bedeutung

Bearbeiten

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 [[Azeotroprektifikation|Azeotrop-] oder Extraktivrektifikation sein.

Aktuelle Parametrisierung

Bearbeiten
 
UNIFAC-Parametermatrix (letzte veröffentlichte Version[7])

Das UNIFAC-Modell erlaubt die Vorhersage von Aktivitätskoeffizienten in Gemischen mit Alkanen, Alkenen, Alkinen, Alkoholen, Aromaten, Estern, Ethern, Aminen, Carbonsäuren, organischen Fluoriden, organischen Chloriden, organischen Bromiden, organischen Iodiden, Thiolen, Sulfonen, Wasser, Furfural, Thiophenen, Pyridinen, Morpholin, Isocyanaten, Silanen, Siloxanen, Amiden und organischen Nitraten.

Für diese funktionellen Gruppen sind Wechselwirkungsparameter mit mindestens einer der anderen Gruppen vorhanden. Es sind jedoch nicht alle Kombinationen vollständig parametrisiert, bspw. kann ein Aktivitätskoeffizient eines Gemisch aus 3-Methylthiophen und n-Hexan berechnet werden, eine Vorhersage eines Aktivitätskoeffizienten in einem Gemisch des Thiophens mit 3-Hexen scheitert jedoch an fehlenden Parametern.

UNIFAC kann nicht verwendet werden, um Aktivitätskoeffizienten in Gemischen zu berechnen, die Salze bzw., allgemeiner, Elektrolyte oder Polymere enthalten. Für diese beiden Stoffklassen gibt es aber Erweiterungen des hier beschriebenen UNIFAC-Modells.

Aktuelle Entwicklungen

Bearbeiten

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

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

Literatur

Bearbeiten
  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. J. Vidal, „Thermodynamique. Méthodes appliquées au raffinage et au Génie Chimique“, Editions Technip, 1997
  5. Gudrun Wienke, „Messung und Vorausberechnung von n-Octanol/Wasser-Verteilungskoeffizienten", Doktorarbeit, Univ. Oldenburg, 1-172, 1993
  6. Weidlich U., Gmehling J., "A Modified UNIFAC Model. 1. Prediction of VLE, hE, and γ", Ind.Eng.Chem.Res., 26(7), p1372-1381, 1987
  7. Wittig R., Lohmann J., Gmehling J., „Vapor-Liquid Equilibria by UNIFAC Group Contribution. 6. Revision and Extension“, Ind.Eng.Chem.Res., 42(1), 183-188, 2003
  8. </Polka H.-M., „Experimentelle Bestimmung und Berechnung von Dampf-Flüssig-Gleichgewichten für Systeme mit starken Elektrolyten", Doktorarbeit, Univ. Oldenburg, 1-144, 1993
  9. 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
  10. Holderbaum T., „Die Vorausberechnung von Dampf-Flüssig-Gleichgewichten mit einer Gruppenbeitragszustandsgleichung“, Fortschrittsber.VDI Reihe 3, 243, 1-154, 1991
  11. 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
  12. 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
Bearbeiten

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