Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inelasticity”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A physics-informed multi-agents model to predict thermo-oxidative/hydrolytic aging of elastomers

This paper introduces a novel physics-informed multi-agents constitutive model to propose prediction in quasi-static constitutive behavior of cross-linked elastomer and the loss of mechanical performance during environmental aging. The presented model is used to simulate the effect of single-mechanism chemical aging (i.e. thermal-inducedor hydrolytic aging) on the behavior of the material in this hybrid framework. Those environmental single-mechanism damages change the polymer matrix over time due to massive chain scission, chain formations, and changing the arrangement of molecules in the polymer matrix. Here we propose a data-driven super-constrained machine-learned engine to represent damage in the polymer matrix and capture the changes in material behavior, including its inelastic features such as Mullins effect and permanent set in the course of aging. We have simplified the 3D stress–strain tensor mapping problem into a small number of super-constrained 1D mapping problems by means of a sequential order reduction. An assembly of multiple replicated conditional neural-network learning-agents (L-agents) is trained to systematically simplify the high-dimensional mapping problem into multiple 1D problems, each represented by a different type of agent. Our hybrid framework is designed to capture the effect of deformation history, aging time, and aging temperature. The model is validated with respect to a comprehensive set of experiments specifically designed to benchmark model capabilities and also against available data in the literature. Thermodynamic consistency and frame independency have been verified. Besides acceptable predictive abilities, a significant reduction of computational cost to predict behavior at multiple states of deformation is the most significant feature of this model.

42 ENGINEERING↗
DOE OSTI2023

Training material models using gradient descent algorithms

High temperature design requires accurate constitutive models to describe material inelastic deformation and failure behavior. Oftentimes, calibrating accurate models devolves into the problem of fitting the model parameters against experimental test data. Here, we present the pyopmat package, an open source framework for calibrating constitutive models against experiment data subjected to various loading conditions using machine learning techniques. The package calculates the exact gradient of the model response with respect to the parameters using a combination of automatic differentiation and the adjoint method. Given this exact gradient, we compare the performance of several gradient-based optimization techniques in fitting realistic constitutive models against data. Here, we demonstrate the efficiency and accuracy of our package through example problems using both synthetic data, generated using known parameter sets, under monotonic and cyclic loading conditions and also with an example applying the techniques developed here to actual high temperature creep-fatigue test data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗
DOE OSTI2020

Activation and surface reactions of CO and H 2 on ZnO powders and nanoplates under CO hydrogenation reaction conditions

Here, activation and surface reactions of CO and H 2 on ZnO powders and nanoplates under CO hydrogenation reaction conditions were (quasi) in situ studied using temperature programmed surface reaction spectra, diffuse reflectance Fourier transform infrared spectroscopy, inelastic neutron scattering spectroscopy and electron paramagnetic resonance. CO undergoes disproportion reaction to produce gaseous CO 2 and surface carbon adatoms, and adsorbs to form surface formate species. H 2 adsorption forms dominant irreversibly-adsorbed surface hydroxyl groups and interstitial H species and very minor surface Zn-H species. Surface formate species and hydroxyl groups react to produce CO 2 and H 2 , while surface carbon adatoms are hydrogenated by surface Zn-H species sequentially to produce CH(a), CH 2 (a), CH 3 (a) and eventually gaseous CH 4 . The ZnO nanoplates, exposing a higher fraction of Zn-ZnO(0001) and O-ZnO(000–1) polar facets, are more active than the ZnO powders to catalyze CO hydrogenation to CH 4 . These results provide fundamental understanding of the reaction mechanisms and structural effects of CO hydrogenation reaction catalyzed by ZnO-based catalysts.

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

High resolution ferromagnetic and antiferromagnetic correlations in La 0.7 Ca 0.3 MnO 3

La 0.7 Ca 0.3 MnO 3 (LCMO30) is a perovskite manganite with a phase transition (T C = 238 K) between a paramagnetic insulator and a ferromagnetic metal. We use single-crystal, inelastic neutron scattering to refine previous work on the low energy (ℏω < 10 meV) magnetic excitations. At 300 K the spectrum consists of a diffusive ferromagnetic component at the pseudo-cubic Brillouin zone center coexisting with a quasielastic excitation located at the pseudo-cubic Brillouin zone edge. Both components are anti-correlated with ferromagnetic long-range order, showing weakened intensity at 215 K in the ferromagnetic phase. Polarized neutron scattering shows both features to be purely magnetic in origin. Furthermore, this indicates the simultaneous presence of ferromagnetic and antiferromagnetic correlations (located, respectively, at the zone center and zone edge) in a three-dimensional perovskite, while this coexistence had been previously observed only in a bilayer compound.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Coupled chemo-mechanical modeling of point-defect diffusion in a crystal plasticity fast Fourier transform framework

Below the yield strength and at moderate-to-high homologous temperatures, the inelastic deformation of metals is mostly governed/rate-controlled by vacancy diffusion-mediated processes. As a function of grain size, stress, temperature and dislocation content, vacancies (or atoms) can adopt preferential diffusion paths across grain interiors, along grain boundaries, or towards and along dislocations, resulting in climb and self-climb. In the steady state and under constant load, grain boundary and grain bulk vacancy diffusion-mediated plasticity have been described in seminal works by Coble and by Nabarro and Herring, respectively. Yet, the interplay between all aforementioned potential diffusion pathways has not been comprehensively mapped. This work presents a thermodynamically-consistent full-field model integrated within a voxel-based elasto-viscoplastic fast Fourier transform framework, which considers the coupling between the diffusion-mediated plasticity mechanisms. In the proposed approach, the kinetics and kinematics of plastic deformation due to vacancy diffusion along grain boundaries and grain bulk, as well as the exchange between grain boundaries and bulk are described explicitly. A homogenization approach at the voxel level is further introduced to simultaneously consider bulk and grain boundary diffusion in a numerically efficient fashion. The new formulation predicts the expected strain rate dependencies and the scaling of the steady-state creep rate with respect to grain size, temperature, and stress. Finally, the model predicts the transition from grain bulk to grain boundary-dominated diffusion with reduction in grain size, a significant step towards capturing transitions in deformation behavior without any phenomenological or ad-hoc adjustments.

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

Spherical Indentation and Implementation of S 3 /P for yield stress determination of brittle materials

Here, a mathematically transparent and robust experimental method has been developed to estimate the yield stress of brittle materials through the analysis of depth-sensing spherical indentation. Employing Hertzian contact mechanics, an elastically invariant ratio based on the simple equation S 3 / P = 6RE r 2 , (where S and P are contact stiffness and indentation load, respectively) has been derived that enables more accurate and confident determination of the transition from elastic to inelastic deformation; a transition that the yield stress dictates and represents. Using two diamond spheres with radii of 3.2 and 8.6 μm, the indentation test method and analyses are applied to two vitreous silicates: Corning's HPFS 7980® fused silica and Vitro's Starphire® soda lime silicate. The estimated yield strengths are 8.15 GPa ± 2.5% for the fused silica and 6.1 GPa ± 3.3% for the soda-lime silicate, and both were independent of indenter radius. Verification of this new experimental method is demonstrated with an as-drawn titanium by showing equivalence of measured yield stress by its spherical indentation and that from uniaxial compression testing. This method will enable easier and more confident estimation of yield stress in brittle materials - a property that historically has been elusive to measure for these materials using common laboratory mechanical test methods.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electronic energy loss and ion velocity correlation effects in track production in swift-ion-irradiated LiNbO3: A quantitative assessment between structural damage morphology and energy deposition

The primary motivation for studying how irradiation modifies the structures and properties of solid materials involves the understanding of undesirable phenomena, including irradiation-induced degradation of components in nuclear reactors and space exploration, and beneficial applications, including material performance tailoring through ion beam modification and defect engineering. In this work, the formation mechanism of latent tracks with different damage morphologies in LiNbO 3 crystals under 0.09–6.17 MeV/u ion irradiation with an electronic energy loss from 2.6–13.2 keV/nm is analyzed by experimental characterizations and numerical calculations. Irradiation-induced damage is preliminarily evaluated via the prism coupling technique to analyze the correlation between the dark-mode spectra and energy loss profiles of irradiated regions. Under the irradiation conditions of different ion velocities and electronic energy losses, different damage morphologies, from individual spherical defects to discontinuous and continuous tracks, are experimentally characterized. During ion penetration process, the ion velocity determines the spatiotemporal distribution of deposited irradiation energy induced by electronic energy loss, meaning that the two essential factors including electronic energy loss and ion velocity co-affect the track damage. The inelastic thermal spike model is used to numerically calculate the spatiotemporal evolutions of energy deposition and the corresponding atomic temperature under different irradiation conditions, and a quantitative relationship is proposed by comparison with corresponding experimentally observed track damage morphologies. Additionally, the obtained quantitative relationship between irradiation conditions and track damage provides deep insight and guidance for understanding the damage behavior of crystal materials in extreme radiation environments and selecting irradiation parameters, including ion species and energies, for ion beam technique application in atomic-level defect manipulation, material modification, and micro/nanofabrication.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A foundation model for non-destructive defect identification from vibrational spectra

