傳热学课程数值计算实验报告.docVIP

传热学课程数值计算实验报告 一 问题重述 有一个用砖砌成的长方形截面的冷气通道，其截面尺寸如下图所示，假设在垂直于纸面方向上冷空气及砖墙的温度变化很小，可以近似忽略。试计算： 砖墙横截面上的温度分布； 垂直于纸面方向上的每米长度上通过砖墙的导热量。 墙壁内、外表面的流体的温度分别为10、30； 内外表面均为第三类边界条件，且已知： =30，=10.33 =10, =3.93 砖墙的导热系数 二 离散 首先考虑到整个墙体截面的对称性，对称地取其四分之一截面进行研究，可以适当简化问题。 为方便与温度场电模拟实验的数据进行对比，对温度场的离散与之保持一致，即网格为0.1m的正方形，且对内、外墙表面的流体温度的模拟各设一排节点进行模拟。在实际的编程过程中，为使程序简洁且更有条理，建立一个网格矩阵，矩阵中每个元素的值根据温度场中对应位置的节点的边界情况或计算特点进行设定从而将温度场节点进行分类。网格矩阵如下。下面会结合该矩阵详细阐述迭代方程的建立过程。 三 方程的建立与求解 在介绍迭代方程之前，需要说明的是：a、为提高迭代速度，采用Gauss-