Wolfgang Schmidt

Berechnung von Stoffdaten und Phasengleichgewichten mit Excel-VBA


Скачать книгу

wie etwa NRTL, Uniquac, Unifac etc. betrachten Wechselwirkungen nur zwischen ungleichen Molekülen. Mit der Lennard-Jones-Methode können Wechselwirkungen zwischen jeder Art von Molekülen berechnet werden.

      Die Elliot-Unifac-Methode enthält neben den bekannten Joback-Gruppen sogenannte Doppelgruppen. Diese wurden der Unifac-Methode zur Berechnung von Wechselwirkungsparametern entliehen. Allerdings werden alle 4 Isomeren, also 2,2-, 2,3-, 2,4- und 2,5-Dimethylhexan mit den gleichen Unifac-Inkrementen beschrieben, nämlich Unifac Sub-Group 1 (CH3) *4, Unifac Sub-Group 2 (CH2) * 2, Unifac Sub-Group 3 (CH) *2. Somit ergeben sich auch für alle dieselben Joback-Unifac-Daten und wir erhalten für den Siedepunkt Tb = 383,658 K. Dieser Wert ist aber besser als der Joback-Wert mit 375,69 K wie der Vergleich mit der Tabelle Second-order-molecules in der Exceldatei Ex_03_07_Group_Contribution zeigt. Eine genauere Strukturanalyse lässt auch Unifac nicht zu.

      Das 2,3- Dimethylhexan hat die Struktur

      Die beiden seitlichen CH3-Gruppen befinden sich am 2. und 3. C-Atom, von links gezählt. Es fällt auf, dass sowohl beide CH3-Gruppen als auch die CH-Gruppen nebeneinander liegen. Beim 2,4-Dimethylhexan

images

      befinden sich die beiden seitlichen CH3-Gruppen am 2. und 4. C-Atom, von links gezählt. Deren Abstand ist damit größer als der Abstand bei 2,3-Dimethylhexan. Auch befindet sich zwischen den beiden CH-Gruppen jetzt eine CH2-Gruppe. Es ist das Ziel weiterer Bemühungen, diese Strukturunterschiede so zu definieren, dass sich daraus bessere Ergebnisse erhalten lassen.

      In der o.g. CHEMCAD-Dokumentation befindet sich eine umfangreiche Beschreibung der Cavett, API, Lee-Kesler, Joback- und der Elliot-Unifac-Methode, sowie Daten für die beiden letztgenannte Methoden. Im neuen VDI-Heat-Atlas ist u.a. die Joback- sowie die Second-Order-Constantinou-/Gani-Methode beschrieben.

      Die DIPPR-Gesellschaft „Design Institute for Physical Property Research“ hat es sich zum Ziel gesetzt, Stoffdatenfunktionen und deren Parameter von Reinstoffen als Datenbank zu speichern und auf dem Markt anzubieten.

      Die Verwendung der DIPPR-Datenbank soll am Beispiel Ethanol in der Prozesssimulationssoftware CHEMCAD 6.5 dargestellt werden. Alle Screenshots sind der Software CHEMCAD entnommen.

      Die Gleichung für die spezifische Wärmekapazität Cp lautet:

      (1.28)image

      Bei näherer Betrachtung der Parameter A bis F fällt auf, dass diese in der genannten Reihenfolge in etwa um den Faktor 1/1000 kleiner werden. Wäre dem nicht so, bestünde die Gefahr einer wellenförmigen Funktion. Dies ist bekanntlich der typische Nachteil der ansonsten einfachen Polynomfunktion. So einfach ein Polynom auch ist, es gibt keinen physikalischen Zusammenhang, der einer Polynomfunktion entspricht. Daher ist ein Polynom zur Beschreibung einer physikalischen Größe nicht geeignet. Die Aussage, mit einem Polynom könne man jeden beliebigen Funktionsverlauf darstellen, stimmt eben nur sehr begrenzt.

      Mit diesen minimalen Daten lassen sich notfalls fehlende Stoffdaten durch Korrelationen abschätzen. Natürlich kann man damit nicht die Genauigkeit gemessener Daten ersetzen.

      Weiterhin findet man die Standardbildungsenthalpie [kJ/kg] (heat of formation) und die Gibbs-Enthalpie of formation [kJ/kg]. Der Zusammenhang beider Enthalpien geht aus der Gleichung hervor.

      (1.29)image

      In diesem Fall sind nur die Daten für die im idealen Zustand bekannt.

      (1.30)image

      Diese Funktion benötigt nur die Parameter A bis D. Diese 4 Parameter liegen für Ethanol vor. Der untere Gültigkeitswert bei Tu = 159,05 K beträgt 19,4128 kmol/m3, bei T0 = 513,92 K beträgt er hingegen 5,96518 kmol/m3.

      Für die Feststoffdichte gilt die DIPPR-Formel Nr. 100. Diese lautet

      (1.31)image

      In dieser Funktion genügt im Minimalfall sogar nur 1 Parameter. Und tatsächlich liegt auch nur ein einziger Parameter, nämlich 22,9 kmol/m3, für Ethanol vor. D.h. die Dichte von festem Ethanol ist also nur mit 22,9 kmol/m3 in der Datenbank enthalten. Vermutlich liegen keine weiteren Messungen vor.

image

      Der Dampfdruck wird nach der DIPPR-Funktion 102 berechnet. Diese lautet

      (1.32)image

      Die Funktion entspricht der Clausius-Clapeyron-Formel und einer empirischen Erweiterung mit den Parametern C, D und E.

      Die Verdampfungsenthalpie wird nach der DIPPR-Funktion