Wie berechnet man den Realgaszustand mit der Peng-Robinson-Gleichung

Berechnen Sie den Realgaszustand mit Peng-Robinson 1976: P = RT/(V-b) - a·α(T)/[V(V+b)+b(V-b)]. VLE-Gleichgewichte, Kohlenwasserstoff-CO₂, NIST REFPROP, AspenTech HYSYS, petrochemische Simulation.

Berechnungsparameter
War dieser Rechner nützlich?

Wie berechnet man den Realgaszustand mit der Peng-Robinson-Gleichung ermöglicht die Schätzung thermodynamischer P-V-T-Eigenschaften von Kohlenwasserstoffen, Erdgas, überkritischem CO₂ und Industriefluiden. Entwickelt von Ding-Yu Peng und Donald B. Robinson an der University of Alberta in 1976 (Veröffentlichung I&EC Fundamentals 15:59-64), verbessert die PR-Gleichung die klassische Soave-Redlich-Kwong (SRK 1972) für Flüssigkeiten und Flüssigkeitsdichten, besonders C5+-Kohlenwasserstoffe und Erdölsysteme. Analytische Form: P = RT/(V-b) - a·α(T)/[V(V+b)+b(V-b)], Parameter: a = 0,45724·R²·T_c²/P_c, b = 0,07780·R·T_c/P_c, α(T) = [1 + κ(1-√(T/T_c))]² mit κ = 0,37464 + 1,54226·ω - 0,26992·ω² (ω = Pitzer 1955 azentrischer Faktor). Kubische Form in Z: Z³ - (1-B)Z² + (A-3B²-2B)Z - (AB-B²-B³) = 0. De-facto-Standard in Aspen HYSYS, Aspen Plus, ProMax, ProSim, Petro-SIM, OLGA, PIPESIM, GERG-2008, Multiflash; für Erdgasverarbeitung, Raffinerien, Flüssiggas, CO₂-Sequestrierung, Gas-to-liquid GTL und Kryotechnik. Grenze: 5-10% Fehler für polare Moleküle (Wasser, Alkohole) — dafür SAFT, CPA oder NRTL+UNIFAC verwenden.

1. Peng-Robinson 1976 Gleichung

Die Peng-Robinson Zustandsgleichung in kubischer Form P(V,T): P = RT/(V-b) - a·α(T)/[V(V+b)+b(V-b)]. Parameter a (Anziehungskräfte) und b (ausgeschlossenes Volumen, Kovolumen) berechnet bei kritischen Bedingungen: a = Ω_a·R²T_c²/P_c mit Ω_a = 0,45724, b = Ω_b·R·T_c/P_c mit Ω_b = 0,07780. Die α(T)-Funktion führt Abhängigkeit von reduzierter Temperatur T_r = T/T_c ein. Azentrischer Faktor ω (Pitzer 1955): ω(Methan)=0,011, ω(Propan)=0,152, ω(n-Decan)=0,489, ω(Wasser)=0,344. Z = PV/RT ist Lösung des kubischen Polynoms — analytisch Cardano oder numerisch Newton-Raphson.

2. Peng-Robinson-Parameter und reine Fluide

  • Methan CH₄: T_c=190,56 K, P_c=4,599 MPa, ω=0,011 — Erdgas.
  • Ethan C₂H₆: T_c=305,32 K, P_c=4,872 MPa, ω=0,099 — Flüssiggas.
  • Propan C₃H₈: T_c=369,83 K, P_c=4,248 MPa, ω=0,152 — Kältemittel R-290.
  • n-Butan C₄H₁₀: T_c=425,12 K, P_c=3,796 MPa, ω=0,200 — Flüssiggas.
  • Kohlendioxid CO₂: T_c=304,13 K, P_c=7,377 MPa, ω=0,225 — Sequestrierung, R-744.
  • Stickstoff N₂: T_c=126,19 K, P_c=3,396 MPa, ω=0,037 — kryogen inert.
  • Wasser H₂O: T_c=647,1 K, P_c=22,064 MPa, ω=0,344 — durch Polarität begrenzt.
  • Wasserstoff H₂: T_c=33,1 K, P_c=1,30 MPa, ω=-0,216 — Quanten, Korrektur nötig.

3. Gültigkeitsbereich und Korrekturen

PR funktioniert optimal für nichtpolare Kohlenwasserstoffe (paraffinisch, olefinisch, aromatisch), CO₂, N₂, H₂S, Erdgas. Fehler <3% für Dampfdichte und binäre VLE. Für reine Flüssigkeiten typischer Dichtefehler 6-12% (Überschätzung); Peneloux 1982 Volumenshift verbessert auf 1-3%. Für polare Moleküle (Wasser, Alkohole, Amine, Glykole) Fehler 10-30%; CPA (Kontogeorgis 1996), SAFT, PC-SAFT (Gross 2001) verwenden. Für Wasserstoff und Helium κ mit Twu 1991 anpassen oder GERG-2008 verwenden. Für überkritische Fluide (scCO₂) ist PR mit Peneloux-Shift hervorragend.

