热流体课程实验报告-二维导热物体温度场的计算机模拟实验.docxVIP

成绩西安交通大学实验报告 成绩 课程： 实验日期 　　　　年　　月　　日 专业班号 组别 交报告日期 　　　　年　　月　　日 姓名 学号 报告退发 　　　　（订正、重做） 同组者 教室审批签字 实验名称二维导热物体温度场的计算机模拟实验 实验名称 一、实验目的 学习电、热类比的原理及边界条件的处理； 通过计算机编程的方式求出墙角导热的离散温度场。 二、实验原理 二维稳态过程，导热方程为 ? 二维稳态导热内部节点的差分方程为 t 于是内部节点的迭代计算式为 t 对于恒温边界条件，除了绝热边界时使用对称性外，只使用上面一个迭代计算式即可。但是对于对流边界，边界上的点，按位置分为内角点、外角点和平直边界，按类型分为对流边界、绝热边界，计算步骤相比恒温边界下更为复杂。 按位置： 内角点：4个方向均有导热热流，有dx2 外角点：2个方向有导热，有dx2 平直边界：3个方向有导热，有dx或dy面积的对流换热 按类型： 绝热边界：该点的绝热一侧没有热流量，基尔霍夫定律中，此方向的热流量代入0计算 对流边界：该点该方向的对流换热量由牛顿冷却公式q=hA 综上所述： 对流边界下的差分方程为： Φ 其中，Φi-1,j,Φi+1,j,Φi,j-1,Φ 代入所有q的计算式，可解得 t 注意： k为实际参与导热的几个方向，对于内角点有4项，外角点有2项，平直边界有3项，绝热边界还要去掉这一方向的那一项 Ak