基于函数展开的连续能量确定论中子输运计算方法技术

本发明专利技术公开了一种基于函数展开的连续能量确定论中子输运计算方法，按照核反应微观截面特点将能量区划分为热中子能区、可辨共振能区、不可辨共振能区和快中子能区；其中每个能量区间又划分为若干个能量段；在在不同的能量区间根据需要选取不同的基函数，并对不同的能量区间内的中子角通量密度进行函数展开，将能量相空间内关于中子角通量密度的方程，转换成频域空间内关于中子角通量密度展开系数矩的方程；在不同的能量区间内求解该方程获得可辨共振能量区间各能量段的中子角通量密度系数矩和中子源强展开系数矩及中子角通量密度函数，进而实现对各种核素的各种核反应率的刻画，定量描述裂变能的释放、核素的消耗与生产、中子束流强度、辐照环境参数等。中子束流强度、辐照环境参数等。中子束流强度、辐照环境参数等。

【技术实现步骤摘要】
基于函数展开的连续能量确定论中子输运计算方法


[0001]本专利技术涉及中子输运过程模拟计算和反应堆物理分析领域，具体涉及连续能量确定论中子输运计算方法。

技术介绍

[0002]裂变核反应堆是进行可控链式裂变反应的装置，可以提供能量、放射性同位素、中子束流、辐照环境等，供工业生产和科学研究使用。在其堆芯活性区内，为了刻画中子在介质中的输运过程，需要求解中子输运方程。该数学物理方程是关于中子角通量密度在三维空间、一维中子能量及二维飞行方向等共六维相空间内的微分
‑
积分方程，属于线性Boltzmann输运方程。该方程的解析解仅在非常简化的情况下才可以获得，具有重要的理论分析价值，但不具备工程实用意义。因此，核反应堆堆芯物理工程计算一般都使用数值计算方法由计算机完成。现阶段中子输运方程的数值求解主要分为两大类方法，一类是基于网格划分与函数展开的确定论方法，用有限个相空间网格内的平均中子通量密度，逼近实际物理问题中连续分布的中子通量密度分布函数；另一类是基于随机抽样的概率论方法，又叫蒙特卡罗方法，把中子通量密度转化成随机过程统计期望的积分形式，通过大量随机抽样和统计估计，获得实际物理问题中子通量密度分布函数的估计值。
[0003]针对蒙特卡洛方法，当且仅当目标物理量的样本贡献数目足够多时，才可以给出可信的结果；尤其是对于分布式物理量，其计算代价将剧增；另外计算大规模问题的精细分布物理量时，该方法从理论上还面临可能出现的伪收敛甚至不收敛问题。因此在结构几何下，不适合进行海量的工程计算。在确定论方法中，对结构几何的计算效率较高，且针对中子输运方程的空间和角度的联合相空间，已经发展了非常成熟的数值离散和迭代求解技术，其中的近似已经越来越少。然而，与空间和角度相空间形成鲜明对比的是，针对中子能量相空间，目前仍在采用传统的多群方法近似处理。在多群近似处理中，整个中子能量范围可以划分为若干(G)个区间，称之为G个能群；在每个能群内认为所有中子与原子核发生相互作用的概率是一样，称之为该能群的平均微观截面；需要先依据近似能谱和评价核数据库获取非共振能群的平均微观截面，再在共振能群内进行共振自屏计算获得共振能群的平均微观截面，才可以进行真正的中子输运模拟计算。该方法通过多群近似方法把原来含有连续能量的方程离散成G个能群方程的求解问题，通过对能群方程的迭代求解即可得到每个能群内的群通量数据，而能群的数目则根据所研究问题的性质和精度要求来决定。
[0004]但是，基于多群近似技术的中子输运过程模拟面临一系列的问题。一方面，多群数据库制作需要预先给定典型中子能谱，无法与具体问题相关；另一方面，共振自屏计算需要考虑复杂的多种影响效应，又包含大量近似修正技术，使得该技术适用范围及修正效果都非常有限，严重限制了核能装置的设计分析工作。
[0005]为了能够彻底消除多群近似带来的问题，并且提高计算精度和计算效率，本专利技术直接从连续能量核数据库出发，提出了基于函数展开的连续能量确定论中子输运计算方法，可以直接从根本上祛除多群数据库制作和共振计算操作，为核能装置的设计分析提供
更加精确的数据支持。

技术实现思路

