SEARCH · Search NASA

Results for “DFT calculations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

250 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗

DOE OSTI2024

Cosmic neutrino decoupling and its observable imprints: insights from entropic-dual transport

Abstract Very different processes characterize the decoupling of neutrinos to form the cosmic neutrino background (CνB) and the much later decoupling of photons from thermal equilibrium to form the cosmic microwave background (CMB). The CνB emerges from the fuzzy, energy-dependent neutrinosphere and encodes the physics operating in the early universe in the temperature rangeT∼ 10 MeV toT∼ 10 keV. This is the epoch where beyond Standard Model (BSM) physics, especially in the neutrino sector, may be influential in setting the light element abundances, the necessarily distorted fossil neutrino energy spectra, and other light particle energy density contributions. Here we use techniques honed in extensive CMB studies to analyze the CνB as calculated in detailed neutrino energy transport and nuclear reaction simulations of the protracted weak decoupling and primordial nucleosynthesis epochs. Our moment method, relative entropy, and differential visibility approach can leverage future high precision CMB and light element primordial abundance measurements to provide new insights into the CνB and any BSM physics it encodes. We demonstrate that the evolution of the energy spectrum of the CνB throughout the weak decoupling epoch is accurately captured in the Standard Model by only three parameters per species, a non-trivial conclusion given the deviation from thermal equilibrium and the impact of the decrease of electron-positron pairs. Furthermore, we can interpret each of the three parameters as physical characteristics of a non-equilibrium system. Though the treatment presented here makes some simplifying assumptions including ignoring neutrino flavor oscillations, the success of our compact description within the Standard Model motivates its use also in BSM scenarios. We further demonstrate how observations of primordial light element abundances can be used to place constraints on the CνB energy spectrum, deriving response functions that can be applied for general deviations from a thermal spectrum. Combined with the description of those deviations that we develop here, our methods provide a convenient and powerful framework to constrain the impact of BSM physics on the CνB.

Astronomy & Astrophysics↗

DOE OSTI2023

Four Years of Type Ia Supernovae Observed by TESS: Early-time Light-curve Shapes and Constraints on Companion Interaction Models

Abstract We present 307 type Ia supernova (SN) light curves from the first 4 yr of the Transiting Exoplanet Survey Satellite mission. We use this sample to characterize the shapes of the early-time light curves, measure the rise times from first light to peak, and search for companion star interactions. Using simulations, we show that light curves must have noise <10% of the peak flux to avoid biases in the early-time light-curve shape, restricting our quantitative analysis to 74 light curves. We find that the mean power-law index t β 1 of the early-time light curves isβ 1 = 1.93 ± 0.57, and the mean rise time to peak is 15.7 ± 3.5 days. The underlying population distribution forβ 1 may instead consist of a Gaussian component with mean 2.29, width 0.34, and a long tail extending to values less than 1.0. We find that the data can rarely distinguish between models with and without companion interaction models. Nevertheless, we find three high-quality light curves that tentatively prefer the addition of a companion interaction model, but the statistical evidence for the companion interactions is not robust. We also find two SNe that disfavor the addition of a companion interaction model to a curved power-law model. Taking the 74 SNe together, we calculate 3σupper limits on the presence of companion signatures to control for orientation effects that can hide companions in individual light curves. Our results rule out common progenitor systems with companions having Roche lobe radii >31R ⊙ (separations >5.7 × 10 12 cm, 99.9% confidence level) and disfavor companions having Roche lobe radii >10R ⊙ (separations >1.9 × 10 12 cm, 95% confidence level). Lastly, we discuss the implications of our results for the intrinsic fraction of single degenerate progenitor systems.

Astronomy & Astrophysics↗

DOE OSTI2026

Single particle – MC-ICP-MS for isotopic analysis of uranium particles

