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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Non-Boltzmann Effects in Chain Branching and Pathway Branching for Diethyl Ether Oxidation

Low-temperature (LT) engine applications have several potential benefits, including reduced emissions and increased efficiency. Attaining these benefits requires accurate kinetic modeling of LT chain branching, which depends heavily on ketohydroperoxide (KHP) decomposition. For diethyl ether (DEE), a promising biofuel, current estimates of the KHP decomposition rate constant are largely based on empirical fits to data. In this study, we investigate the most important KHP isomer in DEE LT oxidation by applying variable reaction coordinate transition state theory to the main pathway for KHP decomposition: OO bond fission to produce •OH and a keto-alkoxy radical, •OQ'O. We also use ab initio kinetics methods to investigate the decomposition of •OQ'O, where we find dominant branching to acetic acid, with the remaining flux going to CH 3 C(O)OCHO. Additionally, new time-resolved measurements of DEE and acetic acid concentrations during LT (450–600 K) DEE oxidation are obtained in a laser photolysis flow reactor coupled with multiplexed photoionization mass spectrometry. These new experimental data, along with jet-stirred reactor data in the literature, are compared with the predictions of a recent DEE mechanism (Tran et al. Proc. Comb. Inst. 2019, 37, 511-519) that was modified with the newly calculated ab initio rate constants for KHP and •OQ'O decomposition. The predictions of the modified mechanism are quite poor when compared to the experimental data; this is primarily due to the new KHP ⇌ •OQ'O + •OH rate constant, which is 1–2 orders of magnitude slower than empirical values employed in recent mechanisms. To reconcile the new KHP rate constant and the experimental data, we explore and quantify the possible role of non-Boltzmann (nB) reaction sequences. The nB reactions have a substantial effect on both the overall mechanism reactivity and the •OQ'O branching to acetic acid. We also provide guidance on the proper implementation of nB reactions in kinetic mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1976

Heat transfer evaluation in a plasma core reactor

Numerical evaluations of heat transfer in a fissioning uranium plasma core reactor cavity, operating with seeded hydrogen propellant, was performed. A two-dimensional analysis is based on an assumed flow pattern and cavity wall heat exchange rate. Various iterative schemes were required by the nature of the radiative field and by the solid seed vaporization. Approximate formulations of the radiative heat flux are generally used, due to the complexity of the solution of a rigorously formulated problem. The present work analyzes the sensitivity of the results with respect to approximations of the radiative field, geometry, seed vaporization coefficients and flow pattern. The results present temperature, heat flux, density and optical depth distributions in the reactor cavity, acceptable simplifying assumptions, and iterative schemes. The present calculations, performed in cartesian and spherical coordinates, are applicable to any most general heat transfer problem.

Smith, D. E.↗
DOE OSTI2022

Layered CAD/CSG geometry for spatially complex radiation transport scenarios

Many spatially complex fission, fusion, and national security Monte Carlo (MC) radiation transport scenarios involve combining computer-aided design (CAD) models with constructive solid geometry (CSG) models. A layered geometry method has been implemented in the Shift MC code to address this need. With layered geometry, multiple CAD and/or CSG models can be clipped, translated, rotated, and placed in overlapping layers to form transport-ready geometries. Here, the utility of this method is demonstrated with two problems: (1) a fixed-source simulation with a layered geometry consisting of a LiDAR-generated CAD model of the Combined Arms Collective Training Facility urban environment overlaid with CSG models of a mock hotel and a detector apparatus, and (2) a k-eigenvalue calculation using a layered geometry model of the Transformational Challenge Reactor consisting of CAD fuel elements placed in a CSG core. Tallied particle flux distributions match expectations, but tracking robustness must be improved prior to general-purpose use.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

H 2 production in the 10 B(n,α) 7 Li reaction in water

