A Three-dimensional Neutron Diffusion Nodal Solver Based on 2D Generic Transverse Integrate Fused with 1D Semi-analytical Nodal Method

doi:10.19596/j.cnki.1001-246x.8377

Abstract

Abstract:

The nodal expansion method (NEM) based on transverse integrated technique is widely used in liquid water reactor core design. It uses an empirical transverse leakage polynomial approximation. Furthermore, it requires pin power reconstruction model for pin power within the assembly. A generic transverse integrate proposed in this work avoids the above shortages. Since the core heterogeneous usually happens in radial direction, the 3D neutron diffusion equation is solved with 2D generic transverse integrate fused with traditional 1D semi-analytical nodal method. With IAEA 2D/3D benchmark problem, the method is demonstrated correct and feasible for core neutron diffusion calculations.

Key words: generic transverse integrate, nodal method, neutron diffusion, semi-analytical

Zhiyong LI. A Three-dimensional Neutron Diffusion Nodal Solver Based on 2D Generic Transverse Integrate Fused with 1D Semi-analytical Nodal Method[J]. Chinese Journal of Computational Physics, 2022, 39(2): 153-158.

Figures/Tables 5

References 7

1	LIAO C K. Study on numerical method for three dimensional nodal space-time neutron kinetics equations and coupled neutronic/thermal-hydraulic core transient analysis [D]. Xi′an: Xi′an Jiaotong University, 2002.
2	ZIMIN V G, NINOKATA H, POGOSBEKYAN L R. Polynomial and semi-analytic nodal methods for nonlinear iteration procedure [C]. Proc Int Conf Phys Nuclear Science and Technology (PHYSOR-98), Long Island, NY, October 5-8, 1998: 994-1002.
3	HU Y M , ZHAO X F . A new nodal Green's function method for cylindrical geometry[J]. Chinese Journal of Computational Physics, 1997, 14 (4): 429- 430.
4	CHENG P D . Partial flux expansion techniques for some advanced nodal methods[J]. Chinese Journal of Computational Physics, 2001, 18 (5): 451- 456.
5	LI Z Y . Full core pin by pin neutron diffusion nodal method based on CMFD acceleration[J]. Chinese Journal of Computational Physics, 2019, 36 (3): 317- 322.
6	DOWNAR T, LEE D, XU Y, et al. PARCS v2.6 US NRC core neutronics simulator theory manual [R]. US: School of Nuclear Engineering Purdue University, 2004.
7	CHAO Y A , PENG L H , TANG C T . A new and rigorous SPN theory, Part Ⅳ: Numerical qualification of GSP30 and the generalized transverse integration nodal method[J]. Annals of Nuclear Energy, 2020, 149, 1- 16.

	方法	k_eff	偏差/pcm	平均偏差ε_avg^a/%	最大偏差ε_max^b/%	计算时间/s
IAEA-2D	参考解	1.029 58
IAEA-2D	本文	1.029 49	-9	0.16	-0.64	0.15
IAEA-3D	参考解	1.029 03
IAEA-3D	本文	1.028 99	-4	0.11	-0.56	7.37

	方法	k_eff	偏差/pcm	平均偏差ε_avg^a/%	最大偏差ε_max^b/%	计算时间/s
IAEA-2D	参考解	1.029 58
IAEA-2D	本文	1.029 49	-9	0.16	-0.64	0.15
IAEA-3D	参考解	1.029 03
IAEA-3D	本文	1.028 99	-4	0.11	-0.56	7.37

	方法	k_eff	偏差/pcm	平均偏差ε_avg^a/%	最大偏差ε_max^b/%	计算时间/s
BIBLIS-2D	参考解	1.025 110
BIBLIS-2D	本文	1.025 095	-1.5	0.30	-0.76	0.21
LMW-3D	参考解	0.999 660
LMW-3D	本文	0.999 719	+5.9	0.03	-0.09	0.87

	方法	k_eff	偏差/pcm	平均偏差ε_avg^a/%	最大偏差ε_max^b/%	计算时间/s
