计算传热学课程设计.docVIP

中国石油大学(华东) 储建学院热能与动力工程系 《计算传热学程序设计》 设计报告 学生姓名： 学 号： 专业班级：热能与动力工程 指导教师： 2011年 7 月 5 日 目录 1 设计题目 2 2 物理与数学模型 2 2.1 物理模型 2 2.2 数学模型 2 3 数值处理与程序设计 3 3.1 数值处理 3 3.1.1 数值处理与方程离散化 3 3.1.2 离散方程的导出 4 3.2 程序设计 6 4 模型与程序验证 7 5 计算结果与分析 9 6 结论 10 7 参考文献 11 8 附录 11 附录一 数据表格 11 附录二 主要程序 12 1 设计题目 有一冷柜初始时处于均匀的温度20℃。开启压缩机后，冷冻室（包括柜门的内表面）以均匀地速度降温，温降速度为18℃/h。柜门尺寸为1.2x1.2㎡，保温材料厚0.08m，导热系数为0.04W/(m·℃)。假设外表面包层很薄，其热阻可以忽略。柜门外受空气自然对流和环境的辐射换热。自然对流可按下式计算， W/(㎡·℃) 其中L为门高。表面发射率为0.8。通过柜门的导热可按一维问题处理。试计算压缩机启动两个小时内通过柜门的冷量损失。 22.1 物理模型 由题可知，柜门左边界温度均匀下降，右边界置于恒温中，并有对流和辐射散热，由于柜门沿高度方向没有散热，所以沿高度方向没有温度的变化，即只有沿厚度方向的温度变化，所以可以视为一维问题处理，同时可以得出左右边界都是第三类边界情况。 2.2以柜门左侧为原点，沿柜门厚度方向为坐标轴正方向建立坐标轴。则有δ=0.08m，基于上述物理模型，只对厚度x方向分析，可以得到该问题的数学模型为 (1)初始条件为： (2) 边界条件为： (3) 图1 柜门模型 33.1 数值处理 3.1.1 数值处理与方程离散化 采用外节点法用均匀网格对求解区域进行离散化，得到的网格系统如图2所示。 图中为了和C语言中数组定义一致，节点标号从0到N-1一共有N个节点。节点间距δx为 图2 宽度x方向的网格划分简化图 采用有限差分法或者有限体积法对控制方程进行离散化，得到的离散方程为 (4a) 式中： (4b) ， , (4c) (4d) 其中，上标“0”表示上一时刻值，λ为节点所在控制容积界面上的导热系数，Δτ为时间步长。由于柜门的材料的λ为常数，在离散方程中取值时，λw=λe=λ=0.04W/(m℃)。 由于左侧界面柜门温度均匀下降，故有: (5) 由题意知右边界为第三类边界条件，故有 (6) 由元体能量平衡法与源项线性化可以得到右边界节点的离散方程，为： (7) 其中 由于柜门的材料的λ为常数，在离散方程中取值时，λw=λe=λ。 3.1.2为了得到离散方程，除了要对空间区域进行离散化之外，还需要对时间坐标进行离散化。对时间进行离散化：从初始时刻开始，取定一个时间步长Δτ构成如下的时间序列即可 ， （p=0，1，2，3，…，） (8) 其中p=0对应于初始时刻，后面为方便起见将τp时刻简称为p时刻。 内部节点P及其与邻点的关系仍如图2所示。求解区域内的任一点P，在p时刻，它仍然满足控制方程，即 (9) 上述方程中，温度对空间坐标的导数项可以表示为 (10) 按照Taylor级数展开法，温度对时间的偏导数采用向后差分格式： (11) 将(12)带入到差分方程(11)中去，并定义 ，， (12) 整理后得到隐式差分格式的离散方程为： (13) 上式即是该题目中一维非稳态导热方程的离散方程。 下面来分别计算左右两边界的边界条件： 对于左边界节点，温度随