郭雪晴 孫中寧 王建軍

（哈爾濱工程大學(xué) 核安全與仿真技術(shù)國(guó)防重點(diǎn)學(xué)科實(shí)驗(yàn)室 哈爾濱 150001）

自然循環(huán)系統(tǒng)由于結(jié)構(gòu)簡(jiǎn)單、固有安全性好等特點(diǎn)在先進(jìn)核反應(yīng)堆的主冷卻劑系統(tǒng)和一些安全系統(tǒng)設(shè)計(jì)中得到了應(yīng)用。自然循環(huán)系統(tǒng)根據(jù)與環(huán)境的開(kāi)放程度分為開(kāi)式自然循環(huán)系統(tǒng)和閉式自然循環(huán)系統(tǒng)[1]。目前關(guān)于自然循環(huán)系統(tǒng)的研究都針對(duì)閉式自然循環(huán)系統(tǒng)，而對(duì)開(kāi)式自然循環(huán)系統(tǒng)的研究很少。不過(guò)由于開(kāi)式自然循環(huán)系統(tǒng)具有回路結(jié)構(gòu)簡(jiǎn)化、系統(tǒng)壓力波動(dòng)小等特點(diǎn)，已開(kāi)始應(yīng)用于一些非能動(dòng)安全系統(tǒng)的設(shè)計(jì)[2?3]。但目前對(duì)開(kāi)式自然循環(huán)系統(tǒng)的瞬態(tài)運(yùn)行行為，特別是其啟動(dòng)特性認(rèn)識(shí)不足。另外，由于自然循環(huán)系統(tǒng)的流動(dòng)特性對(duì)其結(jié)構(gòu)形式比較敏感，不能將閉式自然循環(huán)系統(tǒng)的研究經(jīng)驗(yàn)[4?6]直接應(yīng)用到開(kāi)式自然循環(huán)系統(tǒng)中，因此有必要對(duì)開(kāi)式自然循環(huán)系統(tǒng)的瞬態(tài)行為特性進(jìn)行詳細(xì)的研究。

本文主要對(duì)開(kāi)式自然循環(huán)系統(tǒng)的啟動(dòng)特性進(jìn)行實(shí)驗(yàn)研究，分析系統(tǒng)在啟動(dòng)過(guò)程中出現(xiàn)的流動(dòng)現(xiàn)象和發(fā)生機(jī)理，最后將這些流動(dòng)不穩(wěn)定性進(jìn)行分類，并繪制系統(tǒng)啟動(dòng)過(guò)程中流動(dòng)不穩(wěn)定邊界圖。

1 實(shí)驗(yàn)回路

開(kāi)式自然循環(huán)系統(tǒng)實(shí)驗(yàn)流程如圖1所示，P0?P4為壓力測(cè)點(diǎn)；T1?T6為溫度測(cè)點(diǎn)；F為流量測(cè)點(diǎn)，主要由冷凝罐、循環(huán)管路、水箱和測(cè)量?jī)x表組成，其中循環(huán)管路包括加熱段、下降段和上升段，尺寸和布置如表1所示。在上升段中，管段11和12為可視化的有機(jī)玻璃管，其他管段均為不銹鋼管。有機(jī)玻璃管與不銹鋼管之間采用快速接頭連接，并用耐高溫膠進(jìn)行密封。換熱段由三根相同的傳熱管并列組成，內(nèi)外徑分別為34 mm和38 mm。上升段和下降段中除擴(kuò)張段16的直徑為145 mm外，其余管段直徑均為50 mm。水箱氣空間與大氣環(huán)境聯(lián)通，系統(tǒng)始終處于低壓運(yùn)行狀態(tài)。自然循環(huán)回路內(nèi)布置有若干個(gè)壓力和溫度測(cè)點(diǎn)，用于測(cè)量不同位置流體的壓力和溫度，以監(jiān)測(cè)閃蒸、沸騰等現(xiàn)象的發(fā)生。循環(huán)流量采用電磁流量計(jì)測(cè)量，測(cè)點(diǎn)位于下降段7上。

實(shí)驗(yàn)回路采用純蒸汽加熱，由鍋爐提供。冷凝罐內(nèi)的蒸汽通過(guò)在換熱段9管外冷凝將熱量傳遞給自然循環(huán)系統(tǒng)內(nèi)部流體，形成密度差，從而驅(qū)動(dòng)系統(tǒng)產(chǎn)生循環(huán)流動(dòng)，將熱量導(dǎo)出。由于水箱直接與環(huán)境相通，因此當(dāng)系統(tǒng)產(chǎn)生兩相循環(huán)流動(dòng)時(shí)，從上升段 16流出的汽水混合物中的蒸汽部分將直接排放到環(huán)境中，而液體則依靠重力回落至水箱進(jìn)行再循環(huán)。為了使加熱段始終處于蒸汽環(huán)境中，需要在冷凝罐底部設(shè)置排水閥，適時(shí)地排放冷凝罐底部積累的凝結(jié)水。

