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

1、二维稳态导热实验报告传热学二维导热物体温度场数值模拟作者: 陈振兴学 号: 10037005学院(系): 化工学院专业: 过程装备与控制工程班级: 装备01指导老师: 李增耀试验时间: -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。四: 编程思绪及步骤图 结束否图3-1五、 程序及运行结果第三类边界条件: 1、 试验程序(C语言): / 1.cpp : 定义控制台应用程序入口点。/#include stdafx.h #includ

4、e#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); printf(ntttttt-陈振兴装备n); printf(n题目: 二维导热物体温度场电模拟试验n); printf(n矩形区域, l1=2.2;l2=3;l3=2;l4

5、=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;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

6、;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+) t12j=(t12j-1+t12j+1+2*t11j)/4; for(i=2;i6;i+) ti16=(ti-116+ti+116+2*ti

7、15)/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+) t6j=(dy*h2*t2/lmd+(t6j+1+t6j-1)/2+t5j)/(2+dy*h2/lmd); t11=(h1*dx*t1

8、+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*(t56+t65+t76/2+t67/2)/(3*lmd+h2*dx); for(i=1;i7;i+) for(j=1;j17;j+)

9、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;j7;j+) printf(%3.1f ,tij); l=l+1; if(l=6) printf(n); l=0; fo

10、r(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*(t77-10);q=(q1+q2)/2;e=abs(q2-q1)/q);printf(单位长度上1/4墙体

11、导热量为: %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 27.4 26.9 26.5 26.2 26.0 25.9 25.8 25.7 25.6 25.6

12、 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 26.9 24.8 22.5 19.9 16.7 15.2 14.6 14.3 14.2 14.

13、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 13.6 0 0 0 0 0 0 0 0 0 0;27.2 25.1 22.5 19.5 16.3

14、 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 #include#includeint _tmain(int argc, _TCHAR* argv

15、) 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, 假设区域内无内热源, 导热系数为常熟, 内外表面均为第一类边界条件且已知t1=30;t2=0;LMD=0.5

16、3;求该矩形区域内温度分布及垂直于纸面方向单位长度上经过墙体导热量。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+) aij=tij; for(i=5;i11;i+) for(j=1;j5;j+) ai

17、j=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+abs(tij-aij); dt=dt/(4*14); for(i=5;

18、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+) q1=q1+(30-t1j)*lmd; for(i=1;

19、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）: z=30.0 30.0 30.0

20、 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 12.1 12.0 12.0 12.0

21、 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.0 24.3 18.5 12.6 6

22、.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,yy=meshgrid(y,x);

23、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高级程序设计, 黄河水利出版社, 3、 杨世铭, 陶文铨编著, 传热学, 高等教育出版社,