Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lunar Surface Thermal Modeling”

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.

889 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Feedback, physics, and forecasts: The emerging paradigm of machine learning-driven battery research

Machine learning (ML) is reshaping how we understand, predict, and optimize electrochemical systems. In batteries, ML accelerates discovery across chemistry, design, and operation by transforming massive experimental and simulated datasets into predictive, interpretable models. This review consolidates a decade of progress in ML-driven battery innovation, from early-cycle feature extraction to operando image analysis and physics-informed modeling. We categorize approaches by data domain and physical fidelity, emphasizing interpretable ML for diagnostics, reinforcement learning for control, and multi-objective optimization for lifetime extension strategies. Additionally, we demonstrate how integrated models accelerate discovery, reduce testing time, and guide sustainable design. Economic analyses furthermore illustrate how these advances can lower cost per cycle and improve circularity. Together, these developments chart a path toward self-optimizing, sustainable battery technologies.

artificial intelligence↗
DOE OSTI2026

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗
DOE OSTI2026

Development of Ternary Transition Metal Oxide Catalysts for Oxygen Evolution Reaction

Electrochemical water splitting, a promising method for green hydrogen production, is currently hindered by the high cost of precious-group metal (PGM)-based catalysts for the oxygen evolution reaction (OER). This study addresses this challenge by advancing the development of catalysts for OER, focusing on the development of high-performance PGM-free catalysts using nickel-iron-cobalt (NiFeCo)-based aerogels. The catalysts were synthesized via a sol-gel method and critical point drying to achieve a highly porous structure with an exceptionally high surface area. The designed catalyst structure provides an ideal platform for maximizing catalytic active sites and enhancing mass transport kinetics. Co has been systematically incorporated into the current PGM-free state-of-the-art, NiFeOx catalyst, and the metal ratios have been optimized. In addition to the experimental studies, density functional theory calculations were performed to study the material’s properties of this ternary catalyst and the effect of Co addition on enhancing OER catalysis.

