張恩昊，靖劍平，張春明，楊 軍，*

（1. 華中科技大學(xué)能源與動(dòng)力工程學(xué)院核工程與核技術(shù)系，武漢 430074；2. 生態(tài)環(huán)境部核與輻射安全中心，北京 100082）

目前與核電相關(guān)的各類計(jì)算程序已達(dá)上千個(gè)，包括熱工水力系統(tǒng)程序、堆芯子通道程序、堆芯中子物理程序、安全殼分析程序、燃料性能分析程序等多種類別，其中用于反應(yīng)堆熱工水力分析的計(jì)算程序約占10%［1，2］。在核電廠設(shè)計(jì)和事故研究中，反應(yīng)堆熱工水力系統(tǒng)程序的開發(fā)需要大量的理論模型和工程經(jīng)驗(yàn)，因此，是否具有自主化的熱工水力系統(tǒng)程序在一定程度上代表著國(guó)家在核電領(lǐng)域的技術(shù)水平［3］。在目前的核能數(shù)字化轉(zhuǎn)型潮流中，熱工水力系統(tǒng)程序計(jì)算效率與模擬精度的提高也將為數(shù)字化核電及核電廠數(shù)字孿生的實(shí)現(xiàn)提供重要的內(nèi)核支撐［4］。

在過去的幾十年中，因?yàn)槿狈ψ灾骰姆磻?yīng)堆與核電程序，國(guó)內(nèi)研究人員普遍使用RELAP5、CATHARE等國(guó)外較為成熟的熱工水力系統(tǒng)程序［2］來進(jìn)行相關(guān)的設(shè)計(jì)和分析工作，這影響了我國(guó)的核電“走出去”戰(zhàn)略。2010年前后，眾多單位開始將自主化核電工程軟件的研發(fā)列為重點(diǎn)項(xiàng)目［5］。目前，中核集團(tuán)、中廣核集團(tuán)、國(guó)家電投集團(tuán)、中國(guó)科學(xué)院及相關(guān)高校開發(fā)的第一批反應(yīng)堆熱工水力系統(tǒng)程序已進(jìn)入了初步驗(yàn)證確認(rèn)階段。國(guó)家核安全局頒布了一系列標(biāo)準(zhǔn)與法規(guī)對(duì)核電相關(guān)軟件的驗(yàn)證過程進(jìn)行規(guī)范，如《核電廠設(shè)計(jì)安全規(guī)定》［6］與《核動(dòng)力廠安全評(píng)價(jià)與驗(yàn)證》［7］，這些標(biāo)準(zhǔn)和法規(guī)要求，必須對(duì)核電廠安全分析中使用的計(jì)算機(jī)程序進(jìn)行驗(yàn)證與確認(rèn)，而《核動(dòng)力廠安全分析用計(jì)算機(jī)軟件開發(fā)與應(yīng)用（試行）》［8］則對(duì)安全分析相關(guān)軟件的驗(yàn)證與確認(rèn)做了具體要求。

近年來國(guó)家核安全局已經(jīng)啟動(dòng)了對(duì)國(guó)產(chǎn)自主化軟件的認(rèn)證和評(píng)估工作［9］，截至2021年10月，共有25款程序提交給了國(guó)家核安全局［10］。程序評(píng)估的重點(diǎn)包括功能完整性、相關(guān)關(guān)系式適用性、驗(yàn)證充分性、結(jié)果保守性、輸入敏感性、不確定性、質(zhì)量保證性等。但仍有一些關(guān)鍵問題限制了程序的評(píng)估工作：缺乏可用的實(shí)驗(yàn)數(shù)據(jù)以進(jìn)行程序的驗(yàn)證、不確定性與保守性軟件的評(píng)估需要大量的確認(rèn)工作、程序的認(rèn)證仍需要一段時(shí)間。

一些國(guó)外研究者已經(jīng)對(duì)核電領(lǐng)域的熱工水力系統(tǒng)程序進(jìn)行了系統(tǒng)性的介紹。Alessandro Petruzzi［11］等人對(duì)目前熱工水力系統(tǒng)程序的特點(diǎn)和局限進(jìn)行了介紹，并對(duì)程序的評(píng)估方法與未來發(fā)展進(jìn)行了討論。Glenn A. Roth［12，13］等人對(duì)一些典型熱工水力系統(tǒng)程序的控制方程、流型劃分、數(shù)值方法、閉合關(guān)系式、程序假設(shè)等方面進(jìn)行了比較和討論，并指出了程序未來的改進(jìn)方向。國(guó)內(nèi)的劉志弢［1］等人與靖劍平［2］等人也都對(duì)幾種典型的熱工水力系統(tǒng)程序進(jìn)行了介紹，并分析了我國(guó)熱工水力系統(tǒng)程序的發(fā)展路線。

本文總結(jié)了國(guó)外主流的熱工水力系統(tǒng)程序及其發(fā)展方向，對(duì)我國(guó)目前正在開發(fā)的幾種熱工水力系統(tǒng)程序進(jìn)行了介紹，并對(duì)目前熱工水力系統(tǒng)程序的發(fā)展趨勢(shì)進(jìn)行了討論。具體程序如表1所示。

表1 國(guó)內(nèi)外熱工水力系統(tǒng)程序Table 1 Worldwide thermal hydraulic system analysis codes

1 國(guó)外熱工水力系統(tǒng)程序

國(guó)外的核電技術(shù)及核電軟件的開發(fā)都起步較早，在20世紀(jì)末已經(jīng)有了多個(gè)熱工水力系統(tǒng)程序，如RELAP、CATHARE、ATHLET等。隨著核電技術(shù)與計(jì)算機(jī)技術(shù)的進(jìn)步，這些程序還在不斷完善和擴(kuò)展。

本節(jié)介紹了一些目前仍在使用且不斷更新的國(guó)外典型熱工水力系統(tǒng)程序最新開發(fā)進(jìn)展情況。

1.1 RELAP

RELAP（Reactor Excursion and Leak Analysis Program）程序是由美國(guó)愛達(dá)荷國(guó)家實(shí)驗(yàn)室（INL）開發(fā)的輕水堆瞬態(tài)分析程序，其RELAP5版本也是目前世界上使用最廣泛的熱工水力系統(tǒng)程序。RELAP程序主要有兩個(gè)開發(fā)方向，分別是RELAP5-3D和RELAP-7。

RELAP5-3D［14］程序包含多維流體力學(xué)模型、多維中子動(dòng)力學(xué)模型和BPLU（Border Profiled Lower Upper）矩陣求解器，其主要特點(diǎn)是可以對(duì)反應(yīng)堆組件進(jìn)行三維模擬。此外，程序中包含25種工作流體，可用于第四代堆的安全分析。RELAP5-3D程序的最新版本為2018年6月發(fā)布的RELAP5-3D Version 4.4.2。

RELAP-7是INL于2011年啟動(dòng)開發(fā)的下一代反應(yīng)堆系統(tǒng)安全分析程序。RELAP-7在計(jì)算體系結(jié)構(gòu)、軟件設(shè)計(jì)、數(shù)值方法和物理模型等方面都進(jìn)行了較大的改進(jìn)，預(yù)計(jì)可以將程序的分析能力擴(kuò)展至幾乎所有的反應(yīng)堆系統(tǒng)模擬場(chǎng)景。

相比于RELAP5，RELAP-7的主要改進(jìn)有：（1）建立了七方程兩相流模型（液體、氣體和界面壓力）；（2）數(shù)值近似在空間和時(shí)間上都具有二階精度；（3）對(duì)于長(zhǎng)時(shí)間瞬態(tài)的隱式緊耦合時(shí)間積分；（4）具有與高保真度物理場(chǎng)緊密耦合的能力；（5）易于與正在開發(fā)的多維堆芯模擬器進(jìn)行耦合。

