| 6 | 1/1 | 返回列表 |
| 查看: 3078 | 回復(fù): 5 | ||
popsino新蟲 (初入文壇)
|
[求助]
FORTRAN調(diào)用NIST的REFPROP計(jì)算物性參數(shù),程序只能成功調(diào)用一次,不知道哪里出了問題
|
|
大家好,目前做有機(jī)工質(zhì)的超音速噴嘴設(shè)計(jì),物性參數(shù)計(jì)算通過(guò)調(diào)用NIST的REFPROP實(shí)現(xiàn),初始狀態(tài)給定總溫,總壓,需要編一個(gè)子程序,基于等熵流動(dòng)假設(shè),給出一個(gè)速度V,程序可以計(jì)算靜壓、靜溫等其他參數(shù), h=hs-0.5*velocity**2 ,s=ss,調(diào)用了 HSFLSH 語(yǔ)句, call HSFLSH (h,s,z,t,p,D,Dl,Dv,x,y,q,e,cv,cp,w,ierr,herr),調(diào)用時(shí)第一次可以正確計(jì)算,但是第二次調(diào)用就計(jì)算結(jié)果有問題,參數(shù)傳遞也檢查過(guò)幾遍了,一直不知道是什么原因?若是手動(dòng)查表的話,是可以查出結(jié)果的。請(qǐng)大神幫忙指點(diǎn),可能是那部分的問題?謝謝了,新來(lái)的,只有這么多金幣了。 X Y 速度分量U 速度分量V 馬赫數(shù) 合速度 流動(dòng)角 靜壓 密度 靜溫 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 手動(dòng)查表結(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 以下部分為主要計(jì)算過(guò)程,只使用了一個(gè)函數(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 |
新蟲 (初入文壇)
新蟲 (初入文壇)
新蟲 (初入文壇)
新蟲 (初入文壇)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 +8 | 妮妮ninicgb 2026-03-15 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|