Single particle – multi-collector – inductively coupled plasma – mass spectrometry (SP-MC-ICP-MS) was employed to measure a suspension of 1 µm U3O8 particles (∼1.3 pg total U/particle) to determine their individual isotopic compositions of 234U/238U, 235U/238U, and 236U/238U. The effects of different detector combinations for 235U and 238U, including secondary electron multipliers (SEM) and Faraday detectors (1011 and 1013 Ω amplifiers), were explored for accuracy and precision optimization on the observed 235U/238U. The minor isotopic ratios (i.e., 234U/238U and 236U/238U) were analyzed such that the 234U and 236U were monitored on SEM detectors and the 238U was monitored on a Faraday (1011 Ω) detector. Various integration times (5, 10, 25, and 50 ms) were investigated in all detector configurations to gain a better understanding of their impact on sensitivity, accuracy, and precision. For 235U/238U ratios of 1 µm U3O8 particles, a dual Faraday detector measurement with 1011 Ω was the optimal choice; measurement of 1021 particles yielded an average 235U/238U ratio of 0.00170 (14), a −1.8% relative difference (% RD) from the reference value. The minor isotopic compositions were determined to be 0.0000070 (14) and 0.0000758 (48) for the 234U/238U and 236U/238U, respectively. These measurements correspond to <8% and <1% RD from their reference value for the 234U/238U and 236U/238U, respectively. SP-MC-ICP-MS was also able to provide insight into measurement sensitivity. In these individual particles, merely 15 and 165 atto-grams (ag) of 234U and 236U were present (calculated). Initial limits of detection for SP-MC-ICP-MS were determined to be ∼1.0 ag (when measured via SEM detectors). This valuable approach is applicable to areas including nuclear forensics, nuclear safeguards, and geochemical analysis, which require high-precision measurements of uranium within micron-sized particles.

Manard, Benjamin [ORNL] (ORCID:0000000207400627)↗

DOE OSTI2026

Phonon-mediated ultrafast dynamics in self-assembled monolayers of 4-mercaptobenzoic acid on gold

Non-equilibrium interactions between plasmonic metals and adsorbed molecules lie at the heart of emerging applications such as plasmonic photocatalysis and sensing, though the ultrafast charge and energy transfer mechanisms arising from these interactions are not well understood. Herein, we investigate the ultrafast dynamics of Au nano-islands tethered with a self-assembled monolayer (SAM) of electron-withdrawing 4-mercaptobenzoic acid (4MBA) molecules. Ultrafast UV-visible transient absorption spectroscopy following excitation of the interband transition in Au reveals three well-known, characteristic time constants that quantify electron–electron (el–el), electron–phonon (el–ph) and phonon–phonon (ph–ph) scattering lifetimes. When comparing the dynamics of bare Au and 4MBA-Au, we find that the el–ph and ph–ph scattering lifetimes are notably longer in 4MBA-Au. Density functional perturbation theory calculations ascribe the elongation in el–ph lifetimes in 4MBA-Au to the significant coupling of acoustic phonon modes of Au with certain molecular vibrations of 4MBA, leading to decreased spatial overlap between carrier electronic states and the acoustic modes. We speculate that the elongation of ph–ph scattering lifetimes in 4MBA-Au arises due to poor thermal conductivity of the SAM which disrupts efficient energy dissipation from Au to the environment, thus slowing down the thermalization of phonons. This work provides a glimpse into how molecular adsorbates modify the charge carrier and phonon dynamics of Au and sets the stage for further systematic exploration of plasmonic metal–molecule interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Kondo effect in ferromagnetic quantum critical CeRh 6 ⁢Ge 4

The mechanism of a pressure-induced quantum critical point in the heavy fermion ferromagnet CeRh 6 ⁢Ge 4 has attracted interest, as ferromagnetic quantum criticality in a clean itinerant Ce compound is typically avoided. The localized versus itinerant character of the 4⁢𝑓 electrons is a key aspect for understanding this behavior. We investigated the electronic structure of the 4⁢𝑓 shell in CeRh 6 ⁢Ge 4 using core-level photoelectron and x-ray absorption spectroscopy, demonstrating the hybridization of Ce 4⁢𝑓 with the conduction electrons. Linearly polarized x-ray absorption reveals a temperature-dependent linear dichroism consistent with the crystal-electric-field sequence as inferred from the static susceptibility. This dichroism cannot be described by an ionic full-multiplet model alone, but is reproduced by including the Kondo effect within a single-impurity Anderson model in the noncrossing approximation. The Kondo effect mixes higher-lying crystal-field states into a resulting multiorbital ground state with 4⁢𝑓 occupancy, 𝑛 𝑓 ∼ 0.9. Deviations at low temperatures between the measured linear dichroism and calculated dichroism suggest an orbital-dependent Kondo effect. A scenario in which there is a multiorbital ground state and orbital-dependent Kondo hybridization should be a starting point for a model of pressure-induced criticality in CeRh 6 ⁢Ge 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRSDate not supplied

