粗糙度对大高宽比槽道内流动影响的数值模拟.pdfVIP

第30卷第5期 计 算 物 理 V01．30．NO．5 2013年9月 CHINESE OFCOMPUTATIONALPHYSICS JOURNAL Sep．，2013 文章编号：1001-246X(2013)05-0759-07 粗糙度对大高宽比槽道内流动影响的数值模拟 刘 东1， 蒋 斌1， 刘明侯2 (1．西南科技大学土木工程与建筑学院，绵阳621010；2．中国科学技术大学热科学和能源工程学院，合肥230027) 摘要：基于确定表面粗糙结构形状的PML模型模拟表面粗糙度对大高宽比小槽道内流动特性的影响，并与3．D 模拟结果对比，结果较为一致．基于该模型对不同的流速和不同的粗糙度进行模拟，结果表明：表面粗糙度导致槽 道内出现速度线性分布的流动底层．尺e数相同时，单位长度压降与相对粗糙度成二次方关系．粗糙单元高度相同 时，压降随Re线性增加． 关键词：小槽道；粗糙度；PML模型 中图分类号：TKl24 文献标识码：A 0 引言 随着集成电路封装技术的发展，单位面积上电子器件越来越多，从而导致严重的散热问题¨。．为了解 决这一问题，提出了微槽道散热器、微热管、喷雾冷却等新型冷却方式．而微槽道散热器由于其散热量大、 结构紧凑、低热阻、薄片型、紧凑轻便而受到研究者的广泛重视旧。。．针对微槽道散热器而言，从常规流动 的理论来看，相同流量的流体在流道中层流流动时，压降和槽道当量直径的平方成反比．因此，随着微槽 道尺寸的减小，要保持流量一定，压降会呈现四次方增大．为了减小流动阻力，作者‘41提出了大高宽比小 槽道的设计思路，但是在常规的大高宽比槽道中由于加工技术(线切割，机械加工)的限制，同样也会出现 相对粗糙度要比常规通道大很多的现象，一般其相对粗糙度可达10％一20％，因此很有必要对表面粗糙 度的影响进行研究． 随着计算技术的发展，越来越多的研究者开始使用计算机模拟技术来分析表面粗糙度对流动和换热的 影响．应用较多的方法主要有：确定几何形状法”“1，粗糙度涡旋粘性模型‘7|，多孔介质模型法‘8。91和表面粗 糙度正则摄动法¨¨12]．其中，多孔介质模型法是将流动区域根据不同的性质划分为中心区和边界粗糙层影 响区域两个区域，并采用常规的动量方程和引入阻力向量及经验的壁面渗透率的参数的方法对两个区域分 别处理．在这种方法中，边界粗糙层内阻力向量的确定将直接影响到计算结果的正确与否．现有研究往往通 过实验数据的拟合来确定．即使相同的实验数据，不同的研究者给出的拟合参数也不一致．本文基于确定的 表面粗糙度形状，采用流体力学阻力分析方法，提出确定边界层内阻力向量和经验的壁面渗透率参数的多孔 介质模型(PML)，对大高宽比小槽道内的流动特性进行研究，分析不同的m数情况下表面粗糙度对大高宽 比小槽道内流动特性的影响． 1 几何模型 在大高宽比微结构中，表面粗糙单元的形状非常复杂且各不相同．如果按照表面粗糙度的具体形状来直 接计算，网格将会变得非常庞大，现有计算机还很难满足计算的需求．研究中，通常将表面粗糙单元假设为几 种不同的规则形状来进行模拟”’”1．假设三维表面粗糙度的几何形状为圆柱体，并且周期性排列，粗糙单元 高度为s，两粗糙单元之间的距离为L，粗糙单元的直径为d(图1)．由于大高宽比微结构散热器的高宽比较 大，因此流动简化为三维平板流动，其中平板表面布满了高度为占，直径为d的粗糙单元． 收稿日期：2012—10—08；修回日期：2013一03一16 作者简介：刘东(1984一)，男，博士，讲师，主要从事微尺度流动和传热研究，E-mail：dtldl23@mail．ustc．edu．cn 万方数据 760 计 算 物 理 第30卷 在计算中，相对粗糙度占，定义为粗糙单元高度占与通道宽度 2H的比值，即 占，=8／(2H)． (1) 2基于确定表面形状的多孔介质模型【PML) 采用多孔介质模型模拟表面粗糙度是根据粗糙度对流动的影 响特点，将流动区域划分为流动