IGBT熱管式整體翅片散熱器優化分析

徐鵬程 陶漢中 張紅

(南京工業大學能源學院 南京 211816)

IGBT熱管式整體翅片散熱器優化分析

徐鵬程 陶漢中 張紅

(南京工業大學能源學院 南京 211816)

本文建立了IGBT熱管散熱器模型，并用實驗數據驗證所建立的三維熱管散熱器模型及基于有限體積法的數值模擬結果的可靠性，然后研究翅片厚度t和翅片間距s對散熱器散熱性能的影響。結果表明:數值模型可靠；當t由0.5 mm增至1 mm時整體翅片熱管散熱器的摩擦系數f增加45%，而傳熱因子j最大增加30%；當s由1.2 mm增至2 mm時散熱器的摩擦系數f減小15%，而j的變化較小且s＝1.2 mm時j的平均值最大。

熱管；散熱器；整體翅片；數值模擬

IGBT(Insulated Gate Bipolar Transistor，絕緣柵雙極型晶體管)是一種集高頻、高壓和大電流于一身的功率半導體復合元件，它的出現促進了電力電子技術的發展。但IGBT在工作過程中會有較大的功率損耗，從而產生大量的熱量，使其模塊的溫度快速升高，因此通常需要通過外部器件來輔助散熱。Bo Wang 等[1]研究了結溫對IGBT的電特性的影響，結果顯示低結溫對IGBT的性能影響較小，但超過一定溫度后對IGBT性能產生較大的危害，影響模塊使用壽命。因此良好的散熱是保證模塊長期穩定工作的關鍵因素。宗朝暉[2]概述了通用的電力電子冷卻方法，主要包括空冷散熱、液冷散熱、熱電制冷散熱以及熱管散熱。熱管是一種高效的傳熱元件，特點是結構簡單，散熱效率高，維護方便。熱管的形式有重力熱管、回路熱管、震蕩熱管以及平板熱管，Ya Ping Zhang 等[3－5]對用于電子設備散熱的這幾種形式的熱管進行了研究分析。

熱管式IGBT散熱器的研究主要分為實驗研究和數值模擬。羅愛華等[6]實驗研究了牽引整流器熱管式空氣冷卻器的散熱性能，翅片材料為鋁，熱源采用的是水蒸汽加熱，結果表明熱管散熱器能好的將溫度控制在80℃以下。X.Perpi?à等[7]通過實驗研究了牽引IGBT模塊內部的溫度分布情況，冷卻形式采用的是重力熱管雙面冷卻，結果顯示散熱器安裝在功率模塊上的方式能夠影響散熱效果。李春陽等[8]設計開發的用于CRH2型高速動車組CI用熱管散熱器，著重于散熱器蒸發段熱阻的實驗研究，且開發的散熱器熱阻小于0.005656 K/W。陶漢中等[9]、劉紅等[10]都利用有限元法對熱管散熱器進行了數值模擬，文中對散熱器的模型和邊界條件都進行了詳細的規定。結果證明熱管式散熱器性能明顯優于傳統的翅片散熱器。而熱管冷凝段擴展面的研究文獻較少，大部分文獻都集中在換熱器的翅片管研究上，張瑩等[11]通過實驗和數值模擬研究了翅片中帶有熱管的散熱器的性能，結果表明翅片中帶有熱管的散熱器既可以提高散熱器的傳熱溫差，又可使翅片散熱強化，降低CPU的溫度，但沒有對翅片的厚度和間距進行研究。

在前人研究的基礎上，文中將建立IGBT熱管散熱器的物理數學模型，利用有限體積法對散熱器進行數值模擬，并用實驗數據驗證其可靠性，分析翅片厚度和翅片間距在不同流速下對散熱器的流動和傳熱性能的影響，為用于IGBT的熱管散熱器翅片方面的優化設計提供參考。

1 物理模型與數值計算

1.1 模型建立

