程 坤 陸雅哲 陳宏霞 吳 鵬 初 曉

（核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室，四川 成都 610213）

0 引言

浮動(dòng)核電站是搭載有小型壓水堆的海洋船舶裝備，具有能源輸出高效穩(wěn)定、清潔環(huán)保、機(jī)動(dòng)性好等優(yōu)點(diǎn)，是解決我國(guó)近海島礁與遠(yuǎn)海穩(wěn)定能源供應(yīng)難題的理想途徑[1]。由于船舶載體的內(nèi)部空間有限，現(xiàn)有浮動(dòng)核電站設(shè)計(jì)（如中國(guó)ACPR50S、俄羅斯KLT-40S 等）多采用結(jié)構(gòu)緊湊、換熱高效的螺旋管直流蒸汽發(fā)生器作為反應(yīng)堆一、二回路間的關(guān)鍵換熱設(shè)備[2]，其運(yùn)行穩(wěn)定性將直接影響浮動(dòng)核電站的能源輸出品質(zhì)和反應(yīng)堆安全。

作為直流蒸汽發(fā)生器中實(shí)現(xiàn)能量傳遞的核心部件， 螺旋換熱管內(nèi)存在的汽液兩相流動(dòng)不穩(wěn)定性是直流蒸汽發(fā)生器安全可靠運(yùn)行的重要威脅。 兩相流動(dòng)不穩(wěn)定性會(huì)引起并聯(lián)螺旋管內(nèi)流量、 壓降和空泡份額等熱工參數(shù)的大幅振蕩， 造成螺旋管發(fā)生熱應(yīng)力的交替變化及機(jī)械振動(dòng)， 甚至導(dǎo)致蒸發(fā)器傳熱管破裂（SGTR）等事故的發(fā)生。 國(guó)內(nèi)外研究者針對(duì)螺旋管內(nèi)汽液兩相流動(dòng)不穩(wěn)定性誘發(fā)機(jī)理、影響因素、不穩(wěn)定邊界及預(yù)測(cè)技術(shù)等開展了大量的理論分析、數(shù)值仿真和試驗(yàn)研究。 本文基于公開發(fā)表的文獻(xiàn)，綜述螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定性的國(guó)內(nèi)外研究進(jìn)展和重要成果， 分析過去研究中不足并提出未來研究的相關(guān)建議。

1 靜態(tài)條件下螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定性研究

兩相流動(dòng)不穩(wěn)定性是動(dòng)力裝置熱工水力研究領(lǐng)域長(zhǎng)期關(guān)注的重點(diǎn)問題之一。 Boure 等[3]、Kakac 等[4]、Ruspini 等[5]對(duì)各時(shí)期的重要研究成果進(jìn)行綜述。 可以發(fā)現(xiàn)，過去研究主要關(guān)注靜態(tài)條件下沸騰直管通道內(nèi)存在的兩相流動(dòng)不穩(wěn)定性現(xiàn)象，形成了對(duì)直管通道內(nèi)兩相流動(dòng)不穩(wěn)定性誘發(fā)影響機(jī)理較為清晰的認(rèn)知。與之相比，針對(duì)螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定特性的研究開展相對(duì)較少。

對(duì)于兩相沸騰通道來說，流動(dòng)壓降、沸騰換熱、汽泡動(dòng)力學(xué)行為等熱工水力特性是影響決定其兩相流動(dòng)不穩(wěn)定特性的關(guān)鍵因素。螺旋管內(nèi)流體受到離心力的作用會(huì)產(chǎn)生局部二次流、 汽-液兩相非均勻分布等特殊現(xiàn)象，造成螺旋管的熱工水力特性與直管通道相比存在顯著區(qū)別[6]。 Naphon 等[7]指出二次流現(xiàn)象（見圖1）的存在使得螺旋管內(nèi)的單相流動(dòng)行為相比與直管通道更為復(fù)雜，螺旋管內(nèi)單相摩擦壓降要明顯高于相同流量下的直管通道，且局部二次流的攪混作用能夠強(qiáng)化螺旋管的單相流動(dòng)傳熱能力。 由于汽-液兩相存在較大的密度差異，螺旋管內(nèi)兩相熱工水力特性受離心力的影響更加顯著。 Murai 等[9]和Sun 等[10]分別采用可視化實(shí)驗(yàn)測(cè)量和CFD 模擬的方式研究了離心力對(duì)于螺旋管內(nèi)兩相分布及相界面結(jié)構(gòu)的影響，研究結(jié)果表明受到重力和離心力的作用，液相會(huì)聚集在螺旋外側(cè)和底部， 而汽相則會(huì)聚集在螺旋內(nèi)側(cè)和頂部，造成螺旋管周向壁溫和局部傳熱系數(shù)的分布不均。Santini 等[11]研究了HCSG 傳熱管內(nèi)的單相與兩相壓降特性， 結(jié)果表明兩相摩擦壓降與流量的1.91 次冪成正比關(guān)系，證實(shí)了螺旋管彎曲會(huì)導(dǎo)致強(qiáng)烈的耗散效應(yīng)。 Chen 等[12]實(shí)驗(yàn)發(fā)現(xiàn)與直管存在顯著不同的是，螺旋管中過冷沸騰可以在其壁面平均溫度小于飽和溫度的情況下發(fā)生。 Hardik 等[13]利用紅外熱成像技術(shù)研究了螺旋曲率與汽-液密度比率對(duì)螺旋管兩相沸騰傳熱系數(shù)的影響規(guī)律。

