電子設(shè)備封閉腔內(nèi)自然對(duì)流冷卻效果數(shù)值分析

王 燁，趙興杰，藺虎相，宋榮飛，管?chē)?guó)祥

?

電子設(shè)備封閉腔內(nèi)自然對(duì)流冷卻效果數(shù)值分析

王 燁1,2，趙興杰1，藺虎相1，宋榮飛1，管?chē)?guó)祥1

（1. 蘭州交通大學(xué)環(huán)境與市政工程學(xué)院，蘭州 730070；2. 蘭州交通大學(xué)鐵道車(chē)輛熱工教育部重點(diǎn)實(shí)驗(yàn)室，蘭州 730070）

為了研究安裝于封閉空間內(nèi)的電子設(shè)備散熱元件屬性及空間位置對(duì)腔內(nèi)自然對(duì)流傳熱特性的影響，該文采用FLUENT14.5軟件中的RNG-湍流模型對(duì)流體為空氣、高寬比為1的封閉腔內(nèi)溫度場(chǎng)、流場(chǎng)、壁面?zhèn)鳠崮芰M(jìn)行了數(shù)值分析。結(jié)果表明：在熱壁面1/3高度處布置1個(gè)導(dǎo)熱翅片時(shí)熱壁面的平均數(shù)比相同位置布置絕熱翅片時(shí)提高了9.67%；在熱壁面1/3高度處、冷壁面2/3高度處同時(shí)各布置1個(gè)導(dǎo)熱翅片時(shí)熱壁面平均數(shù)可取得雙翅片工況的最大值39.94，比單翅片的最優(yōu)工況平均數(shù)提高了14.34%。本文研究結(jié)果對(duì)于改善工農(nóng)業(yè)工程中電子元器件的自然對(duì)流冷卻效果、優(yōu)化散熱元器件的空間布置具有一定的理論指導(dǎo)意義。

數(shù)值模擬；溫度；對(duì)流換熱；封閉腔；翅片

0 引 言

封閉腔內(nèi)的自然對(duì)流現(xiàn)象廣泛存在于諸多工業(yè)生產(chǎn)系統(tǒng)中，比如太陽(yáng)能集熱器系統(tǒng)、電子設(shè)備的冷卻、建筑熱工設(shè)計(jì)、核反應(yīng)堆設(shè)計(jì)等。在農(nóng)業(yè)大棚內(nèi)發(fā)生的自然對(duì)流熱質(zhì)交換過(guò)程對(duì)于室內(nèi)溫濕度、氣流速度的控制與產(chǎn)品數(shù)量及質(zhì)量密切相關(guān)。因此，對(duì)封閉空間內(nèi)自然對(duì)流機(jī)理的研究對(duì)于改善這些系統(tǒng)的綜合性能具有重要的工程實(shí)際意義。ButChelor[1]以雙層窗為研究對(duì)象，首次探討了不同瑞利數(shù)下雙層窗玻璃夾層中的自然對(duì)流換熱機(jī)理。后來(lái)人們又從腔體高寬比[2]、瑞利數(shù)[3-8]、普朗特?cái)?shù)[8-10]、腔體傾角[11]等方面進(jìn)行了卓有成效的探索。

在強(qiáng)化封閉腔內(nèi)自然對(duì)流傳熱措施方面，有在腔內(nèi)側(cè)壁面施加正弦變壁溫條件[12]，或在腔內(nèi)設(shè)置一個(gè)旋轉(zhuǎn)圓柱體[13]以及在腔體內(nèi)側(cè)壁上布置水平薄翅片等[14]。Ampofo[15]用實(shí)驗(yàn)方法揭示了內(nèi)置翅片對(duì)封閉腔內(nèi)自然對(duì)流傳熱特性的影響機(jī)理，發(fā)現(xiàn)內(nèi)置翅片強(qiáng)化了熱壁面的傳熱能力，冷壁面邊界層上的速度場(chǎng)、溫度場(chǎng)及湍流參量的峰值均有了顯著提高。目前，內(nèi)置翅片對(duì)封閉腔內(nèi)自然對(duì)流傳熱特性的影響研究重點(diǎn)主要涉及翅片長(zhǎng)度[16-19]、翅片傾斜角[19]、翅片形狀[20-21]、翅片材質(zhì)[22-25]等方面。

