單噴嘴蒸汽射流凝結(jié)引起的壓力振蕩研究

武心壯,邱斌斌,種道彤,嚴(yán)俊杰

(1.國(guó)家核電技術(shù)公司上海核工程研究設(shè)計(jì)院, 200233, 上海;2.西安交通大學(xué)動(dòng)力工程多相流國(guó)家重點(diǎn)實(shí)驗(yàn)室, 710049, 西安)

單噴嘴蒸汽射流凝結(jié)引起的壓力振蕩研究

武心壯1,2,邱斌斌2,種道彤2,嚴(yán)俊杰2

(1.國(guó)家核電技術(shù)公司上海核工程研究設(shè)計(jì)院, 200233, 上海;2.西安交通大學(xué)動(dòng)力工程多相流國(guó)家重點(diǎn)實(shí)驗(yàn)室, 710049, 西安)

為了得到蒸汽在過冷水中浸沒射流凝結(jié)引起的壓力振蕩規(guī)律,針對(duì)不同的過冷水溫度和軸向距離對(duì)壓力振蕩特性進(jìn)行實(shí)驗(yàn)研究。實(shí)驗(yàn)系統(tǒng)主要包括蒸汽發(fā)生器、穩(wěn)壓罐、噴嘴、水箱、儀表、可視化系統(tǒng)和數(shù)據(jù)采集系統(tǒng)。測(cè)量得到了不同軸向位置的壓力振蕩峰值,通過快速傅里葉變換得到了壓力振蕩的主頻,并分析了過冷水溫度和軸向距離對(duì)壓力振蕩特性的影響規(guī)律。結(jié)果表明:壓力振蕩峰值總體上隨溫度升高而增大,但在水溫較高時(shí)出現(xiàn)下降趨勢(shì);主頻隨著水溫的降低而增大,但不隨軸向距離的增大而變化。同時(shí),利用前人提出的公式對(duì)本實(shí)驗(yàn)條件下的主頻進(jìn)行了計(jì)算,計(jì)算值與實(shí)驗(yàn)值具有相同的變化趨勢(shì),且吻合得較好。

壓力振蕩;峰值;主頻

由于蒸汽在過冷水中浸沒射流凝結(jié)具有高效的混合與換熱能力,因此與其相關(guān)的設(shè)備廣泛地用于多種工業(yè)場(chǎng)合,如能源行業(yè)中應(yīng)用的汽水混合加熱器、蒸汽-水噴射器[1-2],特別是核工業(yè)中輕水反應(yīng)堆的卸壓系統(tǒng),例如先進(jìn)反應(yīng)堆AP1000中安全殼內(nèi)換料水箱中的鼓泡器卸壓即采用了蒸汽射流凝結(jié)換熱的機(jī)理。

目前,國(guó)內(nèi)外學(xué)者對(duì)蒸汽浸沒射流流場(chǎng)參數(shù)的研究主要包括壓力、溫度和速度等,其中壓力振蕩特性是相關(guān)工業(yè)設(shè)備設(shè)計(jì)中關(guān)注的內(nèi)容。美國(guó)學(xué)者Simpson在研究亞音速蒸汽浸沒射流的過程中,觀察到了周期性的界面運(yùn)動(dòng),以及氣泡生成、生長(zhǎng)、分離的過程,并測(cè)得了壓力脈動(dòng)特性[3]。韓國(guó)學(xué)者Youn等研究了低蒸汽質(zhì)量流率下間歇流區(qū)域的壓力脈動(dòng)特性,發(fā)現(xiàn)壓力波的產(chǎn)生頻率在10～50 Hz時(shí)隨著蒸汽質(zhì)量流率的增加而增加,但是幾乎不受過冷水溫度的影響[4]。韓國(guó)學(xué)者Cho等實(shí)驗(yàn)研究了多孔蒸汽射流壓力波動(dòng)特性,結(jié)果表明壓力波動(dòng)頻率隨著蒸汽質(zhì)量流率和過冷度的增加而增大,并分析了與氣泡動(dòng)力學(xué)特性的關(guān)系[5]。此外,美國(guó)學(xué)者Eden等、劉光耀、潘冬冬等以及本文作者都先后對(duì)蒸汽射流流場(chǎng)中的壓力分布進(jìn)行了研究和分析[1-2,6-9]。