Defects are ubiquitous in solids and strongly influence materials’ functional properties. However, non-destructive characterization and quantification of defects, especially when multiple types coexist, remain a long-standing challenge. Here, we introduce DefectNet, a foundation machine learning model that predicts the chemical identity and concentration of substitutional point defects with multiple coexisting elements directly from vibrational spectra, specifically phonon density-of-states (PDoS). Trained on over 16,000 simulated spectra from 2,000 semiconductors, DefectNet employs a tailored attention mechanism to identify up to six distinct defect elements at concentrations ranging from 0.2% to 25%. The model generalizes well to unseen crystals across 56 elements and can be fine-tuned on experimental data. Validation using inelastic scattering measurements of SiGe alloys and MgB 2 superconductor demonstrates its accuracy and transferability. Furthermore, our work establishes vibrational spectroscopy as a viable, non-destructive probe for bulk point defect quantification, and highlights the promise of foundation models in data-driven defect engineering.

artificial intelligence↗
DOE OSTI2020

Achieving 5.9% elastic strain in kilograms of metallic glasses: Nanoscopic strain engineering goes macro

The ideal elastic limit is the upper bound of the achievable strength and elastic strain of solids. However, the elastic strains that bulk materials can sustain are usually below 2%, due to the localization of inelastic deformations at the lattice scale. In this study, we achieved >5% elastic strain in bulk quantity of metallic glass, by exploiting the more uniform and smaller-magnitude atomic-scale lattice strains of martensitic transformation as a loading medium in a bulk metallic nanocomposite. The self-limiting nature of martensitic transformation helps to prevent lattice strain transfer that leads to the localization of deformation and damage. This lattice strain egalitarian strategy enables bulk metallic materials in kilogram-quantity to achieve near-ideal elastic limit. This concept is verified in a model in situ bulk amorphous (TiNiFe)-nanocrystalline (TiNi(Fe)) composite, in which the TiNiFe amorphous matrix exhibits a maximum tensile elastic strain of similar to 5.9%, which approaches its theoretical elastic limit. As a result, the model bulk composite possesses a large recoverable strain of similar to 7%, a maximum tensile strength of above 2 GPa, and a large elastic resilience of similar to 79.4 MJ/m 3 . The recoverable strain and elastic resilience are unmatched by known high strength bulk metallic materials. This design concept opens new opportunities for the development of high-performance bulk materials and elastic strain engineering of the physiochemical properties of glasses.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On the effect of strain rate during the cyclic compressive loading of liquid crystal elastomers and their 3D printed lattices