關(guān)于翅片材質(zhì)對(duì)封閉腔內(nèi)傳熱特性的研究比較一致的結(jié)論是：封閉腔內(nèi)傳熱效率隨翅片導(dǎo)熱系數(shù)增大而增大。而Nag等[22]的研究發(fā)現(xiàn)：在豎向熱壁面上布置的導(dǎo)熱翅片對(duì)冷壁面的傳熱有強(qiáng)化作用，在熱壁面上布置的絕熱翅片反而弱化了冷壁面的傳熱。Frederick等[23]對(duì)于有內(nèi)置翅片的封閉腔內(nèi)自然對(duì)流換熱三維數(shù)值研究表明：位于豎向熱壁面上的翅片導(dǎo)熱性越好，翅片的阻擋作用越弱，對(duì)腔內(nèi)的傳熱強(qiáng)化效果越顯著。Xu等[24]對(duì)比研究了封閉腔內(nèi)湍流自然對(duì)流中側(cè)壁所布置單一導(dǎo)熱翅片和絕熱翅片對(duì)流場(chǎng)結(jié)構(gòu)和傳熱效率的影響，但沒(méi)有給出翅片的優(yōu)化配置方案。Elatar等[25]研究了封閉腔內(nèi)層流狀態(tài)下位于側(cè)壁高溫面上的翅片導(dǎo)熱性對(duì)翅片效率和對(duì)流傳熱效率的提高程度，得到了不同翅片導(dǎo)熱率情況下壁面平均努塞爾特?cái)?shù)與數(shù)和翅片高度之間的關(guān)系式。這類(lèi)對(duì)于單一翅片的研究成果已廣泛應(yīng)用于工業(yè)生產(chǎn)中[26-28]。但這些研究只是針對(duì)單一因素進(jìn)行探討的，關(guān)于有內(nèi)熱源情況下封閉腔內(nèi)側(cè)壁翅片位置和數(shù)量對(duì)湍流自然對(duì)流傳熱特性共同影響的研究還未見(jiàn)報(bào)道。而不斷集成化的工業(yè)設(shè)備及系統(tǒng)優(yōu)化設(shè)計(jì)中，在同一封閉空間里往往是多熱源、多元件并存的，如何合理布置這些發(fā)熱元件并控制封閉空間內(nèi)的流動(dòng)冷卻過(guò)程，對(duì)于系統(tǒng)的安全運(yùn)行和可靠性非常重要?；谶@樣的工程背景，本文針對(duì)有內(nèi)熱源的封閉腔內(nèi)側(cè)壁翅片屬性、數(shù)量、位置等因素對(duì)腔內(nèi)自然對(duì)流冷卻的影響規(guī)律及溫度場(chǎng)與流場(chǎng)的耦合作用機(jī)理進(jìn)行了研究，以獲得冷卻性能最優(yōu)的最佳翅片數(shù)量和位置組合，從而為農(nóng)業(yè)及工業(yè)工程中類(lèi)似于封閉空間中電子元件的優(yōu)化布置并提高其使用壽命提供理論參考。

1 封閉腔物理模型和數(shù)學(xué)模型

1.1 物理模型

本文所研究的封閉腔橫斷面結(jié)構(gòu)如圖1所示。以文獻(xiàn)[29]的研究模型為例，寬、高、深分別取=0.75 m、=0.75 m、=1.5 m。

注：Th為熱壁面溫度，℃；Tc為冷壁面溫度，℃；H為腔體高度，m；L為腔體寬度，m；a為熱壁面上翅片與腔體底面之間的距離，m；b為冷壁面上翅片與腔體底面之間的距離，m；c為熱源頂面與腔體底面之間的距離，m。

根據(jù)文獻(xiàn)[29]，當(dāng)封閉腔深寬比大于1.8時(shí)可以不考慮腔體內(nèi)流動(dòng)的三維效應(yīng)。所以本文采用二維數(shù)值模擬。腔體內(nèi)為空氣，流體密度變化采用Boussinesq假設(shè)，=1.58×109。

翅片高度均為80 mm，翅片厚度為4 mm。自然對(duì)流換熱中，流場(chǎng)與溫度場(chǎng)之間存在強(qiáng)烈的耦合關(guān)系，不同翅片位置組合會(huì)對(duì)封閉腔內(nèi)流場(chǎng)和溫度場(chǎng)分析產(chǎn)生不同影響，從而影響自然對(duì)流換熱性能。為了揭示這一特性，本文設(shè)計(jì)了如表1所示的計(jì)算模型。定義無(wú)量綱長(zhǎng)度：=/，=/，，分別為水平和鉛錘方向的實(shí)際距離，m；=0.75 m。

表1 物理模型的尺寸設(shè)置

1.2 控制方程數(shù)學(xué)模型

由于RNG模型在計(jì)算湍流動(dòng)能耗散率方程中產(chǎn)生項(xiàng)系數(shù)時(shí)引入了主流的時(shí)均應(yīng)變率，使得產(chǎn)生項(xiàng)系數(shù)既與流動(dòng)情況有關(guān)，也是空間坐標(biāo)的函數(shù)，從而有效地改善了計(jì)算精度；另一方面，RNG模型也考慮了湍流渦旋結(jié)構(gòu)，而本文模型中氣流通過(guò)翅片的擾流會(huì)產(chǎn)生顯著的渦旋結(jié)構(gòu)。因此，本文采用RNG模型更適合于準(zhǔn)確捕捉翅片與側(cè)壁結(jié)合部位以及腔體內(nèi)部其余位置形成的渦旋結(jié)構(gòu)，從而能更準(zhǔn)確地描述所研究腔體中的流動(dòng)與傳熱過(guò)程。本文計(jì)算中常數(shù)項(xiàng)1ε=1.42，2ε=1.68。

1.3 邊界條件和初始條件

參考文獻(xiàn)[29]，熱壁面溫度為323.15 K，冷壁面溫度為283.15 K，熱源溫度為353.15 K。參考文獻(xiàn)[30]，頂面和底面向外散熱的熱流密度均為-7.0 w/m2，封閉腔內(nèi)初始溫度取冷、熱壁面溫度的算術(shù)平均值，為303.15 K，所有氣固交界面為速度無(wú)滑移條件。

2 數(shù)值求解方法

2.1 數(shù)學(xué)模型驗(yàn)證

