Search NASA⌕ Search

SEARCH · Search NASA

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Simulated Performance of the Micro-Pocket Fission Detector in the Advanced Test Reactor Critical Facility

The Micro-Pocket Fission Detector (MPFD) is a small-form–factor real-time fission chamber. MPFD performance has been simulated in the Advanced Test Reactor Critical Facility (ATRC), located at Idaho National Laboratory (INL). Here, the neutron and gamma-ray flux profiles and magnitudes were simulated using Monte Carlo N-Particle (MCNP) in the near-core B-8 irradiation position. These simulations were performed at 69 discrete axial locations inside the B-8 position for three separate orientations of the nearby hafnium outer shim control cylinders and at a power level of 700 Wth. The resulting neutron and gamma-ray flux values were used to determine the MPFD response for various fissile masses and detector gas pressures. The optimal gas-operating pressure was determined to be between 30 and 60 psig. The required fissile-layer mass was determined to be between 1–2 µg of 235 U. Additionally, the gamma-ray to fission-fragment interaction rate was determined to be 4.42 × 10 -3 with average energy deposition for gamma rays and fission fragments in 30 psig argon gas to be 1 keV and 9.5 MeV, respectively.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
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↗
DOE OSTI2026

Expanded Intercomparison of Nuclear Data Libraries Using Jupiter and Jupiter High-240 Experiments

There is a limited availability of plutonium experiments with sensitivity to lead in the ICSBEP (International Handbook of Evaluated Criticality Safety Benchmark Experiments) Handbook. The Jupiter and Jupiter High-240 experiments were performed at the National Criticality Experiments Research Center as a collaborative effort between Los Alamos National Laboratory and the Japan Atomic Energy Agency to assess lead void coefficients in a plutonium-lead system containing weapons- and reactor-grade plutonium, respectively. Concurrent with benchmark development, an intercomparison of calculations using different nuclear data libraries has been performed to assess the usability of the experimental data for nuclear data adjustment in a “softer-that-fast” neutron energy spectrum. Eigenvalue calculations using MCNP with the ENDF/B-VIII.0 and TENDL-2021 nuclear data libraries calculate closest to the benchmark values for Jupiter. Calculations using JENDL-5 and ENDF/B-VIII.1 match best with the Jupiter High-240 values. Lead void worth calculations using the various nuclear data libraries are all within 3σ of their respective measured values. Perturbation studies between ENDF/B-VIII.0 and ENDF/B-VIII.1 demonstrate an approximate increase in calculated eigenvalues for the Jupiter series experiments by ~240 pcm for plutonium (mostly 239 Pu) and ~120 pcm for lead accompanied by a decrease contributed by ~113 pcm for copper and ~13 pcm for stainless steel. Nuclear data sensitivities and uncertainties investigated using Whisper show slightly lower sensitivity to scatter than a lead-reflected plutonium sphere but greater sensitivity to neutron capture. The sensitivities between Jupiter and Jupiter High-240 for lead are very similar for both ENDF/B-VIII.0 and ENDF/B-VIII.1 nuclear data. These benchmarks are more sensitive to neutron capture in lead than other plutonium benchmark experiments and would be useful for both lead and 240 Pu validation. In conclusion, with the high degree of compensating effects between copper, lead, and plutonium cross sections, additional isolated Pb-Pu and Cu-Pu benchmarks would be beneficial in improving these nuclear data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

MHD Analysis of Dual-Coolant Lead-Lithium Blanket for Spherical Tokamak Advanced Reactor

The tritium breeding blanket is vital for future fusion power plants, with the Spherical Tokamak Advanced Reactor (STAR) project highlighting the dual-coolant lead-lithium (DCLL) design. The DCLL blanket performs shielding, energy exhaust, and tritium breeding using a lead-lithium alloy , with lithium as the breeder and lead as the neutron multiplier. It also serves as the primary coolant, with helium providing supplemental cooling. Reduced-activation ferritic/martensitic steel is used for the blanket structure. Magnetohydrodynamic (MHD) phenomena influence the liquid metal flow in a magnetic field, affecting heat transfer in the breeder affected by energetic neutrons. Understanding key flow parameters in such conditions is critical for efficient DCLL design. This study uses three-dimensional thermofluid MHD analysis with ANSYS CFX software, modified at Princeton Plasma Physics Laboratory, to simulate high Hartmann flows. The neutronics code MCNP, coupled with plasma equilibrium, provides heat source distribution. In conclusion, we examine electromagnetic interactions in adjacent fluid domains and analyze the magnetic field’s impact on flow distribution in the inboard and outboard blanket layout, using detailed mesh generation for accurate results.

