Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MCNP”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An update to δ-ray production in MCNP6

As originally introduced in MCNP6.2, the δ-ray production subroutine had a number of flaws related to treatment of particle or nuclear spin contributions to the production cross section. A recent update by the author to the MCNP6 code corrects these flaws. This report documents those changes for the benefit of MCNP6 users who use the delta ray production capability. Additionally, this report aims to provide a more comprehensive documentation of the delta ray production, building on the prior documentation by C. Anderson and coworkers. On this latter point, the goal is to save future MCNP developers from a visit to the library to dig out ancient texts from 1952, or at least to provide useful points of reference should such a visit be necessary.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2002

Benchmark study for total enery electrons in thick slabs

The total energy deposition profiles when highenergy electrons impinge on a thick slab of elemental aluminum, copper, and tungsten have been computed using representative Monte Carlo codes (NOVICE, TIGER, MCNP), and compared in this paper.

total enery deposition MCNP NOVICE TIGER↗
DOE OSTI2024

Impact of Higher Fidelity Design Iterations on Critical System Criteria

Nuclear criticality experiments are effective at informing the performance of nuclear data libraries across many applications. This work explores the implications of refining critical experiment MCNP models from their low fidelity optimization phase to penultimate neutronic models. Specifically, this work is focused on two series of plutonium fueled experiments funded through internal programs at Los Alamos National Laboratory building off previous efforts under the EUCLID (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) collaboration. Thales, the first of the two collaborations, is a fast spectrum Ta-reflected plutonium experiment to support operations at PF-4. The second experiment are twin configurations designed to target the intermediate energy cross sections in 239 Pu. Motivation for this experiment stems from the PARallel Approach of Differential and InteGral Measurements (PARADIGM) collaboration which hopes to achieve a significant reduction in 239 Pu cross section uncertainties in the intermediate region.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

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↗
DOE OSTI2025

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↗
DOE OSTI2024

Quantifying the Potential of Argon Detection Capabilities for Nuclear Explosion Monitoring

Abstract Current noble gas detection systems for nuclear explosion monitoring are based on the detection of four radioxenon isotopes—Xe-131m, -133, -133m and -135. The data provided by radioxenon detection could be enhanced by other radionuclide signatures such as Ar-37. Activation of Ca-40 in rock by neutrons produces Ar-37, and monitoring for this additional nuclide could help distinguish detections of nuclear explosions from background sources, such as medical isotope production. This work studies the capabilities of a hypothetical argon detection network. A 10 kt explosion was modeled using MCNP and SCALE to determine the inventory of Ar-37 created in a representative granite rock layer, assuming either 0.1, 1 or 10% of the total inventory was released. The Ar-37 inventory was combined with atmospheric transport data from HYSPLIT compiled in a previous study, along with the detection limits of standard Ar-37 detection systems, to determine how many hypothetical monitoring stations would detect Ar-37 from an explosion. This method was repeated for 365 HYSPLIT data sets to create a year’s worth of hypothetical explosions, releases, and detections. The study quantified the average number of detections per release, the number of stations detecting Ar-37, and the possibility of detecting Ar-37 in coincidence with xenon.

37Ar↗
DOE OSTI2025

Quantifying uncertainty in uranium concentration measurements via K-edge densitometry

This study quantifies the uncertainty in uranium concentration predictions of fluoride and chloride-based salts within a steel pipe using K-edge densitometry. Modeling and simulation was conducted with the Monte Carlo N-Particle Transport (MCNP) code. The quality of of this technique’s prediction in a pipe requires proper characterization of the pipe’s thickness, which is dependent on the source size and axial offset from the pipe centerline. The thickness was determined as either the center-line thickness seen by the X-ray source or an average value determined through random sampling. Generally, the predicted concentrations were slightly better at lower offset with the random sampling thickness and using the center-line thickness for the highest offsets. For a line-beam source and varying axial offsets, the relative error of concentration was within 1% of the true value but uncertainty increased by 2 orders of magnitude. Similarly, for no axial offset, the relative error was significantly less than 1% while no trend for uncertainty was found. However, at the largest possible offset for a given source size, the concentrations become erroneous and greater than the allowable 1% relative error. Furthermore, high offsets tended to increase the variance of the transmission spectra by 3 orders of magnitude.

