| 5 | 1/1 | 返回列表 |
| 查看: 3000 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
血未冷鐵桿木蟲 (正式寫手)
|
[求助]
方程求解出錯(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)沸騰抑制因子,無量綱 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 |

鐵桿木蟲 (正式寫手)


鐵桿木蟲 (正式寫手)

鐵桿木蟲 (正式寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +10 | 丹青奶蓋 2026-03-26 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)308分調(diào)劑 +7 | 你好明天你好 2026-03-30 | 7/350 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +6 | @taotao 2026-03-26 | 7/350 |
|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +7 | 26要上岸 2026-03-28 | 7/350 |
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|