鄭海偉 鄂彥志 唐延澤 傅 瑤 王納秀

?

HYSYS軟件在反應(yīng)堆系統(tǒng)仿真中的應(yīng)用探索

鄭海偉1,2鄂彥志1,2唐延澤3傅 瑤1王納秀1

1（中國(guó)科學(xué)院上海應(yīng)用物理研究所 嘉定園區(qū) 上海 201800）2（中國(guó)科學(xué)院大學(xué) 北京 100049）3（中國(guó)科學(xué)院近代物理研究所 蘭州 730000）

系統(tǒng)仿真軟件可以模擬運(yùn)行工況變化對(duì)系統(tǒng)整體運(yùn)行帶來(lái)的影響，在系統(tǒng)瞬態(tài)分析和安全研究中起著重要的作用。Aspen HYSYS軟件是世界知名的油氣過(guò)程仿真和優(yōu)化的系統(tǒng)軟件，具有強(qiáng)大的二次開(kāi)發(fā)功能，可以用于反應(yīng)堆系統(tǒng)仿真。在植入熔鹽物性、修改熔鹽換熱模型的基礎(chǔ)上，建立并調(diào)用點(diǎn)堆模型的動(dòng)態(tài)鏈接庫(kù)，嘗試將HYSYS與點(diǎn)堆耦合起來(lái)，彌補(bǔ)HYSYS無(wú)法對(duì)熔鹽堆等反應(yīng)堆進(jìn)行仿真的缺憾。在此基礎(chǔ)上，對(duì)中國(guó)科學(xué)院上海應(yīng)用物理研究所的熔鹽堆設(shè)計(jì)進(jìn)行了系統(tǒng)仿真，給出了熔鹽堆在不同的運(yùn)行工況下的系統(tǒng)響應(yīng)分析結(jié)果，并與RELAP5仿真結(jié)果進(jìn)行比較。結(jié)果表明，耦合程序有較高的可用性，能夠達(dá)到預(yù)期的效果。

熔鹽堆，HYSYS，耦合，系統(tǒng)仿真

釷基熔鹽堆(Thorium Molten Salt Reactor, TMSR)核能系統(tǒng)項(xiàng)目是中國(guó)科學(xué)院先導(dǎo)科技專(zhuān)項(xiàng)之一，其戰(zhàn)略性目標(biāo)是研發(fā)第四代裂變反應(yīng)堆核能系統(tǒng)[1]。隨著熔鹽堆系統(tǒng)研究的不斷深入，熔鹽堆系統(tǒng)分析的需求愈發(fā)凸顯，現(xiàn)有的反應(yīng)堆系統(tǒng)分析程序大多是針對(duì)壓水堆系統(tǒng)建立的兩相流一維流動(dòng)模型無(wú)法直接應(yīng)用到熔鹽堆的系統(tǒng)模擬當(dāng)中。需要一個(gè)單相一維多工質(zhì)的系統(tǒng)模擬程序用于熔鹽堆的功率升降、流量變化、負(fù)載變化等工況模擬，驗(yàn)證系統(tǒng)設(shè)計(jì)參數(shù)。

Aspen HYSYS是面向油氣生產(chǎn)、氣體處理和煉油工業(yè)的模擬、設(shè)計(jì)和性能監(jiān)測(cè)的流程模擬工具[2]，同時(shí)具有穩(wěn)態(tài)和動(dòng)態(tài)兩種仿真模式。其技術(shù)廣泛應(yīng)用于石油開(kāi)采、儲(chǔ)運(yùn)、天然氣加工、石油化工、精細(xì)化工、制藥、煉制等領(lǐng)域，在世界范圍內(nèi)石油化工模擬、仿真技術(shù)領(lǐng)域占主導(dǎo)地位[3]。該軟件具有豐富的用戶(hù)自定義和軟件通訊接口，具有良好的軟件的擴(kuò)展、二次開(kāi)發(fā)和與其它軟件的鏈接功能[4]。HYSYS具有強(qiáng)大的物性計(jì)算包和物性預(yù)測(cè)能力[5]，可以用作單相和多相一維多工質(zhì)的系統(tǒng)模擬。但是HYSYS中缺少反應(yīng)堆中子動(dòng)力學(xué)計(jì)算模塊，該模塊是反應(yīng)堆與化工系統(tǒng)的關(guān)鍵區(qū)別，要實(shí)現(xiàn)對(duì)反應(yīng)堆的仿真，需要將中子計(jì)算模塊加入到HYSYS中。HYSYS允許用戶(hù)按照自身需求建立物性包、增加新的反應(yīng)方程以及控制參數(shù)的變化等，具有開(kāi)發(fā)成為一款滿(mǎn)足反應(yīng)堆系統(tǒng)仿真需求的高效系統(tǒng)瞬態(tài)仿真軟件的條件和潛力。

