氮氣穩(wěn)壓對主泵運行的影響綜述

摘 "要：隨著小型模塊反應(yīng)堆的發(fā)展，氮氣穩(wěn)壓方式得到廣泛的關(guān)注和重視。通過對氮氣穩(wěn)壓過程中的熱力學(xué)分析，發(fā)現(xiàn)氮氣穩(wěn)壓除具有體積小、結(jié)構(gòu)簡單、響應(yīng)快速等優(yōu)點外，使用氮氣穩(wěn)壓的水循環(huán)在相等時間內(nèi)補水泵啟動補水的頻率更小?；趪鴥?nèi)外常規(guī)泵與核主泵兩相流研究現(xiàn)狀，從氣泡析出、氣泡遷移、氣泡聚集以及對主泵性能影響等方面分析傳統(tǒng)泵兩相流與氮氣穩(wěn)壓下的主泵兩相流的區(qū)別。

關(guān)鍵詞：小型模塊堆；氮氣穩(wěn)壓；主泵；兩相流；熱力學(xué)

中圖分類號：TM623 " " "文獻標志碼：A " " " " "文章編號：2095-2945（2024）20-0096-06

Abstract： With the development of small module reactor， nitrogen pressure stabilization has been paid more and more attention. Through the thermodynamic analysis of the process of nitrogen pressure stabilization， it is found that nitrogen pressure stabilization not only has the advantages of small volume， simple structure and fast response， but also has less frequency of starting water supply pump in the same time. Based on the research status of two-phase flow between conventional pump and nuclear main pump at home and abroad， the difference between conventional pump two-phase flow and main pump two-phase flow under stable nitrogen pressure is analyzed from the aspects of bubble precipitation， bubble migration， bubble accumulation and the influence on the performance of main pump.

Keywords： small module reactor; nitrogen pressure stabilization; main pump; two-phase flow; thermodynamics

近年來，小型模塊化反應(yīng)堆（SMR）得到不斷發(fā)展和研制。在SMR的研究中，對反應(yīng)堆的進一步小型化研究，使SMR可以擁有更廣泛的應(yīng)用場景，是SMR的重要發(fā)展方向之一。小型化的一個研究思路是：使用氮氣穩(wěn)壓裝置替代體積較大的穩(wěn)壓器作為一回路的穩(wěn)壓設(shè)備。對于有氮氣穩(wěn)壓裝置的反應(yīng)堆一回路，當一回路的系統(tǒng)壓力下降時，由于存在氮氣的穩(wěn)定器，其上方的氣體因為壓力減小而發(fā)生膨脹，并逐步將下方的冷卻劑輸送到一回路系統(tǒng)中進行補水。相反地，當一回路的系統(tǒng)壓力上升時，冷卻劑的補充量也會隨之增大，從而使得穩(wěn)壓裝置中的氮氣體積被壓縮，以此達到維持系統(tǒng)壓力穩(wěn)定的效果[1]。

氮氣的化學(xué)性質(zhì)穩(wěn)定且價格便宜，使用其進行壓力調(diào)節(jié)時，不會對一回路系統(tǒng)中的設(shè)備造成腐蝕。然而，引入氮氣穩(wěn)壓也可能導(dǎo)致冷卻劑內(nèi)部溶解某些氮氣，本文將綜合國內(nèi)外相關(guān)研究，探討氮氣穩(wěn)壓如何影響主泵的運行。

1 "氮氣穩(wěn)壓過程的熱力學(xué)分析

對于含有氮氣穩(wěn)壓裝置的高溫?zé)峋W(wǎng)系統(tǒng)，通常包括核反應(yīng)堆、換熱器、主泵和氮氣穩(wěn)壓裝置等設(shè)備，系統(tǒng)圖如圖1所示。對于具有氮氣穩(wěn)壓的高溫?zé)峋W(wǎng)系統(tǒng)中，冷卻劑從核反應(yīng)堆吸收熱量之后流入蒸汽發(fā)生器，然后經(jīng)過熱交換后再通過主泵增壓至核反應(yīng)堆進行下一個循環(huán)。高溫?zé)峋W(wǎng)冷卻劑在循環(huán)過程中，會由于設(shè)備功率變化等原因使系統(tǒng)的壓力變化，并依靠氮氣穩(wěn)壓完成壓力調(diào)節(jié)。