Mizrahi, Michal [Bar-Ilan University, Ramat Gan, I↗
DOE OSTI2025

The phase diagram of quantum chromodynamics in one dimension on a quantum computer

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms.

Quantum information↗
DOE OSTI2024

Factors Influencing Gas Evolution from High‐Nickel Layered Oxide Cathodes in Lithium‐Based Batteries

Abstract Gas evolution from high‐nickel layered oxide cathodes (>90% Ni) remains a major issue for their practical application. Gaseous species, such as CO 2 , O 2 , and CO, that are evolved at high states of charge (SOC) worsen the overall safety of batteries, as pressure build‐up within the cell may lead to cell rupture. Since these gasses are produced during cathode degradation, tracking the formation of gasses is also important in diagnosing cathode failure. Online electrochemical mass spectrometry (OEMS) is a powerful in situ technique to study gas evolution from the cathode during high‐voltage charge. However, the differences in the OEMS experimental setups between different groups make it challenging to compare results between groups. In this perspective, the various factors that influence gas evolution based on the OEMS results collected in this group are presented. The focus is on the conditions that lead to gas release, with a particular emphasis on reactive oxygen formation and subsequent chemical reactions with the electrolyte. Promising strategies, such as electrolytes, compositional tuning, and surface coatings that are effective at suppressing gas evolution from the cathode are highlighted. Critical insights into mitigating cathode degradation and gas evolution are provided to guide the development of safer, high‐energy batteries.

Chemistry↗
DOE OSTI2025

Socio-Techno-Economic Feasibility of Deep Geothermal with EGS for Residential District Heating in the Northeastern United States

In this paper, we present socio-techno-economic simulation results for enhanced geothermal system (EGS) district heating systems in the northeastern United States. We applied the geospatial Distributed Geothermal Market Demand (dGeo) tool coupled with GEOPHIRES to evaluate the feasibility of geothermal deep direct-use, utilizing EGS reservoirs for providing residential heating with district heating systems. The results were combined with JEDI results to assess social impacts. Simulation results assuming doublet EGS reservoirs with target production temperature of 80 degrees C and newly constructed district heating systems indicate a wide range in system levelized cost of heat (LCOH) values ranging from ~$15/MMBtu to over $1,000/MMBtu, with most the attractive regions urban areas with high total thermal demand and high thermal demand density. We estimate about 60 GWth in potential installed capacity with LCOH values under $50/MMBtu and 100 GWth in potential capacity with LCOH values under $100/MMBtu, respectively.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

A Compact Cryogenic Environment for In Situ Neutron Diffraction Under Mechanical Loading

Understanding the deformation mechanisms of materials at cryogenic temperatures is crucial for cryogenic engineering applications. In situ neutron diffraction is a powerful technique for probing such mechanisms under cryogenic conditions. In this study, we present the development of a compact cryogenic environment (CCE) designed to facil-itate in situ neutron diffraction experiments under mechanical loading at tempera-tures as low as 77 K. The CCE features a polystyrene foam cryogenic chamber, alumi-num blocks serving as neutron-transparent cold sinks, a liquid nitrogen dosing system for cryogen delivery, a nitrogen gas flow control system for thermal management, a process controller for temperature control, and a pair of thermally isolated grip adapt-ers for mechanical testing. The CCE achieves reliable temperature control with mini-mal neutron attenuation. Utilizing this setup, we conducted three in situ neutron dif-fraction tensile tests on a 316L stainless steel at 77 , 173 , and 298 K, respectively. The results highlight the pronounced effects of cryogenic temperatures on the material's deformation mechanisms, underscoring both the significance of cryogenic deformation studies and the effectiveness of the CCE.

Yu, Dunji [ORNL] (ORCID:0000000189467851)↗
DOE OSTI2026

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗
DOE OSTI2026

Polymer coatings on zirconia microspheres via rotating flow fluid dynamics

Conventional tristructural isotropic (TRISO) coatings for nuclear fuels require multi-step and expensive formation processes; any breakage of the coatings may lead to fission species release. Polymer derived ceramic (PDC) coatings can be a suitable alternative to address these issues. In this study, allylhydridopolycarbosilane (SMP-10) coatings were created on yttria stabilized zirconia (YSZ) microspheres using a Rotating Flow Fluid Dynamics (RFFD) coating method. The effects of curing temperature, rotation speed, and coating cycle/time on the coating were analyzed. Surface functionalization of YSZ microspheres with NaOH resulted in good adhesion between the polymer precursor and the YSZ kernel particles. Spectroscopy analysis revealed complete curing of SMP-10 coated YSZ at 160 °C. Rotation speed and coating time significantly affect the coating characteristics. After 3 cycles at lower rotation speed (50 rpm) for 10 min each, the coating obtained was uniform and homogeneous. In comparison, the coating showed almost 10 times less eccentricity at a high rotation speed of 300 rpm. Compared with the coatings prepared at 300 rpm condition, the coatings have higher sphericity and lower eccentricity at 50 rpm. Overall, this study provides a novel and effective route for fabricating and curing SMP-10 precursor coatings on YSZ microspheres.

Ravi, Nivetha [University of Alabama, Birmingham]↗
DOE OSTI2025

Topological prethermal strong zero modes on superconducting processors

Abstract Symmetry-protected topological phases 1–4 cannot be described by any local order parameter and are beyond the conventional symmetry-breaking model 5 . They are characterized by topological boundary modes that remain stable under symmetry respecting perturbations 1–4,6–8 . In clean, gapped systems without disorder, the stability of these edge modes is restricted to the zero-temperature manifold; at finite temperatures, interactions with mobile thermal excitations lead to their decay 9–11 . Here we report the observation of a distinct type of topological edge mode 12–14 , which is protected by emergent symmetries and persists across the entire spectrum, in an array of 100 programmable superconducting qubits. Through digital quantum simulation of a one-dimensional disorder-free stabilizer Hamiltonian, we observe robust long-lived topological edge modes over up to 30 cycles for a wide range of initial states. We show that the interaction between these edge modes and bulk excitations can be suppressed by dimerizing the stabilizer strength, leading to an emergent U(1) × U(1) symmetry in the prethermal regime of the system. Furthermore, we exploit these topological edge modes as logical qubits and prepare a logical Bell state, which exhibits persistent coherence, despite the system being disorder-free and at finite temperature. Our results establish a viable digital simulation approach 15–18 to experimentally study topological matter at finite temperature and demonstrate a potential route to construct long-lived, robust boundary qubits in disorder-free systems.

Science & Technology - Other Topics↗
DOE OSTI2025

Hydrogen diffusion induced dislocation transformations in a nickel superalloy

The diffusion of hydrogen in metals and alloys induces embrittlement that can adversely affect the structural properties. We examine the adsorption and diffusion of hydrogen in Inconel-718 (IN-718), and scrutinize the ensuing effects on the dislocation behavior in the alloy to elucidate the fundamental mechanisms of hydrogen-microstructure interactions from classical molecular simulations. Hydrogen adsorption increases with time until the surface saturates, while hydrogen diffusion exhibits strong temperature dependence, with diffusion coefficients converging above 1300 K regardless of the initial hydrogen concentration in the alloy. The diffusion in IN-718 is significantly sluggish than in pure Ni, Fe, or Cr, and is strongly impacted by hydrogen concentrations, resulting in an order of magnitude higher diffusion coefficient for hydrogen (10-14 m2/s relative to 10-15 m2/s) at high concentrations, especially below 600 K. Hydrogen diffusion coefficient varies from 10-12 to 10-15 m2/s in IN-718 depending on temperature (500–1400 K). More critically, our results reveal that increasing hydrogen concentration induces microstructural changes in the alloy, transforming perfect dislocations into stair-rods and Shockley partials, with higher temperatures favoring the latter. The results are significant for hydrogen fuel applications to gain insights into the materials chemistry for designing safer and more efficient propulsion systems, particularly in high-performance environments related to controlled hydrogen combustion applications.

Banerjee, Tanumoy↗
DOE OSTI2026

Draft ASME Code Case to qualify L-PBF 316H material for Section III, Division 5 applications

This report documents the AMMT program’s development and submission of a draft ASME Code Case to qualify Laser Powder Bed Fusion (L PBF) Type 316H stainless steel for Section III, Divi-sion 5 Class A and SM high temperature nuclear applications. It summarizes the technical basis, the comprehensive high temperature mechanical test database assembled between 2023–2026, and the proposed code language and qualification framework submitted to ASME. The work was co-ordinated across multiple national laboratories and leverages prior ASME efforts to integrate additive manufacturing into the Boiler & Pressure Vessel Code. The body of the report describes the experimental database and analysis supporting the Code Case: tensile, creep, fatigue, creep fatigue, and thermal aging tests collected from multiple additive manufacturing sites, machine types, and powder lots, with material processed by a solution anneal heat treatment. The dataset — including both full size and subsized specimens and tests oriented parallel and perpendicular to build direction — shows limited tensile anisotropy, tensile properties comparable to wrought 316H, creep strength within the scatter of wrought material, but markedly reduced creep ductility above about 650 °C associated with rapid σ phase formation in L PBF microstructures. The draft Code Case itself prescribes a staged qualification model (manufacturing process qualification, component qualification, and per build witness testing), treats L PBF components as equivalent to Type 316 weld metal for design and inspection, and requires mechanical, chemical, and metallographic controls tied to ASTM/ISO 52946. Key acceptance criteria include tensile tests within a 90% prediction interval of the AMMT dataset, a creep fatigue screening test adapted from ASME Section III, Division 5, Subsection HB, HBB 2800 but with the cycle acceptance reduced to 100 for L PBF material, and double volumetric inspection of production components. The report concludes that the present data support treating L PBF 316H as analogous to conventional fusion weld metal for Division 5 design and inspection, while highlighting important caveats: the σ phase driven loss of creep ductility above ~650 °C, preliminary indications of enhanced creep fatigue sensitivity in some lots, and remaining gaps in long term aging and additional cyclic testing. Recommended next actions include completing outstanding cyclic and long duration creep/aging tests on the solution annealed condition, supporting inclusion of the 316H chemistry and heat treatment in ASTM/ISO 52946, and continuing engagement with ASME and NRC during balloting and review to enable industry adoption.

Messner, Mark C. (ORCID:0000000200404385)↗
DOE OSTI2025

In‐Plane Anisotropy in van der Waals NiTeSe Ternary Alloy

Abstract The anisotropic properties of materials profoundly influence their electronic, magnetic, optical, and mechanical behaviors and are critical for a wide range of applications. In this study, the anisotropic characteristics of Ni‐based van der Waals materials, specifically NiTe 2 and its alloy NiTeSe, utilizing a combination of comprehensive scanning tunneling microscopy (STM), angle‐resolved photoemission spectroscopy (ARPES), and density functional theory (DFT) calculations, are explored. Unlike 1T‐NiTe 2 , which exhibits trigonal in‐plane symmetry, the substitution of Te with Se in NiTe 2 (resulting in the NiTeSe alloy) induces a pronounced in‐plane anisotropy. This anisotropy is clear in the STM topographs, which reveal a distinct linear order of charge distribution. Corroborating these observations, ARPES measurements and DFT calculations reveal an anisotropic Fermi surface centered at the point, which is notably elongated along the k y direction, leading to directional variations in in‐plane carrier velocities. Consequently, the Fermi velocity is highest along the k x direction where the linear charge distribution aligns in real space and is lowest along the k y direction. These findings offer valuable insights into the tunability of anisotropic properties in ternary transition metal dichalcogenide systems, highlighting their potential applications in the development of anisotropic electronic and optoelectronic devices.

Chemistry↗
DOE OSTI2024

Massive red spiral galaxies in SDSS-IV MaNGA survey

ABSTRACT Massive red spiral galaxies (MRSGs) are supposed to be the possible progenitors of lenticular galaxies (S0s). We select a large sample of MRSGs ($M_*\gt 10^{10.5}\rm {\rm M}_{\odot }$) from Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) DR17 using the g − r colour versus stellar mass diagram, along with control samples of blue spirals and S0s. Our main results are as follows: (1) After comparing the Sérsic index, concentration parameter, asymmetry parameter distribution, size–mass relation, and Σ1 (stellar mass surface density within the central 1 kpc)−mass relation, we find MRSGs are similar to S0s and have more compact and symmetric structures than blue spirals. MRSGs also resemble S0s in Dn4000, metallicity, Mgb/$\rm \left\langle Fe \right\rangle$, and V/σ radial profile. (2) By using MaNGA 2D spectra data, we separate the spatial regions into inner (R < 0.8Re) and outer (0.8 < R < 1.5Re) regions, and detect residual star formation in the outer regions of MRSGs. (3) When we select a sub-sample of MRSGs with NUV − r > 5, we find that they are completely star formation quenched in both inner and outer regions. Compared to optically selected MRSGs, NUV − r selected MRSGs appear to be more concentrated and have more massive dark matter haloes. The similarities between S0s and MRSGs suggest the possible evolutionary trend between MRSGs and S0s.

