Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutron Capture”

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 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 OSTI2021

Tuning spin excitations in magnetic films by confinement

Spin excitations of magnetic thin films are the founding element for novel transport concepts in spintronics, magnonics, and magnetic devices in general. While spin dynamics have been extensively studied in bulk materials, their behaviour in mesoscopic films is less known due to experimental limitations. Here, we employ Resonant Inelastic X-Ray Scattering to investigate the spin excitation spectrum in mesoscopic Fe films, from bulk-like down to 3 unit cells thick. In bulk-like samples, we find isotropic, dispersive ferromagnons consistent with the dispersion observed by neutron scattering in bulk single crystals. As the thickness is reduced, these ferromagnons survive and evolve anisotropically: renormalising to lower energies along the out-of-plane direction while retaining their dispersion in the in-plane direction. This thickness dependence is captured by simple Heisenberg model calculations accounting for the confinement in the out-of-plane direction through the loss of Fe bonds. Overall, our findings highlight the effects of mesoscopic scaling on spin dynamics and identify thickness as a knob for fine-tuning and controlling magnetic properties in films.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Dialing in Direct Air Capture of CO 2 by Crystal Engineering of Bis-iminoguanidines

Direct air capture (DAC) technologies that extract carbon dioxide from the atmosphere via chemical processes have the potential to restore the atmospheric CO 2 concentration to an optimal level. This study elucidates structure-property relationships in DAC by crystallization of bis-iminoguanidine (BIG) carbonate salts. Here, their crystal structures are analyzed by X-ray and neutron diffraction to accurately measure key structural parameters including molecular conformations, hydrogen bonding, and π-stacking. Experimental measurements of key properties, such as aqueous solubilities and regeneration energies and temperatures, are complemented by first-principles calculations of lattice and hydration free energies, as well as free energies of reactions with CO 2 , and BIG regenerations. Minor structural modifications in the molecular structure of the BIGs are found to result in major changes in the crystal structures and the aqueous solubilities within the series, leading to enhanced DAC.

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

Gravitational waves from binary neutron star mergers with a spectral equation of state

In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

First-Principles Studies of Tritium Species Dissociability & Diffusivity Across the Interface of Nickel-Plated Zircaloy-4

Zirconium (Zr) and its alloys (Zircaloy-4) are widely used in nuclear reactors due to their low neutron adsorption cross-section and excellent corrosion resistance. In tritium-producing burnable absorber rods (TPBARs), the metal getter tube located between the cladding and the γ-LiAlO 2 pellets is composed of nickel (Ni)-plated Zircaloy-4, which is used to capture tritium ( 3 H) species (mainly 3 H 2 and 3 H 2 O) generated from γ-LiAlO 2 pellets during irradiation. The 3 H-related products transfer to the surface of metal Ni upon adsorption and dissociation to form new 3 H species and diffuse into the Zircaloy-4 getters to form metal hydrides (Zr 3 H x ). Therefore, exploring 3 H species ( 3 H 2 , 3 H 2 O) dissociation on the surface of Ni and diffusion across the interface of Ni-plated Zircaloy-4 getters can provide a better understanding of 3H species formation and transport from pellets into the getters.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Benchmark Exercise for the Control Rod Swelling Evaluation