本文對(duì)HYSYS進(jìn)行二次開(kāi)發(fā)，將中子動(dòng)力學(xué)計(jì)算模塊與HYSYS進(jìn)行耦合，并對(duì)耦合程序進(jìn)行對(duì)比驗(yàn)證，旨在使HYSYS程序能夠?qū)Ψ磻?yīng)堆系統(tǒng)各種穩(wěn)態(tài)和瞬態(tài)工況進(jìn)行可靠、準(zhǔn)確模擬，為反應(yīng)堆設(shè)計(jì)和安全分析提供幫助。

1 TMSR-SF1介紹

TMSR-SF1是10 MWth固態(tài)燃料釷基熔鹽堆的簡(jiǎn)稱(chēng)，是中國(guó)科學(xué)院TMSR專(zhuān)項(xiàng)的一部分。它具有發(fā)電效率高、安全性好、核燃料有效利用和防止核擴(kuò)散等優(yōu)點(diǎn)[6]。

1.1 TMSR-SF1系統(tǒng)和設(shè)備

TMSR-SF1包括如下主要系統(tǒng)和設(shè)備：堆本體和堆內(nèi)構(gòu)件、反應(yīng)堆冷卻劑系統(tǒng)及設(shè)備、安全相關(guān)系統(tǒng)、核輔助系統(tǒng)、測(cè)量和控制系統(tǒng)等，分為帶走堆芯熱量的一回路、中間換熱的二回路和將熱量帶走的空冷側(cè)三個(gè)部分。

堆本體由堆芯燃料球、石墨反射層、堆芯支撐機(jī)構(gòu)、石墨反射層壓緊機(jī)構(gòu)、堆芯圍筒、主容器、控制棒及控制棒驅(qū)動(dòng)機(jī)構(gòu)、堆內(nèi)測(cè)量機(jī)構(gòu)等組成。采用直徑為6 cm的TRISO (Tri-structural iso-tropic)包覆顆粒燃料球作為燃料。

一回路主要由主循環(huán)泵、中間換熱器(Intermediate Heat Exchanger, IHX)、溢流罐、管路和輔助設(shè)備組成，采用2LiF-BeF2鹽作為冷卻劑。二回路主要由循環(huán)泵、IHX二回路側(cè)、熔鹽-空氣換熱器、溢流罐、管路和輔助設(shè)備組成，采用FLiNaK鹽作為冷卻劑。

正常運(yùn)行時(shí)一回路中熔鹽冷卻劑以強(qiáng)迫循環(huán)的形式從堆芯帶出核裂變能，并通過(guò)雙熔鹽熱交換器把熱能從一回路熔鹽轉(zhuǎn)移到二回路熔鹽；然后二回路熔鹽中的熱量通過(guò)空氣散熱器轉(zhuǎn)移到大氣環(huán)境。

1.2 TMSR-SF1系統(tǒng)主要參數(shù)

釷基熔鹽堆的主要熱工水力設(shè)計(jì)參數(shù)[7]見(jiàn)表1。

表1 TMSR-SF1主要熱工水力參數(shù)
Table 1 Main thermal-hydraulic parameters of TMSR-SF1.

參數(shù) Parameter

數(shù)值 Value

反應(yīng)堆堆芯熱功率

Reactor core thermal power / MW

10

一回路質(zhì)量流量

Flow rate of first loop / kg?s

150

堆芯進(jìn)口溫度

Core inlet temperature / K

873.15

堆芯出口溫度

Core outlet temperature / K