在有氮氣穩(wěn)壓的一回路系統(tǒng)中，主泵的入口與氮氣穩(wěn)壓裝置的補水管連接，當一回路系統(tǒng)壓力降低時，穩(wěn)壓裝置通過氣體膨脹實現(xiàn)對一回路系統(tǒng)的補水，從而確保系統(tǒng)壓力維持穩(wěn)定。氮氣穩(wěn)壓裝置維持系統(tǒng)穩(wěn)壓的過程中，能量轉(zhuǎn)化關(guān)系：當?shù)獨馐艿綁嚎s時，系統(tǒng)需要為它做功，這個能量就儲存在了氮氣之中；反之，當?shù)獨怏w積膨脹時，氣體中的儲能將在體積膨脹過程中轉(zhuǎn)化釋放出來，氣體膨脹過程中可將水從穩(wěn)壓罐中壓入一回路系統(tǒng)中。

當系統(tǒng)運行壓力下降至狀態(tài)點1（p1，V1）最小工作壓力時，如圖2所示，補水管開始向該系統(tǒng)注水，在此過程中，系統(tǒng)狀態(tài)會沿著p-V曲線由狀態(tài)點1變化至狀態(tài)點2（p2，V2），系統(tǒng)壓縮氮氣儲存能量，這部分能量存儲起來將被用于氮氣膨脹和補水，當?shù)竭_狀態(tài)點2（p2，V2）時，這時系統(tǒng)已經(jīng)達到了其工作壓力的最大值，補水支路終止運行。

運行過程中，一回路系統(tǒng)運行壓力和穩(wěn)壓罐氮氣壓力必須保持在一定壓力范圍之內(nèi)，即為圖2中狀態(tài)點1至狀態(tài)點2之間的范圍。狀態(tài)點1受壓轉(zhuǎn)變成狀態(tài)點2的過程中，需壓縮氮氣做功，將能量在穩(wěn)壓罐中儲存起來；狀態(tài)點2膨脹回到狀態(tài)點1時，膨脹過程中儲存在壓縮氮氣中的能量對外釋放，將穩(wěn)壓罐下部的水壓進高溫?zé)峋W(wǎng)系統(tǒng)中，從而實現(xiàn)系統(tǒng)補水的目的。

狀態(tài)點1壓縮轉(zhuǎn)變成狀態(tài)點2的過程是一個等溫過程，在此過程中氮氣的能量變化為

式中：“-”代表了氮氣受壓儲存能力，p指壓力，V指體積，下標1和2對應(yīng)于狀態(tài)點1和狀態(tài)點2。氮氣體積膨脹釋放出來的能力也為這個能量值。穩(wěn)壓罐氮氣在狀態(tài)1和狀態(tài)2之間轉(zhuǎn)變的壓縮和膨脹過程，罐內(nèi)氮氣溫度維持不變，氣體狀態(tài)變化是一個等溫過程，可表示為

式中：ΔVN指氮氣體積變化，“-”表示壓縮，VN為初始體積，即整個過程中穩(wěn)壓罐的最大容量V1，p2為最終壓力。

在高溫系統(tǒng)中，冷卻劑的工作壓力在最大值和最小值之間波動也會導(dǎo)致冷卻劑體積產(chǎn)生一定程度的波動。體積波動可以表示為

式中：“-”表示壓力增加體積減小，KTW為冷卻劑的固定溫度壓縮系數(shù)，VW為總體積，Δp為壓力變化量。

在對循環(huán)系統(tǒng)進行補水時，有氮氣穩(wěn)壓工況下每次補充的水量？駐VW+ΔVN，無氮氣穩(wěn)壓工況下每次補充的水量ΔVW。根據(jù)式（5）和式（6），將系統(tǒng)工作壓力的最大值用Pmax表示，可以得到以下結(jié)論：在相同變化條件下，對于循環(huán)水量為VW的熱網(wǎng)系統(tǒng)，有氮氣穩(wěn)壓和無氮氣穩(wěn)壓工況的補水周期之間有著如下關(guān)聯(lián)

式中：TN為有氮氣穩(wěn)壓的系統(tǒng)補水周期，TW為無氮氣穩(wěn)壓的系統(tǒng)補充周期。根據(jù)公式（7）可以得知，具有氮氣穩(wěn)壓的循環(huán)系統(tǒng)，其補水周期更長，系統(tǒng)運行時補水頻率更低。

2 "氮氣穩(wěn)壓對主泵的影響

在使用氮氣穩(wěn)壓過程中，隨著壓力變化和冷卻劑補充過程，會有一定量的氮氣溶解在冷卻劑中并隨著流體流動進入下游的核主泵，使泵內(nèi)的流動狀態(tài)由純液態(tài)進入氣液兩相流的流動狀態(tài)，從而引發(fā)空化現(xiàn)象?？栈赡芤鸷酥鞅卯a(chǎn)生噪聲和振動從而影響冷卻劑的正常輸送。如果核主泵持續(xù)處于空化故障模式下工作，主泵軸承、水力部件、電機等關(guān)鍵零部件都會增加失效的風(fēng)險。因此，氮氣穩(wěn)壓對主泵影響的研究歸根結(jié)底是主泵的氣液兩相流的研究。

2.1 "泵兩相流技術(shù)研究現(xiàn)狀

國內(nèi)外對泵兩相流技術(shù)的研究始于20世紀80年代初，主要考慮葉輪形狀、安放角、葉頂間隙和泵轉(zhuǎn)速等結(jié)構(gòu)參數(shù)，截面含氣率、兩相速度滑移、氣體的可壓縮性、凝結(jié)效應(yīng)、流體黏性、進口壓力、體積流量和密度變化等介質(zhì)性能，氣泡產(chǎn)生、氣泡的積聚、氣泡的運輸、氣泡的破碎、氣泡的運動軌跡、氣泡的尺寸、氣泡的形狀和氣液兩項分布等發(fā)生機理等因素，采用FLUENT、高速視頻監(jiān)測技術(shù)、高分辨率伽馬射線斷層掃描等技術(shù)，建立多種一維、三維的兩相流分析模型，對泵內(nèi)的兩相流的流動規(guī)律以及對泵的性能的影響開展了詳細的研究，重點對兩相流條件下泵的性能預(yù)測進行的分析。

自1985年開始，F(xiàn)uruya[2]對于雙組份氣體-液體的混合物已經(jīng)展開了深入研究和實驗驗證工作，主要是在考慮泵水力部件的幾何參數(shù)、介質(zhì)含氣率、氣液兩相速度滑移等因素，采用不可壓縮一維兩相流體模型，然后又考慮了氣體的凝結(jié)和可壓特性，眾多學(xué)者在研究工作中參考了該模型。但是這種模型并未考慮流體的黏性導(dǎo)致的水力損失，得到的僅是理論揚程，在實際應(yīng)用中具有一定的局限性。Minemura等[3]提出了一種“泡沫”型態(tài)的三維模型來研究泵內(nèi)的氣液兩相流，并使用有限元作為分析工具去求解液相速度場，通過對氣泡運動微分方程積分求解，得出了流場內(nèi)的含氣量及氣泡運動規(guī)律。但在該模型中，氣泡群的跟蹤計算周期長，適用于含氣量較小的計算模型。Minemura等[4]提出了一種基于固定氣泡的泡狀流動計算模型，在該模型中，將雙相流中積聚和吸附在葉輪邊緣的氣泡作為新的計算邊界進行分析，但分析結(jié)果不夠精確。Akinori等[5]首先通過試驗方法得出，采用開式葉輪、調(diào)整葉輪出口安放角、開設(shè)葉輪回流孔等方式，可提高主泵輸送兩相流的能力，并對上述各種改善方式進行了定性分析。Sato等[6]采用5個不同葉片入口角和出口角的離心泵葉輪進行了試驗研究，通過試驗發(fā)現(xiàn)，在低含氣率、入口和出口角較大的葉輪中，氣泡更容易積聚造成揚程突降；在高含氣率時，入口和出口角度越大揚程越大。Kim等[7]使用離心泵葉輪在不同葉頂間隙條件下開展了單相與氣液兩相流體輸送試驗，采用可視化手段觀察了不同含氣條件下兩相流中的氣泡變化規(guī)律，發(fā)現(xiàn)了開式離心泵葉輪可提高離心泵的兩相流輸送性能。Minemura等[8]在Frurya模型基礎(chǔ)上，增加考慮了氣體可壓縮性和流體黏性，針對葉輪出口位置的流道截面突增、出口角傾斜、兩相均勻混合等因素，采用離心泵模型開展了兩相流分析計算，計算結(jié)果與試驗結(jié)果相一致，該計算方法促進了一維葉輪計算方法的發(fā)展。Poullikkas等[9]利用高速攝影技術(shù)對主泵在失水事故工況中的內(nèi)部流動規(guī)律進行了研究，得出了不同流量時主泵內(nèi)的氣泡分布規(guī)律，促進了主泵基于控制容積法的研究模型的發(fā)展。

