第五嘉瑋,王 燁,2,胡潤鑫
(1.蘭州交通大學(xué) 環(huán)境與市政工程學(xué)院, 蘭州 730070; 2.蘭州交通大學(xué) 鐵道車輛熱工教育部重點(diǎn)實(shí)驗(yàn)室, 蘭州 730070)
自然對流因其自發(fā)性、無額外能耗[1]等優(yōu)點(diǎn)在建筑節(jié)能[2]、變壓器散熱[3]、電子設(shè)備冷卻[4]等諸多工業(yè)生產(chǎn)系統(tǒng)中被廣泛應(yīng)用。陳岳[5]在微尺度下內(nèi)置發(fā)熱電子芯片方腔對流換熱的研究工作中,對電子芯片尺寸、多芯片布局對芯片表面對流冷卻性能的影響進(jìn)行了數(shù)值分析,得到瑞利數(shù)Ra=104、105時,腔內(nèi)對流散熱效率最大的臨界芯片尺寸分別為0.25H、0.39H。由此可見,內(nèi)置熱源封閉腔內(nèi)熱邊界層附近的自然對流換熱是工程應(yīng)用和理論研究的重要問題。自Batchelor[6]最早提出側(cè)加熱封閉腔自然對流模型以來,人們先從腔體高寬比[7]、瑞利數(shù)Ra[8]等方面進(jìn)行相關(guān)探索。后來發(fā)現(xiàn),通過側(cè)壁周期性邊界條件變化[9]、填充納米流體[10]、水平薄翅片側(cè)置[11]等方式可以強(qiáng)化腔內(nèi)對流傳熱。其中,布置翅片操作簡單,節(jié)約能源,更具適用性。有關(guān)翅片材質(zhì)[12]、翅片幾何參數(shù)及定位尺寸[13]、翅片形狀[14]等方面的研究已經(jīng)實(shí)現(xiàn)了技術(shù)成果的多方面突破與應(yīng)用。Patnode[15]數(shù)值分析了翅片對側(cè)加熱方腔內(nèi)低壓氣流的擾動作用,發(fā)現(xiàn)在相同的克努森數(shù)Kn下,熱壁面布置雙翅片相比無翅片空腔,熱壁面平均Nu數(shù)會大幅提升。Bilgen[16]在溫度隨時間周期性變化的封閉腔熱壁面上布置薄固體翅片,數(shù)值分析得出最佳的翅片安裝高度位于熱壁面中心或靠近中心位置。Elatar[17]在封閉腔熱壁面1/2高度處布置單個不同長度的水平導(dǎo)熱翅片,分析了翅片與流體工質(zhì)的導(dǎo)熱系數(shù)比對翅片效率的影響,并得到使腔內(nèi)換熱效果最佳的翅片長度。何騰[18]找到了強(qiáng)化或抑制內(nèi)置熱源封閉腔內(nèi)對流換熱的翅片臨界長度。Dou等[13]數(shù)值分析了不同瑞利數(shù)Ra下,翅片位置、翅片長度及翅片數(shù)量對腔內(nèi)傳熱速率的影響。研究發(fā)現(xiàn),將翅片固定在熱壁中間高度,對腔內(nèi)傳熱效率提升效果最優(yōu)。當(dāng)Ra=3.38×106,翅片數(shù)量大于1,翅片數(shù)量對腔內(nèi)傳熱的影響可忽略不計。當(dāng)Ra=3.38×109,雙翅片工況傳熱速率較單翅片有所降低。為改善工業(yè)裝置封閉空間發(fā)熱電子元件對流冷卻效果,王燁[19]探究了壁面發(fā)射率對內(nèi)置單、雙翅片的封閉腔內(nèi)湍流自然對流的影響,發(fā)現(xiàn)單翅片結(jié)構(gòu)封閉腔的熱壁面局部Nu數(shù)與壁面發(fā)射率呈正相關(guān)變化關(guān)系,壁面發(fā)射率過大不利于雙翅片結(jié)構(gòu)封閉腔內(nèi)的對流散熱。
綜上,內(nèi)置翅片封閉腔相關(guān)研究主要針對翅片單參數(shù)變化對無內(nèi)熱源腔體內(nèi)自然對流換熱產(chǎn)生的影響展開討論,尋找使得內(nèi)置熱源封閉腔內(nèi)對流-輻射耦合傳熱效果最佳的翅片參數(shù)組合的研究,還未見報道。本文中研究了翅片單參數(shù)變化下壁面輻射效應(yīng)對腔內(nèi)對流傳熱特性的影響,并以熱源表面平均Nu數(shù)為評價指標(biāo),引入響應(yīng)面法建立回歸方程,獲得翅片多參數(shù)交互作用下,腔內(nèi)對流-輻射耦合傳熱效率的變化規(guī)律,所得腔內(nèi)對流-輻射耦合傳熱的最優(yōu)翅片參數(shù)組合,為優(yōu)化工業(yè)系統(tǒng)封閉空間內(nèi)散熱元器件的空間布局及對流冷卻效果提供理論參考。
封閉腔物理模型如圖1所示。腔體寬高均為H,左右等溫冷壁面溫度為Tc,底面中心放置邊長0.2H、表面溫度Th的方形熱源。腔內(nèi)初始溫度T0=(Th+Tc)/2,頂面與底面絕熱。對流-輻射耦合傳熱下的腔體壁面及熱源表面發(fā)射率分別取ε1=0.4、ε2=0.8,采用能同時考慮散射氣體與顆粒間輻射換熱,還能考慮非灰體輻射和局部熱源影響的DO輻射模型[20]。腔內(nèi)流體介質(zhì)為空氣,黏性流體物性為常數(shù)且不可壓縮,流動為恒定層流,忽略黏性耗散和體積力[21],瑞利數(shù)Ra=105、普朗特數(shù)Pr=0.71。通過在腔內(nèi)右側(cè)冷壁面布置厚度δ=0.01H的薄翅片來改善腔內(nèi)自然對流傳熱效果。不同翅片參數(shù)變化如表1所示,研究單參數(shù)影響規(guī)律時,其余翅片參數(shù)分別固定為a=0.5H、l=0.16H、θ=90°