Herein we demonstrate a method for measuring the H2 produced in water from the 10 B(n,α) 7 Li fission reaction. Low energy neutrons from the NIST Center for Neutron Research interact with borate-containing water in a temperature-controlled high pressure cell made from titanium. After exposure for one to several hours, the water is extracted and sparged with argon. H 2 entrained in the sparging gas is sampled with a small mass spectrometer. To determine the neutron exposure, a small amount of sodium is included in the borate solution. The water is collected and 24 Na activation is measured in a counting apparatus on the following day. The G-value for H 2 at room temperature is found to be (1.18 ± 0.10) molecules H 2 /100eV, in good agreement with previous estimates and recent modeling calculations.

10B↗
DOE OSTI2020

KRUSTY Experiment: Reactivity Insertion Accident Analysis

The centerpiece of the Kilopower Project, i.e., the Kilowatt Reactor Using Stirling TechnologY (KRUSTY) test, consists of the development and testing of a ground technology demonstration of a small fission power system based on a 1-kW(electric) space science power requirement. The KRUSTY test was authorized by the U.S. Department of Energy’s (DOE’s) National Nuclear Security Administration Nevada Field Office. Authorization was obtained by adding an amendment to the existing regulatory documents for the National Criticality Experiments Research Center to cover the KRUSTY experiment. This amendment was reviewed and approved by the DOE. The most important safety question for the experiment was the addition of over 2 $ of excess reactivity to the reactor system. This amount of excess reactivity meant that the analyst could postulate accidents where the reactor went prompt critical, leading to physical shock or melting of the fuel. This paper analyzes these accidents using computer calculations and examines the controls used to mitigate them. The estimation of the impacts both on accident progression and consequences of reactivity insertion events was a significant part of obtaining approval for the KRUSTY experiment. The regulatory approval of KRUSTY was one of the first to be obtained for a completely new reactor concept in many decades.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

A New Approach for Calculating the Alpha-Decay Half-Life for the Heavy and Super-heavy Elements and an Exact A Priori Result for Beyllium-8

A new general method for calculating the alpha decay half-life is presented. The method predicts an a priori exact value for the beryllium-8 half-life. Beryllium-8 is an exception to the current alpha decay theory captured in the Geiger-Nuttal law. The new method predicts the beryllium-8 alpha half-life using only constants and measured isotopic mass. The method also reliably predicts all the heavier isotope alpha decay half-lives consistent with the Geiger-Nuttal law. With respect to current theory, the inability of the Geiger-Nuttal Law to predict the alpha-decay half-life in the case of beryllium-8 has led to consideration of other decay mechanisms for this isotope, such as fission for example. One result is that given the consistency of the new method presented here for all isotopes including an exact a priori result for beryllium-8, the evidence strongly suggests that the beryllium-8 decay is in fact an alpha decay. A second result is that the method definitively demonstrates that the entire rest mass of the two helium-4 electrons is converted to energy in the decay process and this energy becomes part of the emitted alpha particle kinetic energy.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Application of the bRapid fission matrix burnup methodology to the JSI TRIGA Mark II research reactor

The novel burnup methodology, bRAPID, of the RAPID code system has been applied to the TRIGA Mark II research reactor at the 'Jozef Stefan' Institute. This paper presents the first step in our goal for on-line determination of isotope inventory and burnup, that the RAPID code system with its Fission Matrix (FM) method enables. Burnup analysis of changes in FM coefficients due to TRIGA fuel burnup at different irradiation times and specific power distribution in the reactor core is presented. Our results demonstrate that primary contribution for changes of FM coefficients is the burnup of the destination s, while the spectrum differences, that arise due to fuel's core location, have negligible effect on FM coefficients. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Thermodynamic assessment of lithium halide reciprocal salt systems for energy applications