散熱器是針對耗散功率為1140 W的單個IGBT模塊而設計的。熱管散熱器由熱管，銅基板和翅片組成，散熱器模型結構如圖1所示，模型大小為:0≤X ≤174 mm，0≤Y≤601 mm，0≤Z≤80 mm，采用6根外徑為10 mm的熱管傳導熱量，熱管長度為600 mm，蒸發段長度為129 mm，蒸發段與銅基板相連，銅基板體積為140 mm×130 mm×30 mm，冷凝段長度為450 mm，且在冷凝段上安裝整體式翅片，翅片尺寸為170 mm×60 mm。散熱器空氣進出口布置在整個翅片段，大小為140 mm×130 mm的熱源布置在銅基板上，翅片和基板材料為紫銅，導熱系數為377 W/(m·℃)。

本文首先建立三維熱管散熱器模型，然后通過有限體積法對散熱器的散熱性能進行數值模擬研究。對模型作如下假設:系統處于穩態；流動為不可壓縮湍流流動；空氣的物性為定性溫度下的參數；散熱選用強制對流，因此重力和輻射的影響可以忽略不計；熱管定義為純導熱元件，導熱系數設為20000 W/(m·℃)。

模型邊界條件設置如下，空氣入口設定為速度邊界條件:Z軸正方向速度u＝2m/s，X、Y軸方向速度為v＝w＝0，進口空氣溫度tin＝25℃；出口為自由邊界條件；基板一側設定為面熱源Q＝1140 W，熱流密度為6.26 W/cm2；翅片和熱管采用自身的傳導對流耦合換熱；散熱器其余邊界為無滑移絕熱邊界條件。散熱器的計算模型采用k－ε兩方程模型，采用二次迎風格式，能量方程的計算收斂殘差精度為10－7，其他方程的殘差精度為10－4。

熱管散熱器的雷諾數Re，摩擦系數f，傳熱因子j及換熱系數h定義如下:

式中:u為最小流通截面處的空氣流速，m/s；ν為空氣的動力粘度，m2/s；h為換熱系數，W/(m2· ℃)；d為熱管外徑，m；Δp為空氣進出口壓降，Pa；L是沿氣流方向的翅片長度，m；th為熱源表面的平均溫度，℃，th由數值模擬的結果得出；tin為空氣進口溫度，℃；tm為空氣的定性溫度，℃；Aair為空氣側的換熱面積，m2；Q為IGBT的耗散功率，W；G為空氣質量流量，kg/s。

1.2 網格的獨立性檢驗

為了驗證網格的獨立性，邊界條件設置完成后，對所要計算的模型依次劃分了81萬，98萬，134萬，190萬，252萬等網格，經數值計算后對比發現熱源的平均溫度變化隨網格數量的增加逐漸減小，其驗證結果如圖2所示。網格由187萬增至252萬時熱源平均溫度值變化0.4%，所以這證明了網格具有獨立性。最后考慮計算時長的原因網格劃分選擇在190萬左右。

圖1 散熱器結構簡圖Fig.1 The structure of heat sink

1.3 實驗驗證

為了保證數值模擬的可靠性，采用陶漢中[12]論文中相同幾何尺寸和參數進行數值模擬，表1列出了用于驗證數值模擬的實驗散熱器的參數，實驗中散熱器的結構與圖1一致，銅水熱管的充液量為2.5 mL，實驗時將散熱器安裝在風道中，采用電加熱模擬熱源，輸入功率由變壓器控制，在熱源表面布置12個熱電偶作為測點，風量由ZRQF-J型風速儀測量。在加熱功率為1140 W，風速為2m/s時散熱器穩定工作后測得熱源平均溫度為78℃。同時考慮了銅基塊與熱管間的接觸熱阻R1＝0.0275℃/W，采用上文建立的數值模型模擬得到的熱源平均溫度為76.1℃，與實驗結果相差2.44%，則可以認為文中建立的模型和數值模擬的結果是可靠的。

表1 驗證用散熱器參數Tab.1 Parameters of heat sink for verify

圖2 網格的獨立性檢驗Fig.2 Grid independent study

2 結果與討論分析

文中研究了翅片厚度 t＝1 mm、0.8 mm、0.7 mm、0.6 mm、0.5 mm，翅片間距s＝1 mm，1.2mm、1.4 mm、1.6 mm、1.8 mm、2 mm下，進口風速分別為:u＝2 m/s、4 m/s、6 m/s、8 m/s時對熱管散熱器的摩擦系數f和傳熱因子j的影響。

