Wie Wird Ki Für Überkritische Feedgase Bestimmt? Meistern Sie Verdampfungsberechnungen Für Reaktoren

Für ein überkritisches Feedgas lässt sich der Verdampfungsfaktor (Kᵢ) nicht aus Dampfdrucktabellen ermitteln – er wird mithilfe eines unsymmetrischen thermodynamischen Modells bestimmt. Die maßgebliche Beziehung lautet Kᵢ = (Hᵢ · γᵢ) / (φᵢ · P)*, wobei Hᵢ die Henry-Konstante des Gases im Lösungsmittel ist, γᵢ* der unsymmetrischen Konvention folgt (γᵢ* → 1 für xᵢ → 0), φᵢ der Fugazitätskoeffizient der Dampfphase und P der Gesamtdruck ist. Dieser Ausdruck ersetzt den herkömmlichen K-Wert nach dem Raoultschen Gesetz, da der Rein-Gas-Dampfdruck oberhalb der kritischen Temperatur nicht definiert ist.

Die Bestimmung von Kᵢ für überkritische Feedströme wie Wasserstoff oder Sauerstoff basiert vollständig auf dem Henry-Gesetz mit einem unsymmetrischen Aktivitätskoeffizienten. Aus dem Wert Kᵢ = yᵢ/xᵢ wird (Hᵢ · γᵢ*) / (φᵢ · P), der Löslichkeit, molekulare Wechselwirkungen und Nichtidealität der Gasphase verbindet – eine direkte Folge davon, dass das gelöste Gas als gelöster Stoff behandelt wird, der nie zu einer reinen Flüssigkeit kondensiert.

Die thermodynamische Grundlage für überkritische Gase

Warum das traditionelle Dampf-Flüssig-Gleichgewicht versagt

Der klassische K-Faktor, yᵢ/xᵢ = Pᵢˢᵃᵗ/P, basiert auf dem Rein-Komponenten-Dampfdruck.
Für eine Komponente, deren Temperatur ihre kritische Temperatur übersteigt (z. B. Wasserstoff bei 298 K), existiert Pᵢˢᵃᵗ nicht. Das Gas kann nicht allein durch Druck verflüssigt werden, daher verliert die übliche symmetrische Konvention ihre Bedeutung.

Ingenieure, die Pilotanlagen für Gas-Flüssig-Reaktoren entwickeln, müssen daher zu einem Modell wechseln, bei dem die überkritische Komponente als gelöster Stoff behandelt wird, nicht als kondensierbarer Dampf. Die Gleichgewichtsbedingung fᵢⱽ = fᵢᴸ gilt nach wie vor, aber die Fugazität der flüssigen Phase muss auf das Henry-Gesetz bezogen werden statt auf die Rein-Flüssig-Fugazität.

Die unsymmetrische Konvention und das Henry-Gesetz

Die Flüssigphasen-Fugazität für einen überkritischen gelösten Stoff wird geschrieben als fᵢᴸ = γᵢ · Hᵢ · xᵢ.*
Hierbei ist Hᵢ die Henry-Konstante (eine Eigenschaft des Gelöster-Lösungsmittel-Paares bei der Systemtemperatur), und γᵢ* ist ein Aktivitätskoeffizient, der so normiert ist, dass γᵢ* → 1, wenn der Molenbruch xᵢ gegen Null geht. Diese unsymmetrische Konvention bindet den Standardzustand an unendliche Verdünnung, genau dort, wo das Henry-Gesetz gilt.

Hᵢ erfasst das Grenzverhalten der Löslichkeit – je stärker die Wechselwirkung zwischen Gas und Flüssigkeit ist, desto kleiner ist Hᵢ und desto größer ist die gelöste Konzentration. γᵢ* korrigiert dann Abweichungen vom idealverdünnten Verhalten, wenn xᵢ nicht vernachlässigbar klein ist.

Der vollständige Ausdruck für Kᵢ