4. Numerische Beispiele

Beispiel 1 — Überkritisches CO₂ für scCO₂-Extraktion. T = 350 K, P = 20 MPa. a = 396,2 Pa·m⁶/mol². b = 26,67e-6 m³/mol. α(350K) = 0,888. PR kubisch → Z = 0,57, V = 8,34e-5 m³/mol, ρ = 528 kg/m³. Typische scCO₂-Dichte für Koffein-, Hopfenölextraktion.

Beispiel 2 — Erdgas (95% CH₄, 5% C₂H₆) bei 80 bar, 280 K. Mischungsregeln Van der Waals: k_ij(CH₄-C₂H₆) ≈ -0,003. Berechnetes Z 0,855, Kompressibilitätsfaktor entscheidend für Pipeline-Auslegung ASME B31.8. AGA-8 und GERG-2008 Alternativen.

Beispiel 3 — Reines Propan VLE bei 300 K. PR kubisch liefert 3 Z-Wurzeln. Flüssig (kleineres Z, ρ_L=523 kg/m³, 3% Fehler vs NIST) und Dampf (größeres Z, ρ_V=19 kg/m³). P_sat=1,03 MPa (vs NIST 0,996 MPa, 3,4% Fehler).

5. Praktische Anwendungen der Peng-Robinson

  • Petrochemische Prozesssimulation: Aspen HYSYS, Aspen Plus, ProMax, ProSim, Petro-SIM, OLGA. Standard ExxonMobil, Shell, BP, TotalEnergies, BASF, Linde
  • Erdgas-Pipeline-Auslegung: ASME B31.8, ISO 13443, AGA-8 — reale Z-Berechnung für Pipeline-Dimensionierung (OGE, Gascade, Open Grid Europe)
  • Flüssiggas und LNG: kryogene Zyklen Linde, Air Products APCI, ConocoPhillips Optimized Cascade — Flüssigkeitsdichte und Dampfdruck für Pumpen und BOG
  • CO₂-Sequestrierung (CCS): Injektion in erschöpften Lagerstätten (Sleipner Norwegen, Quest Kanada), überkritischer Pipelinetransport
  • scCO₂-Extraktion: Koffein (Maxwell House), Hopfenöle (Guinness Heineken Bier), Gewürzöle, CBD aus Cannabis
  • Lagerstätten-Engineering: PVT-Simulation (CMG WinProp, Schlumberger Eclipse)
  • Gas-to-liquid GTL Fischer-Tropsch: Sasol Mossel Bay, Shell Pearl Qatar
  • Hydrat-Hemmung: Glykole (MEG, DEG, TEG) zur Vermeidung von Offshore-Hydratbildung

6. Vergleich mit anderen Zustandsgleichungen

Peng-Robinson (PR) 1976 verbessert Soave-Redlich-Kwong (SRK) 1972 für Flüssigkeiten. Van der Waals 1873 nur historisch (30%+ Fehler). Benedict-Webb-Rubin (BWR): 8-11 Parameter, hohe Rechenkosten. Lee-Kesler 1975: hervorragend für Erdgas mit Pitzer-Tabulierungen. SAFT (Chapman 1990): Störungstheorie, hervorragend für Polymere. PC-SAFT (Gross 2001): perturbed-chain. CPA (Kontogeorgis 1996): SRK+Assoziation für Wasser, Glykole. GERG-2008 (Kunz-Wagner): fundamentale 21-Komponenten-Gleichung ISO 20765-2.

7. Peng-Robinson-Entwicklungen und Korrekturen

Moderne Varianten: PRSV 1986 (Stryjek-Vera): modifiziert α(T) für polare. Mathias-Copeman 1983: 3 Parameter κ_1,2,3 für polare. Twu 1991: verallgemeinerte α für Kryotechnik. Volumentranslation Peneloux 1982 und Magoulas-Tassios 1990 korrigiert Flüssigkeitsdichte. Wong-Sandler 1992: GE-basierte Mischungsregeln aus Helmholtz für polare Systeme. Open-Source-Software: CoolProp (Bell-Wronski-Quoilin 2014), Cantera, thermoextensions R/Julia. Kommerziell: NIST REFPROP, Multiflash Infochem-KBC, VMGSim Schlumberger, Symmetry SLB. Für KI/ML entstehen seit 2020 NeuralNetwork+EOS-Hybride (NIST mit MIT) zur Erweiterung der Genauigkeit unter Erhalt thermodynamischer Konsistenz.

