NACH OBEN

Software für die Industrie-Formulation IAPWS-IF97 für Wasser und Wasserdampf

1. Kurze Information zur IAPWS-IF97

Eine detaillierte Beschreibung der IAPWS-IF97 Referenz-Zustandsgleichung für Wasser und Wasserdampf ist hier zu finden. Informationen über den Einsatz unter den verschiedenen Windows Betriebssystemen und Excel Versionen sowie über die Integration der 32-Bit oder 64-Bit Fortran DLL mit Hilfe von Excel Add-Ins und die gleichfalls verfügbare 32-Bit .Net DLL werden im Abschnitt 2.4 der Softwarebeschreibung gegeben.

Dort findet man auch detaillierte Informationen zu dem aktuellen Vertriebsweg über das Ingenieurbüro Dr.-Ing. N. Kurzeja.

Die IAPWS Industrie-Formulation IF97 besteht aus einem Satz von Gleichungen für verschiedene Zustandsbereiche, wobei insgesamt folgender Gültigkeitsbereich abgedeckt wird:

0 °C ≤ t ≤ 800 °C,   p ≤ 1000 bar (100 MPa)

800 °C < t ≤ 2000 °C,   p ≤ 500 bar (50 MPa)

zww_soft_if97_p1

Struktur und Bereiche der IAPWS-IF97.

Die obige Abbildung zeigt die fünf Bereiche, in die der gesamte Gültigkeitsbereich der IAPWS-IF97 aufgeteilt ist. Für die Bereiche 1 und 2 existiert jeweils eine Fundamentalgleichung der spezifischen freien Enthalpie g(p,T), für den Bereich 3 eine Fundamentalgleichung der spezifischen freien Energie f(ρ,T). Die Sättigungslinie, die dem Bereich 4 entspricht, wird durch eine Dampfdruckgleichung ps(T) abgedeckt. Der Hochtemperaturbereich 5 wird ebenfalls durch eine Fundamentalgleichung der Form g(p,T) beschrieben. Diese fünf Gleichungen, die in der Abbildung eingerahmt sind, bilden die sog. Basisgleichungen.

Zusätzlich zu den Basisgleichungen wurden sog. Rückwärtsgleichungen für die Bereiche 1-4 entwickelt. Diese Rückwärtsgleichungen existieren für die folgenden Kombinationen von Zustandsgrößen: Für die Bereiche 1 und 2 als Gleichungen der Form T(p,h), T(p,s), und p(h,s), für den Bereich 3 als Gleichungen der Form T(p,h), v(p,h), T(p,s), v(p,s), p(h,s) und v(p,T). Für den gesamten Bereich 4 ist die Rückwärtsgleichung eine Gleichung für die Sättigungstemperatur in der Form Ts(p) und ps(T). Beide Gleichungen sind mathematisch vollkommen konsistent zu einander. Für den technisch wichtigsten Teil des Bereiches 4 (ss’’ (623.15 K)) gibt es eine Sättigungstemperaturgleichung der Form Ts(h,s). Das obige Bild zeigt die Zuordnung der grau unterlegten Rückwärtsgleichungen zu den einzelnen Bereichen IAPWS-IF97, die Basisgleichungen sind eingerahmt.

Mit Hilfe der Rückwärtsgleichungen lassen sich alle Zustandsgrößen in Abhängigkeit der Variablen (p,h), (p,s) und (h,s) im Bereich 3 auch bei Vorgabe von (p,T) ohne Iteration, und damit sehr schnell berechnen.

Eine ausführliche Beschreibung der IAPWS-IF97 findet man in dem Buch

Kretzschmar, H.-J., Wagner, W., International Steam Tables - Properties of Water and Steam Based on the Industrial Formulation IAPWS-IF97. Springer-Verlag (Berlin), 2019.

In diesem Buch werden auch Gleichungen für die Transportgrößen sowie für die Größen Oberflächenspannung, Dielektrizitätskonstante und Brechungsindex angegeben. Weiterhin enthält das Buch für eine sehr große Anzahl von Zustandsgrößen Druck-Temperatur-Diagramme mit Isolinien aller tabellierten und weiterer Zustandsgrößen einschließlich partieller Differenzialquotienten. Weitere Einzelheiten (Inhaltsverzeichnis, Probeseiten, etc.) findet man hier.