為了保證數(shù)學(xué)模型以及數(shù)值方法的正確性，需要進(jìn)行模型驗(yàn)證。本文在=1.58×109條件下與Tian[29]的試驗(yàn)數(shù)據(jù)進(jìn)行比對(duì)。選擇表1中的物理模型1，網(wǎng)格數(shù)為260×260，時(shí)間步長(zhǎng)為0.01 s，邊界條件與文獻(xiàn)[29]的相同。氣固交界面取速度無(wú)滑移條件。采用RNG-模型對(duì)封閉腔內(nèi)空氣自然對(duì)流與傳熱過(guò)程進(jìn)行數(shù)值模擬，得到腔體側(cè)壁距底面無(wú)量綱長(zhǎng)度=0.5處水平線上的豎向速度值、溫度值與文獻(xiàn)[29]的試驗(yàn)數(shù)據(jù)比對(duì)結(jié)果，如圖2所示。

圖2 本文數(shù)值結(jié)果與文獻(xiàn)[29]試驗(yàn)結(jié)果比對(duì)（距底面無(wú)量綱長(zhǎng)度Y=0.5） Fig. 2 Comparisons between numerical results in this paper and experimental results in reference [29] (non-dimensional length from bottom Y=0.5)

本文模擬結(jié)果的速度值與文獻(xiàn)[29]試驗(yàn)結(jié)果的相對(duì)誤差最大值為4.8%，最小值為0；本文模擬結(jié)果的溫度值與文獻(xiàn)[29]試驗(yàn)結(jié)果的相對(duì)誤差最大值為3.2%，最小值為0。由于近壁區(qū)速度梯度和溫度梯度均很大，導(dǎo)致了較大的測(cè)試誤差，從而使得數(shù)值結(jié)果與試驗(yàn)結(jié)果間存在較大的相對(duì)偏差，但速度及溫度的相對(duì)偏差值已經(jīng)足以證明所采用數(shù)學(xué)模型的正確性和可靠性[10]。模型和數(shù)值方法可用于封閉腔內(nèi)空氣湍流自然對(duì)流與傳熱特性的數(shù)值分析。

2.2 迭代時(shí)間步長(zhǎng)

根據(jù)表1中的物理模型1，分別采用0.02、0.05和0.10 s的時(shí)間步長(zhǎng)進(jìn)行數(shù)值試驗(yàn)，將3種時(shí)間步長(zhǎng)得到的熱壁面距底面無(wú)量綱長(zhǎng)度=0.5處水平線上的豎向速度及溫度分別進(jìn)行對(duì)比。3種時(shí)間步長(zhǎng)對(duì)速度的計(jì)算結(jié)果的最大相對(duì)偏差為0.09%，最小相對(duì)偏差為0。3種時(shí)間步長(zhǎng)對(duì)溫度的計(jì)算結(jié)果的最大相對(duì)偏差為0.02%，最小相對(duì)偏差為0。不同時(shí)間步長(zhǎng)計(jì)算結(jié)果的差異滿足工程計(jì)算精度要求[10]，在保證計(jì)算精度的前提下，為了節(jié)約計(jì)算時(shí)間并確保后續(xù)工況計(jì)算收斂，本文選取時(shí)間步長(zhǎng)為0.05 s。

2.3 網(wǎng)格獨(dú)立性驗(yàn)證

對(duì)于只在熱壁面上布置1個(gè)翅片的情況，選取翅片位于熱壁面距底面以上1/6高度處的物理模型，在相同邊界條件下，采用3套網(wǎng)格（41 358、69 154、98 124）進(jìn)行試算，3套網(wǎng)格所得=0.5處外法線上的豎向速度最大相對(duì)偏差為0.03%、最小相對(duì)偏差為0；溫度最大相對(duì)偏差為0.01%、最小相對(duì)偏差為0。在保證計(jì)算精度的前提下[10]，為了節(jié)約計(jì)算時(shí)間，本文選取69 154為后續(xù)計(jì)算的網(wǎng)格數(shù)。

對(duì)于在熱、冷壁面上同時(shí)布置翅片的情況，選取熱壁面翅片位于距底面以上1/3H處、冷壁面翅片位于距底面以上1/3H處的物理模型為計(jì)算模型，在相同邊界條件下，采用3套網(wǎng)格（45 725、68 067、90 822）進(jìn)行試算，3套網(wǎng)格所得=0.5處外法線上的豎向速度最大相對(duì)偏差為0.8%、最小相對(duì)偏差為0；溫度最大相對(duì)偏差為0.4%、最小相對(duì)偏差為0。在保證計(jì)算精度的前提下[10]，確定計(jì)算網(wǎng)格數(shù)為68 067。

據(jù)此確定適合于本文研究表1中模型1的網(wǎng)格數(shù)為67 081，模型2～模型6的網(wǎng)格數(shù)為69 154，模型7和模型10的網(wǎng)格數(shù)為68 067，模型8和模型9的網(wǎng)格數(shù)為67 967。

2.4 求解計(jì)算

離散格式為：壓力項(xiàng)采用Green-Gauss Cell Based 格式離散[31]，動(dòng)量方程、能量方程、湍流動(dòng)能方程及湍流動(dòng)能耗散率方程均采用二階迎風(fēng)格式離散，采用SIMPLE算法求解速度和壓力耦合問(wèn)題。松弛因子設(shè)置[31]：壓力為0.3，動(dòng)量為0.7，湍流動(dòng)能和湍流動(dòng)能耗散率均為0.8，能量方程、湍流黏度、體積力和密度均為1。湍流模型選用RNG-模型。求解控制方程時(shí)的收斂條件與文獻(xiàn)[31]相同。