Lithium halides are of interest because of their unique applications as fuel or coolants in molten salt reactors (MSRs) and electrolytes in thermally activated batteries. Because iodine is an important fission product for MSRs, the Molten Salt Thermal Properties Database-Thermochemical (MSTDB-TC) is being expanded to include relevant iodide reciprocal salt systems. In this work, we assess the thermochemical properties of lithium halide pseudobinary systems based on reported experimental phase diagrams, eutectic temperatures, and enthalpies of mixing. Within the framework of the modified quasi-chemical model in the quadruplet approximation (MQMQA), extrapolations to the related pseudoternary and pseudoquaternary systems have been performed. The pseudoquaternary representations of phase equilibria, enthalpy, entropy, and heat capacity are represented in the resultant models. Of interest for generating source terms for accident analysis is the ability to compute the vapor pressures of the dominant iodide species for the systems at any composition. Here, an example of such calculations is provided for the LiF-LiI system at x LiI = 0.01 and x LiI = 0.70, in this case indicating that the LiI vapor species are predominant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Comparative study between experimental measurements and model predictions of bubble rise velocity in molten LiCl-KCl

This research investigated the shape and rise velocities of gas bubbles (helium, nitrogen, argon, and krypton) in molten LiCl-KCl at 500 °C using a specialized optical visualization apparatus. The apparatus features a custom quartz rectangular prism cell housed in a furnace, enabling the application of the shadowgraph measurement technique. Bubbles with diameters between 0.45 and 3.2 mm were produced via a capillary, and their equivalent diameters and rise velocities were analyzed using high-resolution imaging and centroid tracking algorithms. Important non-dimensional parameters were calculated to characterize behavior. Bubbles with 0.03 < Eo < 1.32, 17 < Re < 718, and Mo order of 10 -11 were generated in the molten LiCl-KCl, revealing that bubbles with diameters less than 1.46 mm remain spherical, whereas those with larger diameters are ellipsoidal. The transition and oscillation onset for bubble shape occurred at larger diameters than those previously observed in air–water systems. Comparative analysis of the gas-molten salt data against existing air–water bubble rise velocity correlations highlighted inaccuracies in predictions, particularly across the different bubble regimes. An existing correlation was calibrated to the molten salt data to improve the bubble velocity predictions. In conclusion, this study contributes essential experimental data for the design and safety evaluation of molten salt reactors and offers insights for the optimization of sparge systems for fission product removal.

Bubble rise velocity↗
DOE OSTI2021

An Analytic Benchmark for Neutron Boltzmann Transport with Downscattering—Part I: Flux and Eigenvalue Solutions

Computing in the energy dimension is one of the greatest challenges confronting present-day deterministic neutron transport solvers. Accurately resolving the neutron flux as neutrons downscatter across resonances in the nuclear cross sections currently requires considerable computing power and suffers from approximation errors. Flux uncertainty resulting from the uncertainty of the resonance structure is the single-largest cause of reactivity uncertainty. Any additional reference solution for the critical neutron downscattering problem with resonance phenomena would be a boon to verification and validation of neutronics codes. This paper establishes a benchmark to verify the accuracy of neutron transport criticality solvers along the energy dimension. For the first time, the analytic solution of the flux amplitude is derived in the particular case of an infinite homogeneous medium with isotropic scattering in the center of mass and an arbitrary number of no-threshold, neutral particle reaction resonances (e.g., radiative capture, fission, and resonance scattering). Furthermore, original analytic expressions are established to quantify the discrepancy between the ψ k (E) and ψ α (E) flux amplitudes, respective solutions of the multiplication factor k, or the exponential time-evolution frequency α eigenproblems. The physical study of these relations led to analysis of their first-order relative difference near the criticality condition α=0. Finally, numerical solutions are provided to a benchmark problem constituted of the first resonance of 239 Pu, the 6.67-eV resonance of 238 U, and a scattering isotope with a flat cross section, allowing for the computational verification of the energy resolution of current neutron transport criticality codes. Through these novel results, this analytic benchmark can serve as a reference to verify the energy resolution and sensitivity analysis of neutron transport criticality calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Domestic and International Standards for Nuclear Criticality Safety – Overview and Status

