本發(fā)明涉及核反應堆堆芯設計和核反應堆物理計算領(lǐng)域，是一種獲取快中子反應堆少群截面參數(shù)的混合計算方法。

背景技術(shù)：

為了快速、精確地求解堆芯中子學參數(shù)，基于“兩步法”的方法成為快堆工程計算的主要方法。所謂“兩步法”，第一步是對堆芯內(nèi)各種組件材料進行建模并計算，獲得組件內(nèi)的中子通量密度分布從而歸并出少群均勻化群截面；第二步是由前一步計算得到的均勻化參數(shù)對堆芯進行少群中子輸運方程的求解，獲取堆芯有效增殖因子、堆芯功率分布等物理量。

計算少群截面采用的方法可分為蒙特卡洛方法以及確定論方法。蒙特卡洛方法是一種基于統(tǒng)計學的方法，它通過抽樣大量粒子進行中子學特性的模擬，從而獲得中子學相關(guān)的參數(shù)。確定論方法是通過各種理論或數(shù)值方法求解描述中子學特性的數(shù)學方程組，通過獲取方程組的解獲得中子學相關(guān)的參數(shù)。

由于蒙特卡洛方法可以使用連續(xù)能量的數(shù)據(jù)庫，同時對幾何模型可以精確建模，因此其計算結(jié)果具有極高的精度。在蒙特卡洛計算中，少群截面的計算需要對核反應率、中子通量密度分布等信息進行統(tǒng)計。但當這些數(shù)值很小時，會難以獲得準確的結(jié)果，例如高階中子通量矩。由于在高階散射截面的計算中需要用到高階中子通量矩，而高階中子通量矩又難以獲得，因此在蒙特卡洛方法中引入了近似，它使用零階中子通量矩對高階散射截面進行計算。而在確定論方法中，可以利用數(shù)值方法計算出準確的高階中子通量矩，避免了數(shù)值很小時計算結(jié)果不準確的問題。但在確定論方法中，少群截面計算結(jié)果的精確度取決于多群截面的計算結(jié)果，而多群截面的計算又與共振計算方法相關(guān)，因此使用確定論方法時須基于精度較高的共振計算方法。

由于單獨使用蒙特卡洛方法或確定論方法均有難以克服的缺點，因此有必要將這兩種方法進行結(jié)合，發(fā)明一種高精度的快堆少群截面的混合計算方法。

技術(shù)實現(xiàn)要素：

為了克服快堆少群截面計算中蒙特卡洛方法以及確定論方法的缺點，本發(fā)明在計算快堆少群截面的流程中同時使用蒙特卡洛方法和確定論方法，利用蒙特卡洛方法產(chǎn)生核素的微觀總截面、裂變截面、彈性散射截面，利用確定論方法計算其他截面信息以及中子通量密度分布，以此獲得高精度的快堆組件少群截面。

為了實現(xiàn)以上目的，本發(fā)明采取如下的技術(shù)方案予以實施：

一種獲取快中子反應堆少群截面參數(shù)的混合計算方法，包括如下步驟：

步驟1：針對所需計算的任一快堆組件，讀取該組件的幾何信息以及相應材料組分的核素信息；

步驟2：針對步驟1讀取的幾何信息及材料組分信息，讀取所有材料中各個核素的微觀非彈性散射截面值、每次裂變釋放中子數(shù)、裂變譜信息；

步驟3：根據(jù)步驟1讀取的信息，通過指定蒙特卡洛計算的輸入信息建立相應的計算模型并進行計算；在計算過程中統(tǒng)計各個核素在各個區(qū)域中的多群總反應率、裂變反應率、彈性散射反應率以及中子通量密度分布；

步驟4：利用步驟3獲得的各核素的多群總反應率、裂變反應率、彈性散射反應率以及中子通量密度分布，利用公式(1)求得各核素的多群微觀總截面、微觀裂變截面以及微觀彈性散射截面；

式中：

σx,g——第g群的微觀截面，下標x指代微觀總截面、裂變截面或彈性散射截面

rx,g——第g群的微觀反應率，下標x指代微觀總反應率、裂變反應率或彈性散射反應率

φg——第g群的中子通量密度分布

步驟5：:利用公式(2)計算每個核素群到群的散射概率，再根據(jù)步驟4計算得到的各核素的多群微觀彈性散射截面并利用公式(3)計算每個核素的彈性散射矩陣；