Characterization and Analytical Technique↗
DOE OSTI2026

Investigation of irradiation damage and heat deposition: a comparative analysis for HEU-to-LEU conversion in HFIR

The planned conversion of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel requires detailed evaluation of experiment-relevant parameters to ensure continued performance for materials testing and isotope production. Here, this study presents the first comprehensive assessment of displacements per atom (dpa) and heat deposition rates in target materials within the HFIR flux trap with both HEU and candidate LEU core configurations. Seven analyses were conducted to evaluate key performance metrics, including fast neutron flux distribution, cross section response functions, cross section data, and local dpa and heat deposition rates using mesh- and cell-based tallies. Simulations employed Shift, Monte Carlo N-Particle (MCNP), and the HIFR Controller (HFIRCON) tool suite for high-fidelity transport and depletion modeling. The LEU designs—using U 3 Si 2 -Al dispersion fuel and operating at 95 MW—were compared to the current 85 MW HEU configuration. Results show that while the candidate LEU cores exhibit higher dpa rates due to a harder spectrum and extended cycle lengths, they also demonstrate reduced heat deposition rates in irradiation experiments, primarily due to increased gamma self-shielding from higher 238 U content in the core. These findings confirm that LEU conversion can maintain HFIR’s materials irradiation capabilities but may require redesigning existing experimental hardware.

HEU↗
DOE OSTI2025

Porosity in nuclear graphite and its impact on nuclear reactor science and criticality safety applications

Porosity in nuclear-grade graphite significantly influences its low-energy neutron scattering, yet its effect on underlying phonon properties remains debated. Here, this work integrates inelastic and small-angle neutron scattering (INS/SANS) experiments, advanced atomistic simulations with a novel machine-learned potential (DeepMD), total cross-section measurements, and neutronics calculations (SCALE, MCNP, OpenMC) to investigate porosity’s impact on neutron thermalization. INS measurements on diverse graphite grades reveal no discernible porosity effect on phonon spectra, which align with crystalline graphite. Conversely, total cross-section data below ≈10 meV show increased scattering attributable to SANS. Our DeepMD simulations demonstrate that realistic micropores do not distort phonon spectra, challenging the assumptions in current ENDF/B-VIII.1 porosity thermal scattering laws (TSLs). These TSLs, based on random atom removal, produce unphysical phonon spectra and inflate inelastic cross-sections. Augmenting a crystalline TSL with an SANS component accurately captures experimental total cross-sections. Neutronics benchmarks (ICSBEP/IRPhE) show ENDF porosity TSLs unphysically increase neutron multiplication factor, keff. Crucially, incorporating SANS physics (NCrystal/OpenMC) indicates accurately modeled porosity negligibly affects keff, reactor physics, or criticality safety.

Critical benchmarks↗
DOE OSTI2024

Applied radiation tags for warhead dismantlement transparency

Here, this study explores the feasibility of applying radiation tags as a Chain of Custody tool to confirm the dismantlement of individual warheads by investigating time dependent emissions from special nuclear material. Precedents of photon and neutron sources that have been used to irradiate warhead components are presented. Measurements of irradiated Highly Enriched Uranium (HEU) are used to benchmark simulations with the CINDER90 library for MCNP. Delayed gamma emissions are simulated as a function of time since irradiation of moderated HEU and Weapons Grade Plutonium with photon and neutron sources. Models of commercial 9 MV and 15 MV linacs were found to induce ~10 9 fissions/g and ~10 11 fissions/g, respectively, creating 3 to 4 orders of magnitude more delayed gammas than neutron sources, modeled at 10 10 n/s. However, a DT neutron source was found to induce sufficient fissions, ~10 7 fissions/g, to create quantifiable signatures relative to intrinsic emissions. Applied radiation tagging is determined to be technically feasible, but its operational context and practical implementation requires significant development.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2025

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↗
DOE OSTI2025

Optimizing a detection system for fissile material in nuclear disarmament verification