3 計(jì)算結(jié)果與分析

3.1 翅片材質(zhì)對(duì)腔內(nèi)對(duì)流和傳熱的影響

內(nèi)置熱源的封閉腔內(nèi)自然對(duì)流在房間空調(diào)冷卻、電力設(shè)備冷卻等方面應(yīng)用廣泛，但翅片是否導(dǎo)熱對(duì)冷卻效果的影響不同。將翅片分別布置在熱壁面距底面以上1/6、1/3、1/2、2/3、5/6高度處，分析導(dǎo)熱和絕熱翅片位置與內(nèi)置熱源對(duì)封閉腔內(nèi)對(duì)流和傳熱的綜合影響。

絕熱翅片表面熱流密度為0。導(dǎo)熱翅片由于高度（沿壁面外法線方向的尺寸）小、導(dǎo)熱系數(shù)大（銅板=398 W/m·K），計(jì)算中近似認(rèn)為其表面溫度均勻且和熱壁面溫度相同，為323.15 K，熱源除了底面絕熱外，其余各表面溫度均為353.15 K。

圖3為絕熱翅片位于熱壁面不同位置時(shí)沿?zé)岜诿娓叨确较蚓植繑?shù)的變化曲線?？梢钥闯觯?0.1范圍內(nèi)熱壁面局部數(shù)幾乎不受翅片位置影響。但隨著邊界層的不斷發(fā)展，局部數(shù)在翅片的下表面和上表面發(fā)生了階躍性變化，最大值達(dá)到了75，最小值接近0。這是因?yàn)槌崞趥?cè)壁上起到了渦旋發(fā)生器的作用，浮升力作用下不斷上升的氣流首先遇到翅片發(fā)生滯止，隨即繞過(guò)翅片并在翅片上方形成渦旋，從而強(qiáng)化了翅片基部的傳熱過(guò)程。熱壁面局部數(shù)沿流動(dòng)方向因溫度梯度逐漸減小而呈下降趨勢(shì)。所得熱壁面絕熱翅片肋基下部數(shù)趨于0的結(jié)果與文獻(xiàn)[24]的研究結(jié)果一致。

圖3 熱壁面絕熱翅片位置對(duì)熱壁面局部Nu數(shù)的影響

圖4為導(dǎo)熱翅片位于熱壁面不同位置時(shí)局部數(shù)沿?zé)岜诿娓叨确较虻淖兓€。與絕熱翅片情況類(lèi)似，沿邊界層發(fā)展方向局部數(shù)呈下降趨勢(shì)，<0.1范圍內(nèi)熱壁面的局部數(shù)幾乎不受翅片位置影響。這是因?yàn)闊岜诿娴撞繀^(qū)域?yàn)閷恿鳡顟B(tài)，具有較薄的熱層，溫度梯度大，傳熱能力就強(qiáng)。與絕熱翅片不同的是數(shù)只在翅片下方發(fā)生突變，翅片上方增幅較小。這是因?yàn)閷?dǎo)熱翅片下、上表面溫差遠(yuǎn)小于絕熱翅片的情況，上升氣流繞過(guò)翅片后形成的渦旋對(duì)翅片上部傳熱的貢獻(xiàn)較小。

導(dǎo)熱翅片對(duì)于熱壁面平均數(shù)的提高率影響見(jiàn)表 2。由表2可以看出，絕熱翅片位于1/6高度時(shí)熱壁面平均值最大，為31.85，導(dǎo)熱翅片位于1/3高度時(shí)熱壁面平均值最大，為34.93。導(dǎo)熱翅片的平均數(shù)較絕熱翅片提高了9.67%。所以，在封閉空間內(nèi)用于固定散熱元件的基礎(chǔ)盡可能設(shè)計(jì)為導(dǎo)熱系數(shù)較高的材料，并布置在熱壁面距底部約1/3位置處更有利于散熱元件的冷卻。

圖4 熱壁面導(dǎo)熱翅片位置對(duì)熱壁面局部Nu數(shù)影響

表2 翅片位置對(duì)熱壁面平均Nu數(shù)提高率的影響

3.2 翅片數(shù)量對(duì)腔內(nèi)對(duì)流和傳熱的影響

根據(jù)3.1節(jié)的分析結(jié)果，熱壁面布置導(dǎo)熱翅片在強(qiáng)化傳熱效果上優(yōu)于絕熱翅片，因此，下文僅針對(duì)導(dǎo)熱翅片分析翅片位置及數(shù)量對(duì)封閉腔內(nèi)自然對(duì)流傳熱的影響。

當(dāng)封閉腔底面有大功率電子元器件散熱時(shí)，兩側(cè)壁電子元器件的固定位置對(duì)系統(tǒng)的可靠性和使用壽命有一定的影響?；谶@一背景，本節(jié)在熱、冷壁面分別布置1個(gè)翅片，不同翅片位置組合工況設(shè)置如表3所示。其中熱壁面翅片保持恒溫323.15 K，冷壁面翅片保持恒溫283.15 K，熱源位于腔體底面中心，其表面溫度為353.15 K。

