邵一窮　李良星　盛天佑

【摘 要】基于美國桑迪亞國家實驗室（Sandia National Laboratory）利用反應堆熔融材料UO2進行堆內(nèi)實驗獲得的碎片床尺寸分布，采用不同尺寸的玻璃圓球和砂石顆粒按照相似的顆粒尺寸分布，分別構(gòu)建了多尺寸球形顆?；旌隙逊e床和砂石堆積顆粒床，實驗研究了單相和氣-水兩相豎直向上流過顆粒堆積床時的流動阻力特性。研究結(jié)果表明，具有寬廣范圍尺寸分布的顆粒堆積床，其有效直徑的大小受顆粒形狀的影響較小；相似尺寸分布的顆粒堆積床的有效直徑十分接近?；趩蜗鄬嶒灪虴rgun方程獲得的多尺寸顆粒堆積床的有效直徑，其兩相流動阻力壓降與Reed模型的計算值吻合較好。

【關(guān)鍵詞】多孔介質(zhì)；堆積顆粒床；流動特性；壓降

中圖分類號： TK124 文獻標識碼： A 文章編號： 2095-2457（2018）06-0001-004

【Abstract】Based on the size distribution of particles from the in-pile experiments with UO2 in Sandia National Laboratory of USA，the test porous beds packed with multi-size spheres and sands particles are designed and constructed，which have the similar size distribution among them.The single and two phase flow tests are performed and the flow resistance characteristics are studied when single and two phase flow upward the packed beds.The results show that for a porous bed packed with various size particles，the shapes of particles have less influences on the effective diameter，and the approximate effective diameters can be obtained from the packed beds with the similar particles sizes distribution.Employing the effective diameter of bed with multi-size particles derived from the Ergun equation，the calculation results from Reed model predict well the measured pressure drops of the two phase flow.

【Key words】Porous media；Packed particles bed；Flow characteristics；Pressure drops

顆粒堆積多孔介質(zhì)通道內(nèi)的單相/兩相流動現(xiàn)象在許多工程和科學領(lǐng)域中都有著廣泛的應用和研究，涉及包括農(nóng)業(yè)技術(shù)，生物工程，機械工程，石油化工工程、核動力工程等多個領(lǐng)域[1-3]。特別是在輕水反應堆嚴重事故進程中，具有多孔介質(zhì)結(jié)構(gòu)的顆粒堆積碎片床可能會在反應堆內(nèi)不同位置（如壓力容器下封頭等）產(chǎn)生，由于碎片床會產(chǎn)生衰變熱，將嚴重威脅反應堆的安全。因此研究冷卻水在堆積顆粒通道內(nèi)的流動特性，進而分析其對碎片床的冷卻性，對抑制和緩解核反應堆嚴重事故的進一步發(fā)展具有十分重要的工程價值和科學意義[3，4]。

1 理論模型

許多研究者已針對顆粒堆積多孔介質(zhì)通道內(nèi)的流動特性開展了深入的實驗[5-8]和數(shù)值計算研究[9-11]，并提出了不同的流動阻力模型。Liu et al.[12]，閆曉等[1]，Lindholm et al.[13]，Clavier et al.[14]等先后對國內(nèi)外的摩擦阻力預測模型進行了總結(jié)。其中，Ergun方程[15]被普遍接受和廣泛應用于單尺寸球形顆粒堆積床內(nèi)單相流動的摩擦阻力預測分析中

式中，ΔP是壓降；L是多孔介質(zhì)通道的高度；？滋是流體動力粘度；ρ是流體的密度； 150和1.75是基于實驗而得到的經(jīng)驗常數(shù)；ε是多孔介質(zhì)的孔隙率；d是組成多孔介質(zhì)床的顆粒有效直徑；J是多孔介質(zhì)內(nèi)流體的表觀速度；K和η分別被稱為滲透率和穿透率。隨著兩相流動研究的發(fā)展，Lipinski[16]等研究學者在Ergun方程的基礎上，通過引入相對滲透率Kr和相對穿透率ηr，擴展了Ergun方程在兩相流動中的阻力預測的應用。

