Search NASASearch

SEARCH · Search NASA

Results for “MCNP Model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Characterizing the Experiment for Calibration with Uranium (Excalibur) neutron source for use in warhead verification

Neutron sources can play a variety of roles in warhead verification. For transmission radiography, a source of directed high energy neutrons is required, while for applications to detect fissile isotopes, sub-MeV neutrons are preferred. The Excalibur (Experiment for Calibration with Uranium) neutron source has been built and used in a variety of verification-related experiments. Excalibur is based on a commercial deuterium-tritium neutron generator specified and measured to be capable of producing 14 MeV neutrons at rates of up to 8.2 × 10 8 neutrons/s. Here, the generator is enclosed in a carbon-steel 32" diameter, 23.62" high carbon-steel cylinder that moderates the mean neutron energy to under 500 keV. This, in turn, is encased in 5%-borated polyethylene such that the entire assembly is a 48" x 48" box that is 30" tall. For radiographic applications, a narrow, tapered channel in the steel and polyethylene allows 14 MeV neutrons to stream directly from the generator to a test object. Its collimating capability is demonstrated by measuring the neutron flux profile. In the moderated mode of operation, the generator is fully enclosed in the steel, but a large section of the polyethylene is removed, providing a flux of sub-MeV neutrons from a wide range of angles. Neutron angular and spectral measurements using both a nested neutron spectrometer and a commercial liquid scintillator coupled with a 3 He detector show the expected softer neutron spectrum in moderated mode in good agreement with MCNP6 calculations. The gamma-ray spectrum from Excalibur is also in good agreement with MCNP modeling. Based on these findings, the future application of Excalibur in its two configurations is discussed.

Fissile material detection

Ratio of 1,275 and 186 keV Gamma Ray Count Rates for Age Determination of Material in UF 6 Cylinders

In this work, a method for using ratios of passive gamma ray measurements to determine the time since the filling of UF 6 cylinders was explored. This method uses the 186 keV gamma ray emitted by 235 U and the 1,275 keV gamma ray that results from interactions of alpha particles in fluorine. Studies were made using MCNP models of both 30B cylinders and smaller 5A cylinders. Results indicate that for enrichments at or above that of natural uranium, this approach may provide a method for determining age up to 10 years or identifying whether the age is greater than 10 years, even if the uranium enrichment is unknown. The calculated results show the same behavior as the results of a recent measurement campaign, indicating the possibility of using this method for real-world cylinders.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Feasibility of an Active Interrogation System to Classify Waste with He-4 Neutron Spectroscopy

This work investigates a 4 He-detector active interrogation system that leverages neutron spectroscopy to classify nuclear waste streams. MCNP models tested the concept through the simulation of a D-D neutron generator, an array of 4 He detectors, and various waste compositions. The fast-neutron Differential Die-Away signature was augmented with a neutron-energy discrimination signature. This signature isolates induced fission neutrons, the energy of which is greater than that of the D-D monoenergetic spectrum. With the incorporation of this spectroscopic technique, the measurement time decreased by 3–9% (depending on the degree of neutron moderation and absorption presented by the sample), demonstrating how neutron spectroscopy can enhance active interrogation methods. The reduced measurement times would have significant financial and logistical benefits for facilities with large footprints of low-level waste production.

4He detectors

Estimation of General Area Dose Rates During Closure Weld Operations on the DOE Standard Canister Pacific Basin Nuclear Conference 2024

The Department of Energy (DOE) Spent Nuclear Fuel (SNF) Road Ready Demonstration Project is an effort to package DOE-managed SNF for road-ready dry storage. A key operation in this demonstration is performing the closure weld on the DOE Standard Canister after it has been loaded with the SNF. Planning the closure weld operations is ongoing and must consider what parts of the operations (if any) can be done manually and what parts must be done remotely. A key parameter to inform the evolving planning of these operations is the dose rate in the vicinity of the canister and the dose hazard posed to the closure weld operators and collocated workers. This study developed a Monte Carlo N-Particle (MCNP) model of the closure weld operation, which was used to estimate general area dose rates as well as dose rates to the welding location, weld operators, and collocated workers. These estimates provide some insight into the potential radiation exposure hazard for workers during the operations and are a factor to be considered in the ongoing development of designs, technology, processes, and administrative controls for packaging DOE-managed SNF for road-ready dry storage.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Assessment of the Impact of Realistic Sensor Physics and the Integration of Ex-Core Sensors on Reactor Power Synthesis

In the work documented in this report, a weighting function–based core power synthesis method was applied to multiple Monte Carlo N-Particle (MCNP) reactor models, which are informed based on simulated self-powered neutron detector (SPND) responses. The weighting function method used has been coined the point-based iterative (PBI) method. The goal of this application is to assess the impact of considering realistic sensor physics in the generation of the simulated SPND outputs as well as to consider how the synthesis is impacted based on the inclusion of ex-core detectors in the model. The NuScale small modular reactor (SMR) and Westinghouse AP1000 pressurized water reactor (PWR) are the models that served as the testbeds for the assessment of realistic sensor physics; this was achieved by using Geant4 SPND models in comparison with analytical models, such that the effect of electron transport in realistic SPND geometries in the Geant4 model can be understood in terms of synthesis error and convergence time. The comparison was considered for fuel burnup–induced perturbations, for a range of sensor string densities and synthesized power distribution axial fidelities. The Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor MCNP model was used to assess the impact of ex-core sensors; this was done by performing synthesis with and without the ex-core detectors and by quantifying the synthesis error and number of iterations associated with Gaussian-type perturbations in many locations in the core. The TAMU TRIGA model was particularly pertinent for this study because of the interest in future experimental tests with SPNDs in this reactor, as well as the ease of modifying the MCNP model to include ex-core detectors with heterogeneously described response functions. Results from the comparison between the Geant4 and analytical SPND models indicate that similar average and maximum synthesis errors were obtained for burnup-induced perturbations in both the NuScale SMR and the AP1000. This was true for a range of sensor string densities and axial fidelities. However, there were marked differences between both the Geant4 and analytically informed models in terms of the iterations required to converge on the synthesized power distribution. Namely, the Geant4-informed models tended to lead to fewer iterations, except for a few sensor–core configurations that had particularly numerous iterations. Results from the ex-core sensor assessment with the TAMU TRIGA model indicate that the inclusion of ex-core sensors drastically reduces the synthesis error of Gaussian-type perturbations close to the edge of the core, and it slightly reduces synthesis errors for perturbations closer to the center of the core. This was achieved with a minimal increase in computational cost—that is, the number of iterations required for convergence. The errors were identified to be in the same location as the perturbation in the core, indicating that the methodology remains robust for unperturbed regions of the core. A secondary result from this study with the TAMU TRIGA was yielded by analysis of the neutron flux levels in the in-core and ex-core sensor locations of the core; these flux levels indicate that SPNDs could be used as both in-core and ex-core sensors, so long as the emitter material is sensitive to thermal neutrons. The results from these studies provide a quantitative understanding of the importance of considering realistic sensor physics and including ex-core sensors to perform accurate and timely power distribution synthesis of a reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS

NCERC-CL TREAT Free Field Characterization Foil Reaction Rate Results