圖1 封閉腔物理模型

表1 翅片參數(shù)
為進(jìn)一步拓展研究結(jié)論的普適性,需要對腔內(nèi)對流傳熱控制方程進(jìn)行無量綱化處理。引入以下無量綱參數(shù)[22]:

(1)
式(1)中:x、y分別為水平和豎直方向的坐標(biāo)分量,m;u、v分別為x、y方向的速度分量,m/s;p為空氣壓力,Pa;ρf為空氣密度, kg/m3;Tc為左、右側(cè)冷壁面溫度,Tc=297.15 K;Th為熱源表面溫度,Th=323.15 K;為空氣的運(yùn)動黏度系數(shù),m2/s;β為流體的體積膨脹系數(shù),1/K;g為重力加速度,m/s2;Pr為空氣擴(kuò)散系數(shù),W/(mK);Pr為普朗特數(shù);Ra為瑞利數(shù)。
求解封閉腔內(nèi)空氣流動與傳熱過程的無量綱化控制方程如下:
連續(xù)性方程:

(2)
動量方程:
X方向:

(3)
Y方向:

(4)
能量方程:

(5)
式(2)—(5)中:X、Y分別為水平和豎直方向的無量綱坐標(biāo)分量;H為腔體無量綱寬(高)度;U、V分別為X、Y方向無量綱速度分量;P為無量綱壓力;Θ為無量綱溫度。
無量綱邊界條件如下:
左、右側(cè)低溫壁面:U=V=0,Θ=0
熱源表面:U=V=0,Θ=1
冷壁面和熱源表面局部Nu數(shù)、熱源表面平均Nu數(shù)的表達(dá)式分別如下:
冷壁面和熱源表面局部Nu數(shù)為:
熱源表面平均Nu數(shù)為:
利用ANSYS-FLUENT 2020R2軟件求解腔內(nèi)對流換熱過程,采用有限容積法離散控制方程。采用比UPWIND格式精度更高,擁有三階截斷誤差且更具對流穩(wěn)定性的QUICK格式離散對流項(xiàng),用中心差分格式離散擴(kuò)散項(xiàng),用PRESTO算法進(jìn)行壓力方程修正,并使用最具收斂性的SIMPLEC算法處理壓力-速度耦合,求解方程亞松弛因子設(shè)置:壓力為0.3,動量為0.7,能量方程、體積力和密度均為1[23]。通過迭代離散方程得到每個變量,滿足如下準(zhǔn)則,認(rèn)為求解過程達(dá)到收斂。
為確保數(shù)值計算方法的可靠性,與文獻(xiàn)[24]的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對比,得到腔內(nèi)Y=0.5水平線上的豎向無量綱速度分布,如圖2所示,最大相對誤差為5.52%,數(shù)值及趨勢均吻合良好,故認(rèn)為本文中數(shù)學(xué)模型可用于內(nèi)置翅片參數(shù)對含內(nèi)熱源封閉腔內(nèi)純自然對流傳熱及對流-輻射耦合傳熱特性影響的對比研究。

