計(jì)算機(jī)水冷系統(tǒng)散熱器的數(shù)值模擬及結(jié)構(gòu)優(yōu)化研究

1、    計(jì)算機(jī)水冷系統(tǒng)散熱器的數(shù)值模擬及結(jié)構(gòu)優(yōu)化研究    摘 要：本文對(duì)某公司應(yīng)用于計(jì)算機(jī)水冷散熱系統(tǒng)的管帶式交叉流散熱器進(jìn)行性能測(cè)試和研究，采用數(shù)值模擬分析法，利用simulation cfd軟件對(duì)該散熱器3d模型的流場(chǎng)、溫度場(chǎng)、壓力場(chǎng)和不同邊界條件下進(jìn)行數(shù)值模擬，將分析結(jié)果進(jìn)行可視化處理，分析得出散熱器的最佳工況點(diǎn);并據(jù)此利用pro/e三維設(shè)計(jì)軟件對(duì)散熱器3d模型的局部結(jié)構(gòu)進(jìn)行調(diào)整和優(yōu)化，得出散熱器優(yōu)化前、后的性能參數(shù)。比較結(jié)果發(fā)現(xiàn)，優(yōu)化后的散熱器優(yōu)于原設(shè)計(jì)散熱器的散熱性能。關(guān)鍵詞：計(jì)算機(jī)水冷系統(tǒng);管帶式交叉流散熱器;cfd;結(jié)構(gòu)優(yōu)化;數(shù)

2、值模擬computer numerical simulation of water cooling system radiator and structure optimization researchli zhuoxichina energy engineering group guang dong power engineering co.，ltd guangdongguangzhou 510735abstract：in this paper，the performance test and research of the tubular-belt cross-flow

3、 radiator used in the computer water-cooled radiation system of a company are carried out，using the method of numerical simulation，and using the simulation of cfd software to the radiator 3d model of flow field，temperature field，pressure field and different boundary conditions to carry on the numeri

4、cal simulation，the analysis results are visual processing and data analysis;and on the basis of use of pro/e 3d design software for radiator 3d model of the local structure adjustment and optimization，before and after optimization in the radiator performance parameters are obtained.the comparison re

5、sults show that the optimized radiator is superior to the heat dissipation performance of the original radiator.key words：computer water-cooling system;pipe and belt type cross flow radiator;cfd;structure optimization;numerical simulation1 緒論隨著計(jì)算機(jī)和電路集成技術(shù)的不斷發(fā)展，當(dāng)前數(shù)據(jù)中心服務(wù)器和計(jì)算機(jī)的使用性能也不斷增強(qiáng)。散熱問(wèn)題已經(jīng)成為計(jì)算機(jī)芯片開(kāi)發(fā)研

6、制中的瓶頸，它直接影響芯片的性能、使用壽命和生產(chǎn)成本。有效地傳遞芯片發(fā)出的熱量，使芯片在規(guī)定的溫度極限內(nèi)工作，這在當(dāng)今隨著 cpu集成度越來(lái)越大是極為重要的。1管帶式水冷散熱器具有重量輕、結(jié)構(gòu)緊湊、散熱性能好、使用壽命長(zhǎng)等優(yōu)點(diǎn)，已在國(guó)內(nèi)外廣泛應(yīng)用。2通過(guò)三維數(shù)值模擬的方法，對(duì)百葉窗翅片的傳熱和冷卻空氣的對(duì)流換熱進(jìn)行耦合，分析了百葉窗翅片的幾何參數(shù)對(duì)散熱器性能的影響，百葉窗角度為27°，百葉窗間距和翅片間距的比值為1.25時(shí)，具有較好的傳熱性能和綜合性能3;阻力損失主要集中在空氣進(jìn)入百葉窗的入口區(qū)域，百葉窗的前端傳熱效果較好。12 物理模型的建立及數(shù)值模擬分析數(shù)值模擬分析的目的是通過(guò)模

