鐘達(dá)文，孟繼安，李志信

?

朝下溝槽結(jié)構(gòu)表面池沸騰換熱

鐘達(dá)文，孟繼安，李志信

（清華大學(xué)航天航空學(xué)院，工程熱物理研究所，北京 100084）

核電站熔融物堆內(nèi)滯留技術(shù)是一項關(guān)鍵的嚴(yán)重事故應(yīng)對策略，該策略已被核電站廣泛采用。為增強核電站壓力容器下封頭外表面的沸騰換熱能力，實驗研究了常壓下朝下溝槽結(jié)構(gòu)表面的池沸騰換熱，測量了傾角5°、30°、45°、60°和90°下熱通量隨壁面過熱度的變化，獲得了相應(yīng)傾角下的臨界熱通量（CHF）。與光表面相比，朝下溝槽結(jié)構(gòu)表面的CHF可提高65%～90%。實驗觀察發(fā)現(xiàn)，在高熱通量下朝下溝槽結(jié)構(gòu)表面氣泡運動形態(tài)存在蒸汽膜和波浪蒸汽層兩種結(jié)構(gòu)。分析表明，沸騰換熱顯著增強、臨界熱通量大幅提高的原因是溝槽結(jié)構(gòu)大幅增加了換熱面積同時還明顯改善了表面的潤濕性。

傳熱；表面；相變；氣泡；臨界熱通量；傾角

引 言

自美國三哩島核事故發(fā)生后，研究人員發(fā)現(xiàn)，將嚴(yán)重事故條件下的壓力容器浸沒在冷卻水中能夠?qū)崿F(xiàn)堆內(nèi)熔融物滯留壓力容器的目的，由此提出了堆內(nèi)熔融物保持策略（in-vessel retention，IVR）。目前在壓水堆核電站中，IVR已經(jīng)成為緩解嚴(yán)重事故的一項關(guān)鍵措施。在核電站發(fā)生嚴(yán)重事故情況下，通過非能動方式向反應(yīng)堆腔室充水，使壓力容器底部下封頭淹沒，在壓力容器外部發(fā)生池沸騰換熱冷卻壓力容器內(nèi)部熔融物，是實現(xiàn)IVR的一項重要策略[1-4]。保證熱量完全由池沸騰換熱帶走的關(guān)鍵在于使朝下?lián)Q熱表面的熱通量不超過臨界熱通量。因此，提高朝下表面沸騰換熱的臨界熱通量、增強壓力容器外表面的沸騰換熱就顯得格外重要。

近年來，朝下表面池沸騰換熱引起了廣大研究人員的關(guān)注。Vishnev[5]、El-Genk等[6]、Brusstar 等[7]、Kim等[8]、Rouge[9]實驗研究了一維朝下平板表面沸騰換熱狀況；Theofanous等[10-11]研究了一維朝下曲面的沸騰換熱狀況；Chu等[12-13]、Haddad[14]、Yang等[15]和Noh等[16-17]則實驗研究了三維朝下半球下封頭外表面的沸騰換熱狀況。研究表明，臨界熱通量隨傾角的增大而增加（換熱表面豎直布置時的傾角定義為90°，水平朝下布置時的傾角為0°）。Yang等[15]研究了多孔涂層對朝下曲面的沸騰強化，表明涂層能較好地強化沸騰換熱，提高臨界熱通量。然而多孔涂層技術(shù)目前還存在性能不穩(wěn)定、壽命短、可靠性不足等缺陷，對于核電站設(shè)計壽命60年的情況，采用涂層強化朝下表面沸騰換熱存在較大的不確定性。因此，為了增強IVR的安全裕度，迫切需要開發(fā)出性能優(yōu)良、可靠性高、壽命長的強化沸騰換熱的朝下結(jié)構(gòu)表面。本文實驗研究了朝下溝槽結(jié)構(gòu)表面的沸騰換熱性能，測量了不同傾角下沸騰換熱的臨界熱流。

1 實驗裝置和實驗方法

實驗裝置由實驗件系統(tǒng)、操作條件控制系統(tǒng)、數(shù)據(jù)采集系統(tǒng)等組成，圖1為實驗裝置照片。實驗采用電加熱的方法提供實驗件沸騰換熱所需熱量，通過調(diào)節(jié)電加熱器功率控制換熱表面的熱通量，實驗通過觀察換熱表面沸騰時的氣泡運動形態(tài)和測量實驗件表面的溫度變化來判斷沸騰危機是否發(fā)生。

