Search NASA⌕ Search

Engineering topics

Mosher, Scott W.

Publications and source records attributed to Mosher, Scott W..

As-Built Simulation of the High Flux Isotope Reactor

The Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) is an 85 MWt flux trap-type research reactor that supports key research missions, including isotope production, materials irradiation, and neutron scattering. The core consists of an inner and an outer fuel element containing 171 and 369 involute-shaped plates, respectively. The thin fuel plates consist of a U 3 O 8 -Al dispersion fuel (highly enriched), an aluminum-based filler, and aluminum cladding. The fuel meat thickness is varied across the width of the involute plate to reduce thermal flux peaks at the radial edges of the fuel elements. Some deviation from the designed fuel meat shaping is allowed during manufacturing. A homogeneity scan of each fuel plate checks for potential anomalies in the fuel distribution by scanning the surface of the plate and comparing the attenuation of the beam to calibration standards. While typical HFIR simulations use homogenized fuel regions, explicit models of the plates were developed under the Low-Enriched Uranium Conversion Program. These explicit models typically include one inner and one outer fuel plate with nominal fuel distributions, and then the plates are duplicated to fill the space of the corresponding fuel element. Therefore, data extracted from these simulations are limited to azimuthally averaged quantities. To determine the reactivity and physics impacts of an as-built outer fuel element and generate azimuthally dependent data in the element, 369 unique fuel plate models were generated and positioned. This model generates the three-dimensional (i.e., radial–axial–azimuthal) plate power profile, where the azimuthal profile is impacted by features within the adjacent control element region and beryllium reflector. For an as-built model of the outer fuel element, plate-specific homogeneity data, 235 U loading, enrichment, and channel thickness measurements were translated into the model, yielding a much more varied azimuthal power profile encompassed by uncertainty factors in analyses. These models were run with the ORNL-TN and Shift Monte Carlo tools, and they contained upwards of 500,000 cells and 100,000 unique tallies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Application of the Denovo Discrete Ordinates Radiation Transport Code to Large-Scale Fusion Neutronics

Fusion energy systems pose unique challenges to the modeling and simulation community. These challenges must be met to ensure the success of the ITER experimental fusion reactor. ITER’s complex systems require detailed modeling that goes beyond the scale of comparable simulations to date. In this work, the Denovo radiation transport code was used to calculate neutron fluence and kerma for the JET streaming benchmark. This work was performed on the Titan supercomputer at the Oak Ridge Leadership Computing Facility. Denovo is a novel three-dimensional discrete ordinates transport code designed to be highly scalable. Sensitivity studies have been completed to examine the impact of several deterministic parameters. Furthermore, results were compared against experiment as well as the MCNP and Shift Monte Carlo codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Algorithmic Improvements to MCNP5 for High-Resolution Fusion Neutronics Analyses

Neutronics analyses of the ITER experimental fusion reactor rely on increasingly complex geometry models and estimates of energy-dependent neutron flux and radiation dose-rate distributions generated at ever higher resolutions. There are significant practical challenges with applying the Monte Carlo N-Particle (MCNP) continuous-energy transport code to high-resolution analyses. For models consisting of more than 100 000 surfaces and cells, geometry initialization can take several hours, thus slowing down model integration and transport analysis efforts. In multithreaded simulations, the amount of memory consumed by superimposed mesh tally data increases in proportion to the number of threads. This behavior limits either the tally resolution or the number of processor cores that can be utilized in the simulation. This paper describes algorithmic improvements that were implemented in a modified version of MCNP5 to overcome these limitations. These improvements are referred to as the Oak Ridge National Laboratory Transformative Neutronics (ORNL-TN) upgrade. Here, a comparison of the performance and memory usage of both MCNP5 and ORNL-TN on several relevant fusion neutronics models is presented. In these tests and in actual high-resolution neutronics analyses, ORNL-TN reduces geometry processing times from hours to a few seconds and increases in-memory mesh tally capacity from the order of 10 8 to 10 10 space-energy bins.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integration of the Full Tokamak Reference Model with the Complex Model for ITER Neutronic Analysis

The ITER International Organization has developed a number of reference Monte Carlo N-Particle (MCNP) models including the tokamak machine C-model, the Tokamak Complex model, and the neutral beam injection (NBI) systems model. The Tokamak Complex model primarily describes building structures beyond the bioshield. Representation of the tokamak and its systems are not included in this model. The Oak Ridge National Laboratory Radiation Transport Group has conducted two ITER neutronic analysis model integrations: (1) integration of the tokamak C-model with the Tokamak Complex model for shutdown dose rate characterization in Port Cell 16 at level B1, and (2) integration of the NBI model with the Tokamak Complex model for estimating the spatial distribution of biological dose rate at levels L1, L2, and L3 of the Tokamak Complex. The integrated models were further extended to include models of system components that are essential to the neutronic analyses. Furthermore, this paper presents the approach and computer tools used to integrate existing reference models, describes the additional design details implemented in the integrated models, and provides representative neutronic calculations based on the extended models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Estimating consistent biased source probabilities over arbitrary bins - 25621

This paper describes computational algorithms for estimating consistent biased source probabilities for the CADIS and FW-CADIS methods that are compatible with most of the general fixed-source capabilities of the MCNP continuous-energy transport code. The basic algorithm is a novel Monte Carlo sampling technique that enables the straightforward estimation of biased source probabilities over arbitrarily complex distribution trees. In recent work, a method has been developed for generating correct biased spatial source probabilities in the presence of rejection from an arbitrary number of geometry cells. Two numerical test problems that demonstrate the correctness of this new algorithm are presented, including a simplified model of a cooling water system with an N-17 neutron source. (authors)

07 ISOTOPE AND RADIATION SOURCES↗