传热实验报告_2.docVIP

PAGE 6 题目 有一个用砖砌成的长方形截面的冷空气通道，其界面尺寸如图1所示，假设在垂直于纸面方向上冷空气及砖墙的温度变化极小，可以近似的忽略不计。计算: (1).砖墙截面上的温度分布； (2).垂直于纸面方向的每米长度上通过砖墙的导热量。 已知：内外壁分别均匀地维持在0℃及30℃，砖墙导热系数λ=0.53W/m·℃ 图 图 SEQ 图 \* ARABIC 1 实验目的 (1). 学习数值模拟的原理和边界条件的处理； (2).通过数值模拟求出物体的温度场； (3).比较实际值与数值模拟结果之间的异同。 二、实验原理 已知当导热系数为常数时，对均匀的网络，二维稳态导热的差分方程为 Ti+1,j Ti+1,j Fortran90/95可以用于数字模拟，通过一定的控制结构，以一定的循环方式可以得到与实际结果相当接近的结果，这正是我们这里需要使用的。 在进行数字模拟的过程中我们先设定一个温度场的初场，假设均为0℃ ，然后在相应的边界点上建立边界条件。通过（公式1）可以建立中间点的温度场，通过（公式2）可以建立绝热边界上的温度场。 1．等温边界条件是最简单的情况，相应于热系统的等温边界只要是在数字模拟时，让边界上的各个点都等于某个固定值即可。 2.对于绝热条件，可以证明只要 2T(i-1,j)+T(i,j+1)+T(i,j-1)-4T(i,j)=0 即可使边界得到类似。 3.对于热流边界条件则只要 T(i,j+1)+T(i,j-1)+2T(i-1,j)+2hl/λT(i+1)-(4+2hl/λ)T(i,j)=0 经过上述变化之后，我们就能通过数字模拟对二维导热物体温度场进行模拟，从而得到我们想要的结果。 实验过程 问题简化 对于烟囱或者墙角的稳态导热问题，由于其具有对称性，因而可以取1/4部分进行研究。在对称面上a、b上，可认为符合绝热条件。在其内外表面上c、d上，则形成恒壁温条件。 y因此，墙角或烟囱的1/4部分的四个边界条件均为已知，对于角点再具体分析，则可分别列出满足各点要求的差分方程。由联立的差分方程组进行迭代计算，则可得到各点的近似温度值。 y 差分方程 2.1外部角点 2.2内部角点 2.3平直边界上节点 2.4网格节点示意 说明：1）本实验不存在内热源 2）等温边界条件不存在对流换热 3. 程序设计结果 利用上述原理，Fortran90程序如下：（截屏） 数值计算结果： 温度分布: 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 30.00 28.96 27.94 26.95 26.04 25.27 24.72 24.36 24.15 24.02 23.96 23.92 23.90 23.89 23.89 23.89 30.00 27.94 25.90 23.89 22.00 20.40 19.31 18.64 18.26 18.04 17.93 17.87 17.84 17.82 17.82 17.82 30.00 26.96 23.89 20.81 17.78 15.13 13.57 12.73 12.29 12.06 11.93 11.87 11.84 11.82 11.82 11.82 30.00 26.05 22.02 17.80 13.32 8.90 7.23 6.53 6.20 6.04 5.96 5.92 5.90 5.89 5.89 5.89 30.00 25.30 20.44 15.18 8.93 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 30.00 24.74 19.34 13.61 7.25 0.00 30.00 24.38 18.67 12.76 6.54 0.00 30.00 24.17 18.29 12.32 6.22 0.00 30.00 24.05 18.09 12.10 6.07 0.00 30.00 24.00 18.01 12.01 6.01 0.00 与热电比拟结果很接近，说明两种模拟方法均可行。 导热量： 外壁面：59.07681W/m 内壁面：59.06147W/m 等温线绘制： 实验总结 通过本次数值模拟实验我了解了传热学数值模拟的一般过程，对此有了更加直观的认识，也体会到数值模拟在现代热力计算中的方便快捷之处。不仅如此，我还进一步熟悉二维稳态温度场的分布特点和Fortran90在数值模拟方面的应用方法。