程序的输入变量为压力 P
与温度 T，单位分别为 Pa 与 K。输出变量(12 个)如下：
----------------------------------------------------------------------------------------------------------------
比容 v 内能 u 比焓 h 比熵 s 定压比热容 Cp 定容比热容 Cv 声速 w
m^3/kg J/kg J/kg J/kg*K J/kg*K J/kg*K m/s
定压热膨胀率 β 等温压缩率 κ 导热系数 λ 动力粘度 μ 表面张力 σ
1/K 1/Pa W/m*K Pa*s N/m
----------------------------------------------------------------------------------------------------------------
输出变量是以一个数组的形式存储，例如数组为 prop(12)，数组与变量的对应顺序为：
v u h s Cp Cv w β κ λ μ σ
单位如前所述。
1. 使用方法：
将附录 19 个源文件处理成.f90 文件（必须是.f90, 处理成.for 需要排除大量语法错
误），然后加入到所需调用物性的工程中，调试即可。
2. 调用方法：
给定压力 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.9999999*T,prop(12)) 或 饱 和 汽 态 call
h2o(P,1.0000001*T,prop(12)) （原文误差大，修正）
附录：源程序(19个文件)
【1】. const.f90
module const
implicit none
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)
评论2