| 查看: 1230 | 回复: 2 | ||
892226227铜虫 (初入文坛)
|
[求助]
求助修改模拟列管换热器的程序。感激不尽
|
|
求助：帮我修改一下这个模拟列管换热器的程序。感激不尽 %已知 wh=68250/3600; wc=175000/3600; Th1=382; Tc1=275; A=130;%传热面积 Di=0.0225;%管内径 Do=0.025;%管外径 D=0.7;%壳体内径 L=0.6;%管长 AT=0.0438;%管程通道截面积 AS=0.0525;%壳程通道截面积 NB=19;%壳程挡板数 B=0.3;%挡板间距 Ri=0.0005;%管程污垢热阻 Ro=0.0001;%壳程污垢热阻 Ki=12.5;%渣油品特性参数 Ko=12.5;%原油品特性参数 D1=0.919;%渣相对密度 D0=0.850;%原相对密度 Tc20=285;%冷流体的估算出口温度 %计算管中心距 t=Di/(1-AS/(B*D)); %计算壳体当量直径 De=4*(t^2-pi*Do^2/4)/(pi*Do); while 1 %最外层循环求冷热物流出口温度 %计算管内热物流性质 Thm=(Th1+Th20)/2; cpi=((0.7072+0.000551*D1)*Thm-0.318D1)*(0.055*Ki+0.35)*4.18; %计算管外冷物流性质 Tc2=fzero(@wd,295,[],Tc1,Th1,Th20,Thm,D0,cpi,Ko,wh,wc);%求只含Tc2的一元函数的零点 %--------函数定义------- function s =wd(Tc2,Tc1,Th1,Th20,Thm,D0,cpi,Ko,wh,wc); s=(wh*cpi)*(Th1-Th20)-(wc*cpo)*(Tc2-Tc1); Tcm=(Tc1+Tc2)/2; cpo=((0.7072+0.000551*D0)*Tcm-0.318D0)*(0.055*Ko+0.35)*4.18; %计算传热量 Q=wh*(Th1-Th20); %计算管内外流体的相对密度 roui=0.942+0.248*x+0.174*D1^2+0.0841/(x*D1)-0.312*x/D1-0.556*exp(-x); xi=1+Thm/100; rouo=0.942+0.248*x+0.174*D0^2+0.0841/(x*Do0)-0.312*x/D0-0.556*exp(-x); xo=1+Tcm/100; %管内体积流速 ui=wh/(AT*roui); %管外体积流速 uo=wc/(AS*rouo); i=1; fprintf('no=%d\tuo\tui\n',i,uo.ui); %计算管内外流体热导率 lumdi=0.4213*(1-0.00054*Thm)/D1; lumdo=0.4213*(1-0.00054*Tcm)/D0; %计算管内外流体运动粘度 a=log(log(1500+1.22))-blog(50+273); b=log((log(1500+1.22))/log(120+1.22))/log((50+273)/(100+273)); miui=exp(exp(a+b*log(Thm+273)))-1.22; a=log(log(90+1.22))-blog(50+273); b=log((log(90+1.22))/log(13+1.22))/log((50+273)/(100+273)); miuo=exp(exp(a+b*log(Tcm+273)))-1.22; %粘度 ati=miui*roui; ato=miuo*rouo; %计算管内外流体特征数 Reo=Do*Wc/(ato*AS); Rei=Di*Wh/(ati*AT); Pro=cpo*ato/lumdo; Pri=cpi*ati/lumdi; %估算内外壁温 Hi1=0.023*Rei^0.8*Pri^0.33*lumdi/Di; %初值为不受运动粘度影响的表面储热系数 while 1 %迭代求内壁温 Twi=Tmc+Q/(Hi1*A*225/250); a=log(log(1500+1.22))-blog(50+273); b=log((log(1500+1.22))/log(120+1.22))/log((50+273)/(100+273)); miuiw=exp(exp(a+b*log(Twi+273)))-1.22; Hi=Hi1(miui/miuiw)^0.14; %受温度影响的表面传热系数 Tiw=Tmh+Q/(Hi*A*225/250); %替代 if abs(Tiw-Twi)<=0.5 %收敛 break; else Twi=Tiw; end %if语句结束 end %内壁温循环结束 %求管外的表面传热系数的初值 if Reo>2000 Ho1=0.36*Reo^0.55*Pro^0.33*lumdo/Do; else Ho1=0.527*Reo^0.5*Pro^0.33*lumdo/Do; end %if语句结束 while 1 %迭代求外壁温 Two=Tmc+Q/(Ho1*A); a=log(log(1500+1.22))-blog(50+273); b=log((log(1500+1.22))/log(120+1.22))/log((50+273)/(100+273)); miuow=exp(exp(a+b*log(Two+273)))-1.22; Ho=Ho1(miuo/miuow)^0.14; Tow=Tmc+Q/(Ho*A); if abs(Tow-Two)<=0.5 %收敛 break; else Two=Tow;%替代 end %if语句结束 end %内壁温循环结束 %表面传热系数 K=1/((1/Hi+Ri)*250/225+1/Ho+Ro); %计算传热效率 ibcn=2/((1+wh*cpi/wc*cpo)+sqrt(1+wh*cpi/wc*cpo)*(1+exp(-F))/(1-exp(-F))); F=ntu*sqrt(1+(wh*cpi/wc*cpo)^2); ntu=K*A/(wh*cpi); %核算Th2 Th2=Th1-ibcn(Th1-Tc1); Tc2=Tc1+(wh*cpi/wc*cpo)*(Th1-Th2); if abs(Tc20-Tc2)<=0.1 %收敛 break; else T20=Th2;%替代 end %if语句结束 fprintf('no=%d\tuii=%f\tReii=%f\tPrii=%f\tHii=%f\tKii=%f\n',i,ui,Rei,Pri,Hi,Ki); fprintf('no=%d\tuoi=%f\tReoi=%f\tProi=%f\tHoi=%f\tKoi=%f\n',i,uo,Reo,Pro,Ho,Ko); end % 最外层求冷热物流出口温度循环结束 fprintf('计算阻力…'); %计算换热器阻力 %管程阻力 if Rei<100000 fi=0.4513*Rei^(-0.2663); else fi=0.2864*Rei^(-0.2258); end gi=ui*roui; dpt=(fi*L/Di+4)*gi^2*2*1.5/(2000*D1); %壳程阻力 if Reo<150 fo=120*Reo^(-0.993); elseif Reo>1500 fo=0.7664*Reo^(-0.0854); else fo=10*(15.312/(log(Reo))^4.735-0.44); end go=uo*rouo; dps=(D*(NB+1)*fo*0.15/De+10)*go^2/(2000*D0); |
回复此楼» 猜你喜欢
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
面上合作单位盖章
已经有8人回复
-
面上函评意见出来了，像什么等级？
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助啊！求一程序，用matlab程序做，用蒙特卡洛方法模拟
已经有10人回复
- 【求助】对程序设计语言非常喜欢，但是我不会，想参加个培训机构，给点建议，感激不尽
已经有13人回复
- 【求助】翅片式列管换热器换热面积问题
已经有5人回复
- 【求助】关于GAMBIT中质量不好的网格手动修改
已经有5人回复
- 【求助】盘管换热器的K值
已经有4人回复
- 【求助】合成氨过程模拟 结果错误，怎么办？并问下如何修改收敛次数！【已解决】
已经有11人回复
- 【求助】物料属性可以用UDF修改吗
已经有8人回复
- 【求助】换热器计算
已经有13人回复
2楼2011-05-23 23:17:50
只为找工作
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 417.2
- 散金: 20
- 帖子: 77
- 在线: 40.1小时
- 虫号: 3518244
- 注册: 2014-11-04
- 性别: GG
- 专业: 内流流体力学
3楼2015-12-15 22:10:54
投票:
30