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 163 records · Page 9

Assessing High Burnup U-19Pu-10Zr Fuel Performance against Historical and Modeled Behavior

Advancing the deployment of sodium-cooled fast reactors (SFRs) requires thorough testing of metallic fuel pins under accident conditions to establish safe operational limits of high burnup fuel. To conduct transient testing, a comprehensive understanding of steady-state fuel behavior obtained through both experimental characterization and accurate predictive capabilities is needed. This study comparatively assesses the steady-state irradiation performance of two high burnup U-19Pu-10Zr fuel pins, DP-36 and DP-40, irradiated under prototypic fast reactor conditions in preparation for planned safety testing at the Transient Reactor Test Facility. Since DP-40 was designated for use in the test and DP-36 serves as its sibling pin, non-destructive, engineering-scale post-irradiation examinations (PIE) were conducted on both pins while destructive examinations were performed exclusively on DP-36. The results were then assessed against historical performance data from similar fuel pins irradiated in the Experimental Breeder Reactor-II. Additionally, the steady-state irradiation of each pin was modeled using the BISON fuel performance code to assess the accuracy of current modeling capabilities in predicting the baseline irradiation behavior. Non-destructive examinations included neutron radiography to measure fuel column elongation, gamma scanning to verify pin integrity and fission product migration, and profilometry to assess dimensional changes. Benchmarking against existing PIE data revealed consistent patterns in axial fuel column growth and cladding diametral strain, though both pins exhibited longer low-density “fluff” structures, which can have implications for core reactivity and source term calculations. Destructive examinations on DP-36 included fission gas release analysis and sectioning for optical microscopy, which showed more complex constituent redistribution patterns than the traditionally accepted 3-ring model. The axial evolution of fractional areas and porosities of each of the redistributed zones were quantified and presented. Modeling comparisons showed agreement in fractional fission gas release but consistently overestimated axial and radial swelling and disagreed with measured axial porosity patterns. These conservative overpredictions suggested that the pins would appear closer to failure or operational limits at the start of transient tests, potentially leading to higher strain accumulation during the transient. While conservative estimates provide safety margins, they can negatively impact fuel economics. A review of the swelling models identified areas for improvement in the gaseous swelling, solid swelling, and fuel hot-pressing models when applied to ternary fuel. The results of this study highlight the critical importance of conducting pre-test characterization on both test and sibling pins to accurately capture steady-state fuel behavior, providing a precise baseline for post-test evaluations and essential inputs for transient modeling of the planned experiments. The analysis also revealed significant data gaps that require further investigation to enhance the understanding and prediction of fuel swelling and pore dynamics. Collecting comprehensive data across different irradiation conditions, burnup levels, and fuel compositions are essential for refining existing models and developing mechanistic models for both binary and ternary metallic fuels, ultimately improving the integration of modeling and experimental approaches in accident testing.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

AGR-5/6/7 Irradiation Experiment Fission Product Mass Balance

This report presents the fission product mass balance for the AGR-5/6/7 TRISO fuel irradiation experiment. The fission product inventories deposited on capsule components outside of the fuel (e.g., stainless-steel shells, graphite holders, Grafoil disks, and associated hardware) were quantified as part of the post-irradiation examination (PIE) to assess the performance of this fuel. Comparisons were made between these inventories and depletion calculations, non-destructive measurements of fuel fission product inventories in the intact fuel compacts, prior AGR experiments, and results from among each of the five distinct AGR-5/6/7 capsules. The data served as estimates of the condensable fission products released from the fuel during irradiation. Excluding Capsule 1 (which experienced accidental damage during irradiation) and Capsule 3 (which was tested at very high irradiation temperatures of >1300°C), the results indicate that the AGR-5/6/7 fuel performed comparably to fuel from earlier AGR experiments. The mass balance results were also used to estimate the number of particles with in-pile SiC failures. Subject to the assumptions made in these estimates, and excluding Capsules 1 and 3, the in-pile SiC failure rates for AGR-5/6/7 are comparable to those observed in AGR-2.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Distribution functions in plasmas generated by a volume source of fission fragments

The role played by fission fragments and electron distribution functions in nuclear pumped lasers is considered and procedures for their calculations are outlined. The calculations are illustrated for a He-3/Xe mixture where fission is provided by the He-3(n,p)H-3 reaction. Because the dominant ion in the system depends on the Xe fraction, the distribution functions cannot be determined without the simultaneous consideration of a detailed kinetic model. As is the case for wall sources of fission fragments, the resulting plasmas are essentially thermal but the electron distribution functions are non-Maxwellian.

Deese, J. E.↗

The electron Boltzmann equation in a plasma generated by fission fragments

A Boltzmann equation formulation is presented for the determination of the electron distribution function in a plasma generated by fission fragments. The formulation takes into consideration ambipolar diffusion, elastic and inelastic collisions, recombination and ionization, and allows for the fact that the primary electrons are not monoenergetic. Calculations for He in a tube coated with fissionable material show that, over a wide pressure and neutron flux range, the distribution function is non-Maxwellian, but the electrons are essentially thermal. Moreover, about a third of the energy of the primary electrons is transferred into the inelastic levels of He. This fraction of energy transfer is almost independent of pressure and neutron flux but increases sharply in the presence of a sustainer electric field.

Hassan, H. A.↗

Electron distribution function in a plasma generated by fission fragments

A Boltzmann equation formulation is presented for the determination of the electron distribution function in a plasma generated by fission fragments. The formulation takes into consideration ambipolar diffusion, elastic and inelastic collisions, recombination and ionization, and allows for the fact that the primary electrons are not monoenergetic. Calculations for He in a tube coated with fissionable material shows that, over a wide pressure and neutron flux range, the distribution function is non-Maxwellian, but the electrons are essentially thermal. Moreover, about a third of the energy of the primary electrons is transferred into the inelastic levels of He. This fraction of energy transfer is almost independent of pressure and neutron flux.

Hassan, H. A.↗

Efficient Failure Probability Calculations and Modeling Interface Debonding in TRISO Particles with BISON

This document constitutes completion of the NEAMS milestone, which is titled Demonstrate efficiency improvements in TRISO failure probability calculations and the effect of debonding on fission product transport in TRISO fuel particles in Bison . In this report we present the development of the following: (1) the anisotropic elasticity model of pyrolytic carbon and extended material models using the local coordinate system for aspherical particle geometry; (2) a capability of modeling interface debonding in TRISO particles; and (3) an efficient high-fidelity approach to calculate the failure probability of TRISO particles. Benchmark problems and AGR-2 tests are used to verify the models and demonstrate these new capabilities.

36 MATERIALS SCIENCE↗

Spatial Burnout in Water Reactors with Nonuniform Startup Distributions of Uranium and Boron

Spatial burnout calculations have been made of two types of water moderated cylindrical reactor using boron as a burnable poison to increase reactor life. Specific reactors studied were a version of the Submarine Advanced Reactor (sAR) and a supercritical water reactor (SCW) . Burnout characteristics such as reactivity excursion, neutron-flux and heat-generation distributions, and uranium and boron distributions have been determined for core lives corresponding to a burnup of approximately 7 kilograms of fully enriched uranium. All reactivity calculations have been based on the actual nonuniform distribution of absorbers existing during intervals of core life. Spatial burnout of uranium and boron and spatial build-up of fission products and equilibrium xenon have been- considered. Calculations were performed on the NACA nuclear reactor simulator using two-group diff'usion theory. The following reactor burnout characteristics have been demonstrated: 1. A significantly lower excursion in reactivity during core life may be obtained by nonuniform rather than uniform startup distribution of uranium. Results for SCW with uranium distributed to provide constant radial heat generation and a core life corresponding to a uranium burnup of 7 kilograms indicated a maximum excursion in reactivity of 2.5 percent. This compared to a maximum excursion of 4.2 percent obtained for the same core life when w'anium was uniformly distributed at startup. Boron was incorporated uniformly in these cores at startup. 2. It is possible to approach constant radial heat generation during the life of a cylindrical core by means of startup nonuniform radial and axial distributions of uranium and boron. Results for SCW with nonuniform radial distribution of uranium to provide constant radial heat generation at startup and with boron for longevity indicate relatively small departures from the initially constant radial heat generation distribution during core life. Results for SAR with a sinusoidal distribution rather than uniform axial distributions of boron indicate significant improvements in axial heat generation distribution during the greater part of core life. 3. Uranium investments for cylindrical reactors with nonuniform radial uranium distributions which provide constant radial heat generation per unit core volume are somewhat higher than for reactors with uniform uranium concentration at startup. On the other hand, uranium investments for reactors with axial boron distributions which approach constant axial heat generation are somewhat smaller than for reactors with uniform boron distributions at startup.

