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

1、中国石油大学 (华东)储建学院热能与动力工程系 计 算 传 热 学 程 序 设 计 设 计 报 告学生姓名：龚波学 号：08123217专业班级：热能与动力工程 082 班指导教师：黄善波2011 年 7 月 5 日11 设计题目在工程实际中，往往需要增加（对流）传热量，应用比较广泛的较为有效的一种方法就是增加换热面积，即采用肋片在材料消耗量增加较少的条件下能较多地增大换热面积。在一些换热设备中，肋片得到了广泛地应用，如制冷装置的冷凝器、散热器、空气加热器等等。1.1 设计题目某等截面圆柱形直肋，设肋端是绝热的。试分析在一定的金属消耗量下，为使肋片的散热量达到最大时所需要的肋片尺寸，并分析肋片

2、的材料、表面传热系数对该尺寸的影响。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） 。由于圆柱直肋各处受热均匀，再加上肋片通常是由金属材料制成的，导热系数比较大，可以想象肋片内温度将仅沿肋高方向发生明显变化，再直径方

3、向上变化相比很小。因此，假设该圆柱直肋在同一截面上温度相同，则该问题可转化为等截面直肋一维稳态导热问题。 dt/x0Htftfws图 1 圆柱肋片物理模型图22.2 数学模型以肋基为坐标原点，圆柱肋片厚度方向为坐标正方向，建立坐标系如图 2 所示。基于上述物理模型，则该问题的数学模型可描述如下：（1-40cfbsdTAUhTTxa）左右两侧相应的边界条件分别是第一类边界条件和第二类边界条件，分别描述如下：左边界 （1-0wxTb）右边界 （1-0xLdTc） 图 2 圆柱肋片数学模型图3 数值处理与程序设计3.1 数学模型无量纲化为了使数值计算结果具有更普遍的意义，将上述数学模型无量纲化。为此

4、定义， xLfwT（2）控制方程无量纲化后，方程整理为34422 31 0fscbwfwwfTdkdULhTA Txx （3） 定义 ， ， ， 1dkLxUSAchLBi， ， 3bwfRTNffwTsswfT（4）将上述定义带入式（3）中，整理得：（5-22 40LxdcRfsSkBiNkxa）左边界 （5-b)01x右边界 （5-1xdc）3.2 试射法的形式令， 1y12dyx（6）则有试射法形式模型（7-12dyxa）(7-b)42211LxdcRfsSykBiyNx4左边界（7-121PyQWc）其中，P 1=1,Q1=0,W1=1右边界（7-212yd）其中，P 2=0,Q2=1

5、,W2=03.3 程序编写圆柱直肋一维稳态导热数学模型是二阶常微分两点边值问题，可以采用试射法求解。其基本思想是将边值问题转换为初值问题求解。3.3.1 设计特点在主程序外设置全局变量，为使在调用各子程序时，不会因实参与形参的作用范围而无法编译、运行程序。在主程序头部，对参数赋值，对体积和肋高赋值应注意范围和两者的关联性。此处赋值 V=0.00002m3，L=0.5m，保证程序结果为最大传热量，而且保证了足够的计算空间又不至于过分浪费系统资源。利用循环实现计算最大传热量的过程，首先调用肋高函数得到按线性规律递减的肋高，再调用 shoot 函数计算相应肋高时的肋基温度梯度，调用热量函数求解热量Q

6、g，输出各个肋高下的肋基温度梯度和热量，为了便于了解热量随肋高的变化关系。比较各肋高下的热量值，将最大热量值对应下标保留。然后，输出最大热量 Qmax和相应的肋高 LGmax，再根据几何关系求解圆柱肋片的面积 A，半径 r 和此时的最佳长径比 CJB（肋高与半径的比值） 。再次调用 shoot 函数，求解最大传热量时圆柱肋片的温度分布和温度梯度。求出最大传热量使用后，对程序进行验证，用户只需根据实际情况对热量函数RL，用户子程序的相关参数进行设置，不需要对验证程序进行操作，即可对程序结果进行验证。本程序在无辐射和导热率为定值时，即 =C（常数） ，N R=0 时，验证程序自动执行。本程序采用的

7、试射法考虑了物性的变化，辐射的影响，且对模型进行了无量纲化，因此具有普遍的适用性。3.3.2 程序流程5先给程序中相关参数赋值，给定材料体积，利用试射法计算各个肋高是的肋基温度和温度梯度，根据温度梯度求肋片相应肋高的传热量，比较各个传热量值确定最大传热量，最后输出最大传热量对应结果，如果初参数满足验证程序的条件，执行验证程序并输出验证程序的结果，程序结束。程序流程图如下。N开始给程序中相关参数赋值给定材料体积计算各个肋高下肋基温度和温度梯度计算热量、比较并求出最大热量 Y输出最大热量下的结果保存结果至文件判断初参数YN程序结束执行验证程序并输出结果图 3 程序流程图4 模型与程序的验证4.1

