張軍偉,馬兵善,王 剛

(蘭州理工大學土木工程學院,甘肅 蘭州 730050)

多孔介質(zhì)內(nèi)對流傳熱現(xiàn)象廣泛存在于核反應(yīng)堆堆芯元件的設(shè)計與安全運行、換熱器換熱過程、電子設(shè)備冷卻過程及物料干燥等許多工程實際領(lǐng)域中,對其進行研究具有十分重要的意義。因此,近幾十年以來,許多學者已經(jīng)從不同角度對多孔介質(zhì)內(nèi)的對流傳熱問題開展了廣泛深入的研究,取得了大量的研究成果,Nield等[1]對此進行了全面、詳細的回顧。

在電子設(shè)備冷卻換熱的研究中,經(jīng)常會碰到開口的多孔介質(zhì)腔體內(nèi)混合對流換熱問題,腔體內(nèi)按一定方式布置的針翅肋片或堆積的小球體可看作多孔介質(zhì)固體骨架。Bhuiyan等[2]研究了多孔介質(zhì)方形通風腔體內(nèi)混合對流的熱特性,發(fā)現(xiàn)腔體高寬比和流體入口寬度對流體流動和傳熱特性有顯著影響。Mohammed等[3]對開口多孔介質(zhì)腔體內(nèi)混合對流進行了數(shù)值模擬,結(jié)果表明腔體內(nèi)多孔固體骨架對平均傳熱速率的影響主要取決于理查德森數(shù)Ri,即當Ri數(shù)逐漸減小到使混合對流幾乎為強制對流時,平均努塞爾數(shù)增大。Ataei-Dadavi等[4]在開口腔體內(nèi)填充直徑較大的低導(dǎo)熱系數(shù)固體球作為多孔介質(zhì)固體骨架,對其內(nèi)的混合對流換熱進行了實驗研究,實驗結(jié)果表明,根據(jù)Ri數(shù)的大小,腔體內(nèi)存在3種不同的流動及傳熱區(qū)域。Gibanov等[5]研究了部分填充有多孔介質(zhì)且底部壁面具有局部熱源的開口方腔內(nèi)的混合對流,發(fā)現(xiàn)隨Ri數(shù)的增大,混合對流換熱得以增強,而且達西數(shù)Da從10-7增加到10-3,平均努塞爾數(shù)逐漸減小。孫夢楠等[6]對內(nèi)置有熱源的開口方腔內(nèi)的混合對流進行了數(shù)值模擬,計算結(jié)果表明,隨著內(nèi)熱源圓角半徑的增大,平均傳熱速率增大,并且流體進、出口分別位于腔體左側(cè)壁面中部和右側(cè)壁面頂部時,取得最大值。王婷婷等[7]采用格子Blotzmann方法對內(nèi)置有恒溫發(fā)熱圓柱體的開口多孔方腔內(nèi)的混合對流進行了數(shù)值模擬,發(fā)現(xiàn)Ri數(shù)和Da數(shù)對混合對流產(chǎn)生重要的影響。

考慮到開口多孔介質(zhì)方腔中混合對流換熱特性對電子設(shè)備高效冷卻換熱具有重要的意義,因此,對下壁面為高溫壁面的開口多孔介質(zhì)方腔內(nèi)混合對流進行數(shù)值研究,分析不同的Ri數(shù)、Da數(shù)、多孔介質(zhì)孔隙率ε和Ra數(shù)對混合對流換熱特性的影響。

1 物理模型和控制方程

1.1 物理模型

多孔介質(zhì)方腔的物理模型與坐標系統(tǒng)如圖1所示,二維方腔中填充均質(zhì)、各向同性的飽和多孔介質(zhì),腔體的邊長為L,流體進出口長度均為D(D=0.1L),流體入口速度為U0,溫度為Tc,腔體左、右側(cè)壁面及上壁面均為絕熱壁面,下壁面是溫度為Th的恒溫壁面。

圖1 物理模型與坐標系統(tǒng)Fig.1 Physical model and coordinate system

1.2 控制方程

數(shù)值計算中,假設(shè)流體為不可壓縮牛頓流體,多孔介質(zhì)固體骨架與流體之間滿足局部熱平衡,引入Boussinesq假設(shè),滿足無滑移邊界條件,忽略粘性耗散。對于所研究的二維多孔介質(zhì)方腔中的混合對流問題,其無量綱控制方程如下:

(1)

(2)