2. Software für die IAPWS-IF97

Auf Basis der IAPWS-IF97 einschließlich aller Rückwärtsgleichungen sowie der IAPWS-Gleichungen für die Transportgrößen wurde Software für die Berechnung von mehr als 25 Zustandsgrößen entwickelt. Bei der Entwicklung dieser Software wurde ein besonderer Wert auf eine optimale Programmierung im Hinblick auf kurze Rechenzeiten gelegt.

Bei der Anwendung der Software braucht nicht bekannt zu sein, in welchem Bereich der IAPWS-IF97 die zu berechnende Größe liegt. Bei Vorgabe der Eingangsgrößen ermittelt die Software selbsttätig, welche der Bereichsgleichungen anzuwenden ist.

Mit der Software IAPWS-IF97BO zur IAPWS-IF97 lassen sich folgende thermodynamische Zustandsgrößen sowie Transportgrößen und einige andere Zustandsgrößen von Wasser und Wasserdampf berechnen.

2.1 Thermodynamische Zustandsgrößen

Folgende thermodynamische Zustandsgrößen können mit den entsprechenden Gleichungen der IAPWS-IF97 berechnet werden:

Tabelle 1  Berechenbare thermodynamische Größen

Symbol

 Thermodynamische Größen

p

 Druck

T

 Temperatur

v

 spezifisches Volumen

ρ

 Dichte

h

 spezifische Enthalpie

s

 spezifische Entropie

cp

 spezifische isobare Wärmekapazität

cv

 spezifische isochore Wärmekapazität

w

 Schallgeschwindigkeit

u

 spezifische innere Energie

f

 spezifische freie Energie

g

 spezifische freie Energie

f*

 Fugazität

z

 Realgasfaktor

x

 Dampfgehalt

αv

 isobarer Volumenausdehnungskoeffizient, αv = v1(v/T)p

αp

 Relativer Druckkoeffizient, αp = p1(p/T)v

βp

 isochorer Spannungskoeffizient, βp = – p1p/v)T

κ

 Isentropenexponent, κ = – (v/p) (p/v)s

κT

 isotherme Kompressibilität, κT = – v1(v/p)T

µ

 Joule‑Thomson Koeffizient, µ = (T/p)h

(v/p)h

 partielle Ableitung des Volumens nach dem Druck bei konstanter Enthalpie

(ρ/p)h

 partielle Ableitung der Dichte nach dem Druck bei konstanter Enthalpie

(v/h)p

 partielle Ableitung des Volumens nach der Enthalpie bei konstantem Druck

(ρ/h)p

 partielle Ableitung der Dichte nach der Enthalpie bei konstantem Druck

(h/p)T

 partielle Ableitung der Enthalpie nach dem Druck bei konstanter Temperatur

Diese Zustandsgrößen können im gesamten Gültigkeitsbereich der IAPWS-IF97 (siehe Abschn. 1) berechnet werden. In der Region 4 (Zwei-Phasengebiet) können die Zustandsgrößen νρhsufg und x auch innerhalb des Zwei-Phasengebiets berechnet werden. Die anderen Zustandsgrößen können nur an den Phasengrenzen (Siedelinie und Taulinie) berechnet werden, weil sie innerhalb des Zwei-Phasengebiets nicht definiert sind.

2.2 Transportgrößen und weitere Zustandsgrößen

Mit der Software IAPWS-IF97BO lassen sich folgende Transportgrößen und weitere Zustandsgrößen berechnen:

Tabelle 2  Berechenbare Transportgrößen und
weitere Zustandsgrößen

Symbol

 Transportgrößen

η

 dynamische Viskosität

ν

 kinematische Viskosität, ν = η/ρ 

λ

 Wärmeleitfähigkeit

Pr

 Prandtl-Zahl, Pr = η cp/λ 

a

 Temperaturleitfähigkeit, a = λ /(ρcp)

ε

 statische Dielektrizitätskonstante

n

 Brechungsindex

σ

 Oberflächenspannung

Die international vereinbarten Gleichungen zur Berechnung dieser Größen, die nicht zur IAPWS-IF97 gehören, sind in dem im Abschnitt 1 erwähntem Buch International Steam Tables zusammengestellt und beschrieben. Im Einzelnen sind dies das IAPWS Release IAPWS R12-08 für die Viskosität, IAPWS R15-11 für die Wärmeleitfähigkeit, IAPWS R8-97 für die statische Dielektrizitätskonstante, IAPWS R9-97 für den Brechungsindex und IAPWS R1-76(2014) für die Oberflächenspannung, siehe auch http://www.iapws.org/ unter "Releases".