7、擬分析得出散熱器的流阻特性和傳熱特性，為優(yōu)化產(chǎn)品提供數(shù)據(jù)依據(jù)。管帶式散熱器的主要散熱結(jié)構(gòu)為體芯部分，根據(jù)文獻(xiàn)5，帶百葉窗結(jié)構(gòu)的管帶式體芯結(jié)構(gòu)尺寸的影響因素主要有百葉窗開(kāi)口尺寸a和b、開(kāi)口角度、翅片高度h、波距w、波峰數(shù)m、通道數(shù)n等。將這些參數(shù)定為設(shè)計(jì)變量：x=（a，b，h，w，m，n）n（1）式中：n為設(shè)計(jì)變量個(gè)數(shù)。根據(jù)設(shè)計(jì)要求，本文僅對(duì)百葉窗開(kāi)口尺寸a和b、開(kāi)口角度這三個(gè)變量做調(diào)整優(yōu)化，其它邊界條件有：送風(fēng)速度vz（送風(fēng)量）、送風(fēng)溫度tz、進(jìn)水流速vs（進(jìn)水流量）和進(jìn)水溫度ts，下文將對(duì)翅片結(jié)構(gòu)、體芯單元結(jié)構(gòu)和散熱器整體結(jié)構(gòu)進(jìn)行數(shù)學(xué)模型分析。2.1 翅片結(jié)構(gòu)的數(shù)學(xué)模型分析管帶翅片的百葉窗結(jié)

8、構(gòu)復(fù)雜，本文對(duì)散熱翅片單獨(dú)進(jìn)行結(jié)構(gòu)分析并進(jìn)行相應(yīng)的簡(jiǎn)化處理。實(shí)驗(yàn)?zāi)Ｐ鸵?jiàn)圖1。該風(fēng)洞模擬實(shí)驗(yàn)旨在分析百葉窗開(kāi)口角度、a和b尺寸對(duì)翅片的流阻特性和傳熱特性的影響，換熱介質(zhì)為空氣，遵循范德瓦爾方程的變化規(guī)律。（1）分析模型的變量值為開(kāi)窗角度，分別對(duì)10°、15°、20°、25°、30°、35°、40°七個(gè)不同開(kāi)角的模型進(jìn)行模擬分析，選擇k-標(biāo)準(zhǔn)湍流模型，每一組進(jìn)行三次模擬運(yùn)算，一共七組，數(shù)據(jù)結(jié)果取平均值。（2）為分析在不同邊界條件對(duì)百葉窗結(jié)構(gòu)性能的影響，以開(kāi)窗角度25°為例，設(shè)置進(jìn)風(fēng)速度從120 m/s進(jìn)行不同的仿真分析

9、。主要目的是得到不同風(fēng)速對(duì)壓力降的影響。（3）a和b的值分別為1.75mm和6mm，調(diào)整后的值為15mm和5.7mm，翅片結(jié)構(gòu)開(kāi)窗角度為25°。2.2 體芯單元結(jié)構(gòu)的數(shù)學(xué)模型分析圖2為散熱器體芯單元的仿真分析模型，模型中管帶翅片的百葉窗開(kāi)角為22°，結(jié)構(gòu)調(diào)整后為25°，a和b調(diào)整后的值為1.5mm和5.7mm。該模型模擬分析體芯單元部分在結(jié)構(gòu)調(diào)整前后溫度場(chǎng)和速度場(chǎng)的分布情況及傳熱特性和流阻特性，變量部分為模型的結(jié)構(gòu)部分，仿真模擬分為調(diào)整前后兩組，選擇的邊界條件不變。邊界條件：空氣流速為5m/s，溫度為20，水流速度為0.47m/s，進(jìn)水溫度70。網(wǎng)格劃分為混合網(wǎng)格

10、結(jié)構(gòu)，進(jìn)行網(wǎng)格滑移處理。2.3 散熱器整體數(shù)學(xué)模型分析為了方便求解前處理和減少迭代計(jì)算步驟，現(xiàn)將模型中的非關(guān)鍵部位的圓角設(shè)置成直角，將側(cè)板部分剪除，管帶部分將波峰處的圓弧狀改為梯形狀，采用邊界層替換法7簡(jiǎn)化百葉窗結(jié)構(gòu)，如圖3所示。對(duì)散熱器整體的數(shù)值模擬分析，可得到流場(chǎng)和溫度場(chǎng)的分布情況。對(duì)模型不同的邊界條件進(jìn)行模擬仿真，設(shè)置進(jìn)風(fēng)速度和進(jìn)風(fēng)溫度分別為5m/s和25，進(jìn)水溫度為70，進(jìn)水流量分別為0.5l/min、0.75 l/min、1.0 l/min、1.5 l/min、20 l/min、2.5 l/min共六組。3 數(shù)值模擬分析結(jié)果與數(shù)據(jù)處理在上一章中對(duì)數(shù)值模擬分析模型的建立和求解進(jìn)行了介紹