國內(nèi)泵設(shè)備氣液兩相流的研究始于1995年，王旱祥等[10]基于流道流體為無黏性定常流和單個氣泡的假設(shè)，限定了流動過程中僅氣泡大小可發(fā)生改變但氣泡形狀不允許發(fā)生改變，在離心泵葉輪中分析了氣液兩相流的流動狀態(tài)，發(fā)現(xiàn)了葉輪流道內(nèi)部和出口的氣液兩相均存在相對滑動和分離的現(xiàn)象，得出了葉輪結(jié)構(gòu)和工況參數(shù)與氣液兩相的相對滑動系數(shù)的關(guān)系式。黃思等[11]基于無黏性流體和氣相為完全氣體的假設(shè)，提出了一種三維計算模型，對一種軸流泵的氣液兩相流輸送進行了計算分析，發(fā)現(xiàn)了在介質(zhì)流動的橫向截面上壓力梯度較大，易產(chǎn)生兩相分離。但此計算模型被認為接近于單相流模型，需進一步考慮氣液兩相流真實環(huán)境條件。Wu等[12]利用k-e紊流模型分析了兩相流在葉輪內(nèi)的流動規(guī)律，同時對汽油兩相流和單相液油的流動規(guī)律進行比較，發(fā)現(xiàn)在兩相流工況時葉輪出口壓力有所降低且壓力分布不均，但葉輪內(nèi)的流體速度、壓力分布在該分析比較中未得出。班耀濤等[13]基于氣泡不變形、氣體可壓縮液體不可壓縮、介質(zhì)無黏性、氣液相間存在滑移、氣液相間無傳質(zhì)和傳熱等假設(shè)的基礎(chǔ)上，提出了可以考慮徑向流道截面變化的軸流泵多相流動的性能預(yù)測模型，進行了油氣多相流軸流泵的性能預(yù)測分析，通過對比分析和試驗結(jié)果，低含氣率相比于高含氣率時分析結(jié)果較準確。這說明在兩相流中介質(zhì)中含氣率高時，介質(zhì)的流動規(guī)律更復(fù)雜，需進一步模擬真實流動條件改善計算模型。黃思等[14]在已有氣泡的泡狀流模型基礎(chǔ)上，增加考慮氣泡在液體中的阻力、流場壓差、液體慣性力等影響氣泡遷移的因素，優(yōu)化了氣泡運動方程求解，提出了一種適用于兩相流葉片泵的數(shù)值分析方法。盧金鈴等[15]采用歐拉模型，對兩相流在離心泵葉輪內(nèi)部的三維黏性流動規(guī)律進行了研究，發(fā)現(xiàn)在葉輪邊緣含氣率較低，輪蓋側(cè)吸力面邊緣處氣泡會通過凝結(jié)而導(dǎo)致含氣率較高，當增大進氣率時，水力部件流道內(nèi)的壓力降低，泵的揚程下降。金玉珍等[16]在幾條假設(shè)的基礎(chǔ)上，基于一維控制模型，應(yīng)用歐拉公式研究了高揚程小流量離心泵的雙相流計算方法，對比了計算分析結(jié)果和試驗驗證結(jié)果。余志毅等[17]基于雙相流模型構(gòu)建了離心泵雙相流三維湍流流動的模型和分析方法。附加質(zhì)量力和阻力作為相間作用力以及含氣率波動導(dǎo)致的附加源項均為該模型的考慮因素，得出了一種關(guān)于雙相流的修正算法，采用改進模型分析了離心泵在不同含氣率時葉輪內(nèi)部的流場，試驗結(jié)果與計算結(jié)果一致性較高，表明了計算模型的準確性。余志毅等[18]基于RANS和已有雙相流計算模型，采用SIMPLEC算法，對離心泵輸送兩相流的湍流流場進行了計算分析，也開展了試驗驗證工作，結(jié)果表明葉輪在徑向尺寸差較小時能有效降低氣堵的產(chǎn)生；葉輪入口背面的低壓區(qū)容易導(dǎo)致氣體積聚，有必要優(yōu)化葉輪進口的結(jié)構(gòu)設(shè)計。黃思等[19]利用FLUENT采用歐拉方法分析了雙相流軸流泵的流場，分析了泵葉輪內(nèi)部流體速度、壓力分布、氣液兩相滑移等內(nèi)容，研究了泵流道內(nèi)部兩相流的流動規(guī)律。結(jié)果表明：介質(zhì)因受葉輪旋轉(zhuǎn)的離心力作用，兩相流中的液相流動區(qū)域主要集中在葉輪外緣，而氣相的流動區(qū)域主要集中在輪轂附近；兩相流在導(dǎo)葉內(nèi)的分離情況有明顯改善。通過對比試驗驗證結(jié)果和計算分析結(jié)果，驗證了計算方法的有效性。謝鵬等[20]提出了低轉(zhuǎn)速離心泵的設(shè)計方法，研制了由誘導(dǎo)輪、離心葉輪、旋渦泵葉輪等組成的組合式離心旋渦泵，并開展了氣液混輸狀態(tài)下的空化性能試驗研究，以清水作為試驗對照，研究了兩相流流量和兩相流介質(zhì)含氣率2個重要因素對空化性能的影響，結(jié)果表明：泵在兩相流輸送時空化性能陡降，隨著含氣量的增加，空化性能下降越明顯。李咪等[21]通過采用FLUENT軟件模擬了泵內(nèi)部的兩相流場，研究了旋渦泵內(nèi)的氣液泡輸送性能。分析結(jié)果初步得出了旋渦泵內(nèi)部的兩相流的流動規(guī)律，主要是泵葉片根部為氣泡主要集中區(qū)，含氣率越高氣泡聚集程度越高。李雪琴等[22]基于分形理論，通過建立隨機分析模型，得出了泵兩相流的生成過程，研究了兩相流流動規(guī)律與隨機變量之間的關(guān)系，隨機變量包括分形面積、級數(shù)、維數(shù)受兩相流氣相和液相分布特征、氣泡尺寸變化范圍、流道截面含氣量等，研究得出：分形和維數(shù)越大則泵內(nèi)兩相流中混有越多的大尺寸氣泡，分形截面含氣率越高則氣泡幾何尺寸越均勻，分形級數(shù)越大則兩相流中氣泡占比越大。張人會等[23]采用數(shù)值分析的方法對泵內(nèi)部的非穩(wěn)態(tài)氣液兩相流進行研究，形成了兩相流分界面處氣液變化規(guī)律，并探究了氣液兩相流分界面與泵外特性的關(guān)系。李紅等[24]用高速攝影系統(tǒng)對自吸泵的水力部件試驗系統(tǒng)進行觀測，開展了氣液兩相流可視化的研究。李貴東等[25]通過數(shù)值模擬的方法，開展了離心泵的內(nèi)部流場分析計算，研究了泵內(nèi)部氣液兩相流的流動規(guī)律和水力載荷的變化規(guī)律。袁建平等[26]使用數(shù)值模擬的方法，開展了離心泵在氣液兩相流工況分析，得到了氣液兩相流工況下離心泵的內(nèi)部流動特性。

