張輝,趙萬東,高旭
(珠海格力電器股份有限公司,廣東珠海 519070)
電子模塊散熱器基座厚度對換熱性能的影響分析
張輝,趙萬東*,高旭
(珠海格力電器股份有限公司,廣東珠海 519070)
本文應(yīng)用Fluent對不同基座厚度的空調(diào)控制器電子模塊散熱器進(jìn)行換熱性能仿真分析。計算結(jié)果表明,在散熱器總高不變的前提下,基座厚度為5 mm時,散熱器換熱性能最佳。通過功耗元件溫升實驗進(jìn)一步驗證后表明,仿真所得結(jié)論與實驗結(jié)果基本一致。
散熱器;換熱性能;數(shù)值模擬;溫升實驗
變頻空調(diào)運行時,控制器中的整流橋、二極管、絕緣柵雙極晶體管(Insulated Gate Bipolar Transistor,IGBT)及智能功率模塊(Intelligent Power Module,IPM)均會產(chǎn)生較大熱量,導(dǎo)致元器件溫度急劇上升,若熱量無法及時排出,電子元器件的可靠性會急劇下降[1]。有研究指出,單個電子元器件的溫度每升高10 ℃,其失效率將增加1倍[2],致使空調(diào)性能降低,嚴(yán)重時將導(dǎo)致空調(diào)停機(jī)。因此,功耗元器件的散熱器設(shè)計顯得尤為重要。
空調(diào)控制器電子模塊用散熱器形式一般為肋片式散熱器,前人對該類型散熱器的設(shè)計作出了較為廣泛的研究,國外方面,LEUNG等[3]運用實驗測量的方法得出了散熱器翅片間的最佳間距。ANDREW[4]對自然對流散熱器的熱性能進(jìn)行分析,并仿真計算了散熱器處于熱穩(wěn)態(tài)時的相關(guān)熱性能參數(shù),對散熱器直肋尺寸進(jìn)行了優(yōu)化。POULIKAKOS等[5]在強(qiáng)制對流前提下,以熵產(chǎn)率最低為目標(biāo),確定了最佳散熱器結(jié)構(gòu)尺寸。國內(nèi)方面,吳俊鴻等通過分析變頻空調(diào)控制器的智能功率模塊溫度過熱導(dǎo)致壓縮機(jī)頻率受到限制的實例出發(fā),提出了通過合理設(shè)計保護(hù)參數(shù)、散熱器尺寸和形狀合理設(shè)計、板發(fā)熱元器件位置合理設(shè)計、采用散熱性能良好的熱傳導(dǎo)材料等有效手段降低或解決變頻控制器發(fā)熱元器件過熱問題[6]。謝少英等[7]以熱阻最低為優(yōu)化目標(biāo)采用復(fù)合形優(yōu)化算法對散熱器各個幾何參數(shù)對其熱特性的影響進(jìn)行了分析,并用實驗驗證了方法的準(zhǔn)確性。韓寧等[8]采用數(shù)值方法對型材散射器的三維流場及溫度場進(jìn)行了分析計算,定量分析了結(jié)構(gòu)因素對散熱器熱阻的影響。李玉寶等[9]考察了矩形肋片散熱器幾何參數(shù)對散熱效果的影響規(guī)律,認(rèn)為通過調(diào)整合適的肋片高度、長度和間距可有效降低散熱器的熱阻。高一博[10]結(jié)合數(shù)值模擬,在自然對流工況下,運用最少材料法、最小空間體積法、最小熱阻法和層次分析法優(yōu)化直肋片散熱器尺寸,與實驗結(jié)果對比驗證后表明肋片間距與肋片高度的比值對散熱器的換熱性能具有重要影 。
為進(jìn)一步優(yōu)化空調(diào)室外機(jī)控制器電子模塊散熱器的換熱性能,仿真分析了肋片散熱器基座厚度對其換熱性能的影響,并用實驗測試了各基座厚度下整流橋、二極管、IGBT及IPM的溫升效果,結(jié)果表明,仿真所得結(jié)論與實驗結(jié)果較為吻合,在基座厚度為5 mm時,元器件表面溫度降低,散熱器性能得到進(jìn)一步提升。
傳統(tǒng)的制冷空調(diào)產(chǎn)品設(shè)計中主要依賴樣機(jī)的反復(fù)制作與調(diào)試,既增加了開發(fā)費用,又談不上產(chǎn)品的新性能優(yōu)化,不利于節(jié)能節(jié)材。通過計算機(jī)仿真與優(yōu)化等技術(shù)的應(yīng)用,首先在計算機(jī)上面實現(xiàn)產(chǎn)品的優(yōu)化設(shè)計,可以有效地提高優(yōu)化設(shè)計的效率,是制冷空調(diào)產(chǎn)品設(shè)計方法現(xiàn)代化的發(fā)展方向[11]。
對于不可壓縮、常物性、無內(nèi)熱源的流動,對流傳熱的完整微分方程組可描述如下[12-13]。
質(zhì)量守恒方程,又稱連續(xù)性方程:
動量守恒方程:
能量守恒方程:
在固體區(qū)域,F(xiàn)luent中的能量輸運方程如式(6)中所示:
式中:
ρ——固體密度,kg/m3;
h——顯焓,K;
k——熱傳導(dǎo)系數(shù),W/(m·K);
T——固體溫度,K;
Shh— —體熱源產(chǎn)生的能量,J。
式(7)中 Tref的取值取決于求解器的選取,當(dāng)選用壓力基求解器時,其值為298.15 K,而選用密度基求解器時,其值為0 K。
圖1(a)中所示為散熱器在電氣盒中的位置,氣流由圖中左側(cè)通過進(jìn)風(fēng)柵條進(jìn)入散熱器盒蓋,通過散熱器表面后帶走熱量,從盒蓋后方流出,圖1(b)為散熱器模型,屬于直肋型擠塑散熱器,肋片與基座總高為40 mm,基座高度為4mm,肋片數(shù)16片。
圖1 初始模型
圖2為整流橋、二極管、IGBT及IPM在散熱器底部的布置位置示意圖。
由于電器盒結(jié)構(gòu)復(fù)雜,殼體上布滿加強(qiáng)筋及凹槽等復(fù)雜細(xì)小的結(jié)構(gòu),無疑會增加網(wǎng)格數(shù)量,且難以保證網(wǎng)格質(zhì)量。因此,僅截取包裹散熱器的盒蓋部分作為本次仿真的計算域,并忽略空氣入口處的格柵結(jié)構(gòu),簡化后計算模型如圖3所示。
圖2 功耗元器件布置位置
圖3 模型簡化
由于存在溫度傳遞,因此必須開啟能量方程,流動模型選取通過雷諾數(shù) Re計算進(jìn)行判斷,使用熱線風(fēng)速儀對散熱器入口處風(fēng)速進(jìn)行測量,整理后取平均值為1.9 m/s,由雷諾數(shù)計算公式[14]可得入口處Re=9,575.4,屬于湍流范疇,因此,選取k-ε雙方程湍流模型進(jìn)行流動計算。在穩(wěn)態(tài)工況下分析計算研究對象的溫度場,因此所有微分方程中忽略時間項的影響[15]。
為模擬風(fēng)機(jī)的抽吸作用,入口選取自由流入口邊界,入口溫度為 43 ℃,出口為負(fù)壓出口,負(fù)壓值為-10 Pa,散熱器與空氣接觸表面采用耦合壁面,其余壁面為絕熱壁面,為對比散熱器的換熱量及傳熱系數(shù),將各功耗元器件做定壁溫處理,各元器件溫度參數(shù)采用實測值,保證了熱源一致性,邊界條件具體設(shè)置如表1中所示。
表1 數(shù)學(xué)模型及邊界條件設(shè)置
圖4為散熱器方案調(diào)整示意圖,保持散熱器總高H不變,調(diào)整基座厚度d的值分別為4 mm、5 mm、6 mm及7 mm,共計4個計算方案,考察換熱量及換熱系數(shù)的變化情況。
由于散熱器肋片與基座間存在圓弧過渡段,容易引起網(wǎng)格畸變,因此采用四面體網(wǎng)格對空間進(jìn)行離散,三角形網(wǎng)格對面進(jìn)行劃分,并在肋片表面拉升出4層棱柱體網(wǎng)格,最終生成網(wǎng)格總數(shù)達(dá)572萬,節(jié)點數(shù)92萬以上,截面網(wǎng)格如圖5中所示。
圖4 方案調(diào)整示意圖
圖5 肋片表面邊界層網(wǎng)格
圖6為不同基座厚度下散熱器表面的溫度云圖,隨著基座厚度的增加,高溫區(qū)域在肋片表面的分布面積越來越小,基座局部高溫區(qū)域縮小,散熱器表面最高溫度呈下降趨勢。
各方案中,散熱器與各功耗元器件的接觸熱阻均相同,但基座厚度越厚,導(dǎo)熱面積熱阻越大,溫度在散熱器上的傳導(dǎo)阻力越大,因此出現(xiàn)了基座越厚、高溫區(qū)域面積越小的現(xiàn)象。但是換熱器整體的換熱性能優(yōu)劣不能僅 通過表面溫度及基座熱阻衡量,還需綜合考慮散熱器結(jié)構(gòu)不同時肋片在空氣側(cè)換熱面積的影響。理論上,在總高一定的前提下,存在一個最優(yōu)的基座厚度及肋片高度,使得散熱器換熱效果最佳。
圖6 不同基座厚度工況下的散熱器表面溫度云圖
讀取各方案散熱器表面換熱量及表面?zhèn)鳠嵯禂?shù)數(shù)值如表2中所示。從表中可知,基座厚度為5 mm時,散熱器的換熱量最高,表面?zhèn)鳠嵯禂?shù)最大。
表2 換熱性能對比
圖7為換熱量隨基座厚度變化的趨勢曲線,當(dāng)基座厚度增加時,換熱量緩慢上升,當(dāng)基座厚度大于5 mm后,換熱量衰減速率較大。
這是由于在散熱過程中,基座厚度過厚,則元器件產(chǎn)生的熱量積聚在基座處,難以通過肋片將熱量耗散到空氣中,導(dǎo)致整體換熱量下降,基座厚度過薄時,元器件熱量傳導(dǎo)不充分,未能有效利用到每片肋片的換熱面積,同樣會導(dǎo)致?lián)Q熱量的下降。
圖7 換熱量隨基座厚 變化曲線
將外機(jī)與室內(nèi)機(jī)相連,室內(nèi)環(huán)溫 32℃,室內(nèi)機(jī)送風(fēng)溫度23 ℃,室外環(huán)溫43 ℃,實驗方案如表3中所示,方案1為基座厚度4 mm且肋片數(shù)為15片,方案2至方案5分別為基座厚度4 mm至7 mm的散熱器,散熱器總高不變,肋片數(shù)不變。
實驗測試過程中,各功耗元器件溫度波動穩(wěn)定后再記錄溫度數(shù)據(jù),記錄完成后停機(jī),待各元器件溫度下降至常溫后,再進(jìn)行下一方案的測試。
表3 實驗方案
實驗測試數(shù)據(jù)如表4中所示,由實驗數(shù)據(jù)可知,相比于肋片數(shù)量而言,基座厚度對散熱效果影響不大,方案2中,肋片數(shù)減少一片,使得整流橋及二極管溫度大幅上升。方案3中,基座厚度為5 mm時,各元器件溫度均降至 90 ℃以下,散熱效果較佳,均溫性較好?;雍裰? mm厚后,IPM溫度降幅較大,整流橋溫度仍在90 ℃以上,且從成本上考慮,無需采用高基板厚度的散熱器。
表4 實驗數(shù)據(jù)
以方案1中各功耗元器件溫度總和為基準(zhǔn),分別計算方案2至5的元器件總溫,與方案1進(jìn)行對比,得到溫降幅度如表5中所示,方案3與方案5的總溫降均可達(dá)到6%,但方案3溫降幅度更大,成本更低。
表5 溫降效果對比
圖8為基座厚度增加時,各元器件的溫度變化折線圖。由圖中可知,整流橋及二極管溫度呈現(xiàn)先減小再大幅增加,之后緩慢降低的趨勢,IGBT與IPM則在基座厚度增加時,溫度持續(xù)降低,說明基座厚度增加后,對元器件產(chǎn)生的熱量進(jìn)行了儲存,但基座厚度增加,溫度下降速度放緩,效率過低,成本升高。
圖8 元器件溫度隨基座厚度變化圖
實驗所得結(jié)果與仿真結(jié)果基本一致,區(qū)別在于當(dāng)基座厚度增加至7 mm時,實驗結(jié)果顯示元器件溫升情況有一定改善,主要原因是在實驗過程中,為方便方案間的更替,人為剪除了散熱器盒蓋,導(dǎo)致散熱器一部分基座也裸露在空氣中,換熱得到加強(qiáng),而仿真模型中,空氣計算域僅包含肋片區(qū)域及一部分基座區(qū)域,造成了實驗與仿真的偏差,但仿真與實驗結(jié)果均表明,基座厚度為5 mm時散熱器換熱性能最佳。
1)仿真及實驗結(jié)果均表明,基座厚度為5 mm時,所述外機(jī)控制器電子模塊散熱器換熱性能最佳;
2)基座厚度增加至5 mm時,整流橋、二極管、IGBT及IPM的溫升均得到改善,當(dāng)基座厚度進(jìn)一步增加時,整流橋與二極管溫度顯著上升后緩慢下降,IGBT及IPM的溫度則呈現(xiàn)單調(diào)下降的趨勢,但基座加厚使得換熱效率降低,成本大幅升高;
3)肋片數(shù)對散熱器換熱性能的影響極大,后續(xù)設(shè)計散熱器時,可適當(dāng)增加肋片數(shù),肋片高度及寬度應(yīng)盡量設(shè)計到結(jié)構(gòu)允許的最大值;
4)本次仿真及實驗間存在一定差異,后續(xù)仿真方面,可將功耗元器件邊界條件由定壁溫替換為生熱率,更加貼近實際情況,實驗過程中添加散熱器盒蓋,以獲得更精確的實驗值。
[1]胡廣新. 翅片式與微流道式散熱器散熱特性及應(yīng)用研究[D]. 成都: 成都電子科技大學(xué), 2010.
[2]于慈遠(yuǎn), 于湘珍, 楊為民. 電子設(shè)備熱分析[J]. 微電子學(xué), 2000, 30(5): 334-337.
[3]LEUNG C W, ROBERT S D. Heat-exchanger performance: effect of orientation[J]. Applied Energy,1989, 33: 235-252.
[4]ANDREW T M. Optimization of heat sink fin geometry for heat sinks in natural convection[C]. Inter Society Conference on Thermal Phenomena, New York, 1992:145-148.
[5]POULIKAKOS A, BEJAN A. Fin geometry for minimum entropy generation in forced convection[J]. ASME Journal of Heat Transfer, 1982, 104: 616-623.
[6]吳俊鴻, 楊檢群等. 變頻空調(diào)控制器過熱問題的分析及對策[J]. 制冷技術(shù), 2012, 32(4): 51-54.
[7]謝少英, 趙惇殳, 王世萍. 型材散熱器的優(yōu)化設(shè)計[J].電子機(jī)械工程, 2001, 91(3): 28-32.
[8]韓寧, 余墨娟, 趙惇殳, 等. 型材散熱器熱特性分析[J].西安電子科技大學(xué)學(xué)報(自然科學(xué)版), 2002, 29(4):551-555.
[9]李玉寶, 王建萍, 呂召會. 矩形肋片散熱器幾何參數(shù)對散熱的影響分析[J]. 電子機(jī)械工程, 2012, 28(2): 4-7.
[10]高一博. 直肋自然對流散熱器的優(yōu)化研究[D]. 武漢:華中科技大學(xué), 2008.
[11]胡海濤, 丁國良, 王如竹. 第五屆全國制冷空調(diào)新技術(shù)會議反映的研究熱點[J]. 制冷技術(shù), 2008, 28(4):19-25.
[12]王漢青. 暖通空調(diào)流體流動數(shù)值計算方法與應(yīng)用[M].北京: 科學(xué)出版社, 2013.10: 14-20.
[13]楊世銘, 陶文銓. 傳熱學(xué)[M]. 4版. 北京: 高等教育出版社, 2006: 203-206.
[14]莫乃榕. 工程流體力學(xué) [M]. 2版. 武漢: 華中科技大學(xué)出版社, 2009: 74-76.
[15]王棟. 電熱融霜過程中冷風(fēng)機(jī)周圍溫度場分析[J]. 制冷技術(shù), 2012, 32(4):12-14.
Analysis on Influence of Chassis Thickness on Heat Transfer Performance of Electronic Module Radiator
ZHANG Hui, ZHAO Wandong*, GAO Xu
(Gree Electric Appliances, Inc. of Zhuhai, Zhuhai, Guangdong 519070, China)
The heat transfer performance of electronic module radiator of air-condition controller with different foundation thicknesses is simulated with Fluent. Result shows that, as the height is fixed, the radiator has best heat transfer performance at foundation thickness of 5 millimeter. Through temperature rise experiment of power consumption components, the experimental results are basically proved to be the same with the simulation results.
Radiator; Heat transfer performance; Numerical simulation; Temperature rise experiment
10.3969/j.issn.2095-4468.2017.04.201
*趙萬東(1990-),男,助理工程師,碩士。研究方向:傳熱傳質(zhì)。聯(lián)系地址:珠海前山金雞西路,郵編:519070。聯(lián)系電話:0756-8522066。E-mail:zhaowandong115@163.com。