Search NASASearch

SEARCH · Search NASA

Results for “Initialization method”

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 55 records · Page 3

Transport coefficient approach for characterizing nonequilibrium dynamics in soft matter

Nonequilibrium states in soft condensed matter require a systematic approach to characterize and model materials, enhancing predictability and applications. Among the tools, X-ray photon correlation spectroscopy (XPCS) provides exceptional temporal and spatial resolution to extract dynamic insight into the properties of the material. However, existing models might overlook intricate details. We introduce an approach for extracting the transport coefficient, denoted as $J(t)$, from the XPCS studies. This coefficient is a fundamental parameter in nonequilibrium statistical mechanics and is crucial for characterizing transport processes within a system. Our method unifies the Green–Kubo formulas associated with various transport coefficients, including gradient flows, particle–particle interactions, friction matrices, and continuous noise. We achieve this by integrating the collective influence of random and systematic forces acting on the particles within the framework of a Markov chain. We initially validated this method using molecular dynamics simulations of a system subjected to changes in temperatures over time. Subsequently, we conducted further verification using experimental systems reported in the literature and known for their complex nonequilibrium characteristics. The results, including the derived $J(t)$ and other relevant physical parameters, align with the previous observations and reveal detailed dynamical information in nonequilibrium states. This approach represents an advancement in XPCS analysis, addressing the growing demand to extract intricate nonequilibrium dynamics. Further, the methods presented are agnostic to the nature of the material system and can be potentially expanded to hard condensed matter systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Nonunitary Variational Quantum Eigensolver with the Localized Active Space Method and Cost Mitigation

Accurately describing strongly correlated systems with affordable quantum resources remains a central challenge for quantum chemistry applications on near and intermediate term quantum computers. The localized active space self-consistent field (LASSCF) approximates the complete active space self-consistent field (CASSCF) by generating active space-based wave functions within specific fragments while treating interfragment correlation with mean-field approach, hence is computationally less expensive. Hardware-efficient ansatzes (HEA) offer affordable and shallower circuits, yet they often fail to capture the necessary correlation. Previously, Jastrow-factor-inspired nonunitary qubit operators were proposed to use with HEA for variational quantum eigensolver (VQE) calculations (so-called nuVQE), as they do not increase circuit depths and recover correlation beyond the mean-field level for Hartree–Fock initial states. Here, in this study, we explore running nuVQE with LASSCF as the initial state. The method, named LAS-nuVQE, is shown to recover interfragment correlations, reach chemical accuracy with a small number of gates (<70) in both H 4 and square cyclobutadiene (C 4 H 4 ), and produces more accurate energetics than its HEA counterparts at all circuit depths. To further address the inherent symmetry-breaking in HEA, we implemented spin-constrained LAS-nuVQE to extend the capabilities of HEA further and show spin-pure results for square cyclobutadiene. We also mitigate the increased measurement overhead of nuVQE via Pauli grouping and shot-frugal sampling, reducing measurement costs by up to 2 orders of magnitude compared to ungrouped operator, and show that one can achieve better accuracy with a small number of shots (10 3–4 ) per one expectation value calculation compared to noiseless simulations with one or two orders of magnitude more shots. Finally, wall clock time estimates show that, with our measurement mitigation protocols, nuVQE becomes a cheaper and more accurate alternative than vanilla VQE with HEA. Taken together, these developments illustrate a practical pathway toward performing multireference chemical simulations with accuracy and affordable resources on today’s quantum hardware, achieving both accuracy and affordability in challenging correlated systems.

