[發(fā)明專利]基于拋物線插值法的燃料燃耗確定方法及裝置

在審

• 申請(qǐng)?zhí)枺?/td>: 202011209168.5
• 申請(qǐng)日: 2020-11-03
• 公開(kāi)（公告）號(hào): CN112380680A
• 公開(kāi)（公告）日: 2021-02-19
• 發(fā)明（設(shè)計(jì)）人: 郭鳳晨;蘆韡;柴曉明;李慶;劉東;安萍;尹強(qiáng);盧宗健;曾輝;涂曉蘭;馬永強(qiáng);鄭勇;賀濤
• 申請(qǐng)（專利權(quán)）人: 中國(guó)核動(dòng)力研究設(shè)計(jì)院
• 主分類號(hào): G06F30/20
• 分類號(hào): G06F30/20;G21C17/104;G06F119/06
• 代理公司: 核工業(yè)專利中心 11007
• 代理人: 王朋

本公開(kāi)屬于核電技術(shù)領(lǐng)域，具體涉及一種基于拋物線插值法的燃料燃耗確定方法及裝置。本公開(kāi)根據(jù)確定的所述每個(gè)燃耗步末時(shí)刻的各核素的微觀反應(yīng)率的預(yù)測(cè)值和核子密度的第二預(yù)測(cè)值，采用拋物線插值法，確定各核素在當(dāng)前燃耗步末時(shí)刻的微觀反應(yīng)率的修正值；根據(jù)各核素在當(dāng)前燃耗步末時(shí)刻的微觀反應(yīng)率的修正值，各核素在起始時(shí)刻的核子密度，采用燃耗方程，確定各核素在當(dāng)前燃耗步末時(shí)刻的核子密度的修正值。本公開(kāi)通過(guò)使用拋物線插值法，對(duì)燃耗計(jì)算修正步所使用核素的微觀反應(yīng)率進(jìn)行修正，可以減少線性插值修正帶來(lái)的計(jì)算誤差，從而提高計(jì)算精度。

技術(shù)領(lǐng)域

本發(fā)明屬于核電技術(shù)領(lǐng)域，具體涉及一種基于拋物線插值法的燃料燃耗確定方法及裝置。

背景技術(shù)

通常來(lái)說(shuō)，堆芯核設(shè)計(jì)程序系統(tǒng)包括組件中子學(xué)計(jì)算程序和堆芯中子學(xué)計(jì)算程序。無(wú)論組件計(jì)算程序還是堆芯計(jì)算程序，都需要在中子輸運(yùn)方程或中子擴(kuò)散方程的基礎(chǔ)上進(jìn)行燃耗計(jì)算。燃耗計(jì)算是反應(yīng)堆物理計(jì)算分析的重要內(nèi)容，在計(jì)算反應(yīng)堆內(nèi)燃料的增殖與消耗、反應(yīng)性的變化等方面起著至關(guān)重要的作用。但是，相關(guān)技術(shù)中，針對(duì)組件的核子密度的計(jì)算常常會(huì)出現(xiàn)核子密度數(shù)值偏大的現(xiàn)象，嚴(yán)重影響燃耗計(jì)算的準(zhǔn)確性，因此如何準(zhǔn)確的計(jì)算含核素的組件的核子密度成為亟待解決的問(wèn)題。

發(fā)明內(nèi)容

為克服相關(guān)技術(shù)中存在的問(wèn)題，提供了一種基于拋物線插值法的燃料燃耗確定方法及裝置。

根據(jù)本公開(kāi)實(shí)施例的一方面，提供一種基于拋物線插值法的燃料燃耗確定方法，所述方法包括：

針對(duì)多個(gè)核素當(dāng)前燃耗步以及緊鄰當(dāng)前燃耗步的前兩個(gè)燃耗步中的每個(gè)燃耗步進(jìn)行以下操作：

根據(jù)各核素在該燃耗步初時(shí)刻的微觀反應(yīng)率的預(yù)測(cè)值，以及各核素在該燃耗步起始時(shí)刻的核子密度，采用燃耗方程，確定各核素在該燃耗步末時(shí)刻的核子密度的第一預(yù)測(cè)值；

根據(jù)各核素在該燃耗步末時(shí)刻的核子密度的第一預(yù)測(cè)值，采用中子輸運(yùn)方程，確定各核素在該燃耗步末時(shí)刻的微觀反應(yīng)率的預(yù)測(cè)值；

根據(jù)各核素在該燃耗步末時(shí)刻的微觀反應(yīng)率的預(yù)測(cè)值，以及各核素在該燃耗步起始時(shí)刻的核子密度，采用燃耗方程，確定各核素在該燃耗步末時(shí)刻的核子密度的第二預(yù)測(cè)值；

根據(jù)確定的所述每個(gè)燃耗步末時(shí)刻的各核素的微觀反應(yīng)率的預(yù)測(cè)值和核子密度的第二預(yù)測(cè)值，采用拋物線插值法，確定各核素在當(dāng)前燃耗步末時(shí)刻的微觀反應(yīng)率的修正值；

根據(jù)各核素在當(dāng)前燃耗步末時(shí)刻的微觀反應(yīng)率的修正值，各核素在起始時(shí)刻的核子密度，采用燃耗方程，確定各核素在當(dāng)前燃耗步末時(shí)刻的核子密度的修正值；

根據(jù)各核素在當(dāng)前燃耗步末時(shí)刻的核子密度的修正值，以及各核素在當(dāng)前燃耗步末時(shí)刻的核子密度的第一預(yù)測(cè)值，確定各核素在當(dāng)前燃耗步末時(shí)刻的核子密度的最終值。

在一種可能的實(shí)現(xiàn)方式中，根據(jù)確定的所述每個(gè)燃耗步末時(shí)刻的各核素的微觀反應(yīng)率的預(yù)測(cè)值和核子密度的第二預(yù)測(cè)值，采用拋物線插值法，確定各核素在當(dāng)前燃耗步末時(shí)刻的微觀反應(yīng)率的修正值，包括：

針對(duì)每個(gè)核素，進(jìn)行以下操作：

將式一作為該核素的拋物線插值法公式；

通過(guò)式二確定a_i、b_i和c_i的值；

根據(jù)確定的a_i、b_i、c_i，N_i,n，以及式一，確定各核素在當(dāng)前燃耗步末時(shí)刻的微觀反應(yīng)率的修正值；