Error Estimation of k-α Iteration Method for α Eigenvalue Problems

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

Abstract

Abstract: We analyze k-α iteration method for α eigenvalue problems, and point out that parameter in k-α iteration method affects error estimation of alpha eigenvalues, while error estimation has no relationship with the parameter. If the parameter isn't given properly, the estimated error deviates the real error estimation. Estimated error given by MCNP4C may deviates from real error estimation due to improper parameter. A parameter selecting method called average neutron fission lifetime method is proposed to make estimated error and real error estimation equal strictly in mathematics. It was tested with numerical experiments.

Key words: Monte Carlo, α eigenvalue, k-α iteration, error estimation

CLC Number:

TL329

WANG Yi, PAN Liujun, WANG Ruihong. Error Estimation of k-α Iteration Method for α Eigenvalue Problems[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2019, 36(1): 39-46.

References

[1] DENG Li, LI Gang. A summarization on Monte Carlo simulation in particle transport[J].Chinese Journal of Computational Physics, 2010, 27(6):791-798.
[2] 谢仲生,邓力. 中子输运理论数值计算方法[M]. 西安:西北工业大学出版社,2005.
[3] 杜书华. 输运问题的计算机模拟[M]. 湖南:湖南科技出版社,1988.
[4] 裴鹿成,张孝泽. 蒙特卡罗方法及其在粒子输运问题中的应用[M]. 北京:科学出版社,1980.
[5] 李树,田东风,邓力. 中子时间常数的Monte Carlo计算方法[J]. 清华大学学报(自然科学版),2007,47:38-46.
[6] WANG Ruihong, DENG Li, XU Haiyan, et al. Time forced collision methods for multiplication rate calculation[J]. Chinese Journal of Computational Physics,2006, 23(1):1-9.
[7] 王瑞宏,黄正丰,江松,等. 计算α本征值的蒙特卡罗期望估计方法[J]. 核科学与工程,2008,28(3):204-209.
[8] SHIM Hyung Jin, JANG Sang Hoon, KANG Soo Min. Monte Carlo alpha iteration algorithm for a subcritical system analysis[J]. Science and Technology of Nuclear Installations, 2015, 2015(1):1-7.
[9] BROCKWAY D, SORAN P, WHALEN P L. Monte Carlo eigenvalue calculation[R]. La-UR-85-1224. USA:Los Alamos Scientific Laboratory Report, 1985.
[10] DU Jinfeng. Computation of α eigenvalue on neutron transport equation[J]. Chinese Journal of Computational Physics, 2003,20(6):509-513.
[11] YE Tao, CHEN Chaobin, SUN Weili, et al. Prompt time constants of a reflected reactor[C]. Proceedings of the Symposium on Nuclear Data, Japan, 2007.
[12] YAMAMOTO Toshihiro. Higher order mode eigenvalue calculation by Monte Carlo power iteration[J].Nuclear Science and Technology, 2015, 2:826-835.
[13] ZOIA Andrea, BRUN Emeric, MALVAGI Fausto. Alpha eigenvalue calculations with TRIPOLI-4[J].Annals of Nuclear Energy, 2014, 63:276-284.
[14] BRIESMEITER J F L. MCNP-A general Monte Carlo n-particle transport code, version 4C[R]. La-13709-M. USA:Los Alamos Scientific Laboratory Report, 2000.
[15] 邱有恒,李茂生,邓力. 次临界系统中α本征值与k本征值的关系[J]. 清华大学学报(自然科学版),2007,47:967-971.