圖1 螺旋管內(nèi)局部二次流示意圖[8]

鑒于螺旋管與直管通道在流動(dòng)壓降、局部換熱系數(shù)分布、 沸騰起始點(diǎn)等熱工水力特性上存在較大差異，對(duì)于螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定特性的認(rèn)知以及直管通道經(jīng)典流動(dòng)不穩(wěn)定性理論對(duì)螺旋管的適用性成為研究者們關(guān)注的焦點(diǎn)。楊瑞昌等[14，15]以氟利昂12 為介質(zhì)實(shí)驗(yàn)研究了3 根并聯(lián)螺旋蒸發(fā)管內(nèi)的兩相流動(dòng)不穩(wěn)定行為， 獲得了并聯(lián)螺旋管內(nèi)密度波型脈動(dòng)（DWO）的脈動(dòng)特征及不穩(wěn)定邊界預(yù)測(cè)的經(jīng)驗(yàn)關(guān)系式。Guo 等[16]實(shí)驗(yàn)發(fā)現(xiàn)了水平螺旋管在中、低壓條件下存在的三種不穩(wěn)定性類型：密度波型脈動(dòng)（DWO）、壓力降型脈動(dòng)（PDO）和熱力脈動(dòng)（THO），并擬合獲得描述DWO 發(fā)生邊界的關(guān)系式。 馮自平等[17]在此基礎(chǔ)上進(jìn)一步研究發(fā)現(xiàn)低壓下DWO 的發(fā)生會(huì)削弱螺旋管傳熱能力，螺旋管內(nèi)脈動(dòng)時(shí)均傳熱系數(shù)的最低值僅為穩(wěn)定工況下的10%左右。 朱宏曄等[18]通過建立與氦氣對(duì)流換熱的并聯(lián)螺旋管數(shù)值模型，研究了高溫氣冷堆螺旋管蒸汽發(fā)生器內(nèi)發(fā)生的Ledinegg 流量漂移現(xiàn)象。

依托IRIS 小型堆國(guó)際研發(fā)項(xiàng)目，Papini 等[19]利用模擬HCSG 熱工水力行為的全尺寸實(shí)驗(yàn)裝置研究了螺旋管特殊形狀對(duì)兩相不穩(wěn)定性發(fā)生邊界的影響，結(jié)果表明并聯(lián)螺旋管內(nèi)DWO 受過冷度的影響規(guī)律、脈動(dòng)周期與輸運(yùn)時(shí)間的關(guān)系和直管通道經(jīng)典理論存在明顯不同，指出螺旋離心力作用是造成這些差異的重要原因。 Papini 等考慮到螺旋管內(nèi)特殊的空泡分布及兩相壓降特性，認(rèn)為采用與螺旋管等長(zhǎng)度、等傾斜角度的直管模型來預(yù)測(cè)螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定邊界的方法存在較大誤差。

2 海洋條件下兩相流動(dòng)不穩(wěn)定性的研究進(jìn)展

