Search NASA⌕ Search

SEARCH · Search NASA

Results for “lignin-based jet fuel”

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.

687 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Sustainable Production of Biomass‐Derived Graphite and Graphene Conductive Inks from Biochar

Abstract Graphite is a commonly used raw material across many industries and the demand for high‐quality graphite has been increasing in recent years, especially as a primary component for lithium‐ion batteries. However, graphite production is currently limited by production shortages, uneven geographical distribution, and significant environmental impacts incurred from conventional processing. Here, an efficient method of synthesizing biomass‐derived graphite from biochar is presented as a sustainable alternative to natural and synthetic graphite. The resulting bio‐graphite equals or exceeds quantitative quality metrics of spheroidized natural graphite, achieving a RamanI D /I G ratio of 0.051 and crystallite size parallel to the graphene layers (L a ) of 2.08 µm. This bio‐graphite is directly applied as a raw input to liquid‐phase exfoliation of graphene for the scalable production of conductive inks. The spin‐coated films from the bio‐graphene ink exhibit the highest conductivity among all biomass‐derived graphene or carbon materials, reaching 3.58 ± 0.16 × 10 4 S m −1 . Life cycle assessment demonstrates that this bio‐graphite requires less fossil fuel and produces reduced greenhouse gas emissions compared to incumbent methods for natural, synthesized, and other bio‐derived graphitic materials. This work thus offers a sustainable, locally adaptable solution for producing state‐of‐the‐art graphite that is suitable for bio‐graphene and other high‐value products.

Chemistry↗
DOE OSTI2026

Influence of controlled disorder on the dipolar spin-ice state of Ho-based pyrochlores

Pyrochlore magnets of the form 𝑅 2 ⁢𝐵 2 ⁢O 7 , in which rare-earth ions on the 𝑅 site form a three-dimensional network of corner-sharing tetrahedra, provide a canonical setting for geometrical frustration. Ho-based pyrochlores host a dipolar spin-ice ground state, characterized by Ising moments constrained by the ice rules and elementary excitations analogous to magnetic monopoles. Here, in this work, we examine how controlled chemical disorder influences this state by introducing site mixing on the nonmagnetic 𝐵 site in two compounds. Ho 2⁢ GaSbO 7 contains only Ga 3+ /Sb 5+ charge disorder, whereas Ho 2 ⁢ScSbO 7 exhibits both charge and substantial size disorder arising from the large ionic-radius mismatch between Sc 3+ and Sb 5+ . Although both materials retain the pyrochlore structure, neutron-scattering measurements reveal a reduced correlation length for the 𝑅/𝐵-site cation ordering and enhanced local structural distortions in Ho 2 ⁢ScSbO 7 . Despite these structural differences, bulk thermodynamic measurements and magnetic diffuse scattering demonstrate that both systems exhibit the defining signatures of a dipolar spin-ice state. Low-energy inelastic neutron spectroscopy further uncovers broad magnetic excitations that develop within the dipolar spin-ice regime, a feature absent in pristine Ho pyrochlores and indicative of disorder-induced splitting of the non-Kramers ground-state doublet. Together, these results show that controlled disorder generates tunable transverse-field-driven quantum fluctuations in Ho-based pyrochlores, although the dipolar spin-ice state is remarkably robust to this disorder.

magnetic anisotropy↗
DOE OSTI2026

Anion-dependent phase behavior of methylimidazolium-based ionic liquids mixed with water: Correlation between local molecular structure and mesoscale behaviors

The mesoscopic phase behavior of decylmethylimidazolium (C 10 mim) ionic liquids (ILs) bearing three monovalent anions—thiocyanate (SCN − ), nitrate (NO 3 − ), and chloride (Cl − )—mixed with water at relatively high IL contents (50–95 wt%) was investigated. Small-angle and wide-angle X-ray scattering (SAXS/WAXS) were employed to follow the evolution of both local and mesoscale structures across this composition range. In the absence of water, C 10 mimSCN and C 10 mimCl behaved as disordered liquids, whereas C 10 mimNO 3 spontaneously formed a hexagonally ordered cylindrical mesophase and displayed a sticky-solid macroscopic appearance. Upon addition of water, C 10 mimSCN remained a viscous liquid and only weakly ordered lamellar domains were observed. This limited ordering is attributed to the lack of hydrogen bonding and weak interaction energy between SCN − anions and C 10 mim + cations. In contrast, the trigonal-planar NO 3 − and point-like Cl − anions promoted the formation of well-defined hexagonal mesophases up to 35–45 wt% water. The formation of hydrogen bonding of the two anions with imidazolium ring protons likely enabled the creation of compact ion clusters that effectively exclude water molecules from the immediate vicinity of the IL aggregates. These findings demonstrate that the interaction energy between ion pairs dictate IL-water interactions and therefore control the transition from disordered liquids to ordered mesophases in IL/water mixtures with high IL contents. In conclusion, the combined SAXS/WAXS analysis reveals a correlation between local intermolecular structure and the emergence of mesoscopic order, providing a systematic framework for tailoring mesoscale structures in alkylimidazolium-based IL/water systems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Sensitivity-driven Wide Area Protection (SWAP) Coordination Tool for High Penetration of Inverter-based Resources (IBR)

