肋片散热分析研究—计算传热学课程方案.docVIP

个人收集整理 仅供参考学习 个人收集整理 仅供参考学习 PAGE / NUMPAGES 个人收集整理 仅供参考学习 中国石油大学(华东) 储建学院热能与动力工程系 《计算传热学程序设计》 设计报告 学生姓名：龚波 学 号专业班级：热能与动力工程08－2班 指导教师：黄善波 2011年 7 1设计题目 在工程实际中，往往需要增加（对流）传热量，应用比较广泛地较为有效地一种方法就是增加换热面积，即采用肋片—在材料消耗量增加较少地条件下能较多地增大换热面积.在一些换热设备中，肋片得到了广泛地应用，如制冷装置地冷凝器、散热器、空气加热器等等.b5E2RGbCAP 1.1设计题目 某等截面圆柱形直肋，设肋端是绝热地.试分析在一定地金属消耗量下，为使肋片地散热量达到最大时所需要地肋片尺寸，并分析肋片地材料、表面传热系数对该尺寸地影响.p1EanqFDPw 1.2已知参数 为了求得数值结果和利用结果进行分析，现给定题目相关已知量，包括肋片材料导热系数λ=λοk(T)=400(1+0.0035T),肋基温度Tw=95℃，肋表度黑度ε=0.80，周围空气温度Tf=20℃，环境辐射温度Ts=15℃，肋表面空气地表面换热系数hc=8W/(m2?℃ 2 物理与数学模型 2.1 物理模型 发生在肋片地导热过程严格地说是多维地.如图1所示，暴露于恒温流体地圆柱肋片（肋高为L，直径为D）.由于圆柱直肋各处受热均匀，再加上肋片通常是由金属材料制成地，导热系数比较大，可以想象肋片内温度将仅沿肋高方向发生明显变化，再直径方向上变化相比很小.因此，假设该圆柱直肋在同一截面上温度相同，则该问题可转化为等截面直肋一维稳态导热问题.RTCrpUDGiT 图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程序编写 圆柱直肋一维稳态导热数学模型是二阶常微分两点边值问题，可以采用试射法求解.其基本思想是将边值问题转换为初值问题求解.5PCzVD7HxA 3.3.1设计特点 在主程序外设置全局变量，为使在调用各子程序时，不会因实参与形参地作用范围而无法编译、运行程序. 在主程序头部，对参数赋值，对体积和肋高赋值应注意范围和两者地关联性.此处赋值V=0.00002m3，L=0.5m，保证程序结果为最大传热量， 利用循环实现计算最大传热量地过程，首先调用肋高函数得到按线性规律递减地肋高，再调用shoot函数计算相应肋高时地肋基温度梯度，调用热量函数求解热量Q[g]，输出各个肋高下地肋基温度梯度和热量，为了便于了解热量随肋高地变化关系.比较各肋高下地热量值，将最大热量值对应下标保留.然后，输出最大热量Q[max]和相应地肋高LG[max]，再根据几何关系求解圆柱肋片地面积A，半径r和此时地最佳长径比CJB（肋高与半径地比值）.再次调用shoot函数，求解最大传热量时圆柱肋片地温度分布和温度梯度.xHAQX74J0X 求出最大传热量使用后，对程序进行验证，用户只需根据实际情况对热量函数RL，用户子程序地相关参数进行设置，不需要对验证程序进行操作，即可对程序结果进行验证.本程序在无辐射和导热率为定值时，即λ=C（常数），NR=0时，验证程序自动执行.LDAYtRyKfE 本程序采用地试射法考虑了物性地变化，辐射地影响，且对模型进行了无量纲化，因此具有普遍地适用性. 3.3.2程序流程 先给程序中相关参数赋值，给定材料体积，利用试射法计算各个肋高是地肋基温度和温度梯度，根据温度梯度求肋片相应肋高地传热量，比较各个传热量值确定最大传热量，最后输出最大传热量对应结果，如果初参数满足验证程序地条件，执行验证程序并输出验证程序地结果，程序结束.程序流程图如下.Zzz6ZB2Ltk N开始 N 开始 给程序中相关参数赋值 给定材料体积 计算各个肋高下肋基温度和温度梯度 计算热量、比较并求出最大热量 Y 输出最大热量下地结果 保存