圖1 開(kāi)式自然循環(huán)系統(tǒng)實(shí)驗(yàn)流程圖Fig.1 Experimental loop of the open natural circulation system.

表1 管段尺寸和傾斜角Table 1 Dimensions and inclined angles of the pipes.

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

2.1 實(shí)驗(yàn)現(xiàn)象

在較低加熱功率條件下（冷凝罐蒸汽壓力pc為0.2 MPa），開(kāi)式自然循環(huán)系統(tǒng)的啟動(dòng)過(guò)程隨入口過(guò)冷度(Dtsub)的減小會(huì)依次經(jīng)歷單相穩(wěn)定、兩相振蕩和兩相穩(wěn)定幾個(gè)不同的流動(dòng)階段。其中在兩相振蕩階段，不同入口過(guò)冷度條件下系統(tǒng)表現(xiàn)出的流動(dòng)模式也各不相同。

2.1.1 單相穩(wěn)定流動(dòng)

在高入口過(guò)冷度條件下，自然循環(huán)系統(tǒng)始終處于單相流動(dòng)狀態(tài)，回路中流體溫度的變化對(duì)驅(qū)動(dòng)力影響很小，系統(tǒng)質(zhì)量流量穩(wěn)定維持在 0.38 kg·s?1，如圖2所示，圖2中ts_ p1為測(cè)點(diǎn)P1處壓力p1對(duì)應(yīng)的飽和溫度，t1為測(cè)點(diǎn)T1測(cè)得的流體溫度。

圖 2 單相流動(dòng)參數(shù)變化(Dtsub=55.1 oC)Fig.2 Parameter evolutions of the single-phase flow(Dtsub=55.1 oC).

2.1.2 兩相振蕩流動(dòng)

(1) 閃蒸誘發(fā)的流量振蕩

當(dāng)加熱段入口過(guò)冷度Dtsub減小到35.5 oC時(shí)，系統(tǒng)開(kāi)始產(chǎn)生周期性流量振蕩，如圖3所示。從圖3中可以看出，測(cè)點(diǎn)T1、T2、T3處流體溫度t1、t2、t3均未達(dá)到當(dāng)?shù)貕毫ο碌娘柡蜏囟龋醇訜岫魏蜕仙?0?14內(nèi)的流體始終處于單相流動(dòng)狀態(tài)。為了確定系統(tǒng)是否發(fā)生兩相流動(dòng)，利用測(cè)點(diǎn) P3、T3和 P0處的測(cè)量值對(duì)上升段 16出口處的質(zhì)量含氣率進(jìn)行估算(如圖3(b))。圖3中ts_ p1、ts_ p2和ts_ p3分別為測(cè)點(diǎn) P1、P2和P3處壓力p1、p2和 p3對(duì)應(yīng)的飽和溫度，t1、t2和t3分別為測(cè)點(diǎn)T1、T2和T3測(cè)得的流體溫度。結(jié)果發(fā)現(xiàn)該處流體出現(xiàn)周期性的汽化，說(shuō)明此時(shí)流體在上升段間歇地發(fā)生了閃蒸現(xiàn)象，并且閃蒸起始點(diǎn)位于測(cè)點(diǎn)T3和上升段出口之間的區(qū)域。

圖 3 閃蒸引發(fā)的流量振蕩(Dtsub=35.5 oC)Fig.3 Flashing-induced flow oscillations (Dtsub=35.5 oC).

在一個(gè)振蕩周期開(kāi)始時(shí)，系統(tǒng)處于單相流動(dòng)狀態(tài)，質(zhì)量流量較小，傳熱管出口處流體溫度會(huì)升高，但仍然低于當(dāng)?shù)貕毫ο碌娘柡蜏囟?。?dāng)這部分熱流體進(jìn)入上升段時(shí)，隨著當(dāng)?shù)仂o壓的逐漸降低會(huì)誘發(fā)閃蒸現(xiàn)象。閃蒸產(chǎn)生的空泡使得上升管和下降管內(nèi)流體之間的密度差變大，系統(tǒng)驅(qū)動(dòng)力增加，質(zhì)量流量隨之增大。但質(zhì)量流量的增加會(huì)造成自然循環(huán)系統(tǒng)帶走的熱量和傳熱管外通過(guò)凝結(jié)換熱供給的熱量之間的不匹配，反而使得加熱段出口處流體溫度降低，即使在進(jìn)入上升段后也無(wú)法達(dá)到當(dāng)?shù)貕毫ο碌娘柡蜏囟龋W蒸現(xiàn)象逐漸消失，又回到單相流動(dòng)狀態(tài)。在這個(gè)過(guò)程中，由于加熱段入口過(guò)冷度仍相對(duì)較大，因此發(fā)生閃蒸現(xiàn)象的程度并不劇烈，流量振蕩幅度也相對(duì)較小。

