Search NASA⌕ Search

SEARCH · Search NASA

Results for “RESIDUAL STRESS”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Stress engineering for crack and dendrite prevention in solid electrolytes via ion implantation

Solid-state batteries represent a promising technology that offers safer and more densely packed energy storage. A primary cause of failure in solid-state cells is the penetration of metal dendrites through the solid electrolyte. Here, we report that fluorine-ion implantation can enhance the mechanical resistance of solid electrolyte Li 6.5 La 3 Zr 1.5 Ta 0.5 O 12 to dendrite propagation by inducing residual compressive stress in the subsurface of the electrolyte. Ion implantation modifies subsurface residual stress and also alters the electrolyte’s surface, resulting in multifunctional enhancement. The combined chemical and mechanical effects of ion implantation enable reversible lithium metal stripping and plating at elevated current densities while enhancing air stability and mitigating the formation of a harmful carbonate layer at the surface. This study provides new insights into a scalable dendrite-suppression strategy for designing solid-state batteries suitable for cycling at room temperature and low stack pressures.

25 ENERGY STORAGE↗
DOE OSTI2024

Stressful crystal histories recorded around melt inclusions in volcanic quartz

Abstract Magma ascent and eruption are driven by a set of internally and externally generated stresses that act upon the magma. We present microstructural maps around melt inclusions in quartz crystals from six large rhyolitic eruptions using synchrotron Laue X-ray microdiffraction to quantify elastic residual strain and stress. We measure plastic strain using average diffraction peak width and lattice misorientation, highlighting dislocations and subgrain boundaries. Quartz crystals across studied magma systems preserve similar and relatively small magnitudes of elastic residual stress (mean 53–135 MPa, median 46–116 MPa) in comparison to the strength of quartz (~ 10 GPa). However, the distribution of strain in the lattice around inclusions varies between samples. We hypothesize that dislocation and twin systems may be established during compaction of crystal-rich magma, which affects the magnitude and distribution of preserved elastic strains. Given the lack of stress-free haloes around faceted inclusions, we conclude that most residual strain and stress was imparted after inclusion faceting. Fragmentation may be one of the final strain events that superimposes stresses of ~ 100 MPa across all studied crystals. Overall, volcanic quartz crystals preserve complex, overprinted deformation textures indicating that quartz crystals have prolonged deformation histories throughout storage, fragmentation, and eruption.

58 GEOSCIENCES↗
DOE OSTI2024

Phase-field modeling of thermally-grown oxide and damage evolution in environmental barrier coatings

Silicon carbide-based ceramic matrix composites protected by environmental barrier coatings (EBCs) present a promising materials solution for next-generation gas turbines. Developing more robust and efficient EBCs is therefore of significant technological importance. During the service in high-temperature oxidative environments, there is a thermally grown oxide (TGO) layer, spontaneously formed in the EBC system. TGO is recognized as a critical factor for the degradation and failure of EBCs, yet the detailed mechanisms of TGO growth and its effect on EBC failure remain unclear. In this study we develop a comprehensive chemo-mechano-phase-field model to simulate growth of the TGO in EBCs, factoring in creep and deformation, and especially the cracking behaviors. The volume expansion due to TGO growth and the resulting large inelastic deformation are addressed by using our recently developed, so-called incremental realization of inelastic deformation (IRID) algorithm, in combination with an adapted Hu-Chen spectral solver for elasticity. Simulations of TGO growth are performed considering different growth modes of TGOs determined mainly by the ratio of oxidant permeability in the topcoat to that in the TGO itself. Large-scale three-dimensional (3D) simulations are performed to model the formation of interconnecting vertical/channel cracks (often called ‘mud cracks’). The simulated crack morphology are in excellent agreement with the experimental observations from the literature. The simulations also provide insights into the cracking of EBCs and its dependence on the structure and constituent properties of the coating system. Furthermore, these results demonstrate the developed damage model can be a useful tool for design of more durable EBCs.

Cracking↗
DOE OSTI2026

Learning interpretable surface elasticity properties from bulk properties via neural network equation learners