1.2 TRACE

TRACE（TRAC/RELAP Advance Computer Engine）程序是在TRAC-P、TRAC-B、RELAP基礎(chǔ)上由NRC主導(dǎo)開發(fā)的旗艦性熱工水力系統(tǒng)程序［15］。

TRACE程序中使用了兩流體三場(chǎng)模型，來提高對(duì)兩相流動(dòng)的預(yù)測(cè)精度。此外程序可以使用VESSEL組件來進(jìn)行三維幾何模型的建立［16］，該組件主要用于模擬壓力容器，也可以用于其他反應(yīng)堆部件的模擬。

TRACE程序從1997年開始開發(fā)，最新版本為TRACE V5.840［17］。目前TRACE程序的一個(gè)長(zhǎng)期目標(biāo)是將原有的兩流體、兩場(chǎng)模型過渡到兩流體、四場(chǎng)模型。NRC通過程序評(píng)估與維護(hù)計(jì)劃（CAMP）［17］等軟件組織在全世界范圍內(nèi)對(duì)TRACE等程序進(jìn)行驗(yàn)證和確認(rèn)工作。

1.3 CATHARE

CATHARE（Code for Analysis of Thermal-Hydraulics during an Accident of Reactor and Safety Evaluation）程序是由法國(guó)原子能委員會(huì)（CEA）聯(lián)合法國(guó)電力集團(tuán)（EDF）以及法國(guó)核安全防護(hù)研究所（IRSN）共同開發(fā)的熱工水力系統(tǒng)程序［18］。

CATHARE程序使用兩流體六方程模型且組件中包含三維模塊，可以模擬多種熱工水力瞬態(tài)工況［19］。目前程序可以與流體特性數(shù)據(jù)庫(kù)REFPROP結(jié)合，使CATHARE的工作范圍可以擴(kuò)展到許多非核系統(tǒng)。

CATHARE 3的研發(fā)始于2006年。目前的最新版本為2021年12月發(fā)布的CATHARE 3 V2.2.0。CATHARE 3程序目前改進(jìn)的方向包括：更先進(jìn)的兩相流模型，如改為多場(chǎng)模型或增加湍流方程；對(duì)三維模型進(jìn)行改進(jìn)；使用精細(xì)且非均勻的結(jié)構(gòu)網(wǎng)格；與其他尺度的熱工水力系統(tǒng)程序耦合；使程序適用于第四代反應(yīng)堆等。

1.4 ATHLET

ATHLET（Analysis of Thermal-Hydraulic of Leaks and Transients）程序是由德國(guó)核設(shè)施與反應(yīng)堆安全研究中心（GRS）開發(fā)的熱工水力系統(tǒng)程序，主要用于分析核電廠的運(yùn)行工況、異常瞬態(tài)以及各種失水事故［20］。

ATHLET程序提供了五方程模型與兩流體六方程模型來模擬瞬態(tài)運(yùn)行中的流體動(dòng)力學(xué)行為，提供了一個(gè)專門用于穩(wěn)態(tài)計(jì)算的四方程模型，并加入了全范圍漂移通量模型來計(jì)算流體相間的相對(duì)速度。ATHLET程序的工作流體主要為輕水和重水，為方便未來第四代堆的設(shè)計(jì)分析，還可將流體設(shè)置為氦、鈉、液態(tài)鉛、鉛鉍共晶、熔融鹽和用戶自定義的流體等，而程序提供的兩種流體動(dòng)力學(xué)模型也可以支持這些流體的模擬。

ATHLET程序的最新版本為2021年11月發(fā)布的ATHLET 3.3。目前程序的更新方向包括：物理模型的修正、工作流體類型的拓展、對(duì)三維模塊的完善等。

1.5 RETRAN-3D

