Calculation of the Decommissioning Radiation Field of Nuclear
Total Page:16
File Type:pdf, Size:1020Kb
Calculation of the Decommissioning Radiation Field of Nuclear Power Plants Based on the Coupling of MC and Point Kernel Integration GUO Yu-fei1,2, LIU Yan-fang1,2, WANG Li1,2, LIU Shui-qing1, ZHANG Hang-zhou1,2 1Nuclear Power Institute of China,Chengdu,610005,China 2Sichuan Engineering Laboratory for Nuclear Facilities Decommissioning and Radwaste Management,Chengdu,610005,China Table 3.4 Calculation accuracy of the coupling method for different Introduction Single Source Calculation Model grid numbers in the single source model • In the three-dimensional simulation system for FMC-PK/FMC the decommissioning of nuclear power plants, ① ② Point G: 32 G: 64 the radiation field level of the decommissioning A & R: 2.5cm A &R: 1cm site needs to be displayed in real time with the 1 3.177 3.178 decontamination and demolition of the nuclear 2 3.088 3.088 facilities. 3 2.836 2.837 4 3.090 3.091 5 3.502 3.503 6 3.657 3.658 7 3.945 3.946 8 4.128 4.129 Figure 1.1 Visual display of 3D radiation field 9 3.955 3.956 10 4.260 4.262 • It requires the radiation field calculation module Figure 3.1 Single source calculation model (0.662MeV, 1.0E6Bq) Table 3.5 Calculation accuracy of the coupling method for different to calculate the radiation field level quickly and energy group spacings in the single source model accurately. Table 3.1 Result comparison of the coupling method and MC FMC-PK/FMC method in the single source model Point 0.2MeV 0.05MeV • Space problem: geometrically simple regions γ energy flux rates (MeV/(cm2•s)) 1 3.177 2.968 and geometrically complex regions both exist in Ratios 2 3.088 2.877 Coupling in the decommissioning site. Point MCNP5 3 2.836 2.600 method 4 3.090 2.913 1 4.6630E-03 1.4678E-03 3.177 • The three basic methods of radiation field 5 3.502 3.260 2 3.5499E-03 1.1497E-03 3.088 6 3.657 3.443 calculation can not solve the problem and then 7 3.945 3.696 can not meet the requirements of calculation 3 8.3030E-04 2.9278E-04 2.836 8 4.128 3.855 speed and accuracy. 4 6.4362E-03 2.0831E-03 3.090 9 3.955 3.731 10 4.260 4.014 5 2.2223E-03 6.3462E-04 3.502 • The MC-PK coupling method is proposed. 6 4.5223E-03 1.2367E-03 3.657 • γ energy↑ accuracy↑ Table 1.1 The adv. and disad. of three basic methods 7 2.8517E-03 7.2289E-04 3.945 • Grid number↑ calculation time↑ Discrete Point kernel Monte Carlo 8 2.0303E-03 4.9189E-04 4.128 • MC-PK coupling method has good stability. Methods ordinate integration (MC) 9 3.8058E-03 9.6218E-04 3.955 (SN) (PK) Applicable 10 3.1389E-03 7.3675E-04 4.260 Simple Simple Complex Calculation of Qinshan Ⅰ region • FMC-PK/FMC = 3-5 Calculation Long Short Long • Calculation time(CPU: 2.9GHz Memory: 2GB) time • MCNP5: 4h Accuracy High Low High • MC-PK: 30min + 40min MC-PK Coupling Method Table 3.2 Calculation accuracy of the coupling method for • Basic idea: different photon energies in the single source model • To use MC in geometrically complex regions; FMC-PK/FMC • To use PK in geometrically simple regions; • To transform particle parameters on the Point 0.662MeV 0.83MeV 1.17MeV 1.33MeV coupling surface. 1 3.177 2.509 1.792 1.847 2 3.088 2.857 1.947 1.879 3 2.836 2.667 1.913 1.964 4 3.090 2.885 2.017 1.948 5 3.502 3.247 2.270 2.168 Figure 4.1 The simplified model of Qinshan I 6 3.657 3.229 2.232 2.137 Table 4.1 Result comparison of the coupling method and MC 7 3.945 3.464 2.326 2.279 method in Qinshan I 8 4.128 3.685 2.534 2.370 Air-absorbed dose rates (mSv/h) 9 3.955 3.556 2.441 2.298 Coupling Ratios Point MCNP5 10 4.260 3.726 2.524 2.396 method 1 4.3988E+02 7.3010E+01 6.02 2 4.1344E+03 2.7059E+03 1.53 Figure 2.1 Decommissioning site division of a nuclear power plant Table 3.3 Calculation accuracy of the coupling method for different shield thicknesses in the single source model 3 6.3208E+03 6.4018E+03 0.99 4 4.1344E+03 6.5739E+03 0.63 • Assuming that the MC region is a pure emitter, it FMC-PK/FMC can be derived from the Boltzmann equation: • D /D < 7 Point 5cm 20cm 40cm MC-PK MC + jn (rrss,,,, E) = SA ( E ) • Calculation time (CPU: 4.6GHz Memory: 64GB) 1 2.521 3.177 5.362 • MCNP5: 8h 2 2.712 3.088 2.765 • MC-PK: 1h + 5min 3 3.268 2.836 0.897 4 2.672 3.090 3.821 5 3.137 3.502 3.017 Conclusions 6 2.939 3.657 4.106 • The coupling method has reasonable accuracy, 7 3.176 3.945 3.932 8 3.351 4.128 4.023 short calculation time and good stability. 9 3.139 3.955 4.709 • The MC-PK coupling method is adequate for the Figure 2.2 A source on the coupling surface is converted into a requirements of radiation field calculation in sphere source (the length of the arrow in the figure is 10 3.245 4.260 4.912 proportional to the γ energy) decommissioning of nuclear power plants. Contact References [1] GUO Ya-ping, FU Meng-ting, ZHANG Ze-huan, SONG Ying-ming. Study on the Rapid Calculation Method of Radiation Field.[A]. Proceeding of the 2nd Nuclear Power Plant Radioactive Waste Management Symposium of China Nuclear Energy Industry Association[C]: 2018: 69-75. [2] B. Davision. Neutron Transport Theory. Oxford University Press, 174(1957). GUO Yu-fei [3] XI Fu-kun, FANG Bang-cheng. A Monte Carlo Method of the Point Kernel Integration for Gamma-Rays Calculation. Nuclear Power Engineering, 1987,8(2): 66-74. [4] SONG Yu-shou, LIU Hui-lan, LI Wei. Monte Carlo Simulation of Particle Detectors[M]. Chongqing: Chongqing University Press, 2016.05: 4-10. Nuclear Power Institute of China [5] OUYANG Yu. Design and Construction of Qinshan Nuclear Power Plant[J]. Nuclear Power Engineering, 1985, 6(6): 481-489. [6] ZHU Mei-ying, XU Tao-zhong. Modeling Analysis and Estimation Report of Radioactive Retention[R]. Chengdu: Nuclear Power Institute of China, 2018. Sichuan Engineering Laboratory for Nuclear Facilities [7] Frank Herbert Attix. Introduction to Radiological Physics and Radiation Dosimetry[M]. Beijing: Atomic Energy Press, 2013.05: 41. [8] YANG De-feng, XUE Na, CHENG He-ping, MAO Ya-wei. MCNP and DORT Couptation Method[J]. Science and Technology Review, 2013,31(18): 53-56. Decommissioning and Radwaste Management [9] HAN Jing-ru, CHEN Yi-xue, SHI Sheng-chun, YUAN Long-jun, LU Dao-gang. Development of Three-dimensional Coupled Code System Based on Discrete Ordinates and Monte Carlo Method[J]. Chinese Journal of Nuclear Science Email:[email protected] and Engineering, 2012,32(2): 160-164. [10] BI Yuan-jie, YANG Hong-wei, GUO Jin-sen, LUO Zhi-ping, LI Chuan-long, CHEN Ling. Study on Fast Radiation Transport Calculation Methodology in Decommissioning of Research Reactor[J]. Atomic Energy Science and Phone:+86 18380364221 Technology, 2015, 49(S1): 470-474. [11] BI Yuan-jie, YANG Hong-wei, LUO Zhi-ping, GUO Jin-sen, LI Chuan-long. Coupling Method for Fast Calculation of Three-dimensional Radiation Field[J]. Annual Report of China Institute of Atomic Energy. 2013: 246. [12] ZHANG Yong-ling, HU Yi-fei, LIU Meng, ZHANG Kai-yun, DAI Bo, ZHANG Hang-zhou. Study of Nuclear Reactor Decommissioning Simulation Key Technologies[A]. Atomic Energy Press[C]. 2015:64-65-66-67-68-69-70. [13] LUO Shang-geng, ZHANG Zhen-tao, ZHANG Hua. Decommissioning of Nuclear Facilities and Radiation Facilities[M]. Beijing: China Environmental Science Press, 2011. [14] WANG Jun-feng, WANG Shao, Liu Kun-xian, ZHANG Tian-xiang. Decommissioning of Nuclear Facilities[M]. Beijing: Atomic Energy Press, 2013. .