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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2021

High-burnup fuel stress analysis prior to and during a LOCA transient

High-burnup fuel fragmentation and pulverization have been experimentally observed during simulated loss-of-coolant accident (LOCA) transient testing of high-burnup (>62 GWd/tU) nuclear fuel. The cause of this phenomenon is generally accepted to be the formation of stresses arising from the off-normal conditions. However, the primary mechanism(s) contributing to the stress formation is/are still being debated. Leading hypotheses suggest that the concurrence of fission gas bubbles and high burnup structure on the pellet periphery (i.e., the rim) leave the fuel susceptible to pulverization. Other hypotheses suggest that pulverization is artificially inflated during externally heated tests that do not reproduce prototypic fuel stresses and radial temperature profiles as the fuel transitions from full-power nuclear heating to decay heating. This work performs finite element analyses of fuel pellet temperatures and macroscopic stresses prior to and during a LOCA transient to identify (1) operating conditions that may exacerbate the fuel pellet stress state, (2) differences in predicted stresses with nuclear vs. furnace heating, and (3) how fuel cracking impacts the subsequent stress state in the pellet fragments. Analyses predict that macroscopic transient fuel stresses resulting from constrained thermal expansion will remain below the fracture stress of fresh UO 2 once the fragment size reaches ~4 mm or less. Calculated fuel stresses from simulated LOCA tests performed with representative nuclear vs. electrical heating rates are shown to be equal within 5% and essentially unchanged for all times after the first 27 seconds of the transients. Because the fuel stresses do not change once the temperature profile becomes uniform, macroscopic transient stresses would be no different than those following a normal shutdown. Because pulverization of high burnup fuel does not occur prior to a LOCA test, pulverization of high burnup UO 2 cannot be initiated by differential thermal stresses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Sensitivity analysis of homogenized cross section in AP1000 lattices

The two-step method for light water reactor simulation consists of performing lattice-level calculations to determine homogenized properties of the lattice for a variety of configurations and burnups then performing a core-level calculation which relies on those pre-calculated lattice properties. Calculating the homogenized lattice properties for the specific conditions needed in the core-level calculation relies on interpolating from the pre-calculated lattice properties. Although many studies have explored different models that can be used for this interpolation, there is a need to explore how to appropriately parameterize the lattice for those interpolation models. The present study uses linear sensitivity methods to determine the sensitivity of homogenized cross sections to both instantaneous and historical values of state variables. The results from this study can be used to inform parameterizations which more appropriately capture history effects in homogenized cross section interpolation. Two different lattices are used for this analysis from the AP1000 reactor, one includes burnable absorbers. It was found that early in the cycle, the historical values of state variables have little effect on the instantaneous value of the homogenized cross sections explored here - especially in the case of the lattice containing burnable absorbers. However, it was found that later in the cycle, the thermal fission homogenized cross section was sensitive to the historic states of the lattice. In fact, using the average of a state variable as a predictor for interpolation models may be insufficient for capturing historic effects due to more recent states having larger sensitivity measures than less recent states. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Sum-of-Fractions Method

Sum-of-fractions is a method intended to make sure a subcritical margin for aqueous solutions and slurries of fissionable isotopes exists. The method indicates that a system is subcritical if the sum of the ratios of the mass of each isotope (in a mixture) to its individual minimum subcritical mass limit is less than or equal to one. Historically, the basis of the sum-of-fractions has been derived from allowances given in the American National Standards Institute (ANSI)/ American Nuclear Society (ANS)-8.15-1981. However, the allowance was removed in ANSI/ANS-8.15-2014 due to a lack of technical basis. A methodology was developed to assess the validity of using the sum-of-fractions for water- or polyethylene-moderated systems for the following nuclides: 232U, 233U, 234U, 235U, 237Np, 236Pu, 238Pu, 239Pu, 240Pu, 241Pu, 242Pu, 241Am, 242mAm, 243Am, 242Cm, 243Cm, 244Cm, 245Cm, 246Cm, 247Cm, 249Cf, and 251Cf. The methodology uses available benchmark data for mixtures of 233U, 235U, and 239Pu to establish the calculational margin, and a mass limit reduction to establish the margin of subcriticality. Water- or polyethylene-moderated and -reflected mixtures containing the nuclides are evaluated with the code system, SCALE 6.2.4. Including the calculational margin, subcritical mass limits for each nuclide were computed for optimally water- or polyethylene-moderated and fully reflected systems. These masses were used to create nuclide mixtures in which the sum of the mass to subcritical mass limit ratios is one. The various nuclide mixtures were modeled over a range of moderation and demonstrate the keff does not exceed the calculational margin. For additional assurance of subcriticality, a significant mass reduction is applied to each computed minimum critical mass of the nuclides without adequate benchmark data consistent with the method in ANSI/ANS-8.15-2014.