與陸基核電廠不同，海上浮動(dòng)核電站受海面風(fēng)浪及自身機(jī)動(dòng)的影響，會(huì)產(chǎn)生諸如起伏、搖擺等船體運(yùn)動(dòng)，造成反應(yīng)堆系統(tǒng)空間位置的改變，并引入瞬變外力場(chǎng)作用。 一般認(rèn)為，浮動(dòng)核電站等船用核動(dòng)力裝置在實(shí)海況運(yùn)行時(shí)存在如圖2 所示的六自由度空間運(yùn)動(dòng)，這些運(yùn)動(dòng)形式對(duì)核動(dòng)力系統(tǒng)主要產(chǎn)生兩方面的影響：一方面導(dǎo)致系統(tǒng)空間位置的改變，另一方面引入由附加加速度所致的瞬變外力場(chǎng)作用。由于系統(tǒng)空間位置的變化也可近似看作力場(chǎng)作用方向的改變，因此瞬變外力場(chǎng)作用是海洋條件影響核動(dòng)力系統(tǒng)熱工水力特性的關(guān)鍵所在，其中以搖擺運(yùn)動(dòng)所致瞬變外力場(chǎng)的影響最為復(fù)雜典型且最受關(guān)注。

結(jié)合馬健等[20]、Yan 等[21]和程坤等[22]針對(duì)各時(shí)期海洋條件熱工水力特性研究的綜述分析可知，搖擺會(huì)造成加熱通道內(nèi)流動(dòng)阻力、沸騰起始點(diǎn)、換熱系數(shù)、兩相流型、汽泡行為等熱工水力特性的改變，而搖擺瞬變外力場(chǎng)所引起的冷卻劑非穩(wěn)態(tài)流動(dòng)波動(dòng)是導(dǎo)致上述熱工水力特性改變的最主要原因。

圖2 海洋條件下船體六自由度運(yùn)動(dòng)示意圖

在搖擺瞬變外力場(chǎng)作用下的兩相流動(dòng)不穩(wěn)定性研究方面，譚思超等[23，24]實(shí)驗(yàn)發(fā)現(xiàn)搖擺會(huì)造成圓管通道內(nèi)兩相流動(dòng)不穩(wěn)定起始點(diǎn)的提前，搖擺條件下特有的波谷型脈動(dòng)會(huì)與密度波脈動(dòng)會(huì)發(fā)生疊加耦合，加劇系統(tǒng)的不穩(wěn)定。 Guo 等[25]研究發(fā)現(xiàn)搖擺條件下并聯(lián)圓管通道內(nèi)存在兩個(gè)分別位于高含氣率區(qū)和低含氣率區(qū)的不穩(wěn)定區(qū)域。 Wang 等[26]研究了搖擺運(yùn)動(dòng)對(duì)于高壓自然循環(huán)工況下棒束通道流動(dòng)不穩(wěn)定性的影響，發(fā)現(xiàn)系統(tǒng)穩(wěn)定性會(huì)隨著搖擺周期和幅值的增加而減小。Yu 等[27]實(shí)驗(yàn)獲得了搖擺條件下窄矩形通道內(nèi)強(qiáng)迫循環(huán)流動(dòng)不穩(wěn)定性的主要類型及演化規(guī)律。Tang 等[28]發(fā)現(xiàn)搖擺運(yùn)動(dòng)對(duì)于并聯(lián)矩形通道強(qiáng)迫循環(huán)流動(dòng)不穩(wěn)定性界限功率的影響較小。

由于瞬變外力場(chǎng)作用下兩相流動(dòng)不穩(wěn)定性所表現(xiàn)出的強(qiáng)非線性特點(diǎn)，研究者們嘗試將非線性分析方法和混沌理論應(yīng)用到該領(lǐng)域的研究中。 Guo 等[29]發(fā)現(xiàn)采用快速傅里葉變換（FFT）方法可有效評(píng)估海洋條件對(duì)并聯(lián)通道不穩(wěn)定性的影響。 Zhang 等[30]運(yùn)用混沌時(shí)間序列分析的方法研究了搖擺運(yùn)動(dòng)系統(tǒng)的流量脈動(dòng)混沌特性與演變機(jī)理。

由以上綜述可以看出，當(dāng)前該領(lǐng)域相關(guān)研究主要聚焦于搖擺條件下堆芯燃料通道（如直管、矩形通道、棒束通道等）內(nèi)的流動(dòng)不穩(wěn)定特性，尚無(wú)關(guān)于并聯(lián)螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定性受瞬變外力場(chǎng)影響機(jī)制研究的公開文獻(xiàn)報(bào)道。

3 海洋條件潛在影響分析

根據(jù)目前對(duì)于搖擺瞬變外力場(chǎng)影響效應(yīng)與靜態(tài)條件下螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定特性的認(rèn)知，可以定性推測(cè)搖擺瞬變外力場(chǎng)對(duì)并聯(lián)螺旋管內(nèi)的兩相流動(dòng)的潛在影響：

