CN105426342A - 一种获取反应堆物理栅格计算重要参数灵敏度系数的方法 - Google Patents
一种获取反应堆物理栅格计算重要参数灵敏度系数的方法 Download PDFInfo
- Publication number
- CN105426342A CN105426342A CN201510726692.2A CN201510726692A CN105426342A CN 105426342 A CN105426342 A CN 105426342A CN 201510726692 A CN201510726692 A CN 201510726692A CN 105426342 A CN105426342 A CN 105426342A
- Authority
- CN
- China
- Prior art keywords
- phi
- mrow
- formula
- sigma
- subgroup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
- Particle Accelerators (AREA)
Abstract
一种获取反应堆物理栅格计算重要参数灵敏度系数的方法,1、针对具体的反应堆物理计算问题建立计算模型;2、采用子群共振自屏方法进行共振自屏计算,针对快群源项和共振能群源项,以及共振自屏截面,建立对应的广义共轭方程,对广义共轭方程进行求解,得到共振自屏截面灵敏度系数;3、针对需要计算的重要参数,建立并求解与之对应的广义共轭方程,得到该参数的显式灵敏度系数;4、结合共振自屏截面的灵敏度系数和共振自屏截面的显式灵敏度系数,得到隐式灵敏度系数,总的灵敏度系数等于显式灵敏度系数与隐式灵敏度系数之和;能够快速地得到任意能表示成通量或共轭通量泛函的参数的灵敏度系数,同时考虑了隐式影响,使得灵敏度系数更加精确。
Description
技术领域
本发明涉及核反应堆堆芯设计和安全领域,具体涉及一种获取反应堆物理栅格计算重要参数灵敏度系数的方法。
背景技术
传统的核反应堆设计和运行往往采用保守计算,或者设置足够大的安全裕量来保证反应堆的安全。这种做法虽然一定程度上能保证反应堆的安全性,但是丧失了极大的经济性。近年来,国际上在核研究、核工业、核安全和核监管等领域对于提供核反应堆安全参数的“最佳估计值+不确定性”的需求日益增长。
在反应堆物理栅格计算针对核数据的敏感性和不确定性分析中,灵敏度系数计算是获取计算结果不确定度的重要环节,因此灵敏度系数的计算就显得特别重要。根据栅格计算先进行共振计算,再进行多群中子输运计算的特点,灵敏度系数可被分为隐式灵敏度系数和显式灵敏度系数两部分。前者是核数据通过在共振计算中对共振截面的影响,再体现到对中子输运方程计算结果的一种间接影响;后者是核数据通过中子输运方程求解对计算结果的一种直接影响。然而,目前很多敏感性和不确定性分析程序中,作为间接影响的隐式敏感性往往被忽略,而只考虑了作为直接影响的显式敏感性,这势必会引入一定误差。此外,目前隐式敏感性分析往往基于较为简单的共振自屏方法,而针对当前比较主流的共振计算方法——子群方法的隐式敏感性分析并不广泛和深入。另一方面,目前有效增值因子的敏感性和不确定性受到广泛的重视,但对其它重要参数(例如少群常数、功率分布等)的分析计算不够充分。如果要将栅格计算结果的不确定性传递到堆芯计算,诸如少群常数等重要参数的不确定度计算就必不可少。
因此有必要发明一种切实可行的方法,对反应堆物理栅格计算中各种重要参数进行敏感性和不确定性分析,并能够针对目前广泛应用的子群共振自屏方法进行隐式敏感性分析,从而使得敏感性和不确定性的计算更加精确。
发明内容
为了获取反应堆物理栅格计算中重要参数对核数据的灵敏度系数,同时考虑隐式敏感性,本发明提出了一种获取反应堆物理栅格计算重要参数灵敏度系数的方法,是能够精确考虑反应堆物理栅格计算中的灵敏度系数的计算方法。
为了实现上述目的,本发明采用以下技术方案予以实施:
一种获取反应堆物理栅格计算重要参数灵敏度系数的方法,该方法包括以下步骤:
步骤1:基于微扰理论,将多输入参数的反应堆物理计算需要大量扰动计算才能得到灵敏度系数的问题,转化为只需一次前向计算和一次共轭计算就能获取某一计算结果对所有输入参数的灵敏度系数的问题;
设反应堆物理计算中某重要参数为R,以下称为响应,表示成通量或者共轭通量的线性泛函,设为:
式中:
Φ——中子角通量密度;
Φ*——共轭中子角通量密度;
ξ——问题求解的空间;
Σ(ξ)——同求解空间相关的截面数据;
H1[Σ(ξ)],H2[Σ(ξ)]——依赖于截面数据的算子;
建立与响应R的表达式对应的广义共轭输运方程:
式中:
M——输运算子;
Γ——广义通量密度;
M*——输运算子的共轭算子;
Γ*——广义共轭通量密度;
步骤2:采用广泛应用的输运求解方法——组件模块化特征线方法MOC作为二维输运方程的求解方法,对中子输运方程、共轭中子输运方程,以及广义共轭输运方程进行求解;
特征线方法能够实现一般中子输运方程的求解,但是对于共轭中子输运方程以及广义共轭输运方程,需要对求解流程作一点修改;
对于共轭中子输运方程计算,在计算之前,需要如下操作:
(1)将散射矩阵转置;
(2)将材料的裂变产生截面同裂变谱χ向量互换;
(3)将所有截面的能群号按如下方式变化:
完成上述操作之后,输运求解器能够完成基本的共轭中子输运方程的求解;
当求解广义共轭输运方程,还需要如下操作:
(1)将建立共轭源作为问题的外源项;
(2)对于公式(10)将裂变源的更新按如下形式进行:
式中:
FSg——第g群裂变源;
——裂变产生截面;
keff——有效增殖因子;
χ——裂变谱;
Γ*——广义共轭通量密度;
F*——中子输运方程裂变算子的共轭算子;
Φ——前向中子通量密度;
Φ*——共轭中子通量密度;
其计算流程与传统的中子输运方程计算流程基本一致,不同点在于:(1)不需要进行特征值的更新,方程中的特征值使用的是中子输运方程求解得到的特征值;(2)每一次外迭代得到的通量用于裂变源更新时采用公式(18);
对于方程(15)将裂变源的更新按如下形式进行:
式中:
Γ——广义通量密度;
F——中子输运方程裂变算子;
当采用输运求解器完成共轭中子输运方程或者广义共轭输运方程的求解之后,需要将获得的角通量在角度上反转,在能群上颠倒,从而为下一步计算提供正确的通量信息;
步骤3:在步骤1建立的方法的基础上,该步针对子群共振计算方法,建立相应的子群广义共轭方程;采用步骤2建立的求解方法,选取组件模块化特征线方法MOC作为二维输运求解方法,求解子群广义共轭方程,为共振自屏截面的灵敏度系数计算提供广义共轭通量;
子群方法中,共振自屏截面的表达式为:
式中:
g——能群标号;
i——子群标号;
N——子群总数;
σx(E)——x反应的连续能量截面;
φ(E)——权重谱;
σx,g,i——共振截面x的第g共振能群的第i个子群截面;
公式(31)中的子群参数通过帕德近似求解,而子群通量通过求解子群输运方程得到;
设子群输运方程为:
Ω·▽φg,i(r,Ω)+Σt,g,i(r)φg,i(r,Ω)=Qs,g,i(r,Ω)公式(32)
式中:
g——共振能群标号;
i——子群标号;
r——空间位置变量;
Ω——角度变量;
φg,i(r,Ω)——能群g的第i子群的中子通量密度;
Σt,g,i(r)——能群g的第i子群的宏观总截面;
Qs,g,i(r,Ω)——散射源项;
将上述方程写成算子形式为:
Lgφg=Qg公式(33)
式中:
Lg——第g个共振能群的子群输运方程的输运算子;
Qg——第g个共振能群的子群输运方程的源项;
利用输运求解器求解子群输运方程,得到子群通量后,对子群参数加权得到共振自屏截面;
子群共振方法的共振截面灵敏度系数计算公式为:
其中为广义子群共轭方程的解。它由该共振自屏截面的子群输运方程对应的广义子群共轭方程解得,该广义子群共轭方程为:
式中:
——Lg的共轭算子;
——子群广义共轭通量;
——子群广义共轭源;
其中
式中
——第i个子群的源项;
σx,g,i——第i个子群截面;
注意到式(34)中存在子群参数灵敏度系数,考虑采用帕德近似方法求解子群参数的过程中,采用直接扰动方法求解子群参数灵敏度系数;对于各共振核素,逐群扰动其连续能量截面,体现为扰动其共振积分表,设扰动百分比为δ,则根据差商代替微分的方法求得子群参数的灵敏度系数,即
公式(37)
式中为未扰动的子群参数即子群截面σx,g,i或子群概率pg,i,和分别为正向和负向扰动α时的子群参数;δ为α的扰动百分比;
此外,式(34)中存在源项Qg对α的灵敏度系数,这一项的求解同样采用步骤1建立的方法;采用的子群共振自屏计算方法中,源项分成快群散射源项Qf,g和上游共振能群散射源项Qr,g,分别为:
对于共振能群g,分别对上述两种源建立广义共轭方程:
式中:
——快群输运算子;
——广义共轭通量;
g’——快群标号;
Σs,g'→g——快群到共振能群的散射截面;
而对于上游共振能群g’,有
式中:
——子群输运算子的共轭算子;
——广义共轭通量;
g’——上游共振能群标号;
i——第g’群的子群标号;
Σs,g'→g,i——上游共振能群g’到当前共振能群g的第i个子群的散射截面;
采用输运求解器求解方程式(40)和(41),得到对应的广义共轭通量;由下式求得源项的灵敏度系数
式中:
Qg——快群散射源项或上游共振能群散射源项;
Σs——快群或上游共振能群到当前共振能群的散射截面;
Q'——求解源项的输运方程的右端源项;
步骤4:针对反应堆物理栅格计算的重要参数,根据步骤1建立的方法获取与之对应的广义共轭方程,利用步骤2的输运求解器,求解获取通量、共轭通量以及广义通量;从而得到各重要参数的显式灵敏度系数;
设某重要参数为R,表示成通量或者共轭通量的线性泛函,设为:
注意这里R具有一般性,只要是能表示成上述形式的重要参数,都能够进行灵敏度系数的求解;
建立对应的广义共轭源项:
利用输运求解器分别求解:
M*Γ*=Q1公式(46)
MΓ=Q2公式(47)
获取通量、广义共轭通量以及广义通量后,显式灵敏度系数由下式计算得到:
步骤5:根据步骤3得到的共振自屏截面的灵敏度系数和步骤4得到的显式灵敏度系数,计算出考虑了隐式影响的灵敏度系数;
其中j表示共振核素。
与现有技术相比,本发明具有如下优点:
1.由于本发明能够考虑核素截面扰动在子群共振自屏计算过程中对共振自屏截面的影响,因此可以得到核截面数据的考虑了隐式影响的灵敏度系数,使得灵敏度系数的计算更加精确。
2.基于微扰理论,能够建立任意中子通量或共轭通量的泛函的灵敏度系数计算方法,因此可以对反应堆物理计算中多种重要参数进行敏感性计算分析,此外由于微扰理论本身能够快速计算多输入扰动影响的优点,因此本发明可以对存在大量输入参数的反应堆物理计算问题进行快速的敏感性计算,而不要对大量的输入参数进行逐一的扰动,从而提高了计算效率。
附图说明
图1是灵敏度系数计算流程图。
图2是特征线方法求解中子输运方程流程图。
图3是特征线方法求解广义共轭方程流程图。
图4是典型压水堆栅元显式和隐式灵敏度系数对比图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
如图1所示,本发明是一种获取反应堆物理栅格计算重要参数灵敏度系数的方法,包括如下步骤:
步骤1:基于微扰理论,将多输入参数的反应堆物理计算需要大量扰动计算才能得到灵敏度系数的问题,转化为只需一次前向计算和一次共轭计算就能获取某一计算结果对所有输入参数的灵敏度系数的问题。
设中子输运方程为如下算子形式:
MΦ=(L-λF)Φ=0公式(1)
式中:
M——输运算子;
F——裂变源项算子;
L——输运算子中除裂变源项算子以外的其它算子;
Φ——中子角通量密度;
λ——方程的特征值(λ=1/keff),keff为***有效增值因子。
设中子输运方程的共轭方程为:
M*Φ*=(L*-λF*)Φ*=0公式(2)
式中:
M*——输运算子的共轭算子;
F*——裂变源项算子的共轭算子;
L*——输运算子中除裂变源项算子以外的其它算子的共轭算子;
Φ*——共轭中子角通量密度。
设反应堆物理计算中某重要参数为R(以下称为响应),可以表示成通量或者共轭通量的线性泛函,设为:
式中:
ξ——问题求解的空间;
Σ(ξ)——同求解空间相关的截面数据;
H1[Σ(ξ)],H2[Σ(ξ)]——依赖于截面数据的算子。
相对灵敏度系数表示当α相对扰动1%时,R相对改变的百分比,其定义为:
式中:
α——多群微观截面;
SR,α——响应R对α的的灵敏度系数。
则由灵敏度系数定义,有:
式中<>表示在问题求解空间ξ上积分。
整理上式得到:
上式中通量对α的倒数项可以通过微理论建立方法来进行求解。对于中子输运方程两边对α求导:
于是
取与对应的广义共轭函数并作积分:
令
根据式(8)和式(9)有
所以只要广义共轭函数求出来,就能计算出上式左端项,也就得到了式(6)右端第二项。
对于式(6)右端第三项同理,对共轭中子输运方程两边对α求导:
移项有:
对于取其广义共轭通量并积分:
令
根据式(13)和式(14),有:
所以,只要得到了广义共轭通量,就能得到上式左端项,也就得到了式(6)的右端第三项。
从而,可以将灵敏度系数的计算转化为:
根据以上基于微扰理论的推导,对于给定的可以表示成通量或者共轭通量的泛函的任意响应,不需要对输入参数一个一个进行扰动来计算其灵敏度系数,而只需对其建立与之对应的广义(共轭)方程,求解这样带有外源和裂变源的中子输运方程,得到广义(共轭)通量,就能计算出这种响应对所有输入参数的灵敏度系数。与直接扰动方法相比,在响应数量较少时,具有很高的计算效率。
步骤2:由步骤1的方法所得到的方程,需要进行数值求解。选取组件模块化特征线方法(MOC)作为二维输运求解方法,对中子输运方程、共轭中子输运方程以及广义共轭输运方程进行求解。
特征线方法是求解中子输运方程的成熟而有效方法,能够实现一般中子输运方程的求解,其求解流程如图2所示。该求解流程一般称为源迭代法,是一种非常成熟的方法,描述如下:开始任意假定一个初始的裂变源分布,并假定一个初始的特征值,将其带入方程(1)得到一个普通的非齐次方程组,通过组件模块化特征线方法进行内迭代,即对某一群进行求解区域的扫描,在各群内迭代收敛并且能群扫描完成,可以求得更新的中子通量分布和特征值,利用更新的中子通量分布和特征值,进行下一次的区域扫描,直到收敛为止。
但是对于共轭中子输运方程,以及广义共轭方程,需要对求解流程作一点修改。
对于共轭中子输运方程计算,在计算之前,需要如下操作:
(1)将散射矩阵转置;
(2)将材料的(裂变产生截面)同χ(裂变谱)向量互换;
(3)将所有截面的能群号按如下方式变化:
完成上述操作之后,输运求解器可以完成基本的共轭中子输运方程的求解。当求解广义共轭输运方程,还需要如下操作:
(1)将建立共轭源作为问题的外源项;
(2)对于方程(10)将裂变源的更新按如下形式进行:
式中:
FSg——第g群裂变源;
——裂变产生截面;
keff——有效增殖因子;
χ——裂变谱;
Γ*——广义共轭通量密度;
F*——中子输运方程裂变算子的共轭算子;
Φ——前向中子通量密度;
Φ*——共轭中子通量密度。
其求解流程如图3所示。该求解流程同上面描述的图2的求解流程基本一致,但主要有两点区别:(1)不需要进行特征值的更新,方程中的特征值使用的是一般中子输运方程求得的特征值;(2)每一次外迭代得到的通量进行裂变源更新时,都需要采用式(18)的形式进行。
对于方程(15)将裂变源的更新按如下形式进行:
式中:
Γ——广义通量密度;
F——中子输运方程裂变算子。
当采用输运求解器完成共轭中子输运方程或者广义共轭输运方程的求解之后,需要将获得的角通量在角度上反转,在能群上颠倒,从而为下一步计算提供正确的通量信息。
步骤3:在步骤1建立的方法的基础上,该步针对子群共振自屏方法,建立相应的子群广义共轭方程;采用步骤2建立的求解方法,选取组件模块化特征线方法(MOC)作为二维输运求解方法,求解子群广义共轭方程,为共振自屏截面的灵敏度系数计算提供广义共轭通量。
在共振自屏计算过程中,往往需要求解慢化方程或者输运方程来计算问题相关的权重函数φ(E),以得到问题相关的多群自屏截面:
——核素j的x反应的第g群共振截面。
其中φ(E)往往通过求解一个慢化方程或者输运方程得到,设其算子形式为:
Lφ(E)=Q(E)公式(21)
式中:
φ(E)——权重谱;
L——相应的慢化或输运算子;
Q——源项。
各种核素会参与权重谱的求解,其截面扰动会影响到作为权重的通量,从而影响多群自屏截面。如果多群截面为α,共振核素共振自屏截面对非共振核素截面的扰动的灵敏度系数为定义为:
一般地,多群共振自屏截面可以表达成式(20)的形式,为简化表达,省略的核素标志j,α表示某核素某反应的某群截面,将等式两边对α求微分,有:
式(23)的右端第一项称为α扰动的直接影响,该项表示α扰动对该能量段共振截面的直接影响,只有在α为一些基本物理量(共振能量、共振峰宽度等)时该项才非零;右端第二项称为α扰动的间接影响,即通过影响作为权重的通量而间接影响共振截面值。
直接项的计算可以通过具体公式求导或者采用直接扰动方法进行求解。
对于间接项,考虑到通量对α的偏导,将式(21)两边对α求偏导有:
两边同乘以共轭函数在相空间积分得到:
其中<,>代表在求解空间积分,的下标x和g表示广义共轭函数同共振截面和共振能群是相关的。根据共轭方程的性质,有:
其中L*为L的共轭算符。对比式(25)和式(26),不妨取广义共轭方程为:
那么只要得到带入到式(25)右端,就能得到:
下面将上述对一般共振方法的推导结果推广到子群共振计算方法中。
子群共振计算方法中将能群g共振截面从其最小值到最大值划分成若干区间,称为子群。子群截面的定义为:
式中,下标g和i分别表示能群和子群;σx,g,i为第g能群的第i子群的子群截面;φ(E)为中子通量密度。σx,g(E)为能群g的能量相关截面,ΔEi的范围为ΔEi∈{E|σx,i≤σx,g(E)≤σx,i+1}。
本发明中共振计算采用的子群共振自屏方法基于物理概率表,概率表给出了一个给定能群内的一组截面值和相应的权重。基于窄共振近似和物理概率表,可以将子群截面按能群给出如下形式:
式中:N为总子群数,σx,g为第g群的多群截面,σ0为背景截面,σt,g,i为能群g的第i子群的总截面,pg,i为能群g的第i子群的子群概率。
基于帕德近似方法可根据式(30)得到子群截面σx,g,i和子群概率pg,i。
在得到子群截面σx,g,i和子群概率pg,i后,由于各子群对应的能量段必然包含于能群g中且不会重复,根据能群平均截面的定义,可将能群g的平均微观截面写成子群的形式:
Ω·▽φg,i(r,Ω)+Σt,g,i(r)φg,i(r,Ω)=Qs,g,i(r,Ω)公式(32)
式中:
φg,i(r,Ω)——能群g的第i子群的中子通量;
Σt,g,i(r)——能群g的第i子群的宏观总截面;
Qs,g,i(r,Ω)——散射源项;
将上述方程写成算子形式为:
Lgφg=Qg公式(33)
式中:
Lg——第g个共振能群的子群输运方程的输运算子;
Qg——第g个共振能群的子群输运方程的源项。
对比式(21)和式(31)可以发现,二者的区别在于前者是在能量上积分,而后者是在子群上求和,前者中子通量一般是通过空间无关的慢化方程所得,而后者中子通量是通过空间相关的子群输运方程所得,因此很容易将前面一般推导的结果推广到子群方法上面,所以根据式(28),得到子群共振方法的共振截面灵敏度系数为:
其中为广义子群共轭方程的解,广义子群共轭方程为:
式中:
——Lg的共轭算子;
——子群广义共轭通量;
——子群广义共轭源。
其中
式中
——第i个子群的源项;
σx,g,i——第i个子群截面。
注意到式(34)中存在子群参数灵敏度系数,考虑在上面采用帕德近似方法求解子群参数的过程中,采用直接扰动方法求解子群参数灵敏度系数。对于各共振核素,逐群扰动其连续能量截面(体现为扰动其共振积分表),设扰动百分比为δ,则根据差商代替微分的方法求得子群参数的灵敏度系数,即
公式(37)
式中为未扰动的子群参数(子群截面σx,g,i或子群概率pg,i),和分别为正向和负向扰动α时的子群参数;δ为α的扰动百分比。
此外,式(34)中存在源项Qg对α的灵敏度系数,这一项的求解同样采用广义微扰理论。本发明采用的子群共振自屏方法中,源项分成快群散射源项Qf,g和上游共振能群散射源项Qr,g,分别为:
对于共振能群g,分别对上述两种源建立广义共轭方程:
式中:
——快群输运算子;
——广义共轭通量;
g’——快群标号;
Σs,g'→g——快群到共振能群的散射截面。
而对于上游共振能群g’,有
式中:
——子群输运算子的共轭算子;
——广义共轭通量;
g’——上游共振能群标号;
i——第g’群的子群标号;
Σs,g'→g,i——上游共振能群g’到当前共振能群g的第i个子群的散射截面。
采用输运求解器求解方程式(40)和(41),得到对应的广义共轭通量。由下式可以求得源项的灵敏度系数
式中:
Qg——快群散射源项或上游共振能群散射源项;
Σs——快群或上游共振能群到当前共振能群的散射截面;
Q'——求解源项的输运方程的右端源项。
步骤4:针对反应堆物理栅格计算的重要参数,根据步骤1建立的方法获取与之对应的广义(共轭)方程,利用步骤2的输运求解器,求解获取通量、广义共轭通量以及广义通量。从而得到各重要参数的显式灵敏度系数。
设某重要参数为R,可以表示成通量或者共轭通量的线性泛函,设为:
注意这里R具有一般性,只要是能表示成上述形式的重要参数,都可以用本发明的方法进行灵敏度系数的求解。
建立对应的广义共轭源项:
利用步骤2开发的输运求解器分别求解:
M*Γ*=Q1公式(46)
MΓ=Q2公式(47)
获取通量、广义共轭通量以及广义通量后,显式灵敏度系数可由下式计算得到:
步骤5:根据第3步得到的共振自屏截面的灵敏度系数和第4步得到的显式灵敏度系数,计算出考虑了隐式影响的灵敏度系数。
其中j表示共振核素。
如图4所示,给出了典型压水堆栅元问题的一些核素的特征值显式和隐式灵敏度系数比较。传统方法只能得到显式灵敏度系数而往往忽略隐式灵敏度系数,本发明能够计算隐式灵敏度系数,从图中可以看出,某些核素的隐式灵敏度系数相对显式灵敏度系数很大,如果忽略会造成很大的误差。所以考虑隐式灵敏度系数是必要的。
Claims (1)
1.一种获取反应堆物理栅格计算重要参数灵敏度系数的方法,其特征在于:该方法包括以下步骤:
步骤1:基于微扰理论,将多输入参数的反应堆物理计算需要大量扰动计算才能得到灵敏度系数的问题,转化为只需一次前向计算和一次共轭计算就能获取某一计算结果对所有输入参数的灵敏度系数的问题;
设反应堆物理计算中某重要参数为R,以下称为响应,表示成通量或者共轭通量的线性泛函,设为:
式中:
Φ——中子角通量密度;
Φ*——共轭中子角通量密度;
ξ——问题求解的空间;
Σ(ξ)——同求解空间相关的截面数据;
H1[Σ(ξ)],H2[Σ(ξ)]——依赖于截面数据的算子;
建立与响应R的表达式对应的广义共轭输运方程:
式中:
M——输运算子;
Γ——广义通量密度;
M*——输运算子的共轭算子;
Γ*——广义共轭通量密度;
步骤2:采用广泛应用的输运求解方法——组件模块化特征线方法MOC作为二维输运方程的求解方法,对中子输运方程、共轭中子输运方程,以及广义共轭输运方程进行求解;
特征线方法能够实现一般中子输运方程的求解,但是对于共轭中子输运方程以及广义共轭输运方程,需要对求解流程作一点修改;
对于共轭中子输运方程计算,在计算之前,需要如下操作:
(1)将散射矩阵转置;
(2)将材料的裂变产生截面同裂变谱χ向量互换;
(3)将所有截面的能群号按如下方式变化:
完成上述操作之后,输运求解器能够完成基本的共轭中子输运方程的求解;
当求解广义共轭输运方程,还需要如下操作:
(1)将建立共轭源作为问题的外源项;
(2)对于公式(10)将裂变源的更新按如下形式进行:
式中:
FSg——第g群裂变源;
——裂变产生截面;
keff——有效增殖因子;
χ——裂变谱;
Γ*——广义共轭通量密度;
F*——中子输运方程裂变算子的共轭算子;
Φ——前向中子通量密度;
Φ*——共轭中子通量密度;
其计算流程与传统的中子输运方程计算流程基本一致,不同点在于:(1)不需要进行特征值的更新,方程中的特征值使用的是中子输运方程求解得到的特征值;(2)每一次外迭代得到的通量用于裂变源更新时采用公式(18);
对于方程(15)将裂变源的更新按如下形式进行:
式中:
Γ——广义通量密度;
F——中子输运方程裂变算子;
当采用输运求解器完成共轭中子输运方程或者广义共轭输运方程的求解之后,需要将获得的角通量在角度上反转,在能群上颠倒,从而为下一步计算提供正确的通量信息;
步骤3:在步骤1建立的方法的基础上,该步针对子群共振计算方法,建立相应的子群广义共轭方程;采用步骤2建立的求解方法,选取组件模块化特征线方法MOC作为二维输运求解方法,求解子群广义共轭方程,为共振自屏截面的灵敏度系数计算提供广义共轭通量;
子群方法中,共振自屏截面的表达式为:
式中:
g——能群标号;
i——子群标号;
N——子群总数;
σx(E)——x反应的连续能量截面;
φ(E)——权重谱;
σx,g,i——共振截面x的第g共振能群的第i个子群截面;
公式(31)中的子群参数通过帕德近似求解,而子群通量通过求解子群输运方程得到;
设子群输运方程为:
式中:
g——共振能群标号;
i——子群标号;
r——空间位置变量;
Ω——角度变量;
φg,i(r,Ω)——能群g的第i子群的中子通量密度;
Σt,g,i(r)——能群g的第i子群的宏观总截面;
Qs,g,i(r,Ω)——散射源项;
将上述方程写成算子形式为:
Lgφg=Qg公式(33)
式中:
Lg——第g个共振能群的子群输运方程的输运算子;
Qg——第g个共振能群的子群输运方程的源项;
利用输运求解器求解子群输运方程,得到子群通量后,对子群参数加权得到共振自屏截面;
子群共振方法的共振截面灵敏度系数计算公式为:
其中为广义子群共轭方程的解。它由该共振自屏截面的子群输运方程对应的广义子群共轭方程解得,该广义子群共轭方程为:
式中:
——Lg的共轭算子;
——子群广义共轭通量;
——子群广义共轭源;
其中
式中
σx,g,i——第i个子群截面;
注意到式(34)中存在子群参数灵敏度系数,考虑采用帕德近似方法求解子群参数的过程中,采用直接扰动方法求解子群参数灵敏度系数;对于各共振核素,逐群扰动其连续能量截面,体现为扰动其共振积分表,设扰动百分比为δ,则根据差商代替微分的方法求得子群参数的灵敏度系数,即
公式(37)
式中为未扰动的子群参数即子群截面σx,g,i或子群概率pg,i,和分别为正向和负向扰动α时的子群参数;δ为α的扰动百分比;
此外,式(34)中存在源项Qg对α的灵敏度系数,这一项的求解同样采用步骤1建立的方法;采用的子群共振自屏计算方法中,源项分成快群散射源项Qf,g和上游共振能群散射源项Qr,g,分别为:
对于共振能群g,分别对上述两种源建立广义共轭方程:
式中:
——快群输运算子;
——广义共轭通量;
g’——快群标号;
Σs,g'→g——快群到共振能群的散射截面;
而对于上游共振能群g’,有
式中:
——子群输运算子的共轭算子;
——广义共轭通量;
g’——上游共振能群标号;
i——第g’群的子群标号;
Σs,g'→g,i——上游共振能群g’到当前共振能群g的第i个子群的散射截面;
采用输运求解器求解方程式(40)和(41),得到对应的广义共轭通量;由下式求得源项的灵敏度系数
式中:
Qg——快群散射源项或上游共振能群散射源项;
Σs——快群或上游共振能群到当前共振能群的散射截面;
Q'——求解源项的输运方程的右端源项;
步骤4:针对反应堆物理栅格计算的重要参数,根据步骤1建立的方法获取与之对应的广义共轭方程,利用步骤2的输运求解器,求解获取通量、共轭通量以及广义通量;从而得到各重要参数的显式灵敏度系数;
设某重要参数为R,表示成通量或者共轭通量的线性泛函,设为:
注意这里R具有一般性,只要是能表示成上述形式的重要参数,都能够进行灵敏度系数的求解;
建立对应的广义共轭源项:
利用输运求解器分别求解:
M*Γ*=Q1公式(46)
MΓ=Q2公式(47)
获取通量、广义共轭通量以及广义通量后,显式灵敏度系数由下式计算得到:
步骤5:根据步骤3得到的共振自屏截面的灵敏度系数和步骤4得到的显式灵敏度系数,计算出考虑了隐式影响的灵敏度系数;
其中j表示共振核素。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510726692.2A CN105426342B (zh) | 2015-10-30 | 2015-10-30 | 一种获取反应堆物理栅格计算重要参数灵敏度系数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510726692.2A CN105426342B (zh) | 2015-10-30 | 2015-10-30 | 一种获取反应堆物理栅格计算重要参数灵敏度系数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105426342A true CN105426342A (zh) | 2016-03-23 |
CN105426342B CN105426342B (zh) | 2016-10-12 |
Family
ID=55504557
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510726692.2A Active CN105426342B (zh) | 2015-10-30 | 2015-10-30 | 一种获取反应堆物理栅格计算重要参数灵敏度系数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105426342B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106096018A (zh) * | 2016-06-24 | 2016-11-09 | 西安交通大学 | 一种在多群数据库中考虑共振弹性散射的方法 |
CN106096182A (zh) * | 2016-06-24 | 2016-11-09 | 西安交通大学 | 一种反应堆嵌入式共振自屏计算方法 |
CN106126480A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种获取反应堆多群数据库中的多群Pn散射矩阵的方法 |
CN106126926A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种能够处理温度分布效应的子群共振计算方法 |
CN106126818A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种自适应的在线产生子群参数的方法 |
CN106126928A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 获得固态及液态氟盐热中子散射数据库的方法及数据库 |
CN106169019A (zh) * | 2016-06-24 | 2016-11-30 | 西安交通大学 | 一种基于敏感性和不确定分析的目标精度评估方法 |
CN106202862A (zh) * | 2016-06-24 | 2016-12-07 | 西安交通大学 | 一种针对压水堆栅元非均匀共振积分表的制作方法 |
CN106202611A (zh) * | 2016-06-24 | 2016-12-07 | 西安交通大学 | 一种适用于核反应堆组件物理热工耦合的在线计算方法 |
CN107092782A (zh) * | 2017-04-05 | 2017-08-25 | 西安交通大学 | 一种处理共振干涉效应的共振伪核素方法 |
CN107423542A (zh) * | 2017-04-13 | 2017-12-01 | 西安交通大学 | 一种适用于逐棒计算的非均匀泄漏修正方法 |
CN110580935A (zh) * | 2019-08-16 | 2019-12-17 | 西安交通大学 | 一种获取全堆有效共振自屏截面的方法 |
WO2021014820A1 (ja) * | 2019-07-19 | 2021-01-28 | 三菱重工業株式会社 | 核特性の予測方法及び核特性予測装置 |
CN112380744A (zh) * | 2020-11-03 | 2021-02-19 | 中国核动力研究设计院 | 一种加速求解广义共轭中子输运方程的方法 |
CN112632771A (zh) * | 2020-12-21 | 2021-04-09 | 中国核动力研究设计院 | 一种求解燃耗计算响应灵敏度的方法及设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005148079A (ja) * | 2005-01-21 | 2005-06-09 | Toshiba Corp | 炉心性能計算装置 |
-
2015
- 2015-10-30 CN CN201510726692.2A patent/CN105426342B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005148079A (ja) * | 2005-01-21 | 2005-06-09 | Toshiba Corp | 炉心性能計算装置 |
Non-Patent Citations (2)
Title |
---|
卑华等: "SONG多群数据库的设计和制作", 《核动力工程》 * |
杨军等: "中国实验快堆Keff计算值对核数据的灵敏度和不确定性分析", 《核动力工程》 * |
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106202862B (zh) * | 2016-06-24 | 2018-08-21 | 西安交通大学 | 一种针对压水堆栅元非均匀共振积分表的制作方法 |
CN106096018B (zh) * | 2016-06-24 | 2019-04-16 | 西安交通大学 | 一种在多群数据库中考虑共振弹性散射的方法 |
CN106126480A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种获取反应堆多群数据库中的多群Pn散射矩阵的方法 |
CN106126926A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种能够处理温度分布效应的子群共振计算方法 |
CN106126818A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种自适应的在线产生子群参数的方法 |
CN106126928A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 获得固态及液态氟盐热中子散射数据库的方法及数据库 |
CN106169019A (zh) * | 2016-06-24 | 2016-11-30 | 西安交通大学 | 一种基于敏感性和不确定分析的目标精度评估方法 |
CN106202862A (zh) * | 2016-06-24 | 2016-12-07 | 西安交通大学 | 一种针对压水堆栅元非均匀共振积分表的制作方法 |
CN106202611A (zh) * | 2016-06-24 | 2016-12-07 | 西安交通大学 | 一种适用于核反应堆组件物理热工耦合的在线计算方法 |
CN106169019B (zh) * | 2016-06-24 | 2018-07-03 | 西安交通大学 | 一种基于敏感性和不确定性分析的目标精度评估方法 |
CN106126928B (zh) * | 2016-06-24 | 2018-07-03 | 西安交通大学 | 获得固态及液态氟盐热中子散射数据库的方法及数据库 |
CN106126926B (zh) * | 2016-06-24 | 2018-07-03 | 西安交通大学 | 一种能够处理温度分布效应的子群共振计算方法 |
CN106096182A (zh) * | 2016-06-24 | 2016-11-09 | 西安交通大学 | 一种反应堆嵌入式共振自屏计算方法 |
CN106096018A (zh) * | 2016-06-24 | 2016-11-09 | 西安交通大学 | 一种在多群数据库中考虑共振弹性散射的方法 |
CN106096182B (zh) * | 2016-06-24 | 2019-04-09 | 西安交通大学 | 一种反应堆嵌入式共振自屏计算方法 |
CN106126480B (zh) * | 2016-06-24 | 2018-08-28 | 西安交通大学 | 一种获取反应堆多群数据库中的多群Pn散射矩阵的方法 |
CN106126818B (zh) * | 2016-06-24 | 2018-12-18 | 西安交通大学 | 一种自适应的在线产生子群参数的方法 |
CN106202611B (zh) * | 2016-06-24 | 2019-04-09 | 西安交通大学 | 一种适用于核反应堆组件物理热工耦合的在线计算方法 |
CN107092782A (zh) * | 2017-04-05 | 2017-08-25 | 西安交通大学 | 一种处理共振干涉效应的共振伪核素方法 |
CN107423542B (zh) * | 2017-04-13 | 2020-01-17 | 西安交通大学 | 一种适用于逐棒计算的非均匀泄漏修正方法 |
CN107423542A (zh) * | 2017-04-13 | 2017-12-01 | 西安交通大学 | 一种适用于逐棒计算的非均匀泄漏修正方法 |
JP7202984B2 (ja) | 2019-07-19 | 2023-01-12 | 三菱重工業株式会社 | 核特性の予測方法及び核特性予測装置 |
WO2021014820A1 (ja) * | 2019-07-19 | 2021-01-28 | 三菱重工業株式会社 | 核特性の予測方法及び核特性予測装置 |
JP2021018143A (ja) * | 2019-07-19 | 2021-02-15 | 三菱重工業株式会社 | 核特性の予測方法及び核特性予測装置 |
CN110580935A (zh) * | 2019-08-16 | 2019-12-17 | 西安交通大学 | 一种获取全堆有效共振自屏截面的方法 |
CN110580935B (zh) * | 2019-08-16 | 2020-08-28 | 西安交通大学 | 一种获取全堆有效共振自屏截面的方法 |
CN112380744A (zh) * | 2020-11-03 | 2021-02-19 | 中国核动力研究设计院 | 一种加速求解广义共轭中子输运方程的方法 |
CN112632771A (zh) * | 2020-12-21 | 2021-04-09 | 中国核动力研究设计院 | 一种求解燃耗计算响应灵敏度的方法及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN105426342B (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105426342A (zh) | 一种获取反应堆物理栅格计算重要参数灵敏度系数的方法 | |
Lee et al. | MCS–A Monte Carlo particle transport code for large-scale power reactor analysis | |
CN103177154B (zh) | 一种获得核燃料组件共振参数的方法 | |
Pfeifer et al. | Simulation of anyons with tensor network algorithms | |
CN103106301A (zh) | 一种基于蒙特卡罗方法与特征线方法耦合的计算辐射屏蔽的方法 | |
CN105426659B (zh) | 获取不同燃耗下有效增殖因子对截面的灵敏度系数的方法 | |
CN106202866A (zh) | 一种稳定精确的反应堆物理热工耦合计算方法 | |
CN105975977A (zh) | 规则辅助的数据驱动建模方法 | |
CN103116667B (zh) | 一种获取聚变堆实验包层模块中子学参数的方法 | |
CN105808503A (zh) | 针对反应堆逐棒计算中解析求解栅元不连续因子的方法 | |
Li et al. | CAD modeling study on FLUKA and OpenMC for accelerator driven system simulation | |
EP2287853A1 (en) | A computer implemented method for modelling a nuclear reactor core and a corresponding computer program product | |
CN106372763A (zh) | 基于改进的支持向量机的电力负荷预测修正方法 | |
EP4202792A1 (en) | Method and system for determining many-body localization, and storage medium and program product | |
CN106202865A (zh) | 一种计算中子输运离散节块法中任意阶系数的方法 | |
Wang et al. | Exact-to-precision generalized perturbation theory for eigenvalue problems | |
Cornejo et al. | Iteration methods with multigrid in energy for eigenvalue neutron diffusion problems | |
Bao et al. | A data-driven approach to scale bridging in system thermal-hydraulic simulation | |
CN107391443A (zh) | 一种稀疏数据异常检测方法及装置 | |
EP2287854B1 (en) | A computer implemented method for modelling a nuclear reactor core and a corresponding computer program product | |
CN116579205A (zh) | 一种压水堆核热耦合的计算方法 | |
Schlünz et al. | Research reactor in-core fuel management optimisation using the multiobjective cross-entropy method | |
CN105608267A (zh) | 一种多变量全局优化算法 | |
Xu et al. | Control variate multifidelity estimators for the variance and sensitivity analysis of mesostructure–structure systems | |
EP2287855B1 (en) | A computer implemented method for modelling a nuclear reactor core and a corresponding computer program product |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |