-
J Monte Carlo transport code (JMCT), a new three-dimentional (3 D) Monte Carlo transport code, is introduced in this paper. The code is developed on the basis of 3 D geometry infrastructure JCOGIN and composed of multilayer modules. JMCT is capable of simulating the collision of particles with multi-group energy or providing energy data libraries. Two forms of parallelism supported in JMCT are domain decomposition and domain replication. The code has very good expansibility. JMCT Monte Carlo results have been compared with hot zero power (HZP) measurements of BEAVRS benchmark model from the MIT Computational Reactor Physics Group. Included in the comparisons are the eigenvalues, control rod bank worths, isothermal temperature coefficients, axially integrated full core detector measurements, axial detector profiles, etc. The eigenvalues for the HZP condition with different control rods positions and boron concentrations are calculated and the error is less than 0.2% compared with the theoretical error 1.000. The results of JMCT for isothermal temperature coefficients are also listed together with MC21 results and measured data. Each calculation for the eigenvalue is run by 1000 cycles in total, discarding 600 cycles, tracking 4 million neutrons each cycle. It takes 5.3-5.7 hours to run on 200 CPU cores. The JMCT results of axially integrated radial detector relative power distribution (RPD) and axial normalized detector signal are compared with the measured data. Power depressions from grid spacers are clearly seen in the JMCT results and accord with the measured data. The JMCT results of axially integrated assembly RPD power distribution are in good agreement with MC21 results, the maximum difference being 3.173% for 193 assemblies. So is the result of pin power RPD relative power at the axial elevation of peak power; the minimum relative power RPD 0.278 of JMCT is comparable to 0.283 of MC21, and the max relative power RPD 2.422 of JMCT is comparable to 2.452 of MC21. The calculation for RPD is run with 3000 inactive cycles and 5000 active cycles, tracking 4 million particles each cycle. It takes about 4.8 days to run on 200 CPU cores. Shannon entropy is used to demonstrate that the fission source distribution is converged after 3000 inactive cycles. With the development of computers and parallel computing, the Monte Carlo method can be used in reactor design instead of benchmarking other calculated results.
-
Keywords:
- pin-by-pin /
- BEAVRS model /
- Monte Carlo /
- J Monte Carlo transport code
[1] Forrest B, Brian K, David B 2012 Proceedings of PHYSOR 2012-Advances in ReactorPhysics-Linking Research, Industry, and Education, Knoxville, Tennessee, USA, April 15-20, 2012
[2] Du S H, Zhang S F, Feng T G,Wang Y Z, Xing J R 1989 Computer Simulation of Transport Problems (Changsha: Hunan Science and Technology Press) p304 (in Chinese) [杜书华, 张树发, 冯庭桂, 王元璋, 邢静茹 1989 输运问题的计算机模拟(长沙: 湖南科技出版社) 第304页]
[3] Pei L C, Zhang X Z 1980 Monte Carlo Methods and Application in Particle Transportation (Beijing: Science Press) p1 (in Chinese) [裴鹿成,张孝泽 1980 蒙特卡罗方法及其在粒子输运问题中的应用(北京:科学出版社) 第1页]
[4] Deng L, Li G 2010 Chin. J. Comput. Phys. 27 791 (in Chinese) [邓力, 李刚 2010 计算物理 27 791]
[5] Bill M 2007 M&C+SNA, Monterey California, USA, April 15-19, 2007
[6] Forrest B, William M PHYSOR 2010, Pittsburgh, Pennsylvania, USA, May 9-14 2010
[7] Kord S 2003 M&C+SNA, Gatlinburg, Tennessee, USA, April 6-11, 2003
[8] Kord S, Benoit F, 2013 M&C 2013 Sun Valley, Idaho, USA, May 5-9, 2013 p1809
[9] Daniel K, Thomas S, 2012 Proceedings of PHYSOR 2012 Advances in Reactor Physics-Linking Research, Industry, and Education, Knoxville, Tennessee, USA, April 15-20, 2012
[10] Daniel K, Brian A, 2013 M&C 2013 Sun Valley, Idaho, USA, May 5-9, 2013 p2962
[11] Daniel K, Brian A, Paul R, 2014 PHYSOR 2014, Kyoto, Japan, September 28-October 3, 2014
[12] Li G, Zhang B Y, Deng L, 2013 High Power Laser and Particle Beams 25 158 (in Chinese) [李刚, 张宝印, 邓力 2013 强激光与粒子数 25 158]
[13] Li G, Zhang B Y, Deng L 2013 ANS Transactions 109 1425
[14] Zhang B Y, Li G, Deng L 2013 High Power Laser and Particle Beams 25 173 (in Chinese) [张宝印, 李刚, 邓力 2013 强激光与粒子数 25 173]
[15] Taro U, Forrest B 2005 Nuclear Science and Engineering 149 38
-
[1] Forrest B, Brian K, David B 2012 Proceedings of PHYSOR 2012-Advances in ReactorPhysics-Linking Research, Industry, and Education, Knoxville, Tennessee, USA, April 15-20, 2012
[2] Du S H, Zhang S F, Feng T G,Wang Y Z, Xing J R 1989 Computer Simulation of Transport Problems (Changsha: Hunan Science and Technology Press) p304 (in Chinese) [杜书华, 张树发, 冯庭桂, 王元璋, 邢静茹 1989 输运问题的计算机模拟(长沙: 湖南科技出版社) 第304页]
[3] Pei L C, Zhang X Z 1980 Monte Carlo Methods and Application in Particle Transportation (Beijing: Science Press) p1 (in Chinese) [裴鹿成,张孝泽 1980 蒙特卡罗方法及其在粒子输运问题中的应用(北京:科学出版社) 第1页]
[4] Deng L, Li G 2010 Chin. J. Comput. Phys. 27 791 (in Chinese) [邓力, 李刚 2010 计算物理 27 791]
[5] Bill M 2007 M&C+SNA, Monterey California, USA, April 15-19, 2007
[6] Forrest B, William M PHYSOR 2010, Pittsburgh, Pennsylvania, USA, May 9-14 2010
[7] Kord S 2003 M&C+SNA, Gatlinburg, Tennessee, USA, April 6-11, 2003
[8] Kord S, Benoit F, 2013 M&C 2013 Sun Valley, Idaho, USA, May 5-9, 2013 p1809
[9] Daniel K, Thomas S, 2012 Proceedings of PHYSOR 2012 Advances in Reactor Physics-Linking Research, Industry, and Education, Knoxville, Tennessee, USA, April 15-20, 2012
[10] Daniel K, Brian A, 2013 M&C 2013 Sun Valley, Idaho, USA, May 5-9, 2013 p2962
[11] Daniel K, Brian A, Paul R, 2014 PHYSOR 2014, Kyoto, Japan, September 28-October 3, 2014
[12] Li G, Zhang B Y, Deng L, 2013 High Power Laser and Particle Beams 25 158 (in Chinese) [李刚, 张宝印, 邓力 2013 强激光与粒子数 25 158]
[13] Li G, Zhang B Y, Deng L 2013 ANS Transactions 109 1425
[14] Zhang B Y, Li G, Deng L 2013 High Power Laser and Particle Beams 25 173 (in Chinese) [张宝印, 李刚, 邓力 2013 强激光与粒子数 25 173]
[15] Taro U, Forrest B 2005 Nuclear Science and Engineering 149 38
Catalog
Metrics
- Abstract views: 5866
- PDF Downloads: 748
- Cited By: 0
下载:
