传热学课程数值计算实验报告.docxVIP

最新精品文档，知识共享! 最新精品文档，知识共享! 传热学课程数值计算实验报告 一 问题重述 有一个用砖砌成的长方形截面的冷气通道， 其截面尺寸如下图所示， 假设在垂直于纸面 方向上冷空气及砖墙的温度变化很小，可以近似忽略。试计算： （1） 砖墙横截面上的温度分布； （2） 垂直于纸面方向上的每米长度上通过砖墙的导热量。 墙壁内、外表面的流体的温度分别为 10 °C、30 °C ； 内外表面均为第三类边界条件，且已知： tj=30°c, hh =10.33 W/m2 °c t：2=10°C, 0=3.93 W/m2 °C 砖墙的导热系数lamda = 0.53 W / mL°c 离散 离散 首先考虑到整个墙体截面的对称性，对称地取其四分之一截面进行研究，可以适当 简化问题。  为方便与温度场电模拟实验的数据进行对比，对温度场的离散与之保持一致，即网 格为0.1m的正方形，且对内、外墙表面的流体温度的模拟各设一排节点进行模拟。在 实际的编程过程中，为使程序简洁且更有条理，建立一个网格矩阵，矩阵中每个元素的 值根据温度场中对应位置的节点的边界情况或计算特点进行设定从而将温度场节点进 行分类。网格矩阵如下。下面会结合该矩阵详细阐述迭代方程的建立过程。 GOOOOOODOOOOOD0OO 00 0 0 0 0 0 0 0 0 G 0 0 0 ；5 5 3 *555555333^ ! 1 1 1 1 11111111114 2 1111 1111111111^ ￡1111 I 1 1 1 1 1 1 1 1 1 J s 1 1 1 i I 1 1 1 1 1 1 1 1 1 J E 1 1 1 1 1 t~5~5~5~55~5~5S~5~a ￡1111^ 000000000( E 1 1 1 1 4 DDOODODOO( ￡1111^ 0000000601 ￡ 1 1 1 i ￡1111^ aOGGDDOBD( I 0 0—0 D Q—Q EJ 售 0 1 三方程的建立与求解 在介绍迭代方程之前，需要说明的是： a、为提高迭代速度，采用 Gauss-Seidel迭代方 法，即总是将最新得出来的节点数据用到迭代过程中去； b、每次迭代的初始温度场记为 tfi， 正在迭代的新温度场记为 tft ； c、墙壁的内节点初值设为 20。 网格矩阵将温度场的节点按照边界情况或计算条件进行了分类、标记。每个不同数值含 义及对应的迭代方程列举如下： 0：温度场中流体的节点，迭代过程中值保持不变： tft(i,j)=tfi(i,j); :墙壁的内节点，控制方程可有热平衡法或泰勒级数展开得出： tft(i,j)=0.25*(tfi(i,j+1)+tfi(i+1,j)+tft(i-1,j)+tft(i,j-1)); :外墙面边界点(角点除外)，由热平衡法整理可有： tft(i,j)=(h1*d*tft(i,j-1)+0.5*lamda*(tft(i-1,j)+tfi(i+1,j)+2*tfi(i,j+1)))/(2*lamda+h1*d); 3:外墙边界点(角点除外)，与2相似，可由热平衡法整理得出： tft(i,j)=(h1*d*tft(i-1,j)+0.5*lamda*(tft(i,j-1)+tfi(i,j+1)+2*tfi(i+1,j)))/(2*lamda+h1*d); 4：墙壁对称面，绝热边界条件： tft(2:7,17)=tft(2:7,16); 5:内边界点(角点除外)，可由热平衡法整理得出： tft(i,j)=(h2*d*tfi(i+1,j)+0.5*lamda*(tft(i,j-1)+tfi(i,j+1)+2*tft(i-1,j)))/(2*lamda+h2*d); 6:内墙边界点(角点除外)，与5相似，由热平衡法整理得出： tft(i,j)=(h2*d*tfi(i,j+1)+0.5*lamda*(tft(i-1,j)+tfi(i+1,j)+2*tft(i,j-1)))/(2*lamda+h2*d); 7:墙壁对称面，绝热边界条件： tft(13,2:7)=tft(12,2:7); 内墙壁角点，由于其形状位置的特殊性，方程需要单独考虑，仍然由热平衡法整理得出： tft(7,7)=(tfi(8,8)*h2*d+0.5*lamda*(2*tft(6,7)+2*tft(7,6)+tfi(7,8)+tfi(8,7)))/(h2*d+3*lamda); 9:外墙壁角点，由于其形状位置的特殊性，方程需要单独考虑，由热平衡法整理可得： tft(2,2)=(h1*d*tft(1,1)+0.5*lamda*(tfi(2,3)+tfi(3,2)))/(lamda+h1*d); 值得指出的是，在实际的迭代程序编写中， 对于绝热边界条件即 4和7点的计算是要