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

關(guān)鍵詞：計算機(jī)水冷系統(tǒng);管帶式交叉流散熱器;CFD;結(jié)構(gòu)優(yōu)化;數(shù)值模擬

Computer Numerical Simulation of Water Cooling System Radiator and Structure Optimization Research

Li Zhuoxi

CHINA ENERGY ENGINEERING GROUP GUANG DONG POWER ENGINEERING CO.，LTD?GuangdongGuangzhou?510735

Abstract：In this paper，the performance test and research of the Tubular-Belt cross-flow 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 numerical 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 results 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 simulation

1 緒論

隨著計算機(jī)和電路集成技術(shù)的不斷發(fā)展，當(dāng)前數(shù)據(jù)中心服務(wù)器和計算機(jī)的使用性能也不斷增強(qiáng)。散熱問題已經(jīng)成為計算機(jī)芯片開發(fā)研制中的瓶頸，它直接影響芯片的性能、使用壽命和生產(chǎn)成本。有效地傳遞芯片發(fā)出的熱量，使芯片在規(guī)定的溫度極限內(nèi)工作，這在當(dāng)今隨著 CPU集成度越來越大是極為重要的。[1]管帶式水冷散熱器具有重量輕、結(jié)構(gòu)緊湊、散熱性能好、使用壽命長等優(yōu)點，已在國內(nèi)外廣泛應(yīng)用。[2]

通過三維數(shù)值模擬的方法，對百葉窗翅片的傳熱和冷卻空氣的對流換熱進(jìn)行耦合，分析了百葉窗翅片的幾何參數(shù)對散熱器性能的影響，百葉窗角度為27°，百葉窗間距和翅片間距的比值為1.25時，具有較好的傳熱性能和綜合性能[3];阻力損失主要集中在空氣進(jìn)入百葉窗的入口區(qū)域，百葉窗的前端傳熱效果較好。[1]

2 物理模型的建立及數(shù)值模擬分析

數(shù)值模擬分析的目的是通過模擬分析得出散熱器的流阻特性和傳熱特性，為優(yōu)化產(chǎn)品提供數(shù)據(jù)依據(jù)。管帶式散熱器的主要散熱結(jié)構(gòu)為體芯部分，根據(jù)文獻(xiàn)[5]，帶百葉窗結(jié)構(gòu)的管帶式體芯結(jié)構(gòu)尺寸的影響因素主要有百葉窗開口尺寸A和B、開口角度β、翅片高度H、波距W、波峰數(shù)M、通道數(shù)N等。將這些參數(shù)定為設(shè)計變量：

X→=（A，B，β，H，W，M，N）n（1）

式中：n為設(shè)計變量個數(shù)。

根據(jù)設(shè)計要求，本文僅對百葉窗開口尺寸A和B、開口角度β這三個變量做調(diào)整優(yōu)化，其它邊界條件有：送風(fēng)速度Vz（送風(fēng)量）、送風(fēng)溫度Tz、進(jìn)水流速Vs（進(jìn)水流量）和進(jìn)水溫度Ts，下文將對翅片結(jié)構(gòu)、體芯單元結(jié)構(gòu)和散熱器整體結(jié)構(gòu)進(jìn)行數(shù)學(xué)模型分析。

2.1 翅片結(jié)構(gòu)的數(shù)學(xué)模型分析

管帶翅片的百葉窗結(jié)構(gòu)復(fù)雜，本文對散熱翅片單獨進(jìn)行結(jié)構(gòu)分析并進(jìn)行相應(yīng)的簡化處理。實驗?zāi)Ｐ鸵妶D1。

該風(fēng)洞模擬實驗旨在分析百葉窗開口角度β、A和B尺寸對翅片的流阻特性和傳熱特性的影響，換熱介質(zhì)為空氣，遵循范德瓦爾方程的變化規(guī)律。

（1）分析模型的變量值為開窗角度β，分別對10°、15°、20°、25°、30°、35°、40°七個不同開角的模型進(jìn)行模擬分析，選擇k-ε標(biāo)準(zhǔn)湍流模型，每一組進(jìn)行三次模擬運(yùn)算，一共七組，數(shù)據(jù)結(jié)果取平均值。

（2）為分析在不同邊界條件對百葉窗結(jié)構(gòu)性能的影響，以開窗角度25°為例，設(shè)置進(jìn)風(fēng)速度從1～20 m/s進(jìn)行不同的仿真分析。主要目的是得到不同風(fēng)速對壓力降的影響。

（3）A和B的值分別為1.75mm和6mm，調(diào)整后的值為15mm和5.7mm，翅片結(jié)構(gòu)開窗角度為25°。

2.2 體芯單元結(jié)構(gòu)的數(shù)學(xué)模型分析

圖2為散熱器體芯單元的仿真分析模型，模型中管帶翅片的百葉窗開角為22°，結(jié)構(gòu)調(diào)整后為25°，A和B調(diào)整后的值為1.5mm和5.7mm。該模型模擬分析體芯單元部分在結(jié)構(gòu)調(diào)整前后溫度場和速度場的分布情況及傳熱特性和流阻特性，變量部分為模型的結(jié)構(gòu)部分，仿真模擬分為調(diào)整前后兩組，選擇的邊界條件不變。邊界條件：空氣流速為5m/s，溫度為20℃，水流速度為0.47m/s，進(jìn)水溫度70℃。網(wǎng)格劃分為混合網(wǎng)格結(jié)構(gòu)，進(jìn)行網(wǎng)格滑移處理。

2.3 散熱器整體數(shù)學(xué)模型分析

為了方便求解前處理和減少迭代計算步驟，現(xiàn)將模型中的非關(guān)鍵部位的圓角設(shè)置成直角，將側(cè)板部分剪除，管帶部分將波峰處的圓弧狀改為梯形狀，采用邊界層替換法[7]簡化百葉窗結(jié)構(gòu)，如圖3所示。對散熱器整體的數(shù)值模擬分析，可得到流場和溫度場的分布情況。對模型不同的邊界條件進(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ù)處理

在上一章中對數(shù)值模擬分析模型的建立和求解進(jìn)行了介紹，經(jīng)過計算機(jī)中迭代運(yùn)算，得出收斂的仿真結(jié)果。結(jié)果為云圖、截圖、表格和曲線形式，經(jīng)分析對比得出結(jié)論，最后與實驗結(jié)果相比較，驗證數(shù)值模擬的準(zhǔn)確性和可靠性。

3.1 散熱器管帶百葉窗結(jié)構(gòu)與空氣流速的影響

（1）如圖5為在空氣流速為5m/s下，不同開窗角的溫度場分布，分別得出10°～40°七個數(shù)值點的分析結(jié)果。隨著開角增大，對空氣的擾流作用也隨之增強(qiáng)，對邊界層的破壞能力也強(qiáng)。

在對不同百葉窗開角的模型風(fēng)洞模擬分析中，改變開窗角度對模型的面積的影響可忽略不計。根據(jù)比熱容公式及牛頓冷卻公式可知，風(fēng)洞進(jìn)出口的溫差大小，即表現(xiàn)出模型的散熱狀況。已知進(jìn)口溫度都為20℃，在模型中測量空氣出口的溫度即可得出傳熱溫差，表1為不同開窗角度下的空氣出口溫度和模擬風(fēng)洞的壓力降數(shù)值。

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

綜上所述，散熱器管帶翅片百葉結(jié)構(gòu)的開窗角度的最佳值為25°，但考慮到加工工藝條件、技術(shù)的限制以及壓力損失和傳熱效果的綜合影響，開窗角度為23°～27°之間的綜合改善效果基本一致。

（3）以開窗角25°為例，對其進(jìn)行不同空氣流速的仿真，得出壓力降隨空氣流速的數(shù)值，根據(jù)動能方程式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)壓力降和理論比功耗可能有著一定的比值關(guān)系。

上述表明，隨著空氣流速的增大，散熱器換熱量也就越大，但是換熱效率隨之下降，在選擇滿足換熱器的功率要求的空氣流速時，應(yīng)考慮節(jié)能的問題。

（4）圖6為調(diào)整百葉窗結(jié)構(gòu)A和B前后的溫度場分布截面圖，由兩圖對比可以看出，結(jié)構(gòu)調(diào)整后翅片的溫度場分布要比調(diào)整前更加趨于均勻，當(dāng)量的有效的散熱面積更大，數(shù)據(jù)結(jié)果表明管帶翅片得到更好的冷卻散熱效果。