圖4 沸騰和閃蒸共同引發(fā)流量振蕩(Dtsub=26.1 oC)Fig.4 Flow oscillations induced by flashing and boiling (Dtsub=26.1 oC).

(2) 沸騰和閃蒸共同誘發(fā)的流量振蕩

隨著入口過(guò)冷度的繼續(xù)降低，流體在加熱段內(nèi)即可發(fā)生沸騰，加熱段內(nèi)產(chǎn)生的汽水混合物在流經(jīng)上升段的同時(shí)能夠誘發(fā)大規(guī)模的閃蒸，流動(dòng)進(jìn)入沸騰噴發(fā)伴隨閃蒸引發(fā)周期性流量振蕩的模式，如圖4所示。

在低壓蒸汽加熱條件下，過(guò)冷度較低的流體進(jìn)入加熱段后幾乎不會(huì)產(chǎn)生壁面過(guò)冷沸騰，而是積累一定過(guò)熱度后在加熱段內(nèi)突然發(fā)生沸騰，迅速產(chǎn)生蒸汽并急劇膨脹。蒸汽受到加熱段空間的限制會(huì)形成局部高壓，將加熱段內(nèi)的流體向兩端推擠，使得下降段測(cè)得的流量小幅降低甚至出現(xiàn)負(fù)值（如圖4(a)所示），此時(shí)在水平可視化段11內(nèi)能觀察到大氣彈流過(guò)。從加熱段內(nèi)流出的汽水混合物在經(jīng)過(guò)上升段時(shí)由于重位壓頭的降低會(huì)誘發(fā)大量的閃蒸，汽水混合物產(chǎn)生的同時(shí)也會(huì)降低整個(gè)系統(tǒng)的壓力，從而誘發(fā)更大規(guī)模的閃蒸現(xiàn)象，系統(tǒng)質(zhì)量流量快速上升。加熱段出口處的流體溫度隨質(zhì)量流量的快速增加而急劇下降，沸騰過(guò)程結(jié)束，加熱段和上升段自下而上逐漸被單相液體填充，閃蒸過(guò)程也隨之減弱最后消失，系統(tǒng)逐漸恢復(fù)小流量單相流動(dòng)狀態(tài)。在閃蒸逐漸減弱的過(guò)程中，上升段內(nèi)各點(diǎn)的靜壓均逐漸回升，原本能在較低靜壓下發(fā)生閃蒸的流體此時(shí)變?yōu)閱蜗酄顟B(tài)，溫度會(huì)相應(yīng)升高（如t2）。而當(dāng)這部分流體繼續(xù)向上流動(dòng)時(shí)，重位壓降的減小又使其重新發(fā)生程度較弱的閃蒸現(xiàn)象，因此從圖 4(a)中可以看到系統(tǒng)在經(jīng)過(guò)短時(shí)間的單相流動(dòng)后流量會(huì)出現(xiàn)較小的峰值。

隨著入口過(guò)冷度的進(jìn)一步減小，流體更容易在加熱段內(nèi)發(fā)生沸騰現(xiàn)象，且程度更加劇烈，沸騰噴發(fā)對(duì)系統(tǒng)流動(dòng)的影響程度變大。對(duì)比圖 4(a)和圖 5可知，入口過(guò)冷度的減小使振蕩周期逐漸變短，一個(gè)周期內(nèi)加熱段沸騰所占的時(shí)間份額增大，沸騰強(qiáng)度逐漸加強(qiáng)，沸騰噴發(fā)引起的下降段逆流流動(dòng)加劇，圖4中出現(xiàn)的閃蒸引起的小幅流量振蕩消失。隨著加熱段入口過(guò)冷度的減小，單相孕育時(shí)間變短，當(dāng)入口過(guò)冷度為15.9 oC時(shí)，振蕩周期內(nèi)的單相流動(dòng)過(guò)程消失。

