丁杰，唐玉兔

(南車株洲電力機(jī)車研究所有限公司南車電氣技術(shù)與材料工程研究院，湖南株洲 412001)

為化解強(qiáng)化換熱與流動(dòng)阻力的矛盾，優(yōu)化水冷散熱器的設(shè)計(jì)，人們?cè)诹鞯澜Y(jié)構(gòu)、形式、微通道等方面做出了許多努力［1-11］。其中，文獻(xiàn) ［11］就水冷散熱器入口總流量60 L/min時(shí)帶小突起棗核狀、光滑棗核狀、菱形、圓形和方形5種翅柱結(jié)構(gòu)進(jìn)行了對(duì)比研究，得到帶小突起棗核狀的換熱效果最佳、流動(dòng)阻力較大的結(jié)論。由于該翅柱式IGBT水冷散熱器分為3個(gè)并聯(lián)的支路，采用泄壓槽的方式進(jìn)行均流，其均流效果和整個(gè)散熱器的散熱效果有待研究。

為此，文中利用FLUENT軟件對(duì)泄壓槽對(duì)流量分配的影響進(jìn)行研究，對(duì)水冷散熱器的散熱能力進(jìn)行數(shù)值模擬。研究結(jié)果可為水冷散熱器的設(shè)計(jì)工作提供指導(dǎo)。

1 水冷散熱器的結(jié)構(gòu)

圖1是某翅柱式IGBT水冷散熱器結(jié)構(gòu)示意圖，其內(nèi)部槽道采用了并聯(lián)型和串聯(lián)型相結(jié)合的方式。

圖1 水冷散熱器結(jié)構(gòu)示意圖

冷卻介質(zhì)首先從入口流入，均流給3個(gè)支路，每個(gè)支路依次流經(jīng)2個(gè)IGBT元件下方的槽道帶走熱量，然后匯流到出口流出。匯流到出口的槽道設(shè)置了一段細(xì)小彎曲的泄壓槽，其作用是平衡3個(gè)支路的壓降，保證每個(gè)支路的流量相同。每個(gè)支路采用帶小突起棗核狀翅柱結(jié)構(gòu)進(jìn)行叉排布置，以提高換熱效果。

2 仿真策略與仿真建模

由于翅柱的細(xì)節(jié)特征尺寸很小，小突起可以對(duì)流動(dòng)和傳熱產(chǎn)生很大的影響，幾何結(jié)構(gòu)特征限制了網(wǎng)格的尺寸，不能對(duì)整個(gè)水冷散熱器進(jìn)行仿真，否則會(huì)大大超出計(jì)算資源的限制。鑒于此，可以先對(duì)槽道內(nèi)的流體區(qū)域 (見圖2，分為有泄壓槽和無泄壓槽2個(gè)模型，用于分析泄壓槽對(duì)流量分配的影響)劃分0.2 mm尺寸的網(wǎng)格，通過單獨(dú)計(jì)算流場(chǎng) (求解質(zhì)量守恒方程和動(dòng)量守恒方程)得到流向各支路的流量，分析泄壓槽對(duì)3個(gè)并聯(lián)支路均流的影響。然后參考文獻(xiàn)［11］的方法建立如圖3所示的對(duì)稱熱仿真模型，利用前面計(jì)算出的流量作為輸入條件，進(jìn)行流場(chǎng)和溫度場(chǎng)的耦合計(jì)算 (求解質(zhì)量守恒方程、動(dòng)量守恒方程和能量守恒方程)，從而獲得準(zhǔn)確的溫度場(chǎng)分布。

圖2 槽道流動(dòng)仿真模型

圖3 熱仿真模型

冷卻介質(zhì)使用50℃的40%乙二醇+60%水，入口流量為30 L/min。通過計(jì)算可知水冷散熱器入口的Re數(shù)為21 010，大于10 000，屬于湍流狀態(tài)。假定3個(gè)支路入口處的流量相同，計(jì)算可知支路入口處的平均流速為0.527 m/s，Re數(shù)為1 479，該值小于2 300?？紤]到翅柱使流道截面積減小，流速會(huì)增加，且小突起會(huì)起到強(qiáng)化擾流的作用，使得翅柱之間的流動(dòng)狀態(tài)難以界定，因此選擇層流和標(biāo)準(zhǔn)k-ε模型2種最具代表性的仿真模型來分析水冷散熱器內(nèi)部的流動(dòng)狀態(tài)［12］。