實驗件系統(tǒng)包括可旋轉(zhuǎn)水箱和實驗組件，其中可旋轉(zhuǎn)水箱設(shè)置有5個觀察窗，可以方便地在任意傾角觀察實驗件表面的沸騰狀況。實驗組件密閉安裝在旋轉(zhuǎn)水箱中并可360°自由旋轉(zhuǎn)，實驗組件包含實驗件和實驗件隔熱固定組件，實驗件為帶電加熱棒的方形紫銅塊（100 mm×100 mm×75 mm），如圖2所示，內(nèi)插24根長136 mm、直徑9.22 mm、總功率為24 kW的電加熱器，24根電加熱器分4 組均勻布置以保證熱通量均勻。通過晶閘管整流器整流電源，調(diào)節(jié)電加熱器的電壓，從而控制實驗件換熱表面的熱通量。實驗件的一個平面即實驗表面浸沒在水中，而其他表面密封在密閉空間并與水隔熱，為確保實驗中沸騰傳熱僅發(fā)生在實驗表面，在實驗件周邊用隔熱密封件包裹。同時將實驗銅塊、隔熱密封件和其他密封件一起固定，通過旋轉(zhuǎn)轉(zhuǎn)動機構(gòu)就能自由控制平板的傾角。在離換熱表面垂直距離6.5 mm處，銅塊兩側(cè)面各布置5個直徑0.5 mm的K型鎧裝熱電偶，如圖3所示。通過傅里葉導(dǎo)熱定律計算獲得換熱表面的平均壁溫；2個直徑1 mm的熱電偶布置在換熱表面附近的去離子水中，用于測量水溫；2個直徑1 mm的熱電偶布置在銅塊頂部，用于監(jiān)測實驗件頂部的溫度，以免實驗件和電加熱器超溫?zé)龤А?/p>

實驗過程中通過操作控制系統(tǒng)調(diào)節(jié)熱水箱內(nèi)的水溫和壓力，使水箱內(nèi)去離子水處于飽和態(tài)；通過控制水冷機組將實驗過程中多余的熱量排入環(huán)境。數(shù)據(jù)采集系統(tǒng)主要包括熱電偶、壓力變送器以及高速攝影儀、功率變送器等數(shù)據(jù)采集和監(jiān)控設(shè)備。所有測點的溫度、壓力、電壓和電流信號都經(jīng)過數(shù)據(jù)采集系統(tǒng)輸入計算機實時顯示并存儲。在給定傾角下，逐步提升加熱功率，使實驗件在設(shè)定功率下的換熱達(dá)到穩(wěn)態(tài)工況，從而獲得熱通量和壁面過熱度隨時間的變化曲線。當(dāng)發(fā)生沸騰危機時，壁面溫度飛升，為保證整個實驗件的安全，切斷加熱器的電源，記錄相關(guān)數(shù)據(jù)。

電加熱器尺寸、熱電偶插孔孔徑以及插孔距換熱表面的距離使用游標(biāo)卡尺測量，不確定度為±0.02 mm，熱電偶測溫不確定度為±0.5℃，實驗過程中的誤差主要來自于晶閘管整流器調(diào)節(jié)功率所造成的波動和熱電偶測溫誤差。在接近臨界熱通量時，熱通量的增加幅度控制在3%以內(nèi)，沸騰危機發(fā)生的判據(jù)為溫度出現(xiàn)5～15℃·s?1的急劇躍升，臨界熱通量為溫度急劇躍升時的熱通量與前一個穩(wěn)定狀態(tài)時的熱通量的算術(shù)平均值。換熱表面的熱通量由式(1)計算

表面溫度由式(2)計算

熱通量的不確定度為4.0%。表面溫度w的不確定度為5.0%。

2 實驗表面

實驗表面-溝槽結(jié)構(gòu)表面如圖4所示，其中，溝槽結(jié)構(gòu)表面的溝槽寬度為1 mm，深度為1 mm，間距為2 mm，溝槽結(jié)構(gòu)表面相對于光表面的面積增加率為98%。在實驗之前，提前使用丙酮擦拭換熱表面，去除換熱表面的污漬。

3 實驗結(jié)果

