| 查看: 3101 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
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):  ressureREAL,INTENT(OUT):  ensityREAL,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 |
回复此楼» 猜你喜欢
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有168人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有2人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有9人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何用C++6.0调用refprop物性查询软件
已经有9人回复
- 求Nist软件：Refprop9
已经有23人回复
woaiwojia456
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 7.9
- 红花: 1
- 帖子: 42
- 在线: 27.6小时
- 虫号: 3621572
- 注册: 2014-12-31
- 性别: MM
- 专业: 热力学
4楼2016-04-22 10:44:54
2楼2015-07-21 16:40:20
woaiwojia456
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 7.9
- 红花: 1
- 帖子: 42
- 在线: 27.6小时
- 虫号: 3621572
- 注册: 2014-12-31
- 性别: MM
- 专业: 热力学
3楼2016-04-22 10:43:56
woaiwojia456
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 7.9
- 红花: 1
- 帖子: 42
- 在线: 27.6小时
- 虫号: 3621572
- 注册: 2014-12-31
- 性别: MM
- 专业: 热力学
5楼2016-05-25 11:26:46
