局部強內(nèi)熱源移動球床溫度分布特性

王永偉，淮秀蘭，李勛鋒，蔡軍，席文宣

（中國科學(xué)院工程熱物理研究所，北京 100190）

引言

移動球床裝置在核電站、化工、生物等很多工業(yè)生產(chǎn)過程中具有廣泛應(yīng)用，如用于反應(yīng)器、干燥器、過濾器和換熱器等。在這些應(yīng)用中通常伴隨著如核裂變和嬗變、電磁輻射、化學(xué)或生物反應(yīng)等物理和化學(xué)變化過程。這些變化過程往往會放出熱量，在移動球床裝置中形成局部內(nèi)熱源，影響設(shè)備的運行性能。因此，含有局部內(nèi)熱源移動球床的溫度分布特性研究對提高設(shè)備性能具有重要的意義。

許多學(xué)者對含有內(nèi)熱源的填充床傳熱性能進行了實驗和數(shù)值研究[1-4]。Meng等[5]對含有內(nèi)熱源的球床通道中單相流體對流換熱性能進行了實驗研究，獲得了與實驗數(shù)據(jù)吻合良好的量綱1平均傳熱系數(shù)關(guān)系式。Zeisberger等[6]對含有內(nèi)熱源的多孔介質(zhì)填充床的兩相沸騰傳熱狀態(tài)進行了研究。結(jié)果表明，在高于最小加熱速度時增加冷卻液流量會減小蒸氣組分，而系統(tǒng)壓力對沸騰沒有影響。楊劍等[7]采用Forchheimer-Brinkman拓展Darcy模型和局部非熱平衡模型，在層流范圍內(nèi)對骨架發(fā)熱多孔介質(zhì)豎直通道內(nèi)的非達西強制對流換熱進行了數(shù)值研究。研究表明Reynolds數(shù)、有效熱導(dǎo)率比和Darcy數(shù)對流道內(nèi)流動換熱影響顯著。此外，近年來也有很多關(guān)于移動床傳熱性能的研究[8]。Collier等[9]實驗研究了高溫銅球與填充床或流化床之間的傳熱系數(shù)。基于銅球和氣體之間的導(dǎo)熱率獲得了 Nussult數(shù)實驗關(guān)聯(lián)式。Marcelo等[10]構(gòu)建了移動床湍流傳熱方程，研究表明，增加固體顆粒運動速度可以減小界面阻力，當(dāng)流固相速度相同時，湍流傳熱被抑制，而流固相間能量傳遞主要依靠熱傳導(dǎo)作用。Ana等[11]運用雙能量方程對順流移動多孔填充床反應(yīng)器層流流動時固液相間傳熱的數(shù)值研究表明，提高固相和液相之間的相對速度可以減小相間的能量傳遞。總之，以往的研究或者關(guān)注于含有內(nèi)熱源的填充床而沒有考慮固體顆粒的運動問題，或者關(guān)注于移動床而沒有考慮內(nèi)熱源的問題，關(guān)于含有局部內(nèi)熱源移動球床通道內(nèi)的固體顆粒運動和溫度分布等換熱特性研究存在不足。

對于含有局部內(nèi)熱源移動球床來說，顆粒的運動速度是影響溫度分布和傳熱特性的關(guān)鍵因素。顆粒的運動速度與通道內(nèi)溫度分布有著必然聯(lián)系，研究二者的關(guān)系對分析局部強內(nèi)熱源移動球床傳熱特性有著重要的意義。本文運用理論分析和數(shù)值計算相結(jié)合的方法研究含有強內(nèi)熱源移動球床中顆粒的運動速度與出口溫度、中心溫度，以及整體溫度分布特點之間的關(guān)系。

1 物理模型

