Crate seuif97

Expand description

SEUIF97

docs.rs PyPI

This is the Rust implementation of the high-speed IAPWS-IF97 package seuif97 with C and Python binding. It is suitable for computation-intensive calculations，such as heat cycle calculations, simulations of non-stationary processes, real-time process monitoring and optimizations.

Through the high-speed package, the results of the IAPWS-IF97 are accurately produced at about 5-20x speed-up compared to using the powi() of the Rust standard library in the forloop directly when computing the basic equations of Region 1,2,3.

The Fast Methods

The multi-step method unleashes the full power of the compiler optimizations while using powi() with the for loop
The recursive method computes the polynomial values of the base variable and its derivatives

In the package, 36 thermodynamic, transport and further properties can be calculated.

The following 12 input pairs are implemented:

(p,t) (p,h) (p,s) (p,v) 

(t,h) (t,s) (t,v) 

(p,x) (t,x) (h,x) (s,x) 

(h,s)

Usage

Install the crate

cargo add seuif97

The type of functions are provided in the package:

struct o_id_region_args {
   o_id: i32,
   region: i32,
}

fn<R>(f64,f64,R) -> f64
where
    R: Into<o_id_region_args>,

the first,second input parameters(f64) : the input propertry pairs
the third and fourth input parametes
:
- the third : the property ID of the calculated property - o_id
- the fourth option parameter: the region of IAPWS-IF97
the return(f64): the calculated property value of o_id

pt<R>(p:f64,t:f64,o_id_region:R)->f64
ph<R>(p:f64,h:f64,o_id_region:R)->f64
ps<R>(p:f64,s:f64,o_id_region:R)->f64
pv<R>(p:f64,v:f64,o_id_region:R)->f64

th<R>(t:f64,h:f64,o_id_region:R)->f64
ts<R>(t:f64,s:f64,o_id_region:R)->f64
tv<R>(t:f64,v:f64,o_id_region:R)->f64

hs<R>(h:f64,s:f64,o_id_region:R)->f64

px(p:f64,x:f64,o_id:i32)->f64
tx(p:f64,x:f64,o_id:i32)->f64
hx(h:f64,x:f64,o_id:i32)->f64
sx(s:f64,x:f64,o_id:i32)->f64

Example

use seuif97::*;
fn main() {
    
    let p:f64 = 3.0;
    let t:f64= 300.0-273.15;
   
    let h=pt(p,t,OH);
    let s=pt(p,t,OS);
    // set the region
    let v=pt(p,t,(OV,1));
    println!("p={p:.6} t={t:.6} h={t:.6} s={s:.6} v={v:.6}");   
}

The C binding

Building the dynamic link library

cdecl

cargo build -r --features cdecl

stdcall: Win32 API functions

cargo build -r --features stdcall

The convenient compiled dynamic link libraries are provided in the ./dynamic_lib/

seuif97.dll: Windows64 and Windows32
libseuif97.so: Linux64

The functions in C

double pt(double p,double t,short o_id);
double ph(double p,double h,short o_id);
double ps(double p,double s,short o_id);
double pv(double p,double v,short o_id);

double tv(double t,double v,short o_id);
double th(double t,double h,short o_id);
double ts(double t,double s,short o_id);

double hs(double h,double s,short o_id);

double px(double p,double x,short o_id);
double tx(double t,double x,short o_id);
double hx(double h,double x,short o_id);
double sx(double s,double x,short o_id);

Example

./demo_using_lib/: C, Python, C#, Excel VBA, Java, Fortran

#include <stdlib.h>
#include <stdio.h>
#include <string.h>

#define OH 4
#define OS 5

extern double pt(double p,double t,short o_id);

int main(void)
{
    double p = 16.0;
    double t = 530.0;
    double h = pt(p, t, OH);
    double s = pt(p, t, OS);
    printf("p,t %f,%f h= %f s= %f\n", p, t, h, s);
    return EXIT_SUCCESS;
}

The Example of Rankine Cycle Analysis

The Rankine Cycle Steady-state Simulator in Python，C++ and Modelica

The Python binding

Install

pip install seuif97

Example

from seuif97 import *

OH=4

p=16.0
t=535.1
# ??(in1,in2,o_id)
h=pt(p,t,OH)
# ??2?(in1,in2)
s=pt2s(p,t)
print(f"p={p}, t={t} h={h:.3f} s={s:.3f}")

Examples

T-S Diagram

Properties

