Adjoint Neutron Flux Calculation Technique Based on Improved Variational Nodal Method
-
摘要: 共轭中子通量密度对于核安全和压水堆(PWR)中的探测器计算有着重要的意义,为了消除现有节块方法在处理由于控制棒移动带来的非均匀节块(包括非均匀的截面和不连续因子)时所造成的较大误差,本文提出一种改进的变分节块法(VNM)。确定了不同于前向方程的共轭节块方法的连续条件,不同于传统VNM在全局建立泛函,本文方法为每一个节块建立泛函;构建了含非均匀不连续因子的乘子项,以显式处理表面不连续的共轭中子通量密度;除共轭体中子通量密度、截面和表面分中子流密度外,将表面不连续因子展开为分段正交多项式来构造响应矩阵。含有非均匀节块的BEAVRS基准题数值结果证明,同传统VNM相比,改进的VNM可以将非均匀问题的有效共轭增殖系数和燃料区共轭中子通量密度偏差降低2个量级,有利于实现前向与共轭中子通量密度的高精度内积计算。Abstract: The adjoint neutron flux is of great significance for nuclear safety and detector calculation in pressurized water reactor (PWR). However, existing nodal methods would cause a big error due to heterogeneous nodes, including heterogeneous cross sections and discontinuity factors, which will appear frequently with the control rod moving. In this paper, an improved variational nodal method (VNM) is proposed to reduce the error. It determines the continuous conditions for adjoint nodal methods that are different from forward equation. Unlike traditional VNM, which establishes a functional method globally, this paper establishes a functional method for each node. It constructs a multiplier term with a heterogeneous discontinuity factor to explicitly deal with the adjoint neutron flux with surface discontinuity. Apart from the expansions of adjoint neutron flux, cross section and surface partial neutron current densities, the surface discontinuity factor (DF) is also expanded into pieces-wise orthogonal polynomials to construct the nodal response matrixes. The numerical results of the BEAVRS benchmark problem with heterogeneous nodes existing demonstrate that compared with the traditional VNM, the improved VNM can reduce the error by two orders of magnitude for the adjoint neutron flux in fuel area and the adjoint effective multiplication factor, which can help realize high accuracy calculation for the inner product of forward and adjoint neutron flux.
-
Key words:
- PWR /
- Heterogeneous node /
- VNM /
- Adjoint neutron flux
-
表 1 共轭有效增殖因子计算偏差
Table 1. Calculation Error of Adjoint Effective Multiplication Factor
验证算例 $k_{_{{\rm{eff}}}}^* $参考解 计算方法 $k_{_{{\rm{eff}}}}^*$偏差/pcm B62H 0.967789 HOM −30.5 HET 0.1 B62 0.967811 HOM −34.1 HET −0.1 表 2 节块平均共轭中子通量密度计算偏差
Table 2. Calculation Error of Nodal Average Adjoint Neutron Flux
验证算例 计算方法 能群 燃料区节块平均共轭中子
通量密度偏差/%最大值 均方根 B62H HOM 快群 7.229 0.595 热群 6.043 0.591 HET 快群 0.042 0.002 热群 0.028 0.001 B62 HOM 快群 7.979 0.658 热群 6.430 0.653 HET 快群 0.072 0.005 热群 0.050 0.004 -
[1] 李富,罗征培. 不连续因子理论下扩散差分方程的共轭计算[J]. 核动力工程,1998, 19(6): 485-489. [2] DULLA S, MUND E H, RAVETTO P. The quasi-static method revisited[J]. Progress in Nuclear Energy, 2008, 50(8): 908-920. doi: 10.1016/j.pnucene.2008.04.009 [3] CARON D, DULLA S, RAVETTO P. New aspects in the implementation of the quasi-static method for the solution of neutron diffusion problems in the framework of a nodal method[J]. Annals of Nuclear Energy, 2016, 87: 34-48. doi: 10.1016/j.anucene.2015.02.035 [4] 杜柄. 三维时空扩散中子动力学快速计算方法及外源驱动次临界系统中子学动力研究[D]. 西安: 西安交通大学, 2012. [5] ROHDE U, KLIEM S, GRUNDMANN U, et al. The reactor dynamics code DYN3D – models, validation and applications[J]. Progress in Nuclear Energy, 2016, 89: 170-190. doi: 10.1016/j.pnucene.2016.02.013 [6] 李富,王亚奇,罗征培. 节块格林函数法的微扰计算[J]. 核动力工程,1999, 20(2): 103-105. [7] DOWNAR T, XU Y, SEKER V. Theory manual for the PARCS kinetics core simulator module[M]. Michigan: University of Michig, 2010: 89-92. [8] PALMIOTTI G, LEWIS E E, CARRICO C B. VARIANT: VARIational Anisotropic Nodal Transport for multidimensional Cartesian and hexagonal geometry calculation: ANL-95/40[R]. Argonne, IL, USA: Argonne National Laboratory, 1995. [9] ZHANG T F, WU H C, CAO L Z, et al. An improved variational nodal method for the solution of the three-dimensional steady-state multi-group neutron transport equation[J]. Nuclear Engineering and Design, 2018, 337: 419-427. doi: 10.1016/j.nucengdes.2018.07.009 [10] ZHANG T F, ZHU L, WU H C, et al. Development and validation of high flux engineering test reactor fuel management platform HEFT[J]. Atomic Energy Science and Technology, 2013, 47(S1): 467-471. [11] KOBAYASHI K. Boundary condition for the adjoint flux corresponding to the flux discontinuity of the nodal diffusion equation in the perturbation theory[J]. Journal of Nuclear Science and Technology, 1998, 35(1): 20-25. doi: 10.1080/18811248.1998.9733815 [12] LIANG B N, WU H C, LI Y Z. Adaptive expansion order for diffusion variational nodal method[J]. Annals of Nuclear Energy, 2018, 117: 114-130. doi: 10.1016/j.anucene.2018.03.019 [13] 梁博宁,吴宏春,李云召. 控制棒尖齿效应中非均匀不连续因子处理技术[J]. 核动力工程,2020, 41(6): 31-35. doi: 10.13832/j.jnpe.2020.06.0031