First N 2 Profile for Venus’ Deep Lower Atmosphere

We present the first N 2 profile for Venus’ deep lower atmosphere (<15 km) and a constrained isotopic composition for cloud N 2 [1]. This work directly addresses unresolved questions for Venus using legacy observations. Prior to this work, there were no reported measurements for N 2 abundances at <15 km and the isotopic composition remained unconstrained [2]. The N 2 parameters are critical to understanding the evolution and thermal properties of the atmosphere [3-7]. Our N 2 results were obtained by re-analysis of data acquired in 1978 by the Pioneer Venus Large Probe Neutral Mass Spectrometer (LNMS) [8]. The archived mass spectra from 64.2 to 0.2 km were treated using the analytical procedures specifically developed for the LNMS [9-11]. Judicious peak fitting permitted disambiguation of (A) N 2 + and CO + at 28 u and (B) 14 N 15 N+, 13 CO + , and C 2 H 5 + at 29 u. Quality controls included comparing the (A) LNMS CO + /CO 2 + ratios to the NIST database and literature and (B) fitted counts for CO + to the expected counts of CO + calculated from C 18 O + and 13 CO + using the LNMS 16 O/ 18 O and 12 C/ 13 C ratios (obtained from CO 2 ). Our results show that N 2 is uniformly mixed across the deep lower atmosphere between ~0.2 and 15 km (2.49 ± 0.10 v%). In contrast, N 2 is non-uniformly mixed across the sub-cloud atmosphere and clouds (~15–59 km), where N 2 abundances increase by ~ 2-fold between ~15 km (2.45 ± 0.32 v%) and ~59-51 km (5.21 ± 0.18 v%). Using the cloud data, we also obtained a constrained 15 N/ 14 N ratio (2.93×10 -3 ± 0.13×10 -3 ) and δ 15 N value (-204 ± 35‰). Thus, the LNMS results [1] suggest that (A) the atmosphere is unstable at <15 km, (B) N 2 is not well-mixed >15 km, and (C) the cloud δ 15 N falls between Earth and the solar wind [12, 13]. Comparisons of the N 2 abundances and isotopic compositions for nitrogen, carbon, and oxygen to other Venus measurements will be discussed.

deep lower atmosphere↗

NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

DOE OSTI2026

Mapping the Microhardness of Al Matrix in U-7Mo/Al Dispersion Fuels at Medium and High Burn-up by in situ Nanoindentation

The mechanical degradation of the Al matrix in U-7Mo dispersion fuels was detected after irradiation in previous studies, but more details of the degradation mechanism are not clear yet to clarify the root cause. In this study, microhardness mapping in fresh and irradiated U-7Mo dispersion fuels coated with ZrN, with burn-up of 3.35 × 10 21 and 6.28 × 10 21 f/cm 3 , are measured to reveal the hardening of the Al matrix as a function of distance to fuel particles. Micro cubes were pulled by plasma focused ion beam scanning electron microscopy and tested by in situ nanoindenter with Berkovitch tip. Size effects are measured in medium and high burn-up samples. The higher burn-up specimen exhibits a more pronounced size effect when indent load is lower than 20 mN. Size effect correction models are also calculated by plots fitting, which predicts the true microhardness without size effect. Results shows that the hardness of the Al matrix near the ZrN coating has the highest value, then decreases moving away from the coating and become stable when the distance reaches ∼10 µm. The hardness starts increasing again getting closer to the next fuel particle, which makes the hardness distribution between two fuel particles to be “U” shape.

36 - MATERIALS SCIENCE↗

DOE OSTI2026

An integrated in-situ coordination strategy enabling high-performance layered cathodes for sodium-ion batteries

