二维稳态导热实验报告

《二维稳态导热实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、传热学二维导热物体温度场的数值模拟作者:陈振兴学号:10037005化工学院学院（系）:专业:过程装备与控制工程装备01班级:指导教师：李增耀二维导热物体温度场的数值模拟一、物理描述有一个用砖砌成的长方形截面的冷空气通道，其截面尺寸和示意图如图1-1所示，假设在垂直纸面方向上冷空气及砖墙的温度变化很小，可以近似地予以忽略。在以下情况下试计算：(1)砖墙横截面上的温度分布；(2)垂直于纸面方向的每米长度上通过砖墙的导热量。1、内外表面均为第三类边界条件，iL己知：=30oC，尽=10.33W/tn2-°Ct2=l(rq/z2

2、=3.93W/m2-°C砖墙的导热系数A=0.53W/m•°C2、内外壁分布均匀地维持在0°C及30°C;3,]nm-i二、数学描述该结构的导热问题可以作为二维问题处理，并且其截面如图1-1所示，由于对称性，仅研究其］/4部分即可。其网络节点划分如图2-1;上述问题为二维矩形域内的稳态、无内热源、常物性的导热问题，对于这样的物理问题，我们知道，描写其的微分方程即控制方程，就是异热微分方程：第三类边界条件：内外表面均为第三类边界条件，且已知:C,=30°C,h}=10.33W/m2-°C=10°C,/?2=3.93W/m2^

3、C碎墙的导热系数A=0.53W化•a图2-1三：方程的离散如上图2-1所示，用一系列与坐标轴平行的网络线把求解区域划分成许多子区域，以网格线的交点作为需要确定温度值的空间位置，即节点，节点的位置己该点在两个方向上的标号m、n来表示。每一个节点都可以看成是以它为中心的小区域的代表，如上（m，n）:对于（m,n）为内节点时：由级数展开法或热平衡法都可以得到，当Ax=Ay时.•W-l，"+，'卜1）对于（m，n）为边界节点吋：位于平直边界上的节点：t":."=-（tw+.,+2t"卜i，"+1'，'.'卜1）外部角点：如图2-1

4、中a、b、d、e、f点，内部角点：如图2-1屮c点，t〃'，w-+2tw_

5、，M+2tw,w+i+tw,n_j）由已知条件有，当m=l或n=13时的节点的温度衡为G=30°C,当（m=6且n<9）和（n=8H6

6、#include#includcinttmain(intargc,_TCHAR*argv口){inti,j,1;floatdt=l.0,dx=0.1,dy=0.1;floatt[13][17],a[13][17];floatql=0,q2=0,q=0,e;floatlmd=0.53,hl=10.33,h2=3.93,tl=30,t2=10,ep=l.0e~7;A打印出题目Vprintf(〃ttt二维稳态导热问题tt〃)；printf(〃tttttt陈振兴装备〃

7、)；printf("题目：二维导热物体温度场的电模拟实验");printf("矩形区域，11=2.2;12=3;13=2;14=1.2,假设区域内无内热源，导热系数为常熟，内外表面均为第三类边界条件且已知11=30;t2=10;hl=10.33;h2=3.93;LMD=0.53;求该矩形区域内的温度分布及垂直于纸面方向的单位长度上通过墙体的导热量。");A各节点上的温度值吟{for(j=0;j<17;j++){t[0][j]:30.0;}for(i=l;i<13;i++){t[i][0]=30.0;}for

8、(i=7;i<13;i++){t[i][7]=10.0;}for(j=8;j<17;j++){t[7][j]:10.0;}for(i=l;i<7;i++)for(j=l;j<17;j++)t[i][j]=20,a[i][j]=0;for(i=7;i<13;i++)for(j=l;j<7;j++)t[i][j]=20，a[i][j]=0;}while(dt>=ep){{for(i=l;i<7;i++)for(j=l;j<17;j++)for(i=7;i〈13;i++)for(j:l;j<7;j++)for(i=6;i<12;

9、i++)for(j=2;j<6;j++)t[i]=[j]+t[i+l][j]+t[i][j-l]+t[i][j+l])/4;}{for(i=2;i<6;i++)for(j:2;j<16;j++)t[i][j]=(t[i-l][j]+t[i+l][j]+t[i][j-l]+t[i][j+l])/4;}{f