| 5 | 1/1 | 返回列表 |
| 查看: 2999 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +10 | 顧九笙要謙虛 2026-03-24 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 考研調(diào)劑 +6 | 來好運(yùn)來來來 2026-03-21 | 7/350 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
[考研] 321求調(diào)劑 +3 | 璞玉~~ 2026-03-25 | 3/150 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
[考研] 271求調(diào)劑 +4 | 生如夏花… 2026-03-22 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-24 | 3/150 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠心求調(diào)劑 +15 | 吃吃吃才有意義 2026-03-19 | 16/800 |
|
|
[考研]
|
pk3725069 2026-03-19 | 17/850 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +7 | holy J 2026-03-21 | 7/350 |
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
[考研] 材料292調(diào)劑 +8 | 橘頌思美人 2026-03-23 | 8/400 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|