Astronomy & Astrophysics↗
DOE OSTI2026

The impact of aerosol mixing state on immersion freezing: insights from classical nucleation theory and particle-resolved simulations

Immersion freezing, initiated by ice-nucleating particles (INPs) in supercooled aqueous droplets, plays an important role in the formation of ice crystals within clouds. The efficiency of immersion freezing depends strongly on INP composition and, crucially, on the mixing state – how chemical species are distributed across the particle population. Here, we quantify the impact of aerosol mixing state on immersion freezing using a combined theoretical and particle-resolved modeling approach. We derive analytical expressions for the frozen fraction of internally and externally mixed INP populations based on classical nucleation theory, showing that the frozen fraction is sensitive to whether ice-active species are present in all particles or only in a subset of the population. We introduce a multi-species immersion freezing scheme into the particle-resolved model PartMC, using the water activity-based immersion freezing model (ABIFM) to compute freezing probabilities for mixed-composition particles. To improve computational efficiency, we implement a Binned Tau-Leaping algorithm and demonstrate an order-of-magnitude speedup with minimal accuracy loss. Simulations reproduce the analytical trends in limiting cases and extend the analysis to more general aerosol populations, where mixing state continues to exert a substantial control on frozen fraction. Sensitivity analyses across particle size, species type, and cooling condition reveal that the mixing state effect is most pronounced when small amounts of highly efficient INPs are mixed with less efficient materials. These findings underscore the need to represent aerosol mixing state explicitly in models of heterogeneous ice nucleation to reduce uncertainty in cloud-phase partitioning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Physical thickness characterization of the FRIB production targets

