含內熱源多孔介質通道內流動沸騰兩相壓降的預測

張小寧，孫中寧，孟現珂，徐廣展

（哈爾濱工程大學 核安全與仿真技術國防重點學科實驗室，黑龍江 哈爾濱 150001）

含內熱源多孔介質通道內流動沸騰兩相壓降的預測

張小寧，孫中寧，孟現珂，徐廣展

（哈爾濱工程大學 核安全與仿真技術國防重點學科實驗室，黑龍江 哈爾濱 150001）

為預測含內熱源多孔介質通道內流動沸騰的兩相阻力壓降，以Ergun方程為基礎，定義了多孔介質通道的Chisholm參數Y和全液相折算因子Φl0。通過理論分析和實驗數據測量，明確了出口質量含氣率xe、質量流速G和小球直徑d等對參數Y和Φl0的影響，并提出1個Lockhart-Martinelli（L-M）類型的兩相阻力壓降關系式。與文獻中的其他公式相比，本工作提出的公式對實驗結果能做出更好的預測。

兩相壓降；流動沸騰；多孔介質；內熱源

多孔介質通道內氣液兩相的流動特性受到石油化學工業、核能開發利用、地熱資源開采和電子組件冷卻等領域的廣泛關注，對流動過程中的兩相阻力壓降進行預測是這些領域中重要的研究課題之一。閆曉等［1］和廖永浩［2］對多孔介質內兩相流動的研究進行了較全面的探討，對兩相流動阻力壓降的計算模型和預測公式進行了總結。已有的多孔介質內兩相流動的理論模型和計算公式絕大部分是基于絕熱條件下空氣和水的實驗數據建立的，不涉及相態的變化。張楠等［3-4］對絕熱多孔介質通道內氣液兩相流的阻力特性進行了實驗研究，并對兩相流型做了可視化觀察，建立了兩相壓降與流型之間的對應關系。目前，關于多孔介質通道內流動沸騰過程中兩相阻力壓降的研究并不充分。Naik等［5］針對水在直徑為0.59～0.79mm的鋼球形成的多孔介質床內受熱沸騰的情形，提出了一計算兩相流動總壓降的經驗關系式。孟現珂等［6-8］對含內熱源球床通道內的單相對流換熱特性進行了系統地實驗研究。

本文在文獻［6-9］基礎上，利用直徑為5mm和12mm的表面氧化碳鋼球填充形成多孔介質通道，采用電磁感應的方式對球床整體加熱，在常壓下獲取含內熱源多孔介質通道內水的流動沸騰壓降數據，以此為基礎提出1個Lockhart-Martinelli（L-M）類型的兩相阻力壓降預測關系式。

1 阻力壓降的基礎預測模型

多孔介質內的單相阻力壓降通常使用Ergun［10］方程進行預測，該方程綜合考慮了流動過程中的黏性效應和慣性效應對阻力壓降的影響：

其中：Δp為阻力壓降，MPa；ΔL為測壓段長度，m；ε為孔隙率；d為小球直徑，m；μ為流體的動力黏度，Pa·s；ρ為流體的密度，kg·m-3；u為表觀流速，m·s-1；A和B為實驗確定的經驗常數。

多孔介質通道的結構參數列于表1。

表1 多孔介質通道的參數Table 1 Parameters of porous media channel

對式（1）做如下變形：

其中：G為表觀質量流速，kg·m-2·s-1；Rem為修正雷諾數。

令：

由式（2）～（4）可得到：

兩相流中氣相或液相單獨流過多孔介質通道的阻力壓降梯度為：

式中，x為兩相流動系統的流動質量含氣率。

與兩相總質量流速相等的等效氣相或液相通過多孔介質通道的阻力壓降梯度為：

2 Chisholm參數Y

Larkins等［11］和Midoux等［12］采用與Lockhart-Martinelli處理常規通道中兩相摩擦壓降數據類似的方法，在預測絕熱條件下多孔介質通道內兩相流動阻力壓降時引入分液相折算因子Φl和Martinelli參數X，建立了Φl與X的關系式。Φl和X分別定義為：

式中，（Δp／ΔL）tp為兩相阻力壓降梯度。

聯立式（7）、（8）、（12）可得：

X體現了流動質量含氣率、流體物性和多孔介質結構參數等對兩相阻力壓降的影響。在液體受熱沸騰的兩相流動系統中，因x沿流動方向不斷變化，不便于得到確定的x來計算X。因此，Chisholm針對常規通道中的流動沸騰定義了新的參數Y［13］：

多孔介質通道內流動沸騰的Chisholm參數Y可通過聯立式（9）、（10）、（14）得到：

其中：

