壓水堆核電廠穩(wěn)壓器波動管熱分層現(xiàn)象國內(nèi)外研究概述

陶舒暢 葉竹 秦婧 蘇桐

摘 要

在壓水堆核電廠，波動管中由于溫差會出現(xiàn)熱分層現(xiàn)象。熱分層會對波動管產(chǎn)生不利影響，例如疲勞等，特別是對焊縫的安全構(gòu)成嚴(yán)重威脅。對于熱分層問題，國內(nèi)外開展了廣泛的研究，研究方法通常可分成數(shù)值方法和試驗方法。國外起步較早，取得了很多重要成果;國內(nèi)起步較晚，但國家核安全局對該問題已經(jīng)非常重視，推動了國內(nèi)熱分層問題的研究，相關(guān)研究正在逐步開展，取得了一系列進(jìn)展。本文論述了波動管熱分層產(chǎn)生機(jī)理，總結(jié)了國內(nèi)外該領(lǐng)域的研究進(jìn)展。

關(guān)鍵詞

壓水堆核電廠;波動管;熱分層;研究概述

中圖分類號： TL334? ? ? ?文獻(xiàn)標(biāo)識碼： A

DOI：10.19694/j.cnki.issn2095-2457.2020.19.005

0 前言

在壓水堆核電廠中，反應(yīng)堆冷卻劑系統(tǒng)（RCS）通常有多個環(huán)路，其中一條環(huán)路通過波動管與穩(wěn)壓器相連，從而維持RCS系統(tǒng)壓力穩(wěn)定。波動管是穩(wěn)壓器與RCS系統(tǒng)連接的關(guān)鍵設(shè)備，是一回路壓力邊界的重要組成之一，屬于核安全一級設(shè)備，波動管的安全關(guān)系著整個一回路系統(tǒng)的安全。穩(wěn)壓器與主管道通常存在溫差，熱分層會引起管壁溫差，溫差的出現(xiàn)進(jìn)而引起波動管截面產(chǎn)生彎曲及局部熱應(yīng)力。核電廠在運行過程中，冷卻劑流量不可能完全穩(wěn)定，會存在一定范圍的波動，因此波動管中的熱分層也并不是固定不變，可能發(fā)生位置的偏移或交替的消失再出現(xiàn)等，從而造成波動管中頻繁承受冷熱應(yīng)力沖擊，故熱分層對波動管的安全性構(gòu)成嚴(yán)重威脅。

本文簡單介紹了波動管熱分層的基本原理，總結(jié)了與波動管熱分層相關(guān)的政策，最后對波動管熱分層國內(nèi)外相關(guān)研究進(jìn)展進(jìn)行總結(jié)。

1 熱分層產(chǎn)生機(jī)理及相關(guān)政策

1.1 熱分層產(chǎn)生機(jī)理

反應(yīng)堆冷卻劑系統(tǒng)將堆芯裂變放出的熱能帶出反應(yīng)堆并傳遞給二回路工質(zhì)以產(chǎn)生蒸汽?，F(xiàn)代商用壓水堆核電廠反應(yīng)堆冷卻劑系統(tǒng)一般有二至四條并聯(lián)在反應(yīng)堆壓力容器上的封閉環(huán)路。每一條環(huán)路由一臺蒸汽發(fā)生器、一臺反應(yīng)堆冷卻劑泵及相應(yīng)的管道組成，在其中一條環(huán)路的熱管段上，通過波動管與一臺穩(wěn)壓器相連。穩(wěn)壓器的基本功能是建立并維持反應(yīng)堆冷卻劑系統(tǒng)的壓力，避免冷卻劑在反應(yīng)堆內(nèi)發(fā)生容積沸騰。在電廠穩(wěn)態(tài)運行時，穩(wěn)壓器將反應(yīng)堆冷卻劑系統(tǒng)壓力維持在恒定壓力下;在反應(yīng)堆冷卻劑系統(tǒng)瞬態(tài)運行時，將壓力變化限定在允許值以內(nèi);在事故時，防止反應(yīng)堆冷卻劑系統(tǒng)超壓，維護(hù)反應(yīng)堆冷卻劑系統(tǒng)的完整性。此外，穩(wěn)壓器作為反應(yīng)堆冷卻劑系統(tǒng)的緩沖容器，吸收反應(yīng)堆冷卻劑系統(tǒng)水容積的迅速變化。

穩(wěn)壓器波動管一端接反應(yīng)堆冷卻劑系統(tǒng)一條環(huán)路的熱段，另一端接穩(wěn)壓器下封頭，在高溫、高壓和強(qiáng)放射性條件下輸送反應(yīng)堆冷卻劑，并與其他設(shè)備和管道一起構(gòu)成壓力邊界，是反應(yīng)堆冷卻劑系統(tǒng)中最重要的管道之一。波動管共有5條焊縫，分別為A、B、C、D、E，圖2所示為三環(huán)路壓水堆核電廠的穩(wěn)壓器波動管。該波動管布置采用從一條環(huán)路主管道熱段引出，波動管接管嘴與水平面存在17′的夾角，波動管以17′的坡度連續(xù)向上傾斜，最后豎直向上與穩(wěn)壓器底部相接。波動管焊縫A到焊縫E之間的近似水平的管段被稱為“準(zhǔn)水平（quasi horizontal）段”，熱分層現(xiàn)象就在這段管段中出現(xiàn)。熱分層示意圖如下圖3所示。

波動管內(nèi)冷卻劑溫度分層引起管壁溫度差，在波動管截面產(chǎn)生總體彎曲熱應(yīng)力和局部熱應(yīng)力，并產(chǎn)生波動管非預(yù)期大的位移和支撐載荷。核電廠在運行過程中，冷卻劑流量不可能完全穩(wěn)定，會存在一定范圍的波動，因此波動管中的熱分層也并不是固定不變，可能發(fā)生位置的偏移或交替的消失再出現(xiàn)等，從而造成波動管中頻繁承受冷熱應(yīng)力沖擊，故熱分層對波動管的安全性構(gòu)成嚴(yán)重威脅。

1.2 熱分層相關(guān)政策

1979年，NRC發(fā)布了79-13公告。在這個公告中敘述了電廠在無負(fù)荷或低功率運行模式下，蒸汽發(fā)生器給水管中出現(xiàn)緩慢流動，存在熱分層流型。這種流型的兩種不同溫度層之間缺少充分混合，在水平（或準(zhǔn)水平）管道的橫截面上產(chǎn)生不均勻的溫度分布。這種溫度分布的重復(fù)出現(xiàn)增加了蒸汽發(fā)生器給水管或接管嘴處出現(xiàn)疲勞失效（貫穿管壁裂紋）的可能性。

1988年6月22日，NRC發(fā)表了題為“反應(yīng)堆冷卻劑系統(tǒng)相連管道的熱應(yīng)力”的88-08公告。該公告描述了1987年由西屋公司設(shè)計的Alabama電力公司擁有的壓水堆Farlay 2號機(jī)組在安全殼內(nèi)出現(xiàn)了比平時高的濕度和放射劑量。最后查明是由于ECCS（應(yīng)急堆芯冷卻系統(tǒng)）管系中冷水通過球閥泄漏引起熱分層，由此引起熱疲勞并出現(xiàn)了環(huán)向穿壁裂紋，冷卻劑通過環(huán)向穿壁裂紋泄漏到安全殼內(nèi)，從而使安全殼內(nèi)的濕度和放射性比平時高。公告要求所有核電廠評估與反應(yīng)堆冷卻劑系統(tǒng)相連接的非隔離管道溫度承載情況，確保這些管道不會受到不可接受的熱應(yīng)力

