1、肋片散热数值计算干货分享肋片散热数值计算016年12月一、题目-3二、数值计算-4(1)网格划分-4(2)节点方程-(3)计算方式-(4)计算结果-6(5)温度分布云图-7(6)误差分析-10三、结论-10四、程序-11五、参考文献-15一、题目肋片优化问题考虑三种不同形状的肋片,如图所示。材料均为硬铝,热导率为,肋根半厚度为4m,肋高为25mm。对于梯形肋和圆弧边肋,最右端的平面部分半厚度为1m,且圆弧在最右端的切线为水平线。肋根温度即22,肋外流体温度即2,表面对流换热系数为.试编程求解每种肋片的温度分布及散热量,并讨论肋片形状对散热量、材料需求量的优劣。.文档交流 仅供参考.散热量 17
2、15 W 05 147.文档交流 仅供参考.截面面积 00 mm2 25 mm2 101mm.文档交流 仅供参考.qA 85。 124. 15.文档交流 仅供参考.肋顶温度 367 K(94) 355 K(82) 342K(9)二、数值计算鉴于肋片对称,因此只研究上半部分即可,肋半厚度处可以按绝热处理。(1)网格划分如图所示,用方格逼近曲面边界(三种形状的肋片都这样处理),即用图中蓝色网格边界来替代红色实际边界。网格在边界处的取法用以下规则:对于第i列节点,如果其中的第j个节点在红线上方或者恰好在红线上,并且第1个节点在红线下方或者正好在红线上,那么就取第j个节点作为第i列上的边界节点,最后把
3、所有边界节点连成锯齿状得到用方格逼近的边界(蓝色边界)。.文档交流 仅供参考.网格间距delta在运行程序时输入。(建议值0.mm,如果输入的值过大,那么计算结果不精确;输入值过小则程序需要很长时间才能算出结果).文档交流 仅供参考.(2)、节点方程内部节点肋根换热面传热量Q 按照肋根部的导热量计算,忽略根部y方向上的温度梯度(根据程序的运行结果来看,这种近似是合理的,对结果的影响很小),由求得的温度数据,在根部用(t(j,1)t(j,2))/eta作为温度梯度,(由于程序的原因,这里t(j,1)表示第列第个温度值)再用温度梯度乘上微元面积delta1(为计算方便,肋宽取为1m),然后把所有的
4、值求和,就得到从肋根部导出的热量.因为只研究上半部分,求出的散热量只有一半,再乘2即可。.文档交流 仅供参考.具体到程序中所使用的节点方程,一共有一下1类(3)、计算方式简单迭代法,矩形肋片允许误差取为1(-6),梯形肋片和圆弧形肋片的允许误差取为10(4)。()、计算结果肋片形式允许误差e网格间距散热量肋端温度截面面积AQ/A矩形肋片10(6)。1mm7394.22008.97梯形肋片1()。1mm15107。62519993圆弧形肋片10(4)0.1m586511014。12(5)、温度分布云图矩形肋片温度分布梯形肋片温度分布圆弧形肋片温度分布()、误差分析由于用矩形边界代替曲面边界,导致
5、截面周长增加,所以在程序中引入修正项L0,其中L为实际对流换热边界长度,L0为网格对流换热边界长度,将求得的热量乘上该修正项,能够减小误差。.文档交流 仅供参考.修正步骤已经写入程序。三、结论经过反复测试,发现对于矩形肋片,只要设置的允许误差e足够小,迭代求得的肋端温度值基本不受网格间距大小的影响,但是求得的散热量对网格间距比较敏感,当网格间距设为0。1mm,允许误差e=0.000001时,结果已经基本准确。.文档交流 仅供参考.对于变截面肋片,即梯形和圆弧形肋片,测试发现求得的散热量Q和肋端温度对e和间距delt都很敏感;发现设置同样的值,即e=0.01,da0.1m时,求得肋端温度57;d
6、ela0.05mm时,求得的肋端温度反而变小了,68.文档交流 仅供参考.但是可以预见,当e和e取的都足够小时,求得的散热量和肋端温度将会越来越接近,但是鉴于将e和deta同时取很小,程序将运行很长时间,所以不再尝试。.文档交流 仅供参考.对比计算结果可见,其实三种截面的肋片散热量差别不是很大,但是其QA值却差很大,也就是说矩形截面肋片用料比较多,而圆弧形、梯形肋片比较省材.但是也看到,虽然圆弧形肋片省材,但其肋端温度较低,也就是说肋端的散热温差较小,这其实也是一种材料利用不充分的体现.文档交流 仅供参考.本题目属于给定换热系数h和导热系数的情况,而我们知道,肋片效率fth(ml)m,m=,因
7、此肋片的效率取决于肋片的周长与面积比.四、程序矩形肋片程序juxingleipan。mfnctinjuxingleipinh=280;对流换热系数la=7;导热系数sprin(输入网格间距deta,单位为mm,建议值0。或0。(其中0.2运行时间较短).文档交流 仅供参考.dla=input(delta=);输入网格间距,单位mdeta=.00*delta;xnum=0.025/dela;%x轴划分数ynum=0.00/delta;%轴划分数tf0-273;流体温度t0=500-73;%肋根温度x=e(ynum+1,num+1);tzeros(ynm1,xum+1);wzos(num+1,nu
8、m+1);c=0;.文档交流 仅供参考.srintf(输入迭代允许误差,(若elt输入0。,则的建议值为10-6)当两次迭代之间的误差小于该数时,停止迭代).文档交流 仅供参考.e=iu(e=);%输入允许误差for=1:ynum+1;=1:xn+1;t(,)=10;end;任意假定一组初始温度值.文档交流 仅供参考.y=1;whil y=1 for =2:xnu;(1,j)(2t(2,)+(,j+)t(1,j-1)+2etf/lmd)(4+2*hdel/amda);nd;%节点方程.文档交流 仅供参考. foj2:xm;x(+1,j)(2*t(num,j)(ynu+1,j+1)+t(ynum
9、+1,1)+2delta*tf/lmda)/(4+2*delt/ad);e;节点方程.文档交流 仅供参考. for =1:ynum+;(,)=t;nd;节点方程 fori=2:num;x(i,m+)(*(i,xu)(i+1,xnm+)+t(,xnum+1)+2h*deltaflamda)/(42*hdelta/lama);end节点方程.文档交流 仅供参考. fr i=:ynum;j=:xnm;x(i,j)(i,1)+(i,j+1)t(i1,j)(+,j))/4;nd;%节点方程.文档交流 仅供参考. x(1,xnm+1)=(2*hdettflamat(1,xnum)+t(,xnum1)/(2
10、2*ea/lada);节点方程.文档交流 仅供参考. x(num+,xu+1)=(2*hdeltatflama+(nm,num+)+(ynum+1,xnm))/(22heltalamda);节点方程.文档交流 仅供参考. for=1:yum+1;j:num+1;w(i,)=abs(x(i,j)t(i,j));nd.文档交流 仅供参考. if (max(max(w)=e)判断两次迭代的误差是否小于允许值 =0; end =x;=c+1;nd for i:ynum+1;tid(i)(t(i,1)-t(,)delta;n;%求肋根部温度梯度.文档交流 仅供参考.Qlmdaum(td(1,:))dta
11、温度梯度与微元面积乘积求和,得到散热量 m=t(num/2+1,nu)a=linspace(0,25,num+1);=linspace(0,8,y+1);aa,bb=mshrd(,b);figremeh(a,,t);温度分布图fiureconturf(a,b,5);shading fland梯形肋片&圆弧形肋片程序 dderhaped。m注:本程序改变红色字的公式即可分别计算梯形肋片和圆弧形肋片functnladerhaped=2800;%对流换热系数r=4/6;%半径r=0;lad187;%导热系数t=30273;流体温度t0=5023;肋根温度spnf(输入网格间距delt,单位为mm)d
12、eta=ipt(deta=);%输入网格间距,单位mmsprinf(输入设定的误差值,当两次迭代的误差小于该值时停止迭代);uc=input(ucha=);输入设定的误差值dela=0001delta;xnum=005/delta;x轴划分数ynu0。008la;轴划分数tzrs(yum/2+,n+1);%设初始温度场为100zrs(ynum/21,num+1);%下面要确定第i列对应的行数i=1;hile (i=xnum+1)求出每一列的数据个数y()(3*(i1)*delt1000+10)/25;%梯形肋片;肋片形状可以任意更改,只需写出不同截面的方程即可.文档交流 仅供参考.(i)=( 3(i1)*ela1000 -