DCLL↗
DOE OSTI2024

Transfer functions for Q A /Q B international regulatory limits for the safe transport of radioactive materials

This paper presents a proposed revision of the International Atomic Energy Agency transport regulations, related to the A 1 and A 2 limit values used to determine the radioactive transport classification. Based on the 'Q system', a novel methodology was introduced to derive Q A and Q B values related to scenarios involving external exposure from a distant source. These values are key parameters that respectively represent the total effective dose and total equivalent dose to the skin, from all primary and secondary particles contributing to radiation exposure. The International Working Group (WG A 1 /A 2 ) is established and associated with the TRANSSC Technical Expert Group on Radiation Protection. A review of the A 1 and A 2 values is performed in response to identified limitations within the existing Q system. The followed approach is based on Monte Carlo simulations that enabled the development of transfer functions aimed at reducing computational time and increasing the flexibility of dose evaluations for any radionuclide with known particle emission spectra. This method allows updating the Q A and Q B values to account for future data evolutions (decay data, fluence-to-dose conversion coefficients) and standardizing the calculation of regulation limits across all referenced radionuclides and scenarios related to external exposure. The transfer functions are established using three Monte Carlo simulation codes—FLUKA, Geant4, and MCNP—and address the previous limitations of the 'Q system', reflecting the latest International Commission for Radiation Protection recommendations and improvements in calculation techniques. The results of the WG show consistent agreement across the codes, with minor discrepancies observed at low primary energies due to statistical uncertainties and different handling of stopping power for electrons/positrons in the codes. This revised approach aligns with current standards and recommendations, ensuring that the radiological consequences of transport accidents are acceptable for the new A 1 and A 2 limits from a radiological protection perspective.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Experimental validation of an integrated modelling approach to neutron emission studies at JET

An integrated modelling methodology for the calculation of realistic plasma neutron sources for the JET tokamak has been developed. The computational chain comprises TRANSP plasma transport and DRESS neutron spectrum calculations, and their coupling to the MCNP neutron transport code, bridging plasma physics and neutronics. Here, we apply the developed methodology to the analysis of neutron emission properties of deuterium and helium plasmas at JET, and validate individual modelling steps against neutron diagnostic measurements. Two types of JET discharges are modelled—baseline-like and three-ion radio-frequency scenarios—due to their diversity in plasma heating, characteristics of the induced fast ion population, and the imprint of these on neutron emission properties. The neutron emission modelling results are quantitatively compared to the total neutron yield from fission chambers, neutron emissivity profiles from the neutron camera, neutron spectra from the time-of-flight spectrometer, and neutron activation measurements. The agreement between measured and calculated quantities is found to be satisfactory for all four diagnostic systems within the estimated experimental and computational uncertainties. Additionally, the effect of neutrons not originating from the dominating D(D, n) 3 He reactions is studied through modelling of triton burnup DT neutrons, and, in mixed D- 3 He plasmas, neutrons produced in the 9 Be(D, nγ) 10 B reaction on impurities. It is found that these reactions can contribute up to several percent to the total neutron yield and dominate the neutron activation of samples. Finally, the effect of MeV-range fast ions on the neutron activation of 115 In and 27 Al samples is measured and computationally validated.

integrated modelling↗
DOE OSTI2022

AR reactor neutronics multi-code validation*

The affordable, robust, compact (ARC) reactor is a tokamak fusion reactor concept currently under development by Commonwealth Fusion Systems and Massachusetts Institute of Technology. There are three important neutronics considerations for the operation of the ARC reactor: (1) breeding of enough tritium in the blanket to sustain the D–T reaction in the plasma; (2) ensuring low fluence on the superconducting toroidal field coils; and (3) assessing neutron volumetric heating in structural components. Here, this work aims to perform a validation of the neutronics analysis approach by code-to-code comparison. State-of-the-art software stacks are employed for the neutronics analysis of the ARC reactor, and a computer-aided design (CAD) model is used directly for Monte Carlo (MC) neutron transport calculations. Three software stacks, Attila-MCNP, OpenMC-DAGMC, and Shift-DAGMC, are used to perform neutronic analyses of a 90° sector CAD model of the ARC reactor. In conclusion, results show that the flux tallies calculated by the three software stacks are very close. Also, the volumetric heating and tritium breeding values have less than 0.6% relative difference between codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulations of stand-off runaway electron beam termination by tungsten particulates for tokamak disruption mitigation

