分-合式微通道流动传热性能数值模拟_蒋静智.pdf

第4 4卷第3期河 北 科 技 大 学 学 报V o l.4 4,N o.32 0 2 3年6月J o u r n a l o fH e b e iU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g yJ u n e2 0 2 3 文章编号:1 0 0 8-1 5 4 2(2 0 2 3)0 3-0 2 2 8-0 9分-合式微通道流动传热性能数值模拟蒋静智,刘华峰,安京州(河北科技大学机械工程学院,河北石家庄 0 5 0 0 1 8)摘 要:针对矩形微通道进出口压降大、温度分布不均匀,以及分形微通道受到分形维数和分支数限制适用范围较窄的问题,结合矩形微通道和分形微通道的优势设计一种分-合式微通道散热器。使用F l u e n t软件对散热过程进行数值模拟,研究微通道内分支倾斜角度变化对流动和传热性能的影响。结果表明,在1 0 0W/c m2的热流密度下,R e为9 7 0、分支倾斜角度为9 0 时,分-合式微通道平均温度降低了1 1.9K,最高温度降低了1 4.2K,N u增加了8 5.7%,整体传热性能(P E C)也最佳,达到1.4 4。分支的引入可以增加微通道内部换热面积,同时形成新的边界层,在分支内侧产生漩涡,有效提高了微通道散热器的传热性能,为微通道的优化设计提供了新的理论依据。关键词:传热学;微通道;微型设备;数值模拟;传热中图分类号:T K 1 7 2 文献标识码:A D O I:1 0.7 5 3 5/h b k d.2 0 2 3 y x 0 3 0 0 3 收稿日期:2 0 2 3-0 4-2 7;修回日期:2 0 2 3-0 5-0 5;责任编辑:冯 民基金项目:国家自然科学基金(5 1 7 0 6 0 5 8);河北省自然科学基金(E 2 0 1 9 2 0 8 3 4 5)第一作者简介:蒋静智(1 9 7 2),女,河北阜城人,教授,博士,主要从事强化传热及储能方面的研究。E-m a i l:j j z h i 2 0 0 01 2 6.c o m蒋静智,刘华峰,安京州.分-合式微通道流动传热性能数值模拟J.河北科技大学学报,2 0 2 3,4 4(3):2 2 8-2 3 6.J I AN GJ i n g z h i,L I U H u a f e n g,A nJ i n g z h o u.N u m e r i c a ls i m u l a t i o no f f l o wa n dh e a tt r a n s f e rp e r f o r m a n c ei nb r a n c ha n dc o m b i n e d-b r a n c hm i c r o c h a n n e lJ.J o u r n a l o fH e b e iU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y,2 0 2 3,4 4(3):2 2 8-2 3 6.N u m e r i c a l s i m u l a t i o no f f l o wa n dh e a t t r a n s f e rp e r f o r m a n c e i nb r a n c ha n dc o m b i n e d-b r a n c hm i c r o c h a n n e lJ I ANGJ i n g z h i,L I U H u a f e n g,ANJ i n g z h o u(S c h o o l o fM e c h a n i c a lE n g i n e e r i n g,H e b e iU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y,S h i j i a z h u a n g,H e b e i 0 5 0 0 1 8,C h i n a)A b s t r a c t:T h ep r e s s u r ed r o pa t t h e i n l e t a n do u t l e to f r e c t a n g u l a rm i c r o c h a n n e l i s l a r g e,a n dt h et e m p e r a t u r ed i s t r i b u t i o ni su n e v e n;F r a c t a lm i c r o c h a n n e l i s l i m i t e db y f r a c t a l d i m e n s i o na n dt h en u m b e ro fb r a n c h e s,s o i t sa p p l i c a t i o ns c o p e i sn a r r o w.Ab r a n c ha n dc o m b i n e d-b r a n c hm i c r o c h a n n e l h e a t s i n kw a sd e s i g n e db yc o m b i n i n g t h e a d v a n t a g e so f r e c t a n g u l a rm i c r o c h a n n e l a n df r a c t a lm i c r o c h a n n e l,a n dF l u e n t s o f t w a r ew a su s e dt on u m e r i c a l l ys i m u l a t e i t sh e a td i s s i p a t i o np r o c e s st os t u d yt h ef l o wa n dh e a t t r a n s f e rp e r f o r m a n c ew h e nt h eb r a n c ht i l t a n g l e i n s i d e t h em i c r o c h a n n e l c h a n g e s.T h e r e s u l t s s h o wt h a t a t t h eh e a t f l o wo f1 0 0W/c m2,w h e nt h eR ei s9 7 0a n dt h eb r a n c ht i l t a n g l e i s9 0,t h ea v e r a g et e m p e r a t u r eo f t h eb r a n c ha n dc o m b i n e d-b r a n c hm i c r o c h a n n e l d e c r e a s e sb y1 1.9K,t h em a x i m u mt e m p e r a t u r ed e c r e a s e sb y1 4.2K,t h eN ui n c r e a s e sb y8 5.7%,a n dt h ep e r f o r m a n c ee v a l u a t i o nc r i t e r i o n(P E C)i sa l s ot h eb e s tb yr e a c h i n g1.4 4.T h e i n t r o d u c t i o no fb r a n c h e sc a ni n c r e a s et h eh e a tt r a n s f e ra r e ai n s i d et h e m i c r o c h a n n e l,f o r man e wb o u n d a r yl a y e r,a n dg e n e r a t ev o r t i c e sa tt h ei n n e rs i d eo ft h eb r a n c h,第3期蒋静智,等:分-合式微通道流动传热性能数值模拟e f f e c t i v e l y i m p r o v i n gt h eh e a tt r a n s f e rp e r f o r m a n c eo ft h em i c r o c h a n n e lh e a ts i n k s,w h i c hp r o v i d e sn e wt h e o r e t i c a lb a s i sf o ro p t i m i z i n gt h ed e s i g no fm i c r o c h a n n e l s.K e y w o r d s:h e a t t r a n s f e r;m i c r o c h a n n e l;m i c r od e v i c e;n u m e r i c a l s i m u l a t i o n;h e a t t r a n s f e r 随着电子工业技术的迅速崛起,电子元器件正在向集成化、高频化方向发展,电路板热流密度也在不断突破新高。当电子元器件工作时,如果热量无法被及时散出,内部温度就会逐渐上升,从而严重影响自身的性能和可靠性1-2。研究表明:电子元器件的可靠性对温度十分敏感,当工作温度达到7 08 0后,每上升1,可靠性就会下降5%,温度过高已成为现在大多数设备失效的主要原因3。微通道散热器具有体积小、散热快、反应时间短等优势,在解决微型设备器件的高热负荷问题方面极具应用前景4。近年来,矩形微通道5-8和分形树状微通道9-1 1的结构优化研究领域出现了很多优秀研究成果,为微通道的设计提供了良好的理论指导。TU C K E RMAN等1 2首次提出了微通道散热器,发现微尺度下流体的散热效果要比宏观尺度好得多。陈孝根等1 3研究了不同高宽比矩形微通道的散热性能,结果表明,微通道的高宽比为24时最好。GUNNA S E GA R AN等1 4对不同截面形状微通道的传热性能进行了实验验证,结果表明,截面形状对微通道的传热性能有较大的影响。朱崎峰等1 5对不同凹槽形状微通道进行实验验证,结果表明,微通道凹槽为扇形且水力直径较小时换热系数高,凹槽为三角形且水力直径较大时换热系数高。CHA I等1 6在传统平直微通道结构的侧壁上添加扇形肋,结果表明,这种结构的微通道散热器综合性能有了很大提升。P E NG等1 7验证了分形树状微通道的流动性能,结果表明,树状结构可以有效降低流动阻力。HUANG等1 8对变截面分形树状微通道散热器的传热性能进行研究,发现变截面结构可以进一步提高散热器的传热性能。WANG等1 9使用3 D打印技术制造了Y型和H型树状结构散热器,发现H型结构散热器的传热性能较好,Y型结构散热器的压降较小。P E N C E2 0分析了分形树状微通道散热器转角结构和交叉结构对流动和传热性能的影响,结果表明,在转角结构和交叉结构区域,流体流动扰动增强并破坏了流动边界层,从而提高微通散热器的换热能力。矩形微通道散热器结构的改进可以提高传热性能,但是进出口压降也在增加,并且沿流动方向流体的温度会逐渐上升。分形树状微通道散热器在具有更高换热性能和较为均匀的温度分布的同时,具有较低的压降和泵功消耗,但是分形树状微通道散热器受到分形维数和分支数的限制,其入口段换热面积小,换热效率较低,且结构形状被限制,导致其适用性下降。针对上述问题,设计一种分-合式微通道散热器2 1,并对微通道内的流体流动和传热性能进行模拟研究。1 微通道数值模拟HW图1 矩形微通道结构示意图F i g.1 S c h e m a t i cd i a g r a mo f r e c t a n g u l a rm i c r o c h a n n e l s t r u c t u r e1.1 物理模型图1为矩形微通道散热器示意图。该散热器由4 0根截面为矩形的微通道组成,热源基板的尺寸为2 0mm 2 0mm,热流密度为1 0 0 W/c m2。微通道的流道高度H为4 0 0m,宽W为3 0 0m,微通道上下壁厚1均为5 0m,微通道之间的壁厚2为2 0 0m。图2为分-合式微通道散热器示意图。微通道内分支倾斜角度为;分支水力直径表现为宽度变化,倾斜分支的宽度相等为W1,通道分流后宽度减小为W2,合流后宽度增加为W3,故微通道壁厚略有变化;微组分支与合支的距离为h。引入矩形微通道是为了以矩形微通道作为基准微通道散热器,与分-合式微通道散热器进行对比,由于微通道散热器中微通道的排列具有周期性,其传热性能也呈现周期性分布,为了减少模拟仿真的计算时间和计算成本,以及能够完整反映分-合式微通道结构的周期性,选取5根微通道作为仿真模型。MUR R AY2 2通过研究 血 管 的 流 阻 特 性,得 出 父 血 管 半 径 的 立 方 和 等 于 子 血 管 半 径 的 立 方 总922河 北 科 技 大 学 学 报2 0 2 3年hW2W3W1图2 分-合式微通道结构示意图F i g.2 S c h e m a t i cd i a g r a mo fb r a n c ha n dc o m b i n e d-b r a n c hm i c r o-c h a n n e l s t r u c t u r e和,即:D3=D31+D32+,(1)其中,D为血管的水力直径。当微通道的分支宽度W1和分流后主流道宽度W2相等时,根据式(1)计算可得各流道的宽度,分支和合支均匀排布,计算可得W1=1 7 0mm,W2=1 7 0mm,W3=3 6 0mm,h=1 1 5mm。以分支的倾斜角度为变量,设计了5种微通道结构,分支倾斜角度编号及角度如表1所示。表1 微通道结构参数表T a b.1 T a b l eo fm i c r o s t r u c t u r es t r u c t u r ep a r a m e t e r编号1#2#3#4#5#6#/()3 04 56 07 59 01.2 数学模型及控制方程为简化模型,在模拟过程中作如下假设:1)微通道散热器的传热和流体在流道中的流动为稳态过程,且流体流动过程为充分发展的层流流动(微尺度下R e小于10 0 0时,流动为层流);2)假设不受热辐射和重力的影响;3)流体流动为充分发展的层流流动,故不考虑流道内轴向传导和黏性耗散;4)微通道内固体材料为硅,其物性参数如表2所示,流体材料为去离子水。表2 硅的热物性参数表T a b.2 T a b l eo f t h e r m o p h y s i c a l p a r a m e t e r so f s i l i c o n材料密度/(k gm-3)比热容/(Jk g-1K-1)导热系数/(Wm-1K-1)硅23 2 8.37 0 01 4 8根据上述假设并依据传热学2 3,本模型采用的控制方程如下。连续性方程:xi(ui)=0。(2)动量方程:xi(fuiuj)=-pxj+xifujxi+uixj。(3)能量方程:xi(fuicpfT)=xifTxi+f2uixi2+(ujxi+uixj)2。(4)对于微通道的固体区域,能量方程为xi(sTxi)=0,(5)式中:i,j=1,2,3;下标f表示流体,s表示固体;为密度,k g/m3;u为速度,m/s;p为压力,P a;为流体动力黏度,P as;cp为定压比热容,J/(k gK);T为温度,K;为导热率,W/(mK)。模型的边界条件是:散热器入口为速度入口边界条件,流体入口温度为2 9 3K;出口设置为压力出口边界条件,出口面压力为1.0 1 31 05P a,即表示压力为0;流体与固体接触面为耦合边界条件,无速度滑移;散热器底面为恒定热流密度边界,1 0 0W/c m2;其余面为绝热边界条件。流动和传热的评价指标如下。通过压降的计算结果以及模拟结果可以求得流体的摩擦阻力系数f,摩擦阻力系数的计算公式如下:f=2 pu2mDcL,(6)式中:Dc为通道当量直径,m;um为平均流速,m/s;L为通道的长度,m;p为通道压降,P a。对流换热系数反映了流体与固体换热能力的强弱,平均对流传热系数h的计算公式为032第3期蒋静智,等:分-合式微通道流动传热性能数值模拟h=q AAb(Tw-Tf,a v e),(7)式中:A为微通道基底面面积,m2;Ab为微通道固液交界面面积,m2;Tf,a v e为流体平均温度,K。N u表征了微通道的换热能力,N u计算公式为N u=h Dcf。(8)为 综 合 考 虑 微 通 道 散 热 器 的 传 热 特 性 与 摩 擦 阻 力,采 用 整 体 传 热 性 能(p e r f o r m a n c ee v a l u a t i o nc r i t e r i o n,P E C)来验证微通道散热器的综合能力,其计算公式为P E C=N u/N u0(f/f0)1/3,(9)式中:N u0和f0为基准微通道散热器计算值;N u和f为分-合式微通道计算值。1.3 网格设置采用基于有限元方法进行模拟计算的F l u e n t软件进行模拟仿真,采用S I MP L E C算法用来求解连续性、动量和能量方程,当能量方程的残差值小于1 0-8,其他方程的残差值小于1 0-5时,模拟收敛,计算完成。有限元方法的核心是通过网格划分的方式将复杂的几何实体离散成有限个规则的几何单元,网格的质量以及数量直接决定结果的准确性。使用微通道分支倾斜角度为4 5 时的模型进行网格无关性验证,采用4种不同网格数目进行对比,验证结果如表3所示。网格数目最多时的值F0和其他网格数目所得值F之间的相对误差e定义为e=|F0-FF0|1 0 0%。(1 0)表3 微通道网格无关性验证T a b.3 V e r i f i c a t i o no fm i c r o c h a n n e l g r i d i n d e p e n d e n c e网格数目N ueN u/%fef/%ReR/%18 8 88 2 01 8.3 81.6 20.3 2 843.2 91.0 3 413.0 429 3 04 7 61 8.9 31.3 20.3 3 670.8 61.0 1 310.4 353 3 46 0 81 8.7 60.4 00.3 3 920.1 11.0 0 440.0 866 1 52 8 71 8.6 80.3 3 961.0 0 36如表3所示,随着网格数目的不断增加,计算结果之间的误差也在不断减小,当网格数目为53 3 46 0 8时,N u、摩擦阻力系数和热阻的误差分别为0.4 0%,0.1 1%和0.0 8%。为了兼顾模拟计算的准确性与计算机性能,采用网格数目为53 3 46 0 8时的网格划分方式进行后续的模拟计算。1.4 数值计算结果有效性验证S T E I NK E等2 4对微通道中的单相流体摩擦系数进行研究,微通道截面为矩形、流体流动为层流时的对图3 微通道散热器摩擦系数对比曲线F i g.3 C o m p a r i s o nc u r v eo f f r i c t i o nc o e f f i c i e n to fm i c r o c h a n n e lh e a t s i n k流传热摩擦系数的理论计算公式为fR e=9 6(1-1.3 5 53c+1.9 4 672c-1.7 0 123c+0.9 5 644c-0.2 5 375c),(1 1)其中,c为通道的宽高比。图3为微通道散热器摩擦阻力系数对比曲线。从图3可以看出,当流体的物性参数为常数时,摩擦阻力系数的模拟值与理论值相差较大;当流体为变物性时,其物性参数根据参考文献2 3 选取,从图3可以看出,流体物性参数随温度变化时,摩擦系数的模拟值与理论值很接近,最大差值为2.8%,因此微通道内的流体采用变物性。由于变物性模拟结果与理论公式计算结果相差很小,从而验证了模拟结果的可靠性。2 结果与讨论132河 北 科 技 大 学 学 报2 0 2 3年对6种不同模型进行模拟计算后,从流动特性、压力特性以及温度特性3方面对不同结构微通道散热器进行分析,以验证微组的引入对传热性能的影响。2.1 流动特性当R e=3 2 2时,微通道x-y截面的速度云图如图4所示(z=0.2 5mm,0mmx2.5mm,7.5mmy1 2.5mm)。图4 不同微通道x-y截面的速度云图F i g.4 V e l o c i t yc l o u dm a p so fd i f f e r e n tm i c r o c h a n n e lx-yc r o s s-s e c t i o n s从图4可以看出,微通道贴近壁面处流体的流速几乎为0,且越靠近流道中心,其流速越大。矩形微通道由于流道形状没有变化,流动状态一致,对流速影响较小。对于分-合式微通道散热器,当流体没有流经分支时,流体在微通道的流动状态与矩形微通道内的流动状态相似;当流体流经分支时,由于局部流道变宽,流体的流速有明显的下降趋势,且在分支处有新的边界层形成;当流体流经合支时,由于局部流道变窄,流体的流速有明显的上升趋势,在合支处也有新的边界层形成。流体每流经一处分支或合支时,其流动状态的变化相似。同时从图4可以看出,分支倾斜角度的变化对微通道分支处的流体流动有着较大的影响,对主流道区流体的流动影响相对较小。随着分支倾斜角度的不断增加,流道内的最大流速在不断变大,并且随着分支倾斜角度的不断增加,分支处的流速也出现小幅的上升,并且都形成新的边界层。当分支倾斜角度为9 0 时,截面处的最大流速与其他角度时的最大流速相比最大,为2.2 4m/s。当R e=3 2 2时,微通道x-y截面的速度矢量图如图5所示(z=0.2 5mm,0.2 5mmx0.7 5mm,1 2.9mmy1 3.6 5mm)。从图5可以看出,矩形微通道的流体速度呈线性分布,几乎无扰动。分-合式微通道的流体在分支内侧产生了漩涡,并且随着分支倾斜角度的增加,漩涡的大小明显增加,流体的扰动程度不断增加,这主要是因为流体流动存在惯性,当流经流道弯曲段时流体不能突然变向,导致流体在转弯内壁面开始产生分离现象,弯度越大,分离现象越明显,甚至出现与主流区流速相反的回流区域,即漩涡。当分支倾斜角度为9 0 时,微通道分支处的漩涡最大。2.2 压力特性不同R e下,微通道总压降随分支倾斜角度的变化曲线如图6所示。从图6可以看出,随着R e的增加,不同结构分-合式微通道的总压降均在不断增加,且增加幅度也随着分支倾斜角度的增加而不断变大。在不同R e下,分-合式微通道的总压降均高于矩形微通道的总压降,这主要是因为在分支处产生了漩涡,引起流体的碰撞与摩擦,使分支处流体的局部阻力损失增加,从而增加了流体流动的能量损失。并且微通道的总压降也随着分支倾斜角度的增加而增加,这是由于随着分支倾斜角度的增加,分支处的漩涡不断变大,导致232第3期蒋静智,等:分-合式微通道流动传热性能数值模拟图5 不同微通道x-y截面的速度矢量图F i g.5 V e l o c i t yv e c t o rp l o t so fd i f f e r e n tm i c r o c h a n n e lx-yc r o s s-s e c t i o n s局部阻力损失增加。当R e为9 7 0、分支倾斜角度为9 0 时,微通道的总压降最大。不同R e下,微通道摩擦阻力系数随分支倾斜角度的变化曲线如图7所示。从图7可以看出,在不同R e下,分-合式微通道的摩擦阻力系数都大于矩形微通道摩擦阻力系数。随着分支倾斜角度的增加,微通道的摩擦阻力系数也在增加,这主要是因为在微通道内分支的存在产生了新的壁面,增加了微通道的流通面积,并且分支处产生漩涡,能量损失增加,使微通道的总压降增加,摩擦阻力系数变大,因此随着分支倾斜角度的增加,摩擦阻力系数出现逐渐上升的趋势。图6 微通道总压降随分支倾斜角度的变化曲线F i g.6 C u r v eo f t o t a lp r e s s u r ed r o p i nm i c r o c h a n n e la sa f u n c t i o no fb r a n c ht i l t a n g l e图7 微通道摩擦阻力系数随分支倾斜角度的变化曲线F i g.7 V a r i a t i o nc u r v eo fm i c r o c h a n n e l f r i c t i o nr e s i s t a n c ec o e f f i c i e n tw i t hb r a n c ht i l t a n g l e2.3 温度特性当R e=3 2 2时,微通道x-y截面的温度云图如图8所示(z=0.2 5mm,0mmx2.5mm,7.5mmy1 2.5mm)。从图8可以看出,微通道内流体沿流动方向,其温度分布都呈现出中心处温度低、壁面温度高的趋势。矩形微通道内由于流速分布一致,沿着流动方向,微通道内的温度在不断上升。分-合式微通道内,分支倾斜角度的变化对微通道内流体温度分布有较大的影响,当分支倾斜角度较小时,分支处的漩涡较332河 北 科 技 大 学 学 报2 0 2 3年图8 不同微通道x-y截面的温度云图F i g.8 T e m p e r a t u r ec l o u dm a p so fd i f f e r e n tm i c r o c h a n n e lx-yc r o s s-s e c t i o n s小,流体扰动程度小,因此分支处的局部换热系数较小,分支处的流体温度较低;当分支倾斜角度较大时,分支处的漩涡变大,流体扰动程度增加,局部换热系数增加,分支处的流体温度升高。因此,随着分支倾斜角度的增加,微通道内的温度分布更加均匀,高温区域也在不断减少。不同R e下,微通道平均温度和最高温度随分支倾斜角度的变化曲线分别如图9和图1 0所示(横坐标为R处代表矩形微通道散热器,与下文情况一致)。图9 微通道平均温度随分支倾斜角度的变化曲线F i g.9 V a r i a t i o nc u r v eo f t h ea v e r a g e t e m p e r a t u r eo fm i c r o c h a n n e lw i t hb r a n c ht i l t a n g l e图1 0 微通道最高温度随分支倾斜角度的变化曲线F i g.1 0 V a r i a t i o nc u r v eo f t h em a x i m u mt e m p e r a t u r eo fm i c r o c h a n n e lw i t hb r a n c ht i l t a n g l e 从图9可以看出,随着R e的增加,各个结构微通道的平均温度均不断降低,但是其变化幅度在不断减小,这说明增加流体的流速可以提高微通道的换热能力,但当流速增加到一定值后,换热能力的提升效果趋于平缓。在不同R e下,分-合式微通道的平均温度均低于矩形微通道的平均温度,这主要是因为分支处的流体扰动程度大,微通道换热能力提升。随着分支倾斜角度的不断增加,微通道内的平均温度在不断下降,这是由于随着分支倾斜角度的增加,分支处的流体扰动程度在不断增加。当R e为9 7 0、分支倾斜角度为9 0 时,微通道的平均温度最低,为3 0 8.9K,与矩形微通道相比降低了1 1.9K。从图1 0可以看出,在不同R e下,分-合式微通道的最高温度均低于矩形微通道的最高温度,且随着分支倾斜角度的增加,微通道内的最高温度在不断下降。当R e为9 7 0、分支倾斜角度为9 0 时,微通道的最高432第3期蒋静智,等:分-合式微通道流动传热性能数值模拟温度最低,为3 1 6.6K,与矩形微通道相比降低了1 4.2K。不同R e下,微通道N u随分支倾斜角度的变化曲线如图1 1所示。从图1 1可以看出,随着R e的增加,各个结构微通道的N u均不断增加,但其增加幅度却越来越小。在不同R e下,分-合式微通道的N u均高于矩形微通道的N u,并且N u也随着分支倾斜角度的增加而不断增加。当R e为9 7 0、分支倾斜角度为9 0 时,N u取得最大值,表示其具有最佳的传热性能,与矩形微通道相比,在该条件下,N u最大,为3 5.8,增加幅度最大,为8 5.7%。不同R e下,微通道P E C随分支倾斜角度的变化曲线如图1 2所示。从图1 2可以看出,在不同R e下,分-合式微通道的P E C都随着分支倾斜角度的增加而增加,当分支倾斜角度小于4 5 时,增加幅度较大,大于4 5 时增加幅度较小。分支倾斜角度为4 5,R e为6 4 6时,数据变化趋势明显不同,分析认为,当分支倾斜角度为4 5 时,R e由4 8 4增加到6 4 6与R e由6 4 6依次增加到9 7 0时相比,摩擦阻力系数的下降趋势明显要大(见图7),而N u下降趋势相差不大(见图1 1)。P E C与摩擦阻力系数成反比,所以当分支倾斜角度为4 5,R e为6 4 6时,P E C增加幅度要大。在本文模拟计算范围内,当R e为9 7 0、分支倾斜角度为9 0 时,P E C取得最大值,为1.4 4,表示其具有最佳的综合传热性能。图1 1 微通道N u随分支倾斜角度的变化曲线F i g.1 1 C u r v eo fm i c r o c h a n n e lN uc h a n g i n gw i t hb r a n c ht i l t a n g l e图1 2 微通道P E C值随分支倾斜角度的变化曲线F i g.1 2 C u r v eo fm i c r o c h a n n e lP E Cv a l u e sv a r i a t i o nw i t hb r a n c ht i l t a n g l e3 结 论通过对分-合式微通道散热器的散热过程进行模拟计算,研究了分支倾斜角度变化对微通道内流体流动和传热性能的影响,主要结论如下。1)随着R e的增加,矩形微通道与分-合式微通道的总压降不断增加,但微通道的换热能力都得到提升,且提升的幅度越来越小。2)分支结构打破了微通道的流动边界层,并且流体在分支处产生分离现象,形成漩涡,增加流体扰动程度,冷热流体混合更加充分,高温区域不断减少,微通道内的温度分布更加均匀。3)随着分支倾斜角度的增加,微通道的流速不断增加,分支处的流体漩涡不断增大,微通道的总压降和摩擦阻力系数都呈现不断增加的趋势;最高温度也在不断下降,微通道的N u和P E C也在不断增加。当分-合式微通道的分支倾斜角度为9 0 时,微通道的平均温度为3 0 8.9K,降低了1 1.9K;最高温度为3 1 6.6K,降低了1 4.2K;N u为3 5.8,增加了8 5.7%;P E C最佳,为1.4 4。研究发现,分支的引入可以明显提高微通道散热器的传热性能,并且分支处流体扰动程度增加,改善了微通道的温度分布,但是沿流动方向,流体的温度逐渐升高,换热能力减少,微通道后端温度比前端温度高出很多。为解决这一问题,后续研究可从改进微通道内分支的排布方式出发,设计更合理的结构,以改善微通道散热器温度分布不均的问题,为微通道的发展提供新的思路和方法。参考文献/R e f e r e n c e s:1 宁璐璐,刘清江,杨凤,等.电子器件冷却技术研究进展J.制冷与空调,2 0 2 1,2 1(1 2):1-7.532河 北 科 技 大 学 学 报2 0 2 3年N I N GL u l u,L I U Q i n g j i a n g,YAN GF e n g,e t a l.R e s e a r c hp r o g r e s so fe l e c t r o n i cd e v i c ec o o l i n gt e c h n o l o g yJ.R e f r i g e r a t i o na n dA i r-c o n d i-t i o n i n g,2 0 2 1,2 1(1 2):1-7.2 孙亚松,刘红敏.微通道强化换热研究进展J.应用化工,2 0 2 2,5 1(4):1 1 1 6-1 1 1 8.S UNY a s o n g,L I U H o n g m i n.A d v a n c e s i nm i c r o c h a n n e l e n h a n c e dh e a t t r a n s f e rJ.A p p l i e dC h e m i c a l I n d u s t r y,2 0 2 2,5 1(4):1 1 1 6-1 1 1 8.3 P A T I LNG,HO T T ATK.Ar e v i e wo nc o o l i n go f d i s c r e t eh e a t e dm o d u l e su s i n g l i q u i d j e t i m p i n g e m e n tJ.F r o n t i e r s i nH e a t a n dM a s sT r a n s f e r,2 0 1 8,1 1.D O I:1 0.5 0 9 8/h m t.1 1.1 6.4 赵文博,应展烽,宋新阳,等.基于热约束的功率器件散热器成本最小化方法J.科学技术与工程,2 0 2 1,2 1(3 0):1 2 9 4 6-1 2 9 5 1.Z HAO W e n b o,Y I N GZ h a n f e n g,S ON G X i n y a n g,e ta l.C o s t-m i n i m i z a t i o nm e t h o df o rh e a t s i n ko fp o w e rd e v i c e sb a s e do nt h e r m a lc o n-s t r a i n t sJ.S c i e n c eT e c h n o l o g ya n dE n g i n e e r i n g,2 0 2 1,2 1(3 0):1 2 9 4 6-1 2 9 5 1.5 D A T T AA,S HA RMAV,S ANYA LD,e t a l.Ac o n j u g a t eh e a t t r a n s f e r a n a l u s i s o f p e r f o r m a n c e f o r r e c t a n g u l a rm i c r o c h a n n e lw i t h t r a p r z o i-d a l c a v i t i e sa n dr i b sJ.I n t e r n a t i o n a l J o u r n a l o fT h e r m a lS c i e n c e s,2 0 1 9,1 3 8:4 2 5-4 4 6.6 肖润锋,侯予,吕坤鹏,等.大深宽比微通道热沉流动换热特性数值模拟及优化J.低温工程,2 0 1 8(6):1-7.X I AOR u n f e n g,HOUY u,L YUK u n p e n g,e t a l.N u m e r i c a l s i m u l a t i o na n do p t i m i z a t i o no f f l o wa n dh e a t t r a n s f e r c h a r a c t e r i s t i c s o f l a r g e-a s-