Fox, Thomas A.↗

Neutron fluences and dose equivalents measured with passive detectors on LDEF

Neutron fluences were measured on LDEF in the low energy (< 1 MeV) and high energy (> 1 MeV) ranges. The low energy detectors used the 6Li(n,alpha)T reaction with Gd foil absorbers to separate thermal (< 0.2 eV) and resonance (0.2 eV-1 MeV) neutron response. High energy detectors contained sets of fission foils (181Ta, 209Bi, 232Th, 238U) with different neutron energy thresholds. The measured neutron fluences together with predicted spectral shapes were used to estimate neutron dose equivalents. The detectors were located in the A0015 and P0006 experiments at the west and Earth sides of LDEF under shielding varying from 1 to 19 g/cm2. Dose equivalent rates varied from 0.8 to 3.3 microSv/d for the low energy neutrons and from 160 to 390 microSv/d for the high energy neutrons. This compares with TLD measured absorbed dose rates in the range of 1000-3000 microGy/d near these locations and demonstrates that high energy neutrons contribute a significant fraction of the total dose equivalent in LEO. Comparisons between measurements and calculations were made for high energy neutrons based on fission fragment tracks generated by fission foils at different shielding depths. A simple 1-D slab geometry was used in the calculations. Agreement between measurements and calculations depended on both shielding depth and threshold energy of the fission foils. Differences increased as both shielding and threshold energy increased. The modeled proton/neutron spectra appeared deficient at high energies. A 3-D model of the experiments is needed to help resolve the differences.

LDEF Project↗

Complete β-decay patterns of 142 Cs, 142 Ba, and 142 La determined using total absorption spectroscopy

Background: The β decays of fission products produced in nuclear fuel are important for nuclear energy applications and fundamental science of reactor antineutrinos. In particular, nuclear reactor safety is related to the decay modes of radioactive neutron-rich nuclei, primarily via the emission of γ rays, neutrons, and electrons. Additionally nuclear reactors are the most powerful man-made source of antineutrinos emitted during the β decay of fission products. These antineutrinos are used to inspect fundamental properties of leptons as well as informing reactor operation. However, the majority of data on complex decays of fission products collected in the evaluated nuclear data repositories like Evaluated Nuclear Structure Data File (ENSDF) and Evaluated Nuclear Data Files (ENDF) are based on low-efficiency and often incomplete measurements resulting in questionable reference reactor antineutrino flux predictions, see the analysis by [Nichols, J. Nucl. Sci. Technol. 52, 17 (2015)]. Various assessments like the one done under the auspices of the [Yoshida et al., Assessment of Fission Product Decay Data for Decay Heat Calculations: A report by the Working Party on International Evaluation Co-operation of the Nuclear Energy Agency Nuclear Science Committee (Nuclear Energy Agency, Organization for Economic Co-operation and Development, Paris, France, 2007), Vol. 25], as well as by [Sonzogni, Johnson, and McCutchan, Phys. Rev. C 91, 011301(R) (2015)] and [Dwyer and Langford, Phys. Rev. Lett. 114, 012502 (2015)], list the A = 142 isobars with high cumulative fission yield among the important nuclei where data for reactor decay heat and/or antineutrino production should be verified and/or improved. Purpose: Here, our goal is to improve the quality of β -decay measurements and evaluate the impact of modified decay schemes on reactor decay heat and antineutrino energy spectra, for fission products along the A = 142 isobaric chain. This work is an in depth follow-up on [Rasco et al., Phys. Rev. Lett. 117, 092501 (2016)]. which presented briefly the impact of the corrected decay scheme of 142 Cs . Here, we extend the data to full isobaric decay chain including the daughter nuclei, 142 Ba and 142 La, and present more details on the 142 Cs results. Method: The decays of neutron-rich isobars of mass A = 142 produced by means of proton-induced fission of 238 U were measured using the Modular Total Absorption Spectrometer (MTAS) array on-line at the mass separator and Tandem accelerator at Oak Ridge National Laboratory. Results: The β -decay schemes for 142 Cs and 142 La were modified with respect to the nuclear data repositories. A small β-delayed neutron branching ratio for 142 Cs emitter was remeasured as $0.10^{+5}_{–3}% %. Improved precision on the measured half-lives is reported. Small corrections to the low-energy decay of 142 Ba are made. The β-decay patterns for 142 La and 142 Cs are presented. The decay heat release and cross section for the detection of reactor antineutrinos are deduced and compared to earlier results. Conclusions: The β-feeding pattern for 142 Cs having decay energy value $Q_β$ of over 7 MeV was substantially modified with respect to the current ENSDF entry. Smaller changes were encountered for 142 La, but since this A = 142 isobar also has a large cumulative yield in fission, the changes influence both decay heat and the antineutrino spectra. The previously known β intensities for 142 Ba decay ($Q_β$ value of 2.2 MeV) were verified and slightly modified. Overall, increased decay heat values and lower flux of antineutrinos interacting with matter are presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Uncertainty Quantification of Calculated Temperatures for the AGR 5/6/7 Experiment