2.1 翅片厚度對散熱器的影響

為了研究整體翅片厚度對散熱器性能的影響，保持翅片間距s＝2 mm，比較了不同雷諾數下5種翅片厚度的散熱情況。由圖3(a)可以看出隨著Re的增加，摩擦系數f均逐漸減小。摩擦系數f隨著翅片厚度的增加而增加。t由0.5mm增至1 mm時f平均增大了45%。圖3(b)展示了隨著Re的增加傳熱因子j都逐漸減小。j隨著t的增加而增加，但是在高Re下增加的較小。ln(Re)＜10.5時，t由0.5 mm增至1 mm時j平均增加了30%，ln(Re)＞10.5時，j平均增加了25%。可以看出f增加的幅度大于j增加的幅度，雖然翅片厚度的增加可以增加散熱效果，但是帶來的流動阻力大于帶來的換熱效果。

圖3 翅片厚度對散熱器換熱和流動性能的影響Fig.3 The effect of fin thickness on heat transfer and flow

表2 2 m/s下不同翅片厚度的數值模擬結果Tab.2 Numerical simulation results of different fin thickness at 2 m/s

由表2可以看出，在進口風速為2 m/s下，不同翅片厚度的散熱器均能保證熱源最大溫度低125℃，且翅片厚度t＝1 mm時熱源平均溫度最低。

2.2翅片間距對散熱器的影響

為了研究整體翅片間距對散熱器性能的影響，保持翅片厚度t＝1 mm，比較了不同雷諾數下5種不同翅片間距的散熱影響。由圖4(a)可以看出隨著Re的增加，摩擦系數f均逐漸減小。摩擦系數f隨著翅片間距s的增加而減小。s由1.2 mm增至2 mm時f平均減小了15%。圖3(b)展示了隨著Re的增加傳熱因子j都逐漸減小。但j隨著s的變化不大，s＝1.2 mm時換熱因子j最大。可以看出f和j隨著s的變化均不大，選擇s＝1.2 mm較好。

表3 2 m/s下不同翅片間距的數值模擬結果Tab.3 Numerical simulation results of different fin spacing at 2 m/s

圖4 翅片間距對散熱器換熱和流動性能的影響Fig.4 The effect of spacing on heat transfer and flow

由表3可以看出，在進口風速為2 m/s下，不同翅片間距的散熱器均能保證熱源最大溫度低125℃，且翅片間距s＝1.2 mm時熱源平均溫度最低。

3 結論

文中首先利用IGBT熱管散熱器現有的實驗數據結果驗證建立的三維熱管散熱器模型和數值模擬結果的可靠性。然后比較了不同翅片厚度，翅片間距對整體翅片熱管散熱器的性能影響，得到以下結論:

1)整體翅片熱管散熱器的摩擦系數f隨著t的增加而增加，隨著s的增大而減小，其中受t的影響大于受s的影響。j因子隨著t或s的增大而增加；

2)整體翅片熱管散熱器的傳熱因子j隨著t的增加而增加，j因子受s的影響較小，s＝1.2 mm是平均j值最大；

3)本文中t＝0.5 mm，s＝1.2 mm時散熱器的散熱效果最好。

因此，不同翅片厚度和翅片間距對熱管散熱器的傳熱與流動性能的影響結果的獲得對散熱器的設計及應用提供了參考。

[1] Bo Wang，Yong Tang，Ming Chen.Study on Electric Characteristic of IGBT at Different Junction Temperature [C]//Power and Energy Engineering Conference，2011: 1-4.

[2] 宗朝輝.現代電力電子的冷卻技術[J].變流技術與電力牽引，2007(4):6-12.(Zong Chaohui.Cooling Technology of Modern Power Electronics[J].Converter Technology＆Electric Traction，2007(4):6-12.)

[3] Yaping Zhang，Xiaoling Yu，Quanke Feng.An experimental investigation of a heat pipe integrated power module [C]//Industrial Electronics and Applications，3rd IEEE Conference，2008:820-823.