圖5顯示了溝槽結(jié)構(gòu)表面在5個傾角下的熱通量隨壁面過熱度的變化。從圖中可以看出，熱通量隨壁面過熱度增加而增加，CHF隨著傾角的增大而增大。圖6為溝槽結(jié)構(gòu)表面的沸騰傳熱系數(shù)隨壁面過熱度的變化，表明傳熱系數(shù)隨壁面過熱度的增加而增大。圖7(a)為傾角90°、熱通量約130 kW·m-2時換熱表面的沸騰狀況，溝槽結(jié)構(gòu)表面的溝槽使得在較小的過熱度下汽化核心就能被激活。由圖7可以看到氣泡大部分從溝槽內(nèi)產(chǎn)生，并形成穩(wěn)定的汽化核心。隨著熱通量的增加，換熱表面的氣泡逐漸增多，同時氣泡相互聚合形成較大的氣泡，如圖7(b)、(c)所示，大氣泡逐漸覆蓋整個換熱表面，從圖中可以觀察到溝槽結(jié)構(gòu)能夠更好地減少氣泡聚合，有利于打破氣液兩相流的邊界層，利于提高換熱壁面的潤濕性。在熱通量較高時，溝槽結(jié)構(gòu)能改善換熱表面的潤濕性能，從而使沸騰換熱顯著增強，延遲沸騰危機的到來，CHF大幅提高。圖7(d)為傾角90°、換熱表面發(fā)生臨界危機時的沸騰狀況，圖中顯示氣泡在換熱表面形成一層很薄的氣膜覆蓋整個換熱表面，換熱表面和水之間的換熱主要依靠輻射換熱進(jìn)行，氣膜表面脫離的氣泡數(shù)量很少，同時相比于核態(tài)沸騰，膜態(tài)沸騰的氣泡變小，并且伴隨巨大噪聲。

前期實驗測量了朝下光表面的CHF隨傾角的變化規(guī)律[18]，其中，光表面和溝槽結(jié)構(gòu)表面在不同傾角下的CHF及對應(yīng)的壁溫見表1。從表可以看出，相比于光表面，溝槽結(jié)構(gòu)表面CHF的增幅可達(dá)65%～90%，而且，傾角越大，沸騰換熱強化效果越明顯。在傾角0°時測得的CHF為956 kW·m-2，增加65%；在傾角90°時，CHF最高達(dá)到2075 kW·m-2，增加90%。圖8顯示了實驗中的溝槽結(jié)構(gòu)表面、光表面的CHF和Yang等[15]CHF性能的對比結(jié)果，從圖中可以看出溝槽結(jié)構(gòu)表面和涂層表面的CHF性能明顯好于光表面，同時溝槽結(jié)構(gòu)表面的CHF隨著傾角的增大而增大，但涂層表面的CHF隨著傾角的增大，先減小后增大，在小傾角下，溝槽結(jié)構(gòu)表面明顯高于涂層表面，但隨著傾角的增大，兩者的CHF性能逐漸相當(dāng)。實驗表明溝槽結(jié)構(gòu)表面相比于涂層表面CHF性能更為優(yōu)良，同時考慮到結(jié)構(gòu)可靠性和壽命，溝槽結(jié)構(gòu)表面能更為有效可靠地增強IVR安全裕度。

表1 不同傾角下光表面和溝槽結(jié)構(gòu)表面的CHF及對應(yīng)的壁溫（飽和溫度100℃） Table 1 Local CHF limits at various orientation for plain surface and structured surface (Tsat=100℃)

4 沸騰換熱機理分析

實驗觀察記錄的朝下溝槽結(jié)構(gòu)表面高熱通量時在不同傾角下的氣泡運動形態(tài)，如圖9所示，此時的氣泡結(jié)構(gòu)可以分為蒸汽膜結(jié)構(gòu)和波浪蒸汽層結(jié)構(gòu)。如圖9(a)所示，在傾角為5°時，氣泡持續(xù)長大，和相鄰氣泡聚集成蘑菇狀氣泡群并覆蓋換熱表面，由于換熱表面朝下布置造成氣泡難于逃逸并不斷聚集成大氣泡群直到氣泡群到達(dá)換熱表面的邊界，氣泡開始以大氣泡的形態(tài)逃逸換熱表面，冷卻液體迅速潤濕換熱表面，開啟新的一輪氣泡生長周期。圖9(b)為蒸汽膜結(jié)構(gòu)示意圖，換熱表面的邊界總是被冷卻液體反復(fù)潤濕，但冷卻液體難以進(jìn)入到換熱表面的中心區(qū)域，以致中心區(qū)域容易產(chǎn)生干涸點，隨著熱通量的增加，干涸區(qū)域擴大并相互連成片，直至覆蓋整個換熱表面，從而引發(fā)沸騰危機。