O3-type layered transition metal oxide cathodes hold tremendous potential in sodium-ion batteries (SIBs) due to their low cost and high energy density. However, the structure instability associated with detrimental phase transitions and severe interface parasitic reactions exacerbate the material's electrochemical performance degradation. Herein, we develop an integrated in-situ coordination strategy via heteroatomic modulation inducing coherent epitaxial layer to collaboratively enhance the overall framework robustness from surface to bulk. The theoretical calculation and multiple in/ex-situ characterizations demonstrate the charge density around oxygen is redistributed, which promotes the electron localization, thus widening the NaO 2 lattice space and accelerating the Na + transport dynamics. Furthermore, the formed strengthened oxygen bond energy effectively distributes the long-range coordination of Mn 3+ O 6 octahedron, thereby alleviating Jahn-Teller distortion and local stress. Importantly, the in-situ formed conformal buffer layer dramatically relieves the adverse interface side reactions, facilitating the construction of robust cathode-electrolyte interface, which ameliorate the whole structure stability of designed materials. Consequently, the optimized NFMZ@NZO-1.0 exhibits the excellent cycling stability with 80.2% capacity retention after 300 cycles at 1C, and delivers a high discharge capacity of 107.1 mAh g −1 at 10C. In conclusion, this distinctive coupling strategy provides valuable insights for developing high-performance layered cathode materials in SIBs.

Coherent epitaxial layer↗

DOE OSTI2026

Block-Type Antiferromagnetism in Single Chain Quasi-One-Dimensional K 3 ⁢Fe 2 ⁢Se 4

One-dimensional (1D) structures provide a unique platform to study the correlated quantum interactions and phase transitions such as unconventional magnetism and superconducting states. Here, we report that iron chalcogenide K 3 ⁢Fe 2 ⁢Se 4 exhibits an unusual block-type canted antiferromagnetic (AFM) order with a clear single chain quasi-1D structure, which is structurally different from the two-leg ladder BaFe 2 ⁢Se 3 , through both experimental measurements and density matrix renormalization group (DMRG) calculations. The narrow bandgap semiconductor K 3 ⁢Fe 2 ⁢Se 4 has a quasi-1D edge-shared FeSe4 tetrahedra chain structure and orders antiferromagnetically below 110 K. The magnetic moments couple antiferromagnetically along the quasi-1D chain direction of the 𝑏 axis and form an up-down-down-up (↑−↓−↓−↑)–like spin structure with a commensurate propagation vector 𝒌=⁢(0,0,0), where block-type spin ↑−↑ or ↓−↓ coupling are between the longer Fe-Fe bonds of the quasi-1D chain. DMRG results show that block antiferromagnetic state is stable in K 3 ⁢Fe 2 ⁢Se 4 and reveal that the block-ordered arrangement of Fe 2.5+ ions spins arise from the competition between ferromagnetic and AFM interaction in the presence of strong electronic correlation. Our research results not only report the discovery of a clear block-type canted antiferromagnetic structure in a real quasi-1D chain material but also provide a theoretical approach to understand the block-type antiferromagnetism in quasi-1D iron chalcogenides.

antiferromagnetism↗

NASA NTRS2026

Final Report Document: Microgravity Medical Eyewash B

Senior Capstone Design Team 15 has been commissioned by the National Aeronautics and Space Administration (NASA) to redesign the current microgravity eye wash station. The current design exhibits four primary limitations: single-use operation, excessive mass and volume, operational complexity requiring coordination with external systems, and dependence on frequent Earth resupply missions. Our updated eyewash design will be deployed on long-range space missions that could last multiple years. Based on these and other requirements from our customer, we have developed our updated eyewash design. Our compact, single-eye system uses capillary-driven flow and a free-floating water ball. This design allows us to eliminate pumps and bulky tubing, decreasing the mass and volume of the overall system. The containment eyecup is held to the astronaut's face with a strap, and a silicone seal prevents leakage while improving user comfort. Contaminated water is contained within a disposal bag containing wicking material to pull it back out of the eyecup after washing the eye. This system is reusable, efficient, and easy to activate quickly. We evaluated our design through calculations, microfluidic testing, and user try-on testing to verify our requirements. Our work demonstrated that steady-state flow was capable with our capillary framework, achieving 1.31 L/min. The prototype construction confirmed mass and volume reductions, and try-on-testing evaluated the ability to put on the system quickly. Almost all major design requirements were achieved, and water loss could be validated by testing in microgravity. Our design went through a number of iterations to reach the final solution presented in this report. These changes were driven by our testing and collaboration from each member of the team. Updated models optimized tube placement and geometry of the eyecup to better direct flow as it pools across the eye. Changes to our tubing bends improved capillary efficiency and increased the flow rate we were capable of achieving. And improved bracket placements adjusted the fit, comfort, and seal of the eyecup to the astronaut's face. Our capillary-based eyewash system is technically feasible and has been theoretically validated to operate effectively in microgravity conditions. Testing and analysis indicate that it can meet or exceed the required performance metrics, including flow rate and safety constraints. Overall, the design represents a viable alternative to current ISS-dependent eyewash systems for future space missions.