RETRAN（RELAP4 TRANSIENT）程序是由美國(guó)電力研究院（EPRI）開發(fā)的系統(tǒng)性瞬態(tài)熱工水力系統(tǒng)程序，用于分析輕水反應(yīng)堆的運(yùn)行瞬態(tài)、小破口失水事故等事件［21］。

RETRAN-3D程序由RETRAN-02程序發(fā)展而來。程序中的流動(dòng)模型以一維均勻平衡態(tài)（HEM）模型（三方程）為基礎(chǔ)，增加兩種可以用于提高兩相流模擬精度的可選模型：通過增加滑移流模型選項(xiàng)來模擬兩相流動(dòng)過程中速度的差異（四方程），通過增加滑移流模型和蒸汽連續(xù)性方程來計(jì)算兩相速度的不均勻性和溫度的不平衡性（五方程）。程序中的堆芯物理模型可以選擇一維或三維中子動(dòng)力模型，以及點(diǎn)堆模型。

1.6 COBRA/TRAC

COBRA/TRAC（Coolant Boiling in Rod Arrays Code/Transient Reactor Analysis Code）［22，23］程序是由美國(guó)西北太平洋國(guó)家實(shí)驗(yàn)室（PNNL）開發(fā)的最佳估算系統(tǒng)性熱工水力系統(tǒng)程序，由COBRA-TF與TRAC-PD2組合而成，主要用于分析反應(yīng)堆大破口失水事故。

TRAC-PD2程序由美國(guó)洛斯阿拉莫斯國(guó)家實(shí)驗(yàn)室開發(fā)，用于模擬反應(yīng)堆一次側(cè)系統(tǒng)，求解一維五守恒方程漂移通量模型。

COBRA-TF［24］程序最初在1980年由PNNL在NRC贊助下開發(fā)，也被賓夕法尼亞州立大學(xué)（PSU）的反應(yīng)堆動(dòng)力學(xué)和燃料建模小組（RDFMG）改進(jìn)和更新過，用于對(duì)反應(yīng)堆容器進(jìn)行三維建模，求解兩流體三場(chǎng)的八個(gè)三維守恒方程。在COBRA/TRAC程序中COBRA-TF作為TRAC-PD2的堆芯模塊，兩者互為邊界條件。

COBRA/TRAC于1983年開發(fā)，后由西屋公司發(fā)展為WCOBRA/TRAC，用于AP系列反應(yīng)堆的大破口失水事故最佳估算。WCOBRA/TRAC程序也被用于對(duì)我國(guó)CAP1400堆型的大破口失水事故進(jìn)行評(píng)估［25］。

1.7 MARS

MARS［26］（Multi-dimensional Analysis of Reactor Safety）是由韓國(guó)原子能研究院（KAERI）通過整合與重組RELAP5/MOD 3.2.1.2和COBRATF開發(fā)的多維熱工水力系統(tǒng)程序，用于對(duì)輕水堆瞬態(tài)進(jìn)行多維、多用途的真實(shí)熱工水力系統(tǒng)分析。

MARS程序于1996年開始開發(fā)，目前的最新版本為MARS -KS1.5。MARS程序最初可以對(duì)輕水堆一維和三維熱工水力系統(tǒng)及燃料的瞬態(tài)響應(yīng)進(jìn)行分析。在之后的開發(fā)中，其模擬能力逐漸擴(kuò)展到重水反應(yīng)堆、研究型反應(yīng)堆和許多先進(jìn)反應(yīng)堆類型［27］。

在開發(fā)出MARS后，KAERI開始開發(fā)CUPID程序，用于分析核反應(yīng)堆中組件級(jí)的瞬態(tài)多維兩相流現(xiàn)象，以滿足多維分析的需求［28］。CUPID可以與系統(tǒng)程序耦合，在系統(tǒng)的特定部分進(jìn)行CUPID模擬，對(duì)其他部分進(jìn)行系統(tǒng)程序模擬，從而為系統(tǒng)程序提供局部縮放的能力。