（1）并聯(lián)螺旋管內(nèi)兩相流體受瞬變外力場(chǎng)的驅(qū)迫作用會(huì)產(chǎn)生不同程度的非穩(wěn)態(tài)流量脈動(dòng)，造成流動(dòng)阻力、沸騰傳熱、兩相流型轉(zhuǎn)換等熱工水力特性改變，進(jìn)而影響兩相流動(dòng)不穩(wěn)定起始點(diǎn)位置和發(fā)生邊界。

（2）螺旋管內(nèi)不同位置處兩相流體所受搖擺瞬變外力場(chǎng)的切向分力會(huì)周期性增強(qiáng)/削弱螺旋離心力的作用效果，如圖3 所示,，造成管內(nèi)局部二次流方向和強(qiáng)度的瞬態(tài)變化；特別對(duì)于兩相流體來說，由于汽、液相間存在質(zhì)量流速、密度等差異，兩者所受的離心力與搖擺附加外力的作用程度不同， 可能會(huì)造成管內(nèi)復(fù)雜的相場(chǎng)時(shí)空分布特性，引發(fā)局部沸騰加劇等現(xiàn)象發(fā)生，對(duì)螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定誘發(fā)機(jī)理造成潛在影響。

圖3 搖擺條件下螺旋管與直管通道受力分析對(duì)比

（3）并聯(lián)螺旋管中搖擺瞬變外力場(chǎng)驅(qū)迫造成的同相流量脈動(dòng)與兩相熱致流動(dòng)不穩(wěn)定性造成的反相管間脈動(dòng)之間會(huì)發(fā)生非線性的疊加耦合，引發(fā)更為復(fù)雜劇烈的熱工參數(shù)脈動(dòng)行為，進(jìn)一步增加螺旋管結(jié)構(gòu)損傷的潛在風(fēng)險(xiǎn)。

基于以上定性推測(cè)，可梳理形成如圖4 所示瞬變外力場(chǎng)作用下并聯(lián)螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定性的各潛在影響因素間的內(nèi)在關(guān)系。 可以看出，搖擺瞬變外力場(chǎng)和螺旋管的離心力是影響搖擺運(yùn)動(dòng)下并聯(lián)螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定特性的兩個(gè)核心因素，研究明確這兩種力場(chǎng)的耦合作用關(guān)系是揭示搖擺運(yùn)動(dòng)對(duì)螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定性影響機(jī)制、實(shí)現(xiàn)不穩(wěn)定發(fā)生邊界準(zhǔn)確預(yù)測(cè)的關(guān)鍵所在。

圖4 瞬變外力場(chǎng)對(duì)螺旋管兩相流動(dòng)不穩(wěn)定性的潛在影響關(guān)系

4 結(jié)論

本文基于國(guó)內(nèi)外公開文獻(xiàn)，梳理總結(jié)了靜態(tài)條件下螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定性的研究進(jìn)展，定性分析了海洋條件對(duì)于螺旋管內(nèi)兩相流動(dòng)穩(wěn)定性的潛在影響。 對(duì)于研究現(xiàn)狀的總結(jié)和未來研究的建議如下：

（1）螺旋管所特有的離心力作用是影響其流動(dòng)穩(wěn)定性的一個(gè)不可忽略的重要因素?，F(xiàn)有研究形成了對(duì)于靜態(tài)條件下螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定宏觀現(xiàn)象特征和關(guān)鍵熱工參數(shù)影響規(guī)律的基本認(rèn)知，但對(duì)于螺旋離心力影響兩相流動(dòng)不穩(wěn)定性的具體程度和內(nèi)在機(jī)理的認(rèn)識(shí)仍存在較大不足。未來需針對(duì)螺旋離心力對(duì)兩相流動(dòng)不穩(wěn)定性的影響程度及影響機(jī)理開展深入的研究。

（2）海洋條件所引入的外力場(chǎng)效應(yīng)對(duì)于船用螺旋管直流蒸汽發(fā)生器的潛在影響不可忽略，當(dāng)前相關(guān)研究的開展極為有限。 基于定性分析，瞬變外力場(chǎng)和螺旋離心力是影響海洋條件下螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定特性的兩個(gè)重要因素，建議未來可從研究這兩種力場(chǎng)耦合作用關(guān)系的角度出發(fā)，揭示海洋條件對(duì)螺旋管內(nèi)兩相流動(dòng)不穩(wěn)定性影響機(jī)制。