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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Westinghouse Accident Tolerant Fuel Phase 2B with Higher Enriched and Higher Burnup Add-On Project Final Technical Report Deliverable Volume 2

The Westinghouse Electric Company LLC (Westinghouse) accident tolerant fuel (ATF) program funded by Westinghouse and the US Department of Energy (DOE) utilized chromium (Cr) coated zirconium alloy cladding with doped UO 2 (ADOPT(TM)) and uranium silicide (U 3 Si 2 ) high density/high thermal conductivity fuel for its lead test rod (LTR) program with irradiation beginning in 2019. ADOPT is a Cr 2 O 3 +Al 2 O 3 doped UO 2 pellet with increased oxidation resistance and density, and increased resistance to fission gas release. Due to the issues with U 3 Si 2 fuel reaction in pressurized water reactor (PWR) environments, uranium nitride (U 15 N) has been substituted for U 3 Si 2 as a long-term fuel option. Cr coated cladding with ADOPT fuel is the near-term Westinghouse EnCore ® fuel product. The long-term EnCore fuel product is SiGA ® SiC/SiC composite cladding with high density/high thermal conductivity UN fuel. In 2020, the higher burnup and higher enriched and (HBHE) program was integrated into the ATF program. The objective of this expanded program is to extend ATF burnups to at least 75 MWd/kgU to economically facilitate 24-month cycles in PWRs. The ATF program now includes ADOPT fuel with enrichments >5% 235 U. Over the past several years, Westinghouse has tested the Cr coated zirconium (Zr) and silicon carbide (SiC) claddings in Westinghouse Churchill autoclaves and the Massachusetts Institute of Technology (MIT) reactor. Additionally, the Cr coated cladding has been tested in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) and in LTR programs at the Doel 4 and Bryon 2 commercial reactors. High temperature tests at the state-of-the-art Westinghouse facilities in Churchill, PA and at Karlsruhe Institute of Technology (KIT) have been carried out to determine the time and temperature limits for the Cr coated zirconium claddings. These tests indicate that Cr coated Zr cladding can take temperatures up to about 1500°C for short periods of time without becoming totally oxidized. The main issue is the formation of the Cr-Zr eutectic at 1333°C resulting in the migration of this eutectic inward with the formation of ZrCr 2 which rapidly oxidizes to ZrO 2 on the outside of the tube. The Cr coated Zr cladding also delays the bursting of the tube, reduces the burst area, and reduces hydriding and embrittlement of the Zr which is a major benefit for reducing fuel fragmentation, redistribution, and dispersal (FFRD), the major licensing issue faced by HBHE. The manufacturing parameters for the SiC have been found to have a significant effect on the corrosion rate of the SiC in light water reactor (LWR) conditions and significant improvements have been made in the yield, quality, and oxidation resistance of the SiC cladding by identifying and tightening manufacturing process parameters. Current autoclave results for SiC composite claddings indicate that a corrosion rate of fewer than 2 micrometers per year can be achieved which meets corrosion requirements under normal operating conditions. High temperature steam oxidation tests at the state-of-the-art Westinghouse facilities in Churchill, PA and at Karlsruhe Institute of Technology (KIT) have been carried and determined that SiC cladding can withstand temperatures up to about 1800°C to 1900°C before excessive corrosion reactions begin, and that the SiC will not balloon and burst. ADOPT fuel pellets were incorporated into the Westinghouse ATF program because of their increased density and oxidation resistance. A topical report supporting the use of ADOPT pellets has been submitted to the NRC and will be approved shortly. ADOPT pellets have been a product in Europe for over 15 years. The additives work to make larger grain sizes which appear to reduce fission gas release during transients as well as increase the density of the pellet. Fuel rod and assembly design in preparation for the lead test rod (LTR) and lead test assembly (LTA) programs is underway as well as licensing efforts with the Nuclear Regulatory Commission (NRC). Finally, accident analyses coupled with economic evaluations for both operating savings as well as fuel savings have been initiated. Modular Accident Analysis Program 5 (MAAP5) calculations indicate that both Cr coated cladding and SiC composite cladding can increase the time to fuel rod loss of geometry for up to two hours longer than current Zr based cladding in a station blackout scenario. This additional time is due to the much higher oxidation resistance of the SiC and Cr coated cladding. These two hours can be used to implement additional response options instituted through the Diverse and Flexible Mitigation Capability (FLEX) program by reactor operators. Both ATF cladding options also reduce hydrogen production which dramatically reduces primary system and containment pressure and the risk of fission product release beyond containment in the unlikely event of an accident. The lower pressure in the system allows more time to feed cooling water to the core, resulting in the avoidance of fuel melting. The coping time is extended indefinitely if the modest water flow provided by FLEX continues. Experimental work on methods to rapidly fabricate SiC composite structures with high density and reduce the fabrication price of SiC fibers while maintaining a high level of performance is needed. Methods for mitigating or stopping the Cr-Zr eutectic are also needed to further improve the oxidation resistance of Cr coated Zr. Minimal (<1%) swelling of U 3 Si 2 and subsequent fission gas release has been demonstrated up to 20 MWd/kgU. Irradiation experiments with U 3 Si 2 fuel in ATR to determine these properties at 40 MWd/kgU were completed but post irradiation examinations were not done since U 3 Si 2 has been replaced by UN. UN has three issues that are being addressed. The first is increasing the oxidation resistance so that there is not excessive reaction up to ~1500°C. While UN increases the pellet density and additives to the UN postponed the temperature at which rapid oxidation occurred by ~100°C to 200°C, this is not enough, and pellet coating options are now being pursued. The second is developing a method to manufacture that does not require the multi-step process of UF 6 to UO 2 to UC to UN. Efforts are underway looking at UF 6 to UN 2 to UN and UF 6 to UF 4 to UN 2 to UN reactions. Finally, an economically acceptable process for potentially enriching the 15 N content of the nitrogen used to make UN to >95% 15 N was identified though not experimentally validated. In addition to the work supported by the DOE, research and testing activities are being carried on in a world-wide effort funded by many countries such as Sweden, United Kingdom, Belgium, Netherlands, Spain, Germany, Japan, and France. This work is being facilitated through the Westinghouse led Collaboration for Advanced Research on Accident Tolerant Fuel (CARAT) program. Annual meetings were organized in 2018 and 2019 by Westinghouse as a venue for presentation of this work and to provide for the cross-fertilization of ideas among the many outstanding researchers in the ATF area. No meetings were held in 2020 and 2021 due to the COVID-19 related restrictions on traveling and gatherings. Since SiC, coated cladding, and high-density fuel options are not currently used in the nuclear industry; support from the government and industry members is needed to further the significant effort of setting new standards. The same is true for the Nuclear Regulatory Commission (NRC) which must review and approve the commercial use of these new fuels since all current regulations are oriented toward Zr/UO 2 fuel. Several meetings have been held with the NRC to generate a fast-track approach to licensing ATF using a combination of atomic modeling, in-rod sensors, and in-reactor testing. This approach was memorialized in an “Accelerated Fuel Qualification White Paper” prepared by the Accelerated Fuel Qualification Working Group lead by General Atomics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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↗
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↗
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↗
DOE OSTI2021

