Date of Award


Degree Type


Degree Name

Doctor of Philosophy


Nuclear Engineering

Major Professor

G. Ivan Maldonado

Committee Members

Ronald E. Pevey, Arthur E. Ruggles, James D. Freels, David H. Cook, David Chandler (courtesy member)


The computational ability to accurately predict the conditions in an experiment under irradiation is a valuable tool in the operation of a research reactor whose scientific mission includes isotope production, materials irradiation, and neutron activation analysis. Understanding of different governing physics is required to ascertain satisfactory conditions within the experiment: the neutron transport behavior throughout the reactor and the coupled behavior of heat transfer, structural mechanics and fluid flow. Computational methods and tools were developed for robust numerical analysis of experiment behavior at the Oak Ridge National Laboratory (ORNL) High Flux Isotope Reactor (HFIR), including fully-coupled thermo-mechanics in three plutonium-238 (238Pu) production targets. In addition, a new computational tool was developed that solves neutron transport using the discrete ordinates method on a finite element mesh and offers multiphysics coupling.

The thermo-mechanical models of the 238Pu targets are solved using the COMSOL heat transfer and solid mechanics modules with irradiation behavior and thermophysical properties taken from measurements performed at ORNL. The experiments, placed in the permanent beryllium (PB) reflector, consist of neptunium dioxide/aluminum (NpO2/Al) pellets in Al containment, the model taking advantage of axisymmetry in two-dimensional R-Z cylindrical geometry. At times, extended analysis was needed for incomplete data sets and time schedules; however, the thermal-structure models ensured progression through three project phases of target qualification.

The neutron transport equation was solved in COMSOL, using the discrete ordinates formulation in the weak form partial differential equation (PDE) interface. Validation studies were performed for the dimensions developed (one-, two- and three- dimensional Cartesian as well as two-dimensional R-Z cylindrical/axi-symmetric) and compared to external deterministic and stochastic codes. The method was then applied to a beginning-of-cycle (BOC) simplified HFIR core, with good comparison to other static solutions of the HFIR, and a time-dependent extension to this tool was created and exhibited for a benchmark problem.

The research presented in this dissertation is the continued progress towards creating a comprehensive multiphysics methodology for studying the dynamic behavior of the HFIR core, and shows the capabilities of detailed space-time reactor physics studies and of multiphysics analyses for experiment qualification and safety analyses at a research reactor.

Available for download on Friday, December 15, 2017

Files over 3MB may be slow to open. For best results, right-click and select "save as..."