criticality safety, Actinide↗

DOE OSTI2024

Nuclear Uncertainties Associated with the Nucleosynthesis in Ejecta of a Black Hole Accretion Disk

Abstract The simulation of heavy element nucleosynthesis requires input from yet-to-be-measured nuclear properties. The uncertainty in the values of these off-stability nuclear properties propagates to uncertainties in the predictions of elemental and isotopic abundances. However, for any given astrophysical explosion, there are many different trajectories, i.e., temperature and density histories, experienced by outflowing material, and thus different nuclear properties can come into play. We consider combined nucleosynthesis results from 460,000 trajectories from a black hole accretion disk and find the spread in elemental predictions due solely to unknown nuclear properties to be a factor of a few. We analyze this relative spread in model predictions due to nuclear variations and conclude that the uncertainties can be attributed to a combination of properties in a given region of the abundance pattern. We calculate a cross-correlation between mass changes and abundance changes to show how variations among the properties of participating nuclei may be explored. Our results provide further impetus for measurements of multiple quantities on individual short-lived neutron-rich isotopes at modern experimental facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗

DOE OSTI2024

Sodium-Cooled Fast Reactor Reference Plant Model

This report details the progress of Idaho National Laboratory (INL) in creating a reference plant multiphysics model for the Advanced Burner Test Reactor (ABTR). This model was developed under Task 13 of the U.S. Nuclear Regulatory Commission project “Development and Modeling Support for Advanced Non-Light Water Reactors,” and is an extension of the reference plant model developed in Task 4b, which was improved upon in the following ways. (1) The discrete ordinates method was used in lieu of the super-homogenization (SPH)-corrected diffusion approximation in order to better capture the anisotropic scattering contribution and the neutron leakage change due to thermal expansion. (2) The novel neutronic spatial discretization approach, termed the ring-heterogeneous (RH) approximation, was conceptualized and introduced to capture the differential expansion of the materials in the core. This new technique proved capable of preserving fission rates and maintaining the eigenvalue within 2.5% and 266 pcm with 9 neutron energy groups, respectively. Separating the different materials in the core enables the differential expansion of materials to be explicitly accounted for, eliminating the need for problem-specific cross-section functionalization techniques. (3) The SAM model for the core and system thermal-hydraulics analysis was updated to include 61 channels instead of just four representative ones. This enables users to obtain improved spatial resolution for sodium temperature and density scalar fields. (4) All the mesh files were created via the Multiphysics Object-Oriented Simulation Environment (MOOSE) Reactor module, eliminating all reliance on external tools for mesh creation. (5) Finally, the fuel axial expansion now leverages the HT9 and UPuZr material properties that have been validated against experimental data. The reference plant model was used to perform a full-core unprotected loss of flow (ULOF) transient calculation, including neutronics, thermal and mechanical feedback mechanisms. Future work will be devoted to further enhancements of the model. Potential improvements to the model include the addition of the control rod driveline expansion feedback and the upgrading of the support plate model so as to explicitly include 3D effects. Additionally, a Nuclear Energy Advanced Modeling and Simulation funded parallel effort has completely automated the creation of the ring-heterogeneous (RH) mesh from the fully heterogeneous (FH) geometry, thus maximizing user friendliness for the sodium fast reactor sodium-cooled fast reactor (SFR) workflow and will be incorporated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