8、模型验证为了方便利用解析解验证程序，将本题简化为常物性、无辐射等截面直肋一维稳态导热模型。已知肋片材料导热系数 =100 W/(m ),肋基温度 Tw=95 ，周围空气温度 Tf=20 ，肋表面空气的表面换热系数 hc=8W/(m2 )。建立坐标系，列出其控制方程式及定解条件：(8-a)2dx(8-b)0,;wfxt,0dHx其中过余温度 ， 为一常量。ft/cmhpA6式（8-a）是一个二阶线性微分方程，由两边界条件可求出精确解为0chmx（9）4.2 程序验证将式（9）中参数换算成无量纲形式，然后编程，计算出每个节点温度的解析解（验证程序见附录）和数值解（验证源程序见附录） ，进行比较，如

9、表格 1。表 1 =100、无辐射圆柱直肋无量纲温度值数值解和分析解x 数值解 理论解 百分误差/%0 1 1 00.05 0.966899 0.966899 00.1 0.935906 0.935905 0.000106850.15 0.906951 0.906951 00.2 0.879973 0.879973 00.25 0.854913 0.854912 0.000116970.3 0.831715 0.831715 00.35 0.81033 0.810329 0.000123410.4 0.79071 0.79071 00.45 0.772814 0.772813 0.000129

10、40.5 0.756602 0.756601 0.000132170.55 0.742038 0.742037 0.000134760.6 0.729091 0.72909 0.000137160.65 0.717733 0.717732 0.000139330.7 0.707939 0.707938 0.000141260.75 0.699687 0.699686 0.000142920.8 0.69296 0.692959 0.000144310.85 0.687744 0.687742 0.000290810.9 0.684026 0.684024 0.000292390.95 0.68

11、1798 0.681797 0.000146671 0.681056 0.681055 0.00014683700.20.40.60.811.20 0.3 0.6 0.9 1.2无 量 纲 位 置无量纲温度数 值 解理 论 解图 4 圆柱直肋无量纲温度分布曲线 由上述图表可知圆柱肋片分析解和数值解相差不大，二者吻合较好，可以说明所编制的数值解法的程序是正确的。5 计算结果与分析5.1 肋高与热量的关系材料的导热率 =400(1+0.0035T)，肋基温度 Tw=95 ，肋表度黑度 =0.80，周围空气温度 Tf=20 ，环境辐射温度 Ts=15 ，肋表面空气的表面换热系数 hc=8W/(m2

12、)时，圆柱肋片肋基无量纲温度梯度和传热量见表 2。表 2 不同长度下肋片的传热量 任意长度 L/m 肋基温度梯度 热量 /W0.5 10 -54.92460.49 10 -57.18930.48 -1.98619 11.83710.47 -1.92857 11.9880.46 -1.87095 12.1410.45 -1.81331 12.29570.4 -1.52509 13.08830.39 -1.46747 13.24790.36 -1.29491 13.71960.35 -1.23757 13.87210.34 -1.18038 14.02080.33 -1.12339 14.1647

13、0.32 -1.06663 14.30280.27 -0.78888 14.85980.26 -0.73514 14.93240.25 -0.68225 14.98890.24 -0.63034 15.02660.23 -0.57955 15.04320.22 -0.53001 15.03660.21 -0.48189 15.00440.2 -0.43533 14.94410.19 -0.39049 14.85290.18 -0.34756 14.72950.17 -0.30666 14.57040.1 -0.08934 12.26740.07 -0.03744 10.49210.06 -0.

14、0256 9.76310.05 -0.01629 8.94650.04 -0.00935 8.02070.03 -0.00456 6.95490.02 -0.00166 5.69630.01 -0.00029 3.9285从表中结果易看出在肋高为 0.50m 和 0.49m 的时候，肋基温度梯度为正值，传热量为负值，与实际情况不符。这是因为肋高增加，一定的耗材下，肋片直径变小，对流换热量处理成广义热源已不合适，即不能作为一维稳态导热模型看待。此时，增大肋片体积、增加导热率或减小对流换热系数，又能满足模型使用调节。忽略表中前两行的数据绘图见图 5。02468101214160 0.2 0.4 0

15、.6实 际 肋 高 /m传热量/W图 5 不同长度下肋片传热量曲线由图可知传热量随着肋高先增后减，传热量最大在肋高 L=0.23m 取得。因为 =-Ay0 1，在肋基，温度始终为 tw，即导热系数不变，肋基温度梯度 y01为负且和圆柱截面积 A 随肋高增加而变小，所以存在最佳肋高使传热量最大。95.2 表面换热系数的影响材料的导热率 =400(1+0.0035T)，肋表度黑度 =0.80，圆柱肋片在不同表面换热系数 h 时，为使传热量最大，相应最佳肋高 LGmax和最佳长径比 CJB（肋高与半径的比值）见表 3。表 3 不同表面传热量下的 LGmax和 CJBh/(W/m2) 4 5 7 9

16、10 12LGmax/m 0.15 0.14 0.13 0.13 0.12 0.12CJB 23.019014 20.755915 18.572264 18.572264 16.471066 16.471066由表易知，随着表面传热系数的增加，最佳肋高是逐渐减小的。表面传热系数的增加，传热量增加，由 =-Ay0 1知，需要增大肋基导热面积，所以最佳肋高减小。5.3 材料导热率的影响无辐射，肋表面空气的表面换热系数 hc=8W/(m2 )时，不同导热系数 时，为使传热量最大，相应最佳肋高 LGmax和最佳长径比 CJB（肋高与半径的比值）见表4。表 4 不同传热系数下的 LGmax和 CJB/(

17、W/mk) 100 300 400 600 700 800LGmax/m 0.18 0.28 0.31 0.36 0.39 0.41CJB 30.25928 58.706592 68.389963 85.586167 96.504316 104.022099由表易知，随着导热系数的增加，最佳肋高是逐渐增大的。因为导热系数变大，传热量增加，由 =hcA(t-tf)知，需增加圆柱侧面积以加强换热，所以最佳肋高增加。6 结论在肋基，温度始终维持不变，即导热系数不变，肋基温度梯度为负且和圆柱截面积 A 随肋高增加而变小，由傅里叶公式可知存在最佳肋高使肋片传热量最大，在题目已知条件下，当肋高 L=0.2

18、3m 时取得最大散热量 =15.0432W；表面传热系数的增加，传热量增加，由傅里叶公式知，需要增大肋基导热面积，所以最佳肋高减小；导热系数变大，传热量增加，由对流换热公式知，需增加圆柱侧面积以加强换热，所以最佳肋高增加。参考文献1 黄善波，刘中良.计算传热学基础.中国石油大学（华东）热能与动力工程系，2009102 杨世铭，陶文铨.传热学（第四版）.高等教育出版社，2007附录 1 主要程序表 5 程序列表序号 程序名称 程序功能 对应图表1 全功能程序 按题目要求，输出各个肋高下的传热量，求解最大传热量以及此时肋片的尺寸，输出温度和温度梯度的分布，验证程序图 3、 4、5，表 1、2 、3

19、、4已知参数赋值/输入圆柱肋体积 V、任意给定肋高 L 及肋高变化步长 bcV=0.00002;L=0.5; bc=0.01;void fct(int N,double x,double y,double f)/函数子程序，用户根据具体条件进行修改double Kd,Slx,Bic;Kd=0.2625/(2.026025+0.2625*y0);Slx=2*LGZ*sqrt(3.14*LGZ/V);Bic=0.08*LGZ/4.0;Nr=0.000191*LGZ/4.0;f0=y1; /f0=dy1/dxf1=-Kd*y1*y1+Slx*(Bic*y0+Nr*(pow(y0+3.908667),

20、4)-pow(3.868667,4);return;void pqw1(double Y,double *P,double *Q,double *W)/左边界处的第三类边界条件(x=xa)/P1*y1+Q1*y2=W1-用户应根据具体条件进行修改11*P=1.0;*Q=0.0;*W=1;void pqw2(double Y,double *P,double *Q,double *W)/右边界处的第三类边界条件(x=xb)/P2*y1+Q2*y2=W2-用户应根据具体条件进行修改*P=0;*Q=1;*W=0;return;求最大传热量/调用各个函数求最大传热量 for(g=0;L-g*bc0;g

21、+)/用肋高函数求肋高LGk=LeiGao(L,bc,k);LGZ=LGk;x=xa;/利用试射法确定 mshoot(N,NS,x,h,M1,M2,Eps,y0,y);/求热量Qk=RL(V, y0);LGg=LeiGao(L,bc,g);LGZ=LGg;x=xa;/利用试射法确定 mshoot(N,NS,x,h,M1,M2,Eps,y0,y);12Qg=RL(V, y0);/输出任意长度及所对应的热量 printf(“%f “,LGg); /输出长度fprintf(fp,“%f “,LGg); printf(“%f “,y01); /输出肋基出温度梯度fprintf(fp,“%f “,y01

22、);printf(“%6.4f n“,Qg); /输出对应传热量fprintf(fp,“%6.4f n“,Qg); /保存最大传热量if(QkQg)k=g;/输出最大传热量，并输出对应圆肋的尺寸max=k;LGmax=LeiGao(L,bc,max);LGZ=LGmax;A=V/LGmax;r=sqrt(A/3.14);CJB=LGmax/r;printf(“最佳肋高 LGmax=%fn“,LGZ);fprintf(fp,“最佳肋高 LGmax=%fn“,LGZ);printf(“最佳面积 A=%fn“,A);fprintf(fp,“最佳面积 A=%fn“,A);printf(“最佳半径 r=

23、%fn“,r);fprintf(fp,“最佳半径 r=%fn“,r);printf(“最佳长径比 CJB=%fn“,CJB);fprintf(fp,“最佳长径比 CJB=%fn“,CJB);x=xa;/利用试射法确定 mshoot(N,NS,x,h,M1,M2,Eps,y0,y);Qmax=RL(V, y0);/求最大传热量13printf(“最大热量 Qmax=%6.4f “,Qmax);fprintf(fp,“最大热量 Qmax=%6.4f “,Qmax);printf(“nnn“);fprintf(fp,“nnn“);/输出最大传热量时的温度分分布printf(“输出最大传热量时的温度分

24、布n“);/显示在屏幕上fprintf(fp,“输出最大传热量时的温度分布n“);/保存到文件中/输出表头printf(“ x y1 y2n“);/显示在屏幕上fprintf(fp,“ x y1 y2 n“);/保存到文件中/输出 x=a 时的结果printf(“%6.4f “,xa);fprintf(fp,“%6.4f “,xa);for(i=0;iN;i+)printf(“%10.6f “,y0i);fprintf(fp,“%10.6f “,y0i);printf(“n“);fprintf(fp,“n“);x=xa;/调用 R-K 方法计算并输出后续各点的值for(i=0;iN;i+)yi

25、=y0i;for(j=0;jNS;j+)rungek(N, /根据求出的 m 解决问题printf(“%6.4f “,x);fprintf(fp,“%6.4f “,x);for(i=0;iN;i+)14printf(“%9.6f “,yi);fprintf(fp,“%9.6f “,yi);printf(“n“);fprintf(fp,“n“);验证程序/令 lmd=常数，无辐射，验证程序的正确性if(kd=0.0fprintf(fp,“nn 令 lmd=100，无辐射，验证程序理论温度n“);MM=sqrt(0.16/r);MMH=MM*LGZ;printf(“ x y1 y2n“);/显示在

26、屏幕上fprintf(fp,“ x y1 y2 n“);/保存到文件中x=0.0000;for(i=0;h*i=1.0;i+) x=h*i;printf(“%6.4f “,x);fprintf(fp,“%6.4f “,x);LLWD=(exp(MMH*(x-1)+exp(MMH*(1-x)/(exp(MMH)+exp(-1.0*MMH);/温度分布printf(“%9.6f “,LLWD);fprintf(fp,“%9.6f “,LLWD);printf(“n“);fprintf(fp,“n“);15附录 2 数学模型的无量纲化过程推导针对式（1）进行无量纲化处理，为此定义， xLfwT（10

27、）其中 Tf、Tw 均为常数，假定 =k(T)，则无量纲化过程如下。（112dTdTAAxxx 22ddTAx） （12wffwfL-a） (12 2fTdTdxx2-b) （12-cwfwffdkdkxTT）22wfwfwfAddkdALTLxx(13)40cfbffsUhT整理(14)4422 31 0fscbwfwwfTdkdLhATxx 定义 ， ， ， 1dkLxUSAchLBi， ， (15)3bwfRTNffwTsswfT16将上述定义带入式（3）中，整理得：(16-a) 22 40LxdcRfsSkBiNkx左边界 (16-b)01x右边界 (16c)1xd热能与动力工程系计算传热学程序设计成绩考核表指标 考核内容 分值 得分1模型、方法和计算结果的可靠性（含程序考核）302 讨论、分析的充分性、详实性 253 报告格式的规范性 204 报告内容的阐述、问答问题的情况 105 平时的表现 15176 合计 100教师签字： 此页单独占一页！