北科大传热传质学二维墙角传热大作业VIP免费

传热大作业二维墙角稳态导热问题的数值模拟目录1.问题描述..............................................................................................32.数学建模..............................................................................................33.MATLAB程序源代码............................................................................64.证明网格无关性................................................................................105.数值模拟结果....................................................................................116.温度场的数据文件（另附excel表格在压缩包中）........................157.实验测量数据绘图与模拟等温线对比.............................................158.实验感悟............................................................................................171.问题描述等温边界条件：t130C,t20C对流换热边界条件：t130C,h19.30W/(m2C)t210C,h13.93W/(m2C)2.数学建模外部角点内部角点边界点内部点外部角点ktm-1,ntm,n2ktm,n12h(△x△y)(tftm,n)022h=h1△x=△y=d内部角点ktm-1,ntm,n△x△yktm,n1tm,n△y△xktm,n1tm,n2△y△xktm1,ntm,n2△x△yh(△x△y)(tftm,n)022h=h2边界角点ktm-1,ntm,n△x△yktm,n1tm,n2△y△xktm,n1tm,n2△y△xh△y(tftm,n)0内边界tf10h=h2外边界tf30h=h1内部节点ktm-1,ntm,n△x△yktm,n1tm,n△y△xktm,n1tm,n△y△xktm1,ntm,n△x△y0内部节点区域划分方式3.MATLAB程序源代码当对流换热系数h无限大时（h>10^6），对流换热边界条件可近似认为定壁温边界条件，这样就能实现一个程序两种计算的实现。除此之外，我们组的程序可通过改变A,B,C,D四个数的数值来实现不同墙壁长度的代码复用。详细程序如下：源程序：clcL1=2.2;L2=3.0;L3=2.0;L4=1.2;k=0.53;%设置初始条件begin=1;while(begin==1)c=input('输入1，计算定壁温边界条件；输入2，计算对流换热边界条件');switchccase1T1=30;T2=0;h1=10^6;h2=10^6;%对流换热系数很大时可近似认为是定壁温条件case2T1=30;T2=10.0;h1=9.30;h2=3.93;endd=input('输入节点间距离d=');%定义相邻节点的距离M=L1/d;N=L2/d;p=N*M;%计算节点个数以及节点总数A=0.5;%A=(L1-L4)/2;可计算出A,B,C,D的数值，然后更改数值来实现程序的复用。B=0.5;%B=(L2-L3)/2;C=1.7;%C=(L1+L4)/2;D=2.5;%D=(L2+L3)/2;disp('所取节点个数为');disp(p);T=zeros(M+1,N+1);%定义数组Tm=zeros(M+1,N+1);%定义数组n=1;e=1;o=input('请输入收敛精度要求');while(e>o)m=0;while(m<=n)forj=1:M+1fori=1:N+1Tm(j,i)=T(j,i);%满足收敛条件时将迭代值赋给对应的结点endendT(1,1)=(2*h1*d*T1/k+T(1,2)+T(2,1))/(2*h1*d/k+2);%外左上结点T(M+1,1)=(2*h1*d*T1/k+T(M,1)+T(M+1,2))/(2*h1*d/k+2);%外左下结点T(M+1,N+1)=(2*h1*d*T1/k+T(M,N+1)+T(M+1,N))/(2*h1*d/k+2);%外右下结点T(1,N+1)=(2*h1*d*T1/k+T(1,N)+T(2,N+1))/(2*h1*d/k+2);%外右上结点%内边界的四个角节点T(1+A/d,1+B/d)=(2*T(1+A/d,B/d)+2*T(A/d,1+B/d)+T(2+A/d,1+B/d)+T(1+A/d,2+B/d)+2*h2*d*T2/k)/(6+2*h2*d/k);%内左上节点T(1+C/d,1+B/d)=(2*T(1+C/d,B/d)+2*T(C/d,1+B/d)+T(2+C/d,1+B/d)+T(1+C/d,2+B/d)+2*h2*d*T2/k)/(6+2*h2*d/k);%内左下节点T(1+C/d,1+D/d)=(2*T(1+C/d,D/d)+2*T(C/d,1+D/d)+T(2+C/d,1+D/d)+T(1+C/d,2+D/d)+2*h2*d*T2/k)/(6+2*h2*d/k);%内右下节点T(1+A/d,1+D/d)=(2*T(1+A/d,D/d)+2*T(A/d,1+D/d)+T(2+A/d,1+D/d)+T(1+A/d,2+D/d)+2*h2*d*T2/k)/(6+2*h2*d/k);%内右上节点%外边界的节点（非角点）fori=2:NT(1,i)=(0.5*T(1,i-1)+0.5*T(1,i+1)+T(2,i)+h1*d*T1/k)/(h1*d/k+2);%外上边界T(M+1,i)=(0.5*T(M+1,i-1)+0.5*T(M+1,i+1)+T(M,i)+h1*d*T1/k)/(h1*d/k+2);%外下边界endforj=2:MT(j,1)=(0.5*T(j-1,1)+0.5*T(j+1,1)+T(j,2)+h1*d*T1/k)/(h1*d/k+2);%外左边界T(j,N+1)=(0.5...