哈尔滨工程大学传热学大作业数值计算matlab程序内容

1传热学作业数值计算2数值计算matlab程序内容：tw1=10;%赋初值tw2=20;c=1.5;p2=20;p1=c*p2;L2=40;L1=c*L2;deltaX=L2/p2;a=p2+1;b=p1+1;ti=ones(a,b)*5;m1=ones(a,b);m1(a,2:b-1)=zeros(1,b-2);m1(2:a,1)=zeros(a-1,1);m1(2:a,b)=zeros(a-1,1);m1(1,:)=ones(1,b)*2;k=0;max1=1.0;tn=ti;while(max11e-6)max1=0;k=k+1;fori=1:1:aforj=1:1:b3m=m1(i,j);n=ti(i,j);switchmcase0tn(i,j)=tw1;case1tn(i,j)=0.25*(tn(i,j+1)+tn(i,j-1)+tn(i+1,j)+tn(i-1,j));case2tn(i,j)=tw1+tw2*sin(pi*(j-1)/(b-1));ender=abs(tn(i,j)-n);ifermax1max1=er;endendendti=tn;endktimax1t2=ones(a,b);%求解析温度场fori=a:-1:1forj=1:1:by=deltaX*(a-i);x=deltaX*(j-1);t2(i,j)=tw1+tw2*sin(pi*x/L1)*(sinh(pi*y/L1))/(sinh(pi*L2/L1));endendt2迭代次数k=706数值解温度场ti4数值解每次迭代的最大误差max1=9.8531e-07解析温度场t2取第11行的解析解和数值解的点曲线为第11行的解析解的直线，散点为其数值解的点第11行的误差=[数值解(11行)–解析解(11行)]/解析解5数值温度场图像解析温度场图像数值解与解析解的误差数值计算matlab程序内容：tw1=10;tw2=20;c=1.5;p2=20;p1=c*p2;L2=20;deltaX=L2/p2;L1=c*L2;a=p2+1;b=p1+1;ti=ones(a,b)*5;6m1=ones(a,b);m1(a,2:b-1)=zeros(1,b-2);m1(2:a,1)=zeros(a-1,1);m1(2:a,b)=zeros(a-1,1);m1(1,:)=ones(1,b)*2;k=0;max1=1.0;tn=ti;while(max11e-6)max1=0;k=k+1;fori=1:1:aforj=1:1:bm=m1(i,j);n=ti(i,j);switchmcase0tn(i,j)=tw1;case1tn(i,j)=0.25*(tn(i,j+1)+tn(i,j-1)+tn(i+1,j)+tn(i-1,j));case2tn(i,j)=tw2;ender=abs(tn(i,j)-n);ifermax1max1=er;endendendti=tn;endktimax1tx=ones(a,b);fori=1:1:aforj=1:1:by=(a-i)*deltaX;x=(j-1)*deltaX;m=sym('m');g=(((-1)^(m+1)+1)/m)*sin(m*pi*x/L1)*sinh(m*pi*y/L1)/sinh(m*pi*L2/L1);h=symsum(g,m,1,100);tx(i,j)=2*h*(tw2-tw1)/pi+tw1;end7endtx迭代次数k=695数值解温度场ti数值解每次迭代的最大误差max1=9.8243e-07解析温度场tx=取第11行的解析解和数值解的点曲线为第11行的解析解的直线，散点为其数值解的点第11行的误差=[数值解(11行)–解析解(11行)]/解析解8数值温度场ti图像：解析温度场tx图像：数值解与解析解的误差9数值计算matlab程序内容：t0=90;tf=10;L=10;c=0.25;p2=20;p1=p2/c;B=c*L;d=0.5*B;h=10;a=p2+1;b=p1+1;deltaX=B/p2;lambda=160;Bi=h*deltaX/lambda;ti=ones(a,b)*10;m1=ones(a,b)*3;m1(2:a-1,1)=zeros(a-2,1);m1(a,2:b-1)=ones(1,b-2);m1(1,2:b-1)=ones(1,b-2)*6;m1(2:a-1,b)=ones(a-2,1)*2;m1(1,b)=ones(1,1)*4;m1(a,b)=ones(1,1)*5;m1(1,1)=7;m1(a,1)=8;tn=ti;max1=1.0;k=0;while(max11e-6)k=k+1;max1=0;fori=1:1:a10forj=1:1:bm=m1(i,j);n=tn(i,j);switchmcase0tn(i,j)=t0;case1tn(i,j)=(2*tn(i-1,j)+tn(i,j-1)+tn(i,j+1)-4*tf)/(4+2*Bi)+tf;case2tn(i,j)=(2*tn(i,j-1)+tn(i-1,j)+tn(i+1,j)-4*tf)/(4+2*Bi)+tf;case3tn(i,j)=0.25*(tn(i,j-1)+tn(i,j+1)+tn(i-1,j)+tn(i+1,j));case4tn(i,j)=(tn(i,j-1)+tn(i+1,j)-2*tf)/(2*Bi+2)+tf;case5tn(i,j)=(tn(i,j-1)+tn(i-1,j)-2*tf)/(2*Bi+2)+tf;case6tn(i,j)=(2*tn(i+1,j)+tn(i,j-1)+tn(i,j+1)-4*tf)/(4+2*Bi)+tf;case7tn(i,j)=t0;case8tn(i,j)=t0;ender=abs(tn(i,j)-n);ifermax1max1=er;endendendti=tn;endktita=ones(a,b);Bi1=h*d/lambda;sbi=sqrt(Bi1);fori=1:1:aforj=1:1:bifi(a+1)/2y=-(i-(a+1)/2)*deltaX;elsey=((a+1)/2-i)*deltaX;endx=deltaX*(j-1);ta(i,j)=(cosh(sbi*(L-x)/d)+sbi*sinh(sbi*(L-x)/d))*(t0-tf)/(cosh(sbi*L/d)+sbi*sinh(sbi*L/d))+tf;11endendta迭代次数k=1461数值解温度场ti解析温度场ta取第11行的解析解和数值解的点曲线为第11行的解析解的直线，散点为其数值解的点第11行的误差=[数值解(11行)–解析解(11行)]/解析解图像如下12数值温度场图像解析温度场图像数值解与解析解的误差13数值计算matlab程序内容：tw=10;L2=15;c=0.75;L1=L2/c;p2=24;p1=p2/c;deltaX=2*L2/p2;a=p2+1;b=p1+1;lambda=16;qv0=24;ti=ones(a,b)*5;m1=ones(a,b);m1(1,:)=zeros(1,b);m1(2:a,b)=zeros(a-1,1);m1(2:a,1)=zeros(a-1,1);m1(a,2:b-1)=zeros(1,b-2);tn=ti;max1=1.0;k=0;while(max11e-6)max1=0;k=k+1;fori=1:1:aforj=1:1:bm=m1(i,j);n=tn(i,j);switchmcase0tn(i,j)=tw;case1tn(i,j)=0.25*(tn(i-1,j)+tn(i+1,j)+tn(i,j-1)+tn(i,j+1)+qv0*(deltaX^2)/lambda);14ender=abs(tn(i,j)-n);ifermax1max1=er;endendendti=tn;endtik;tx=ones(a,b);fori=1:1:aforj=1:1:bifi(a+1)/2y=-(i-(a+1)/2)*deltaX;elsey=((a+1)/2-i)*deltaX;endifj(b+1)/2x=(j-(b+1)/2)*deltaX;elsex=-((b+1)/2-j)*deltaX;endm=sym('m');xi=(2*m-1)*pi/2;g=((-1)^m)/(xi^3)*(cosh(xi*y/L1)/cosh(xi*L2/L1))*cos(xi*x/L1);h=symsum(g,m,1,100);tx(i,j)=2*qv0*L1^2/lambda*h+qv0*(L1^2-x^2)/(2*lambda)+tw;endendtx数值温度场ti解析温度场tx15取第13行的解析解和数值解的点曲线为第13行的解析解的直线，散点为其数值解的点第13行的误差=[数值解(13行)–解析解(13行)]/解析解16数值温度场图像解析温度场图像数值解与解析解的误差