Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear Graphite”

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

A method to predict texture effect on ion beam channeling analysis of polycrystals and the application to study the mosaic spreading effect in highly oriented pyrolytic graphite

We propose a method to convert the channeling Rutherford backscattering spectrum yield map of a single crystal to a polycrystal through a matrix rotation technique. The rotation matrix is determined by the deviation of the crystal axial direction from the original z axis. The final yield map is created after averaging the rotated yields using the texture function as the weight factor. For highly oriented pyrolytic graphite (HOPG) exhibiting mosaic spread, the method leads to a Gaussian kernel averaging of the map obtained from a single crystal. The yield map of a single crystal is obtained by a simulation of ion trajectories in a potential field described by Moliere screened Coulomb potentials. Yield maps are calculated under various σ values (standard deviations of mosaic spread). The simulated results are compared with experimental results obtained using 1.2 MeV alpha particle. σ is extracted through the best fitting, demonstrating that the method can be used to obtain texture details. The effects of mosaic spread on minimum yield χ min and the half-width at half maximum of angular scans ψ ½ are systematically modeled and compared with previous theoretical equations. The study also shows that previous theoretical equations are valid only at small σ values. The proposed method can be applied to any type of polycrystal and is not limited to HOPG. It provides near-surface mosaic spread and crystallography information with a longitudinal depth resolution of tens of nanometers and is not influenced by grain shapes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

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

Improved Cross Section Generation Capability of Griffin in FY22

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of deterministic steady-state transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication as well as transient solvers for spatial dynamics with the improved quasi-static method. Griffin uses cross section data in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. In recent years, the MC 2 -3 modules have been added to Griffin for fast reactor cross section generation, and the self-shielding application programming interface (SSAPI) was implemented in the ISOXML module for thermal reactor cross section generation. The on-the-fly slowing down method and double-heterogeneity treatment have been implemented to SSAPI and verified against particulate fuel-bearing graphite-moderated thermal reactor problems with high accuracy. This year, work has been focused on improving the cross section generation capability of ISOXML and streamlining the cross section generation procedures. In addition, the form function data were added to ISOXML in order to support the pin power reconstruction capability that was newly implemented in Griffin in this fiscal year. To facilitate the cross section generation using MC 2 -3 and SSAPI in Griffin, the cross section generation workflows have been set up for both fast and thermal spectrum reactors. The MOOSE action system tool was devised for fast spectrum problems, and the MOOSE stochastic tool was adopted to the branch calculation procedure for thermal spectrum problems. Meanwhile, to ensure the accuracy of group-constants, the thermal up-scattering kernel calculator accounting for resonance scattering was implemented in ISOXML, demonstrating the accurate computation of a Doppler-broadened scattering kernel of any Legendre order within a reasonable timescale. Other aspects of ISOXML, such as deletion solver and data, documentation, ISOXML file management, and the interface for Mixture, were improved as well.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Preparation of a uranium monocarbide anode and electrochemical characterization in molten LiCl-KCl-UCl 3

Porous uranium carbide (UC) pellets possessing moderate electrical conductivity were synthesized by reaction of UO 2 with graphite at temperatures up to 1550°C under rough vacuum. Conversions as high as 98% were achieved at soak times of 2-4 hours. The electrochemistry of the UC pellets in molten LiCl-KCl-6.5 wt% UCl 3 was explored using a variety of techniques including DC polarization methods, cyclic voltammetry, chronopotentiometry and bulk electrolysis. Here, the electrode reaction for anodic dissolution was found to be kinetically controlled by dissociation of UC to a transition state complex that was hypothesized to consist of a uranium atom partially complexed by chloride ions. Precise measurements of current efficiencies using chronopotentiometry indicated upper limits of 90.9 ± 3.4% and 98.3 +1.7/-3.7% for anode and cathode, respectively, when operating at anodic overpotentials near +300 mV. Bulk electrolysis of a UC pellet performed by passing 98% of the theoretical charge resulted in nearly complete recovery of its uranium content as highly pure metal at the cathode.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On the product phases and the reaction kinetics of carbothermic reduction of UO 2 +C at relatively low temperatures

