肋片散热分析—计算传热学课程设计肋片散热分析—计算传热学课程设计.docVIP

中国石油大学(华东) 储建学院热能与动力工程系 《计算传热学程序设计》 设计报告 学生姓名：龚波 学 号专业班级：热能与动力工程08－2班 指导教师：黄善波 2011年 7 月 5 日 1 设计题目 在工程实际中，往往需要增加（对流）传热量，应用比较广泛的较为有效的一种方法就是增加换热面积，即采用肋片—在材料消耗量增加较少的条件下能较多地增大换热面积。在一些换热设备中如制冷装置的冷凝器、散热器、空气加热器等等。某λ=λοk(T)=400(1+0.0035T),肋基温度Tw=95℃，肋表度黑度ε=0.80，周围空气温度Tf=20℃，环境辐射温度Ts=15℃，肋表面空气的表面换热系数hc=8W/(m2?℃)。 2 物理与数学模型 2.1 物理模型 发生在肋片的导热过程严格地说是多维的。如图1所示，暴露于恒温流体的圆柱肋片（肋高为L，直径为D）。由于圆柱直肋各处受热均匀，再加上肋片通常是由金属材料制成的，导热系数比较大，可以想象肋片内温度将仅沿肋高方向发生明显变化，再直径方向上变化相比很小。因此，假设该圆柱直肋在同一截面上温度相同，则该问题可转化为等截面直肋一维稳态导热问题。 图1 圆柱肋片物理模型图 2.2 数学模型 以肋基为坐标原点，圆柱肋片厚度方向为坐标正方向，建立坐标系如图2所示。 基于上述物理模型，则该问题的数学模型可描述如下： （1-a） 左右两侧相应的边界条件分别是第一类边界条件和第二类边界条件，分别描述如下： 左边界 （1-b） 右边界 （1-c） 图2 圆柱肋片数学模型图 3 数值处理与程序设计 3.1数学模型无量纲化 为了使数值计算结果具有更普遍的意义，将上述数学模型无量纲化。为此定义 ， （2） 控制方程无量纲化后，方程整理为 （3） 定义 ，， ， ，， （4）将上述定义带入式（3） （5-a） 左边界 （5-b) 右边界 （5-c） 3.2 试射法的形式 令 ， （6） 则有试射法形式模型 （7-a） (7-b) 左边界 （7-c） 其中，P1=1,Q1=0,W1=1 右边界 （7-d） 其中，P2=0,Q2=1,W2=0 3.3 程序编写 圆柱直肋一维稳态导热数学模型是二阶常微分两点边值问题，可以采用试射法求解。其基本思想是将边值问题转换为初值问题求解。 3.3.1 设计特点 在主程序外设置全局变量，为使在调用各子程序时，不会因实参与形参的作用范围而无法编译、运行程序。 在主程序头部，对参数赋值，对体积和肋高赋值应注意范围和两者的关联性。此处赋值V=0.00002m3，L=0.5m，保证程序结果为最大传热量，而且保证了足够的计算空间又不至于过分浪费系统资源。 利用循环实现计算最大传热量的过程，首先调用肋高函数得到按线性规律递减的肋高，再调用shoot函数计算相应肋高时的肋基温度梯度，调用热量函数求解热量Q[g]，输出各个肋高下的肋基温度梯度和热量，为了便于了解热量随肋高的变化关系。比较各肋高下的热量值，将最大热量值对应下标保留。然后，输出最大热量Q[max]和相应的肋高L