Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ORIGEN”

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

New capabilities for rapid depletion analysis of pebble-bed reactors in SCALE

The SCALE Leap-In for Cores at Equilibrium (SLICE) method leverages capabilities available in the SCALE modeling and simulation suite to facilitate rapid estimation of equilibrium core inventories for flowing-pebble reactor systems in support of a variety of fuel cycle applications. New capabilities that compliment the SLICE method have been developed for the ORIGAMI interface to ORIGEN in SCALE to facilitate rapid depletion calculations for flowing-pebble systems. New features include a more generalized and flexible means of specifying interpolation dimensions, the ability to evaluate pebbles moving through user-defined “axial zones” in each of the pebble’s multiple passes through the core, and the treatment of differential velocities across radial channels (representing both pebble-to-pebble and wall-to-pebble friction effects). ORIGAMI thus provides an efficient user interface to define a pebble’s path and irradiation histories of its multiple passes through the core to calculate the pebble’s time-dependent inventories, which can be useful in various applications domains such as safeguards, criticality safety, and disposal analyses. The consistency of new ORIGAMI approach was verified against the solution obtained using the SLICE method with ORIGEN-ARP, showing excellent agreement.

Skutnik, Steven E. [Oak Ridge National Laboratory ↗

The Relationship Between Dose Rate and Decay Heat for Spent Nuclear Fuel Casks

Decay heat and dose rate are two important limits used for determining the allowable contents of spent fuel (SNF) in dry storage systems and transportation packages. While the decay heat limit is used to maintain fuel cladding integrity and ensure retrievability, dose rates are used to demonstrate compliance with regulatory requirements on radiation protection. Because both dose rate and decay heat result from decay of radioisotopes in SNF, this study is an attempt to examine the relationship between dose rate and decay heat for a given cask design. Dose rates were evaluated for 198 cask configurations, that include various SNF system designs (e.g., storage, transfer, transport), SNF characteristics (e.g., fuel types, burnup, cooling time), and loading maps (e.g., uniform loading, zone loading), while a constant decay heat was maintained. The decay heat was calculated using US Nuclear Regulatory Commission (NRC) Regulatory Guide (RG) 3.54, Revision 2, and verified using ORIGEN sequence within SCALE code system. The ORIGEN outputs were used as source terms in the dose analysis using 198 different configurations. A computer script was developed to calculate the cooling time necessary to achieve a given decay heat for a given enrichment, assembly average burnup, assembly mass, and in-core history using a rootfinder algorithm. Initially the computer script was developed to provide cooling time and burnup calculations directly to the analysis of dose rates and decay heats, so a comparison between Used Nuclear Fuel-Storage Transportation and Disposal Analysis Resource Data System (UNF-ST&DARDS) results and RG3.54 data was made; results are included in the appendix to this document. However, an iterative approach was used to compute cooling time, and its accuracy did not depend on the results of the RG3.54r2 algorithm, although the algorithm was still used. The results of the evaluation presented herein clearly demonstrate that a given decay heat does not correspond to a unique dose rate for a variety of cask and package designs. There is no clear pattern to develop a correlation between decay heat and the source terms. Depending on burnup, enrichment, cask type, and loading pattern, dose rates varied for the exact same decay heat—in some cases by 400% for a given cask. For cases in which decay heat was held constant through selection of the appropriate cooling time, dose rates would decrease with increasing burnup, and in other cases, dose rates would increase. The large variation in dose rates for a constant decay heat indicates that casks loaded based on decay heat—that is allowing any burnup, cooling time, and enrichment combinations that yield the qualified decay heat limit(s) —cannot ensure that an Independent Spent Fuel Storage Installation or a spent fuel transportation package will meet the regulatory limits set forth by the respective regulations, i.e., 10 CFR 72 or 10 CFR 71.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

HFIR Activity Workbook Generator (HAWK) User Guide

The HFIR Activity WorkbooK generator (HAWK) is a Python code that automates and streamlines the activity calculation of samples after irradiation in the High Flux Isotope Reactor (HFIR). HAWK’s results provide estimates of the activity and nuclide inventory of irradiated specimens before they are moved to hot cell facilities, where they undergo post-irradiation examination. The samples’ activity results guide the packing of shipping containers and inform the accountable inventories for the hot cell facilities. The toolkit was originally developed by Charles Daily, a former R&D staff member at Oak Ridge National Laboratory (ORNL). As of May 2025, HAWK is developed by the Radiation Transport & HPC Methods Group (Nuclear Energy and Fuel Cycle Division) at ORNL. Figure 1 presents HAWK’s workflow. To use HAWK, users need to: 1. Develop an Excel input workbook (i.e., XLSX extension) containing data from the experiment’s materials, irradiation history (cycles), and irradiation positions. 2. Make minor edits to an existing template JSON file (i.e., auxiliary_data.JSON) and to the Python driver. The driver sets the necessary environment variables, defines the material compositions, and ultimately calls HAWK. Once configured, HAWK runs the Oak Ridge Isotope Generation code (ORIGEN) to calculate the masses, activities, and heat load at the end of irradiation for each isotope in the specimen. ORIGEN is part of SCALE, ORNL’s in-house computational tool for performing nuclear safety and design calculations. Following this step, HAWK postprocesses the results and generates three output workbooks summarizing the activity calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Modeling and Simulation of Fuel Burnup in Pebble Bed Reactors

Modeling and simulation of fuel burnup plays important roles in reactor design, operation, safety, and security as well as nuclear material control and accounting (MC&A) [1]. This task is uniquely challenging for pebble bed reactors (PBR) because the pebbles are continuously added and recycled into the reactor core, and their paths through the core are random. To address this problem, we present two simulation models in this paper. Brookhaven National Laboratory (BNL) developed a simple lattice model of a PBR in Serpent software to generate used pebble isotopic concentrations. The benefit of using Serpent software in this specific application is that it helps streamline the data generation process without having to use too many independent software codes in combination to achieve a simple task. For example, transport, burnup and zero power decay can be implemented in a single pass. Three-dimensional core models were developed using Serpent to simulate the burnup process of 5 subject pebbles starting from fresh till they reach nearly target burnup, with each pebble placed in one of the five artificially designated radial channels to capture the changing neutron spectra along the core radius. Equilibrium isotopic concentrations were assumed in all other pebbles in the core. To provide a verification for the Serpent isotope transmutation and decay results, Oak Ridge National Laboratory (ORNL) performed simple SCALE/ORIGEN calculations using the average neutron spectra calculated by Serpent for each of the 5 pebbles. The 252-group neutron spectra from Serpent were then used by ORIGEN to produce the one-group library for depletion and decay calculations. The isotopic concentrations of a few nuclides of interest and neutron and gamma source terms produced from the ORIGEN calculations were compared with the ones from the Serpent calculations. The model simulated in this work was based on the Pebble Bed Modular Reactor (PBMR)-400 design because many data needed for the simulation such as core power profiles, fuel and reflector temperatures, and equilibrium core composition are publicly available. In this paper, we will compare the results between these two approaches and benchmark the results against a set of well-established simulation results for PBMR-400.

Dim, Odera↗

Modeling and Simulation of Fuel Burnup in Pebble Bed Reactors

Modeling and simulation of fuel burnup plays important roles in reactor design, operation, safety, and security as well as nuclear material control and accounting (MC&A) [1]. This task is uniquely challenging for pebble bed reactors (PBR) because the pebbles are continuously added and recycled into the reactor core, and their paths through the core are random. To address this problem, we present two simulation models in this paper. Brookhaven National Laboratory (BNL) developed a simple lattice model of a PBR in Serpent software to generate used pebble isotopic concentrations. The benefit of using Serpent software in this specific application is that it helps streamline the data generation process without having to use too many independent software codes in combination to achieve a simple task. For example, transport, burnup and zero power decay can be implemented in a single pass. Three-dimensional core models were developed using Serpent to simulate the burnup process of five subject pebbles starting from fresh till they reach nearly target burnup, with each pebble placed in one of the five artificially designated radial channels to capture the changing neutron spectra along the core radius. Equilibrium isotopic concentrations were assumed in all other pebbles in the core. To provide a verification for the Serpent isotope transmutation and decay results, Oak Ridge National Laboratory (ORNL) performed simple SCALE/ORIGEN calculations using the average neutron spectra calculated by Serpent for each of the five pebbles. The 252-group neutron spectra from Serpent were then used by ORIGEN to produce the one-group library for depletion and decay calculations. The isotopic concentrations of a few nuclides of interest and neutron and gamma source terms produced from the ORIGEN calculations were compared with the ones from the Serpent calculations. The model simulated in this work was based on the Pebble Bed Modular Reactor (PBMR)-400 design because many data needed for the simulation such as core power profiles, fuel and reflector temperatures, and equilibrium core composition are publicly available. In this paper, we will compare the results between these two approaches and benchmark the results against a set of well-established simulation results for PBMR-400.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Assessment of alternative radionuclides for use in a radioisotope thermoelectric generator

238 Pu is the current fuel source for MMRTGs. This is due to favorable properties, such as high decay heat and low gamma emissions. 238 Pu supplies are dwindling rapidly, necessitating the review of potential replacement isotopes for deep space missions. This work considers nine potential replacements, Ac-227, Am-241, Cf-252, Cm-244, Cs-137, Po-210, Sr-90, Th-228, and Th-170. We have assessed performance of alternatives to 238 Pu by modeling the specific decay-heat and gamma spectra of several potential radioisotopes in SCALE-ORIGEN, and comparing these with 238 Pu. An analysis of heat output, heat stability, safety, availability and required mass, baselined against 238 Pu demonstrates that of those radioisotopes studied, 241 Am has properties that make it most conducive for 238 Pu replacement.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Comparison of Irradiated TRISO Fuel Radioactivity from Multiple Advanced Reactor Designs

The SCALE code system was used to model, deplete, and compare several different TRISO-fueled reactor designs: a helium-cooled prismatic reactor, a helium-cooled pebble-bed reactor (PBR), and a Fluoride-Lithium-Beryllium (FLIBE) molten-salt-cooled PBR. The purpose of this comparison was to understand how differences in the reactor designs affect the radioactivity of the fuel after discharge and whether those differences are significant. First, the various reactor designs were build and depleted in the TRITON module for each design and fuel enrichment. Then, the TRITON outputs were used to create burn-up dependent reactor libraries. These libraries were then used by ORIGEN to determine the activities of discharged fuel for each reactor, which were compared to generic Westinghouse 17x17 fuel. Overall, the results showed that short term decays are dominated by reactors with higher operating powers, and the reactor type, initial fuel enrichment, and maximum burn-up are of only secondary importance. Although this analysis only focuses on activities in becquerels, these dependencies are consistent with the expected behavior of decay heat. However, analysis of long- term time periods post-irradiation shows that the reactor type and maximum burn-up have strong impacts on the activities; initial fuel enrichment has a secondary impact while operating power is inconsequential. These results would be useful for analyses, such as dose assessment and modeling in post- release scenarios; normal fuel handling operations; and spent fuel transport, storage and disposal. Of particular interest, the results in this report show that analyses that focus on spent nuclear fuel of advanced reactors need to consider each parameter carefully. Unsurprisingly, if the correct operating power is not used in short term analyses, the results will not be correct. Perhaps unexpectedly, however, if the correct reactor type is not used, then the long term results will also be incorrect, especially for areas such as permanent disposal. Even though this report focuses on the total activity of nuclear fuel, it provides initial results on the effects of various input parameters and also provides a framework to extend the work into other analyses of spent fuel from advanced reactors, especially those employing TRISO fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Advanced User Interface Capabilities [Slides]

This tutorial will review the data plotting and geometry visualization capabilities in the Fulcrum user interface. This tutorial will help you become familiar with Fulcrum’s 2D plot, and 2D and 3D geometry visualization features. You will learn how to identify plottable data items, compose and export plot and plot data for SCALE plot formats (SDF, Ampx MG/CE, PLT, F71, PTP, SPF, ORIGEN Gamma data, etc.) and visualize, navigate, cut, hide, and export the geometry and spatial data (fission-, dose-map, etc.) overlays in 2D and 3D. No prior experience with SCALE is required. Attendees can follow along using 6.3.0-beta.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

AMPX and SCALE Nuclear Data Libraries for Depletion [Slides]

This presentation touches on AMPX and SCALE nuclear data libraries for depletion. AMPX is a nuclear data processing code suite that generates continuous-energy and multigroup neutron and gamma libraries used in SCALE and VERA, covariance data used in S/U calculations, and depletion, activation and decay origin data for Oak Ridge Isotope Generation code (ORIGEN).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fission Product Yields: SCALE/ORNL Perspective [Slides]

This presentation covers a few key points. First and foremost, that ORIGEN is a consumer of FPY data. Secondly, TRITON and POLARIS are assets which can be used to validate new FPY data, e.g., from existing radiochemical benchmarks and new irradiation experiments. Lastly, this presentation finds that new capabilities with covariance data should be developed, and V and V methods and analysis are also needed for the covariance data.

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↗

Radioxenon signatures of molten salt reactors

Developments of molten salt reactor (MSR) technologies are making rapid progress across the globe. These reactor designs involve properly controlling radionuclides through off-gas systems. This work examines radioxenon emissions and assesses treaty detection technology to be used for monitoring. This work also includes a sensitivity study with ORIGEN on the nuclear forensic signals possible from a MSR. Multiple Isotope Ratio Comparison (MIRC) plots are used to compare the ratios of radioxenon isotopes in an MSR to a highly enriched uranium (HEU) reactor. Finally, the results show that MSR emissions may significantly overlap with signals produced in a HEU reactor and that changing reactor operations or retaining emissions in an off-gas system can significantly shift the radioxenon signature away from the HEU pulse.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Californium-252 production at the High Flux Isotope Reactor - I: Validation study using campaign data

This paper presents a series of 252 Cf production validation and code-to-code comparison studies performed based on data from the production campaigns at the High Flux Isotope Reactor (HFIR). These studies support efforts to convert HFIR from using highly enriched uranium (HEU) fuel to low-enriched uranium (LEU) fuel. HFIR must maintain its world-class performance and missions following this conversion, and because 252 Cf is a vital neutron-emitting radioisotope used for a variety of high-impact applications (e.g., reactor startup, cancer treatment), the ability to efficiently produce 252 Cf must be preserved. In this work, the HFIRCON, Shift, ORIGEN, and TCOMP codes were deployed, and several sets of data libraries were investigated to better understand the calculation codes and the data biases. As-loaded target composition data, as-run irradiation history data, and post-irradiation measurements from recent multi-cycle irradiation campaigns of the HEU core were used to validate and determine methodology biases. Further, the findings demonstrated a good agreement, with results falling within 3 standard deviations of measurements. This paper lays the ground work for the second paper, which evaluates and compares 252 Cf production and safety metrics with the HEU core and a proposed LEU core.

07 ISOTOPE AND RADIATION SOURCES↗

Modeling and simulation of oxygen transport in high burnup LWR fuel

We have developed a formulation for oxygen transport in uranium dioxide nuclear fuel that accounts for the effects of irradiation. The overall simulation combines the evolving isotopic composition, thermochemistry, and oxygen transport in irradiated fuel. The driving forces for oxygen transport are computed from local thermodynamic equilibrium calculations and account for the effects of temperature gradients and composition, including fission products. The proposed method provides a mechanism for including complex thermodynamic models of nuclear fuel in modeling of mass redistribution, and alleviates difficulties associated with the common thermodiffusion formulation. Here, the transport model has been implemented within the nuclear fuel performance code BISON utilizing the thermochemistry code Thermochimica, with burnup calculations provided by the ORIGEN isotopic transmutation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fission product distribution in irradiated safety-tested and as-irradiated AGR-2 TRISO particles

In this work, two tristructural isotropic (TRISO)-coated nuclear fuel particles were examined by electron probe microanalysis (EPMA) as part of the Advanced Gas Reactor program. The compacts’ average irradiation temperatures ranged from approximately 1260 to 1290 °C. One particle was examined in the as-irradiated condition, while the other was subject to 1600 °C post-irradiation safety testing. This study was undertaken to test a newly-developed EPMA technique to determine fission product masses in TRISO particles on a layer-by-layer basis, and to compare fission product distributions between an as-irradiated and safety-tested particle. Fission product concentration profiles were collected along two radii in each particle, with measured concentrations used to compute the fission product mass in each TRISO particle layer. These measured masses were then compared to those predicted from ORIGEN modeling calculations. Data collected from these measurements show that for these two particles, masses determined via EPMA were within ± 20% of the calculated masses for the rare-earth elements, Mo, Zr, Cs, I, and Pd. Elements that tend to be less homogeneously distributed include Sr, Te, Eu, Ag, and possibly Ba. Measured Ag masses differed by more than 40% from the calculated mass. Lanthanides other than Eu remain primarily within the fuel kernel in the as-irradiated particle but in the safety-tested particle these element masses were divided approximately equally between the kernel and kernel periphery. In both particles, the majority of Sr and Eu accumulated in the carbon-rich kernel periphery, although in the safety-tested particle, Sr and Eu accumulated farther from the fuel kernel than occurred with irradiation alone. A greater mass fraction of mobile elements, such as Cs and I accumulated in the buffer and IPyC in the safety-tested particle as compared to the as-irradiated particle. When fully developed and tested, this mass balance approach to TRISO particle analysis has the potential to provide insight into fuel behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling of Am-241 as an alternative fuel source in a radioisotope thermoelectric generator

Radioisotope thermoelectric generators (RTGs) have been utilized in the USA to power satellites and space exploration equipment for over half a century. In our preceding paper, Dustin and Borrelli (2021), we assessed the performance of 241 Am as an alternative to 238 Pu by modeling the specific decay-heat and gamma spectra of several potential radioisotopes in SCALE-ORIGEN, and comparing these with 238 Pu An analysis of heat output, heat stability, safety, availability and required mass showed that of the isotopes reviewed, 241 Am had the greatest potential for replacing 238 Pu, with additional design considerations. In our current work, we have assessed radiation transport for two configurations of an 241 Am powered Multi-Mission Radioisotope Thermoelectric Generator (MMRTG), due to the potential health effects for operators, fabricators, and others with the potential to be exposed to neutrons and lower energy gamma rays. We have modeled the use of americium as a fuel source replacement for plutonium in current MMRTG designs. The assessment has been performed using MCNP6, in terms of neutron and gamma emissions that might interfere with sensitive equipment and power system production. The radiation transport behavior in terms of dose rate of the 241 Am fueled models is compared with that of a 238 Pu fueled model. Finally, additional work was performed to validate modeling simplifications with regards to thermocouple material configurations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Secondary-Source Core Reload Modeling with VERA

The CASL reactor simulation package VERA has been adapted to provide high-fidelity simulation capabilities for modeling source range detector response during subcritical reactor configurations. New features include the activation and shuffling of secondary-source assemblies, use of burned fuel neutron emission data from the ORIGEN depletion solver to the MPACT deterministic neutron transport solver, allowance of user-defined sources in MPACT based on material composition, ability to solve the subcritical source-driven system with neutron multiplication using the MPACT diffusion solver, and transfer of the calculated fission source from MPACT to the continuous-energy Monte Carlo solver Shift for final detector response evaluation using the CADIS methodology for variance reduction. These new capabilities were benchmarked against Watts Bar Unit 1 plant operating data for the first few fuel loading steps and were found to demonstrate excellent agreement with the measured data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Cross-Code Verification of Neutronics Analysis Tools at INL Applied for 238 Pu Production in the Advanced Test Reactor

Here, analyses are completed for experiments prior to experiment irradiation in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). Various codes are used to qualify all experiments planned for insertion in the reactor, thereby ensuring that all safety and programmatic requirements are satisfied preirradiation. Among the common experiment analysis tools at INL are MCNP5 coupled to ORIGEN2 (MOPY) and MC21. MOPY uses MCNP5 for transport calculations along with calculations for fluxes and select reaction rates, and then ORIGEN2 handles the step-by-step and postirradiation depletion. MC21 handles all in-reactor transport and step-by-step, during-irradiation, depletion calculations, and then ORIGEN (SCALE 6.2.3) is used for decay and dose calculations postirradiation. The MOPY results, along with those obtained via two variations of the MC21 model, were compared in terms of 238 Pu production in the ATR’s H10 position. For the MOPY model, the MC21 model utilizing the HELIOS-based fission product (FP) library, and the MC21 model utilizing the expanded 1300 FP library, the during-cycle irradiation in-core heating results were sufficiently equivalent; however, the MOPY model and the MC21 model with the HELIOS library showed some differences relating to the respective FP libraries. Ultimately, the MC21 model with a 1300 FP library produced the most consistent results throughout the cycle, whereas the MC21 model that utilized the (smaller) HELIOS library was able to handle during-irradiation analysis but lacked certain short-lived FPs that significantly contributed to the total decay heat at shutdown. MOPY, on the other hand, was found to overpredict fission gas production, as a result of limitations in the ORIGEN2 code.

ATR↗