The synthesis of UC using carbothermic reduction of UO 2 and C mixtures has been well studied at high temperatures. However, the product phase behavior of carbothermic reduction at low temperatures (≤1773 K) is not well studied. Such a study is important as low temperatures permit single phase UC synthesis without forming secondary higher carbides, and it further supports the knowledge base of the process that needs to be used for transuranic elements such as plutonium that have high vapor pressures at elevated temperatures. Therefore, a low temperature carbothermic reduction of two different C/UO 2 molar ratios under inert and reducing environments have been studied here. Two different sample holding crucibles, alumina (Al 2 O 3 ) and graphite, were also used here to differentiate the hypostoichiometric (UC 1-a ) and oxygen dissolved (UC 1-x O x ) uranium monocarbide phases adding more details on the two systems. Also, the reaction kinetics involved in the formation of UC via the carbothermic reduction of UO 2 +C using product phases instead of evolved gases such as carbon monoxide is reported here. Under inert atmospheres but with significant oxygen partial pressures, the low temperature carbothermic reduction of UO 2 +C produced up to 90 wt.% UC 1-x O x type oxycarbides as was confirmed by Xray powder diffraction. Reducing Ar-4%H 2 environments at these temperatures were not successful in synthesizing UC as it reduces the amount of C required for the carbothermic reduction, leaving UC phase at a non-equilibrium state. Inert atmospheres with low or negligible oxygen partial pressures on the other hand produced near stoichiometric UC at high phase purity, especially at 1673 – 1773 K temperature range. An activation energy of 377±75 kJmol -1 was also calculated using product phase concentrations of the carbothermic reduction of UO 2 +C under these inert Ar (g) atmospheres.

36 MATERIALS SCIENCE↗
DOE OSTI2022

FY22 Status Report on the ART-GCR CMVB and CNWG International Collaborations

This work presents the numerical model for the high temperature test reactor (HTTR) loss of forced cooling (LOFC) experiment with the INL codes Griffin, BISON, and RELAP-7, based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Promising results were obtained, with the overall behavior of the reactor successfully captured. Changes in the heat transfer coupling, as compared to the fiscal year (FY)-21 model, enabled drastic improvement of the steady-state solution, both in terms of computational time and global energy discrepancies. The former was reduced by a factor of roughly 60, whereas the latter decreased from 11% to less than 2%. Furthermore, the discrepancy in the steady-state multiplication factor was improved from +2,300 and +2,900 pcm (for the 30 and 9 MW cases) to -700 and +1,200 pcm, respectively, and now falls well within the large measurement uncertainties stemming from graphite impurities. Validation of the Monte Carlo model used to generate cross sections was also performed against available measurements. Though significant, the discrepancies remain acceptable overall in light of the large uncertainty. Numerous improvements are still needed to better compare with the experiment involving the 9 MW case and to instill greater confidence in the model’s ability to accurately predict the 30 MW behavior. Specifically, the power levels predicted by the 9 MW transient simulation following re-criticality remain low, pointing to an underestimation of the passive cooling of the core. A key aspect of future work will be to better understand the flow pattern during the LOFC event, particularly to determine if natural or forced convection is occurring inside the reactor pressure vessel (RPV). More generally, additional validation data would be immensely useful for further enhancing the numerical model and better matching the experiments. In addition, a more sophisticated thermal-hydraulics model that simulates all the channels as a single system model should be considered to take into account the rest of the primary loop. Finally, even if the results are in better agreement with the experiments, sensitivity analysis and uncertainty quantification will be necessary to evaluate the model uncertainty.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

AGR-3/4 TRISO Fuel Compact Ceramography