其中，Jl和Jg分別代表兩相流動中液相和氣相的表觀速度；α是流道內(nèi)的截面含氣率，F(xiàn)i是相間摩擦力。在實驗基礎上，Reed[17]，Schulenberg and Mller[18]，Hu and Theofanous[19]等研究人員先后提出了各自的兩相流動阻力模型基本參數(shù)。表1列出了上述研究人員提出的基本參數(shù)，其中，s表示多孔介質(zhì)結(jié)構(gòu)單元內(nèi)的飽和度，s=1-α。

由上述模型可以看出，堆積顆粒的有效直徑d是多孔介質(zhì)通道內(nèi)單相、兩相流動壓降計算的重要參數(shù)。對單尺寸球形顆粒組成的多孔介質(zhì)通道，顆粒直徑即為其有效直徑。然而反應堆嚴重事故進程中形生的碎片床是由具有一定尺寸分布的多種碎片顆粒混合堆積而成[20，21]。因此，獲得多尺寸顆粒混合堆積結(jié)構(gòu)的有效直徑，進而驗證研究通道內(nèi)流動阻力分析模型，對反應堆嚴重事故碎片床冷卻性研究具有重要的學術(shù)意義。

表1 多孔介質(zhì)通道內(nèi)兩相流動摩擦阻力模型參數(shù)

Table.1 Parameters in Frictional Resistance models of two-phase flow in a porous media

為研究多尺寸顆?；旌隙逊e多孔介質(zhì)通道內(nèi)單相、兩相流動特性，本文基于美國桑迪亞國家實驗室（Sandia National Laboratory）利用反應堆熔融材料UO2進行堆內(nèi)實驗獲得的碎片床尺寸分布，采用不同尺寸的玻璃圓球和砂石按照相似的顆粒尺寸分布，分別構(gòu)建了多尺寸球形顆粒混合堆積床和砂石堆積顆粒床，實驗研究了單相和氣-水兩相豎直向上流過顆粒堆積床時的流動阻力特性。

2 實驗系統(tǒng)及實驗描述

2.1 實驗裝置

為研究顆粒堆積多孔介質(zhì)通道內(nèi)的單相/兩相流動特性，西安交通大學動力工程多相流國家重點實驗室設計和搭建了顆粒堆積床流動特性實驗系統(tǒng)，如圖1所示。

整個系統(tǒng)由供水系統(tǒng)、供氣系統(tǒng)、汽水混合段、實驗測試段及實驗數(shù)據(jù)測量和采集系統(tǒng)構(gòu)成。其中實驗測試段內(nèi)徑120 mm，高度600 mm。在測試段的底部和頂部的法蘭連接處，分別放置不銹鋼網(wǎng)柵，用于支撐和固定堆積顆粒床。實驗過程中的溫度、壓力等信號由NI公司生產(chǎn)的數(shù)據(jù)采集系統(tǒng)進行收集和整理。

2.2 實驗顆粒床

在本研究中，多種不同尺寸的玻璃球和砂石分別被填充到實驗測試段中，進行單相和兩相豎直向上流動實驗，進而研究分析流動摩擦壓降及其變化規(guī)律。表2列出了實驗顆粒床信息，顆粒尺寸分布如圖2所示。其中，圖2中DCC-2實驗[21]顆粒尺寸分布是美國桑迪亞國家實驗室（Sandia National Laboratory）基于反應堆熔融材料UO2進行堆內(nèi)實驗（In-Pile Exp.）獲得的碎片床尺寸分布，而本次實驗使用的球形顆粒堆積床和砂石顆粒堆積床的尺寸分布均是按照DCC-2實驗[21]顆粒尺寸分布進行均勻混合堆積。

其中，Mj是裝入實驗段顆粒j的總重量，ρj是其密度，V0是實驗測試段的體積。

3 實驗結(jié)果與分析

3.1 單相流動實驗結(jié)果及分析