1.8 SPACE

SPACE（Safety and Performance Analysis Code for Nuclear Power Plants）［29-31］程序是由韓國(guó)核工業(yè)的多家單位共同開發(fā)的多維熱工水力系統(tǒng)程序。SPACE程序從2006年開始開發(fā)，參與者包括韓國(guó)電力研究院（KEPRI）、韓國(guó)電力工程公司（KOPEC）、韓國(guó)核燃料有限公司（KNF）。

SPACE程序采用了多維兩流體三場(chǎng)模型，可以模擬壓水堆各種瞬態(tài)與假定事故中遇到的兩相流現(xiàn)象。程序中的網(wǎng)格系統(tǒng)包括結(jié)構(gòu)化/交錯(cuò)網(wǎng)格和非結(jié)構(gòu)化/同位網(wǎng)格，可以對(duì)反應(yīng)堆系統(tǒng)中的復(fù)雜幾何模型進(jìn)行建模。

目前已通過兩相空泡流與兩相環(huán)狀流等驗(yàn)證實(shí)驗(yàn)，驗(yàn)證了SPCAE程序中的多維兩流體三場(chǎng)模型。

2 國(guó)內(nèi)熱工水力系統(tǒng)程序

為突破核電軟件這一核電國(guó)際化發(fā)展的“卡脖子”問題，近年來國(guó)內(nèi)各單位大力推動(dòng)自主化核電軟件的開發(fā)。本節(jié)介紹了國(guó)內(nèi)目前一些自主化反應(yīng)堆熱工水力系統(tǒng)程序的開發(fā)現(xiàn)狀。（如圖1所示）

圖1 國(guó)內(nèi)熱工水力系統(tǒng)程序總體發(fā)展Fig.1 Development of domestic system thermalhydraulic code

2.1 LOCUST

LOCUST程序?yàn)橹袊?guó)廣核集團(tuán)有限公司（CGN）自主研發(fā)的兩相流熱工水力系統(tǒng)程序，主要用于壓水堆失水事故分析計(jì)算，目前版本為1.2版［32］。LOCUST程序使用的兩相流模型為兩流體六方程模型，可以模擬反應(yīng)堆運(yùn)行及事故過程中的熱工水力現(xiàn)象。同時(shí)為滿足保守分析的需求，程序設(shè)置有滿足美國(guó)聯(lián)邦法規(guī)10 CFR 50附錄K要求的保守模型。

目前中廣核已經(jīng)開展了11個(gè)典型驗(yàn)證算例的分析計(jì)算，結(jié)果表明程序能正確地反映物理現(xiàn)象的響應(yīng)過程，并針對(duì)華龍一號(hào)開展了分離效應(yīng)臺(tái)架和整體效應(yīng)臺(tái)架試驗(yàn)確認(rèn)，結(jié)果表明軟件可以合理地預(yù)測(cè)有關(guān)事故的進(jìn)程及關(guān)鍵現(xiàn)象。

2.2 ARSAC

ARSAC（Advanced Reactor System Analysis Code）為中國(guó)核動(dòng)力研究設(shè)計(jì)院（NPIC）開發(fā)的自主化熱工水力系統(tǒng)程序［33］，主要用于對(duì)壓水堆熱工水力響應(yīng)特性進(jìn)行分析。該程序從2010年開始規(guī)劃，2014年正式開展程序研發(fā)，目前已完成程序的編制、測(cè)試和部分驗(yàn)證工作。