The domestic and international consensus standards for nuclear criticality safety (NCS) have provided guidance for staff performing hands-on work in operations with fissionable materials. These consensus standards have contributed directly to the significant reduction in the rate of criticality accidents in process facilities since the 1940s. The last known criticality accident inside the United States was in 1978 (nearly 43 years ago) at the Idaho Chemical Processing Plant, and outside the United States, an accident occurred at Tokai-mura, Japan, in 1999 (23 years ago). The domestic consensus standards for NCS include the American Nuclear Society (ANS) standards. The ANS Standards Board, the NCS Consensus Committee, and the ANS-8 Subcommittee oversee the development and maintenance of these standards. There are currently eighteen standards in the ANS-8 series. Currently, there are six ANS-8 standards in revision mode and eleven in a maintenance mode with one new standard under development. The international consensus standards for NCS calculations, procedures, and practices are maintained and developed within the International Organization for Standardization (ISO), Technical Committee 85 on Nuclear Energy, Subcommittee 5 on Nuclear Fuel Technology, and Working Group 8, “Nuclear Criticality Safety.” Eleven standards are currently available, three standards are in revision mode, and two standards are development. This paper provides the NCS community with an overview and status report of domestic and international NCS consensus standards to stimulate interest and to support their continued development.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

X-ray absorption spectroscopy of trivalent Eu, Gd, Tb, and Dy chlorides and oxychlorides

Growing interest in next generation molten salt nuclear reactors necessitates a fundamental understanding of the chemical and physico-chemical properties of fission products, including elements of the lanthanide series. In this study, we probed coordination environment of Eu, Gd, Tb, and Dy in their chloride salts with both laboratory and synchrotron scale X-ray absorption spectroscopy (XAS). Here the aerobic melting was carried out with LnCl 3 ·nH 2 O, where Ln is Eu, Gd, Tb, and Dy, and yielded predominant formation of oxychlorides (LnOCl), that was confirmed by extended X-ray absorption fine structure (EXAFS) analysis and X-ray diffraction (XRD). Density functional theory (DFT) and time-dependent density functional theory (TDDFT) computations were performed to obtain ground state structures, lattice parameters, and to simulate L 3 -edge XANES (X-ray absorption near edge structure) spectra of the LnCl 3 ·nH 2 O and LnOCl. A blue shift, initially predicted by our calculations, was also observed for the experimental Ln L 3 -edge XANES spectra of the hydrated trichlorides vs oxychlorides, highlighting the power of combining the predictive ability of theory with experiment to elucidate the properties of the compounds of interest. The blue shift is indicative of structural stabilization of LnCl 3 ·nH 2 O due to higher symmetry compared with LnOCl. To the best of knowledge, this is a first study providing the framework for fundamental insights into the structure of Ln trichlorides and oxychlorides in support of future characterization of lanthanides in the molten salt environment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Domestic and International Standards for Nuclear Criticality Safety - Overview & Status [Abstract]

The domestic and international consensus standards for nuclear criticality safety (NCS) have provided guidance for staff performing hands-on work in operations with fissionable materials. These consensus standards have contributed directly to the significant reduction in the rate of criticality accidents in process facilities since the 1940s. The last known criticality accident inside the United States was in 1978 (nearly 43 years ago) at the Idaho Chemical Processing Plant, and outside the United States, an accident occurred at Tokai-mura, Japan, in 1999 (22 years ago). The domestic consensus standards for NCS include the American Nuclear Society (ANS) standards. The ANS Standards Board, the NCS Consensus Committee, and the ANS-8 Subcommittee oversee the development and maintenance of these standards. There are currently eighteen standards in the ANS-8 series. Currently, there are six ANS-8 standards in re- vision mode and eleven in a maintenance mode with one new standard under development. The international consensus standards for NCS calculations, procedures, and practices are maintained and developed within the International Organization for Standardization, Technical Committee 85 on Nuclear Energy, Subcommittee 5 on Nuclear Fuel Technology, and Working Group 8, “Nuclear Criticality Safety.” Eleven standards are currently available, three standards are in revision mode, and two standards are development. This paper provides the NCS community with an overview and status report of domestic and international NCS consensus standards to stimulate interest and to support their continued development.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2022

