盧 迪,夏榜樣
(中國核動(dòng)力研究設(shè)計(jì)院 核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室,四川 成都 610041)
超臨界水冷堆(SCWR)是最具發(fā)展前景的第4代核能系統(tǒng)之一,不僅具有機(jī)組熱效率高、系統(tǒng)簡化等優(yōu)點(diǎn)[1],且充分借鑒和繼承了傳統(tǒng)壓水堆、沸水堆和超臨界火電設(shè)計(jì)技術(shù)。目前,各核電大國均將SCWR 作為輕水堆的后續(xù)發(fā)展方向,制定了相應(yīng)的研究計(jì)劃,提出多種設(shè)計(jì)概念。目前,SCWR 正處于初步概念設(shè)計(jì)階段,數(shù)值分析模擬對(duì)于堆芯設(shè)計(jì)及物理特性分析至關(guān)重要。由于SCWR 具有強(qiáng)烈的物理-熱工耦合特性,堆芯設(shè)計(jì)更加復(fù)雜,導(dǎo)致傳統(tǒng)的壓水堆堆芯燃料管理程序無法滿足SCWR 堆芯設(shè)計(jì)要求。因此,國內(nèi)外雖已提出了多種SCWR概念設(shè)計(jì)方案,但適用于SCWR 的基于物理-熱工耦合的堆芯燃料管理程序相對(duì)較少。
本文針對(duì)SCWR 堆芯特點(diǎn),開發(fā)基于節(jié)塊法的超臨界堆芯燃料管理程序Xpack。
Xpack程序由3個(gè)模塊構(gòu)成,分別為中子學(xué)程序模塊、熱工水力程序模塊及物理-熱工耦合計(jì)算流程構(gòu)建模塊。
1)中子學(xué)程序模塊
中子學(xué)程序模塊由數(shù)據(jù)庫模塊、組件計(jì)算模塊、少群截面擬合模塊及臨界-燃耗計(jì)算模塊4個(gè)子模塊構(gòu)成。它的作用是實(shí)現(xiàn)程序系統(tǒng)的中子學(xué)計(jì)算功能,并為熱工水力分析模塊提供所需的輸入?yún)?shù)。其中,數(shù)據(jù)庫模塊由3個(gè)截面數(shù)據(jù)庫組成,包括快群數(shù)據(jù)庫、熱群數(shù)據(jù)庫及MCROSS數(shù)據(jù)庫,可提供1×10-5eV~10 MeV能量范圍內(nèi)300種核素的107群截面數(shù)據(jù)。組件計(jì)算模塊是基于碰撞概率法的燃料組件均勻化程序,可直接進(jìn)行全組件輸運(yùn)-燃耗均勻化計(jì)算,將107群截面歸并為2群,為少群截面擬合模塊提供必要的少群截面常數(shù)。少群截面擬合模塊首先將組件均勻化計(jì)算得到少群截面的離散數(shù)值制成關(guān)于慢化劑溫度、冷卻劑溫度及燃耗深度的列表,然后通過對(duì)列表進(jìn)行線性插值,得到相應(yīng)慢化劑溫度、冷卻劑溫度和燃耗深度下的均勻化截面。臨界-燃耗計(jì)算模塊可進(jìn)行多維堆芯臨界-燃耗計(jì)算,它的臨界計(jì)算基于擴(kuò)散理論,在求解堆芯臨界-擴(kuò)散方程時(shí)采用的數(shù)值算法是第二類邊界條件格林函數(shù)節(jié)塊法,并在計(jì)算給出堆芯組件平均功率后進(jìn)行組件精細(xì)功率重構(gòu)。第二類邊界條件格林函數(shù)節(jié)塊法通過引入第二類邊界格林函數(shù)將微分?jǐn)U散方程轉(zhuǎn)化為積分方程,在較寬的網(wǎng)格劃分下即可獲得足夠精確的擴(kuò)散方程數(shù)值解,從而大幅減少程序計(jì)算所需時(shí)間,使Xpack程序具有較高的計(jì)算效率。
2)熱工水力程序模塊
熱工水力程序模塊是基于單通道模型的熱工水力計(jì)算程序,其作用是實(shí)現(xiàn)Xpack程序熱工水力計(jì)算分析功能。由于SCWR 具有強(qiáng)烈的物理-熱工耦合效應(yīng),因此熱工水力計(jì)算對(duì)于SCWR 堆芯設(shè)計(jì)及燃料管理非常重要,該模塊的作用主要包括評(píng)價(jià)熱工水力反饋對(duì)中子學(xué)計(jì)算的影響及獲得堆芯設(shè)計(jì)所需的相關(guān)熱工水力特性參數(shù)。熱工水力分析程序模塊能實(shí)現(xiàn)對(duì)每個(gè)組件的冷卻劑通道進(jìn)行獨(dú)立的流量搜索與分配,它的目的在于:(1)搜索出最小給水流量;(2)計(jì)算堆芯各通道的冷卻劑、慢化劑溫度分布;(3)計(jì)算最大包殼溫度,確保滿足設(shè)計(jì)限值。程序計(jì)算流程分為3 步:(1)利用最大功率分布,按照給定的包殼溫度限值搜索最小給水流量,該流量在后面兩步的計(jì)算中作為已知量處理;(2)利用搜索出的給水流量及組件平均功率分布計(jì)算冷卻劑、慢化劑的溫度分布;(3)利用搜索出的給水流量及組件最大功率分布計(jì)算最大包殼溫度,以檢驗(yàn)是否符合設(shè)計(jì)限值。
3)物理-熱工耦合計(jì)算流程構(gòu)建模塊
物理-熱工耦合計(jì)算流程構(gòu)建模塊的作用是將Xpack程序系統(tǒng)的各子程序模塊按計(jì)算流程連接起來,并設(shè)定各子程序所需的輸入?yún)?shù)及運(yùn)行環(huán)境。通過各模塊的交互運(yùn)行,實(shí)現(xiàn)整個(gè)壽期內(nèi)的物理-熱工耦合計(jì)算,最終完成平衡循環(huán)堆芯方案的搜索及分析計(jì)算。該模塊能實(shí)現(xiàn)每燃耗步內(nèi)的物理-熱工耦合,以模擬SCWR強(qiáng)烈的物理-熱工耦合效應(yīng),從而滿足超臨界水冷堆堆芯設(shè)計(jì)的要求。程序?qū)崿F(xiàn)平衡循環(huán)堆芯設(shè)計(jì)的流程如圖1所示。
圖1 Xpack程序平衡循環(huán)堆芯設(shè)計(jì)流程Fig.1 Flow chart for equilibrium core design of Xpack code
中國核動(dòng)力研究設(shè)計(jì)院提出了百萬千瓦級(jí)SCWR 概念設(shè)計(jì)方案CSR1000[2],其堆芯活性區(qū)由157個(gè)燃料組件構(gòu)成,組件尺寸為23.9cm×23.9cm,活性區(qū)高度為420cm。詳細(xì)堆芯參數(shù)列于表1。采用Xpack 程序?qū)SR1000 堆芯方案在整個(gè)壽期內(nèi)進(jìn)行物理-熱工耦合計(jì)算,計(jì)算結(jié)果與OKUMURA 等[3]開發(fā)的SRAC2K6/SPROD 程序進(jìn)行比較分析。在程序求解臨界-擴(kuò)散方程時(shí),Xpack程序采用節(jié)塊法,計(jì)算時(shí)將每一組件劃分為4(2×2)個(gè)節(jié)塊。SRAC2K6/SPROD 程序采用細(xì)網(wǎng)有限差分法,為獲得足夠精確的數(shù)值解,需進(jìn)行細(xì)網(wǎng)劃分,將每個(gè)組件劃分為324(18×18)個(gè)均勻網(wǎng)格。
堆芯全提棒下keff的計(jì)算結(jié)果列于表2。圖2示出壽期初、壽期末堆芯平均功率的比較。
表1 CSR1000堆芯總體參數(shù)Table 1 General parameters of CSR1000core
表2 Xpack和SRAC2K6/SPROD程序計(jì)算得到的keff對(duì)比Table 2 Contrast of keffcalculated by Xpack and SRAC2K6/SPROD codes
從表2、圖2 可看出,在整個(gè)燃耗壽期內(nèi),兩種計(jì)算方法的keff計(jì)算結(jié)果相對(duì)偏差最大值僅為0.049 0%,組件平均功率相對(duì)偏差最大值為1.121%,計(jì)算結(jié)果滿足偏差限值要求。由于Xpack程序采用了第二類邊界條件格林函數(shù)節(jié)塊法求解臨界擴(kuò)散方程,相對(duì)于采用了細(xì)網(wǎng)有限差分法的SRAC2K6/SPROD 程序,計(jì)算效率提高了60倍左右。綜上所述,Xpack程序不僅能夠滿足超臨界堆芯物理計(jì)算精度的要求,而且具有較高的計(jì)算效率。
圖2 平均功率分布對(duì)比Fig.2 Contrast of average power distribution
堆芯中子學(xué)計(jì)算結(jié)束后,Xpack程序運(yùn)行熱工水力程序模塊進(jìn)行相關(guān)的熱工水力分析計(jì)算,最終給出堆芯第Ⅱ流程冷卻劑流量分配方案、冷卻劑出口溫度分布等。圖3 示出Xpack程序和SRAC2K6/SPROD 程序計(jì)算得到的堆芯第Ⅱ流程冷卻劑相對(duì)流量分布對(duì)比。圖4示出計(jì)算得到的堆芯冷卻劑出口溫度分布對(duì)比。從圖3、4 可看出,與SRAC2K6/SPROD 程 序相比,Xpack程序計(jì)算得到的堆芯冷卻劑出口流量最大相對(duì)偏差為3.794%,出口溫度最大相對(duì)偏差為3.384%,較組件的平均功率偏差要大。原因是Xpack 程序進(jìn)行熱工水力計(jì)算時(shí),需要輸入組件精細(xì)功率分布,而組件精細(xì)功率的最大相對(duì)偏差要大于組件平均功率的相對(duì)偏差,因此,熱工水力計(jì)算結(jié)果較堆芯功率計(jì)算結(jié)果而言,偏差較大,但仍在允許范圍內(nèi),可認(rèn)為Xpack程序物理-熱工耦合的計(jì)算精度滿足超臨界堆芯計(jì)算要求。
圖3 冷卻劑相對(duì)流量分布對(duì)比Fig.3 Contrast of relative coolant flow rate distribution
圖4 冷卻劑出口溫度分布對(duì)比Fig.4 Contrast of coolant outlet temperature distribution
SCWR 作為第4 代核能系統(tǒng)的候選堆型之一,由于采用了超臨界水作為堆芯冷卻劑,使堆芯具有非常強(qiáng)烈的物理-熱工耦合特性,導(dǎo)致無法直接應(yīng)用傳統(tǒng)的壓水堆程序進(jìn)行超臨界堆芯燃料管理計(jì)算。針對(duì)SCWR 堆芯特點(diǎn)開發(fā)了Xpack 程序,程序包含中子學(xué)程序模塊、熱工水力程序模塊及物理-熱工耦合計(jì)算流程構(gòu)建模塊,能實(shí)現(xiàn)壽期內(nèi)每個(gè)燃耗步的物理-熱工耦合,滿足SCWR 堆芯計(jì)算的要求。數(shù)值計(jì)算結(jié)果表明:Xpack程序是高效和可靠的,適用于超臨界堆芯設(shè)計(jì)及性能分析,可為CSR1000研究提供高效的數(shù)值分析工具。
[1] TIM A,SUE I.Generation-Ⅳnuclear power:A review of the state of the science[J].Energy Policy,2008,36:4 323-4 330.
[2] 李翔,李慶,夏榜樣,等.超臨界水冷堆CSR1000總體設(shè)計(jì)研究[J].核動(dòng)力工程,3013,34(1):5-8.LI Xiang,LI Qing,XIA Bangyang,et al.Overview on overall design of CSR1000[J].Nuclear Power Engineering,2013,34(1):5-8(in Chinese).
[3] OKUMURA K.COREBN:A core burn-up calculation module for SRAC2006[M].Japan:Japan Atomic Energy Research Institute,2007.