3 槽道流動(dòng)情況的仿真與分析

圖4是30 L/min入口流量時(shí)，有、無泄壓槽兩種情況分別采用標(biāo)準(zhǔn)k-ε模型計(jì)算得到的流速跡線圖，從圖中可以大致看出冷卻介質(zhì)在槽道內(nèi)的流動(dòng)情況，最大流速主要出現(xiàn)在轉(zhuǎn)折區(qū)域。

圖4 流速跡線圖

有泄壓槽的最大流速為2.02 m/s，比無泄壓槽的最大流速略小。利用FLUENT軟件的Report功能，可以得到3個(gè)支路的平均流速，再通過平均流速與截面積的乘積得到平均流量。有泄壓槽時(shí)，第1、2、3支路的平流流量分別為:10.021、10.028和9.953 L/min。無泄壓槽時(shí)，第1、2、3支路的平流流量分別為:10.154、9.657和10.19 L/min。說明了泄壓槽對(duì)3個(gè)支路的均流有明顯的效果。

圖5是30 L/min入口流量時(shí)，有、無泄壓槽兩種情況分別采用標(biāo)準(zhǔn)k-ε模型計(jì)算得到的截面壓力分布圖。有泄壓槽的壓降為32.23 kPa，比無泄壓槽的壓降 (35.18 kPa)稍小。

圖5 壓力分布結(jié)果

為了確定泄壓槽在不同入口流量和流動(dòng)狀態(tài)下的均流效果，分別用層流模型和湍流模型對(duì)10～60 L/min入口流量進(jìn)行了計(jì)算。表1列出了有、無泄壓槽的流量分配情況，可以看出:有泄壓槽時(shí)支路1、支路2和支路3的平均流速相差不足0.5%;無泄壓槽時(shí)支路1、支路2和支路3的平均流速可以相差約5%。這說明了泄壓槽可以調(diào)節(jié)各支路的壓力，使得有泄壓槽的均流效果要明顯優(yōu)于無泄壓槽的。

表1 各支路的平均流速 m/s

圖6是對(duì)入口流量為10～60 L/min、有泄壓槽水冷散熱器進(jìn)行仿真得到的壓力特性曲線。入口流量逐漸增加時(shí)，壓降迅速增大。層流模型的壓降比湍流模型要小，兩種流動(dòng)狀態(tài)模型的壓降差會(huì)隨著入口流量的增加而增大。

圖6 有泄壓槽水冷散熱器的壓力特性曲線

4 散熱性能的仿真與分析

從前面的分析結(jié)果可知，泄壓槽保證了3個(gè)支路流量的均勻分配，整體水冷散熱器的性能基本上可以以一個(gè)支路的散熱性能來代表。圖7—10是水冷散熱器在入口總流量為30 L/min(折算成熱仿真模型入口流速為0.527 m/s)采用標(biāo)準(zhǔn)k-ε模型進(jìn)行仿真得到的計(jì)算結(jié)果。

圖7表示中間截面的流速分布，最大流速為2.03 m/s?？梢钥闯隼鋮s介質(zhì)非常均勻地對(duì)叉排的翅柱進(jìn)行沖刷。

圖7 中間截面流速分布

圖8表示中間截面壓力分布，壓降為30.55 kPa，與整個(gè)水冷散熱器的流阻 (32.23 kPa)相比，要小一些，這是由于水冷散熱器入口段、出口段和泄壓槽也會(huì)產(chǎn)生一定的流動(dòng)阻力。

圖8 中間截面壓力分布

圖9 IGBT元件安裝面溫度

圖10 表面溫度