Propertry	Unit	Symbol	o_id	o_id(i32)
Pressure	MPa	p	OP	0
Temperature	°C	t	OT	1
Density	kg/m³	ρ	OD	2
Specific Volume	m³/kg	v	OV	3
Specific enthalpy	kJ/kg	h	OH	4
Specific entropy	kJ/(kg·K)	s	OS	5
Specific exergy	kJ/kg	e	OE	6
Specific internal energy	kJ/kg	u	OU	7
Specific isobaric heat capacity	kJ/(kg·K)	cp	OCP	8
Specific isochoric heat capacity	kJ/(kg·K)	cv	OCV	9
Speed of sound	m/s	w	OW	10
Isentropic exponent		k	OKS	11
Specific Helmholtz free energy	kJ/kg	f	OF	12
Specific Gibbs free energy	kJ/kg	g	OG	13
Compressibility factor		z	OZ	14
Steam quality		x	OX	15
Region		r	OR	16
Isobari cubic expansion coefficient	1/K	ɑv	OEC	17
Isothermal compressibility	1/MPa	kT	OKT	18
Partial derivative (∂V/∂T)p	m³/(kg·K)	(∂V/∂T)p	ODVDT	19
Partial derivative (∂V/∂p)T	m³/(kg·MPa)	(∂v/∂p)t	ODVDP	20
Partial derivative (∂P/∂T)v	MPa/K	(∂p/∂t)v	ODPDT	21
Isothermal throttling coefficient	kJ/(kg·MPa)	δt	OIJTC	22
Joule-Thomson coefficient	K/MPa	μ	OJTC	23
Dynamic viscosity	Pa·s	η	ODV	24
Kinematic viscosity	m²/s	ν	OKV	25
Thermal conductivity	W/(m.K)	λ	OTC	26
Thermal diffusivity	m²/s	a	OTD	27
Prandtl number		Pr	OPR	28
Surface tension	N/m	σ	OST	29
Static Dielectric Constant		ε	OSDC	30
Isochoric pressure coefficient	1/K	β	OPC	31
Isothermal stress coefficient	kg/m³	βp	OBETAP	32
Fugacity coefficient		fi	OFI	33
Fugacity	MPa	f*	OFU	34
Relative pressure coefficient	1/K	αp	OAFLAP	35

Cite as

Cheng Maohua. (2023). The Rust implementation of the high-speed IAPWS-IF97 package: SEUIF97 (1.0.9). Zenodo. https://doi.org/10.5281/zenodo.8246380

Structs

o_id_region_args
the paramters:
o_id: the propertry of id;
region: the region in the IAPWS-IF97,optional

Constants

OALFAP
αp - Relative pressure coefficient 1/K
OBETAP
βp - Isothermal stress coefficient, kg/m³
OCP
cp - Specific isobaric heat capacity kJ/(kg·K)
OCV
cv - Specific isochoric heat capacity kJ/(kg·K)
OD
ρ - Density kg/m³
ODPDT
(∂p/∂t)v - Partial derivative MPa/K
ODV
η - Dynamic viscosity Pa.s dv
ODVDP
(∂V/∂P)T - Partial derivative m³/(kg·MPa)
ODVDT
(∂V/∂T)p - Partial derivative m³/(kg·K)
OE
e - Specific exergy kJ/kg
OEC
δt - Isobaric cubic expansion coefficient 1/K
OF
f - Specific Helmholtz free energy kJ/kg
OFI
fi- Fugacity coefficient
OFU
f* - Fugacity MPa
OG
g - Specific Gibbs free energy kJ/kg
OH
h - Specific enthalpy kJ/kg
OIJTC
δt - Isothermal throttling coefficient kJ/(kg·MPa)
OJTC
μ - Joule-Thomson coefficient K/MPa joule
OKS
k - Isentropic exponent
OKT
kT - Isothermal compressibility 1/MPa
OKV
ν - Kinematic viscosity m²/s
OP
p - Pressure MPa
OPC
β - Isochoric pressure coefficient 1/K
OPR
Pr - Prandtl number
OR
r - Region
OS
s - Specific entropy kJ/(kg·K) )
OSDC
ε - Static Dielectric Constant
OST
σ - Surface tension N/m
OT
t - Temperature °C
OTC
λ - Thermal conductivity W/(m.K) tc
OTD
a - Thermal diffusivity m²/s
OU
u - Specific internal energy kJ/kg
OV
v - Specific Volume m³/kg
OW
w - Speed of sound m/s
OX
x - Steam quality
OZ
z - Compressibility factor

Functions

hs
hs(h,s,o_id) - the propertry of o_id (thermodynamic,transport,etc)
hs(h,s,(o_id,reg)) - the propertry of o_id in the region of reg(thermodynamic,transport,etc)
hx
hx(h,x,o_id) - the propertry of o_id(thermodynamic)
ph
ph(p,h,o_id) - the propertry of o_id (thermodynamic,transport,etc)
ph(p,h,(o_id,reg)) - the propertry of o_id in the region of reg(thermodynamic,transport,etc)
ps
ps(p,s,o_id) - the propertry of o_id (thermodynamic,transport,etc)
ps(p,s,(o_id,reg)) - the propertry of o_id in the region of reg(thermodynamic,transport,etc)
pt
pt(p,t,o_id) - the propertry of o_id (thermodynamic,transport,etc)
pt(p,t,(o_id,reg)) - the propertry of o_id in the region of reg(thermodynamic,transport,etc)
pv
pv(p,v,o_id) - the propertry of o_id(thermodynamic,transport,etc)
pv(p,v,(o_id,reg)) - the propertry of o_id in the region of reg(thermodynamic,transport,etc)
px
px(p,x,o_id) - the propertry of o_id (thermodynamic)
sx
sx(s,x,o_id) - the propertry of o_id(thermodynamic)
th
th(t,h,o_id) - the propertry of o_id(thermodynamic,transport,etc)
th(t,h,(o_id,reg)) - the propertry of o_id in the region of reg(thermodynamic,transport,etc)
ts
ts(t,s,o_id) - the propertry of o_id(thermodynamic,transport,etc)
ts(t,s,(o_id,reg)) - the propertry of o_id in the region of reg(thermodynamic,transport,etc)
tv
tv(t,v,o_id) - the propertry of o_id(thermodynamic,transport,etc)
tv(t,v,(o_id,reg)) - the propertry of o_id in the region of reg(thermodynamic,transport,etc)
tx
tx(t,x,o_id) - the propertry of o_id (thermodynamic)