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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
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 OSTI2023

X-ray Inspection Model Validation with Physical Dosimetry

A non-functional printed circuit board assembly was developed using ULTIboard and a simple computational transport model was created for use with the Monte Carlo software MCNP. EBT 3 film radiography was used to compare with Monte Carlo simulation to validate the board representation. The computational transport model was then revised to include connection pins, solder balls, highly attenuating internal structures, and copper trace distribution. Results from the revised model were compared to the EBT 3 films to observe improvements to dose profiles. It was found that this method was useful in verifying the placement of components, as the dose profiles were observed to follow the same trends. The experiment was then repeated using XRQA 2 films to achieve the same level of contrast with 1% the dose of EBT 3 films. It was found that high contrast may be achieved using these films to identify major issues with the model geometry, at a cost of dose profile accuracy. A second validation method was applied to the model using 37 CaF 2 thermoluminescent dosimeters (TLDs). TLD measurements were compared with the simplified and complex transport models to identify the features that have the greatest impact on simulation accuracy. The TLD calibration to CaF 2 was found to be accurate within 5.6%, while calibration to dose in Si was found to be accurate within 4.7%. Finally, it was observed that the accurate representation of solder balls and proper modeling of highly attenuating internal structures had the greatest impact on simulation accuracy.

61 RADIATION PROTECTION AND DOSIMETRY↗
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 OSTI2021

Neutronic Performance of a Compact 10 MWe Nuclear Reactor with Low Enrichment (ThxU1−x)N Fuel

Abstract Recent interest in compact nuclear reactors for applications in space or in remote locations drives innovation in nuclear fuel design, especially non-oxide ceramic nuclear fuels. This work details neutronic modeling designed to support the development of a new nuclear fuel concept based on a mixture of thorium and uranium nitride. A Monte Carlo N-Particle Version 6.2 (MCNP-6) model of a compact 10 MWe reactor design which incorporates (Th x U 1− x )N fuel is presented. In this context, a “compact” reactor is a completely assembled reactor which may be emptied of coolant and transported by specialized commercial vehicle, deployed by a C130J aircraft, or launched into space. Core geometry, reflector barrels, and the heat exchange zones are designed to support reduction of overall reactor volume of core components while maintaining criticality with a fixed total fuel mass of 4500 kg. Dense mixed nitrides of thorium nitride (ThN) additions in uranium nitride (UN) in 5 wt.% increments between $$0.05 \le x \le 0.5$$ 0.05 ≤ x ≤ 0.5 have been considered for calculation of $$k_{\infty }$$ k ∞ and $$k_{{{\text{effective}}}}$$ k effective . ThN additions in UN results in a slight increase in the magnitude of the temperature coefficient of reactivity, which is negative by design. The isotopic distribution of the principal actinide inventory as a function of burnup, time, and initial fuel composition is presented and discussed within the context of the proliferation risk of this core design.

Parker, S. S. (ORCID:000000018597440X)↗
DOE OSTI2020

SCALE capabilities for high temperature gas-cooled reactor analysis

The SCALE code system’s ability to address stochastic distributions of fuel particles within a graphite matrix has been revisited in both multigroup (MG) features and continuous-energy (CE) Monte Carlo methods. Furthermore, this paper describes current and emergent SCALE capabilities within the CSAS sequence to address double-heterogeneous systems and presents verification and validation studies of these methods and data. Good agreement was obtained for a high temperature gas-cooled reactor (HTGR) fuel pebble model between CSAS MG eigenvalue calculations and corresponding CE reference solutions. Code-to-code comparisons for this HTGR pebble model showed good agreement of CSAS-KENO and CSAS-Shift CE calculations and the Serpent and MCNP codes in terms of eigenvalues and reaction rate ratios. Validation studies based on two HTGR experiments resulted in good agreement between MG and CE results, as well as between experiment and calculation, although the level of agreement was significantly influenced by the applied ENDF/B nuclear data library.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Improved MPACT energy deposition and explicit heat generation coupling with CTF