901.15

堆芯平均流速

Core mean velocity / m?s

0.06

二回路冷卻劑流量

Coolant flow rate of second loop / kg?s

161.0

IHX二次側(cè)進(jìn)口溫度

Secondary side inlet temperature of IHX / K

803.15

IHX二次側(cè)出口溫度

Secondary side outlet temperature of IHX / K

836.15

2 仿真模型

該耦合程序以Aspen HYSYS和中子動(dòng)力學(xué)計(jì)算程序?yàn)榛A(chǔ)，模擬反應(yīng)堆系統(tǒng)在不同運(yùn)行工況下的響應(yīng)情況，重點(diǎn)觀測(cè)系統(tǒng)溫度等參數(shù)隨時(shí)間變化。在中子動(dòng)力學(xué)計(jì)算中，采用6群中子動(dòng)力學(xué)方程模型，并考慮燃料溫度變化對(duì)堆芯溫度反饋的影響，具有較高的計(jì)算精度和效率，可以用于各種堆型的系統(tǒng)仿真。

2.1 HYSYS計(jì)算模型

采用HYSYS程序進(jìn)行熱工水力計(jì)算，在HYSYS中，所有單元操作均采用集總參數(shù)模型。這里主要對(duì)仿真建模中涉及的數(shù)學(xué)計(jì)算模型進(jìn)行簡(jiǎn)要的介紹，而不再對(duì)用到的各個(gè)元件和設(shè)備的具體的計(jì)算原理進(jìn)行詳細(xì)的描述。守恒關(guān)系是HYSYS中數(shù)學(xué)計(jì)算的基礎(chǔ)。下面將給出質(zhì)量、組分、能量等的動(dòng)態(tài)守恒關(guān)系式，并簡(jiǎn)述軟件的求解策略。

2.1.1 物料守恒

質(zhì)量守恒的基本關(guān)系式為：

式中：i、o分別為進(jìn)、出口流量；i、o分別為進(jìn)、出口的流體密度；為體積。

式(1)只是Aspen HYSYS中所使用方程的簡(jiǎn)化形式，在軟件中考慮了反應(yīng)、蒸發(fā)密度變化等更復(fù)雜的情況。

2.1.2 組分守恒

組分守恒的基本關(guān)系為：系統(tǒng)中某一組分的變化等于流入系統(tǒng)的該組分減去流出該系統(tǒng)的組分再加上該組分發(fā)生反應(yīng)的量。

對(duì)于多組分物流，組分的平衡關(guān)系式如下：

式中：Ci為入口物流中組分的濃度；Co為出口物流中的濃度；R為組分的反應(yīng)率。

在反應(yīng)堆中，中子的影響對(duì)物質(zhì)的改變非常小，而堆中物質(zhì)相容性非常好，化學(xué)反應(yīng)可以忽略不計(jì)，認(rèn)為發(fā)生反應(yīng)的量為零。

2.1.3 能量守恒

能量流入或流出系統(tǒng)是通過(guò)對(duì)流或者傳導(dǎo)，通過(guò)邊界向系統(tǒng)中增加能量是通過(guò)熱傳導(dǎo)或輻射。

對(duì)于有散熱的連續(xù)反應(yīng)器，可以采用如下方程：

式中：、、分別為單位質(zhì)量的內(nèi)能、動(dòng)能和勢(shì)能；為系統(tǒng)單位時(shí)間內(nèi)所做的軸功；o為容器的壓力；i為進(jìn)料的壓力；為通過(guò)邊界增加的熱量；r為反應(yīng)熱。

2.1.4 求解策略

HYSYS中采用歐拉隱式求解方法，對(duì)體積（壓力-流量）、能量和組分的求解并不同步，以減少計(jì)算機(jī)求解工作量，他們分別對(duì)應(yīng)的默認(rèn)求解頻率分別為1、2、10個(gè)時(shí)間步。

2.2 中子計(jì)算模型

2.2.1 點(diǎn)堆模型

考慮到計(jì)算的速度、精確性以及計(jì)算模型的復(fù)雜程度，選擇廣泛采用的具有6組緩發(fā)中子的反應(yīng)堆中子動(dòng)力學(xué)方程模型：