圖9是溫度場(chǎng)分布結(jié)果。圖9表示IGBT元件安裝面的溫度場(chǎng)分布，溫度標(biāo)尺采用了所見即所得方式，最高溫度為73.4℃，溫升為23.4℃，可以說明2個(gè)IGBT元件安裝面的溫度分布并不一致，且不均勻，最高溫度主要出現(xiàn)在芯片對(duì)應(yīng)位置的正下方。圖10表示整體溫度場(chǎng)分布，IGBT元件芯片上的溫度最高，為105.96℃，溫升為55.96℃。分別用層流模型和標(biāo)準(zhǔn)k-ε湍流模型對(duì)0.175 7～1.054 m/s的入口流速 (對(duì)應(yīng)水冷散熱器的入口流量為10～60 L/min)進(jìn)行計(jì)算，可以得到IGBT元件芯片最高溫度和IGBT元件安裝面最高溫度隨入口流量變化的曲線，分別見圖11和圖12。隨著入口流量由10 L/min增加至30 L/min，溫度下降速度較快，入口流量繼續(xù)由30 L/min往上增加時(shí)，溫度下降的速度變得非常緩慢。隨著入口流量由10 L/min增加至60 L/min，標(biāo)準(zhǔn)k-ε湍流模型得到的溫度與層流模型的差別越來越大，湍流效應(yīng)的影響越來越強(qiáng)。此外，還可看出入口流量為10 L/min時(shí)，芯片最高溫度會(huì)超出允許的125℃結(jié)溫限制，過小的入口流量很容易造成IGBT元件因溫度過高而損壞。

圖11 芯片最高溫度隨入口流量變化曲線

圖12 安裝面最高溫度隨入口流量變化曲線

5 結(jié)論

由于翅柱的細(xì)節(jié)特征尺寸很小，限制了網(wǎng)格的尺寸，難以對(duì)整個(gè)水冷散熱器建立完整模型進(jìn)行仿真計(jì)算，因此，文中采取了先對(duì)流體區(qū)域進(jìn)行計(jì)算，對(duì)比分析有、無泄壓槽對(duì)3個(gè)支路流量分配的情況，發(fā)現(xiàn)泄壓槽方式可以有效調(diào)節(jié)3個(gè)支路的壓力和流量，使得3個(gè)支路的流量基本上相等?；诹髁肯嗟鹊挠?jì)算結(jié)果，可以充分利用對(duì)稱性來建立熱仿真模型進(jìn)行散熱性能的仿真，得到了水冷散熱器在不同入口流量下的芯片最高溫度和IGBT元件安裝面最高溫度變化曲線。該方法的應(yīng)用，大大降低了對(duì)計(jì)算機(jī)資源的需求，有效解決和實(shí)現(xiàn)了復(fù)雜問題的仿真計(jì)算。

［1］HAN Minsub，LEE Su-Dong，HONG Chanook，et al.Development of Water-cooled Heat Sink for High-Power IGBT Inverter［C］//The 7th International Conference on Power E-lectronics，Korea:Daegu，2008:295-299.

［2］SALEM T E，BAYNE S B，PORSEHET D，et al.Thermal Performance of Water-cooled Heat Sinks［C］//The 20th Annual IEEE Applied Power Electronics Conference and Exposition，Texas，2005:129-133.

［3］唐亞男.填充式液冷散熱器的傳熱強(qiáng)化及流阻性能研究［D］.廣州:華南理工大學(xué)，2010.

［4］張建剛.航空電子設(shè)備新型結(jié)構(gòu)功能模塊的設(shè)計(jì)開發(fā)與研究［D］.西安:電子科技大學(xué)，2009.

［5］程勁嘉.綜合模塊化航空電子的液冷設(shè)計(jì)［J］.電訊技術(shù)，2011，51(6):151-155.

［6］張平，唐良寶.矩形微通道冷板的實(shí)驗(yàn)研究［J］.武漢理工大學(xué)學(xué)報(bào)，2009，31(15):91-94.

［7］李文頂.中壓礦用變頻器主電路損耗分析及散熱設(shè)計(jì)［D］.上海:上海交通大學(xué)，2009.

［8］余小玲.電力電子集成模塊及新型翅柱復(fù)合型散熱器的傳熱性能研究［D］.西安:西安交通大學(xué)，2005.

［9］葛浩.新型微通道熱沉的設(shè)計(jì)和數(shù)值模擬［D］.上海:上海交通大學(xué)，2007.

［10］揭貴生，孫馳，汪光森，等.大容量電力電子裝置中板式水冷散熱器的優(yōu)化設(shè)計(jì)［J］.機(jī)械工程學(xué)報(bào)，2010，46(2):99-105.

［11］丁杰，何多昌，唐玉兔.翅柱式水冷散熱器的性能研究［J］.大功率變流技術(shù)，2013(2):32-36.

［12］丁杰，李江紅，陳燕平，等.流動(dòng)狀態(tài)與熱源簡(jiǎn)化方式對(duì)IGBT水冷板仿真結(jié)果的影響［J］.機(jī)車電傳動(dòng)，2011(5):21-25.