式中：

f^l(g→g')——第l階由第g群散射到第g’群的散射概率

α——(a‐1)²/(a+1)²，a為核素的原子質(zhì)量

μc——質(zhì)心坐標系下的散射角余弦值

f(e,μc)——散射概率函數(shù)

pl(μs)——關(guān)于μs的l階勒讓德多項式

μs——實驗室坐標系下的散射角余弦值

δeg——第g群的能群間隔

式中：

——第l階由第g群散射到第g’群的彈性散射截面

σs,g——第g群的彈性散射截面

步驟6：根據(jù)由步驟2至步驟5獲得的各個核素的多群微觀截面信息，結(jié)合步驟1讀取的組件幾何信息及材料信息，進行基于確定論方法的中子輸運方程的求解，以此獲得多群的各階中子通量矩分布；

步驟7：基于步驟2至步驟5獲得的多群微觀截面信息以及步驟6求解得到的多群的各階中子通量矩分布，并利用公式(4)、(5)對多群微觀截面進行能群、空間的歸并，從而得到組件的少群截面；

式中：

σx,g,i——第g群、屬于第i區(qū)的微觀截面，下標x指代微觀總截面、微觀裂變截面

σx,g,i——第g群、屬于第i區(qū)的微觀截面，下標x指代微觀總截面、微觀裂變截面

φg,i——第g群、屬于第i區(qū)的中子通量密度分布

vi——第i子區(qū)域的體積

式中：

——第l階由第g群散射到第g’群、屬于第i區(qū)的微觀散射截面

——第l階由第g群散射到第g’群、屬于第i區(qū)的微觀散射截面

——第g群、屬于第i區(qū)的第l階的中子通量密度分布。

與現(xiàn)有技術(shù)相比，本發(fā)明有如下突出的優(yōu)點：

1.在混合計算方法中，多群微觀總截面、裂變截面、彈性散射截面信息由蒙特卡洛方法進行計算，保證了高精度；

2.在混合計算方法中，各階的彈性散射矩陣由確定論方法進行計算，避免了蒙特卡洛方法中的近似，保證了散射截面精度；

3.在混合計算方法中，各階的中子通量矩由確定論方法進行計算，避免了蒙特卡洛方法中的近似；同時利用高階中子通量矩對高階散射截面進行能群及區(qū)域的歸并，獲得更加準確的少群截面。

附圖說明

圖1為快堆少群截面計算的混合方法流程圖。

圖2為均勻問題各個核素24群微觀截面的誤差。

圖3為一維圓柱問題24群宏觀總截面的誤差。

具體實施方式

下面結(jié)合附圖和具體實施方式對本發(fā)明做進一步詳細說明：

本發(fā)明通過將蒙特卡洛方法與確定論方法相結(jié)合，利用蒙特卡洛方法精細考慮快堆組件的共振效應，計算各個核素精確的多群微觀總截面、裂變截面、彈性散射截面，同時利用確定論方法求解快堆組件的各階彈性散射截面以及各階中子通量矩，以此將多群截面歸并為少群截面，如圖1所示，該發(fā)明包括如下步驟：

步驟1：針對所需計算的任一快堆組件，讀取該組件的幾何信息以及相應材料組分的核素信息；

步驟2：針對步驟1讀取的幾何信息及材料組分信息，讀取所有材料中各個核素的微觀非彈性散射截面值、每次裂變釋放中子數(shù)、裂變譜信息；

步驟3：根據(jù)步驟1讀取的信息，通過指定蒙特卡洛計算的輸入信息建立相應的計算模型并進行計算；在計算過程中統(tǒng)計各個核素在各個區(qū)域中的多群總反應率、裂變反應率、彈性散射反應率以及中子通量密度分布；

步驟4：利用步驟3獲得的各核素的多群總反應率、裂變反應率、彈性散射反應率以及中子通量密度分布，利用公式(1)求得各核素的多群微觀總截面、微觀裂變截面以及微觀彈性散射截面；

式中：

σx,g——第g群的微觀截面，下標x指代微觀總截面、微觀裂變截面或微觀彈性散射截面

rx,g——第g群的微觀反應率，下標x指代微觀總反應率、微觀裂變反應率或微觀彈性散射反應率

φg——第g群的中子通量密度分布

步驟5：:利用公式(2)計算每個核素群到群的散射概率，再根據(jù)步驟4計算得到的各核素的多群微觀彈性散射截面并利用公式(3)計算每個核素的彈性散射矩陣；

式中：

f^l(g→g')——第l階由第g群散射到第g’群的散射概率

α——(a‐1)²/(a+1)²，a為核素的原子質(zhì)量

μc——質(zhì)心坐標系下的散射角余弦值

f(e,μc)——散射概率函數(shù)

pl(μs)——關(guān)于μs的l階勒讓德多項式

μs——實驗室坐標系下的散射角余弦值

δeg——第g群的能群間隔

式中：

——第l階由第g群散射到第g’群的彈性散射截面

σs,g——第g群的彈性散射截面

步驟6：根據(jù)由步驟2至步驟5獲得的各個核素的多群微觀截面信息，結(jié)合步驟1讀取的組件幾何信息及材料信息，進行基于確定論方法的中子輸運方程的求解，以此獲得多群的各階中子通量矩分布；

步驟7：基于步驟2至步驟5獲得的多群微觀截面信息以及步驟6求解得到的多群的各階中子通量矩分布，并利用公式(4)、(5)對多群微觀截面進行能群、空間的歸并，從而得到組件的少群截面。

式中：

σx,g,i——第g群、屬于第i區(qū)的微觀截面，下標x指代微觀總截面、微觀裂變截面

σx,g,i——第g群、屬于第i區(qū)的微觀截面，下標x指代微觀總截面、微觀裂變截面

φg,i——第g群、屬于第i區(qū)的中子通量密度分布

vi——第i子區(qū)域的體積

式中：

——第l階由第g群散射到第g’群、屬于第i區(qū)的微觀散射截面

——第l階由第g群散射到第g’群、屬于第i區(qū)的微觀散射截面

——第g群、屬于第i區(qū)的第l階的中子通量密度分布

在本發(fā)明中，任一快堆組件的幾何及材料信息由步驟1讀入，根據(jù)步驟1讀入的信息須通過步驟2在多群數(shù)據(jù)庫中讀取各個核素的微觀非彈性散射截面值、每次裂變釋放中子數(shù)、裂變譜信息。適用于反應堆物理計算的多群數(shù)據(jù)庫均包含上述信息，本發(fā)明對多群數(shù)據(jù)庫的選擇沒有限制。

步驟3的蒙特卡洛計算中，須對統(tǒng)計區(qū)域做出定義，即設置統(tǒng)計區(qū)域的空間所包含的范圍、能量區(qū)間所處的上下限以及具體統(tǒng)計核素的名稱信息，統(tǒng)計時可采用任意的蒙特卡洛統(tǒng)計方法，如碰撞計數(shù)法或徑跡長度計數(shù)法，本發(fā)明對反應率及通量的統(tǒng)計不受統(tǒng)計方法的限制。

步驟6中各階中子通量矩分布的求解可采用不同空間維度的中子輸運方程求解方法，如碰撞概率方法、離散縱標方法、特征線方法，空間維度包含均勻問題、一維問題。本發(fā)明對于各階中子通量矩分布以及少群截面歸并的求解過程不受中子輸運方程求解方法的限制。

為驗證本發(fā)明的有效性，圖2展示了利用本發(fā)明計算得到的一均勻快堆組件問題各個核素的24微觀總截面與參考解的相對誤差。計算結(jié)果表明，利用本發(fā)明計算得到的少群截面與參考解相比，少群截面具有極高的精度，各個核素各個能群的截面誤差均在1％以內(nèi)。為驗證本發(fā)明計算真實組件的能力，針對一維圓柱幾何快堆組件，分別采用蒙特卡洛方法、確定論方法以及混合計算方法進行計算，以蒙特卡洛方法為參考解，比較其它兩種方法的24群少群宏觀總截面的誤差。如圖3所示，相比于確定論方法，混合計算方法所產(chǎn)生的截面具有更高的精度，截面誤差均在1％以內(nèi)。

利用準確的少群截面，就可以進行快堆堆芯的各種中子學計算，如穩(wěn)態(tài)計算、燃耗計算、瞬態(tài)計算等。本發(fā)明可以得到較高精度的快堆組件少群截面，可應用于實際的工程計算當中。