 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3094 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
popsino新蟲 (初入文壇)
|
[求助]
FORTRAN調(diào)用NIST的REFPROP計算物性參數(shù)，程序只能成功調(diào)用一次，不知道哪里出了問題
|
|
|
大家好，目前做有機(jī)工質(zhì)的超音速噴嘴設(shè)計，物性參數(shù)計算通過調(diào)用NIST的REFPROP實現(xiàn),初始狀態(tài)給定總溫，總壓，需要編一個子程序，基于等熵流動假設(shè)，給出一個速度V，程序可以計算靜壓、靜溫等其他參數(shù)， h=hs-0.5*velocity**2 ，s=ss，調(diào)用了 HSFLSH 語句， call HSFLSH (h,s,z,t,p,D,Dl,Dv,x,y,q,e,cv,cp,w,ierr,herr)，調(diào)用時第一次可以正確計算，但是第二次調(diào)用就計算結(jié)果有問題，參數(shù)傳遞也檢查過幾遍了，一直不知道是什么原因？若是手動查表的話，是可以查出結(jié)果的。請大神幫忙指點，可能是那部分的問題？謝謝了，新來的，只有這么多金幣了。 X Y 速度分量U 速度分量V 馬赫數(shù) 合速度 流動角 靜壓 密度 靜溫 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 手動查表結(jié)果： 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 以下部分為主要計算過程，只使用了一個函數(shù)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 |
回復(fù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有154人回復(fù)
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 如何用C++6.0調(diào)用refprop物性查詢軟件
已經(jīng)有9人回復(fù)
- 求Nist軟件：Refprop9
已經(jīng)有23人回復(fù)
woaiwojia456
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 7.9
- 紅花: 1
- 帖子: 42
- 在線: 27.6小時
- 蟲號: 3621572
- 注冊: 2014-12-31
- 性別: MM
- 專業(yè): 熱力學(xué)
3樓2016-04-22 10:43:56
2樓2015-07-21 16:40:20
woaiwojia456
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 7.9
- 紅花: 1
- 帖子: 42
- 在線: 27.6小時
- 蟲號: 3621572
- 注冊: 2014-12-31
- 性別: MM
- 專業(yè): 熱力學(xué)
4樓2016-04-22 10:44:54
woaiwojia456
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 7.9
- 紅花: 1
- 帖子: 42
- 在線: 27.6小時
- 蟲號: 3621572
- 注冊: 2014-12-31
- 性別: MM
- 專業(yè): 熱力學(xué)
5樓2016-05-25 11:26:46
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
 10+4
|
chibby 2026-03-25 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 322求調(diào)劑 +4 | 舊吢 2026-03-24 | 4/200 |
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +7 | 曼111 2026-03-24 | 8/400 |
|
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
|
[考研] 316求調(diào)劑 +5 | 江辭666 2026-03-26 | 5/250 |
|
|
|
[考研] 329求調(diào)劑 +7 | 鈕恩雪 2026-03-25 | 7/350 |
|
|
|
[考研] 333求調(diào)劑 +7 | 87639 2026-03-21 | 12/600 |
|
|
|
[考研] 求調(diào)劑 +5 | 蘆lty 2026-03-25 | 6/300 |
|
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
|
[考研] 材料277求調(diào)劑 +5 | min3 2026-03-24 | 5/250 |
|
|
|
[考研] 0856求調(diào)劑 +8 | zhn03 2026-03-25 | 9/450 |
|
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap－ 2026-03-22 | 5/250 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
|
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|