| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2995 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||
血未冷鐵桿木蟲(chóng) (正式寫(xiě)手)
|
[求助]
方程求解出錯(cuò),報(bào)錯(cuò)“未定義函數(shù)或變量 'hi'”求助! 已有1人參與
|
||
|
function htp=calhtpbyGW1(Pcr,Pevap,Molmass,fi,Gwf,uf,ufg,xdry,di,T,deltaT,densf,densnp,densg,conductf,conductnp,Cpnp,Cpf,rp,ifg) htpxmin=1.0e+1; htpxmax=3.0e+4; opt=optimset( 'fzero'); opt=optimset( opt,'TolX', 1.0e-3); htp=fzero(@(htpx)calerror(htpx,Pcr,Pevap,Molmass,fi,Gwf,uf,ufg,xdry,di,T,deltaT,densf,densnp,densg,conductf,conductnp,Cpnp,Cpf,rp,ifg),[htpxmin htpxmax],opt); end function error=calerror(htpx,Pcr,Pevap,Molmass,fi,Gwf,uf,ufg,xdry,di,T,deltaT,densf,densnp,densg,conductf,conductnp,Cpnp,Cpf,rp,ifg) unf=(1+2.5*fi+6.25*fi)*uf; %運(yùn)用Batchelor公式計(jì)算納米流體動(dòng)力黏度 fiw=1/((1+(1-fi))*(densf/densnp));% 體積分?jǐn)?shù)換算成質(zhì)量分?jǐn)?shù) kB=1.3806505*10^-23; %玻爾茲曼常數(shù),單位J/K densnf=fiw*densnp+(1-fiw)*densf; %計(jì)算納米流體密度 Cpnf=(fi*densnp*Cpnp+(1-fi)*densf*Cpf)/densnf;%計(jì)算納米流體的比熱 conductnf=(conductnp+2*conductf-2*fi*(conductf-conductnp))/(conductnp+2*conductf+fi*(conductf-conductnp))+((fi*(densnp*Cpnp)/2*conductf))*(kB*T/3*3.14*rp*unf)^0.5;%運(yùn)用宣益民模型計(jì)算導(dǎo)熱系數(shù) Prnf=unf*Cpnf/conductnf; %計(jì)算納米流體的普朗特?cái)?shù) Rel=(Gwf*(1-xdry)*di/unf); %計(jì)算液相Re Frl=Gwf^2/(densf^2*9.81*di); %計(jì)算Fr數(shù) Xtt=((1-xdry)/xdry)^0.9*(densg/densnf)^0.5*(unf/ufg)^0.1; %計(jì)算Xtt數(shù) F=1+24000*(hi*deltaT/Gwf*ifg)^1.16+1.37*(1/Xtt)^0.86; S=1/(1+1.15*10^(-6)*F^2*Rel^1.17); %%計(jì)算核態(tài)沸騰抑制因子,無(wú)量綱 hl=0.023*Rel^0.8*Prnf^0.4*conductnf/di; %%計(jì)算液相對(duì)流換熱系數(shù),W/(m2.K) Prevapwf=Pevap/Pcr; q=htpx*deltaT; hnb=55*Prevapwf^0.12*(-log10(Prevapwf))^(-0.55)*Molmass^(-0.5)*q^0.67; %% cooper公式計(jì)算核態(tài)沸騰換熱系數(shù),W/(m2.K)(COOPER公式單位必須用W/(m2.K) hi0=Frl^(0.1-2*Frl)*F*hl+Frl^0.5*S*hnb; %%計(jì)算對(duì)流沸騰換熱系數(shù)hi,W/(m2.K) hi1=hl*F+hnb*S; %%計(jì)算對(duì)流沸騰換熱系數(shù)hi,W/(m2.K) if Frl<=0.05 hi=hi0; else hi=hi1; end error=hi-htpx; end |

鐵桿木蟲(chóng) (正式寫(xiě)手)
|
那個(gè)地方打錯(cuò)了。。。不過(guò)還是沒(méi)能運(yùn)行出來(lái)。。初值htp=calhtpbyGW1(3661800,100000,152.93,0.0001,1600,0.00040588,0.000010826,0.3,0.00545,300.61,22.54,1457.6,0.00397,6.3918,0.075734,29,800,1026,10^(-8),2000)輸入后,出現(xiàn)了“錯(cuò)誤使用 fzero (line 274)區(qū)間端點(diǎn)處的函數(shù)值必須具有不同的符號(hào)!鼻蠼!! |


鐵桿木蟲(chóng) (正式寫(xiě)手)

鐵桿木蟲(chóng) (正式寫(xiě)手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +5 | 生物工程調(diào)劑 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 ,六級(jí)已過(guò),有科研經(jīng)歷 +14 | 曦熙兮 2026-03-15 | 14/700 |
|
|
[考研] 326求調(diào)劑 +4 | mlpqaz03 2026-03-15 | 4/200 |
|
|
[考研] 求調(diào)劑 +5 | 十三加油 2026-03-21 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過(guò) 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 一志愿中國(guó)石油大學(xué)(華東) 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 3/150 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|