等温边界与对流边界换热完整matlab代码.pdf

《等温边界与对流边界换热完整matlab代码.pdf》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、西安交大能动专业传热学上机实验等温边界：clear;clc;T=zeros(12,16);temp=zeros(12,16);%赋初值forj=1:16T(1,j)=30;endfori=1:12T(i,1)=30;endforj=6:16fori=6:12T(i,j)=0;endendtemp=T;t=max(max(T));%单次迭代最高温度eps=ones(12,16);k=1;%迭代次数%热平衡方程组while(max(max(eps))>1.0e-6)%“正常”点fori=2:5forj=2:15T(i,j)=(T

2、(i-1,j)+T(i+1,j)+T(i,j+1)+T(i,j-1))/4;endendfori=6:11forj=2:5T(i,j)=(T(i-1,j)+T(i+1,j)+T(i,j+1)+T(i,j-1))/4;endend%热对称边界点fori=2:5T(i,16)=T(i-1,16)/4+T(i+1,16)/4+T(i,15)/2;endforj=2:5T(12,j)=T(11,j)/2+T(12,j-1)/4+T(12,j+1)/4;endt=max(max(T));%单次迭代最高温度eps=abs(T-temp)

3、/t;%精度temp=T;k=k+1;%迭代次数endA=flipud(T);%将矩阵上下颠倒pcolor(A);shadinginterp%figure(1),contour(A,1000);%等温线集合figure,contour(A,[12,18,24])%作出三条等温线%计算内外壁导热量，偏差dx=0.1;dy=0.1;lamda=0.53;temp1=zeros(13,1);temp2=zeros(13,1);temp3=zeros(13,1);temp4=zeros(13,1);forj=3:15temp1(j-

4、2)=(T(1,j)-T(2,j))*lamda*dx*1/dy;endfori=3:11temp2(i-2)=(T(i,1)-T(i,2))*lamda*dy*1/dx;endheat_out=sum(temp1)+sum(temp2)+(T(2,1)-T(2,2))*lamda*dy*1/dx+(T(1,2)-T(2,2))*lamda*dx*1/dy+(T(1,16)-T(2,16))*lamda*dx/2*1/dy+(T(12,1)-T(12,2))*lamda*dy/2*1/dx;forj=6:15temp3(j-

5、5)=(T(5,j)-T(6,j))*lamda*dx*1/dy;endfori=6:11temp4(i-5)=(T(i,5)-T(i,6))*lamda*dy*1/dx;endheat_in=sum(temp3)+sum(temp4)+(T(5,16)-T(6,16))*lamda*dx/2*1/dy+(T(12,5)-T(12,6))*lamda*dy/2*1/dx;heat=(heat_in+heat_out)/2;e=abs(heat_in-heat_out)/heat;fprintf('外壁传热量=%2.3f'

6、,heat_out)fprintf('内壁传热量=%2.3f',heat_in)fprintf('平均传热量=%2.3f',heat)fprintf('相对偏差=%d',e)fprintf('墙角总传热量=%f',4*heat)对流边界：clear;clc;T=zeros(13,17);temp=zeros(13,17);dx=0.1;dy=dx;lamda=0.53;h_out=10.6;h_in=3.975;tf_out=30;tf_in=10;%参数%赋初值forj=1:17T(1,j)=30;endf

7、ori=1:13T(i,1)=30;endforj=8:17fori=8:13T(i,j)=10;endendtemp=T;t=max(max(T));%温度最高值eps=ones(13,17);%精度k=1;%热平衡方程组while(max(max(eps))>1.0e-6)%“正常”点fori=3:6forj=3:16T(i,j)=(T(i-1,j)+T(i+1,j)+T(i,j+1)+T(i,j-1))/4;endendfori=7:12forj=3:6T(i,j)=(T(i-1,j)+T(i+1,j)+T(i,j+1

8、)+T(i,j-1))/4;endend%绝热边界fori=3:6T(i,17)=T(i-1,17)/4+T(i+1,17)/4+T(i,16)/2;endforj=3:6T(13,j)=T(12,j)/2+T(13,j-1)/4+T(13,j+1)/4;end%外平直边界forj=3:1