IAPWS-IF97水物性计算Fortran90程序.doc

IAPWS-IF97水物性计算Fortran90程序 本文给出了基于IAPWS-IF97的水物性计算Fortran90语言程序代码，程序编写为动态链接库(dll)形式，以便于调用。下面具体讲述。 如上图所示，水与水蒸气可分为5个区，但本程序只适用于4个区，即图中的第1、2、3、4区，不包括第5区。压力与温度范围为： 程序的输入变量为压力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 β κ λ μ σ 单位如前所述。 现给出具体的调用方法： 给定压力P(Pa)与温度T(K)，计算物性，程序自己判断所属区域： call h2o(P,T,prop(12)) 若想得到给定压力P(Pa)或温度T(K)下的饱和态物性，则可如下调用： 给定压力call PsatW(P,T)或给定温度call TsatW(P,T) 饱和液态call h2o(P,1.001*T,prop(12))或饱和汽态call h2o(P,0.999*T,prop(12)) 当然还有其他方法，读者详细了解程序内容后，请自行发现。 这里给出获得动态链接库的方法： 新建一个动态链接库工程 加入附录所示源程序文件 编译，获得最重要的下面两个文件：.lib与.dll 上面两个文件就可用于其他程序的物性调用 Wolfgang Wagner, Hans-Joachim Kretzschmar. Properties ofWater and Steam Based on the Industrial formulation IAPWS-IF97, Second Edition.附录：源程序(19个文件) 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 reg1.f90 module reg1 implicit none integer :: I(34),J(34) real*8 :: n(34) data I / 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 , 1 , 1 , 1 , 1 , 1 , 2 , 2 , 2 , 2 , 2 , 3 , 3 , 3 , 4 , 4 , 4 , 5 , 8 , 8 , 21 , 23 , 29 , 30 , 31 , 32 / data J / -2 , -1 , 0 , 1 , 2 , 3 , 4 , 5 , -9 , -7 , -1 , 0 , 1 , 3 , -3 , 0 , 1 , 3 , 17 , -4 , 0 , 6 , -5 , -2 , 10 , -8 , -11 , -6 , -29 , -31 , -38 , -39 , -40 , -41 / data n / 0.14632971213167 , -0.84548187169114 , -3.756360367204 , 3.3855169168385 , -0.95791963387872 , 0.15772038513228 , -0.016616417199501 , 8.1214629983