圖3和圖4分別顯示了單相水豎直向上流過球形顆粒堆積床和砂石顆粒堆積床時測量得到的阻力壓降，圖中以“Δ”符號表示。由圖中可以看出，隨著流速升高，兩種堆積顆粒床的單相阻力壓降均逐漸增大。

基于Ergun方程和測量的阻力壓降，砂石顆粒床的有效直徑可以計算得到，為1.44 mm。圖3和圖4中也給出了使用1.5 mm和1.44 mm為堆積顆粒床的有效直徑時Ergun方程的計算值，以黑實線表示。由圖中可以看出，在實驗工況范圍內(nèi)計算值和實驗值吻合的很好。DCC-2實驗[21]中獲得的碎片床的有效直徑和孔隙率分別為1.42 mm和0.41，Ergun方程也計算預測了單相水流過UO2顆粒床時的阻力壓降，如圖4中虛線所示。由圖4可以看出，砂石堆積顆粒床的單相阻力壓降與UO2顆粒堆積床的預測計算值很接近，最大處的相對偏差僅為7.6%。

表3列出了上述三種顆粒堆積床的孔隙率和有效直徑，包括DCC-2實驗[21]形成的顆粒床及本文的球形顆粒堆積床和砂石顆粒堆積床?？傮w上，三種顆粒床具有相似的顆粒尺寸分布（如圖2所示），分別以具有不同形狀的球形顆粒和砂石顆粒進行模擬反應堆材料堆積碎片床時，獲得的有效直徑與反應堆材料堆積碎片床的有效直徑差值很小（球形顆粒床差值最大，為0.8 mm），相對偏差僅為5.6%。由表3對比結(jié)果可以說明，具有寬廣范圍尺寸分布的顆粒堆積床，其有效直徑的大小受顆粒形狀的影響較?。幌嗨频某叽绶植碱w粒床具有接近的有效直徑。

表3 具有相似尺寸分布的三種顆粒床的孔隙率和有效直徑

Table 3 Porosities and effective diameters of three particles beds with the similar size distribution

進一步對比圖3和圖4中測量獲得的阻力壓降可以發(fā)現(xiàn)，相同條件下，圖3中球形顆粒堆積床的單相阻力壓降遠高于圖4中砂石顆粒堆積床，這主要因為砂石顆粒具有不規(guī)則形狀，受其影響，不同尺寸砂石顆?；旌隙逊e床的孔隙率較高（表3中為0.4），而球形顆?；旌隙逊e床的空隙率較低（表3中為0.33）。顯然，有效直徑差別不大的情況下，較低的空隙率將導致更高的阻力壓降。

3.2 兩相流動實驗結(jié)果及分析

在單相流動特性研究的基礎上，本文也進行了氣-水兩相流動實驗。圖5和圖6分別顯示了氣-水兩相豎直向上流過球形顆粒堆積床和砂石堆積顆粒床時測量的壓降梯度，圖中以“Δ”符號表示。為驗證兩相流模型的準確性，圖中也分別給出了表1中不同兩相模型的兩相阻力壓降計算預測值，其中顆粒直徑采用單相流動實驗獲得的有效直徑。圖中，實線表示Lipinski模型[16]計算預測值，虛線表示Reed模型[17]，點劃線表示Schulenberg & Muller 模型[18]，雙點劃線表示Hu & Theofanous 模型[19]的計算預測值。

由圖5和圖6可以看出，當液速保持恒定時，隨氣速增大，堆積顆粒床的兩相阻力壓降逐漸升高?？傮w上，利用單相流動實驗和Ergun方程獲得的有效直徑，Reed模型的計算預測值與實驗值吻合的較好，尤其是流速較高的時候。

