IAPWS-IF97水物性计算Fortran90程序

IAPWS-IF97 水物性计算 Fortran90 程序

代码开源，版权所有，引用请标明出处（2013/3/27 修改）

本文给出了基于 IAPWS-IF97 的水物性计算 Fortran90 语言程序代码，程序编写为动态

链接库(dll)形式，以便于调用。下面具体讲述。

如上图所示，水与水蒸气可分为 5 个区，但本程序只适用于 4 个区，即图中的第

1、2、3、4 区，不包括第 5 区。压力与温度范围为：

程序的输入变量为压力 P 与温度 T，单位分别为 Pa 与 K。输出变量(12 个)如下：

----------------------------------------------------------------------------------------------------------------

v u h s Cp Cv w β κ λ μ σ

单位如前所述。

现给出具体的调用方法：

给定压力 P(Pa)与温度 T(K)，计算物性，程序自己判断所属区域：

call h2o(P,T,prop)

若想得到给定压力 P(Pa)或温度 T(K)下的饱和态物性，则可如下调用：

a. 给定压力 call PsatW(P,T)或给定温度 call TsatW(P,T)

b. 饱和液态 call h2o(P,0.999*T,prop(12))或饱和汽态 call h2o(P,1.001*T,prop(12))

当然还有其他方法，读者详细了解程序内容后，请自行发现。

这里给出获得动态链接库的方法：

1. 新建一个动态链接库工程

2. 加入附录所示源程序文件

附录：源程序(19个文件)

【1】. const.f90

real*8 :: Rg=0.461526E3 !J/kg/K

real*8 :: Tc=647.096 !K

real*8 :: Pc=22.064E6 !Pa

real*8 :: rhoc=322.0 !kg/m^3

end module const

【2】. reg1.f90

module reg1

implicit none

integer :: I(34),J(34)

real*8 :: n(34)

0 ,&

0 ,&

0 ,&

1 ,&

1 ,&

1 ,&

1 ,&

1 ,&

1 ,&

2 ,&

3 ,&

4 ,&

4 ,&

4 ,&

5 ,&

8 ,&

8 ,&

21 ,&

30 ,&

31 ,&

32 &

/

data J /&

-2 ,&

-1 ,&

1 ,&

2 ,&

3 ,&

4 ,&

3 ,&

-3 ,&

0 ,&

1 ,&

3 ,&

17 ,&

-4 ,&

0 ,&

6 ,&

-5 ,&

-31 ,&

-38 ,&

-39 ,&

-40 ,&

-41 &

/

data n /&

0.14632971213167 ,&

3.3855169168385 ,&

-0.95791963387872 ,&

0.15772038513228 ,&

-0.016616417199501 ,&

8.1214629983568E-04 ,&

2.8319080123804E-04 ,&

-6.0706301565874E-04 ,&

-0.032529748770505 ,&

-0.021841717175414 ,&

-5.283835796993E-05 ,&

-4.7184321073267E-04 ,&

-8.5205128120103E-10 ,&

-2.2425281908E-06 ,&

-6.5171222895601E-07 ,&

-1.4341729937924E-13 ,&

-4.0516996860117E-07 ,&

-1.2734301741641E-09 ,&

-1.7424871230634E-10 ,&

-6.8762131295531E-19 ,&

1.4478307828521E-20 ,&

2.6335781662795E-23 ,&

module reg2

implicit none