二维稳态导热实验报告.doc

《二维稳态导热实验报告.doc》由会员分享，可在线阅读，更多相关《二维稳态导热实验报告.doc（19页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、,传热学二维导热物体温度场的数值模拟作者：陈振兴学 号：10037005学院(系)：化工学院专业：过程装备与控制工程班级：装备01指导教师：李增耀实验时间：2012-10二维导热物体温度场的数值模拟 一、物理描述有一个用砖砌成的长方形截面的冷空气通道，其截面尺寸和示意图如图1-1所示，假设在垂直纸面方向上冷空气及砖墙的温度变化很小，可以近似地予以忽略。在以下情况下试计算：（1）砖墙横截面上的温度分布；（2）垂直于纸面方向的每米长度上通过砖墙的导热量。1、内外表面均为第三类边界条件，且已知：.33 砖墙的导热系数2、内外壁分布均匀地维持在0及30； 图1-1 二、数学描述该结构的导热问题可以作为

2、二维问题处理，并且其截面如图1-1所示，由于对称性，仅研究其1/4部分即可。其网络节点划分如图2-1；上述问题为二维矩形域内的稳态、无内热源、常物性的导热问题，对于这样的物理问题，我们知道，描写其的微分方程即控制方程，就是导热微分方程： 第三类边界条件：内外表面均为第三类边界条件，且已知： 砖墙的导热系数 a f（m，n）c b = n e m d 图2-1三：方程的离散如上图2-1所示，用一系列与坐标轴平行的网络线把求解区域划分成许多子区域，以网格线的交点作为需要确定温度值的空间位置，即节点，节点的位置已该点在两个方向上的标号m、n来表示。每一个节点都可以看成是以它为中心的小区域的代表，如上

3、（m，n）：对于（m，n）为内节点时：由级数展开法或热平衡法都可以得到，当=时： 对于（m，n）为边界节点时：位于平直边界上的节点：外部角点：如图2-1中a、b、d、e、f点，内部角点：如图2-1中c点，由已知条件有，当m=1或n=13时的节点的温度衡为=30，当（m=6且n9）和（n=8且6m17）时的节点的温度为=10。四：编程思路及流程图 开始 输入已知参数 结束否说明边界条件取定初始试探值TA（i, j）=0计算新的内节点和边界点温度T(i, j)T=TA比较所有节点|TA(i, j)- T(i, j)|EPSILO ？计算内外边界散热量及热平衡偏差输出内外边界散热量及热平衡偏差结束图

4、3-1五、程序及运行结果第三类边界条件：1、实验程序(C语言)：/ 1.cpp : 定义控制台应用程序的入口点。/#include stdafx.h #include#includeint _tmain(int argc, _TCHAR* argv) int i,j,l; float dt=1.0,dx=0.1,dy=0.1; float t1317,a1317; float q1=0,q2=0,q=0,e; float lmd=0.53,h1=10.33,h2=3.93,t1=30,t2=10,ep=1.0e-7; /*打印出题目*/ printf(ttt二维稳态导热问题tt); print

5、f(ntttttt-陈振兴装备n); printf(n题目：二维导热物体温度场的电模拟实验n); printf(n矩形区域，l1=2.2;l2=3;l3=2;l4=1.2，假设区域内无内热源，导热系数为常熟，内外表面均为第三类边界条件且已知t1=30;t2=10;h1=10.33;h2=3.93;LMD=0.53;求该矩形区域内的温度分布及垂直于纸面方向的单位长度上通过墙体的导热量。n); /*各节点上的温度值*/ for(j=0;j17;j+) t0j=30.0; for(i=1;i13;i+) ti0=30.0; for(i=7;i13;i+) ti7=10.0; for(j=8;j17;

6、j+) t7j=10.0; for(i=1;i7;i+)for(j=1;j17;j+)tij=20,aij=0; for(i=7;i13;i+)for(j=1;j=ep) for(i=1;i7;i+)for(j=1;j17;j+)aij=tij; for(i=7;i13;i+)for(j=1;j7;j+)aij=tij; for(i=6;i12;i+)for(j=2;j6;j+)tij=(ti-1j+ti+1j+tij-1+tij+1)/4; for(i=2;i6;i+) for(j=2;j16;j+) tij=(ti-1j+ti+1j+tij-1+tij+1)/4; for(j=2;j6;j

7、+)t12j=(t12j-1+t12j+1+2*t11j)/4; for(i=2;i6;i+)ti16=(ti-116+ti+116+2*ti15)/4; for(i=2;i12;i+)ti1=(dx*h1*t1/lmd+(ti+11+ti-11)/2+ti2)/(2+dx*h1/lmd); for(j=2;j16;j+)t1j=(dy*h1*t1/lmd+(t1j+1+t1j-1)/2+t2j)/(2+dy*h1/lmd); for(i=7;i12;i+)ti6=(dx*h2*t2/lmd+(ti+16+ti-16)/2+ti5)/(2+dx*h2/lmd); for(j=7;j16;j+)