含有局部強內(nèi)熱源移動球床模型如圖1所示。圓筒形移動球床通道長為L，半徑為r2。移動球床通道的中上部有一恒定功率的內(nèi)熱源，加熱區(qū)域也為圓筒形，長為L1，半徑為r1，內(nèi)熱源區(qū)加熱功率在兆瓦級別屬于強熱源。移動球床中隨機填滿了等徑球形金屬顆粒，球徑遠遠小于通道尺寸，顆粒間的空隙中充滿了氣體介質(zhì)，并隨顆粒一起運動。初始溫度為T0的球形顆粒以恒定速度v自通道上部入口進入，通過內(nèi)熱源區(qū)加熱后溫度升高至T1，然后自通道下部出口離開移動球床，出口溫度為T2。

圖1 含有局部強內(nèi)熱源移動球床Fig.1 Schematic diagram of moving pebble bed with local internal heat generation

本文的研究基于以下條件展開：① 內(nèi)熱源加熱區(qū)中球形金屬顆粒被均勻加熱，加熱功率完全轉(zhuǎn)化為顆粒內(nèi)能；② 球形顆粒之間接觸緊密，沿軸向勻速運動；③ 內(nèi)熱源加熱功率恒定，移動球床通道壁面保持恒定溫度；④ 直徑相同的金屬顆粒在通道中隨機無序填充，孔隙率與顆粒尺寸無關(guān)。

2 顆粒運動速度與溫度分布的關(guān)系

2.1 理論分析

基于球形金屬顆粒在內(nèi)熱源加熱區(qū)中被均勻加熱，且顆粒沿軸向勻速運動，無徑向摻混，以及壁面恒溫等假設(shè)條件，可以將移動球床內(nèi)熱源加熱區(qū)至通道壁面之間的熱量傳遞簡化為通過圓筒壁的導(dǎo)熱問題。在圖1中，含有持續(xù)恒定功率的內(nèi)熱源的移動球床通道沿軸向分為上下兩個部分，上部長L1且含有內(nèi)熱源，下部長L2。上部通道沿徑向可分為兩個部分，中心為內(nèi)熱源加熱區(qū)，半徑為r1，邊界溫度為Tw1；加熱區(qū)至通道壁面之間為壁面導(dǎo)熱區(qū)（r1與r2之間），由于加熱區(qū)與通道壁面之間存在溫差，使得熱量沿徑向由內(nèi)向外傳導(dǎo)。與上部對應(yīng)的可將下部通道沿徑向也分為兩個部分，中心為高溫區(qū)，半徑為r1，邊界溫度為Tw2；同樣，下部r1與r2之間為導(dǎo)熱區(qū)。Tw1、Tw2分別取各自進出口溫度的平均值。在整個過程中，通道壁面溫度Tw保持恒定。加熱區(qū)至通道壁面之間的壁面導(dǎo)熱區(qū)為圓筒形，可視為隨機填充的球形顆粒與氣體組成的多孔介質(zhì)，其徑向有效導(dǎo)熱率為[12]。

在內(nèi)熱源的恒定加熱功率為P的移動球床通道中，溫度為T0的固體球形顆粒自通道上部進入，沿軸向運動距離L1至加熱區(qū)出口，加熱區(qū)中的顆粒受到內(nèi)熱源的均勻加熱使得出口溫度升高為T1。假設(shè)上部沿徑向由內(nèi)向外傳導(dǎo)的熱流量為Φ1，上部壁面導(dǎo)熱區(qū)中顆粒吸收的熱流量為Φ2，則內(nèi)熱源加熱區(qū)中顆粒自入口運動至出口吸收的熱流量為P-Φ1-Φ2；然后，固體球形顆粒繼續(xù)向下運動距離L2至通道出口，通道出口平均溫度為T2。假設(shè)下部高溫區(qū)沿徑向由內(nèi)向外傳導(dǎo)的熱流量為Φ3，下部高溫區(qū)出口處平均溫度為T′2，則下部高溫區(qū)中顆粒自入口運動至出口，因徑向熱傳導(dǎo)使得高溫區(qū)顆粒的溫度由T1下降致T′2。由以上分析可知，顆粒在移動球床通道中自上而下勻速運動，顆粒溫度先升高后降低，因此，對于含有強內(nèi)熱源移動球床來說，最高溫度出現(xiàn)在加熱區(qū)出口處。而加熱區(qū)出口溫度的高低取決于顆粒在加熱區(qū)停留時間的長短，停留時間又與球形顆粒運動速度呈反比。

將含有恒定功率內(nèi)熱源的移動球床通道沿軸向分為上下兩個部分，經(jīng)過理論推導(dǎo)，在移動球床通道上部（內(nèi)熱源加熱區(qū)），加熱區(qū)出口溫度T1（最高溫度）與球形顆粒運動速度v之間的關(guān)系式為

在移動球床通道下部，下部高溫區(qū)出口處平均溫度T′2與球形顆粒運動速度v之間的關(guān)系式為

2.2 理論計算結(jié)果

在本研究中，取移動球床通道尺寸為L=1.4 m，r2= 0.1 m；內(nèi)熱源的恒定加熱功率為P=10 MW，內(nèi)熱源加熱區(qū)尺寸為L1= 0.8 m，r1= 0.05 m；通道入口溫度T0= 473 K，通道壁面溫度Tw= 473K，固體球形顆粒材質(zhì)為鎢，空隙中充滿了氦氣。通道中顆粒孔隙率為0.359～0.44，計算時分別取0.359、0.4和0.44，由式(1)和式(3)可以計算得到加熱區(qū)出口溫度T1、運動速度v，通道出口溫度T2等。

加熱區(qū)出口溫度及通道出口溫度隨顆粒運動速度的變化情況如圖2所示。從圖中可以看出，加熱區(qū)出口溫度T1及通道出口溫度T2隨著顆粒運動速度的增大而降低。這是由于隨著運動速度增大使得顆粒在加熱區(qū)停留時間縮短，在相同加熱功率下，顆粒獲取的熱量減少，因而溫度降低。同時，從圖中還可以看出孔隙率較大的通道中顆粒溫度較高。這是由于通道孔隙率越大即單位體積內(nèi)的顆粒越少，在加熱功率和運動速度相同的條件下，孔隙率較大時單位體積顆粒獲取的熱量較多，因而加熱區(qū)出口溫度較高，通道出口溫度也隨之高于孔隙率較小時的出口溫度。

3 移動球床通道溫度分布數(shù)值模擬

3.1 數(shù)值方法

圖2 理論計算出口溫度隨顆粒運動速度的變化情況Fig.2 Outlet temperature at various particles velocities using theoretical correlation

在Fluent軟件平臺上，對三維通道模型中顆粒與氣體混合流動與傳熱過程進行了數(shù)值計算。物理模型如圖1所示，模型尺寸與理論計算時完全相同。以通道出口截面中心為坐標原點建立坐標系，顆粒由上部入口進入通道，沿z軸向下運動。內(nèi)熱源總加熱功率P為10 MW，通道壁面為恒溫邊界，以入口速度v和溫度T0作為初始條件，選用標準k-ε模型分別對0.426、0.605、0.7 m·s-1，3組入口速度條件下的通道中溫度場的動態(tài)變化過程進行了數(shù)值預(yù)測。經(jīng)計算，通道中溫度場自初始條件起經(jīng) 8 s后趨于穩(wěn)定。數(shù)值計算邊界條件和初始條件與理論計算時相同，通道入口溫度T0= 473 K，通道壁面溫度Tw= 473 K，球形顆粒材質(zhì)為鎢，顆粒空隙間為氦氣。通道中顆粒孔隙率取0.4。

為考察局部內(nèi)熱源移動球床通道中顆粒運動速度與溫度分布，分別提取通道軸向和徑向的關(guān)鍵位置的溫度數(shù)據(jù)進行了對比。在沿顆粒運動方向上，T1為內(nèi)熱源加熱區(qū)出口面平均溫度（z= 0.6 m），T2為通道出口面平均溫度（z= 0）。在沿徑向方向上，Tcentre為移動床通道中心溫度（r= 0）。此外，還獲得了通道總體溫度分布云圖。

3.2 數(shù)值結(jié)果及討論

