Application Research on VITAS—a General-purpose Neutron Transport Code
-
摘要: 为提高确定论全堆芯中子输运程序的适用性,开发了通用型中子输运程序VITAS。针对TAKEDA3基准题(矩形组件)、TAKEDA4基准题(六角形组件)、Dodds基准题(R-Z几何)和C5G7-TD5基准题(压水堆高保真计算)的验证结果表明,高阶的空间和角度基函数能够使结果稳定地向参考解渐进收敛,达到与多群蒙卡相当的计算精度水平。与参考解相比,TAKEDA3基准题有效增殖系数(keff)偏差低于60pcm(1pcm=10−5),控制棒价值偏差为−3pcm,中子通量密度分布均方根(RMS)偏差为1.03%;TAKEDA4基准题keff偏差低于20pcm,控制棒价值偏差为32pcm,中子通量密度分布RMS偏差为0.70%;Dodds基准题的功率最大偏差低于1%;C5G7-TD5基准题的功率偏差低于0.9%。本文研究表明VITAS有望成为一套精确求解中子输运问题的通用型计算工具。Abstract: In order to improve the applicability of deterministic whole-core neutron transport code, a general-purpose neutron transport code VITAS is developed. The verification results of TAKEDA3 benchmark problem (rectangular assembly), the TAKEDA4 benchmark problem (hexagonal assembly), the Dodds benchmark problem (R-Z geometry) and the C5G7-TD5 benchmark problem (PWR high fidelity calculation) show that higher-order spatial and angular basis functions can make the results converge asymptotically and steadily to the reference solution, reaching the calculation accuracy level equivalent to that of multi-group Monte Carlo method. Compared with the reference solution, the deviations are less than 60pcm (1pcm = 10−5) in effective multiplication coefficient (keff), −3pcm in control rod worth and 1.03% in neutron flux distribution root mean square (RMS) for TAKEDA3 benchmark problem; less than 20pcm in keff, 32pcm in control rod worth and 0.70% in neutron flux distribution RMS for TAKEDA4 benchmark problem; less than 1% (maximum) in power for Dodds benchmark problem; and less than 0.9% in power for C5G7-TD5 benchmark problem. The research in this paper shows that VITAS has the potential to become a general-purpose calculation tool for accurately solving the neutron transport problems.
-
表 1 计算基准题介绍
Table 1. Introduction to Calculation Benchmark Problems
基准题名称 问题类型 计算所用节块类型 TAKEDA3 三维、矩形、稳态 Cartesian TAKEDA4 三维、六角形、稳态 Hexagonal-z Dodds 三维、R-Z、瞬态 Triangular-z C5G7-TD5 三维、燃料棒精细描述、瞬态 Cartesian-FE 表 2 TAKEDA3基准题的keff偏差对比
Table 2. Comparison of keff Deviations for TAKEDA3 Benchmark Problem
计算对象 角度阶数 keff偏差/pcm 4/0阶 5/1阶 6/2阶 7/2阶 7/3阶 8/4阶 CR-In P1 2081 281 246 247 244 243 P3 1952 122 91 92 88 88 P5 1931 98 66 66 63 63 P7 1925 90 58 58 / / CR-Out P1 1637 216 199 199 198 197 P3 1547 99 85 85 83 83 P5 1531 80 65 65 63 63 P7 1526 74 59 59 / / 控制棒价值 P1 −551 −82 −62 −63 −61 −61 P3 −507 −30 −12 −13 −10 −10 P5 −501 −24 −5 −5 −4 −4 P7 −500 −21 −3 −3 / / “/”—计算超出了内存的限制,计算未启动,其余同 表 3 TAKEDA4基准题的keff偏差对比
Table 3. Comparison of keff Deviations for TAKEDA4 Benchmark Problem
计算对象 角度阶数 keff偏差/pcm 4/0阶 5/1阶 6/2阶 7/3阶 8/4阶 CRi P1 6783 821 518 544 560 P3 6336 353 52 59 67 P5 6278 289 −14 −17 / P7 6261 277 −29 / / CRh P1 6048 706 504 516 527 P3 5658 278 73 76 79 P5 5615 222 10 6 / P7 5603 210 −7 / / CRo P1 4841 460 376 377 381 P3 4549 163 70 70 69 P5 4519 120 17 12 / P7 4510 109 2 / / 控制棒价值 P1 −4268 −669 −353 −385 −402 P3 −3992 −318 −9 −18 −29 P5 −3950 −272 32 32 / P7 −3938 −266 39 / / 表 4 二维C5G7问题的计算结果对比
Table 4. Comparison of Calculation Results for Two-dimensional C5G7 Problem
角度阶数 keff偏差/pcm 棒功率偏差/% 最大值 RMS 32FE/2阶 96FE/4阶 144FE/6阶 32FE/2阶 96FE/4阶 144FE/6阶 32FE/2阶 96FE/4阶 144FE/6阶 P3_3 477 526 551 2.03 1.79 1.83 0.87 0.76 0.79 P5_3 261 313 325 1.48 1.27 1.28 0.57 0.48 0.48 P7_3 199 239 252 1.31 1.08 1.09 0.48 0.37 0.37 P9_3 132 160 168 1.14 0.88 0.88 0.39 0.27 0.27 P11_3 102 123 132 1.07 0.80 0.80 0.36 0.24 0.23 P13_3 71 86 93 1.00 0.71 0.70 0.32 0.21 0.20 P15_3 56 69 75 0.97 0.67 0.65 0.31 0.19 0.19 P17_3 41 51 56 0.93 0.62 0.60 0.29 0.18 0.18 P19_3 33 42 46 0.92 0.60 0.60 0.29 0.18 0.18 P21_3 24 30 33 0.89 0.61 0.64 0.28 0.18 0.17 P23_3 18 23 26 0.88 0.62 0.66 0.27 0.17 0.17 -
[1] 曹良志, 谢仲生, 李云召. 近代核反应堆物理分析[M]. 北京: 中国原子能出版社, 2017:4-7. [2] DOWNAR T, STAUDENMIER J. Theory Manual for the PARCS Neutronics core simulator: PARCS v2.6 U. S. NRC Core Neutronics Simulator, theory manual[Z]. Indiana, United States: Purdue University, 2004. [3] LAWRENCE R D. The DIF3D nodal neutronics option for two- and three-dimensional diffusion theory calculations in hexagonal geometry: ANL-83-1[R]. Argonne, Illinois, United States: Argonne National Lab, 1983. [4] VER PLANCK D M, COBB W R, BORLAND R S, et al. SIMULATE-E: a nodal core-analysis program for light-water reactors. Computer code user's manual: EPRI-NP-2792-CCM[R]. Framingham: Yankee Atomic Electric Co, 1983. [5] YANG W, WU H C, LI Y Z, et al. Development and verification of PWR-core fuel management calculation code system NECP-Bamboo: part Ⅱ Bamboo-Core[J]. Nuclear Engineering and Design, 2018, 337: 279-290. doi: 10.1016/j.nucengdes.2018.07.017 [6] ZHANG T F, LI Z P. Variational nodal methods for neutron transport: 40 years in review[J]. Nuclear Engineering and Technology, 2022, 54(9): 3181-3204. doi: 10.1016/j.net.2022.04.012 [7] 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 [8] ZHANG T F, XIONG J B, LIU L G, et al. Development and implementation of an integral variational nodal method to the hexagonal geometry nuclear reactors[J]. Annals of Nuclear Energy, 2019, 131: 210-220. doi: 10.1016/j.anucene.2019.03.031 [9] ZHUANG K, SHANG W, LI T, et al. Variational nodal method for three-dimensional multigroup neutron diffusion equation based on arbitrary triangular prism[J]. Annals of Nuclear Energy, 2021, 158: 108285. doi: 10.1016/j.anucene.2021.108285 [10] ZHANG T F, LEWIS E E, SMITH M A, et al. A variational nodal approach to 2D/1D pin resolved neutron transport for pressurized water reactors[J]. Nuclear Science and Engineering, 2017, 186(2): 120-133. doi: 10.1080/00295639.2016.1273023 [11] ZHANG T F, WANG Y P, LEWIS E E, et al. A three-dimensional variational nodal method for pin-resolved neutron transport analysis of pressurized water reactors[J]. Nuclear Science and Engineering, 2017, 188(2): 160-174. doi: 10.1080/00295639.2017.1350002 [12] XIAO W, LI X Y, LI P J, et al. High-fidelity multi-physics coupling study on advanced heat pipe reactor[J]. Computer Physics Communications, 2022, 270: 108152. doi: 10.1016/j.cpc.2021.108152 [13] GEUZAINE C, REMACLE J F. Gmsh: a 3-D finite element mesh generator with built-in pre- and post-processing facilities[J]. International Journal for Numerical Methods in Engineering, 2009, 79(11): 1309-1331. doi: 10.1002/nme.2579 [14] LOELIGER J, MCCULLOUGH M. Version control with Git: powerful tools and techniques for collaborative software development[M]. 2nd ed. 1005 Gravenstein Highway North, Sebastopol: O'Reilly Media, 2012:336. [15] BONNEY G E, KISSLING G E. Gram‐Schmidt orthogonalization of multinormal variates: applications in genetics[J]. Biometrical Journal, 1986, 28(4): 417-425. doi: 10.1002/bimj.4710280407 [16] TAKEDA T, IKEDA H. 3-D neutron transport benchmarks[J]. Journal of Nuclear Science and Technology, 1991, 28(7): 656-669. doi: 10.1080/18811248.1991.9731408 [17] ROMANO P K, FORGET B. The OpenMC Monte Carlo particle transport code[J]. Annals of Nuclear Energy, 2013, 51: 274-281. doi: 10.1016/j.anucene.2012.06.040 [18] SEUBERT A, SUREDA A, BADER J, et al. The 3-D time-dependent transport code TORT-TD and its coupling with the 3D thermal-hydraulic code ATTICA3D for HTGR applications[J]. Nuclear Engineering and Design, 2012, 251: 173-180. doi: 10.1016/j.nucengdes.2011.09.067 [19] GOLUOGLU S, BENTLEY C, DEMEGLIO R, et al. A deterministic method for transient, three-dimensional neutron transport: MOL. 19980504.0231[R]. Las Vegas: Yucca Mountain Project, 1998. [20] Argonne Code Center. Argonne Code Center: benchmark problem book: ANL-7416[R]. Argonne, Illinois, United States: American Nuclear Society, 1977. [21] HOU J, IVANOV K N, BOYARINOV V F, et al. OECD/NEA benchmark for time-dependent neutron transport calculations without spatial homogenization[J]. Nuclear Engineering and Design, 2017, 317: 177-189. doi: 10.1016/j.nucengdes.2017.02.008 [22] SHEN Q C, WANG Y R, JABAAY D, et al. Transient analysis of C5G7-TD benchmark with MPACT[J]. Annals of Nuclear Energy, 2019, 125: 107-120. doi: 10.1016/j.anucene.2018.10.049