This report documents the quantification of uncertainty of the calculated temperature data for the Advanced Gas Reactor (AGR) 5/6/7 fuel irradiation experiment conducted in the Advanced Test Reactor at Idaho National Laboratory in support of the Advanced Reactor Technologies? research and development program. Recognizing uncertainties inherent in physics and thermal simulations of the AGR 5/6/7 capsules, the results of the numerical simulations are used in combination with statistical analysis methods to improve qualification of measured data. The calculated fuel temperatures for AGR tests are also used for validation of the fission product transport and fuel performance simulation models. These crucial roles of the calculated fuel temperatures in ensuring achievement of the AGR experimental program objectives require accurate determination of the model temperature uncertainties. This report covers temperature uncertainty results for each of the five AGR 5/6/7 capsules. To quantify the uncertainty of calculated temperatures determined using the ABAQUS finite element heat transfer code, this study identifies and analyzes model parameters of potential importance to the calculated temperatures of fuel compacts and thermocouples. The selection of input parameters for uncertainty quantification is based on the ranking of their influences upon temperature predictions. Thus, selected input parameters include those with high sensitivity and those with the largest uncertainty. Propagation of model parameter uncertainty and sensitivity is then used to quantify the overall uncertainty of calculated temperatures. Measurement uncertainty, analysis of modeling assumptions, and expert judgment are used as the basis to quantify the uncertainty range for selected input parameters. The input uncertainties are dynamic, accounting for the effect of unplanned events and changes in thermal properties of capsule components over extended exposure to high temperatures and fast neutron irradiation. The sensitivity analysis performed in this work went beyond the traditional local sensitivity. Using experimental design, analysis of pairwise interactions of model parameters was performed to establish sufficiency of the time dependent first order (linear) expansion terms in constructing the temperature response surface. To achieve completeness, uncertainty propagation made use of pairwise noise correlations of model parameters. Furthermore, using an interpolation scheme over the input parameter domain, the analysis obtains time dependent sensitivity over the test campaign duration. This allows computation of uncertainty for the calculated fuel temperatures and the calculated graphite temperatures at thermocouple locations during the entire irradiation period.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Mechanistic nuclear fuel performance modeling of uranium nitride