2.2 "核主泵兩相流技術(shù)研究現(xiàn)狀

核主泵因重要的功能被國內(nèi)外學(xué)者單獨展開研究，基于水泵的研究成果、核主泵的結(jié)構(gòu)以及運行工況等要求，對主泵內(nèi)部氣液兩相流對主泵的相關(guān)影響進行預(yù)測，主要表現(xiàn)為水力性能的研究和安全分析。

自1986年開始，Narabayshi等[27]開展了穩(wěn)態(tài)和瞬態(tài)條件下的核主泵兩相流研究，利用RRAC計算程序分析了不同轉(zhuǎn)速、流量、入口壓力和含氣率時的氣液兩相流流場，計算結(jié)果與試驗結(jié)果較吻合。Lee等[28]基于邊界層理論分析了主泵氣液兩相流，研究了泵葉輪內(nèi)部的氣液滑移規(guī)律，并使用RELAP5程序計算了兩相流條件下的水力性能，程序計算與試驗結(jié)果一致性較高，進一步證明了該程序的準確性。Lee等[29]在高壓工況下進行了全流量主泵兩相流試驗，研究了泵入口局部含氣率與泵流量之間的關(guān)系，得出了泵進口溫度對主泵性能影響較大的結(jié)論。Poullikkas等[30]將泵內(nèi)兩相流氣液分離、密度、空間尺寸、氣液轉(zhuǎn)化、可壓縮性等聯(lián)系起來，分析了氣液兩相流中氣相存在會造成泵揚程損失，提出了一種主泵失水事故條件下的兩相流水力性能的計算模型。Poullikkas等[9]基于高速攝影技術(shù)對主泵失水工況的內(nèi)部流動規(guī)律進行了研究，研究了在失水事故中氣泡流入葉片的遷移過程，得到了高、中、小流量條件下的氣體分布規(guī)律。付強等[31]應(yīng)用開展了主泵兩相流在瞬態(tài)流動情況下的模擬分析，重點研究了主泵流道中兩相流的流動不穩(wěn)定性規(guī)律。同年，王秀禮[32]采用CFD計算的方法對核主泵在兩相流中的流動特性開展了計算分析，得到了主泵流量和含氣率對器水力性能的影響規(guī)律。Bin等[33]針對ACP100核主泵開展了四象限外特性和氣液兩相降級函數(shù)預(yù)測，同時采用葉輪歐拉揚程的曲線分布來評估葉輪流道內(nèi)流向位置和展向位置的液體流動狀況，并引入單位體積熵生成率的方式對葉輪流道內(nèi)流動損失進行定量分析，最終確認兩相流工況含氣率為70%時揚程的降級程度最大。