Surface elasticity is central to understanding the mechanics and stability of surfaces and interfaces. It is characterized by quantities such as surface tension, residual surface stress, and surface stiffness. However their analytical expressions are typically difficult to derive from atomistic data, and depend strongly on modeling choices. This work presents a neural network-based equation learner which combines customized activation functions and connection-based pruning to discover parsimonious, closed-form equations for surface elasticity from atomistic simulations. Applying the method to seven face-centered cubic (FCC) metals, our equation learner uncovers interpretable equations that describe both low-Miller index and high-Miller index surface properties, capturing long-tail property distributions accurately. The discovered expressions are decoupled into two components: a universal, geometry-driven orientation function, and material-specific baseline coefficients. We find that lower-order properties such as surface tension are fundamentally geometry dependent, while higher-order properties such as surface stress and elasticity show more complex geometry and material dependence. We also relate material dependent coefficients to bulk properties, forming a clear map from bulk material properties to surface elasticity. Overall, this approach demonstrates that interpretable neurosymbolic machine learning can bridge the gap between atomistic simulations and physical laws, enabling the discovery of generalizable structure–property relationships for materials science phenomena such as surface elasticity.

Equation learning↗
DOE OSTI2025

Ion Size Effects on the Thermodynamic, Kinetic, and Mechanical Properties during Ion Exchange in Solid-State Electrolytes

Ion exchange offers a pathway to impose residual compressive stresses to mitigate the electro-chemo-mechanical cracking of solid-state electrolytes such as lithium lanthanum zirconium oxide. This study uses a coupled multiscale framework (integrating density functional theory (DFT), molecular dynamics (MD), and continuum modeling) to examine how exchange ion size influences stress, diffusion, fracture toughness, and electronic properties. Larger isovalent ions (Na + , Ag + , K + ) were exchanged with Li + , with DFT confirming their preference for octahedral 96h sites and a linear relationship between ion size and chemical free expansion coefficient. MD simulations reveal stress and concentration effects on exchange ion diffusivity at elevated temperatures, with Na + and Ag + maintaining favorable mobility while K + showing concentration-dependent clustering. Continuum modeling predicts the range of fracture strength improvements and the required ion exchange concentration profile. It was shown that a 5% surface exchange concentration can induce ∼0.6 GPa of surface compressive stress using Na + and ∼1.0 GPa of surface compressive stress using Ag + . On the other hand, larger ion exchange species may penalize Li + diffusivity by increasing the activation volume and activation energy. Interestingly, Na + has a negligible penalty on Li-ion diffusivity. The room temperature Li + ion diffusivity is reduced by ∼40% with Ag + ion exchange. Electronic band structure analysis shows no size-dependent change in the bandgap, though Ag + introduces localized defect states near the valence band maximum. This study highlights ion size as a key factor in optimizing LLZO properties, offering a framework to improve the solid-state battery performance.