This report summarizes the analysis and preliminary MCNP reaction rate modeling of free field characterization measurements performed in the Transient Reaction Test Facility (TREAT) at Idaho National Laboratory (INL) between 12/18/2023 and 2/22/2024. A series of seven irradiations were performed using the Big-BUSTER (Broad Use Specimen Transient Experiment Rig) core configuration. Foil sets including 19.5% LEU-Zr alloy wire, S, Au, Fe, Ni, Co, Ti, and Zr were deployed in each irradiation. Identical foil sets were supplied to Los Alamos National Laboratory (LANL), Lawrence Livermore National Laboratory (LLNL), and INL. All foils were supplied by INL. Additional details on the experiment can be found in [2]. The foil sets were shipped to the National Criticality Experiments Research Center Counting Laboratory (NCERC-CL) branch at LANL (NCERC-CL NISC) and were received on 3/25/2024. Additional details on receipt and counting are reported in [3]. All reactor dosimetry measurements were performed in adherence to the ASTM standards applicable to reactor dosimetry. Measured reaction rates from Au, Co, Fe, Ni, and Ti are reported. Select reaction rate ratios and simulated reaction rates and ratios are discussed. Reaction rates for all measured reaction products were modeled with MCNP 6.3.1 and a TREAT input deck supplied by Edward Lum. Neutron emission estimates calculated using the measured and simulated reaction rates show excellent agreement between the three fast-threshold n-p reactions with the percent differences being < 9% between the three reactions. The capture reactions had poorer, but still reasonable, agreement with < 20% percent differences between the neutron emission estimates of the three capture reactions. The percent difference between the fast-threshold and capture reactions is ∼ 100% for each reaction. This supports the TREAT model is not accurately modeling the neutron spectrum and additional measurements are required to characterize TREAT and match measurements to simulation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Evaluation of simulated HPGe detector efficiencies in OpenMC compared to MCNP

In late 2023, the open-source radiation transport code OpenMC introduced a pulse-height tally (PHT) feature, enabling users to track the total energy deposited by individual photons in OpenMC cells. This function represents an important improvement of the OpenMC code because it allows users to simulate the response function of a gamma detector without relying on closed-source alternatives. Despite this, limited work has been published evaluating OpenMC as a radiation transport code for simulating gamma spectroscopy experiments. This study attempts to demonstrate the usefulness of OpenMC in this space by directly comparing its PHT output to MCNP (a trusted industry-standard Monte Carlo code). In the first half of this study, a radiation transport experiment is described in which a detector with a complex internal geometry is exposed to various gamma-emitting isotopes ( 133 Ba, 137 Cs, and 60 Co) over a range of distances. The setup of this experiment was modeled in MCNP with sufficient detail to capture the efficiency characteristics of a high purity germanium (HPGe) detector for the primary gammas of each isotope. After good agreement between the MCNP model and experiment was achieved, an identical model was produced in OpenMC to allow direct comparison between the full energy peak (FEP) values produced in OpenMC and MCNP. The results show strong agreement between OpenMC and MCNP across the full range of tested energies, with each model’s FEP values typically within 2% of each other and most FEP areas within 4% of experimental data. Worse agreement was seen between the Monte Carlo codes and experiment below 300 keV (an expected result). For the 662 keV line of 137 Cs, both codes were found in poor agreement with experiment and each other over the full range of distances tested (possibly indicating an error with the 137 Cs experimental data). Ignoring the anomalous results of the 662 keV line, all other data show good qualitative and quantitative agreement between MCNP and OpenMC. This result demonstrates the accuracy of OpenMC’s PHT feature for spectroscopic applications in which detector efficiency is a primary concern.

07 - ISOTOPES AND RADIATION SOURCES

Developing a Nuclear Quality Assurance Compliant Design Methodology for Neutronic Analysis of Xe-100 Design

The primary objective of this work is to develop a design methodology compliant with nuclear quality assurance standards for the Xe-100 neutronic design verification studies. To achieve this, a Monte Carlo model of the Xe-100 reactor was constructed using the exclusion principle, transformation technique, and universe-based level specification following Idaho National Laboratory (INL) NQA level-1 compliant standards and an NQA-1 compliant version of MCNP6. The model encompasses the entire reactor core structures, including the upper plenum, core region, and lower plenum sections, along with all sub-components. The active core section was represented using the spectral regions, each comprising a particular fuel composition and temperature averaged over the considered zone, calculated by X-energy using Very Superior Old Programs (VSOP). Additionally, a component-wise temperature map was implemented into the model, not only for the core region but also for the structural components. Temperature-dependent cross-section libraries, along with thermal scattering law libraries, generated using INL NQA-1 compliant version of NJOY21, were utilized for each isotope in the burnt fuel and the structural materials. Furthermore, the volume of each modeled component was estimated using a stochastic approach with the ray tracing method in MCNP and criticality calculations were performed.

