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

中国石油大学(华东)储建学院热能与动力工程系计算传热学程序设计设计报告学生姓名：龚波学 号：08123217专业班级：热能与动力工程082班指导教师：黄善波 2011年 7 月 5 日1 设计题目在工程实际中，往往需要增加（对流）传热量，应用比较广泛的较为有效的一种方法就是增加换热面积，即采用肋片在材料消耗量增加较少的条件下能较多地增大换热面积。在一些换热设备中，肋片得到了广泛地应用，如制冷装置的冷凝器、散热器、空气加热器等等。1.1 设计题目某等截面圆柱形直肋，设肋端是绝热的。试分析在一定的金属消耗量下，为使肋片的散热量达到最大时所需要的肋片尺寸，并分析肋片的材料、表面传热系数对该尺寸的影响。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）。由于圆柱直肋各处受热均匀，再加上肋片通常是由金属材料制成的，导热系数比较大，可以想象肋片内温度将仅沿肋高方向发生明显变化，再直径方向上变化相比很小。因此，假设该圆柱直肋在同一截面上温度相同，则该问题可转化为等截面直肋一维稳态导热问题。图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=03.3 程序编写圆柱直肋一维稳态导热数学模型是二阶常微分两点边值问题，可以采用试射法求解。其基本思想是将边值问题转换为初值问题求解。3.3.1 设计特点在主程序外设置全局变量，为使在调用各子程序时，不会因实参与形参的作用范围而无法编译、运行程序。在主程序头部，对参数赋值，对体积和肋高赋值应注意范围和两者的关联性。此处赋值V=0.00002m3，L=0.5m，保证程序结果为最大传热量，而且保证了足够的计算空间又不至于过分浪费系统资源。利用循环实现计算最大传热量的过程，首先调用肋高函数得到按线性规律递减的肋高，再调用shoot函数计算相应肋高时的肋基温度梯度，调用热量函数求解热量Qg，输出各个肋高下的肋基温度梯度和热量，为了便于了解热量随肋高的变化关系。比较各肋高下的热量值，将最大热量值对应下标保留。然后，输出最大热量Qmax和相应的肋高LGmax，再根据几何关系求解圆柱肋片的面积A，半径r和此时的最佳长径比CJB（肋高与半径的比值）。再次调用shoot函数，求解最大传热量时圆柱肋片的温度分布和温度梯度。求出最大传热量使用后，对程序进行验证，用户只需根据实际情况对热量函数RL，用户子程序的相关参数进行设置，不需要对验证程序进行操作，即可对程序结果进行验证。本程序在无辐射和导热率为定值时，即=C（常数），NR=0时，验证程序自动执行。本程序采用的试射法考虑了物性的变化，辐射的影响，且对模型进行了无量纲化，因此具有普遍的适用性。3.3.2 程序流程先给程序中相关参数赋值，给定材料体积，利用试射法计算各个肋高是的肋基温度和温度梯度，根据温度梯度求肋片相应肋高的传热量，比较各个传热量值确定最大传热量，最后输出最大传热量对应结果，如果初参数满足验证程序的条件，执行验证程序并输出验证程序的结果，程序结束。程序流程图如下。N开始给程序中相关参数赋值给定材料体积计算各个肋高下肋基温度和温度梯度计算热量、比较并求出最大热量Y输出最大热量下的结果保存结果至文件判断初参数YN程序结束执行验证程序并输出结果图3 程序流程图4 模型与程序的验证4.1 模型验证为了方便利用解析解验证程序，将本题简化为常物性、无辐射等截面直肋一维稳态导热模型。已知肋片材料导热系数=100 W/(m),肋基温度Tw=95，周围空气温度Tf=20，肋表面空气的表面换热系数hc=8W/(m2)。建立坐标系，列出其控制方程式及定解条件： (8-a) (8-b)其中过余温度，为一常量。 式（8-a）是一个二阶线性微分方程，由两边界条件可求出精确解为 （9）4.2 程序验证将式（9）中参数换算成无量纲形式，然后编程，计算出每个节点温度的解析解（验证程序见附录）和数值解（验证源程序见附录），进行比较，如表格1。表1 =100、无辐射圆柱直肋无量纲温度值数值解和分析解x数值解理论解百分误差/%01100.050.9668990.96689900.10.9359060.9359050.000106850.150.9069510.90695100.20.8799730.87997300.250.8549130.8549120.000116970.30.8317150.83171500.350.810330.8103290.000123410.40.790710.7907100.450.7728140.7728130.00012940.50.7566020.7566010.000132170.550.7420380.7420370.000134760.60.7290910.729090.000137160.650.7177330.7177320.000139330.70.7079390.7079380.000141260.750.6996870.6996860.000142920.80.692960.6929590.000144310.850.6877440.6877420.000290810.90.6840260.6840240.000292390.950.6817980.6817970.0001466710.6810560.6810550.00014683 图4 圆柱直肋无量纲温度分布曲线 由上述图表可知圆柱肋片分析解和数值解相差不大，二者吻合较好，可以说明所编制的数值解法的程序是正确的。5 计算结果与分析5.1 肋高与热量的关系材料的导热率=400(1+0.0035T)，肋基温度Tw=95，肋表度黑度=0.80，周围空气温度Tf=20，环境辐射温度Ts=15，肋表面空气的表面换热系数hc=8W/(m2)时，圆柱肋片肋基无量纲温度梯度和传热量见表2。表2 不同长度下肋片的传热量 任意长度L/m肋基温度梯度热量/W0.510-54.92460.4910-57.18930.48-1.9861911.83710.47-1.9285711.9880.46-1.8709512.1410.45-1.8133112.29570.4-1.5250913.08830.39-1.4674713.24790.36-1.2949113.71960.35-1.2375713.87210.34-1.1803814.02080.33-1.1233914.16470.32-1.0666314.30280.27-0.7888814.8590.26-0.7351414.93240.25-0.6822514.98890.24-0.6303415.02660.23-0.5795515.04320.22-0.5300115.03660.21-0.4818915.00440.2-0.4353314.94410.19-0.3904914.85290.18-0.3475614.72950.17-0.3066614.57040.1-0.0893412.26740.07-0.0374410.49210.06-0.02569.76310.05-0.016298.94650.04-0.009358.02070.03-0.004566.95490.02-0.001665.69630.01-0.000293.9285 从表中结果易看出在肋高为0.50m和0.49m的时候，肋基温度梯度为正值，传热量为负值，与实际情况不符。这是因为肋高增加，一定的耗材下，肋片直径变小，对流换热量处理成广义热源已不合适，即不能作为一维稳态导热模型看待。此时，增大肋片体积、增加导热率或减小对流换热系数，又能满足模型使用调节。忽略表中前两行的数据绘图见图5。图5 不同长度下肋片传热量曲线由图可知传热量随着肋高先增后减，传热量最大在肋高L=0.23m取得。因为=-Ay0 1，在肋基，温度始终为tw，即导热系数不变，肋基温度梯度y01为负且和圆柱截面积A随肋高增加而变小，所以存在最佳肋高使传热量最大。5.2 表面换热系数的影响材料的导热率=400(1+0.0035T)，肋表度黑度=0.80，圆柱肋片在不同表面换热系数h时，为使传热量最大，相应最佳肋高LGmax和最佳长径比CJB（肋高与半径的比值）见表3。表3 不同表面传热量下的LGmax和CJBh/(W/m2)45791012LGmax/m0.150.140.130.130.120.12CJB23.01901420.75591518.57226418.57226416.47106616.471066由表易知，随着表面传热系数的增加，最佳肋高是逐渐减小的。表面传热系数的增加，传热量增加，由=-Ay0 1知，需要增大肋基导热面积，所以最佳肋高减小。5.3 材料导热率的影响无辐射，肋表面空气的表面换热系数hc=8W/(m2)时，不同导热系数时，为使传热量最大，相应最佳肋高LGmax和最佳长径比CJB（肋高与半径的比值）见表4。表4 不同传热系数下的LGmax和CJB/(W/mk)100300400600700800LGmax/m0.180.280.310.360.390.41CJB30.2592858.70659268.38996385.58616796.504316104.022099由表易知，随着导热系数的增加，最佳肋高是逐渐增大的。因为导热系数变大，传热量增加，由=hcA(t-tf)知，需增加圆柱侧面积以加强换热，所以最佳肋高增加。6 结论在肋基，温度始终维持不变，即导热系数不变，肋基温度梯度为负且和圆柱截面积A随肋高增加而变小，由傅里叶公式可知存在最佳肋高使肋片传热量最大，在题目已知条件下，当肋高L=0.23m时取得最大散热量=15.0432W；表面传热系数的增加，传热量增加，由傅里叶公式知，需要增大肋基导热面积，所以最佳肋高减小；导热系数变大，传热量增加，由对流换热公式知，需增加圆柱侧面积以加强换热，所以最佳肋高增加。参考文献1 黄善波，刘中良.计算传热学基础.中国石油大学（华东）热能与动力工程系，20092 杨世铭，陶文铨.传热学（第四版）.高等教育出版社，2007附录1 主要程序表5 程序列表序号程序名称程序功能对应图表1全功能程序按题目要求，输出各个肋高下的传热量，求解最大传热量以及此时肋片的尺寸，输出温度和温度梯度的分布，验证程序图3、4、5，表1、2、3、4已知参数赋值 /输入圆柱肋体积V、任意给定肋高L及肋高变化步长bc V=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/dx f1=-Kd*y1*y1+Slx*(Bic*y0+Nr*(pow(y0+3.908667),4)-pow(3.868667,4); return;void pqw1(double Y,double *P,double *Q,double *W)/左边界处的第三类边界条件(x=xa)/P1*y1+Q1*y2=W1-用户应根据具体条件进行修改 *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+) /用肋高函数求肋高 LGk=LeiGao(L,bc,k); LGZ=LGk; x=xa; /利用试射法确定m shoot(N,NS,x,h,M1,M2,Eps,y0,y); /求热量 Qk=RL(V, y0); LGg=LeiGao(L,bc,g); LGZ=LGg; x=xa; /利用试射法确定m shoot(N,NS,x,h,M1,M2,Eps,y0,y); Qg=RL(V, y0); /输出任意长度及所对应的热量 printf(%f ,LGg); /输出长度 fprintf(fp,%f ,LGg); printf(%f ,y01); /输出肋基出温度梯度 fprintf(fp,%f ,y01); 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=%fn,r); fprintf(fp,最佳半径r=%fn,r); printf(最佳长径比CJB=%fn,CJB); fprintf(fp,最佳长径比CJB=%fn,CJB); x=xa; /利用试射法确定m shoot(N,NS,x,h,M1,M2,Eps,y0,y); Qmax=RL(V, y0);/求最大传热量 printf(最大热量Qmax=%6.4f ,Qmax); fprintf(fp,最大热量Qmax=%6.4f ,Qmax); printf(nnn); fprintf(fp,nnn); /输出最大传热量时的温度分分布 printf(输出最大传热量时的温度分布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=y0i; for(j=0;jNS;j+) rungek(N,&x,h,y); /根据求出的m解决问题 printf(%6.4f ,x); fprintf(fp,%6.4f ,x); for(i=0;iN;i+) printf(%9.6f ,yi); fprintf(fp,%9.6