Investigation of improved methods for assessing convergence of models in MCNP using Shannon entropy

Monte Carlo computationals methods are widely used in academia to analyze nuclear systems design and operation because of their high accuracy and the relative ease of use in comparison to deterministic methods. However, current Monte Carlo codes require an extensive knowledge of the physics of a problem as well as the computational methods being used in order to ensure accuracy. This investigation aims to provide better on-the-fly diagnostics for convergence using Shannon entropy and statistical checks for tally undersampling in order to reduce the burden on the code user, hopfully increasing the use and accuracy of Monte Carlo codes. These methods were tested by simulating the OECD/NEA benchmark #1 problem in MCNP. It was found that Shannon entropy does accurately predict the number of batches required for a source distribution to converge, though only when when the Shannon entropy mesh was the size of the tally mesh. The investigation of undersampling showed evidence of methods to predict undersampling on-the-fly using Shannon entropy as well as laying out where future work should lead.

Department dc:contributor.department

Massachusetts Institute of Technology. Dept. of Nuclear Science and Engineering.

Grantor dc:publisher

Massachusetts Institute of Technology

Year dc:date.issued

2012

Author dc:creator

• Macdonald, Ruaridh (Ruaridh R.)

Advisor dc:contributor.advisor

• Benoit Forget.

• Nuclear Science and Engineering.

Statement dc:rights

• M.I.T. theses are protected by copyright. They may be viewed from this source for any purpose, but reproduction or distribution in any format is prohibited without written permission. See provided URL for inquiries about permission.

Licence dc:rights.uri

• http://dspace.mit.edu/handle/1721.1/7582

Language dc:language.iso

eng

Handle dc:identifier.uri

http://hdl.handle.net/1721.1/76954

OAI identifier oai:identifier

oai:dspace.mit.edu:1721.1/76954