Diese Gleichungen sind Funktionen der Temperatur und der Dichte (nicht Druck). Das hat zur Folge, dass bei der Vorgabe von Temperatur und Druck als Eingangsgrößen Iterationen notwendig sind, die in der Software IAPWS-IF97BO durchgeführt werden. Die Oberflächenspannung ist eine reine Temperaturfunktion. Bis auf die Dielektrizitätskonstante ε können alle übrigen Eigenschaften im gesamten Gültigkeitsbereich der IAPWS-IF97 berechnet werden.

Die Gleichung für ε ist nur bis 873,15 K gültig, kann aber bis 1073,15 K sinnvoll extrapoliert werden. Die Gleichung für σ(T) gilt nur für den Bereich 4 der IAWPS-IF97.

2.3 Mögliche Kombination der Eingangsgrößen für die Berechnung der verschiedenen Zustandsgrößen

2.3.1 Regionen 1-5

Für die Regionen 1-3 und 5 (Ein-Phasengebiet) und Region 4 (Zwei-Phasengebiet) können alle in den Tabellen 1 und 2 aufgeführten Größen als Funktion der in Tabelle 3 aufgelisteten Kombinationen von Eingangsgrößen berechnet werden. Die Berechnung der Zustandsgrößen für andere Kombinationen von Eingangsgrößen als (p,T) erfolgt über entsprechende Iterationen; für die Region 5 gibt es keine Rückwärtsgleichungen. In der Region 4 lassen sich keine Berechnungen als Funktion von (p,T) durchführen, weil in dieser Region p und T nicht unabhängig voneinander sind.

Tabelle 3  Mögliche Kombinationen
der Eingangsgrößen für die Regionen
1 bis 5; für Region 4 kann (p,T) nicht
als Eingangsgrößen benutzt werden,
siehe den Text.

(p,T)

(T,h)

(ν,h)

(h,s)

(p,h)

(T,s)

(ν,s)

(p,s)

(T,ν)

 

(p,ν)

(T,ρ)

 

Für die in Tabelle 3 aufgeführten Paare von Eingangsgrößen können die wichtigsten Größen direkt berechnet werden. Die Größen, für die es für die gewählten Kombinationen keine direkten Funktionen in der Software gibt, können durch Kombinationen von in der Software vorhandenen Funktionen ermittelt werden.

2.3.2 Berechnung von Zustandsgrößen nur in der Region 4 (Zwei-Phasengebiet)

