劉曉麗，席明智，謝琰

（內(nèi)蒙古工業(yè)大學(xué)能源與動力工程學(xué)院，內(nèi)蒙古呼和浩特010051）

新型礦用車管芯式散熱器阻力性能的數(shù)值模擬

劉曉麗，席明智，謝琰

（內(nèi)蒙古工業(yè)大學(xué)能源與動力工程學(xué)院，內(nèi)蒙古呼和浩特010051）

利用ANSYS Flotran模塊對不同速度下的管芯式散熱器的空氣流動速度和壓降進(jìn)行數(shù)值模擬，模擬值與試驗值吻合，表明利用ANSYS對管芯式散熱器阻力性能進(jìn)行模擬的方法是可行的。

管芯式散熱器阻力性能數(shù)值模擬

1 前言

隨著計算機(jī)的發(fā)展，對流動的數(shù)值計算成為可能，并成為有效的研究手段及現(xiàn)代設(shè)計工具。ANSYS中的Flotran模塊可以模擬不同類型、不同結(jié)構(gòu)的散熱器在不同環(huán)境溫度下的空氣流動阻力，避免了試驗繁雜、周期長、費(fèi)用高、個別數(shù)據(jù)難測量等問題。

管芯式散熱器維修性能佳，是近年來礦用汽車發(fā)動機(jī)散熱器首選的管芯型式[1]。為使其應(yīng)用更為廣泛，要對其性能進(jìn)行深入研究，而國內(nèi)外對散熱器的性能研究多集中于管片式和管帶式[2-4]，并且多為試驗研究，而對管芯式散熱器的研究尚不多見。

本文以空氣流動作為對象，利用Pro/Engineer軟件建模，ANSYS Flotran模塊對模型進(jìn)行分析，模擬新型管芯式散熱器空氣流動速度變化和空氣阻力。

2 試驗系統(tǒng)

試驗系統(tǒng)按文獻(xiàn)[5]，由內(nèi)蒙古工業(yè)大學(xué)能源與動力工程學(xué)院搭建而成。試驗臺架主要由風(fēng)機(jī)、水泵、風(fēng)筒、試件、循環(huán)水路、加熱水箱、測量儀器和控制裝置組成，以模擬散熱器的真實工作環(huán)境。系統(tǒng)中風(fēng)速測量采用畢托管，并配用最小刻度為1 mm的傾斜式微壓計，進(jìn)風(fēng)溫度和出風(fēng)溫度選用銅－康銅熱電偶進(jìn)行測量。四排管散熱器的參數(shù)如表1。

表1 散熱器的結(jié)構(gòu)參數(shù)

3 基本假設(shè)與數(shù)學(xué)模型

認(rèn)為空氣流動為定常流動，氣體為理想氣體，不可壓縮，物性參數(shù)取為常數(shù)。不考慮散熱片的散熱，不考慮冷卻水的溫度影響，認(rèn)為空氣速度在散熱器的迎風(fēng)面上均勻分布，忽略重力的影響。

計算在不同的流速下進(jìn)行，可根據(jù)散熱器入口處空氣流速來估算雷諾數(shù)。因所確定流動區(qū)域內(nèi)空氣流動為外流，空氣與實體的相對運(yùn)動對流體邊界的影響可以不計，雷諾數(shù)的計算公式如下：

或

其中，為流體密度，為流體的動力粘度，

為流體的運(yùn)動粘度，為流場的特征速度，即進(jìn)口風(fēng)速；為流場的特征長度，取氣流方向上一段翅片的長度。

實際計算時，由于雷諾數(shù)大于4 000，故為湍流流動，計算采用標(biāo)準(zhǔn)湍流模型。根據(jù)文獻(xiàn)[6]，其質(zhì)量守恒方程和動量守恒方程分別如下：

補(bǔ)充湍流動能方程和湍流動能耗散率方程如下：

其中，為湍動粘度，為流體密度，是由于平均速度梯度引起的湍流動能的產(chǎn)生項，是由于浮力引起的湍流動能的產(chǎn)生項，代表可壓湍流中脈動擴(kuò)張的貢獻(xiàn)，和為經(jīng)驗常數(shù)，和分別是與湍流動能k和耗散率對應(yīng)的Prandtl數(shù)，和是用戶定義的源項。

4 管芯式散熱器阻力性能分析

4.1 散熱器幾何模型的建立

管芯式散熱器每個冷卻水管及其上面的散熱帶如同一個散熱單元，散熱器芯子型式如圖1所示。根據(jù)散熱器的幾何結(jié)構(gòu)，空氣沿軸負(fù)向流動，冷卻水沿y軸負(fù)向流動，水管和翅片沿軸錯排布置。管芯式散熱器局部幾何模型示意圖如圖2所示。

圖1 管芯式芯子

圖2 管芯式散熱器局部幾何模型示意圖

4.2 計算區(qū)域的確定及邊界條件的設(shè)置

為節(jié)省計算機(jī)資源，根據(jù)散熱器的實體結(jié)構(gòu)，確定計算區(qū)域如下：沿向（空氣流動方向）取四排管散熱器芯體的總長；沿向（水流動的方向）取三段翅片的長度；沿向取能夠涵蓋散熱器水管錯列位置變化為一個周期的區(qū)域。將散熱器實體部分減掉后便得到計算散熱器空氣阻力的空氣流動區(qū)域，即計算區(qū)域。以此為周期即可得到流過管排的空氣的壓降。