Abstract Stand-off runaway electron termination by injected tungsten particulates offers a plausible option in the toolbox of disruption mitigation. Tungsten is an attractive material choice for this application due to large electron stopping power and high melting point. To assess the feasibility of this scheme, we simulate runaway collisions with tungsten particulates using the MCNP program for incident runaway energies ranging from 1 to 10 MeV. We assess runaway termination from energetics and collisional kinematics perspectives. Energetically, the simulations show that 99% of runaway beam energy is removed by tungsten particulates on a timescale of 4–9 µ s. Kinematically, the simulations show that 99% of runaways are terminated by absorption or backscattering on a timescale of 3–4 µ s. By either metric, the runaway beam is effectively terminated before the onset of particulate melting. Furthermore, the simulations show that secondary radiation emission by tungsten particulates does not significantly impact the runaway termination efficacy of this scheme. Secondary radiation is emitted at lower particle energies than the incident runaways and with a broad angular distribution such that the majority of secondary electrons emitted will not experience efficient runaway re-acceleration. Overall, the stand-off runaway termination scheme is a promising concept as a last line of defense against runaway damage in ITER, SPARC, and other future burning-plasma tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Trabecula: A Random Generator of Computational Phantoms for Bone Marrow Dosimetry

This study was motivated by the efforts on evaluation of radiation risk for leukemia incidence in the Techa River Cohort, where the main bone marrow dose contributors were 89,90Sr (bone-seeking beta emitters). Energy deposition in bone marrow targets was evaluated by simulating radiation particle transport using computational phantoms. The present paper describes the computer program “Trabecula” implementing an algorithm for parametric generation of computational phantoms, which serve as the basis for calculating bone marrow doses. “Trabecula” is a user-friendly tool that automatically converts analytical models into voxelized representations that are directly compatible as input to MCNP.

Shishkina, Elena A.↗
DOE OSTI2022

Applying ALARA Principles in the Design of New Radiological Facilities

The application of ALARA (As Low As Reasonably Achievable) principles to the design of new radiological facilities at Argonne National Laboratory provides a consistent radiation safety basis for future facility operations. The Department of Energy Rule 10 CFR 835 specifies design objectives to be adopted during the design of new facilities for controlling personal radiation exposure. One is to keep exposure levels below 20% of the applicable standards in 10 CFR 835.202. For a radiation worker 20% of the standard corresponds to 5 µSv hr -1 for a 2,000 hour work year. For a member of the public the DOE design objective corresponds to 0.2 mSv in a calendar year. For a new facility worst case radionuclides and their source strengths are chosen. Local shielding is specified to reduce exposure rates to less than 50 µSv hr -1 at 30 cm from the shielding. The current version MCNP6 of the Los Alamos radiation shielding computer program MCNP is then used to calculate the exposure rates elsewhere. Design modifications are made to meet the criteria. The calculations and resulting facility design modifications are discussed for two new radiological facilities.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Neutron-Neutron Correlations in the Photofission of 238U

Background: In the fission of actinides, the nearly back-to-back motion of the fission fragments has a strong effect on the kinematics of fission neutrons. This leads to a favoring of opening angles near 0? and 180? in the neutron-neutron (n-n) opening angle distributions of correlated neutron pairs from the same fission event. Purpose: To measure the n-n opening angle and energy correlations in the photofission of 238U. As of this writing, measurements of correlated n-n opening angle distributions have been reported only for the spontaneous and neutron-induced fission of actinides. This work is the first to report such a measurement using photofission and will provide useful experimental input for photofission models used in codes such as MCNP and FREYA. Method: Fission is induced using bremsstrahlung photons produced via a low duty factor, pulsed, linear electron accelerator. The bremsstrahlung photon beam impinges upon a 238U target that is surrounded by a large neutron scintillation detection system capable of measuring particle position and time of flight, from which n-n opening angle and energy are measured. Neutron-neutron angular correlations are determined by taking the ratio between a correlated neutron distribution and an uncorrelated neutron distribution formed by the pairing of neutrons produced during different beam pulses. This analysis technique greatly diminishes effects due to detector efficiencies, acceptance, and experimental drifts. Results: The angular correlation of neutrons from the photofission of 238U shows a high dependence on neutron energy as well as a dependence on the angle of the emitted neutrons with respect to the incoming photon beam. Angular correlations were also measured using neutrons from the spontaneous fission of 252Cf, showing good agreement with past measurements. Conclusions: The measured angular correlations reflect the underlying back-to-back nature of the fission fragments. An anomalous decline in n-n yield was observed for opening angles near 180? for 238U.