[0006]为了克服上述现有技术存在的问题，本专利技术的目的在于提供一种基于函数展开的连续能量确定论中子输运计算方法，可以直接从根本上祛除多群数据库制作和共振计算操作，为核能装置的设计分析提供更加精确的数据支持。
[0007]为了达到上述目的，本专利技术采取了以下计算方案予以实施：
[0008]一种基于函数展开的连续能量确定论中子输运计算方法，步骤如下：
[0009]步骤1：按照核反应微观截面特点将能量区划分为：热中子能区、可辨共振能区、不可辨共振能区和快中子能区；其中每个能量区间又可划分为若干个能量段，能量段的划分越密计算精度越高；
[0010]步骤2：在不同的能量区间根据需要选取不同的基函数R
n
(u)，n＝0,1,2,
…
,N，并对不同的能量区间内的中子角通量密度进行函数展开，将中子角通量密度展开成未知系数与不同阶基函数的乘积之和的形式，如下所示：
[0011][0012]φ(r,Ω,u)
‑‑
在空间r处沿角度Ω飞行且飞行能量的对数能量为u的中子角通量密度(cm
‑2s
‑1)；其中u＝lnE，E为出射中子飞行能量，单位为：MeV；
[0013]ψ
n
(r,Ω)
‑‑
待求的中子角通量密度未知展开系数矩；
[0014]不同能量区间内选取的基函数R
n
(u)，可精确刻画该能量区间内中子角通量密度随能量的波动变化；基函数展开阶数N越大，越能满足对不同能区的刻画精度；
[0015]步骤3：对应到快中子能区，核反应微观截面随入射中子能量连续光滑变化，使相应的中子角通量密度也呈现连续光滑变化的趋势；针对该特点，选择光滑变化的基函数，仅考虑本能区自散射源项(无较高能区散射到该能区下散射源项)，通过对快中子能区的中子角通量密度进行基函数展开，将能量相空间内关于中子角通量密度的方程，转换成频域空间内关于中子角通量密度展开系数矩的方程，该能区函数展开后的中子输运方程格式为公式(1)，求解该方程即获得快中子能量区间各能量段的中子角通量密度系数矩及中子角通量密度函数；避免了多群方法中计算需要预设能谱的情况，使计算更精准；
[0016][0017]式中：
[0018]Ω
‑‑
空间中出射中子角度；
[0019]Ω'
‑‑
空间中入射中子角度；
[0020]r
‑‑
空间位置；
[0021]u
–
出射中子能量的对数能量；
[0022]u
’–
入射中子能量的对数能量；
[0023]N
k
(r)
‑‑
空间r处第k种核素的原子核密度(cm
‑3)；
[0024]σ
t,k
(u)
‑‑
第k种核素与对数能量为u的中子发生核反应的微观总截面(cm2)；σ
s,k
(u'
→
u,Ω'
→
Ω)
‑‑
第k种核素与对数能量为u'、飞行方向为Ω'的中子发生散射核反应并产生对数能量为u、飞行方向为Ω中子的微观散射截面(cm2)；
[0025]q
n
(r,Ω)
‑‑
待求的中子源强未知展开系数矩；
[0026]R
1n
(u)
‑‑
第I1区的基函数；
[0027]I1‑‑
快中子能量区间；
[0028]步骤4：利用步骤3中计算得到的快中子能量区间各能量段的中子角通量密度系数矩，可得到快中子能区到不可辨共振能区的下散射源项为，
[0029][0030]R
1n
(u)
...

【技术保护点】

【技术特征摘要】
1.一种基于函数展开的连续能量确定论中子输运计算方法，其特征在于：步骤如下：步骤1：按照核反应微观截面特点将能量区划分为：热中子能区、可辨共振能区、不可辨共振能区和快中子能区；其中每个能量区间又可划分为若干个能量段；步骤2：在不同的能量区间根据需要选取不同的基函数R
n
(u)，n＝0,1,2,
…
,N，并对不同的能量区间内的中子角通量密度进行函数展开，将中子角通量密度展开成未知系数与不同阶基函数的乘积之和的形式，如下所示：φ(r,Ω,u)
‑‑
在空间r处沿角度Ω飞行且飞行能量的对数能量为u的中子角通量密度，单位为：cm
‑2s
‑1；其中u＝lnE，E为出射中子飞行能量，单位为：MeV；ψ
n
(r,Ω)
‑‑
待求的中子角通量密度未知展开系数矩；不同能量区间内选取的基函数R
n
(u)，能精确刻画该能量区间内中子角通量密度随能量的波动变化；基函数展开阶数N越大，越能满足对不同能区的刻画精度；步骤3：选择光滑变化的基函数，仅考虑本能区自散射源项，通过对快中子能区的中子角通量密度进行基函数展开，将能量相空间内关于中子角通量密度的方程，转换成频域空间内关于中子角通量密度展开系数矩的方程，该能区函数展开后的中子输运方程格式为公式(1)，求解该方程即获得快中子能量区间各能量段的中子角通量密度系数矩及中子角通量密度函数；式中：Ω
‑‑
空间中出射中子角度；Ω'
‑‑
空间中入射中子角度；r
‑‑
空间位置；u
–
出射中子能量的对数能量；u
’–
入射中子能量的对数能量；N
k
(r)
‑‑
空间r处第k种核素的原子核密度，单位为：cm
‑3；σ
t,k
(u)
‑‑
第k种核素与对数能量为u的中子发生核反应的微观总截面，单位为：cm2；σ
s,k
(u'
→
u,Ω'
→
Ω)
‑‑
第k种核素与对数能量为u'、飞行方向为Ω'的中子发生散射核反应并产生对数能量为u、飞行方向为Ω中子的微观散射截面，单位为：cm2；q
n
(r,Ω)
‑‑
待求的中子源强未知展开系数矩；R
1n
(u)
‑‑
第I1区的基函数；I1‑‑
快中子能量区间；步骤4：利用步骤3中计算得到的快中子能量区间各能量段的中子角通量密度系数矩，
得到快中子能区到不可辨共振能区的下散射源项为，R
1n
(u)
‑‑
第I1区的基函数；I1‑‑
快中子能量区间；I2‑‑
不可辨共振能量区间；然后需要考虑不可辨共振能量区间内的自散射源项，在不可辨共振能区，微观截面随入射中子能量变化的激发曲线上具有过度密集的共振峰，实验上无法测定指...

【专利技术属性】
技术研发人员：李云召，刘浩泼，黄星，吴宏春，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：