由于蒸汽浸沒射流凝結(jié)引起的壓力振蕩是設(shè)備設(shè)計(jì)應(yīng)關(guān)注的問題,但目前國(guó)內(nèi)外對(duì)此相關(guān)的研究卻比較匱乏,故本文針對(duì)不同的過冷水溫度和軸向距離,對(duì)蒸汽浸沒射流凝結(jié)引起的壓力振蕩特性進(jìn)行實(shí)驗(yàn)研究,得到了上述參數(shù)對(duì)壓力振蕩特性的影響規(guī)律,為豐富汽液兩相流理論積累了實(shí)驗(yàn)依據(jù),并為相關(guān)的設(shè)備設(shè)計(jì)提供了一定的參考。

1 實(shí)驗(yàn)系統(tǒng)與方法

(a)實(shí)驗(yàn)系統(tǒng)示意圖

圖1所示為蒸汽浸沒射流凝結(jié)引起的壓力振蕩實(shí)驗(yàn)研究的系統(tǒng)示意圖和測(cè)點(diǎn)布置示意圖。實(shí)驗(yàn)系統(tǒng)主要包括蒸汽發(fā)生器、穩(wěn)壓罐、噴嘴、水箱、三維支架、壓力傳感器、熱電偶、可視化系統(tǒng)和數(shù)據(jù)采集系統(tǒng)。本實(shí)驗(yàn)中直接測(cè)量的參數(shù)有飽和蒸汽壓力、水箱中過冷水的溫度,以及不同測(cè)點(diǎn)處的壓力振蕩特性。蒸汽參數(shù)由噴嘴入口的傳感器測(cè)量得到,壓力振蕩特性參數(shù)由裝在三維支架上的高頻動(dòng)態(tài)壓力傳感器測(cè)量,過冷水溫度由裝在水箱中的熱電偶測(cè)量。

(b)測(cè)點(diǎn)布置示意圖

本實(shí)驗(yàn)在軸向位置布置了3個(gè)測(cè)點(diǎn),以研究不同軸向位置的壓力振蕩特性,表1所示為詳細(xì)的測(cè)點(diǎn)位置和汽水參數(shù)。實(shí)驗(yàn)中當(dāng)汽水參數(shù)達(dá)到設(shè)定值后,移動(dòng)三維支架上的壓力探針測(cè)量設(shè)定軸向位置的壓力振蕩特性,本實(shí)驗(yàn)采用NI采集系統(tǒng),采樣頻率為5 kHz;高頻動(dòng)態(tài)壓力傳感器響應(yīng)頻率為4 kHz,量程為-100～100 kPa,滿量程精度為0.25%;K型熱電偶精度為1 ℃。由于實(shí)驗(yàn)研究中存在不可避免的測(cè)量誤差和系統(tǒng)誤差,本文使用Moffat方法[10]和文獻(xiàn)[11]中介紹的對(duì)快速傅里葉變換不確定度的計(jì)算方法來(lái)分別分析壓力振蕩幅值和振蕩主頻的不確定度。

表1 實(shí)驗(yàn)條件

注：R和X分別為徑向距離和軸向距離與噴嘴出口直徑的比值。

對(duì)于直接測(cè)量量的不確定度

(1)

(2)

ub=Δ/31/2

(3)

對(duì)于快速傅里葉變換引入的不確定度

(4)

urel1=(8m+3)1/2σ

(5)

σ2=2-2B/12

(6)

urel2=ε/81/2

(7)

(8)

式中:urel1和urel2分別為疊形運(yùn)算舍入誤差引入的不確定度分量和窗函數(shù)以及柵欄效應(yīng)引入的不確定度分量;σ2為相對(duì)誤差的方差;N為傅里葉變換長(zhǎng)度;m=lbN為級(jí)數(shù);B為計(jì)算機(jī)計(jì)算長(zhǎng)度;ε為窗函數(shù)以及柵欄效應(yīng)引入的最大誤差。根據(jù)實(shí)驗(yàn)中使用的壓力傳感器等的精度以及上述不確定度的計(jì)算方法,可以得到實(shí)驗(yàn)中壓力振蕩幅值和振蕩主頻的最大不確定度分別為12.6%和12.8%。

2 實(shí)驗(yàn)結(jié)果與分析

2.1 幅值

(a)tw=20 ℃

(b)tw=40 ℃

(c)tw=60 ℃