SAM Code Development for Source Term Modeling in Fluoride-salt-cooled High-temperature Reactors

The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. Many concepts leverage higher operating temperatures for improved efficiency with a variety of coolants and structural materials to support those needs. Such design choices may also pose unique radiological source term risks. There are continuous efforts to incorporate new physics and phenomena relevant to advanced reactor concepts, including capabilities to enable system-level source term modeling. The Pebble-Bed Fluoride-salt-cooled High-temperature Reactor (PB-FHR) is a promising candidate among advanced nuclear reactor concepts with its improved passive safety characteristics and high thermal efficiency. In addition to past efforts to support the development and utilization of SAM for PB-FHR safety analysis, the species transport modeling capabilities in SAM have been extended to simulate source term phenomena in concepts like the PB-FHR. The PB-FHR concept utilizes pebble-form TRISO fuel and FLiBe salt coolant to provide robust barriers to the release of almost all fission products and radiological source terms. However, tritium poses a unique radiological risk due to its significant production from neutron interactions with 6 Li and 9 Be in the FLiBe as well as its high mobility at elevated temperatures where it can permeate through structural metals. Secondly, the use of graphite at elevated temperatures poses the risk of oxidation damage to any structural components or fuel pebbles that are exposed to accidental air ingress. This report summarizes progress made in modeling these phenomena in SAM, which leverages the system-level multi-physics thermal hydraulic simulation to support an effective engineering-scale source term modeling capability. A tritium transport model is developed to simulate the various transport pathways in FHR and MSR concepts and is presented through verification and validation examples as well as demonstrations simulating experimental test loops and reference plant FHR models. A graphite oxidation model is also introduced to calculate local oxidation rates and is presented with initial validation results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Effect of Nuclear Data Covariances on Integral Experiment Design with Sensitivity and Uncertainty Analysis

Washington River Protection Solutions (WRPS) uses MCNP6.2 and the Whisper code for criticality safety analyses of the Hanford Tank Farm. Together the codes derive baseline upper subcritical limits (USLs) for the waste models using experimental benchmarks. Whisper returns higher USLs, i.e. , has less of a conservative penalty, when the neutronic similarity of the experimental benchmarks to the application is high. Unfortunately, few critical benchmarks have high similarity to the Hanford tanks. The waste in the tanks is highly dilute in plutonium and contains large masses of weakly neutron-absorbing elements like iron and manganese. Experimental benchmarks typically have low sensitivity to these absorbers because they are present as structural materials. Lacking similar benchmarks, new Thermal Epithermal eXperiment (TEX) configurations with high Pu content and interstitial iron absorbers have been designed for the criticality safety validation. The features of the design have been iterated upon to maximize the similarity between the experiment and different Hanford waste models. The similarity is quantified with sensitivity analysis and uncertainty quantification using the representativity coefficient, or c k . The representativity calculation requires nuclear data covariances, which may differ between nuclear data libraries and between library versions. Because of these variations, the optimal design may depend on the nuclear data covariances library. A scenario can be envisioned where an experiment is designed, and c k is maximized, with one set of covariance data. However, when the covariance data is changed, say from ENDF/B-VII.1 to ENDF/B-VIII.0, and the benchmark is used in a criticality safety evaluation, the experiment becomes suboptimal with respect to c k . In this paper, we present how the optimal design of the new TEX experiments varied depending on the nuclear data covariances used to calculate c k . We compare ENDF/B-VII.1 and ENDF/B-VIII.0, as if the library had been updated since the design of the experiment. Additionally, we use JEFF3.3 to simulate if the covariance data of a different library had been used. The results show that the covariances do have an important effect on the designs, less so for thermal systems (where the data are more consistent between evaluations) and more so for epithermal systems where more differences exist.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Godiva IV Simulated Radiation Field Characterization and Variance Reduction

