随机中子动力学及其数值模拟研究进展

摘要/Abstract

摘要： 随机中子动力学是核动力设计和核反应堆安全中的重要课题，本文从随机中子动力学的基础概念和研究方法出发，介绍随机中子动力学研究的历史发展和研究现状.裂变中子与光子的多重性是反应堆零功率中子噪声主要来源，对中子涨落的方程描述及其求解，演化出零功率中子噪声与功率反应堆噪声的随机理论.随机中子动力学的重要应用包括反应性微观测量、功率反应堆噪声测量和分析、核临界漂移分析和核材料识别与检测等.在半个多世纪的研究中，以脉冲堆点火过程的脉冲爆发等待时间分布为代表的随机性，一直缺乏定量分析方法和工具.直到近几年，模拟随机中子动力学过程的广义半马尔科夫过程模拟方法取得了重要进展，很好地揭示了脉冲堆实验中子点火规律.最后讨论随机中子动力学研究中有待解决的研究课题.

关键词: 中子动力学, 反应堆噪声, 随机过程

Abstract: Stochastic neutron dynamics is important task in nuclear power plant design and nuclear reactor safety. In this paper, basic concepts and research methods for stochastic neutron kinetics are introduced and historical development and research status in the field are sketched out. Multiplicities of fission neutrons and photons are main sources of zero power reactor noise. Based on description equation of neutron fluctuation and its solution, stochastic theory of zero power and power reactor noise is evolved. Stochastic neutron dynamics are applied in important areas such as reactivity microscopic measurement, power reactor noise measurement and analysis, nuclear criticality excursion analysis, nuclear material detection and identification and so on. In the past half century, however, it lacks available method and tool to implement quantitative analysis of problems such as probability distribution of burst waiting time of neutron initiation in pulse reactor. In recent years, important progress is made in generalized semi-Markov process simulation method, which is applied to stochastic neutron kinetic process simulation. It reveals inherent law of neutron initiation experiments conducted in pulse reactor. At last, research topics to be solved in stochastic neutron dynamics are discussed.

Key words: neutron dynamics, reactor noise, stochastic process

中图分类号:

O571.33

杨俊云, 应阳君, 肖刚. 随机中子动力学及其数值模拟研究进展[J]. 计算物理, 2017, 34(2): 127-141.