圖2 數(shù)學(xué)模型驗(yàn)證
圖3為3組近壁面局部加密的結(jié)構(gòu)化網(wǎng)格計算所得的腔內(nèi)Y=0.5水平線上的無量綱溫度,吻合良好,后續(xù)計算選取130×130的網(wǎng)格數(shù)。

圖3 網(wǎng)格獨(dú)立性檢驗(yàn)
研究發(fā)現(xiàn),翅片單參數(shù)變化條件下,翅片長度對腔內(nèi)物理場邊界層發(fā)展?fàn)顩r影響最大,故圖4與圖5分別展示了不考慮壁面輻射(左)和考慮壁面輻射(右)時,不同翅片長度對腔內(nèi)溫度場和流場結(jié)構(gòu)的影響對比。

圖4 翅片長度對腔內(nèi)溫度場的影響對比(a=0.5H,θ=90°)

圖5 翅片長度對腔內(nèi)流場的影響對比(a=0.5H,θ=90°)
由圖4、圖5可知,無論是否考慮輻射,當(dāng)翅片長度l≤0.02H,溫度場與流場結(jié)構(gòu)均呈近似對稱分布,影響不顯著。但隨著翅片長度增大,腔體右側(cè)自頂部沿冷壁面向下不斷形成的溫度邊界層與速度邊界層逐漸變厚且遭受到擠壓破壞,存在1個翅片臨界長度(0.16H≤l≤023H),使得腔體右半?yún)^(qū)的等溫線與順時針渦流由被擠壓分別到發(fā)生“錯位斷裂”與“分割”。隨著翅片長度增加,翅片上表面與腔體右壁面交角處會出現(xiàn)冷滯留現(xiàn)象,翅片周圍低溫區(qū)域持續(xù)擴(kuò)大。當(dāng)翅片長度l>0.16H,由于腔體右壁面附近向下流動的空氣受翅片阻礙作用加大,導(dǎo)致翅片上下表面附近空氣流動方向不同,翅片上下區(qū)域溫度梯度變大,腔體底部區(qū)域等溫線向翅片附近發(fā)生顯著偏移。值得注意的是,當(dāng)翅片長度l≥0.23H時,腔內(nèi)流函數(shù)值整體會減小。
與不考慮壁面輻射相比,相同翅片長度下,壁面輻射使得腔內(nèi)整體區(qū)域的溫度值及流函數(shù)值均變化不大,但頂部區(qū)域溫度梯度明顯減小,腔內(nèi)靠近腔體頂部區(qū)域的同水平等溫線及溫度區(qū)域分布結(jié)構(gòu)發(fā)生了一定的變化,且前述冷滯留現(xiàn)象也有所加強(qiáng),翅片上表面與腔體右壁面交角處的低溫區(qū)域較不考慮輻射工況有所擴(kuò)大。由此可見,輻射效應(yīng)會通過改變腔內(nèi)空氣的整體溫度分布結(jié)構(gòu)進(jìn)而影響腔內(nèi)的對流傳熱特性。相比不考慮輻射,同一翅片長度下,腔內(nèi)各區(qū)域的流函數(shù)值及左、右兩側(cè)逆時針與順時針對稱分布的流場結(jié)構(gòu)變化不大。
圖6(a)—圖6(c)左側(cè)為不同翅片單參數(shù)變化對腔內(nèi)熱源表面局部Nu數(shù)的影響對比。無論是否考慮輻射,不同翅片單參數(shù)工況下熱源表面局部Nu數(shù)關(guān)于腔體半寬處(X=0.5)呈對稱分布。腔體底部的低溫空氣經(jīng)熱源AD和BC面加熱后受浮升力作用向腔體頂部遷移,經(jīng)過熱源頂點(diǎn)A和B時,空氣在慣性與熱浮升力作用下運(yùn)動方向依然保持向上趨勢,經(jīng)過A、B點(diǎn)后與熱源頂面AB分離,導(dǎo)致熱源AB面局部Nu數(shù)自A和B向中點(diǎn)迅速減少,呈“U”型分布。而AD和BC面的局部Nu數(shù)則分別自D向A和自C向B逐漸增大,這是由于熱源頂點(diǎn)A與B附近的空氣受熱程度較靠近腔體底部絕熱壁面的D與C點(diǎn)更高,導(dǎo)致受熱浮升力沿D到A與沿C到B向上運(yùn)動的空氣對流換熱強(qiáng)度逐漸增強(qiáng),而熱源所有端點(diǎn)的局部Nu數(shù)均會隨流動邊界層的發(fā)展發(fā)生急劇變化。