Godiva IV is a system comprised of highly enriched uranium alloyed with molybdenum in the form of fuel plate rings. The reactor, along with its predecessors, was designed with the unique ability to satisfy interests in the super-prompt-critical reactor operation space. Originally, the reactor was part of the Los Alamos Critical Experiments Facility (LACEF) at Technical Area-18 (TA-18). The radiation field around Godiva at this facility was well characterized and understood. As a fast neutron system, the neutron spectrum in and around Godiva was close to a Watt Fission spectrum. The Kiva where Godiva IV was located at LACEF was made of thin, sheet metal walls which did not contribute significantly to the neutron spectrum. Following the transition of LACEF to the National Critical Experiments and Research Center (NCERC) in Nevada, Godiva-IV was moved from TA-18 to the Device Assembly Facility (DAF) at the Nevada National Security Site (NNSS). Part of this move brought renewed interest in radiation field characterization. The new facility introduced significant changes to the environment surrounding Godiva, and preliminary foil irradiation results suggested that the room contribution to the neutron spectrum was significant. Unlike at TA-18, a large thermal neutron signature was added to the fast spectrum from Godiva due to significant room return. A primary goal due to the additional complexity that the room return adds to the Godiva IV radiation emission spectrum was the development of an efficient Monte Carlo N-Particle (MCNP) calculation capable of characterizing the neutron spectrum anywhere in the room around Godiva. A campaign of activation foil irradiations and analysis were completed to support the validation of the MCNP model. The modeling of these foils in MCNP can be easily done with a standard volumetric neutron flux tally. However, given the multitude of locations and reaction rates to be modeled, further steps must be taken to increase the efficiency of these calculations in MCNP. During this study, a benchmark model currently under development for Godiva IV was used. A qualitative assessment of the thermal neutron contributors was performed using spatial neutron distribution plots. Additional detail was added to the model based on the qualitative results showing the thermal spectrum’s large sensitivity to hydrogenous material. Neutron energy spectra was evaluated at discrete locations in the room around Godiva to quantify the relative contribution of various components. It was discovered that the concrete walls are the largest contributor to the thermal signature, with minor contributions from plastic components surrounding Godiva. Following these results, two different variance reduction techniques were implemented to improve the problem efficiency in these calculations. In the first approach, an F5 point detector tally was implemented in the standard Godiva IV criticality problem. The second approach involved a weight-window generator implementation with an F5 point detector tally in a fixed source problem. The weight window implementation reduced the runtime from 42739.55 minutes to 1803.34 minutes (computer time), compared to the F5 KCODE implementation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Nuclear level density and γ -decay strength of Sr 93

This work presents the first experimentally determined nuclear level density and γ-ray strength function of the short-lived fission product 93 Sr, accomplished using the β-Oslo method. Direct measurement of the 92 Sr(n, γ) 93 Sr cross section is not currently possible, as the half-life of 2.66 hours is too short; instead, 93 Sr was formed through β decay of 93 Rb to excitation energies around the neutron separation energy. The γ-ray spectra were measured using a total absorption spectrometer at the National Superconducting Cyclotron Laboratory (NSCL) at Michigan State University (MSU). The statistical properties of the 93 Sr nucleus were experimentally determined, including the γ-ray strength function and nuclear level density. At low energies, the γ-ray strength function exhibits a constant γ-decay strength, rather than a slightly increasing strength with decreasing γ-ray energy as had been previously observed for several nuclei in this mid-mass region. Finally, these statistical properties were then implemented in the reaction code TALYS1.95 to calculate the 92 Sr(n, γ) 93 Sr cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Scoping Study for Fast Flux Testing in the Advanced Test Reactor