Burggraf, Jeff↗
DOE OSTI2020

Cross section of neutrons from the H 2 ( n , 2 n ) reaction at E n = 15 MeV

In this work, the double-differential cross section of the deuteron breakup reaction 2 H(n, 2n) has been studied experimentally with a neutron beam energy of 15 MeV. Special attention has been devoted to estimation of background condition and multiple scattering effect in the scattering sample. Experimental data have been compared with models based on phase-space approximation used in the ENDF/B-VIII.0 data library and in the MCNP neutron transport code, as well as with rigorous model based on Faddeev equations used for cross section evaluations in JENDL data library. It was found that experimental data are better reproduced by Faddeev model, however, the model overestimates data in the low-energy region of the neutron spectrum (<4 MeV)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Differential cross sections of the 16 O (𝑛,𝛼) reaction at neutron energies from 3.8 to 15 MeV

The importance of studying the 16 O (𝑛,𝛼) reaction is motivated by multiple nuclear applications which rely on an accurate nuclear reaction data library for oxygen. Discrepancies between past experimental data on the 16 O (𝑛,𝛼) 13 C reaction and its time-reverse 13 C ⁡(𝛼,𝑛) 16 O reaction have led to various different nuclear data evaluations. Here, we have measured 16 O (𝑛,𝛼) reaction cross sections using the LENZ instrument with the unmoderated white neutron source at LANSCE. Results from 2016/2017 data are discussed and used to benchmark the mcnp and geant simulations of the LENZ experimental setup. We report partial differential cross sections of 16 O (𝑛,𝛼 0 ) at 𝐸 𝑛 = 3.8–15 MeV and 16 O ⁡(𝑛,𝛼 1 +𝛼 2 +𝛼 3 ) at 𝐸 𝑛 =9–15 MeV, based on the new measurement in 2021. The resonances that we observed are in good agreement with the levels in 17 O that were previously measured. The LENZ cross sections are in better agreement with the ENDF/B-VIII.0 evaluation than the with reduced cross section found in ENDF/B-VII.1 up to 6 MeV. However, the current results for 16 O ⁡(𝑛,𝛼 0 ) 13 C g.s. appear to be in the best agreement over the entire energy range with that of JENDL/AN-2005 (ENDF/B-VI.0).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Development of nuclear data processing code FRENDY version 2

Nuclear data processing is an important interface between an evaluated nuclear data library and neutronics calculation codes. JAEA has been developed the new nuclear data processing code FRENDY from 2013. FRENDY version 1 generates the ACE files used for the continuous-energy Monte Carlo codes including PHITS, Solomon, Serpent,and MCNP and it was released as the open-source software under the 2-clause BSD license in 2019.After we released FRENDY version 1, many functions, e.g., the multi-group neutron cross-section library generation, the statistical uncertainty quantification of the probability tables for unresolved resonance cross-section, the 'perturbation' of the ACE file for uncertainty quantification, and the modification of the ENDF-6 formatted nuclear data file, have been developed. We released FRENDY version 2 including these functions. The present paper gives an overview of FRENDY and presents features of new functions implemented in FRENDY version 2. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Neutron Yield of Thermo Scientific P385 D-T Neutron Generator vs . Current and Voltage

The Thermo Scientific P385 Neutron Generator is a compact neutron source, producing 14 MeV neutrons through the deuterium-tritium (DT) fusion reaction. It is important to measure and understand the dependence of the neutron production rate on the accelerator current and voltage. In this study we evaluated neutron production with an absolutely calibrated liquid scintillator neutron spectrometer (BTI N-Probe), an absolutely calibrated He-3 detector surrounded by HDPE shells (Detec Nested Neutron Spectrometer, NNS), and two uncalibrated ZnS fast neutron scintillators (EJ-410), for both A3082 and A3083 sealed tubes. Here we also modeled the neutron yield using the TRIM code, which calculates the trajectory and the energy loss of deuterons and tritons within the target. Experimental results showed an essentially linear dependence on beam current, as expected. A 3.59 ±0.08 power law dependence on the operating voltage was measured, in effective agreement with the modeled value of 3.5. A series of absolute NNS and N-Probe measurements, matched against MCNP calculations, showed that the A3083 and A3082 tubes provide a maximum neutron yield of 8.2 × 10 8 n/s and 4.7 × 10 8 n/s respectively, with estimated uncertainty of ±10%.We showed, through modeling, that tritium decay is not a significant consideration for tubes, such as these, with lifetimes of less than 10 years.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