總反應(yīng)性表示為：

式中：0為初始反應(yīng)性；exp為控制機(jī)構(gòu)動(dòng)作等原因所引起的反應(yīng)性變化，通過(guò)該項(xiàng)可以引入反應(yīng)性瞬變；f、r、c和c0分別表示燃料、反射層、下降管和上升管中冷卻劑按質(zhì)量取的平均溫度；f、m、r和c分別表示燃料、慢化劑、反射層、冷卻劑的溫度系數(shù)。

程序中采用gear算法求解中子動(dòng)力學(xué)方程，該算法可用于求解剛性微分方程組。程序模型中將反射層、燃料、冷卻劑等劃分成多個(gè)節(jié)點(diǎn)，反射層、燃料、冷卻劑等每個(gè)節(jié)點(diǎn)相關(guān)的量，在程序中默認(rèn)為均勻劃分。

2.2.2 耦合方法

選擇Automation方式來(lái)實(shí)現(xiàn)相應(yīng)的耦合功能。將中子動(dòng)力學(xué)計(jì)算程序改編成動(dòng)態(tài)鏈接庫(kù)文件(Dynamic Link Library, DLL)，以供其他程序在需要的時(shí)候調(diào)用，該DLL可以返回修正的堆芯功率。進(jìn)行HYSYS中建立反應(yīng)堆系統(tǒng)模型時(shí)，用加熱器代替堆芯對(duì)整個(gè)系統(tǒng)進(jìn)行加熱，堆芯部分軸向分成多段，加載不同的反應(yīng)堆功率，熱量通過(guò)換熱器傳遞到二回路中。讀取時(shí)刻HYSYS 中流經(jīng)堆芯的溫度()及添加的反應(yīng)性()后，傳遞到DLL中；DLL中進(jìn)行燃料的傳熱計(jì)算和堆芯燃料、流體的溫度反饋，返回修正后的堆芯功率到HYSYS中，經(jīng)重新分配作為下一時(shí)刻+1的堆芯功率(,1)，再利用HYSYS進(jìn)行熱工水力計(jì)算。

3 仿真驗(yàn)證

為驗(yàn)證仿真計(jì)算結(jié)果的準(zhǔn)確性，本文針對(duì)熔鹽堆系統(tǒng)進(jìn)行仿真，并將仿真計(jì)算結(jié)果與RELAP5的計(jì)算結(jié)果進(jìn)行比較。RELAP5程序是由美國(guó)愛(ài)達(dá)荷國(guó)家工程實(shí)驗(yàn)室開(kāi)發(fā)的熱工水力系統(tǒng)分析程序，被廣泛應(yīng)用于事故分析和安全評(píng)審等方面[8]。在空氣冷卻器流量增加和流經(jīng)空氣冷卻器的冷空氣溫度升高兩種情況下，對(duì)流經(jīng)堆芯的熔鹽進(jìn)出口溫度和反應(yīng)堆的反應(yīng)性進(jìn)行分析，模擬時(shí)間均為2000 s。

建立的系統(tǒng)如圖1所示。初始時(shí)刻反應(yīng)堆均處于滿(mǎn)功率運(yùn)行狀態(tài)，堆芯額定功率為10 MW，HYSYS軟件模擬的反應(yīng)堆滿(mǎn)功率運(yùn)行的進(jìn)出口溫度分別為875.6 K、904.5 K。

圖1 HYSYS模擬系統(tǒng)圖Fig.1 System diagram of HYSYS simulation.

3.1 空氣冷卻器流量增加

具體工況描述如下：在0時(shí)刻，熔鹽空氣冷卻器空氣側(cè)流量增加了3%，并維持在新流量水平上運(yùn)行，其他條件保持不變?？諝饬髁康脑黾油ㄟ^(guò)增大空氣側(cè)進(jìn)口壓力使進(jìn)出口壓差增大的方法實(shí)現(xiàn)。

