On the use of shannon entropy of the fission distribution for assessing convergence of Monte Carlo criticality calculations
- Los Alamos National Laboratory, MS F663, PO Box 1663, Los Alamos, NM 87544 (United States)
Monte Carlo calculations of k-eigenvalue problems are based on a power iteration procedure. To obtain correct results free of contamination from the initial guess for the fission distribution, it is imperative to determine when the iteration procedure has converged, so that a sufficient number of the initial batches are discarded prior to beginning the Monte Carlo tallies. In this paper, we examine the convergence behavior using both theory and numerical testing, demonstrating that k{sub eff} may converge before the fission distribution for problems with a high dominance ratio. Thus, it is necessary to assess convergence of both k{sub eff} and the fission distribution to obtain correct results. To this end, the Shannon entropy of the fission distribution has been found to be a highly effective means of characterizing convergence of the fission distribution. The latest version of MCNP5 includes new capabilities for computing and plotting the Shannon entropy of the fission distribution as an important new tool for assessing problem convergence. Examples of the application of this new tool are presented for a variety of practical criticality problems. (authors)
- Research Organization:
- American Nuclear Society, 555 North Kensington Avenue, La Grange Park, IL 60526 (United States)
- OSTI ID:
- 22039478
- Resource Relation:
- Conference: PHYSOR-2006: American Nuclear Society's Topical Meeting on Reactor Physics - Advances in Nuclear Analysis and Simulation, Vancouver, BC (Canada), 10-14 Sep 2006; Other Information: Country of input: France; 8 refs.
- Country of Publication:
- United States
- Language:
- English
Similar Records
Automated convergence detection in Monte Carlo criticality calculations using student's bridge statistics based on K{sub eff} or Shannon entropy
"K-Effective of the World" and Other Concerns for Monte Carlo Eigenvalue Calculations