Uranium mononitride (UN) is a nuclear fuel candidate for advanced reactor designs and an alternative being considered for light water reactors due to its higher thermal conductivity and uranium density than UO 2 . As with any nuclear fuel, swelling and fission gas release are important factors for safety, while also being some of the hardest phenomena to predict with a high degree of confidence. Getting a grasp on the gas swelling behavior and release is crucial to lower the barrier for UN utilization. An accelerated swelling rate at high temperatures observed experimentally, sometimes referred to as “breakaway swelling,” further complicates the prediction of fuel performance of UN. A mechanistic model has been developed using a multiscale approach to describe the intragranular and intergranular fission gas behavior. Lower-length-scale calculations have been employed to inform models of the gas and self-diffusion behavior, resolution rate, and bubble shape. Leveraging previous work on high burnup UO 2 , two populations of intragranular bubbles are considered; small bulk bubbles and larger bubbles located along dislocations. The dislocation bubbles were found to be crucial to the overall swelling behavior, and the breakaway swelling transition was associated with the transition in the gas atom diffusion mechanism from an irradiation-induced athermal diffusion regime at lower temperatures to an intrinsic thermal equilibrium regime at higher temperatures, accelerating the growth of the dislocation bubbles. Similarly, the threshold for fission gas release was associated with the grain boundary vacancy diffusivity surpassing the gas atom diffusivity at sufficiently high temperatures, allowing the over-pressurized grain boundary bubble to grow in size and interconnect. Using thermo-mechanical models with the fission gas model, two integral fuel pin assessment cases were simulated. Finally, this work demonstrates the ability of a multiscale approach to accelerate the understanding of advanced fuel forms when experimental data is limited.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Parametric analyses of planned flowing uranium hexafluoride critical experiments

Analytical investigations were conducted to determine preliminary design and operating characteristics of flowing uranium hexafluoride (UF6) gaseous nuclear reactor experiments in which a hybrid core configuration comprised of UF6 gas and a region of solid fuel will be employed. The investigations are part of a planned program to perform a series of experiments of increasing performance, culminating in an approximately 5 MW fissioning uranium plasma experiment. A preliminary design is described for an argon buffer gas confined, UF6 flow loop system for future use in flowing critical experiments. Initial calculations to estimate the operating characteristics of the gaseous fissioning UF6 in a confined flow test at a pressure of 4 atm, indicate temperature increases of approximately 100 and 1000 K in the UF6 may be obtained for total test power levels of 100 kW and 1 MW for test times of 320 and 32 sec, respectively.

Rodgers, R. J.↗

Machine learning guided prediction of solute segregation at coherent and semi-coherent metal/oxide interfaces

Investigation of semi-coherent metal/oxide interfaces with misfit dislocations using density functional theory (DFT) is computationally intensive to the point of being prohibitive, as it involves several hundreds to many thousands of atoms. In this study, we examined the solute segregation behavior at the Fe/Y 2 O 3 interface—a model interface for cladding applications in nuclear fission reactors—using a combination of DFT calculations and machine learning (ML) approaches. Both coherent and semi-coherent interfaces were considered. ML models were trained on DFT-calculated segregation energies to identify the key chemical, geometric and strain energy related features that govern solute segregation behavior at coherent Fe/Y 2 O 3 interfaces. Furthermore, it was found that ML models when trained on DFT calculated segregation energy of elements at a coherent interface, comprising of about a hundred-atom supercell, can predict the segregation energy of elements at a semi-coherent Fe/Y 2 O 3 interface (with multiple hundreds of atoms) at a fraction of computational cost (1/35th), with an accuracy comparable to DFT calculations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Neutronic Safety Analysis of Pu-238 Production at Idaho National Laboratory

This analysis was completed to support the irradiation of plutonium fuel services (PFS) targets in the NEFT in the Advanced Test Reactor (ATR) as a part of the campaign to restart domestic production of plutonium-238 used in radioisotope power systems (RPS) by the National Aeronautical and Space Administration (NASA) and Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Nuclear Infrastructure Program (NE-3). Referred to as the PFS-ATR-GEN1-NEFT experiment, the assembly was designed to hold 46 PFS targets in the NEFT. The scope of this paper is to outline the MOPY method used to calculate the heat generation rates (HGRs), flux, fission density, and quantify the viability of the target design for Pu-238 production in ATR.

