發(fā)布時間：2021-08-18 02:16

　　在核反應(yīng)堆設(shè)計中通常采用有限差分方法對擴散方程求解,但由于其邊界條件反照率的選取帶有很多經(jīng)驗因子,無法得到精確的解析解。針對以上問題,采用節(jié)塊法的一種邊界條件,利用Matlab編程數(shù)值求解反應(yīng)堆堆芯二群中子擴散方程。分析了1/4反應(yīng)堆堆芯內(nèi)數(shù)值計算中子通量密度的劃分節(jié)塊依據(jù),綜合邊界條件研究了堆芯不同位置處節(jié)塊的具體差分方程形式以及相關(guān)系數(shù)計算表達(dá)式,討論了擴散方程組系數(shù)矩陣的具體構(gòu)造形式,通過編程構(gòu)造方程組系數(shù)矩陣結(jié)合SOR迭代法以及乘冪法進(jìn)行內(nèi)外迭代對堆芯中子通量密度進(jìn)行數(shù)值模擬。分析表明,計算結(jié)果基本合理,程序收斂迅速。 

【文章來源】：兵器裝備工程學(xué)報. 2020,41(02)北大核心

【文章頁數(shù)】：6 頁

【部分圖文】：

第4次計算Keff收斂過程

表1中界面AB等按照圖1對應(yīng)，其中，BCA,BCB分別為利用第一種和第二種方法得到的反照率數(shù)值解。在編制計算程序的過程中，對反照率邊界條件進(jìn)行進(jìn)一步簡化，綜合考慮補充邊界節(jié)點代數(shù)方程的方法和附加源項法[9]，假設(shè)在反射層中與堆芯緊鄰處添加節(jié)點，且該節(jié)塊內(nèi)的系數(shù)與堆芯相同，同時假設(shè)其中子通量密度Φr與該緊鄰堆芯節(jié)塊內(nèi)中子通量密度Φr關(guān)系為:

式(3)為普遍情況下節(jié)塊的差分方程。對于元胞1，即1～8號組件第1層節(jié)塊，上表面需要考慮全反射邊界條件，第150個節(jié)塊即(1,150)號節(jié)塊，則還要考慮反照率邊界條件。式(3)中各項系數(shù)如下所示:

【參考文獻(xiàn)】：
期刊論文
[1]AP1000核反應(yīng)堆控制棒價值特性的MC模擬[J]. 謝明亮,于雷,陳玉清.  兵器裝備工程學(xué)報. 2016(03)

博士論文
[1]中子輸運方程特征線解法及嵌入式組件均勻化方法的研究[D]. 湯春桃.上海交通大學(xué) 2009



本文編號：3348990