The CASL neutronics code MPACT assumes local energy deposition with equilibrium delayed energy for steady-state and transient calculations. These approximations limit the accurate representation of the heat generation in space and its variations over time, both of which are essential for power distribution and thermal–hydraulic coupling analyses. This paper presents an explicit energy deposition model considering the neutron and gamma heating of all regions and the time-dependent delayed energy. The MPACT/CTF interfaces are updated to transfer the explicit heat sources. The new model agrees well with MCNP for problems without thermal–hydraulic feedback. The MPACT/CTF coupled calculation for a hot-full-power quarter-core case using the new model exhibits a reduction of peak pin power by 2.3% and fuel centerline temperature by 17 K. The new model also shows a keff difference up to -100 pcm in depletion and a 30 K peak fuel temperature reduction in a reactivity-initiated-accident problem.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Measurements of effective actinides neutron capture cross sections in a cadmium-filtered epithermal neutron spectrum

In this work, a total of 20 actinide samples were irradiated in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL): 233 U (1), 235 U (1), 236 U (1), 238 U (1), 237 Np (2), 239 Pu (1), 240 Pu (1), 242 Pu (2), 244 Pu (2), 241 Am (2), 243 Am (2), 244 Cm (2), and 248 Cm (2). The samples were loaded in a basket containing a 1.1-mm thick cadmium filter. The presence of the filter drastically changes the neutron spectrum in the samples and accentuates the reaction rates in the resonance region. Isotopic ratios were determined using Multi-Collector Inductively Coupled Plasma Mass Spectrometry (MC-ICP-MS) that allowed inferring effective capture cross sections. Detailed best-estimate MCNP models of the experiment using ENDF/B-VII.0 nuclear data were developed to calculate the effective neutron cross sections and, subsequently, compare them with the cross sections inferred from experimental results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Serpent transient analyses of GIACINT geometrical change experiments

This paper simulates seven transient experiments with geometrical movements performed at the GIACINT facility of Belarus. These experiments include three slow control rods insertion in more than 25 s, one fast control rod insertion in less than 1 s, and three water moderator draining transients. In the fast control rod insertion experiment, the control rods are first pushed by a mechanical spring and then let fall by gravity in the assembly. In this experiment, the control rods speed varies because of the initial spring force, the gravitational acceleration, and the water buoyancy force. In all other transient experiments, the geometry change occurs with a constant speed. The Serpent computer program was utilized to simulate these transient experiments with geometrical changes by performing two separate computations. One time-independent computation is performed in criticality mode and the other time-dependent computation is performed in dynamic-source mode. The first computation writes the neutron population and delayed neutron precursors for starting the transient simulation. The second computation divides the transient time into several time bins, reads the two files from the first computation, and updates them at the end of each time bin. The Serpent simulation results obtained for this set of experiments are in good agreement with the time dependent experimental measurements. In addition, Serpent and MCNP simulations results show an excellent agreement for the water moderator transient

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

High-level fuel fabrication facility designs from discrete-event simulation

Like other industrial processes, the production of metallic nuclear fuels (MNF) requires that fabrication facilities be able to reliably meet production demands, operate efficiently, and adhere to federal and local safety regulations. In turn, the set of design variables employed by such a facility, such as operations policies, infrastructure and machinery purchased, and the type and number of staff hired, directly impact a facility’s ability to satisfy these goals. Therefore, facility designers must carefully determine which set of design variable values optimally satisfies these constraints. In this paper, we explore how values for these high-level design variables, namely hiring requirements, can be determined in the context of nuclear fuel manufacturing through the coupling of physics-based and discrete-event simulation technologies. Using the Versatile Test Reactor (VTR) program as a case study, we demonstrate how SCALE and MCNP nuclear physics model outputs can be integrated into ExtendSim discrete-event simulation (DES) models of the fuel fabrication process to determine the optimal number of staff hired to ensure fuel production goals are met, operations comply with effective dose limit regulations, and overall project costs are reduced.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Utilization of ACE nuclear data file toolkit ACEtk to calculate relative sensitivity coefficients of point-kinetics parameters

