DOI QR코드

DOI QR Code

Analysis of inconsistent source sampling in monte carlo weight-window variance reduction methods

  • Griesheimer, David P. (Naval Nuclear Laboratory) ;
  • Sandhu, Virinder S. (Naval Nuclear Laboratory)
  • Received : 2017.05.31
  • Accepted : 2017.07.27
  • Published : 2017.09.25

Abstract

The application of Monte Carlo (MC) to large-scale fixed-source problems has recently become possible with new hybrid methods that automate generation of parameters for variance reduction techniques. Two common variance reduction techniques, weight windows and source biasing, have been automated and popularized by the consistent adjoint-driven importance sampling (CADIS) method. This method uses the adjoint solution from an inexpensive deterministic calculation to define a consistent set of weight windows and source particles for a subsequent MC calculation. One of the motivations for source consistency is to avoid the splitting or rouletting of particles at birth, which requires computational resources. However, it is not always possible or desirable to implement such consistency, which results in inconsistent source biasing. This paper develops an original framework that mathematically expresses the coupling of the weight window and source biasing techniques, allowing the authors to explore the impact of inconsistent source sampling on the variance of MC results. A numerical experiment supports this new framework and suggests that certain classes of problems may be relatively insensitive to inconsistent source sampling schemes with moderate levels of splitting and rouletting.

Keywords

References

  1. G.I. Bell, S. Glasstone, Nuclear Reactor Theory, Van Nostrand and Reinhold, New York, 1967.
  2. J.C. Wagner, A. Haghighat, Automated variance reduction of Monte Carlo shielding calculations using the discrete ordinates adjoint function, Nucl. Sci. Eng. 128 (1998) 186-208. https://doi.org/10.13182/NSE98-2
  3. J.C. Wagner, D.E. Peplow, S.W. Mosher, FW-CADIS method for global and regional variance reduction of Monte Carlo radiation transport calculations, Nucl. Sci. Eng. 176 (2014) 37-57. https://doi.org/10.13182/NSE12-33
  4. S.C. Wilson, D.P. Griesheimer, Rejection sampling from a spherical harmonics angular flux functional expansion for a discrete ordinates to Monte Carlo splice, Proc. International Conference on Mathematics, Computational Methods & Reactor Physics (M&C 2009), CD-ROM, Saratoga Springs, NY, 2009.
  5. X-5 Monte Carlo Team, MCNP - A General Monte Carlo N-Particle Transport Code, Version 5, LA-CP-03-0284, Los Alamos National Laboratory, 2003.
  6. T.E. Booth, J.S. Hendricks, Importance estimation in forward Monte Carlo calculations, Nucl. Technol. Fusion 5 (1984) 90-100. https://doi.org/10.13182/FST84-A23082
  7. J. Spanier, E.M. Gelbard, Monte Carlo Principles and Neutron Transport Problems, Dover Publications, New York, 2008.