The VTR core has six reactivity control assemblies and three safety assemblies. The control assemblies or primary control rods are adjusted during the normal operation to balance the core reactivity and to control the reactor power. A typical control assembly radial layout is presented in Figure 1. The figure shows the swelled absorber (B 4 C) rod. Initially, helium gas fills the gap between the pin and the cladding before irradiation swelling takes place. For VTR, HT9 steel was selected as the cladding and duct material. The main neutron absorbing material used in the VTR is B 4 C. When residing in the core, the neutronics, thermophysical, and mechanical properties of the materials used in a control assembly will degrade due to accumulated neutron damage. Material degradation limits how long a control assembly can reside in the core. Many phenomena affect the control assembly lifetime, such as the loss of reactivity worth due to B 4 C depletion, the mechanical interaction of the absorber rod and the cladding due to B 4 C swelling, the helium gas buildup in the pin due to B-10 capture, etc. B 4 C swelling, which causes closure of the gap between the absorber rod and the cladding, is usually considered as the main limiting factor from past experience. An initial study was conducted at PNNL to evaluate the irradiation behavior of a VTR control assembly. The evaluation was performed using the CNRD2 code that was initially developed for the FFTF. The study also included an assessment of the VTR control assembly and focused on a 61-pin control assembly design, which is different from that used (37-pin design) in the core design study. The study conducted by PNNL was reviewed independently by ANL. A Python script referred to as the Control Assembly Evaluation Script (CAES) was developed for the independent review and additional assessment of 37-pin control assembly design. The script has focused on the assessment of the absorber rod swelling for its importance in determining the control assembly lifetime. CAES uses geometry, neutronics, materials data as input to predict the swelling of the absorber rod during its residence in the reactor core. The results from CAES showed some non-negligible differences against the PNNL results. Some of the differences can be attributed to the different interpretation of the control rod assembly dimensions. To resolve this issue, a benchmark exercise was proposed. The benchmark specification was developed by PNNL. The benchmark exercise was performed independently at PNNL and ANL using different codes/scripts (CRND2 and CAES). This memo documents the results calculated using the different codes. However, this report is limited to presenting the results obtained. Further investigation of the cause of the observed difference will be performed as part of future activities, pending continuation of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

How effective is the Brink–Axel hypothesis for astrophysical weak rates?

We explore the effectiveness of the Brink–Axel hypothesis (BAH) for the computation of stellar electron capture (EC) and β-decay (BD) rates, namely that the transition strength function depends only upon the transition energy and not upon the details of the initial state. For this purpose, we calculated Gamow–Teller (GT) strength distributions for a selection of sd-shell nuclides, using two different microscopic models, namely the proton–neutron quasiparticle random phase approximation and the full configuration-interaction shell model, taking into account the first 100 states of both the initial and final nuclides. The GT transition strengths among these levels evolve with initial state energy. These transition strength functions we folded into weak-interaction mediated rates in stellar matter, specifically EC and BD rates, for a range of densities 10 g cm –3 ≤ ρ ≤ 10 11 g cm –3 and range of temperatures 1 GK ≤ T ≤ 30 GK. When transitions from excited states were approximated using the BAH, augmented by so-called 'back-resonance' transitions, the rates were affected by up to three orders of magnitude or more at high temperatures and densities. Furthermore the BAH is not a reliable approximation for the calculation of stellar rates, especially in high temperature–density environments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Ab initio translationally invariant nucleon-nucleus optical potentials

We combine the ab initio symmetry-adapted no-core shell model (SA-NCSM) with the single-particle Green's function approach to construct optical potentials rooted in first principles. Specifically, we show that total cross sections and phase shifts for neutron elastic scattering from a 4 He target with projectile energies between 0.5 and 10 MeV closely reproduce the experiment. In addition, we discuss an important new development that resolves a long-standing issue with spurious center-of-mass motion in the Green's function formalism for many-body approaches. Furthermore, the new development opens a path for first-principle predictions of cross sections for elastic scattering of single-nucleon projectiles, nucleon capture, and deuteron breakup reactions, feasible for a broad range of open-shell spherical and deformed nuclei in the SA-NCSM approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Methodology to Generate the Decay Gamma Source for a Second Target Station Target Wedge