The FRIB heavy-ion accelerator, commissioned in 2022, is a leading facility for producing rare isotope beams (RIBs) and exploring nuclei beyond the limits of stability. These RIBs are produced via reactions between stable primary beams and a graphite target. Approximately 20–40% of the primary beam power is deposited in the target, requiring efficient thermal dissipation. Currently, FRIB operates with a primary beam power of up to 20 kW. To enhance thermal dissipation efficiency, a single-slice rotating graphite target with a diameter of approximately 30 cm is employed. The effective target region is a 1 cm-wide outer rim of the graphite disc. To achieve high RIB production rates, the areal thickness variation must be constrained within 2%. This paper presents physical thickness characterizations of FRIB production targets with various nominal thicknesses, measured using a custom-built non-contact thickness measurement apparatus.

Instrumentation for radioactive beams (fragmentati↗
DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dark Energy Survey year 6 results: Magnification modeling and its impact on galaxy clustering and galaxy-galaxy lensing cosmology

Gravitational lensing magnification alters the observed spatial distribution of galaxies and must be accounted for to prevent biases in cosmological probes of the large-scale structure. We investigate its effects on the Dark Energy Survey Year 6 galaxy clustering and galaxy-galaxy lensing analyses using the fiducial lens (position tracer) sample M ag L im++. Magnification bias is parameterized by a coefficient that describes the response of the number of selected objects per unlensed area element to a change in the lensing convergence. We quantify this coefficient using the BALROG synthetic source injection catalog to account for the complexity of the selection function, and compare these results with simplified estimates. The resulting values of the magnification coefficients for each redshift bin are [3.16 ± 0.08, 2.76 ± 0.21, 4.09 ± 0.15, 4.42 ± 0.16, 4.90 ± 0.29, 4.83 ± 0.25]. Relative to Year 3, this analysis provides more precise and accurate magnification bias estimates through a larger BALROG area and reweighting to better match the data properties. Here, the cosmological results are robust when tested against various magnification parameter prior choices and also when adding cross-clustering between lens redshift bins. Neglecting magnification, however, introduces significant systematic shifts: relative to the fiducial analysis with Gaussian priors centered on the BALROG -derived estimates, we observe shifts of 1.37σ in S 8 and -0.84σ in Ω m (with cosmic shear included: -0.61σ in S 8 and -0.71σ in Ω m ), in agreement with findings from simulated data, demonstrating that magnification must be modeled to avoid biases. Freeing the magnification bias in lens bin 2 leads to unphysical negative values, further justifying its exclusion from the fiducial Year 6 analysis.

Cosmological parameters↗