The thermophysical properties of water and steam are very important to researchers and engineers in various domains (energy, power, chemistry, materials ...).
Usually, researchers and engineers rely on the NIST Refprop to calculate the properties. However, Refprop is commercial and not open-source. Therefore, I wrote this code to provide an alternative to the Refprop.
Program Name: Ultra-high-resolution Water and Steam Thermophysical Property Calculation
Purpose: Calculate the very-high-resolution thermophysical properties at given point(s)
License: MIT
Technical Reference: IAPWS IF97
For library use, CMake builds the thermox_if97 shared library and its small,
namespaced public API:
cmake -S . -B build
cmake --build build
ctest --test-dir buildInclude thermox_if97.h and call thermox_if97_state_pt,
thermox_if97_state_ph, thermox_if97_state_ps, or
thermox_if97_saturation_p. The saturation query returns the coexisting
saturated-liquid and saturated-vapor states at a pressure from the minimum
IF97 saturation pressure through the Region 1/2 boundary (approximately
16.5 MPa). The API uses SI units and explicit status and phase values. The
optional file-based CLI can be built with
-DTHERMOX_IF97_BUILD_LEGACY_CLI=ON.
You need a C compiler to build.
- For Microsoft Windows users, mingw-w64 is required.
- For GNU/Linux Distro or other *nix users, the GNU Compiler Collections, known as gcc, is recommended.
- For macOS users, clang is easy to install and use (brew is not needed to install clang on macOS).
- Use
gitto clone this code:git clone https://github.com/zhenrong-wang/water-steam-if97.git - Build command example:
gcc *.c -o my-if97.exe -lm
Note: the -lm may not be valid for Windows or macOS. It is necessary for GNU/Linux distros.
Suppose you have followed the steps above and built the my-if97.exe in the source code folder.
cd test../my-if97.exe(UNIX-like OS) or..\my-if97.exe(Microsoft Windows)
If you see the output below, congrats, you've built the binary successfully.
...
1 25000000.00000000 873.00000000
1 23960000.00000000 811.00000000
1 6000000.00000000 811.00000000
1 6000000.00000000 573.00000000
# Calculation finished. Please check _properties.dat for results.
# Press any key to exit.
@ Any problems found, please contact the author.
@ Zhenrong Wang, zhenrongwang@live.com.
Suppose the working directory is /home/amber/water-steam/
Suppose the absolute path of the executable is /home/amber/bin/my-if97.exe
$ cd /home/amber/water-steam/$ vim _input.dat# Create an input file named '_input.dat' with some input data (See Section 3.3 below)$ /home/amber/bin/my-if97.exe
Suppose the working directory is C:\Users\amber\water-steam\
Suppose the absolute path of the executable is C:\Users\amber\bin\my-if97.exe
- Open a Command Prompt Window
cd C:\Users\amber\water-steam\notepad _input.datC:\Users\amber\bin\my-if97.exe
For the _input.dat file, you need to input the data in lines. A single line refers to a point. Strict format:
POINT_TYPE(int),PARAM1(float/double),PARAM2 (float/double)
Please separate the params by a comma ','
POINT_TYPE:
- 1 Pressure and Temperature
- 2 Pressure and Density
- 3 Pressure and Specific Internal Energy
- 4 Pressure and Specific Enthalpy
- 5 Pressure and Specific Entropy
- 6 Temperature and Density
- 7 Temperature and Specific Internal Energy
- 8 Temperature and Specific Enthalpy
- 9 Temperature and Specific Entropy
- 10 Specific Enthalpy and Specific Entropy
- 11 Pressure and Steam Dryness
- 12 Temperature and Steam Dryness
Example: 1,100000,300
Example: 1,1e5,3e2
Units:
- Pressure(p): Pa
- Temperature(t): K
- Density(r): kg/m3
- Specific Internal Energy(u): J/kg
- Specific Enthalpy(h): J/kg
- Specific Entropy(s): J/(kg.K)
The program writes out a file named '_properties.dat' in the working directory.
This code was written years ago, and I didn't continuously develop it. Therefore, bugs may occur.
If you are interested in this project, please submit issues to this repo. I'd be glad to communicate on any issues.
Or, you can also email me: zhenrongwang@live.com