Efficient Monte Carlo algorithm of time-dependent particle transport problem in multi-physics coupling calculation

被引:3
|
作者
Shangguan Dan-Hua [1 ]
Yan Wei-Hua [1 ]
Wei Jun-Xia [1 ]
Gao Zhi-Ming [1 ]
Chen Yi-Bing [1 ]
Ji Zhi-Cheng [1 ]
机构
[1] Inst Appl Phys & Computat Math, Beijing 100094, Peoples R China
基金
中国国家自然科学基金;
关键词
multi-physics calculation; time-dependent transport problem; Monte Carlo method; DIAGNOSTICS;
D O I
10.7498/aps.71.20211474
中图分类号
O4 [物理学];
学科分类号
0702 ;
摘要
Multi-physics coupling calculation has applications in many important research fields. If particle transportprocess is included in this calculation, Monte Carlo method is often used to simulate this process and usually alarge amount of calculation time is needed. So, efficient Monte Carlo algorithm for time-dependent particletransport problem is important for an efficiently coupling calculation, which inevitably relies on large-scaleparallel calculation. Based on the characteristic of time-dependent particle transport problem, two methods areproposed in this paper to achieve high- efficiency calculation. One is a tally-reducing algorithm which is used inthe coupling of transport simulation and burnup calculation. By reducing the quantity of data which should bereduced necessarily, this method can reduce the calculation time largely. It can be seen that a new couplingmode for these two processes in MPI environment has a larger value when model scale is larger than the samplesize. The other method is an adaptive method of setting the sample size of Monte Carlo simulation. The law oflarge number assures that the Monte Carlo method will obtain an exact solution when the sample scale tends toinfinity. But generally, no one knows which sample scale is big enough for obtaining a solution with targetprecision in advance. So, the common strategy is to set a huge-enough sample scale by experience and conductthe posterior check for all results. Apparently, this way cannot be efficient because the calculation will go onafter the precision of solution has reached an object value. Another popular method is to set the sample size torely on the relative error of some single calculation. The sample size is enlarged without a break until therelative error is less than some presetting value. This method is not suitable either, because Monte Carloparticle transport simulation will gives feedbacks to other process which is composed of many tallies. It isinappropriate to adjust the sample size according to the relative error of any calculation. Relying on thegeneralization of the Shannon entropy concept and an on-the-fly diagnosis rule for a entropy value sequence, theadaptive method proposed in this paper can reduce the original huge sample scale to a reasonable level. Bynumerically testing some non-trivial examples, both algorithms can reduce the calculation time largely, with theresults kept almost unchanged, so the efficiency is high in these cases
引用
收藏
页数:8
相关论文
共 15 条
  • [1] Domain decomposition models for parallel Monte Carlo transport
    Alme, HJ
    Rodrigue, GH
    Zimmerman, GB
    [J]. JOURNAL OF SUPERCOMPUTING, 2001, 18 (01): : 5 - 23
  • [2] Brown F B, 2004, LAUR044532
  • [3] Danhua SG, 2016, NUCL SCI ENG, V182, P555
  • [4] Du S H, 1989, COMPUTER SIMULATION, pp47
  • [5] Kelly D J, 2013, P M C 2013 SUN VALL, P2962
  • [6] Kelly D J, 2012, P PHYSOR 2012 KNOXV, P1
  • [7] Li G, 2014, NUCL POWER ENG, VS2, P228
  • [8] The sandwich method for determining source convergence in Monte Carlo calculation
    Naito, Y
    Yang, JA
    [J]. JOURNAL OF NUCLEAR SCIENCE AND TECHNOLOGY, 2004, 41 (05) : 559 - 568
  • [9] Romano P.K., 2011, Progress in Nuclear Science and Technology, V2, P670
  • [10] Shangguan D H, 2018, HIGH POWER LASER PAR, V30