HYSYS仿真結(jié)果與RELAP5仿真結(jié)果對(duì)比如圖2所示。無(wú)耦合情況下模擬堆芯功率沒(méi)有發(fā)生變化。從圖2中可以看出，有無(wú)耦合對(duì)HYSYS模擬結(jié)果影響較大，無(wú)耦合情況下HYSYS模擬所得到的堆芯進(jìn)出口溫度均呈現(xiàn)一直下降的狀態(tài)，而耦合后HYSYS模擬的堆芯進(jìn)出口溫度則有溫度的短時(shí)間小幅上升情況出現(xiàn)，且無(wú)耦合模擬結(jié)果溫度變化速率和幅度明顯快于耦合后模擬結(jié)果，這是因?yàn)轳詈虾蠹尤肓素?fù)溫度反饋特性，對(duì)溫度變化幅度和速率均產(chǎn)生了抑制變化的效果，使得溫度的降低速率和幅度變慢，逐漸趨于新的平衡狀態(tài)。耦合后HYSYS模擬結(jié)果與RELAP5模擬結(jié)果的總體趨勢(shì)是一致的：初始時(shí)均有小幅的溫度降低再升高過(guò)程，然后溫度緩慢降低。出現(xiàn)該現(xiàn)象的原因是：熔鹽空氣冷卻器的流量增大使得該換熱器的空氣側(cè)換熱能力增強(qiáng)，從二回路帶走的熱量增多，從而導(dǎo)致流經(jīng)該換熱器的二回路FLiNaK熔鹽溫度降低；二回路與一回路通過(guò)雙熔鹽換熱器進(jìn)行熱量傳遞，二回路熔鹽溫度的降低引起一回路流經(jīng)反應(yīng)堆堆芯的FLiBe熔鹽溫度降低；由于反應(yīng)堆的負(fù)溫度反饋特性，溫度的降低引入正的反應(yīng)性，堆芯功率增加，使得短時(shí)間內(nèi)對(duì)一回路熔鹽加熱功率大于冷卻功率，熔鹽溫度升高；后由于溫度的變化量減小，正反應(yīng)性逐漸減小，冷卻功率大于加熱功率，溫度逐漸降低。

但是由于初始時(shí)刻，HYSYS進(jìn)出口溫度較RELAP5略低，在整個(gè)模擬時(shí)間內(nèi)，HYSYS模擬得到的堆芯進(jìn)出口溫度均低于RELAP5的進(jìn)出口溫度。而且圖2中HYSYS溫度變化曲線斜率略大于RELAP5，即溫度變化略快于RELAP5，HYSYS耦合的負(fù)溫度反應(yīng)性系數(shù)大于RELAP5的反應(yīng)性系數(shù)。這種情況的出現(xiàn)與熱工計(jì)算及溫度權(quán)重等因素有關(guān)。

圖2 空氣流量增加時(shí)堆芯進(jìn)出口溫度隨時(shí)間的變化Fig.2 Temperature change of the core inlet and outletvs.time with air flow rate increasing.

3.2 空氣冷卻器空氣溫度升高

具體工況描述如下：在0時(shí)刻，熔鹽空氣冷卻器空氣側(cè)溫度增加了10 K，空氣流量保持不變。分析空氣進(jìn)口溫度變化對(duì)反應(yīng)堆運(yùn)行的影響，溫度增加可以直接在HYSYS中通過(guò)設(shè)定進(jìn)口溫度實(shí)現(xiàn)。

