Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fission Calculations”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗
DOE OSTI2024

Three-Dimensional Full-Core BEAVRS Using OpenMOC with Transport Equivalence

Using an optimized implementation of the three-dimensional (3D) method of characteristics for neutron transport, along with a novel equivalence method for transport calculations that was designed to correct self-shielding errors from neglecting the angular dependence of resonant group absorption, a 3D full-core light water reactor hybrid stochastic-deterministic eigenvalue calculation was achieved. This paper presents the optimizations developed and compares the transport solutions obtained. For the statepoint, run times near 10 000 CPU hours are achieved—improving on previous works by an order of magnitude—with near 1% error on pin fission to 238 U capture ratios and a few dozen pcms on the eigenvalue.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

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

Validation of KENO Delayed Neutron Fraction Capabilities

The delayed neutron fraction, which is necessary for the safe operation of reactors, varies by fissionable material and neutron spectrum. A prompt fission neutron spectrum has been included in the SCALE 6.2 release to account for the removal of delayed neutrons from precursor materials from the neutron balance. This paper describes how experimental measurements from criticality safety benchmarks were used to validate the SCALE calculational procedure for the delayed neutron fraction. Additionally, the effect of nuclear data uncertainty on the KENO calculation of β eff was examined with the SAMPLER sequence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Microstructure, Thermal, and Mechanical Properties Relationships in U and UZr Alloys (Final Report)

Uranium-zirconium (U-Zr) alloys are candidate fuel systems for transmutation based reactors that can be used to burn long-lived minor actinides and fission products in fast spectrum reactors. Metallic fuels have also been gaining more recent attention for applications as accident tolerant fuels, as well as for use in small modular reactors. This research focused on a “science-based” approach to capture the connections between U and U-Zr alloys’ three-dimensional (3-D) microstructure, thermal properties, and mechanical properties through closely coordinated experiments and modeling efforts from the unirradiated to the irradiated fuels. Advanced characterization and modeling techniques were used to understand irradiation-induced microstructural evolution and its direct impact on the thermal and mechanical properties of U and U-Zr fuel. Closely coordinated experiments and modeling were performed to provide crucial data that does not currently exist. Overall, this research spanned multiple length and time scales within the models and experiments. The scope of the research encompassed the understanding of the irradiation effects in U and various U-Zr alloys with particular attention paid to three task areas: (1) microstructural evolution, (2) in-situ/ex-situ thermal and mechanical properties, and (3) multiscale modeling connections to microstructure, thermal, and mechanical properties. This research resulted in (1) the 3-D characterization of neutron irradiated U-Zr fuel in multiple phase regions to better understand fission gas swelling and constituent redistribution, (2) development of a microstructural model linking thermal and mechanical properties via in situ Raman and nanoindentation, (3) and mesoscale phase field modeling was coupled with the AEH method in the MOOSE framework was used to calculate the effective thermal conductivities of U-Zr fuels consisting of α-U and δ-UZr 2 heterogeneous microstructures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1971

Uranium-bearing minerals of lunar rock 12013.

The U distribution in rock 12013 was studied by fission track and elemental mapping techniques. Major U-bearing phases are whitlockite, apatite, zircon, and phase beta, which is a Zr-Ti mineral rich in Fe, Nb, Y, REE, and containing up to 3.6% UO2, 4.7% ThO2 and 4.2% PbO. Calculated microprobe ages for phase beta average 4.0 b.y. and are in reasonable agreement with Rb-Sr and K-Ar ages. Phase beta plays a significant role in the U-Th-Pb systematics of rock 12012 and may play a similar role in the model ages of lunar soil.

Haines, E. L.↗
NASA NTRS1976

Plasma core reactor applications