The combined third and fourth irradiation in the Advanced Gas Reactor (AGR) program (AGR-3/4) contained tristructural isotropic (TRISO)-coated particle fuel and designed-to-fail (DTF) fuel particles. The DTF particles were only coated with highly anisotropic pyrocarbon (PyC) so they would purposely fail during the AGR-3/4 irradiation and provide a source of fission products for measurement. To observe the post-irradiation morphology of these DTF particles and the TRISO-coated particles, three AGR-3/4 compacts were mounted in epoxy, sectioned above their centerlines, and ground/polished. Three rounds of grinding/polishing and optical microscopy were performed so that the particles could be observed at multiple planes. Each compact contained approximately 1,872 TRISO particles and exactly 20 DTF particles. The three AGR-3/4 compacts examined covered a wide temperature range and featured both the hottest average irradiation temperature (1375°C) and the coldest/lowest burnup (872°C and 5.5% fissions per initial metal atom) of any compacts to undergo post-irradiation examination (PIE) in the AGR program to date. A total of 29 DTF particles were located and observed via microscopy across the three compacts. All the DTF particles observed via microscopy had completely failed PyC coatings. Some different DTF kernel and PyC morphologies were observed that have been attributed to the differences in irradiation temperature and/or burnup. At low and medium irradiation temperatures (roughly 850 to 1050°C), it appears irradiation-induced dimensional changes in the DTF PyC caused it to fracture and fold in on itself, and the kernel deformed to accommodate this PyC deformation. The extent of DTF kernel deformation to accommodate DTF PyC buckling tended to be greater for the medium-temperature fuel compared to the cold fuel. In the high-irradiation-temperature compact (1375°C), the DTF PyC layer appears to have completely reacted chemically with the kernel material. The only material surrounding the DTF particles in this hot compact resembles the compact matrix material. Observations of the TRISO-coated driver fuel particles and fuel compact graphitic matrix were also made. The TRISO particle kernels in the medium-temperature compact (1047°C), and many in the high-temperature compact (1375°C), showed morphologies consistent with what was seen in AGR-1 and AGR-2. In other TRISO particles from the high-temperature compact, a spatial gradient in the kernel appearance was evident. Here the cool side of the kernel (away from the center of the compact) had a much darker appearance (like that of the buffer and pyrocarbon). This is the first observation of this kind of kernel spatial gradient in AGR fuel. It is believed that the high irradiation temperature of that compact (an average of 1375°C) activates or accelerates an unidentified chemical reaction, and the temperature gradient within the fuel causes a sharp spatial gradient. The low-temperature/low-burnup compact also showed unique kernel morphologies where the TRISO-coated particles had clearly distinguishable oxide and carbide phases at the center of the kernel, an oxide rind surrounding this, and remnants of a carbide skin surrounding the oxide rind. These are features observed in the as-fabricated fuel that are no longer present in higher burnup fuel. Occasional gaps between the outer pyrolytic carbon (OPyC) and graphitic matrix material were observed in all compacts. Gaps were most often found in the small, matrix-filled spaces between adjacent particles. Finally, buffer morphologies in AGR-3/4 TRISO particles were like those observed in AGR-1 and AGR-2, and the buffer fracture frequencies in AGR-3/4 and AGR-2 TRISO particles were plotted versus fast neutron fluence and irradiation temperature. In AGR-3/4 the buffer fracture rate was highest (23%) for an irradiation temperature of 1047°C and a fast fluence of 5.18E25 n/m 2 . Increasing the irradiation temperature to 1375°C reduced the buffer fracture rate to 14%. Reducing the irradiation temperature to 872°C and the fast fluence by a factor of 3 reduced the buffer fracture rate to 7%. The temperature and fluence dependencies observed for AGR-3/4 buffer fracture were consistent with those observed in AGR-2. Lower fluences and/or higher temperatures significantly reduced the buffer fracture frequencies. Lower buffer fracture frequencies were also observed at lower temperatures as long as the fluence was also significantly reduced. This is believed to be due to lower fluence resulting in less irradiation-induced dimensional change (shrinkage) of the buffer, and a higher temperature promoting enhanced creep relaxation of stresses within the buffer, leading to less buffer fracture.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Hydrogen Burning of 29 Si and Its Impact on Presolar Stardust Grains from Classical Novae

Abstract Presolar stardust grains found in primitive meteorites are believed to retain the isotopic composition of stellar outflows at the time of grain condensation. Therefore, laboratory measurements of their isotopic ratios represent sensitive probes for investigating open questions related to stellar evolution, stellar explosions, nucleosynthesis, mixing mechanisms, dust formation, and galactic chemical evolution. For a few selected presolar grains, classical novae have been discussed as a potential source. For SiC, silicate, and graphite presolar grains, the association is based on the observation of small N ( 12 C)/ N ( 13 C) and N ( 14 N)/ N ( 15 N) number abundance ratios compared to solar values, and abundance excesses in 30 Si relative to 29 Si, as previously predicted by models of classical novae. We report on a direct measurement of the 29 Si(p, γ ) 30 P reaction, which strongly impacts simulated δ 29 Si values from classical novae. Our new experimental 29 Si(p, γ ) 30 P thermonuclear reaction rate differs from previous results by up to 50% in the classical nova temperature range ( T = 100–400 MK), while the rate uncertainty is reduced by up to a factor of 3. Using our new reaction rate in Monte Carlo reaction network and hydrodynamic simulations of classical novae, we estimate δ 29 Si values with much reduced uncertainties. Our results establish δ 29 Si values measured in presolar grains as a sensitive probe for assessing their classical nova paternity. We also demonstrate that δ 30 Si values from nova simulations are currently not a useful diagnostic tool unless the large uncertainty of the 30 P(p, γ ) 31 S reaction rate can be significantly reduced.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Design of a prototypical natural circulation water-based reactor cavity cooling system (RCCS) for a pebble-bed generic FHR