Nematic liquid crystal elastomers (LCEs) are a unique class of network polymers with the potential for enhanced mechanical energy absorption and dissipation capacity over conventional network polymers because they exhibit both conventional viscoelastic behavior and soft-elastic behavior (nematic director changes under shear loading). This additional inelastic mechanism makes them appealing as candidate damping materials in a variety of applications from vibration to impact. The lattice structures made from the LCEs provide further mechanical energy absorption and dissipation capacity associated with packing out the porosity under compressive loading. Understanding the extent of mechanical energy absorption, which is the work per unit mass (or volume) absorbed during loading, versus dissipation, which is the work per unit mass (or volume) dissipated during a loading cycle, requires measurement of both loading and unloading response. Here, in this study, a bench-top linear actuator was employed to characterize the loading-unloading compressive response of polydomain and monodomain LCE polymers and polydomain LCE lattice structures with two different porosities (nominally, 62% and 85%) at both low and intermediate strain rates at room temperature. As a reference material, a bisphenol-A (BPA) polymer with a similar glass transition temperature (9 °C) as the nematic LCE (4 °C) was also characterized at the same conditions for comparing to the LCE polymers. Based on the loading-unloading stress-strain curves, the energy absorption and dissipation for each material at different strain rates (0.001, 0.1, 1, 10 and 90 s -1 ) were calculated with considerations of maximum stress and material mass/density. The strain-rate effect on the mechanical response and energy absorption and dissipation behaviors was determined. The energy dissipation ratio was also calculated from the resultant loading and unloading stress-strain curves. All five materials showed significant but different strain rate effects on energy dissipation ratio. The solid LCE and BPA materials showed greater energy dissipation capabilities at both low (0.001 s -1 ) and high (above 1 s -1 ) strain rates, but not at the strain rates in between. The polydomain LCE lattice structure showed superior energy dissipation performance compared with the solid polymers especially at high strain rates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A finite-strain rate- and pressure-dependent constitutive framework for analyzing shock compression behavior of cemented tungsten carbides to 100 GPa

In the present study a thermodynamically-consistent finite-strain rate-and-pressure-dependent constitutive framework is implemented to analyze the shock-compression behavior of cemented tungsten carbides to 100 GPa. Central to this framework is the use of logarithmic strain with a set of invariant basis that allow the Cauchy stress tensor to be expressed as a sum of three response terms that are mutually orthogonal, thus permitting a complete separation of the deviatoric and volumetric (pressure) response. An overstress viscoplasticity model that includes strain and strain rate hardening along with thermal softening is used to represent the deviatoric response, while a complete Mie-Grüneisen equation of state (EoS) is used to obtain the pressure response. Using this formulation, the shock-induced compression behavior of cemented tungsten carbide - obtained from planar plate impact experiments using a 30 mm powder gun to peak stresses of up to ~100 GPa - is analyzed to better understand the structure of the measured shock wave profiles and the associated in-material shock quantities. Of particular interest is the evolution of material inelasticity and strength, and temperature in the tungsten carbide samples during the shock compression process.

Cemented tungsten carbide↗
DOE OSTI2019

In-situ, microscale characterization of heterogeneous deformation around notch in martensitic Shape Memory Alloy

Deformation characterization of low-symmetry phases such as martensite in SMAs is difficult due to a fine-scale hierarchical microstructure. With X-ray microLaue diffraction, in-situ deformation of martensite is examined in a notched NiTi specimen. The local deformation is influenced by the notch stress field, initial martensite microstructure, and interaction between notch stress field and the external load. The microstructure evolves heterogeneously and inelastically, with detwinning and twin nucleation occurring simultaneously through loading. These results contrast with the traditional view of martensite deformation that is partitioned in three distinct regimes: elasticity, reorientation and de-twinning.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Room temperature sodium diffusion in sodium iron phosphate investigated with neutron scattering

Sodium iron phosphate, Na x FePO 4 , is currently one of the most promising cathode materials for sodium ion batteries. Here we employ a systematic approach using three neutron scattering techniques to investigate temperature dependent properties with a focus on the charge transfer processes. Two distinct Na + diffusive motions are observed with quasielastic neutron scattering, including room temperature motion. These processes correspond to rapid double jumps through extrinsic vacancy sites of ∼6 Å along with the more fundamental nearest neighbor jumps of 3.1 Å. Experimental observations of room temperature ionic diffusion are highly desirable, as the diffusion pathways are naturally well suited for efficient transport. Inelastic neutron scattering was used to probe the vibrational density of states that assist in the ionic and electronic (polaronic) transport processes, which are further supported with density functional theory. Neutron diffraction examines the disordering of the sodium sublattice upon transition to the solid-solution phase, which is believed to enhance the overall kinetics during sodium insertion/extraction. Finally, 57 Fe Mössbauer spectroscopy measurements probes polaron hopping and the local structural evolution that modifies the electron density surrounding the redox active iron sites. This study sheds light on the temperature induced atomic scale dynamics that occur in olivine Na x FePO 4 .