Analysis of Approximations in Modeling of BWR Bundle Void Distributions

In boiling water reactors, complex heterogeneous bundle designs, control blades adjacent to the corner of bundles, and the presence of boiling can lead to complex internal void distributions. A few approximations exist to model these void distributions. They could be modeled using a 1D axial solver in which each axial node is assumed to be at an average void, or each pin cell could be modeled with its own void concentration. In the latter case, the void could be discretized in pin-centered or coolant-centered channels. The goal of this project was to quantify the effect of using the different approximations for modeling internal void distributions on neutronics calculations. Using 3D void distributions calculated with CTF, Monte Carlo Neutral Particle (MCNP) transport code models were created for GE-9 and GE-14 lattices. For each model, the internal void distribution from CTF at a given axial node was selected, and a lattice calculation was carried out with MCNP. Comparisons between models using a lattice-averaged void, or using a void distribution in coolant-centered channels, showed large differences in reactivity which in some cases were well above 1,000 pcm, and it also showed differences in normalized fission rates greater than 20%. It was also found that using a lattice average void can lead to a significant difference in the estimation of the worth of a control blade. The differences found when comparing results from models using pin-centered and coolant-centered channels were up to 200 pcm in reactivity and up to 1.4% in the normalized fission rates. In addition to these two sets of comparisons, MCNP models were set up so that each subchannel had a saturated liquid component around the fuel pins and a saturated vapor component in the center to approximate annular flow. In comparison to the models using coolant-centered subchannels, up to 1–3% differences in normalized fission rates could be found.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Niowave Neutron Source Converter: Lead-Bismuth-Eutectic (LBE) Windowless Target Design and Evaluation

Los Alamos National Laboratory (LANL) is working with Niowave on the design and evaluation of their lead-bismuth-eutectic (LBE) windowless target (i.e., neutron source converter). Niowave plans to use 200 kW electron beam at 40 MeV beam energy to produce neutrons by photonuclear reaction with LBE. Then, the neutrons undergo fission at the surrounding uranium target assembly (UTA) to produce Molybdenum 99 (Mo-99) as a fission product, which eventually decays to Technetium-99 (Tc-99m). Tc-99m is one of important radioisotopes that is used for medical diagnostics. LANL conducted 3D multiphysics analysis for the Niowave neutron converter design and provided design assessment in thermohydraulic aspects. LANL conducted radiation transport calculations using Monte-Carlo N-Particle (MCNP) code with unstructured meshing scheme. The 3D volumetric heating profiles in the LBE and Stainless-Steel (SS) housing were imported into multiphase computational fluid dynamics (CFD) to obtain 3D temperature profiles of LBE and SS through conjugate heat transfer (CHT) analysis. The key findings are: LBE film thickness at the center of the beam is approximately 1.6 cm with a maximum LBE velocity of approximately 1.8 m/s, which is below a 2 m/s limit to avoid erosion issues on supporting structures; Heat deposition in the LBE peaks at ~1 cm depth from the LBE free-surface because of the forward interactions of electron, photon, and neutron with LBE; LBE maximum temperature is ~360 °C, which is below LBE evaporation initiative temperature, ~450 °C; LBE-SS interface temperature is ~350 °C, which is below the safety thermal limit to prevent severe corrosion on SS. The results indicate that Niowave’s neutron converter design satisfies both hydraulic and thermal criteria for safe operation. By virtue of such computational analysis, Niowave can move toward establishing an experimental setup to experimentally test their LBE neutron converter. The following sections describe the detailed work done by LANL.

43 PARTICLE ACCELERATORS↗