圖2所示為X=0,水溫為20、40和60 ℃時(shí),在時(shí)間0.05 s內(nèi)的壓力振蕩特性。從圖中可以明顯看出壓力振蕩規(guī)律,當(dāng)水溫比較低時(shí),壓力振蕩幅值比較小,隨著水溫的升高,壓力振蕩的幅值逐漸增大。這是由蒸汽凝結(jié)特性所決定的,當(dāng)水溫比較低時(shí),過冷水的冷凝作用強(qiáng),蒸汽迅速被凝結(jié),汽液界面小且穩(wěn)定,振蕩幅值小;當(dāng)水溫升高時(shí),冷凝作用減弱,汽液界面增大并且不穩(wěn)定,振蕩幅值大。為了便于分析,圖3給出了該蒸汽質(zhì)量流率下,不同水溫對(duì)應(yīng)的凝結(jié)形態(tài)。低水溫時(shí),蒸汽射流形成了穩(wěn)定的汽羽,水溫從40 ℃到50 ℃時(shí),汽羽逐漸不穩(wěn)定,特別是在高水溫時(shí),汽羽逐漸發(fā)散,所以蒸汽射流的壓力振蕩特性與凝結(jié)特性和凝結(jié)形態(tài)是一致的。

圖3 不同水溫下的凝結(jié)形態(tài)

為了分析壓力振蕩幅值的規(guī)律,本文定義壓力振蕩峰值pmax和均方根值p′如下

pmax=max(p1…pN)

(9)

(10)

本文對(duì)2 s內(nèi)的壓力振蕩特性進(jìn)行分析,圖4所示為不同軸向距離時(shí),壓力振蕩峰值和均方根值隨水溫的變化規(guī)律。圖4a中,在本實(shí)驗(yàn)選取的近距離測(cè)點(diǎn)內(nèi),壓力振蕩峰值在高水溫時(shí)可達(dá)到20 kPa,當(dāng)測(cè)點(diǎn)較遠(yuǎn)時(shí),壓力振蕩峰值約為5 kPa。圖4b反映了壓力振蕩均方根值的變化趨勢(shì)?？傮w上,壓力隨著溫度的升高而升高,但是當(dāng)水溫達(dá)到60 ℃左右時(shí),壓力開始出現(xiàn)下降的趨勢(shì)。韓國(guó)學(xué)者Cho在介紹前人的研究結(jié)果時(shí),報(bào)告過類似的規(guī)律,壓力峰值出現(xiàn)在水溫60～80 ℃左右[5]。本文上述討論中對(duì)壓力振蕩幅值隨水溫升高而增大進(jìn)行了分析,實(shí)際上,當(dāng)在水溫逐漸升高趨于飽和溫度的過程中,冷凝能力逐漸降低,蒸汽射流凝結(jié)現(xiàn)象越來(lái)越不明顯,從而由凝結(jié)產(chǎn)生的壓力振蕩會(huì)逐漸消失。因此,隨著水溫上升到一定程度,壓力振蕩幅值出現(xiàn)下降趨勢(shì)。

(b)均方根值隨過冷水溫度的變化

2.2 頻率

圖5 不同水溫下壓力振蕩的主頻(X=0)

圖2中描述的是蒸汽射流凝結(jié)引起的時(shí)間域壓力振蕩信號(hào),對(duì)于壓力振蕩頻率可以用快速傅里葉變換(FFT)方法分析,并用主頻表征。傅里葉變換把復(fù)雜的信號(hào)考慮成由一定振幅、相位、頻率的基本正弦(余弦)信號(hào)組合而成,將函數(shù)向一組正交的正弦、余弦函數(shù)展開,傅里葉變換的目的就是找出這些基本正弦(余弦)信號(hào)中振幅較大信號(hào)對(duì)應(yīng)的頻率,從而找出復(fù)雜信號(hào)中的主要振動(dòng)頻率特點(diǎn),快速傅里葉變換是傅里葉變換的快速算法。圖5所示為X=0時(shí),不同水溫下壓力振蕩的頻域特性。從圖中可以看出,主頻在200～600 Hz之間,且主頻隨著水溫的降低而增大,這也可以直觀地從圖2的時(shí)間域壓力振蕩特性中看出。這是由于隨著水溫降低,冷水的過冷度增大,冷凝作用增強(qiáng),蒸汽迅速被凝結(jié),汽液界面處氣泡存在的時(shí)間短暫,所以引起的壓力振蕩頻率大。圖6所示為在水溫為40 ℃時(shí),不同歸一化軸向距離下壓力振蕩的頻域特性。從圖中可以看出,雖然隨著軸向距離的增大,幅值逐漸減小,但是主頻相同。這也說(shuō)明在壓力振蕩的傳播中,能量逐漸減小,但是頻率不變。

圖6 不同歸一化軸向距離下壓力振蕩的主頻(tw=40 ℃)

Fukuda在對(duì)蒸汽凝結(jié)引起的壓力振蕩的頻率研究中得到主頻正比于冷水過冷度Δt,反比于噴嘴直徑de,并提出了計(jì)算主頻的簡(jiǎn)單關(guān)系式[12]

f=60Δt/de

(11)