ARSAC程序的流動(dòng)模型以氣液兩相非均勻流和非平衡態(tài)流體動(dòng)力學(xué)模型為基礎(chǔ)，可以對(duì)核電廠的熱工水力瞬態(tài)和事故過程進(jìn)行模擬分析。程序的主要特點(diǎn)為較先進(jìn)的矩陣求解算法和物性分析模塊、較精細(xì)化的壁面換熱模型和再淹沒分析模型。ARSAC程序的驗(yàn)證需要獲取大量可靠的實(shí)驗(yàn)數(shù)據(jù)［34］。對(duì)于大破口失水事故過程，動(dòng)力院計(jì)劃采用LOFT、UPTF、PKL、LOBI等臺(tái)架的相關(guān)工況開展驗(yàn)證工作；對(duì)于小破口失水事故，計(jì)劃采用SEMISCALE、BETHSY、PKL、ATLAS等臺(tái)架的相關(guān)工況開展驗(yàn)證。除此之外，已經(jīng)使用ARSAC程序?qū)θA龍一號(hào)大破口失水事故進(jìn)行了模擬，驗(yàn)證了該程序?qū)θ穗娬镜乃矐B(tài)熱工水力分析及事故分析的能力。

2.3 TRANTH

TRANTH（TRansient ANalysis code of Thermal Hydraulics）是由NPIC自主研發(fā)的具有完全自主知識(shí)產(chǎn)權(quán)的壓水堆瞬態(tài)分析軟件［35］，是由中國(guó)核工業(yè)集團(tuán)有限公司（CNNC）開發(fā)的應(yīng)用于“華龍一號(hào)”等三代核電機(jī)組的核電設(shè)計(jì)與分析軟件包NESTOR的一部分［36］。

TRANTH程序主要適用于二代和三代先進(jìn)壓水堆設(shè)計(jì)與安全分析。程序可模擬多環(huán)路壓水堆冷卻劑系統(tǒng)邊界完整情況下的所有瞬態(tài)，即非失水事故工況和蒸汽發(fā)生器傳熱管破裂的情況。

TRANTH已完成V2.0版本的開發(fā)，目前的軟件研發(fā)工作主要圍繞華龍系列堆型的非能動(dòng)余熱排出功能開展，根據(jù)非能動(dòng)系統(tǒng)的要求為程序增加相應(yīng)的模塊。

2.4 cosFlow

2010年，國(guó)家電力投資集團(tuán)中央研究院（SPICRI）開始聯(lián)合國(guó)內(nèi)10余家單位，研發(fā)具有完全自主知識(shí)產(chǎn)權(quán)的核電廠安全分析與工程設(shè)計(jì)一體化軟件包COSINE（COre and System INtegrated Engine for design and analysis）［37］。

cosFlow是COSINE軟件包中的熱工水力與安全分析程序［38，39］，包含子通道分析程序cosSUBC［40］、反應(yīng)堆系統(tǒng)分析程序cosSYST［41］、安全殼程序cosCONT［42］等。程序采用一維瞬態(tài)兩流體模型，計(jì)算范圍覆蓋從正常運(yùn)行到設(shè)計(jì)基準(zhǔn)事故的各類工況。

cosFlow程序目前已大量用于全范圍模擬機(jī)的熱工水力建模工作。除此之外，cosFlow還可以與嚴(yán)重事故堆芯機(jī)理分析程序cosSA耦合來對(duì)嚴(yán)重事故過程進(jìn)行計(jì)算。

2.5 GINKGO

GINKGO［43-45］程序是由CGN自主開發(fā)的熱工水力系統(tǒng)程序，主要用于研究壓水堆核電廠運(yùn)行或非失水事故瞬態(tài)過程中一二回路系統(tǒng)響應(yīng)。

GINKGO程序使用了點(diǎn)堆模型以及兩相流模型、控制系統(tǒng)模型等模擬壓水堆瞬態(tài)過程，可以對(duì)反應(yīng)堆系統(tǒng)熱工水力相應(yīng)、燃料-冷卻劑傳熱、蒸汽發(fā)生器內(nèi)相變過程等進(jìn)行模擬。程序可用于壓水堆核電廠的系統(tǒng)設(shè)計(jì)、安全分析與設(shè)備校驗(yàn)。

2.6 其他