Aktualisiert im Juli 2026

Häufig gestellte Fragen

Was ist die Peng-Robinson-Gleichung?

Die Peng-Robinson-Gleichung (1976, Univ. Alberta) ist eine kubische Zustandsgleichung für Realgase: P = RT/(V-b) - a·α(T)/[V(V+b)+b(V-b)]. Sie ist der Standard für Kohlenwasserstoffe, Erdgas, CO₂ in petrochemischer Prozesssimulation (Aspen HYSYS). Verbessert SRK 1972 für Flüssigkeitsdichten. Typischer Fehler <3% für nichtpolare VLE, 5-10% für polare (Wasser, Alkohole).

Wie werden die Parameter a und b berechnet?

a = 0,45724·R²·T_c²/P_c (Anziehungskräfte); b = 0,07780·R·T_c/P_c (Kovolumen); α(T) = [1 + κ(1-√(T/T_c))]² mit κ = 0,37464 + 1,54226·ω - 0,26992·ω². ω ist der Pitzer 1955 azentrische Faktor: ω(CH₄)=0,011, ω(C₃H₈)=0,152, ω(CO₂)=0,225, ω(H₂O)=0,344. Kritische Konstanten T_c, P_c, ω tabelliert in NIST REFPROP und Reid-Prausnitz-Poling.

Was ist der Pitzer-azentrische Faktor ω?

Der azentrische Faktor ω (Kenneth Pitzer 1955) misst, wie stark ein Molekül von der sphärischen Symmetrie abweicht. Definiert ω = -log₁₀(P_sat/P_c)|_T_r=0.7 - 1. Für sphärische Moleküle (Ar, Kr, Xe): ω≈0. Für Kettenkohlenwasserstoffe: wächst mit n (n-Decan ω=0,489). Für polare: H₂O 0,344, NH₃ 0,253, Methanol 0,556. ω ist der 3. Corresponding-States-Parameter nach T_c, P_c.

Wie löst man die Peng-Robinson-Kubische?

Man schreibt sie in Z: Z³ - (1-B)Z² + (A-3B²-2B)Z - (AB-B²-B³) = 0, mit A=aP/R²T², B=bP/RT. Lösung mit Cardano (analytisch) oder Newton-Raphson. 3 reelle Wurzeln unter T_c: die größte ist Dampf, die kleinste Flüssig, die mittlere unphysikalisch. Über T_c: nur 1 reelle Wurzel (Gas oder überkritisch). Software: CoolProp, REFPROP, Aspen HYSYS Property Method PR.

Wann PR vs SRK vs SAFT verwenden?

PR für C5+-Kohlenwasserstoffe, Erdgas, CO₂, nichtpolare Mischungen — petrochemisches Workhorse. SRK für leichte Kohlenwasserstoffe (C1-C4) und Süßgas. CPA (SRK+Assoziation) für Wasser, Alkohole, Glykole (Statoil, Linde). PC-SAFT für Polymere, Tenside, fortgeschrittene Petrochemie. GERG-2008 für Erdgas mit hoher Genauigkeit ISO 20765-2 (OGE, Gascade). Für Wasserdampf H₂O: IAPWS-IF97 spezifisch.

Was sind Mixing Rules in PR?

Für Mischungen werden Van-der-Waals-Regeln verwendet: a_mix = ΣΣ x_i·x_j·(a_i·a_j)^0,5·(1-k_ij), b_mix = Σ x_i·b_i. Parameter k_ij (binäre Wechselwirkung) wird aus experimentellen VLE-Daten regressiert. Tabelliert in HYSYS, ProMax, Multiflash. Beispiele: k_ij(CH₄-C₂H₆)≈-0,003, k_ij(CO₂-CH₄)≈0,1. Wong-Sandler 1992 und Huron-Vidal 1979 sind erweiterte Regeln für polare.

Welcher Peneloux-Volumenshift?

Die Peneloux-Korrektur (1982) löst den 6-12%-Flüssigkeitsdichtefehler von reinem PR: V_korrigiert = V_PR - c, wobei c ein tabellierter Parameter pro Substanz ist. Flüssigfehler sinkt auf 1-3%. Magoulas-Tassios 1990 erweitert c(T). Implementiert in Aspen HYSYS (PR-VC method) und ProMax. Entscheidend für LNG, Flüssiggas, überkritische CO₂-Sequestrierung, wo Flüssigkeitsdichte Schlüssel ist.

Kommentare (0)

Noch keine Kommentare. Sei der Erste!