如圖9(a)所示，當(dāng)傾角大于30°時，氣泡在浮力作用下容易沿著換熱表面朝上逃逸，由于氣液界面不穩(wěn)定從而導(dǎo)致氣液界面呈波浪形，波動的氣液界面有利于液體潤濕換熱表面，增強沸騰換熱。Mudawar等[19]和Kim等[20]曾報道過類似的氣液界面不穩(wěn)定周期現(xiàn)象，同時他們提出了轉(zhuǎn)折角的概念，認(rèn)為這個傾角在15°附近，本文實驗結(jié)果發(fā)現(xiàn)轉(zhuǎn)折角在30°附近，造成這種差異的原因可能在于實驗件尺寸、實驗條件等因素的影響。圖9(c)顯示了一個傾斜表面上的波浪蒸汽層結(jié)構(gòu)，氣泡的動態(tài)行為顯示潤濕界面形成在氣液界面接近換熱表面時，換熱表面的邊界也同樣被潤濕。冷卻液體在壓力和表面張力的作用下流向干涸表面，溝槽結(jié)構(gòu)的存在大大增強了換熱表面的潤濕性，從而能夠大幅提高沸騰傳熱系數(shù)，延緩沸騰危機的到來。當(dāng)熱通量增加到足夠大后，潤濕界面的潤濕速率低于冷卻液體的蒸干速率，沸騰危機被觸發(fā)。

5 結(jié) 論

相比于光表面，朝下溝槽結(jié)構(gòu)表面能大幅增強沸騰換熱強度，CHF可提高65%～90%，而且，傾角越大，沸騰換熱強化效果越明顯。實驗觀察發(fā)現(xiàn)了朝下表面高熱通量下氣泡運動形態(tài)存在蒸汽膜結(jié)構(gòu)和波浪蒸汽層兩種結(jié)構(gòu)，轉(zhuǎn)折角在30°附近。朝下溝槽結(jié)構(gòu)表面CHF大幅增加的原因在于，溝槽結(jié)構(gòu)在大幅增加換熱面積的同時顯著改善了換熱表面的潤濕性。

符 號 說 明

A——換熱面積，m2 B——溝槽寬度，m g——重力加速度，m·s?2 I——電流，A K——溝槽深度，m L——正方形換熱表面的長度，m Q——電加熱器功率，W q——熱通量，kW·m?2 S——溝槽間距，m T——溫度，℃ V——電壓，V x——熱電偶測點和換熱表面的距離 λ——熱導(dǎo)率，W·m-1·K?1 η——考慮實驗件散熱損失的保溫系數(shù) 下角標(biāo) m——平均 w——壁面 sat——飽和

References

[1] HENRY R E, FAUSKE H K. External cooling of a reactor vessel under severe accident conditions [J]. Nucl. Eng. Des., 1993, 139(1): 31-43.

[2] REMPE J L, SUH K Y, CHEUNG F B,. In-vessel retention-recent efforts and future needs [C]// Proceedings of the Sixth International Topical Meeting on Nuclear Thermal Hydraulics, Operations and Safety. Nara, Japan, 2004, Paper ID N6P045.

[3] REMPE J L, SUH K Y, CHEUNG F B,. Insights from investigations of in-vessel retention for high powered reactors [C]// Proceedings, 11th International Topical Meeting on Nuclear Reactor Thermal-Hydraulics. Avignon, France, 2005, Paper No. 271.

[4] REMPE J L, SUH K Y, CHEUNG F B,. In-vessel retention strategy for high power reactors final report [R]. INEEL/EXT-04- 02561, Idaho National Laboratory Report, 2005.

[5] VISHNEV I P. Effect of orienting the hot surface with respect to the gravitational field on the critical nucleate boiling of a liquid [J]. Journal of Engineering Physics, 1973, 24(1): 43-48.

[6] EL-GENK M S, GUO Z. Transient boiling from inclined and downward-facing surfaces in a saturated pool [J]. Int. J. Ref., 1993, 16(6): 414-422.

[7] BRUSSTAR M J, MERTE J H, KELLER R B,. Effects of heater surface orientation on the critical heat flux (Ⅰ): An experimental evaluation of models for subcooled pool boiling [J]. Int. J. Heat Mass Trans., 1997, 40(17): 4007-4019.

[8] KIM J, SUH K. One-dimensional critical heat flux concerning surface orientation and gap size effects [J]. Nucl. Eng. Des., 2003, 226(3): 277-292.

[9] ROUGE S. SULTAN test facility for large-scale vessel coolability in natural convection at low pressure [J]. Nucl. Eng. Des.,1997, 169(1/2/3): 185-195.