A prototypical natural circulation water-based reactor cavity cooling system (RCCS) for the UC Berkeley Mark-1 pebble-bed gFHR is designed based on one-dimensional thermal hydraulics modeling and optimization implemented in an in-house MATLAB code. The model employs a lumped core consisting of fuel pebble and graphite pebble regions, while the graphite reflector, vessel, and RCCS are represented as separate but energy-coupled regions. The model is derived based on steady state energy balance equations accounting for conductive heat transfer from the lumped core to the vessel and subsequent radiative heat transfer to the RCCS and convective heat transfer in the water. Further, mass flow rate is calculated based on momentum and thermal energy balance in the RCCS. Starting with ANL’s water NSTF as baseline, effects of RCCS design parameters such as size of the plate between pipes, pipe diameter, source-to-sink distance, surface emissivity, and pitch from the core are examined in an integral effects framework which accounts for subsequent variation of the mass flow rate, number of pipes, convection coefficient, mean region and surface temperatures, radiative view factors, and water outlet temperature. The importance of including a prototypical reactor design in RCCS design calculations is emphasized as the calculations show that temperature drops from the core to the RCCS fluid to transfer the heat tend to be the limitation not the heat removal capacity of the fluid. A closed-loop design is obtained with capability to safely remove up to 0.72% of the nominal reactor power at maximum estimated peak conditions and 0.40% at shutdown. The design is based on physics calculations and does not account for economic optimization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

2023 Critical Materials Strategy