圖5為4種工況在準(zhǔn)定常階段時(shí)腔體內(nèi)部的等溫線圖。由圖5可知，腔內(nèi)熱壁面附近流體在受到熱壁加熱和熱源輻射的共同作用下沿?zé)岜诿嫦蛏狭鲃?dòng)，遇到翅片時(shí)原有熱邊界層和速度邊界層被破壞，隨即形成新的發(fā)展邊界層和渦旋結(jié)構(gòu)并強(qiáng)化了熱壁面與空氣之間的換熱過(guò)程。熱氣流在浮升力作用下繼續(xù)向上爬升，沿著腔體頂部向冷壁面一側(cè)運(yùn)動(dòng)，在腔體頂部與冷壁面交界處開(kāi)始形成沿冷壁面向下發(fā)展的邊界層，熱壁面翅片對(duì)熱壁面邊界層結(jié)構(gòu)的影響與冷壁面翅片對(duì)冷壁面邊界層結(jié)構(gòu)的影響產(chǎn)生了強(qiáng)烈的耦合作用。如圖 5a所示，當(dāng)熱壁面翅片距離底面較近時(shí)沿?zé)岜诿娴倪吔鐚幽茉谙掠蔚玫匠浞职l(fā)展，當(dāng)冷壁面翅片距離底面較近時(shí)沿冷壁面的邊界層在下游的發(fā)展得到了抑制。隨著冷壁面翅片上移，腔體頂部高溫區(qū)變小，底部附近溫度梯度變小，如圖5b所示。另外，從圖5c及5d也可看出，冷壁面翅片位置對(duì)溫度場(chǎng)的影響要比熱壁面翅片位置的影響顯著。

表3 冷熱壁面翅片相對(duì)位置

Table 3 Relative fin position on cool and hot walls

工況Working conditionab 11/3H1/3H 21/3H2/3H 32/3H1/3H 42/3H2/3H

圖5 冷熱壁面翅片位置對(duì)溫度場(chǎng)的影響

Fig. 5 Influence of fins position of hot and cold walls on temperature profile

圖6為準(zhǔn)定常階段時(shí)4種工況下的腔體內(nèi)部流函數(shù)圖。由圖6可知，兩側(cè)壁翅片位置的不同組合對(duì)空腔內(nèi)流動(dòng)結(jié)構(gòu)的影響存在較大差異。工況1和工況2中流體繞過(guò)熱壁面翅片時(shí)被加速，但工況1冷壁面邊界層起始段對(duì)流較弱，而工況2冷壁面邊界層下游段對(duì)流較強(qiáng)，這與文獻(xiàn)[32]的研究結(jié)論一致。由圖6c可知，冷壁面邊界層在下游得到了抑制，流體在熱壁面附近浮升力誘導(dǎo)下沿?zé)嵩错敳砍療岜诿娣较蜻\(yùn)動(dòng)。工況4中熱壁面邊界層上游和冷壁面邊界層下游均得到了充分發(fā)展，流體被加速，這一區(qū)段的傳熱被得到顯著強(qiáng)化。

a. 工況1a. Working condition 1b. 工況2b. Working condition 2c. 工況3c. Working condition 3d. 工況4d. Working condition 4

圖7為熱壁面1/3處布置翅片與熱壁面1/3、冷壁面2/3處同時(shí)布置翅片在有無(wú)內(nèi)熱源情況下沿封閉腔底面中心鉛錘線上的水平速度變化曲線。由于空腔側(cè)壁上翅片的存在，腔內(nèi)的對(duì)流速度分布在單翅片、雙翅片情況下表現(xiàn)出了較大的差異?？梢钥闯?，單翅片與雙翅片對(duì)空腔頂部附近水平速度的影響微弱，對(duì)腔體底部附近區(qū)域水平速度場(chǎng)結(jié)構(gòu)的擾動(dòng)激烈，且影響趨勢(shì)很接近，這一現(xiàn)象與無(wú)熱源時(shí)并不一致。翅片對(duì)腔體核心區(qū)水平速度的影響存在較大差異：?jiǎn)纬崞瑢?duì)于腔體核心區(qū)較低位置的擾動(dòng)較明顯，而雙翅片則對(duì)于上部區(qū)域的流場(chǎng)擾動(dòng)更為強(qiáng)烈。不同數(shù)量翅片通過(guò)對(duì)流場(chǎng)結(jié)構(gòu)的改變導(dǎo)致了封閉腔內(nèi)溫度場(chǎng)的不同分布特征。

圖7 翅片數(shù)量和位置對(duì)水平速度的影響

圖8a為不同工況下局部數(shù)沿?zé)岜诿娓叨确较虻淖兓€?？梢钥闯?，冷熱壁面導(dǎo)熱翅片位置對(duì)=0.1以下區(qū)域局部傳熱能力幾乎沒(méi)有影響，熱壁面底部區(qū)域局部值最大，這是因?yàn)闊岜诿嫠俣冗吔鐚蛹皽囟冗吔鐚泳鹗加跓岜诿娴撞浚鵁岜诿娴撞繉儆趯恿鳡顟B(tài)，存在很薄的熱層，從而使得這一區(qū)域表現(xiàn)出了較強(qiáng)的傳熱能力。工況1和工況2由于其熱壁面上翅片位置相同從而表現(xiàn)出類(lèi)似的局部傳熱特性，工況3和工況4的局部傳熱特性較一致也是源于這一因素。

