基于MC的辐射场快速计算程序设计与验证
Program Design and Verification of Rapid Calculation of Radiation Field Based on MC
DOI: 10.12677/NST.2022.103019, PDF,  被引量   
作者: 谢明亮*:中核武汉核电运行技术股份有限公司,湖北 武汉;海军工程大学核科学技术学院,湖北 武汉;彭 波, 魏 巍, 李 青, 谢政权:中核武汉核电运行技术股份有限公司,湖北 武汉;陈玉清, 于 雷:海军工程大学核科学技术学院,湖北 武汉;沈华韵, ZhongBin :北京应用物理与计算数学研究所,北京
关键词: 蒙特卡罗MCNP5NPTS能谱Monte Carlo MCNP5 NPTS Spectrum
摘要: 针对当前核电厂辐射场剂量计算需求,基于蒙特卡罗方法,跟踪模拟粒子输运与碰撞过程,对辐射场快速计算方法进行研究,采用分裂、轮盘赌、偏倚抽样等三类减方差技巧加速收敛;同时采用MPI进行并行加速,设计并验证辐射场计算程序NPTS;选取15个MCNP5自带的临界基准题,较为全面地考察该程序在临界计算方面的正确性,选取4个固定源计算例题,验证固定源问题处理的正确性,同时对并行效率进行验证分析。结果表明:选取的临界基准题Keff相对偏差小于5‰,4个固定源基准题模型的中子伽马能谱与MCNP5程序基本一致,验证了粒子输运的精准性;而且计算效率优于MCNP5程序,计算选取模型百核并行效率高于95%,缩短了程序计算时间,为建立核与辐射防护剂量快速计算体系提供参考和依据。
Abstract: For the current demand of radiation field dose calculation of nuclear power plant, based on the Monte Carlo method, the rapid calculation method of radiation field is studied by tracking and simulating the particle transport and collision process, and three kinds of subtraction variance reduction such as splitting, roulette and bias sampling are adopted. At the same time, MPI was used for parallel acceleration, and the radiation field calculation program NPTS was designed and verified. 15 MCNP5 critical reference questions were selected to examine the correctness of the program in the critical calculation more comprehensively, 4 fixed source calculation questions were selected to verify the correctness of the fixed source problem, and verify the parallel efficiency. The results showed that the relative deviation of Keff was less than 5‰, and the neutron gamma spectrum of 4 fixed source models was basically the same as that of MCNP5, which verified the accuracy of particle transport. The computational efficiency is better than that of the MCNP5 program, and the parallel efficiency of hundreds cores for model calculation is higher than 95%, which significantly shortens calculating time of the program and provides a reference and basis for the establishment of a re-al-time three-dimensional radiation field calculation system.
文章引用:谢明亮, 彭波, 魏巍, 李青, 谢政权, 陈玉清, 于雷, 沈华韵, ZhongBin. 基于MC的辐射场快速计算程序设计与验证[J]. 核科学与技术, 2022, 10(3): 183-194. https://doi.org/10.12677/NST.2022.103019

参考文献