圖7對比了相同的液速和氣速條件下，球形顆粒堆積床和砂石顆粒堆積床的兩相阻力壓降。由圖7可以看出，同等條件下，砂石顆粒堆積床的兩相阻力壓降遠低于球形顆粒堆積床；且隨著氣速增大，兩者的差值越大。研究表明[22]，顆粒直徑越小，孔隙率越低，流體經(jīng)過堆積顆粒床的阻力壓降越高。由表3可知，相對于球形顆粒床，砂石顆粒床具有更小的有效直徑（1.44 mm）和更高的孔隙率（0.4）。分析認為，對寬廣尺寸分布的顆粒堆積床，其有效直徑的大小受顆粒形狀的影響較小，與具有相似尺寸分布的球形顆粒床相比，雖然砂石顆粒具有不規(guī)則形狀，但砂石堆積顆粒床的有效直徑僅僅減小了4%；而其孔隙率卻增加了21.2%。因此堆積顆粒的不規(guī)則形狀主要提高了砂石顆粒堆積床的孔隙率，進而產(chǎn)生了較低的兩相流動阻力。

4 結(jié)論

本文基于美國桑迪亞國家實驗室（Sandia National Laboratory）利用反應堆熔融材料UO2進行堆內(nèi)實驗獲得的碎片床尺寸分布，采用不同尺寸的玻璃圓球和砂石按照相似的顆粒尺寸分布，分別構(gòu)建了多尺寸球形顆?；旌隙逊e床和砂石堆積顆粒床，實驗研究了單相和氣-水兩相豎直向上流過顆粒堆積床時的流動阻力特性。研究結(jié)果表明：

（1） 具有寬廣范圍尺寸分布的顆粒堆積床，其有效直徑的大小受顆粒形狀的影響較??；相似尺寸分布的顆粒堆積床的有效直徑十分接近。

（2） 基于單相實驗和Ergun方程獲得的多尺寸顆粒堆積床的有效直徑，其兩相流動阻力壓降與Reed模型的計算值吻合較好。

【參考文獻】

[1]閆曉，肖澤軍，黃彥平，等.多孔介質(zhì)中流動換熱特性的研究進展[J].核動力工程，2006，27（1）：77-82.YAN Xiao， XIAO Zejun， HUANG Yanping， et. al. Research Progress on Flow and Heat Transfer in Porous Media[J]. Nuclear Power Engineering， 2006， 27（1）： 77-82. （in Chinese）.

[2]Jamialahadi M， Müller-Steinhagen H， Izadpanah M. Pressure Drop， Gas Hold-up and Heat Transfer during Single and Two-phase Flow through Porous Media[J]. International Journal of Heat and Fluid Flow， 2005， 26（1）： 156-172.

[3]Sehgal B R.， Stabilization and termination of severe accidents in LWRs [J]. Nuclear Engineering & Design， 2006， 236（19–21）：1941-1952.

[4]巴爾·拉吉·塞加爾主編. 輕水堆核安全[M]. 馬衛(wèi)民，趙博等譯. 輕水堆核安全：嚴重事故現(xiàn)象學[M] 北京，中國原子能出版社，2015， 206-234. Sehgal B R. Nuclear safety in light water reactor[M]. Translated by Ma WM， Zhao B. Nuclear safety in light water reactor： severe accident phenomenology [M]. Beijing： China atomic energy press， 2015， 206-234. （in Chinese）.

[5]Chikhi N， Clavier R， Laurent JP， etal. Pressure drop and average void fraction measurements for two-phase flow through highly permeable porous media， Annals of Nuclear Energy. 2016， 94（8）： 422-432.

[6]Takasuo E. An experimental study of the coolability of debris beds with geometry variations， Annals of Nuclear Energy. 2016， 92（6）：251-261.

[7]李良星，鄒旭毛，孔劉波.不同堆積結(jié)構(gòu)顆粒床內(nèi)的流動特性研究[J].西安交通大學學報，2016，50（9）：49-54.Li， Liangxing， Zou Xumao， Kong Luibo. Investigation on the Flow Characteristics in Particulate Beds with Different Packing Structures [J]， Journal of XiAn Jiaotong University， 2016， 50（9）： 49-54. （in Chinese）