Das Gleichsetzen der Fugazitäten von Dampf- und Flüssigphase ergibt φᵢ · yᵢ · P = γᵢ · Hᵢ · xᵢ.*
Durch Umstellen erhält man direkt den Verdampfungsfaktor:

[ K_i = \frac{y_i}{x_i} = \frac{H_i \cdot \gamma_i^*}{\phi_i \cdot P} ]

Jeder Term ist zugänglich: Hᵢ stammt aus Literaturkorrelationen oder einer begrenzten Anzahl von Löslichkeitsmessungen in der Pilotanlage; γᵢ* aus Modellen für überschüssige Gibbs-Energie (z. B. NRTL, UNIQUAC), die an die unsymmetrische Konvention angepasst sind; φᵢ aus einer kubischen Zustandsgleichung (Peng-Robinson, Soave-Redlich-Kwong); und P ist einfach der Betriebsdruck des Reaktors. Dieser Ausdruck ist der logische und notwendige Ersatz für das Analogon nach dem Raoultschen Gesetz, wenn ein Feedstrom oberhalb seiner kritischen Temperatur zugeführt wird.

Anwendung der Gleichung im Pilotanlagenbetrieb

Ermittlung der Henry-Konstante

Hᵢ ist der Grundpfeiger der Berechnung.
Sie ist stark temperaturabhängig und muss für mittlere bis hohe Drücke über die Krichevsky-Kasarnovsky-Gleichung [ \ln H_i(P) = \ln H_i(P_{\text{sat}}) + \frac{\bar{v}i^\infty}{RT}(P - P{\text{sat}}) ] korrigiert werden, wobei (\bar{v}_i^\infty) das partielle molare Volumen des gelösten Stoffes bei unendlicher Verdünnung ist.

Für viele in Pilotanlagen verwendete Gase (H₂, O₂, CO) in gängigen organischen Lösungsmitteln sind zuverlässige Hᵢ-Werte tabelliert oder können über Gruppenbeitragsmethoden vorhergesagt werden. Fehlen Daten, lässt sich die Henry-Konstante unter den Temperaturbedingungen des Reaktors schnell durch wenige Löslichkeitsexperimente im Labor bestimmen.

Abschätzung des Aktivitätskoeffizienten (γᵢ*)

Selbst in relativ verdünnten Lösungen kann γᵢ von Eins abweichen, wenn das Lösungsmittel komplex ist oder mehrere gelöste Stoffe wechselwirken.*
Bleibt der Molenbruch des überkritischen Gases in der Flüssigphase unter etwa 0,01, führt die Annahme γᵢ* = 1 oft zu einem Fehler von weniger als 5 % – was für die Dimensionierung von Pilotanlagen in der frühen Phase akzeptabel ist.

Für stärker konzentrierte Systeme können Modelle wie UNIQUAC oder NRTL für die unsymmetrische Konvention neu parametrisiert werden. Die notwendigen binären Wechselwirkungsparameter werden typischerweise aus Gaslöslichkeitsdaten regressiert oder über molekulare Ähnlichkeit abgeschätzt. Der Entscheidungsbaum für den Pilotanlageningenieur teilt sich also auf: Verwenden Sie γᵢ* = 1 für eine schnelle Vorauslegung, oder investieren Sie in ein vollständiges Aktivitätsmodell, wenn eine genaue Vorhersage des Stofftransfers entscheidend ist.

Berechnung des Dampfphasen-Fugazitätskoeffizienten (φᵢ)

Bei den mäßigen Drücken, die in vielen Pilotreaktoren üblich sind (5–20 bar), weicht φᵢ selten um mehr als wenige Prozent von Eins ab.
Hochdruckhydrierungen oder Oxidationen arbeiten jedoch oft bei 50–200 bar, wo die Nichtidealität der Gasphase signifikant wird. Eine kubische Zustandsgleichung, passend zur Gasmischung und Temperatur, liefert φᵢ auf einfache Weise.

