连续能量共轭加权蒙特卡罗动态参数计算

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

计算物理 ›› 2018, Vol. 35 ›› Issue (1): 87-94.DOI: 10.19596/j.cnki.1001-246x.7570

连续能量共轭加权蒙特卡罗动态参数计算

彭钢

中国核动力研究设计院核反应堆系统设计技术国家级重点实验室, 四川成都 610041

收稿日期:2016-11-01 修回日期:2017-02-27 出版日期:2018-01-25 发布日期:2018-01-25
作者简介:彭钢(1973-),男,上海,高级工程师,硕士,主要从事反应堆物理研究工作,E-mail:penggang1973@vip.sina.com
基金资助:
中国核动力研究设计院探索基金（JJXM-TSQZ-2009-12-01）资助项目

Calculation of Kinetic Parameters with Continuous Energy Adjoint Weighted Monte Carlo Method

PENG Gang

National Key Laboratory of Reactor System Design Technology, Nuclear Power Institute of China, Chengdu 610041, China

Received:2016-11-01 Revised:2017-02-27 Online:2018-01-25 Published:2018-01-25

摘要/Abstract

摘要： 本文开展基于连续能量共轭加权蒙特卡罗的堆芯动态参数计算研究，这些参数主要包括缓发中子有效份额、瞬发中子代时间和瞬发中子衰减常数，在目前普遍采用的迭代裂变概率（IFP）的基础上，扩展原有IFP方法中共轭通量的选择，比较径迹长度估计、碰撞估计、吸收估计和tally记数估计的差别，及协方差和方差权重计算.同时，给出多个概率区间的动态参数分布范围，对迭代裂变概率法进行深入的研究，对比并给出下一代事件估计计数和不同中间代迭代裂变概率估计计数对计算结果的影响，从中得出合适的迭代代数.相关程序在计算完成后自动输出各种粒子分布状况统计.

关键词: 连续能量, 蒙特卡罗, 共轭权重, 动态参数

Abstract: A study of continuous energy adjoint weighted Monte Carlo method is carried out in several kinetic parameters, including effective delayed neutron fraction, prompt neutron generation time and prompt neutron decay constant, which is based on Iterated Fission Probability (IFP) method and selection of adjoint flux is extented. Difference among tracklength, collision, absorption and tally estimator are compared. Covariance and variance weight are also adopted based on these estimators. Meanwhile, different distribution domain of kinetic parameters are provided. IFP method is studied in more detail. Next event estimate and IFP method with different latent generation number are compared and influence on result are obtained. Fitness number of latent generation number is obtained. Particle kind and distribution information are automatically listed in output file.

Key words: continuous energy, Monte Carlo, adjoint weighted, kinetic parameters

中图分类号:

TL327

彭钢. 连续能量共轭加权蒙特卡罗动态参数计算[J]. 计算物理, 2018, 35(1): 87-94.

