针对三维瞬态多群中子扩散方程的非结构网格变分节块方法

1.本发明涉及的是一种核反应堆建造领域的技术，具体是一种针对三维瞬态多群中子扩散方程的非结构网格变分节块方法。


背景技术：

2.反应堆物理的核心任务是准确高效地求解中子输运/扩散方程，多群瞬态中子输运/扩散方程的数值求解需要对时间和空间进行离散。现有的变分节块法(vnm)通过对每个典型节块构造一系列响应矩阵，可以较好地处理四边形、正六边形及正三角形的空间离散，然而随着核科学与技术的不断发展，各种具有强烈的非结构特点的新概念反应堆设计相继提出，现有的vnm已不能满足新的工程设计需求。另外在瞬态计算中，vnm的时间推进需要对响应矩阵进行多次重构，增大计算代价，这些因素限制vnm在工程实践中的应用性。


技术实现要素：

3.本发明针对现有vnm无法应用于复杂几何的新概念反应堆设计以及瞬态非结构网格vnm计算代价相对较大的问题，提出一种针对三维瞬态多群中子扩散方程的非结构网格变分节块方法，提高瞬态计算下变分节块法的几何适应性，并通过刚性限制法(scm)减少非结构网格vnm在时间推进中的响应矩阵重构次数以提高瞬态计算效率，能够准确模拟实际不规则几何的复杂问题，从而提高vnm在工程实践中的应用价值。
4.本发明是通过以下技术方案实现的：
5.本发明涉及一种针对三维瞬态多群中子扩散方程的非结构网格变分节块方法，首先使用非结构vnm对稳态中子扩散方程计算后初始化中子通量密度、缓发中子先驱核以及动力学截面，然后在瞬态计算的每个时间步内通过迭代更新反应截面以及更新动力学频率并重构响应矩阵后，使用非结构vnm求解刚性限制法(scm)方程，实现任意几何形状下的中子的瞬态过程的模拟。技术效果
6.本发明通过网格的空间映射，将任意形状的三角形网格映射到标准节块中，从而使得vnm可以处理任意几何的实际问题；本发明利用刚性限制法(scm)处理瞬态中子扩散方程的时间项，消除方程的刚性，减少非结构vnm的计算代价。
7.与现有的瞬态中子扩散求解器相比，本发明在实现较高的计算精度的同时，具有更好的几何兼容性，可以计算任意几何的实际问题。同时可在大时间步长下进行计算，具有较好的数值稳定性，避免响应矩阵的多次重构。在以dodds基准题为代表的复杂几何瞬态中子扩散问题中与参考解误差在1％以内；时间步由传统的0.01s减少至0.1s。
附图说明
8.图1为本发明流程图；
9.图2为实施例坐标变换示意图；
10.图3为dodds基准题几何示意图；
11.图中：(a)基准题轴向布置，(b)基准题径向布置；
12.图4为实施例效果示意图。
具体实施方式
13.如图1所示，为本实施例涉及一种针对三维瞬态多群中子扩散方程的非结构网格变分节块方法，首先使用非结构vnm对稳态中子扩散方程计算后初始化中子通量密度、缓发中子先驱核以及动力学截面，然后在瞬态计算的每个时间步内通过迭代更新反应截面以及更新动力学频率并重构响应矩阵后，使用非结构vnm求解刚性限制法(scm)方程。
14.本发明是通过以下技术方案实现的：
15.步骤1、将特征值问题(evp)的区域划分为若干个三角形非结构网格，应用变分原理，建立包含瞬态中子扩散方程和边界条件在内的泛函fv[φ，j]，即等价为寻找一组φ，j使得泛函fv[φ，j]取得极小值，具体为：[φ，j]取得极小值，具体为：其中：dv为三角形节块体积微元，dγ为三角形的表面微元。
[0016]
步骤2、对变量φ，s，j
γ
在空间上采用一系列标准正交多项式展开，得到：其中：f
t
(r)，分别为三角形节块内部以及表面完备正交的多项式，通过gramm-schmidt正交化得到；分别为上述三个变量的展开矩。
[0017]
在实际的响应矩阵计算中，通过构建标准三角形节块进行gramm-schmidt正交化得到正交基函数，并利用坐标映射关系求解实际节块下的响应矩阵，具体如图2所示，实际节块坐标系为(ξ，η，τ)，标准节块坐标系为(x，y，z)，根据映射关系推导出任意三角形与标准三角形的积分关系满足：其中：vv，vs分别为实际结块与标准节块的体积；|jv|为坐标变换的雅各比行列式，v
ii
′
为克罗内克函数；进一步计算得到离散形式的泛函：其中：a，m为系数矩阵，分别为：
[0018]
标准节块为一等边三角形，如图2所示，其体积vs恒为1，便于数值计算的实现。
[0019]
步骤3、引入经推导后得到响应矩阵方程为：(i-rλπ)j