Im Peng-Robinson-Modell ist φᵢ beispielsweise eine Funktion von reduzierter Temperatur, Druck und azentrischem Faktor. Für reinen Wasserstoff kann die φᵢ-Korrektur Kᵢ bei 100 bar um 10–30 % gegenüber der Idealgasannahme erhöhen, sodass die Ignorierung dieses Terms zu einer Unterschätzung des Gasphasenmolenbruchs führt und die Berechnung der Stofftransportantriebskraft verzerrt.

Die Rolle des Drucks

Der Druck erscheint explizit im Nenner des Kᵢ-Ausdrucks, sodass Kᵢ umgekehrt proportional zu P ist, wenn alle anderen Terme konstant bleiben.
Das bedeutet: Wenn der Reaktordruck steigt, sinkt der Gleichgewichts-Gasphasenmolenbruch bei gegebener flüssigphasiger Konzentration – eine direkte Reflektion des Henry-Gesetzes: Mehr Gas löst sich, also weniger verbleibt im Kopfraum oder in der Blasenphase.

Bei der Auslegung von Pilotanlagen muss dieser Druckeffekt mit der oft übersehenen Druckabhängigkeit von Hᵢ selbst kombiniert werden. Bei hohen Drücken steigt Hᵢ, was den P-Effekt teilweise kompensiert. Ein vollständiges Kᵢ-Modell koppelt daher die Krichevsky-Kasarnovsky-Korrektur mit dem Nenner und stellt sicher, dass der Verdampfungsfaktor die Löslichkeitserhöhung unter tatsächlichen Betriebsbedingungen korrekt widerspiegelt.

Verständnis von Kompromissen und Grenzen

Die Verdünnungsannahme

Die Behandlung von γᵢ* = 1 ist verlockend, aber nicht immer sicher. In hochkonzentrierten Hydrierungsgemischen oder wenn der Anteil des gelösten Gases einige Molprozent übersteigt, kann der Aktivitätskoeffizient um 10–20 % von Eins abweichen. Der resultierende Kᵢ-Fehler überträgt sich auf Schätzungen der flüssigphasigen Konzentrationen und letztendlich auf Vorhersagen des Reaktorvolumens. Für eine frühe Machbarkeitsstudie ist der Fehler oft tolerabel; für die endgültige Scale-up ist er es nicht.

Temperatur- und Druckempfindlichkeit

Henry-Konstanten für leichte Gase verdoppeln oder verdreifachen sich typischerweise über eine Temperaturänderung von 50 K, und der Fugazitätskoeffizient ändert sich mit dem Druck. Ein Kᵢ-Wert, der an einem Betriebspunkt gültig ist, kann ungenau werden, wenn sich das Temperaturprofil des Reaktors ändert – was häufig bei Anfahrvorgängen oder Katalysatordeaktivierung der Fall ist. Das Pilotanlagenteam muss Kᵢ als dynamische Variable behandeln, nicht als statische Zahl, wenn es den Reaktor über seinen gesamten Betriebsbereich modelliert.

Mehrkomponentengemische

Wenn der Gasfeed mehrere überkritische Komponenten enthält (z. B. H₂ und CO im Synthesegas), wird der Fugazitätskoeffizient φᵢ zu einer Eigenschaft der Mischung, und Kreuzwechselwirkungsparameter in der Zustandsgleichung spielen eine Rolle. Die Verwendung von Rein-Komponenten-φᵢ-Schätzungen kann zu systematischen Verzerrungen führen. Wenn möglich, sollten Pilotanlageningenieure Kᵢ-Vorhersagen anhand einiger Messungen von Gesamtdruck und Dampfzusammensetzung validieren und die EOS-Mischungsregeln bei Bedarf verfeinern.

Die richtige Wahl für Ihr Pilotanlagenziel

Der passende Detaillierungsgrad bei der Bestimmung von Kᵢ hängt vollständig davon ab, wofür Sie den Wert benötigen.

Wenn Ihr Hauptaugenmerk auf einer schnellen Machbarkeitsprüfung liegt: Verwenden Sie die vereinfachte Form Kᵢ ≈ Hᵢ/P mit Literaturwerten für Hᵢ und gehen Sie von γᵢ* = φᵢ = 1 aus. Damit erhalten Sie Größenordnungsabschätzungen für die Gas-Flüssig-Teilung und Mindestlösungsmittelflüsse an einem einzigen Nachmittag.
Wenn Ihr Hauptaugenmerk auf einer genauen Stofftransfermodellierung liegt: Verwenden Sie den vollständigen Ausdruck mit einer druckkorrigierten Hᵢ (über die Krichevsky-Kasarnovsky-Beziehung), einer modernen EOS für φᵢ und mindestens einem erstgradigen Aktivitätsmodell für γᵢ*. Ihre Simulationen von Blasensäulen oder Rührkesseln werden Sie mit KLa-Werten belohnen, die auch bei der Maßstabsvergrößerung Bestand haben.
Wenn Ihr Hauptaugenmerk auf dem Hochdruckbetrieb (oberhalb 30 bar) liegt: Vernachlässigen Sie niemals die Fugazitätskoeffizientenkorrektur. Eine kubische EOS – angepasst an veröffentlichte Binärdaten oder gemessene Taupunkte – ist die Mindestanforderung, um eine systematische Unterschätzung von Gasphasenkonzentrationen zu vermeiden.
Wenn Ihr Hauptaugenmerk auf dem Lösungsmittelscreening oder mehrkomponentigen Feedströmen liegt: Investieren Sie in experimentelle Löslichkeitsläufe, um Hᵢ- und γᵢ*-Parameter zu regressieren. Der Aufwand amortisiert sich durch eine drastische Reduzierung der Unsicherheit bei der Reaktordimensionierung und vermeidet Überraschungen in nachgeschalteten Trenneinheiten.

Die genaue Bestimmung von Kᵢ für überkritische Feedströme ist keine Übung aus dem Lehrbuch – sie ist der Dreh- und Angelpunkt, der Löslichkeit, Stofftransport und Reaktordimensionierung zu einer vertrauenswürdigen Pilotanlagenauslegung verbindet.

Zusammenfassungstabelle:

Parameter	Symbol	Rolle in der $K_i$-Gleichung	Wesentliche Abschätzmethode
Verdampfungsfaktor	$K_i$	Misst das Gas-Flüssig-Gleichgewicht ($y_i/x_i$)	Berechnet als $(H_i \cdot \gamma_i^*) / (\phi_i \cdot P)$
Henry-Konstante	$H_i$	Definiert das Grenzverhalten der Löslichkeit	Korrektur nach der Krichevsky-Kasarnovsky-Gleichung
Aktivitätskoeffizient	$\gamma_i^*$	Korrigiert Nichtidealität in der Flüssigphase	NRTL / UNIQUAC (unsymmetrische Konvention)
Fugazitätskoeffizient	$\phi_i$	Korrigiert Nichtidealität in der Dampfphase	Kubische Zustandsgleichung (z. B. Peng-Robinson)
Systemdruck	$P$	Stellt den Reaktorbetriebsdruck dar	Expliziter Nenner; umgekehrt proportional zu $K_i$

Die Optimierung thermodynamischer Modelle für Gas-Flüssig-Reaktoren erfordert präzise Versuche im Pilotmaßstab. LABPARK bietet hochwertige Bildungs- und Berufsbildungs-Pilotanlagen für Einheitsoperationen in Chemieingenieurwesen, Bioprozess- und Biotechnologie sowie Umwelt- und Wasseraufbereitung für Universitäten, Forschungsinstitute und Unternehmen. Stärken Sie Ihre Forschung und Ausbildung mit zuverlässigen, für die Maßstabsvergrößerung geeigneten Systemen – kontaktieren Sie LABPARK noch heute, um Ihr Projekt zu besprechen!