The global effort to curb carbon emissions is accelerating demand for clean energy technologies and the materials they rely on. Demand for these materials will only continue to grow, especially as some nations aim to achieve net zero emissions by 2050. While some major materials like steel, copper, and aluminum are already powering the fossil fuel economy, others are more minor materials with potential supply risks. These risks could jeopardize the ability to reduce greenhouse gas emissions within the desirable timeframe to avoid significant climate change. In some cases, it may be necessary to take action to improve the resilience of material supply chains and mitigate supply risks. Understanding the importance of individual materials to clean energy and the supply risks associated with them is necessary to identify which materials may serve as potential roadblocks to a clean energy future. The U.S. Department of Energy (DOE) issued a series of 13 supply chain deep dive assessment reports on various energy technologies in 2022 in response to President Biden’s Executive Order on America’s Supply Chains (E.O. 14017). These reports emphasized that supply chain bottlenecks can occur at any stage of the value chain from mining and refining to component and even sub-system manufacturing. The bottlenecks are a combination of factors such as material availability, equipment availability, work force availability and quality, logistics, regulatory framework, and market conditions. These bottlenecks were worsened during the global Covid-19 pandemic. Its lingering impacts have hindered capacity expansion for material supply chains and prevented product lead-time recovery. One approach to reduce supply chain risks for the United States is to have a strong domestic manufacturing sector with a diverse set of producers. Boosting responsible domestic production would require leveraging the latest science not only in material extraction but also in developing substitutes, recycling, reuse, and remanufacturing. This report is an updated analysis of previous Critical Materials Strategy (CMS) reports published by the DOE in 2010, 2011, and 2019 based on national and global priorities, technology advancement, and technology adoption trends. Like the CMS reports, this analysis presents the results of a formal material criticality assessment to identify which materials are critical to the continued deployment of clean energy technologies globally. The analysis in this report leveraged the DOE supply chain deep dive assessments to develop the initial list of materials to evaluate. This DOE Critical Materials Assessment (CMA) is conducted independently of criticality assessments performed by other U.S. government agencies, such as that conducted by the U.S. Geological Survey (USGS). This analysis complements the USGS critical minerals determination in three aspects. First, the DOE assessment is performed from a global perspective, while the USGS analysis focusses on the importance of minerals to the U.S. economy. Second, this report focuses on the importance of materials to clean energy technologies, rather than to the economy in general. Lastly, this study is forward looking to 2035 based on clean energy deployment scenarios, whereas the USGS assessment is retrospective. Materials evaluated in this report that do not appear in the USGS Critical Minerals List include copper, uranium, electrical steel, and SiC. A draft version of this report received ~80 public comments related to supporting data and methodological improvement. Those comments have been incorporated as much as possible where appropriate. Highlights of findings from this 2023 CMA include: Rare earth materials (neodymium, praseodymium, dysprosium, and terbium) used in magnets in electric vehicle (EV) motors and wind turbine generators continue to be critical. While dysprosium (Dy) and terbium (Tb) are both heavy rare earth elements that serve the same function in magnets, the criticality of Tb is slightly lower than that for Dy in the short term due to the widespread use of Dy in high-grade magnets and Tb’s present role as a substitute. Similarly, praseodymium (Pr) is critical in the medium term but only near critical in the short term because it is more substitutable in magnets than neodymium (Nd); Materials used in batteries for EVs and stationary storage are now considered to be critical. While cobalt (Co) was found to be critical in this and previous reports, lithium (Li) becomes critical in the medium term due to its broader use in various battery chemistries and the rampant growth of the EV industry. Natural graphite is a new addition in this assessment and is also found to be critical; Platinum group metals used in hydrogen electrolyzers, such as platinum (Pr) and iridium (Ir), are critical due to an increased focus on hydrogen technologies to achieve net zero carbon emissions, while those used in catalytic converters, such as rhodium (Rh) and palladium (Pd), were screened out due to the decreased importance of catalytic converters in the medium term; Gallium (Ga) continues to be critical due to its use in light-emitting diodes (LEDs). In addition, the use of Ga has increased in magnet manufacturing and in semiconductor in forms such as gallium arsenide (GaAs) or gallium nitride (GaN); Major materials like Aluminum (Al), copper (Cu), nickel (Ni), and silicon (Si) move from noncritical in the short term to near critical in the medium term due to their importance in electrification; Electrical steel is near critical due to its use in transformers for the grid and electric motors in EVs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Investigating the AGN-201M research reactor's unique dominance ratio

The AGN-201M reactor is a small, thermal spectrum reactor located at the University of New Mexico (UNM). It is moderated by polyethylene, reflected by graphite, with fuel comprised of uranium microspheres embedded in polyethylene plates that are separated by an aluminum baffle. The dominance ratio of a multiplying system is the ratio of the first higher-mode eigenvalue to the fundamental, k{sub 1}/k{sub 0}. It is a key parameter in characterizing the rate of convergence, analyzing system stability, and determining how tightly coupled the system's neutron fission is. The purpose of this paper is to investigate the system's unusually low dominance ratio of 0.632. This investigation varied the reactor's physical parameters, such as the fuel geometry, fuel density, and reflector thickness, and examined their impact on the reactor's dominance ratio. Additionally, this study examined neutronically similar systems in hopes of identifying common causes for systems with low dominance ratios. The reason why the dominance ratio of the reactor is small compared to larger thermal reactors was determined to be because of its significantly smaller geometry and the composition of the fuel plates. The reflector's effect on the dominance ratio is small in comparison to the other factors, but was found to have a non-zero effect. Furthermore, the AGN-201M was found to have a significantly lower dominance ratio than systems in which it shares a very high (c{sub k} >95%) degree of neutronic similarity.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1997

Extraterrestrial Helium (He@C60) Trapped in Fullerenes in the Sudbury Impact Structure