Traditionally, power system generation sources have been composed of synchronous generators, of which the fault current behavior is understood with minimal differences between generation size and types due to the physics of their construction. Present protection schemes and modeling methods are based upon these understood characteristics. Most renewable generation is composed of inverter-based resources (IBR), in which fault current is determined by switching control software and hardware limitations, each of which can vary between manufacturers and even between models of the same manufacturer. The resulting fault current is low in magnitude, low in negative-sequence current, unpredictable phase angles, and is a challenge to model. These characteristics also result in a challenge to traditional protection schemes and fault simulation software. To address several of these concerns, the project has the following goals: 1. Improve IBR models: Improve IBR models used in short circuit (SC) programs to accurately capture the response of IBRs at the bulk power system (BPS) level for fault and protection studies. 2. Develop automation tool: Develop an automation tool that allows engineers to identify protection coordination and sensitivity issues by performing SC and protection coordination studies in a high IBR-penetrated grid by applying variations to the IBR models, faults, contingencies, etc. 3. Develop schemes: Develop new protection mitigation solution schemes that complement the existing protection systems to ensure safe operation of the BPS with higher IBR penetration levels. The project team did not achieve this final goal, as the Department of Energy (DOE) stopped the project early due to changes in DOE funding priorities. The termination notice came at the beginning of the final project phase, while the team was identifying and beginning to investigate protection issues. It should be noted that the team discussed a 100% penetration scenario. However, this scenario would require the use of grid-forming IBR models that are not presently available. Since developing these models requires additional effort, the 100% penetration scenario was not pursued during this project. In the future, developing the methodology and models for the 100% scenario could benefit the industry.

14 SOLAR ENERGY↗
DOE OSTI2026

Universal Relationship between Limiting Current and Electrochemical Transport Properties in Malonate-Based Polymer Electrolytes

There is considerable interest in developing high-performance electrolytes for rechargeable lithium batteries. For practical applications, the electrolyte must support large dc currents. However, the parameters most often reported in the literature, conductivity, κ, and current fraction, ρ + , reflect ion transport in the limit of infinitesimal currents. In this limit, the efficacy of an electrolyte is given by the product κρ + . The limiting current density, i lim , is the maximum current density that can be applied across an electrolyte; the cell voltage diverges if the applied current density exceeds ilim. This parameter reflects ion transport in the limit of large dc currents and is therefore of practical interest. It would therefore be convenient if i lim could be predicted from measurements of κρ + . In order to explore this possibility, we studied six malonate-based polymers and PEO at a fixed salt concentration (r = 0.08) and temperature (90°C) using symmetric cells with planar electrodes. Unfortunately, there is no correlation between ilim and κρ + . When the applied current density, i, is less than ilim, the cell voltage approaches a stable plateau, ϕ plateau . Here, we found a linear dependence between i and thickness-normalized plateau potential, ϕ plateau L –1 , irrespective of the magnitude of the applied current. In all seven polymer electrolytes, we found a linear correlation between ilim and the slopes of these lines, σ. In other words, measurements of σ can be used to predict the limiting current.