[8]張楠，孫中寧，趙忠南.球床通道內(nèi)氣液兩相豎直向上流動流型試驗研究[J].原子能科學技術(shù)，2011，45（9），1040-1044.Zhang Nan， SUN Zhongning， ZHAO Zhongnan. Experimental investigation on Flow Patterns of Gas-liquid Two-phase upward Flow through Packed Channel with spheres [J]. Atomic Energy Science and Technology， 2011， 45 （9）：1040-1044. （in Chinese）

[9] 于立章， 孫立成， 孫中寧. 多孔介質(zhì)通道單相流阻力特性數(shù)值模擬[J]. 原子能科學技術(shù)， 2010， 44 （12）：1441-1444.

YU Lizhang， SUN Licheng， SUN Zhongning. Numerical Simulation of Single-Flow Resistance in Porous Media Channel [J]， Atomic Energy Science and Technology， 2010， 44 （12）：1441-1444. （in Chinese）

[10] 李健，宋曉明，魯劍超，等. 基于有序堆積球床的單相阻力特性實驗研究與數(shù)值模擬[J]. 原子能科學技術(shù)， 2012， 46 （Suppl.）： 807-810.LI Jian， SONG Xiaoming， LU Jiangchao， et.al. Experimental Research and Numerical Simulation on Single Phase Flow Resistance Characteristics of Orderly Packed Bed [J]， Atomic Energy Science and Technology， 2012， 46 （Suppl.）： 807-810. （in Chinese）

[11]Tung VX， Dhir， VK. A hydrodynamic model for two-phase flow through porous media. International Journal of Multiphase Flow， 1988， 14（1）： 47-65.

[12]Liu， S.， Afacan， A.， Masliyah， J.H. Steady incompressible laminar flow in porous media [J]. Chemical Engineering Science， 1994， 49（21）： 3565-3586.

[13]Lindholm， I.， Holmstr？m， S.， Miettinen， J.， et al. Dryout heat flux experiments with deep heterogeneous particle bed [J]. Nuclear Engineer and Design， 2006， 236（19-21）： 2060-2074.

[14]Clavier R.， Chikhi N.， Fichot F.， et al. Modeling of Inertial Multi-phase Flows through High Permeability Porous Media： Friction Closure Laws [J]. International Journal of Multiphase Flow， 2017， 91： 243-261.

[15]Ergun， S. Fluid flow through packed columns. Chemical Engineering Progress， 1952， 48（2）： 89-94.

[16]Lipinski， R.J. A one dimensional particle bed dryout model. ANS Transactions， 1981， 38： 386-387.

[17]Reed， A.W. The effect of channeling on the dryout of heated particulate beds immersed in a liquid pool. PhD thesis， Massachusetts Institute of Technology， Cambridge， USA， 1982.

[18]Schulenberg， T. and M？ller， U. An improved model for two-phase flow through beds of course particles[J]. International Journal of Multiphase flow， 1987， 13 （1）： 87-97.

[19]Hu， K. and Theofanous， T.G. On the measurement and mechanism of dryout in volumetrically heated coarse particle beds[J]. International Journal of Multiphase Flow， 1991， 17 （4）： 519-532.

[20]Ma，W. and Dinh， T.N.. The effects of debris bed's prototypical characteristics on corium coolability in a LWR severe accident [J]. Nuclear Engineering and Design， 2010， 240：598-608.

[21]Reed， A.W.， Bergeron， E.D.， Boldt， K.R. and Schmidt， T.R. Coolability of UO2 debris beds in pressurized water pools： DCC-1 and DCC-2 experiment results [J]. Nuclear Engineering and Design， 1986， 97： 81-88.

[22]李良星，李會雄. 多尺寸顆粒堆積多孔介質(zhì)床有效直徑研究[J]. 工程熱物理學報，2014， 35（9）： 1785-1788.LI Liangxing， LI Huixiong. Study on the effective diameter of Porous bed packed with multi-size particles [J]. Journal of Engineering Thermophysics， 2014， 35（9）： 1785-1788. （in Chinese）.