圖8b為冷熱壁面導(dǎo)熱翅片不同位置時(shí)熱壁面平均數(shù)比較?？梢钥闯觯?dāng)熱壁面導(dǎo)熱翅片位于1/3高度處、冷壁面導(dǎo)熱翅片位于2/3高度處時(shí)，熱壁面平均數(shù)達(dá)到了39.94，壁面與腔體內(nèi)部的平均傳熱效率最高，腔體內(nèi)部的對(duì)流換熱最強(qiáng)。因此，在散熱空腔內(nèi)雙側(cè)布置翅片時(shí)采取工況2的布置方式可以最大程度地對(duì)散熱元件進(jìn)行自然對(duì)流冷卻降溫，以保證散熱元件的安全使用。

圖8 翅片位置對(duì)熱壁面局部Nu數(shù)及平均Nu數(shù)的影響

與表2中單翅片各工況所得熱壁面的平均數(shù)相比，圖8b所示雙翅片各工況的強(qiáng)化傳熱效果均有明顯改善。例如對(duì)于單導(dǎo)熱翅片位于1/3H的最優(yōu)工況及雙導(dǎo)熱翅片的最優(yōu)工況2，對(duì)應(yīng)的熱壁面平均數(shù)分別為34.93和39.94，雙翅片較單翅片提高了14.34%。這說(shuō)明在封閉腔內(nèi)利用自然對(duì)流中流場(chǎng)與溫度場(chǎng)間的耦合關(guān)系，采用在熱、冷壁面上合理位置處同時(shí)布置翅片可以有效提高換熱效率，對(duì)于存在散熱元件冷卻問(wèn)題的封閉腔內(nèi)結(jié)構(gòu)設(shè)計(jì)具有重要的理論指導(dǎo)意義。

4 結(jié) 論

通過(guò)在有內(nèi)熱源的封閉方腔內(nèi)側(cè)壁上布置導(dǎo)熱翅片或絕熱翅片，數(shù)值分析了不同翅片位置、數(shù)量及不同翅片材質(zhì)對(duì)封閉腔內(nèi)湍流自然對(duì)流傳熱特性的影響。得到了如下主要結(jié)論：

1）在封閉腔側(cè)壁布置翅片，可以通過(guò)改變速度邊界層與溫度邊界層結(jié)構(gòu)來(lái)強(qiáng)化封閉腔內(nèi)對(duì)流換熱過(guò)程。在單側(cè)壁面布置單個(gè)翅片與在兩側(cè)壁面同時(shí)布置翅片對(duì)于封閉腔側(cè)壁傳熱能力的提升幅度存在差異。

2）絕熱翅片位于熱壁面距底面以上1/6高度處時(shí)熱壁面的平均傳熱能力最強(qiáng)。導(dǎo)熱翅片位于熱壁面距底面以上1/3高度處時(shí)熱壁面的平均傳熱能力最強(qiáng)。熱壁面翅片位于熱壁面距底面以上1/3高度與冷壁面翅片位于冷壁面距底面以上2/3高度組合使熱壁面的平均傳熱能力最強(qiáng)。

3）在熱壁面距底面以上1/3高度處布置導(dǎo)熱翅片比在同位置布置絕熱翅片時(shí)熱壁面平均數(shù)提高了9.67%；在熱壁面距底面以上1/3高度處布置導(dǎo)熱翅片、冷壁面距底面以上2/3處布置導(dǎo)熱翅片時(shí)熱壁面平均數(shù)比只在熱壁面距底面以上1/3處布置導(dǎo)熱翅片的平均數(shù)提高了14.34%。

封閉腔內(nèi)側(cè)壁布置翅片在強(qiáng)化對(duì)流傳熱的同時(shí)，翅片附近漩渦的存在會(huì)導(dǎo)致局部熱滯留現(xiàn)象，在實(shí)際應(yīng)用中對(duì)應(yīng)電子元器件與母材的連接處局部散熱效果下降，從而降低元器件的使用可靠性。因此，針對(duì)特定的封閉空間內(nèi)各表面不同發(fā)射率時(shí)如何確定翅片結(jié)構(gòu)參數(shù)、安裝方式等還需要進(jìn)一步的探索。

[1] Bachelor G K. Heat transfer by free convection across a closed cavity between vertical boundaries at different temperatures[J]. Quarterly of Applied Mathematics, 1954, 12(3): 209－233.

[2] 王迪. 物理幾何條件對(duì)多熱源小空間自然對(duì)流換熱的影響[D].長(zhǎng)春：吉林建筑大學(xué)，2016. Wang Di. Influence of Physical Geometry on Natural Convection Heat Transfer in a Small Space with Multiple Heat Sources[D]. Changchun: Jilin Jianzhu University, 2016. (in Chinese with English abstract)

[3] Cheikh N B, Beya B B, Lili T. Influence of thermal boundary conditions on natural convection in a square enclosure partially heated from below[J]. International communications in Heat and Mass Transfer, 2007, 34(3): 369－379.

[4] 李娜，過(guò)增元，李志信，方腔自然對(duì)流中力的尺度效應(yīng)[J]. 清華大學(xué)學(xué)報(bào)，2002，42(11)：1508－1514. Li Na, Guo Zengyuan, Li Zhixin. Size effect of forces on natural convection in a square cavity[J]. Journal of Tsinghua University, 2002, 42(11): 1508－1514. (in Chinese with English abstract)

[5] Sharma A K, Velusamy K, Balaji C. Turbulent natural convection in an enclosure with localized heating from below[J]. International Journal of Thermal Sciences, 2007, 46(12): 1232－1241.