1982年，美國Trojan核電廠開始監(jiān)測穩(wěn)壓器波動管的運動特性，在換料期間發(fā)現(xiàn)波動管出現(xiàn)意外的運動位移，波動管甩擊限制器兩個間隙閉合，然后根據(jù)不同的假設(shè)條件所做出的分析結(jié)果反復(fù)調(diào)整墊塊和間隙的尺寸，問題仍然得不到妥善解決。直到1987年，各種分析證實在原始設(shè)計中沒有考慮的熱分層現(xiàn)象是導(dǎo)致波動管意外運動的根源。此問題引起了NRC的高度重視。1988年10月7日，NRC發(fā)布了88-80信息通告，公布了Trojan核電廠由于熱分層引起的波動管意外運動的研究過程。1988年12月20日，NRC發(fā)布88-11公告“穩(wěn)壓器波動管熱分層”，這份報告要求業(yè)主考慮在熱分層的影響下，按照美國最新ASME規(guī)范，重新對穩(wěn)壓器波動管的完整性進(jìn)行評定。公告的發(fā)表促使了穩(wěn)壓器波動管熱分層研究的廣泛開展。

波動管熱分層現(xiàn)象成為國際核能界重點研究的熱門課題，并引起了我國核安全當(dāng)局的高度重視，在嶺澳二期核電工程的初步安全分析報告中，提出了非常關(guān)注的問題，并要求在壓力邊界的評定中，充分考慮波動管熱分層現(xiàn)象的影響。

2 國內(nèi)外相關(guān)研究概述

2.1 國外研究概述[1-3]

國外對波動管熱分層現(xiàn)象的研究起步較早，關(guān)注度也較高。早在1989年，V.N.Shah和D.A.Conley（美國愛達(dá)荷國家實驗室）就提出對壓水堆核電廠穩(wěn)壓器波動管和噴淋管線由熱分層引起的疲勞損傷進(jìn)行研究。對波動管熱分層現(xiàn)象研究比較深入的國家是美國和法國，最近幾年韓國和日本也做了相應(yīng)的研究;研究方法主要有兩種：數(shù)值方法和試驗方法。

2.1.1 數(shù)值方法

對核電廠系統(tǒng)中穩(wěn)壓器波動管熱分層現(xiàn)象進(jìn)行數(shù)值方法研究是近年來國際上反應(yīng)堆熱工分析的研究熱點。應(yīng)用CFD方法分析波動管熱分層現(xiàn)象，研究對象多針對波動管內(nèi)部流體而忽略管壁，通過數(shù)值計算，得到波動管內(nèi)部流體的分層流動形式和溫度分布。在這方面，國外很多學(xué)者取得了大量有意義的研究成果。V.K.DHIR（美國，1988）等人以水平管道為研究對象，給出了預(yù)測水平管熱分層流場的一維模型，發(fā)現(xiàn)在水平管道中，考慮冷熱流體之間的熱交換得到的熱分層現(xiàn)象比不考慮熱交換時減弱18%;法瑪通公司（法國，1993）在EXPRESS試驗研究的基礎(chǔ)上利用自己開發(fā)的TRIO程序模擬了二維波動管熱分層流場，該程序利用有限元方法求解控制方程，湍流模型采用k-ε模型，模擬結(jié)果與試驗結(jié)果符合良好;F.C.CHANG（美國，1998）在自己開發(fā)的COMMIX程序中研究了雷諾應(yīng)力模型（RSM）在熱分層流場計算中的應(yīng)用，認(rèn)為在穩(wěn)態(tài)情況下RSM模型比k-ε模型更適用于熱分層流場的計算，但在瞬態(tài)情況下RSM模型在該研究領(lǐng)域的應(yīng)用還需要做進(jìn)一步的研究;Ahn Jang Sun（美國，2007）采用低雷諾數(shù)k-ε模型用數(shù)值方法模擬了二維坐標(biāo)系下的穩(wěn)壓器波動管在反應(yīng)堆啟動時的熱分層流場;Peniguel Christophe（美國，2012）采用k-ε模型數(shù)值模擬了三維坐標(biāo)系下的穩(wěn)壓器波動管熱分層流場，這是國際上首次對對穩(wěn)壓器波動管熱分層全模型實現(xiàn)數(shù)值模擬;Hirofumi Hattori等人（日本，2014）用DNS（直接數(shù)值模擬）方法計算了定常和非定常的熱邊界層流場。DNS的最大好處是無需對湍流流動做任何簡化或近似，理論上可以得到相對準(zhǔn)確的計算結(jié)果。這在熱分層流場數(shù)值計算方法上是一個突破，為波動管熱分層的數(shù)值方法研究提供了一個新的途徑。Ildik′o Boros，Attila Asz′odi（匈牙利，2018）針對VVER-440機(jī)組，利用CFX軟件計算了穩(wěn)壓器波動管和高壓安注管的熱分層流場，VVER- 440機(jī)組的波動管為兩條對稱布置的管道組成，在計算中，作者假設(shè)管內(nèi)流體流動方式為層流，溫差為100℃，波動流速為0.03m/s。計算結(jié)果顯示兩條對稱布置的管道中的熱分層現(xiàn)象并不對稱。