對于本實驗的Ergun常數A和B，函數f（Rem）＝lg（ARe-1m＋B）／lg Rem在Rem＞0的范圍內是單調變化的，即只有當μl與μg相等時才會使kl0等于kg0。在低于臨界壓力的情況下，μg不等于μl，則kl0和kg0不相等，于是由式（15）定義的Chisholm參數Y不僅表征了流體的物性，其同時也受到質量流速G及多孔介質的結構參數d和ε的影響。圖1示出質量流速G和小球直徑d對Chisholm參數Y的影響。由圖1可見，Chisholm參數Y隨質量流速的增加而增大，隨小球直徑的增大而減小。

圖1 G和d對Chisholm參數Y的影響Fig.1 Effects of Gand don Chisholm parameter Y

3 兩相流動阻力壓降預測關系式

在液體受熱沸騰的兩相流動系統中，質量含氣率x沿流動方向不斷變化，絕熱狀態下定義的分液相折算因子Φl的值不便于確定，此時使用全液相折算因子Φl0更為方便。參照Martinelli-Nelson［14］、Thom［15］及Chisholm［13］有關預測常規通道中流動沸騰兩相摩擦壓降的研究，對Φl0定義如下：

利用實驗結果計算得到的Φl0示于圖2。由圖2可看出，多孔介質通道內流動沸騰的全液相折算因子Φl0受到通道出口質量含氣率xe、流體的質量流速G及形成多孔介質的小球直徑d的影響。通道出口質量含氣率增加時，全液相折算因子增大；質量流速增加時，全液相折算因子減小。球徑對全液相折算因子的影響與出口質量含氣率有一定的關系。在圖2b所示的出口質量含氣率范圍內，球徑增加，全液相折算因子增大。

圖2 出口質量含氣率與全液相折算因子的關系Fig.2 Φl0vs.xe

綜合上述各影響因素，全液相折算因子Φl0可表示為xe、G和d的函數，即：

考慮到式（14）定義的Chisholm參數Y不僅表征了流體的物性，同時也體現了質量流速G及多孔介質的結構參數d和ε的影響，參照Chisholm兩相壓降預測關系式［13］，將式（19）具體表示為如下形式：

其中，Γ為參數。

圖3示出不同Y時Г與xe的關系。由圖3可見，參數Г與xe之間呈線性變化關系，且不同Chisholm參數Y對應不同的直線，故將Г表示成Y和xe的冪指數的形式，即：

圖3 不同Y時Г與xe的關系Fig.3 Гvs.xefor different Chisholm parameter Y

其中，a、m、n為通過實驗數據確定的經驗常數。

對實驗數據回歸擬合得到，a＝1.628× 106，m＝-5.832，n＝-1.022。式（21）的適用范圍為：12.3kg·m-2·s-1＜G＜38.1kg· m-2·s-1，0＜xe＜0.53。

利用式（20）和（21）計算得到的全液相折算因子Φl0與本文實驗值的比較示于圖4，其中小球直徑d＝12mm，質量流速G＝37.9kg· m-2·s-1，圖4中同時示出Naik關系式［5］和植田辰洋關系式［16］的計算結果。Naik關系式是基于顆粒直徑為0.59～0.79mm的球體堆積床在質量流速小于2.78kg·m-2·s-1情況下的沸騰壓降數據得到的，其計算值比本文實驗值大1～2個數量級。植田辰洋關系式是計算常規通道中流動沸騰摩擦壓降常用的公式，其計算結果與本文實驗值最大相差160%。本工作提出的式（20）、（21）對實驗結果做出了很好的預測，相對偏差在±20%以內。

圖4 公式預測值與本文實驗值的比較Fig.4 Comparison of evaluated and experimental values

4 結論

1）多孔介質通道的Chisholm參數Y不僅表征了流體的物性，同時也受到質量流速G及多孔介質的結構參數d和ε的影響。

2）全液相折算因子Φl0隨出口質量含氣率的增加而增大，隨質量流速的增加而減?。磺驈皆黾?，全液相折算因子增大。

3）利用本工作提出的Lockhart-Martinelli類型公式能對含內熱源多孔介質通道內流動沸騰的兩相阻力壓降做出很好的預測，與實驗值的相對偏差在±20%以內。

［1］ 閆曉，肖澤軍，黃彥平，等.多孔介質中流動換熱特性的研究進展［J］.核動力工程，2006，27（1）：77-78.

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-78（in Chinese）.

［2］ 廖永浩.微球床多孔介質通道氣液兩相流特性研究［D］.哈爾濱：哈爾濱工程大學，2009.

［3］ 張楠，孫中寧.多孔介質通道內氣-液兩相流動阻力特性實驗［J］.核動力工程，2011，32（3）：106-110，126.