[6] Paroncini M, Corvaro F, Montucchiari A, et al. A numerical and experimental analysis on natural convective heat transfer in a square enclosure with partially active side walls[J]. Experimental Thermal and Fluid Science, 2012(36): 118－125.

[7] Xu Feng, Saha S C. Transition to an unsteady flow induced by a fin on the sidewall of a differentially heated air-filled square cavity and heat transfer[J]. International Journal of Heat and Mass Transfer, 2014(71): 236－244.

[8] Xia Keqing, Lam S, Zhou Shengqi. Heat-flux measurement in high-Prandtl-number turbulent Rayleigh-Bénard convection [J]. Physical Review Letters, 2002, 88(6): 1501-1504.

[9] 王燁，王良壁. 一種用于分析封閉腔內(nèi)湍流自然對(duì)流換熱的新模型[J]. 應(yīng)用力學(xué)學(xué)報(bào)，2014，31(5)：814－818. Wang Ye, Wang Liangbi. A new-model for the turbulent natural convection heat transfer in enclosure[J]. Chinese Journal of Applied Mechanics, 2014, 31(5): 814－818. (in Chinese with English abstract)

[10] 陽(yáng)祥，陶文銓. 高瑞利數(shù)下封閉腔內(nèi)自然對(duì)流的數(shù)值模擬[J].西安交通大學(xué)學(xué)報(bào)，2014，48(5)：27－31. Yang Xiang, Tao Wenquan. Numerical simulations for natural convection with high Rayleigh number in a tall rectangular cavity[J]. Journal of Xi’an Jiaotong University, 2014, 48(5): 27－31. (in Chinese with English abstract)

[11] Heiland H G, Rubes D, Wozniak G. Convection experiments in an inclined narrow cavity[J]. Heat Mass Transfer, 2012, 48(12): 2007－2012.

[12] Kwak H S, Kuwahara K, Hyun J M. Resonant enhancement of natural convection heat transfer in a square enclosure[J].International Journal of Heat and Mass Transfer, 1998, 41(18): 2837－2846.

[13] Fu Wushung, Cheng Chaosheng, Shieh W J. Enhancement of natural convection heat transfer of an enclosure by a rotating circular cylinder[J]. International Journal of Heat and Mass Transfer, 1994, 37(13): 1885－1897.

[14] Xu Feng, Patterson J C, Lei Chengwang. Temperature oscillations in a differentially heated cavity with and without a fin on the sidewall[J]. International Communications in Heat and Mass Transfer, 2010, 37(4): 350－359.

[15] Ampofo F. Turbulent natural convection of air in a non-partitionedor partitioned cavity with differentially heated verticaland conducting horizontal walls[J]. Experimental Thermal and Fluid Science, 2005, 29(2): 137－157.

[16] Liu Yang, Lei Chengwang, Patterson J C. Plume separation from an adiabatic horizontal thinfin placed at different heights on the sidewall of a differentially heated cavity[J]. International Communications in Heat and Mass Transfer, 2015, 61: 162－169

[17] Jani S,Mahmoodi M,Amini M, et al. Numerical investigation of natural convection heat transfer in a symmetrically cooled square cavity with a thin fin on its bottom wall[J]. Thermal Science, 2014, 18(4): 1119－1132.

[18] Xu Feng, Patterson J C, Lei Chengwang. Effect of the fin length on natural convection flow transition in a cavity[J]. International Journal of Thermal Sciences,2013, 70: 92－101.

[19] Khanafer K, AlAmiri A, Bull J. Laminar natural convection heat transfer in a differentially heated cavity with a thin porous fin attached to the hot wall[J]. International Journal of Heat and Mass Transfer, 2015, 87: 59－70.

[20] Lorenzini G, Machado B S, Isoldi L A, et al. Structure design of rectangular fin intruded into mixed convective lid-driven cavity flows[J]. Journal of Heat Transfer, 2016, 138(10): 2501-2512.

[21] 黃格永. 受限空間電子元件自然對(duì)流散熱特性研究[D]. 重慶：重慶大學(xué)，2016. Huang Geyong. Characteristics of Natural Convection Heat Dissipation of Electronic Components in Confined Space[D]. Chongqing: Chongqing University, 2016. (in Chinese with English abstract)

[22] Nag A, Sarkar A, Sastri V M K. Natural convection in a differentially heated square cavity with horizontal partition plate on the hot wall[J]. Computer Methods in Applied Mechanics and Engineering, 1993, 110(1-2): 143－156.

[23] Frederick R L, Valencia A. Heat transfer in a square cavity with a conducting partition on its hot wall[J]. International communications in Heat and Transfer, 1989, 16(3): 347－354.

[24] Xu Feng, Patterson J C, Lei Chengwang. Unsteady flow and heat transfer adjacent to the sidewall of a differentially heated cavity with a conducting and an adiabatic fin[J]. International Journal of Heat and Fluid Flow, 2011, 32(3): 680－687.

[25] Elatar A, Teamah M A, Hassab M A. Numerical study of laminar natural convection inside squareenclosure with single horizontal fin[J].International Journal of Thermal Sciences, 2016 (99): 41－51.

[26] Peng Hao, Ling Xiang, Li Juan. Numerical simulation and experimental verification onthermal performance of a novel fin-plate thermosyphon[J]. Applied Thermal Engineering, 2012(40): 181－188.

[27] Frederick R L.Heat transfer enhancement in cubical enclosures with verticalfins[J]. Applied Thermal Engineering,2007, 27(8-9): 1585－1592.