Jagad, Harsh D. [Brown Univ., Providence, RI (Unit↗
DOE OSTI2025

Material Needs and Measurement Challenges for Advanced Semiconductor Packaging: Understanding the Soft Side of Science

This Perspective builds upon insights from the National Institute of Standards and Technology (NIST)-organized workshop, “Materials and Metrology Needs for Advanced Semiconductor Packaging Strategies,” held at the 35th annual Electronics Packaging Symposium in Binghamton, NY, on September 5, 2024. It outlines critical challenges and opportunities related to polymer-based “soft” materials in advanced semiconductor packaging, with emphasis on polymer science, measurement science (metrology), and the strategic development of Research-Grade Test Materials (RGTMs). These efforts, led by the NIST CHIPS team, aim to advance the fundamental understanding of structure-property-processing relationships, promote standardized guidelines and innovative methods for material characterization, and accelerate the development, qualification, and adoption of next-generation packaging materials. The Perspective also distills key insights from the panel discussion with industry experts, emphasizing the need for close collaboration among materials scientists, process engineers, and metrology experts to enable a holistic strategy, further highlighting the importance of cross-sector partnerships among industry, academia, and government to address pressing challenges in packaging materials and processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Evaluation of In-Situ AM Process Monitoring Techniques and Potential for Detecting Process Anomalies and Undesirable Microstructures

The US Department of Energy’s Advanced Materials and Manufacturing Technologies (AMMT) program is pursuing rapid qualification of new materials for fabrication of nuclear relevant components using advanced manufacturing techniques. Particular interest is placed on code-qualifying stainless steel (SS) 316H processed by laser powder bed fusion (LPBF) additive manufacturing. A paradigm that incorporates data from in-situ sensing during the printing, ex-situ characterization, and advanced artificial intelligence–based models was established under the Transformation Challenge Reactor (TCR) program to develop a pedigree for each fabricated component that could be tracked from the feedstock to the component’s release for application. Under the TCR program, the Peregrine software was developed as a tool for incorporating the vast amounts of in-situ and ex-situ characterization data collected; all data stored on a rapidly growing digital platform. The digital platform allows for users to link site-specific process anomalies to the macro- and microstructure. The platform will eventually be able to predict component performance, which will be crucial to qualifying materials and components in risk-averse industries such as those supporting and building nuclear reactors. Current in-situ process monitoring techniques that are already integrated with software like Peregrine are advantageous for identifying process anomalies including powder spatter, component edge swelling, recoating-build interactions, and so on. However, additional data are required to fully predict the resulting microstructures needed for identifying relationships to component performance. The rapid cooling rates observed in LPBF are some of the highest of any bulk manufacturing process, resulting in heterogenous microstructures and typically causing anisotropy in mechanical properties. Moreover, evolved residual thermal stresses are high, which can cause severe defects such as delamination or cracking. Therefore, other in-situ monitoring methods are warranted for exploration to measure and map the thermal history, and potentially the stress state, of each build. This report summarizes different in-situ monitoring strategies proposed for LPBF with a focus on the more developed sensor systems. Novel capabilities for measuring melt pool temperatures are also addressed to better inform modeling efforts.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Analysis of Molten Salt Thermal Energy Storage Tank Sizing and Preheating

Molten Salt Thermal energy storage (TES) is essential for concentrating solar power (CSP) plants, enhancing their capacity factor and dispatchability to ultimately reduce electricity costs. This technology is becoming more important for other technologies including nuclear power, thermal-electrical storage, and concentrating solar thermal (CST). While molten salt TES tanks have been successfully deployed in lower-temperature parabolic-trough plants, the current state-of-the-art central receiver systems operating at temperatures up to 565 degrees Celsius have experienced costly and frequent failures. These failures are associated with the immaturity of the technology and problems in tank design, fabrication, commissioning and, sometimes, aggressive operation. A primary contributing factor is the absence of a dedicated technical standard for high-temperature molten nitrate TES tanks; designers must currently rely on existing codes like API 650 and ASME Section II, which are insufficient because they fail to account for the unique demands of high temperatures, thermal cycling, and transient conditions.

25 ENERGY STORAGE↗
DOE OSTI2024

Environmentally Assisted Fatigue in Light Water Reactor Environment

This report summarizes the Environmentally Assisted Fatigue (EAF) research conducted at ANL under the US DOE Light Water Reactor Sustainability (LWRS) program. Starting from a rich background in theoretical and experimental EAF, ANL previously developed an approach to evaluate fatigue performance of reactor materials in light water reactor environments with the correction factor F en . The approach was based on a large body of experimental work performed at ANL and elsewhere, and was consistent with American Society of Mechanical Engineers (ASME)’s methodology governing the design and construction of reactor components. In recent years, the program was focused on component fatigue prediction and made several major and fundamental contributions in this area. These accomplishments help meet the needs identified by the industry concerning component level fatigue predictions in complex, transient conditions. The main contribution of the ANL program involved the development of a system-level model for estimating residual strain and life of nuclear reactor coolant system components under connected-system-thermal-mechanical boundary conditions. The goal was to predict the stress hotspots, strain residuals, strain amplitudes and the resulting fatigue lives. Thermal-mechanical stress analysis was performed considering thermal stratification and a design-basis reactor loading cycle. Based on the finite element (FE) model results, the strain residuals, strain amplitudes and resulting fatigue lives of reactor coolant system (RCS) components were predicted. The results show that some of the RCS components can have significantly different strain amplitudes, residual strain, and fatigue lives, despite having similar geometry and material. In addition, the simulated component-level strain profile can guide the selection of appropriate test inputs for conducting laboratory-scale EAF tests. Building upon the system-level model, ANL developed a digital twin (DT) framework to predict the structural states and associated fatigue life of components in real-time. This framework is a comprehensive system designed to predict the structural states and fatigue lives of reactor components. It includes multiple models and integrates artificial intelligence (AI), machine learning (ML), and FE based modeling tools to evaluate the structural states and fatigue lives.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Uncovering hidden bias in neutron diffraction residual strain measurements

When calculating residual strain via neutron or X-ray diffraction, uncertainties propagated from the peak fit are often inadequate to describe the true scatter of measurements about a singular strain state, such as one that should describe a macroscopic continuum. Because diffraction is inherently a selective process, orientation-dependent scatter arises from the sub-sampling of strong microstructure and strain gradients. This paper investigates the appropriateness of propagated uncertainties with reference to their original intention, i.e. noise about a mean value. Thirty-six unique orientations of strain measurements are taken at multiple locations within an additive friction-stir deposition component with fine-scale gradients (∼200 µm) of plastic strain, texture and residual elastic strain. Multiple strain and stress calculation pathways are compared: direct substitution of three measurements into Hooke's law, direct inversion of any six unique orientations into the strain state tensor and thirty-six measurement least-squares estimation. For the last two cases, the appropriateness of the uncertainty interval is statistically evaluated on the basis of a physical constraint: common agreement under the strain transformation law. For this sample, the direct inversion of six measurements retains a conservative estimate of the uncertainty. However, propagated uncertainties in the least-squares solution greatly underestimate the true experimental scatter. A simple pathway to estimate appropriate uncertainty intervals is suggested. These results demonstrate that the interpretation of uncertainty in residual strain is strongly dependent on intrinsic sample-dependent effects, and that oversampling orientations and statistical analysis can give more accurate results with realistic uncertainties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Characterization of kerogen nanopores using 2D NMR relaxation and MD simulations

The characterization of kerogen nanopores is crucial for predicting the geostorage capacity and recoverability of natural gas in unconventional gas shale reservoirs. Towards this end, a powerful technique is presented which integrates 2D NMR T 1 -T 2 relaxation measurements with molecular dynamics (MD) simulations of hydrocarbons confined in the nanopores of kerogen. The integrated NMR-MD technique is demonstrated using T 1 -T 2 measurements of kerogen isolates and organic-rich chalks saturated with heptane, together with MD simulations of heptane completely dissolved in a realistic kerogen model. The NMR-MD results are used to extract the swelling ratio and nanopore size distribution of kerogen as a function of depth in the reservoir. The effects of organic nanoconfinement on the T 1 relaxation dispersion and T 2 residual dipolar coupling of heptane are investigated, as well as the effect of downhole effective stress on the kerogen nanopore size as a function of depth and compaction. Potential applications in partially depleted gas shale reservoirs are discussed, including CO 2 utilization/geostorage, geostorage of green H 2 , and integration of the NMR-MD technique with thermodynamic models for predicting the competitive sorption of gas mixtures in kerogen.

Compaction↗
DOE OSTI2024

Informative and non-informative decomposition of turbulent flow fields

Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative Φ I (x, t) and residual Φ R (x, t) components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in Φ I (x, t), whereas Φ R (x, t) offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Attribution of heterogeneous stress distributions in low-grain polycrystals under conditions leading to damage

In high-purity polycrystalline metallic materials, voids tend to favor grain boundaries as nucleation sites due to the elevated stress states produced by granular interactions and the weakened grain boundary from the relative atomic disorder. To quantify the key factors of this elevated stress state, simple compression of a small multi-grain cylinder of body-centered cubic tantalum was simulated using a single crystal plasticity model that incorporates non-Schmid effects. Four increasingly complex synthetic microstructures were created to tractably incorporate grain boundary interactions, and a statistically significant number of combinations were performed by varying the initial crystallographic orientations of the microstructure. Most of these simulations produce the maximum von Mises stress on a grain boundary and less frequently at the multi-grain junctions. To build a statistical model for the maximum von Mises stress at the grain boundary, physically based features that could contribute to the elevated stress state were selected. Then, a learning algorithm based on information theory was used to identify which of these features contributed the most information to the data set. The identified features include a grain’s propensity to accommodate both elastic and plastic deformations and their directional components. The misalignment of the direction of each grain’s mechanical response was found to be strongly correlated to the magnitude of the stress near the grain boundary. For all of the synthetic microstructures, the statistical models produce a residual distribution that is nearly Gaussian with a variance of, at most, 10% of the prior distribution. The successful performance of the statistical model implies the correct identification of the physical features that cause severe stress localization in polycrystalline materials. The statistical models constructed here can be used to formulate a physically motivated void nucleation model which is sensitive to a microstructure’s propensity to produce elevated stress states. As a result, these statistical models also enable the design of material microstructures, in which the crystallographic orientation is chosen to resist void nucleation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Glass production rate in an electric melter: Melting rate correlation and primary foam stability

Waste loading and melting rate are major factors determining the lifecycle of nuclear waste disposal by vitrification in electric melters in which the heat to the cold cap is delivered from the melt pool. One of the crucial parameters that affect the rate of melting is the temperature at which the primary foam collapses at the cold cap bottom. Further, apart from the transient glass-forming melt viscosity, primary foam stability is affected by the shear stress imposed by the flow in the melt pool and the presence of residual solid particles.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Utah FORGE: Triaxial Direct Shear Results - February 2025

This dataset contains results from nine triaxial direct shear tests conducted by Los Alamos National Laboratory on samples from FORGE Well 16A(78)-32. The primary objectives of this work were to determine the shear strength in both intact and residual states, evaluate dilation against displacement, assess permeability in relation to displacement, time, and normal stress, understand the relationship between aperture and normal stress, and monitor the effluent chemistry as a function of time. The data includes time-series measurements of stress, displacement, permeability, and effluent chemistry, with and without experimental dilution corrections. Additional materials include profilometry data, photographic documentation of the experimental setups and apparatus, and test notes. The dataset is organized into folders corresponding to each test, containing hydromechanical data, effluent chemistry measurements, and images. The hydromechanical data consists of detailed time-series records capturing parameters such as shear force, confining pressure, permeability, and temperature. Effluent chemistry data tracks fluid composition changes over time. Also included are conference papers, presentation slides, and a summary document outlining the experiments.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Neural network potentials with effective charge separation for non-equilibrium dynamics of ionic solids: a ZnO case study

Developing neural network potentials (NNPs) accurate under non-equilibrium dynamics is challenging, as such systems require extensive sampling beyond equilibrium phases. Here we construct high-fidelity NNPs for zinc oxide (ZnO), a polymorphic ionic solid, using density functional theory (DFT) reference data. To efficiently capture transitional configurations, we combine enhanced-sampling molecular dynamics with empirical potentials, data distillation, and pretraining on short-range atomic energies (A-Train), followed by transfer learning with DFT-relabeled datasets. This hierarchical approach improves transferability across polymorphs and stress states. We further introduce effective charge separation, treating long-range Coulombic terms analytically while short-range residual interactions are learned by the NNP. The optimal effective charges fall in the range 0.5–1.0 q e , consistent with dielectric-screened values derived from formal charges but distinct from Bader estimates. Motivated by this observation, we propose a simple data-driven protocol in which effective charges are optimized by comparing DFT reference energies with explicit Coulomb calculations, without additional NNP training. This strategy improves accuracy and transferability in DFT-level predictions of energies, forces, and stress. Together, these results provide a practical charge-selection framework for robust NNP development in ionic solids, enabling reliable simulation of polymorphic phase transformations and non-equilibrium dynamics.

Chemistry↗