圖5 較低過(guò)冷度時(shí)沸騰和閃蒸共同引發(fā)的流量振蕩Fig.5 Flow oscillations induced by flashing and boiling under lower inlet subcooling conditions.

2.1.3 兩相穩(wěn)定流動(dòng)

由于開(kāi)式自然循環(huán)系統(tǒng)的啟動(dòng)過(guò)程中兩相振蕩時(shí)間較短，水箱橫截面積很大，啟動(dòng)過(guò)程中水箱內(nèi)的水位變化可以忽略，因此當(dāng)水箱溫度達(dá)到當(dāng)?shù)卮髿鈮合碌娘柡蜏囟葧r(shí)，流體的入口過(guò)冷度幾乎保持不變(Dtsub=12.1 oC)，系統(tǒng)進(jìn)入兩相穩(wěn)定流動(dòng)階段，如圖6所示。在兩相穩(wěn)定流動(dòng)過(guò)程中，系統(tǒng)質(zhì)量流量較大(1.36 kg·s?1)，加熱段出口處流體的溫度達(dá)不到當(dāng)?shù)貕毫?duì)應(yīng)的飽和溫度而處于單相狀態(tài)。從圖6中可以看到，上升段中測(cè)點(diǎn)T2處的流體也沒(méi)有達(dá)到當(dāng)?shù)仫柡蜏囟龋鴾y(cè)點(diǎn)T3處流體溫度t3卻超過(guò)了當(dāng)?shù)仫柡蜏囟?，這說(shuō)明閃蒸起始點(diǎn)位于上升段 15內(nèi)，并且發(fā)生閃蒸現(xiàn)象時(shí)流體處于過(guò)熱狀態(tài)。

圖6 兩相穩(wěn)定流動(dòng)Fig.6 Stable two-phase flow.

2.2 流動(dòng)不穩(wěn)定性分類

開(kāi)式自然循環(huán)系統(tǒng)在啟動(dòng)過(guò)程中產(chǎn)生的流量振蕩均是由加熱段沸騰或上升段閃蒸引起的系統(tǒng)流量、空泡份額和壓降之間的多重反饋造成的。在兩相自然循環(huán)系統(tǒng)中，密度波不穩(wěn)定性是最常的流動(dòng)不穩(wěn)定。密度波的基本特征是流動(dòng)振蕩周期(tfo)是流體流經(jīng)上升段所需時(shí)間(t)的 1.5?2 倍[6?7]，其中t可以通過(guò)上升段流道容積與平均體積流量的比獲得。圖7給出了流量振蕩周期tfo與t的關(guān)系，從圖7中看出，不論是閃蒸引起的流量振蕩還是沸騰和閃蒸共同誘發(fā)的流量振蕩，周期均在流體流過(guò)上升段所需時(shí)間的1.5?2倍之間，說(shuō)明這兩種不同的振蕩模式均屬于密度波振蕩。

圖7 振蕩周期與流體流經(jīng)上升段所需時(shí)間的關(guān)系Fig.7 Relationship of flow oscillation period and the time required to pass through the chimney.

2.3 流動(dòng)不穩(wěn)定邊界

圖8 為冷凝罐內(nèi)蒸汽壓力不同時(shí)所對(duì)應(yīng)的開(kāi)式自然循環(huán)系統(tǒng)的啟動(dòng)過(guò)程。當(dāng)加熱蒸汽壓力較低時(shí)(pc=0.2 MPa、pc=0.3 MPa)，系統(tǒng)的啟動(dòng)過(guò)程隨入口過(guò)冷度的減小會(huì)依次經(jīng)歷上述所描述的單相穩(wěn)定流動(dòng)、兩相振蕩和兩相穩(wěn)定流動(dòng)三個(gè)階段。隨加熱蒸汽壓力增大(pc=0.4 MPa)，加熱段內(nèi)的沸騰噴發(fā)作用權(quán)重增大，使得系統(tǒng)即使在低入口過(guò)冷度條件下也無(wú)法形成穩(wěn)定的兩相流動(dòng)，而系統(tǒng)啟動(dòng)過(guò)程中的流動(dòng)模式只有單相穩(wěn)定流動(dòng)和周期性振蕩流動(dòng)兩種。

圖8 不同蒸汽壓力下的啟動(dòng)過(guò)程Fig.8 Start-up processes with different steam pressures.