Jana, Rounak [Lawrence Berkeley National Laborator↗
DOE OSTI2026

High-Q superconducting lumped-element resonators for low-mass axion searches

Low-frequency superconducting lumped-element resonators have recently attracted significant attention in the context of axion dark matter searches. Here, we present the design and implementation of a fixed-frequency superconducting resonator operating near 250 kHz, possessing an inductor volume of ∼1 liter and achieving an unloaded quality factor Q ≈ 2.1 × 10 6 . This resonator represents a significant improvement over the state of the art and informs the design of searches for low-mass axions.

Bosons↗
DOE OSTI2026

Cave-based chronology of valley incision in the glaciated Northern Carpathians

While glaciers leave a mark on mountain landscapes, the importance of glacial erosion in mountains at a global scale remains insufficiently constrained despite over a century of research aimed at quantifying the rates and magnitudes of such erosion. Measuring rates of landscape change across the Quaternary in repeatedly glaciated mountain landscapes requires the ability to accurately partition the relative role of glacial versus fluvial processes in forming relief; however, the modification of landforms by surface erosion hinders quantification of process rates over the last few million years. This study presents 32 paired 26 Al and 10 Be cosmogenic nuclide burial ages from allochthonous, quartz-bearing sediment deposited within a paleophreatic karst system in the northern Carpathians. These data quantify Pliocene-to-present timing of valley incision and karst system reactivation, as well as erosion rates over time. Post-burial, muon production-correction yields 24 fully resolved ages, six minimum ages, and two unresolvable ages. Burial ages are interpreted to document the tempo of valley incision via phreatic passage formation and reveal that ∼80% of current valley relief (∼700 m) had developed by approximately 1.6 Ma. Paleoerosion rates derived from this dataset suggest little change in the processes generating sediment washed into the karst networks since the Pliocene. These results suggest that river incision dominated relief creation in Carpathian valleys that bear the imprint of glaciation. Our findings suggest the importance of investigating other mid-latitude mountainous regions worldwide to determine if they have similar erosion and topographic histories.

Carpathians↗
DOE OSTI2024

A butterfly-shaped acceptor with rigid skeleton and unique assembly enables both efficient organic photovoltaics and high-speed organic photodetectors

ABSTRACT It remains challenging to design efficient bifunctional semiconductor materials in organic photovoltaic and photodetector devices. Here, we report a butterfly-shaped molecule, named WD-6, which exhibits low energy disorder and small reorganization energy due to its enhanced molecular rigidity and unique assembly with strong intermolecular interaction. The binary photovoltaic device based on PM6:WD-6 achieved an efficiency of 18.41%. Notably, an efficiency of 19.42% was achieved for the ternary device based on PM6:BTP-eC9:WD-6. Moreover, the photodetection device based on WD-6 demonstrated an ultrafast response speed (205 ns response time at λ of 820 nm) and a high cutoff frequency of −3 dB (2.45 MHz), surpassing the values of most commercial Si photodiodes. Based on these findings, we showcased an application of the WD-6-based photodetection device in high-speed optical communication. These results offer valuable insights into the design of organic semiconductor materials capable of simultaneously exhibiting high photovoltaic and photodetective performance.

Science & Technology - Other Topics↗
DOE OSTI2025

Photosynthetic capacity is reduced by warming but unaffected by elevated CO2 in seedlings of five boreal tree species

Abstract Increasing atmospheric CO2 concentrations fuel global warming, with boreal regions warming at a faster rate than many other areas. Boreal forests are an important component of the global carbon cycle, yet we have little data on photosynthetic responses of boreal trees to elevated CO2 (EC) and warming. We grew seedlings of 5 widespread North American boreal tree species (from Betula, Larix, Picea, and Pinus) under current (410 ppm) or elevated (750 ppm) CO2 and either ambient (+0 °C) or increased (+4 °C or +8 °C) temperature, then measured photosynthetic traits over a range of leaf temperatures. Our results were generally consistent across species: photosynthetic capacity (maximum rates of Rubisco carboxylation, Vcmax, and electron transport, Jmax) was unaffected by EC but decreased under +8 °C warming. Accordingly, net photosynthesis measured at the growth CO2 concentration (Agrowth) was reduced under warming and increased under EC. The thermal optimum for Agrowth (ToptA) increased by ∼1.8 °C with EC but increased with warming in only two species. In contrast, the activation energies and thermal optima for Vcmax and Jmax, which are used to estimate photosynthesis in Earth System Models, were unaffected by growth environment. There were a few interactions between growth, CO2, and warming. These results suggest increased photosynthesis of widespread boreal tree species under EC may be offset by future reductions in photosynthetic capacity related to warming. We also show that the temperature sensitivities of parameters used to estimate global photosynthesis in large-scale models are generally unaffected by simulated climate change in these species.

Plant Sciences↗
DOE OSTI2026

Electron/hadron separation with the TRD based on GaAs pixel detectors

Transition Radiation Detectors (TRDs) are widely used for particle identification in both high-energy physics and astroparticle physics. They are typically equipped with gaseous detectors. The main limitation of these types of detectors is that TR photons and ionization losses cannot be decoupled, which significantly reduces the particle separation power. Recent advancements in the development of pixel detectors based on GaAs sensors offer a unique opportunity to effectively detect TR photons and separate them from ionization losses. Such detectors represent novel devices that combine precise tracking capabilities with particle identification (PID) properties. The present work is dedicated to an experimental study of particle identification properties of the TRD prototype based on 500 μ m-thick GaAs sensor bonded to a Timepix3 chip. Studies were performed at the CERN SPS and it was shown that at a particle momentum of 20 GeV/c, the probability of misidentifying a hadron as an electron is below 10 −2 for an electron detection efficiency of 98%–99%, and it is 1 . 6⋅10 −4 for electron detection efficiency of 90%. This performance surpasses the electron/hadron rejection power of all known TRDs by an order of magnitude for the same detector length.

GaAs detectors↗
DOE OSTI2025

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗
DOE OSTI2025

Generative learning of densities on manifolds

A generative modeling framework is proposed that combines diffusion models and manifold learning to efficiently sample data densities on manifolds. The approach utilizes Diffusion Maps to uncover possible low-dimensional underlying (latent) spaces in the high-dimensional data (ambient) space. Two approaches for sampling from the latent data density are described. The first is a score-based diffusion model, which is trained to map a standard normal distribution to the latent data distribution using a neural network. The second one involves solving an Itô stochastic differential equation in the latent space. Additional realizations of the data are generated by lifting the samples back to the ambient space using Double Diffusion Maps , a recently introduced technique typically employed in studying dynamical system reduction; here the focus lies in sampling densities rather than system dynamics. The proposed approaches enable sampling high dimensional data densities restricted to low-dimensional, a priori unknown manifolds. The efficacy of the proposed framework is demonstrated through a benchmark problem and a material with multiscale structure.

Double diffusion maps↗
DOE OSTI2025

Gas-mediated defect engineering in earth-abundant Mn-rich layered oxides for non-aqueous sodium-based batteries

Gases are often by-products of battery materials during cell formation and degradation, affecting the cycle life and safety of rechargeable batteries. However, understanding gas-mediated (electro)-chemical reactions and nanoscale structural transformations during the synthesis of battery electrode materials remains challenging because of the lack of suitable characterization routes and the complexity of the interplay between thermodynamics and kinetics. Here, in this study, we use operando synchrotron X-ray diffraction, in situ transmission X-ray microscopy and multiscale modelling to elucidate the reaction pathways and microstructural defect development of earth-abundant Mn-rich layered oxides as positive electrode materials for sodium-based batteries. In particular, we demonstrate the dominant role of CO 2 over O 2 and H 2 O (g) in modulating the competition between entropy and enthalpy during solid-state synthesis. Using Ni 0.25 Mn 0.75 CO 3 as a model precursor, we reveal that CO 2 generation favours the formation of entropy-driven metastable intermediates, suppresses closed pore/nanovoids formation and decreases chemical heterogeneity and residual lattice strain of Mn-rich layered oxides during the synthesis. This result motivates a fast-sintering strategy to promote CO 2 release, which ultimately leads to improved chemo-mechanical and electrochemical stability of the Mn-rich positive electrodes when tested in non-aqueous Na metal coin cells.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Influence of functional additives, fillers, and pigments on thermal and catalytic pyrolysis of polyethylene for waste plastic upcycling

Pyrolysis offers a relatively green and economical method to convert waste plastics into valuable chemicals and fuels without the need for harmful solvents, toxic chemicals, or costly high-pressure reactors. Despite its popularity among chemical upcycling technologies, industrial adoption suffers from feedstock heterogeneity, low-quality products, and catalyst deactivation. Most plastics in our daily lives are formulated with functional additives, fillers, and colorants. These additives remaining in end-of-life waste streams increase feedstock heterogeneity, creating a challenging issue in recycling plastics. Still, the potential impacts of additives on the chemical upcycling of plastics have been poorly understood. In this study, polyethylene compounded with a range of widely used additives (antioxidants, stabilizers, pigments, fillers, slip agents, and flame retardants) was subjected to both thermal pyrolysis and catalytic pyrolysis in different catalyst-to-feedstock contact modes. It showed that many inorganic additives, such as talc, kaolin, CaCO 3 , TiO 2 , carbon black, and zinc stearate, facilitated polymer decomposition during pyrolysis, increasing light hydrocarbons while also promoting aromatic and carbon residue formation. Conversely, antioxidants and stabilizers inhibited depolymerization, favoring heavier hydrocarbons. During catalytic pyrolysis with HZSM-5 zeolite, additives strongly enhanced aromatic and catalytic coke formation, especially when there was direct contact between plastics and catalysts. Although certain additives seem beneficial in the short term by promoting polymer cracking and improving the selectivity of aromatics, the transport of the additives and their degradation products and increased carbon coking can contaminate products, deactivate or modify catalysts, and foul reactors. These findings address a critical knowledge gap in effectively converting waste plastics via a greener route.

42 ENGINEERING↗
DOE OSTI2026

Detecting the full photoemission cone from laser-based ARPES experiments by leveraging deflector technology

Angle-resolved photoemission spectroscopy (ARPES) provides a direct access to the electronic band structure of solid and molecular systems. The momentum range accessible by this technique depends directly on the photon energy used, and low-photon-energy sources are insufficient to photoemit electrons over the full Brillouin zone of most quantum materials. In addition, while electrons are emitted over a 2π solid angle, conventional hemispherical analyzers only collect a small subset of those electrons. A previous work [Gauthier et al., Rev. Sci. Instrum. 92, 123907 (2021)] demonstrated that electrons emitted over a larger field-of-view can be acquired in one fixed configuration by accelerating them toward the analyzer with a bias voltage. Here, in this study, we extend this work by leveraging the deflector technology of novel ARPES hemispherical analyzers. We demonstrate the ability to detect all 2π photoemitted electrons in a fixed configuration for various materials, such as gold, cuprates, and transition-metal dichalcogenides. This approach is especially advantageous for time-resolved ARPES, as electron dynamics over a large momentum range can be accessed with identical measurement conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Mechanistic origin of solvent-dependent thermal stability in sodiated Sn anodes for sodium-ion batteries

Understanding the thermal stability of high-energy density alloy anodes is critical for the safe deployment of sodium-ion batteries (SIBs). Here, accelerating rate calorimetry (ARC), post-mortem characterizations, and density functional theory (DFT) calculations are combined to understand the thermal reactivity of fully sodiated Sn, Sn-hard carbon (HC) blends, and HC anodes in carbonate- and ether-based electrolytes. ARC measurements show that propylene carbonate (PC) causes earlier self-heating rate (SHR) onset and higher reactivity than tetraethylene glycol dimethyl ether (TEGDME), indicating inferior thermal stability. Sodiated Sn exhibits better thermal stability than sodiated HC, while Sn-HC blends show intermediate behavior that improves with increasing Sn content. Post-ARC analyses reveal desodiation of Na15Sn4 to metallic Sn with particle coalescence, whereas Sn-HC blends and HC retain finer morphologies. PC promotes Sn oxidation to SnO, while TEGDME suppresses oxide formation; NaPF6-containing electrolytes additionally form NaF. DFT calculations show that PC adsorption lowers Na extraction energy and enhances interfacial electronic interactions, facilitating Na release and reductive decomposition. These results establish a direct correlation between solvent-dependent reaction pathways and thermal stability in SIB alloy anodes.

Accelerating rate calorimetry↗
DOE OSTI2026

Quantifying System Strength From Grid-Forming Resources Using Frequency Scan Approach: Preprint

Current industry practices for quantifying the system strength contribution from grid-forming (GFM) resources to ensure stability of power systems dominated by inverter-based resources (IBRs) are primarily based on iterative electromagnetic transient (EMT) time-domain simulation studies. While feasible, these approaches are resource-intensive, lack scalability and intuition, and might not evaluate the system strength contribution over the entire frequency range of interest. This paper introduces a novel, frequency-domain approach to quantify system strength support provided by a GFM resource using frequency scans. The proposed method uses transfer functions from the grid voltage magnitude (V) and phase (?), respectively, to the reactive (Q) and active power (P) output of a GFM resource for quantifying its contribution to system strength. These transfer functions provide a direct measure of the ability of a GFM resource to behave as a stiff voltage source behind a reactance over a specified frequency range, enabling robust quantification of its system strength contribution. The key innovation of this work is the development of a frequency domain system strength metric called the dynamic short-circuit ratio (dSCR) that is suitable for IBR-dominated power systems and is directly related with the familiar short circuit ratio (SCR) metric. The new metric, dSCR, enables the assessment of system strength contributions from both synchronous machines and converter-based GFM resources using a unified benchmark, which is not possible with the traditional SCR metric. The paper also demonstrates how impedance scans could identify if an unstable condition observed during weak grid conditions is a result of the lack active or reactive power support or both. By leveraging the proposed frequency-domain dSCR metric for quantifying system strength contribution from GFM IBRs, the paper demonstrates targeted mitigation strategies for weak grid instabilities without resorting to repeated, time-consuming time-domain simulations. The result is a scalable and efficient approach to remediate stability challenges in power systems with high shares of IBRs and accelerating the integration of GFM technologies for system strength support in power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue↗