Fullerenes (C60 and C70) have recently been identified in a shock-produced breccia (Onaping Formation) associated with the 1.85-Ga Sudbury Impact Crater. The presence of parts-per-million levels of fullerenes in this impact structure raises interesting questions about the processes that led to the formation of fullerenes and the potential for delivery of intact organic material to the Earth by a large bolide (e.g., asteroid or comet). Two possible scenarios for the presence of fullerenes in the Sudbury impact deposits are that (1) fullerenes are synthesized within the impact plume from the C contained in the bolide; or (2) fullerenes are already present in the bolide and survived the impact event. The correlation of C and trapped noble gas atoms in meteorites is well established. Primitive meteorites contain several trapped noble gas components that have anomalous isotopic compositions, some of which may have a presolar origin. Several C-bearing phases, including SiC, graphite, and diamond, have been recognized as carriers of trapped noble gases. It has also been suggested that fullerenes (C60 and C70) might be a carrier of noble gas components in carbonaceous chondrites. Recently, fullerenes have been detected in separate samples in the Allende meteorite. Carbon-60 is large enough to enclose the noble gases He, Ne, Ar, Kr, and Xe, but it is too small to contain diatomic gases such as N2 or triatomic gases such as CO2. Recent experimental work has demonstrated that noble gases of a specific isotopic composition can be introduced into synthetic fullerenes at high temperatures and pressures; these encapsulated gases can then be released by the breaking of one or more C bonds during step-heating under vacuum. These thermal-release patterns for He encapsulated within the C60 molecule (He@C60) are similar to the patterns for acid residues of carbonaceous chondrites, suggesting that fullerenes could be an additional carrier of trapped noble gases in acid residues of meteorites. Analysis and Results: In order to characterize the noble gas compositions of the Sudbury fullerenes, we undertook a systematic study of acid-resistant residues throughout the C-rich layer (Black member) of the Onaping Formation. Samples were demineralized and extracted using standard techniques. The Onaping extracts were analyzed using several techniques, including UV-Vis adsorption, electro spray mass spectrometry, and laser desorption (linear and reflectron) time-of-flight (TOF) mass spectrometry (LDMS). The Sudbury fullerenes were then separated and purified using HPLC coupled with a photo diode array detector. The HPLC extracts containing the purified fullerenes were loaded into a metal tube furnace within a glove box under a N atmosphere in preparation for noble gas analyses. The 3-He and 4-He content of the fullerene extracts was measured using previously reported standard techniques . Discussion: Fullerenes (C60 and C70) in the Sudbury Impact Structure have been found to contain trapped He with a 3-He/4-He ratio greater than 5 x 10(exp -4). The 3-He/4-He ratio exceeds the accepted solar value by more than 30% and is more than 10x higher than the maximum reported mantle value. Terrestrial nuclear reactions or cosmic-my bombardment are not sufficient to generate such a high ratio. The 3-He/4-He ratios in the Sudbury fullerenes are similar to those determined for interplanetary dust particles. The greater-than-solar ratios of 3-He/4-He in the Sudbury fullerenes may indicate a presolar origin, although alternative mechanisms occurring in the ISM to explain these high ratios (e.g., spallation reactions, selective He implantation, etc.) cannot be entirely ruled out. We are currently attempting to isolate enough fullerene material to measure anomalous Ne (or Kr or Xe) contained within the C60 (e.g., the "pure" 22-Ne component) and thus determine whether the Sudbury fullerenes are indeed presolar in origin.

Becker, L.↗
DOE OSTI2026

The stellar origins of 96 Zr excesses in presolar graphites from the Murchison meteorite

Context. Zirconium-96 is a stable isotope that can be synthesized under different neutron-rich nucleosynthetic conditions. Astrophysical models predict its production to occur in various stellar environments: from low-to-intermediate-mass asymptotic giant branch (AGB) stars to massive stars and core-collapse supernovae. Aims. Detections of 96Zr excesses, in combination with other isotopic measurements from presolar grains can provide unique constraints on its stellar origin. Presolar grains are microscopic particles found in primitive Solar System materials, which formed in stellar winds and supernova ejecta. The isotopic composition of each grain can provide us a snapshot of the nucleosynthetic processes that took place during the parent star’s lifetime. Methods. In this study, we measured the stable isotopes of C, N, O, Mo, Zr, and Ru in high-density presolar graphite grains from the Murchison meteorite and found four grains that contain positive isotopic anomalies in 96 Zr carried by their internal subgrains. We analyzed multi-element isotopic datasets from each grain to explore the source of the observed 96 Zr excesses. Results. Comparisons with stellar models indicate that two grains likely condensed in an intermediate-mass AGB star with initial metallicity of Z ≤ 0.014. Their 96 Zr/ 94 Zr ratios also match those predicted for born-again AGB stars undergoing a very late thermal pulse and rapidly accreting white dwarfs. After comparing the relative populations of the aforementioned dust-producing stars, we propose rapidly accreting white dwarfs as a new, and more likely, stellar source for one of the presolar grains. The remaining two grains could have originated in the supernova ejecta of massive stars, due to correlated excesses in the p-nuclides, 92, 94 Mo. Thus, grains with 96 Zr anomalies can have a variety of stellar origins, in agreement with theoretical studies. Conclusions. Our study highlights the importance of multi-element analysis in constraining the types of stars where presolar grains have condensed. These data will help improve our understanding of various nucleosynthesis processes in different stellar phases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Development of the TREAT M2 experiment as a transient benchmark

