離散坐標(biāo)方法在反射層參數(shù)計(jì)算中的應(yīng)用研究

陳長(zhǎng)　肖鋒

【摘 要】本文介紹離散坐標(biāo)方法在堆芯反射層參數(shù)計(jì)算中的應(yīng)用研究。針對(duì)使用TPFAP程序近似計(jì)算反射層參數(shù)的情況，利用二維離散坐標(biāo)法程序DOT的輸運(yùn)計(jì)算模塊替換了TPFAP程序內(nèi)的穿透幾率法子程序TPXY，開(kāi)發(fā)了反射層參數(shù)計(jì)算程序TPREF。TPREF程序計(jì)算得到的反射層參數(shù)與TPFAP程序較為接近，表明程序開(kāi)發(fā)是正確的。

【關(guān)鍵詞】離散坐標(biāo)方法；反射層參數(shù)

【Abstract】This paper introduces research on the reflector parameters calculation using the discrete ordinates method. The reflector parameters calculated by TPFAP. In this instance， the code （TPREF） using the two-dimensional discrete ordinates code DOT which replaces the transmission probability method （TPM） code TPXY is encoded. The results of reflector parameters demonstrate that this method is correct.

【Key words】Discrete ordinates method；Reflector parameters

0 引言

反射層是指包圍在堆芯外面用以反射從堆芯泄漏出來(lái)的中子的材料。反射層的主要作用有兩個(gè)方面，一是可以減少中子泄漏，提高中子、燃料利用率；二是使得堆芯通量分布更為平坦。良好的慢化劑材料通常也是良好的反射層材料，如：輕水、重水、石墨和鈹?shù)取Ｔ趬核阎校ǔＲ暂p水作為反射層。為了約束燃料組件在堆芯中的徑向位置，在堆芯活性區(qū)邊界布置有2～4cm的圍板。圍板的材料通常是不銹鋼，由于不銹鋼對(duì)中子有較強(qiáng)的吸收能力，堆芯功率分布有較大的影響。圍板的存在也使得圍板/反射層參數(shù)計(jì)算變得困難，國(guó)內(nèi)一些學(xué)者對(duì)此進(jìn)行了較為深入的研究[1-4]。

在堆芯核設(shè)計(jì)中，圍板/水反射層參數(shù)計(jì)算使用TPFAP程序[5]，構(gòu)造計(jì)算幾何模型是固定的，通過(guò)柵元、組件的輸運(yùn)計(jì)算，分別給出圍板及水反射層的兩群宏觀參數(shù)。圍板和水的兩群宏觀參數(shù)直接提供給堆芯程序進(jìn)行堆芯計(jì)算，其中圍板和水是作為兩種材料。反射層參數(shù)對(duì)堆芯計(jì)算尤其是堆芯功率分布影響較大，改進(jìn)反射層參數(shù)計(jì)算方法具有一定的研究意義。

本文介紹了利用DOT程序的輸運(yùn)計(jì)算模塊替換TPFAP程序中的組件兩維穿透幾率計(jì)算模塊，形成程序TPREF的研究過(guò)程，以及對(duì)反射層參數(shù)的影響分析。

1 研究思路

總的研究思路是在TPFAP程序內(nèi)增加接口，在柵元多群輸運(yùn)計(jì)算之后，調(diào)用DOT程序進(jìn)行組件輸運(yùn)計(jì)算并屏蔽TPFAP該部分計(jì)算，計(jì)算輸出通量分布，再繼續(xù)完成原有的TPFAP計(jì)算，輸出圍板和水的兩群宏觀參數(shù)。據(jù)此形成TPREF程序計(jì)算流程如圖1所示。

根據(jù)總體研究思路，DOT程序只需被調(diào)用而不必改動(dòng)，因此程序分析及修改主要針對(duì)TPFAP程序。

2 程序分析

2.1 TPFAP程序分析

TPFAP程序是一個(gè)組件參數(shù)計(jì)算程序。它可以計(jì)算燃料組件、控制組件、可燃毒物組件以及反射層等的少群截面參數(shù)。該程序采用碰撞幾率方法對(duì)各類(lèi)柵元進(jìn)行多群輸運(yùn)計(jì)算，給出組件內(nèi)各類(lèi)柵元的均勻化截面參數(shù)。然后用穿透幾率方法對(duì)組件進(jìn)行寬群二維輸運(yùn)計(jì)算，并群后得到組件的少群常數(shù)。

TPFAP程序的核心是組件二維輸運(yùn)計(jì)算，包括柵元均勻化截面讀取跟蹤、組件二維輸運(yùn)計(jì)算調(diào)用子程序等。

在應(yīng)用穿透幾率法程序TPXY解兩維組件問(wèn)題時(shí)，將組件幾何劃分為NI×NJ個(gè)矩形子區(qū)。通常將一個(gè)柵元視為一個(gè)子區(qū)，且認(rèn)為每個(gè)子區(qū)的核性質(zhì)是均勻的。對(duì)每個(gè)子區(qū)寫(xiě)出分群近似中子積分輸運(yùn)方程為（略去下標(biāo)（i， j））：

（1）

式中：

TPXY的包括四個(gè)子程序，其中SETUP功能主要是初始化源分布、柵元通量分布等，ITERA功能主要是通量計(jì)算及給出組件特征值。FLUXY1的功能是寫(xiě)通量分布。由此可以看出，以上子程序的分析和修改是本研究的核心內(nèi)容。

2.2 DOT程序分析

DOT程序是解二維幾何空間的BOLTZMANN輸運(yùn)方程的程序，主要用于中子或光子輸運(yùn)問(wèn)題。對(duì)反應(yīng)堆的齊次問(wèn)題和外獨(dú)立源問(wèn)題均可求解。求解方法是用離散坐標(biāo)法來(lái)近似求解，即將空間方向（x，y），（r，z），（r，θ）變量離散化，只考慮每個(gè)變量的離散值。能量變量用多群近似，每一群對(duì)應(yīng)一能量間隔，平均出多群截面，在所給的能量范圍內(nèi)考慮所有粒子的運(yùn)動(dòng)及其與介質(zhì)之間的相互作用。粒子運(yùn)動(dòng)方向是通過(guò)選擇一組有序的離散方向及與其有關(guān)的權(quán)重來(lái)處理，粒子沿著離散方向運(yùn)動(dòng)且要滿(mǎn)足平衡條件。

圖1 TPREF程序計(jì)算流程圖

Fig.1 Calculation process diagram of TPREF

該程序可以計(jì)算核反應(yīng)堆的本征值問(wèn)題，如有效增殖系數(shù)keff。也可計(jì)算獨(dú)立外源問(wèn)題，如固定體分布源、固定邊界源及首次解析碰撞源。

其基本方程為積分微分形式的BOLTZMANN輸運(yùn)方程：

（2）

式中：

d■為方向微分立體角元；

dE為能量微分元；

■為位置矢量；

Ψ（■，E，■）為角通量：即定義在■附近的體元dV，■附近的立體角元d■，能量E附近dE內(nèi)的粒子數(shù)為ΨdVdEd■；

σs為從（E'，■'）到（E，■）的散射截面；endprint

σt為總截面；

σf為裂變截面；

v為每次裂變產(chǎn)生的中子數(shù)；

х為裂變譜。

程序中所作簡(jiǎn)化物理假設(shè)為：

（1）方程與時(shí)間無(wú)關(guān)；

（2）σt與■無(wú)關(guān)；

（3）裂變?cè)粗慌c能群有關(guān)而與方向無(wú)關(guān)。

由于本研究不涉及DOT程序的修改，因此對(duì)DOT程序的理論模型及程序結(jié)構(gòu)不再詳細(xì)介紹。

3 程序開(kāi)發(fā)

根據(jù)前述TPXY有關(guān)子程序的邏輯順序，首先在SETUP子程序中增加DOT輸運(yùn)計(jì)算所需初始化的內(nèi)容（包括計(jì)算區(qū)域幾何信息）；其次是增加輸運(yùn)計(jì)算子程序SNDOT，并在ITERA子程序中進(jìn)行調(diào)用，并寫(xiě)入通量分布。

為了形成便于堆芯計(jì)算使用的圍板和水的宏觀參數(shù)，對(duì)后處理子程序INPUT和EDITGC作了少量的修改，目的是對(duì)參數(shù)進(jìn)行重新排序。

在原TPFAP程序基礎(chǔ)上進(jìn)行開(kāi)發(fā)修改之后，重新編譯獲得新的計(jì)算程序，與原版程序的主要區(qū)別是修改了兩維組件輸運(yùn)計(jì)算。將新的計(jì)算程序命名為T(mén)PREF。

4 程序檢驗(yàn)

利用TPREF計(jì)算堆芯反射層參數(shù)，見(jiàn)表1。

表1 反射層參數(shù)計(jì)算結(jié)果比較

Table1 Comparison of reflector parameters calculation

將以上反射層參數(shù)用于堆芯計(jì)算，功率分布計(jì)算結(jié)果基本相同，通量分布計(jì)算結(jié)果差別較小。

分析其原因，主要是由于維持了原有反射層計(jì)算幾何模型不變，然后利用DOT輸運(yùn)計(jì)算模塊替代了穿透幾率計(jì)算模塊，對(duì)程序計(jì)算幾何模型沒(méi)有影響。因此計(jì)算結(jié)果相符合，這表明程序修改是正確的，實(shí)現(xiàn)了兩維組件輸運(yùn)計(jì)算模塊的更換。

5 結(jié)論

通過(guò)本研究結(jié)果表明，將二維離散坐標(biāo)程序DOT程序的輸運(yùn)計(jì)算模塊替代TPFAP程序中的穿透幾率輸運(yùn)計(jì)算模塊是可行的，程序的開(kāi)發(fā)和編制是正確的。作為兩種不同的輸運(yùn)計(jì)算方法，穿透幾率法和離散坐標(biāo)法在二維組件輸運(yùn)計(jì)算中，其計(jì)算精度是相當(dāng)?shù)摹?/p>

【參考文獻(xiàn)】

[1]成紅武，謝仲生，章宗耀.圍板和反射層的等效均勻化常數(shù)的計(jì)算[J].核科學(xué)與工程，1996，16（2）：97-103.

[2]徐曉勤.壓水堆圍板對(duì)堆芯物理特性的影響[J].核動(dòng)力工程，1999，20（3）：205-208.

[3]司勝義.壓水堆堆芯反射層的截面參數(shù)化及優(yōu)化設(shè)計(jì)[J].核電工程與技術(shù)，2004（2）.

[4]吳宏春，劉海峰.反射層參數(shù)的輸運(yùn)等效計(jì)算程序[J].核動(dòng)力工程，2005，26（3）：224-227.

[5]章宗耀，李大圖，姚棟，等.輕水堆燃料組件計(jì)算程序包TPFAP[J].核動(dòng)力工程，1993，14（2）：117-121，192.endprint