Wang, Qiaohong [Univ. of Chicago, IL (United State

Insights Into Thermal Runaway Mechanisms: Fast Tomography Analysis of Metal Agglomerates in Lithium-Ion Batteries

Thermal Runaway (TR) in lithium-ion batteries (LIB) is a critical technological and social concern. Whilst such events are rare, TR is characterized by uncontrollable heating leading to catastrophic failures. To deepen the understanding of the failure process and subsequently develop more accurate TR prediction models and as a result safer battery systems, we present in this work high-speed X-ray tomography for in-depth investigations of the copper current collector melting and agglomeration during TR. The melting process presents valuable real-time internal information about heat evolution during TR, previously challenging to access but crucially important for validating TR models. In this work, controlled failure studies combined with high-speed X-ray tomography were performed on two different commercial LIB models, subjecting them to both external heating and nail penetration to induce TR. Through real-time observation via high-speed tomography, followed by segmentation, rendering, and analysis, the formation of copper agglomerates was qualitatively and quantitatively characterized and visualized for the first time. Agglomerates tended to form either from the battery's outermost layers or centrally, depending on the method of TR initiation, and gives an indirect insight into the internal temperature evolution and distribution. Moreover, an initial comparative analysis between the battery models also revealed differences in agglomerate size, which has been linked to the thicker copper current collectors of one of the cell models. We further discuss the impact of larger copper agglomerates on heat distribution and safety. This study not only sheds light on the intricate dynamics of TR in LIBs but also underscores the pivotal role of 'gold-standard' imaging techniques in advancing battery safety, crucial for the robust modeling of TR and the future design of electric vehicle safety systems.

25 ENERGY STORAGE

A generalizable machine learning-assisted fast Fourier transform algorithm to simulate the large strain phenomena in polycrystalline materials

Machine learning methods have shown initial promise in constitutive modeling for single crystals or homogenized polycrystals, delivering notable computational efficiency. However, existing machine learning-based constitutive models often lack generalizability, limiting their application across diverse boundary value problems. This study introduces a thermodynamics-informed artificial neural network model to accelerate rate-tangent crystal plasticity fast Fourier transform simulations for cross-scale deformation behaviors of polycrystals under complex loading. Our model integrates microstructural variability and local interactions effectively. To address local effects in each grain, we employ K-means clustering to group Gauss points within the microstructure into clusters assumed to be in similar mechanical states. This approach, based on self-clustering analysis, extends model scope from macroscopic stress response to the granular level, capturing mechanical responses and orientation evolution across grains. This reduces the number of nonlinear problems to solve, with cluster responses propagated throughout each group. The thermodynamics-based artificial neural network-extracted features are further processed using local material state clusters to account for history-dependent deformation and evolving microstructures. Additionally, representative volume element simulations with rate-tangent crystal plasticity fast Fourier transform provide reliable datasets for model training. The proposed model demonstrates high efficiency, accuracy, self-consistency, and enhanced generalizability in predicting strain–stress responses and orientation evolution at both individual grain and aggregate scales under complex loading conditions, such as biaxial tension and arbitrary loading scenarios.

36 MATERIALS SCIENCE

Fast-RF-Shimming: Accelerate RF shimming in 7T MRI using deep learning

Ultrahigh field (UHF) Magnetic Resonance Imaging (MRI) offers an elevated signal-to-noise ratio (SNR), enabling exceptionally high spatial resolution that benefits both clinical diagnostics and advanced research. However, the jump to higher fields introduces complications, particularly transmit radiofrequency (RF) field ($B^{+}_{1}$) inhomogeneities, manifesting as uneven flip angles and image intensity irregularities. These artifacts can degrade image quality and impede broader clinical adoption. Traditional RF shimming methods, such as Magnitude Least Squares (MLS) optimization, effectively mitigate $B^{+}_{1}$ inhomogeneity, but remain time-consuming. Recent machine learning approaches, including RF Shim Prediction by Iteratively Projected Ridge Regression and other deep learning architectures, suggest alternative pathways. Although these approaches show promise, challenges such as extensive training periods, limited network complexity, and practical data requirements persist. In this paper, we introduce a holistic learning-based framework called Fast-RF-Shimming, which achieves a 5000 ​× ​speed-up compared to the traditional MLS method. In the initial phase, we employ random-initialized Adaptive Moment Estimation (Adam) to derive the desired reference shimming weights from multi-channel $B^{+}_{1}$ fields. Next, we train a Residual Network (ResNet) to map $B^{+}_{1}$ fields directly to the ultimate RF shimming outputs, incorporating the confidence parameter into its loss function. Finally, we design Non-uniformity Field Detector (NFD), an optional post-processing step, to ensure the extreme non-uniform outcomes are identified. Comparative evaluations with standard MLS optimization underscore notable gains in both processing speed and predictive accuracy, which indicates that our technique shows a promising solution for addressing persistent inhomogeneity challenges.

Deep learning

Chemical Activation of Frontal Ring-Opening Metathesis Polymerization

In this study, we present an innovative approach to frontal ring-opening metathesis polymerization (FROMP) through chemical activation. By locally concentrating Grubbs-type initiators, we accelerate exothermic polymerization at the activation site. Sufficient thermal energy initiates a polymerization front that travels through the remainder of the sample even with reduced initiator concentrations in the bulk resin. Initiation location and timing are controlled by adjusting the initiator reactivity, the properties of the transport solvent, and the air–resin interface. Our technique enables multipoint front activation, either through synchronized injection or by timing the relative reactivity of various initiators. This accessible method of chemical activation enables on-demand production of polymer and composite materials using only the thermal energy generated by the reaction, removing the need for external power sources and enabling deployment in remote or infrastructure-limited environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Study of e + e − → π + π − π 0 at s from 2.00 to 3.08 GeV at BESIII

With the data samples taken at center-of-mass energies from 2.00 to 3.08 GeV with the BESIII detector at the BEPCII collider, a partial wave analysis on the e + e − → π + π − π 0 process is performed. The Born cross sections for e + e − → π + π − π 0 and its intermediate processes e + e − → ρ π and ρ ( 1450 ) π are measured as functions of s . The results for e + e − → π + π − π 0 are consistent with previous results measured with the initial state radiation method within one standard deviation, and improve the uncertainty by a factor of ten. By fitting the line shapes of the Born cross sections for the e + e − → ρ π and e + e − → ρ ( 1450 ) π , a structure with mass M = 2119 ± 11 ± 15 MeV / c 2 and width Γ = 69 ± 30 ± 5 MeV is observed with a significance of 5.9 σ , where the first uncertainties are statistical and the second ones are systematic. This structure can be interpreted as an excited ω state. Published by the American Physical Society 2024

Astronomy & Astrophysics

Physics-constrained superresolution diffusion for six-dimensional phase space diagnostics

Adaptive physics-constrained superresolution diffusion is developed for noninvasive virtual diagnostics of the six-dimensional (6D) phase space density of charged particle beams. An adaptive variational autoencoder embeds initial beam condition images and scalar measurements to a low-dimensional latent space from which a 32 6 pixel 6D tensor representation of the beam's 6D phase space density is generated. Projecting from a 6D tensor generates physically consistent two-dimensional projections. Physics-guided superresolution diffusion transforms low-resolution images of the 6D density to high resolution 256 × 256 pixel images. Unsupervised adaptive latent space tuning enables tracking of time-varying beams without knowledge of time-varying initial conditions. The method is demonstrated with experimental data and multiparticle simulations at the HiRES UED. The general approach is applicable to a wide range of complex dynamic systems evolving in high-dimensional phase space. The method is shown to be robust to distribution shift without retraining. Published by the American Physical Society 2025

43 PARTICLE ACCELERATORS

Capturing Infrastructure Interdependencies for Power Outages Prediction During Extreme Events

As extreme weather events such as hurricanes, severe thunderstorms, and floods grow in frequency and intensity, the disruption of power grid systems poses significant challenges, including widespread electrical outages, economic losses, and threats to public safety. This paper presents a forward-looking approach that leverages geographical graph-based machine learning models to predict county-level maximum power outages during such events. By capturing the intricate interdependencies within power system networks, our approach aims to provide precise and actionable predictions that can optimize emergency response efforts and enhance grid resilience. Through the integration of real-world data, including hurricane advisories and power outage records, we have trained and benchmarked multiple machine learning models, demonstrating the feasibility and potential of this method. While our initial results are promising, this paper also charts a course for advancing these models, addressing the remaining challenges, and ultimately transforming how we anticipate and respond to the impacts of extreme weather on power systems.

Lee, Sangkeun (Matt) [ORNL] (ORCID:000000021317511

Economic Evaluation of Modernization Expenditures for Electric Utility Distribution Systems: A Guide for Utility Regulators

Cost-effectiveness evaluation of potential grid modernization solutions is integral to integrated distribution system planning (IDSP). This report synthesizes and updates prior cost-effectiveness approaches identified by the U.S. Department of Energy for grid modernization investments, incorporating multi-objective decision-making. The structure of the report follows the IDSP process, from identifying objectives and priorities to prioritizing grid solutions. The report details two methods for conducting initial cost-effectiveness screening of individual and interdependent grid components: (1) Lowest Reasonable Cost and (2) Benefit-Cost Analysis. Drawing on solutions that passed the initial cost-effectiveness screen, the utility prioritizes grid solutions based on policy objectives and priorities, regulatory compliance, operational efficiencies, and other factors to develop a portfolio of distribution solutions. The utility uses multi-objective decision-making to assess each proposed expenditure against each objective and respective metric. The utility applies a weighting factor, reflecting the priority ranking of the objective, to the total numerical “score” based on the contribution of the proposed solution to addressing each objective. Ultimately, the utility ranks each expenditure from highest to lowest priority using the final score for each solution as well as its cost. The report includes examples of emerging best practices for states and utilities and a cost-effectiveness process checklist.

24 POWER TRANSMISSION AND DISTRIBUTION

Improved Weld Residual Stress Modeling System in BlackBear

This report presents enhancements to the MOOSE-based BlackBear application aimed at improving its capability to simulate welding and other thermo-mechanical manufacturing processes. Two primary avenues of improvement are pursued. First, to enhance user accessibility, we introduce a centralized default block restriction mechanism that ensures coverage checks are performed within user-specified default blocks. This default setting is applied consistently to all block-describable objects, such as variables, kernels, and more. In addition, we develop a modular action for moving heat source simulations, which integrates path file parsing, subdomain modification, and heat source kernel enforcement into a single, streamlined configuration. Second, to improve solver robustness, we implement an alternative method for assigning initial conditions to the updated active domain during the simulation, thereby enhancing convergence behavior. To validate the framework, we design and conduct several benchmark simulations, including heat conduction with progressive material addition, linear elasticity with time-dependent material deposition, and viscoplasticity model with isotropic hardening under similar conditions. Finally, we demonstrate the effectiveness of the proposed framework through large-scale thermo-mechanical welding simulations in both two and three dimensions.

42 ENGINEERING

Direct Processing of Lithium Chloride Based Waste Salt in a Ceramic Waste Form

Ceramic waste form materials were fabricated with a LiCl based salt mixture to demonstrate the suitability of direct processing of waste salt from an all LiCl pyroprocessing flow sheet on the formation, microstructure, and durability of ceramic waste forms. Materials were successfully synthesized by using the simplified direct processing technique with salt loadings of 5, 7.5, and 10 wt% at the laboratory scale. Materials fabricated with LiCl based salt indicated full conversion of zeolite to sodalite occurred to form products with low open porosity. This indicates that the Na 2 O content of NBS4 glass is suitable for processing salt chemistries which do not contain NaCl. Generated sodalite was microencapsulated by the glass binder, halite occlusions were detected within sodalite domains at all salt loadings, and a Cs-rich phase believed to be Cs-pollucite was detected in the material made with the lowest salt loading (5 wt%). Salt inclusions were encapsulated within the binder glass at all loadings, and found to consist primarily of NaCl with Cs present. A small amount of salt was found at the surface of the wasteform made with the highest salt loading of 10 wt%, which may indicate an upper limit to the amount of salt that can be accommodated or incomplete mechanical mixing of the reagents. Durability testing by using the ASTM C1308 method indicated that initial rapid dissolution of exposed halite inclusion phases was directly correlated to salt loading. Wasteform degradation during all tests occurred primarily due to dissolution of the sodalite phase. Durability of the binder glass was found to increase as the waste salt loading increased. Results for the release of cations from the salt indicate the changes in wasteform durability and halite inclusion chemistry are strongly influenced by ion exchange phenomena between the salt and the glass. This appears to be the result of the relatively stronger affinity for the glass phase of Li, compared to that of other cations such as Na and Cs. The effects from cation exchange were diminished at lower salt loadings.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W

FY26 MID-YEAR REPORT 238 Pu Measurement in Bulk Environmental Samples by Thermal Ionization Mass Spectrometry (TIMS)

Funding was received in December 2026 to optimize a nascent measurement capability for mass spectrometry of 238 Pu, as it would apply to relevant environmental collections in safeguards applications, i.e. NWAL swipes or other bulk collections such as soil, sediment, waters, or air filters. The measurement relies on differing ionization temperatures for U and Pu during thermal ionization and therefore leverages TIMS to partially separate and measure these elements. Collected counts of 238 Pu can be mathematically corrected for interfering counts of background 238 U, by addition of a U tracer (e.g. 235 U) to account for the 238 U present. The method has undergone initial characterization at LANL regarding detection limit, precision, and accuracy but requires further applicability-testing for possible use in safeguards.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Synthesizing, Compounding, and Characterizing a Heat Labile Polyurethane Foam

ABSTRACT A need exists for a packaging foam material that can be converted from solid to gaseous degradation products at reasonably low energy levels or temperatures, such as 100O C. This paper will primarily discuss the approaches currently being used to synthesize and characterize such a material. These approaches include the incorporation of novel polyols such as azo containing diols, polycarbonate diols, and polypropylene carbonate polyols into polyurethane foams. Characterization methods include NMR, FTIR, TGA, finite element analysis, impact strength, and others. This project will be funded for a duration of three years. Year 1 focused on developing the proposed test methods and producing an initial rigid polyurethane foam. Year 2 focuses on refining the materials and test methods. If appropriate, design of experiment (doe) techniques will be used to optimize components, component levels, density and other variables to attain required final material properties (TGA weight loss, impact strength, etc.). Year 3 focuses on scaling up to larger engineering quantities. The application for this material is in load securement for transportation of low level radioactive waste materials within the US Department of Energy (DOE) complex. Foam in place process equipment and operators will be shielded using this novel method over current practice. Current practice can involve time consuming methods of load securement in low level radiation environments. This new technique would eliminate exposure time securing the load and greatly improve the As Low As Reasonably Attainable (ALARA) conditions. The objective is to progress to higher Technical Readiness Levels (TRL) and larger pilot scale quantities. This paper discusses methodologies and presents current results to date.

Kranjc, Mark D.

Distribution Feeder Characteristics and Their Resiliency to Natural Hazards

This paper introduces a method and approach for initially screening actions to improve the reliability of distribution systems exposed to three natural hazards (wildlife, weather, and vegetation). The method is well-suited to needs and capabilities of smaller utilities because it relies on readily available information and does not depend on simulations involving detailed engineering models. The method involves first aggregating feeders into like types and then correlating them, separately by hazard, to standard measures of reliability, namely, System Average Interruption Frequency Index (SAIFI) and System Average Interruption Duration Index (SAIDI), and the individual constituents of these measures. The correlations provide both new insights into the efficacy of accepted reliability management practices while and re-confirm long-understood, accepted practices. Importantly, the method does not seek to replace the need for detailed engineering analysis. Instead, in view of the significant costs involved in conducting these more involved analyses, the method is intended to help smaller utilities prioritize their more limited resources to maximize the efficacy of actions they take to improve reliability.

24 POWER TRANSMISSION AND DISTRIBUTION

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING

Gallium oxide (Ga2O3) as a radiation detector distinguishing neutrons and gammas

Gallium oxide (Ga2O3) is a promising ultrawide bandgap semiconductor for radiation detection with the potential of integrating electronic and scintillation functions within a single crystal device. This study establishes the scintillation response of β-Ga2O3 gamma irradiation from yttrium-88 (88Y). Then, californium-252 (252Cf) is used as a spontaneous fission source of mixed neutron and gamma radiation field to measure scintillation signals. Pulse shape discrimination and constant fraction discrimination techniques were used to separate neutron and gamma interaction events. Further investigation indicates that the prompt temporal responses of β-Ga2O3 for gammas and neutrons may enable discrimination of the two by prompt pulse fitting methods, focused around the initial peak. For gamma irradiation, we observed a rise time (τr) of 2.1 ns, decay time (τd) of 9.5 ns, and a full width at half maximum (FWHM) of 6.2 ns. For neutrons, it showed a τr of 2.3 ns, a τd of 12.1 ns, 9.4 ns FWHM, and reduced peak intensity. A diamond detector exhibited a more symmetrical τr and τd for both gamma and neutron signals and therefore is less effective at discriminating between the two by this method. This draws attention to β-Ga2O3’s ability to distinguish neutron and gamma particles. These findings showcase Ga2O3’s potential as a next-generation semiconductor for applications in nuclear safety and medical imaging, where precise discrimination between neutron and gamma interactions is essential.

Valdes, D. J. (ORCID:0000000304373720)

Toward engineering lattice structures with the material point method (MPM)

This study examines the potential of two variants of the material point method—the generalized interpolation material point (GIMP) and dual domain material point (DDMP) methods—in developing a robust computational framework for engineering lattice structures under different loading conditions. The study begins with assessing the ability of the two methods in predicting elastic buckling phenomena using column geometries with and without initial geometric imperfections. The results indicate that both methods effectively capture buckling phenomena when initial geometric imperfections are introduced. After this verification step, we create several models of tetrahedral lattice structures with varying strut diameter and orientation and subject them to quasi-static loading. We then validate the numerical results using laboratory test results. The results show that, while both methods accurately predict load–displacement curves in the pre-buckling regime, their predictive capabilities diminish in the post-buckling regime. Through visual comparison between the numerical and experimental deformed shapes, it appears that the discrepancies between model and experimental results are attributed to initial geometric imperfections in the lattices that occurred during 3D printing. We then establish a second set of lattice models where different types of initial geometric imperfections are considered. The results from these models show that imperfections have a negligible influence in the pre-buckling regime but affect the behavior considerably in the post-buckling regime. As a final step in this work, we subject the lattice models to impact loading and employ hypothetical soft and stiff materials. These results show that the lattice stiffness, which depends on material stiffness, strut diameter, and orientation, significantly influences the ability of a lattice structure to resist impact. In particular, we find that a stiffer lattice (i.e., one made with a stiff material and thicker struts) is capable of absorbing more energy than a softer one during impact. Although material nonlinearities, inelasticity, and detailed contact formulations are not considered in this study, the findings obtained herein lay the groundwork for engineering lattice structures under extreme loading conditions through a simulation-driven framework based on particle-based methods.

97 MATHEMATICS AND COMPUTING