| 查看: 3006 | 回复: 6 | ||
血未冷铁杆木虫 (正式写手)
|
[求助]
方程求解出错，报错“未定义函数或变量 '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; %运用Batchelor公式计算纳米流体动力黏度 fiw=1/((1+(1-fi))*(densf/densnp));% 体积分数换算成质量分数 kB=1.3806505*10^-23; %玻尔兹曼常数，单位J/K densnf=fiw*densnp+(1-fiw)*densf; %计算纳米流体密度 Cpnf=(fi*densnp*Cpnp+(1-fi)*densf*Cpf)/densnf;%计算纳米流体的比热 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;%运用宣益民模型计算导热系数 Prnf=unf*Cpnf/conductnf; %计算纳米流体的普朗特数 Rel=(Gwf*(1-xdry)*di/unf); %计算液相Re Frl=Gwf^2/(densf^2*9.81*di); %计算Fr数 Xtt=((1-xdry)/xdry)^0.9*(densg/densnf)^0.5*(unf/ufg)^0.1; %计算Xtt数 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); %%计算核态沸腾抑制因子，无量纲 hl=0.023*Rel^0.8*Prnf^0.4*conductnf/di; %%计算液相对流换热系数，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公式计算核态沸腾换热系数，W/(m2.K)（COOPER公式单位必须用W/(m2.K) hi0=Frl^(0.1-2*Frl)*F*hl+Frl^0.5*S*hnb; %%计算对流沸腾换热系数hi,W/(m2.K) hi1=hl*F+hnb*S; %%计算对流沸腾换热系数hi,W/(m2.K) if Frl<=0.05 hi=hi0; else hi=hi1; end error=hi-htpx; end |
回复此楼» 猜你喜欢
348分环境工程·调剂
已经有5人回复
-
考研调剂
已经有4人回复
-
283分材料与化工求调剂
已经有25人回复
-
298求调剂
已经有8人回复
-
266求调剂
已经有17人回复
-
265求调剂
已经有19人回复
-
考研调剂
已经有3人回复
-
材料295
已经有6人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有39人回复
-
285求调剂
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB 基础问题求教
已经有6人回复
2楼2015-05-22 11:14:04
血未冷
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 7512.9
- 散金: 76
- 红花: 1
- 帖子: 326
- 在线: 38.7小时
- 虫号: 2637318
- 注册: 2013-09-05
- 性别: GG
- 专业: 工程热物理与能源利用
3楼2015-05-22 11:25:33
血未冷
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 7512.9
- 散金: 76
- 红花: 1
- 帖子: 326
- 在线: 38.7小时
- 虫号: 2637318
- 注册: 2013-09-05
- 性别: GG
- 专业: 工程热物理与能源利用
4楼2015-05-22 11:28:40
血未冷
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 7512.9
- 散金: 76
- 红花: 1
- 帖子: 326
- 在线: 38.7小时
- 虫号: 2637318
- 注册: 2013-09-05
- 性别: GG
- 专业: 工程热物理与能源利用
5楼2015-05-22 11:30:02
6楼2015-05-22 11:45:45
 |
7楼2015-05-22 21:40:51