本文利用上述關(guān)系式計(jì)算得到的頻率與本實(shí)驗(yàn)壓力振蕩FFT分析得到的主頻比較結(jié)果如圖7所示,計(jì)算值與實(shí)驗(yàn)值的變化趨勢(shì)一致,并且吻合得較好,誤差在±8%以內(nèi)。

圖7 壓力振蕩頻率計(jì)算值和實(shí)驗(yàn)值的比較

3 結(jié) 論

本文針對(duì)不同位置和過冷水溫度,對(duì)飽和蒸汽在過冷水中浸沒射流凝結(jié)引起的壓力振蕩特性進(jìn)行了實(shí)驗(yàn)研究,主要結(jié)論如下。

(1)蒸汽射流的壓力振蕩特性與凝結(jié)特性和凝結(jié)形態(tài)是一致的。當(dāng)水溫較低時(shí),壓力振蕩幅值較小,隨著水溫的升高,壓力振蕩的幅值逐漸增大。

(2)壓力振蕩幅值隨著溫度的升高而增大,但是當(dāng)水溫達(dá)到60 ℃左右時(shí),壓力開始出現(xiàn)下降的趨勢(shì)。這是由于在水溫逐漸升高趨于飽和溫度的過程中,蒸汽射流凝結(jié)現(xiàn)象越來(lái)越不明顯,由凝結(jié)產(chǎn)生的壓力振蕩會(huì)逐漸消失,從而隨著水溫上升到一定程度,壓力振蕩幅值出現(xiàn)下降趨勢(shì)。

(3)主頻隨著水溫的降低而增大,這是由于水溫降低時(shí),冷水的過冷度大,蒸汽迅速被凝結(jié),汽液界面處氣泡存在的時(shí)間短暫,所以引起的壓力振蕩頻率大。當(dāng)水溫相同時(shí),主頻不隨軸向距離的增大而變化,這也說(shuō)明在壓力振蕩的傳播中,能量逐漸減小,但是頻率不變。

(4)利用前人給出的關(guān)系式計(jì)算得到的頻率與本文得到的主頻變化趨勢(shì)一致,并且吻合得較好,誤差在±8%以內(nèi)。

[1] 武心壯, 潘如東, 嚴(yán)俊杰. 蒸汽浸沒射流凝結(jié)引起的沖擊特性研究 [J]. 工程熱物理學(xué)報(bào), 2012, 33(5): 805-808. WU Xinzhuang, PAN Rudong, YAN Junjie. Research on impact characteristic of steam jet in water [J]. Journal of Engineering Thermophysics, 2012, 33(5): 805-808.

[2] 武心壯, 李文軍, 嚴(yán)俊杰. 音速蒸汽浸沒射流凝結(jié)軸向全壓分布特性研究 [J]. 核動(dòng)力工程, 2012, 33(6): 76-80. WU Xinzhuang, LI Wenjun, YAN Junjie. Research on axial total pressure distributions of sonic steam jet in subcooled water [J]. Nuclear Power Engineering, 2012, 33(6): 76-80.

[3] SIMPSON M E, CHAN C K. Hydrodynamics of a subsonic vapor jet in subcooled liquid [J]. ASME Journal of Heat Transfer, 1982, 104(2): 271-278.

[4] YOUN D H, KO K B, LEE Y Y, et al. The direct contact condensation of steam in a pool at low mass flux [J]. Journal of Nuclear and Technology, 2003, 40(10): 881-885.

[5] CHO S, CHUN S Y, BAEK W P, et al. Effect of multiple holes on the performance of sparger during direct contact condensation of steam [J]. Experimental Thermal and Fluid Science, 2004, 28(6): 629-638.

[6] EDEN T J, MILLER T F, JACOBS H R. The centerline pressure and cavity shape of horizontal plane choked vapor jets with low condensation potential [J]. ASME Journal of Heat Transfer, 1998, 120(4): 999-1007.

[7] 劉光耀, 嚴(yán)俊杰, 潘冬冬, 等. 超音速蒸汽浸沒射流汽羽形狀及壓力分布的實(shí)驗(yàn)研究 [J]. 工程熱物理學(xué)報(bào), 2010, 31(5): 781-784.

LIU Guangyao, YAN Junjie, PAN Dongdong, et al. Research on the steam plume shape and pressure fields of supersonic steam jet in subcooled water [J]. Journal of Engineering Thermophysics, 2010, 31(5): 781-784.