3 "結(jié)論

通過上述調(diào)研，傳統(tǒng)主泵兩相流的研究已經(jīng)取得一些成果，但與氮氣穩(wěn)壓下核主泵的兩相流技術(shù)研究相比，兩者還存在一定的區(qū)別，具體表現(xiàn)為氣泡的產(chǎn)生析出機理、氣泡的遷移、氣泡積聚，以及對泵整體性能的影響，包括水力性能、介質(zhì)潤滑軸承性能、電機換熱性能等。

從析出的角度來說，傳統(tǒng)主泵兩相流的研究中主要考慮主泵內(nèi)部發(fā)生汽蝕或在入口處混有氣體進入，即由于泵內(nèi)壓力降低至飽和蒸汽壓致使冷卻劑發(fā)生汽化，或假定主泵入口處混有不溶解氣體，從而使得主泵內(nèi)部發(fā)生兩相流。而氮氣穩(wěn)壓下的兩相流是指氮氣在冷卻劑中，由于溶解度變化而導(dǎo)致的析出，其析出涉及氮氣的飽和度、化學(xué)勢、臨界體積濃度、臨界尺寸等概念。因此，兩者存在本質(zhì)的區(qū)別。

從遷移積聚的角度來說，傳統(tǒng)主泵兩相流的研究中無論是相變的氣泡還是假定混入的氣泡均是大尺寸的氣泡，而氮氣穩(wěn)壓條件下析出的氣泡是納米級的微氣泡，兩者尺度上的不同導(dǎo)致其遷移和積聚的機理也不同。