In arms control treaties, verification plays a crucial role in detecting non-compliance, deterring future violations, and building trust between state parties. Neutron interrogation that induces fission reactions in fissile isotopes and measures the resulting fission neutrons, could be employed for this purpose. This study aims to develop a system which can determine the presence of fissile material while intrinsically protecting information. In this paper, we focus on optimizing the system for discriminating between an enriched uranium block from a depleted uranium (DU) block. The system was built and we report on benchmark measurements with DU and 16% enriched uranium blocks. Furthermore, the Excalibur (Experiment for Calibration with Uranium) neutron source, a neutron spectrometer (redBubble Technology Industries (BTI) N-Probe), and superheated droplet detectors were used for these measurements. MCNP simulations provided insights into detector responses to fissile materials with varying isotopic compositions, confirming that the system functioned as designed.

Active neutron interrogation↗
DOE OSTI2026

Fuel Fabrication Specification Impact Analysis for NBSR LEU Conversion

As part of a national initiative to enhance nuclear security and reduce proliferation risks, significant efforts have been undertaken by the National Nuclear Security Administration Material Management and Minimization Office of Reactor Conversion Program to convert U.S. high performance research reactors (USHPRRs) from the use of highly enriched uranium (HEU) to low-enriched uranium (LEU), including the National Bureau of Standards Reactor (NBSR). The current plan is to procure LEU fuel assemblies from commercial fabricators according to fuel specifications tailored for each USHPRR. The analysis conducted at Brookhaven National Laboratory was part of an effort to identify the sources of uncertainty in the fuel specifications that may impact the performance of the NBSR core after its conversion and, in particular, to assess the range of acceptable tolerance limits from the perspective of core safety and reactor performance. Using the stochastic neutronics code MCNP 6.2, the variations in important NBSR neutronics characteristics were analyzed as a function of the specification parameters independently and in combination. The important NBSR specification parameters analyzed were the fuel isotopic composition, the amount of impurity content in cladding, the fuel plate thickness, and the fuel element 235U mass loading. The range of variation of each specification parameter was based on the technical specification limit or available as-fabricated assay data and uncertainties. The NBSR neutronics characteristics selected for analysis were the reactor reactivity characteristics at equilibrium and the equilibrium fuel cycle length. Results show that with variations in the fabrication parameters of the as-fabricated U-10Mo fuel within the specification limitations, the excess reactivity of the NBSR LEU core remains well below the 15% Δk/k technical specification limit, and the shutdown margin is always significantly greater than the required 0.68% Δk/k. This ensures that the NBSR can be operated safely and reliably shut down for all analyzed cases within the specified fabrication limits after the LEU conversion. In the prototypic case, the fuel cycle length was 1.5 days longer than the targeted 38.5 days. In a credible worst-case scenario, where all low-reactivity parameters were combined, the fuel cycle length was reduced to 35.5 days, which is still considered manageable for reactor operations. Variations in cycle length are primarily driven by changes in 235U loading, with other parameters having secondary effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

GEAR-MC and Differential-Operator Methods Applied to Electron-Photon Transport in the Integrated TIGER Series

The sensitivity analysis algorithms that have been developed by the radiation transport community in multiple neutron transport codes, such as MCNP and SCALE, are extensively used by fields such as the nuclear criticality community. However, these techniques have seldom been considered for electron transport applications. In the past, the differential-operator method with the single scatter capability has been implemented in Sandia National Laboratories’ Integrated TIGER Series (ITS) coupled electron-photon transport code. This work is meant to extend the available sensitivity estimation techniques in ITS by implementing an adjoint-based sensitivity method, GEAR-MC, to strengthen its sensitivity analysis capabilities. To ensure the accuracy of this method being extended to coupled electron-photon transport, it is compared against the central-difference and differential-operator methodologies to estimate sensitivity coefficients for an experiment performed by McLaughlin and Hussman. Energy deposition sensitivities were calculated using all three methods, and the comparison between them has provided confidence in the accuracy of the newly implemented method. Unlike the current implementation of the differential-operator method in ITS, the GEAR-MC method was implemented with the option to calculate the energy-dependent energy deposition sensitivities, which are the sensitivity coefficients for energy deposition tallies to energy-dependent cross sections. The energy-dependent cross sections could be the cross sections for the material, elements in the material, or reactions of interest for the element. Further, these sensitivities were compared to the energy-integrated sensitivity coefficients and exhibited a maximum percentage difference of 2.15%.

42 ENGINEERING↗
DOE OSTI2025

Demonstration of TOFFEE: A Response Uncertainty Quantification Tool