The value of fast spectrum reactors remains prominent in the nuclear technology portfolio. The performance of these reactors can be maximized with advancements in nuclear fuel technologies, but development of these technologies is currently held back by lack of fast spectrum test reactors available to the United States. Spectral modification of experiment positions in the thermal spectrum Advanced Test Reactor (ATR) has long been used to support fast reactor fuel development, but these methods have not been progressed to their full potential. This study investigated the use of concentric rings of aluminum-clad fuel plates in ATR flux traps and thermal neutron absorbing filters to increase fast neutron flux on test specimens. This concept was termed the Boosted Energy Advanced Spectrum Test (BEAST). This approach will enable irradiation of advanced fuel designs in prototypic-length fuel pins and representative flux environment to support post irradiation exams, enable transient testing, and produce the type of data that will permit lead test assembly irradiations in true Sodium Fast Reactors (SFRs) when they become available. Neutronic predictions were performed to investigate BEAST design options and thermal hydraulic models were produced to ensure feasibility of BEAST. Two versions were considered based on the geometric limitations of ATR’s small and large flux traps. The small version was found to be preferable due to slightly higher fast flux and fast-to-thermal neutron ratio. Perhaps more influentially, the small flux trap option was also preferred to avoid conflict with ongoing very high temperature reactor fuel irradiation programs in ATR’s large northeast flux trap. The small flux trap option provided less than half the test volume of the large version, but still had adequate volume for seven SFR pins in cross section which could be stacked two-high in ATR’s 1.2m long core to accommodate up to 14 EBR-II size pins. The preference for the small flux trap configuration should be revisited if additional collaborative test programs emerge with the need to irradiate a significant volume of additional specimens. Calculations were performed regarding a lithium deuteride ring to convert thermal neutrons into 14 MeV fusion neutrons. At the time this report was written these calculations were partially complete and it remains to be seen whether the concept would be worth including in BEAST. Given the preference for the small flux trap option, which does not afford enough room for the 14 MeV ring, it was concluded to defer future work on the lithium deuteride ring. This decision could be revisited if fusion material research programs emerge for collaborative testing in BEAST. A cadmium-lined specimen holder design was found to be adequate in filtering thermal neutrons and preferred over other neutron absorbers based on past experience with cadmium baskets. It was acknowledged that cadmium-bearing hardware would become depleted and need to be replaced occasionally, which appeared feasible from a mechanical design perspective. Neutronic studies investigated different enrichment levels in the booster fuel using uranium-molybdenum alloy dispersion fuel which has performed well in past ATR irradiations. Both options were able to drive fuel pins to SFR-like fission heating rates. The high enriched booster fuel option outperformed the low enriched option by ~20% on key metrics including fast flux and fast-to-thermal ratio, but the low enriched option was favored in order to broaden options for potential fuel suppliers. The preferred BEAST design options including cadmium filter with low enriched booster fuel in the small flux trap configuration was predicted to achieve 6.2E14 n/cm2sec fast flux (>0.1 MeV) with a fast-to-thermal ratio of 44.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Evaluation of the Sum-of-Fractions Methodology for Water and Polyethylene Moderated Systems

Sum-of-Fractions is a method intended to assure a subcritical margin for aqueous solutions and slurries of fissionable isotopes. The method indicates that a system is subcritical if the sum of the ratios of the mass of each isotope in a mixture to its individual minimum subcritical mass limit is less than or equal to one. The basis of the Sum-of Fractions has historically been derived from allowances given in ANSI/ANS-8.15-1981. However, the allowance was removed in ANSI/ANS-8.15-2014 due to a lack of technical basis. A methodology was developed to assess the validity of using the Sum-of-Fractions for water or polyethylene moderated systems for the following nuclides: 232 U, 233 U, 234 U, 235 U, 237 Np, 236 Pu, 238 Pu, 239 Pu, 240 Pu, 241 Pu, 242 Pu, 241 Am, 242 mAm, 243 Am, 242 Cm, 243 Cm, 244 Cm, 245 Cm, 246 Cm, 247 Cm, 249 Cf, and 251 Cf. The methodology uses available benchmark data for mixtures of 233 U, 235 U, and 239 Pu to establish the calculational margin, and a mass limit reduction to establish the margin of subcriticality. Water or polyethylene moderated and reflected mixtures containing the nuclides are evaluated with SCALE 6.2.4. Including the calculational margin, subcritical mass limits for each nuclide were computed for optimally water or polyethylene moderated and fully reflected systems. These masses were used to create nuclide mixtures in which the sum of the mass to subcritical mass limit ratios is one. The various nuclide mixtures were modeled over a range of moderation and demonstrate the k eff does not exceed the calculational margin. For additional assurance of subcriticality, a significant mass reduction is applied to each computed minimum critical mass of the nuclides without adequate benchmark data consistent with the method in ANSI/ANS-8.15-2014.