除上述幾個(gè)針對(duì)傳統(tǒng)輕水反應(yīng)堆開發(fā)的大型熱工水力系統(tǒng)程序外，國(guó)內(nèi)也開發(fā)了一些針對(duì)某一堆型或某類堆型的微機(jī)型熱工水力程序。（如表2所示）

表2 國(guó)內(nèi)其他熱工水力程序開發(fā)情況Table 2 Development of other domestic thermal hydraulic analysis code

西安交通大學(xué)的秋穗正等人［46］基于模塊化建模思想開發(fā)了鈉冷快堆瞬態(tài)熱工水力及安全分析軟件THACS（Transient Thermal-Hydraulic Analysis Code for Sodiumcooled fast Reactors），程序中建立了鈉冷快堆系統(tǒng)關(guān)鍵部件的熱工水力模型和輔助模型，并通過了國(guó)際基準(zhǔn)題EBR-Ⅱ SHRT-17實(shí)驗(yàn)的初步驗(yàn)證。魏詩(shī)穎等人［47］建立了用于鉛鉍冷卻沸水快堆（PBWFR）主回路系統(tǒng)熱工水力分析的數(shù)學(xué)物理模型，開發(fā)了適用于PBWFR的熱工水力系統(tǒng)程序SACOL。王偉偉等人［48］針對(duì)AP1000的具體結(jié)構(gòu)和運(yùn)行特點(diǎn)，使用FORTRAN語(yǔ)言開發(fā)了AP1000非失水事故瞬態(tài)熱工水力程序RETAC（REactor Transient Analysis Code）。

NPIC的劉偉［49］等人針對(duì)板型燃料反應(yīng)堆的堆芯結(jié)構(gòu)及系統(tǒng)特點(diǎn)，開發(fā)了板型燃料反應(yīng)堆系統(tǒng)瞬態(tài)熱工水力程序SYSTRAN。

西北核技術(shù)研究所的田曉艷［50］等人建立了熱管冷卻雙模式空間堆（HP-BSNR）的熱工水力模型，采用FORTRAN語(yǔ)言開發(fā)了用于HPBSNR堆芯概念設(shè)計(jì)和安全分析的穩(wěn)態(tài)程序STHA_HPBSNR。

3 熱工水力系統(tǒng)程序發(fā)展趨勢(shì)

目前的熱工水力系統(tǒng)程序仍存在一些局限性［51］，限制了系統(tǒng)程序計(jì)算精度與范圍的提升，主要的限制可以分為三類：（1）部分物理模型存在精度不足或適用范圍小等缺陷，程序無法適用于新的反應(yīng)堆型或事故工況；（2）目前程序已經(jīng)達(dá)到了兩流體六方程模型的內(nèi)在局限性；（3）由于計(jì)算機(jī)性能限制，反應(yīng)堆回路和組件在建模時(shí)必須被簡(jiǎn)化，降低了模擬的準(zhǔn)確性。

程序物理模型中缺陷存在的原因主要有兩種：一是程序開發(fā)時(shí)僅重點(diǎn)關(guān)注了部分較為敏感的現(xiàn)象，對(duì)于重要性較低的現(xiàn)象模擬精度較差；二是部分復(fù)雜物理現(xiàn)象的模型建立主要依靠實(shí)驗(yàn)數(shù)據(jù)，對(duì)于超出實(shí)驗(yàn)數(shù)據(jù)范圍的工況計(jì)算結(jié)果的可靠性較差。對(duì)于這些問題，需要建立實(shí)驗(yàn)臺(tái)架收集大量實(shí)驗(yàn)數(shù)據(jù)來驗(yàn)證并完善系統(tǒng)程序中的物理模型，如我國(guó)的ACME［52］臺(tái)架、韓國(guó)的ATLAS［53］臺(tái)架、德國(guó)的PKL［54］臺(tái)架等。此外，許多國(guó)外系統(tǒng)程序中還加入了三維模型，以更方便地研究壓力容器、蒸汽發(fā)生器等部件的復(fù)雜熱工水力現(xiàn)象，而國(guó)內(nèi)熱工水力系統(tǒng)程序的開發(fā)則對(duì)三維模型關(guān)注較少。