PENG Gang. Calculation of Kinetic Parameters with Continuous Energy Adjoint Weighted Monte Carlo Method[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2018, 35(1): 87-94.

参考文献

[1] KEEPIN G R. Physics of nuclear kinetics[M]. Addison-Wesley,1965.
[2] X-5 MONTE CARLO TEAM. MCNP-A General N-particle transport code[CP/OL]. Version 5.LA-UR-03-1987, Los Alamos National Laboratory[Apr 2003]. https://www.lanl.gov.
[3] BRETSCHER E E. Evalution of reactor kinetic parameters without the need for perturbation codes[C]. Proc Int Mtg Reduced Enrichment for Research and Test Reactors,JacksonHole,Wyoming,Oct 7-10,1997.
[4] STANKOVSKIY E Y. Calculation of the effective delayed neutron fraction with MCNP[J]. Trans Am Nucl Soc, 2009, 100:634-635.
[5] MEULEKAMP R K, STEVEN C, van der MARCK. Calculating the effective delayed neutron fraction with Monte Carlo[J]. Nucl Sci Eng,2006, 152:142-148.
[6] van der MARCK SC, MEULEKAMP R K. Calculating the effective delayed neutron fraction using Monte Carlo techniques[C]. PHYSOR 2004,The Physics of Fuel Cycles and Advanced Nuclear Systems:Global Developments, Chicago, Illinois, April 25-29, 2004.
[7] YASUSHI Nauchi,TAKAMORI Kameyama. Proposal of direct calculation of kinetic parameters β_eff and Λ based on continuous energy Monte Carlo method[J]. J Nucl Sci Techn,2005, 42:503-514.
[8] BRIAN C, Kiedrowski, et al.Calculating kinetics parameters and reactivity changes with continuous energy Monte Carlo[C]. Proc PHYSOR2010, Pittsburgh,Pennsylvania,USA,May 9-14,2010.
[9] JUDITH F, Briesmeister. MCNP^TM-A general Monte Carlo N-particle transporrt code[CP/OL]. Version 4C, LA-13709-M,March, 2000. https://mcnp.lanl.gov.
[10] BRIAN Kiedrowski. Adjoint-weighted kinetics parameters with continuous energy Monte Carlo[J]. Trans Am Nucl Soc,2009, 100:297-299.
[11] YASUNOBU Nagaya,KEN Nakajima. Approximate estimation of effective delayed neutron fraction with correated sampling method[C]. Proc PHYSOR2006,Interlaken,Switzerland,Sep 14-19,2008.
[12] STEPHANIE C Frankle. A suite of criticality benchmarks for validating nuclear data[EB/OL]. LA-13594,April 1999, https://www.lanl.gov.
[13] 徐传效,等.高通量工程试验堆(HFETR)运行十年论文集[M]. 四川科学技术出版社,1990.
[14] RAUCK S. Science V2 nuclear code package-qualification report[EB/OL]. FRAMATOME ANP,Ver.A,2004. https://www.areva.com.
[15] HALPERIN M. Almost linearly optimum combination of unbiased estimates[J]. Amer Stat Ass J, 1961, 56:36-43.
[16] SCHEFFH. On solutions of the Behrens Fisher problem based on the t-distribution[J]. Annals of Mathematical Statistics,1943,14:35-44.

连续能量共轭加权蒙特卡罗动态参数计算

Calculation of Kinetic Parameters with Continuous Energy Adjoint Weighted Monte Carlo Method

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	王敏, 申玉清, 陈震宇, 徐鹏. 随机多孔介质的蒙特卡罗重构与渗流特性模拟[J]. 计算物理, 2021, 38(5): 623-630.
[2]	杨波, 衷斌, 徐琪, 成立, 潘流俊, 沈华韵. (α, n)反应中子产额与能谱计算方法[J]. 计算物理, 2021, 38(4): 393-400.
[3]	邓楚楚, 张鑫源, 陈水源, 吴宇晗, 叶晴莹, 黄志高. 小缺陷铁纳米环的磁特性[J]. 计算物理, 2020, 37(4): 473-478.
[4]	任娟, 张宁超, 刘萍萍. 碳气凝胶储氢性能的理论研究[J]. 计算物理, 2019, 36(6): 749-756.
[5]	吴言, 刘思琪, 李胜强. 芯片表面的冷极性分子静电阱与微阱阵列[J]. 计算物理, 2019, 36(4): 483-490.
[6]	王一, 潘流俊, 王瑞宏. 求解α本征值的k-α迭代方法误差估计[J]. 计算物理, 2019, 36(1): 39-46.
[7]	刘聪, 张斌, 张亮, 郑君萧, 陈义学. S_N共轭函数用于蒙特卡罗粒子输运自动减方差的研究[J]. 计算物理, 2018, 35(5): 535-544.
[8]	马先林, 周德胜. 基于水平集的两步MCMC方法对河道形态的识别[J]. 计算物理, 2018, 35(3): 321-329.
[9]	徐小东, 周彬, 郭金川, 易明皓, 肖飞虎. 微结构X射线源中电子与阳极作用的蒙特卡罗模拟[J]. 计算物理, 2018, 35(1): 95-102.
[10]	李胜强. 一种针对弱场搜寻态冷极性分子的可控静电表面囚禁方案[J]. 计算物理, 2017, 34(6): 731-739.
[11]	李刚, 张宝印, 邓力, 上官丹骅, 李瑞, 马彦, 付元光, 胡小利. 蒙特卡罗区域分解并行计算中确保串并行结果一致的伪随机数应用[J]. 计算物理, 2017, 34(1): 67-72.
[12]	杨师俨, 何曼丽, 蒋丹枫, 王国辉, 苑超, 戴耀东. 反应堆²³⁵U裂变源辐射防护材料的优化设计[J]. 计算物理, 2017, 34(1): 73-81.
[13]	上官丹骅, 许海燕. 非定常输运模拟中基于粒子标识分类的源偏倚算法[J]. 计算物理, 2016, 33(6): 639-644.
[14]	刘雄国, 邓力, 胡泽华, 李瑞, 付元光, 李刚. 基于VENUS-III国际基准模型的JMCT程序验证[J]. 计算物理, 2016, 33(5): 570-580.
[15]	付元光, 郑建华, 上官丹骅, 李瑞, 李刚, 马彦, 邓力. 蒙特卡罗粒子输运软件JMCT的网格计数功能设计与实现[J]. 计算物理, 2016, 33(5): 581-586.