圖3 含有局部內(nèi)熱源移動球床通道溫度分布Fig.3 Temperature distributions in moving pebble bed with local internal heat generation

3種顆粒速度條件下含有局部內(nèi)熱源移動球床通道溫度穩(wěn)定分布如圖3所示，圖中的虛線表示的是加熱區(qū)的邊界位置。從圖中可以看出，顆粒進入移動床通道后，在強內(nèi)熱源的加熱作用下溫度迅速升高，位于通道中心的加熱區(qū)溫度高于外圍的壁面導(dǎo)熱區(qū)。在通道上部，隨著顆粒在移動球床通道加熱區(qū)中自上而下勻速運動，顆粒溫度不斷升高并在加熱區(qū)出口處出現(xiàn)溫度的最高值。在通道下部，球形顆粒離開內(nèi)熱源加熱區(qū)后中心高溫區(qū)溫度不再升高。隨著顆粒的向下運動，中心高溫區(qū)不斷向外圍傳導(dǎo)熱量溫度逐漸降低，中心與外圍區(qū)域球形顆粒溫度逐漸趨于一致。這與理論分析獲得的溫度分布趨勢相同，在此基礎(chǔ)上，進一步對3種顆粒速度下移動球床通道出口平均溫度的數(shù)值模擬結(jié)果和理論計算結(jié)果進行了對比。結(jié)果發(fā)現(xiàn)，出口平均溫度T2的數(shù)值模擬結(jié)果和理論計算結(jié)果差值ΔT2分別為55.3、20.7、14.5 K，相對誤差分別為6.1%、2.6%、1.9%。由此可以說明，兩者結(jié)果比較吻合，數(shù)值計算所得結(jié)果能夠比較準確地預(yù)測含有局部內(nèi)熱源移動球床通道的溫度場分布。

在含有局部內(nèi)熱源的移動球床中，顆粒的運動速度與通道內(nèi)溫度分布有著必然聯(lián)系。通過數(shù)值模擬計算，對0.426、0.605、0.7 m·s-13組入口速度條件下的通道中溫度分布進行了對比，分析了運動速度與溫度分布之間的關(guān)系。不同顆粒運動速度下加熱區(qū)出口溫度T1及通道出口溫度T2的對比情況如圖4所示，圖中的虛線表示的是加熱區(qū)的邊界位置。從圖4(a)中可以看出，加熱區(qū)出口溫度T1沿通道徑向呈拋物線狀分布，位于通道中部的加熱區(qū)的溫度明顯高于外圍的壁面導(dǎo)熱區(qū)。最高溫度出現(xiàn)在通道徑向中心即通道的軸線上，沿通道徑向自中心至通道壁面的溫度梯度較大。隨著顆粒運動速度的增大，加熱區(qū)出口溫度逐漸降低，徑向溫度梯度也隨之減小。從圖4(b)中可以看出，通道出口溫度T2沿通道徑向同樣呈中部高外圍低的分布，而且隨著顆粒運動速度的增大，通道出口溫度逐漸降低。但是，通道出口溫度分布較為平緩且通道出口最高溫度及徑向溫度梯度明顯低于加熱區(qū)出口。

圖4 數(shù)值計算出口溫度隨顆粒運動速度的變化情況Fig.4 Outlet temperature at various particle velocities using numerical simulation

3種顆粒運動速度下通道中心溫度Tcentre的對比情況如圖5所示，圖中的虛線表示的是加熱區(qū)的出口位置。從圖中可以看出，以加熱區(qū)出口為分界，通道中心溫度Tcentre沿通道軸線方向（顆粒運動方向）先升高后降低，最高溫度出現(xiàn)在加熱區(qū)出口處。由于內(nèi)熱源位于通道上部中心，因此通道中心溫度高于加熱區(qū)邊緣溫度，而且通道上部軸向溫度梯度高于通道下部。隨著球形顆粒運動速度的增大，通道中心溫度隨之降低。

圖5 通道中心溫度隨速度的變化情況Fig.5 Channel centre temperature at various particle velocities

通過對比不同入口速度條件下的通道中溫度分布情況可以得出，對于移動球床來說，在內(nèi)熱源加熱功率一定的條件下提高球形顆粒的運動速度可以減少在加熱區(qū)的停留時間，從而有效減少固體顆粒獲得的熱量，降低顆粒溫度；反之，減小顆粒運動速度可以增加在加熱區(qū)的停留時間，從而增加球形顆粒獲得的熱量，提高顆粒溫度。

4 結(jié) 論

本文對含有強內(nèi)熱源移動球床中顆粒運動速度與溫度分布特性進行了理論分析，結(jié)論如下：含有局部內(nèi)熱源移動球床通道中心溫度沿通道軸線方向（顆粒運動方向）先升高后降低，最高溫度出現(xiàn)在加熱區(qū)出口處。隨著顆粒運動速度的增大，加熱區(qū)出口溫度、通道出口溫度和通道中心溫度逐漸降低。與加熱區(qū)出口溫度分布比較，通道出口溫度分布較為平緩且通道出口最高溫度明顯低于加熱區(qū)出口。孔隙率較大的移動球床通道中，加熱區(qū)出口和通道出口溫度也較高。

[1] Akiyama Tomohiro, Takahashi Reijiro, Yagi Jun-ichiro.Measurements of heat transfer coefficients between gas and particles for a single sphere and for moving beds [J].ISIJ International, 1993,33(6)：703-710.

[2] Amiri A,Vafai K. Analysis of dispersion effects and non-thermal equilibrium, non-Darcian, variable porosity incompressible flow through porous media [J].International Journal of Heat and Mass Transfer, 1994, 37(6)：939-954.

[3] Jiang Peixue,Lu Xiaochen. Numerical simulation of fluid flow and convection heat transfer in sintered porous plate channels [J].International Journal of Heat and Mass Transfer, 2006, 49(9/10)：1685-1695.

[4] Jiang Peixue,Lu Xiaochen. Numerical simulation and theoretical analysis of thermal boundary characteristics of convection heat transfer in porous media [J].International Journal of Heat and Fluid Flow, 2007, 28(5)：1144-1156.

[5] Meng Xianke, Sun Zhongning, Xu Guangzhan. Single-phase convection heat transfer characteristics of pebble-bed channels with internal heat generation [J].Nuclear Engineering and Design, 2012,252：121-127.

[6] Zeisberger Andreas, Mayinger Franz. Void fraction and heat transport in two-dimensional mixed size particle beds with internal heat sources [J].Nuclear Engineering and Design, 2005, 235(20)：2209-2218.

[7] Yang Jian(楊劍 ), Zeng Min(曾敏 ), Yan Xiao(閆 曉), Wang Qiuwang(王秋旺). Numerical study of forced convection in a vertical channel filled with heat-generating porous medium [J].Nuclear Power Engineering(核動力工程), 2009, 30(2)：16-20.

[8] Henda Redhouane, Falcioni Daniel J. Modeling of heat transfer in a moving packed bed：case of the preheater in nickel carbonyl process[J].Journal of Applied Mechanics, 2005, 73(1)：47-53.

[9] Collier A P, Hayhurst A N, Richardson J L, Scott S A. The heat transfer coefficient between a particle and a bed (packed or fluidised) of much larger particles [J].Chemical Engineering Science, 2004, 59(21)：4613-4620.

[10] Lemos Marcelo J S de, Saito Marcelo B. Computation of turbulent heat transfer in a moving porous bed using a macroscopic two-energy equation model [J].International Communications in Heat and Mass Transfer, 2008, 35(10)：1262-1266.

[11] Pivem Ana C, Lemos Marcelo J S de. Laminar heat transfer in a moving porous bed reactor simulated with a macroscopic two-energy equation model [J].International Journal of Heat and Mass Transfer,2012, 55(7/8)：1922-1930.

[12] Lin Ruitai(林瑞泰). Heat and Mass Transfer in Porous(多孔介質(zhì)傳熱傳質(zhì)引論) [M]. Beijing：Science Press, 1994.