BISON validation of FeCrAl cladding mechanical failure during simulated reactivity-initiated accident conditions*

Here, a reactivity-initiated accident (RIA) is a postulated design basis accident in light water reactors (LWRs) in which a rapid reactivity insertion induces a fission rate increase and a fuel pellet temperature rise. During RIA, the fuel pellet thermally expands and may cause pellet-cladding mechanical interaction (PCMI). Separate effects PCMI tests were performed on C26M FeCrAl cladding tube samples, introducing biaxial stress via a well understood modified burst test (MBT) system. A high-speed camera in the MBT system captured the projections, covering a 360° view of the cladding deformation and enabling a digital image correlation (DIC) method to quantify the surface strains with high fidelity. Representative hot zero-power RIA mechanical loading conditions were applied to the sample, and the test duration ranged from 20 to 500 ms at an average temperature of 573 K. BISON finite element–based fuel performance code modeling was performed against the high-fidelity DIC data produced from the MBTs. Validation calculations were conducted with 2D models and systematically compared with test data of the burst time, burst pressure, burst hoop strains, and hoop strain rates. Based on the behaviors from the separate effects test, BISON calculations satisfactorily predicted the cladding deformation behaviors. Sensitivity analysis was conducted to identify highly influential mechanical properties responsible for the cladding failure behavior during the MBT experiments. The results highlight the significance of the cladding's mechanical strength in governing cladding strain, followed by the significance of the mechanical interactions between the pellet and the cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Subcritical Californium Source Drive Noise Analysis Measurements With Unreflected Uranium (93.15) Hydride

On July 6, 1989, subcritical californium source-driven noise analysis (CSDNA) measurements were performed with bare uranium hydride cylindrical assemblies at the Los Alamos National Laboratory Critical Experiments Facility in KIVA 1. Three configurations of ~7.5 cm radius-enriched (93.15 wt. % 235 U) uranium hydride with a density of ~10 g/cm 3 cylinders were assembled with uranium hydride heights of ~11, ~14, and ~16 cm. The uranium hydride (~2 and ~3 cm high) was in thin-welded stainless steel cans. The neutron multiplication factors obtained on-line from the measured ratios of spectral densities were 0.942 ± 0.002, 0.917 ± 0.003, and 0.867 ± 0.004 for uranium hydride heights of ~16, ~14, and ~11 cm., respectively. Neutron multiplication factors for these three configurations calculated by Los Alamos National Laboratory Monte Carlo methods were 0.952 ± 0.005, 0.922 ± 0.005, and 0.860 ± 0.005, all of which are in good agreement with the measurements. The break frequency noise analysis data could be fitted to obtain the prompt neutron decay constant at all subcritical states. Some data presented in this report are from notes that are not in the logbooks. A discussion of the two different point kinetics theories of the measurements that illustrated the deficiencies of the rigorous theory at low neutron multiplication factors (k eff < 0.80) is given in an Appendix. These measurements may be the only configuration of unreflected, highly enriched uranium hydride configuration assembled. The presence of hydrogen diluted the uranium density, decreasing the reactivity, but this is offset by slowing down the neutron energy to the energy ranges of higher-fission cross sections. Extrapolations to delayed critical indicate that adding one more 3 cm can and one more 2 cm can would result in the unreflected cylindrical system being slightly above critical with a critical mass of ~37 kg, which is considerably less than the critical mass of an unreflected and unmoderated highly enriched uranium metal sphere. An unreflected, highly enriched uranium hydride sphere would have a critical mass less than 37 kg.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Gamma Rays from Kilonova: A Potential Probe of r-process Nucleosynthesis