07 ISOTOPE AND RADIATION SOURCES↗

NASA NTRS1974

Roche limit of a solid body

Tidal fission of both impacting and orbiting linear elastic solid bodies based on Kelvin's theory of earth tides is considered. It is shown that there can be more than one mutually exclusive modes of fracture - the particular mode in which a body fractures depending on its size and strength. The analysis gives a vivid picture of the propagation of the fracture with a decreasing distance from the planet. Expressions for the initiation and completion of fracture are obtained which are displayed graphically for a rigid body. The effect of elasticity on the breakup altitude is discussed. For orbiting solid bodies, the study gives the upper limit of the breakup altitude as 0.38R (where R is the radius of planet), which is much less than the value 1.44R used for such bodies in the past. The results presented include a previously given theory by Sekiguchi as a part. For the special case of a liquid body, comparison is made with Roche's calculation and the difference explained.

Aggarwal, H. R.↗

NASA NTRSDate not supplied

An Extremely High Isp Spacecraft Propulsion System

Specific Impulse, Isp, is a measure of a rocket engine’s efficiency. It is calculated relative to the Earth’s gravitational field, where Isp = ve/ go, with go= 9.8 m/s2 and the escape velocity of the propellant, ve, in m/s. Chemical rockets have ve < 4.4x103 m/s and Isp < 450 seconds. As an alternative, the NASA Glenn Research Center developed multiple generations of Solar Electric Propulsion (SEP), high Isp, ion engines using Xe gas as a propellant. Consequently, over 100 SEP Ion Thrusters provide geo-synchronous station keeping along with deep space probes like Deep Space One and Dawn. These have ve » 2.9 x 104 m/s with Isp = 3x103 seconds. These thrusters have continuous operating lifetimes of thousands of hours allowing continuous acceleration making up for the very low thrust. Chang-Diaz’ Variable Specific Impulse Magnetoplasma Rocket (VASIMR) has the potential for four times the propellant escape velocity and four times the specific impulse1. Unfortunately, this comes at a tremendous electrical power cost, estimated at 200 kWe for maintaining the International Space Station in Low Earth Orbit (LEO). Although nuclear fission and fusion reactors2 have been suggested for powering nuclear thermal propulsion (NTP) it only doubles the Isp over chemical rockets but with comparable thrust. Instead, we propose using Lattice Confinement Fusion (LCF) reactions.

High Isp Propulsion↗

DOE OSTI2022

Numerical study of multi-component flow and mixing in a scaled fission product venting system

Numerical simulation of isothermal multi-component gas flow was conducted to study fluid flow and mixing in a scaled low-pressure, low-temperature test facility for the fission product venting system (FPVS) of a gas cooled fast reactor (GFR). The geometry of FPVS test facility was an open loop including gradual expansion coupling and two 90°pipe elbows. First, the validation of large eddy simulation (LES) wall-adapting local eddy-viscosity (WALE) turbulent model was performed for transitional flow in a circular pipe with gradual expansion. The reactingFOAM solver in OpenFOAM v8 was employed. By introducing appropriate turbulence disturbance at the flow inlet, such as turbulence intensity, the numerical results showed good agreement in (1) the velocity profile downstream of the pipe expansion measured by particle image velocimetry (PIV) in this study, and (2) the reattachment length reported in the literature. Using this validated model, the numerical simulation of flow and component mixing was performed for the FPVS with an inlet flowrate corresponding to a Reynolds number of 2,400 to investigate the flow behavior and component mixing. Standard deviation of the mass fraction of component, referred to as absolute mixing index (AMI), was calculated to quantify the mixing of components, showing that the component mixing length determined by AMI is related to the reattachment length downstream of the expansion. Additionally, this study was able to identify the best location in the FPVS test facility to measure the velocity and concentration profies where multicomponent flow is well mixed.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