A transient benchmark based on the Transient Reactor Test Facility (TREAT) M2 Calibration experiment (M2-CAL) is under development. TREAT, at Idaho National Laboratory, is a graphite moderated air-cooled research reactor which has been used extensively for fuel material testing under extreme and accident conditions. Accurate benchmark models are a beneficial component in the operations, experimental planning, and development of TREAT. In this work, we present a transient benchmark model for the M2-CAL experiment core loading. The benchmark model incorporates the coupling between the Monte Carlo Code SERPENT and the computational fluid dynamics code OpenFOAM to capture the temperature feedback mechanism. The M2-CAL transient 2580 was simulated in this work. A pre-transient analysis was performed to determine the optimum core composition and conditions before the beginning of the transient. The analysis was validated against historic TREAT kinetic measurements and the worth of the transient rod T-2. The axial power distribution in the flux wire for the M2-CAL experiment was determined and contrasted with the experiment. The model was then used to simulate the M2-CAL transient 2580 experiment based on the reported pre-transient and transient conditions. Several transient observables were calculated and compared to the experiment. The model shows good agreement with the experimental power traces, period, and average increase of power at the power ramp (time > 7.8 s). Inverse point kinetics analysis was performed during the period of the power ramp for the model and the experiment based on that, the temperature feedback component was isolated and contrasted with the experimental feedback. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Thermomechanical Processing of Uranium Alloys with 10 and 50 Weight Percent Zirconium

This work focuses on casting and thermomechanical processing for moderate sized (500-1000g) 10% by weight and 50% by weight uranium-zirconium alloys. Three castings of uranium – 10wt% zirconium and a fourth casting of uranium – 50wt% zirconium were completed using different castings processes. The first casting process relied on only induction mixing within the VIM system thereby eliminating the need for an arc-melting feedstock preparation step. The second casting process used the arc-melting process prior for feedstock to the same VIM system. In this study, the induction mixing process alone with a graphite crucible appears insufficient for the fabrication of UZr alloys. Graphite is suitable with the use of master alloy UZr feedstock produced by arc-melting. There was a moderate increase of approximately 100 ppm of carbon within the U10Zr alloys. However, there was a near 1000 ppm increase of carbon in the U50Zr alloy. This carbon content is considered too high for forming processes, and therefore, warrants the use of a ceramic crucible for UZr castings near 50wt% zirconium. Zirconia does not appear suitable for UZr castings near 50wt% zirconium. The use of Beryllia or Thoria may be a more suitable option. However, due to the significantly lower liquidus and lower zirconium content, Zirconia crucibles do appear viable for carbon reduction in the U10Zr system. Two thermomechanical processes are also studied in this work - hot extrusion of the U50Zr system and heat treatment in an argon atmosphere. This work demonstrated that extrusion of the U50Zr alloy near 680°C is a viable choice even with +1000 ppm carbon. Extrusion of the U50Zr system near 630°C was unsuccessful and does not appear viable unless significant force and tooling strength is available. There appears to be a moderate increase in extrudate hardness when compared to as-cast. Heat-treatment reduced the hardness of as-cast alloys in the U10Zr castings as well as resolve some variance in the microstructures. However, for the U50Zr billet that did not successfully extrude at 630°C, heat treatment at 900°C for 96 hours did not have a significant impact on hardness or microstructure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Molten Salt Reactor Experiment Simulation using Shift/Griffin