Sensitivity and uncertainty methods are quintessential for nuclear criticality safety and experiment design. This type of analysis relies on calculations of sensitivity coefficients; sensitivity coefficients of the effective neutron multiplication factor with respect to some nuclear data are predominantly calculated and used. As a part of the Laboratory Directed Research & Development project EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data) at Los Alamos National Laboratory, sensitivity coefficients of many radiation detector measurement responses with respect to nuclear data were investigated. Specifically, this paper outlines a method to calculate point-kinetics parameters relative sensitivity coefficients with respect to nuclear data. Point-kinetics parameters such as the prompt neutron decay constant, effective delayed neutron fraction, and neutron generation time are especially important to experimenters and reactor operators designing systems with dynamic neutron populations. This method couples capabilities of the ACE (A Compact ENDF) nuclear data file toolkit, ACEtk, with the ability to load cross sections into the radiation transport code Monte Carlo N-Particle (MCNP). In conclusion, key aspects of optimizing this method for a particular application and sensitivity profiles of the Jezebel criticality experiment are examined and discussed.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
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 OSTI2022

Demonstration of RELAP5-3D for transient analysis of a dual coolant lead lithium fusion blanket concept

Our research aims to build a multiphysics framework for transient analysis of the Dual- Coolant Lead-Lithium (DCLL) blanket design of the proposed Fusion Nuclear Science Facility (FNSF) for all Onticipated Operational Occurrences (AOOs) using Reduced Order Models (ROMs). This paper presents an assessment of using RELAP5-3D for transient thermal-hydraulic blanket analysis of the DCLL design. Preliminary validation studies for the RELAP5-3D properties were conducted using a representative vertical flow loop model and a comparison study using experimental data from the thermal convection corrosion loops at ORNL. These studies inspected the thermal hydraulic response of Lead-Lithium Eutectic (PbLi) systems within the RELAP5-3D code; mainly pressure drop and heat transfer, as compared to analytical calculations and experimental data. The Magnetohydrodynamic (MHD) pressure drop effect for liquid metal under the effects of a uniform magnetic field was implemented into our models using an equivalent forms loss coefficient based on the phasic momentum equations. We verified our implementation of MHD pressure drop within RELAP5-3D using a systematic approach consisting of two studies. Models were developed to perform comparison studies based on experimental data from the Argonne Liquid metal EXperiment (ALEX) facility and the Magnetohydrodynamic PbLi Experiment (MaPLE) facility. Using the developed validation basis, we built a simplified DCLL channel model within RELAP5-3D, inspected MHD pressure drop and performed transient temperature analysis. The DCLL channel model used radial heating data from the International Thermonuclear Experimental Reactor (ITER) which was normalized to the FNSF design. Here, a representative startup transient was performed using a relative power curve from the initial ramp up of a long-term power pulse of the Tore Supra reactor. This study demonstrated that RELAP5-3D ROMs are capable of accurate transient analysis of the DCLL blanket. We have since utilized MCNP to develop heating profiles based on the FNSF design and are working toward implementing them into our DCLL model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Virtual prototyping of liquid metal blanket performance in fusion pilot plant