對(duì)于兩流體六方程模型的限制，需要對(duì)現(xiàn)有的兩相流方程進(jìn)行改進(jìn)。目前已經(jīng)有一些系統(tǒng)程序在這方面做出了嘗試，如RELAP-7程序中使用了七方程兩相流模型；普渡大學(xué)的Ishii［55］等人建立了界面面積輸運(yùn)方程來為兩流體模型提供界面面積濃度，該方法已經(jīng)被用于TRACE的兩流體三場(chǎng)模型中［56］；SPCAE程序中也使用了兩流體三場(chǎng)模型；此外TRACE程序還正在開發(fā)兩流體四場(chǎng)模型。

關(guān)于回路與組件在建模時(shí)被簡(jiǎn)化所導(dǎo)致的問題，目前的解決方法仍是對(duì)程序中的物理模型進(jìn)行改進(jìn)。如RELAP5程序中設(shè)置了許多特殊過程模型，這些模型可以由使用者決定是否在某些組件中開啟，來更精確地模擬這些組件中由于節(jié)點(diǎn)簡(jiǎn)化可能被忽略的物理現(xiàn)象。

除了對(duì)系統(tǒng)程序本身進(jìn)行改進(jìn)外，一些熱工水力系統(tǒng)程序也開始與嚴(yán)重事故分析程序、子通道程序、計(jì)算流體力學(xué)（CFD）程序、中子物理計(jì)算程序等耦合來形成具備多尺度、多物理場(chǎng)模擬能力的計(jì)算平臺(tái)，提高程序的計(jì)算能力。例如，將RELAP5程序與SCDAP程序耦合來將RELAP5的計(jì)算范圍擴(kuò)展至嚴(yán)重事故過程；將MARS程序與CUPID程序耦合來對(duì)局部現(xiàn)象進(jìn)行更精確的模擬；將GINKGO程序與堆芯程序COCO耦合來對(duì)堆芯進(jìn)行更準(zhǔn)確的模擬與建模等。

4 結(jié)論

本文介紹了目前國(guó)內(nèi)外熱工水力系統(tǒng)程序的最新發(fā)展現(xiàn)狀，并對(duì)未來熱工水力系統(tǒng)程序的發(fā)展方向進(jìn)行了探討，主要結(jié)論如下：

（1）國(guó)外一些成熟的熱工水力系統(tǒng)程序已經(jīng)可以較全面地對(duì)核電廠的各種系統(tǒng)瞬態(tài)及事故工況進(jìn)行計(jì)算。其中部分程序中加入了三維模型來模擬堆芯中的復(fù)雜現(xiàn)象，并設(shè)置了多種工作流體以滿足新型反應(yīng)堆設(shè)計(jì)與分析的需求。

（2）國(guó)內(nèi)的熱工水力系統(tǒng)程序開發(fā)和驗(yàn)證已取得了較大的進(jìn)展，部分熱工水力系統(tǒng)程序在計(jì)算能力、適用范圍等方面已經(jīng)可以比肩國(guó)外的主流程序。但這些程序仍需要大量的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行驗(yàn)證，且經(jīng)過核安全局評(píng)估后才能正式應(yīng)用到核電廠的安全分析中。

（3）目前熱工水力系統(tǒng)程序的發(fā)展主要分為兩方面，一是對(duì)系統(tǒng)程序本身進(jìn)行改進(jìn)，通過修改程序中的物理模型、改進(jìn)兩相流方程等方式突破以往程序的限制，提升程序的計(jì)算能力；二是與其他類型的程序進(jìn)行耦合來擴(kuò)展程序的計(jì)算范圍，或?qū)植楷F(xiàn)象進(jìn)行更精確的模擬。