[4] Manfred Groll，Marcus Schneider，Valérie Sartre，et al. Thermal control of electronic equipment by heat pipes[J]. Revue Générale de Thermique，1998，37(5):323-352.

[5] 田金穎，諸凱，劉建林，等.冷卻電子芯片的平板熱管散熱器傳熱性能研究[J].制冷學報，2007，28(6) (6):18-22.(Tian Jinying，Zhu Kai，Liu Jianlin，et al. Experiments Investigation on Thermal Conductivity of Flat Heat Pipes Used to Cool Electric Die[J].Journal of Refrigeration，2007，28(6):18-22.)

[6] 羅愛華，金鷹，岳良，等.牽引整流器熱管式空氣冷卻器散熱性能實驗研究[J].制冷學報，2013，34 (5):90-94.(Luo Aihua，Jin Ying，Yue Liang，et al.

Experimental Study on Heat Pipe Cooler of Electric Locomotive Traction Converters[J].Journal of Refrigeration，2013，34(5):90-94.)

[7] X Perpi?à，M Piton，M Mermet-Guyennet，et al.Local thermal cycles determination in thermosyphon-cooled traction IGBT modules reproducing mission profiles[J].Microelectronics Reliability，2007，47(9/10/11):1701-1706.

[8] 李春陽，肖寧，楊閎盛.CRH2型高速動車組CI用熱管散熱器的研制[C]//第十三屆全國熱管學術會議，2012.

[9] 陶漢中，張紅，莊駿.高速芯片模塊熱管散熱器的數值傳熱分析[J].南京工業大學學報，2004，26(1):68-71.(Tao Hanzhong，Zhang Hong，Zhuang Jun.Thermal simulating analysis on heat pipe heat sink for high2speed chip module[J].Journal of Nanjing University of Technology，2004，26(1):68-71.)

[10]劉紅，童思成，蔣蘭芳.熱管散熱器數值仿真模型[J].半導體光電，2012，33(2):194-200.(Liu Hong，Tong Sicheng，Jiang Lanfang.Numerical Simulation Model for Radiator with Heat Pipe[J].Semiconductor optoelectronics，2012，33(2):194-200.)

[11]張瑩，諸凱.翅片中帶有熱管的散熱器數值模擬和實驗研究[J].低溫與超導，2013，41(2):56-61.(Zhang Ying，Zhu Kai.Numerical analysis and experimental research on finned heat sinks with heat pipes[J].Cryo.＆Supercond，2013，41(2):56-61.)

[12]陶漢中.熱管式大功率電子模塊散熱器的研究與開發[D].南京:南京工業大學，2003.

Optimization and Analysis of IGBT Heat Pipe Heat Sink with Integral Fin

Xu Pengcheng Tao Hanzhong Zhang Hong

(College of Energy，Nanjing University of Technology，Nanjing，211816，China)

The model of IGBT heat pipe heat sink was established.The established three－dimensional model and the reliability of numerical simulation results based on the finite volume method were verified by experimental data.Then the effects of fin thickness(t)and fin spacing(s)was investigated on the radiator cooling performance.The results show that the numerical model is reliable.When t is increased from 0.5 mm to 1 mm，the friction coefficient(f)of integral fin tube radiator is increased by 45%，and j factor is greatly increased by 30%.When s is increased from 1.2 mm to 2 mm，the friction coefficient(f)of radiator is decrease by 15%，but the change of j is small and its average reaches the maximum when s＝1.2 mm.

heat pipe；heat sink；integral fin；numerical simulation

TB657.5；TK172

A

0253－4339(2014)05－0101－05

10.3969/j.issn.0253－4339.2014.05.101

陶漢中，男(1978－)，副教授，南京工業大學能源學院，13776668774，E-mail:taohanzhong＠njut.edu.cn。研究方向:電子設備熱控制、太陽能熱利用。

2013年12月30日

About the corresponding author

Tao Hanzhong(1978－)，male，Associate Professor，College of Energy，Nanjing University of Technology，13776668774，E-mail:taohanzhong＠njut.edu.cn.Research fields:Thermal control of electronic equipment，Solar energy thermal application.