[10] THEOFANOUS T G, SYRI S. The coolability limits of a reactor pressure vessel lower head [J]. Nucl. Eng. Des., 1997, 169(1): 59-76.

[11] THEOFANOUS T G, TU J P, SALMASSI T,. Quantification of limits to coolability in ULPU-2000 configuration: CRSS-02.05.3 [R/OL]. Santa Barbara: Center for Risk Studies and Safety University of California, 2002. http://pbadupws.nrc.gov/docs/ML0216/ ML021620559. pdf.

[12] CHU T Y, BAINBRIDGE B L, SIMPSON R B,. Ex-vessel boiling experiments: laboratory- and reactor-scale testing of the flooded cavity concept for in-vessel core retention (Ⅰ): Observation of quenching of downward-facing surfaces [J]. Nucl. Eng. Des.,1997, 169(1/2/3): 77-88.

[13] CHU T Y, BENTZ J H, SLEZAK S E,. Ex-vessel boiling experiments: laboratory-and reactor-scale testing of the flooded cavity concept for in-vessel core retention (Ⅱ): Reactor-scale boiling experiments of the flooded cavity concept for in-vessel core retention [J]. Nucl. Eng. Des., 1997, 169(1/2/3): 89-99.

[14] HADDAD K H. An experimental and theoretical study of two-phase boundary layer ?ow on the outside of curved downward-facing surfaces [D]. PA: Pennsylvania State University, 1996.

[15] YANG J, DIZON M B, CHEUNG F B,. CHF enhancement by vessel coating for external reactor vessel cooling [J]. Nucl. Eng. Des., 2006, 236(10): 1089-1098.

[16] NOH S W, SUH K Y. Critical heat flux in various inclined rectangular straight surface channels [J]. Exp. Therm. Fluid Sci., 2014, 52: 1-11.

[17] NOH S W, SUH K Y. Critical heat flux for APR1400 lower head vessel during a severe accident [J]. Nucl. Eng. Des., 2013, 258: 116-129.

[18] ZHONG D W, MENG J A, LI Z X,. Critical heat flux for downward-facing saturated pool boiling on pin fin surfaces [J]. Int. J. Heat Mass Trans., 2015, 87: 201-211.

[19] HOWARD A H, MUDAWAR I. Orientation effects on pool boiling critical heat flux (CHF) and modeling of CHF for near-vertical surfaces [J]. Int. J. Heat Mass Trans., 1999, 42(9): 1665-1688.

[20] KIM Y H, KIM S J, KIM J J,. Visualization of boiling phenomena in inclined rectangular gap [J]. Int. J. Mul. Flow, 2005, 31(5): 618-642.

Saturated pool boiling from downward facing structured surfaces with grooves

ZHONG Dawen, MENG Ji’an, LI Zhixin

(Institutes of Engineering Thermal Physics, School of Aerospace Engineering, Tsinghua University, Beijing 100084, China)

In-vessel retention is a key severe accident management strategy now widely adopted by some nuclear power plants. The saturated pool boiling heat transfer coefficients and critical heat flux (CHF) were measured from downward facing structured surface with grooves in deionized water to enhance the CHF. The orientations were 5°, 30°, 45°, 60° and 90° (vertical). The results showed that the nucleate boiling heat transfer coefficients and the local CHF for structured surface with grooves were consistently higher than those for plain surface. Compared with plain surface, the CHF increase could reach 65%—90% for structured surface. The vapor blanket and wavy vapor layer bubble structures with different CHF trigger mechanisms on the downward facing surface at high heat fluxes were observed on structured surface. The enhancement of the local CHF and the nucleate boiling heat transfer coefficients was mainly due to the significantly increase of heat transfer area and the grooves that effectively improve surface wettability.

heat transfer; surface; phase change; bubble; critical heat flux; orientation

supported by the National Science and Technology Major Project of China (2011ZX06004-008) and the National Basic Research Program of China (2013CB228301).

date: 2016-01-06.

MENG Ji’an, mja@tsinghua.edu.cn

TQ 026.4

A

0438—1157（2016）09—3559—07

10.11949/j.issn.0438-1157.20160015

大型先進(jìn)壓水堆核電站重大專項子課題項目（2011ZX06004-008）；國家重點基礎(chǔ)研究發(fā)展計劃項目（2013CB228301）。

2016-01-06收到初稿，2016-05-31收到修改稿。

聯(lián)系人：孟繼安。第一作者：鐘達(dá)文（1989—），男，博士研究生。