BIBLIS-2D	参考解	1.025 110
BIBLIS-2D	本文	1.025 095	-1.5	0.30	-0.76	0.21
LMW-3D	参考解	0.999 660
LMW-3D	本文	0.999 719	+5.9	0.03	-0.09	0.87

[1]	LI Zhiyong. Full Core Pin by Pin Neutron Diffusion Nodal Method Based on CMFD Acceleration [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2019, 36(3): 317-322.
[2]	CHEN Gong, WANG Yizheng, WANG Ye, ZHANG Chunyu. Reduced Basis Finite Element Method for Fast Solution of Parameterized Partial Differential Equations [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2018, 35(5): 515-524.
[3]	ZHANG Yifei, ZHANG Zhiquan, CHEN Yuefeng. Semi-analytical Integration for Coefficients of Charge on Plane Elements in Surface Charge Method [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2010, 27(4): 567-572.
[4]	PENG Jianshe, LIU Yan, YANG Jie. Convolution Semi-analytical Differential Quadrature Method for Structural Dynamic Response [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2009, 26(3): 409-414.
[5]	GUO Yonghui, DUAN Yaoyong, KUAI Bin. Hydrogen-like Wave Functions in a Non-relativistic Average Atom Model for Plasmas [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2007, 24(4): 487-493.
[6]	YANG Yin-ling, GUO Xiu-lan, MA Yi-chen, ZHANG Zhi-bin. MIXED AND DISCONTINUOUS FINITE ELEMENT METHOD USED TO NEUTRON DIFFUSION EQUATION [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2001, 18(6): 563-568.
[7]	FU Xue-dong, LIU Cheng-an. A NUMERICAL STABILIZATION METHOD FOR NONLINEAR ANALYTIC NODAL METHOD [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2001, 18(2): 111-114.
[8]	Hu Yongming, Zhao Xianfeng. A NEW NODAL GREEN'S FUNCTION METHOD FOR CYLINDRICAL GEOMETRY [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 1997, 14(S1): 429-430.
[9]	Xu Yunlin, Jing Xingqing, Li Junli. THE RESEARCH OF IMPROVED QUASI-STATIC METHOD IN REACTOR KINETICS [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 1997, 14(S1): 524-526.
[10]	Chen Chao, Zhu Guangning, Wang Kan. RECONSTRUCTION OF PIN POWER IN HEXAGONAL GEOMETRY [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 1997, 14(S1): 617-618,616.
[11]	Peng Jianshe, Zhang Jingyu,. Yang Jie. A SEMI-ANALYTICAL METHOD BASED ON CONVOLUTIONTYPE VARIATIONAL PRINCIPLE TO SOLVE THE TRANSIENT HEAT TRANSFER PROBLEMS [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 1996, 13(2): 237-242.
[12]	Dong Hefeng, Wang Kexie, Dong Qingde. NUMERICAL SIMULATION OF ACOUSTIC FIELD IN A FLUID-FILLED BOREHOLE SURROUNDED BY HORIZONTICAL THIN-LAYERED MEDIA WITH THE SEMI-ANALYTICAL METHOD [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 1995, 12(4): 449-456.
[13]	Peng Jianshe, Zhang Jingyu. A SEMI-ANALYTICAL METHOD BASED ON GURTIN VARIATIONAL PRINCIPLE TO SOLVE 1-D INITIAL-VALUE PROBLEMS OF DYNAMICS UNDER ALTERNATE FORCES [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 1995, 12(4): 571-575.
[14]	Cheng Pingdong, Zhao Huimin, Shao Tianwei, Zhou Qianfu, Shen Sen. SNENM/RTSP:A THREE-DIMENSIONAL NODAL AND SUBCHANNEL ANALYSIS CODE PACKAGE FOR THE CORE TRANSIENT IN PWPS [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 1992, 9(S1): 579-582.
[15]	Zhang Zongyao. THE STATUS AND PROSPECT OF REACTOR ANALYSIS METHODS FOR LWR [J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 1992, 9(S1): 593-596.