在该包中，可以计算 36 个热力学、输运和其他属性。

实现了以下 12 组输入

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

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

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

(h,s)  

用法

安装包

cargo add seuif97

包中提供了函数类型

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

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

• 第一个、第二个输入参数(f64)：输入属性对
• 第三个和第四个输入参数
  • 第三个：计算属性的属性 ID - o_id
  • 第四个 option 参数：IAPWS-IF97 的区域
• 返回(f64)：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

示例

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}");   
}

C 绑定

构建动态链接库

• cdecl

cargo build -r --features cdecl

• stdcall：Win32 API 函数

cargo build -r --features stdcall

在 ./dynamic_lib/ 提供了方便编译的动态链接库

• seuif97.dll：Windows64 和 Windows32

• libseuif97.so：Linux64

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);

示例

• ./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;
}

Rankine 循环分析示例

• Python，C++ 和 Modelica 中的 Rankine 循环稳态仿真器

Python 绑定

安装

pip install seuif97

示例

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}")

示例

• T-S 图

• H-S 图

• 蒸汽轮机膨胀的 H-S 图

• Python 中的 Rankine 循环稳态仿真器的混合版

属性

引用方式

• 程茂华. (2023). 高速 IAPWS-IF97 包的 Rust 实现：SEUIF97 (1.0.9). Zenodo. https://doi.org/10.5281/zenodo.8246380