Für die Berechnung von Zustandsgrößen in Region 4 (Zwei-Phasengebiet (0 < x < 1) [Sättigungsdruck ps, Sättigungstemperatur Ts, Siedelinie [('), x = 0], Taulinie [("), x = 1]), die in der linken Spalte der Tabelle 4 aufgelistet sind, existieren in der Software direkte Funktionen für die in der zweiten Spalte aufgeführten Eingangsgrößen.

Tabelle 4..Berechenbare Zustandsgrößen für den gesamten Bereich 4
(Siedelinie, Taulinie und innerhalb des Zweiphasengebiets) mit den in
der Software vorhandenen Funktionen

Berechenbare

Größen

Eingangsgrößen1

Erläuterung

ps

Thsρx

1 Zusätzlich zu den aufgeführten Eingangsgrößen können die Berechnungen auch für die in Tabelle 3 aufgeführten Paare von Eingangsgrößen [außer (p,T)] erfolgen; die Software enthält die entsprechenden Funktionen. Die Software ermittelt, ob der Zustandspunkt in Region 4 liegt.

Ts

phsρ

vh, s

(T,x), (p,x)

νf, g

(T,x), (p,x)

(p,h), (p,s), (p,v), (T,ρ)

Die Zustandsgrößen, die in den Tabellen 1 und 2 aufgeführt und nicht in der linken Spalte von Tabelle 4 enthalten sind, können nicht innerhalb von Region 4 (Zwei-Phasengebiet, 0 < x < 1) berechnet werden, sondern nur auf den Phasengrenzen Siedelinie (x = 0) und Taulinie (x = 1) als Funktion von T oder von p; diese Größen sind innerhalb des Zwei-Phasengebiets nicht definiert.

2.4 IAPWS-IF97 DLL and Excel Add-In für die Berechnung thermopyhsikalischer Eigenschaften

Zur Einbindung der IAPWS-IF97 in Windows-basierte Standardprogramme enthält das Softwarepaket IAPWS-IF97BO die Dynamic Link Library (DLL) IF97BO.dll, die mit der Intel Parallel Studio XE 2015 Composer Edition for Fortran (Version 15.07.287) optimiert und kompiliert worden ist. Die DLL enthält zahlreiche Funktionen, die die Berechnung aller aufgelisteten Zustandsgrößen für alle weiter oben genannten Kombinationen der Eingangsgrößen ermöglichen.

In der Software kann gewählt werden, ob die gewünschten Zustandsgrößen mit Hilfe der in der IAPWS-IF97 enthaltenen Rückwärtsgleichungen berechnet werden sollen oder nur die Basisgleichungen mit Iterationen benutzt werden sollen (falls Iterationen für die Berechnung der entsprechenden Zustandsgröße notwendig sind). In diesem Fall werden die Rückwärtsgleichungen zur Berechnung der Startwerte benutzt, so dass diese Iterationen deutlich schneller sind als ohne diese Startwerte. Der Aufruf erfolgt jeweils über einfache Funktionsnamen, die sich aus Abkürzungen für die zu berechnenden Zustandsgröße und die notwendigen Eingangsgrößen ergeben. So steht zum Beispiel zur Berechnung der spezifischen Enthalpie h für gegebene Werte des Drucks p und der Temperatur T mit Verwendung der Rückwärtsgleichungen die Funktion HBPT und ohne Verwendung der Rückwärtsgleichungen die Funktion HOPT zur Verfügung.

Zusätzlich enthält das in 32-Bit oder 64-Bit verfügbare IAPWS-IF97BO Softwarepaket ein Microsoft Excel-Add-In, das eine einfache Integration der DLL in Microsoft Excel ermöglicht. Da durch die Installation des IAPWS-IF97BO Excel Add-Ins alle mit der DLL exportierten Funktionen dem Standard-Funktionsumfang von Microsoft Excel hinzugefügt werden, können die gewünschten thermodynamischen Zustandsgrößen und Transportgrößen entweder durch die direkte Eingabe der im Manual aufgeführten Formel in ein Excel-Tabellenblatt oder sehr komfortabel mit Hilfe des Funktionsassistenten von Microsoft Excel aufgerufen werden.

zww_soft_if97_p3

Aufruf der IAPWS-IF97 Funktion VOPT mit Hilfe des Excel Funktionsassistenten

zww_soft_if97_p4

Eingabe der Argumente für den Funktionsaufruf VOPT

Ein einfacher Zugriff auf die vom IAPWS-IF97BO Softwarepaket angebotenen Funktionsaufrufe für die thermodynamischen Zustandsgrößen und die Transportgrößen ist natürlich auch innerhalb Excel über für eigene Zwecke geschriebene VBA-Skripte möglich, die dann als Makro aufgerufen werden können.

Darüber hinaus enthält das IAPWS-IF97BO Softwarepaket die Import-Bibliothek IF97BO.lib, mit deren Hilfe die DLL in anwenderspezifische Fortran- und unmanaged C++-Anwendungen integriert werden kann.

Alle in der Software IAPWS-IF97BO enthaltenen Funktionen sind in dem mitgelieferten MANUAL IAPWS-IF97BO_x86/x64.PDF detailliert beschrieben.

Als Beispiel zeigt der folgende Screenshot des Referenztabellenblatts IF97BO.xlsm die Ergebnisse für die Berechnung der Zustandsgrößen für die Temperatur T = 300 K und den Druck p = 1 MPa. Für T = 300 K wurden auch alle Größen auf der Siede- und Taulinie berechnet. Für den Dampfgehalt von x = 0,5 wurden außerdem alle Zustandsgrößen berechnet, die im Zweiphasengebiet definiert sind. Dieses Excel Referenztabellenblatt enthält alle Funktionsaufrufe, die das Softwarepaket IAPWS-IF97BO zur Verfügung stellt, und kann daher als schnelles und einfaches Berechnungstool für die thermodynamischen Zustandsgrößen und Transportgrößen von Wasser und Wasserdampf verwendet werden. Natürlich kann auch ein benutzerspezifisches Excel Tabellenblatt mit den gewünschten Funktionen erstellt werden, oder jede gewünschte Funktion der IAPWS-IF97BO Software kann auch innerhalb eines VBA-Skripts aufgerufen werden. https://www.thermo.ruhr-uni-bochum.de/thermo/wagner.html.de

zww_soft_if97_p2

Screenshot des Excel Tabellenblatts IF97BO.xlsm als Referenz für die mit IAPWS-IF97BO berechenbaren Zustandsgrößen.

Die 32-Bit Version der IAPWS-IF97BO DLL kann mit Hilfe des mitgelieferten Excel Add-Ins IF97BO.xla (für Excel 2003) bzw. IF97BO.xlam (für Excel 2007 – 2021/Microsoft 365) unter den verschiedenen Windows Betriebssystemen [Windows 2000 oder XP (32 Bit), Windows Vista, 7, 8 und 10 (32 Bit/64 Bit)] in den 32-Bit Versionen von Microsoft Excel (2003 – 2021/Microsoft 365) benutzt werden.

Die IAPWS-IF97BO Software ist auch in einer 64-Bit Version erhältlich, die auf den 64-Bit Betriebssystemen von Windows XP bis Windows 11 mit Hilfe des mitgelieferten 64-Bit Excel Add-Ins IF97BO.xlam in die 64-Bit Versionen von Excel 2010 – Excel 2021/Microsoft 365 integriert werden kann. In einer 32-Bit Version von Excel entwickelte Arbeitsblätter mit IAPWS-IF97BO Funktionsaufrufen können damit ohne Änderungen in einer 64-Bit  Version von Excel weiterverwendet werden.

Zusätzlich ist auch eine 32-Bit .NET DLL von IAPWS-IF97BO erhältlich, die in die .NET Programmiersprachen (Visual C#, Visual C++ oder Visual Basic) von Microsoft Visual Studio eingebunden werden kann.

Die IAPWS-IF97BO Software ist nicht kostenlos und wird vertrieben über das Ingenieurbüro Dr.-Ing. N. Kurzeja.

 

Bitte beachten Sie, dass ein entsprechendes Softwarepaket auch für die Referenzzustandsgleichung GERG-2008 und für über 95 reine Fluide über die Software FluidCal verfügbar ist. Weitergehende Informationen zum gesamten Softwareangebot finden Sie unter https://www.thermo.ruhr-uni-bochum.de/thermo/wagner.html.de.

Kontakt:

Ingenieurbüro Dr.-Ing. N. Kurzeja
feos-software
Hevener Mark 26
D-58455 Witten
Germany
phone: +49 (0)2302 279849
mobile: +49 (0)175 224 8598

Norbert.Kurzeja@ruhr-uni-bochum.de

Das Ingenieurbüro arbeitet  auf dem Gebiet der Entwicklung und des Vertriebs thermodynamischer Software sowie in der Forschung zu kalorimetrischen Messungen des Brennwerts und zum Verhalten reinen fluider Stoffe im sog. kritischen Gebiet. Beispiele dafür sind die Publikationen

Kurzeja, N., Span, R. Accurate measurements of the gross calorific value of methane by the renewed GERG calorimeter.  Publiziert in der Featured Collection of J. Phys. Chem. Ref. Data 48, 043103 (2019). https://doi.org/10.1063/1.5110054

Wagner, W., Kurzeja, N. Letter to the editor. Int. J. Thermophys. 32 (2011), 549-552. https://doi.org/10.1007/s10765-011-0945-7

Kurzeja, N., Tielkes, Th., Wagner, W. The nearly classical behavior of a pure fluid on the critical isochore very near the critical point under influence of gravity. Int. J. Thermophysics 20 (1999), 531-561. https://doi.org/10.1023/A:1022657121329

Wagner, W., Kurzeja, N., Pieperbeck, B. The thermal behavior of fluid substances in the critical region - experiences from recent pρT measurements on SF6 with a multi-cell apparatus. Fluid Phase Equilibria 79 (1992), 151-174. https://doi.org/10.1016/0378-3812(92)85127-T