传热学大作业报告

第一题假设1流体表面传热系数为常数2一维稳态导热3肋片物性为常数4肋片顶端绝热工况1分为8*4网格，则节点数为9*5m=2~9,n=1时θm=2~9,n=2~5θm=2~9.n=5θm=9,n=2~4θm=9,n=5θm=9,n=1θ编程如下Bi=50*(0.02/4)/100;a=ones(9,5);a(1,:)=80;b=a;c=a;C=1;whileC>1.0E-6form=2:8b(m,1)=(a(m-1,1)+a(m+1,1)+2*a(m,2))/4;endform=2:8forn=2:4b(m,n)=(a(m+1,n)+a(m-1,n)+a(m,n+1)+a(m,n-1))/4;endendform=2:8b(m,5)=(a(m-1,5)+a(m+1,5)+2*a(m,4))/(4+2*Bi);endforn=2:4b(9,n)=(a(9,n+1)+a(9,n-1)+2*a(8,n))/(4+2*Bi);endb(9,5)=(a(9,4)+a(8,5))/(2+2*Bi);b(9,1)=(a(8,1)+a(9,2))/(2+1*Bi);fori=1:9forj=1:5c(i,j)=abs((a(i,j)-b(i,j))/max(max(a)));endendC=max(max(c));a=b;endcontour(b',50)N=(0.5*(b(1,5)+b(9,1))+sum(b(:,5))-b(1,5)+sum(b(9,:))-b(9,1)-b(9,5))/(12*80)工况2类似于工况1把网格的节点数改为17*5Bi=400*(0.02/4)/8;a=ones(17,5);a(1,:)=80;b=a;c=a;C=1;whileC>1.0E-6form=2:16b(m,1)=(a(m-1,1)+a(m+1,1)+2*a(m,2))/4;endform=2:16forn=2:4b(m,n)=(a(m+1,n)+a(m-1,n)+a(m,n+1)+a(m,n-1))/4;endendform=2:16b(m,5)=(a(m-1,5)+a(m+1,5)+2*a(m,4))/(4+2*Bi);endforn=2:4b(17,n)=(a(17,n+1)+a(17,n-1)+2*a(16,n))/(4+2*Bi);endb(17,5)=(a(17,4)+a(16,5))/(2+2*Bi);b(17,1)=(a(16,1)+a(17,2))/(2+1*Bi);fori=1:17forj=1:5c(i,j)=abs((a(i,j)-b(i,j))/max(max(a)));endendC=max(max(c));a=b;endcontour(b',50)N=(0.5*(b(1,5)+b(17,1))+sum(b(:,5))-b(1,5)+sum(b(17,:))-b(17,1)-b(17,5))/(20*80)第二题此题不需迭代是一维非稳态平板导热的显示格式a=zeros(10001,11);%10000时间层10等分a(1,:)=80;%t0F0=1.39E-5*0.18/(0.1/10)^2;Bi=1163*(0.1/10)/50;tf=300;fori=1:10000forn=2:10a(i+1,n)=F0*(a(i,n-1)+a(i,n+1))+(1-2*F0)*a(i,n);a(i+1,11)=(1-2*F0)*a(i,11)+2*F0*a(i,10);a(i+1,1)=a(i,1)*(1-2*F0*Bi-2*F0)+2*F0*a(i,2)+2*F0*Bi*tf;endendx=1:11;plot(x,a(201,:),'o-',x,a(1001,:),'o-',x,a(2001,:),'o-',x,a(3001,:),x,a(4001,:),'o-',x,a(5001,:),'o-',x,a(6001,:),'o-',x,a(7001,:),'o-',x,a(8001,:),'o-',x,a(10001