(3)

(4)

方程(1)～(4)中涉及的無量綱參數(shù)如下:

(5)

所研究問題的無量綱邊界條件如下:

流體入口截面:U=1,V=0,θ=0；

流體出口截面:?U/?X=0,V=0, ?θ/?X=0；

上部絕熱壁面:U=V=0, ?θ/?Y=0；

豎直絕熱壁面:U=V=0, ?θ/?X=0；

下部高溫壁面:U=V=0,θ=1。

(6)

腔體熱壁面的平均努塞爾數(shù)Nuav的計算式為

(7)

2 數(shù)值計算方法

采用SIMPLE算法對離散后的無量綱控制方程進行計算求解,控制方程對流項采用QUICK格式進行離散,使計算結(jié)果具有較高的精度[9]。為了得到更準確的結(jié)果,計算之前對文獻[10]中所研究的多孔介質(zhì)腔體內(nèi)自然對流問題進行了驗證,平均努塞爾數(shù)Nuav的計算結(jié)果如表1所列。從表1可以看出,平均努塞爾數(shù)Nuav計算值與文獻[10]中的相應(yīng)值吻合較好,保證了計算程序的可靠性。同時,采用60×60、70×70、80×80、90×90、100×100 5套非均分網(wǎng)格進行了網(wǎng)格獨立性驗證,結(jié)果表明,當網(wǎng)格數(shù)為90×90時,可得到與網(wǎng)格無關(guān)的解,故數(shù)值計算時均采用90×90的非均分網(wǎng)格。

表1 平均努塞爾數(shù)Nuav計算值與文獻值的比較Table 1 Comparison of average Nusselt number with that of the previous results in reference

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

對有空氣流入和流出的開口多孔介質(zhì)方腔內(nèi)的混合對流換熱進行了數(shù)值計算,分析研究了達西數(shù)Da、理查德森數(shù)Ri、瑞利數(shù)Ra和多孔介質(zhì)孔隙率ε對腔體內(nèi)混合對流換熱特性的影響。計算中取Pr=0.703、Ri=0.01～100、Da=10-5～10-2、Ra=103～106、ε=0.5～0.9。

3.1 Ri數(shù)對混合對流換熱的影響

Da=10-3,Ra=104,ε=0.6時,不同Ri數(shù)時多孔介質(zhì)方腔中流線和等溫線的分布如圖2所示。從圖2中可以看出,當Ri=0.01時,在流體入口的上部出現(xiàn)了漩渦,隨著Ri數(shù)的不斷增大,漩渦會逐漸減小,在Ri=0.1時漩渦消失;在Ri≥1時,隨著Ri數(shù)的增大,流線的分布幾乎不再受其影響。從等溫線分布圖中可看出,當Ri≤0.1時,等溫線集中在腔體的右下部,且靠近入口處分布較密集,而靠近出口處的右上部分布較為稀疏;當Ri≥1時,隨著Ri數(shù)的不斷增大,等溫線分布逐漸向左側(cè)移動,最終主要集中分布于流體入口處。

圖2 不同Ri數(shù)下流線和等溫線分布 (Da=10-3,Ra=104,ε=0.6)Fig.2 Distributions of streamlines and isotherms at different Ri number (Da=10-3,Ra=104,ε=0.6)

Ra=103,ε=0.5及不同Da數(shù)時,Nu數(shù)隨Ri數(shù)變化的曲線如圖3所示。從圖3可以看出,當Ra數(shù)、ε和Da數(shù)相同時,Nu數(shù)隨Ri數(shù)的增大而減小,Ri=100相對于Ri=0.01時Nu數(shù)減小了將近60%;當Ri≤1時,Ri數(shù)對Nu數(shù)的影響較為明顯,這是由于強制對流占主導(dǎo)地位,自然對流的影響幾乎可忽略；Ri<1時,對于相同的Ri數(shù),Nu數(shù)隨Da數(shù)的增大而增大,但當Da數(shù)為10-4和10-5時,其對Nu數(shù)幾乎無影響；在Ri=0.01時,Da=10-2的Nu數(shù)比Da=10-5時增加了18.7%;在Ri=0.1時,Da=10-2的Nu數(shù)比Da=10-5時增加了8.7%;當Ri≥1時,Da數(shù)對Nu數(shù)幾乎沒有影響。

3.2 Da數(shù)對混合對流換熱的影響