A Fast Reactor Irradiation Experiment Design in the ATR

Modern modeling techniques were used to investigate a proposed method for fast neutron irradiations in an existing thermal-spectrum reactor, the Advanced Test Reactor (ATR). This method builds upon pre-existing ideas, where fast flux is increased by surrounding the specimens with fissionable “booster fuel” but diverges from historical approaches by using an already developed fuel element design used in the Belgian Reactor 2 (BR2) as the booster fuel while leveraging modern 3-D modeling and simulation techniques. Design evaluations and neutronics simulations were performed to evaluate the performance of a BR2 fuel element irradiated in an ATR flux trap with test pins in the central channel of the BR2 fuel element. These efforts have yielded promising results. Adding a BR2 fuel element in the northeast (NE) flux trap of ATR was predicted to result in a 150% increase to the incident fast neutron flux with a fast (>0.1 MeV) to thermal (<0.625 eV) neutron flux ratio ranging from approximately 50 to 150, dependent on the material used for thermal neutron filtering and volume of moderator within the central channel of the BR2 element. The predicted annual fast neutron fluence (>0.1 MeV) ranges from 7.9 × 1021 to 9.1 × 1021 n/cm2. Given the relatively large fast to thermal neutron flux ratio, the calculated radial power profiles within 4.3 mm outer diameter U10Zr fueled specimens irradiated within the BR2 booster fuel element are adequate representations of those within fast neutron reactors. The predicted radial power profiles are not flat, but they are more prototypic than those seen in advanced fuels tests which began in ATR in 2003. Another distinct advantage of this experiment design is that full-scale test pins can be irradiated to augment the ongoing series of reduced scale advanced fuels tests. The proposed experiment design irradiated within a BR2 fuel element in a flux trap of ATR offers an improved alternative to the current testing of advanced reactor fuels in ATR. Selected results from this design evaluation are presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Thermodynamic evaluation and modeling of the CeF 3 -LaF 3 -LiF-NaF system

The thermodynamic properties of the LiF-NaF-CeF 3 -LaF 3 pseudo-quaternary system have been evaluated, and models for the relevant phases have been revised or newly developed. The solid solubility of NaF in LaF 3 and CeF 3 was described in the model based on reported phase-equilibria observations and results from density functional theory calculations. The Gibbs energy functions for the NaCeF 4 and NaLaF 4 phases have been updated with recently reported thermochemical measurements. The Gibbs energy functions and model relations for the LiF-NaF-LaF 3 and LiF-NaF-CeF 3 pseudo-ternary systems were successfully optimized based on thermochemical properties and phase equilibria observations.

CALPHAD↗

DOE OSTI2025

Re-evaluating the prompt fission neutron spectrum of spontaneously fissioning 252 Cf

The prompt fission neutron spectrum (PFNS) of spontaneously fissioning 252 Cf is a Neutron Data Standards observable. Nearly all fission spectra of actinides were measured relative to it, using efficiencies derived from it, or analyzed with simulations validated by it. The current Standards evaluation was published by W. Mannhart in 1987. It could not be updated because the evaluation input, experimental mean values and covariances, were lost. First, we attempt to reproduce it. However, Mannhart’s evaluation can only be reproduced within its one-σ uncertainties as some of its aspects (e.g., experimental covariances, rejected data points) remain unknown. Therefore, a new evaluation is presented: We revisit all existing experimental 252 Cf(sf) PFNS data, including those published after the release of the current Standards evaluation, and re-estimate associated covariances. The newly evaluated 252 Cf(sf) PFNS differs distinctly from Mannhart’s below 300 keV and extends it to lower and higher outgoing neutron energies (500 eV–25 MeV). The new evaluated uncertainties are larger from 3–9 MeV and smaller otherwise. Spectrum averaged cross sections of importance to the International Reactor Dosimetry and Fusion File community calculated with the new spectrum are close to those calculated with Mannhart’s evaluation and agree with experimental values well within their uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