Novak, Eric [Swarthmore College, Swarthmore, PA (U↗
DOE OSTI2021

Effects of recoil spectra and electronic energy dissipation on defect survival in 3C-SiC

The coincidence of electronic and damage energy dissipation from energetic ions to an atomic lattice can significantly affect damage production along the ion trajectory due to spatial overlap of inelastic and elastic processes. Damage production and disordering in single crystal 3C-SiC from 5 MeV Si and 10 MeV Au ions is investigated using ion-channeling experiments. While defects are created by damage energy dissipation via elastic scattering, electronic energy dissipation via electron-phonon coupling decreases defect survival along the ion trajectory for Si ions. Finally, the more energetic recoil spectrum for 10 MeV Au ions leads to weaker spatial coupling of electronic and damage energy dissipation processes, and damage production is only weakly affected.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Vibrational dynamics in the undercooled liquid of ultra-fragile metallic glasses

Here, the lattice vibrational dynamics in ultra-fragile Pt 57 Cu 23 P 20 and Pt 60 Cu 20 P 20 glasses were investigated using in situ inelastic neutron scattering to measure the evolution of the phonon density of states (DOS) through the glass transition, into the undercooled liquid phase, and through crystallization. The phonon DOS curves are remarkably similar for the undercooled liquids and the amorphous solids. The vibrational entropy was obtained with these phonon DOS curves, and changes in vibrational entropy caused by the glass transition are less than 0.05 k B /atom. For these highly fragile metal/metalloid glasses, the excess vibrational entropy of the undercooled liquid is about 0.5–1.5 Jmol -1 K -1 , compared to 0.1 Jmol -1 K -1 reported earlier for stronger Cu–Zr glasses. Our work suggests that, consistent with our previous report for less fragile metallic glasses, the fraction of the total excess entropy that is vibrational is small, and the discontinuity in excess entropy at the glass transition is mainly due to configurational entropy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Anomalous neutron nuclear-magnetic interference spectroscopy

The electron–phonon interaction plays a critical role in materials’ electrical, thermal, optical, and superconducting properties. However, measuring the phonon mode-resolved electron–phonon interaction has been challenging. Here we propose neutron-scattering-based Anomalous Neutron nUclear-Magnetic Interference Spectroscopy (ANUBIS), where the co-existence of neutron nuclear scattering and magnetic scattering leads to anomalous dynamical structure factor under the presence of the electron–phonon interaction. Such anomalous structure factor is linear in electron–phonon coupling constant at the phonon wavevector, and is directly proportional to the momentum and energy-resolved dielectric function. The experimental configuration can be achieved using existing polarized inelastic neutron scattering setup, and an order-of-magnitude estimate shows the viability to observe the anomalous scattering signal is around to 10 –4 to 10 –3 relative to phonon scattering, which is examined through the preliminary experiments and achievable at emerging neutron facilities. As a result, our proposal offers an alternative neutron-based metrology to probe the crucial electronic properties.

Electron-phonon interactions↗
DOE OSTI2026

Disorder-induced proximate quantum spin ice phase in Pr2Sn2O7

We report a comprehensive bulk characterization and neutron scattering investigation of single-crystalline Pr2Sn2O7, a magnetic pyrochlore synthesized via a flux-growth method. Unpolarized neutron diffuse scattering reveals the emergence of spin-ice correlations below T ∼ 1 K, evidenced by the development of anisotropic pinch-point features that are consistent with quantum-spin-ice (QSI) behavior. A.C. susceptibility measurements indicate a progressive slowing of spin dynamics in this regime, culminating in complete spin freezing below Tf ≈ 0.15 K. Inelastic neutron scattering at T = 0.5 K reveals a broad spectrum of quasi-elastic magnetic excitations, with intensity in the low-energy range [0, 0.2] meV significantly suppressed below Tf. Meanwhile, an incipient (100)-type magnetic order begins to nucleate, and a gapped excitation centered at ℏω = 0.23 meV persists. We further identify two distinct dynamical timescales above Tf , a slow component τ_slow ∼ 10^−5 s and a fast component τ_fast ∼ 10^−10 s, in quantitative agreement with theoretical predictions for QSI systems. Taken together, these results indicate that Pr2Sn2O7 enters a disorder-induced spin-frozen phase below Tf , lying in close proximity to a U(1) quantum spin liquid

Luo, Yi [ORNL] (ORCID:000900048710428X)↗