24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

popsino

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 8
帖子: 5
在线: 8.2小时
虫号: 2656537
注册: 2013-09-15
专业: 工程热物理相关交叉领域

[求助] FORTRAN调用NIST的REFPROP计算物性参数，程序只能成功调用一次，不知道哪里出了问题

大家好，目前做有机工质的超音速喷嘴设计，物性参数计算通过调用NIST的REFPROP实现,初始状态给定总温，总压，需要编一个子程序，基于等熵流动假设，给出一个速度V，程序可以计算静压、静温等其他参数， h=hs-0.5*velocity**2  ，s=ss，调用了 HSFLSH 语句，
call HSFLSH (h,s,z,t,p,D,Dl,Dv,x,y,q,e,cv,cp,w,ierr,herr)，调用时第一次可以正确计算，但是第二次调用就计算结果有问题，参数传递也检查过几遍了，一直不知道是什么原因？若是手动查表的话，是可以查出结果的。请大神帮忙指点，可能是那部分的问题？谢谢了，新来的，只有这么多金币了。

X Y 速度分量U 速度分量V 马赫数合速度流动角静压密度静温
0.106 0 137.9 0 0.9998 137.9 0 46.1823 0.0193 293.2
0.1058 0.0526 137.9 0 0.3097 137.9 0 0 0 3094.2
0.1049 0.1053 138 0 0.3099 138 0 0 0 3094.2
0.1034 0.1579 138.1 0 0.3102 138.1 0 0 0 3094.2
0.1013 0.2105 138.3 0 0.3106 138.3 0 0 0 3094.2
0.0987 0.2632 138.5 0 0.3111 138.5 0 0 0 3094.2
0.0955 0.3158 138.8 0 0.3118 138.8 0 0 0 3094.2

手动查表结果：

　 X Y 合速度静焓(kJ/kg) 熵(kJ/kg-K) 静压(MPa) 密度(kg/m?) 静温(K) 音速(m/s)
1 0.106 0 137.9 422.22 1.8197 0.046041 2.5837 293.15 137.88
2 0.1058 0.0526 137.9 422.22 1.8197 0.046041 2.5837 293.15 137.88
3 0.1049 0.1053 138 422.12 1.8197 0.045775 2.5697 293.02 137.86
4 0.1034 0.1579 138.1 422.02 1.8197 0.045521 2.5563 292.89 137.84
5 0.1013 0.2105 138.3 421.81 1.8197 0.044988 2.5282 292.63 137.8
6 0.0987 0.2632 138.5 421.60 1.8197 0.04446 2.5004 292.37 137.76
7 0.0955 0.3158 138.8 421.29 1.8197 0.043691 2.4599 291.98 137.7

SUBROUTINE THERMO2 (velocity,hs,ss,pressure,density,temperature,sonicspeed,mach)
   implicit double precision (a-h,o-z)
   implicit integer (i-k,m,n)
   parameter (ncmax=20) !max number of components in mixture
   dimension x(ncmax),xliq(ncmax),xvap(ncmax),f(ncmax)
   character hrf*3, herr*255
   character*255 hf(ncmax),hfmix
   REAL,INTENT(IN):: velocity
   REAL,INTENT(IN):: hs
   REAL,INTENT(IN):: ss
   REAL,INTENT(OUT):

ressure
REAL,INTENT(OUT):

ensity
   REAL,INTENT(OUT)::temperature
   REAL,INTENT(OUT)::sonicspeed
   REAL,INTENT(OUT)::mach
黄色底纹为初始化固定格式
   i=1
   hf(1)=\'R245FA.fld\'
   hfmix=\'hmx.bnc\'
   hrf=\'DEF\'
   call SETUP (i,hf,hfmix,hrf,ierr,herr)
   if (ierr.ne.0) write (*,*) herr
以下部分为主要计算过程，只使用了一个函数HSFLSH ，
   h=hs-0.5*velocity**2
   s=ss
   call HSFLSH (h,s,z,t,p,D,Dl,Dv,x,y,q,e,cv,cp,w,ierr,herr)
   temperature  = t
   sonicspeed = w
   mach       = velocity/sonicspeed
   pressure    = p
   Density    = D
   RETURN
   END

回复此楼