07 ISOTOPE AND RADIATION SOURCES↗

Improving Fission Products at CARIBU: Near Field Detection (Q1/FY22 Quarterly Progress Report)

We worked with our collaborators at ANL to perform a measurement of the masses (and therefore energy difference) of the isomer and ground state for several fission products. The results of that measurement will have an impact on fission yields as well as on astrophysical calculations of r-process nucleosynthesis. We are currently analyzing the results of these measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Revisiting Source Convergence Diagnostics in the KENO Monte Carlo Neutron Transport Codes [Abstract]

Monte Carlo criticality transport codes, which rely on the power iteration procedure, are a fundamental tool for nuclear criticality safety practitioners in assessing the neutron multiplication factor (k eff ) for problems involving fissile material. In these calculations, ensuring the convergence of both the fission source distributions and the k eff estimate for accurate results is crucial. However, a converged k eff estimate does not necessarily mean the fission source distribution is also converged because the fission source and flux distribution may continue to evolve even after k eff convergence. Therefore, most Monte Carlo transport criticality codes now offer various diagnostic tests to assess fission source convergence in addition to the k eff convergence by analyzing the trends of these quantities over multiple generations.

AZURE↗

Comparison of NMC estimates with trans-stilbene, EJ-309, and He-3 detection systems

Neutron multiplicity counting (NMC) is a technique for the assay of fissile material. In this work, three detection systems are utilized for active interrogation assay of shells of the Rocky Flats shells (highly enriched uranium, 93% 235U) stacked from 13.25-54.92 kg assemblies. The singles ($R_1$) and doubles ($R_2$) rates are calculated with each system to estimate two sample parameters: $M_L$- the leakage multiplication and F - the sample fission rate. The estimated mass, m, is found by dividing F by a constant activity per unit mass. Since we are interrogating HEU, the α ratio of (α; n) to fission neutrons is taken to be zero. The system of equations to calculate these quantities was originally derived. The Neutron Multiplicity 3 He Array Detector (NOMAD) consists of 15 3 He tubes inside polyethylene and is the traditional, capture-based detection system for NMC. The polyethylene moderates incoming neutrons for thermal capture in individual tubes. The low gamma background, discrete capture signals, and high efficiency of the NOMAD are beneficial for NMC. However, the time to slow down neutrons to thermal energies leads to a system die away time on the scale of microseconds. Comparatively, the Rossi-alpha Measurements – Rapid Organic (n, γ) Discrimination Detector (RAMRODD) and the Organic Scintillator Array (OSCAR) are scatter-based detection systems. RAMRODD consists of 8, 5.08 by 5.08 cm EJ-309 liquid scintillators in 4 pairs, evenly-spaced with 90 degrees separation about the center of each assembly. OSCAR consists of a single array of 12, 5.08 by 5.08 cm trans-stilbene crystals aligned with the center of the assembly. Scatter-based systems detect fast neutrons without any moderation, leading to system die away times on the scale of tens of nanoseconds. This allows much shorter timing gate widths compared to thermal systems, thus increasing counting statistics of true correlated fission events. However, scatter-based systems are susceptible to neutron cross-talk when an incident neutron scatters off one detector and interacts in an adjacent detector, causing two seemingly correlated detection signals. The equations developed adjust for cross-talk to conduct NMC with scatter-based systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Gamma-ray spectrometry analysis methods for radioisotope dissolution and mixing for nuclear forensics applications

We report the production of surrogate nuclear explosive debris involves the dissolution of irradiated fissile material, elemental isolation, and recombination based on required fission product ratios. Gamma-ray spectrometry is used following each step to assess the radioactive fission product inventory. A high-purity germanium detector is used to measure each sample, and the detector efficiencies are determined prior to sample measurements. Through radiochemical techniques, elements are isolated and recombined into samples with specific fission product ratios. As a result, a matrix-based approach can account for imperfect chemical separations to calculate the proper mixture volumes to obtain intended volatile-to-refractory fission product ratios, or R-values.

07 ISOTOPE AND RADIATION SOURCES↗