UNIFAC

Dampf-Flüssig-Gleichgewicht mit Sattelazeotrop
berechnet mit UNIFAC bei P=1 atm. Die roten Linien sind Dampfzusammensetzungen, die blauen die Zusammensetzung der Flüssigkeit.

UNIFAC^[1][2] (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 Technischen Chemie häufig eingesetzt wird.

UNIFAC ist von UNIQUAC abgeleitet. UNIQUAC verwendet stoffspezifische Parameter, die aus experimentell ermittelten Daten bestimmt werden müssen. Diese Parameter werden in UNIFAC vorhergesagt, so dass keine gemessenen Daten mehr benötigt werden.

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

Beispiel einer Gruppenzuordnung in UNIFAC

Das Prinzip der Gruppenbeitragsmethoden besteht im Wesentlichen darin, 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^[3] 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.

Aktuelle Parametrisierung

UNIFAC-Parametermatrix (letzte veröffentlichte Version^[4])

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.

Gleichungen

Ausführliche Formulierungen können im Originalartikel gefunden werden oder auch in Lehrbüchern^[5][6].

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 und wird aus van der Waalsschen Oberflächen (F) und Volumina (V) sowie Molenbrüchen (x) berechnet.

$ln \; \gamma_i^C = 1 - V_i + ln \; V_i - 5 q_i \left( 1 - { {V_i} \over {F_i} } + ln \; { {V_i} \over {F_i} } \right)$

mit

$V_i = { {r_i} \over { \sum_j {r_j x_j} } }$

$F_i = { {q_i} \over { \sum_j {q_j x_j} } }$

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

$q_i = \sum { \nu_k^{(i)} Q_k}$

$r_i = \sum { \nu_k^{(i)} R_k}$

ν_k⁽ⁱ⁾ ist die Häufigkeit der Gruppe k im Molekül i.

Restanteil

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

Der Restanteil wird aus Gruppenaktivitätskoeffizienten Γ_k berechnet.

$ln \; \gamma_i^R = \sum { \nu_k^{(i)} \left( ln \; \Gamma_k - ln \; \Gamma_k^{(i)} \right)}$

Die Gruppenaktivitätskoeffizienten in der Mischung Γ_k und im Reinstoff Γ_k⁽ⁱ⁾ werden über die Beziehung

$ln \; \Gamma_k = Q_k \left[ 1 - ln \left( \sum_m{\Theta_m \Psi_{mk}} \right) - \sum_m{{\Theta_m \Psi_{km}} \over {\sum_m {\Theta_n \Psi_{mn} }}} \right]$

berechnet.

Θ ist der Oberflächenanteil

$\Theta={{Q_m X_m} \over {\sum_n{Q_n X_n}}}$

und X der Gruppenmolenbruch.

$X_m={ {\sum_j {\nu_m^{(j)}x_j}} \over {\sum_j \sum_m {\nu_n^{(j)} x_j}} }$

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

$\Psi = exp { \frac{-a_{nm}}{T}}$ mit a_nm als Wechselwirkungsparameter zwischen den Gruppen.

Beispielrechnung

Die Bilder zeigen das Dampf-Flüssiggleichgewicht des Gemischs aus Toluol und Methanol. Die roten Linien zeigen die mit (mod.) UNIFAC berechneten Daten. Für dieses System zeigt die Berechnung eine sehr gute Übereinstimmung mit den experimentell ermittelten Daten.

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^[7]. 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^[8] 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 Prozesssimulation 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^[9]
• die Abschätzung von Viskositäten^[10]
• die Integration von UNIFAC in Mischungsregeln für Zustandsgleichungen (PSRK)^[11]
• die Vorhersage von UNIFAC-Wechselwirkungsparametern^[12]
• die Erweiterung des UNIFAC-Modells auf spezielle Stoffgruppen^[13]
• die Ableitung des UNIFAC-Modells zur Vorhersage von Verdampfungsenthalpien und Exzessenthalpien: Gruppenbeitragsmodelle UNIVAP & EBGCM^[14]
• die Überarbeitung von Wechselwirkungsparametern für bestehende Gruppen und die Ergänzung von Wechselwirkungsparametern für neue funktionelle Gruppen, um die Anwendbarkeit und die Qualität des UNIFAC-Modells zu verbessern (UNIFAC-Konsortium).

Literatur

1. ↑ Fredenslund A., Jones R.L., Prausnitz J.M., „Group-Contribution Estimation of Activity Coefficients in Nonideal Liquid Mixtures“, AIChE J., 21(6), 1086-1099, 1975
2. ↑ Å. Fredenslund, J. Gmehling, P. Rasmussen; „Vapor-Liquid Equilibria Using UNIFAC - A Group Contribution Method;“ Elsevier, Amsterdam (1977)
3. ↑ Bondi A., "van der Waals Volumes and Radii", J.Phys.Chem., 68(3), 441-451, 1964
4. ↑ 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
5. ↑ J. Gmehling, B. Kolbe; „Thermodynamik“, 2. Auflage; VCH-Verlag, Weinheim (1992)
6. ↑ J. Vidal, „Thermodynamique. Méthodes appliquées au raffinage et au Génie Chimique“, Editions Technip, 1997
7. ↑ Gudrun Wienke, „Messung und Vorausberechnung von n-Octanol/Wasser-Verteilungskoeffizienten", Doktorarbeit, Univ. Oldenburg, 1-172, 1993
8. ↑ Weidlich U., Gmehling J., „A Modified UNIFAC Model. 1. Prediction of VLE, hE, and γ^∞“", Ind.Eng.Chem.Res., 26(7), 1372-1381, 1987
9. ↑ Polka H.-M., „Experimentelle Bestimmung und Berechnung von Dampf-Flüssig-Gleichgewichten für Systeme mit starken Elektrolyten", Doktorarbeit, Univ. Oldenburg, 1-144, 1993
10. ↑ 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
11. ↑ Holderbaum T., „Die Vorausberechnung von Dampf-Flüssig-Gleichgewichten mit einer Gruppenbeitragszustandsgleichung“, Fortschrittsber.VDI Reihe 3, 243, 1-154, 1991
12. ↑ 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
13. ↑ 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
14. ↑ Ulbig, P., „Gruppenbeitragsmodelle UNIVAP & EBGCM Entwicklung der Gruppenbeitragsmodelle UNIVAP und EBGCM zur Vorhersage thermodynamischer Größen sowie Bestimmung der Modellparameter unter Verwendung evolutionärer Algorithmen“, Hannemann-Verlag, 1996

Weblinks

• UNIFAC-Industriekonsortium an der Carl von Ossietzky Universität Oldenburg
• Berechnung von Aktivitätskoeffizienten in einigen binären Gemischen
• UNIFAC-Strukturgruppen und -Parameter