《传热学模拟大实验-二维导热温度场的数值计算》-毕业设计（论文）.docxVIP

传热大模拟实验 PAGE1 / NUMPAGES20 传热学模拟大实验 —二维导热温度场的数值计算 学院：能动学院 名字： 班级： 专业：能源动力系统及自动化 学号： 2012.12.30 一、实验题目： 有一个用砖砌成的长方形截面的冷空气空道，其截面尺寸如下图所示，假设在垂直于纸面方向上冷空气及砖墙的温度变化很小，可以近似地予以忽略。在下列两种情况下试计算：（1）砖墙横截面上的温度分布；（2）垂直于纸面方向的每米长度上通过砖墙的导热量。外矩形长为3.0m，宽为2.2m；内矩形长为2.0m，宽为1.2m。 第一种情况：内外壁分别均匀地维持在及（等温边界条件）； 第二种情况：内外表面均为第三类边界条件（对流边界条件），且已知： 外壁： , 内壁： , 砖墙的导热系数； 由于对称性，仅研究1/4部分即可。 二、数学描述 1.基本方程 此问题为无内热源的二维稳态导热问题，所以可得出方程： 2.边界条件 对边界条件进行分析，对于边界恒温的第一种情况有： 边界绝热；q=0，且， 对于第三类边界条件，可得： 三、方程离散 将上图进行区域离散，且，并建立内外节点的离散方程 第一种情况： 内点： 边界点： ab边界（a、b除外）i=2，3，4，5； j=1 de边界（d、e除外）i=16; j=8，9，10，11 af、fe边界(包括a、f、e点） i=1；j=1，2，3，…，12 j=12；i=1，2，3，…，16 bc、cd边界（包括b、c、d点） i=6；j=1，2，…，7 j=7；i=6，8，…，16 第二种情况： 内点： 边界点： ab边界（a、b除外）i=2，3，4，5； j=1 de边界（d、e除外）i=16; j=8，9，10，11 af、fe边界(不包括a、f、e点） i=1；j=2，3，…，11 j=12；i=2，3，…，15 bc、cd边界（不包括b、c、d点） i=6；j=2，3，…，6 j=7；i=7，8，…，15 a点： b点： c点： d点： e点： f点： 四、编程方程式 采用高斯-塞德尔迭代法 设立初值,通过程序计算与上一次的值进行对比,设立需要的最大差值,进行多次迭代,最后得出结果.（具体编程式见后附） 五、实验结果 1.第一种情况:(式中k为迭代次数) k = 151 s（温度场分布） = Columns 1 through 14 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 30.0000 29.0345 28.0690 27.1174 26.2214 25.4556 24.8940 24.5255 24.3001 24.1684 24.0935 24.0517 24.0289 24.0167 30.0000 28.0691 26.1241 24.1790 22.3129 20.7068 19.5951 18.9077 18.5067 18.2799 18.1538 18.0846 18.0470 18.0272 30.0000 27.1176 24.1793 21.1619 18.1442 15.4636 13.8718 13.0037 12.5391 12.2907 12.1574 12.0858 12.0474 12.0273 30.0000 26.2220 22.3137 18.1449 13.6384 9.1317 7.4249 6.6960 6.3554 6.1863 6.0994 6.0536 6.0295 6.0169 30.0000 25.4568 20.7088 15.4655 9.1328 0 0 0 0 0 0 0 0 0 30.0000 24.8964 19.5990 13.8757 7.4272 0 0 0 0 0 0