Liquid metal blanket is a dominant design option for the next step fusion devices responsible for harvesting energy from fusion reaction, and simultaneously producing fuel for the same reaction through tritium breeding. Liquid metal blankets introduce additional complexity to the design due to fluid motion, fluid structure interaction, and magnetohydrodynamic (MHD) effects arising from the motion of the conducting fluid through the magnetic field. They are also directly affected by the plasma heat flux and neutronic fluence. PPPL is currently developing a virtual prototyping system for numerical analysis of the liquid metal blankets for future fusion devices. The system has a customized 3D computational fluid dynamics (CFD) code in its core, allowing MHD flow and conjugate heat transfer analysis in blankets fluids and solids. The code was successfully used before for dual coolant blanket analysis [A. Khodak et al., Fusion Eng. and Des. 137 (2018)]. Recently the same code was modified to allow verified simulation of MHD flows at high Hartmann numbers of several thousand typical for blanket applications. CFD code receives volumetric heat source distribution from the neutronic analysis based on MCNP code. In addition, direct tritium breeding simulation will be performed allowing optimization of the blanket performance. 2D axisymmetric version of neutronics code will be used for rapid optimization, with 3D version employed for detailed analysis. The surface heat distribution on the plasma facing wall will be defined by the software HEAT allowing 3D modeling of the heat flux based on the magnetic field distribution including gyro-orbit effects. Results of thermal analysis are imported into structural analysis code also included in the system. Finally, direct import of CAD geometry will be used for analyzing all components and as a result design option can be efficiently optimized.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

FAST irradiations and initial post irradiation examinations – Part I

The Advanced Fuels Campaign Fission Accelerated Steady-state Test (FAST) at Idaho National Laboratory (INL) completed its first irradiation cycle within the Advanced Test Reactor (ATR). The test focused on the irradiation of alloy fuel forms for use in sodium fast reactors. Tests included the use of fuel additives (Sn, Sb, and Pd) for reducing the gettering of lanthanide fission products, zirconium liners to mitigate fuel-cladding chemical interaction, and annular fuel geometry for sodium-free fuel designs. The first cycle of FAST testing was completed and a series of rodlets were removed for the initial post irradiation examination (PIE) of FAST rodlets. The rodlet irradiation conditions was evaluated using Monte Carlo N-Particle (MCNP) for as-run power history and COMSOL for temperature analysis. These rodlets include a set of low burnups (~2.5 % heavy metal [%FIMA]) control rodlets and a helium bonded annular rodlet (4.7 %FIMA). Initial non-destructive PIE has been completed and includes, visual inspection, neutron radiography and gamma scanning of the FAST rodlets. Radiography confirmed the integrity of the experiments as well as showing the complete filling of the annulus in the annular fuel at a modest burnup (4.7 %FIMA) and potential slumping of the cooler rodlets at lower burnup. Precision gamma scanning indicated mostly usual fission product behavior with the exception of the ceasium in the He-bonded annular fuel. Future destructive PIE will be necessary in order to fully understand the effects of accelerated irradiation on U-Zr metallic fuel behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

New evaluation on angular distributions and energy spectra for neutron-induced charged-particle measurements

Recent development on neutron-induced charged particle reaction measurements can provide high-quality double differential cross sections at Los Alamos Neutron Science Center. With newly measured cross section data, we found lacking nuclear outputs for neutron-induced charged particle reactions in the latest ENDF/B-VIII.0 library, when compared with Monte Carlo simulations. Here, this paper presents how we improved the missing angular distributions and energy spectra for these reactions and the validation of our new evaluations with the measurements performed with the LENZ instrument at LANSCE. The impact of this additional evaluation is demonstrated by comparing the MCNP simulations with the LENZ experimental data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Confirming the absence of nuclear warheads via passive gamma-ray measurements

Arms-control agreements between the United States and Russia negotiated after the end of the Cold War have imposed limits on the number of deployed strategic nuclear weapons. Verification of these agreements has relied on onsite inspections, sometimes supported by radiation detection techniques to confirm the absence of a nuclear warhead when ambiguities arise. So far, these measurements have sought to detect neutron emissions associated with the presence of plutonium, but they would be inadequate for uranium devices. In an effort to offer instruments that could be used to confirm the absence of both plutonium and uranium weapons, here we propose an inspection system that uses only passive gamma radiation detection techniques. Such a system would be particularly valuable for next-generation arms-control agreements that limit total numbers of weapons and would involve containerized items in storage. We conducted extensive Monte Carlo simulations to support the development of a verification protocol and detection algorithm. Here, we demonstrate the viability of the technique using standard laboratory check sources and MCNP simulations for simplified configurations of special nuclear material.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