A key characteristic in neutron transport is nuclear data. Cross-section uncertainty is not used in MCNP6.3 to propagate response uncertainty without external analysis. Here, the TOol For Fast Error Estimation (TOFFEE) is a Python-based code developed to automate the propagation of cross-section uncertainty for MCNP evaluations. TOFFEE implements the sandwich rule to calculate the uncertainty from cross sections with sensitivity coefficients from MCNP6.3 and ENDF/B covariance data. In this paper, TOFFEE has been tested with benchmark experiments, and it has been compared to the uncertainty quantification capabilities of Sampler and TSUNAMI, within SCALE, to verify the application’s capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Challenges and Opportunities Specific to Microreactors: Compact and Light-Weight Shielding During Operation and Transportation

The problem of designing a shield for a mobile microreactor has been addressed. The principal challenges encountered in this endeavor are the dual goals of increasing radiation protection and decreasing weight. It is recognized that these two goals are contradictory. Increasing shielding naturally increases weight. Conversely, decreasing weight naturally calls for decreasing material inventory, including possibly decreasing shield mass. The operational and regulatory drivers for shielding and transportation constraints are identified and summarized. The main constraints pertain to volume, mass, radiation dose, material feasibility and availability, and decommissioning expectations. The latter are limits on residual contamination and dose rates after the reactor is removed. Progress toward solving the challenge of shielding a transportable microreactor is made by considering a tungsten tetraboride (WB4) shell surrounding the reactor vessel. The shield and a generic microreactor were modeled using the MCNP code. It was determined that the presence of a 16-cm-thick WB4 shield reduces the dose equivalent rate during operations to 3.6 mrem/h just outside of a concrete vault and dome and soil berm surrounding the reactor, which is very near the occupational dose rate limits. Without the WB4 layer, the dose rate at the same location outside the berm would be 120 mrem/h. Most importantly, the presence of the WB4 shield substantially protects the surrounding concrete and soil from activation, thus decreasing the duration of administrative control for the site following removal of the reactor from 16 to 3.4 years. At that point, the site may be decommissioned and returned to unrestricted public access. A limitation on the incorporation of a WB4 shield as an intrinsic component of a microreactor is that it results in a design that is transportable only by using the largest road conveyances and that the load would require special permits for transit through the U.S. highway system. New directions for further improvements and possibly more effective solutions to the shielding problem of a mobile microreactor are discussed in the conclusion section.

Microreactor↗
DOE OSTI2025

Computational Optimization of 133m Xe Production via Neutron Irradiation in a TRIGA Reactor

Here, the Comprehensive Nuclear-Test-Ban Treaty bans all nuclear tests worldwide. As part of treaty compliance, the concentration of radioactive nuclides in the atmosphere is monitored to detect nuclear weapons tests. Radioactive noble gas fission products, specifically radioxenon, can vent into the atmosphere after a nuclear weapons test, even if the test is well contained underground or underwater. Radioxenon thus serves as a signal for nuclear weapons tests. All atmospheric monitoring systems require samples of radioxenon isotopes for detector calibration, quality control, and certification. Here, we present a novel, improved method for creating samples of 133m Xe via neutron irradiation of 132 Xe in the Washington State University TRIGA reactor. 132 Xe neutron absorption results in either 133 Xe or 133m Xe—thermal neutron absorption results in 133m Xe 12% of the time, while fast neutron absorption (above ~1 MeV) results in 133m Xe ~50% of the time. To optimize the production of 133m Xe via neutron absorption in 132 Xe in the thermal TRIGA reactor, spectral tuning using an irradiation chamber is required to maximize the fraction of fast neutrons being absorbed and minimize the number of thermal neutrons interacting with the 132 Xe. We used MCNP to tally 132 Xe absorptions with the isotopic tally function, flux tallies and neutron attenuation to estimate the number of neutrons reaching the 132 Xe through the irradiation chamber, and the adjoint importance function to improve the source strength estimate. Additionally, we performed a heat transfer analysis for safety considerations. It was determined that the use of a 96% enriched 10 B boron carbide chamber, placed next to the fuel elements in reactor position D8, increases the 133m Xe/ 133 Xe activity ratio from a baseline value of 0.3 to 1.0, a 233% increase. Additionally, it was determined that the alpha heating produced in the boron does not become an unmanageable problem in the Washington State University reactor.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