The mergers of compact binaries with at least one neutron star component are the potential leading sites of the production and ejection of r-process elements. Discoveries of galactic binary pulsars, short gamma-ray bursts, and gravitational-wave detections have all been constraining the rate of these events, while the gravitational wave plus broadband electromagnetic coverage of binary neutron star merger (GW170817) has also placed constraints on the properties (mass and composition) of the merger ejecta. But uncertainties and ambiguities in modeling the optical and infrared emission make it difficult to definitively measure the distribution of heavy isotopes in these mergers. In contrast, gamma rays emitted in the decay of these neutron-rich ejecta may provide a more direct measurement of the yields. We calculate the gamma production in remnants of neutron star mergers, considering two epochs: a kilonova epoch, lasting about two weeks, and a much later epoch of tens and hundreds of thousands of years after the merger. For the kilonova epoch, when the expanding ejecta is still only partially transparent to gamma radiation, we use 3D radiative transport simulations to produce the spectra. Futhermore, we show that the gamma-ray spectra associated with beta- and alpha-decay provide a fingerprint of the ejecta properties and, for a sufficiently nearby remnant, may be detectable, even for old remnants. We compare our gamma spectra with the potential detection limits of next generation detectors, including the Lunar Occultation Explorer (LOX), the All-sky Medium Energy Gamma-ray Observatory (AMEGO), and the Compton Spectrometer and Imager (COSI). We show that fission models can be discriminated via the presence of short-lived fission fragments in the remnant spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Advanced Porous Materials as Designer Platforms for Sequestering Radionuclide Pertechnetate

Technetium-99 ( 99 Tc), predominantly present as pertechnetate ( 99 TcO 4 – ), is a challenging contaminant in nuclear waste from artificial nuclear fission. The selective removal of 99 TcO 4 – from nuclear waste and contaminated groundwater is complex due to (i) the acidic and intricate nature of high-level liquid wastes; (ii) the highly alkaline environment in low-activity level tank wastes, such as those at Hanford, and in high-level wastes at locations like Savannah River; and (iii) the potential for 99 TcO 4 – to leak into groundwater, risking severe water pollution due to its high mobility. This Review focuses on recent developments in advanced porous materials, including metal–organic frameworks (MOFs), covalent organic frameworks (COFs), and their amorphous counterparts, porous organic polymers (POPs). These materials have demonstrated exceptional effectiveness in adsorbing 99 TcO 4 – and similar oxyanions. We comprehensively review the adsorption mechanisms of these anions with the adsorbents, employing macroscopic batch/column experiments, microscopic spectroscopic analyses, and theoretical calculations. In conclusion, we present our perspectives on potential future research directions, aiming to overcome current challenges and explore new opportunities in this area. Our goal is to encourage further research into the development of advanced porous materials for efficient 99 TcO 4 – management.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Preliminary Investigation of Target keff Search for Fuel Feed in MSRs

Fluid fuel molten salt reactor (MSR) concepts may achieve continuous-at-power operation with continuous or batch fresh fuel feed to counter reactivity loss from fuel burnup and fission product buildup. Accurate modeling and simulation of MSRs requires consideration of online fuel feed during depletion as well as a capability to determine the required feed amount to maintain a target multiplication factor (keff). As a first step supporting a full implementation of a new search capability into the SCALE code system’s TRITON reactor physics and depletion sequence, a preliminary investigation was performed to identify the best approach for a critical feed search capability considering computation time and accuracy. Using a simple iterative approach in which initial feed amounts were guessed based on previous depletion steps, the main drivers of computation time were found as the Monte Carlo neutron transport settings and the applied tolerance for the target keff. In the applied model, with tight convergence criteria (keff convergence to 8 pcm and a tolerance of ±15 pcm), the critical search calculation required twice the number of neutron transport calculations compared to a depletion calculation without critical search. In contrast, with a relaxed tolerance of ±50 pcm, the overall computation time for the critical search for merely increased by 30-40%. Further studies will be performed before implementing the critical search capability into TRITON.

Faulkner, Jonathon [ORNL]↗