從對泵的影響角度來說，傳統(tǒng)主泵兩相流的研究中，主要針對軸封泵，重點關(guān)注發(fā)生兩相流后氣泡在葉輪中的分布情況、氣泡在泵腔內(nèi)的積聚情況以及對主泵水力性能的影響，著重預(yù)測不同含氣量對主泵揚程、振動噪聲等性能。而氮氣穩(wěn)壓條件下的主泵為屏蔽泵，泵腔與電機腔相通，泵腔和電機腔內(nèi)部溫差較大，會導(dǎo)致氣泡的析出，此外氮氣穩(wěn)壓下的氣泡為微納氣泡，氣泡的遷移和積聚不僅僅發(fā)生在葉輪當中，因此若發(fā)生兩相流，不僅對水力性能造成影響，而且對介質(zhì)潤滑軸承、電機換熱等性能產(chǎn)生影響。

參考文獻：

[1] 陳智，廖龍濤，張英，等.ACP100S浮動核電站核蒸汽供應(yīng)系統(tǒng)控制系統(tǒng)分析與仿真研究[J].南華大學(xué)學(xué)報（自然科學(xué)版），2018，32（6）：9-14.

[2] FURUYA O. An analytical model for prediction of two-phase（noncondensable） flow pump performance[J]. Journal of Fluids Engincering， 1985，107（3）：139-147.

[3] MINEMURA K， UCHIYAMA T. Three-dimension calculation of air-water two-phase flow in centrifugal pump impeller based on a bubbly flow model[J]. ASME Journal of Fluids Engineering， 1993，115（4）：781-783.

[4] MINEMURA K， UCHIYAMA T. Three-dimension calculation of air-water two-phase flow in centrifugal pump impeller based on a bubbly flow model with fixed cavity[J]. JSME Internal Journal， Series B， 1994，37（4）：726-735.

[5] AKINORI F， SHIN-ICHIROU S， SATO. Experiments on air-water flow pump impeller with rotating- stationary circular cascades and recirculating flow holes[J]. JSME International Journal， Series B， 1996，39（3）：575-582.

[6] SATO S， FURUKAWA A， TAKATSU Y. Air-water two-phase performance of centrifugal pump impellers with various blade angles[J]. JSME International Journal， Series B，1996，39（2）：223.220.

[7] KIM Y， TANAKA K. Influence of tip clearance on pump performance of two phase flow in a screw centrifugal pump， Proceedings of The Third Internal Conference on Pump and Fan[J].Tsinghua University Press， 1998：346-357.

[8] MINEMURA K， UCHIYAMA T. Prediction of air-water two-phase flow performance of a centrifugal pump based on one-dimensional two-fluid model[J]. ASME Journal of Fluids Engineering， 1998，120（3）：327-334.

[9] POULLIKKAS A. Effects of two-phase liquid-gas flow on the performance of nuclear reactor coolmgpumps[J].Progress in Nuclear Energy， 2003，42（1）：3-10.

[10] 王旱祥，萬邦烈.電動潛油離心泵中氣液兩相混合物的流動分析[J].石油大學(xué)學(xué)報然科學(xué)版，1995，19（2）：52-58.

[11] 黃思，李汗強，班耀濤，等.葉片式氣液混輸泵水力設(shè)計的一種選橋方法[J].石油機械，1999（12）：9-11.

[12] WU Y L. Three-Dimension calculation of oil-bubble flows through a centrifugal pump impeller[C]//Proceedings of The Third Internal Conference On Pump and Fan， 'Tsinghua University Prese October.

[13] 班耀濤，趙宏，薛敦松.螺旋軸流式油氣多相泵的性能預(yù)測模型[J].工程熱物理學(xué)報，2000，21（2）：187-190.