總體來說，對熱分層流動的數(shù)值模擬所采用的湍流模型還是比較單一，一般是k-ε模型，而k-ε模型計算熱浮力項是采用傳統(tǒng)的模化方法，其結(jié)果不是非常精確。因此，對三維坐標(biāo)系中波動管的熱分層流場采用何種湍流模型更加有效，特別是切應(yīng)力輸運模型（SST）和雷諾應(yīng)力模型（RSM）在該領(lǐng)域的應(yīng)用還有很廣的探討空間。

2.1.2 試驗方法

除了數(shù)值方法以外，試驗測定和經(jīng)驗研究也是國外穩(wěn)壓器波動管熱分層現(xiàn)象研究的主要方法。Paul Hirschberg，Arthur F.Deardorff（西班牙），John Carey（英國）給出了多座核電廠的運行數(shù)據(jù)，提供了核電廠中與反應(yīng)堆冷卻劑系統(tǒng)相聯(lián)的管道由于發(fā)生泄漏或破裂而產(chǎn)生熱疲勞損傷的經(jīng)驗 。

美國西屋公司在各電廠業(yè)主的支持下，對美國國內(nèi)大量壓水堆核電廠穩(wěn)壓器波動管實施瞬態(tài)溫度監(jiān)測，獲得了這些電廠在正常工況以及升/降溫過程中，在考慮熱分層情形下的穩(wěn)壓器波動管瞬態(tài)，建立了升/降溫過程的瞬態(tài)數(shù)據(jù)庫，此外，還開發(fā)了熱分層分析方法，并按照NRC 88-11公報的要求對核電廠穩(wěn)壓器波動管的完整性重新進(jìn)行了評定。

法國FRAMETOME和EDF公司自1981年開始對穩(wěn)壓器波動管熱分層現(xiàn)象進(jìn)行了大量的研究工作。在大量的實驗和理論分析的基礎(chǔ)上，針對CPY機(jī)組核電廠形成了一套系統(tǒng)的熱分層分析方法。另外，在N4、EPR等電廠的設(shè)計過程中，進(jìn)行大量的試驗和分析，對波動管進(jìn)行改進(jìn)設(shè)計，減小或消除熱分層現(xiàn)象對波動管的疲勞破壞。

FRAMATOME、EDF和CEA聯(lián)合進(jìn)行了EXPRESS試驗，通過布置在試驗裝置管壁上的熱電偶測量波動管的瞬態(tài)溫度分布。試驗結(jié)果指出，流體的溫度分布主要與無量綱數(shù)Froude數(shù)相關(guān)。Froude數(shù)定義為：