The Department of Energy (DOE)’s NEAMS focuses its efforts on the development of advanced modeling and simulation (M&S) tools for light-water reactors (LWRs) and non–LWRs (i.e., molten salt reactors, high-temperature gas reactors, microreactors, and fast reactors). In the previous fiscal year, the Multiphysics Applications Driver Technical Area funded molten salt reactor (MSR) M&S at Oak Ridge National Laboratory (ORNL) to generate multigroup macroscopic cross sections with Shift for a MSRE 2D lattice model in Griffin. In addition, Shift’s capability to calculate gamma dose rates from activated components in the primary exchangers in a molten salt breeder reactor was also demonstrated. In fiscal year 2023, ORNL generated multigroup macroscopic cross sections using Shift for a 3D MSRE core model. MSRE depletion calculations using Griffin were also demonstrated in this fiscal year. For the depletion calculation, one-group microscopic cross sections for the 3D MSRE core were generated using Shift, and the decay transmutation library from ORIGEN was converted to an ISOXML file, which is required as input in Griffin. Several Monte Carlo codes, such as OpenMC and Serpent, were also used to benchmark and supplement multigroup cross sections generated by Shift. Multigroup libraries were generated with 8 and 20 group structures, and the study found the 8-group structure to be more accurate when comparing Griffin results to continuous energy (CE) Monte Carlo results. The average flux from CE Shift calculations is up to ~6% higher than the CE Serpent calculations because of different values applied for the energy released per fission (κ values). The average flux in the fuel salt calculated by Griffin using cross sections generated with Shift agrees well with the reference CE Shift solution; the same is valid for the corresponding Serpent results. The maximum relative error is ~6% and ~2% compared to the CE Shift and Serpent reference solutions, respectively. Meanwhile, the average flux calculated by Griffin in the graphite moderator shows a higher difference in the thermal range when compared to both reference Monte Carlo solutions; this result suggests a need for improvement in cross section generation for the graphite moderator in the thermal range in both Monte Carlo codes. Griffin depletion calculations using cross sections from Shift and Serpent were performed and compared against ORIGEN calculations, and the nuclide densities calculated by Griffin were found to be generally in agreement with those of ORIGEN. Because a different approach was taken to calculate the energy released per fission ( κ values) in Shift and Serpent, a difference in nuclide densities from differences in the average flux was observed between Griffin using Serpent and Shift cross sections. Griffin calculations with Shift cross sections produced higher average flux in the salt than with Serpent cross sections, leading to higher consumption of 235 U and higher production of 135 Xe. For time-dependent depletion calculations, cross sections were generated with Serpent, and Griffin’s results using these cross sections were compared to CE Serpent depletion results, demonstrating good agreement. The average difference in keff between Serpent and Griffin as a function of burnup is about 155 pcm. Similarly, good agreement with small differences up to ~0.5% was also noticed in the nuclide density of 235 U and 135 Xe. More details regarding the methodologies invoked to generate the cross section to make code-to-code comparisons are discussed further in this report. User feedback on Griffin and Shift capabilities that will enhance these calculations is provided in this report for future consideration. The work performed this fiscal year can be extended further for multiphysics coupling of Griffin-Pronghorn/SAM with Mole to study precursor flow and salt chemistry.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Experimental validation of a high fidelity Monte Carlo neutron transport model of the MIT graphite exponential pile

High-fidelity modeling and simulation were performed for the MIT graphite exponential pile (MGEP) using Monte Carlo neutron transport codes OpenMC and MCNP, and the results were validated by experimental data. The MGEP is being used as the test bed for the design of an autonomous control system for the pile's neutron flux distribution. The main contribution of this work is to generate the training data sets of neutron flux distributions with different locations of control rods that perturb the neutron flux profiles. First, code -to-code cross verification between OpenMC and MCNP was performed to ensure consistency of the numerical modeling within statistical uncertainties. To validate the accuracy of this high-fidelity model, a series of neutron flux measurements were conducted using a Helium-3 (He-3) neutron detector on a mobile platform that is placed inside the pile. Second, the neutron flux profiles were measured in four vertical layers of interest, and compared to the corresponding simulation results. The comparison results shows that the root mean square error is less than 2.5% in the two upper layers, and less than 4.5% in all four measured layers. Here the results validated the accuracy of the modeling and simulation. Finally, the relative change of the neutron flux profiles from moving control rods was analyzed, which identified the layer that has the best sensitivity regarding the control rods movements. Thus, this work identified and provided training data sets of both simulated and experimental neutron flux profiles in the most sensitive layer, paving the path forward to the real-time experimental demonstration of the autonomous control system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