11、，經(jīng)過(guò)計(jì)算機(jī)中迭代運(yùn)算，得出收斂的仿真結(jié)果。結(jié)果為云圖、截圖、表格和曲線(xiàn)形式，經(jīng)分析對(duì)比得出結(jié)論，最后與實(shí)驗(yàn)結(jié)果相比較，驗(yàn)證數(shù)值模擬的準(zhǔn)確性和可靠性。3.1 散熱器管帶百葉窗結(jié)構(gòu)與空氣流速的影響（1）如圖5為在空氣流速為5m/s下，不同開(kāi)窗角的溫度場(chǎng)分布，分別得出10°40°七個(gè)數(shù)值點(diǎn)的分析結(jié)果。隨著開(kāi)角增大，對(duì)空氣的擾流作用也隨之增強(qiáng)，對(duì)邊界層的破壞能力也強(qiáng)。在對(duì)不同百葉窗開(kāi)角的模型風(fēng)洞模擬分析中，改變開(kāi)窗角度對(duì)模型的面積的影響可忽略不計(jì)。根據(jù)比熱容公式及牛頓冷卻公式可知，風(fēng)洞進(jìn)出口的溫差大小，即表現(xiàn)出模型的散熱狀況。已知進(jìn)口溫度都為20，在模型中測(cè)量空氣出口的溫度即可得

12、出傳熱溫差，表1為不同開(kāi)窗角度下的空氣出口溫度和模擬風(fēng)洞的壓力降數(shù)值。（2）將表1中的出口溫度與壓降的平均值用坐標(biāo)圖表示（圖5），空氣的出口溫度隨開(kāi)窗角度的增大而增大，擬合曲線(xiàn)的拐點(diǎn)出現(xiàn)在25°，在25°以下，曲線(xiàn)為下凹形狀，在這階段中增大開(kāi)角總是有利于改善傳熱效果的;當(dāng)開(kāi)角大于25°之后，曲線(xiàn)為上凹形狀，隨著開(kāi)角的繼續(xù)增大，影響作用力減小，傳熱效果的改善作用逐漸減弱，最后趨于一個(gè)穩(wěn)定的狀態(tài)。壓降隨開(kāi)窗角度的增大而增大，接近25°點(diǎn)時(shí)趨于平緩，曲線(xiàn)為上凹形狀。當(dāng)開(kāi)角繼續(xù)增大后，壓力降大幅上升，曲線(xiàn)為下凹形狀，說(shuō)明壓力降隨著開(kāi)角增大而急劇上升，能量損失則越

13、大。綜上所述，散熱器管帶翅片百葉結(jié)構(gòu)的開(kāi)窗角度的最佳值為25°，但考慮到加工工藝條件、技術(shù)的限制以及壓力損失和傳熱效果的綜合影響，開(kāi)窗角度為23°27°之間的綜合改善效果基本一致。（3）以開(kāi)窗角25°為例，對(duì)其進(jìn)行不同空氣流速的仿真，得出壓力降隨空氣流速的數(shù)值，根據(jù)動(dòng)能方程式ek=12mv2，消耗的理比機(jī)械功與速度的平方的二分之一成正比，如表2所示：利用graph數(shù)值函數(shù)軟件將表2中的壓力降和理論比功耗的數(shù)值擬合成函數(shù)，得到壓力降和空氣流速的擬合函數(shù)為f（x）=0.4x2，擬合因數(shù)r2=1，即與二次函數(shù)的相似度為99.07%。比較結(jié)果發(fā)現(xiàn)壓力降和理論比功

