求解α本征值的k-α迭代方法误差估计

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

计算物理 ›› 2019, Vol. 36 ›› Issue (1): 39-46.DOI: 10.19596/j.cnki.1001-246x.7794

求解α本征值的k-α迭代方法误差估计

王一¹, 潘流俊², 王瑞宏²

1. 中国工程物理研究院研究生院, 北京 100088;
2. 北京应用物理与计算数学研究所, 北京 100094

收稿日期:2017-11-03 修回日期:2018-01-11 出版日期:2019-01-25 发布日期:2019-01-25
通讯作者: 王瑞宏,E-mail:wang_ruihong@iapcm.ac.cn
作者简介:王一(1994-),研究生,研究方向:蒙特卡罗粒子输运方法及应用
基金资助:
国家重点基础研究发展计划（2011CB309705）及国家自然科学基金（11771051）资助项目

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

WANG Yi¹, PAN Liujun², WANG Ruihong²

1. Graduate School of China Academy of Engineering Physics, Beijing 100088, China;
2. Institute of Applied Physics and Computational Mathematics of Beijing, Beijing 100094, China

Received:2017-11-03 Revised:2018-01-11 Online:2019-01-25 Published:2019-01-25

摘要/Abstract

摘要： 分析求解α本征值的k-α迭代方法，指出迭代法中参数的选择影响使用标准差的平均值估计出来的α误差.如果参数选取不当，实际估计误差可能和真实估计误差相差很大.MCNP4C程序中使用的参数给出的误差估计通常和真实估计误差有一定的偏差.本文给出一个新的参数选择办法，能够让实际估计误差几乎等于真实估计误差，并用数值实验进行验证.

关键词: 蒙特卡罗α本征值, k-α迭代, 误差估计

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

中图分类号:

TL329

王一, 潘流俊, 王瑞宏. 求解α本征值的k-α迭代方法误差估计[J]. 计算物理, 2019, 36(1): 39-46.

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.

参考文献

[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.

求解α本征值的k-α迭代方法误差估计

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

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 13

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	马岩, 胡健伟. 计算圆柱绕流拉力系数的自适应迎风有限元方法[J]. 计算物理, 2009, 26(3): 378-388.
[2]	聂存云, 舒适, 盛志强. 非结构四边形网格下的一类保对称有限体元格式[J]. 计算物理, 2009, 26(2): 175-183.
[3]	郭双冰, 张志跃. 一类非线性发展方程的特征中心差分法[J]. 计算物理, 2007, 24(6): 637-646.
[4]	张志娟. 非饱和水流问题的迎风有限体积元方法及数值模拟[J]. 计算物理, 2007, 24(6): 725-732.
[5]	王兆清, 李淑萍. 多边形单元平均值插值的误差估计及应用[J]. 计算物理, 2007, 24(2): 217-221.
[6]	李长峰, 袁益让. 抛物方程显隐修正迎风区域分裂差分法[J]. 计算物理, 2007, 24(2): 239-246.
[7]	刘伟, 袁益让. 二维半导体问题在复合网格上的有限差分格式[J]. 计算物理, 2006, 23(6): 721-730.
[8]	袁益让, 杜宁, 韩玉笈, 冯国祥, 杨成顺, 左东华. 油资源运移聚集并行数值模拟和分析[J]. 计算物理, 2005, 22(1): 27-37.
[9]	徐茂林, 邱钧, 范惠荣, 张兆田, 李兴东. 用一种新滤波函数作CT图像局部重建[J]. 计算物理, 2004, 21(3): 362-368.
[10]	王同科. 一维非线性对流占优扩散方程的变网格特征差分方法[J]. 计算物理, 2003, 20(6): 493-497.
[11]	张瑞红, 王鹏林, 郭长海, 张敬宇. 铸造系统热传导有限元法自适应分析[J]. 计算物理, 2000, 17(S1): 59-64.
[12]	蒋勇, 戴煜. 任意阶导函数的数值逼近[J]. 计算物理, 1993, 10(2): 191-196.
[13]	吴新元, 吴宏伟. 一个新的高精度二重数值积分公式[J]. 计算物理, 1991, 8(4): 437-441.