Analytical and experimental investigations were conducted to demonstrate the feasibility of fissioning uranium plasma core reactors and to characterize space and terrestrial applications for such reactors. Uranium hexafluoride fuel is injected into core cavities and confined away from the surface by argon buffer gas injected tangentially from the peripheral walls. Radiant heat transfer calculations were performed for a six-cavity reactor configuration. Axial working fluid channels are located along a fraction of each cavity peripheral wall. Results of calculations for outward-directed radiant energy fluxes corresponding to radiating temperatures of 2000 to 5000 K indicate total operating pressures from 80 to 650 atm, centerline temperatures from 6900 to 30,000 K, and total radiated powers from 25 to 2500 MW, respectively. Applications are described for this type of reactor such as (1) high-thrust, high specific impulse space propulsion, (2) highly efficient systems for generation of electricity, and (3) hydrogen or synthetic fuel production systems using the intense radiant energy fluxes.

Latham, T. S.↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
DOE OSTI2020

Self-assembly of solid nanoclusters in molybdenum under gas ion implantation

The self-assembly of defects in materials could create unique physical and chemical properties. In nuclear fuels, the formation of self-assembled fission products mitigates fission-gas-release and swelling. In this work, we report the formation of a solid nanocluster superlattice in molybdenum, an important constituent in metallic fuels, under krypton gas ion implantation. X-ray absorption near edge structure measurements coupled with density functional theory calculations show that the pressure of krypton nanoclusters in molybdenum is ~2 GPa and molecular dynamic simulations reveal that the krypton atoms are solid in the nanoclusters with a krypton-to-vacancy ratio of ~0.28 at 300 K.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Yrast structures of 98 Nb and 99 Mo

Neutron-rich nuclei in the A ≈ 100 mass region are interesting due to a rapid shape transition, especially pronounced in the Zr isotopes, and more recently observed in Nb isotopes. 98 Nb, with only one proton and one neutron outside the subshell closure nucleus of 96 Zr, is amenable to shell model calculations. To further examine the rapid shape transition the yrast structure of 98 Nb was established in this work. This was the only yrast structure missing from all immediate neighbors to 96 Zr. The yrast structure of 98 Nb was studied in the fission of the compound systems formed in three heavy-ion induced reactions, 24 Mg (134.5 MeV) + 173 Yb, 23 Na (129 MeV) + 176 Yb, and 18 O (91 MeV) + 208 Pb. Prompt γ-ray spectroscopy was accomplished using the Gammasphere array. Excitation energies up to 3 MeV were observed for the first time in 98 Nb. The yrast structure above the previously known (5) + isomer was established. In the process of studying 98 Nb the yrast structure of positive-parity states in 99 Mo was extended to 3.7 MeV excitation energy, the previously-known 99 Nb level scheme was enriched, and two new levels were added in the level scheme of 97 Zr. The coupling of the odd proton occupying the g 9/2 orbital to the yrast states in the core nucleus of 97 Zr can account for all observed states in 98 Nb. The yrast structure for the positive-parity states of 99 Mo is compared to the deformed ground-state bands of the 101 Ru isotone and of 98,100 Mo.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1972

Post impact behavior of mobile reactor core containment systems

The reactor core containment vessel temperatures after impact, and the design variables that affect the post impact survival of the system are analyzed. The heat transfer analysis includes conduction, radiation, and convection in addition to the core material heats of fusion and vaporization under partially burial conditions. Also, included is the fact that fission products vaporize and transport radially outward and condense outward and condense on cooler surfaces, resulting in a moving heat source. A computer program entitled Executive Subroutines for Afterheat Temperature Analysis (ESATA) was written to consider this complex heat transfer analysis. Seven cases were calculated of a reactor power system capable of delivering up to 300 MW of thermal power to a nuclear airplane.

Puthoff, R. L.↗
DOE OSTI2022

New capabilities of the MORET 6 Monte Carlo neutron transport code

The MORET code is a simulation tool that solves the transport equation for neutrons using the Monte Carlo method. It allows users to model complex three-dimensional geometrical configurations in a user-friendly way. New features have been introduced to extend the application field of MORET beyond the usual criticality calculations for which it has been initially designed. The most important change is the addition of an analog fixed source mode which allows studies of systems of any reactivity combined with very flexible outputs. Other useful improvements have been added concerning the geometric part, the fission matrix, the multigroup sensitivity coefficients and the outputs. This paper presents an overview of these new features. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