[14] 黃思，吳玉林.葉片式泵內(nèi)氣液兩相泡狀流的三維數(shù)值計算[J].水利學(xué)報，2001（6）：57-61.

[15] 盧金鈴，席光，祁火同.離心泵葉輪內(nèi)氣液兩相井維流動數(shù)值研究[J].工程熱物理學(xué)報，2003，24（2）：237-240.

[16] 金玉珍，謝腐，脫旭東.小流量高揚程離心旋涌泵氣液混輸揚程的分析[J].浙江理工大學(xué)學(xué)報，2007，24（4）：420-423：.

[17] 余志毅，曹樹良，王國玉.葉片泵內(nèi)氣液兩相流的三維流動數(shù)值模擬[J].北京理工大學(xué)學(xué)報，2007，27（12）：1057-1060，1064.

[18] 余志毅，曹樹良，王國玉.葉片式混輸泵內(nèi)氣液兩相流的數(shù)值計算[J].工程熱物理學(xué)報，2007，28（1）：46-48.

[19] 黃思，王宏君，鄭茂溪.葉片式混輸泵氣液兩相流及性能的數(shù)值分析[J].華南理工大學(xué)學(xué)報：自然科學(xué)版，2007，35（12）：11-16.

[20] 謝鵬，朱祖超.低比轉(zhuǎn)速旋渦泵的汽液混輸汽蝕試驗分析[J].水利學(xué)報，2009，40（12）：1506-1511.

[21] 李咪，盛英英.汽液泡輸旋渦泵內(nèi)部氣相分布與性能預(yù)測研究[J].機電工程，2009，26（7）：32-33，56.

[22] 李雪琴，二君，趙鵬，等.多相泳輸泵內(nèi)汽液兩相分布的隨機分形模型[J].廣西大學(xué)學(xué)報：自然科學(xué)版，2010，35（5）：756-761.

[23] 張人會，郭廣強，楊軍虎，等.液環(huán)泵內(nèi)部氣液兩相流動及其性能分析[J].農(nóng)業(yè)機械學(xué)報，2014，45（12）：99-103.

[24] 李紅，姜波，陸天橋.泵自吸過程氣液兩相流的可視化試驗[J].農(nóng)業(yè)機械學(xué)報，2015，46（8）：59-65.

[25] 李貴東，王洋，鄭意，等.氣液兩相條件下離心泵內(nèi)部流態(tài)及受力分析[J].排灌機械工程學(xué)報，2016，34（5）：369-374.

[26] 袁建平，張克玉，司喬瑞，等.基于非均相流模型的離心泵氣液兩相流動數(shù)值研究[J].農(nóng)業(yè)機械學(xué)報，2017（1）：89-95.

[27] NARABAYSHI T， ARAI K， KUBOKOYA T， et al. Centrifugal pump behavior in steady and transient two-phase flow[J]. Journal of Nuclear Science Technology， 1986，23（2）：136-150.

[28] LEE S， BANG Y S， KIM H J. Prediction of reactor coolant pump performance under two-phase flow conditions[J]. Nuclear Enginecring and Design， 1994，26（2）：179-189.

[29] LEE S， BANG Y S， KIM H J. Experimental study of two-phase pump performance using a full size nuclear reactor pump[J]. Nuclear Engineering and Design， 1999，193：159-172.

[30] POULLIKKAS A. Two phase flow performance of nuclear reactor cooling pumps[J]. Journal of the Korean Nuclear Society， 2000，36（2）：123-130.

[31] 付強，袁壽其，朱榮生，等.失水事故工況下核主泵氣液兩相瞬態(tài)流動特性[J].華中科技大學(xué)學(xué)報（自然科學(xué)版），2013，41（9）：112-116.

[32] 王秀禮.核主泵內(nèi)多相流動瞬態(tài)水力特性研究[D].鎮(zhèn)江：江蘇大學(xué)，2013.

[33] BIN H， MEIMEI Z， KEXIN P， et al. Study On the Four-Quadrant Homologous Characteristic and Two-Phase Flow Head Degradation of a Reactor Coolant Pump[J]. Journal of Pressure Vessel Technology， 2021，143（4）：041404.