[8] 潘冬冬, 武心壯, 嚴(yán)俊杰, 等. 超音速蒸汽浸沒射流尾部速度的實(shí)驗(yàn)研究 [J]. 工程熱物理學(xué)報(bào), 2010, 31(8): 1324-1326. PAN Dongdong, WU Xinzhuang, YAN Junjie, et al. Experimental study on velocity distribution for the tail of supersonic steam jet [J]. Journal of Engineering Thermophysics, 2010, 31(8): 1324-1326.

[9] WU X Z, YAN J J, LI W J, et al. Experimental study on a steam-driven turbulent jet in subcooled water [J]. Nuclear Engineering and Design, 2010, 240(10): 3259-3266.

[10]MOFFAT R J. Describing the uncertainties in experimental results [J]. Experimental Thermal and Fluid Science, 1988, 1(1): 3-17.

[11]宋瓊, 朱長(zhǎng)春, 牛寶良. 基2快速傅里葉變換引入的不確定度估算 [J]. 計(jì)量學(xué)報(bào), 2004, 3(25): 281-283. SONG Qiong, ZHU Changchun, NIU Baoliang. The uncertainty evaluation of base 2 fast Fourier transfer [J]. Acta Metrologica Sinica, 2004, 3(25): 281-283.

[12]FUKUDA S. Pressure variations due to vapor condensation in liquid: II Phenomena at larger vapor mass flow flux [J]. Journal of the Atomic Energy Society of Japan, 1982, 24(6): 466-474.

[本刊相關(guān)文獻(xiàn)鏈接]

韓輝,何雅玲,唐桂華.硫酸蒸氣在換熱器表面沉積特性的數(shù)值研究.2013,47(11):7-13.[doi:10.7652/xjtuxb201311 002]

王斯民,文鍵.無(wú)短路區(qū)新型螺旋折流板換熱器換熱性能的實(shí)驗(yàn)研究.2012,46(9):12-15.[doi:10.7652/xjtuxb201209 003]

李瑜,李亮,鐘剛云,等.末級(jí)透平對(duì)低壓缸氣動(dòng)和凝結(jié)特性的影響.2012,46(7):16-20.[doi:10.7652/xjtuxb201207 004]

張僑禹,趙紅利,侯予.間壁式篩網(wǎng)網(wǎng)格換熱器的性能分析.2012,46(1):24-29.[doi:10.7652/xjtuxb201201005]

郭憲民,王善云,汪偉華,等.環(huán)境參數(shù)對(duì)空氣源熱泵蒸發(fā)器表面霜層影響研究.2011,45(3):30-34.[doi:10.7652/xjtuxb201103005]

何緯峰,戴義平,馬慶中,等.環(huán)境風(fēng)影響下直接空冷單元背壓預(yù)測(cè)研究.2011,45(1):15-20.[doi:10.7652/xjtuxb2011 01004]

(編輯 荊樹蓉)

PressureOscillationInducedbySteamJetCondensationinWaterThroughaNozzle

WU Xinzhuang1,2,QIU Binbin2,CHONG Daotong2,YAN Junjie2

(1. State Nuclear Power Corporation Ltd., Shanghai Nuclear Engineering Research and Design Institute, Shanghai 200233, China;2. State Key Laboratory of Multiphase Flow in Power Engineering, Xi’an Jiaotong University, Xi’an 710019, China)

In order to obtain the pressure oscillation characteristics due to steam jet in water, an experiment was carried out to investigate the pressure oscillation characteristics at various axial positions and different water temperatures. The experimental system mainly consists of a steam generator, a surge tank, a nozzle, a water tank, instrumentations, a camera and a data acquisition system. The pressure oscillation amplitudes were measured and the main frequency was gained by the FFT method. The effects of the water temperature and position on the oscillation characteristic were analyzed. The results show that pressure oscillation amplitudes increased with water temperature, but a decreasing trend was found at high water temperatures. The main frequency decreased as water temperature increased, while the axial position had no effect on the main frequency. In addition, an existing correlation was used to predict the main frequency, and a good agreement was achieved between the predictions and the experiment data.

pressure oscillation; peak magnitude; main frequency

10.7652/xjtuxb201401009

2013-03-29。 作者簡(jiǎn)介: 武心壯(1982—),男,博士,高級(jí)工程師。 基金項(xiàng)目: 國(guó)家核電技術(shù)公司員工自主創(chuàng)新項(xiàng)目專項(xiàng)資金資助項(xiàng)目(SNP-KJ-CX-2011-0008);國(guó)家自然科學(xué)基金資助項(xiàng)目(51125027)。

時(shí)間: 2013-10-17 網(wǎng)絡(luò)出版地址: http:∥www.cnki.net/kcms/detail/61.1069.T.20131017.0827.011.html

TK2

:A

:0253-987X(2014)01-0048-05