＝bs，其中：r＝[g i]-1
[g-i]，π为出射入射关系的转移矩阵，λ为由于三角形映射引起的基函数正负关系矩阵。通过裂变源迭代、多群迭代以及群内迭代对响应矩阵方程进行求解可对稳态evp问题进行求解。
[0020]
步骤4、利用scm对三维瞬态多群中子扩散方程进行变形后得到evp形式的瞬态scm方程，具体包括：
[0021]
①
三维瞬态多群中子扩散方程为：三维瞬态多群中子扩散方程为：其中：其中：下标g为第g群，d为扩散系数，φ为中子通量密度，vg为中子速度，∑
x
为各类反应截面，v为一次裂变释放中子数，χg为中子裂变谱，ci(r，t)为缓发中子先驱核，β为缓发中子份额；
[0022]
②
将动力学频率引入三维瞬态多群中子扩散方程以消除方程的刚性，并将动力学频率进一步拆分为通量形状频率ω
s，g
(r，t)以及通量幅值频率ω
t
(t)；
[0023]
③
对缓发中子先驱核ci(r，t)拆分为缓发中子先驱核频率将动力学有效增殖系数kd引入扩散方程得到变形后的三维瞬态多群中子扩散方程：以及以及其中：动力学总截面动力学裂变谱
[0024]
步骤5、通过求解稳态中子扩散方程的解得到瞬态计算的中子通量密度的初值φ(r，t0)，缓发中子先驱核浓度并初始化新动力学频率为0。
[0025]
步骤6、利用步骤1-步骤3计算瞬态scm方程，得到动力学有效增殖系数kd，若kd不满足收敛条件，更新动力学频率ω
s，g
(tn)、动力学总截面∑
′
t，g
(r，t)和动力学裂变谱χ
′g(r，t)，再据此进一步重新构造步骤3中的响应矩阵；
[0026]
所述的动力学有效增殖系数kd收敛条件具体为：当|k
d-1|＜ε时则认为其收敛，结束当前时间步计算。其中ε为预设的收敛限，可取为1e-6。
[0027]
所述的更新动力学频率所述的更新动力学频率其中：vv代表节块体积微元，下标n代表时间步，上标m代表迭代步。
[0028]
步骤7、重复步骤6直至kd收敛，结束当前时刻计算。
[0029]
步骤8、根据已知条件更新反应截面，重复步骤7进行下一时刻的计算，直至达到最大预设时间，从而可以得到待求问题所有时间点下的中子通量密度分布。
[0030]
经过具体实际实验，在intel i7-8700 cpu的pc机上采用上述方法单核运行如图2所示的圆柱瞬态基准题dodds，以测试本方法在复杂非结构几何问题中性能与应用。
[0031]
所有外边界采用真空边界条件；节块内部以及表面的空间展开阶数分别为4阶与2
阶；特征值、裂变源以及中子通量密度的收敛限分别为1e-6、1e-5、1e-5；时间步长设置为0.1s。瞬态计算归一化功率随时间变化的数值结果见表1。
[0032]
如表1所示，在0.1s的时间步长下数值计算结果已趋于收敛。图3为不同方法的dodds基准题的计算结果，其中vitas-t应用本方法提出的三维瞬态多群中子扩散方程的非结构网格变分节块法。计算结果显示vitas-t在0.1s的时间步下与torttd-grs参考解的最大相对误差不超过1％，可以验证该方法在非结构网格复杂几何的dodds瞬态基准题中取得较高的计算精度。
[0033]
表1
[0034]
如表1可见，该计算所采用的时间步长为0.1s，而现有瞬态计算一般采用的时间步长为0.01s；相对现有技术，本方法可采用更大的时间步长，具有更好的数值稳定性。
[0035]
本方法对复杂非结构网格具有良好的适应性，可在大时间步长下准确高效地进行非结构多维瞬态中子扩散方程的求解，从而提高vnm在反应堆设计以及安全分析中的工程实用性。
[0036]
与现有技术相比，本方法将瞬态变分节块法的适用范围推广到非结构几何，可处理任意几何形状的问题；可以采用较大的时间步长以减少响应矩阵重构次数，具有较好的数值稳定性，提高计算效率。
[0037]
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整，本发明的保护范围以权利要求书为准且不由上述具体实施所限，在其范围内的各个实现方案均受本发明之约束。