3.2 散熱器體芯單元的數(shù)值模擬結(jié)果

圖7模型的溫度場分布情況。該體芯單元數(shù)值模擬結(jié)果數(shù)據(jù)表明，其他條件不變時，結(jié)構(gòu)優(yōu)化前的水側(cè)進(jìn)出口溫差為0.743℃，散熱功率為23.7W，結(jié)構(gòu)優(yōu)化后的水側(cè)進(jìn)出口溫差為0.917℃，散熱功率為29.3W，散熱性能提高23.6%。

3.3 整體散熱器的數(shù)值模擬分析結(jié)果

（1）圖8為散熱器的穩(wěn)態(tài)分析的溫度場切片圖，從圖中可分析得出，由于中間的空氣擾流情況復(fù)雜，空氣經(jīng)散熱器換熱后中間部分的空氣溫度明顯高于兩側(cè)的溫度，從中間向兩側(cè)的溫度逐漸降低。在整個空氣流域中溫度的等值面應(yīng)為四面雙曲線錐形結(jié)構(gòu)，散熱器前側(cè)的溫度低于后側(cè)的溫度。

（2）圖9為仿真分析的流場分布切片圖，從圖中可直觀得出空氣流經(jīng)散熱器時的流動情況。空氣流經(jīng)散熱器管帶翅片時受到阻力作用而減小，同時空氣的流動受到強(qiáng)烈干擾，增強(qiáng)了換熱效果。

（3）進(jìn)水流量對散熱器散熱性能的影響。在定風(fēng)量的情況下，選擇不同水流量對模型進(jìn)行仿真分析，設(shè)定進(jìn)水溫度為70℃恒溫進(jìn)水，進(jìn)水流量和模擬結(jié)果見表3。

將表3中數(shù)據(jù)整理到坐標(biāo)圖中，如圖10為散熱器的進(jìn)水流量和散熱量的關(guān)系，隨著水流量的增加，散熱器散熱量先增后減，在進(jìn)水流量為1.0L/min時有極值，表明該散熱器在進(jìn)水流量約為1.0L/min（溫度為70℃）的散熱性能為最佳工況值。圖中顯示了結(jié)果的模擬值和實驗值的對比關(guān)系，兩者的整體變化情況基本一致，說明數(shù)值模擬與實驗的吻合性。

4 結(jié)語

通過利用simulation CFD軟件對管帶式散熱器流場和溫度場進(jìn)行了數(shù)值模擬計算，并對散熱器的結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計，得到了散熱單元流場的內(nèi)部流動細(xì)節(jié)與散熱器整體的換熱分析計算數(shù)據(jù)，最后和實驗結(jié)果對比，各項數(shù)據(jù)偏差范圍合理，驗證計算方法可行。

總結(jié)全文：

（1）散熱器散熱管帶的百葉窗開角在25±2°區(qū)間的散熱效果最佳;

（2）調(diào)整百葉窗結(jié)構(gòu)的開窗尺寸A和B，可提高散熱性能;

（3）當(dāng)選擇進(jìn)水流量約為1.0L/min（水溫為70℃）時，散熱器的散熱量為最大值。

參考文獻(xiàn)：

[1]李世偉，吳國榮，袁津震，等.百葉窗結(jié)構(gòu)參數(shù)對管帶式散熱器流動傳熱特性影響的研究[J].內(nèi)燃機(jī)與動力裝置，2012，1：28-31.

[2]徐國濤.臺式計算機(jī)水冷散熱系統(tǒng)優(yōu)化研究[D].北京：華北電力大學(xué)，2008.

[3]漆波，田曉虎，李隆鍵，等.車用百葉窗翅片式熱交換器空氣側(cè)性能的CFD研究[J].天津理工大學(xué)學(xué)報，2007，32（2）：64-66.

[4]李重?zé)?，崔淑琴，玄哲?管帶式散熱器優(yōu)化設(shè)計實例[J].汽車工程學(xué)報，1995，4：18-21.

[5]付衛(wèi)東，吳金星.管帶式換熱器動態(tài)調(diào)節(jié)性能研究[J].制冷與空調(diào)，2006，6（6）：17-19.

作者簡介：李卓喜（1993-），男，工學(xué)本科，助理工程師，研究方向：熱能與動力工程。