圖 9給出了不同熱流密度和不同入口過(guò)冷數(shù)(Nsub=(hf?hin)/hfg×(rf?rg)/rg)時(shí)系統(tǒng)的流動(dòng)不穩(wěn)定性邊界圖。從圖9中看出，系統(tǒng)的流動(dòng)不穩(wěn)定性只發(fā)生在一定的入口過(guò)冷度和加熱功率范圍內(nèi)；隨著加熱功率的增大，單相穩(wěn)定流動(dòng)區(qū)域逐漸減小，流動(dòng)不穩(wěn)定區(qū)域會(huì)增大；當(dāng)熱流密度超過(guò)120 kW·m?2·K?1時(shí)，系統(tǒng)很難再出現(xiàn)穩(wěn)定的兩相流動(dòng)。在兩相振蕩區(qū)域內(nèi)，靠近高過(guò)冷度分界線的條件更容易使系統(tǒng)產(chǎn)生閃蒸不穩(wěn)定性，而過(guò)冷度較低時(shí)的流動(dòng)振蕩則是由沸騰噴發(fā)和閃蒸共同引起的，并且沸騰噴發(fā)的作用程度隨過(guò)冷度降低而加劇。

圖9 啟動(dòng)過(guò)程流動(dòng)不穩(wěn)定性邊界Fig.9 Flow instability boundary during startup process.

3 結(jié)語(yǔ)

本文對(duì)開(kāi)式自然循環(huán)系統(tǒng)的啟動(dòng)特性進(jìn)行了實(shí)驗(yàn)研究，并分析了不同加熱功率和入口過(guò)冷度對(duì)系統(tǒng)流動(dòng)模式和流動(dòng)特性的影響，最后給出了系統(tǒng)的不穩(wěn)定流動(dòng)邊界。結(jié)論如下：

(1) 開(kāi)式自然循環(huán)系統(tǒng)在低加熱功率條件下隨入口過(guò)冷度的降低會(huì)依次出現(xiàn)單相穩(wěn)定流動(dòng)、兩相振蕩和兩相穩(wěn)定流動(dòng)。在兩相振蕩流動(dòng)過(guò)程中，高入口過(guò)冷度條件會(huì)使系統(tǒng)產(chǎn)生閃蒸不穩(wěn)定性，而較低入口過(guò)冷度下的流動(dòng)振蕩則是由加熱段沸騰噴發(fā)和上升段閃蒸共同引發(fā)的，并且隨入口過(guò)冷度的降低沸騰噴發(fā)的作用越強(qiáng)烈。

(2) 在兩相振蕩流動(dòng)中，閃蒸誘發(fā)的流量振蕩以及沸騰噴發(fā)和閃蒸共同誘發(fā)的流量振蕩均屬于密度波流動(dòng)不穩(wěn)定性。

(3) 系統(tǒng)在啟動(dòng)過(guò)程中，兩相流動(dòng)不穩(wěn)定區(qū)域隨加熱功率的增加會(huì)擴(kuò)大。當(dāng)加熱功率超過(guò)一定值時(shí)，系統(tǒng)很難再形成穩(wěn)定的兩相流動(dòng)。

1 姚偉. 沸騰兩相自然循環(huán)系統(tǒng)穩(wěn)定性的實(shí)驗(yàn)與理論研究[D]. 上海: 上海交通大學(xué), 2000 YAO Wei. Experimental and theoretical research on the instability mechanism of boiling two-phase natural circulation system[D]. Shanghai: Shanghai Jiao Tong University, 2000

2 Byun C S, Jerng D W, Todreas N E, et al. Conceptural design and analysis of a semi-passive containment cooling system for a large concrete containment[J]. Nuclear Engineering and Design, 2000, 199: 227?242

3 Ick Soo Kyung, Sang Yong Lee. Periodic flow excursion in an open two-phase natural circulation loop[J]. Nuclear Engineering and Design, 1996, 162: 233?244

4 Aritomi M, Chiang J H, Mori M. Fundamental studies on safety-related thermo-hydraulics of natural circulation boiling parallel channel flow systems under start-up conditions (mechanism of geysering in parallel channels)[J]. Nuclear Safety, 1992, 33(2): 170?182

5 Chiang J H, Aritomi M, Mori M. Fundamental studies on thermo-hydraulics during start-up in natural circulation boiling water reactors (natural circulation oscillation induced by hydrostatic head fluctuation)[J]. Journal of Nuclear Science and Technology, 1993, 30(3): 203?211

6 Furuya M, Inada F, van der Hagen T H J J.Flashing-induced density wave oscillation in a natural circulation BWR-mechanism of instability and stability map[J]. Nuclear Engineering and Design, 2005, 235:1557?1569

7 Boure J A, Bergles J E, Tong L S. Review of two-phase flow instability[J]. Nuclear Engineering and Design, 1973,25: 165?192