Ri=0.1,ε=0.5及不同Ra數(shù)時,Nu數(shù)隨Da數(shù)的變化曲線如圖4所示。從圖4中可以看出,隨著Da數(shù)的增大,Nu數(shù)有增大趨勢但不明顯。當Da>10-4時,隨Da數(shù)的增大,Nu數(shù)逐漸增大,Da=10-2相較Da=10-3,Nu數(shù)最大增大了12.4%,這是由于隨著Da數(shù)的增大,多孔介質(zhì)阻力對混合對流的影響逐漸減小,熱壁面的Nu數(shù)受Da數(shù)的影響逐漸增大。當Ra數(shù)較小時,Da數(shù)的變化對Nu數(shù)幾乎沒影響,這是由于Ra較小時,腔體內(nèi)流體的流動強度相對較弱,同時多孔介質(zhì)阻力對腔體內(nèi)混合對流有一定的影響。當Da數(shù)相同時,Nu數(shù)隨Ra數(shù)的增大而增大,這是因為隨著Ra數(shù)的不斷增大,腔體內(nèi)流體的流動強度在逐漸加強。Da=10-2時,Ra=105的Nu數(shù)較Ra=104增大了59.5%,Ra=106的Nu數(shù)相較Ra=105增大了63.6%。

圖4 不同Ra數(shù)時Nu數(shù)隨Da數(shù)的變化 (Ri=0.1,ε=0.5)Fig.4 Variation of Nu number with Da number at different Ra number (Ri=0.1,ε=0.5)

3.3 孔隙率對混合對流換熱的影響

Ra=106,Ri=100及不同Da數(shù)時,Nu數(shù)隨孔隙率ε的變化如圖5所示。從圖5中可以看出,當Da≤10-4時,ε的變化對Nu數(shù)幾乎沒有影響;當Da>10-4時,隨ε的增大,Nu數(shù)逐漸增大,ε=0.9相較ε=0.5最大增加了5%,這是由于孔隙率增大,流體更易流動。對于相同ε而言,Nu數(shù)隨Da數(shù)的增加顯著增大,當ε=0.8時,Da=10-2相較于Da=10-5的Nu數(shù)增加了39.2%。

圖5 不同Da數(shù)時Nu數(shù)隨ε的變化 (Ra=106,Ri=100)Fig.5 Variation of Nu number with ε at differentDa number (Ra=106,Ri=100)

3.4 Ra數(shù)對混合對流換熱的影響

Da=10-3,ε=0.8及不同Ri數(shù)時,Nu數(shù)隨Ra數(shù)的變化曲線如圖6所示。從圖6中可以看出,Nu數(shù)隨Ra數(shù)的增大而增大,這是因為在Ra數(shù)較小時,腔體內(nèi)流體流動強度較弱,換熱以導(dǎo)熱為主,隨著Ra數(shù)的增大,腔體內(nèi)流體流動強度逐漸增強,對流換熱強度逐漸增強。對于相同的Ra數(shù),Nu數(shù)隨Ri數(shù)的減小而增大,Ri=0.01的Nu數(shù)幾乎是Ri=100的6倍,這是由于隨著Ri數(shù)的減小,強制對流在混合對流中逐漸占據(jù)主導(dǎo)地位,而自然對流逐漸減弱。

圖6 不同Ri數(shù)時Nu數(shù)隨Ra數(shù)的變化 (Da=10-3,ε=0.8)Fig.6 Variation of Nu number with Ra number at different Ri number (Da=10-3,ε=0.8)

4 結(jié)論

對開口多孔介質(zhì)方腔內(nèi)的混合對流換熱進行了數(shù)值模擬,分析了不同Ri數(shù)、Da數(shù)、ε和Ra數(shù)對腔體內(nèi)混合對流換熱的影響,得出如下主要結(jié)論:

(1) 隨著Ri數(shù)增大,等溫線從主要分布于方腔右下部向左側(cè)移動,最終主要集中分布于流體入口處;隨著Ri數(shù)增大,入口處上部的漩渦逐漸減小,在Ri=0.1時消失,之后流線分布幾乎不受Ri數(shù)的影響。

(2) 腔體熱壁面Nu數(shù)隨Ri數(shù)的增大而減小,隨Ra數(shù)的增大而增大,當Da<10-4時,Nu數(shù)變化不再受Ri數(shù)影響;ε一定時,Nu數(shù)隨Da數(shù)的增大而增大,而ε變化對Nu數(shù)影響不明顯。