22 GENERAL STUDIES OF NUCLEAR REACTORS

STARTR: An Open-Source MARVEL model for the NRIC Virtual Test Bed [Poster]

The National Reactor Innovation Center (NRIC) seeks to improve the understanding of microreactor physics in industry and academia through the development of a Microreactor Applications Research Validation and Evaluation (MARVEL) reactor-based model, published on the Virtual Test Bed (VTB). To achieve this goal, the Sodium-cooled Thermal-spectrum Advanced Research Test Reactor (STARTR) model was built using publicly available MARVEL specifications where possible and approximations where applicable, and was optimized for fast runtimes for researchers to receive rapid simulation feedback. STARTR will fill a gap between stakeholder interest and available models, as the first Sodium-cooled Thermal Reactor (STR) hosted on the VTB with baseline performance sanctioned by INL. This project involved the definition of all materials used in the reactor, geometry and all reactor subcomponents, and assertion of tallies and simulation settings within OpenMC 0.13.3. This poster details a small subset of the overall reactor physics testing: the two-dimensional power peaking factors and the flux energy spectrum, as well as plots of the created geometry. Future work includes code-to-code verification between the OpenMC-based model and a separately designed MCNP 6.2-based model.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN

Benchmarking Monte Carlo codes for the modelling of low-energy neutron production target reactions

The increasing adoption of accelerator-based neutron sources (ABNS) for applications including neutron capture therapy (NCT) research has highlighted the need for accurate simulation tools. Precise modelling of the neutron production target is crucial to ensure that simulated predictions of neutron beam characteristics used for subsequent beam shaping assembly design are reliable. This work presents a comprehensive benchmarking of four widely-used Monte Carlo codes - Geant4, PHITS, FLUKA (CERN), and MCNP - for modelling low-energy neutron production target reactions. Using their recommended physics models and cross-section libraries, we evaluate each code’s performance in simulating four beam-target reactions: 7 Li(p,n) 7 Be, 9 Be(p,n) 9 B, 9 Be(d,n) 10 B, and C(d,n)N. Predictions of neutron yield, angular distributions, and energy spectra are compared against available thick target experimental data. Results show varying levels of agreement between the codes depending on the reaction type, energy range, and beam characteristics. Geant4, MCNP and PHITS are the overall best performing codes for the simulation of total neutron yield and yield in the forward direction across most reactions. Across energies where experimental benchmarks exist, inter-code discrepancies in total and forward-directed yield are typically 10 to 30%, with larger deviations at near-threshold incident ion energies. PHITS provides the best overall reproduction of experimental spectra, particularly for the 9 Be(p,n) 9 B reaction. Additionally, PHITS demonstrates superior computational performance for most reactions. These findings provide valuable guidance for ABNS design, highlighting the strengths and limitations of each code for the simulation of low-energy neutron production reactions.

43 PARTICLE ACCELERATORS

Analysis of SPND Material Candidates for Real-Time Neutron Spectrum Monitoring Poster

The future of the nuclear industry will be shaped by advanced reactor designs requiring innovative instrumentation for safe operation. A key area of innovation is the ability to accurately and reliably monitor the neutron energy spectrum during operations, which would also enhance the experimental capabilities of existing research reactors. To address this need, our research focuses on the potential use of Self-Powered Neutron Detectors (SPNDs) for real-time neutron spectrum monitoring. We investigated various candidate materials for SPNDs by evaluating their energy-dependent reaction rates. Promising materials were further modeled using MCNP to determine the free electron production per neutron across the neutron energy spectrum, approximating their energy-dependent signal sensitivity. Our findings identified Samarium (Sm) and Gadolinium (Gd) as top candidates for thermal neutrons, Hafnium (Hf) and Thulium (Tm) for epithermal neutrons, and Tantalum (Ta) for fast neutrons. These materials will undergo future experimental validation to confirm the modeled responses and assess their feasibility for real-time neutron spectrum measurements.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN

Concrete Compositions Used in Neutronics Analyses for Design of the Second Target Station Project

The purpose of this letter is to document the concrete compositions used by the Second Target Station (STS) Project Neutronics Group. Members of the Neutronics Group all have access to these compositions in our MCNP Master Model, which is stored on the ORNL GitLab server. It is expected that all STS Neutronics staff members will use these compositions in their analyses, unless all stakeholders agree to an exception. It is impossible to know the detailed composition of concrete before its constituents have been purchased, mixed, and a chemical analysis performed. Therefore, the STS Neutronics Group decided to use the same concrete compositions used by the Neutronics Group at the First Target Station (FTS). Furthermore, this leads to consistency between the analyses performed by the FTS and STS Neutronics Groups. The most important parameter to match between design simulations and construction regarding these concrete compositions is the density. Next in importance is the fraction of iron and hydrogen in the mixture.

43 PARTICLE ACCELERATORS

Benchmarking a Defeatured Geant4 NIF Model with HTOAD and HNED Neutron Data

This study benchmarks a simplified Geant4 NIF Target Chamber (TC) model against foil measurements in two instruments: the HTOAD and HNED Snout. The number of product atoms per source neutron per gram ( N 0 /n/g ) is compared across multiple locations and material configurations. The model reproduces high energy threshold reactions dominated by 14 MeV neutrons to within a few percent, depending on location. Discrepancies occur where scattering and moderation are significant for low energy threshold reactions. The defeatured Geant4 model achieves comparable statistical precision with ∼400 times less CPU time than a full fidelity TC model in MCNP. The results can be produced locally on a laptop without the need for high performance computing. Low energy biasing can be improved by preserving room return pathways and implementing IRDFF cross section evaluations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Processing MCNP Elemental Edit Outputs

The Monte Carlo N-Particle (MCNP) transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. A UM geometry is a collection of elements representing a solid geometry. The first step of MCNP UM modeling is using other software packages to create a finite element mesh representation of a solid 3D geometry. Computer-aided design (CAD) or computer-aided manufacturing (CAM) software is typically used to create a solid geometry model, which is later imported into mesh generation software to create a UM model. The MCNP UM feature was originally designed for models generated by the Abaqus/CAE software. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. MCNP can process a UM model consisting of several different element types including linear tetrahedral or hexahedral elements and calculate quantities of interest such as flux and energy deposition at elements. An MCNP UM simulation provides high-fidelity elemental edit (i.e., tally) outputs, which can be further used in multiphysics calculations. The MCNP UM feature was used for multiphysics simulations where quantities of interest calculated by MCNP are used as inputs for heat transfer calculations in Abaqus. MCNP6.3 can produce two types of elemental edit output (EEOUT) file formats: ASCII and HDF5. An EEOUT file type must be requested on an EMBED card while output type (flux or energy deposition) must be requested on an EMBEE card. We wrote Python3 scripts to extract energy deposition values in an ASCII or HDF5 EEOUT file and compute a heat flux profile for an Abaqus heat transfer calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

AGR 1 and 2 Isotopic Depletion Validation

This presentation includes work performed developing a computational benchmark for the AGR 1 and 2 experiments. These experiments involve irradiation of TRISO particle compact in the Advanced Test Reactor (ATR). Many post-irradiation measurements were made from these experiments, including isotopic concentrations of fission products. The benchmark involves a MCNP-ORIGEN coupled model of the ATR and AGR experiment capsules to calculate a comprehensive inventory of fission product and actinide concentrations in the compacts. This presentation describes the experiments, the measurements obtained, the benchmark model, and results of the model evaluation.

97 - MATHEMATICS AND COMPUTING