ZHANG Nan，SUN Zhongning.Experimental investigation on resistance characteristics of twophase flow through porous media channel［J］.Nuclear Power Engineering，2011，32（3）：106-110，126（in Chinese）.

［4］ 張楠，孫中寧，趙忠南.球床通道內氣液兩相豎直向上流動流型實驗研究［J］.原子能科學技術，2011，45（9）：1 040-1 044.

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）：1 040-1 044（in Chinese）.

［5］ NAIK A S，DHIR V K.Forced flow evaporative cooling of a volumetrically heated porous layer［J］.Int J Heat Mass Transfer，1982，25（4）：541-552.

［6］ 孟現珂，孫中寧，周平，等.強內熱源球床通道單相對流換熱特性實驗研究［J］.原子能科學技術，2012，46（9）：1 061-1 066.

MENG Xianke，SUN Zhongning，ZHOU Ping，et al.Experimental study on single-phase convection heat transfer characteristics of pebble bed channels with internal heat generation［J］.Atomic Energy Science and Technology，2012，46（9）：1 061-1 066（in Chinese）.

［7］ 孟現珂，孫中寧，徐廣展，等.含內熱源堆積球床對流換熱特性的實驗研究［J］.哈爾濱工程大學學報，2012，33（9）：1 122-1 126.

MENG Xianke，SUN Zhongning，XU Guangzhan，et al.Experimental study on convection heat transfer characteristics of stacked pebble-bed channels with an internal heat source［J］.Journal of Harbin Engineering University，2012，33（9）：1 122-1 126（in Chinese）.

［8］ 孟現珂，孫中寧，鄧振國，等.含內熱源球床通道換熱特性實驗研究［J］.核動力工程，2012，33（2）：71-74.

MENG Xianke，SUN Zhongning，DENG Zhenguo，et al.Experimental study on heat transfer characteristics of pebble-bed channels with inter-nal heat generation［J］.Nuclear Power Engineering，2012，33（2）：71-74（in Chinese）.

［9］ 徐廣展，孫中寧，孟現珂，等.含內熱源多孔介質通道內ONB特性［J］.化工學報，2012，63（10）：3 080-3 084.

XU Guangzhan，SUN Zhongning，MENG Xianke，et al.ONB characteristics in porous media with internal heat source［J］.CIESC Journal，2012，63（10）：3 080-3 084（in Chinese）.

［10］ERGUN S.Fluid flow through packed columns［J］.Chemical Engineering Progress，1952，48（2）：89-94.

［11］LARKINS R P，WHITE R R.Two-phase concurrent flow in packed beds［J］.ALChE Journal，1961，7（2）：231-239.

［12］MIDOUX N，FAVIER M，CHARPENTTER J C.Flow pattern，pressure loss and liquid hold up data in gas-liquid down flow in packed-beds with foaming and non forming hydrocarbons［J］.Journal of Chemical Engineering of Japan，1976，9（5）：350-356.

［13］CHISHOLM D.Pressure gradients due to friction during the flow of evaporating two-phase mixtures in smooth tubes and channels［J］.Int J Heat Mass Transfer，1973，16：347-358.

［14］MARTINELLI R C，NELSON D B.Prediction of pressure drop during forced-circulation boiling of water［J］.Trans Amer Soc Mech Engrs，1948，70：695-702.

［15］THOM J R S.Prediction of pressure drop during forced circulation boiling of water［J］.Int J Heat Mass Transfer，1964，7：709-724.

［16］閻昌琪.氣液兩相流［M］.哈爾濱：哈爾濱工程大學出版社，2010：74.

Prediction of Two-phase Pressure Drop During Flow Boiling in Porous Media Channel With Internal Heat Source

ZHANG Xiao-ning，SUN Zhong-ning，MENG Xian-ke，XU Guang-zhan
（Fundamental Science on Nuclear Safety and Simulation Technology Laboratory，Harbin Engineering University，Harbin150001，China）

Based on Ergun equation，the Chisholm parameter Yand two-phase multiplier Φl0for porous media channel were defined to predict two-phase pressure drop during flow boiling in porous media channel with internal heat source.Effects of exit quality xe，mass flux Gand particle diameter d on parameters Y andΦl0were identified by theoretical analysis and experimental measurements.A Lockhart-Martinelli（L-M）type correlation for predicting two-phase pressure drop was proposed and it can make more reasonable predictions to experimental results than other correlations available in literature.

two-phase pressure drop；flow boiling；porous media；internal heat source

TL333

A

1000-6931（2014）02-0241-05

10.7538／yzk.2014.48.02.0241

2012-11-12；

2013-01-06

國家自然科學基金資助項目（11075042）

張小寧（1989—），男，山西大同人，碩士研究生，核能科學與工程專業