式中Ti表示主管道流體溫度，Tp表示穩(wěn)壓器底部流體溫度，T表示某一截面位置的流體溫度。Fr-τ數(shù)據(jù)庫如表1所示，表中下標(biāo)“-”表示截面下部，下標(biāo)“+”表示截面上部。該數(shù)據(jù)庫用于確定熱分層在瞬態(tài)工況下是否消除，并用于確定部分熱分層。當(dāng)τ=1時，說明波動管內(nèi)流體充分?jǐn)嚮旎蛘呋旌希瑹岱謱蝇F(xiàn)象消失，當(dāng)τ<1時，說明當(dāng)時的瞬態(tài)運行工況產(chǎn)生的波動管流動不充分，不足以使熱分層完全消失。但是截面下部的溫度狀態(tài)也會受到波動流量的影響，其影響程度取決于τ值的大小。以焊縫E為例，由表1可以看到，隨著Fr的逐漸增大，焊縫E上部截面的混合函數(shù)值始終等于1，而下部截面的混合函數(shù)τ<1且隨著Fr而變化，焊縫E處存在熱分層現(xiàn)象;當(dāng)Fr=0.13時，焊縫E處下部截面τ=1，由此可以判定，此時焊縫E下部截面的溫度等于來自穩(wěn)壓器的高溫流體溫度，熱分層現(xiàn)象消失。

上述試驗還表明，當(dāng)波動管內(nèi)流量由主管道流向穩(wěn)壓器時，熱分層現(xiàn)象消失，波動管內(nèi)流體溫度等于主管道流體溫度，截面上部承受冷沖擊，沖擊幅度為穩(wěn)壓器和主管道的溫差。

在上述試驗結(jié)果的基礎(chǔ)上，F(xiàn)ROMATOME公司形成了一種分析CPY機(jī)組波動管熱分層熱工水力狀態(tài)的一般方法，但僅適用于CPY機(jī)組。

2.2 國內(nèi)研究概述[4-6]

我國對波動管熱分層現(xiàn)象研究起步較晚。近年來，國家核安全局對該問題也非常重視。因此，國內(nèi)對波動管熱分層現(xiàn)象的研究正在逐步展開，而且也取得了一定的成果。

上海交通大學(xué)的周美伍和徐濟(jì)鋆于2001年描述了核電廠穩(wěn)壓器波動管熱分層現(xiàn)象以及發(fā)生的原因;清華大學(xué)陳麗星對壓水堆穩(wěn)壓器波動管熱分層現(xiàn)象作了比較多的研究，提出了適合熱分層數(shù)值計算的方程RSM模型，并于2003對壓水堆穩(wěn)壓器波動管熱分層現(xiàn)象做了計算，由于缺乏工程背景，作者對波動管幾何模型做了不切實際的近似（用矩形管代替圓管），所得的計算結(jié)果價值不大，但是該文提出的湍流模型和計算方法值得后續(xù)研究者借鑒;2006年，中國核動力研究設(shè)計院的張毅雄等人從理論上分析了穩(wěn)壓器波動管熱分層發(fā)生的條件，以百萬千瓦級三環(huán)路壓水堆核電廠啟堆為例，建立了熱分層瞬態(tài)，研究了熱分層應(yīng)力計算方法。除此之外，席時桐等人對熱分層數(shù)值模擬的算法做了比較多的研究，提出了很多可行的計算方法。

近年來國內(nèi)也開展了試驗方法的研究。中國核動力研究設(shè)計院空泡物理和自然循環(huán)重點實驗室于2003年對秦山二期2號堆穩(wěn)壓器波動管外壁瞬態(tài)溫度和壓力容器安注管外壁瞬態(tài)溫度進(jìn)行測量，得到溫度測試數(shù)據(jù)，這是國內(nèi)首次針對熱分層現(xiàn)象對核電廠實施溫度測量。測試數(shù)據(jù)表明在運行過程中波動管與壓力容器安注管均存在明顯的熱分層現(xiàn)象，但測試過程中沒有記錄波動流量和穩(wěn)壓器及主管道內(nèi)冷卻劑的溫度，無法得到波動管出現(xiàn)熱分層現(xiàn)象的規(guī)律。圖4所示為升溫過程波動管某截面的溫度分布情況。