Jeffrey Allen↗

DOE OSTI2026

Modeling Flow in a Microreactor Core: From High-Fidelity CFD to Subchannel Analysis

The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

NASA NTRS2020

Evaluation of Fatigue Damage Accumulation Functions for Delamination Initiation and Propagation

The present report follows on the cohesive fatigue damage model methodology proposed in NASA-TP-2018-219838. In that report, an empirical function describing the incremental damage due to cyclic loading was used to calculate fatigue damage within a cohesive formulation. The form of the function was developed such that, when integrated at a constant stress amplitude from no damage to failure, it produces a life versus load response that is consistent with an S-N diagram. Therefore, the parameters of the damage model could be obtained by fitting the model predictions to an S-N diagram. The finite element analyses performed demonstrate that the cohesive fatigue accumulation function provides a link between the S-N diagram that describes crack initiation, and the Paris law that characterizes the rate of crack propagation. However, when the model was proposed, it was not known whether the form of the damage accumulation function associated with a desired S-N diagram is unique and, if not, if the link between S-N and the Paris law is unique and independent of the fatigue function selected. In the effort described herein, several alternative forms of the damage function that reproduce the desired features of S-N diagrams were found and evaluated. The effects of each of these functions on the predicted parameters of the Paris law and the propagation threshold are discussed. The results indicate that the predicted exponent m of the Paris law is indeed independent of the damage accumulation function. However, different functions predict different values for the pre-factor C of the Paris law. Therefore, the proper damage accumulation function must be selected by comparison with experiments. One of the new damage accumulation functions proposed herein was found to be particularly useful for analysis because of the ease with which the model parameters can be determined with a minimal amount of experimental information. The effectiveness of the proposed methodology and damage function was demonstrated by conducting analyses of a double cantilever beam test, a mixed-mode bending test, and a three-point bending test of a skin/doubler specimen. The results indicate that the same set of model parameters can provide accurate predictions of the rate of fatigue crack propagation for a variety of material interfaces, mode mixities, load levels, and stress ratios.

DKIN/Stiffener Debonding↗

DOE OSTI2026

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2004

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda↗

NASA NTRS2026

X-Hab 2026: LiDAR-Powered Autonomous Charging Service Capability for Surface Rovers and Systems

This document details the timeline of the LiDAR-Powered Autonomous Charging Service Capability for Surface Rovers and Systems project, initiated by the Fall 2025 semester class and completed by the Spring 2026 semester class. This project focuses on developing a fully autonomous system composed of a mobile surface rover and an induction charging station with a robotic arm, both controlled by their own NVIDIA Jetson Orin Nano. Structural improvements to the rover suspension system and body eliminated excessive camber, reduced stress and strain on the plexiglass body, and improved maneuverability and durability of the rover. The charging station robotic arm was fully redesigned to increase reach while minimizing weight and increasing misalignment tolerance during docking on uneven terrain. Electrical system improvements addressed previous torque and power limitations of both the rover and charging station arm. High-torque servo motors were selected based on updated calculations which incorporated terrain slope and Factor of Safety, enabling zero-point turning for the rover and increased payload capacity of the charging station arm. Significant progress was made in autonomy and perception. The rover now employs 3D LiDAR and SLAM mapping for localization, mapping, and path planning. The Battery Monitoring System (BMS) was created to coordinate battery management between the rover and charging station. The BMS provides continuous monitoring of battery state of charge, temperature, current, and will enable the autonomous initiation, execution, and termination of the charging cycle via Bluetooth communication. Testing of the WIBOTIC induction charging system demonstrated reliable power transfer under both aligned and misaligned conditions. This project demonstrated the ability of an autonomously navigating surface rover to independently plan a path to the charging station, dock, and the charging station to autonomously deploy a robotic charging arm and initiate charging of the rover. This work details the progress made to demonstrate the feasibility of autonomous surface rover navigation and recharging systems.

Megan Steele↗