TRIPOLI-4 neutronics calculations for IAEA-CRP benchmark of CEFR start-up tests using new libraries JEFF-3.3 and ENDF/B-VIII

Sodium-cooled Fast Reactors (SFR) are one of the promising Generation IV fast reactors. The China Experimental Fast Reactor (CEFR) is a 65 MWth pool-type SFR with a high neutron leakage core using high enriched uranium oxide fuel. The CEFR start-up tests in 2010 consist of series of neutronics experiments. Essential experimental data are available from the 2018 IAEA-CEFR-CRP benchmark document and useful for the validation of neutron transport codes and nuclear data libraries. The TRIPOLI-4 Monte Carlo transport code is a general-purpose neutronics code using continuous-energy nuclear data libraries. It has a rich validation database covering different computational and experimental benchmark data sets to assure the accuracy and credibility of numerical studies. Previous TRIPOLI-4 SFR core physics calculations used mainly MOX fuels. SFR control rod worth studies with TRIPOLI-4 were mainly on big-size cores. To test the modeling capability of TRIPOLI-4 on fuel loading patterns and on variable control rods positions for SFRs and to check recent nuclear data libraries, it is interesting to investigate the CEFR start-up tests in this work, including core states from subcritical to supercritical, control rods and rod group worth, point kinetics parameters, and radial foil activation measurements. TRIPOLI-4 calculation results using new data libraries JEFF-3.3, ENDF/B-VIII, and those of ENDF/B-VII.1 for different CEFR core states are reported here. By means of different options of TRIPOLI-4, control-rod SAs reactivity worth, point kinetics parameters, and radial {sup 237}Np(n,f) fission rate distribution were successfully evaluated. Both ENDF/B-VIII and JEFF-3.3 nuclear data libraries provided reliable results with a difference of 170 +/- 11 pcm in k{sub eff}.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Multi-scale modeling of wastage layer formation in metallic fuel cladding

Fuel-cladding chemical interaction (FCCI) is a major concern for U-Zr metallic fuels' performance, primarily due to the formation of a brittle layer (wastage) in the cladding. This brittle layer, resulting from intermetallic compounds between cladding constituents Fe, Cr, and lanthanide fission products, significantly impacts the cladding's mechanical integrity. Recent efforts focus on developing a mechanistic modeling framework to understand lanthanide production, transport to the fuel-cladding interface, and phase transformation to intermetallic phases. A multi-scale computational approach has been used to calculate lanthanide transport rates, with atomistic calculations determining Nd diffusivities through the solid fuel matrix and along pore surfaces. These diffusivities inform a mesoscale model to determine an effective diffusion coefficient, accounting for porosity and infiltration with bond sodium. This effective diffusivity is used in engineering-scale simulations via the BISON fuel performance code, which has been validated against EBR-II and FFTF reactor experiments.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

On the Statistical Uncertainty of Monte Carlo-Calculated Scattering Sensitivities

Sensitivity coefficients calculated with Monte Carlo codes are widely used for nuclear data uncertainty quantification in the modeling and simulation of complex 3D reactor systems. This study systematically compares sensitivity coefficients and associated statistical uncertainties for the multiplication factor and fuel temperature reactivity across multiple Monte Carlo codes (SCALE/KENO, SCALE/Shift, MCNP, and Serpent) using simple models representing light-water reactors and advanced reactor concepts. For multiplication factor sensitivities, statistical uncertainties are generally acceptable, although scattering sensitivities show significantly larger statistical uncertainties than, for example, fission and capture reactions. Fuel temperature reactivity sensitivities show significantly larger statistical uncertainties across all reactions. Elastic scattering sensitivities are the most problematic: all Monte Carlo codes fail to resolve energy-dependent coefficients, and they produce dramatically different energy-collapsed values. Critically, the use of these sensitivity coefficients in nuclear data uncertainty propagation leads to reduced statistical uncertainties in individual uncertainty contributions. This can lead to the masking of unusable sensitivity coefficients and producing misleading uncertainty results. The findings of this study show that new or enhanced methods are needed to improve Monte Carlo elastic scattering sensitivity calculations. Additionally, this study shows the relevance of verifying sensitivity coefficients through direct perturbation calculations for individual nuclide reactions, instead of only for total cross sections as commonly done.

