二维稳态导热实验报告

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且n<9）和（n=8且6<m<17）时的节点的温度为=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<stdio.h>#include<math.h>int _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

5、,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=1.2，假设区域内无内热源，导热系数为常熟，内外表面均为第三类边界条件且已知t1=30;t2=10;h1=10.33;h2=3.93;LMD=0.53;求该矩形区域内的温度分布及垂直于纸面方向的单位长度上通过墙体的导热量。n");

6、 /*各节点上的温度值*/ for(j=0;j<17;j+) t0j=30.0; for(i=1;i<13;i+) ti0=30.0; for(i=7;i<13;i+) ti7=10.0; for(j=8;j<17;j+) t7j=10.0; for(i=1;i<7;i+)for(j=1;j<17;j+)tij=20,aij=0; for(i=7;i<13;i+)for(j=1;j<7;j+)tij=20,aij=0; while(dt>=ep) for(i=1;i<7;i+)for(j=1;j<17;j+)aij=tij; f

7、or(i=7;i<13;i+)for(j=1;j<7;j+)aij=tij; for(i=6;i<12;i+)for(j=2;j<6;j+)tij=(ti-1j+ti+1j+tij-1+tij+1)/4; for(i=2;i<6;i+) for(j=2;j<16;j+) tij=(ti-1j+ti+1j+tij-1+tij+1)/4; for(j=2;j<6;j+)t12j=(t12j-1+t12j+1+2*t11j)/4; for(i=2;i<6;i+)ti16=(ti-116+ti+116+2*ti15)/4; for(i=2;i<12

8、;i+)ti1=(dx*h1*t1/lmd+(ti+11+ti-11)/2+ti2)/(2+dx*h1/lmd); for(j=2;j<16;j+)t1j=(dy*h1*t1/lmd+(t1j+1+t1j-1)/2+t2j)/(2+dy*h1/lmd); for(i=7;i<12;i+)ti6=(dx*h2*t2/lmd+(ti+16+ti-16)/2+ti5)/(2+dx*h2/lmd); for(j=7;j<16;j+)t6j=(dy*h2*t2/lmd+(t6j+1+t6j-1)/2+t5j)/(2+dy*h2/lmd); t11=(h1*dx*t1+lmd*(t21+

9、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;i<7;i+)for(j=1;j<17;j+)dt=d

10、t+abs(tij-aij); dt=dt/(6*16);for(i=7;i<13;i+)for(j=1;j<7;j+)dt=dt+abs(tij-aij); dt=dt/(6*6); printf("温度分布为：ttttttttt");l=0; for(i=1;i<7;i+) for(j=1;j<17;j+) printf("%3.1f ",tij); l=l+1; if(l=16) printf("n"); l=0; l=0; for(i=7;i<13;i+) for(j=1;j<7;j+)pr

11、intf("%3.1f ",tij); l=l+1; if(l=6) printf("n"); l=0; for(j=2;j<16;j+)q1=q1+(30-t1j)*h1*dx; for(i=2;i<12;i+)q1=q1+(30-ti1)*h1*dy;for(j=7;j<17;j+)q2=q2+(t6j-10)*h2*dx;for(i=7;i<12;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/

12、2*(t616-10)+dy/2*(t126-10)+dx*(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 2

13、8.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 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

14、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.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

15、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 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,contou

16、r(xx,yy,z,v);grid on数值模拟图：图3-3图3-4第一类边界条件1、 实验程序(C语言)：#include "stdafx.h" #include<stdio.h>#include<math.h>int _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(&

17、quot;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.53;求该矩形区域内的温度分布及垂直于纸面方向的单位长度上通过墙体的导热量。n"); /*各节点上的温度值*/ for(j=0;j<16;j+) t0j=30.0; for(i=1;i&

18、lt;12;i+) ti0=30.0; for(i=5;i<12;i+) ti5=0.0; for(j=6;j<16;j+) t5j=0.0; for(i=1;i<5;i+)for(j=1;j<15;j+)tij=20,aij=0; for(i=5;i<11;i+)for(j=1;j<5;j+)tij=20,aij=0; while(dt>=ep) for(i=1;i<5;i+)for(j=1;j<15;j+)aij=tij; for(i=5;i<11;i+)for(j=1;j<5;j+)aij=tij; for(i=5;i&l

19、t;11;i+)for(j=1;j<5;j+)tij=(ti-1j+ti+1j+tij-1+tij+1)/4; for(i=1;i<5;i+)for(j=1;j<15;j+) tij=(ti-1j+ti+1j+tij-1+tij+1)/4; for(j=1;j<5;j+)t11j=(t11j-1+t11j+1+2*t10j)/4; for(i=1;i<5;i+)ti15=(ti-115+ti+115+2*ti14)/4; for(i=1;i<5;i+)for(j=1;j<15;j+)dt=dt+abs(tij-aij); dt=dt/(4*14);fo

20、r(i=5;i<11;i+)for(j=1;j<5;j+)dt=dt+abs(tij-aij); dt=dt/(4*6); printf("温度分布为：ttttttttt");l=0; for(i=0;i<6;i+) for(j=0;j<16;j+) printf("%4.1f ",tij); l=l+1; if(l=16) printf("n"); l=0; l=0; for(i=6;i<12;i+) for(j=0;j<6;j+)printf("%4.1f ",tij); l

21、=l+1; if(l=6) printf("n"); l=0; for(j=1;j<15;j+)q1=q1+(30-t1j)*lmd; for(i=1;i<11;i+)q1=q1+(30-ti1)*lmd;for(j=6;j<15;j+)q2=q2+t4j*lmd;for(i=6;i<11;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("单位长度上墙体的

22、导热量为：%4.2fW,偏差为：%3.2f",4*q,e);getchar();getchar();return 0;运行结果图：2、数值模拟程序（matlab）：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 1

23、8.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 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.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;