This report details the specifics of applying the Position-Averaged Method methodology for the calculation of the rotating lasagna-style target wedge activation and the resulting decay gamma source terms. The methodology uses MCNP® Code Version 6.2.0 with the RNUCS patch coupled with CINDER2008 from the AARE V1.0 package. The process begins with model preparations that include ensuring impurities are included in the materials in the regions of interest along with segmentation of the geometry in order to capture the decay gamma source gradients. This report also details the physics models and other MCNP®, CINDER2008, and ADVANTG 3.2.0 related options used to calculate the decay gamma source terms. ADVANTG 3.2.0 is used to generate a weight window set to aid in convergence of the neutron fluxes in the initial MCNP® protons-on-target particle transport calculation to calculate the neutron fluxes and spallation products in the regions of interest. The neutron fluxes and spallation products are then used in the CINDER2008 transmutation calculation, and the decay gamma sources are generated by the AARE_GAMMA_SOURCE_SCRIPT as a part of the AARE V1.0 package. These decay gamma sources include the energy distributions as well as the necessary distributions to sample the location of the source using the cell rejection method. The gamma sources that are produced using the methodology may be used for a wide variety of analysis involving the target wedge. This report includes the material specifications, physics options, proton source term description, geometrical configuration, and post-processing details.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Supramolecular Control of Redox-Mediated Interactions for Ion-Selective Electroseparations

This report summarizes efforts from the Su Group at the University of Illinois at Urbana-Champaign in researching and developing materials and processes for the recovery of rare earth elements and other metal elements that are critical to the United States economy and national security, under DOE Award No. DE-SC0021409. This report highlights efforts in the recovery of: • Transition metal elements (including arsenic, chromium, molybdenum, rhenium, vanadium): Redox-active metallopolymers are able to capture these elements in their oxyanion form. Direct control of target ion capture/release through electrochemical-mediated switching of metallocene binding sites reduces chemical input and lowers generation of waste in element recovery processes. We study the effects of metallopolymer structure on targeted ion binding through a combination of experimental screening of ion selectivity from multi-component mixtures, computational efforts which elucidate the charge-transfer origins of metallocene binding and provide insight on how binding site structure influences selectivity, and state-of-the-art neutron reflectivity of metallocene polymer films to understand the interplay of polymer-ion (charge-transfer), polymer-solvent (polymer hydrophilicity/hydrophobicity), and ion-solvent (ion hydration) interactions with spatial resolution of polymer films. • Rare earth elements (including yttrium, neodymium, europium, gadolinium, dysprosium, cerium): Here we employ redox-copolymers for recovery of rare earth elements by electrochemically regenerated ion-exchange. Rare earth elements require an alternative route to capture/release than directly with metallopolymers since rare earth elements typically exist in cationic form in water. To address this, copolymerization of traditional ion-exchange moieties with metallocene moieties allowed for capture of rare earth elements with the ion-exchange moieties and electrochemical-mediated release of rare earth elements through oxidation of metallocene moieties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗
DOE OSTI2024

Application of research reactor environments for validation of the IRDFF-II dosimetry cross section library

Activation data from seven different reactor-based reference neutron fields are examined to provide enhanced validation evidence for the newly released IRDFF-II library. A least-squares based spectrum adjustment methodology is used and rigorous statistical metrics demonstrate the consistency of the set of IRDFF-II dosimetry cross sections. The use of these reactor-based neutron fields provides validation evidence for nine more reactions than have been addressed in the 252 Cf(sf) and 235 U(th) benchmark validation testing. The use of covers to shift their energy response increases the power of the response validation but makes it challenging to properly capture response correlations.

Griffin, Patrick↗
DOE OSTI2021

Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation

Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Binary-stripped Stars as Core-collapse Supernovae Progenitors

Most massive stars experience binary interactions in their lifetimes that can alter both the surface and core structure of the stripped star with significant effects on their ultimate fate as core-collapse supernovae. However, core-collapse supernovae simulations to date have focused almost exclusively on the evolution of single stars. We present a systematic simulation study of single and binary-stripped stars with the same initial mass as candidates for core-collapse supernovae (11–21 M⊙). Generally, we find that binary-stripped stars core tend to have a smaller compactness parameter, with a more prominent, deeper silicon/oxygen interface, and explode preferentially to the corresponding single stars of the same initial mass. Such a dichotomy of behavior between these two modes of evolution would have important implications for supernovae statistics, including the final neutron star masses, explosion energies, and nucleosynthetic yields. Binary-stripped remnants are also well poised to populate the possible mass gap between the heaviest neutron stars and the lightest black holes. Our work presents an improvement along two fronts, as we self-consistently account for the pre-collapse stellar evolution and the subsequent explosion outcome. Even so, our results emphasize the need for more detailed stellar evolutionary models to capture the sensitive nature of explosion outcome.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Development for Integrated System-Level Analysis Capabilities in SAM for Molten Salt Reactors

