Development of high fidelity methods for 3D Monte Carlo transient analysis of nuclear reactors

Thesis: Ph. D., Massachusetts Institute of Technology, Department of Nuclear Science and Engineering, 2018.

Bibliographic Details
Main Author: Shaner, Samuel Christopher
Other Authors: Kord S. Smith and Benoit Forget.
Format: Thesis
Language:eng
Published: Massachusetts Institute of Technology 2018
Subjects:
Online Access:http://hdl.handle.net/1721.1/119034
_version_ 1826188080122429440
author Shaner, Samuel Christopher
author2 Kord S. Smith and Benoit Forget.
author_facet Kord S. Smith and Benoit Forget.
Shaner, Samuel Christopher
author_sort Shaner, Samuel Christopher
collection MIT
description Thesis: Ph. D., Massachusetts Institute of Technology, Department of Nuclear Science and Engineering, 2018.
first_indexed 2024-09-23T07:54:17Z
format Thesis
id mit-1721.1/119034
institution Massachusetts Institute of Technology
language eng
last_indexed 2024-09-23T07:54:17Z
publishDate 2018
publisher Massachusetts Institute of Technology
record_format dspace
spelling mit-1721.1/1190342019-04-09T15:27:08Z Development of high fidelity methods for 3D Monte Carlo transient analysis of nuclear reactors Shaner, Samuel Christopher Kord S. Smith and Benoit Forget. Massachusetts Institute of Technology. Department of Nuclear Science and Engineering. Massachusetts Institute of Technology. Department of Nuclear Science and Engineering. Nuclear Science and Engineering. Thesis: Ph. D., Massachusetts Institute of Technology, Department of Nuclear Science and Engineering, 2018. This electronic version was submitted by the student author. The certified thesis is available in the Institute Archives and Special Collections. Cataloged from student-submitted PDF version of thesis. Includes bibliographical references (pages 136-140). Monte Carlo is increasingly being used to perform high-fidelity, steady-state neutronics analysis of power reactor geometries on today's leadership class supercomputers. Extending Monte Carlo to time dependent problems has proven to be a formidable challenge due to the significant computational resource and data processing requirements. In this thesis, a transient methodology is proposed and implemented to enable accurate and computationally tractable time dependent Monte Carlo analysis. The frequency transform method has been described and implemented in Monte Carlo for the first time. The attractiveness of this method lies in its ability to accurately capture the space and time dependent distribution of the delayed neutron source throughout a transient. Nuances to the algorithmic implementation are described and validated through a series of simple analytical test problems. Comparison with the adiabatic method currently employed for Monte Carlo transient analysis shows significant improvement in the spatial distribution and magnitude of the power for a negative reactivity insertion transient in the 2D and 3D C5G7 geometry. To aid in understanding the effect of statistical uncertainty in the tallied quantities on the time dependent flux solution, a simplified point kinetics model was developed and used for insightful analysis on simple transient test problems. This revealed how the time dependent flux profiles for a series of independent trials can be approximated by a normal distribution at low uncertainties in the tallied reactivity, but deviates from a normal distribution at relatively modest uncertainties in reactivity. Given the compuational constraints of solving large problems, having a simple model that can provide insight on the expected behavior and flux distribution can be very valuable. The frequency transform methodology belongs to a class of indirect space-time factorization methods that perform high-order calculations (e.g. Monte Carlo) over long time steps and low-order, computationally-efficient calculations (e.g. Point Kinetics) over short time steps as an approach to balance performance and accuracy. The coarse mesh finite difference (CMFD) diffusion operator is employed as the low-order solver in Monte Carlo transient analysis for the first time. The CMFD diffusion operator is attractive due to its potential to increase the time step size between the computationally expensive high-order solves. Implementing this methodology is important as continuous energy Monte Carlo is reactor-agnostic and able to treat complex geometries without difficulty, opening up the possibility of solving transients on new experimental geometries for which there is little data. by Samuel Christopher Shaner. Ph. D. 2018-11-15T15:51:54Z 2018-11-15T15:51:54Z 2018 2018 Thesis http://hdl.handle.net/1721.1/119034 1059513405 eng MIT theses are protected by copyright. They may be viewed, downloaded, or printed from this source but further reproduction or distribution in any format is prohibited without written permission. http://dspace.mit.edu/handle/1721.1/7582 140 pages application/pdf Massachusetts Institute of Technology
spellingShingle Nuclear Science and Engineering.
Shaner, Samuel Christopher
Development of high fidelity methods for 3D Monte Carlo transient analysis of nuclear reactors
title Development of high fidelity methods for 3D Monte Carlo transient analysis of nuclear reactors
title_full Development of high fidelity methods for 3D Monte Carlo transient analysis of nuclear reactors
title_fullStr Development of high fidelity methods for 3D Monte Carlo transient analysis of nuclear reactors
title_full_unstemmed Development of high fidelity methods for 3D Monte Carlo transient analysis of nuclear reactors
title_short Development of high fidelity methods for 3D Monte Carlo transient analysis of nuclear reactors
title_sort development of high fidelity methods for 3d monte carlo transient analysis of nuclear reactors
topic Nuclear Science and Engineering.
url http://hdl.handle.net/1721.1/119034
work_keys_str_mv AT shanersamuelchristopher developmentofhighfidelitymethodsfor3dmontecarlotransientanalysisofnuclearreactors