散热器传热性能场协同仿真剖析.docVIP

散热器传热性能场协同仿真剖析 　　摘要管带式散热器外表面的空气流动与换热状况关系到发动机的动力性及燃油经济性。笔者根据场协同理论，采用数值模拟的方法对流经散热器翅片周围空气的速度场和温度场进行模拟计算，探讨了翅片的结构参数对流场的场协同角的影响。该方法从场协同角的角度为散热器翅片的结构优化设计提供了参考。 　　关键词散热器数值模拟场协同角优化设计 　　基金项目：衢州职业技术学院科研计划项目“建筑房顶集中式太阳能吸收制冷技术应用研究”（项目编号：QZYZ1408）。 　　作者简介：毛方（1984-），男，湖北咸宁人，衢州职业技术学院硕士研究生，研究方向：动力机械运动学仿真分析、动力机械传热分析。E-mail：maohuan2211@sinacom。 　　0引言 　　对于对流换热问题存在两个基本的矢量场，一个是流体流动时的速度场，另外一个是由于流体的温度不均匀时产生的温度梯度场，由于温度梯度具有方向性，因此温度梯度场也是一个矢量场。场协同原理由过增元教授及其合作者于1998年提出，其内容可以表述为：对流换热的性能除了和流体的速度和物性以及流体与固壁的温差有关外，还与流体速度场与流体热流场间协同的程度相关。在相同的速度和温度边界条件下，它们的协同程度愈好，换热强度愈高＼[1＼]。 　　1二维流场场协同理论 　　从二维平板层流边界层问题的能量守恒方程出发，将能量守恒方程和导热问题的能量守恒方程进行对比，可得到对流换热问题的场协同角度的物理机制＼[2＼]。 　　ρcpuTx+vTy=yλTy（1） 　　将式（1）在热边界层内积分得到 　　∫δt，x0ρcpuTx+vTy-=-λTy|w=qw（x）（2） 　　将式（2）左边的对流项改写成矢量点积形式U―?△T。 　　其中： 　　U―=（u，v） 　　△T=Tx+Ty 　　U―?△T=uTx+vTy=U―?△Tcosβ（3） 　　式中 β―速度矢量和温度梯度矢量的夹角，也叫场协同角。 　　从式（3）可以看出，速度矢量和温度梯度矢量的夹角，即场协同角对积分结果有直接关系，当该角度为0°时，积分值最大，该角度为90°时，积分值为零。角度越小，对传热越有利，反之亦然。 　　2翅片结构对散热器空气侧强化换热的场协同仿真分析 　　在管带式散热器表面布置百叶窗结构，一方面是为了增加散热表面积，另一方面是为了增强对空气的扰动，削弱层流底层的影响。扰动作用的结果是改变速度矢量和温度梯度矢量之间的协同程度，即改变整个或者局部流场区域内的场协同角＼[3＼]。使得场协同角度尽可能小，积分值尽可能大，获得最佳的对流换热系数。 　　图1所示为散热器翅片及其周围空气的有限元网格模型。流场计算方法为SIMPLE算法。SIMPLE算法是一种主要用于求解不可压流场的数值方法。 　　从图2和图3可以看出，在无百叶窗翅片部分，流体速度矢量一致，流场分布如图2左边部分所示。等温线分布均匀，其温度梯度方向基本垂直于翅片表面，未开百叶窗部分，翅片周围区域场协同角分布均匀，场协同角度在85°～91°之间，在离翅片较远的区域，场协同角稍小，在70°～85°之间。根据场协同理论，场协同角较大时不利于传热，因此，无百叶窗翅片部分传热效果不佳。 　　图1管带式散热器空气侧传热分析有限元模型 　　图2无百叶窗翅片部分的温度场与流场分布图 　　图3无百叶窗翅片部分翅片周围的场协同角分布 　　图4有百叶窗部分的温度场与流场分布图 　　图5有百叶窗部分翅片周围的场协同角分布特征 　　从图4、图5可知，在有百叶窗部分，百叶窗的存在加强了对流场的扰动作用，空气流经百叶窗时流动方向发生改变，沿百叶窗表面，由于边界层开始被间断的百叶窗破坏，温度场和流场分布也相应发生变化，速度矢量方向与温度梯度方向夹角在这一区域开始减小（90°），流场与温度场协同性得到改善，传热得到强化。 　　在开百叶窗部分，流场和温度场的协同角分布不均匀，在翅片周围以及离翅片较远的区域，场协同角分布均呈现明显的分布不均。但总的来说，场协同角小于90°的区域远大于未开百叶窗部分。这表明，百叶窗在改善空气和翅片之间的换热方面具有强化作用，而这一强化机制可以从改善速度场和温度场的协同性方面进行阐述。 　　3空气流速对散热器空气侧强化换热的场协同仿真分析 　　空气流速对换热的影响和雷诺数相关，在临界雷诺数前后，流场将发生层流到紊流的转变。从场协同的角度来看，实质是流场的场协同性由不一致向一致的转变。从整个空气流场来看，随着气流速度的不断增大，流场由于扰动增强发生显著变化，温度也随之发生改变。两者之间的协同性将空气的流场和温度场对于换热的影响统一起来。从翅片通道的局部区域来看，由于百叶窗的存在，空气流速的增大，百叶窗周围的流场和温度场的