該工況下，HYSYS仿真結(jié)果與RELAP5仿真結(jié)果對(duì)比如圖3所示，無(wú)耦合情況下模擬堆芯功率沒(méi)有發(fā)生變化。從圖3中可以看出，有無(wú)耦合對(duì)HYSYS模擬結(jié)果影響較大，無(wú)耦合情況下HYSYS模擬所得到的堆芯進(jìn)出口溫度均呈現(xiàn)不斷升高的狀態(tài)，而耦合后HYSYS模擬的堆芯進(jìn)出口溫度則有溫度的短時(shí)間小幅下降情況出現(xiàn)，且無(wú)耦合模擬結(jié)果溫度變化速率和幅度明顯快于耦合后模擬結(jié)果，這是因?yàn)轳詈虾蠹尤肓素?fù)溫度反饋特性，對(duì)溫度變化幅度和速率均產(chǎn)生了抑制變化的效果，使得溫度的升高速率和幅度變慢，逐漸趨于新的平衡狀態(tài)。耦合后HYSYS模擬結(jié)果與RELAP5模擬結(jié)果的總體趨勢(shì)是一致的：兩者總體變化趨勢(shì)相同，均出現(xiàn)初始時(shí)刻溫度短暫升高后降低、再不斷升高的情況，且溫度升高速度不斷降低，溫度變化曲線趨于平緩。出現(xiàn)該現(xiàn)象的原因?yàn)椋涸谙嗤目諝饬髁壳闆r下，突然升高熔鹽-空氣換熱器進(jìn)口溫度，使得換熱器的冷卻能力降低；二次側(cè)FliNaK熔鹽溫度升高，進(jìn)而導(dǎo)致一回路總FliBe熔鹽溫度升高；由于該反應(yīng)堆具有負(fù)溫度系數(shù)，溫度的升高引入負(fù)反應(yīng)性，堆芯加熱功率減小，初始時(shí)功率降低量大于冷卻能力降低量，即發(fā)熱量小于冷卻量，故一回路熔鹽溫度降低，進(jìn)而導(dǎo)致一回路溫度升高量降低，引入反應(yīng)性減小，溫度降低速度變慢；溫度降至最小值時(shí)，冷卻能力小于堆芯功率，一回路溫度身高，由于負(fù)反饋的影響，溫度升高逐漸變慢，但加熱功率仍占主導(dǎo)地位。

圖3 入口空氣溫度升高時(shí)堆芯進(jìn)出口溫度隨時(shí)間的變化Fig.3 Temperature change of the core inlet and outletvs.time with inlet temperature of air increasing.

但在溫度從最低值開(kāi)始升高時(shí)，RELAP5的升溫速率大于HYSYS升溫速率，可以得出，HYSYS耦合的負(fù)溫度反應(yīng)性系數(shù)大于RELAP5的反應(yīng)性系數(shù)，這與熱工計(jì)算及溫度權(quán)重等因素有關(guān)。

3.3 階躍引入反應(yīng)性

通常情況下反應(yīng)堆測(cè)量系統(tǒng)和反應(yīng)性控制系統(tǒng)存在一定的誤差，導(dǎo)致實(shí)際運(yùn)行時(shí)可能出現(xiàn)棒位略高于臨界棒位而引入少量反應(yīng)性。這里就針對(duì)反應(yīng)堆引入正反應(yīng)性的情況進(jìn)行仿真。假設(shè)反應(yīng)堆處于額定功率運(yùn)行狀態(tài)，在0時(shí)刻突然引入5 pcm的反應(yīng)性，分析反應(yīng)性小幅度擾動(dòng)對(duì)反應(yīng)堆運(yùn)行的影響。

引入5 pcm反應(yīng)性后HYSYS耦合仿真結(jié)果與RELAP5仿真結(jié)果的對(duì)比如圖4所示。

圖4 引入正反應(yīng)性時(shí)堆芯進(jìn)出口溫度隨時(shí)間的變化Fig.4 Temperature change of the core inlet and outletvs.time with additional positive reactivity introduced.

以上幾個(gè)工況計(jì)算結(jié)果的對(duì)比表明，HYSYS計(jì)算的結(jié)果與RELAP5的結(jié)果趨勢(shì)相同，計(jì)算結(jié)果相差不大，證明了通過(guò)HYSYS耦合進(jìn)行系統(tǒng)仿真的可行性。

4 結(jié)語(yǔ)

基于熔鹽堆系統(tǒng)與普通化工系統(tǒng)的不同，本文對(duì)HYSYS軟件進(jìn)行深度擴(kuò)展，在加入熔鹽物性和修正換熱系數(shù)計(jì)算公式的基礎(chǔ)上，將點(diǎn)堆模塊與HYSYS軟件進(jìn)行耦合，引入堆芯反饋特性。通過(guò)軟件模擬和RELAP5仿真計(jì)算的對(duì)比結(jié)果來(lái)看，擴(kuò)展后的HYSYS 軟件成功地模擬了熔鹽堆系統(tǒng)的瞬態(tài)過(guò)程。軟件計(jì)算的回路溫度隨時(shí)間的變化趨勢(shì)、溫度峰值大小等均與RELAP5計(jì)算結(jié)果符合較好。

