-
在多计算步的非定常输运问题的蒙特卡罗模拟中, 为自动调整每一步的样本数以获得较高的计算效率, 可以有多种准则. 一种可选的方法是在每一步每隔若干样本监测一次系统中未死亡粒子属性分布对应的香农熵的收敛情况以决定何时停止追加样本, 此种方法需要在每一步频繁计算香农熵值. 由于在MPI消息传递并行编程环境下香农熵的经典计算方法必须广播大量的数据, 导致每一步的计算时间随香农熵计算频率的提高而快速增大, 这显然是不能满足实际需求的. 本文提出了一种适应于消息传递并行编程环境的香农熵计算新方法, 该方法计算得到的香农熵值并不等价于经典方法, 但二者之间的差别会随着样本数的增加而趋于零. 新方法的最大优势是高频计算香农熵值的时间代价大为降低, 为最终实现基于香农熵收敛判断的每步样本数的自动调整奠定了必要的基础.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
计量
- 文章访问数: 6553
- PDF下载量: 168
- 被引次数: 0