8、t6j=(dy*h2*t2/lmd+(t6j+1+t6j-1)/2+t5j)/(2+dy*h2/lmd); t11=(h1*dx*t1+lmd*(t21+t12)/2)/(lmd+h1*dy);t116=(h1*dx*t1+lmd*(t115+t216)/2)/(lmd+h1*dy);t616=(h2*dx*t2+lmd*(t615+t516)/2)/(lmd+h2*dy);t121=(h1*dx*t1+lmd*(t122+t111)/2)/(lmd+h1*dy);t126=(h2*dx*t2+lmd*(t125+t116)/2)/(lmd+h2*dy);t66=(h2*dy*t2+lmd*(

9、t56+t65+t76/2+t67/2)/(3*lmd+h2*dx); for(i=1;i7;i+)for(j=1;j17;j+)dt=dt+abs(tij-aij); dt=dt/(6*16);for(i=7;i13;i+)for(j=1;j7;j+)dt=dt+abs(tij-aij); dt=dt/(6*6); printf(温度分布为：ttttttttt);l=0; for(i=1;i7;i+) for(j=1;j17;j+) printf(%3.1f ,tij); l=l+1; if(l=16) printf(n); l=0; l=0; for(i=7;i13;i+) for(j=1

10、;j7;j+)printf(%3.1f ,tij); l=l+1; if(l=6) printf(n); l=0; for(j=2;j16;j+)q1=q1+(30-t1j)*h1*dx; for(i=2;i12;i+)q1=q1+(30-ti1)*h1*dy;for(j=7;j17;j+)q2=q2+(t6j-10)*h2*dx;for(i=7;i12;i+)q2=q2+(ti6-10)*h2*dy;q1=q1+h1*(dx/2*(30-t116)+dy/2*(30-t121)+dx*(30-t11); q2=q2+h2*(dx/2*(t616-10)+dy/2*(t126-10)+dx*(

11、t77-10);q=(q1+q2)/2;e=abs(q2-q1)/q);printf(单位长度上1/4墙体的导热量为：%4.2fW,偏差为：%3.2f,q,e);getchar();getchar();return 0;运行结果图：图3-2实验算得导热量为97.62W，与数值模拟的偏差为(26.73*4W-97.62W)/（26.73*4）W*100%=8.7%2、数值模拟程序（matlab）：z=29.9 29.7 29.5 29.3 29.1 29.0 28.8 28.7 28.7 28.6 28.6 28.5 28.5 28.5 28.5 28.229.7 29.1 28.5 28.0

12、27.4 26.9 26.5 26.2 26.0 25.9 25.8 25.7 25.6 25.6 25.6 25.529.5 28.5 27.6 26.6 25.7 24.8 24.1 23.6 23.3 23.1 22.9 22.8 22.8 22.7 22.7 22.729.3 28.0 26.6 25.2 23.8 22.5 21.6 20.9 20.5 20.2 20.0 19.9 19.9 19.8 19.8 19.829.1 27.4 25.7 23.8 21.9 19.9 18.7 17.9 17.5 17.2 17.1 17.0 16.9 16.9 16.9 16.929.0

13、 26.9 24.8 22.5 19.9 16.7 15.2 14.6 14.3 14.2 14.1 14.0 14.0 14.0 14.0 13.928.8 26.5 24.1 21.5 18.6 15.2 0 0 0 0 0 0 0 0 0 0;28.7 26.2 23.5 20.8 17.8 14.6 0 0 0 0 0 0 0 0 0 0;28.6 25.9 23.2 20.3 17.3 14.2 0 0 0 0 0 0 0 0 0 0;28.5 25.7 22.8 20.0 17.0 14.0 0 0 0 0 0 0 0 0 0 0;28.3 25.4 22.6 19.7 16.7

14、13.6 0 0 0 0 0 0 0 0 0 0;27.2 25.1 22.5 19.5 16.3 12.5 0 0 0 0 0 0 0 0 0 0;v=18 22 26;xx,yy=meshgrid(y,x);surf(xx,yy,z);colorbar;xlabel(x);ylabel(y);zlabel(z);az = 0;el =- 90;view(az, el);shading interp;axis tight;figure,contour(xx,yy,z,v);grid on数值模拟图：图3-3图3-4第一类边界条件1、 实验程序(C语言)：#include stdafx.h #

15、include#includeint _tmain(int argc, _TCHAR* argv) int i,j,l; float dt=1.0,dx=0.1,dy=0.1; float t1216,a1216; float q1=0,q2=0,q=0,e; float lmd=0.53,t1=30,t2=0,ep=1.0e-7; /*打印出题目*/ printf(ttt二维稳态导热问题tt); printf(ntttttt-陈振兴装备n); printf(n题目：二维导热物体温度场的电模拟实验n); printf(n矩形区域，l1=2.2;l2=3;l3=2;l4=1.2，假设区域内无内热

16、源，导热系数为常熟，内外表面均为第一类边界条件且已知t1=30;t2=0;LMD=0.53;求该矩形区域内的温度分布及垂直于纸面方向的单位长度上通过墙体的导热量。n); /*各节点上的温度值*/ for(j=0;j16;j+) t0j=30.0; for(i=1;i12;i+) ti0=30.0; for(i=5;i12;i+) ti5=0.0; for(j=6;j16;j+) t5j=0.0; for(i=1;i5;i+)for(j=1;j15;j+)tij=20,aij=0; for(i=5;i11;i+)for(j=1;j=ep) for(i=1;i5;i+)for(j=1;j15;j+

17、)aij=tij; for(i=5;i11;i+)for(j=1;j5;j+)aij=tij; for(i=5;i11;i+)for(j=1;j5;j+)tij=(ti-1j+ti+1j+tij-1+tij+1)/4; for(i=1;i5;i+)for(j=1;j15;j+) tij=(ti-1j+ti+1j+tij-1+tij+1)/4; for(j=1;j5;j+)t11j=(t11j-1+t11j+1+2*t10j)/4; for(i=1;i5;i+)ti15=(ti-115+ti+115+2*ti14)/4; for(i=1;i5;i+)for(j=1;j15;j+)dt=dt+ab

18、s(tij-aij); dt=dt/(4*14);for(i=5;i11;i+)for(j=1;j5;j+)dt=dt+abs(tij-aij); dt=dt/(4*6); printf(温度分布为：ttttttttt);l=0; for(i=0;i6;i+) for(j=0;j16;j+) printf(%4.1f ,tij); l=l+1; if(l=16) printf(n); l=0; l=0; for(i=6;i12;i+) for(j=0;j6;j+)printf(%4.1f ,tij); l=l+1; if(l=6) printf(n); l=0; for(j=1;j15;j+)

19、q1=q1+(30-t1j)*lmd; for(i=1;i11;i+)q1=q1+(30-ti1)*lmd;for(j=6;j15;j+)q2=q2+t4j*lmd;for(i=6;i11;i+)q2=q2+ti4*lmd;q1=q1+lmd*(30-t115)/2+(30-t111)/2); q2=q2+lmd*(t415/2+t114/2+t44);q=(q1+q2)/2;e=abs(q2-q1)/q);printf(单位长度上墙体的导热量为：%4.2fW,偏差为：%3.2f,4*q,e);getchar();getchar();return 0;运行结果图：2、数值模拟程序（matlab

20、）：z=30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0;30.0 29.0 28.1 27.1 26.2 25.5 24.9 24.5 24.3 24.2 24.1 24.1 24.0 24.0 24.0 24.0;30.0 28.1 26.1 24.2 22.3 20.7 19.6 18.9 18.5 18.3 18.2 18.1 18.0 18.0 18.0 18.0;30.0 27.1 24.2 21.2 18.1 15.5 13.9 13.0 12.5 12.3 12.2 1

21、2.1 12.0 12.0 12.0 12.0;30.0 26.2 22.3 18.1 13.6 9.1 7.4 6.7 6.4 6.2 6.1 6.1 6.0 6.0 6.0 6.0;30.0 25.5 20.7 15.5 9.1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;30.0 24.9 19.6 13.9 7.4 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;30.0 24.5 18.9 13.0 6.7 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;30.

22、0 24.3 18.5 12.6 6.4 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;30.0 24.2 18.3 12.3 6.2 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;30.0 24.1 18.2 12.2 6.1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;30.0 24.1 18.2 12.2 6.1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;v=12 18 24;x=1:1:12;y=1:1:16;xx

23、,yy=meshgrid(y,x);surf(xx,yy,z);colorbar;xlabel(x);ylabel(y);zlabel(z);az = 0;el =- 90;view(az, el);shading interp;axis tight;figure,contour(xx,yy,z,v);grid on六、实验感想这次花了很长时间来完成传热学数值模拟实验，总的来说还是收获很大的。首先从实验本身来看，结合书本上关于数值求解的讲解，以及在网上搜索的一些例题，我对数值计算中起决定性作用的划分网格、求各节点的相应温度值的计算方法有了比较清晰的认识。数值求解在解决实际问题中确实能够发挥很大作用，能最接近真实情况的将实际问题的解用数学模型表示出来，从而加深对这些问题的理解，同时也能加强对二维稳态导热问题等相关知识的掌握。其次从具体实施的角度来看，我在编改程序的过程中也再次详细认识和掌握了很多C语言以及matlab的相关知识，对迭代求解的思路有了比较清晰地认识 。而且也认识到编写程序所需要的严谨与认真是又快又好解决问题的关键。否则浪费了时间精力还做了无用功，得不偿失。七、参考文献1、西安交通大学等编，传热学实验指导书，热与流体实验中心2、周振红等主编，Fortran90/95高级程序设计，黄河水利出版社，20053、杨世铭，陶文铨编著，传热学，高等教育出版社，2006