圖6 不同翅片單參數(shù)對冷壁面和熱源表面局部Nu數(shù)的影響
不同翅片單參數(shù)改變,均會使得熱源BC面附近的氣流渦旋尺度及被擠壓程度發(fā)生變化,對熱源附近溫度場及流場結(jié)構(gòu)產(chǎn)生不同程度的影響,從而導(dǎo)致不同翅片單參數(shù)工況對應(yīng)的熱源局部Nu數(shù)數(shù)值分布及最大值獲取位置存在一定差異。同時發(fā)現(xiàn),考慮壁面輻射時,翅片單參數(shù)變化對于不同熱源表面的局部Nu數(shù)的提升程度不盡相同。其中,翅片安裝角度對熱源表面局部Nu數(shù)影響程度為:AD面>BC面>AB面,翅片長度與翅片安裝高度則均為:AD面>AB面>BC面。由數(shù)值試驗(yàn)可知,與無翅片工況相比,熱源AD、AB、BC面平均Nu數(shù)分別最高提升了15.14%、12.86%、13.51%,單參數(shù)改變對熱源AD面對流傳熱效果的提升最為顯著,任一翅片單參數(shù)變化相比無翅片工況均提升了熱源AD、AB和BC面的散熱能力。而不考慮壁面輻射,發(fā)現(xiàn)只有當(dāng)翅片安裝角度θ≥90°,翅片長度l≤0.23H,翅片安裝高度a≥0.5H時,3個熱源表面的散熱能力才都會得到提升。表明一定強(qiáng)度的壁面輻射會有利于熱源表面的對流散熱。
由前述腔體物理場分析可得,翅片附近空氣對流擾動狀況的變化最為顯著,分析放置翅片一側(cè)的冷壁面局部Nu數(shù)變化有助于翅片附近位置壁面局部邊界層發(fā)展規(guī)律的得出。圖6(a)—圖6(c)右側(cè)為不同翅片單參數(shù)變化下,冷壁面(X=1.0)局部Nu數(shù)的變化曲線。無論是否考慮輻射,翅片附近的冷壁面局部邊界層發(fā)展?fàn)顩r較其余區(qū)域受影響更大,冷壁面局部Nu數(shù)會由于翅片對附近空氣產(chǎn)生不同程度的擾動而發(fā)生有差異的波動。壁面局部邊界層發(fā)展?fàn)顩r受翅片擾動,冷壁面局部Nu數(shù)會在不同單參數(shù)變化下的翅片附近產(chǎn)生急劇變化,且考慮輻射之后波動更為劇烈。由于腔內(nèi)空氣經(jīng)熱源加熱向腔體上部流動,冷壁面附近空氣沿Y軸正方向的對流換熱強(qiáng)度不斷加強(qiáng),因此,翅片附近的冷壁面局部Nu數(shù),雖然存在不同程度波動,但總體呈增大趨勢。不同單參數(shù)改變,均會使得來流空氣通過翅片附近時的冷滯留區(qū)與渦流結(jié)構(gòu)發(fā)生變化,從而對腔內(nèi)溫度場及流場結(jié)構(gòu)產(chǎn)生不同程度的影響,導(dǎo)致不同翅片單參數(shù)工況對應(yīng)的冷壁面局部Nu數(shù)數(shù)值分布及最大值獲取位置存在一定差異。
輻射效應(yīng)對于內(nèi)置相同單參數(shù)翅片的冷壁面局部Nu數(shù)的總體變化趨勢影響不大,但翅片附近的冷壁面局部Nu數(shù)波動程度較不考慮輻射顯著增強(qiáng),壁面局部邊界層發(fā)展?fàn)顩r更為劇烈,不同翅片單參數(shù)變化對冷壁面局部區(qū)域的對流換熱強(qiáng)化效果較不考慮輻射更為顯著。數(shù)值計算得知,考慮壁面輻射,當(dāng)θ=120°(a=0.375H,l=0.16H),冷壁面平均Nu數(shù)提升最顯著,相比無翅片工況提升了16.74%。
圖7為僅考慮翅片單參數(shù)變化腔內(nèi)熱源表面平均Nu數(shù)及其提升率η的變化曲線。由圖7可得,無論是否考慮輻射,隨任一翅片單參數(shù)變化,熱源表面平均Nu數(shù)與無翅片工況相比均有不同程度的提升。熱源表面平均Nu數(shù)及其提升率η,隨翅片角度的增加呈先減小后增大的趨勢,當(dāng)60°<θ<90°,變化幅度最大。考慮輻射與不考慮輻射的情況下,熱源表面平均Nu數(shù)最大值分別為25.80、11.14,相比無翅片工況分別提升了10.15%、7.00%。

圖7 不同翅片單參數(shù)對熱源表面平均Nu數(shù)及其提升率的影響