[28] Varol Y, Oztop H F.Control of buoyancy-induced temperature and flow fieldswith an embedded adiabatic thin plate in porous triangular cavities[J]. Applied Thermal Engineering,2009, 29(2-3): 558－566.

[29] Tian Y S, Karayiannis T G. Low turbulence natural convection in an air filled square cavity-Part I: the thermal and fluid flow fields[J]. International Journal of Heat and Mass Transfer, 2000, 43 (6): 849－866.

[30] Omri M, Galanis N. Numerical analysis of turbulent buoyant flows in enclosures[J]. Influence of grid and boundary conditions, 2007(46): 727－738.

[31] 王燁，徐燃. 青藏線空調(diào)列車(chē)室內(nèi)氣流組織的沿線變化數(shù)值分析[J]. 重慶大學(xué)學(xué)報(bào)，2014，37(9)：59－67. Wang Ye, Xu Ran. Numerical study of airflow distribution in air-conditioned train in Qinghai-Tibet Railway[J]. Journal of Chongqing University, 2014, 37(9): 59－67. (in Chinese with English abstract)

[32] Xu Feng, Patterson J C, Lei Chengwang. An experimental study of the unsteady thermal flow around a thin fin on a sidewall of a differentially heated cavity[J]. International Journal of Heat and Fluid Flow, 2008, 29(4): 1139－1153.

Numerical analysis of natural convection cooling effect in closed cavity of electronic equipment

Wang Ye1,2, Zhao Xingjie1, Lin Huxiang1, Song Rongfei1, Guan Guoxiang1

(1.,,, 730070,; 2.,,,, 730070,)

Natural convection in closed cavity has been received considerable attention due to its widely applications in industry, such as in solar energy collectors design, cooling of electronic instruments, energy saving of building and nuclear reactor design, the study of natural convection mechanism in enclosed cavity is of great practical significance for improving the comprehensive performance of these systems. With the rapid development of science and technology, more and more heat sources and components are located in the same enclosed cavity, which will lead to large number of heat generation. How to arrange the heat source elements is very important to the cooling and reliability of the system by the natural convection heat transfer. Lots of studies have shown that the fins attached to the vertical walls inside the cavity can enhance the heat transfer rate in the cavity. In order to study the effect of the material and space position of heat radiation in enclosed space on the heat transfer characteristics of natural convection, the RNG-model of FLUENT14.5 was adopted to analyze the temperature field, flow field, local Nusselt numbers and the average Nusselt numbers of the vertical wall of the air-filled cavity (left side was hot wall, right side was cold wall and the top side was the heat source) with aspect ratio of 1. The results showed that the heat transfer enhancement of the cavity was strongly dependent on the position and material of the fins attached to the vertical walls. The influence of single fin and double fin on the horizontal velocity near the top of the cavity was weak, and the disturbance to the horizontal velocity field near the bottom of the cavity was intense, and the influence trend was very close. This phenomenon was not consistent with that of without heat source. The effect of fins on the horizontal velocity in the core region of the cavity was quite different: the disturbance of single fin on the lower position of the core region of the cavity was more obvious, while the disturbance of double fins on the upper region was more intense. Different number of fins lead to different distribution characteristics of temperature field in the enclosed cavity by changing the flow field structure. When the adiabatic fin was located at 1/6 height above the bottom of the hot wall, the average heat transfer capacity of the hot wall was the strongest. When the heat conducting fin was located at 1/3 height above the bottom of the hot wall, the average heat transfer capacity of the hot wall was the strongest. The combination of 1/3 height above the bottom of the hot wall and 2/3 height above the bottom of the cold wall made the average heat transfer capacity of the hot wall the strongest. The averagenumber of thermal conductive fins arranged at 1/3 height above the bottom of the hot wall was maximum, which was 34.93 and 9.67% higher than that of adiabatic fin arranged at the same location. The averagenumber of thermal conductive fins arranged at 1/3 height above the bottom of the hot wall and 2/3 above the bottom of the cold wall was maximum, which was 14.3% and higher than that of thermal conductive fins arranged only at 1/3 above the bottom of the hot wall. The study has a certain theoretical significance for improving the natural convection cooling effect of electronic components in industrial and agricultural engineering and optimizing the spatial location of heat dissipation components.

numerical simulation; temperature; heat convection; enclosed cavity; fin

2018-10-18

2018-11-28

國(guó)家自然科學(xué)基金（51476073，51266004）；甘肅省建設(shè)科技攻關(guān)項(xiàng)目（KJ2016-2）

王燁，教授，博士，博士生導(dǎo)師，主要從事自然對(duì)流及暖通空調(diào)研究。Email：wangye@mail.lzjtu.cn

10.11975/j.issn.1002-6819.2019.06.026

TK124

A

1002-6819(2019)-06-0214-08

王 燁，趙興杰，藺虎相，宋榮飛，管?chē)?guó)祥. 電子設(shè)備封閉腔內(nèi)自然對(duì)流冷卻效果數(shù)值分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2019，35(6)：214－221. doi：10.11975/j.issn.1002-6819.2019.06.026 http://www.tcsae.org

Wang Ye, Zhao Xingjie, Lin Huxiang, Song Rongfei, Guan Guoxiang. Numerical analysis of natural convection cooling effect in closed cavity of electronic equipment[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(6): 214－221. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2019.06.026 http://www.tcsae.org