IF97

IF97 is the high-speed package of IAPWS-IF97 in Rust. 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

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

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

The following 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)  

Release Notes

• 1.2.0 - The multi-step method unleashes the full power of the compiler optimizations while using powi() with the for loop

  • Add the optional parameter of region to computer the properties of the specified region quickly

• 1.0.9 - The shortest addition chain computes integer powers of a number to provide the stable speed performance under the various hardware and software platforms

Usage

Install the crate

cargo add if97

The type of functions are provided in the if97 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 if97::*;
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}");   
}

Properties