[1] Shang, X.T., Yu, G.L., Song, J. and Wang, K. (2018) A Direct Calculation Method for Subcritical Multiplication Factor in Reactor Monte Carlo Code RMC. Annals of Nuclear Energy, 118, 81-89. [Google Scholar] [CrossRef
[2] Wang, Z.Y., Wu, B., Hao, L.J., Liu, H.F. and Song, J. (2018) Validation of SuperMC with BEAVRS Benchmark at Hot Zero Power Condition. Annals of Nuclear Energy, 111, 709-714. [Google Scholar] [CrossRef
[3] Deng, L. (2018) JMCT Monte Carlo Code with Ca-pability of Integrating Nuclear System Feedback. Proceedings of 2018 2nd International Conference on Applied Mathe-matics, Modelling and Statistics Application (AMMSA 2018), Sanya, 27-28 May 2018, 58-64. [Google Scholar] [CrossRef
[4] Xie, M.-L., Xie, F., Shan, F.-C., et al. (2019) Estimation of the Dose Rate of Spent Fuel Related Components of Lingao Nuclear Power Plant Using ORIGEN2 and MCNP5. Wuhan University Journal of Natural Sciences, 24, 64-70. [Google Scholar] [CrossRef
[5] 谢明亮, 于雷, 陈玉清. AP1000核反应堆控制棒价值特性的MC模拟[J]. 兵器装备工程学报, 2016, 37(3): 121-125.
[6] Lee, Y.-K. (2018) Radiation Shielding Calculations for a 3D ITER Benchmark Model Using TRIPOLI-4® Monte Carlo Code. Fusion Engineering and Design, 136, 602-607. [Google Scholar] [CrossRef
[7] Ellis, J.A., Evans, T.M., Hamilton, S.P., et al. (2019) Opti-mization of Processor Allocation for Domain Decomposed Monte Carlo Calculations. Parallel Computing, 87, 1-12.
[8] Tracheva, N. and Ukhinov, S. (2019) A New Monte Carlo Method for Estimation of Time Asymptotic Pa-rameters of Polarized Radiation. Mathematics and Computers in Simulation, 161, 84-92. [Google Scholar] [CrossRef
[9] Amato, E., Italiano, A. and Auditore, L. (2018) Radiation Pro-tection from External Exposure to Radionuclides: A Monte Carlo Data Handbook. Physica Medica, 46, 160-167. [Google Scholar] [CrossRef] [PubMed]
[10] Sevastyanov, V.D. and Trykov, L.A. (2008) Neutron and γ-Radiation Field Characteristics for 14 MeV Neutron Generators Used with a Stilbene-Crystal Spectrometer. Measure-ment Techniques, 5, 541-549. [Google Scholar] [CrossRef
[11] Cunsolo, S., Baillis, D. and Bianco, N. (2019) Improved Monte Carlo Methods for Computational Modelling of Thermal Radiation Applied to Porous Cellular Materials. International Journal of Thermal Sciences, 137, 161-179. [Google Scholar] [CrossRef
[12] Demir, N., Kuluozturk, Z.N. and Akkurt, I. (2019) FLUKA Monte Carlo Calculations for Angular Distribution of Bremsstrahlung Photons from Thin Targets. Nuclear Instruments and Methods in Physics Research Section B, 443, 19-24. [Google Scholar] [CrossRef
[13] Dupont, C., Baert, G. and Mordon, S. (2019) Parallelized Mon-te-Carlo Dosimetry Using Graphics Processing Units to Model Cylindrical Diffusers Used in Photodynamic Therapy: From Implementation to Validation. Photodiagnosis and Photodynamic Therapy, 26, 351-360. [Google Scholar] [CrossRef] [PubMed]
[14] 聂星辰, 李佳, 赵平辉, 等. 深穿透屏蔽计算中MCNP减方差技巧应用及比较[J]. 核电子学与探测技术, 2016, 36(7): 729-733.
[15] 张磊, 贾铭椿, 龚军军, 等. 几何分裂与轮盘赌减方差技巧中分裂比对计算效率影响研究[J]. 核电子学与探测技术, 2017, 37(8): 815-818.
[16] Mol, A.C.A., Pereira, C.M.N.A. and Freitas, V.G.G. (2011) Radiation Dose Rate Map Interpolation in Nuclear Plants Using Neural Networks and Virtual Reality Techniques. Annals of Nuclear Energy, 38, 705-712. [Google Scholar] [CrossRef
[17] Kahani, M., Kamali-Asl, A. and Tabrizi, S.H. (2019) Proposi-tion of a Practical Protocol for Obtaining a Valid Radiology Image Using Radiography Tally of MCNPX Monte Carlo Code. Applied Radiation and Isotopes, 149, 114-122. [Google Scholar] [CrossRef] [PubMed]