2007年，中國核動力研究設(shè)計院和西安交通大學(xué)動力工程多相流國家重點實驗室聯(lián)合進(jìn)行了一項研究主管道波動接管嘴傳熱特性的試驗。該試驗根據(jù)相似理論，采用1：9的比例模型，分別在低溫低壓和高溫高壓條件下研究了不同運行工況下近壁流體溫度分布以及流體和壁面的傳熱性能。試驗結(jié)果表明，支流和主流的流速比是影響結(jié)構(gòu)流體溫度分布和傳熱的決定性因素。

2011年，中國核動力研究設(shè)計院賴建永通過研究，發(fā)現(xiàn)了熱分層現(xiàn)象出現(xiàn)的主要原因，建立了焊縫處熱分層現(xiàn)象的消失與Fr的關(guān)系，為改善熱分層現(xiàn)象對波動管的影響提出了可行的建議;2019年，中國核動力研究設(shè)計院唐鵬等以概率論為基礎(chǔ)，利用ANSYS 程序中的可靠性模塊對波動管模型進(jìn)行隨機(jī)抽樣分析，求出在一定置信度下的可靠度曲線，并對輸出隨機(jī)變量的靈敏度和抽樣過程進(jìn)行了分析，求得對結(jié)果影響最大的因素。結(jié)果表明，計算模型可以有效地反映波動管熱分層的實際情況，為波動管結(jié)構(gòu)可靠性分析提供參考。

綜上所述，數(shù)十年來國內(nèi)外諸多學(xué)者通過理論研究、實地測量、數(shù)值模擬和試驗研究等不同的方法對波動管熱分層現(xiàn)象進(jìn)行分析，對其產(chǎn)生的原因以及對波動管結(jié)構(gòu)完整性的影響有了一定程度的認(rèn)識。在對波動管熱分層現(xiàn)象進(jìn)行數(shù)值分析的研究中，國內(nèi)外學(xué)者在湍流模型的選擇上做了大量的工作，但尚未給出統(tǒng)一的結(jié)論，k-ε模型，k-ω模型，RSM模型等常見的湍流模型均有應(yīng)用實例。此外，在改善波動管熱分層現(xiàn)象的方法上，國際上尚缺少深入的研究，還需要繼續(xù)進(jìn)行研究。

3 結(jié)論

本文先分析了波動管熱分層產(chǎn)生機(jī)理，隨后對波動管熱分層在核電領(lǐng)域的政策沿革進(jìn)行了總結(jié)，最后闡述了國內(nèi)外波動管熱分層研究情況概述，并論述了目前國內(nèi)外在該領(lǐng)域的研究所存在的問題。本文對了解穩(wěn)壓器波動管熱分層問題機(jī)理、選擇研究方向具有參考意義。

參考文獻(xiàn)

[1]Wilhelm Kleinder，Hans-Jürgen Golembiewski.Monitoring for Fatigue Examples for Unexpected Component LoadingFramatome ANP GmbH，F(xiàn)ANP NDM2，freyesleberstr.1，91058 Erlangen， Germany，2001.

[2]US NRC.Pressurizer Surge Line Thermal Stratification.NRC Bulletin 88-11， 1988.

[3]Hirofumi Hattori，Tomoya Houra，Yasutaka NaganoDirect numerical simulation of stable and unstable turbulent thermal boundary layers.Department of Mechanical Engineering，Nago ya Institute of Technology，Gokiso-cho，Showa-ku，Japan. Department of Environmental Technology，Nagoya Institute of Technology，Nagoya， Japan，2014.

[4]賴建永，黃偉.布置方式對波動管熱分層現(xiàn)象的影響分析[J].核動力工程，2011年12月第32卷第6期.

[5]賴建永.壓水堆核電廠穩(wěn)壓器波動管熱分層現(xiàn)象研究[D].中國核動力研究設(shè)計院碩士學(xué)位論文，2009.

[6]唐鵬，等.穩(wěn)壓器波動管熱分層應(yīng)力強(qiáng)度可靠性分析[J].核動力工程，2019年6月第40卷增刊1.