HYSYS 軟件能夠?qū)崿F(xiàn)對(duì)熔鹽堆系統(tǒng)的瞬態(tài)過(guò)程比較準(zhǔn)確地模擬，證明軟件適用于反應(yīng)堆回路瞬態(tài)仿真，體現(xiàn)了其強(qiáng)大的可二次開(kāi)發(fā)能力，顯示了較強(qiáng)的可靠性。在缺乏熔鹽堆系統(tǒng)專(zhuān)用仿真軟件的現(xiàn)狀下，對(duì)HYSYS軟件功能進(jìn)行深度挖掘，開(kāi)發(fā)出具有較完善功能和模塊的熔鹽堆系統(tǒng)仿真軟件是熔鹽堆研究的一種有益探索，是對(duì)熔鹽堆系統(tǒng)仿真軟件的一個(gè)有效補(bǔ)充，對(duì)現(xiàn)階段熔鹽堆的熱工水力和安全研究具有重要的意義。

HYSYS軟件的耦合不僅利用了HYSYS軟件物性計(jì)算能力強(qiáng)大和模擬計(jì)算準(zhǔn)確等優(yōu)勢(shì)，而且彌補(bǔ)了通用熱流體計(jì)算程序不能實(shí)現(xiàn)反應(yīng)性反饋計(jì)算的不足，還能避免許多熱工分析程序無(wú)法實(shí)現(xiàn)不同物性冷卻劑混合、流體接觸空氣等情況的缺點(diǎn)，可以進(jìn)行一些特殊現(xiàn)象的模擬。后續(xù)可以利用HYSYS軟件進(jìn)行更多典型工況和特殊物理現(xiàn)象的仿真，開(kāi)發(fā)熔鹽物性計(jì)算方法，豐富和完善其在熔鹽堆中的應(yīng)用。

致謝 該工作的完成得到了中國(guó)科學(xué)院上海應(yīng)用物理研究所TMSR中心反應(yīng)堆物理部計(jì)算組職工何龍的幫助，在此表示衷心感謝。

1 江綿恒, 徐洪杰, 戴志敏. 未來(lái)先進(jìn)核裂變能——TMSR核能系統(tǒng)[J]. 中國(guó)科學(xué)院院刊, 2012, 27(3): 366?374. DOI: 10.3969/j.issn.1000-3045.2012.03.016
JIANG Mianheng, XU Hongjie, DAI Zhimin. Advanced fission energy program-TMSR nuclear energy system[J]. Bulletin of Chinese Academy of Sciences, 2012, 27(3): 366?374. DOI: 10.3969/j.issn.1000-3045.2012.03.016

2 蔡飛云, 劉江楠, 張強(qiáng). 流程模擬在焦化吸收穩(wěn)定系統(tǒng)中的應(yīng)用[J]. 化工自動(dòng)化及儀表, 2014, 41(3): 311?314, 325. DOI: 10.3969/j.issn.1000-3932.2014.03.019
CAI Feiyun, LIU Jiangnan, ZHANG Qiang. Application of process simulation in stable system for coking absorption[J]. Control and Instruments in Chemical Industry, 2014, 41(3): 311?314, 325. DOI: 10.3969/ j.issn.1000-3932.2014.03.019

3 金大天, 曹義鳴, 王麗娜, 等. HYSYS的二次開(kāi)發(fā)及其在多組分氣體膜分離過(guò)程模擬中的應(yīng)用[J]. 膜科學(xué)與技術(shù), 2012, 32(5): 87?91. DOI: 10.3969/j.issn.1007-8924. 2012.05.016
KIM Daechun, CAO Yiming, WANG Lina,. Secondary development of HYSYS and its application in the simulation of multi-component gas separation by membrane[J]. Membrane Science and Technology, 2012, 32(5): 87?91. DOI: 10.3969/j.issn.1007-8924.2012.05. 016