Bostelmann, Rike [ORNL] (ORCID:0000000165968088)↗
DOE OSTI2020

Calculations of Astrophysical Reaction Rates Using ENDF/B-VIII.0 Library

Recent observations of neutron stars merger (GW170817) renewed interest in stellar nucleosynthesis. Stellar nucleosynthesis modeling requires fully traceable, unbiased, high fidelity nuclear data. The Evaluated Nuclear Data File (ENDF) libraries contain complete collections of reaction data sets over the nuclear industry standard 10 -5 eV - 20 MeV energy span. For the first time, the (n,γ), (n,fission), (n,p), and (n,α) astrophysical reaction rates were computed within the 0.01-10 GK range of temperatures using the ENDF/B-VIII.0 library, and nuclear Reaction rate Library for astrophysics (REACLIB) fit parameters were deduced. The present results were used to estimate the slow neutron capture timescale. These findings demonstrate the potential astrophysical applications of the ENDF libraries and the complementary relations between the nuclear industry and astrophysics data developments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Implementation of a Drift-Flux Model in SAM for Modeling of Passively Transported GAS in Molten Salt Reactors

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is further developing the Multiphysics Object Oriented Simulation Environment (MOOSE) system thermal hydraulics (T/H) code, System Analysis Module (SAM), to include modeling of molten salt reactor (MSR) designs. MSR designs in which the fuel is dissolved in the salt coolant poses unique modeling and simulation (M&S) challenges because fission products—some of which will be noncondensable gasses—travel with the coolant, thus allowing for deposition and interaction throughout the entire system flow loop. This behavior must be modeled to successfully understand MSR system behavior during normal and off-normal operation events. In past work, SAM has been extended to include a passive species transport capability for modeling of MSR fission products that travel with the salt coolant. This work extends the model to include a drift-flux model to capture the tendency of the gas phase to move at a different velocity than the bulk liquid. Additional constitutive models have also been added for calculation of gas volume fraction, bubble radius, and interfacial area, which will be important for calculation of species phase migration in future coupling activities. Testing has been performed to ensure that the new model agrees with expected values for several simple models, but additional future work will be required to test the model for more complex geometries that will be needed for full-scale MSR simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Application of SCALE to Molten Salt Fueled Reactor Physics in Support of Severe Accident Analyses

As part of a US Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with SCALE and MELCOR, SCALE was used for the modeling and simulation of a molten salt-fueled reactor (MSR). SCALE capabilities for the modeling of MSR physics were demonstrated based on the Molten Salt Reactor Experiment (MSRE). Of primary interest were the determination of the system’s nuclide inventory, as well as the inventories in the various regions of the loop, considering that the fuel is continuously pumped through the system. This report contains discussions on the following: 1. Determination of the system-average fuel salt inventory considering fission gas removal in the off-gas system and noble metal removal through plating out at the heat exchanger using recent enhancements in SCALE’s depletion sequence TRITON, 2. Assessment of the nuclide spatial distribution throughout the loop using SCALE’s depletion solver ORIGEN, 3. Calculation of the core’s power profile, flux profile, temperature reactivity coefficients, and xenon reactivity using full-core calculations with SCALE’s Monte Carlo code KENO-VI. The results obtained with SCALE were post-processed to provide the MELCOR team with the core inventory and decay heat of the system, as well as the inventory and decay heat of individual regions in the loop, a zone-wise power profile, temperature feedback coefficients, and the xenon worth.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