-
For the Monte Carlo simulation of the non-static transport problem, there must be many calculation steps. Because some particles cannot finish their transport in the last step, they are naturally used as the source particles of the present step. These particles are called undied particles. It is difficult to adjust the history number of each step to obtain higher efficiency because the adjusting rule is hard to find. The most direct method is to set a large enough history number for all steps. But evidently, it is unnecessary for some steps. Among all possible rules, one candidate of adjusting the history number is to check the convergence situation of Shannon entropy (corresponding to the distribution of some undied particle attributes) every some samples in each step to determine whether or not to simulate more particles. So, this method needs to calculate the Shannon entropy frequently. Because the classical method of calculating Shannon entropy in message passing parallel programming environment must reduce massive data, it is unpractical to be used in this situation for the great increasing of computation time with the high frequency of entropy calculation. In this paper, we propose an efficient method of calculating the entropy in the message passing parallel programming environment by letting each process calculate its entropy value based on the local data in each processer and calculating the final entropy by averaging all the entropy values gotten by all processes. The entropy value calculated by this method is not the same as that by the classical method when using finite history number, but the difference goes to zero when the history number goes to infinity. The most remarkable advantage of this method is the small increasing of computation time when calculating the entropy frequently. It is a suitable method of calculating Shannon entropy when adjusting the history number automatically based on the judgment of the convergence situation of Shannon entropy.
-
Keywords:
- non-static transport problem /
- Monte Carlo method /
- Shannon entropy
[1] Martin W 2012 Nucl. Eng. Technol. 44 2
[2] Li G, Deng L, Li S, Mo Z Y 2011 Acta Phys. Sin. 60 022401 (in Chinese) [李刚, 邓力, 李树, 莫则尧 2011 物理学报 60 022401]
[3] Li S, Li G, Tian D F, Deng L 2013 Acta Phys. Sin. 62 249501 (in Chinese) [李树, 李刚, 田东风, 邓力 2013 物理学报 62 249501]
[4] Horelik N, Herman B R, Forget B, Smith K 2013 Proceedings of MC 2013 Sun Valley, Idaho, USA, May 5-9, 2013 p2986
[5] Smith K S, Forget B 2013 Proceedings of MC 2013, Sun Valley, Idaho, USA, May 5-9, 2013 p1809
[6] Lan K, He X T, Lai D X, Li S G 2006 Acta Phys. Sin. 55 3789 (in Chinese) [蓝可, 贺贤土, 赖东显, 李双贵 2006 物理学报 55 3789]
[7] Brown F B 2006 Proceedings of PHYSOR 2006 Vancouver, BC, Canada, September 10-14, 2006C081
[8] Ueki T, Brown F B 2005 Nucl. Sci. Eng. 149 38
-
[1] Martin W 2012 Nucl. Eng. Technol. 44 2
[2] Li G, Deng L, Li S, Mo Z Y 2011 Acta Phys. Sin. 60 022401 (in Chinese) [李刚, 邓力, 李树, 莫则尧 2011 物理学报 60 022401]
[3] Li S, Li G, Tian D F, Deng L 2013 Acta Phys. Sin. 62 249501 (in Chinese) [李树, 李刚, 田东风, 邓力 2013 物理学报 62 249501]
[4] Horelik N, Herman B R, Forget B, Smith K 2013 Proceedings of MC 2013 Sun Valley, Idaho, USA, May 5-9, 2013 p2986
[5] Smith K S, Forget B 2013 Proceedings of MC 2013, Sun Valley, Idaho, USA, May 5-9, 2013 p1809
[6] Lan K, He X T, Lai D X, Li S G 2006 Acta Phys. Sin. 55 3789 (in Chinese) [蓝可, 贺贤土, 赖东显, 李双贵 2006 物理学报 55 3789]
[7] Brown F B 2006 Proceedings of PHYSOR 2006 Vancouver, BC, Canada, September 10-14, 2006C081
[8] Ueki T, Brown F B 2005 Nucl. Sci. Eng. 149 38
Catalog
Metrics
- Abstract views: 6532
- PDF Downloads: 167
- Cited By: 0
下载:
