墙角传热上机报告(共5页).doc

西 安 交 通 大 学 实 验 报 告课 程 传热学上机实验实验名称 墙角导热问题的数值计算 实验时间 学 院 能源与动力工程学院 姓 名 学 号 班 级 一、问题描述有一墙角模型，尺寸如图1所示，导热系数0.53W/(mK)，墙角内外壁为第一类边界条件。求解该模型的温度分布及导热量。图1二、计算原理 根据热平衡法列出节点方程，各方向导入单元体的热量之和为零。内节点和绝热边界点（图1点划线上的点）的方程形式不同。图2图2所示的内节点和绝热边界节点方程如下：内节点： 绝热边界点：三、计算过程等温边界条件：内壁0，外壁30说明：1.计算内部节点时， t(i,j)=0.25*(t(i-1,j)+t(i+1,j)+t(i,j-1)+t(i,j+1)，计算边界节点时此式稍作变动。2.计算导热量时，先从外表面计算，再从内表面计算，最终结果为二者的平均值。程序：program xxqiangjiao