二维导热物体温度场地数值模拟

《二维导热物体温度场地数值模拟》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文案二维导热物体温度场的数值模拟作者：学号：学院(系)：能源与动力工程学院专业：能源动力系统及自动化班级：二维导热物体温度场的数值模拟精彩文档实用标准文案一：物理问题有一个用砖砌成的长方形截面的冷空气通道，其截面尺寸和示意图如图1-1所示，假设在垂直纸面方向上冷空气及砖墙的温度变化很小，可以近似地予以忽略。在下列两种情况下试计算：（1）砖墙横截面上的温度分布；（2）垂直于纸面方向的每米长度上通过砖墙的导热量。第一种情况：内外壁分布均匀地维持在0及30；第二种情况：内外表面均为第三类边界条件，且已知：砖墙的

2、导热系数二：数学描述该结构的导热问题可以作为二维问题处理，并且其截面如图1-1所示，由于对称性，仅研究其1/4部分即可。精彩文档实用标准文案其网络节点划分如图fac（m，n）b=nemd上述问题为二维矩形域内的稳态、无内热源、常物性的导热问题，对于这样的物理问题，我们知道，描写其的微分方程即控制方程，就是导热微分方程：第一类边界条件：内外壁分布均匀地维持在0及30；=30=0第三类边界条件：内外表面均为第三类边界条件，且已知：砖墙的导热系数精彩文档实用标准文案三：方程的离散如上图所示，用一系列与坐标轴平行的网络线

3、把求解区域划分成许多子区域，以网格线的交点作为需要确定温度值的空间位置，即节点，节点的位置已该点在两个方向上的标号m、n来表示。每一个节点都可以看成是以它为中心的小区域的代表，如上（m，n）：对于（m，n）为内节点时：由热平衡法可以得到，当=时：①　对于（m，n）为边界节点时：l恒温边界只需特殊考虑位于绝热平直边界上的节点：l对流边界分为角点、绝热边界点和对流边界点。1.绝热边界点：2.对流边界点：3.外角点：4.内角点：四：编程思路及流程图开始输入已知参数精彩文档实用标准文案说明边界条件取定初始试探值TA（i,

4、j）=0计算新的内节点和边界点温度T(i,j)T=TA比较所有节点

5、TA(i,j)-T(i,j)

6、

7、edai1(t);fori=2:5forj=2:16whilet(i,j)-t0(i,j)>=0.00001;t0=t;t=diedai1(t);endendendfori=6:12forj=2:5whilet(i,j)-t0(i,j)>=0.00001;t0=t;t=diedai1(t);endendendt0=tfori=2:11Q1x=Q1x+0.53*(t(i,1)-t(i,2));endQ1x=Q1x+0.53*(t(12,1)-t(12,2))/2;forj=2:15Q1y=Q1y+0.53*(t(1

8、,j)-t(2,j));endQ1y=Q1y+0.53*(t(1,16)-t(2,16))/2;Q1=(Q1x+Q1y)*4fori=6:11Q2x=Q2x+0.53*(t(i,5)-t(i,6));endQ2x=Q2x+0.53*(t(12,5)-t(12,6))/2;forj=6:15Q2y=Q2y+0.53*(t(5,j)-t(6,j));精彩文档实用标准文案endQ2y=Q2y+0.53*(t(5,16)-t(6,16))/2;Q2=(Q2x+Q2y)*4n=2*abs(Q1-Q2)/(Q1+Q2)t0=

9、tf;tf=diedai2(tf);fori=2:5forj=2:16whiletf(i,j)-t0(i,j)>=0.00001;t0=tf;tf=diedai2(tf);endendendfori=6:12forj=2:5whiletf(i,j)-t0(i,j)>=0.00001;t0=tf;tf=diedai2(tf);endendendt0=tfQ1x=0;Q1y=0;Q1=0;Q2x=0;Q2y=0;Q2=0;n=0;fori=1:11Q1x=Q1x+10*0.1*(30-tf(i,1));endQ1x=

10、Q1x+10*0.05*(30-tf(12,1));forj=2:15Q1y=Q1y+10*0.1*(30-tf(1,j));endQ1y=Q1y+10*0.05*(30-tf(1,16));fori=6:11Q2x=Q2x+4*0.1*(tf(i,6)-10);endQ2x=Q2x+4*0.05*(tf(12,6)-10);forj=7:15Q2y=Q2y+4*0.1*