In recent years, there has been renewed interest in Molten Salt Reactors (MSRs) for their potential advantages compared to reactors that rely on solid fuel. In response to such interest, many methods and codes have been developed to capture the unique features of MSRs. Among them, the System Analysis Module (SAM) is a modern system analysis tool that provides fast-running, modest-fidelity, whole-plant transient analysis capabilities, essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. For liquid-fuel MSRs, the complex physics and chemistry involved in MSR operation—such as reactor kinetics, fluid flow, heat transfer, and salt composition dynamics—pose significant challenges for system-level modeling. Specific modeling capabilities are needed for system-level transient simulation. This paper presents recent advancements in SAM capability enhancements for system-level modeling of MSRs, focusing on improved simulation fidelity, computational efficiency, and multi-physics integration. Key enhancements include the development of species transport, Delayed Neutron Precursor (DNP) drift, modified Point Kinetics Equations (PKE), decay heat modeling, key fission product behavior, salt corrosion, and thermal-hydraulic coupling, as well as code robustness and performance enhancements for MSR applications. The code enhancement allows for better predictive accuracy in safety analysis, transient behavior, and operational optimization, thus supporting the design and licensing of next-generation MSRs. Results from case studies are presented to demonstrate the benefits of these enhancements in accurately capturing key reactor transient behaviors.

Hu, Rui (ORCID:0000000237712920)↗
DOE OSTI2025

Stellar population and metal production in AGN discs

ABSTRACT As gravitational wave detections increase the number of observed compact binaries (consisting of neutron stars or blacks), we begin to probe the different conditions producing these binaries. Most studies of compact remnant formation focus either on stellar collapse from the evolution of field binary stars in gas-free environments or on the formation of stars in clusters where dynamical interactions capture the compact objects, forming binaries. But a third scenario exists. In this paper, we study the fate of massive stars formed, accrete gas, and evolve in the dense discs surrounding supermassive black holes. We calculate the explosions produced and compact objects formed by the collapse of these massive stars. Nucleosynthetic yields may provide an ideal, directly observable, diagnostic of the formation and fate of these stars in active galactic nuclei. We present a first study of the explosive yields from these stars, comparing these yields with the observed nucleosynthetic signatures in the discs around supermassive stars with quasars. We show that, even though these stars tend to form black holes, their rapid rotation leads to discs that can eject a considerable amount of iron during the collapse of the star. The nucleosynthetic yields from these stars can produce constraints on the number of systems formed in this manner, but further work is needed to exploit variations from the initial models presented in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Report on initial development of a database of nuclear graphite characteristics based on microstructural characterization

This report outlines the current and future efforts to produce a comprehensive library of microstructures of nuclear graphite and carbon-based materials that are candidate materials for nuclear applications in the United States. This library must contain extensive characterizations of unirradiated graphite materials, a guide to some of the techniques used to characterize graphite, a compendium of characterization data of neutron-irradiated or oxidized material, and a compendium of microstructural information of carbon-based materials. These characterization efforts are being conducted at various length scales to understand these complex materials’ local structure and property relationships. Achieving this goal requires further developing or adapting advanced characterization techniques that capture graphite’s most relevant characteristics. Some of the general objectives of this project are to aid with the material selection, licensing, management, and core assessments of a graphite core by documenting the unirradiated microstructure of relevant grades or by characterizing the evolution of the microstructure under the reactor environment. Moreover, this project aims to provide additional information, guidelines for characterizing graphite, and a protocol to assess a nuclear graphite grade. This report also summarizes some of the initial results and some of the techniques commonly used to characterize nuclear graphite.

36 MATERIALS SCIENCE↗