散熱器計算區(qū)域及邊界設(shè)置[7]，如圖3所示。固體壁面采用無滑移邊界，入口邊界為速度，出口邊界壓力值設(shè)置為0。

圖3 計算區(qū)域及邊界設(shè)置示意圖

4.3 計算模型的建立及有限元網(wǎng)格的劃分

采用三維建模軟件Pro/E建立幾何計算模型，并直接生成ANSYS模型，選定Flotran142單元進(jìn)行分析。對幾何模型（圖4）采用MeshTool進(jìn)行自由網(wǎng)格劃分，單元尺寸設(shè)置為2 mm，并進(jìn)行智能化控制。網(wǎng)格劃分的后節(jié)點(diǎn)數(shù)為111 174，單元數(shù)為503 332，網(wǎng)格劃分后的模型如圖5所示。

圖4 Pro/E幾何模型

圖5 ANSYS中劃分網(wǎng)格結(jié)果圖

4.4 加載、Flotran參數(shù)設(shè)置及求解

入口處施加速度和溫度邊界，出口處為壓力邊界。試驗中最大風(fēng)速為25 m/s，在最大風(fēng)量范圍內(nèi)選取5個工況點(diǎn)，風(fēng)速分別為7.5 m/s、12.6 m/s、17.9 m/s、21.9 m/s和23.3 m/s，出口相對壓力設(shè)置為0。翅片結(jié)構(gòu)共316個面，空氣接觸翅片和管子壁面施加無滑移邊界，水管排與排之間的空氣流通面采用對稱邊界，其余兩面施加周期性邊界。

在Fluid Properties對話框中選擇【air-SI】選項。設(shè)置參考壓力為101.35 kPa，空氣的物性參數(shù)取為定性溫度下的物性參數(shù)，其定性溫度取進(jìn)出口的平均溫度。選擇【Run Flotran】，進(jìn)行求解。通過模擬求出管芯式散熱器在不同速度下的風(fēng)阻值，及速度變化情況（圖6），空氣的流動跡線如圖7，圖8為湍流動能分布，動能最大值出現(xiàn)在翅片周圍區(qū)。下列圖示是速度為12.6 m/s的云圖顯示。

圖6 速度變化

圖7 空氣流動跡線圖

圖8 湍流動能分布

4.5 試驗結(jié)果與模擬結(jié)果對比

風(fēng)阻值的模擬結(jié)果與試驗數(shù)據(jù)擬合曲線對比如圖9所示，二者的最大相對誤差為6.88%，吻合度較好。

圖9 風(fēng)阻特性試驗值與計算值比較

5 結(jié)論

（1）能量損失主要在翅片周圍區(qū)域，隨著入口速度變大，即風(fēng)量增加，風(fēng)阻隨之增加。

（2）根據(jù)散熱器的幾何結(jié)構(gòu)，以散熱器中的空氣流體作為分析對象來實現(xiàn)對散熱器空氣阻力性能的分析，利用周期性邊界對計算區(qū)域進(jìn)行簡化，只取一個周期的空氣區(qū)域為計算區(qū)域，簡化了計算過程，節(jié)省計算資源。

（3）數(shù)值模擬結(jié)果與試驗數(shù)據(jù)的最大相對誤差小于7%，吻合較好，表明所建模型及算法對管芯式散熱器的風(fēng)阻分析具有實用價值。

1楊家騏.汽車散熱器（第一版）[M].人民交通出版社，1982.

2朱恂廖強(qiáng).管帶式汽車散熱器流動阻力與傳熱性能分析[J].重慶大學(xué)學(xué)報，2002(8).

3李岳林，張志沛，郭曉汾.管片式散熱器在汽車發(fā)動機(jī)冷卻系匹配設(shè)計中的數(shù)值模擬[J].汽車工程，2001(1).

4 Aoki H.Shinagawa T.Suga K.An Experimental study of the Local Heat Transfer Characteristics in Automotive Louvered Fin[J].Experimental Thermal and Fluid Science 1989,2(3):293-300.

5中華人民共和國第一機(jī)械工業(yè)部.汽車、拖拉機(jī)風(fēng)筒試驗方法[S].JB2293-1978，1978.

6王福軍.計算流體動力學(xué)分析[M].北京：清華大學(xué)出版社，2004.

7索文超，畢小平，李賀佳.車用散熱器空氣流動阻力預(yù)測研究[J].汽車工程，2008，30(9).

Numerical Simulation of Resistance Performance for New-Style Tube-Core Radiator Used in Mine Truck

Liu Xiaoli, Xi Mingzhi, Xie Yan
（College of Energy and Power Engineering, Inner Mongolia University of Technology, Hohhot 010051, China)

By using ANSYS Flotran module，numerical simulation at various velocity and pressure drop is conducted. The simulation result is in good agreement with the experimental data, which shows ANSYS can be used in resistance performance simulation of tube-core radiator.

tube-core radiator, resistance performance, numerical simulation

10.3969/j.issn.1671-0614.2010.02.006

來稿日期：2009-12-01

劉曉麗（1983－），女，助教，在讀碩士，主要研究方向為動力機(jī)械性能分析和研究。