dependences

• Support for both python branch:

  * python 2.7
  * python 3.4 or later

• Numpy-scipy: library with mathematic and scientific tools

install

In debian you can find in official repositories in jessie, testing and sid. In ubuntu it’s in official repositories from ubuntu saucy (13.10). In other system you can install using pip:

pip install iapws

or directly from the github repository:

pip install git+https://github.com/jjgomera/iapws.git

This is the recommended option to have the latest version.

documentation

To see the full documentation of package, see readthedocs

For a rapid usage demostration, see this examples

IAPWS-IF97 (see full documentation)

from iapws import IAPWS97
sat_steam = IAPWS97(P=1,x=1)                #saturated steam with known P
sat_liquid = IAPWS97(T=370, x=0)            #saturated liquid with known T
steam = IAPWS97(P=2.5, T=500)               #steam with known P and T
print(sat_steam.h, sat_liquid.h, steam.h) #calculated enthalpies

IAPWS-95 (see full documentation)

from iapws import IAPWS95
sat_steam = IAPWS95(P=1,x=1)                #saturated steam with known P
sat_liquid = IAPWS95(T=370, x=0)            #saturated liquid with known T
steam = IAPWS95(P=2.5, T=500)               #steam with known P and T
print(sat_steam.h, sat_liquid.h, steam.h) #calculated enthalpies

For calculation of multiple states is possible use multiprocessing to speed up calculation, 6x in my laptop. Method valid too for D2O

from iapws import IAPWS95
from numpy import arange
from time import time

x = arange(0, 1.01, 0.01)

def fi(x):
    return IAPWS95(P=20.8, x=x)

start = time()
for xi in x:
    fi(xi)
print(f'Without multiprocessing: {time() - start}')

start = time()
states = IAPWS95.from_list("P", 20.8, "x", x)
print(f'With multiprocessing: {time() - start}')

IAPWS-17 for Heavy water (see full documentation)

from iapws import D2O
sat_liquid = D2O(T=370, x=0)            #saturated liquid with known T
print(sat_liquid.h)                   #calculated enthalpy

IAPWS-06 for Ice Ih (see full documentation)

from iapws import _Ice
ice = _Ice(273.15, 0.101325)            #Ice at normal melting point
print(ice["rho"])                     #Calculated density

IAPWS-08 for seawater (see full documentation)

from iapws import SeaWater
state = SeaWater(T=300, P=0.101325, S=0.001)    #Seawater with 0.1% Salinity
print(state.cp)                                 # Get cp