4 唐延澤, 鄭海偉, 陳玉爽, 等. HYSYS軟件應(yīng)用于熔鹽冷卻系統(tǒng)瞬態(tài)分析的可行性研究[J]. 核技術(shù), 2015, 38(4): 040601. DOI: 10.11889/j.0253-3219.2015.hjs.38. 040601
TANG Yanze, ZHENG Haiwei, CHEN Yushuang,. Feasibility study of HYSYS code for molten salt cooled system in dynamics analysis[J]. Nuclear Techniques, 2015, 38(4): 040601. DOI: 10.11889/j.0253-3219.2015.hjs.38. 040601

5 郭廣智. 石油化工動(dòng)態(tài)模擬軟件HYSYS[J]. 石油化工設(shè)計(jì), 1997, 14(3): 29?33
GUO Guangzhi. Petrochemical dynamic simulation software HYSYS[J]. Petrochemical Design, 1997, 14(3): 29?33

6 秋穗正, 張大林, 蘇光輝, 等. 新概念熔鹽堆的固有安全性及相關(guān)關(guān)鍵問(wèn)題研究[J]. 原子能科學(xué)技術(shù), 2009, 43(S1): 64?75
QIU Suizheng, ZHANG Dalin, SU Guanghui,. Research on inherent safety and relative key issues of a molten salt reactor[J]. Atomic Energy Science and Technology, 2009, 43(S1): 64?75

7 TMSR卓越中心. 10 MWt固態(tài)燃料釷基熔鹽實(shí)驗(yàn)堆概念設(shè)計(jì)報(bào)告[R]. 上海: 中國(guó)科學(xué)院上海應(yīng)用物理研究所, 2014
Center for Excellence in TMSR. The conceptual design report on 10 MWt solid thorium-based molten salt experiment reactor[R]. Shanghai: Shanghai Institute of Applied Physics, Chinese Academy of Sciences, 2014

8 劉余, 張虹. RELAP5程序耦合接口的開(kāi)發(fā)[J]. 核動(dòng)力工程, 2009, 30(6): 38?40, 45
LIU Yu, ZHANG Hong. Development of RELAP5 coupling interface[J]. Nuclear Power Engineering, 2009, 30(6): 38?40, 45

國(guó)家自然科學(xué)基金(No.91326201)、中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專(zhuān)項(xiàng)(No.XDA02001002)資助

Supported by National Natural Science Foundation of China (No.91326201), Strategic Priority Program of Chinese Academy of Sciences

(No.XDA02001002)

Applied exploration of HYSYS code in reactor system simulation

ZHENG Haiwei1,2E Yanzhi1,2TANG Yanze3FU Yao1WANG Naxiu1

1(Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Jiading Campus, Shanghai 201800, China)2(University of Chinese Academy of Sciences, Beijing 100049, China)3(Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China)

Background: The system code can be used to simulate the effects of the whole system with different operation conditions, which plays an important role in dynamic analysis and safety study. As a world-known modeling tool for process simulation and optimization of oil and gas, Aspen HYSYS has the ability of secondary development for nuclear reactor simulation. Purpose: The aim is to make the coupling of HYSYS and point reactor model and verify its efficiency. Methods: Based on implanted and modified molten salt properties, HYSYS can show the features of molten salt accurately. A dynamic link library is built based on the point reactor model. Then the DLL (Dynamic Link Library) is called by HYSYS to calculate the reactivity which is caused by changes of operation conditions. Some operating conditions of thorium molten salt reactor are simulated and then compared with RELAP5. Results: The results of HYSYS is in good agreement with that of RELAP5. Conclusion: The coupled code is available in simulating molten salt-cooled system.

Molten salt reactor, HYSYS, Coupling, System simulation

ZHENG Haiwei, male, born in 1989, graduated from Harbin Engineering University in 2013, master student, focusing on reactor thermal hydraulic

WANG Naxiu, E-mail: wangnaxiu@sinap.ac.cn

TL333

10.11889/j.0253-3219.2016.hjs.39.060602

鄭海偉，男，1989年出生，2013年畢業(yè)于哈爾濱工程大學(xué)，現(xiàn)為碩士研究生，研究領(lǐng)域?yàn)榉磻?yīng)堆熱工水力學(xué)

王納秀，E-mail: wangnaxiu@sinap.ac.cn

2016-03-07，

2016-03-22