14、耗可能有著一定的比值關(guān)系。上述表明，隨著空氣流速的增大，散熱器換熱量也就越大，但是換熱效率隨之下降，在選擇滿(mǎn)足換熱器的功率要求的空氣流速時(shí)，應(yīng)考慮節(jié)能的問(wèn)題。（4）圖6為調(diào)整百葉窗結(jié)構(gòu)a和b前后的溫度場(chǎng)分布截面圖，由兩圖對(duì)比可以看出，結(jié)構(gòu)調(diào)整后翅片的溫度場(chǎng)分布要比調(diào)整前更加趨于均勻，當(dāng)量的有效的散熱面積更大，數(shù)據(jù)結(jié)果表明管帶翅片得到更好的冷卻散熱效果。3.2 散熱器體芯單元的數(shù)值模擬結(jié)果圖7模型的溫度場(chǎng)分布情況。該體芯單元數(shù)值模擬結(jié)果數(shù)據(jù)表明，其他條件不變時(shí)，結(jié)構(gòu)優(yōu)化前的水側(cè)進(jìn)出口溫差為0.743，散熱功率為23.7w，結(jié)構(gòu)優(yōu)化后的水側(cè)進(jìn)出口溫差為0.917，散熱功率為29.3w，散熱性能提高

15、23.6%。3.3 整體散熱器的數(shù)值模擬分析結(jié)果（1）圖8為散熱器的穩(wěn)態(tài)分析的溫度場(chǎng)切片圖，從圖中可分析得出，由于中間的空氣擾流情況復(fù)雜，空氣經(jīng)散熱器換熱后中間部分的空氣溫度明顯高于兩側(cè)的溫度，從中間向兩側(cè)的溫度逐漸降低。在整個(gè)空氣流域中溫度的等值面應(yīng)為四面雙曲線(xiàn)錐形結(jié)構(gòu)，散熱器前側(cè)的溫度低于后側(cè)的溫度。（2）圖9為仿真分析的流場(chǎng)分布切片圖，從圖中可直觀得出空氣流經(jīng)散熱器時(shí)的流動(dòng)情況?？諝饬鹘?jīng)散熱器管帶翅片時(shí)受到阻力作用而減小，同時(shí)空氣的流動(dòng)受到強(qiáng)烈干擾，增強(qiáng)了換熱效果。（3）進(jìn)水流量對(duì)散熱器散熱性能的影響。在定風(fēng)量的情況下，選擇不同水流量對(duì)模型進(jìn)行仿真分析，設(shè)定進(jìn)水溫度為70恒溫進(jìn)水，進(jìn)水流

16、量和模擬結(jié)果見(jiàn)表3。將表3中數(shù)據(jù)整理到坐標(biāo)圖中，如圖10為散熱器的進(jìn)水流量和散熱量的關(guān)系，隨著水流量的增加，散熱器散熱量先增后減，在進(jìn)水流量為1.0l/min時(shí)有極值，表明該散熱器在進(jìn)水流量約為1.0l/min（溫度為70）的散熱性能為最佳工況值。圖中顯示了結(jié)果的模擬值和實(shí)驗(yàn)值的對(duì)比關(guān)系，兩者的整體變化情況基本一致，說(shuō)明數(shù)值模擬與實(shí)驗(yàn)的吻合性。4 結(jié)語(yǔ)通過(guò)利用simulation cfd軟件對(duì)管帶式散熱器流場(chǎng)和溫度場(chǎng)進(jìn)行了數(shù)值模擬計(jì)算，并對(duì)散熱器的結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計(jì)，得到了散熱單元流場(chǎng)的內(nèi)部流動(dòng)細(xì)節(jié)與散熱器整體的換熱分析計(jì)算數(shù)據(jù)，最后和實(shí)驗(yàn)結(jié)果對(duì)比，各項(xiàng)數(shù)據(jù)偏差范圍合理，驗(yàn)證計(jì)算方法可行?？偨Y(jié)全文：（1）散熱器散熱管帶的百葉窗開(kāi)角在25±2°區(qū)間的散熱效果最佳;（2）調(diào)整百葉窗結(jié)構(gòu)的開(kāi)窗尺寸a和b，可提高散熱性能;（3）當(dāng)選擇進(jìn)水流量約為1.0l/min（水溫為70）時(shí)，散熱器的散熱量為最大值。參考文獻(xiàn)：1李世偉，吳國(guó)榮，袁津震，