YANG Junyun, YING Yangjun, XIAO Gang. Progress in Research on Stochastic Neutron Dynamics and Its Numerical Simulation[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2017, 34(2): 127-141.

参考文献

[1] GANG X. The simulation of probability distribution of the burst waiting time of neutron initiation[J]. Nucl Sci Eng, 2011, 169(1): 56-67.
[2] LIU Xiaobo, PENG Xianjue, LEI Jiarong, et al. Experimental study of neutron initiation probability at CFBR-II with pulsed neutrons[J]. Nucl Sci Eng, 2015, 181(1): 1-10.
[3] YANG Junyun, XIAO Gang, YING Yangjun. CriSAA: A simulation code for evaluation of criticality transients involving low source startup[J]. Nucl Sci Eng, 2016, 182(4):538-554.
[4] 杨俊云, 肖刚, 应阳君. 临界漂移分析[J].核动力工程, 2016, 37(2): 52-55.
[5] VALENTINE T E. Review of subcritical source-driven noise analysis measurements[R]. ORNL/TM-1999/288, Lockheed Martin Energy Research Corp, Oak Ridge National Laboratory, 1999.
[6] MIHALCZO J T, MULLENS J A, MATTINGLY J K, et al. Physical description of nuclear materials identification system (NMIS) signatures[J]. Nucl Instrum Meth A, 2000, 450(2): 531-555.
[7] MEDOWS J W. ²⁵²Cf fission neutron spectrum from 0.003 to 15.0 MeV[J]. Physical Review, 1967, 157(4): 1076.
[8] ETHVIGNOT T, DEVLIN M, DUARTE H, et al. Neutron multiplicity in the fission of ²³⁸U and ²³⁵U with neutrons up to 200 MeV[J]. Phys Rev Lett, 2005, 94(5): 052701.
[9] DIVEN B C, MARTIN H C, TASCHEK R F, et al. Multiplicities of fission neutrons[J]. Physical Review, 1956, 101(3): 1012-1015.
[10] TERRELL J. Distributions of fission neutron numbers[J]. Physical Review, 1957, 108(3): 783-789.
[11] 黄祖洽. 核反应堆动力学基础[M]. 北京:北京大学出版社, 2007.
[12] GWIN R, SPENCER R R, INGLE R W. Measurements of the energy dependence of prompt neutron emission from ²³³U, ²³⁵U, ²³⁹Pu, and ²⁴¹Pu for E_n=0.005 to 10 eV relative to emission from spontaneous fission of ²⁵²Cf[J]. Nucl Sci Eng, 1984, 87(4): 381-404.
[13] GWIN R, SPENCER R R, INGLE R W. Measurements of the energy dependence of prompt neutron emission from ²³³U, ²³⁵U, and ²³⁹Pu for E_n=0.000 5 to 10 MeV relative to emission from spontaneous fission of ²⁵²Cf[J]. Nucl Sci Eng, 1986, 94(4): 365-379.
[14] BOLDEMAN J W. The energy dependence of ν_p [J]. Physics of Neutron Emission in Fission, 1989: 21-45.
[15] MADLAND D G, NIX J R. New calculation of prompt fission neutron spectra and average prompt neutron multiplicities[J]. Nucl Sci Eng, 1982, 81(2): 213-271.
[16] FORT E, FREHAUT J, TELLIER H, et al. Evaluation of ν_p for ²³⁹Pu impact for applications of the fluctuations at low energy[J]. Nucl Sci Eng, 1988, 99(4): 375-389.
[17] ZUCKER M S, HOLDEN N E. Energy dependence of neutron multiplicity P_ν in fast-neutron-induced fission for ^235,238U and ²³⁹Pu[J]. BNL-38491, 1986.
[18] FREHAUT J. Neutron multiplicity distribution in fast neutron-induced fission[J]. Physics of Neutron Emission in Fissioin, 1989: 81-92.
[19] BOLDEMAN J W, HINES M G. Prompt neutron emission probabilities following spontaneous and thermal neutron fission[J]. Nucl Sci Eng, 1985, 91(1): 114-116.
[20] VALENTINE T E, MIHALCZO J T. MCNP-DSP: A neutron and gamma ray Monte Carlo calculation of source-driven noise-measured parameters[J]. Ann Nucl Energy, 1996, 23(16): 1271-1287.
[21] PLEASONTON F, FERGUSON R L, SCHMITT H W. Prompt gamma rays emitted in the thermal-neutron-induced fission of U²³⁵[J]. Phys Rev C, 1972, 6(3): 1023.
[22] FREHAUT J. Neutron gamma competition in fast fission[J]. Physics of Neutron Emission in Fission, 1989: 99.
[23] HOLDEN N E, ZUCKER M S. A reevaluation of the average prompt neutron emission multiplicity (nubar) values from fission of uranium and transuranium nuclides[M]. National Nuclear Data Center, Brookhaven National Laboratory, 1984.
[24] CULLEN D E. Sampling ENDL Watt fission spectra[R]. Lawrence Livermore National Lab, Livermore, CA (US), 2004.
[25] BRUNSON G S. Multiplicity and correlated energy of gamma rays emitted in the spontaneous fission of californium-252[D]. Dept of Mechanical and Industrial Engineering, University of Utah, 1982.
[26] NIFENECKER H, SIGNARBIEUX C, RIBRAG M, et al. Gamma-neutron competition in the de-excitation mechanism of the fission fragments of ²⁵²Cf[J]. Nucl Phys A, 1972, 189(2): 285-304.
[27] WAGEMANS C. The nuclear fission process[M]. CRC Press, 1991.
[28] 黄祖洽, 丁鄂江. 输运理论[M]. 北京:科学出版社, 1987.
[29] PAL L. On the theory of stochastic processes in nuclear reactors[J]. II Nuovo Cimento, 1958, 7: 25-42.
[30] BELL G I. On the stochastic theory of neutron transport[J]. Nucl Sci Eng, 1965, 21(3): 390-401.
[31] PAZSIT I, PAL L. Neutron fluctuations: A treatise on the physics of branching processes[M]. Elsevier, 2007.
[32] HURWITZ Jr H. Kinetics of low source reactor startups: Part I[J]. Nucl Sci Eng, 1963, 15(2): 166-186.
[33] HURWITZ Jr H. Kinetics of low source reactor startups: Part II[J]. Nucl Sci Eng, 1963, 15(2): 187-196.
[34] MUNOZ-COBO J L, PEREZ R B, VERDU G. Stochastic neutron transport theory: Neutron counting statistics in nuclear assemblies[J]. Nucl Sci Eng, 1987, 95(2): 83-105.
[35] HUMBERT P. Stochastic neutronics with PANDA deterministic code[J]. Nucl Math and Compu Sciences: A Century in Review, A Century Anew, 2003: 6-11.
[36] BAKER R S. Deterministic methods for time-dependent stochastic neutron transport[C].International Conference on Mathematics, Computational Methods, and Reactor Physics, (Saratoga Springs, NY), pp. CD-ROM, American Nuclear Society, 2009.
[37] GREENMAN G M, PROCASSINI R J, CLOUSE C J. A Monte Carlo method for calculating initiation probability[C]//Proceedings of the Joint International Topical Meeting on Mathematics & Computation and Supercomputing in Nuclear Applications, 2007: 16-19.
[38] MIHALCZO J T. Use of californium-252 as a randomly pulsed neutron source for prompt neutron decay measurements[J]. Nucl Sci Eng, 1970, 41(2): 296-298.
[39] FICARO E P, WEHE D K. KENO-NR: A Monte Carlo code for simulating ²⁵²Cf-source-driven noise measurements to determine subcriticality[J]. Nucl Sci Eng, 1994, 117(3): 158-176.
[40] VALENTINE T E, MIHALCZO J T. MCNP-DSP: A neutron and gamma ray monte carlo calculation of source-driven noise-measured parameters[J]. Ann Nucl Energy, 1996,23(16): 1271-1287.
[41] POZZI S A, PADOVANI E, MARSEGUERRA M. MCNP-PoliMi: A Monte-Carlo code for correlation measurements[J]. Nucl Instrum Meth A, 2003, 513(3): 550-558.
[42] 肖刚, 李天柁. 系统可靠性分析中的蒙特卡罗方法[M]. 北京: 科学出版社, 2003.
[43] 杨俊云, 肖刚, 应阳君. 点堆中子动力学的蒙特卡罗方法[J]. 核动力工程, 2015, 36(6): 4-9.
[44] 杨俊云, 应阳君, 肖刚. 随机中子动力学应用[D].2016 (未发表).
[45] 杨俊云, 肖刚, 胡泽华, 应阳君. 三维瞬态输运程序在中子点火中的应用[J].原子能科学技术, 2016,50(7)1320-1323.
[46] KUANG Z F, PAZSIT I. The general backward theory of neutron fluctuations in subcritical systems with multiple emission sources[J]. II Nuovo Cimento A, 1999, 112(10): 1067-1092.
[47] YAMAMOTO T. Higher order mode analyses in Feynman-α method[J]. Ann Nucl Energy, 2011, 38(6): 1231-1237.
[48] ORNDOFF J D. Prompt neutron periods of metal critical assemblies[J]. Nucl Sci Eng, 1957, 2(4): 450-460.
[49] ZOLOTUKHIN V G, MOGIL'NER A I. Distribution of the number of counts on a neutron detector placed in a reactor[J]. Atomic Energy, 1964, 15(1): 664-670.
[50] SZELESS A, RUBY L. The exact probability distribution of reactor neutron noise[J]. Nucl Sci Eng, 1971, 45(1): 7-13.
[51] BABALA D. Neutron counting statistics in nuclear reactors[R]. Institutt for Atomenergi, Kjeller (Norway), 1966.
[52] 白云, 杨波, 李金鸿, 等. 脉冲源法测量系统次临界度的进一步分析[J]. 原子能科学技术, 2010, 44(z1): 325-328.
[53] PAZSIT I, KITAMURA Y, WRIGHT J, et al. Calculation of the pulsed Feynman-alpha formulae and their experimental verification[J]. Ann Nucl Energy, 2005, 32(9): 986-1007.
[54] 宋凌莉, 周浩军, 金宇, 等. ²⁵²Cf随机脉冲源法测量深次临界瞬发中子衰减常数[J]. 原子能科学技术, 2007, 40(6): 714-717.
[55] KITAMURA Y, TAGUCHI K, MISAWA T, et al. Calculation of the stochastic pulsed Rossi-alpha formula and its experimental verification[J]. Prog Nucl Energ, 2006, 48(1): 37-50.
[56] DEGWEKER S B. Some variants of the Feynman alpha method in critical and accelerator driven sub critical systems[J]. Ann Nucl Energy, 2000, 27(14): 1245-1257.
[57] KITAMURA Y, PAZSIT I, WRIGHT J, et al. Calculation of the pulsed Feynman-and Rossi-alpha formulae with delayed neutrons[J]. Ann Nucl Energy, 2005, 32(7): 671-692.
[58] KOSALY G, VALKO J. Remarks on the use of well-known reactivity measuring techniques by the pulsed-source method[J]. Ann Nucl Energy, 1975, 2(6): 477-482.
[59] BAJPAI S, ROY T, RAY N, et al. Evaluation of spatial correction factors for BRAHMMA subcritical assembly[J]. Nucl Sci Eng, 2015, 181(3): 1-7.
[60] SPRIGGS G D, BUSCH R D, WILLIAMS J G. Two-region kinetic model for reflected reactors[J]. Ann Nucl Energy, 1997, 24(3): 205-250.
[61] MUNOZ-COBO J L, BERGLOF C, PENA J, et al. Feynman-α and Rossi-α formulas with spatial and modal effects[J]. Ann Nucl Energy, 2011, 38(2): 590-600.
[62] SPRIGGS G D. The reactor noise threshold[J]. Nucl Sci Eng, 1994, 116(1): 67-72.
[63] SPRIGGS G D, BUSCH R D, SAKURAI T, et al. The equivalent fundamental-mode source[J]. Ann Nucl Energy, 1999, 26(3): 237-264.
[64] WIMETT T F, WHITE R H, STRATTON W R, et al. Godiva II—An unmoderated pulse-irradiation reactor[J]. Nucl Sci Eng, 1960, 8(6): 691-708.
[65] HANSEN G E. Assembly of fissionable material in the presence of a weak neutron source[J]. Nucl Sci Eng, 1960, 8(6): 709-719.
[66] AUTHIER N, RICHARD B, HUMBERT P. Initiation of persistent fission chains in the fast burst reactor Caliban[J]. Nucl Sci Eng, 2014, 177(2): 169-183.
[67] BELL G I. Probability distribution of neutrons and precursors in a multiplying assembly[J]. Annals of Physics, 1963, 21(2): 243-283.
[68] BELL G I, ANDERSON W A, GALBRAITH D. Probability distribution of neutrons and precursors in multiplying medium[J]. Nucl Sci Eng, 1963, 16(1): 118-123.
[69] 杨俊云, 肖刚, 应阳君. 脉冲堆弱源点火实验的蒙特卡罗数值模拟[J].核科学与工程, 2015, 35(4): 36-44.
[70] MCLAUGHLIN T P, MONAHAN S P, PRUVOST N L, et al. A review of criticality accidents 2000 revision[R]. Los Alamos National Lab, NM (US), 2000.
[71] MCLAUGLIN T P. Process criticality accident likelihoods, consequences and emergency planning[J]. Nucl Energy, 1992, 31(2): 143-147.
[72] TUCK G. Simplified methods of estimating the results of accidental solution excursions[J]. Nucl Technol, 1974, 23(2): 177-188.
[73] MIYOSHI Y, YAMANE Y, OKUBO K, et al. Inter-code comparison exercise for criticality excursion analysis[C]. NEA, 2009.
[74] NAKAJIMA K, YAMAMOTO T, MIYOSHI Y. Modified quasi-steady-state method to evaluate the mean power profiles of nuclear excursions in fissile solution[J]. J Nucl Sci Technol, 2002, 39(11): 1162-1168.
[75] PAIN C C, DE OLIVEIRA C R E, GODDARD A J H, et al. Transient criticality in fissile solutions-compressibility effects[J]. Nucl Sci Eng, 2001, 138(1): 78-95.
[76] YAMAMOTO T, OKUNO H, NOMURA Y. Development of a new simulation code for evaluation of criticality transients Involving fissile solution boiling[M]. JAERI, 1998.
[77] HANKINS D E. Effect of reactivity addition rate and of weak neutron source on the fission yield of uranium solutions[J]. Nucl Sci Eng, 1966, 26(1): 110-116.
[78] RAMSEY S D, HUTCHENS G J. Deterministic and stochastic evaluation of criticality excursion power bursts[J]. Nucl Sci Eng, 2011, 168(3): 265-277.
[79] NAKAJIMA K, YAMANE Y, OGAWA K, et al. Tracy transient experiment databook 1: Pulse withdrawal experiment[M]. 2002.
[80] YAMANE Y, OGAWA K, AIIZAWA E, et al. TRACY Transient experiment databook 3: Ramp feed experiment[M]. JAERI, 2002.
[81] REILLY D, ENSSLIN N, SMITH Jr H, et al. Passive nondestructive assay of nuclear materials[R]. Nuclear Regulatory Commission, Washington, DC (United States); Office of Nuclear Regulatory Research; Los Alamos National Lab, NM (United States), 1991.
[82] GOZANI T. Active nondestructive assay of nuclear materials-principles and applications[R]. US Nuclear Regulatory Commission, NUREG/CR-0602, 1981.
[83] MATTINGLY J K, NEAL J S, MIHALCZO J T. Time-dependent coincidence method to measure plutonium mass and multiplication[R]. Y/LB-16121, 2001.
[84] LANGNER D G, STEWART J E, PICKRELL M M, et al. Application guide to neutron multiplicity counting[R]. Los Alamos National Laboratory, Los Alamos, NM, 1998.
[85] CRYE J M. Enrichment determination of uranium metal in shielded configurations without calibration standarts[D]. University of Tennessee, 2013.
[86] 李建胜, 叶岑明, 谢文雄, 等. 定时DT中子源驱动铀部件的相关函数测量[J]. 核技术, 2013, 36(9): 090204.
[87] XIE W, LI J, ZHU J. Uranium mass and neutron multiplication factor estimates from time-correlation coincidence counts[J]. Nucl Instrum Meth A, 2015, 797: 182-187.
[88] 刘晓波, 胡倩, 肖建国. 铀材料的中子多重性测量分析方法[J]. 核动力工程, 2009, 30(1): 74-77.
[89] PAZSIT I, POZZI S A. Calculation of gamma multiplicities in a multiplying sample for the assay of nuclear materials[J]. Nucl Instrum Meth A, 2005, 555(1): 340-346.